主权项 |
一种GAO‑YONG理性湍流模型的处理方法,其特征在于:所述GAO‑YONG理性湍流模型的机械能方程为<maths num="0001"><math><![CDATA[<mrow><mo>|</mo><mi>L</mi><mo>|</mo><mo>=</mo><mfrac><mrow><mi>A</mi><mo>-</mo><mrow><mo>(</mo><mover><mi>ρ</mi><mo>‾</mo></mover><mover><mi>U</mi><mo>~</mo></mover><mo>·</mo><mo>(</mo><mrow><mi>L</mi><mi>D</mi><mi>i</mi><mi>r</mi><mi>e</mi><mi>c</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi><mo>·</mo><mo>▿</mo><mover><mi>U</mi><mo>‾</mo></mover></mrow><mo>)</mo><mo>)</mo></mrow></mrow><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><mover><mi>ρ</mi><mo>‾</mo></mover><mover><mi>U</mi><mo>~</mo></mover><mo>·</mo><mrow><mo>(</mo><mi>L</mi><mi>D</mi><mi>i</mi><mi>r</mi><mi>e</mi><mi>c</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi><mo>·</mo><mo>(</mo><mrow><mi>L</mi><mi>D</mi><mi>i</mi><mi>r</mi><mi>e</mi><mi>c</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi><mo>·</mo><mo>▿</mo><mrow><mo>(</mo><mrow><mo>▿</mo><mover><mi>U</mi><mo>‾</mo></mover></mrow><mo>)</mo></mrow></mrow><mo>)</mo><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000895693320000011.GIF" wi="1083" he="192" /></maths>其中,<maths num="0002"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mrow><mo>(</mo><mi>L</mi><mi>D</mi><mi>i</mi><mi>r</mi><mi>e</mi><mi>c</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi><mo>·</mo><mo>(</mo><mrow><mo>-</mo><mo>▿</mo><mrow><mo>(</mo><mrow><mover><mi>U</mi><mo>‾</mo></mover><mover><mi>U</mi><mo>~</mo></mover><mo>+</mo><mover><mi>p</mi><mo>~</mo></mover><mi>I</mi></mrow><mo>)</mo></mrow><mo>+</mo><mo>▿</mo><mo>·</mo><mrow><mo>(</mo><mrow><mover><mi>σ</mi><mo>~</mo></mover><mo>+</mo><msup><mover><mi>σ</mi><mo>~</mo></mover><mrow><mi>t</mi><mi>u</mi><mi>r</mi><mi>b</mi></mrow></msup><mo>-</mo><mi>σ</mi><mo>-</mo><msup><mi>σ</mi><mrow><mi>t</mi><mi>u</mi><mi>r</mi><mi>b</mi></mrow></msup></mrow><mo>)</mo></mrow></mrow><mo>)</mo><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000895693320000012.GIF" wi="1244" he="79" /></maths>L是特征长度,运算符·表示取一阶内积,运算符<img file="FDA0000895693320000013.GIF" wi="52" he="58" />表示矢量<img file="FDA0000895693320000014.GIF" wi="306" he="160" />其中<img file="FDA0000895693320000015.GIF" wi="74" he="131" />表示对某参数求x方向偏导数,<img file="FDA0000895693320000016.GIF" wi="77" he="139" />表示对某参数求y方向偏导数,<img file="FDA0000895693320000017.GIF" wi="74" he="131" />表示对某参数求z方向偏导数,<img file="FDA0000895693320000018.GIF" wi="59" he="72" />是气流平均速度,<img file="FDA0000895693320000019.GIF" wi="58" he="89" />是气流平均密度,<img file="FDA00008956933200000110.GIF" wi="54" he="74" />是气流脉动速度,<img file="FDA00008956933200000111.GIF" wi="49" he="74" />是气流脉动压力,I表示单位矩阵<maths num="0003"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA00008956933200000112.GIF" wi="254" he="235" /></maths><img file="FDA00008956933200000113.GIF" wi="49" he="59" />表示由气流粘性引起的脉动流雷诺应力项,<img file="FDA00008956933200000114.GIF" wi="98" he="67" />表示由计算得到的涡粘性引起的脉动流雷诺应力项,<img file="FDA00008956933200000115.GIF" wi="54" he="56" />表示由气流粘性引起的平均流雷诺应力项,<img file="FDA00008956933200000116.GIF" wi="101" he="66" />表示由计算得到的涡粘性引起的平均流雷诺应力项;所述求解过程包含以下步骤:S1,对计算域进行网格划分;S2,对边界条件进行设置;S3,利用GAO‑YONG理性湍流模型的机械能方程对划分好网格、设置完边界条件的计算域进行CFD模拟。 |