首页/文章/ 详情

一文深入了解8类湍流模型,是8类,不是8个

1月前浏览1724

在湍流的计算建模中,一个共同的目标是获得一个能够预测感兴趣量(如流体速度)的模型,以便在模拟系统工程设计中使用。对于湍流流动,涉及湍流现象的长度尺度范围和复杂性使得大多数建模方法成本过高;解决湍流中涉及的所有尺度所需的分辨率超出了计算能力。在这种情况下,主要的方法是创建数值模型来近似未解决的现象。这里列出了一些常用的湍流流动计算模型。

湍流模型可以根据计算成本进行分类,这与被模拟与被解析的尺度范围相对应(解析的湍流尺度越多,模拟的分辨率越细,因此计算成本越高)。如果大多数或全部湍流尺度都没有被模拟,计算成本非常低,但代价是准确性降低。

除了广泛的长度和时间尺度以及相关的计算成本外,流体动力学的控制方程包含一个非线性对流项和一个非线性及非局部的压力梯度项。这些非线性方程必须在适当的边界和初始条件下用数值方法求解。


1. 雷诺平均纳维-斯托克斯方程

   
   

Reynolds-averaged Navier–Stokes 雷诺平均纳维-斯托克斯(RANS)方程是湍流建模最古老的方法。解决的是控制方程的集 合版本,引入了新的表观应力,称为雷诺应力。这增加了一个二阶张量的未知数,各种模型可以提供不同程度的封闭。一个常见的误解是,RANS方程不适用于具有时间变化平均流的流动,因为这些方程是“时间平均”的。实际上,统计上非稳态(或非定常)流动也可以同样处理。有时这被称为URANS。雷诺平均没有固有的东西可以阻止这一点,但用于封闭方程的湍流模型只有在这些平均变化的时间与包含大部分能量的湍流运动的时间尺度相比很大时才有效。

RANS模型可以分为两种广泛的途径:

a)Boussinesq假设

这种方法涉及使用雷诺应力的代数方程,包括确定湍流粘度,并根据模型的复杂程度,求解传输方程以确定湍流动能和耗散。模型包括k-ε(Launder和Spalding),混合长度模型(Prandtl)和零方程模型(Cebeci和Smith)。这种方法中可用的模型通常根据与方法相关联的传输方程数量来引用。例如,混合长度模型是“零方程”模型,因为没有求解传输方程;k-ε是“二方程”模型,因为求解了两个传输方程(一个用于k,一个用于ε)。

b)雷诺应力模型(RSM)

这种方法尝试实际解决雷诺应力的传输方程。这意味着引入了几个传输方程,用于所有雷诺应力,因此这种方法在CPU工作上的成本要高得多。


2. 大涡模拟

   
   

Large eddy simulation大涡模拟(LES)是一种技术,在这种技术中,通过过滤操作去除流动的最小尺度,并使用子网格尺度模型来建模它们的效果。这允许解析湍流的最大和最重要的尺度,同时大大减少了最小尺度引起的计算成本。这种方法需要比RANS方法更多的计算资源,但比DNS要便宜得多。


3. 分离涡模拟

   
   

Detached eddy simulation分离涡模拟(DES)是RANS模型的修改,在足够细分的区域中,模型切换到子网格尺度表述,以便进行LES计算。靠近固体边界和湍流长度尺度小于网格最大尺寸的区域被分配RANS模式的解决方案。当湍流长度尺度超过网格尺寸时,使用LES模式解决这些区域。因此,DES的网格分辨率不像纯LES那样苛刻,从而大幅降低了计算成本。尽管DES最初是为Spalart-Allmaras模型(Philippe R. Spalart等人,1997年)制定的,但它也可以通过适当修改RANS模型中明确或隐含涉及的长度尺度,与其他RANS模型(Strelets,2001年)一起实现。因此,基于Spalart-Allmaras模型的DES就像带壁面模型的LES,而基于其他模型(如二方程模型)的DES则表现为混合RANS-LES模型。由于RANS-LES切换,网格生成比简单的RANS或LES情况更复杂。DES是一种非区域性方法,为解决方案的RANS和LES区域提供了单一平滑的速度场。


4. 直接数值模拟

   
   

Direct numerical simulation直接数值模拟(DNS)解析了湍流长度尺度的全部范围。这边缘化了模型的效果,但极其昂贵。计算成本与Re^3成正比。对于具有复杂几何形状或流动配置的流动,DNS是不可行的。


5. 相干涡旋模拟

   
   

Coherent vortex simulation相干涡旋模拟方法将湍流场分解为相干部分,包括有组织的涡旋运动,以及不相干部分,即随机背景流动。这种分解是通过小波滤波完成的。这种方法与LES有很多共同之处,因为它使用分解并只解析过滤部分,但不同之处在于它不使用线性低通滤波器。相反,过滤操作基于小波,并且过滤器可以根据流动场的演变进行调整。Farge和Schneider使用两种流动配置测试了CVS方法,并表明流动的相干部分表现出与总流动的能量谱-40/39相对应的相干结构(涡旋管),而流动的不相干部分组成了均匀背景噪声,没有表现出有组织的结构。Goldstein和Vasilyev将FDV模型应用于大涡模拟,但没有假设小波滤波器从子滤波尺度中消除了所有相干运动。通过同时使用LES和CVS过滤,他们表明SFS耗散由SFS流动场的相干部分主导。


