CN115357849B - 笛卡尔网格下的壁面距离的计算方法及装置 - Google Patents
笛卡尔网格下的壁面距离的计算方法及装置 Download PDFInfo
- Publication number
- CN115357849B CN115357849B CN202211300472.XA CN202211300472A CN115357849B CN 115357849 B CN115357849 B CN 115357849B CN 202211300472 A CN202211300472 A CN 202211300472A CN 115357849 B CN115357849 B CN 115357849B
- Authority
- CN
- China
- Prior art keywords
- target
- cartesian grid
- determining
- object plane
- distance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/901—Indexing; Data structures therefor; Storage structures
- G06F16/9027—Trees
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/903—Querying
- G06F16/90335—Query processing
- G06F16/90348—Query processing by searching ordered data, e.g. alpha-numerically ordered data
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- Image Analysis (AREA)
Abstract
本申请公开了笛卡尔网格下的壁面距离的计算方法及装置,通过获取物面几何文件;根据物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构;确定目标笛卡尔网格点;根据基于嵌套包围盒概念的数据索引结构和目标笛卡尔网格点。通过优化基于嵌套包围盒概念的数据索引结构的最近邻搜索算法中距物面较远数据点回溯过程,实现最小距离对应的三角形的快速定位,确定目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形,精确高效地计算目标笛卡尔网格中心到定位的三角形之间的最小距离;将最小距离确定为壁面距离。
Description
技术领域
本申请属于流体力学技术领域,尤其涉及笛卡尔网格下的壁面距离的计算方法及装置。
背景技术
笛卡尔网格下实现壁面距离的精确高效计算具有重要意义:一方面,由于笛卡尔网格的非贴体特性,通常采用虚拟单元法处理物面边界,需要根据物面确定虚拟单元的壁面距离,寻找参考点并插值获得其物理量,壁面距离精确计算对虚拟单元物面处理精准度有较大影响;另一方面,对于非定常流动或运动物体,笛卡尔网格会根据流场的变化或物体的运动进行自适应,自适应后的网格壁面距离需重新计算,壁面距离计算效率成为制约流动计算效率的关键因素之一。
目前,壁面距离的计算方式主要分为求解偏微分方程和采用几何方法直接计算两类:第一类方法,将最小距离转换为关于波传播问题的偏微分方程数值求解,需要额外的计算花费,壁面距离计算精度受数值离散精度的影响,且对于复杂外形适应性较差。第二类方法是根据空间点与物面离散网格的几何关系计算壁面距离。通常,为简单起见,在物面网格足够密的情况下,求空间网格点的壁面距离时一般采用空间点到物面网格中心的距离代替最小距离,但是这种近似距离往往引起计算误差,不仅影响计算结果的精度,而且对计算稳定性造成影响。为提高计算精度,可采用根据空间点关于物面网格所在平面投影点与物面网格的几何关系计算壁面距离的3D投影法,以及通过坐标变换转化为二维问题的2D法,但这些方法涉及大量向量计算和关系判断,计算量较大。
在几何方法求解壁面距离时,由于物面网格数目较大,计算效率的核心问题是如何快速定位到最短距离相关的物面网格,物面网格数据结构是关键。
最常见的方法是将物面网格顺序存入数组,采用遍历法搜索物面点得到壁面距离。网格规模较小时,计算量尚在承受范围之内。但是对于三维复杂几何外形的流场计算,物面网格数量可达到O(105)~O(106)量级,空间网格点可达到O(108)~O(109)量级,采用遍历法计算量达到O(1013) ~O(1015)量级。此种规模的计算量对整个计算周期的影响是巨大的。采用ADT(Alternating Digital Tree,交替数字二叉树)叉树数据结构存储物面网格,计算效率可大幅提高,成为目前最常用的物面网格数据存储方式,但对复杂外形,ADT平衡性下降,计算效率仍然差强人意。采用传统的KDT存储物面网格,虽然解决了ADT平衡性问题,计算效率进一步提升,但对于离物面较远的空间点,在运用KDT(KDT树,一种多维空间数据索引结构)最近邻搜索算法时,超球面与分割面位置关系判断失准,需要反复回溯,计算量大。
发明内容
本发明意在提供笛卡尔网格下的壁面距离的计算方法及装置,以解决现有技术中存在的不足,本发明要解决的技术问题通过以下技术方案来实现。
第一个方面,本发明实施例提供笛卡尔网格下的壁面距离的计算方法,所述方法包括:
获取物面几何文件;
根据所述物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构;
确定目标笛卡尔网格点;
根据所述基于嵌套包围盒概念的数据索引结构和所述目标笛卡尔网格点,确定所述目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形;
计算所述目标笛卡尔网格点到所述空间三角形的最小距离;
将所述最小距离确定为壁面距离。
可选地,所述根据所述物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构,包括:
确定物面几何文件的目标数据集G中三角形元素的数量N;
根据N个三角形元素,确定所述目标数据集G的最小包围盒的两个顶点;
计算N个三角形元素的中心点在各个维度的方差;
将最大的方差对应的维度确定为划分维;
将所述N个三角形元素的中心点在方差最大的维度上按从小到大的顺序进行排序;
将所述N个三角形元素的中心点加入到所述基于嵌套包围盒概念的数据索引结构中。
可选地,所述根据所述基于嵌套包围盒概念的数据索引结构和所述目标笛卡尔网格点,确定所述目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形,包括:
从所述基于嵌套包围盒概念的数据索引结构的根节点开始查询,根据所述目标笛卡尔网格点与根节点剖分面的位置关系,确定与所述目标笛卡尔网格点距离最近点的目标子树;
对所述目标子树进行递归处理;
当递归访问至叶子节点时,计算所述目标笛卡尔网格点与当前节点之间的距离,并将所述距离作为当前最近距离,将所述当前节点作为与目标笛卡尔网格点最近的节点确定为目标节点;
对所述目标节点进行回溯处理,计算目标笛卡尔网格点与所述目标节点的父节点剖分轴之间的关系;
若所述目标笛卡尔网格点与所述目标节点的剖分轴的距离小于当前最小距离,则对父节点的另一子节点进行查找潜在最近节点,直至根节点;
根据所述目标节点,确定所述目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形。
可选地,所述对所述目标节点进行回溯处理,计算目标笛卡尔网格点与所述目标节点的父节点剖分轴之间的关系,包括:
判断以目标笛卡尔网格点为圆心,目标笛卡尔网格点和目标数据集之间的距离为半径的圆是否与预设线段相交,其中,所述预设线段至少包括剖分轴和当前包围盒;
若相交,则进行正常的回溯操作;
若不相交,则直接排除对应节点另外子树存在最小距离的可能性。可选地,
所述计算所述目标笛卡尔网格点到所述空间三角形的最小距离,包括:
可选地,所述方法还包括:
在存储所述物面几何文件的物面网格时,按逆时针顺序存储物面网格顶点信息。
第二个方面,本发明实施例提供笛卡尔网格下的壁面距离的计算装置,所述装置包括:
获取模块,用于获取物面几何文件;
第一确定模块,用于根据所述物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构;
第二确定模块,用于确定目标笛卡尔网格点;
第一计算模块,用于根据所述基于嵌套包围盒概念的数据索引结构和所述目标笛卡尔网格点,确定所述目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形;
第二计算模块,用于计算所述目标笛卡尔网格点到所述空间三角形的最小距离;
第三确定模块,用于将所述最小距离确定为壁面距离。
可选地,所述第一确定模块,用于:
确定物面几何文件的目标数据集G中三角形元素的数量N;
根据N个三角形元素,确定所述目标数据集G的最小包围盒的两个顶点;
计算N个三角形元素的中心点在各个维度的方差;
将最大的方差对应的维度确定为划分维;
将所述N个三角形元素的中心点在方差最大的维度上按从小到大的顺序进行排序;
将所述N个三角形元素的中心点加入到所述基于嵌套包围盒概念的数据索引结构中。
可选地,所述第一计算模块用于:
从所述基于嵌套包围盒概念的数据索引结构的根节点开始查询,根据所述目标笛卡尔网格点与根节点剖分面的位置关系,确定与所述目标笛卡尔网格点距离最近点的目标子树;
对所述目标子树进行递归处理;
当递归访问至叶子节点时,计算所述目标笛卡尔网格点与当前节点之间的距离,并将所述距离作为当前最近距离,将所述当前节点作为与目标笛卡尔网格点最近的节点确定为目标节点;
对所述目标节点进行回溯处理,计算目标笛卡尔网格点与所述目标节点的父节点剖分轴之间的关系;
若所述目标笛卡尔网格点与所述目标节点的父节点剖分轴的距离小于当前最小距离,则对父节点的另一子节点进行查找潜在最近节点,直至根节点;
根据所述目标节点,确定所述目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形。
可选地,所述第一计算模块用于:
判断以目标笛卡尔网格点为圆心,目标笛卡尔网格点和目标数据集之间的距离为半径的圆是否与预设线段相交,其中,所述预设线段至少包括剖分轴和当前包围盒;
若相交,则进行正常的回溯操作;
若不相交,则直接排除对应节点另外子树存在最小距离的可能性。
可选地,所述第二计算模块用于:
可选地,所述获取模块用于:
在存储所述物面几何文件的物面网格时,按逆时针顺序存储物面网格顶点信息。
第三个方面,本发明实施例提供一种终端设备,包括:至少一个处理器和存储器;
所述存储器存储计算机程序;所述至少一个处理器执行所述存储器存储的计算机程序,以实现第一个方面提供的笛卡尔网格下的壁面距离的计算方法。
第四个方面,本发明实施例提供一种计算机可读存储介质,该计算机可读存储介质中存储有计算机程序,所述计算机程序被执行时实现第一个方面提供的笛卡尔网格下的壁面距离的计算方法。
本发明实施例包括以下优点:
本发明实施例提供的笛卡尔网格下的壁面距离的计算方法及装置,通过获取物面几何文件;根据物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构;确定目标笛卡尔网格点;根据基于嵌套包围盒概念的数据索引结构和目标笛卡尔网格点,确定目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形;计算目标笛卡尔网格点到空间三角形的最小距离;将最小距离确定为壁面距离。通过优化基于嵌套包围盒概念的数据索引结构的最近邻搜索算法中距物面较远数据点回溯过程,实现最小距离对应的三角形的快速定位,进一步,基于三角形参数化方法,通过符号判断和加减乘运算,精确高效地计算目标笛卡尔网格中心到定位的三角形之间的最小距离。
附图说明
为了更清楚地说明本申请实施例或现有的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本申请一实施例中一种笛卡尔网格下的壁面距离的计算方法的流程图;
图2为本申请一实施例中又一种笛卡尔网格下的壁面距离的计算方法的流程图;
图3为本申请一实施例中剖分面相交的情况;
图4为本申请一实施例中目标点与KDT节点之间的位置关系;
图5为本申请一实施例中回溯方法改进前后的回溯过程;
图6是本发明的一种笛卡尔网格下的壁面距离的计算装置实施例的结构框图;
图7是本发明的一种终端设备的结构示意图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,下面将结合具体实施例及相应的附图对本申请的技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
本发明一实施例提供一种笛卡尔网格下的壁面距离的计算方法,用于计算壁面距离。本实施例的执行主体为笛卡尔网格下的壁面距离的计算装置,设置在终端设备上,例如,终端设备至少包括计算机终端等。
参照图1,示出了本发明的一种笛卡尔网格下的壁面距离的计算方法实施例的步骤流程图,该方法具体可以包括如下步骤:
S101、获取物面几何文件;
具体地,终端设备获取物面几何文件。
S102、根据物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构;
具体地,终端设备根据物面三角形顶点信息建立KDT。
S103、确定目标笛卡尔网格点;
具体地,输入目标笛卡尔网格点。
S104、根据基于嵌套包围盒概念的数据索引结构和目标笛卡尔网格点,确定目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形;
终端设备确定目标笛卡尔网格点到物面最小距离的点对应的三角形。
S105、计算目标笛卡尔网格点到空间三角形的最小距离;
S106、将最小距离确定为壁面距离。
本发明实施例增加壁面距离查询模块,壁面距离查询模块是一个单独的模块,在计算壁面距离时,仅需知道空间网格点的位置信息以及物面三角形信息。因此,本文发展的壁面距离计算方法不仅对笛卡尔网格,对结构网格、非结构网格以及重叠网格等也可较好兼容,可方便程序的移植。
本发明实施例提供的笛卡尔网格下的壁面距离的计算方法,通过获取物面几何文件;根据物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构;确定目标笛卡尔网格点;根据基于嵌套包围盒概念的数据索引结构和目标笛卡尔网格点,确定目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形;计算目标笛卡尔网格点到空间三角形的最小距离;将最小距离确定为壁面距离。通过优化基于嵌套包围盒概念的数据索引结构的最近邻搜索算法中距物面较远数据点回溯过程,实现最小距离对应的三角形的快速定位,进一步,基于三角形参数化方法,通过符号判断和加减乘运算,精确高效地计算目标笛卡尔网格中心到定位的三角形之间的最小距离。
本发明又一实施例对上述实施例提供的笛卡尔网格下的壁面距离的计算方法做进一步补充说明。
图2为本申请一实施例中又一种笛卡尔网格下的壁面距离的计算方法的流程图,可选地,根据物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构,包括:
确定物面几何文件的目标数据集G中三角形元素的数量N;
根据N个三角形元素,确定目标数据集G的最小包围盒的两个顶点;
计算N个三角形元素的中心点在各个维度的方差;
将最大的方差对应的维度确定为划分维;
将N个三角形元素的中心点在方差最大的维度上按从小到大的顺序进行排序;
将N个三角形元素的中心点加入到基于嵌套包围盒概念的数据索引结构中。
具体地,根据物面三角形顶点信息建立基于嵌套包围盒概念的KDT。具体分为两个过程:①首先需要确定目标数据集G中三角形元素的数量N。然后,需要根据这N个元素确定G的最小包围盒的两个顶点。②确定划分维。为了保证KDT的平衡,划分维的确定是计算N个三角形的中心点Tcenter在各个维度的方差,然后找出方差最大的维度即为划分维。然后将此N个元素中心点在方差最大的维度上按从小到大的顺序进行排序。将位于中间的元素插入到当前KDT的节点中去。小于此元素的将全部位于当前节点的左子树Left,大于此元素的将全部位于当前节点的右子树Right。对当前节点左子树Left和右子树Right的元素递归执行②,直至所有元素全部插入到KDT中。值得注意的是,在向KDT中插入节点的同时,节点对应的包围盒也在进行同步划分。由于节点包围盒表示此空间区域内所有三角形元素集合的最小包围盒,而在进行空间区域剖分时的依据是三角形中心点确定的平面,所以节点包围盒与包围盒之间,剖分面与剖分面之间会存在相交的部分,如图3所示。
可选地,根据基于嵌套包围盒概念的数据索引结构和目标笛卡尔网格点,确定目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形,包括:
从基于嵌套包围盒概念的数据索引结构的根节点开始查询,根据目标笛卡尔网格点与根节点剖分面的位置关系,确定与目标笛卡尔网格点距离最近点的目标子树;
对目标子树进行递归处理;
当递归访问至叶子节点时,计算目标笛卡尔网格点与当前节点之间的距离,并将距离作为当前最近距离,将当前节点作为与目标笛卡尔网格点最近的节点确定为目标节点;
对目标节点进行回溯处理,计算目标笛卡尔网格点与目标节点的父节点剖分轴之间的关系;
若目标笛卡尔网格点与目标节点的父节点剖分轴的距离小于当前最小距离,则对父节点的另一子节点进行查找潜在最近节点,直至根节点;
根据目标节点,确定目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形。
具体地,确定目标笛卡尔网格点到物面最小距离的点对应的三角形。具体分为三个过程:
①首先从KDT的根节点开始查询,根据目标点与根节点剖分面的位置关系,确定与目标点T潜在的最近点是在根节点的左子树还是右子树。然后对确定的子树递归过程①的操作。
②当递归访问至叶子节点时,记录T与当前节点之间的距离作为当前最近距离dismin,当前节点作为与目标点T最近的节点记为Knearest。
③进行回溯操作。回溯操作从当前最近的点Knearest开始向上进一步查找。构造以T为球心,dismin为半径的球。若球面与Knearest父节点的剖分面相交,则查找Knearest的兄弟节点即Knearest父节点的另外子树的节点,同时计算T与Knearest父节点之间的距离,若小于当前最小距离,则Knearest和dismin均进行更新;若球面与剖分面不相交,则说明Knearest的兄弟节点中不存在比Knearest与T更近的节点。递归此过程直至根节点,则可获得最小距离的点所在的物面三角形。
值得注意的是,当空间网格点与物面距离较远的情况时,上述方法会面临回溯操作效率低下的问题,本发明对回溯过程进行优化。为了直观说明问题,以二维数据点为例进行说明,如图4所示。
首先判断以目标点T为圆心,T和G之间的距离DisTG为半径的圆是否与线段ab相交(ab由剖分轴和当前包围盒确定)。若相交,则进行正常的回溯操作;若不相交,则直接排除对应节点另外子树存在最小距离的可能性。
改进后的回溯方法可极大减少不必要的回溯过程,进而提高最小距离查询效率。改进前后的回溯过程如图5中实线箭头所示。
可选地,对目标节点进行回溯处理,计算目标笛卡尔网格点与目标节点的父节点剖分轴之间的关系,包括:
判断以目标笛卡尔网格点为圆心,目标笛卡尔网格点和目标数据集之间的距离为半径的圆是否与预设线段相交,其中,预设线段至少包括剖分轴和当前包围盒;
若相交,则进行正常的回溯操作;
若不相交,则直接排除对应节点另外子树存在最小距离的可能性。
可选地,计算目标笛卡尔网格点到空间三角形的最小距离,包括:
具体地 ,计算目标点到空间三角形最小距离。此问题可描述为点P和三角形之间的最小距离,其中/>,B为三角形的一个顶点,E0和E1分别为此顶点对应的三角形两条边。点P和三角形内任一点/>距离的平方为椭圆函数:
最小距离问题转换为在D内求连续可微函数的极值问题。令:,当/>且/>时,有。如果/>,最小距离为点P到/>的距离;否则,最小距离出现在D边界上。其中,D是一块二维空间的区域,(s,t)是坐标位于这块区域的点;
可选地,该方法还包括:
在存储物面几何文件的物面网格时,按逆时针顺序存储物面网格顶点信息。
具体地,在提高壁面距离计算的准确性方面,本发明将物面三角形进行参数化(s和t)表示,根据s和t的取值范围将st平面分为7个区域,对不同区域进行分情况讨论,从而将求最小距离问题转化为求约束条件下一维极值问题,计算精确度提高的同时计算量大幅减少。
在提高计算效率方面,采用平衡的KDT存储物面三角形所有顶点信息并建立相应的嵌套包围盒,使得在最小距离查询时尽可能快速排除在目标范围之外的单元;同时优化了当前超球面与剖分平面的位置判断法则,解决了在距离物面较远处因无效回溯次数过多导致查询效率下降的问题。
采用不同几何外形对本发明进行精确性及效率测试,结果如表1所示。从表1可以看出,本发明采用的基于嵌套包围盒概念的KDT方法效率较遍历法以及ADT方法2~3个量级的提升。
存储物面网格时,按逆时针顺序存储物面网格顶点信息(方便笛卡尔网格物面处理),还同时包含KDT节点对应所有物面网格的最小包围盒,KDT节点从上到下形成嵌套包围盒形式。嵌套包围盒的引入保证KDT每个节点包含当前节点对应三角形的所有信息,且在最小距离查询时,仅需要通过与剖分面(剖分面的大小由包围盒确定)坐标对比,就可以快速排除与目标点距离大于当前最小距离的三角形,效率大幅提升。
表1 不同方法壁面距离查询时间
引入三角形参数化方法,将空间点到三角形物面离散网格的最小距离问题转换为约束条件下一维极值问题,仅需通过符号判断和少量加减乘运算,即可确定最小距离,计算精度和效率大幅提高。
需要说明的是,对于方法实施例,为了简单描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本发明实施例并不受所描述的动作顺序的限制,因为依据本发明实施例,某些步骤可以采用其他顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作并不一定是本发明实施例所必须的。
本发明实施例提供的笛卡尔网格下的壁面距离的计算方法,通过获取物面几何文件;根据物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构;确定目标笛卡尔网格点;根据基于嵌套包围盒概念的数据索引结构和目标笛卡尔网格点,确定目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形;计算目标笛卡尔网格点到空间三角形的最小距离;将最小距离确定为壁面距离。通过优化基于嵌套包围盒概念的数据索引结构的最近邻搜索算法中距物面较远数据点回溯过程,实现最小距离对应的三角形的快速定位,进一步,基于三角形参数化方法,通过符号判断和加减乘运算,精确高效地计算目标笛卡尔网格中心到定位的三角形之间的最小距离。
本发明另一实施例提供一种笛卡尔网格下的壁面距离的计算装置,用于执行上述实施例提供的笛卡尔网格下的壁面距离的计算方法。
参照图6,示出了本发明的一种笛卡尔网格下的壁面距离的计算装置实施例的结构框图,该装置具体可以包括如下模块:获取模块601、第一确定模块602、第二确定模块603、第一计算模块604、第二计算模块605和第三确定模块606,其中:
获取模块601用于获取物面几何文件;
第一确定模块602用于根据物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构;
第二确定模块603用于确定目标笛卡尔网格点;
第一计算模块604用于根据基于嵌套包围盒概念的数据索引结构和目标笛卡尔网格点,确定目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形;
第二计算模块605用于计算目标笛卡尔网格点到空间三角形的最小距离;
第三确定模块606用于将最小距离确定为壁面距离。
本发明实施例提供的笛卡尔网格下的壁面距离的计算装置,通过获取物面几何文件;根据物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构;确定目标笛卡尔网格点;根据基于嵌套包围盒概念的数据索引结构和目标笛卡尔网格点,确定目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形;计算目标笛卡尔网格点到空间三角形的最小距离;将最小距离确定为壁面距离。通过优化基于嵌套包围盒概念的数据索引结构的最近邻搜索算法中距物面较远数据点回溯过程,实现最小距离对应的三角形的快速定位,进一步,基于三角形参数化方法,通过符号判断和加减乘运算,精确高效地计算目标笛卡尔网格中心到定位的三角形之间的最小距离。
本发明又一实施例对上述实施例提供的笛卡尔网格下的壁面距离的计算装置做进一步补充说明。
可选地,第一确定模块,用于:
确定物面几何文件的目标数据集G中三角形元素的数量N;
根据N个三角形元素,确定目标数据集G的最小包围盒的两个顶点;
计算N个三角形元素的中心点在各个维度的方差;
将最大的方差对应的维度确定为划分维;
将N个三角形元素的中心点在方差最大的维度上按从小到大的顺序进行排序;
将N个三角形元素的中心点加入到基于嵌套包围盒概念的数据索引结构中。
可选地,第一计算模块用于:
从基于嵌套包围盒概念的数据索引结构的根节点开始查询,根据目标笛卡尔网格点与根节点剖分面的位置关系,确定与目标笛卡尔网格点距离最近点的目标子树;
对目标子树进行递归处理;
当递归访问至叶子节点时,计算目标笛卡尔网格点与当前节点之间的距离,并将距离作为当前最近距离,将当前节点作为与目标笛卡尔网格点最近的节点确定为目标节点;
对目标节点进行回溯处理,计算目标笛卡尔网格点与目标节点的父节点剖分轴之间的关系;
若目标笛卡尔网格点与目标节点的父节点剖分轴的距离小于当前最小距离,则对父节点的另一子节点进行查找潜在最近节点,直至根节点;
根据目标节点,确定目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形。
可选地,第一计算模块用于:
判断以目标笛卡尔网格点为圆心,目标笛卡尔网格点和目标数据集之间的距离为半径的圆是否与预设线段相交,其中,预设线段至少包括剖分轴和当前包围盒;
若相交,则进行正常的回溯操作;
若不相交,则直接排除对应节点另外子树存在最小距离的可能性。
可选地,第二计算模块用于:
可选地,获取模块用于:
在存储物面几何文件的物面网格时,按逆时针顺序存储物面网格顶点信息。
对于装置实施例而言,由于其与方法实施例基本相似,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本发明实施例提供的笛卡尔网格下的壁面距离的计算装置,通过获取物面几何文件;根据物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构;确定目标笛卡尔网格点;根据基于嵌套包围盒概念的数据索引结构和目标笛卡尔网格点,确定目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形;计算目标笛卡尔网格点到空间三角形的最小距离;将最小距离确定为壁面距离。通过优化基于嵌套包围盒概念的数据索引结构的最近邻搜索算法中距物面较远数据点回溯过程,实现最小距离对应的三角形的快速定位,进一步,基于三角形参数化方法,通过符号判断和加减乘运算,精确高效地计算目标笛卡尔网格中心到定位的三角形之间的最小距离。
本发明再一实施例提供一种终端设备,用于执行上述实施例提供的笛卡尔网格下的壁面距离的计算方法。
图7是本发明的一种终端设备的结构示意图,如图7所示,该终端设备包括:至少一个处理器701和存储器702;
存储器存储计算机程序;至少一个处理器执行存储器存储的计算机程序,以实现上述实施例提供的笛卡尔网格下的壁面距离的计算方法。
本实施例提供的终端设备,通过获取物面几何文件;根据物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构;确定目标笛卡尔网格点;根据基于嵌套包围盒概念的数据索引结构和目标笛卡尔网格点,确定目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形;计算目标笛卡尔网格点到空间三角形的最小距离;将最小距离确定为壁面距离。通过优化基于嵌套包围盒概念的数据索引结构的最近邻搜索算法中距物面较远数据点回溯过程,实现最小距离对应的三角形的快速定位,进一步,基于三角形参数化方法,通过符号判断和加减乘运算,精确高效地计算目标笛卡尔网格中心到定位的三角形之间的最小距离。
本申请又一实施例提供一种计算机可读存储介质,该计算机可读存储介质中存储有计算机程序,计算机程序被执行时实现上述任一实施例提供的笛卡尔网格下的壁面距离的计算方法。
根据本实施例的计算机可读存储介质,通过获取物面几何文件;根据物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构;确定目标笛卡尔网格点;根据基于嵌套包围盒概念的数据索引结构和目标笛卡尔网格点,确定目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形;计算目标笛卡尔网格点到空间三角形的最小距离;将最小距离确定为壁面距离。通过优化基于嵌套包围盒概念的数据索引结构的最近邻搜索算法中距物面较远数据点回溯过程,实现最小距离对应的三角形的快速定位,进一步,基于三角形参数化方法,通过符号判断和加减乘运算,精确高效地计算目标笛卡尔网格中心到定位的三角形之间的最小距离。
应该指出,上述详细说明都是示例性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语均具有与本申请所属技术领域的普通技术人员的通常理解所相同的含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式。此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
需要说明的是,本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的术语在适当情况下可以互换,以便这里描述的本申请的实施方式能够以除了在这里图示或描述的那些以外的顺序实施。
此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含。例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
为了便于描述,在这里可以使用空间相对术语,如“在……之上”、“在……上方”、“在……上表面”、“上面的”等,用来描述如在图中所示的一个器件或特征与其他器件或特征的空间位置关系。应当理解的是,空间相对术语旨在包含除了器件在图中所描述的方位之外的在使用或操作中的不同方位。例如,如果附图中的器件被倒置,则描述为“在其他器件或构造上方”或“在其他器件或构造之上”的器件之后将被定位为“在其他器件或构造下方”或“在其他器件或构造之下”。因而,示例性术语“在……上方”可以包括“在……上方”和“在……下方”两种方位。该器件也可以其他不同方式定位,如旋转90度或处于其他方位,并且对这里所使用的空间相对描述作出相应解释。
在上面详细的说明中,参考了附图,附图形成本文的一部分。在附图中,类似的符号典型地确定类似的部件,除非上下文以其他方式指明。在详细的说明书、附图及权利要求书中所描述的图示说明的实施方案不意味是限制性的。在不脱离本文所呈现的主题的精神或范围下,其他实施方案可以被使用,并且可以作其他改变。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.笛卡尔网格下的壁面距离的计算方法,其特征在于,所述方法包括:
获取物面几何文件;
根据所述物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构;
所述根据所述物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构,包括:
确定物面几何文件的目标数据集G中三角形元素的数量N;
根据N个三角形元素,确定所述目标数据集G的最小包围盒的两个顶点;
计算N个三角形元素的中心点在各个维度的方差;
将最大的方差对应的维度确定为划分维;
将所述N个三角形元素的中心点在方差最大的维度上按从小到大的顺序进行排序;
将所述N个三角形元素的中心点加入到所述基于嵌套包围盒概念的数据索引结构中;
确定目标笛卡尔网格点;
根据所述基于嵌套包围盒概念的数据索引结构和所述目标笛卡尔网格点,确定所述目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形;
计算所述目标笛卡尔网格点到所述空间三角形的最小距离;
将所述最小距离确定为壁面距离。
2.根据权利要求1所述的笛卡尔网格下的壁面距离的计算方法,其特征在于,所述根据所述基于嵌套包围盒概念的数据索引结构和所述目标笛卡尔网格点,确定所述目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形,包括:
从所述基于嵌套包围盒概念的数据索引结构的根节点开始查询,根据所述目标笛卡尔网格点与根节点剖分面的位置关系,确定与所述目标笛卡尔网格点距离最近点的目标子树;
对所述目标子树进行递归处理;
当递归访问至叶子节点时,计算所述目标笛卡尔网格点与当前节点之间的距离,并将所述距离作为当前最近距离,将所述当前节点作为与目标笛卡尔网格点最近的节点确定为目标节点;
对所述目标节点进行回溯处理,计算目标笛卡尔网格点与所述目标节点的父节点剖分轴之间的关系;
若所述目标笛卡尔网格点与所述目标节点的父节点剖分轴的距离小于当前最小距离,则对父节点的另一子节点进行查找潜在最近节点,直至根节点;
根据所述目标节点,确定所述目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形。
3.根据权利要求2所述的笛卡尔网格下的壁面距离的计算方法,其特征在于,所述对所述目标节点进行回溯处理,计算目标笛卡尔网格点与所述目标节点的父节点剖分轴之间的关系,包括:
判断以目标笛卡尔网格点为圆心,目标笛卡尔网格点和目标数据集之间的距离为半径的圆是否与预设线段相交,其中,所述预设线段至少包括剖分轴和当前包围盒;
若相交,则进行正常的回溯操作;
若不相交,则直接排除对应节点另外子树存在最小距离的可能性。
5.根据权利要求1所述的笛卡尔网格下的壁面距离的计算方法,其特征在于,所述方法还包括:
在存储所述物面几何文件的物面网格时,按逆时针顺序存储物面网格顶点信息。
6.笛卡尔网格下的壁面距离的计算装置,其特征在于,所述装置包括:
获取模块,用于获取物面几何文件;
第一确定模块,用于根据所述物面几何文件的物面三角形顶点信息,确定基于嵌套包围盒概念的数据索引结构; 所述第一确定模块,用于:
确定物面几何文件的目标数据集G中三角形元素的数量N;
根据N个三角形元素,确定所述目标数据集G的最小包围盒的两个顶点;
计算N个三角形元素的中心点在各个维度的方差;
将最大的方差对应的维度确定为划分维;
将所述N个三角形元素的中心点在方差最大的维度上按从小到大的顺序进行排序;
将所述N个三角形元素的中心点加入到所述基于嵌套包围盒概念的数据索引结构中;
第二确定模块,用于确定目标笛卡尔网格点;
第一计算模块,用于根据所述基于嵌套包围盒概念的数据索引结构和所述目标笛卡尔网格点,确定所述目标笛卡尔网格点到物面几何文件的最小距离的点对应的空间三角形;
第二计算模块,用于计算所述目标笛卡尔网格点到所述空间三角形的最小距离;
第三确定模块,用于将所述最小距离确定为壁面距离。
7.一种终端设备,其特征在于,包括:至少一个处理器和存储器;
所述存储器存储计算机程序;所述至少一个处理器执行所述存储器存储的计算机程序,以实现权利要求1-5中任一项所述的笛卡尔网格下的壁面距离的计算方法。
8.一种计算机可读存储介质,其特征在于,该计算机可读存储介质中存储有计算机程序,所述计算机程序被执行时实现权利要求1-5中任一项所述的笛卡尔网格下的壁面距离的计算方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211300472.XA CN115357849B (zh) | 2022-10-24 | 2022-10-24 | 笛卡尔网格下的壁面距离的计算方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211300472.XA CN115357849B (zh) | 2022-10-24 | 2022-10-24 | 笛卡尔网格下的壁面距离的计算方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115357849A CN115357849A (zh) | 2022-11-18 |
CN115357849B true CN115357849B (zh) | 2023-03-28 |
Family
ID=84008511
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211300472.XA Active CN115357849B (zh) | 2022-10-24 | 2022-10-24 | 笛卡尔网格下的壁面距离的计算方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115357849B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115564925B (zh) * | 2022-11-21 | 2024-04-16 | 广州中望龙腾软件股份有限公司 | 基于B-rep模型和笛卡尔网格切片的网格生成方法 |
CN116502567B (zh) * | 2023-06-28 | 2023-09-12 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种非结构网格流场的插值求解方法、装置、设备、介质 |
CN116861822B (zh) * | 2023-09-05 | 2024-03-08 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于笛卡尔网格的物面边界的处理方法及装置 |
CN117390775B (zh) * | 2023-12-13 | 2024-02-27 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种粘性笛卡尔网格生成方法、设备、介质及系统 |
CN117540507A (zh) * | 2024-01-08 | 2024-02-09 | 中国空气动力研究与发展中心计算空气动力研究所 | 全环网格壁面距离的计算方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111783343A (zh) * | 2020-07-08 | 2020-10-16 | 中国人民解放军国防科技大学 | 非结构网格方向模板选择方法与流场求解方法 |
CN113673186A (zh) * | 2021-08-26 | 2021-11-19 | 苏州流场信息技术有限公司 | 一种基于stl文件的笛卡尔网格快速生成方法 |
CN113689556A (zh) * | 2021-10-25 | 2021-11-23 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种块自适应型笛卡尔网格快速图映射方法及系统 |
CN114943167A (zh) * | 2022-05-31 | 2022-08-26 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种结构网格壁面距离的计算方法、系统、介质和设备 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10002109B2 (en) * | 2016-06-24 | 2018-06-19 | Freeport-Mcmoran Inc. | Systems and methods of correlating satellite position data with terrestrial features |
CN107391950A (zh) * | 2017-08-03 | 2017-11-24 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种湍流壁面距离快速计算方法 |
CN108764282A (zh) * | 2018-04-19 | 2018-11-06 | 中国科学院计算技术研究所 | 一种类别增量行为识别方法和系统 |
CN109948002A (zh) * | 2019-03-14 | 2019-06-28 | 空气动力学国家重点实验室 | 基于平衡kd树的非结构网格最近壁面距离求解方法 |
CN110110446B (zh) * | 2019-05-09 | 2023-01-13 | 中国航空工业集团公司西安航空计算技术研究所 | 一种可实现重叠区快速优化的并行嵌套网格方法及装置 |
CN113449208B (zh) * | 2020-03-26 | 2022-09-02 | 阿里巴巴集团控股有限公司 | 空间查询方法、设备、系统及存储介质 |
FR3110974A1 (fr) * | 2020-05-28 | 2021-12-03 | Orange | Procédé et dispositif de détermination de positions géographiques d’un traceur de localisation géographique |
US11282291B1 (en) * | 2021-02-09 | 2022-03-22 | URC Ventures, Inc. | Determining object structure using fixed-location cameras with only partial view of object |
CN112927364B (zh) * | 2021-04-08 | 2022-05-17 | 北京航空航天大学 | 一种面向地质结构力学分析的复合场网格构造方法 |
CN114357625A (zh) * | 2022-01-07 | 2022-04-15 | 江苏科技大学 | 一种适用于自由变形参数化的几何约束计算技术 |
CN114722745B (zh) * | 2022-06-10 | 2022-08-26 | 中国空气动力研究与发展中心计算空气动力研究所 | 湍流壁面距离计算方法、装置、计算机设备和存储介质 |
CN114996858B (zh) * | 2022-07-14 | 2022-10-25 | 中国空气动力研究与发展中心计算空气动力研究所 | 飞行器仿真方法、装置、终端设备和存储介质 |
-
2022
- 2022-10-24 CN CN202211300472.XA patent/CN115357849B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111783343A (zh) * | 2020-07-08 | 2020-10-16 | 中国人民解放军国防科技大学 | 非结构网格方向模板选择方法与流场求解方法 |
CN113673186A (zh) * | 2021-08-26 | 2021-11-19 | 苏州流场信息技术有限公司 | 一种基于stl文件的笛卡尔网格快速生成方法 |
CN113689556A (zh) * | 2021-10-25 | 2021-11-23 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种块自适应型笛卡尔网格快速图映射方法及系统 |
CN114943167A (zh) * | 2022-05-31 | 2022-08-26 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种结构网格壁面距离的计算方法、系统、介质和设备 |
Also Published As
Publication number | Publication date |
---|---|
CN115357849A (zh) | 2022-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115357849B (zh) | 笛卡尔网格下的壁面距离的计算方法及装置 | |
Shi et al. | Adaptive simplification of point cloud using k-means clustering | |
JP4325877B2 (ja) | 行列の高速高精度特異値分解方法、プログラムおよび装置 | |
WO2022011851A1 (zh) | 度量空间划分方式评价方法、装置、计算机设备及存储介质 | |
US20040233191A1 (en) | Robust tetrahedralization and triangulation method with applications in VLSI layout design and manufacturability | |
US10580114B2 (en) | Methods and systems for real time 3D-space search and point-cloud registration using a dimension-shuffle transform | |
KR20060047436A (ko) | 2차원 및 3차원 도형의 데이터를 컴퓨터의 메모리에기록하는 데이터 구조, 프로그램 및 기록 매체 | |
Doraiswamy et al. | Output-sensitive construction of Reeb graphs | |
JP2010505170A (ja) | 数値的にシミュレートされた2つのオブジェクト間の衝突を検出する方法および装置 | |
de Magalhães et al. | An efficient and exact parallel algorithm for intersecting large 3-D triangular meshes using arithmetic filters | |
KR20230043958A (ko) | 3 차원 재구성 중 이미지 그루핑 방법 및 장치, 전자 기기 및 컴퓨터 판독 가능한 저장매체 | |
CN114820979B (zh) | 三维网格模型的处理方法、装置和存储介质 | |
US6865479B2 (en) | Programmatically calculating paths from a spatially-enabled database | |
CN111340100B (zh) | Bim模型的相似度计算方法 | |
CN115346005B (zh) | 基于嵌套包围盒概念用于物面网格的数据结构构建方法 | |
Olliff et al. | Efficient searching in meshfree methods | |
CN116484702B (zh) | 一种适用于任意单元类型的非结构嵌套网格计算方法 | |
CN117115393A (zh) | 一种基于gpu的nurbs曲面并行求交方法、设备及存储介质 | |
CN110570511B (zh) | 点云数据的处理方法、装置、系统和存储介质 | |
Khamayseh et al. | Use of the spatial kD-tree in computational physics applications | |
CN114818309A (zh) | 一种基于离散点表示的二维标量场插值方法、装置及计算机设备 | |
Li et al. | Contracting medial surfaces isotropically for fast extraction of centred curve skeletons | |
CN115145930A (zh) | 基于树状层次索引的gis矢量数据层次编码方法及装置 | |
Haunert et al. | Partitioning polygons via graph augmentation | |
JP2007193552A (ja) | 面モデルの作成装置と作成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |