[发明专利]一种脱体涡模拟模型构造方法有效
申请号: | 202111239727.1 | 申请日: | 2021-10-25 |
公开(公告)号: | CN114139465B | 公开(公告)日: | 2023-04-07 |
发明(设计)人: | 郭启龙;李辰;刘朋欣;孙东;袁先旭;陈坚强;黄文锋 | 申请(专利权)人: | 中国空气动力研究与发展中心计算空气动力研究所 |
主分类号: | G06F30/28 | 分类号: | G06F30/28;G06F111/10;G06F113/08;G06F119/14 |
代理公司: | 成都九鼎天元知识产权代理有限公司 51214 | 代理人: | 孙杰 |
地址: | 621052 四*** | 国省代码: | 四川;51 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 脱体涡 模拟 模型 构造 方法 | ||
1.一种脱体涡模拟模型构造方法,其特征在于,通过改进RANS模型长度实现,包括以下步骤:
步骤1、获取RANS模型控制方程中的湍动能破坏项,确定湍动能破坏项中的长度尺度;
步骤2、对长度尺度进行修正,完成新的替代破坏项构建;
步骤3、将新的替代破坏项替代原控制方程中的湍动能破坏项,并对新的替代破坏项进行网格尺度修正,得到新的改进脱体涡模型的控制方程;
步骤4、基于改进脱体涡模型控制方程构建计算机程序,读取流场数据,利用数值方法离散求解改进脱体涡模型的控制方程,得到关于湍流粘性系数的时间空间演变的数据;
所述步骤2中,修正后的长度尺度为:
其中,κ=0.40~0.41为卡门常数,CL=1.35为经验系数,lm为普朗特混合长度,定义为:
其中,vt为湍流涡粘性系数,为流场的涡量,为由速度场梯度得到的反对称张量,xi,xj分别表示空间中i方向和j方向的坐标,i与j取值为1、2、3,表示空间的三个方向;ui为流体的速度矢量在xi方向的分量,uj为流体的速度矢量在xj方向的分量;
所述步骤2中替代破坏项为:
其中,φ为控制方程的输送变量,fk表示联系φ与湍动能k之间的因子。
2.根据权利要求1所述的脱体涡模拟模型构造方法,其特征在于,所述网格尺度修正的具体方法为:采用替代修正后的长度尺度L*;
其中,Δ=max(Δx,Δy,Δz)为当地最大网格尺度,
3.根据权利要求1所述的脱体涡模拟模型构造方法,其特征在于所述RANS模型为S-A模型、k-ωSST模型、k-ω模型或者k-kL模型。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国空气动力研究与发展中心计算空气动力研究所,未经中国空气动力研究与发展中心计算空气动力研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202111239727.1/1.html,转载请声明来源钻瓜专利网。