一类非线性扩散方程的间断有限元方法研究

来源 :中国工程物理研究院 | 被引量 : 10次 | 上传用户:ningmengpan
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
非线性扩散方程的求解是流体力学、传热学、生物学、化学和斑图动力学等领域的一个热门研究课题。由于非线性扩散问题的复杂性,大多数定解问题还是不能精确求解的,因此对于非线性扩散方程数值方法的研究就显得非常重要。间断Galerkin(DG)有限元方法具有局部守恒,容易实现自适应和并行化的特点,因此在众多领域得到了广泛地应用。本文围绕间断Galerkin方法展开讨论,我们发展了新型的DG格式,给出相应的收敛性和稳定性分析,并应用DG方法对一类非线性扩散方程进行了研究。这些方程包括:Burgers类方程,N-载体系统,带有间断系数的扩散方程,非线性Schrodinger方程,反应扩散方程以及辐射扩散方程。本文的主要工作有以下几点:首先在第一章,我们对求解扩散方程的各种DG方法进行了回顾,这些方法包括:局部间断Galerkin方法(LDG), Baumann-Oden DG方法,Gassner广义黎曼问题DG方法,直接间断Galerkin方法(DDG), Cheng-Shu提出的间断Galerkin方法和内惩罚类方法。我们分析这些DG方法之间的联系,并比较了各自特点。在时间离散方面,我们着重介绍保持强稳定性质的Runge-Kutta方法,指数时间差分方法和隐积分因子方法,其中隐积分因子方法既可以避免显式方法苛刻的时间步长限制,又可以保证DG方法在每个单元上局部求解的特点,适用于DG格式离散后得到的半离散形式。第二章应用局部间断Galerkin方法求解了一类Burgers方程和N-载体系统。数值算例得到的结果与精确解和其他结果进行了比较,验证了该方法的高精度和可靠性。第三章应用间断Galerkin方法求解带有间断系数的扩散方程。对于带有间断系数的抛物方程,针对扩散系数间断的特点,我们提出了一种满足能流连续的加权数值流量,发展了一种新型的DG格式,并证明了该格式的稳定性和收敛性,数值算例验证了收敛性分析。对于带有间断系数的椭圆方程,我们构造一种新的加权内惩罚方法,证明了双线性形式的连续性和强制性,并给出收敛性证明,数值算例表明我们的DG方法对于求解强间断系数问题十分有效。第四章讨论非线性Schrodinger方程(NLS)的间断Galerkin方法。我们应用直接间断Galerkin方法求解了一维、二维以及耦合的非线性Schrodinger方程,证明了直接间断Galerkin方法得到的数值解满足质量守恒律。与徐岩等人的局部间断Galerkin方法相比,直接间断Galerkin方法不用引入辅助变量,从而减少了计算量。通过对众多数值算例的求解,我们得到了非常好的数值模拟效果。第五章讨论反应-扩散方程的间断Galerkin方法。空间离散采用直接间断Galerkin方法,时间离散应用隐式积分因子方法求解。隐式积分因子方法在得到高精度的同时又具有良好的稳定性,避免了显式时间离散方法对时间步长的苛刻限制。更为重要的是该方法无需求解大规模的非线性代数方程组,可以应用Newton方法逐个单元求解非线性代数方程组,保持了DG方法局部求解的特性。第六章将隐式积分因子方法与间断Galerkin方法结合求解辐射扩散方程。非平衡辐射扩散方程在不同的物质界面处扩散系数可以相差量级,我们将第三章构造的加权数值流量推广至此,得到一种新的加权DG方法。对于一维非平衡辐射扩散方程,由于自由度数不是很大,我们应用全隐格式进行时间离散;对于二维非平衡辐射扩散方程,我们首先利用线性化对非线性扩散系数进行显式处理,然后再利用积分因子方法求解非线性常微分方程组,得到一种新的隐-显积分因子DG方法。本文首次将DG方法应用到强耦合、多介质、强非线性的辐射扩散方程组,这是间断Galerkin方法一个新的应用拓展。第七章对全文做了总结,并对下一步的工作进行了展望。
其他文献
一、经济全球化的正面与负面都影响世界全球化已经成为当代世界的发展潮流和客观趋势。正确认识经济全球化和政治多极化,有助于我们在融入世界的同时,更好地维护国家安全。党的
目的分析慢阻肺患者应用自拟慢阻肺汤合康复训练的临床疗效。方法选取86例慢阻肺患者,经单双病号分组,2组各43例,常规组应用康复训练治疗,研究组应用自拟慢阻肺汤合康复训练
目的分析医院肾内科辅助用药超说明书用药的合理性,为临床合理用药提供科学依据,并进一步规范医院辅助用药超说明书用药行为。方法通过HIS系统查询医院肾内科2015年1月-12月
<正> 我们的国家,是以工人阶级为领导的、以工农联盟为基础的人民民主专政的、国家。对人民内部的民主和对反动派的专政,这两方面互相结合起来,就是人民民主专政。这个专政是
针对目前图像融合算法基于像素级实现、未充分考虑图像纹理特征、融合效果不理想的现状,本文通过结构张量特征值构造特征模板,在梯度域通过特征加权获得融合梯度场,充分考虑
在早期阅读倍受关注的今天,儿童读物的选择成为很多家长和老师非常关心的问题。市面上涌现出大量的儿童读物,其内容、质量、价格等方面差异较大,导致儿童读物的选择存在很多
<正>能随着光线调整的百叶窗,因其所创造的美式休闲风格而备受世人喜爱,同时也让空间的视觉感受随光影挪移更显丰富。百叶窗虽小,但功能多样:夏天遮阳隔热、冬天保温隔冷、平
据英国媒体2018年11月28日报道,英国巴克莱银行(Barclays)与英国连锁咖啡品牌Costa Coffee合作推出了一款可以重复使用的咖啡杯Clever Cup,该咖啡杯底部装有非接触式芯片,顾
本文主要研究对流扩散问题非协调有限元逼近的残量型后验误差估计.针对一系列的稳定化有限元方法,我们在统一的框架下推导了半健壮的和健壮的后验误差估计,并将此理论结果推
西方发达国家凭借其在经济和信息全球化中的主导地位,不断对中国进行文化扩张和渗透,使得我国主流意识形态和民族文化受到较大冲击,面临文化生态保护形势严峻、文化创新能力