【技术实现步骤摘要】
一种应用于飞行器外形的壁面距离计算方法及装置
本专利技术属于计算流体动力学
，更具体地，涉及一种应用于飞行器外形的壁面距离计算方法及装置。
技术介绍
目前，现有应用于飞行器外形的壁面距离计算方法主要分为两类，第一类是基于搜索比较的近似计算方法，第二类是基于求解偏微分方程的计算方法。其中，基于求解偏微分方程的计算方法需要迭代隐式求解Eikonal和Hamilton-Jacobi方程，不仅需要另外的编程来实现，而且如何保证对于各种复杂的工程问题也能达到合理快速的收敛性尚不得知。而对于基于搜索比较的近似计算方法而言，较为合理的还是直接计算方法，即采用遍历法计算每个离散点到所有壁面离散点的最小距离，该方法的优点是方法简单、程序实现容易以及不存在迭代收敛问题等，但当壁面离散点密集时，所需的计算量非常巨大，计算耗时长。
技术实现思路
本专利技术的目的在于解决现有应用于飞行器外形的壁面距离计算方法的通用性差或者效率低的问题。为了实现上述目的，本专利技术提供一种应用于飞行器外形的壁面距离计算方法及装置。
【技术保护点】
1.一种应用于飞行器外形的壁面距离计算方法，其特征在于，包括：/n获取包含所述飞行器外形的所有壁面网格点的外接方盒；/n采用二等分法对所述外接方盒进行循环切割，得到多个小外接方盒；/n检测目标多块结构网格中的目标网格单元的相邻网格单元是否计算过壁面距离；/n响应于所述相邻网格单元均未计算过壁面距离的检测结果，采用循环盒子法计算所述目标网格单元的壁面距离；/n响应于至少一个所述相邻网格单元计算过壁面距离的检测结果，采用阵面推进法计算所述目标网格单元及下一网格层中对应网格单元的壁面距离。/n
【技术特征摘要】
1.一种应用于飞行器外形的壁面距离计算方法，其特征在于，包括：
获取包含所述飞行器外形的所有壁面网格点的外接方盒；
采用二等分法对所述外接方盒进行循环切割，得到多个小外接方盒；
检测目标多块结构网格中的目标网格单元的相邻网格单元是否计算过壁面距离；
响应于所述相邻网格单元均未计算过壁面距离的检测结果，采用循环盒子法计算所述目标网格单元的壁面距离；
响应于至少一个所述相邻网格单元计算过壁面距离的检测结果，采用阵面推进法计算所述目标网格单元及下一网格层中对应网格单元的壁面距离。
2.根据权利要求1所述的壁面距离计算方法，其特征在于，所述获取包含所述飞行器外形的所有壁面网格点的外接方盒的步骤包括：
收集所述所有壁面网格点，并存放于一个数组中；
对所述所有壁面网格点按照x坐标进行排序；
在笛卡儿坐标系中计算包括所述所有壁面网格点的外接长方体，得到所述外接方盒。
3.根据权利要求2所述的壁面距离计算方法，其特征在于，所述对所述所有壁面网格点按照x坐标进行排序的步骤采用堆排序的方式实现。
4.根据权利要求1所述的壁面距离计算方法，其特征在于，所述采用二等分法对所述外接方盒进行循环切割，得到多个小外接方盒的步骤包括：
沿所述外接方盒的最长边对所述外接方盒进行二等分，得到两个中间外接方盒；
按照以上切割方式持续对得到的每个中间外接方盒进行二等分，直至切割层次达到预定的切割层次阈值或者得到的中间外接方盒所包含的壁面网格点的数量达到预定的壁面网格点数量阈值；
如此得到的多个中间外接方盒为所述多个小外接方盒。
5.根据权利要求4所述的壁面距离计算方法，其特征在于，所述切割层次阈值为10层；
和/或，
所述壁面网格点数量阈值为50个。
6.根据权利要求1所述的壁面距离计算方法，其特征在于，所述目标网格单元的相邻网格单元为所述目标网格单元在I、J和K矢量方向上的六个邻居单元。
7.根据权利要求1所述的壁面距离计算方法，其特征在于，所述响应于所述相邻网格单元均未计算过壁面距离的检测结果，采用循环盒子法计算所述目标网格单元的壁面距离的步骤包括：
按照距离所述目标网格单元由近及远的顺序对所述多个...
【专利技术属性】
技术研发人员：刘松，马继魁，陈培，刘耀峰，
申请(专利权)人：中国航天空气动力技术研究院，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。