[知识分享]WENO格式保姆级教程(1):Euler方程的通量表达式

绿颊鹦鹉

于 2024-12-28 16:59:18 发布

阅读量473

点赞数 6

分类专栏：计算流体力学文章标签：算法学习

本文链接：https://blog.csdn.net/m0_51867190/article/details/144791348

版权

计算流体力学专栏收录该内容

1 篇文章

订阅专栏

写这篇文章，一方面是对我的学习过程做个总结，另一方面是想帮助后来的同学快速入门。我将从初学者的视角分享WENO-JS格式，力求详细、易懂，把各种背景知识介绍清楚。各位大佬珠玉在前，学生只是拾人牙慧，水平有限，还请各位不吝赐教。

Euler方程的通量表达式

WENO格式的作用是对N-S方程中的对流项做离散化。为了便于研究，我们去掉N-S方程中的粘性项和源项，保留瞬态项和对流项。这样N-S方程就退化为Euler方程了：

$\frac{\partial{U}}{\partial{t}}+\frac{\partial{F_{1}(U)}}{\partial{x}}+\frac{\partial{F_{2}(U)}}{\partial{y}}+\frac{\partial{F_{3}(U)}}{\partial{z}}=0$

守恒量

$U=\begin{pmatrix} \rho \\ \rho u \\ \rho v \\ \rho w \\ E \end{pmatrix}$ ,

包含质量密度、动量密度和能量密度。流通矢量

$F_{1}(U)=\begin{pmatrix} \rho u\\ \rho u^2+p \\ \rho uv \\ \rho uw \\ u(E+P) \end{pmatrix}$ , $F_{2}(U)=\begin{pmatrix} \rho v\\ \rho vu \\ \rho v^2+p \\ \rho vw \\ v(E+P) \end{pmatrix}$ , $F_{3}(U)=\begin{pmatrix} \rho w\\ \rho wu \\ \rho wv \\ \rho w^2+p \\ w(E+P) \end{pmatrix}$ ,

包含质量、动量和能量的通量。有兴趣可以验证一下。接下来我将以有限差分法为例进行说明。

假如我们将计算区域划分成很多个小的正方体单元，单元中心点的守恒量 $U|_{(i,j,k)}$ 已知，而单元边界的信息未知(有限差分法的性质)。将单元

$x\in{(i-\frac{1}{2},i+\frac{1}{2})}, y\in{(j-\frac{1}{2},j+\frac{1}{2})}, z\in{(k-\frac{1}{2},k+\frac{1}{2})}$

视作一个控制体，如图1所示。

图1

将Euler方程写成半离散的形式：

$\frac{\partial{U}}{\partial{t}}|_{(i,j,k)} + \frac {\hat{F}_{1,i+1/2} - \hat{F}_{1,i-1/2}} {\Delta{x}}+ \frac {\hat{F}_{2,j+1/2} - \hat{F}_{2,j-1/2}} {\Delta{y}}+ \frac {\hat{F}_{3,k+1/2} - \hat{F}_{3,k-1/2}} {\Delta{z}} =0$

可以看出守恒量的更新完全依赖于单元边界的数值通量。WENO格式就是一种构造数值通量的方法。

注释：1.式(1.2)采用了守恒格式，将流通矢量的导数表达为 $\frac{\partial{F_{2}}} {\partial{y}}|_{j} = \frac {\hat{F}_{2,j+1/2} - \hat{F}_{2,j-1/2}} {\Delta{y}}$ ,其中 $\hat{F}_{2,j+1/2}$ 是单元边界处的数值通量，上标"^"表示估计值。