有限体积方在守恒律中的应用.pdfVIP

中英文摘要 §o．1 中文摘要 来自自然科学与工程领域中的大多数微分方程在数学上表现为守恒形式 J二0 r可，，(M z∈R“．t∈R+三(0，。c)，¨∈R“ “J．11 z∈R“ 0j--=}l【}(r) 它是自然界中的守恒定律在数学上的直接反映，对流体力学方程组而言，它就是质量，动 1 量，能量守恒得到的方程。由于双曲守恒律(o1)没有其它项，如色散(dispersion)，扩 流体的流动引起的输运)时，守恒律(1】11)的解失去光滑性(这里不特殊说明守恒律就 指该意义下)，甚至即使光滑的初始数据，解随着时间的发展会变成不连续，这在物理上 表现为激波的形成。从流体力学的角度上看，(o1．1)事实上就是粘性很小的近似。当考 1 虑粘性后，即在数学上反映为(01)中多了扩散项(二阶导数项)．即使很粗糙的初始数 据，解在瞬间内变的很光滑．这由于流体的粘性扩散引起，这种对流一扩散问题可用古 典的微分方程来研究。自然的想法就是当粘性趋于零时，带粘性的对流一扩散问题的解 I 在某意义下趋于无粘性问题(01)的解，这就是正则化方法。另一办法从离散(数值)角 1 度上研究仅有对流项的守恒律(O1)．如构造它的差分格式，甚至更一般的有限体积格 式，有限元及谱方法等，从这些格式构造近似解(常表现为分片多项式)来逼近原守恒律 的解。正则化方法与数值角度研究都在于估计解的构造，其次考虑估计解的稳定性，即 在适当度量下估计解大小的上界，而且着重研究与正剐化系数【如粘性大小)与数值的离 散参数(如时间，空间步长大小)无关的上界。在利用紧框架(如BV紧，Ll紧，补偿紧 等)得到估计解有子列在某种模意义下收敛到某一个函数，再利用估计格式(估计解所满 足的方程)逼近守恒律来证明所得到的函数就是在某种意义下守恒律的解．由于大时间 1 范围内守恒律(01)的解表现为很差的正则性，它不能在古典意义下定义，即在每一点 下的导数无意义．使得古典办法研究遇到很大困难，它只能在弱意义下定义弱解，但往 往这种弱解不唯一，需要某条件限制确保解的唯一性，在数学上称为熵条件，满足该条 件的弱解称为熵解。从流体力学来看，它事实上是热力学第二定理的反映，即熵越过激 波f一种间断1要增加．各种估计格式构造的估计解应反映这一事宴，即满足熵不等式． 从实际计算来看，总是通过离散化求解，不考虑计算的积累误差，它的稳定性与计算精 度都依赖与真解的光滑性，一般说，在解较光滑的区域有较好的稳定性与计算精度，而 l 1)的求解来说，一般保持一定精度的数值格式在 在较粗糙的区域则相反．对守恒律(0 光滑性较差的区域计算效果不大好，首先，捕捉解的光滑性较差的区域(间断处)较差， 其次，在该区域附近逼近不大好，如不能模拟物理上的震荡。所有这些都反映了解的正 则性的缺乏给数值研究带来困难。 varial?le lj1,Izkov[48】在7()年代用双变量技巧(Doubletechnique)解决了多维单个双 l 曲守恒律(01)的熵解的适定性问题，即熵解的存在性，唯一性及正规性结果，特别熵解 物正规(粘性方法)或用f}“上的矩形网格构造差分格式，并得到了相应的￡1收敛速度估 I 与P(,rthame{51推广到用非线形退化守恒的抛物问题逼近双曲守恒律(0I)，并得到两 21把双变量技巧得到非线形 者熵解的误差的显式表示，1999年Cockimrn与Gripenbeg[I 退化守恒的抛物问题半群解在L1意义