水动力学计算 薛雷平.pptVIP

SIP: 要求 M? ? A ? (A?)P = 且设 对比系数有 4大型稀疏方程组迭代求解 NW N W P E S SE 4大型稀疏方程组迭代求解 代入 4大型稀疏方程组迭代求解 可得 ?=0 时，为普通的非完全三角分裂法 4大型稀疏方程组迭代求解 SIP方法的迭代矩阵M，可结合共轭梯度类方法结合，作为CG法（对称系数矩阵）以及BiCGSTAB（非对称系数矩阵）的前置优化矩阵使用。尤其对于多块结构化网格，可以取得更快的迭代收敛速度。 * * 谢谢！ * * BICGSTAB主要思想 * * * * BICGSTAB主要思想 BICGSTAB计算流程 * * BICGSTAB计算流程 * * CG 计算流程 * * 现代密码学理论与实践之五 水动力学计算若干问题2012年研究生暑期学校上海交通大学 薛雷平 提纲 数值计算要点 有限体积法 压力修正法 大型稀疏方程组迭代求解 * * 1 数值计算要点 基本模型的选择 计算方法的选择 边界条件和计算流程的选择 网格生成及其收敛性研究 结果分析和总结 * * 1 数值计算要点 基本模型的选择 可压缩非定场粘性流动方程 不可压缩NS方程 + 湍流模型（RANS） 势流方程（理想流体） 边界层方程 2D简化、几何简化、定解条件简化 … 适用范围 ? 简化模型的误差估计? 可通过典型算例不同层次的简化计算 得到合适的基本计算模型 * * 1 数值计算要点 计算方法的选择 离散方法：有限体积，有限差分，有限元… 计算格式： 伪可压缩模型的时间推进法（显式）， 压力修正法（隐式）… 湍流模型、多相流及空化模型、数值造波\消波… * * 1 数值计算要点 边界条件和计算流程的选择 边界位置及给出的边界条件并不一定要与实际问题完全一致 通过物理意义的考虑和试算得到能够简单实现， 但又不使实际问题失真的边界条件 注意点：来流速度分布和湍流度 非远场出口边界条件 对称几何的非对称流动 空化流动时的压力条件的实现 计算流程：对于复杂问题，通常可以先计算物理、几何简化的初始解，逐步加入复杂性 * * 1 数值计算要点 网格生成 截断误差取决与于 网格间距 和 解的梯度，合理分布网格，通过初步计算重新分布网格通常是必要的 粘性流动计算时边界层网格（最好是结构化网格）， 不同湍流模型对y+的要求 网格收敛性研究 网格收敛性研究指采用不断加密的网格使解收敛 不是指采用某一密度的网格时解与实验值一致 这是获得可靠数值计算结果必不可少的步骤 * * 1 数值计算要点 结果分析和总结 数值计算结果分析与实验结果分析类似 尽量采用无量刚表达式 CFD不是colour fluid dynamics，少用云图，多用曲线 分析总结的层次： 1：初步的定性规律（相关性和大小） 2：进一步的定性规律（线性、幂次、指数…，先用量纲分析和理论储备获得基本规律，多变量时考虑采用变量分离法） 3：定量规律（系数的拟合，加入特征点） 4：分析定量规律背后的物理机理 * * * * 2 有限体积法 流体力学基本方程的一般形式 对流扩散方程 随体变化=时间变化+对流 = 扩散项+体积源项 +面积源余项 * * 2 有限体积法 对流扩散方程的体积分形式 V为任意空间体积（控制体），F为体积V的表面 该输运方程可由 - 微分方程体积分后直接由GAUSS积分公式得到， - 也是控制体中物理量变化规律的描述 * * 2 有限体积法 将积分方程应用于一个微元有限体积 * * 2 有限体积法 对流项的离散 迎风格式 + 高阶修正 * * 2 有限体积法 扩散项的离散——梯度格式 * * 2 有限体积法 扩散项的离散——梯度格式 * * 2 有限体积法 扩散项的离散——梯度格式 * * 2 有限体积法 有限体梯度项的数值计算 * * 2 有限体积法 时间变化项的离散 2 有限体积法 时间变化项的离散 与对流项合并（对任意f: