CN114492224A - 喷流干扰流场混合网格计算方法、系统、设备及介质 - Google Patents

喷流干扰流场混合网格计算方法、系统、设备及介质 Download PDF

Info

Publication number
CN114492224A
CN114492224A CN202111614125.XA CN202111614125A CN114492224A CN 114492224 A CN114492224 A CN 114492224A CN 202111614125 A CN202111614125 A CN 202111614125A CN 114492224 A CN114492224 A CN 114492224A
Authority
CN
China
Prior art keywords
grid
structural
flow field
sub
transition layer
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
CN202111614125.XA
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.)
China Academy of Aerospace Aerodynamics CAAA
Original Assignee
China Academy of Aerospace Aerodynamics CAAA
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 China Academy of Aerospace Aerodynamics CAAA filed Critical China Academy of Aerospace Aerodynamics CAAA
Priority to CN202111614125.XA priority Critical patent/CN114492224A/zh
Publication of CN114492224A publication Critical patent/CN114492224A/zh
Pending legal-status Critical Current

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/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • 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]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种喷流干扰流场混合网格计算方法、系统、设备及介质。方法包括：S1:将计算域划分为结构求解区域和非结构求解区域；S2:生成结构子网格和非结构子网格；S3:选取交界处的部分结构子网格作为过渡层；S4:标记过渡层的网格单元在结构子网格和非结构子网格之间的对应关系；S5:将过渡层的结构网格与非结构子网格一并输出为非结构网格；S6:调用结构求解器和非结构求解器分别进行CFD求解计算；S7:根据S4标记的对应关系，将过渡层的结构网格的流场信息与过渡层的非结构网格的流场信息进行数据交换，并更新结构网格的流场信息和非结构网格的流场信息；S8:判断是否收敛，若是，则计算完成，否则，返回S6。实现提高喷流干扰流场计算的精度及效率。

Description

