CN107403466A - 基于全局加密的超大规模非结构网格生成方法 - Google Patents

基于全局加密的超大规模非结构网格生成方法 Download PDF

Info

Publication number
CN107403466A
CN107403466A CN201710655494.0A CN201710655494A CN107403466A CN 107403466 A CN107403466 A CN 107403466A CN 201710655494 A CN201710655494 A CN 201710655494A CN 107403466 A CN107403466 A CN 107403466A
Authority
CN
China
Prior art keywords
object plane
point
ultra
overall situation
encryption
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.)
Pending
Application number
CN201710655494.0A
Other languages
English (en)
Inventor
赵钟
何磊
张扬
张来平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN201710655494.0A priority Critical patent/CN107403466A/zh
Publication of CN107403466A publication Critical patent/CN107403466A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明公开了一种基于全局加密的超大规模非结构网格生成方法，包括如下步骤：步骤一、在初始网格单元的基础上，将所有网格单元剖分加密；步骤二、对剖分加密后的网格单元进行物面保形；步骤三、采用并行化的径向基函数方法实现空间点变形。与现有技术相比，本发明的积极效果是：由于采用了基于全局初始网格加密的方式，因此可以并行化生成任意类型的网格单元；只要有足够的计算机资源，就可以不断加密，生成任意超大规模的网格；本发明的方法是在初始网格上不断“全局加密‑物面投影‑空间变形”生成，因此效率极高。

Description

