CN114970403A - 一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法 - Google Patents

一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法 Download PDF

Info

Publication number
CN114970403A
CN114970403A CN202210816721.4A CN202210816721A CN114970403A CN 114970403 A CN114970403 A CN 114970403A CN 202210816721 A CN202210816721 A CN 202210816721A CN 114970403 A CN114970403 A CN 114970403A
Authority
CN
China
Prior art keywords
turbulence
series model
function
gas
grid
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
Application number
CN202210816721.4A
Other languages
English (en)
Inventor
柳阳威
汪光宇
唐雨萌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN202210816721.4A priority Critical patent/CN114970403A/zh
Publication of CN114970403A publication Critical patent/CN114970403A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • 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
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • 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)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于湍流能谱耦合v2‑f系列模型的网格自适应湍流模拟方法，其实现的主要过程包括：步骤一，判断是否应用屏蔽函数；步骤二，识别当地网格尺度；步骤三，基于湍流能谱积分耦合v2‑f系列模型构造尺度相关的调节函数；步骤四，使用调节函数重构v2‑f系列模型的湍流粘性；步骤五，基于v2‑f系列模型，使用重构的湍流粘性进行湍流模拟。本发明通过识别当地网格大小，确定当地网格长度尺度Δ*，进而通过湍流能谱积分构造尺度相关函数对湍流粘性进行重构，实现网格自适应模拟，有效克服了现有RANS/LES混合模型对网格的经验依赖度高的问题，大幅减少计算耗费，显著加快湍流模拟进程，为解决复杂工程流动问题的快速高精度模拟提供了高效的数值模拟方法。

Description