喷流干扰流场混合网格计算方法、系统、设备及介质
技术领域
本发明属于计算流体动力学领域，更具体地，涉及一种喷流干扰流场混合网格计算方法、系统、设备及介质。
背景技术
反作用控制系统(RCS)利用发动机喷流产生反作用力快速改变飞行器运动姿态或轨道的直接力控制系统，其作用在于补充气动舵面效率不足和快速改变飞行状态，具有广阔的工程应用前景。RCS喷流进入超/高超声速外流时将产生相互作用，形成包含边界层分离与再附、激波、膨胀波、马赫盘、剪切层等多种流动现象的复杂的激波/边界层干扰流场，产生随喷流参数、飞行条件、布局形式等强非线性变化的气动力/热干扰。
工程应用中，要求尽可能准确模拟这种复杂干扰流场，并提取因干扰产生的气动力/热变化量作为控制系统设计的输入，保证设计精度。因此，对于喷流干扰流场一般采用结构网格进行模拟，结构网格具有贴体性，可以精确、高效的模拟边界层、间断等，计算精度高。然而结构网格对复杂几何构型的适应能力较差，网格生成难度大、耗时长，大大限制了计算效率，甚至因某些过于复杂的构型无法直接生成结构化网格或网格质量太差，而不得不进行简化处理。
非结构网格最大的优点是几何适应能力强、网格生成效率高，与结构网格混合起来可以用于弥补结构网格在处理前述复杂构型问题上的不足。
对于复杂外形飞行器的数值模拟，综合结构网格和非结构网格优势的混合网格技术代表了未来网格技术的发展趋势。然而，当前绝大多数基于混合网格的CFD软件一般将结构网格视为非结构网格的特例，即将结构网格转化为非结构网格的数据结构，进而统一运用非结构网格求解器进行数值计算，这实际上没有充分发挥结构网格的优势，这种仅网格层面的混合方法不适用于精度要求高的喷流干扰流场计算。
发明内容
本发明的目的是提出一种喷流干扰流场混合网格计算方法、系统、设备及介质，实现提高喷流干扰流场计算的精度及效率。
第一方面，本发明提出一种喷流干扰流场混合网格计算方法，包括：
步骤S1:将喷流干扰流场的计算域划分为结构求解区域和非结构求解区域；
步骤S2:将所述结构求解区域生成结构子网格，将所述非结构求解区域生成非结构子网格；
步骤S3:选取位于所述结构子网格与所述非结构子网格交界处的部分结构子网格作为过渡层；
步骤S4:标记所述过渡层的网格单元在所述结构子网格和所述非结构子网格之间的对应关系；
步骤S5:将所述过渡层的结构网格与所述非结构子网格一并输出为非结构网格，将所述结构子网格输出为结构网格；
步骤S6:调用结构求解器对结构网格进行CFD求解计算，获得结构网格的流场信息，并调用非结构求解器对非结构网格进行CFD求解计算，获得非结构网格的流场信息；
步骤S7:根据步骤S4标记的对应关系，将所述结构解算器计算得到的过渡层的结构网格的流场信息与所述非结构解算器计算得到的过渡层的非结构网格的流场信息进行数据交换，并更新所述结构网格的流场信息和所述非结构网格的流场信息；
步骤S8:判断所述结构网格及所述非结构网格的CFD计算结果是否收敛，若是，则计算完成，否则，返回步骤S6。
可选地，所述步骤S1具体包括：
根据喷流干扰流场的计算域外形，将所述计算域划分为非结构求解区域和结构求解区域。
可选地，所述步骤S3具体包括：
以所述结构子网格与所述非结构子网格的交接面为起始，向所述结构子网格内扩展设定层数的结构网格作为过渡层。
可选地，所述步骤S7具体包括:
根据步骤S4标记的对应关系，将所述结构解算器计算得到的过渡层的结构网格的流场信息赋值于过渡层对应的非结构网格，并更新所述非结构网格的流场信息；
根据步骤S4标记的对应关系，将所述非结构解算器计算得到的过渡层的非结构网格的流场信息赋值于过渡层对应的结构网格，并更新所述结构网格的流场信息。
第二方面，本发明提出一种电子设备，所述电子设备包括：
至少一个处理器；以及，
与所述至少一个处理器通信连接的存储器；其中，
所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行第一方面所述的喷流干扰流场混合网格计算方法。
第三方面，本发明提出一种非暂态计算机可读存储介质，该非暂态计算机可读存储介质存储计算机指令，该计算机指令用于使计算机执行第一方面所述的的喷流干扰流场混合网格计算方法。
第四方面，本发明提出一种喷流干扰流场混合网格计算系统，包括：
网格生成模块，用于将喷流干扰流场的计算域划分为结构求解区域和非结构求解区域，并将所述结构求解区域生成结构子网格，将所述非结构求解区域生成非结构子网格；
过渡层选取模块，用于选取位于所述结构子网格与所述非结构子网格交界处的部分结构子网格作为过渡层；
标记模块，用于标记所述过渡层的网格单元在所述结构子网格和所述非结构子网格之间的对应关系；
网格输出模块，用于将所述过渡层的结构网格与所述非结构子网格一并输出为非结构网格，将所述结构子网格输出为结构网格；
结构求解器，用于对结构网格进行CFD求解计算；
非结构求解器，用于对非结构网格进行CFD求解计算；
数据交换模块，用于根据所述过渡层选取模块标记的对应关系，将所述结构解算器计算得到的过渡层的结构网格的流场信息与所述非结构解算器计算得到的过渡层的非结构网格的流场信息进行数据交换；
流场信息计算模块，用于执行以下步骤：
步骤S11:调用结构求解器对结构网格进行CFD求解计算，获得结构网格的流场信息，并调用非结构求解器对非结构网格进行CFD求解计算，获得非结构网格的流场信息；
步骤S12:根据所述标记模块标记的对应关系，通过所述数据交换模块将所述结构解算器计算得到的过渡层的结构网格的流场信息与所述非结构解算器计算得到的过渡层的非结构网格的流场信息进行数据交换，并更新所述结构网格的流场信息和所述非结构网格的流场信息；
步骤S13:判断所述结构网格及所述非结构网格的CFD计算结果是否收敛，若是，则计算完成，否则，返回步骤S11。
可选地，所述网格生成模块具体用于：
根据喷流干扰流场的计算域外形，将所述计算域划分为非结构求解区域和结构求解区域，并将所述结构求解区域生成结构子网格，将所述非结构求解区域生成非结构子网格。
可选地，所述过渡层选取模块具体用于：
以所述结构子网格与所述非结构子网格的交接面为起始，向所述结构子网格内扩展并选取设定层数的结构网格作为过渡层。
可选地，所述数据交换模块具体用于:
根据所述标记模块标记的对应关系，将所述结构解算器计算得到的过渡层的结构网格的流场信息赋值于过渡层对应的非结构网格；
根据所述标记模块标记的对应关系，将所述非结构解算器计算得到的过渡层的非结构网格的流场信息赋值于过渡层对应的结构网格。
本发明的有益效果在于：
本发明通过选取位于结构子网格与非结构子网格交界处的部分结构子网格作为过渡层，并标记过渡层的网格单元在结构子网格和非结构子网格之间的对应关系，然后将过渡层的结构网格与非结构子网格一并输出为非结构网格，并将结构子网格输出为结构网格，之后调用结构求解器和非结构求解器分别对结构网格及非结构网格进行CFD求解计算，获得结构网格的流场信息和非结构网格的流场信息，然后基于之前标记的对应关系，将结构解算器计算得到的过渡层的结构网格的流场信息与非结构解算器计算得到的过渡层的非结构网格的流场信息进行数据交换，并更新结构网格的流场信息和非结构网格的流场信息，反复迭代计算直到计算结果收敛，从而实现结构/非结构求解器层面的耦合计算，相较于现有技术，本发明的方法不需要将结构求解器与非结构求解器进行统一架构，两者可独立，程序改动少，易实现，扩展性好，结构求解器与非结构求解器各自具有很大独立性，求解器内部可根据需求自由选择离散格式、湍流模型等算法组合，或对现有算法进行更新替换，而不影响耦合计算，同时兼顾了计算精度与效率，适用于喷流干扰流场高精度模拟。
本发明的系统具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据本发明实施例1的一种喷流干扰流场混合网格计算方法的流程图。
图2示出了根据本发明实施例1的一种喷流干扰流场混合网格计算方法中过渡层网格的示意图。
图3示出了根据本发明实施例1的一种喷流干扰流场混合网格计算方法中结构网格和非结构网格的流场信息数据交换示意图。
图4示出了根据本发明实施例1的一具体应用示例中旋成体/舵面外形侧向喷流干扰流场结构/非结构混合计算网格图。
图5示出了根据本发明实施例1的一具体应用示例中结构/非结构混合网格计算结果与纯结构网格计算结果的对比示意图。
具体实施方式
针对现有技术存在的问题，一种最佳的选择是在结构网格上运行结构解算器，而在非结构网格上运行非结构解算器，在算法层次上实现真正的结构/非结构网格耦合计算。
对此，本发明提供一种用于喷流干扰流场的结构/非结构算法层次的混合计算方法，非结构网格仅用于处理计算域中局部复杂部件，除此之外的计算域全部用结构网格处理，最大范围利用结构网格以保证计算精度，并严格保证结构网格、非结构网格交界面的通量守恒,避免插值误差。
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
实施例1
图1示出了根据本发明实施例1的一种喷流干扰流场混合网格计算方法的流程图。
如图1所示，一种喷流干扰流场混合网格计算方法，包括：
步骤S101:将喷流干扰流场的计算域划分为结构求解区域和非结构求解区域；
在一具体应用场景中，根据喷流干扰流场的计算域外形，将计算域划分为非结构求解区域和结构求解区域。
步骤S102:将结构求解区域生成结构子网格，将非结构求解区域生成非结构子网格；
在上述具体应用场景中，根据计算域外形，针对复杂部件生成非结构网格，其余生成结构网格，将结构与非结构网格交接处网格面定义为交接面。
步骤S103:选取位于结构子网格与非结构子网格交界处的部分结构子网格作为过渡层；
在上述具体应用场景中，如图2所示，以结构子网格与非结构子网格交接面为起始，向结构子网格内层扩展一定层数的网格作为过渡层，用于数据交换，过渡层作为非结构子网格的一部分。
步骤S104:标记过渡层的网格单元在结构子网格和非结构子网格之间的对应关系；
在上述具体应用场景中，结构和非结构子网格在过渡层重合，根据重合的网格单元中心坐标相等的条件，找出并标记过渡层网格单元在结构和非结构子网格之间的对应关系。
步骤S105:将过渡层的结构网格与非结构子网格一并输出为非结构网格，将结构子网格输出为结构网格；
在上述具体应用场景中，将过渡层所在的结构网格与原非结构网格部分一并按非结构网格输出，用于非结构求解器读取与计算，将原结构网格部分输出用于结构求解器读取与计算。
步骤S106:调用结构求解器对结构网格进行CFD求解计算，获得结构网格的流场信息，并调用非结构求解器对非结构网格进行CFD求解计算，获得非结构网格的流场信息；
在上述具体应用场景中，开始迭代计算，调用结构求解器和非结构求解器分别对结构网格和非结构网格同步进行CFD求解。
步骤S107:根据步骤S4标记的对应关系，将结构解算器计算得到的过渡层的结构网格的流场信息与非结构解算器计算得到的过渡层的非结构网格的流场信息进行数据交换，并更新结构网格的流场信息和非结构网格的流场信息；
在上述具体应用场景中，如图3所示，结构网格和非结构网格同步处理边界条件，分别获取结构子网格和非结构子网格在过渡层的流场信息，然后，根据步骤S104标记的对应关系，将结构解算器计算得到的过渡层内层结构网格的流场信息赋值于过渡层对应的非结构网格，并更新非结构网格的流场信息；根据步骤S104标记的对应关系，将非结构解算器计算得到的过渡层内层非结构网格的流场信息赋值于过渡层对应的结构网格，并更新结构网格的流场信息。
具体实施过程中，可以在结构求解器和非结构求解器外部添加数据交换模块，用于过渡层数据的相互交换，从而实现结构/非结构求解器层面的耦合计算。
步骤S108:判断结构网格及非结构网格的CFD计算结果是否收敛，若是，则计算完成，否则，返回步骤S106。
下面通过一具体应用示例对本发明的方法做进一步的解释。
应用示例1的混合网格计算方法，针对旋成体外形侧向喷流干扰流场开展了结构/非结构混合网格计算测试，来流马赫数为4.5，来流总压为1491000Pa，来流总温为306K，喷管出口马赫数为1，喷流总压为3877000Pa，喷流总温为306K，旋成体/舵面外形侧向喷流干扰流场结构/非结构混合计算网格如图4所示。
计算结果与采用纯结构网格的计算结果进行了对比，如图5所示，图中，P/P∞为压强与来流静压之比；X/L为x轴坐标与弹长之比，L的长度为290mm；实线为纯结构网格计算结果，虚线为混合网格计算结果。可以看出混合网格计算结果与纯结构网格计算结果基本一致，从而验证了本发明的混合网格计算方法用于计算喷流干扰流场的有效性。
综上，本发明的喷流干扰流场混合网格计算方法具有以下优点：
(1)不需要将结构求解器与非结构求解器进行统一架构，两者可独立；
(2)程序改动少，易实现，只需添加数据交换模块；
(3)扩展性好，结构求解器与非结构求解器各自具有很大独立性，求解器内部可根据需求自由选择离散格式、湍流模型等算法组合，或对现有算法进行更新替换，而不影响耦合计算；
(4)兼顾了计算精度与效率，适用于喷流干扰流场高精度模拟。
实施例2
一种喷流干扰流场混合网格计算系统，包括：
网格生成模块，用于将喷流干扰流场的计算域划分为结构求解区域和非结构求解区域，并将结构求解区域生成结构子网格，将非结构求解区域生成非结构子网格；
过渡层选取模块，用于选取位于结构子网格与非结构子网格交界处的部分结构子网格作为过渡层；
标记模块，用于标记过渡层的网格单元在结构子网格和非结构子网格之间的对应关系；
网格输出模块，用于将过渡层的结构网格与非结构子网格一并输出为非结构网格，将结构子网格输出为结构网格；
结构求解器，用于对结构网格进行CFD求解计算；
非结构求解器，用于对非结构网格进行CFD求解计算；
数据交换模块，用于根据过渡层选取模块标记的对应关系，将结构解算器计算得到的过渡层的结构网格的流场信息与非结构解算器计算得到的过渡层的非结构网格的流场信息进行数据交换；
流场信息计算模块，用于执行以下步骤：
步骤S11:调用结构求解器对结构网格进行CFD求解计算，获得结构网格的流场信息，并调用非结构求解器对非结构网格进行CFD求解计算，获得非结构网格的流场信息；
步骤S12:根据标记模块标记的对应关系，通过数据交换模块将结构解算器计算得到的过渡层的结构网格的流场信息与非结构解算器计算得到的过渡层的非结构网格的流场信息进行数据交换，并更新结构网格的流场信息和非结构网格的流场信息；
步骤S13:判断结构网格及非结构网格的CFD计算结果是否收敛，若是，则计算完成，否则，返回步骤S11。
本实施例中，网格生成模块具体用于：
根据喷流干扰流场的计算域外形，将计算域划分为非结构求解区域和结构求解区域，并将结构求解区域生成结构子网格，将非结构求解区域生成非结构子网格。
本实施例中，过渡层选取模块具体用于：
以结构子网格与非结构子网格的交接面为起始，向结构子网格内扩展并选取设定层数的结构网格作为过渡层。
本实施例中，数据交换模块具体用于:
根据标记模块标记的对应关系，将结构解算器计算得到的过渡层的结构网格的流场信息赋值于过渡层对应的非结构网格；
根据标记模块标记的对应关系，将非结构解算器计算得到的过渡层的非结构网格的流场信息赋值于过渡层对应的结构网格。
实施例3
一种电子设备，所述电子设备包括：
至少一个处理器；以及，
与所述至少一个处理器通信连接的存储器；其中，
所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行实施例1所述的喷流干扰流场混合网格计算方法。
根据本公开实施例的电子设备包括存储器和处理器，该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(RAM)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(ROM)、硬盘、闪存等。
该处理器可以是中央处理单元(CPU)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。
本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。
有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。
实施例4
一种非暂态计算机可读存储介质，该非暂态计算机可读存储介质存储计算机指令，该计算机指令用于使计算机执行实施例1所述的的喷流干扰流场混合网格计算方法。
根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。
上述计算机可读存储介质包括但不限于：光存储介质(例如：CD－ROM和DVD)、磁光存储介质(例如：MO)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置ROM的媒体(例如：ROM盒)。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。

