偏微分方程的数值离散方法.ppt

插值*当积分点的函数是线性插值时Secondorder插值*QUICK(quadraticupwindinterpolationforconvectivekinematics)插值三阶精度，但积分（差分）往往只有二阶精度。插值*高精度：N阶精度的quadrture需要N-1阶多项式插值公式。界面上导数可以用插值公式的微分求出。030102用边界条件替代面积分有限体积法的边界条件*入口：通常给定对流通量(mass,momentum,energy,etc.)壁面和对称面：通量为零边界上函数值给定：和内部CV的值共同构建边界上的导数FV例子*3.2.6守恒律的有限体积方法

Godunov格式3.2.6.1Godunov方法的思想一阶迎风格式(CIR格式)*用Godunov思想

说明CIR格式=Godunov格式*（三）偏微分方程的数值离散方法3.1有限差分法3.2有限体积法(有限元，谱方法，谱元，无网格，有限解析，边界元，特征线）3.1有限差分法*模型方程的差分逼近差分格式的构造差分方程的修正方程守恒型差分格式差分方法的理论基础偏微分方程的全离散方法0102030405063.1.1模型方程的差分逼近*3.1.2差分格式的构造*3.1.3差分方程的修正方程*04030102差分方程所精确逼近的微分方程称为修正方程对于时间发展方程，利用展开的方程逐步消去带时间的高阶导数，只留空间导数。Warming-Hyett方法：差分方程(2)写成算子的形式：3.1.3差分方程的修正方程(续)3.1.3差分方程的修正方程(续)3.1.4差分方法的理论基础*相容性，稳定性，收敛性01等价性定理02Fourier稳定性分析033.1.4差分方法的理论基础（续）*Fourier(VonNeumann)稳定性分析差分方法的理论基础（续）Fourier(VonNeumann)稳定性分（续）称为CFL条件(Courant,Friedrichs,Levy)123.1.5守恒型差分格式*1流体力学方程组描述物理量的守恒性；守恒律组：2定义3.1.5守恒型差分格式（续）*守恒性质：非守恒的差分格式一般没有对应于原始守恒律的“离散守恒律”。3.1.5守恒型差分格式（续）*守恒型差分格式的Lax-Wendroff定理：如果守恒型差分格式是和守恒律相容的，且当时间和空间步长趋于零时，差分解一致有界，几乎处处收敛于分片连续可微的函数，则这个收敛的函数就是守恒律的一个弱解。推论：守恒型差分各式的收敛解能自动满足间断关系。用途:(加上熵条件）可以得到正确的激波，研究中大量使用例如：Lax-Friedrichs格式，Lax-Wendroff格式，MacCormack格式3.1.6偏微分方程的全离散方法*对差分格式的一般要求：有精度、格式稳定、求解效率高特殊要求物理定律（守恒性）、物理特征（激波、湍流、旋涡、多介质、化学反应等）、有界性(正密度、正温度、正湍动能、正组分浓度等)主要指非定常方程的时间离散010302偏微分方程的全离散方法(续）*Crank-Nicolson格式、P-C格式、Lax-Wendroff格式、MacCormack格式Runge-Kutta方法时空全守恒：如Godunov格式、central-upwind格式、CESE方法两层格式Leap-Frog格式、Adams-Bashforth格式、后三点隐格式多层格式3.1.6.1两层格式*0102030405Crank-Nicolson格式Predictor-Corrector格式Lax-Wendroff格式MacCormack格式Runge-Kutta方法3.1.6.1两层格式(cont.）Lax-Wendroff格式一步LW格式3.1.6.1两层格式(cont.）两步LW格式常系数Jacobian时与单步LW等价。但计算更简单，不涉及矩阵相乘。Lax-Wendroff格式3.1.6.1两层格式(cont.）MacCormack格式(1969)两步格式比LW更简单，不需要计算函数在半点上的值。LW两步格式和MC各式的缺点：定常解的误差依赖于时间步长。01030204MacCormack格式的构造*3.1.6.2三层格式*Leap-Frog格式Adams-Bashforth格式第二课