第八讲 平流方案有限差分课件.pptx

第八讲 平流方案有限差分课程内容平流方程的物理背景平流方案的基本知识有限差分方案举例平流方程的物理背景平流方程：在流场给定的情况下，计算示踪物密度的演变。当边界无流入流出，以及无源汇情况时，示踪物总质量守恒尽管是被动的（示踪物不直接改变流场），但是与其它方程相耦合，属于复杂地球流体力学过程的一部分平流过程无处不在复杂的气候系统地球上的各种成分（示踪物）空间分布不连续性高大气中的平流方程?最重要的示踪物：水物质（水汽、云水、云冰）大气化学传输模式需要计算大量化学示踪物（、、气溶胶等）火山灰、人为污染物等上百种！计算效率成问题海洋中的平流方程平流位温、盐度密度由状态方程诊断得到平流生物地球化学、辐射污染物等示踪物海洋模式平流辐射污染物http://www.nersc.no/data/ocean-tracer-simulations平流方程通量形式：利于保持总质量守恒平流形式：质量混合比?平流方程是双曲型发展方程波动传播速度（即传播矩阵的特征值）是实数且有限，在平流方程中的传播速度即是流速波动沿特征线传播平流方程的波动性?代入波动解因此波速，群速度特征线流体质点沿特征线运动由于流速存在空间差异，流体微团不仅做平移运动，同时其形状也会发生改变描述平流方程的两种观点欧拉1707~1783拉格朗日1736~1813欧拉观点?适于将流体离散为固定网格拉格朗日观点?轨迹方程：?适于将流体离散为移动微团两种观点比较欧拉方案在气候模式中占主导地位，发展最为悠久和成熟，理论分析完善拉格朗日方案能够避免一些欧拉方案存在的缺陷，但需要面对流体形变的影响课程内容平流方程的物理背景平流方案的基本知识有限差分方案举例平流方案平流方案是求解平流方程的数值计算方案平流方案是动力框架的基础算法需要保持多种物理性质不被破坏需要保持的物理性质风场-质量场协调性示踪物间的相关关系保持性总质量守恒性计算精度计算稳定性单调性、正定性需要保持的物理性质总质量守恒性最基本的要求不满足守恒性的方案是无法进行长期气候模拟的需要保持的物理性质计算精度不仅对于光滑连续的示踪物要有高的精度（截断误差的阶数），对于不连续的示踪物也要能够有好的模拟效果Gibbs现象（Fourier展开对间断函数的振荡性）→ 高精度不代表好效果需要保持的物理性质计算稳定性在具有好效果的前提下，稳定性好的方案可以取较长时间步长，从而节省计算时间需要保持的物理性质单调性不会产生新的极值（避免小波动出现）正定性不会产生负的示踪物密度需要保持的物理性质质量-风场协调性动力框架求解的连续方程与平流方案计算存在差异只要两种计算方案不同就会有差异需要保持的物理性质随着更多系统（化学、生物）纳入气候模式，需要进行更多示踪物的平流计算（提高了对多示踪物计算效率的要求）多种示踪物存在一定相互关系（correlations），这种相互关系可能是由于化学反应或平流导致的（Thuburn and McIntyre 1997；Plumb 2007）具有化学活性的示踪物间的相互关系变化会激活化学反应（提高对示踪物间的相互关系保持的要求）摘自Plumb (2007)平流方案与模式动力框架平流方案平流方案分类欧拉方案有限差分方案有限体积方案有限元方案谱方案半拉格朗日方案追踪点的半拉格朗日方案追踪面积的半拉格朗日方案全拉格朗日方案微团无形状的全拉格朗日方案微团有形状的全拉格朗日方案有限差分方案有限差分方法就是将连续微分转变为离散差分（前面课程介绍过）将采用修正方程法（或Hirt启示法）获得差分格式对应的近似微分方程将采用von Neumann线性稳定性分析法考察查对光滑谐波的保持情况“Attempts to solve partial differential equations by finite-difference methods often end in disaster. A finite-difference equation may have a rapidly growing and oscillating solution that in no way resembles the solution expected from the differential equation. In such a situation the difference equation is said to be computationally unstable. Obviously, it is desirable to avoid these disasters.”~ C. W. Hirt (1968)修正方程分析方法Modified equation method（Hirt 1968；Warming and Hyett 1974）原理：代表方