一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟 方法
技术领域
本发明涉及工程流体力学计算领域，尤其涉及一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法。
背景技术
湍流现象在自然界和工程领域中普遍存在，对湍流的精准预测是研究复杂流动问题的一大难题。以流体机械内部流动为例，其流动结构十分复杂，广泛存在多尺度、非线性、非定常等复杂湍流流动现象，这些复杂湍流对流体机械的性能有很大影响。因此，在工程设计中，迫切的需要发展预测精度高、计算效率高的湍流模拟方法，使之能准确模拟以流体机械内部流动为例的复杂工程问题中的湍流流动。
现有工程常用的湍流模拟方法以求解雷诺平均NS方程(RANS)方法为主，虽然计算量较小，但由于其基于简单基本流动发展而来，对工程问题中多尺度、非定常、大分离等复杂湍流流动预测准确性不佳，从而难以实现对以流体机械内部流动为例的工程问题中复杂流动机理的精细化研究，严重制约流体机械等领域设计水平的提升。大涡模拟(LES)方法作为高精度数值模拟方法，对网格数的要求很高，使其相比于RANS方法计算耗费呈指数级增长，远高于工程应用中可以承受的计算耗费水平。对于流体机械内部流动等工程流动问题，雷诺数往往较高，在现有计算能力下，LES方法还难以在工程设计层面上应用于复杂工程领域的流动预测。
近二十年来诞生了RANS-LES混合模拟方法，通过在近壁区域采用RANS方法，主流区域采用LES方法，来平衡计算精度和计算效率，对解决复杂流动的高精度模拟计算耗费大的问题，提供了很好的解决策略。但目前经典的RANS-LES混合模型对网格有较为严格的要求，需要使用者具有较为丰富的高精度数值模拟经验，因而对以流体机械内部流动为例的复杂工程流动的预测能力有限。因此，降低RANS-LES混合模型对网格的经验依赖程度，对更好地实现复杂工程流动问题中多尺度、非线性、非定常等湍流流动现象的准确、高效的高精度模拟具有重要意义。
发明内容
(一)要解决的技术问题
本发明的目的在于提出一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法，有效克服现有RANS-LES混合模型对网格的经验依赖度高的问题，在提升计算准确度的同时，大幅减少计算耗费，显著加快湍流模拟进程，为解决复杂工程流动问题中多尺度、非线性、非定常等湍流流动的快速高精度模拟提供高效的数值模拟方法。
(二)技术方案
为了解决上述技术问题，本发明提供一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法，包括以下步骤：
步骤一，判断是否应用屏蔽函数；
步骤二，识别当地网格尺度；
步骤三，基于湍流能谱积分耦合v2-f系列模型构造尺度相关的调节函数；
步骤四，使用调节函数重构v2-f系列模型的湍流粘性；
步骤五，基于v2-f系列模型，使用重构的湍流粘性进行湍流模拟；
①所述判断是否应用屏蔽函数包括：
结合所模拟的流动状态类型，判断是否采用屏蔽函数FGAS，具体地，当所述流动状态类型为自由剪切流动，则不采用屏蔽函数，此时所述屏蔽函数FGAS＝0；当所述流动状态类型为近壁流动，则采用屏蔽函数，所述屏蔽函数FGAS可使用如：来源于DDES-SST模型中的F1屏蔽函数、F2屏蔽函数和来源于DDES-SA模型中的Fd屏蔽函数；
②所述识别当地网格尺度包括：
结合步骤一中的屏蔽函数FGAS，确定当地网格长度尺度Δ*；
所述当地网格长度尺度Δ*由下式给出：
Δ*＝CGAS[(1-FGAS)Δvol+FGASΔmax]
Δmax＝max(Δx,Δy,Δz)
Figure BDA0003740956500000011
其中，Δx为当地六面体网格的长、Δy为当地六面体网格的宽、Δz为当地六面体网格的高，CGAS为经验系数取0.6；
③所述基于湍流能谱积分耦合v2-f系列模型构造尺度相关的调节函数包括：
按照v2-f系列模型中对湍动能的模化方式，得原始模化的湍动能km，根据步骤二中所述当地网格长度尺度Δ*，基于湍流能谱通过积分得到实际应模化的湍动能ku；
所述实际应模化的湍动能ku由下式得到：
Figure BDA0003740956500000021
其中，Ck为柯尔莫哥洛夫常系数，取1.5，ε为实际的湍流耗散率，κc为可解湍流截断波数，由步骤二中所述当地网格长度尺度Δ*决定：
Figure BDA0003740956500000022
根据所述实际应模化的湍动能ku、所述原始模化的湍动能km和步骤一中所述屏蔽函数FGAS，构造动态尺度相关的调节函数Df；定义尺度之比为所述实际应模化的湍动能ku和所述原始模化的湍动能km的比值，所述动态尺度相关的调节函数Df为所述尺度之比相关的函数，由下式得出：
Figure BDA0003740956500000023
lGAS＝(1-FGAS)lu+FGASlm
Figure BDA0003740956500000024
其中，lu为网格相关尺度，lm为v2-f系列模型给出的湍流长度尺度；
④所述使用调节函数重构v2-f系列模型的湍流粘性包括：
采用步骤三所述动态尺度相关的调节函数Df对v2-f系列模型的湍流粘性νt进行调控，得到重构的湍流粘性νsfs，由下式得出：
νsfs＝Df·νt
⑤所述基于v2-f系列模型，使用重构的湍流粘性进行湍流模拟包括：
采用步骤四中所述重构的湍流粘性νsfs计算雷诺应力，并对v2-f系列模型的输运方程进行更新，替代v2-f系列模型的原始输运方程中的湍流粘性νt，与v2-f系列模型结合得到所述基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法。
(三)有益效果
本发明提供了一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法，具有以下有益效果：通过识别当地网格大小，确定当地网格长度尺度Δ*，进而通过湍流能谱积分构造尺度相关函数对湍流粘性进行重构，实现网格自适应模拟，有效克服了现有RANS-LES混合模型对网格的经验依赖度高的问题，在提升计算准确度的同时，大幅减少计算耗费，显著加快湍流模拟进程，为解决复杂工程流动问题中多尺度、非线性、非定常等湍流流动的快速高精度模拟提供高效的数值模拟方法；本发明方法形式简明，可移植性强，能够很好地与v2-f系列模型进行结合，便于植入现有CFD代码进行应用和拓展，具有广阔的学术、工程应用前景。
附图说明
图1为本发明基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法的流程图；
图2本发明基于湍流能谱构造动态尺度相关的调节函数Df的示意图。
图3为v2-f湍流模型计算的周期山流动算例的湍流涡结构图；
图4为本发明基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法计算的周期山流动算例的湍流涡结构图。
具体实施方式
以下结合附图，以v2-f湍流模型为例，对本发明的具体实施方式进行进一步详细说明，以下实施例仅用于说明本发明，但不用于限制本发明的范围。
选取周期山流动算例作为计算算例，采用高质量六面体网格对计算域进行空间离散，总网格数约为22万。
本发明提供一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法，包括如下步骤：
步骤一，判断是否应用屏蔽函数；
在这一步骤中，结合所模拟的流动状态类型，判断是否采用屏蔽函数FGAS，具体地，当流动状态类型为自由剪切流动，则不采用屏蔽函数，此时屏蔽函数FGAS＝0；当流动状态类型为近壁流动，则采用屏蔽函数，所述屏蔽函数FGAS可使用如：来源于DDES-SST模型中的F1屏蔽函数、F2屏蔽函数和来源于DDES-SA模型中的Fd屏蔽函数。
在本实施例中，使用DDES-SA模型中的Fd屏蔽函数，则有：
FGAS＝Fd
Figure BDA0003740956500000031
Figure BDA0003740956500000032
其中，Ui,j为速度梯度张量；νt为v2-f湍流模型所给出的湍流粘性，ν为流体的粘性系数，d为当地网格到壁面的距离；κ为卡门常数，取值为0.41；Cd1、Cd2为经验系数，分别取20和3；
步骤二，识别当地网格尺度；
在这一步骤中，结合步骤一中的屏蔽函数FGAS，确定当地网格长度尺度Δ*；
当地网格长度尺度Δ*由下式给出：
Δ*＝CGAS[(1-FGAS)Δvol+FGASΔmax]
Δmax＝max(Δx,Δy,Δz)
Figure BDA0003740956500000033
其中，Δx为当地六面体网格的长、Δy为当地六面体网格的宽、Δz为当地六面体网格的高，CGAS为经验系数，取0.6；
步骤三，基于湍流能谱积分耦合v2-f系列模型构造尺度相关的调节函数：
在这一步骤中，按照v2-f湍流模型中对湍动能的模化方式，得原始模化的湍动能km；
根据步骤二中的当地网格长度尺度Δ*，基于湍流能谱通过积分得到实际应模化的湍动能ku，由下式给出：
Figure BDA0003740956500000034
其中，Ck为柯尔莫哥洛夫常系数，取1.5；ε为实际的湍流耗散率；κc为可解湍流截断波数，由步骤二中的当地网格长度尺度Δ*决定：
Figure BDA0003740956500000035
根据实际应模化的湍动能ku、原始模化的湍动能km和步骤一中的屏蔽函数FGAS，构造动态尺度相关的调节函数Df，如图2所示；定义尺度之比为实际应模化的湍动能ku和原始模化的湍动能km的比值，动态尺度相关的调节函数Df为尺度之比相关的函数，由下式得出：
Figure BDA0003740956500000036
lGAS＝(1-FGAS)lu+FGASlm
Figure BDA0003740956500000037
其中，lu为网格相关尺度，lm为v2-f湍流模型给出的湍流长度尺度，由下式得出：
Figure BDA0003740956500000038
其中，km为v2-f模型给出的原始模化的湍动能，εm为v2-f模型给出的原始模化的耗散率，ν为流体的粘性系数，CL和Cη为v2-f湍流模型中的常系数，其中，CL取0.25，Cη取85；
步骤四，使用调节函数重构v2-f湍流模型的湍流粘性：
在这一步骤中，采用步骤三中动态尺度相关的调节函数Df对v2-f湍流模型所得湍流粘性νt进行调控，得到重构的湍流粘性νsfs，如下式所示：
Figure BDA0003740956500000041
其中，T为v2-f湍流模型中的湍流时间尺度，其表达式为：
Figure BDA0003740956500000042
特别地，Cμ和CT为v2-f湍流模型中的常系数，Cμ取0.22，CT取6；
步骤五，基于v2-f湍流模型使用重构的湍流粘性进行湍流模拟：
在这一步骤中，采用步骤四中重构的湍流粘性νsfs计算雷诺应力，并对v2-f湍流模型的输运方程进行更新，替代v2-f湍流模型的原始输运方程中的湍流粘性νt，所得到的新的输运方程如下式所示：
Figure BDA0003740956500000043
Figure BDA0003740956500000044
Figure BDA0003740956500000045
Figure BDA0003740956500000046
Figure BDA0003740956500000047
Figure BDA0003740956500000048
特别地，常系数C1取1.4，常系数C2取0.45，常系数
Figure BDA0003740956500000049
取1；
将所得到的新的输运方程与v2-f湍流模型结合得到基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法，并用于周期山流动算例的数值模拟。
采用全隐式耦合求解技术进行瞬态计算，时间步长满足工程计算流体力学中的CFL条件。另外，为更好地显示本发明实施例的优点，同时选取v2-f湍流模型对周期山流动算例进行数值模拟，与基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法的数值模拟结果进行对比。
图3为v2-f湍流模型计算的周期山流动算例的湍流涡结构图，采用无量纲流向速度进行着色。
图4为本发明基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法计算的周期山流动算例的湍流涡结构图，采用无量纲流向速度进行着色。
图3、图4的对比分析表明，本发明提出的基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法计算所得到的湍流涡结构解析能力强于v2-f湍流模型计算所得到的结果，在相同的网格数目下能够捕捉更为丰富的湍流结构，能够提供更加精确的流场细节。
以上所述仅为本发明专利的较佳实施例，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替代、改进等，均应包含在本发明的保护范围之内。
综上而言，本发明实施例通过识别当地网格大小，确定当地网格长度尺度Δ*，进而通过湍流能谱积分构造尺度相关函数对湍流粘性进行重构，实现网格自适应模拟，有效克服了现有湍流混合模型对网格的经验依赖度高的问题，在提升计算准确度的同时，可大幅减少计算耗费，显著加快湍流模拟进程，为复杂工程领域的流动预测提供新方法。

