计算流体力学CFD教案.ppt

方程的一般变换 从物理平面(x,y,t) 计算平面(?,?,?) 方程的一般变换 方程的一般变换 方程的一般变换 度量和雅可比行列式 度量和雅可比行列式 从物理平面(x,y,t) 计算平面(?,?,?) 涉及网格几何性质的项，如??/?x,??/?y,??/?x,??/?y称为度量。 度量和雅可比行列式 从物理平面(x,y,t) 计算平面(?,?,?) 如果上述变换式用解析形式给出，则度量也能得到解析值。 度量和雅可比行列式 从物理平面(x,y,t) 计算平面(?,?,?) 大部分情况下，上述关系式是用数值形式给出的，则度量用有限差分计算。 度量和雅可比行列式 从物理平面(x,y,t) 计算平面(?,?,?) 逆变换 下面推导用逆变换来表示导数 从计算平面(?,?,?) 物理平面(x,y,t) 度量和雅可比行列式 令 u的全微分为： 则： 度量和雅可比行列式 由 得： 度量和雅可比行列式 分母上的行列式称为雅可比行列式，记作 度量和雅可比行列式 由 得： 度量和雅可比行列式 度量和雅可比行列式 写成更一般的形式： 度量和雅可比行列式 用逆变换来表示导数（含J）： 度量和雅可比行列式 用直接变换来表示导数(不含J)： 下面根据逆度量和直接度量之间的关系式来推导怎么用逆变换来表示导数。 度量和雅可比行列式 考虑二维的直接变换： 有： 度量和雅可比行列式 即： 度量和雅可比行列式 再考虑二维的逆变换： 有： 误差与稳定性分析 由于误差?也满足差分方程，故有 而 则 误差与稳定性分析 解得 放大因子 误差与稳定性分析 要使 放大因子 必须满足 误差与稳定性分析 上式就是差分方程 的稳定性条件。 对于给定的?x，?t的值必须足够小，才能满足上述稳定性条件，以保证计算过程中误差不会放大。 误差与稳定性分析 稳定性条件的具体形式取决于差分方程的形式。 的差分方程 是无条件不稳定的。 比如，一阶波动方程： 误差与稳定性分析 但如果用 则 （Lax方法） 误差与稳定性分析 令误差 则放大因子 式中 误差与稳定性分析 则放大因子 稳定性要求 则 误差与稳定性分析 稳定性要求 式中的C称为柯朗(Courant)数。 误差与稳定性分析 稳定性要求 上式称为柯朗－弗里德里奇－列维(Courant-Friedrichs-Lewy)条件，一般写成CFL条件。 误差与稳定性分析 下面来看CFL条件的物理意义。 CFL条件： 也是二阶波动方程： 的稳定性条件。 误差与稳定性分析 下面来看CFL条件的物理意义。 二阶波动方程： 的特征线为 CFL条件的物理意义：要保证稳定性，数值解的依赖区域必须全部包含解析解的依赖区域。 误差与稳定性分析 CFL条件的物理意义：要保证稳定性，数值解的依赖区域必须全部包含解析解的依赖区域。 误差与稳定性分析 计算流体力学CFD(3) 网格生成与坐标变换 引言 引言 在CFD里，确定适当的网格是一件非常重要的事情。 网格会影响数值计算的成功与失败。 引言 标准的有限差分方法需要均匀网格。 如果在流场内生成了非均匀网格，也需要将它变换成均匀分布的矩形网格。 引言 采用均匀网格计算翼型绕流有如下问题： 1）有些网格点落入翼型内部，而不是在流场里，如何给定这些点上的流动参量？ 引言 采用均匀网格计算翼型绕流有如下问题： 2）只有少量的网格点落在翼型表面上，这也不好。因为翼型表面是极其重要的边界，翼型表面上的边界条件确定了整个流动。 引言 1）翼型内部没有网格点 2）网格点落在翼型表面上 引言 网格既不是矩形的，也不是均匀分布的。通常的差分很难应用，必须转换。 引言 控制方程从物理平面(x,y)转换到计算平面(?,?) 物理平面 计算平面 贴体网格 方程的一般变换 方程的一般变换 考虑二维非定常流场， 从物理平面(x,y,t) 计算平面(?,?,?) 物理平面 计算平面 方程的一般变换 从物理平面(x,y,t) 计算平面(?,?,?) 下标表示求偏导数过程中保持不变的量 方程的一般变换 从物理平面(x,y,t) 计算平面(?,?,?) 下标表示求偏导数过程中保持不变的量 方程的一般变换 从物理平面(x,y,t) 计算平面(?,?,?) 方程的一般变换 方程的一般变换 方程的一般变换 显式方法 边界条件已知 显式方法 显式方法中每一个差分方程只包含一个未知数，从而这个未知数可以用直接计算的方法显式地求解。显式方法是最简单的方法。 隐式方法 隐式方法 克兰克－尼科尔森格式 隐式方法 对于排列在同一时间层所有网格点上的未知量，必须将它们联立起来同时求解，才能求出这些未知量，这种方法就定义为隐式方法。 隐式方法 由于需要求解联立的代数方程组，隐式方法通常涉及大型矩阵的运算。隐式方法比显式方法需要更多、更复杂的计算。 隐式方法 隐式方