[发明专利]一种飞行器流场特性的重映方法有效
申请号: | 201910792599.X | 申请日: | 2019-08-26 |
公开(公告)号: | CN110705016B | 公开(公告)日: | 2022-10-28 |
发明(设计)人: | 黄飞;靳旭红;陈智;程晓丽;俞继军;沈清 | 申请(专利权)人: | 中国航天空气动力技术研究院 |
主分类号: | G06F30/20 | 分类号: | G06F30/20 |
代理公司: | 中国航天科技专利中心 11009 | 代理人: | 庞静 |
地址: | 100074 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 飞行器 特性 重映 方法 | ||
1.一种飞行器流场特性的重映方法,其特征在于:包括以下步骤:
(1)读入源网格单元的已知流场变量Fi及源网格单元与目标网格单元各自的位置坐标记为原始坐标,规定两个单元网格的顶点的编号顺序标号,所述的源网格单元与目标网格单元均为二维网格;
(2)针对当前目标网格单元,依次对与当前目标网格单元相交的每个源网格单元进行等参变换,在变换后的网格单元内均匀布置n个随机点,将每个随机点ci的位置坐标重新变换至原始坐标(Xi,Yi),i=1~n;
(3)判断并统计n个随机坐标点ci中位于目标网格单元的个数m,并求解每个源网格单元与当前目标网格单元的相交面积Si;
(4)根据获得的相交面积确定当前目标网格单元内的流场变量F;
(5)对每一个目标网格单元重复(2)~(4)步骤,则可得到飞行器整个目标网格单元下的流场变量;
所述的目标网格单元为待求解流场信息量的新网格单元,源网格单元为已知流场信息的网格单元。
2.一种飞行器流场特性的重映方法,其特征在于:包括以下步骤:
(1)读入源网格单元的已知流场变量Fi及源网格单元与目标网格单元各自的位置坐标记为原始坐标,所述的源网格单元与目标网格单元均为三维网格;
(2)针对当前目标网格单元,依次对与当前目标网格单元相交的每个源网格单元进行等参变换,在变换后的网格单元内均匀布置n个随机点,将每个随机点ci的位置坐标重新变换至原始坐标(Xi,Yi,Zi),i=1~n;
(3)判断并统计n个随机坐标点ci中位于目标网格单元的个数m,并求解每个源网格单元与当前目标网格单元的相交体积Vi;
(4)根据获得的相交体积确定当前目标网格单元内的流场变量F;
(5)对每一个目标网格单元重复(2)~(4)步骤,则可得到飞行器整个目标网格单元下的流场变量;
所述的目标网格单元为待求解流场信息量的新网格单元,源网格单元为已知流场信息的网格单元。
3.根据权利要求1或2所述的一种飞行器流场特性的重映方法,其特征在于:所述步骤(1)中流场变量F包括密度ρ、动量密度ρU、能量密度ρE。
4.根据权利要求1所述的一种飞行器流场特性的重映方法,其特征在于:所述步骤(1)中源网格单元与目标网格单元为三角形单元或者四边形单元,源网格单元与目标网格单元的网格形状可以相同,也可以不同。
5.根据权利要求1所述的一种飞行器流场特性的重映方法,其特征在于:所述的等参变换过程中,源网格单元内变换前后的随机点坐标转换式为其中,b为源网格单元的顶点数;
当b=3时,λ1=rf1,λ2=rf2(1-rf1),λ3=1-rf1-rf2,rf1、rf2为均匀随机数;
当b=4时,λ1=(1-rf1)(1-rf2),λ2=rf1(1-rf2),λ3=rf2rf1,λ4=rf2(1-rf1),rf1、rf2为均匀随机数。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国航天空气动力技术研究院,未经中国航天空气动力技术研究院许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910792599.X/1.html,转载请声明来源钻瓜专利网。