计算流体CFD2004-2.pptVIP

第一章 引 言 1.3 一维问题的有限体积法FVM 1.4 Neumann边界条件 2.有限体积法 1.5 几个例子 1.5.2 NEUMANN边界条件 1.6总结 第二章 控制方程 影响域 一般二维偏微分方程 特征线种类 二、例：模型方程及其性质 2、模型方程的数学性质： 特征线上物理量连续而导数可以不连续(弱解) 波动方程 表：方程特征 三、方程组的数学性质： 2、方程组特征值 3、三维流体力学控制方程组的性质列表如下： 2.2 Navier-Stokes 2、N-S方程组 定常N-S方程组(二维标量式) ②守恒形式： 守恒变量 有限体积法 2.3 边界条件 举例 2.4 总结 * 华中科技大学能源与动力工程学院 1.2 有限差分法的一维计算 0x1 当x=0时 u=0 ; x=1时 u=0 Dirichlet边界条件 由此可得精确解： 差分解？ u 图1.2.1 有限差分近似 ui-ui-1 Δx ui+1-ui-1 2Δx ui+1-ui Δx Δx Δx i-1 i i+1 x x=1/2( )处 精确解： (1.2.4) 差分解： CS1(控制面1) CS2(控制面2) · × · × · i-1 i-1/2 i i+1/2 i+1 Δx/2 Δx=h Δx/2 CV1 CV2(控制体2) CV3 图1.3.1. 有限体积近似 (1.3.4)=(1.2.4) 1.有限差分法——①利用虚节点法 边界点3处应用努曼边界条件和有限差分方程有 i-1 i i+1 o o o o 1 2 3 4 u(0)=0 虚点 求解得到 ②节点3用二阶精度的差分方程 2点： 3点： 联立求解得到 1.5.1 Dirichlet边界条件：x=0时u=0；x=1时u=-1 精确解： u2=1/9; u3= -2/9 1、有限差分方法: O O O O 2 3 4 三元素四节点(FDM用) 节点2： 节点2： 节点3： 解： u2=1/9; u3= -2/9; 与精确解相同 2、有限体积法： 与FDM相同，该例中也能够给出精确解。 边界条件：x＝0时u＝0， x＝1时 ＝-3 精确解：u＝－2x2＋x; u2=1/9, u3= -2/9 有限差分方法: 三元素 不利用虚拟节点有：u2=-0.018, u3= -0.51 利用虚拟节点有：u1=0, u2=1/9, u3= -2/9, u4= -1 使用节点4上du/dx的前向二阶精度公式:可得精确解 使用节点4上d2u/dx2的二阶精度公式:也可得精确解 如果使用两元素，解的误差很大 有限体积法: 两元素 节点3利用FDM O O × O 1 2 5/2 3 0 h=1/2 1 节点2利用FDM 解得 u2=0; u3= -1; 与精确解u＝-2x2+x一致 本章的意图是通过非常简单的一维线性二阶微分方程组使读者了解有限差分，有限体积法，提供了彼此相同的最终形式的代数方程，给出了Dirichlet问题的同样结果。 Neumann边界条件在FDM中近似满足，但是在FVM中精确实现。他们在推导公式的过程中“自然”地显示出来。由于这个原因，Neumann边界条件被称作自然边界条件。对于FDM则不是这种情况，尽管对于简单的问题得到了精确的解。 尽管在