Claims (10)

1.一种喷流干扰流场混合网格计算方法，其特征在于，包括：
步骤S1:将喷流干扰流场的计算域划分为结构求解区域和非结构求解区域；
步骤S2:将所述结构求解区域生成结构子网格，将所述非结构求解区域生成非结构子网格；
步骤S3:选取位于所述结构子网格与所述非结构子网格交界处的部分结构子网格作为过渡层；
步骤S4:标记所述过渡层的网格单元在所述结构子网格和所述非结构子网格之间的对应关系；
步骤S5:将所述过渡层的结构网格与所述非结构子网格一并输出为非结构网格，将所述结构子网格输出为结构网格；
步骤S6:调用结构求解器对结构网格进行CFD求解计算，获得结构网格的流场信息，并调用非结构求解器对非结构网格进行CFD求解计算，获得非结构网格的流场信息；
步骤S7:根据步骤S4标记的对应关系，将所述结构解算器计算得到的过渡层的结构网格的流场信息与所述非结构解算器计算得到的过渡层的非结构网格的流场信息进行数据交换，并更新所述结构网格的流场信息和所述非结构网格的流场信息；
步骤S8:判断所述结构网格及所述非结构网格的CFD计算结果是否收敛，若是，则计算完成，否则，返回步骤S6。
2.根据权利要求1所述的喷流干扰流场混合网格计算方法，其特征在于，所述步骤S1具体包括：
根据喷流干扰流场的计算域外形，将所述计算域划分为非结构求解区域和结构求解区域。
3.根据权利要求1所述的喷流干扰流场混合网格计算方法，其特征在于，所述步骤S3具体包括：
以所述结构子网格与所述非结构子网格的交接面为起始，向所述结构子网格内扩展设定层数的结构网格作为过渡层。
4.根据权利要求1所述的喷流干扰流场混合网格计算方法，其特征在于，所述步骤S7具体包括:
根据步骤S4标记的对应关系，将所述结构解算器计算得到的过渡层的结构网格的流场信息赋值于过渡层对应的非结构网格，并更新所述非结构网格的流场信息；
根据步骤S4标记的对应关系，将所述非结构解算器计算得到的过渡层的非结构网格的流场信息赋值于过渡层对应的结构网格，并更新所述结构网格的流场信息。
5.一种电子设备，其特征在于，所述电子设备包括：
至少一个处理器；以及，
与所述至少一个处理器通信连接的存储器；其中，
所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1-4任一所述的喷流干扰流场混合网格计算方法。
6.一种非暂态计算机可读存储介质，其特征在于，该非暂态计算机可读存储介质存储计算机指令，该计算机指令用于使计算机执行权利要求1-4任一所述的的喷流干扰流场混合网格计算方法。
7.一种喷流干扰流场混合网格计算系统，其特征在于，包括：
网格生成模块，用于将喷流干扰流场的计算域划分为结构求解区域和非结构求解区域，并将所述结构求解区域生成结构子网格，将所述非结构求解区域生成非结构子网格；
过渡层选取模块，用于选取位于所述结构子网格与所述非结构子网格交界处的部分结构子网格作为过渡层；
标记模块，用于标记所述过渡层的网格单元在所述结构子网格和所述非结构子网格之间的对应关系；
网格输出模块，用于将所述过渡层的结构网格与所述非结构子网格一并输出为非结构网格，将所述结构子网格输出为结构网格；
结构求解器，用于对结构网格进行CFD求解计算；
非结构求解器，用于对非结构网格进行CFD求解计算；
数据交换模块，用于根据所述过渡层选取模块标记的对应关系，将所述结构解算器计算得到的过渡层的结构网格的流场信息与所述非结构解算器计算得到的过渡层的非结构网格的流场信息进行数据交换；
流场信息计算模块，用于执行以下步骤：
步骤S11:调用结构求解器对结构网格进行CFD求解计算，获得结构网格的流场信息，并调用非结构求解器对非结构网格进行CFD求解计算，获得非结构网格的流场信息；
步骤S12:根据所述标记模块标记的对应关系，通过所述数据交换模块将所述结构解算器计算得到的过渡层的结构网格的流场信息与所述非结构解算器计算得到的过渡层的非结构网格的流场信息进行数据交换，并更新所述结构网格的流场信息和所述非结构网格的流场信息；
步骤S13:判断所述结构网格及所述非结构网格的CFD计算结果是否收敛，若是，则计算完成，否则，返回步骤S11。
8.根据权利要求1所述的喷流干扰流场混合网格计算系统，其特征在于，所述网格生成模块具体用于：
根据喷流干扰流场的计算域外形，将所述计算域划分为非结构求解区域和结构求解区域，并将所述结构求解区域生成结构子网格，将所述非结构求解区域生成非结构子网格。
9.根据权利要求1所述的喷流干扰流场混合网格计算系统，其特征在于，所述过渡层选取模块具体用于：
以所述结构子网格与所述非结构子网格的交接面为起始，向所述结构子网格内扩展并选取设定层数的结构网格作为过渡层。
10.根据权利要求1所述的喷流干扰流场混合网格计算系统，其特征在于，所述数据交换模块具体用于:
根据所述标记模块标记的对应关系，将所述结构解算器计算得到的过渡层的结构网格的流场信息赋值于过渡层对应的非结构网格；
根据所述标记模块标记的对应关系，将所述非结构解算器计算得到的过渡层的非结构网格的流场信息赋值于过渡层对应的结构网格。
CN202111614125.XA 2021-12-27 2021-12-27 喷流干扰流场混合网格计算方法、系统、设备及介质 Pending CN114492224A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111614125.XA CN114492224A (zh) 2021-12-27 2021-12-27 喷流干扰流场混合网格计算方法、系统、设备及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111614125.XA CN114492224A (zh) 2021-12-27 2021-12-27 喷流干扰流场混合网格计算方法、系统、设备及介质

