异质偏微分方程的区域分解和深度学习算法研究

来源 :东北师范大学 | 被引量 : 0次 | 上传用户:hbbzy
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
偏微分方程是科学家对自然物理现象建模的一类重要手段,例如Poisson方程用于描述势场与电荷或质量密度分布的关系;Navier-Stokes方程用于描述粘性流体物质运动过程中速度以及压强的变化;Helmholtz方程用于描述物理波的传播现象等等.但是在现实世界中存在着大量的、复杂的物理现象无法通过单一偏微分方程进行刻画.当物理区域中不同位置上的介质具有不同的特性时,常用异质偏微分方程来描述异质区域中相关物理量的状态和变化,此时方程具体表现为系数函数是不连续的.异质偏微分方程的数值求解通常比求解单一方程困难许多,直接应用经典的数值方法往往会遇到一些棘手的问题.传统的有限元方法或有限差分法离散异质偏微分方程得到的线性代数系统通常是病态的,另外基于物理定律的界面约束条件需要被特殊地处理以保证算法收敛.本文关心异质偏微分方程的数值求解问题,以流体-流体耦合模型,异质Poisson方程,异质Helmholtz方程为研究对象,以区域分解方法为主要解耦和算法工具,提出求解模型问题的稳定、高效、可扩展性强的数值算法.流体-流体耦合模型描述的是在关心的物理区域上两种不相容的流体相互作用的物理过程,数学上使用具有摩擦型界面条件的耦合Navier-Stokes方程进行刻画.本文通过对流体-流体耦合模型进行线性化和时间离散得到Stokes-Stokes耦合方程,进一步提出求解该系统的Schwarz交替算法.我们通过能量分析法得到Schwarz交替算法对于任意大于0的牵引系数(traction coeffcient)都是收敛的.进一步借助Fourier分析,我们首次得到算法的收敛因子分别与流体粘性系数、牵引系数、时间步长参数的关系.特别地,理论分析证明了算法收敛因子与空间离散步长无关.为了得到收敛更快的区域分解算法,本文提出针对于模型问题的优化Schwarz方法.不同于使用基于物理界面传输条件的Schwarz交替算法,优化Schwarz方法构造极限状态下满足物理定律的人工界面传输条件,基于优化的方法得到最优的传输参数.具体而言,本文利用Fourier分析得知最优界面算子是伪微分算子,其对应的界面传输条件在数值计算中难以应用,因此构造近似的Robin型和双边Robin型界面传输条件.我们通过求解最小最大问题得到两类近似界面条件的最优传输参数,并在空间步长很小的渐近意义下分析最优传输参数和算法收敛因子.一系列的数值实验说明了提出的算法是稳定且高效的,同时验证了理论的分析结果.本文首次提出深度区域分解方法(deep domain decomposition method),其联结深度学习和区域分解的思想,可求解定义在复杂区域上的异质偏微分方程.基于网格离散的传统区域分解方法在求解复杂区域上的偏微分方程时,计算精度会受到网格离散以直代曲误差的影响.而深度学习在求解偏微分方程问题上是无网格的,因此可适用于复杂区域上的模型问题.但深度学习无法直接处理异质偏微分方程中界面条件的异性性质和解的不连续性,不可直接用于求解异质偏微分方程.深度区域分解方法继承了深度学习和区域分解方法的优点,具有无网格、可并行计算的特性.该方法首先将异质偏微分方程根据区域性质分解成多个子问题,然后使用深度学习方法分别求解子问题,最后基于界面传输条件迭代更新数值解.本文在Tensor Flow框架下使用随机梯度下降法训练神经网络参数,得到以神经网络为模型的偏微分方程数值解.在使用训练数据优化神经网络的损失函数时,可以将大量的内部训练点分组,并分别将每一组数据和边界训练点合并成多个批次的训练集.事实上,深度区域分解算法在求解不同的方程时,我们根据方程的性质选择界面传输条件和神经网络模型.在求解异质Poisson方程的数值实验中,发现深度区域分解方法继承了传统区域分解方法的一些性质.算法达到收敛的所需要的迭代次数只与子区域的个数和相邻区域的重叠量有关,与函数空间大小无关.在多种实验设定下,深度区域分解算法的数值收敛率与传统区域分解算法理论收敛率相吻合.面对具有曲界面异质Helmholtz方程的数值求解问题,本文采用激活函数为平面波基函数的复值神经网络作为模型,使得在神经网络层数少和(或)在内部区域训练点少的情况下仍然能达到相对误差为10-~3的精度.异性介质上齐次和非齐次的Helmholtz方程在多种情况下的数值实验结果验证了深度区域分解算法保留了传统区域分解的性质的同时,还可以扩展应用到更复杂的偏微分方程模型问题.本文的研究丰富了Schwarz交替算法和优化Schwarz方法在流体-流体耦合模型中的理论,我们首次提出深度区域分解算法以探索深度学习和区域分解方法在数值求解异质偏微分方程问题上的更多可能性.在应用层面,面对需要数值模拟流体-流体耦合系统的工业应用中,本文提供了高效且稳定的优化Schwarz算法.而面对异质Helmholtz方程或是复杂区域上异质偏微分方程的数值求解问题时,本文提出的深度区域分解方法是一类可行且富有竞争力的算法.
其他文献
本文以德国“通学派”(Sachphilologie)古典学者A.伯克(August Boeckh,1785-1856)与K.O.穆勒(Karl Otfried Müller,1797-1840)通信为研究对象,通过分析书信中两位学者围绕古典学范畴内涉及历史、地理、铭文、考古、神话、建筑、钱币等诸多分支领域的具体问题展开的学术交流与讨论,窥见引领德国古典学走向辉煌的“百科全书式”研究范式具体如何在1
学位
本论文主要研究了李超代数的理想、导子和导子对.主要内容包括理想、c-理想、拟理想、导子、型心、上同调、形变和中心扩张.具体内容分为以下六章:第一章的内容是引言.本章介绍了研究背景与发展现状以及我们所做的工作和本论文结构.第二章的内容是李超代数的三类理想.首先,定义了李超代数的完美理想和导序列,通过极大完美理想和导序列长度刻画了李超代数的可解性.其次,定义了李超代数的几乎完美理想和降中心序列,利用极
学位
在临床实践中,不同的患者对同一种治疗的反应存在很大差异,这使得生物医学逐步从传统的一刀切治疗转变为个性化医疗的新模式.所谓个性化医疗,其核心目标是根据每个个体的独特特征开发个性化药物并确定最优的治疗方案.在许多临床试验中收集到的患者疾病相关数据都是高维的,因此研究者们希望从中筛选出对不同治疗方案有显著影响的因素.此外,慢性病的治疗通常是长期的,多阶段的,而现有的方法大多数考虑的是单阶段高维数据下的
学位
场景能够基于价值主导逻辑演变与企业动态能力进阶交互赋能商业模式创新,但其机理不清、路径不明。综合运用案例研究与程序化扎根分析方法,发掘和提炼场景基于价值主导逻辑演变与企业动态能力进阶交互协同演化赋能商业模式创新的机理和路径。研究发现:(1)企业从“产品功能—服务效用—场景体验”的价值主导逻辑演变出发,通过“感知能力—整合能力—学习能力—创新能力”的企业动态能力进阶,识别市场机遇把握、个性定制驱动和
期刊
传统的煤泥浮选只出精煤、尾煤两种产品,往往存在保证了精煤灰分,尾煤灰分可能偏低的问题。某选煤厂煤泥灰分为55%,采用汽车运输销售,随着环保要求越来越严格,煤泥的运输销售已严重影响该厂正常生产。经过增加二次浮选系统,优化煤泥水处理系统等技术改造后,从煤泥中回收部分中煤,将煤泥灰分提高至75%以上,与水洗矸石一起进行综合利用,彻底解决煤泥滞销问题的同时,提高了经济效益。
期刊
由于时滞现象广泛存在于客观世界,泛函微分方程作为数学模型被应用于物理学、生物学、医学、自动化控制等领域.自然界中的系统不可避免地受到来自系统内部或外部的随机干扰,因此利用随机泛函微分方程建模更加现实有意义.由于解是取值于无穷维空间的随机过程,几乎很难得到其闭形式,因此构造可靠且易于操作的数值格式在理论研究和实际应用中具有重要意义.本文主要研究受布朗运动驱动的非线性随机泛函微分方程的显式逼近理论,包
学位
本篇论文研究了几类具有奇异性质的积分/微分方程的反馈控制问题.对于一类有限维和一类无限维奇异线性Volterra积分方程,分别讨论了它们的线性二次最优控制问题(LQ问题),给出了最优控制的因果状态反馈表示.对于一类常微分方程,考虑了其LQ理论在反馈爆破能控性问题中的应用.方程的解在有限时刻发生爆破也是方程所具有的一种奇异性质.本篇论文的内容分为以下三部分.本篇论文的第一章介绍了反馈控制问题的背景和
学位
<正>消费是联结国内大循环和国际国内双循环的核心,也是拉动经济增长的驱动力。随着北京国际消费中心城市建设有序推进,优质消费资源加快集聚,消费供给持续优化,符合居民消费升级需求的多样化、个性化商品和服务供给增加,这些都有助于释放消费潜力。受各类因素影响,当前消费市场既出现了数字化、智能化、信息化等新动态,也面临着消费能力受挫、消费信心不足等局限。为完善北京消费环境,全力促进国际消费中心城市建设,建议
期刊
回复运动是动力系统的重要研究课题,在研究动力系统稳定性讨论中占有重要地位.回复性的相关概念非常之多,如:周期性、概周期性、几乎自守性、Birkhoff回复性、Poisson稳定性、仿射周期等等.关于微分方程周期解存在性的最早结果可以追溯到Ⅰ.Newton对于开普勒第一定律的证明.自微分方程定性理论建立之后,随着周期解研究重要性的提升以及对各类周期现象的深入研究,人们发现了许多近似周期的自然规律.现
学位
报纸