基于全局加密的超大规模非结构网格生成方法
技术领域
本发明涉及一种基于全局加密的超大规模非结构网格生成方法。
背景技术
在航空航天飞行器设计、流体机械、高速列车/汽车、风能和风工程等工业领域，流体力学是进行产品设计的主干学科。计算流体力学(CFD)是流体力学、计算数学、计算机的交叉学科，采用计算机对流体动力学方程进行仿真，从而得到空间流动，为产品设计提供依据。随着计算机的发展，CFD在工业产品的设计过程中正发挥越来越重要的作用。CFD计算的第一步是要进行网格生成，即将现实的三维空间离散成为四面体、六面体等单元。自上世纪50年代以来，网格生成技术得到了飞速的发展，从一开始的只能生成二维简单外形的网格，到如今能生成全尺寸飞行器的千万、甚至上亿数量网格。近年来，随着以“天河二号”、“神威.太湖之光”等为代表的超级高性能计算机的发展，国内外的计算机性能显著提高，在为CFD技术的发展提供充足动力的同时也带来了挑战，即这些超级计算机无一例外都要求计算软件具有并行化处理能力。遗憾的是，网格生成技术一直受限于其算法的特殊性，迄今为止，不管是商业软件还是国外的研究型代码，绝大部分都是串行程序。换句话说，目前在国内外都没有成熟的并行化网格生成技术。
对网格生成技术的并行化技术，目前还处于探索阶段，国外已经开始有一些相关文献可查。然而，这些网格生成技术都仅仅是只能并行生成某些特定类型的单元。例如，一些并行化算法通过对传统的Denaulay方法进行改造，只能并行生成四面体单元；另外一些是对传统的层推进(阵面推进)法进行改造后，只能生成棱柱单元。总之，目前还没有一种能适应于任意网格类型的并行化网格生成技术。此外，这些并行化算法还面临一个巨大的问题，即在给定边界条件下，难以生成上百亿量级的超大规模网格。
发明内容
为了克服现有技术的上述缺点，本发明提出了一种基于全局加密的超大规模非结构网格生成方法，可以利用千核量级的计算机资源，生成上百亿量级的任意类型网格单元(如四面体/三棱柱/六面体/金字塔等)。
本发明解决其技术问题所采用的技术方案是：一种基于全局加密的超大规模非结构网格生成方法，包括如下步骤：
步骤一、在初始网格单元的基础上，将所有网格单元剖分加密；
步骤二、对剖分加密后的网格单元进行物面保形；
步骤三、采用并行化的径向基函数方法实现空间点变形。
与现有技术相比，本发明的积极效果是：
1、由于采用了基于全局初始网格加密的方式，因此可以并行化生成任意类型的网格单元；
2、只要有足够的计算机资源，就可以不断加密，生成任意超大规模的网格；
3、本发明的方法是在初始网格上不断“全局加密-物面投影-空间变形”生成，因此效率极高。
附图说明
本发明将通过例子并参照附图的方式说明，其中：
图1为本发明方法的流程图；
图2为全局网格剖分加密示意图；
图3为物面保形示意图；
图4为空间点变形示意图；
图5为并行化的径向基函数方法的流程图。
具体实施方式
一种基于全局加密的超大规模非结构网格生成方法，如图1所示，包括如下步骤：
步骤一、全局网格剖分加密：在初始网格单元的基础上，将所有网格单元剖分加密。
如图2所示，图中的加粗线条表示并行计算过程中的分区的边界，其中初始网格是分区后被分配至每个进程的，每个进程上有1个或多个子区；剖分加密按照CGNS数据结构中给定的剖分方式进行，例如，将1个六面体剖分为8个六面体，1个三棱柱剖分为6个三棱柱。这样的剖分方式能基于任意的初始网格实现大规模网格的并行生成。
步骤二、物面保形：
进行物面保形前后的图如图3所示，图中以圆形代表物体外形，加粗线条表示并行计算过程中的分区的边界，保形前物面是由线段构成的，保形后物面是完全真实的外形。将步骤一中物面边界上新添加的不处于CAD几何上的点pi做保形处理，使点pi投影到CAD几何上得到新投影的物面点pi’。并且物面点的投影是并行处理的，每个进程分别读入数模，并将各自的子区上的物面点投影。物面保形调用了开源软件OpenCascade将物面上新添加的点投影到几何上，使得新的网格的物面边界能完全真实地代表实际的物体外形。
添加点的方式是按照步骤一中的方法，根据CGNS数据结构(行业标准)，直接在每条线段中点上添加新的点。
步骤三、采用并行化的径向基函数方法实现空间点变形：
如图4所示，图中的加粗线条表示并行计算过程中的分区的边界，将与新投影的物面点pi’相交的空间点变形。
径向基函数方法(Radial Basis Function，以下简称为“RBF方法”)，是一类利用距离基函数进行数据插值的方法，在计算数学、数值模拟等工程领域有广泛的应用。例如，在进行气动/热/结构多场耦合模拟中，用于在空间单元和物面单元之间的流场变量插值；在弹性变形、多体分离等动边界问题中，通过将物面上的点的变形插值到空间，实现空间网格变形。
RBF方法的基本原理为：
这里，N是控制点子集中点的数目(控制点子集由离空间点最近的多个物面变形点构成)，f(r)为空间中点Pfield,i的位移量，r为Pfield,i的位置，ri为第i个物面控制点Psurface,i的初始位置，||r-ri||表示Pfield,i和Psurface,i之间的欧式距离，表示所选用的径向基函数，wi表示第i个物面控制点Psurface,i处的权重系数,以x方向为例，该系数满足如下关系式：
其中，△xi为第i个物面点Psurface,i的变形量，基函数一般选用C2函数定义为：
通过式(1)，将物面点的位移插值到空间，得到空间点的变形量。
虽然径向基函数方法已经在工程领域得到了大量使用，但是迄今为止，几乎绝大部分的应用都是基于串行的方法，这是由于RBF方法中涉及到式(2)中的矩阵运算，而该式右端项中的矩阵由全局计算域中的物面点构成，使得RBF方法难以并行化。虽然目前在学术界已经有一些初步的并行化算法，但是一般都是采用PETs、LAPACK等第三方数学库，对式(2)中的矩阵进行并行化地求逆操作，存在健壮性差、效率低、通用性差的问题。
如图5所示，并行化的径向基函数方法包括如下内容：
(1)建立并行架构：
1)使用常规的分区算法(如调用METIS)，将计算域的网格划分为n个分区块Z＝{Z1、Z2、Z3······Zn}，每个区块对应一个进程P＝{P1、P2、P3······Pn}；这里的分区块数n由初始网格规模和计算机内存决定，一般使得分区后每个分区内的网格单元数少于100万。
2)将每个分区块Zi上的物面单元发送至对应的服务器进程Pi；
3)服务器进程将收集到的所有物面单元S＝{s1、s2、s3······sl组装为整个物面S；这里，l<＝n。
4)服务器进程将包含所有物面单元S＝{s1、s2、s3······sl}的整个物面发送给每一个进程Pi，使得每个进程Pi都有一个全局物面单元S的备份；
(2)遍历空间的所有点，进行并行RBF变形。对于每个进程Pi中的每一个空间点Pfield,i变形：
1)在整个物面单元集S上的物面控制点中，查找离空间点Pfield,i最近的N个控制点(N一般取10～30)，构成Pfield,i的控制点子集{ci}；；
2)对空间点Pfield,i，基于其拥有的控制点子集{ci}和这些控制点的变形量，根据式(1)、式(2)，求解得到空间点Pfield,i的变形。
下表为本实施例与其他现有方法的比较结果：

Claims (7)

