交坐标变换将物理空间不规则的计算域转换为计算空间规则的计算域坐标变换
如下
t x
x
x
式中
u
d
u
d
在新坐标下上下壁面分别位于
和
新旧坐标微分之间的关系为
t
t
t
x
x
x
x
x
x
x
x
i
i
i
其中
i
i
i
并且
i
x
i
i
i
i
i
t
t
需要指出的是Laplace 算子可表示为
x
j
x
j
j
j
S
式中
S
j
j
j
j
j
j
j
是由于坐标变换所产生的附加项在文献的公式中上式右端最后 项被遗漏
数 值 方 法
21时间推进
坐标变换后NS 方程变为
u
i
t
u
i
F
i
x
i
Re
u
i
x
j
x
j
和连续方程一起 控制方程利用 Karniadakis 等
所介绍的时间分裂法来进行时间推进该方
法曾被用于常规槽道湍流的直接数值模拟
在时间分裂法中每一个完整的时间步可以分
为 个子步非线性步压力步和粘性步在非线性步中引入速度场 u
s
i
满足
u
s
i
J
i
q
q
u
s q
i
J
e
q
q
t
u
i
F
i
s q
在上式中显式计算非线性项精度为 J
i
阶其他线性项采用隐式处理精度为 J
e
阶s 代表时间
运动壁面槽道流动的直接数值模拟