(七)不可压缩流的数值方法 7.1 MAC方法 7.2 投影法 7.3 人工压缩性方法 7.4 SIMPLE方法 7.5 其他方法: 谱方法、谱元法、有限元法
不可压N-S方程 不可压N-S方程 MAC法 (1965), 和投影法Chorin (1968), Temam(1969) 投影法(分数步法),解不可压非定常流使用最多的方法
不可压方程的求解困难 不是时间发展型方程,而是每个时间步带散度为零约束条件。 半交错网格或者一般同位网格上,中心差分使得压力没有唯一解,或振荡,或奇偶失联。
7.1 MAC法 (Marker and Cell) Harlow and Welch (1965), 属于后来的投影法(1968,1969) 右边的称为MAC网格 优点:符合有限体积法概念和守恒特性,边界条件可以得到正确实施。
MAC法:方程
p u v MAC法: 离散点其他值的定义
MAC法的计算步骤
MAC法:稳定性要求 (14a),(14b)显式格式,基于简化分析:
边界条件 速度 p u v
p u v 边界条件 压力
7.2 投影法(1):基础 Theorem of orthogonal decomposition (an instance of the general Hodge orthogonal decomposition) 思想:先计算一个不满足散度为零的中间速度场u*(通常不计压力的效应),然后将它分解为满足散度为零的速度un+1和正比于压力的标量场的梯度之和。
投影法(2): 算法 显式投影法由三个步骤组成
投影法 (3):分析 分析:(3)+(4)
其他类型的投影法 CN (Kim and Moin)
其他类型的投影法 CNPC (pressure correction, Van Kan 1986)
其他类型的投影法 Backward 3-point second order PC (Shen et al): assume its stability is better than CN
精度 速度 L2-norm:second order 压力 L2-norm: 阶 或 阶(rotational) 数值边界层