发明名称 一种GAO-YONG理性湍流模型的处理方法
摘要 本发明公开了一种GAO-YONG理性湍流模型的处理方法,涉及CFD计算技术领域。所述GAO-YONG理性湍流模型的处理方法的求解过程包含以下步骤:S1,对计算域进行网格划分;S2,对边界条件进行设置;S3,利用GAO-YONG理性湍流模型机械能方程计算分析;该模型采用了新的计算方法,对划分好网格、设置完边界条件的计算域进行CFD模拟。本发明的有益效果在于:通过采用新提出的对机械能方程的处理方法,可以降低方程求解过程中的耗散,进而得到更准确的湍流长度的分布情况,更加准确地进行数值模拟。
申请公布号 CN105701272A 申请公布日期 2016.06.22
申请号 CN201511020017.4 申请日期 2015.12.30
申请人 中国航空工业集团公司沈阳发动机设计研究所 发明人 李孝堂;郑文涛
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京航信高科知识产权代理事务所(普通合伙) 11526 代理人 周良玉
主权项 一种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>&rho;</mi><mo>&OverBar;</mo></mover><mover><mi>U</mi><mo>~</mo></mover><mo>&CenterDot;</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>&CenterDot;</mo><mo>&dtri;</mo><mover><mi>U</mi><mo>&OverBar;</mo></mover></mrow><mo>)</mo><mo>)</mo></mrow></mrow><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><mover><mi>&rho;</mi><mo>&OverBar;</mo></mover><mover><mi>U</mi><mo>~</mo></mover><mo>&CenterDot;</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>&CenterDot;</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>&CenterDot;</mo><mo>&dtri;</mo><mrow><mo>(</mo><mrow><mo>&dtri;</mo><mover><mi>U</mi><mo>&OverBar;</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>&CenterDot;</mo><mo>(</mo><mrow><mo>-</mo><mo>&dtri;</mo><mrow><mo>(</mo><mrow><mover><mi>U</mi><mo>&OverBar;</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>&dtri;</mo><mo>&CenterDot;</mo><mrow><mo>(</mo><mrow><mover><mi>&sigma;</mi><mo>~</mo></mover><mo>+</mo><msup><mover><mi>&sigma;</mi><mo>~</mo></mover><mrow><mi>t</mi><mi>u</mi><mi>r</mi><mi>b</mi></mrow></msup><mo>-</mo><mi>&sigma;</mi><mo>-</mo><msup><mi>&sigma;</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模拟。
地址 110015 辽宁省沈阳市沈河区万莲路1号