CN114943167B - 一种结构网格壁面距离的计算方法、系统、介质和设备 - Google Patents
一种结构网格壁面距离的计算方法、系统、介质和设备 Download PDFInfo
- Publication number
- CN114943167B CN114943167B CN202210606484.9A CN202210606484A CN114943167B CN 114943167 B CN114943167 B CN 114943167B CN 202210606484 A CN202210606484 A CN 202210606484A CN 114943167 B CN114943167 B CN 114943167B
- Authority
- CN
- China
- Prior art keywords
- grid
- wall surface
- structural
- unstructured
- temporary
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Image Generation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明适用于空气动力学领域,提供了一种结构网格壁面距离的计算方法、系统、介质和设备,结构网格壁面距离的计算方法包括如下步骤:通过将结构网格的拓扑关系转换为非结构网格的拓扑关系,获取临时非结构网格;采用非结构壁面距离的计算方法对所述临时非结构网格的壁面距离进行计算,得到临时非结构网格的壁面距离信息;将临时非结构网格的壁面距离信息映射回所述结构网格中,得到所述结构网格的壁面距离信息。本发明提供的一种结构网格壁面距离的计算方法、系统、介质和设备实现了结构、非结构网格壁面距离同时计算的需求,降低了程序的维护成本和使用难度。
Description
技术领域
本发明涉及空气动力学领域,尤其是涉及一种结构网格壁面距离的计算方法、系统、介质和设备。
背景技术
计算流体力学(Computational Fluid Dynamics,简称为CFD),是一门利用计算机和数值算法,对流体力学问题进行数值仿真分析的交叉学科,其研究的目的是为了在以航空航天为代表的众多工业领域得到良好的应用,解决航空航天飞行器研制中的众多关键气动问题。
真实流体运动有层流和湍流两种性质决然不同的流动状态,自然界中的流动绝大多数是湍流。利用CFD模拟层流流动已经较为成熟,但是面对实际流动中的湍流问题,由于其物理本身的高度复杂性,只能做简化处理。目前湍流数值模拟主要有三类方法:直接数值模拟方法(DNS)、大涡模拟方法(LES)、雷诺平均N-S方程方法(RANS)。由于DNS及LES方法需要极其庞大的计算资源,RANS方法是当前CFD工程应用中主要采用的方法。
目前的工程应用中,绝大部分RANS数值模拟都是采用SA和SST两类湍流模型,这两类模型方程中均涉及壁面距离,即空间点到几何壁面的最近距离。壁面距离的精确性直接影响到CFD求解器的稳定性、收敛性和计算结果的精度。
壁面距离计算方法主要有两类:一类是通过几何方法直接计算,另一类是通过求解偏微分方程的形式求解,目前实际应用中基本都是用几何方法计算壁面距离。几何方法计算壁面距离时,对空间中的任意一个点,搜索壁面上的壁面单元(如三角形、四边形)离其最近的距离。根据壁面单元搜索方法的不同,几何计算方法主要有:直接搜索法、沿网格线搜索法、循环盒子法和基于交替数字树算法(Alternating Digital Tree,简称ADT)数据结构的搜索算法。
目前CFD采用的网格拓扑主要是结构网格(Structured Grid)和非结构网格(Unstructured Grid)。结构网格数据结构简单、存储方便、计算简单快捷、网格质量高、计算结果精度高,缺点是生成复杂外形的结构网格十分困难,需要耗费大量人工。非结构网格能很好地处理复杂几何外形,有很优越的几何灵活性,其随机存储的数据结构有利于进行网格自适应。但是,由于需要记忆单元之间的关联信息,因而非结构网格存储量和计算量大。
在实际的工程应用中,根据计算条件的不同,使用的网格类型也不一样。在不同的应用领域,结构、非结构网格有着各自的优势。例如,在低速、亚跨声速领域,由于外形比较复杂,往往采用非结构网格离散,快速获得流场数据。在高超声速领域,由于热流和摩阻精确预测、激波和转捩等现象精细捕捉的需要,具有较高离散精度的结构网格则应用较多。
由于网格拓扑的不同,单一的壁面距离计算程序无法同时对结构、非结构网格进行处理。现有的处理方法是对结构、非结构网格采用不同的程序进行壁面距离计算。这样的处理方法存在两个问题:一是使用两套不同的壁面距离计算程序处理结构、非结构网格,增加了程序的维护成本和使用难度;二是单一的壁面距离计算程序无法满足结构、非结构网格的复杂应用场景。
发明内容
本发明的目的解决不同网格结构中壁面距离的计算方法无法统一的技术问题,提供了一种结构网格壁面距离的计算方法、系统、介质和设备。
本发明第一方面提供了一种结构网格壁面距离的计算方法,包括如下步骤:
步骤S11:通过将结构网格的拓扑关系转换为非结构网格的拓扑关系,获取临时非结构网格;
步骤S12:采用非结构壁面距离的计算方法对所述临时非结构网格的壁面距离进行计算,得到临时非结构网格的壁面距离信息;
步骤S13:将临时非结构网格的壁面距离信息映射回所述结构网格中,得到所述结构网格的壁面距离信息。
进一步的,步骤S11中,将结构网格的拓扑关系转换为非结构网格的拓扑关系的方法为:对结构网格的结构网格点和结构网格单元进行编号;构建所述结构网格点与临时非结构网格点之间的映射关系,构建所述结构网格单元与临时非结构网格单元之间的映射关系;得到临时非结构网格的拓扑关系。
进一步的,所述临时非结构网格的拓扑关系为临时非结构网格单元与临时非结构网格单元点之间的拓扑连接关系;所述临时非结构网格的拓扑连接关系包括临时非结构网格单元的序号、临时非结构网格点的数量和非结构网格点的序号。
进一步的,所述临时非结构网格的壁面距离信息的计算方法为:遍历临时非结构网格的所有壁面单元,将所述壁面单元的坐标信息存入交替数字树数据结构中;构建空间点的几何坐标盒子,在所述数字树数据结构中查找所述几何坐标盒子中包含的壁面单元,得到盒内壁面单元;计算空间点到所述盒内壁面单元的最近距离,得到临时非结构网格的壁面距离信息。
本发明第二方面提供了一种网格壁面距离的计算系统,包括临时非结构网格获取模块、壁面距离信息获取模块、壁面距离信息转换模块,所述临时非结构网格获取模块用于通过将结构网格映射到非结构网格获取临时非结构网格;所述壁面距离信息获取模块用于采用非结构壁面距离的计算方法对所述临时非结构网格的壁面距离进行计算,得到临时非结构网格的壁面距离信息;所述壁面距离信息转换模块用于通过所述映射关系将临时非结构网格的壁面距离信息映射回所述结构网格中,得到所述结构网格的壁面距离信息。
本发明第三方面提供了一种混合网格壁面距离的计算方法,所述混合网格包括结构网格和非结构网格,包括如下步骤:
步骤S10:采用所述的一种结构网格壁面距离的计算方法,得到结构网格的壁面距离信息;
步骤S20:采用非结构壁面距离计算方法,得到非结构网格的壁面距离信息;
步骤S30:得到混合网格壁面距离信息。
进一步的,非结构壁面距离计算方法包括直接搜索法、沿网格线搜索法和/或循环盒子法和交替数字树算法。
本发明第四方面提供了一种混合网格壁面距离的计算系统,包括结构网格壁面距离信息获取模块、非结构网格壁面距离信息获取模块和混合网格壁面距离信息获取模块,所述结构网格壁面距离信息获取模块用于采用所述的一种结构网格壁面距离的计算方法,得到结构网格的壁面距离信息;所述非结构网格壁面距离信息获取模块用于采用非结构壁面距离计算方法,得到非结构网格的壁面距离信息;所述混合网格壁面距离信息获取模块用于得到混合网格壁面距离信息。
本发明第五方面提供了一种可读存储介质,用于存储程序,所述存储程序被执行时,用于实现所述的网格壁面距离的计算方法,或用于实现所述的混合网格壁面距离的计算方法。
本发明第六方面提供了一种电子设备,包括:一个或多个处理器;存储器,其上存储有一个或多个程序;当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现所述的网格壁面距离的计算方法或实现所述的混合网格壁面距离的计算方法。
综上所述,本发明至少具有如下技术效果:
1. 本发明通过将结构网格的拓扑关系转换为非结构网格的拓扑关系,得到临时非结构网格,从而使得结构网格可临时转换为非结构网格,进而为通过非结构网格壁面距离计算方法的应用提供前提条件;
2. 本发明通过临时非结构网格的壁面距离信息映射回所述结构网格中,得到结构网格的壁面距离信息,进而通过非结构网格壁面距离计算方法得到了结构网格的壁面距离信息,为结构网格壁面距离的计算方法提供了一种新思路;
3. 本发明通过构建可采用非结构壁面距离的计算方法进行壁面距离信息获取的临时非结构网格,避免了现有技术中由于网格拓扑的不同,单一的壁面距离计算程序无法同时对结构、非结构网格进行处理的问题,实现了结构、非结构网格壁面距离同时计算的需求,降低了程序的维护成本和使用难度;
4. 采用本发明中结构网格壁面距离的计算方法,使用一套壁面距离计算程序即可能够满足于混合网格中壁面距离的计算,可同时满足高超声速领域结构网格以及低速、亚跨声速领域非结构网格的壁面距离计算的需求,拓宽了现有技术中壁面距离计算方法的应用场景,大大降低了混合网格中壁面距离运算程序的存储量和运算量,有效提高了运算效率,具有成本低、简单易行的优势。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对本发明实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面所描述的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明中结构网格壁面距离计算方法的示意图;
图2是本发明中二维结构网格的示意图;
图3是本发明中结构网格到临时非结构网格的网格点信息映射关系的示意图;
图4是本发明中结构网格到临时非结构网格的网格单元信息映射关系的示意图;
图5是本发明中临时非结构网格的网格单元信息和网格点信息的示意图;
图6是本发明中临时非结构网格拓扑连接关系的示意图;
图7是本发明中ADT方法查找壁面单元的示意图;
图8是本发明中混合网格壁面距离计算方法的示意图。
具体实施方式
以下的说明提供了许多不同的实施例、或是例子,用来实施本发明的不同特征。以下特定例子所描述的元件和排列方式,仅用来精简的表达本发明,其仅作为例子,而并非用以限制本发明。
下面将结合本发明实施例中附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一:
如图1所示,本发明实施例一提供了一种结构网格壁面距离的计算方法,包括如下步骤:
步骤S11:通过将结构网格的拓扑关系转换为非结构网格的拓扑关系,获取临时非结构网格;
步骤S12:采用非结构壁面距离的计算方法对所述临时非结构网格的壁面距离进行计算,得到临时非结构网格的壁面距离信息;
步骤S13:将临时非结构网格的壁面距离信息映射回所述结构网格中,得到所述结构网格的壁面距离信息。
如图2所示,在处理结构网格时,对结构网格A的网格拓扑进行重构,构建一套与结构网格A对应的非结构网格拓扑关系,得到临时非结构网格B,如图5所示。同时将结构网格A与临时非结构网格B之间几何信息的映射关系保存。然后采用非结构壁面距离计算方法对临时非结构网格B的壁面距离进行计算,得到临时非结构网格B的壁面距离信息。最后利用结构网格A与非结构网格B之间的几何信息的映射关系,将非结构网格B的壁面距离信息映射回原结构网格A中,得到结构网格A的壁面距离信息。
通过将结构网格的拓扑关系转换为非结构网格的拓扑关系,得到临时非结构网格,从而使得结构网格可临时转换为非结构网格,进而为通过非结构网格壁面距离计算方法的应用提供前提条件。
通过临时非结构网格的壁面距离信息映射回所述结构网格中,得到结构网格的壁面距离信息,进而通过非结构网格壁面距离计算方法得到了结构网格的壁面距离信息,为结构网格壁面距离的计算方法提供了一种新思路。
通过构建可采用非结构壁面距离的计算方法进行壁面距离信息获取的临时非结构网格,避免了现有技术中由于网格拓扑的不同,单一的壁面距离计算程序无法同时对结构、非结构网格进行处理的问题,实现了结构、非结构网格壁面距离同时计算的需求,降低了程序的维护成本和使用难度。
进一步的,步骤S11中,将结构网格的拓扑关系转换为非结构网格的拓扑关系的方法为:对结构网格的结构网格点和结构网格单元进行编号;构建所述结构网格点与临时非结构网格点之间的映射关系,构建所述结构网格单元与临时非结构网格单元之间的映射关系;得到临时非结构网格的拓扑关系。
结构网格的网格拓扑按照特定的顺序或规律编号,如在
x/y/z三个方向分别用
i-
j-k对点、面、单元进行编号。非结构网格的网格拓扑则没有特定的结构或规律,常用的非结构网格拓扑连接关系有几类:面—点关系、面—单元关系、单元—点关系。
本实施例中将结构网格转换为非结构网格,就是在不修改网格几何信息的前提下,将结构网格的网格拓扑改为非结构的网格拓扑,具有简单易行的特点。
如图2所示,当结构网格为二维结构网格时(
i=5,
j=5),将结构网格的拓扑关系转换为非结构网格的拓扑关系的具体步骤如下:首先,建立结构网格到临时非结构网格的网格点信息的映射关系,如图3所示;同时,建立结构网格到临时非结构网格的网格单元信息的映射关系,如图4所示;合并临时非结构网格的网格点信息和网格单元信息,如图5所示;最后,构建临时非结构网格的拓扑关系
cell2node(单元—点关系),如表1所示。
表1 二维临时非结构网格的拓扑关系示意表
第1列 | 第2列 | 第3列 | 第4列 |
13:4;16、17、22、21 | 14:4;17、18、23、22 | 15:4;18、19、24、23 | 16:4;19、20、25、24 |
9:4;11、12、17、16 | 10:4;12、13、18、17 | 11:4;13、14、19、18 | 12:4;14、15、20、19 |
5:4;6、7、12、11 | 6:4;7、8、13、12 | 7:4;8、9、14、13 | 8:4;9、10、15、14 |
1:4;1、2、7、6 | 2:4;2、3、8、7 | 3:4;3、4、9、8 | 4:4;4、5、10、9 |
进一步的,所述临时非结构网格的拓扑关系为临时非结构网格单元与临时非结构网格单元点之间的拓扑连接关系;所述临时非结构网格的拓扑连接关系包括临时非结构网格单元的序号、临时非结构网格点的数量和非结构网格点的序号。
当结构网格为二维结构网格时,非结构网格的拓扑关系采用单元—点关系时,网格单元与网格点的拓扑连接关系包含三个要素:每个单元的编号、每个单元包含的点的数量和每个单元包含的点序列。具体的,对于每一个单元的单元—点关系,第1个数字是单元编号,第2个数字是单元的点数(通常为4),后续数字是该单元含的点编号(通常为4组数字)。例如,1号单元的单元—点关系如图6所示。
非结构网格的拓扑关系采用面—点关系时,网格面与网格点的拓扑连接关系包含三个要素:每个面的编号、每个面包含的点的数量和每个面包含的点序列。
当结构网格为三维结构网格时,非结构网格的拓扑关系采用单元—点关系时,网格单元与网格点的拓扑连接关系包含三个要素:每个单元的编号、每个单元包含的点的数量和每个单元包含的点序列。具体的,对于每一个单元的单元—点关系,第1个数字是单元编号,第2个数字是单元的点数(通常为8),后续数字是该单元含的点编号(通常为8组数字),可类比如图6所示得到。
当结构网格为三维结构网格时,非结构网格的拓扑关系也可采用面—点关系或面—单元关系,具体转换方式可参照上述方式进行类比得到。
进一步的,所述临时非结构网格的壁面距离信息的计算方法为:遍历临时非结构网格的所有壁面单元,将所述壁面单元的坐标信息存入交替数字树数据结构中;构建空间点的几何坐标盒子,在所述数字树数据结构中查找所述几何坐标盒子中包含的壁面单元,得到盒内壁面单元;计算空间点到所述盒内壁面单元的最近距离,得到临时非结构网格的壁面距离信息。
本实施例使用的非结构网格壁面距离计算方法为交替数字树数据算法(ADT)方法。ADT方法采用一种高效的ADT数据结构存储壁面单元的数据,其基本搜索原理是:对每一个壁面单元,将其几何坐标的最小盒子,依次按照笛卡尔坐标三个方向插入到ADT数据结构中。搜索时,根据设定的搜索半径,构建出空间点的几何坐标盒子,查找此盒子中包含的壁面单元,最后计算空间点到上述搜索到的壁面单元的最近距离,即为此空间点的壁面距离。
如图7所示,本实施例中临时非结构网格的壁面距离信息的计算方法具体为:首先,遍历计算网格的所有壁面单元,将壁面单元的坐标信息存入ADT数据结构中;其次,在计算空间中任意一点的壁面距离时,根据设定的搜索半径
r,构建空间点的几何坐标盒子,即搜索壁面单元的搜索范围,如图7中的虚线方形;然后,利用构建好的空间点的几何坐标盒子,在ADT数据结构中查找几何坐标盒子中包含的壁面单元,即图7中实线圆形与虚线方形相交的壁面单元;最后,计算空间点到几何坐标盒子中壁面单元的最近距离,以此作为空间点的壁面距离,即得到临时非结构网格的壁面距离信息。
实施例二:
如图1所示,本发明实施例二提供了一种网格壁面距离的计算系统,包括临时非结构网格获取模块、壁面距离信息获取模块、壁面距离信息转换模块,所述临时非结构网格获取模块用于通过将结构网格映射到非结构网格获取临时非结构网格;所述壁面距离信息获取模块用于采用非结构壁面距离的计算方法对所述临时非结构网格的壁面距离进行计算,得到临时非结构网格的壁面距离信息;所述壁面距离信息转换模块用于通过所述映射关系将临时非结构网格的壁面距离信息映射回所述结构网格中,得到所述结构网格的壁面距离信息。
实施例三:
如图8所示,本发明实施例三提供了一种混合网格壁面距离的计算方法,所述混合网格包括结构网格和非结构网格,包括如下步骤:
步骤S10:采用所述的一种结构网格壁面距离的计算方法,得到结构网格的壁面距离信息;
步骤S20:采用非结构壁面距离计算方法,得到非结构网格的壁面距离信息;
步骤S30:得到混合网格壁面距离信息。
进一步的,非结构壁面距离计算方法包括直接搜索法、沿网格线搜索法和/或循环盒子法和交替数字树算法。
采用本实施例中结构网格壁面距离的计算方法,使用一套壁面距离计算程序即可能够满足于混合网格中壁面距离的计算,可同时满足高超声速领域结构网格以及低速、亚跨声速领域非结构网格的壁面距离计算的需求,拓宽了现有技术中壁面距离计算方法的应用场景,大大降低了混合网格中壁面距离运算程序的存储量和运算量,有效提高了运算效率,具有成本低、简单易行的优势。
实施例四:
如图8所示,本发明实施例四提供了一种混合网格壁面距离的计算系统,包括结构网格壁面距离信息获取模块、非结构网格壁面距离信息获取模块和混合网格壁面距离信息获取模块,所述结构网格壁面距离信息获取模块用于采用所述的一种结构网格壁面距离的计算方法,得到结构网格的壁面距离信息;所述非结构网格壁面距离信息获取模块用于采用非结构壁面距离计算方法,得到非结构网格的壁面距离信息;所述混合网格壁面距离信息获取模块用于得到混合网格壁面距离信息。
实施例五:
本发明实施例五提供了一种可读存储介质,用于存储程序,所述存储程序被执行时,用于实现所述的网格壁面距离的计算方法,或用于实现所述的混合网格壁面距离的计算方法。
实施例六:
本发明实施例六提供了一种电子设备,包括:一个或多个处理器;存储器,其上存储有一个或多个程序;当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现所述的网格壁面距离的计算方法或实现所述的混合网格壁面距离的计算方法。
本发明提供的一种结构网格壁面距离的计算方法、系统、介质和设备,降低了程序的维护成本和使用难度,基于结构-非结构网格壁面距离统一计算方法,使用一套壁面距离计算程序即可同时完成结构、非结构网格的壁面距离计算工作。并且,扩宽了应用场景,使用结构-非结构网格壁面距离统一计算方法,可同时满足高超声速领域结构网格以及低速、亚跨声速领域非结构网格的壁面距离计算的需求。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种结构网格壁面距离的计算方法,其特征在于,包括如下步骤:
步骤S11:通过将结构网格的拓扑关系转换为非结构网格的拓扑关系,获取临时非结构网格;
步骤S12:采用非结构壁面距离的计算方法对所述临时非结构网格的壁面距离进行计算,得到临时非结构网格的壁面距离信息;
步骤S13:将临时非结构网格的壁面距离信息映射回所述结构网格中,得到所述结构网格的壁面距离信息;
步骤S11中,将结构网格的拓扑关系转换为非结构网格的拓扑关系的方法为:对结构网格的结构网格点和结构网格单元进行编号;构建所述结构网格点与临时非结构网格点之间的映射关系,构建所述结构网格单元与临时非结构网格单元之间的映射关系;得到临时非结构网格的拓扑关系。
2.如权利要求1所述的一种结构网格壁面距离的计算方法,其特征在于,所述临时非结构网格的拓扑关系为临时非结构网格单元与临时非结构网格单元点之间的拓扑连接关系;所述临时非结构网格的拓扑连接关系包括临时非结构网格单元的序号、临时非结构网格点的数量和非结构网格点的序号。
3.如权利要求1所述的一种结构网格壁面距离的计算方法,其特征在于,所述临时非结构网格的壁面距离信息的计算方法为:遍历临时非结构网格的所有壁面单元,将所述壁面单元的坐标信息存入交替数字树数据结构中;构建空间点的几何坐标盒子,在所述数字树数据结构中查找所述几何坐标盒子中包含的壁面单元,得到盒内壁面单元;计算空间点到所述盒内壁面单元的最近距离,得到临时非结构网格的壁面距离信息。
4.一种采用如权利要求1-3任一所述的一种结构网格壁面距离的计算方法的结构网格壁面距离的计算系统,其特征在于,包括临时非结构网格获取模块、壁面距离信息获取模块、壁面距离信息转换模块,
所述临时非结构网格获取模块,用于通过将结构网格的拓扑关系转换为非结构网格的拓扑关系,获取临时非结构网格;
所述壁面距离信息获取模块,用于采用非结构壁面距离的计算方法对所述临时非结构网格的壁面距离进行计算,得到临时非结构网格的壁面距离信息;
所述壁面距离信息转换模块,将临时非结构网格的壁面距离信息映射回所述结构网格中,得到所述结构网格的壁面距离信息。
5.一种混合网格壁面距离的计算方法,其特征在于,所述混合网格包括结构网格和非结构网格,包括如下步骤:
步骤S10:采用如权利要求1-3之一所述的一种结构网格壁面距离的计算方法,得到结构网格的壁面距离信息;
步骤S20:采用非结构壁面距离计算方法,得到非结构网格的壁面距离信息;
步骤S30:得到混合网格壁面距离信息。
6.如权利要求5所述的一种混合网格壁面距离的计算方法,其特征在于,非结构壁面距离计算方法包括直接搜索法、沿网格线搜索法、循环盒子法和交替数字树算法。
7.一种混合网格壁面距离的计算系统,其特征在于,包括结构网格壁面距离信息获取模块、非结构网格壁面距离信息获取模块和混合网格壁面距离信息获取模块,所述结构网格壁面距离信息获取模块用于采用如权利要求1-3之一所述的一种结构网格壁面距离的计算方法,得到结构网格的壁面距离信息;所述非结构网格壁面距离信息获取模块用于采用非结构壁面距离计算方法,得到非结构网格的壁面距离信息;所述混合网格壁面距离信息获取模块用于得到混合网格壁面距离信息。
8.一种可读存储介质,其特征在于,用于存储程序,所述存储程序被执行时,用于实现如权利要求1-3之一所述的网格壁面距离的计算方法,或用于实现如权利要求5-6之一所述的混合网格壁面距离的计算方法。
9.一种电子设备,其特征在于,包括:一个或多个处理器;存储器,其上存储有一个或多个程序;当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现如权利要求1-3之一所述的网格壁面距离的计算方法或实现如权利要求5-6之一所述的混合网格壁面距离的计算方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210606484.9A CN114943167B (zh) | 2022-05-31 | 2022-05-31 | 一种结构网格壁面距离的计算方法、系统、介质和设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210606484.9A CN114943167B (zh) | 2022-05-31 | 2022-05-31 | 一种结构网格壁面距离的计算方法、系统、介质和设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114943167A CN114943167A (zh) | 2022-08-26 |
CN114943167B true CN114943167B (zh) | 2023-04-11 |
Family
ID=82908464
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210606484.9A Active CN114943167B (zh) | 2022-05-31 | 2022-05-31 | 一种结构网格壁面距离的计算方法、系统、介质和设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114943167B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115357849B (zh) * | 2022-10-24 | 2023-03-28 | 中国空气动力研究与发展中心计算空气动力研究所 | 笛卡尔网格下的壁面距离的计算方法及装置 |
CN117540507B (zh) * | 2024-01-08 | 2024-06-04 | 中国空气动力研究与发展中心计算空气动力研究所 | 全环网格壁面距离的计算方法及装置 |
CN118070621B (zh) * | 2024-04-22 | 2024-06-18 | 中国空气动力研究与发展中心计算空气动力研究所 | 固壁边界的处理方法、装置、终端设备及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103793584A (zh) * | 2014-03-13 | 2014-05-14 | 北京大学工学院南京研究院 | 一种快速查找空间网格到壁面单元最短距离的方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107391950A (zh) * | 2017-08-03 | 2017-11-24 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种湍流壁面距离快速计算方法 |
CN109948002A (zh) * | 2019-03-14 | 2019-06-28 | 空气动力学国家重点实验室 | 基于平衡kd树的非结构网格最近壁面距离求解方法 |
CN112818573B (zh) * | 2021-01-26 | 2022-06-14 | 北京理工大学 | 一种用于非结构网格的获取边界层非当地变量信息的方法 |
CN113723027A (zh) * | 2021-09-06 | 2021-11-30 | 复旦大学 | 一种针对弹性飞机的静气动弹性计算方法 |
-
2022
- 2022-05-31 CN CN202210606484.9A patent/CN114943167B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103793584A (zh) * | 2014-03-13 | 2014-05-14 | 北京大学工学院南京研究院 | 一种快速查找空间网格到壁面单元最短距离的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114943167A (zh) | 2022-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114943167B (zh) | 一种结构网格壁面距离的计算方法、系统、介质和设备 | |
CN111768502A (zh) | 一种基于gpu加速技术的非结构网格二维洪水模拟系统 | |
CN112100937B (zh) | 基于壁面距的高效重叠网格并行装配方法 | |
CN104036095A (zh) | 基于区域分解的耦合高精度复杂外形流场快速算法 | |
Capizzano | Automatic generation of locally refined Cartesian meshes: Data management and algorithms | |
CN103778191A (zh) | 一种顾及空间邻近关系的矢量等高线数据划分方法 | |
CN114676522B (zh) | 融合gan和迁移学习的气动形状优化设计方法及系统及设备 | |
CN114861500A (zh) | 基于三维点云自动生成隧道结构有限元模型的方法及系统 | |
Maack et al. | Parallel Computation of Piecewise Linear Morse-Smale Segmentations | |
Yilmaz et al. | Surface conformed linear mesh and data subdivision technique for large-scale flow simulation and visualization in variable intensity computational environment | |
CN109767492A (zh) | 一种变电站三维模型的间距计算方法 | |
CN110610040B (zh) | 一种多相材料梯度点阵结构设计方法 | |
Athkuri et al. | The mid‐point Green‐Gauss gradient method and its efficient implementation in a 3D unstructured finite volume solver | |
CN107967397B (zh) | 一种基于有限元分析的飞行器结构质心漂移量高精度设计方法 | |
CN107644139B (zh) | 一种从cad模型到cae模型的属性映射方法 | |
Gao et al. | A multi-level parallel tie-dye algorithm for auto-CFD | |
CN115146408A (zh) | 一种基于三叉元结构的机床结构件正向设计方法 | |
CN114119882A (zh) | 飞行器动态流场分析中高效的嵌套网格宿主单元搜索方法 | |
CN109049715B (zh) | 一种基于多相材料拓扑优化的3d打印模型数据处理方法 | |
CN114611423A (zh) | 一种三维多相可压缩流固耦合的快速计算方法 | |
Norgard et al. | Robust computation of Morse–Smale complexes of bilinear functions | |
Jamroz et al. | Asynchronous communication in spectral-element and discontinuous Galerkin methods for atmospheric dynamics–a case study using the High-Order Methods Modeling Environment (HOMME-homme_dg_branch) | |
CN114119929B (zh) | 飞行器流场分析中的物面距阵面推进并行计算方法 | |
Ma et al. | Geometric modeling and mesh generation by radial basis functions and their application to ship flow simulations | |
CN117274534B (zh) | 一种三维结构化网格的表面流线生成方法、设备及介质 |
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 |