1) TLNS equations
薄层Navier-Stokes方程
2) Navier-Stokes equations
Navier-Stokes方程
1.
The boundary treatment of the fourth-order compact finite difference scheme for the incompressible Navier-Stokes equations;
关于不可压流体Navier-Stokes方程的四阶精度有限差分紧致格式的边界处理
2.
A fourth-order finite volume compact method for the incompressible Navier-Stokes equations;
求解不可压流体Navier-Stokes方程的四阶精度有限容积紧致格式
3.
Third-order projection method for solving the incompressible Navier-Stokes equations;
求解不可压Navier-Stokes方程的三阶精度投影方法
3) Navier Stokes equations
Navier-Stokes方程
1.
A cell centered finite volume scheme that applies a multistage time stepping scheme in conjunction with steady state acceleration techniques is used to solve the integral form of mass averaged three dimensional Navier Stokes equations for flows over a swept NACA0012 wing and a thin low aspect ratio wing mounted in a wind tunnel.
应用三维可压、雷诺平均 Navier-Stokes方程数值模拟了机翼半模实验风洞侧壁干扰和三维机翼半模与安装侧壁结合部流场。
2.
In this paper, the harmonic oscillating flows were solved by unsteady Navier Stokes equations, and the calculation method of damping in roll derivatives was presented under the Etkin s theory when body oscilated around a fixed axis.
采用非定常 Navier-Stokes方程描述物体简谐振动流场 ,并在 Etkin理论下给出绕定轴转动时滚转阻尼导数的计算公式 ,定常流场的计算采用空间二阶精度的交替方向隐式分解的 NND格式 ,非定常流场的计算采用时、空二阶精度的 Runge-Kutta多步格式 ,采用代数方法生成物体静、动网格。
3.
This paper is devoted to the study of the Navier Stokes equations with variable viscosity in a noncylindrical domain.
本文研究的是非柱状区域中关于可变粘性的Navier-Stokes方程的解的存在性与唯一性。
4) Navier-Stokes equation
Navier-Stokes方程
1.
Finite Element Analysis on Three-dimensional Navier-Stokes Equation;
三维Navier-Stokes方程的有限元分析
2.
Solving Navier-Stokes equation by mixed interpolation method;
混合插值法求解Navier-Stokes方程(英文)
3.
Global attractor for the penalized 2D Navier-Stokes equation;
带惩罚项的二维Navier-Stokes方程的全局吸引子
5) N-S equations
Navier-Stokes方程
1.
A new method is worked out for numerical simulation of 3-D unsteady incompressible N-S equations.
对三维、非定常、不可压Navier-Stokes方程提出了一种新的数值计算方法。
补充资料:边界层方程数值解法
边界层理论是德国L.普朗特在20世纪初建立起来的。当流体流经物体表面时,靠近壁面边界很薄的一层,粘性效应很重要。利用粘性边界层很薄的特点,可以把流体力学运动方程(即纳维-斯托克斯方程)中量级较小的各项忽略掉,简化成为边界层方程。边界层理论为粘性流体力学的应用开辟了广阔的道路,在近代力学中起着重要的作用。
以平面问题为例:定常二维不可压缩流的边界层方程组,由一个连续性方程和两个动量方程组成,即
式中u、v为沿着x、y方向上的速度分量;p、ρ和v分别表示压力、密度和运动粘性系数。边界条件要求在不渗透的固体表面上,两个速度分量为零。在边界层外缘,u渐近地等于外缘速度ue(x),所以有:
(2)另外,还要给定压力梯度дp/дx。由于式(1c)中的压力p只是x的函数,它与外缘速度之间的关系为:
。方程组(1)是非线性偏微分方程组,求解很困难,一般需用数值方法,这里主要介绍相似性解法和差分解法。
相似性解法 其要点是引进无量纲相似参数,将偏微分方程转换成常微分方程,然后再用数值方法求解。德国Н.布拉西乌斯在1907年首次用此法解压力为常数的平板绕流问题。在连续性方程中引进流函数Ψ,即u=дΨ/дy,v=-дΨ/дx,并定义一个相似参数同时令f(η) 为无量纲的流函数。速度分量u、v及其导数дu/дy和д2u/дy2均可以从Ψ 求出,而且都可以用函数 f(η)及其高阶导数表示。最后,原方程组(1)变成一个三阶常微分方程:
f冺+ff″=0,
(3)对应于边界条件(2), 要求f(0)=f′(0)=0,f′(∞)=1。这是两点边值问题。一般的作法是先假设f″(0)=α, 从η=0的地方对方程(3)进行数值积分。当η→∞时,要求f′(η)→1。如果条件不能满足,必须更改 α的初值,反复迭代到满足 f′(∞)=1的条件为止。但通过变数的转换,也可将这个两点边值问题换成初值问题,求解时不需要反复迭代。令ζ=α1/3η,α仍然代表f″(0);再令f(η)=α1/3F(ζ),则f′(η)=α2/3F′(ζ),f″(η)=αF″(ζ),f冺(η)=α4/3F冺(ζ)。代入方程式(3),得到一个同样形式的方程:
F冺(ζ)+F(ζ)F″(ζ)=0,(4)
但边界条件有些不同,变成F(0)=F′(0)=0,F″(0)=1三个初始条件,正好用数值积分直接求F(ζ),而后利用f′(∞)=1=α2/3F′(∞)求α,即
(5)方程(4)的具体解法, 是把它改为三个一阶常微分方程,令F的一阶导数为G,二阶导数为H,则有:
F′=G,G′=H,H′+FH=0,
(6)
F、G、H为三个未知变数,相应的初始条件为:F(0)=0,G(0)=0,H(0)=1。这组一阶常微分方程可用一般的数值积分法求解。
差分解法 这种解法是将微分算符近似地用差商代替,把微分方程改为差分方程然后再求解。在有压力梯度的流动中,相似条件不能满足。用前面相同的坐标变换,即但此处应令由于相似性假设不适用,流函数f是ξ、η的函数。通过坐标转换,方程(1b)变为: ,
(7)式中 f′、f″、f冺 均为 η 的导数;f 为 ξ 的导数;为压力梯度参数。差分-微分方程是将上式的 ξ导数项改用差分形式,而在η方向仍保持微分形式。这样,方程(7)变成在 η 方向上的常微分方程,具有在η=0,η=∞的两点边界条件,可用迭代法求解。近来,人们直接将边界层方程的所有偏导数均用差分表示。这类差分法的格式很多(见有限差分方法),现以凯勒的差分格式为例。 此法首先将原方程〔如方程(7)〕改写成几个一阶偏微分方程组,而后将所有一阶导数均用中心差分,给出具有二阶精度的差分方法。现将 f(ξ,η)对 η的一阶导数用 g(ξ,η)表示,二阶导数用h(ξ,η)表示。方程(7)可改为:
(8a)
。 (8b)上两式均在点上取值,它们的差分方程为:
(9a)
(9b)方程(8b)则在点上取值,如
在这些式子中,还有一些非线性项,如g卾,(fh)i+1,须进行线性化,如果把gi+1和gi的差值看作小量,并忽略小量二阶以上的项,即得出线性化关系式:
将以上各式代入(8b),即可得出在i+1截面上的线性差分方程。连同(9a)和(9b)一起,并结合相应的边界条件,便可联立求解三个未知量f、g和h。从f即可求流函数Ψ,从而可计算出两个速度分量u和v。
以平面问题为例:定常二维不可压缩流的边界层方程组,由一个连续性方程和两个动量方程组成,即
式中u、v为沿着x、y方向上的速度分量;p、ρ和v分别表示压力、密度和运动粘性系数。边界条件要求在不渗透的固体表面上,两个速度分量为零。在边界层外缘,u渐近地等于外缘速度ue(x),所以有:
(2)另外,还要给定压力梯度дp/дx。由于式(1c)中的压力p只是x的函数,它与外缘速度之间的关系为:
。方程组(1)是非线性偏微分方程组,求解很困难,一般需用数值方法,这里主要介绍相似性解法和差分解法。
相似性解法 其要点是引进无量纲相似参数,将偏微分方程转换成常微分方程,然后再用数值方法求解。德国Н.布拉西乌斯在1907年首次用此法解压力为常数的平板绕流问题。在连续性方程中引进流函数Ψ,即u=дΨ/дy,v=-дΨ/дx,并定义一个相似参数同时令f(η) 为无量纲的流函数。速度分量u、v及其导数дu/дy和д2u/дy2均可以从Ψ 求出,而且都可以用函数 f(η)及其高阶导数表示。最后,原方程组(1)变成一个三阶常微分方程:
f冺+ff″=0,
(3)对应于边界条件(2), 要求f(0)=f′(0)=0,f′(∞)=1。这是两点边值问题。一般的作法是先假设f″(0)=α, 从η=0的地方对方程(3)进行数值积分。当η→∞时,要求f′(η)→1。如果条件不能满足,必须更改 α的初值,反复迭代到满足 f′(∞)=1的条件为止。但通过变数的转换,也可将这个两点边值问题换成初值问题,求解时不需要反复迭代。令ζ=α1/3η,α仍然代表f″(0);再令f(η)=α1/3F(ζ),则f′(η)=α2/3F′(ζ),f″(η)=αF″(ζ),f冺(η)=α4/3F冺(ζ)。代入方程式(3),得到一个同样形式的方程:
F冺(ζ)+F(ζ)F″(ζ)=0,(4)
但边界条件有些不同,变成F(0)=F′(0)=0,F″(0)=1三个初始条件,正好用数值积分直接求F(ζ),而后利用f′(∞)=1=α2/3F′(∞)求α,即
(5)方程(4)的具体解法, 是把它改为三个一阶常微分方程,令F的一阶导数为G,二阶导数为H,则有:
F′=G,G′=H,H′+FH=0,
(6)
F、G、H为三个未知变数,相应的初始条件为:F(0)=0,G(0)=0,H(0)=1。这组一阶常微分方程可用一般的数值积分法求解。
差分解法 这种解法是将微分算符近似地用差商代替,把微分方程改为差分方程然后再求解。在有压力梯度的流动中,相似条件不能满足。用前面相同的坐标变换,即但此处应令由于相似性假设不适用,流函数f是ξ、η的函数。通过坐标转换,方程(1b)变为: ,
(7)式中 f′、f″、f冺 均为 η 的导数;f 为 ξ 的导数;为压力梯度参数。差分-微分方程是将上式的 ξ导数项改用差分形式,而在η方向仍保持微分形式。这样,方程(7)变成在 η 方向上的常微分方程,具有在η=0,η=∞的两点边界条件,可用迭代法求解。近来,人们直接将边界层方程的所有偏导数均用差分表示。这类差分法的格式很多(见有限差分方法),现以凯勒的差分格式为例。 此法首先将原方程〔如方程(7)〕改写成几个一阶偏微分方程组,而后将所有一阶导数均用中心差分,给出具有二阶精度的差分方法。现将 f(ξ,η)对 η的一阶导数用 g(ξ,η)表示,二阶导数用h(ξ,η)表示。方程(7)可改为:
(8a)
。 (8b)上两式均在点上取值,它们的差分方程为:
(9a)
(9b)方程(8b)则在点上取值,如
在这些式子中,还有一些非线性项,如g卾,(fh)i+1,须进行线性化,如果把gi+1和gi的差值看作小量,并忽略小量二阶以上的项,即得出线性化关系式:
将以上各式代入(8b),即可得出在i+1截面上的线性差分方程。连同(9a)和(9b)一起,并结合相应的边界条件,便可联立求解三个未知量f、g和h。从f即可求流函数Ψ,从而可计算出两个速度分量u和v。
说明:补充资料仅用于学习参考,请勿用于其它任何用途。
参考词条