1.一种基于全局加密的超大规模非结构网格生成方法，其特征在于：包括如下步骤：
步骤一、在初始网格单元的基础上，将所有网格单元剖分加密；
步骤二、对剖分加密后的网格单元进行物面保形；
步骤三、采用并行化的径向基函数方法实现空间点变形。
2.根据权利要求1所述的基于全局加密的超大规模非结构网格生成方法，其特征在于：步骤二所述对剖分加密后的网格单元进行物面保形的方法为：将经过剖分加密后物面边界上新添加的不处于CAD几何上的点pi投影到CAD几何上得到新投影的物面点pi’。
3.根据权利要求2所述的基于全局加密的超大规模非结构网格生成方法，其特征在于：所述物面点的投影是并行处理的，由每个进程分别读入数模，并将各自的子区上的物面点投影。
4.根据权利要求1所述的基于全局加密的超大规模非结构网格生成方法，其特征在于：在对初始网格单元进行剖分加密时在物面边界上添加点的方法为直接在每条线段中点或者面心上添加新的点。
5.根据权利要求1所述的基于全局加密的超大规模非结构网格生成方法，其特征在于：步骤三所述采用并行化的径向基函数方法实现空间点变形的方法为：(1)建立并行架构，将计算域的网格进行分区，每个分区块对应一个进程，并使每个进程均有一个全局物面单元点的备份；(2)遍历每个进程中的所有空间点，进行并行RBF变形，其中，各空间点以离其最近的多个物面变形点为控制点子集，执行串行RBF变形，求解得到空间点的变形量。
6.根据权利要求5所述的基于全局加密的超大规模非结构网格生成方法，其特征在于：求解空间点变形量的方法如下：按如下公式计算空间点的位移量：
其中：N是控制点子集中点的数目，f(r)为空间点Pfield,i的位移量，r为空间点Pfield,i的位置，ri为第i个物面控制点Psurface,i的初始位置，||r-ri||表示Pfield,i和Psurface,i之间的欧式距离，表示所选用的径向基函数，wi表示第i个物面控制点Psurface,i处的权重系数。
7.根据权利要求6所述的基于全局加密的超大规模非结构网格生成方法，其特征在于：所述权重系数wi满足如下关系式：
其中，△xi为第i个物面控制点Psurface,i的变形量，
CN201710655494.0A 2017-08-03 2017-08-03 基于全局加密的超大规模非结构网格生成方法 Pending CN107403466A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710655494.0A CN107403466A (zh) 2017-08-03 2017-08-03 基于全局加密的超大规模非结构网格生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710655494.0A CN107403466A (zh) 2017-08-03 2017-08-03 基于全局加密的超大规模非结构网格生成方法

Publications (1)

Publication Number Publication Date
CN107403466A true CN107403466A (zh) 2017-11-28

Family

ID=60402659

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710655494.0A Pending CN107403466A (zh) 2017-08-03 2017-08-03 基于全局加密的超大规模非结构网格生成方法

Country Status (1)

