动翼型边界条件

计算科学 流体动力学 边界条件 有限体积
2021-12-11 08:45:43

我正在尝试以恒定速度(马赫=0.755,aoa=1.25)模拟移动的翼型。我用 Roe 的方法求解欧拉方程。我有两个边界条件:

  1. 远处的荒(田)野
  2. 滑动壁(翼型)

对于所有面孔,我给出的面孔速度为:

vbx = -Mach cos (aoa)

vby = -Mach sin (aoa)

在远场,我设置

(rho,u,v,p) = (1, 0, 0, 1/gamma)

其中 gamma 是比热比。

对于滑动壁,我将速度设置为边界单元附近相邻单元的切向速度。

流在任何地方都用零(u,v)初始化。但是在这些边界条件下,残差为零,因此没有必要进行迭代。残差是通量的总和,其中通量定义为:

Fn0 = rho*(un-vb)

Fn1 = rho*(un-vb)*u + p*nx

Fn2 = rho*(un-vb)*v + p*ny

Fn3 = E*(un-vb) + p*un

那么应该怎样才是合适的边界条件呢?

1个回答

我发现了问题。我在更新翼型表面时忘记添加面速度。它应该是:

u = tangX + vbx
v = tangY + vby

其中,tangX 和 tangY 是相邻内部单元的切向速度的 x 和 y 分量。那么,翼型表面的通量应该是:

Fn0 = 0;
Fn1 = p*nx;
Fn2 = p*ny;
Fn3 = un * p;

请注意,由于面速度,un 不为零,但 (un-vb) 为零。