一种基于湍流能谱耦合SA系列模型的网格自适应湍流模拟方法

专利主权项内容

1.一种基于湍流能谱耦合SA系列模型的网格自适应湍流模拟方法，其特征在于，包括如下步骤：步骤一，判断是否应用屏蔽函数；步骤二，识别当地网格尺度；步骤三，基于湍流能谱积分耦合SA系列模型构造尺度相关的调节函数；步骤四，使用调节函数重构SA系列模型的湍流粘性；步骤五，基于SA系列模型，使用重构的湍流粘性进行湍流模拟；①所述判断是否应用屏蔽函数包括：结合所模拟的流动状态类型，判断是否采用屏蔽函数F，具体地，当所述流动状态类型为自由剪切流动，则不采用屏蔽函数，此时所述屏蔽函数F＝0；当所述流动状态类型为近壁流动，则采用屏蔽函数，所述屏蔽函数F可使用如：来源于DDES-SA模型中的F屏蔽函数，和来源于IDDES-SA模型中的F屏蔽函数；GASGASGAS<<②所述识别当地网格尺度包括：结合步骤一中所述屏蔽函数F，确定当地网格长度尺度Δ<GAS<所述当地网格长度尺度Δ由下式给出：<<＝C[(1-F)Δ+F<]<GASGASvolGASmax<＝max(Δ, Δ, Δ<max<<<其中，Δ为当地六面体网格的长，Δ为当地六面体网格的宽，Δ当地六面体网格的高，C为经验系数取0.016；<<<GAS③所述基于湍流能谱积分耦合SA系列模型构造尺度相关的调节函数包括：根据SA系列模型得到当地模化的湍流粘性ν和速度梯度张量U，由所述当地模化的湍流粘性ν和所述速度梯度张量U构造新的代数经验公式对模化的湍动能k进行表征，如下式所示：<i, j<i, jm根据步骤二中所述当地网格长度尺度Δ，基于湍流能谱通过积分得到实际应模化的湍动能k；*u所述实际应模化的湍动能k由下式得到：u其中，C为柯尔莫哥洛夫常系数，取1.5，ε为实际的湍流耗散率；κ为可解湍流截断波数，由步骤二中所述当地网格长度尺度Δ决定：kc*根据所述实际应模化的湍动能k、所述原始模化的湍动能k和步骤一中所述屏蔽函数F，构造动态尺度相关的调节函数D；定义尺度之比为所述实际应模化的湍动能k和所述原始模化的湍动能k的比值，所述动态尺度相关的调节函数D为所述尺度之比相关的函数，由下式得出：umGASfumfl＝(1-F)l+FlGASGASuGASm其中，l为网格相关尺度，l为湍流长度尺度；um所述湍流长度尺度l的表达式如下式所示：m其中，ν为SA系列模型原始湍流粘性，U为SA系列模型得到的速度梯度张量；a为经验系数，取5.794；ti, j0④所述使用调节函数重构SA系列模型的湍流粘性包括：采用步骤三所述动态尺度相关的调节函数D对SA系列模型中的湍流粘性ν进行调控，得到重构的湍流粘性ν，由下式得出：ftsfsν＝D·νsfsft⑤所述基于SA系列模型，使用重构的湍流粘性进行湍流模拟包括：采用步骤四中所述重构的湍流粘性ν，计算雷诺应力，并更新质量、动量、能量输运方程，与SA系列模型结合得到所述基于湍流能谱耦合SA系列模型的网格自适应湍流模拟方法。sfs