一种流体问题中物理量求解方法技术

本发明专利技术公开了一种流体问题中物理量求解方法，在实际流动问题计算区域内，生成直角坐标系并进行网格剖分；基于剖分后的计算区域，得到描述流体运动的Euler或Navier

【技术实现步骤摘要】
一种流体问题中物理量求解方法


[0001]本专利技术属于流体力学计算
，具体涉及一种流体问题中物理量求解方法。

技术介绍

[0002]高精度数值算法具备较低的数值耗散，相同条件下相比于低阶格式能够得到更为精细的流动结构，所以一直是计算流体力学、计算数学、计算电磁学、计算声学等领域的重点研究课题。在众多的高精度算法中，WENO格式扮演着非常重要的角色。WENO格式于1994年首次被Liu等人提出，并在1996年被Jiang和Shu完善并给出了任意高阶WENO格式的统一框架。基于该框架，众多学者展开了诸多的后续研究。比如：Balsara等人构造了7
‑
13阶精度的保单调WENO格式；Henrick等人设计了WENO
‑
M格式用于恢复传统WENO格式在极值点附近的数值精度；Dumbser等人发展了非结构网格体系下的有限体积WENO方法，并通过Tsoutsanis等人开发的UCNS3D代码实现了非结构WENO格式的工程化应用。上述WENO格式构造过程相对比较复杂，涉及繁琐的线性权计算过程。因此，Zhu等人于2016年设计了一类称之为非等距WENO的格式。该格式的线性权值可以任取一组和为一的正数，极大地简化了实施过程(尤其是非结构网格体系下)，并且对于定常问题具备更加优异的定常收敛性。
[0003]毫无疑问的是WENO格式是一类受欢迎的激波捕捉类高精度格式，具备极好的工程应用前景。但是，由于其构造过程中涉及特征分解以及复杂的非线性权值计算，从而导致其计算效率不高。因此，如何继承WENO格式优良特性的同时降低其计算成本一直是WENO格式的研究热点。杂交WENO格式是一种可行的替代方案，其通过间断识别器对网格单元进行分类，然后在光滑区域使用廉价高效的线性方法，而在间断区域使用昂贵但激波捕捉能力优异的WENO格式。这种结合使得数值算法能够充分发挥线性方法和非线性WENO格式各自的优点。但是正如Pirozzoli所言，在杂交WENO格式的过程中，对网格单元进行分类的间断识别器这一角色显得尤为重要，如何能够快速、准确、高效且无需人为参与的区分光滑流场单元和间断流场单元一直是一个开放的研究课题。

技术实现思路