Country Link
CN (1) CN107403466A (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109783965A (zh) * 2019-01-25 2019-05-21 中国空气动力研究与发展中心计算空气动力研究所 一种结构网格自动分块加密方法
CN110334450A (zh) * 2019-07-09 2019-10-15 中国空气动力研究与发展中心计算空气动力研究所 一种多块结构网格生成中物面投影错误的修复方法
CN111125949A (zh) * 2019-12-06 2020-05-08 北京科技大学 一种有限元分析的大规模并行网格划分系统及方法
CN111400969A (zh) * 2020-03-12 2020-07-10 中国航天空气动力技术研究院 一种非结构直角网格加速生成方法
CN111538487A (zh) * 2020-04-17 2020-08-14 中国空气动力研究与发展中心计算空气动力研究所 分布式并行网格生成软件框架
CN111814246A (zh) * 2020-06-11 2020-10-23 中国空气动力研究与发展中心计算空气动力研究所 一种基于生成对抗网络的翼型反设计方法
CN112463360A (zh) * 2020-10-29 2021-03-09 空气动力学国家重点实验室 一种百亿百gb量级网格数据文件的并行读入方法
CN112632869A (zh) * 2020-12-23 2021-04-09 中国空气动力研究与发展中心计算空气动力研究所 一种基于网格框架的非结构附面层网格生成方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101013454A (zh) * 2007-02-02 2007-08-08 郑州机械研究所 Cae软件系统网格剖分的智能化方法
CN105302951A (zh) * 2015-10-20 2016-02-03 北京工业大学 一种有限元网格曲面剖分方法
CN105550424A (zh) * 2015-12-08 2016-05-04 中国航空工业集团公司西安航空计算技术研究所 一种基于rbf网格变形插值序列的筛选方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101013454A (zh) * 2007-02-02 2007-08-08 郑州机械研究所 Cae软件系统网格剖分的智能化方法
CN105302951A (zh) * 2015-10-20 2016-02-03 北京工业大学 一种有限元网格曲面剖分方法
CN105550424A (zh) * 2015-12-08 2016-05-04 中国航空工业集团公司西安航空计算技术研究所 一种基于rbf网格变形插值序列的筛选方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
PRADEEP GOPALAKRISHNAN ET AL: "A parallel boundary fitted dynamic mesh solver for applications to flapping flight", 《COMPUTERS & FLUIDS》 *
唐静: "非结构网格自适应技术研究与应用", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
孙岩等: "基于径向基函数与混合背景网格的动态网格变形方法", 《航空学报》 *
马文鹏等: "径向基函数在动网格中的应用及可并行性研究", 《科研信息化技术与应用》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109783965B (zh) * 2019-01-25 2022-08-02 中国空气动力研究与发展中心计算空气动力研究所 一种结构网格自动分块加密方法
CN109783965A (zh) * 2019-01-25 2019-05-21 中国空气动力研究与发展中心计算空气动力研究所 一种结构网格自动分块加密方法
CN110334450B (zh) * 2019-07-09 2021-06-15 中国空气动力研究与发展中心计算空气动力研究所 一种多块结构网格生成中物面投影错误的修复方法
CN110334450A (zh) * 2019-07-09 2019-10-15 中国空气动力研究与发展中心计算空气动力研究所 一种多块结构网格生成中物面投影错误的修复方法
CN111125949A (zh) * 2019-12-06 2020-05-08 北京科技大学 一种有限元分析的大规模并行网格划分系统及方法
CN111400969A (zh) * 2020-03-12 2020-07-10 中国航天空气动力技术研究院 一种非结构直角网格加速生成方法
CN111400969B (zh) * 2020-03-12 2022-10-28 中国航天空气动力技术研究院 一种非结构直角网格加速生成方法
CN111538487B (zh) * 2020-04-17 2022-06-14 中国空气动力研究与发展中心计算空气动力研究所 分布式并行网格生成软件框架
CN111538487A (zh) * 2020-04-17 2020-08-14 中国空气动力研究与发展中心计算空气动力研究所 分布式并行网格生成软件框架
CN111814246A (zh) * 2020-06-11 2020-10-23 中国空气动力研究与发展中心计算空气动力研究所 一种基于生成对抗网络的翼型反设计方法
CN111814246B (zh) * 2020-06-11 2022-09-27 中国空气动力研究与发展中心计算空气动力研究所 一种基于生成对抗网络的翼型反设计方法
CN112463360A (zh) * 2020-10-29 2021-03-09 空气动力学国家重点实验室 一种百亿百gb量级网格数据文件的并行读入方法
CN112632869A (zh) * 2020-12-23 2021-04-09 中国空气动力研究与发展中心计算空气动力研究所 一种基于网格框架的非结构附面层网格生成方法

Similar Documents

Publication Publication Date Title
CN107403466A (zh) 基于全局加密的超大规模非结构网格生成方法
Feng et al. Three-dimensional solar wind modeling from the Sun to Earth by a SIP-CESE MHD model with a six-component grid
CN112862972B (zh) 一种表面结构网格生成方法
Cai et al. Adaptive radial-basis-function-based multifidelity metamodeling for expensive black-box problems
Xiong et al. Efficient parallel implementation of the lattice Boltzmann method on large clusters of graphic processing units
CN116151084B (zh) 基于结构网格的模拟方法、装置、终端设备及存储介质
CN108228970B (zh) 结构动力学分析显式异步长并行计算方法
Gang et al. Mesh deformation on 3D complex configurations using multistep radial basis functions interpolation
CN106294975B (zh) 一种基于降阶模型的梁式结构自由振动分析方法
Zhang et al. Adaptive ANCF method and its application in planar flexible cables
CN103345580A (zh) 基于格子Boltzmann方法的并行CFD方法
CN104360896A (zh) 一种基于gpu集群的并行流体仿真加速方法
CN112632874A (zh) 一种直升机流场数值模拟的优化方法及系统
CN110610042A (zh) 基于阵面推进的非结构网格结构化处理方法
CN111125963A (zh) 基于拉格朗日积分点有限元的数值仿真系统及方法
Mueller‐Roemer et al. Ternary sparse matrix representation for volumetric mesh subdivision and processing on GPUs
CN112001109A (zh) 再生核粒子算法实现结构冲击动力学仿真方法
Negi et al. Algorithms for uniform particle initialization in domains with complex boundaries
Deng et al. Multimoment finite volume solver for euler equations on unstructured grids
CN107563080A (zh) 基于gpu的两相介质随机模型并行生成方法、电子设备
CN104992425A (zh) 一种基于gpu加速的dem超分辨率方法
Gao et al. A multi-level parallel tie-dye algorithm for auto-CFD
Xu et al. Parallelizing a high-order CFD software for 3D, multi-block, structural grids on the TianHe-1A supercomputer
CN114925627A (zh) 一种基于图形处理器的直升机流场数值模拟系统及方法
CN114611423A (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20171128

RJ01 Rejection of invention patent application after publication