Publications (1)

Publication Number Publication Date
CN114492224A true CN114492224A (zh) 2022-05-13

Family

ID=81496658

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111614125.XA Pending CN114492224A (zh) 2021-12-27 2021-12-27 喷流干扰流场混合网格计算方法、系统、设备及介质

Country Status (1)

Country Link
CN (1) CN114492224A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116451592A (zh) * 2023-06-14 2023-07-18 北京航天天美科技有限公司 一种可扩展智能方舱布局优化方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111199118A (zh) * 2019-12-13 2020-05-26 同济大学 非齐平双涵喷口喷流噪声计算方法、计算装置及存储介质
US20200410147A1 (en) * 2019-06-28 2020-12-31 Viettel Group Aerodynamic derivatives calculation method for flight vehicle
CN112287615A (zh) * 2020-10-27 2021-01-29 哈尔滨工程大学 一种复杂结构弹簧阀非定常流场多块结构网格划分方法
CN112699498A (zh) * 2021-03-23 2021-04-23 中国空气动力研究与发展中心计算空气动力研究所 基于归一化物理量间断特征的喷流模拟激波快速判别方法
CN112699623A (zh) * 2021-03-24 2021-04-23 中国空气动力研究与发展中心计算空气动力研究所 基于非结构网格规则化重构技术的高精度热流计算方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200410147A1 (en) * 2019-06-28 2020-12-31 Viettel Group Aerodynamic derivatives calculation method for flight vehicle
CN111199118A (zh) * 2019-12-13 2020-05-26 同济大学 非齐平双涵喷口喷流噪声计算方法、计算装置及存储介质
CN112287615A (zh) * 2020-10-27 2021-01-29 哈尔滨工程大学 一种复杂结构弹簧阀非定常流场多块结构网格划分方法
CN112699498A (zh) * 2021-03-23 2021-04-23 中国空气动力研究与发展中心计算空气动力研究所 基于归一化物理量间断特征的喷流模拟激波快速判别方法
CN112699623A (zh) * 2021-03-24 2021-04-23 中国空气动力研究与发展中心计算空气动力研究所 基于非结构网格规则化重构技术的高精度热流计算方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116451592A (zh) * 2023-06-14 2023-07-18 北京航天天美科技有限公司 一种可扩展智能方舱布局优化方法
CN116451592B (zh) * 2023-06-14 2023-08-15 北京航天天美科技有限公司 一种可扩展智能方舱布局优化方法