[0004]本专利技术所要解决的技术问题是针对上述现有技术的不足，提供一种流体问题中物理量求解方法，设计一种无需人为参与且准确高效的间断识别器，然后对计算区域内的网格单元进行分类，分别利用廉价高效的线性方法和激波捕捉能力优异的非等距WENO格式对流体力学控制方程进行离散，最后利用具备TVD性质的三阶Runge
‑
Kutta时间离散方法实现控制方程在时间方向上的推进，进而得到任意时刻计算区域内所有网格点处的物理量，并通过无粘和黏性可压缩流动问题验证方法的有效性。
[0005]为实现上述技术目的，本专利技术采取的技术方案为：
[0006]一种流体问题中物理量求解方法，包括以下步骤：
[0007]步骤一、在实际流动问题计算区域内，生成直角坐标系并进行网格剖分；
[0008]步骤二、基于剖分后的计算区域，得到描述流体运动的Euler或Navier
‑
Stokes方
程组的空间离散化方程；
[0009]步骤三、基于所述空间离散化方程，对计算区域内的网格点进行分类并标记，得到光滑网格点和间断网格点；
[0010]步骤四、分别在光滑网格点处使用线性格式和间断网格点处使用非等距有限差分WENO格式对其空间导数项进行空间离散；
[0011]步骤五、基于步骤四，采用经典的Runge
‑
Kutta时间离散方法对流体力学控制方程组进行时间方向上的推进，从而得到任意时刻网格离散点处的物理量值，进而数值模拟出实际流动问题中各个时刻的状态。
[0012]为优化上述技术方案，采取的具体措施还包括：
[0013]上述的实际流动问题包括前台阶反射问题、Rayleigh
‑
Taylor不稳定问题、粘性激波管问题、三维激波
‑
剪切层相互作用问题和Taylor
‑
Green涡演化问题。
[0014]上述的步骤一中，对计算区域构建笛卡尔直角坐标系，并进行网格剖分，相应的网格点坐标记为I
i,j,k
＝(x
i
,y
j
,z
k
)，其中x，y，z表示空间变量，下标i，j，k表示网格点序号。
[0015]上述的步骤二中，描述流体运动的Euler或Navier
‑
Stokes方程组为：
[0016][0017]其中，t表示时间变量，x＝(x,y,z)表示空间变量向量，U(x,t)＝(ρ,ρu,ρv,ρw,E)
T
表示守恒变量向量，T表示转置；
[0018]表示总能，ρ、u、v、w、e都是关于时间变量t和空间变量x，y，z的函数，分别表示流体密度、x方向速度、y方向速度、z方向速度以及内能；
[0019]p＝(γ
‑
1)ρe表示压力，其中γ表示比热比常数；
[0020]▽
表示散度算子，σ为Euler方程和Navier
‑
Stokes方程的转换符号，σ＝0表示Euler方程，σ＝1表示Navier
‑
Stokes方程；
[0021]Re为Reynolds数，表示对流通量向量，表示黏性通量向量，上标x，y，z表示对应的空间方向，下标c和v表示对流和黏性；
[0022]F
c
(U)和F
v
(U,
▽
U)中各分量和的具体表达式为
[0023][0024][0025][0026][0027][0028][0029]其中，τ
xx
、τ
yy
、τ
zz
分别表示x，y，z方向的正应力，τ
xy
、τ
xz
、τ
yz
分表表示xy、xz、yz平面的切应力，q
x
、q
y
、q
z
表示x，y，z方向的热通量；
[0030]基于空间直角坐标系下的剖分网格，可以得到如下空间离散方程
[0031][0032]其中，U
i,j,k
表示守恒变量U(x,t)在网格点(x
i
,y
j
,z
k
)处的离散值，表示x方向对流通量在网格边界(x
i
±
1/2
,y
j
,z
k
)处的离散值，表示y方向对流通量F
cy
在网格边界(x
i
,y
j
±
1/2
,z
k
)处的离散值，表示z方向对流通量在网格边界(x
i
,y...

【技术保护点】

【技术特征摘要】
1.一种流体问题中物理量求解方法，其特征在于，包括以下步骤：步骤一、在实际流动问题计算区域内，生成直角坐标系并进行网格剖分；步骤二、基于剖分后的计算区域，得到描述流体运动的Euler或Navier
‑
Stokes方程组的空间离散化方程；步骤三、基于所述空间离散化方程，对计算区域内的网格点进行分类并标记，得到光滑网格点和间断网格点；步骤四、分别在光滑网格点处使用线性格式和间断网格点处使用非等距有限差分WENO格式对其空间导数项进行空间离散；步骤五、基于步骤四，采用经典的Runge
‑
Kutta时间离散方法对流体力学控制方程组进行时间方向上的推进，从而得到任意时刻网格离散点处的物理量值，进而数值模拟出实际流动问题中各个时刻的状态。2.根据权利要求1所述的一种流体问题中物理量求解方法，其特征在于，所述实际流动问题包括前台阶反射问题、Rayleigh
‑
Taylor不稳定问题、粘性激波管问题、三维激波
‑
剪切层相互作用问题和Taylor
‑
Green涡演化问题。3.根据权利要求1所述的一种流体问题中物理量求解方法，其特征在于，所述步骤一中，对计算区域构建笛卡尔直角坐标系，并进行网格剖分，相应的网格点坐标记为I
i,j,k
＝(x
i
,y
j
,z
k
)，其中x，y，z表示空间变量，下标i，j，k表示网格点序号。4.根据权利要求3所述的一种流体问题中物理量求解方法，其特征在于，所述步骤二中，描述流体运动的Euler或Navier
‑
Stokes方程组为：其中，t表示时间变量，x＝(x,y,z)表示空间变量向量，U(x,t)＝(ρ,ρu,ρv,ρw,E)
T
表示守恒变量向量，T表示转置；表示总能，ρ、u、v、w、e都是关于时间变量t和空间变量x，y，z的函数，分别表示流体密度、x方向速度、y方向速度、z方向速度以及内能；p＝(γ
‑
1)ρe表示压力，其中γ表示比热比常数；表示散度算子，σ为Euler方程和Navier
‑
Stokes方程的转换符号，σ＝0表示Euler方程，σ＝1表示Navier
‑
Stokes方程；Re为Reynolds数，表示对流通量向量，表示黏性通量向量，上标x，y，z表示对应的空间方向，下标c和v表示对流和黏性；F
c
(U)和中各分量和的具体表达式为的具体表达式为的具体表达式为的具体表达式为的具体表达式为
其中，τ
xx
、τ
yy
、τ
zz
分别表示x，y，z方向的正应力，τ
xy
、τ
xz
、τ
yz
分表表示xy、xz、yz平面的切应力，q
x
、q
y
、q
z
表示x，y，z方向的热通量；基于空间直角坐标系下的剖分网格，可以得到如下空间离散方程其中，U
i,j,k
表示守恒变量U(x,t)在网格点(x
i
,y
j
,z
k
)处的离散值，表示x方向对流通量在网格边界(x
i
±
1/2
,y
j
,z
k
)处的离散值，表示y方向对流通量在网格边界(x
i
,y
j
±
1/2
,z
k
)处的离散值，表示z方向对流通量在网格边界(x
i
,y
j
,z
k
±
1/2
)处的离散值，表示x方向黏性通量在网格边界(x
i
±
1/2
,y
j
,z
k
)处的离散值，表示y方向黏性通量在网格边界(x
i
,y
j
±
1/2
,z
k
)处的离散值，表示z方向黏性通量在网格边界(x
i
,y
j
,z
k
±
1/2
)处的离散值。5.根据权利要求4所述的一种流体问题中物理量求解方法，其特征在于，所述步骤三中，对计算区域内网格点进行分类并标记的具体过程如下：步骤31、对于目标单元，即目标网格点I
i,j,k
，先固定y方向和z方...

【专利技术属性】
技术研发人员：王镇明，朱君，田琳琳，赵宁，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：