6. PDF方法

   
   

Probability density function (PDF) methods,由Lundgren首次引入的概率密度函数(PDF)湍流方法,基于跟踪速度的一点PDF,fV(v;x,t)dv,它给出了在点x的速度在v和v+dv之间的概率。这种方法类似于气体的动力学理论,其中气体的宏观属性由大量粒子描述。PDF方法的独特之处在于它们可以应用于不同湍流模型的框架内;主要差异发生在PDF传输方程的形式上。例如,在大涡模拟的背景下,PDF变为过滤后的PDF。PDF方法也可以用来描述化学反应,并且特别适用于模拟化学反应流动,因为化学源项是封闭的,不需要模型。PDF通常通过使用拉格朗日粒子方法来跟踪;当与大涡模拟结合时,这导致了一个用于子滤波粒子演化的Langevin方程。


7. 涡度约束方法

   
   

Vorticity confinement method,涡度约束(VC)方法是一种在模拟湍流尾迹时使用的欧拉技术。它使用类似孤立波的方法来产生一个没有数值扩散的稳定解。VC可以在仅有2个网格单元内捕获小尺度特征。在这些特征内,求解一个非线性差分方程,而不是有限差分方程。VC类似于捕获冲击波的方法,其中满足守恒定律,以便准确计算基本积分量。


8. 线性涡旋模型

   
   

Linear eddy model,线性涡旋模型是一种用于模拟湍流中发生的对流混合的技术。具体来说,它提供了一种数学方法来描述标量变量在矢量流场内的相互作用。它主要用于湍流的一维表示,因为它可以应用于广泛的长度尺度和雷诺数范围。这个模型通常用作更复杂流动表示的构建块,因为它提供了在大范围流动条件下保持的高分辨率预测。


邀您关注

▽ 纯粹CFD:软件教程、行业应用、专业理论、

基础科普、研究前沿、严选培训广告

▽ 只聊CFD相关的大小事,信手天成,娓娓道来



来源:CFD饭圈
FluentCFX非线性燃烧化学湍流PolyflowUM理论科普控制ParaView
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-09-29
最近编辑:1月前
CFD饭圈
硕士 分享CFD文章,感谢关注
获赞 22粉丝 21文章 376课程 0
点赞
收藏
作者推荐

20个CFD常用方程,各自都有假设前提,你知道几个?

