圆柱周围潜在流动的适当边界条件

计算科学 pde 数值分析 流体动力学 边界条件 计算物理学
2021-12-02 09:10:11

我正在使用一般坐标中的离散化计算围绕圆柱体的静止、不可压缩、无粘性和无旋流动。

我在由下式给出的物理域(左图)中推导出了 PDE 和适当的边界条件G:={(x1,x2:R12x12+x22R22,x20}. 这个偏微分方程是,其中是扰动势。 Δϕ=0ϕ物理和逻辑域

然后任务要求在一般坐标中简化拉普拉斯算子。为此计算协变/逆变基向量和雅可比行列式。接下来我定义映射x(ξ)

但我坚持在逻辑域中推导适当的边界条件(右图)。这些必须根据协变/逆变基向量和度量张量来定义。作为提示,给出:在这个阶段不要使用特定的映射。使用域的边界形式为常数这一事实,因此边界上的单位法向量由ξα(x)=±a(α)/|a(α)|

我看不到如何使用协变/逆变基向量。我们为极坐标系的物理域 G 有了新的基础,那么它们与计算域有什么关系呢?

不变形式的边界条件:

  1. 沿圆柱体(流动无粘性),其中给出的速度场给出的总速度势函数现在方向上的自由流速度,而是扰动势。0=un=Un1uuα=ΦxαΦΦ=Ux1+ϕUx1ϕ
  2. 在物理域的水平边界上,我们有,所以ux2=02ϕx1x2=02ϕx22=0
  3. 沿着域的顶部,我们有一个均匀的狄利克雷边界条件用于扰动势。对于此边界上ϕ(x)=0x

我遗漏了很多细节,因为我不确定这些是否需要。

1个回答

小心 Neumann BC。并不意味着,例如u/x=02u/x2=0u=x2

您在物理域中定义了 BC...并且您有一个与物理坐标相关的表达式x在计算域中具有计算坐标 ,反之亦然。在您的 BC 中替换这些并应用链式规则。如果地图定义明确,那么计算它的逆将不会有任何问题。ξx=x(ξ)