Claims (1)

1.一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法，其特征在于，包括以下步骤：
步骤一，判断是否应用屏蔽函数；
步骤二，识别当地网格尺度；
步骤三，基于湍流能谱积分耦合v2-f系列模型构造尺度相关的调节函数；
步骤四，使用调节函数重构v2-f系列模型的湍流粘性；
步骤五，基于v2-f系列模型，使用重构的湍流粘性进行湍流模拟；
①所述判断是否应用屏蔽函数包括：
结合所模拟的流动状态类型，判断是否采用屏蔽函数FGAS，具体地，当所述流动状态类型为自由剪切流动，则不采用屏蔽函数，此时所述屏蔽函数FGAS＝0；当所述流动状态类型为近壁流动，则采用屏蔽函数，所述屏蔽函数FGAS可使用如：来源于DDES-SST模型中的F1屏蔽函数、F2屏蔽函数和来源于DDES-SA模型中的Fd屏蔽函数；
②所述识别当地网格尺度包括：
结合步骤一中的屏蔽函数FGAS，确定当地网格长度尺度Δ*；
所述当地网格长度尺度Δ*由下式给出：
Δ*＝CGAS[(1-FGAS)Δvol+FGASΔmax]
Δmax＝max(Δx,Δy,Δz)
Figure FDA0003740956490000011
其中，Δx为当地六面体网格的长、Δy为当地六面体网格的宽、Δz为当地六面体网格的高，CGAS为经验系数取0.6；
③所述基于湍流能谱积分耦合v2-f系列模型构造尺度相关的调节函数包括：
按照v2-f系列模型中对湍动能的模化方式，得原始模化的湍动能km，根据步骤二中所述当地网格长度尺度Δ*，基于湍流能谱通过积分得到实际应模化的湍动能ku；
所述实际应模化的湍动能ku由下式得到：
Figure FDA0003740956490000012
其中，Ck为柯尔莫哥洛夫常系数，取1.5，ε为实际的湍流耗散率，κc为可解湍流截断波数，由步骤二中所述当地网格长度尺度Δ*决定：
Figure FDA0003740956490000013
根据所述实际应模化的湍动能ku、所述原始模化的湍动能km和步骤一中所述屏蔽函数FGAS，构造动态尺度相关的调节函数Df；定义尺度之比为所述实际应模化的湍动能ku和所述原始模化的湍动能km的比值，所述动态尺度相关的调节函数Df为所述尺度之比相关的函数，由下式得出：
Figure FDA0003740956490000014
lGAS＝(1-FGAS)lu+FGASlm
Figure FDA0003740956490000015
其中，lu为网格相关尺度，lm为v2-f系列模型给出的湍流长度尺度；
④所述使用调节函数重构v2-f系列模型的湍流粘性包括：
采用步骤三所述动态尺度相关的调节函数Df对v2-f系列模型的湍流粘性νt进行调控，得到重构的湍流粘性νsfs，由下式得出：
νsfs＝Df·νt
⑤所述基于v2-f系列模型，使用重构的湍流粘性进行湍流模拟包括：
采用步骤四中所述重构的湍流粘性νsfs计算雷诺应力，并对v2-f系列模型的输运方程进行更新，替代v2-f系列模型的原始输运方程中的湍流粘性νt，与v2-f系列模型结合得到所述基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法。
CN202210816721.4A 2022-07-12 2022-07-12 一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法 Pending CN114970403A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210816721.4A CN114970403A (zh) 2022-07-12 2022-07-12 一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210816721.4A CN114970403A (zh) 2022-07-12 2022-07-12 一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法

