CN104462829B - 空间推进求解中的复杂多区网格的处理方法 - Google Patents
空间推进求解中的复杂多区网格的处理方法 Download PDFInfo
- Publication number
- CN104462829B CN104462829B CN201410766514.8A CN201410766514A CN104462829B CN 104462829 B CN104462829 B CN 104462829B CN 201410766514 A CN201410766514 A CN 201410766514A CN 104462829 B CN104462829 B CN 104462829B
- Authority
- CN
- China
- Prior art keywords
- grid
- block
- overall situation
- processing method
- region
- 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
Landscapes
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
- Stored Programmes (AREA)
Abstract
本发明公开了一种空间推进求解中的多区网格自动处理方法,其包括:网格拓扑关系带权有向图化、流向网格编号全局化、最短路径搜索算法和基于I向网格全局编号的剖分方式。本发明针对PNS等空间推进程序复杂多区结构网格计算问题,不需手动修改边界信息,可自动生成空间推进程序必需的流向网格分段信息,保证空间推进程序高效求解复杂外形流动,提高空间推进程序的易用性。
Description
技术领域
本发明涉及PNS等空间推进程序复杂多区网格自动处理。
背景技术
PNS或Euler方程空间推进程序可用于高效快速获得超声速飞行器气动性能,是一种理想的飞行器设计初期快速气动性能评估工具。尽管该类工具计算效率较高,但对复杂多块网格适应性较差,使用过程中往往要求对所生成的网格信息进行二次加工,计算准备工作繁琐,极大地抑制了工程人员的使用热情。
如图1所示流向分块网格,空间推进程序计算时需要分两段进行,首先计算块1流场,计算完毕后将块1终止截面流动参数传给块2初始截面作为初场信息,逐层向后推进计算,获得块2流场。为了完成整个流动计算,必须提供给空间推进程序必需的信息,这些信息包括:计算分段包括的网格块、分段范围、分段初场信息如何获得等。对于简单网格可手工生成这些信息,如果网格十分复杂,则手工生成必要信息会十分繁琐,且容易出错,每次网格数量或拓扑关系改变时都要重新计算相关信息,严重影响了空间推进程序的应用。
发明内容
为了解决上述问题,本发明提供了一种空间推进求解中的多区网格自动处理方法,其包括:将网格拓扑关系转化为带权有向图,所述带权有向图以邻接矩阵表示;根据获得的I向网格总数,对尾部任一网格块进行I向网格全局编号;通过I向对接关系对已完成I向网格全局编号的尾部任一网格块的上游直到头部网格块进行I向网格全局编号;通过最短路径搜索算法搜索未进行I向网格全局编号的网格块到已进行全局编号的网格块的最短路径,沿所述最短路径逐一对所述最短路径上的网格块进行I向网格全局编号;采用排序算法对各网格块的全局编号进行排序,并进行I向网格剖分以生成分段信息。
根据如上所述的多区网格自动处理方法中,优选,在将网格拓扑关系转化为带权有向图之前,所述多区网格自动处理方法还包括:通过读取CGNS或Gridgen边界格式获取网格拓扑关系。
根据如上所述的多区网格自动处理方法中,优选,采用递归方法获得I向网格总数。
本发明实施例提供的技术方案带来的有益效果是:
针对空间推进程序复杂多区网格分段问题,本方法能够从Gridgen原生边界格式或CGNS格式中读取网格拓扑信息,将其转化为带权有向图表示,采用最短路径算法对网格I向编号进行全局统一编号,在此基础上完成网格流向自动剖分,生成流向分段信息供空间推进程序使用,从而不需手工对网格信息进行二次处理,避免了人为出错的可能性,避免了繁琐的手工操作,提高了空间推进程序的易用性;以带权有向图表示网格拓扑关系,便于采用成熟的最短路径搜索算法寻找最短路径;I向对接关系以较小权值表示,获得的最短路径倾向于I向对接,极大简化后续的全局编号操作;适应任意复杂网格。
附图说明
图1为本发明实施例提供的一种流向分段推进求解示意图;
图2为本发明实施例提供的一种多区网格拓扑关系带权有向图示意图;
图3为本发明实施例提供的一种空间推进求解中的多区网格自动处理方法的流程示意图;
图4为本发明实施例提供的一种网格块分段示意图。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步详细说明。
本发明实施例提供了一种空间推进求解中的多区网格自动处理方法,参见图1-4,其包括如下步骤:
步骤101,通过读取CGNS或Gridgen边界格式获取网格拓扑关系;
步骤102,将网格拓扑关系转化为带权有向图,带权有向图以邻接矩阵表示;
具体地,首先将多区对接网格的拓扑关系转化为带权有向图。如图2所示,带权有向图以邻接矩阵表示,该矩阵元素(i,j)表示块i指向块j的连接关系,若两块通过I向连接,则认为连接是单向的,由下游块指向上游块并赋以较小权值,若两块通过J向或K向连接,则认为连接是双向的并赋以较大权值,若两块之间没有对接关系,则以无穷大权值表示。在图2中,较小权值以数字1表示,较大权值以数字100表示,在其他的实施中,只需满足较小权值的数值小于较大权值的数值即可。需要说明的是:I,J,K为结构网格的三个方向,对接关系是CFD计算中多块网格之间的一种边界条件类型。由于本发明实施例主要处理对接关系,因此假定附图中网格块之间只要共边就认为存在对接关系。
步骤103,根据获得的I向网格总数,对尾部任一网格块进行I向网格全局编号;
具体地,对全部网格采用递归算法获得I向网格总数,在图1中,I向网格总数为64。由于每个网格都有自己的索引,比如I向有20个网格点,那么I向索引就从1到20。I向全局编号指的是从最上游开始依次向后递增编号。以块2为例,块1的I向共有20个点,那么块2原来本地编号为1的点就成为全局编号为20的点,块2原来本地I向编号为45的点就成为64。网格块的最小和最大I向全局编号的形式分别为Imin,Imax。如图1的块2,其I向网格全局编号为Imin=20,Imax=64。
步骤104,通过I向对接关系对已完成I向网格全局编号的尾部任一网格块的上游直到头部网格块进行I向网格全局编号;
具体地,在对尾部任一网格块进行I向网格全局编号后,如图2中的网格块6,对其I向网格进行全局编号,且通过I向连接对其上游直到头部的网格块进行全局编号,块6的上游直到头部的块只有块5。由于流体流动是有方向的,从哪儿流过来称为上游,头部块(或头部网格快)即最上游的块,尾部块(或尾部网格块)即最下游的块。
步骤105,通过最短路径搜索算法搜索未进行I向网格全局编号的网格块到已进行全局编号的网格块的最短路径,沿最短路径逐一对最短路径上的网格块进行I向网格全局编号;
具体地,对未进行I向网格统一编号的块,采用最短路径算法,通过带权有向图,寻找该块到已进行全局编号块的最短路径,沿此路径逐一对该路径上的块进行I向网格全局编号。以图2为例,假定块6为找到的尾部块并进行全局编号,由步骤104块5也已完成全局编号,现拟对块1进行全局编号,则采用最短路径搜索算法可找到一条最短路径1->5,由块5和块1间的对接信息可对块1进行全局编号,可见随着已进行全局编号的块越来越多,则最短路径查找会越来越快。
步骤106,采用排序算法对各网格块的全局编号进行排序,并进行I向网格剖分以生成分段信息。
具体地,获得所有网格块全局编号后,采用排序算法将各块Imax全局编号由小到大进行排序,形成排序列表。从该列表中依次取值,记所取值为Isect,则任一块只要满足Imin<Isect≤Imax,则该块划入当前段中。若Isect<Imax,则当前块需进行切割,并被划分到前后两个不同的段中;若Isect=Imax,则当前块不需要切割。最后,将分段信息按所需格式输出到文件中供空间推进程序使用。如图4所示,图中有4个网格块,需要分成三个段,段1是Block 1和Block3连接,段2是Block2和Block 3连接,段3是Block2和Block4连接。对于块3来说,如果不切割将横跨两个段:段1和段2。分段信息根据具体的需要(来自程序的需要)而不同,但至少包括:各段的I向范围,每个段包含的网格块列表,以及分段处的边界类型(是同一网格块内部的切割还是与下游的网格块相连)。每个分段首尾都有一个I向全局编号,分别记为Isect_min和Isect_max,若一个段的Isect_max=Isect,则该段为当前段。
即:本发明实施例首先将网格拓扑关系以带权有向图表示,I向对接关系以较小权值表示,横向(J向、K向)对接关系以较大权值表示;然后获得I向网格总数并找到任一尾部网格块,对其进行I向网格统一编号;采用最短路径搜索算法,对未进行全局编号的块进行I向网格全局编号;基于各块网格I向全局编号进行网格剖分,生成分段信息。
综上所述,本发明实施例的有益效果如下:
针对空间推进程序复杂多区网格分段问题,本方法能够从Gridgen原生边界格式或CGNS格式中读取网格拓扑信息,将其转化为带权有向图表示,采用最短路径算法对网格I向编号进行全局统一编号,在此基础上完成网格流向自动剖分,生成流向分段信息供空间推进程序使用,从而不需手工对网格信息进行二次处理,避免了人为出错的可能性,避免了繁琐的手工操作,提高了空间推进程序的易用性;以带权有向图表示网格拓扑关系,便于采用成熟的最短路径搜索算法寻找最短路径;I向对接关系以较小权值表示,获得的最短路径倾向于I向对接,极大简化后续的全局编号操作;适应任意复杂网格。
本发明未详细说明部分属本领域技术人员公知常识。
Claims (3)
1.一种空间推进求解中的多区网格自动处理方法,其用于PNS或Euler方程空间推进程序中,其特征在于,所述多区网格自动处理方法包括:
将网格拓扑关系转化为带权有向图,所述带权有向图以邻接矩阵表示,该矩阵元素(i,j)表示块i指向块j的连接关系,若两块通过I向连接,则认为连接是单向的,由下游块指向上游块并赋以较小权值,若两块通过J向或K向连接,则认为连接是双向的并赋以较大权值,若两块之间没有对接关系,则以无穷大权值表示;
根据获得的I向网格总数,对尾部任一网格块进行I向网格全局编号;
通过I向对接关系对已完成I向网格全局编号的尾部任一网格块的上游直到头部网格块进行I向网格全局编号;
通过最短路径搜索算法搜索未进行I向网格全局编号的网格块到已进行全局编号的网格块的最短路径,沿所述最短路径逐一对所述最短路径上的网格块进行I向网格全局编号;
采用排序算法对各网格块的全局编号进行排序,并进行I向网格剖分以生成分段信息。
2.根据权利要求1所述的多区网格自动处理方法,其特征在于,在将网格拓扑关系转化为带权有向图之前,所述多区网格自动处理方法还包括:
通过读取CGNS或Gridgen边界格式获取网格拓扑关系。
3.根据权利要求1所述的多区网格自动处理方法,其特征在于,采用递归方法获得I向网格总数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410766514.8A CN104462829B (zh) | 2014-12-11 | 2014-12-11 | 空间推进求解中的复杂多区网格的处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410766514.8A CN104462829B (zh) | 2014-12-11 | 2014-12-11 | 空间推进求解中的复杂多区网格的处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104462829A CN104462829A (zh) | 2015-03-25 |
CN104462829B true CN104462829B (zh) | 2019-01-18 |
Family
ID=52908858
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410766514.8A Active CN104462829B (zh) | 2014-12-11 | 2014-12-11 | 空间推进求解中的复杂多区网格的处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104462829B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106021665B (zh) * | 2016-05-10 | 2019-05-24 | 中国航天空气动力技术研究院 | 高超声速飞行器多批次气动数据快速预测方法 |
CN115374542B (zh) * | 2022-10-24 | 2023-03-24 | 中汽研(天津)汽车工程研究院有限公司 | 基于三角形网格位置的网格调整方法、设备和存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009035737A2 (en) * | 2007-06-13 | 2009-03-19 | United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Adaptive refinement tools for tetrahedral unstructured grids |
CN103226540A (zh) * | 2013-05-21 | 2013-07-31 | 中国人民解放军国防科学技术大学 | 基于分组多流的gpu上多区结构网格cfd加速方法 |
CN103345580A (zh) * | 2013-07-02 | 2013-10-09 | 上海大学 | 基于格子Boltzmann方法的并行CFD方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8259104B2 (en) * | 2008-08-28 | 2012-09-04 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Domain decomposition by the advancing-partition method for parallel unstructured grid generation |
-
2014
- 2014-12-11 CN CN201410766514.8A patent/CN104462829B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009035737A2 (en) * | 2007-06-13 | 2009-03-19 | United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Adaptive refinement tools for tetrahedral unstructured grids |
CN103226540A (zh) * | 2013-05-21 | 2013-07-31 | 中国人民解放军国防科学技术大学 | 基于分组多流的gpu上多区结构网格cfd加速方法 |
CN103345580A (zh) * | 2013-07-02 | 2013-10-09 | 上海大学 | 基于格子Boltzmann方法的并行CFD方法 |
Non-Patent Citations (4)
Title |
---|
CFD并行计算中的多区结构网格二次剖分方法与实现;王勇献等;《计算机研究与发展》;2013083;第50卷(第8期);正文第1节、第2.1.2节和第3节 |
基于多区域多代码问题的自动分块算法;张娟等;《计算机工程》;20100531;第36卷(第9期);第73-76页及第79页 |
大规模CFD多区结构网格任务负载平衡算法;唐波等;《计算机工程与科学》;20140731;第36卷(第7期);第1213-1220页 |
经纬网格寻址中最短路径算法优化;吕燕等;《电脑编程技巧与维护》;20120731;第1节和第4节 |
Also Published As
Publication number | Publication date |
---|---|
CN104462829A (zh) | 2015-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Picheny | A stepwise uncertainty reduction approach to constrained global optimization | |
CN106776963B (zh) | 轻量化的bim大数据在线可视化方法和系统 | |
CN102810113B (zh) | 一种针对复杂网络的混合型聚类方法 | |
CN106202335B (zh) | 一种基于云计算框架的交通大数据清洗方法 | |
CN103324796B (zh) | 一种大规模集成电路设计中的绕障直角斯坦纳树构造方法 | |
CN105976048A (zh) | 一种基于改进人工蜂群算法的输电网扩展规划方法 | |
CN112699623A (zh) | 基于非结构网格规则化重构技术的高精度热流计算方法 | |
CN104462829B (zh) | 空间推进求解中的复杂多区网格的处理方法 | |
CN102880760B (zh) | 基于多目标粒子群及拓扑扩张的微波电路设计方法 | |
CN104361185A (zh) | 电缆虚拟设计用布线空间自动生成方法 | |
WO2023241423A1 (zh) | 一种基于流量预测的路由决策方法及系统 | |
Adacher et al. | Decentralized spatial decomposition for traffic signal synchronization | |
CN104935476A (zh) | 一种基于sdn的网络流量矩阵测量方法 | |
CN109214596A (zh) | 求取具有方向约束和障碍限制的栅格最短路径afw算法 | |
CN109584371A (zh) | 空间曲线覆盖三角网格曲面的方法 | |
CN103149840B (zh) | 一种基于动态规划的语义服务组合方法 | |
CN106951615B (zh) | 基于密母算法的电力运输网络拓扑结构设计方法 | |
CN103595652B (zh) | 一种电力通信网络中QoS能效的分级方法 | |
CN104407567B (zh) | 基于特征的蒙皮碎面重构方法 | |
CN109766578A (zh) | 一种三维机翼结冰后表面重构方法 | |
CN107276093B (zh) | 基于场景削减的电力系统概率潮流计算方法 | |
CN104615679A (zh) | 一种基于人工免疫网络的多智能体数据挖掘方法 | |
Bekkerman et al. | Heuristic approximation method for a random flow of events by an MC-flow with arbitrary number of states | |
CN104866678A (zh) | Fpga时序约束布局方法 | |
CN113258567B (zh) | 基于混合聚类算法的电网关键输电断面确定方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |