CN106682398A - 一种基于湍动能和湍流长度尺度截断的rans/les方法 - Google Patents
一种基于湍动能和湍流长度尺度截断的rans/les方法 Download PDFInfo
- Publication number
- CN106682398A CN106682398A CN201611136614.8A CN201611136614A CN106682398A CN 106682398 A CN106682398 A CN 106682398A CN 201611136614 A CN201611136614 A CN 201611136614A CN 106682398 A CN106682398 A CN 106682398A
- Authority
- CN
- China
- Prior art keywords
- rans
- les
- tubulence energy
- turbulent flow
- energy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于湍动能和湍流长度尺度截断的RANS/LES方法。其主要是建立一种新的湍流粘性系数的截断公式Fr,得到截断湍流粘性系数。Fr=1,计算呈现RANS,0<Fr<1,呈现类LES。根据分子粘性系数的定义和湍流尺度的能谱满足:E(κ)∝ε2/3κ‑5/3推导出Fr计算公式。结合标准K‑ε模型,可实现RANS/LES方法。通过上述方式,本发明可以在较少网格条件下,计算出速度脉动,并且得到的平均速度与DNS结果较吻合。
Description
技术领域
本发明涉及流体力学(Computational Fluid Dynamics,CFD)领域,尤其涉及一种基于湍动能和湍流长度尺度截断的RANS/LES方法。是在近壁面区域求解RANS方程,在远壁面区域呈现类LES计算的湍流模型。
背景技术
湍流流动广泛存在于大自然和工程领域中,例如大气、海洋和瀑布等流动,飞机机翼和船体附近流体的流动,流体机械内的流动都属于湍流流动。湍流是是十分复杂的多尺度随机不规则的流动,而研究湍流对于大气环境预测,飞机机翼和船体周围的受力分析以及流体机械的机械强度分析有着重要的意义。由于描述流体流动动量守恒的Navier-Stokes方程的解析解只在极少数情况下才存在,故理论分析对于湍流研究极其局限。湍流研究方法主要集中在实验和数值模拟。由于实验费用昂贵且部分实验条件难以重现,加上计算机的迅速发展,因而采用数值模拟方法研究湍流日益得到重视。
湍流数值模拟方法有三种:直接数值模拟(Direct Numerical Simulation,DNS)、大涡模拟(Large Eddy Simulation,LES)和雷诺平均模拟(Reynolds Averaged NavierStokes,RANS)。DNS不需要对湍流建立模型,采用数值计算方法直接求解控制方程。由于湍流尺度的复杂性,要获得所有湍流尺度的信息,对于空间和时间的分辨率要求极高,需要巨大的计算内存和时间,网格数量大概为Re2.7,对于较小Re=1000,网格数量级就将到达108,所以DNS只能求解简单湍流问题,而对于复杂湍流问题难以求解。LES和RANS对控制方程进行过滤得到过滤方程:
式中ui和uj是速度,ρ和p是密度和压力,μ是粘性系数。对于RANS,上横线表示时间平均,湍流粘性系数为K是湍动能,ε是湍动能耗散率,Cμ=0.09。对于LES,上横线表示空间平均,湍流粘性系数为μt=ρ(CsΔ)2S,Δ是网格体积,S是应变率张量大小,Cs=0.1。K和ε的控制方程为:
式中Pk=2μtSijSij,σk=1.0,σε=1.3,Cε1=1.44,Cε2=1.92。
RANS方法在工程模拟领域占据主导位置,应用广泛,但是对于复杂流动,例如分离与再附着流动,RANS表现比较差,难以模拟出主导作用的流动结构。LES方法的基本思想是模式化小涡,求解大涡,能够捕捉到比RANS更多的湍流结构,但是LES对于网格精度要求比RANS高很多,网格数量到达Re1.76,这主要是因为LES在靠近壁面处对于网格的要求极其高。为结合RANS和LES方法的优点,研究人员发展了混合RANS/LES方法,其基本思想:在靠近壁面区域求解呈现RANS方法,远离壁面区域呈现类LES方法。应用比较广泛的混合RANS/LES有分离涡模拟(Detached Eddy Simulation,DES,部分平均NS模拟(Partially AveragedNavier Stokes,PANS)以及超大涡模拟(Very Large Eddy Simulation,VLES),DES的思想是对湍流长度尺度进行截断,PANS对湍动能进行截断,VLES对湍动能进行截断。
发明内容
本发明的目的在于克服上述现有技术的缺点和不足,提供一种基于湍动能和湍流长度尺度截断的RANS/LES方法,以解决湍流模型问题,能够综合RANS和LES的优点。
本发明通过下述技术方案实现:
一种基于湍动能和湍流长度尺度截断的RANS/LES方法,建立新的湍流粘性系数的截断公式Fr,得到截断湍流粘性系数:
其中
Fr=1,计算呈现RANS,0<Fr<1,呈现类LES;其中υ是粘性系数,K是湍动能,ε是湍动能耗散率,ρ是密度,Cμ=0.09。
步骤一:在RANS标准K-ε的基础上,对涡粘系数修正;针对平板槽道内充分发展湍流计算,建立连续性、动量、湍动能和湍动能耗散率方程;
步骤二,对平板通道网格划分;
步骤三,利用步骤一和步骤二的网格,通过有限体积法计算,得到计算区域的平均速度场和瞬时速度场,与DNS结果进行对比。
所需的网格数量小于DNS所需的网格,并且计算结果与DNS结果吻合。
本发明的核心是得出截断湍流粘性系数。
对于RANS,湍流粘性系数类比分子粘性系数,正比于湍流速度尺度u′和长度尺度L,有以下关系:
μt∝u′L u′∝K1/2→μt∝K1/2Li (5)
分别对于湍动能和湍流长度尺度截断后,湍流粘性系数依然符合式(5)
μtc∝Kc 1/2Lc (6)
定义截断系数,Fr=1时,呈现RANS,0<Fr<1,呈现类LES。
当截断尺度处于湍流惯性子区时,湍动能能谱有以下关系:
E(κ)∝ε2/3κ-5/3 (8)
可得
κK,κc,κi分别是Kolmogorov尺度波数,截断波数,积分尺度波数,与各自长度尺度关系如下:
式(10)代入(9)可得
为保证Fr在0到1范围内,可做以下限制:
可得截断湍流系数为:
结果式(1)(2)(3)(4)(13)可实现本发明的一种基于湍动能和湍流长度尺度截断的RANS/LES模拟。
本发明通过截断湍动能和湍流长度尺度来实现混合RANS/LES。
本发明可以在较少网格条件下,计算出速度脉动,并且得到的平均速度与DNS结果较吻合。
附图说明
图1计算区域。
图2速度单等值面图。
图3速度云图一
图4速度云图二
图5平均流向速度曲线图。
具体实施方式
下面结合图1至5及实施例对本发明作进一步具体详细描述。
平板通道内充分发展湍流是验证湍流模型的基准算例。
为验证本发明的正确性,接下来将湍流模型结果来对比DNS的模拟结果进行对比。具体如下:
图1是计算区域,即平板通道的示意图,尺寸为x=3.2δ,y=2δ,z=1.6δ,x方向为流向,y为垂直壁面方向,z为展向。网格数量为35×70×35,流向和展向平均离散,垂直壁面方向,靠近壁面的第一个网格点设置在y+=1.0处。流体流动的特征参数表1。
表1流动参数
为求解方程(1)(2)(3)(4)(13),采用有限体积法(FVM)对方程进行离散。动量方程采用二阶中心差分格式,湍动能和湍动能耗散率方程擦用二阶迎风格式,时间步进采用二阶隐性格式,压力速度耦合采用SIMPLEC方法,壁面模型采用增强壁面函数模型,时间步长设置为0.002s,保证CFL<1,总时间为100倍流过时间,以减小统计结果误差。流向和展向设置为周期性边界条件,流向压力梯度为1Pa/m。
通过数值计算足够时间,统计结果稳定。图2、图3和图4分别是计算结果的速度等值面图和不同y值面速度云图,可以看到,本发明可以有效地计算出壁面附近的速度脉动。对比DNS,如图5所示,利用本发明得到的平均流向速度与DNS结果吻合较好。
如上所述,便可较好地实现本发明。
本发明的实施方式并不受上述实施例的限制,其他任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (3)
1.一种基于湍动能和湍流长度尺度截断的RANS/LES方法,其特征在于:建立湍流粘性系数的截断公式Fr,得到截断湍流粘性系数:
其中
Fr=1,计算呈现RANS,0<Fr<1,呈现类LES;其中υ是粘性系数,K是湍动能,ε是湍动能耗散率,ρ是密度,Cμ=0.09。
2.根据权利要求1所述基于湍动能和湍流长度尺度截断的RANS/LES方法,其特征在于:
步骤一:在RANS标准K-ε的基础上,对涡粘系数修正;针对平板槽道内充分发展湍流计算,建立连续性、动量、湍动能和湍动能耗散率方程;
步骤二,对平板通道网格划分;
步骤三,利用步骤一和步骤二的网格,通过有限体积法计算,得到计算区域的平均速度场和瞬时速度场,与DNS结果进行对比。
3.根据权利要求2所述基于湍动能和湍流长度尺度截断的RANS/LES方法,其特征在于:所需的网格数量小于DNS所需的网格,并且计算结果与DNS结果吻合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611136614.8A CN106682398A (zh) | 2016-12-12 | 2016-12-12 | 一种基于湍动能和湍流长度尺度截断的rans/les方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611136614.8A CN106682398A (zh) | 2016-12-12 | 2016-12-12 | 一种基于湍动能和湍流长度尺度截断的rans/les方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106682398A true CN106682398A (zh) | 2017-05-17 |
Family
ID=58868819
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611136614.8A Pending CN106682398A (zh) | 2016-12-12 | 2016-12-12 | 一种基于湍动能和湍流长度尺度截断的rans/les方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106682398A (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108319793A (zh) * | 2018-02-12 | 2018-07-24 | 大连理工大学 | 一种基于cfd的低速倾斜入水自由运动的分析方法 |
CN109858148A (zh) * | 2019-01-30 | 2019-06-07 | 南京航空航天大学 | 一种基于部分滤波的湍流计算方法 |
CN111191394A (zh) * | 2019-12-23 | 2020-05-22 | 国家海洋环境预报中心 | 一种气泡羽流湍流闭合的方法 |
CN111751396A (zh) * | 2020-07-17 | 2020-10-09 | 北京唯实兴邦科技有限公司 | 一种机械结构件微结构损伤失效检测分析方法 |
CN112464581A (zh) * | 2020-10-23 | 2021-03-09 | 黄河勘测规划设计研究院有限公司 | 一种基于流速滑移与能耗耦合湍流模型的流动控制方法 |
CN113111610A (zh) * | 2021-05-10 | 2021-07-13 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种新型亚格子尺度模型建立方法 |
CN113158340A (zh) * | 2021-04-16 | 2021-07-23 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种针对k-epsilon湍流模型的湍流长度尺度修正方法 |
CN113378492A (zh) * | 2021-06-10 | 2021-09-10 | 杭州电子科技大学 | 一种基于电磁修正的低磁雷诺数磁流体湍流数值计算方法 |
CN114970403A (zh) * | 2022-07-12 | 2022-08-30 | 北京航空航天大学 | 一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法 |
CN115034162A (zh) * | 2022-07-12 | 2022-09-09 | 北京航空航天大学 | 一种基于湍流能谱耦合k-ε系列模型的网格自适应湍流模拟方法 |
CN115438601A (zh) * | 2022-09-19 | 2022-12-06 | 中国航空发动机研究院 | 基于les自适应网格滤波尺度的混合rans-les方法和相关设备 |
CN116070071A (zh) * | 2023-04-04 | 2023-05-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 壁面函数的选择方法、装置、终端设备及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103742893A (zh) * | 2013-12-24 | 2014-04-23 | 华电重工股份有限公司 | 余热锅炉内部噪声的控制方法及装置 |
CN103968362A (zh) * | 2013-12-24 | 2014-08-06 | 华电重工股份有限公司 | 用于余热锅炉内部的穿孔吸声体的设计方法、结构及应用 |
CN104699947A (zh) * | 2015-01-14 | 2015-06-10 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种采用rans/les混合技术模拟飞行器摇滚运动的方法 |
CN104794293A (zh) * | 2015-04-24 | 2015-07-22 | 南京航空航天大学 | 风力机尾流计算方法 |
CN105183975A (zh) * | 2015-09-01 | 2015-12-23 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种基于尾涡流场传递的多飞行器飞行编队数值模拟方法 |
CN105527383A (zh) * | 2016-01-13 | 2016-04-27 | 长春工程学院 | 一种复杂环境下多工况汽车尾气测量方法 |
-
2016
- 2016-12-12 CN CN201611136614.8A patent/CN106682398A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103742893A (zh) * | 2013-12-24 | 2014-04-23 | 华电重工股份有限公司 | 余热锅炉内部噪声的控制方法及装置 |
CN103968362A (zh) * | 2013-12-24 | 2014-08-06 | 华电重工股份有限公司 | 用于余热锅炉内部的穿孔吸声体的设计方法、结构及应用 |
CN104699947A (zh) * | 2015-01-14 | 2015-06-10 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种采用rans/les混合技术模拟飞行器摇滚运动的方法 |
CN104794293A (zh) * | 2015-04-24 | 2015-07-22 | 南京航空航天大学 | 风力机尾流计算方法 |
CN105183975A (zh) * | 2015-09-01 | 2015-12-23 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种基于尾涡流场传递的多飞行器飞行编队数值模拟方法 |
CN105527383A (zh) * | 2016-01-13 | 2016-04-27 | 长春工程学院 | 一种复杂环境下多工况汽车尾气测量方法 |
Non-Patent Citations (2)
Title |
---|
宋科等: "多段翼型大迎角分离流动的 DELAYEDRANS/LES混合算法", 《航空计算基数》 * |
王乐等: "基于可压缩修正 k-ε 模型的弹丸气动计算", 《弹箭与制导学报》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108319793A (zh) * | 2018-02-12 | 2018-07-24 | 大连理工大学 | 一种基于cfd的低速倾斜入水自由运动的分析方法 |
CN109858148B (zh) * | 2019-01-30 | 2023-06-09 | 南京航空航天大学 | 一种基于部分滤波的湍流计算方法 |
CN109858148A (zh) * | 2019-01-30 | 2019-06-07 | 南京航空航天大学 | 一种基于部分滤波的湍流计算方法 |
CN111191394A (zh) * | 2019-12-23 | 2020-05-22 | 国家海洋环境预报中心 | 一种气泡羽流湍流闭合的方法 |
CN111751396A (zh) * | 2020-07-17 | 2020-10-09 | 北京唯实兴邦科技有限公司 | 一种机械结构件微结构损伤失效检测分析方法 |
CN112464581A (zh) * | 2020-10-23 | 2021-03-09 | 黄河勘测规划设计研究院有限公司 | 一种基于流速滑移与能耗耦合湍流模型的流动控制方法 |
CN112464581B (zh) * | 2020-10-23 | 2023-06-20 | 黄河勘测规划设计研究院有限公司 | 一种基于流速滑移与能耗耦合湍流模型的流动控制方法 |
CN113158340A (zh) * | 2021-04-16 | 2021-07-23 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种针对k-epsilon湍流模型的湍流长度尺度修正方法 |
CN113111610A (zh) * | 2021-05-10 | 2021-07-13 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种新型亚格子尺度模型建立方法 |
CN113111610B (zh) * | 2021-05-10 | 2022-10-14 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种亚格子尺度模型建立方法 |
CN113378492A (zh) * | 2021-06-10 | 2021-09-10 | 杭州电子科技大学 | 一种基于电磁修正的低磁雷诺数磁流体湍流数值计算方法 |
CN113378492B (zh) * | 2021-06-10 | 2024-03-12 | 杭州电子科技大学 | 一种基于电磁修正的低磁雷诺数磁流体湍流数值计算方法 |
CN114970403A (zh) * | 2022-07-12 | 2022-08-30 | 北京航空航天大学 | 一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法 |
CN115034162A (zh) * | 2022-07-12 | 2022-09-09 | 北京航空航天大学 | 一种基于湍流能谱耦合k-ε系列模型的网格自适应湍流模拟方法 |
CN115034162B (zh) * | 2022-07-12 | 2024-06-07 | 北京航空航天大学 | 一种基于湍流能谱耦合k-ε系列模型的网格自适应湍流模拟方法 |
CN115438601A (zh) * | 2022-09-19 | 2022-12-06 | 中国航空发动机研究院 | 基于les自适应网格滤波尺度的混合rans-les方法和相关设备 |
CN115438601B (zh) * | 2022-09-19 | 2023-06-20 | 中国航空发动机研究院 | 基于les自适应网格滤波尺度的混合rans-les方法和相关设备 |
CN116070071A (zh) * | 2023-04-04 | 2023-05-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 壁面函数的选择方法、装置、终端设备及存储介质 |
CN116070071B (zh) * | 2023-04-04 | 2023-06-27 | 中国空气动力研究与发展中心计算空气动力研究所 | 壁面函数的选择方法、装置、终端设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106682398A (zh) | 一种基于湍动能和湍流长度尺度截断的rans/les方法 | |
Vasaturo et al. | Large eddy simulation of the neutral atmospheric boundary layer: performance evaluation of three inflow methods for terrains with different roughness | |
Li et al. | Large-eddy simulation of flow and pollutant dispersion in high-aspect-ratio urban street canyons with wall model | |
Liu et al. | On the prediction of air and pollutant exchange rates in street canyons of different aspect ratios using large-eddy simulation | |
Liu et al. | Large-eddy simulation of flow and scalar transport in a modeled street canyon | |
Schumann et al. | Numerical study of pressure and velocity fluctuations in nearly isotropic turbulence | |
CN115525978B (zh) | 用于飞行器波浪情况水动力分析的速度入口数值造波方法 | |
Weerasuriya | Computational Fluid Dynamic (CFD) simulation of flow around tall buildings | |
CN103631992A (zh) | 一种自吸泵自吸过程流动模拟的计算方法 | |
CN110276090A (zh) | 一种基于相干结构的湍流大涡模拟方法 | |
Lee et al. | Performance of different turbulence models in predicting flow kinematics around an open offshore intake | |
Omari | CFD simulations of lid driven cavity flow at moderate Reynolds number | |
Zadravec et al. | The influence of rotating domain size in a rotating frame of reference approach for simulation of rotating impeller in a mixing vessel | |
Vosper et al. | Flow separation and rotor formation beneath two‐dimensional trapped lee waves | |
Bicer et al. | Numerical simulation of cavitation phenomena in diesel injector nozzles | |
Atzori et al. | High-resolution large-eddy simulations of simplified urban flows | |
Beg et al. | Investigation of the flow field inside a drainage system: Gully-Pipe-Manhole | |
Gillissen | Turbulent drag reduction using fluid spheres | |
Liu et al. | Study of unsteady cavitating flow around Clark-Y hydrofoil using nonlinear PANS model with near-wall correction | |
Agafonova et al. | Limits of wasp modelling in comparison with сfd for wind flow over two-dimensional hills | |
El-Okda et al. | Large-eddy simulation of flow over a surface-mounted prism using a high-order finite-difference scheme | |
Panopoulos et al. | CFD simulation of street canyon phenomenon as gas-liquid two-phase flow separation mechanism | |
Liu et al. | Numerical study of turbulent flow over complex topography by LES model | |
da Silva | CFD Analysis of the Air Flow Through Horizontal Axis Wind Turbines | |
Wittwer et al. | Statistical analysis of wind tunnel and atmospheric boundary layer turbulent flows |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170517 |
|
RJ01 | Rejection of invention patent application after publication |