Publications (1)

Publication Number Publication Date
CN114970403A true CN114970403A (zh) 2022-08-30

Family

ID=82970175

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210816721.4A Pending CN114970403A (zh) 2022-07-12 2022-07-12 一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法

Country Status (1)

Country Link
CN (1) CN114970403A (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017004103A (ja) * 2015-06-05 2017-01-05 株式会社東芝 乱流のシミュレーション方法
CN106682398A (zh) * 2016-12-12 2017-05-17 华南理工大学 一种基于湍动能和湍流长度尺度截断的rans/les方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017004103A (ja) * 2015-06-05 2017-01-05 株式会社東芝 乱流のシミュレーション方法
CN106682398A (zh) * 2016-12-12 2017-05-17 华南理工大学 一种基于湍动能和湍流长度尺度截断的rans/les方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DONG, YIDAO 等: "An Enhanced Version of Delayed Detached-Eddy Simulation Based on the v2 − f Model", FLOW TURBULENCE COMBUST (2019), 3 July 2018 (2018-07-03), pages 167 - 188, XP036742476, DOI: 10.1007/s10494-018-9957-8 *
吴佳诺 等: "Sajben扩压器复杂流动中湍流模型的性能评估", 航空动力学报, vol. 32, no. 5, 31 May 2017 (2017-05-31), pages 1216 - 1223 *

Similar Documents

Publication Publication Date Title
CN113111608B (zh) 一种新型局部湍流脉动生成方法
Ghia et al. Analysis of incompressible massively separated viscous flows using unsteady Navier–Stokes equations
CN113158340B (zh) 一种针对k-epsilon湍流模型的湍流长度尺度修正方法
Chen et al. A machine learning based solver for pressure Poisson equations
CN112417785B (zh) 基于微纳米沟槽壁面滑移效应的跨尺度数值模拟方法
CN114970403A (zh) 一种基于湍流能谱耦合v2-f系列模型的网格自适应湍流模拟方法
Yadav et al. Analysis of retaining walls along the Beas River’s meandering curve
CN115186608A (zh) 一种基于湍流能谱耦合rsm模型的网格自适应湍流模拟方法
CN115034162B (zh) 一种基于湍流能谱耦合k-ε系列模型的网格自适应湍流模拟方法
CN115186607A (zh) 一种基于湍流能谱耦合k-ω系列模型的网格自适应湍流模拟方法
CN113158339B (zh) 一种针对sst湍流模型的湍流长度尺度修正方法
CN115186608B (zh) 一种基于湍流能谱耦合rsm模型的网格自适应湍流模拟方法
Lulekar et al. Cfd-based analysis and surrogate-based optimization of bio-inspired surface riblets for aerodynamic efficiency
CN117787130A (zh) 一种基于Vreman动态系数耦合v2-f系列模型的网格自适应湍流模拟方法
Evje et al. CFL-violating numerical schemes for a two-fluid model
CN118133701A (zh) 一种基于Vreman动态系数耦合k-ε系列模型的网格自适应湍流模拟方法
CN117763838A (zh) 一种基于湍流能谱耦合sa系列模型的网格自适应湍流模拟方法
CN117763996A (zh) 一种基于Vreman动态系数耦合RSM模型的网格自适应湍流模拟方法
Vengadesan et al. Evaluation of LES models for flow over bluff body from engineering application perspective
CN117763846A (zh) 一种基于Vreman动态系数耦合k-ω系列模型的网格自适应湍流模拟方法
CN117763995A (zh) 一种基于Vreman动态系数耦合SA系列模型的网格自适应湍流模拟方法
Skarolek et al. Transitional flow over a SD7003 wing using flux reconstruction scheme
Khokhar et al. A numerical analysis of mixing and separating of Newtonian fluids in a channel filled with porous materials
Sherar et al. An Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm.
CN113111610B (zh) 一种亚格子尺度模型建立方法

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