Similar Documents

Publication Publication Date Title
Nakahashi et al. Intergrid-boundary definition method for overset unstructured grid approach
CN110309543B (zh) 一种多介质流体界面运动的仿真流程设计方法
US20150127301A1 (en) Updating A CAD Model To Reflect Global Or Local Shape Changes
CN116245049B (zh) 节点式非结构网格的边界修正方法、装置、设备及介质
Chen et al. Tetrahedral mesh improvement by shell transformation
CN114492224A (zh) 喷流干扰流场混合网格计算方法、系统、设备及介质
Meng et al. Shape optimization of axisymmetric solids with the finite cell method using a fixed grid
CN116227043A (zh) 一种飞行器数值模拟方法、系统、设备及计算机存储介质
Ando et al. Automotive aerodynamic design exploration employing new optimization methodology based on CFD
CN112765871B (zh) 一种基于曲线坐标的并行粒子追踪方法和装置
KR102131305B1 (ko) 충돌 사건의 수치상 시뮬레이션에서 저차 및 2차 유한 요소들을 수반하는 접촉면 정의 생성
US8762917B2 (en) Automatically modifying a circuit layout to perform electromagnetic simulation
US11176739B1 (en) Enhancing tetrahedral meshing for self-intersecting boundaries
CN112182771B (zh) 基于数值模拟的数据处理方法、存储介质、电子装置
CN113515915B (zh) 插入填充单元的方法、装置、设备以及存储介质
CN110955380A (zh) 访存数据生成方法、存储介质、计算机设备和装置
Nakahashi et al. An automatic Euler solver using the unstructured upwind method
CN117763735A (zh) 一种结构-非结构混合数据交换方法、系统及设备
Zhou et al. Research on parametric modeling optimization of trajectory correction projectile based on fvm
CN105975683B (zh) 有限元网格模型更新的快速查找方法
Dokken et al. Requirements from isogeometric analysis for changes in product design ontologies
Bedregal et al. A study on size-optimal longest edge refinement algorithms
CN118313314B (zh) 粒子类模拟自由来流边界处理方法、装置、设备及介质
US12094059B1 (en) Anisotropic polyhedra boundary layer adaptation method
CN115544893A (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