CN112052512A - 湍流边界层分层判据的方法 - Google Patents
湍流边界层分层判据的方法 Download PDFInfo
- Publication number
- CN112052512A CN112052512A CN202010715245.8A CN202010715245A CN112052512A CN 112052512 A CN112052512 A CN 112052512A CN 202010715245 A CN202010715245 A CN 202010715245A CN 112052512 A CN112052512 A CN 112052512A
- Authority
- CN
- China
- Prior art keywords
- boundary layer
- dimensionless
- turbulent
- turbulent boundary
- formula
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
Description
技术领域
本发明属于湍流边界分析技术领域,特别涉及一种湍流边界层分层判据的方法。
背景技术
物理风洞试验是预估各类航天航空飞行器气动性能的重要手段,随着计算流体力学的发展和计算机资源的日益丰富,数值模拟技术在飞行器气动设计中发挥着越来越重要的作用,尤其是涉及到时空多尺度流动结构的湍流流动中。尽管近几十年来计算机资源飞速增长,巨大的计算资源消耗仍然是限制直接数值模拟DNS和大涡模拟LES在工程复杂湍流问题上应用的一个根本因素。即使是采用计算资源消耗最少的雷诺平均湍流模型模拟RANS方法,要获得能够满足工程设计需求精度的比较准确的壁面摩擦阻力和壁面热流,通常需要在湍流边界层中最靠近固体壁面的粘性底层中布置相当数量的网格分布,从而使得壁面附近的网格十分细密,不仅会极大地增加计算过程迭代收敛的步数,还会因壁面附近较大的网格长细比而带来较为严重的数值刚性问题,导致计算过程的稳定性大幅下降。对湍流边界层采用壁面函数技术来模化处理,可以大幅放宽壁面法向第一层网格的尺度,不再要求第一层网格位于粘性底层内,可以使其位于湍流脉动已充分发展的对数律层内。湍流边界层壁面函数技术的应用一般情况下可以将壁面法向第一层网格的尺度放大两到三个数量级。一方面,最细密网格尺度的增加能够加速计算迭代收敛的过程;另一方面,壁面附近网格长细比的降低能明显改善数值刚性问题,有利于增强计算过程的稳定性。
采用壁面函数技术模拟湍流边界层时,首先需要对湍流边界层进行分层,然后对各层内的流动分别采用不同的函数进行近似模化处理。湍流边界层分层的主要依据是不同位置的流动特征各不相同。湍流边界层分为内层和外层,通常情况下,内层又可细分为粘性底层、过渡层和对数律层,外层一般可分为速度亏损层和粘性上层。湍流边界层的壁面函数技术主要研究湍流边界层内层中的流动分布情况。
湍流边界层分层的主要依据是不同层中流动特征各不相同,其具体表现为:粘性底层紧挨于固体壁面,由于受到固体壁面的约束限制,湍流脉动小到可以忽略,流动粘性中分子粘性占主导,湍流涡粘性可以忽略;在对数律层中,湍流脉动充分发展,流动粘性由湍流涡粘性主导,分子粘性效应可以忽略;在粘性底层和对数律层之间的过渡层,分子粘性和湍流涡粘性同等重要,二者共同影响流动分布。
现有技术对湍流边界层内层进行分层的主要依据是无量纲的壁面距离+的取值范围,且一般认为y+≤5为粘性底层,5<y+<30为过渡层,y+≥30则为对数律层
尽管目前上述依据无量纲壁面距离y+对湍流边界层进行分层的处理方法得到较为广泛的应用,但其存在的问题也是显而易见的,主要有:
1、现有技术的无量纲壁面距离+与局部流动特征无直接关系,不能完全体现湍流边界层内的湍流脉动特征;
2、现有技术的基于无量纲壁面距离+的分层判据的确定主要依赖于经验,缺乏严格的理论支撑,其分层标准在具体实施过程中依赖于使用者的经验,变化范围较大,影响了对湍流边界层的模拟精度;
3、在流向压力梯度不可忽略时,湍流边界层在壁面法向上的流动分布会发生显著变化,而现有无量纲壁面距离y+判据并不能反映这些变化特征。
这些问题的存在增加了壁面函数方法在湍流边界层中应用的不确定性,也使得目前壁面函数方法的应用效果严重依赖于使用者的经验,亟待改进分层判据以适应不同流动特征的湍流边界层,以更好地应用壁面函数方法。
因此,目前迫切需要本领域技术人员解决的一个技术问题就是:如何能够创新的提出一种有效的湍流边界层分层判据的方法,以完全体现湍流边界层内的湍流脉动特征,有效提高湍流边界层的模拟精度,克服现有技术不适用于存在流向压力梯度的流动问题。
发明内容
本申请实施例的目的是提出一种湍流边界层分层判据的方法,以解决现有技术的不能完全体现湍流边界层内的湍流脉动特征、对湍流边界层的模拟精度不高、不适用于存在流向压力梯度的流动的问题。
为解决上述技术问题,本申请实施例是这样实现的:
根据本申请实施例的一个方面,提出了一种湍流边界层分层判据的方法,包括:
依据无量纲粘性系数的取值范围对湍流边界层进行分层。
基于本发明的湍流边界层分层判据的方法的另一个实施例中,所述依据无量纲粘性系数的取值范围对湍流边界层进行分层包括:
当无量纲粘性系数满足μ+<0.1时,为粘性底层;
当无量纲粘性系数满足0.1≤μ+≤10时,为过渡层;
当无量纲粘性系数满足μ+>10时,为对数律层。
对于不可压缩流动,忽略密度变化的影响,则:
式中,uc=uτ+up,其中,ρ为流体密度,y为壁面距离,μL为流动分子粘性系数,uc为uτ和up之和,为基于速度uc的无量纲壁面距离,uτ为摩擦速度,up为压力梯度引起的速度量,AT是相对于的倒数相关量;Ap是相对于的倒数相关量,和分别为基于uτ和up的无量纲壁面距离。
基于本发明的湍流边界层分层判据的方法的另一个实施例中,所述摩擦速度uτ、压力梯度引起的速度量up的算法公式为:
基于本发明的湍流边界层分层判据的方法的另一个实施例中,所述分段函数f1的算法公式为:
式中,a1、a2、a3、b0、b1、b2、b3、c0、c1、c2、c3、c4为多项式拟合系数,为常数,k、B为常数。
基于本发明的湍流边界层分层判据的方法的另一个实施例中,所述分段函数f2的算法公式为:
式中,a2、a3、b0、b1、b2、b3、b4、c0、c1、c2、c3、c4为多项式拟合系数,为常数,α、β为常数。
本发明给出的无量纲粘性系数与当地的湍流脉动特征直接相关,能够完全体现湍流边界层内的湍流脉动特征,分层标准确定严格,实施过程不依赖于使用者的经验,能够有效提高湍流边界层的模拟精度,可直接应用于存在压力梯度的流动中,能够有效克服现有技术的湍流边界层分层判据不适用于存在流向压力梯度的流动问题。
附图说明
图1示出了根据本申请的一个实施例提供的湍流边界层分层判据的方法的流程图。
具体实施方式
以下结合附图对本发明的具体实施方式作出详细说明。
图1示出了根据本申请的一些实施例提供的湍流边界层分层判据的方法的流程图,结合图1对示例实施例中的湍流边界层分层判据的方法进行详细的说明。所述湍流边界层分层判据的方法包括:
20,依据无量纲粘性系数的取值范围对湍流边界层进行分层。
所述依据无量纲粘性系数的取值范围对湍流边界层进行分层包括:
当无量纲粘性系数满足μ+<0.1时,为粘性底层;
当无量纲粘性系数满足0.1≤μ+≤10时,为过渡层;
当无量纲粘性系数满足μ+>10时,为对数律层。
对于不可压缩流动,忽略密度变化的影响,则:
式中,uc=uτ+up,其中,ρ为流体密度,y为壁面距离,μL为流动分子粘性系数,uc为uτ和up之和,为基于速度uc的无量纲壁面距离,uτ为摩擦速度,up为压力梯度引起的速度量。AT是相对于的倒数相关量;Ap是相对于的倒数相关量,和分别为基于uτ和up的无量纲壁面距离。
所述摩擦速度uτ、压力梯度引起的速度量up的算法公式为:
所述分段函数f1的算法公式为:
式中,a1、a2、a3、b0、b1、b2、b3、c0、c1、c2、c3、c4为多项式拟合系数,为常数,k、B为常数。
多项式系数a1、a2、a3、b0、b1、b2、b3、c0、c1、c2、c3、c4的取值如下表所示:
a<sub>1</sub> | a<sub>2</sub> | a<sub>3</sub> | ||
1.0 | 1.0E-02 | -2.9E-03 | ||
b<sub>0</sub> | b<sub>1</sub> | b<sub>2</sub> | b<sub>3</sub> | b<sub>4</sub> |
-0.872 | 1.465 | -7.02E-02 | 1.66E-03 | -1.495E-05 |
c<sub>0</sub> | c<sub>1</sub> | c<sub>2</sub> | c<sub>3</sub> | c<sub>4</sub> |
8.6 | 0.1864 | -2.006E-03 | 1.144E-05 | -2.551E-08 |
所述分段函数f2的算法公式为:
式中,a2、a3、b0、b1、b2、b3、b4、c0、c1、c2、c3、c4为多项式拟合系数,为常数,α、β为常数。
多项式系数a2、a3、b0、b1、b2、b3、b4、c0、c1、c2、c3、c4的取值如下表所示:
a<sub>2</sub> | a<sub>3</sub> | |||
0.5 | -7.31E-03 | |||
b<sub>0</sub> | b<sub>1</sub> | b<sub>2</sub> | b<sub>3</sub> | b<sub>4</sub> |
-15.138 | 8.4688 | -0.81976 | 3.7292E-02 | -6.3866E-04 |
c<sub>0</sub> | c<sub>1</sub> | c<sub>2</sub> | c<sub>3</sub> | c<sub>4</sub> |
11.925 | 0.934 | -2.7805E-02 | 4.6262E-04 | -3.1442E-06 |
对于本领域技术人员而言,显然本发明实施例不限于上述示范性实施例的细节,而且在不背离本发明实施例的精神或基本特征的情况下,能够以其他的具体形式实现本发明实施例。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明实施例的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化涵括在本发明实施例内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。此外,显然“包括”一词不排除其他单元或步骤,单数不排除复数。系统、装置或终端权利要求中陈述的多个单元、模块或装置也可以由同一个单元、模块或装置通过软件或者硬件来实现。第一,第二等词语用来表示名称,而并不表示任何特定的顺序。
最后应说明的是,以上实施方式仅用以说明本发明实施例的技术方案而非限制,尽管参照以上较佳实施方式对本发明实施例进行了详细说明,本领域的普通技术人员应当理解,可以对本发明实施例的技术方案进行修改或等同替换都不应脱离本发明实施例的技术方案的精神和范围。
Claims (8)
2.根据权利要求1所述的湍流边界层分层判据的方法,其特征在于,所述依据无量纲粘性系数的取值范围对湍流边界层进行分层包括:
当无量纲粘性系数满足μ+<0.1时,为粘性底层;
当无量纲粘性系数满足0.1≤μ+≤10时,为过渡层;
当无量纲粘性系数满足μ+>10时,为对数律层。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010715245.8A CN112052512B (zh) | 2020-07-23 | 2020-07-23 | 湍流边界层分层判据的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010715245.8A CN112052512B (zh) | 2020-07-23 | 2020-07-23 | 湍流边界层分层判据的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112052512A true CN112052512A (zh) | 2020-12-08 |
CN112052512B CN112052512B (zh) | 2023-01-10 |
Family
ID=73601028
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010715245.8A Active CN112052512B (zh) | 2020-07-23 | 2020-07-23 | 湍流边界层分层判据的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112052512B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114004175A (zh) * | 2021-10-29 | 2022-02-01 | 中船重工奥蓝托无锡软件技术有限公司 | 一种快速查找全域壁面距离和无量纲壁面距离的方法 |
CN116070071A (zh) * | 2023-04-04 | 2023-05-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 壁面函数的选择方法、装置、终端设备及存储介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050246111A1 (en) * | 2004-03-10 | 2005-11-03 | Gysling Daniel L | Method and apparatus for measuring parameters of a stratified flow |
CN102831276A (zh) * | 2012-08-31 | 2012-12-19 | 中国能源建设集团广东省电力设计研究院 | 环境风对超大型自然通风冷却塔热力性能影响的计算方法 |
CN103884485A (zh) * | 2014-02-25 | 2014-06-25 | 国家电网公司 | 一种基于多尾流模型的风机尾流分析方法 |
US20150168600A1 (en) * | 2012-05-24 | 2015-06-18 | Lanzhou University | Monitoring system for turbulence of atmospheric boundary layer under wind drift sand flow or sand dust storm environment |
CN109446691A (zh) * | 2018-11-08 | 2019-03-08 | 南京林业大学 | 基于激光点云与空气动力学的活立木抗风性能分析方法 |
CN110309552A (zh) * | 2019-06-10 | 2019-10-08 | 中国航天空气动力技术研究院 | 一种考虑质量引射效应的飞行器湍流预测方法及系统 |
CN110489709A (zh) * | 2019-08-01 | 2019-11-22 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于可压缩流动的解析壁面函数的数值模拟方法 |
CN111079310A (zh) * | 2019-12-30 | 2020-04-28 | 空气动力学国家重点实验室 | 一种湍流区域识别方法 |
-
2020
- 2020-07-23 CN CN202010715245.8A patent/CN112052512B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050246111A1 (en) * | 2004-03-10 | 2005-11-03 | Gysling Daniel L | Method and apparatus for measuring parameters of a stratified flow |
US20150168600A1 (en) * | 2012-05-24 | 2015-06-18 | Lanzhou University | Monitoring system for turbulence of atmospheric boundary layer under wind drift sand flow or sand dust storm environment |
CN102831276A (zh) * | 2012-08-31 | 2012-12-19 | 中国能源建设集团广东省电力设计研究院 | 环境风对超大型自然通风冷却塔热力性能影响的计算方法 |
CN103884485A (zh) * | 2014-02-25 | 2014-06-25 | 国家电网公司 | 一种基于多尾流模型的风机尾流分析方法 |
CN109446691A (zh) * | 2018-11-08 | 2019-03-08 | 南京林业大学 | 基于激光点云与空气动力学的活立木抗风性能分析方法 |
CN110309552A (zh) * | 2019-06-10 | 2019-10-08 | 中国航天空气动力技术研究院 | 一种考虑质量引射效应的飞行器湍流预测方法及系统 |
CN110489709A (zh) * | 2019-08-01 | 2019-11-22 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于可压缩流动的解析壁面函数的数值模拟方法 |
CN111079310A (zh) * | 2019-12-30 | 2020-04-28 | 空气动力学国家重点实验室 | 一种湍流区域识别方法 |
Non-Patent Citations (4)
Title |
---|
SINA KHANI等: "Effective eddy viscosity in stratified turbulence", 《JOURNAL OF TURBULENCE》 * |
周根水等: "分层流体中尾流混合效应数值模拟", 《水动力学研究与进展(A辑)》 * |
李兆辉: "无平均剪切稳定分层二层流非湍流/湍流密度界面处的湍流", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》 * |
邵文毅: "非均匀湍流路径下大气分层特性研究", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114004175A (zh) * | 2021-10-29 | 2022-02-01 | 中船重工奥蓝托无锡软件技术有限公司 | 一种快速查找全域壁面距离和无量纲壁面距离的方法 |
CN114004175B (zh) * | 2021-10-29 | 2023-09-22 | 中船奥蓝托无锡软件技术有限公司 | 一种快速查找全域壁面距离和无量纲壁面距离的方法 |
CN116070071A (zh) * | 2023-04-04 | 2023-05-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 壁面函数的选择方法、装置、终端设备及存储介质 |
CN116070071B (zh) * | 2023-04-04 | 2023-06-27 | 中国空气动力研究与发展中心计算空气动力研究所 | 壁面函数的选择方法、装置、终端设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112052512B (zh) | 2023-01-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Frink | Tetrahedral unstructured Navier-Stokes method for turbulent flows | |
Coles et al. | Flying-hot-wire study of flow past an NACA 4412 airfoil at maximum lift | |
Coder et al. | Comparisons of theoretical methods for predicting airfoil aerodynamic characteristics | |
CN112052512B (zh) | 湍流边界层分层判据的方法 | |
US20120232860A1 (en) | Generating a simulated fluid flow over an aircraft surface using anisotropic diffusion | |
Lörstad et al. | Assessment of volume of fluid and immersed boundary methods for droplet computations | |
Rumsey et al. | HLPW-4/GMGW-3: overview and workshop summary | |
Ranzenbach et al. | Cambered airfoil in ground effect-an experimental and computational study | |
Shrivastava et al. | A novel level set-based immersed-boundary method for CFD simulation of moving-boundary problems | |
Araújo et al. | Experimental evaluation of 2-D entropy model for open-channel flow | |
Jin et al. | Improvement of fast fluid dynamics with a conservative semi-Lagrangian scheme | |
Ramnefors et al. | Accuracy of drag predictions on cars using CFD-effect of grid refinement and turbulence models | |
Ursachi et al. | Output-based adaptive RANS solutions using higher-order FEM on a multi-element airfoil | |
Goc et al. | Wind tunnel and grid resolution effects in large-eddy simulations of the high-lift common research model | |
Li et al. | Numerical and experimental investigations on drag-reducing effects of riblets | |
Abdolmaleki et al. | Performance analysis of elliptic-profile airfoil cascade for designing reversible axial flow fans | |
CN105372579A (zh) | 一种快速有效的电路单元重要性测度方法 | |
Chang et al. | Effects of Reynolds number and protuberance amplitude on twin-protuberance airfoil performance | |
CN112417785A (zh) | 基于微纳米沟槽壁面滑移效应的跨尺度数值模拟方法 | |
Sclafani et al. | OVERFLOW analysis of the NASA trap wing model from the first high lift prediction workshop | |
De La Calzada et al. | Heat transfer in separated flows on the pressure side of turbine blades | |
Liu et al. | Performance and inter-blade flow of axial flow fans with different blade angles of attack | |
Leishman | Unsteady aerodynamics of airfoils encountering traveling gusts and vortices | |
Xu et al. | On the development of turbomachine blade aerodynamic design system | |
Zastawny et al. | Validation of intermittency-based transition prediction models |
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 |