CN110489709B - 基于可压缩流动的解析壁面函数的数值模拟方法 - Google Patents
基于可压缩流动的解析壁面函数的数值模拟方法 Download PDFInfo
- Publication number
- CN110489709B CN110489709B CN201910705306.XA CN201910705306A CN110489709B CN 110489709 B CN110489709 B CN 110489709B CN 201910705306 A CN201910705306 A CN 201910705306A CN 110489709 B CN110489709 B CN 110489709B
- Authority
- CN
- China
- Prior art keywords
- equation
- bottom layer
- term
- temperature
- flow
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了一种基于可压缩流动的解析壁面函数的数值模拟方法,包括:根据可压缩流动特点建立Navier‑Stockes方程,并进行简化,然后解析得到全湍流区速度方程、粘性底层速度方程、全湍流区温度方程和粘性底层温度方程,并进一步定义应力方程和热流方程;给定粘性底层温度Tv的初值,并根据Tv计算得到应力和热流;利用计算得到的应力和热流更新应力项和热流项;再通过计算湍动能生成项和湍动能耗散项的平均量更新湍动能方程中的生成项和耗散项的值;最后利用粘性底层温度方程重新计算Tv,在下个时间步循环中重复进行更新。本发明基于可压缩流动特点,特别适合于高超声速流动,可以更加准确的预测壁面热流。
Description
技术领域
本发明涉及计算流体力学领域,尤其是一种基于可压缩流动的解析壁面函数的数值模拟方法。
背景技术
目前大多数工程应用的壁面函数基于局部平衡假设,在壁面附近速度和温度满足对数定律,因此在非平衡流条件下,标准壁面函数只能得到很低的预测精度。然而工业应用中的复杂流动通常包含流动分离和再附,这些流动都不满足局部平衡条件,因此将标准壁面函数应用到复杂流动问题上是不合适的。除了标准壁面函数,解析壁面函数由于在壁面处不涉及太多假设,在粗网格上的预测精度可以接近低雷诺数模型的结果,而计算时间比低雷诺数低一到两个数量级。由于壁面函数在实际编程,计算中的鲁棒性和可操作性等优点,都使得其工程应用得到广泛关注。而解析壁面函数在可压缩流体中的应用,由于需要考虑流体可压缩性,动量方程和能量方程的耦合等,其方法研究对于先进壁面函数在复杂超声速和高超声速流动中的应用也具有实际意义。
发明内容
本发明所要解决的技术问题是:针对上述存在的问题,提供一种基于可压缩流动的解析壁面函数的数值模拟方法。
本发明采用的技术方案如下:
一种基于可压缩流动的解析壁面函数的数值模拟方法,包括如下步骤:
步骤1,在壁面网格内根据可压缩流动特点建立Navier-Stockes方程;
步骤2,简化步骤1建立的Navier-Stockes方程,得到壁面动量方程和壁面能量方程;
步骤3,通过对壁面动量方程两次积分得到全湍流区速度方程和粘性底层速度方程;并通过对壁面能量方程两次积分得到全湍流区温度方程和粘性底层温度方程;
步骤4,根据全湍流区速度方程和粘性底层速度方程,以及全湍流区温度方程和粘性底层温度方程,定义应力方程和热流方程;
步骤5,给定粘性底层温度Tv的初值;
步骤6,根据粘性底层温度Tv计算壁面网格中粘性底层的密度和层流粘性系数;
步骤7,将步骤6得到的粘性底层的密度和粘性系数代入应力方程和热流方程,计算得到应力和热流;并利用计算得到的应力和热流更新步骤1建立的Navier-Stockes方程中的应力项和热流项;
步骤8,通过解析全湍流区速度方程,计算湍动能生成项和湍动能耗散项的平均量;并利用湍动能生成项和湍动能耗散项的平均量更新湍动能方程中的生成项和耗散项的值;
步骤9,利用粘性底层温度方程重新计算粘性底层温度Tv,在下个时间步循环中重复执行步骤6~9,形成闭环的数值模拟回路。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
本发明的基于可压缩流动的解析壁面函数的数值模拟方法,基于可压缩流动特点,特别适合于高超声速流动,由于其壁面动量方程和壁面能量方程之间的耦合更加紧密,可以更加准确的预测壁面热流。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明基于可压缩流动的解析壁面函数的数值模拟方法的流程框图。
图2为壁面网格内Navier-Stockes方程各项对比图。
图3为壁面网格内层流粘性系数示意图。
图4为壁面网格内湍流粘性系数示意图。
图5为壁面网格内湍动能破坏项示意图。
图6a为在马赫数为5,入射角分别为10°的壁面压力对比图。
图6b为在马赫数为5,入射角分别为14°的壁面压力对比图。
图7a为在马赫数为5,入射角分别为10°的摩擦系数对比图。
图7b为在马赫数为5,入射角分别为14°的摩擦系数对比图。
图8a为在马赫数为5,入射角分别为10°的壁面热流对比图。
图8b为在马赫数为5,入射角分别为14°的壁面热流对比图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明,即所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供的一种基于可压缩流动的解析壁面函数的数值模拟方法,包括如下步骤:
步骤1,在壁面网格内根据可压缩流动特点建立Navier-Stockes方程;
如图2所示,在壁面网格内有如下假设成立:
(1)动量和能量的耗散在壁面法向远大于流向,因此流向的耗散可以忽略;
(2)压力梯度在平行于壁面方向是常数;
(3)在可压缩流动中粘性耗散项相较于其它项不可忽略;
因此,所述步骤1中,在壁面网格内根据可压缩流动特点建立的Navier-Stockes方程为:
步骤2,简化步骤1建立的Navier-Stockes方程,得到壁面动量方程和壁面能量方程;
通过边界层假设,可得到所述步骤2中简化步骤1建立的Navier-Stockes方程,得到的壁面动量方程和壁面能量方程:
1)壁面动量方程:
2)壁面能量方程:
其中,U、V分别表示沿x和y方向速度,p表示压力,μ表示层流粘性系数,μt表示湍流粘性系数;T表示温度。
其中,可以通过DNS数值实验观测总结,得到所述层流粘性系数μ和湍流粘性系数μt:
1)如图3所示,层流粘性系数μ:
当yn>y>yv时,μ=μv;
2)如图4所示,湍流粘性系数μt:
当y<yv时,μt=0;
步骤3,通过对壁面动量方程两次积分得到全湍流区速度方程和粘性底层速度方程;并通过对壁面能量方程两次积分得到全湍流区温度方程和粘性底层温度方程;其中,在进行该积分时,对流项和压力梯度使用网格点P的值作为常数;
其中的主要参数为:
其中的主要参数为:
步骤4,根据全湍流区速度方程和粘性底层速度方程,以及全湍流区温度方程和粘性底层温度方程,定义应力方程和热流方程;
1)应力方程为:
2)热流方程为:
步骤5,给定粘性底层温度Tv的初值;具体地,给定的粘性底层温度Tv的初值为壁面温度TW。
步骤6,根据粘性底层温度Tv计算壁面网格中粘性底层的密度和层流粘性系数;
其中,气体常数R=287.06m2S-2K-1。
步骤7,将步骤6得到的粘性底层的密度和粘性系数代入应力方程和热流方程,计算得到应力和热流;并利用计算得到的应力和热流更新步骤1建立的Navier-Stockes方程中的应力项和热流项;具体地,将计算得到的应力和热流作为源项,通过添加源项的方式更新步骤1建立的Navier-Stockes方程中的应力项和热流项;其中,
1)计算得到应力的源项为:
2)计算得到热流的源项为:
步骤8,通过解析全湍流区速度方程,计算湍动能生成项和湍动能耗散项的平均量;并利用湍动能生成项和湍动能耗散项的平均量更新湍动能方程中的生成项和耗散项的值;
1)湍动能生成项的平均量为:
2)如图5所示,在壁面网格内重新选取位置yd,则壁面网格内的湍动能耗散项分布为:
则湍动能耗散项的平均量为:
所述湍动能方程为:
其中,σk=1。
在步骤8中,可以通过将湍动能生成项和湍动能耗散项的平均量直接替换湍动能方程中的生成项和耗散项的值的方式进行更新。
步骤9,利用粘性底层温度方程重新计算粘性底层温度Tv,在下个时间步循环中重复执行步骤6~9,形成闭环的数值模拟回路。
分别比较本发明基于可压缩流动的解析壁面函数的数值模拟方法(AWF),现有的标准避免函数(SWF),低雷诺数湍流模型(LS,Launder and Sharma k-ε模型),以及风洞实验值(EXP)在马赫数为5,入射角分别为10°和14°的斜激波边界层干扰的壁面压力、摩擦系数和壁面热流。其中,本发明(AWF)和现有的标准避免函数(SWF)使用的网格为120×45,壁面的第一层网格y+=20;低雷诺数湍流模型(LS,Launder and Sharma k-ε模型)使用的网格为240×80。则有:
如图6a所示的在马赫数为5,入射角分别为10°的壁面压力对比图;
如图6b所示的在马赫数为5,入射角分别为14°的壁面压力对比图;
如图7a所示的在马赫数为5,入射角分别为10°的摩擦系数对比图;
如图7b所示的在马赫数为5,入射角分别为14°的摩擦系数对比图;
如图8a所示的在马赫数为5,入射角分别为10°的壁面热流对比图;
如图8b所示的在马赫数为5,入射角分别为14°的壁面热流对比图;
通过以上对比可知,使用粗网格的本发明的基于可压缩流动的解析壁面函数的数值模拟方法,可以准确的模拟壁面流动,壁面附近的压力、摩擦系数和壁面热流的结果接近用使用密网格的低雷诺数湍流模型,而使用粗网格的现有的标准壁面函数却不能准确预测壁面热流。因此,本发明的基于可压缩流动的解析壁面函数的数值模拟方法,基于可压缩流动特点,特别适合于高超声速流动,由于其壁面动量方程和壁面能量方程之间的耦合更加紧密,可以更加准确的预测壁面热流。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种基于可压缩流动的解析壁面函数的数值模拟方法,其特征在于,包括如下步骤:
步骤1,在壁面网格内根据可压缩流动特点建立Navier-Stockes方程;
步骤2,简化步骤1建立的Navier-Stockes方程,得到壁面动量方程和壁面能量方程;
步骤3,通过对壁面动量方程两次积分得到全湍流区速度方程和粘性底层速度方程;并通过对壁面能量方程两次积分得到全湍流区温度方程和粘性底层温度方程;
步骤4,根据全湍流区速度方程和粘性底层速度方程,以及全湍流区温度方程和粘性底层温度方程,定义应力方程和热流方程;
步骤5,给定粘性底层温度Tv的初值;
步骤6,根据粘性底层温度Tv计算壁面网格中粘性底层的密度和层流粘性系数;
步骤7,将步骤6得到的粘性底层的密度和粘性系数代入应力方程和热流方程,计算得到应力和热流;并利用计算得到的应力和热流更新步骤1建立的Navier-Stockes方程中的应力项和热流项;
步骤8,通过解析全湍流区速度方程,计算湍动能生成项和湍动能耗散项的平均量;并利用湍动能生成项和湍动能耗散项的平均量更新湍动能方程中的生成项和耗散项的值;
步骤9,利用粘性底层温度方程重新计算粘性底层温度Tv,在下个时间步循环中重复执行步骤6~9,形成闭环的数值模拟回路;
所述步骤1建立的Navier-Stockes方程为:
所述步骤2中,简化步骤1建立的Navier-Stockes方程,得到的壁面动量方程和壁面能量方程分别为:
1)壁面动量方程:
2)壁面能量方程:
其中,U、V分别表示沿x和y方向速度,p表示压力,μ表示层流粘性系数,μt表示湍流粘性系数;T表示温度;
所述层流粘性系数μ和湍流粘性系数μt分别为:
1)层流粘性系数μ:
当yn>y>yv时,μ=μv;
2)湍流粘性系数μt:
当y<yv时,μt=0;
4.根据权利要求3所述的基于可压缩流动的解析壁面函数的数值模拟方法,其特征在于,所述步骤5中给定的粘性底层温度Tv的初值为壁面温度TW。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910705306.XA CN110489709B (zh) | 2019-08-01 | 2019-08-01 | 基于可压缩流动的解析壁面函数的数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910705306.XA CN110489709B (zh) | 2019-08-01 | 2019-08-01 | 基于可压缩流动的解析壁面函数的数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110489709A CN110489709A (zh) | 2019-11-22 |
CN110489709B true CN110489709B (zh) | 2022-06-10 |
Family
ID=68548959
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910705306.XA Active CN110489709B (zh) | 2019-08-01 | 2019-08-01 | 基于可压缩流动的解析壁面函数的数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110489709B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111523220B (zh) * | 2020-04-17 | 2023-03-21 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种考虑流动影响的风扇、压气机转静干涉宽频噪声快速预测方法 |
CN112052512B (zh) * | 2020-07-23 | 2023-01-10 | 中国空气动力研究与发展中心计算空气动力研究所 | 湍流边界层分层判据的方法 |
CN112069463B (zh) * | 2020-08-28 | 2022-06-03 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种风洞群高压空气资源消耗的计算方法 |
CN112560358B (zh) * | 2020-12-10 | 2022-06-14 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于径向基函数的dsmc计算结果降噪处理方法 |
CN113051846B (zh) * | 2021-03-15 | 2022-05-17 | 中国空气动力研究与发展中心计算空气动力研究所 | 考虑可压缩和热传导效应的壁面第一层网格厚度估算方法 |
CN113158338B (zh) * | 2021-04-13 | 2022-05-03 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于粗网格快速湍流壁面函数气动力预测方法 |
CN113158340B (zh) * | 2021-04-16 | 2022-10-18 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种针对k-epsilon湍流模型的湍流长度尺度修正方法 |
CN113065201B (zh) * | 2021-05-08 | 2022-03-18 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种考虑滑移修正的辐射平衡温度计算方法 |
CN113343596B (zh) * | 2021-05-25 | 2022-07-08 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于Launder-Sharma k-epsilon模型的可压缩修正方法 |
CN114139465B (zh) * | 2021-10-25 | 2023-04-07 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种脱体涡模拟模型构造方法 |
CN116070071B (zh) * | 2023-04-04 | 2023-06-27 | 中国空气动力研究与发展中心计算空气动力研究所 | 壁面函数的选择方法、装置、终端设备及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012177625A (ja) * | 2011-02-25 | 2012-09-13 | Honda Motor Co Ltd | 圧縮性流体拡散予測装置及び方法 |
CN107292001A (zh) * | 2017-06-06 | 2017-10-24 | 北京航空航天大学 | 考虑边界层燃烧放热效应的可压缩壁函数计算方法 |
CN107742050A (zh) * | 2017-11-14 | 2018-02-27 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种面向混合网格高精度阻力预测的熵修正方法 |
CN108763800A (zh) * | 2018-06-04 | 2018-11-06 | 北京理工大学 | 一种空化可压缩流动激波动力学数值模拟方法 |
CN108984874A (zh) * | 2018-07-02 | 2018-12-11 | 中国航空发动机研究院 | 获取不可压缩流动的流场的数值模拟方法 |
CN109859311A (zh) * | 2019-01-29 | 2019-06-07 | 河海大学 | 一种基于Liutex-Omega涡识别理论的空化流动数值模拟方法 |
-
2019
- 2019-08-01 CN CN201910705306.XA patent/CN110489709B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012177625A (ja) * | 2011-02-25 | 2012-09-13 | Honda Motor Co Ltd | 圧縮性流体拡散予測装置及び方法 |
CN107292001A (zh) * | 2017-06-06 | 2017-10-24 | 北京航空航天大学 | 考虑边界层燃烧放热效应的可压缩壁函数计算方法 |
CN107742050A (zh) * | 2017-11-14 | 2018-02-27 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种面向混合网格高精度阻力预测的熵修正方法 |
CN108763800A (zh) * | 2018-06-04 | 2018-11-06 | 北京理工大学 | 一种空化可压缩流动激波动力学数值模拟方法 |
CN108984874A (zh) * | 2018-07-02 | 2018-12-11 | 中国航空发动机研究院 | 获取不可压缩流动的流场的数值模拟方法 |
CN109859311A (zh) * | 2019-01-29 | 2019-06-07 | 河海大学 | 一种基于Liutex-Omega涡识别理论的空化流动数值模拟方法 |
Non-Patent Citations (2)
Title |
---|
Generalized wall function and its application to compressible turbulent boundary layer over a flat plate;J Liu 等;《Journal of Physics: Conference Series》;20170131;第1-7页 * |
高速可压缩流动壁函数边界条件的改进与应用;高振勋 等;《中国科学》;20131231;第43卷(第12期);第1359-1373页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110489709A (zh) | 2019-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110489709B (zh) | 基于可压缩流动的解析壁面函数的数值模拟方法 | |
Dewan | Tackling turbulent flows in engineering | |
Aupoix | Roughness corrections for the k–ω shear stress transport model: Status and proposals | |
Walters et al. | A three-equation eddy-viscosity model for Reynolds-averaged Navier–Stokes simulations of transitional flow | |
Tafti | Evaluating the role of subgrid stress modeling in a ribbed duct for the internal cooling of turbine blades | |
Alinot et al. | k‐ϵ model for the atmospheric boundary layer under various thermal stratifications | |
Chang et al. | Unsteady numerical simulations of turbulence and coherent structures in axial flow near a narrow gap | |
Abdel-Wahab et al. | Large eddy simulation of flow and heat transfer in a 90 deg ribbed duct with rotation: Effect of Coriolis and centrifugal buoyancy forces | |
CN107832260B (zh) | 一种平板冲击射流传热问题的数值模拟方法 | |
Issa et al. | Pressure-based compressible calculation method utilizing total variation diminishing schemes | |
Sun et al. | Evaluation of macroscopic properties in the direct simulation Monte Carlo method | |
Jammalamadaka et al. | Numerical investigations of shock wave interactions with a supersonic turbulent boundary layer | |
Raman et al. | Effect of axis ratio on fluid flow around an elliptic cylinder—a numerical study | |
Kanani et al. | Simulations of slot film-cooling with freestream acceleration and turbulence | |
Ganatra et al. | Comparison of various RANS models for impinging round jet cooling from a cylinder | |
Sewall et al. | Large Eddy Simulation of Flow and Heat Transfer in the 180‐Deg Bend Region of a Stationary Gas Turbine Blade Ribbed Internal Cooling Duct | |
Han et al. | A wall function approach in lattice Boltzmann method: algorithm and validation using turbulent channel flow | |
Zangeneh | Development of a new algorithm for modeling viscous transonic flow on unstructured grids at high Reynolds numbers | |
Zhalehrajabi et al. | Effects of mesh grid and turbulence models on heat transfer coefficient in a convergent–divergent nozzle | |
Musa et al. | Verification study of a cfd-rans code for turbulent flow at high reynolds numbers | |
Goldberg et al. | Rough surface wall function treatment with single equation turbulence models | |
AMERI et al. | Navier-Stokes turbine heat transfer predictions using two-equation turbulence closures | |
Ilinca et al. | Adaptive finite element solution of compressible turbulent flows | |
Kumar Singh et al. | Analysis of Turbulent Natural and Mixed Convection Flows Using the v 2–f Model | |
Atik | Estimation of Discretization Uncertainty Using the γ− Re θ Transition Model for Transitional Flows on 6: 1 Spheroid |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |