[发明专利]多孔介质LBM计算网格生成方法有效
申请号: | 202110416484.8 | 申请日: | 2021-04-19 |
公开(公告)号: | CN112992294B | 公开(公告)日: | 2021-08-10 |
发明(设计)人: | 肖光明;杜雁霞;刘磊;张超;魏东;杨肖峰;向静;邱芷薇 | 申请(专利权)人: | 中国空气动力研究与发展中心计算空气动力研究所 |
主分类号: | G16C60/00 | 分类号: | G16C60/00 |
代理公司: | 成都九鼎天元知识产权代理有限公司 51214 | 代理人: | 贾年龙 |
地址: | 621052 四*** | 国省代码: | 四川;51 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 多孔 介质 lbm 计算 网格 生成 方法 | ||
1.多孔介质LBM计算网格生成方法,其特征在于,包括步骤:
S1,获取多孔介质几何数模文件;
S2,根据多孔介质几何数模文件的所有表面单元各顶点的坐标来判断几何空间域的范围;
S3,根据LBM计算的网格规模和几何空间域的范围,计算LBM网格空间分辨率、物理求解域范围及各方向的网格节点数;
S4,根据多孔介质几何数模文件和LBM网格空间分辨率设置,计算多孔介质几何数模文件的各表面单元在LBM计算空间的区域范围;
S5,根据各表面单元在LBM计算空间j-k平面的自适应区域范围,然后将该区域范围内的方向射线还原为真实几何空间的方向射线,利用射线相交来判断与射线对应表面单元是否存在交点并标记;
S6,根据LBM计算空间方向射线的交点标记结果进行流/固求解域与边界的判断;
S7,根据交点排列的奇偶性对流/固区域进行区分;
S8,根据各表面单元在LBM计算空间k-i平面与i-j平面的自适应区域范围,分别再从j方向和k方向对LBM计算空间内的流/固区域与边界进行标识,进而完成LBM计算网格的生成。
2.根据权利要求1所述的多孔介质LBM计算网格生成方法,其特征在于,在步骤S1中,通过X-CT扫描设备获得多孔介质几何数模文件;所述多孔介质几何数模文件具有所有表面单元各顶点的空间坐标信息;其中,为表面单元编号,为表面单元顶点编号,下标和分别表示第个表面单元及该表面单元对应的第个顶点。
3.根据权利要求1或2任一所述的多孔介质LBM计算网格生成方法,其特征在于,在步骤S2中,通过寻找多孔介质几何数模文件的所有表面单元各顶点的x,y,z坐标的最大、最小值 来判断几何空间域的范围;所述几何空间域的范围指的是x, y, z各方向的最长距离,采用如下公式计算:
。
4.根据权利要求3所述的多孔介质LBM计算网格生成方法,其特征在于,在步骤S3中,根据LBM计算的网格规模以及步骤S2中的几何空间域的范围,计算LBM物理求解域范围及各方向的网格节点数,在LBM网格空间分辨率满足
的基础上,几何空间域的范围包含于LBM物理求解域范围,即
并使LBM计算网格总节点数接近于设置的网格规模,即
此时,对应LBM计算空间的区域范围为、及;其中,、分别为LBM物理求解域x,y,z方向的总长度,、、分别为LBM计算空间i,j,k方向的网格节点数。
5.根据权利要求4所述的多孔介质LBM计算网格生成方法,其特征在于,
在步骤S4中,根据多孔介质几何数模文件及LBM网格空间分辨率设置,依次计算各表面单元在LBM计算空间的区域范围、及;其中,
同时,使得表面单元的各顶点包含于对应的LBM物理求解域,即
从而实现复杂三维几何模型在LBM计算空间的区域自适应;其中,为表面单元编号,为表面单元顶点编号,下标和分别表示第个表面单元及该表面单元对应的第个顶点,、分别为表面单元对应i方向最小及最大网格节点编号,、分别为表面单元对应j方向最小及最大网格节点编号,、分别为表面单元对应k方向最小及最大网格节点编号,、、分别表示表面单元对应顶点的x,y,z坐标。
6.根据权利要求5所述的多孔介质LBM计算网格生成方法,其特征在于,
在步骤S5中,根据各表面单元在LBM计算空间j-k平面的自适应区域范围
将该区域范围内的方向射线还原为真实几何空间的方向射线,对应的y、z坐标分别为
并利用射线相交法判断与对应表面单元是否存在交点;若存在交点,则计算相关交点的x坐标位置,并对LBM计算空间位置的i方向射线进行交点标记,即交点总数;若不存在交点,,交点坐标;按照上述方式依次完成所有表面单元适应区域范围内i方向射线的交点标记;其中,、分别表示在j-k平面选取i方向射线与表面单元进行交点判断的j编号与k编号,、分别表示在j-k平面选取i方向射线j、k编号对应的实际y,z坐标。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国空气动力研究与发展中心计算空气动力研究所,未经中国空气动力研究与发展中心计算空气动力研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202110416484.8/1.html,转载请声明来源钻瓜专利网。