极坐标中气体动力学的可压缩欧拉方程与二维人工扩散的确切公式是什么?

计算科学 pde 计算几何 双曲-pde
2021-12-06 08:38:07

感兴趣的方程是对流-扩散方程。典型的例子之一是 Navier-Stokes 方程。但是,我想让扩散系数常数变为零,ϵ0,而网格尺寸,Δx,Δdy0, 和时间步长,Δt0,趋近于零。然后,我可以得到无粘可压缩方程,可以作为气体动力学欧拉方程的一个例子。以下方程是二维笛卡尔坐标中的标量平流-扩散方程。

dfdt+Udfdx+Vdfdy=ϵ(d2fdx2+d2fdy2)

主要兴趣是在圆柱坐标系中找到气体动力学欧拉方程的精确公式(r,θ)带有人工扩散部分。

编辑:我对二维气体动力学欧拉方程感兴趣的主要方程。笛卡尔向量形式可以表示为

(ρρuρvE)t+(ρuρu2+pρuvu(E+p))x+(ρvρuvρv2+pv(E+p))y=0

在这里,质量守恒、-x 和 -y 方向的动量和能量方程可以以矩阵形式看到。但是,我想在 Navier-Stokes 方程中添加具有小常数因子的人工粘度,以平滑数值模拟期间可能发生的振荡。出于这个原因,我需要知道 Navier-Stokes 方程的柱坐标中的向量拉普拉斯形式。然后我可以应用上述理论。

1个回答

看看Manuel Baumann 的学士论文“Modellierung und Simulation von Dispersionen in turbulenten Strömungen”的第 3 章。

它是德文的,但可以理解方程式和图片。

在第 3.1.1 章中,解释了特别影响对流的坐标变换。

在第 3.2.1 章中,公式化了变换后的RANS方程。将正确的项设置为零,您将得到欧拉方程(第 29 页的方程 (19))。

对于欧拉方程,主要工作是对流项,在极坐标中写为

在此处输入图像描述

请注意,在极坐标中,坐标单位向量eθ,er取决于角坐标 θ. 这意味着,例如,

θ(Uθeθ)=θ(Uθ)eθUθer
是什么使附加项出现在上述等式的最后一个等式中。

作为参考,我可以推荐 Longwell 的“流体流动力学”。

免责声明:我正在监督上述论文。