CFD可以看作是一组计算方法(下面讨论),用于解决控制流体流动的方程。在应用CFD时,一个关键步骤是决定需要使用哪些物理假设和相关方程来解决手头的问题。为了说明这一步,以下总结了20个常用方程中所采取的物理假设/简化。注意,以下一些方程可能可以通过多种方式推导。1. 守恒定律:Conservation laws (CL),这些是CFD中考虑的最基本方程,例如,所有以下方程都可以从它们推导出来。对于单相、单种、可压缩流,考虑质量守恒、线性动量守恒和能量守恒。2. 连续介质守恒定律:Continuum conservation laws (CCL),从CL开始。假设质量、动量和能量是局部守恒的:这些量是守恒的,不能“瞬间移动”到另一个地方,只能通过连续流动移动(见连续性方程)。另一种解释是,从CL开始,并假设一个连续介质(见连续介质力学)。得到的方程组是未封闭的,因为要解决它需要进一步的关系/方程:(a) 粘性应力张量的本构关系;(b) 扩散热通量的本构关系;(c) 状态方程(EOS),如理想气体定律;以及,(d) 将温度与焓或内能等量联系起来的热状态方程。3. 可压缩纳维-斯托克斯方程:Compressible Navier-Stokes equations (C-NS),从CCL开始。假设牛顿粘性应力张量(见牛顿流体)和傅里叶热通量(见热通量)。C-NS需要增加一个EOS和一个热状态方程,才能有一个封闭的方程组。4. 不可压缩纳维-斯托克斯方程:Incompressible Navier-Stokes equations (I-NS),从C-NS开始。假设密度始终且无处不在是常数。另一种获得I-NS的方法是假设马赫数非常小,并且流体中的温度差异也非常小。结果是,质量守恒和动量守恒方程与能量守恒方程解耦,因此只需要解决前两个方程。5. 可压缩欧拉方程:Compressible Euler equations (EE),从C-NS开始。假设无摩擦流动且无扩散热通量。6. 弱可压缩纳维-斯托克斯方程:Weakly compressible Navier-Stokes equations (WC-NS),从C-NS开始。假设密度变化仅依赖于温度而不依赖于压力。例如,对于理想气体,使用 𝜌 = 𝑝0 /(𝑅𝑇),其中𝑝0是方便定义的参考压力,始终且无处不在是常数,𝜌是密度,𝑅是特定气体常数,𝑇是温度。结果是,WC-NS不捕捉声波。在WC-NS中也常见忽略能量守恒方程中的压强功和粘性加热项。WC-NS也称为低马赫数近似的C-NS。7. 布辛涅斯克方程:Boussinesq equations,从C-NS开始。假设密度变化始终且无处不在可以忽略,除了动量守恒方程中的重力项(其中密度乘以重力加速度)。还假设各种流体属性如粘度、热导率和热容始终且无处不在是常数。布辛涅斯克方程在微观气象学中广泛使用。8. 可压缩雷诺平均纳维-斯托克斯方程和可压缩Favre平均纳维-斯托克斯方程:Compressible Reynolds-averaged Navier–Stokes equations和compressible Favre-averaged Navier-Stokes equations (C-RANS and C-FANS),从C-NS开始。假设任何流动变量𝑓,如密度、速度和压力,可以表示为𝑓=𝐹+𝑓′′,其中𝐹是任何流动变量的集 合平均,𝑓′′是从这个平均的扰动或波动。𝑓′′不一定很小。如果𝐹是经典的集 合平均(见雷诺分解),则获得雷诺平均纳维-斯托克斯方程。如果𝐹是密度加权的集 合平均,则获得Favre-平均纳维-斯托克斯方程。结果是,根据雷诺数,运动的尺度范围大大减少,这导致与解决C-NS相比,解决方案更快。然而,信息丢失了,得到的方程组需要封闭各种未封闭的项,特别是雷诺应力。9. 理想流动或势流方程:Ideal flow or potential flow equations,从EE开始。假设流体粒子无旋转(零涡度)和无流动膨胀(零散度)。结果的流场完全由几何边界决定。在现代CFD中,理想流动可以用于初始化模拟。10. 线性化可压缩欧拉方程:Linearized compressible Euler equations (LEE),从EE开始。假设任何流动变量𝑓,如密度、速度和压力,可以表示为𝑓=𝑓0+𝑓′,其中𝑓0 是某个参考或基本状态的流动变量值,𝑓′是从这种状态的扰动或波动。此外,假设这种扰动𝑓′与某个参考值相比非常小。最后,假设𝑓0满足“它自己的”方程,如EE。LEE及其许多变体在计算气动声学中广泛使用。11. 声波或声波方程:Sound wave or acoustic wave equation,从LEE开始。忽略𝑓0和𝑓′的所有梯度,并假设参考或基本状态的马赫数非常小。结果的密度、动量和能量方程可以操作成一个压力方程,给出著名的声波方程。12. 浅水方程:Shallow water equations (SW),考虑在墙附近的流动,其中与墙平行的长度尺度远大于与墙垂直的长度尺度。从EE开始。假设密度始终且无处不在是常数,忽略垂直于墙的速度分量,并考虑与墙平行的速度在空间上是常数。13. 边界层方程:Boundary layer equations (BL),从可压缩(不可压缩)边界层的C-NS(I-NS)开始。假设在墙附近有薄区域,垂直于墙的空间梯度远大于平行于墙的梯度。14. 伯努利方程:Bernoulli equation,从EE开始。假设密度变化仅依赖于压力变化。见伯努利原理。15. 稳态伯努利方程:Steady Bernoulli equation,从伯努利方程开始,假设稳态流动。或从EE开始,假设流动是稳态的,并沿流线积分结果方程。16. 斯托克斯流动或蠕动流动方程:Stokes Flow或creeping flow equations,从C-NS或I-NS开始。忽略流动的惯性。这种假设在雷诺数非常低时是合理的。结果,得到的方程组是线性的,这大大简化了它们的解决。17. 二维通道流动方程:Two-dimensional channel flow equation ,考虑在两个无限平行板之间的流动。从C-NS开始。假设流动是稳态的、二维的、充分发展的(即,速度剖面沿流向方向不变)。注意,这种广泛使用的充分发展假设在某些情况下可能不充分,如某些可压缩、微通道流动,这种情况下可以被局部充分发展假设所取代。18. 一维欧拉方程或一维气体动力学方程:One-dimensional Euler equations or one-dimensional gas-dynamic equations (1D-EE),从EE开始,假设所有流动量仅依赖于一个空间维度。19. 范诺流动方程:Fanno flow equation,考虑在具有恒定面积和绝热壁的管道内的流动。从1D-EE开始。假设稳态流动,无重力效应,并在动量守恒方程中引入经验项以恢复壁摩擦效应(在EE中被忽略)。为了封闭范诺流动方程,需要这种摩擦项的模型。这种封闭涉及问题依赖的假设。20. 瑞利流动方程:Rayleigh flow equation,考虑在具有恒定面积的管道内的流动,管道要么有非绝热壁且无体积热源,要么有绝热壁但有体积热源。从1D-EE开始。假设稳态流动,无重力效应,并在能量守恒方程中引入经验项以恢复壁热传递或热源效应(在EE中被忽略)。邀您关注▽ 纯粹CFD:软件教程、行业应用、专业理论、基础科普、研究前沿、严选培训广告▽ 只聊CFD相关的大小事,信手天成,娓娓道来来源:CFD饭圈

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