CN115841063A - 用于动导数数值模拟的滑移动网格生成方法和装置 - Google Patents
用于动导数数值模拟的滑移动网格生成方法和装置 Download PDFInfo
- Publication number
- CN115841063A CN115841063A CN202310160122.6A CN202310160122A CN115841063A CN 115841063 A CN115841063 A CN 115841063A CN 202310160122 A CN202310160122 A CN 202310160122A CN 115841063 A CN115841063 A CN 115841063A
- Authority
- CN
- China
- Prior art keywords
- grid
- dynamic
- interface
- flow field
- area
- 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.)
- Granted
Links
Images
Classifications
-
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请提供一种用于动导数数值模拟的滑移动网格生成方法和装置,所述方法包括:基于预设数值模拟要求将目标计算区域划分为动网格区域和静网格区域,并分别生成所述动网格区域和所述静网格区域的基础网格单元;所述动网格区域为以待研究物体的质心为球心的球形区域,所述动网格区域的外边界为所述静网格区域的内边界;在所述动网格区域运动的情况下,确定所述动网格区域与所述静网格区域交界面上各基础网格单元的实时位置,并基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数;基于所述交界面的流场参数确定所述交界面对应的通量值,能够兼顾计算效率、计算精度和通用性。
Description
技术领域
本申请涉及数据处理技术领域,尤其涉及一种用于动导数数值模拟的滑移动网格生成方法和装置。
背景技术
动导数是飞行器动态特性研究的重要内容,数值模拟是研究动导数的重要手段,为此,动导数研究工作者对动网格生成方法进行了广泛研究,形成各种动网格生成方法。但现有的动网格生成方法根据所求解问题的不同,具有不同的侧重点,无法同时兼顾计算效率、计算精度和通用性。
发明内容
本申请提供一种用于动导数数值模拟的滑移动网格生成方法和装置,以用于解决现有的动网格生成方法无法同时兼顾计算效率、计算精度和通用性的问题。
本申请提供一种用于动导数数值模拟的滑移动网格生成方法,所述方法包括:
基于预设数值模拟要求将目标计算区域划分为动网格区域和静网格区域,并分别生成所述动网格区域和所述静网格区域的基础网格单元;所述动网格区域为以待研究物体的质心为球心的球形区域,所述动网格区域的外边界为所述静网格区域的内边界;
在所述动网格区域运动的情况下,确定所述动网格区域与所述静网格区域交界面上各基础网格单元的实时位置,并基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数;
基于所述交界面的流场参数确定所述交界面对应的通量值。
根据本申请提供的一种用于动导数数值模拟的滑移动网格生成方法,所述基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数,具体包括:
基于所述交界面上各基础网格单元的实时位置,确定各基础网格单元对应的流场参数平均值;
对各基础网格单元对应的流场参数平均值进行重构,以得到各基础网格单元对应的流场参数节点值。
根据本申请提供的一种用于动导数数值模拟的滑移动网格生成方法,所述基于所述交界面的流场参数确定所述交界面对应的通量值,具体包括:
基于所述各基础网格单元对应的流场参数节点值,利用通量向量分裂方法计算得到各基础网格单元对应的通量值以及所述通量值对应的正特征值和负特征值。
根据本申请提供的一种用于动导数数值模拟的滑移动网格生成方法,所述基于所述交界面上各基础网格单元的实时位置,确定各基础网格单元对应的流场参数平均值,具体包括:
基于所述交界面上各基础网格单元的实时位置,利用单元中心的有限体积方法求解积分形式的流体力学方程组以得到各基础网格单元对应的流场参数平均值。
本申请还提供一种用于动导数数值模拟的滑移动网格生成装置,所述装置包括:
基础网格生成模块,用于基于预设数值模拟要求将目标计算区域划分为动网格区域和静网格区域,并分别生成所述动网格区域和所述静网格区域的基础网格单元;所述动网格区域为以待研究物体的质心为球心的球形区域,所述动网格区域的外边界为所述静网格区域的内边界;
流场参数重构模块,用于在所述动网格区域运动的情况下,确定所述动网格区域与所述静网格区域交界面上各基础网格单元的实时位置,并基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数;
通量确定模块,用于基于所述交界面的流场参数确定所述交界面对应的通量值。
根据本申请提供的一种用于动导数数值模拟的滑移动网格生成装置,所述基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数,具体包括:
基于所述交界面上各基础网格单元的实时位置,确定各基础网格单元对应的流场参数平均值;
对各基础网格单元对应的流场参数平均值进行重构,以得到各基础网格单元对应的流场参数节点值。
根据本申请提供的一种用于动导数数值模拟的滑移动网格生成装置,所述基于所述交界面的流场参数确定所述交界面对应的通量值,具体包括:
基于所述各基础网格单元对应的流场参数节点值,利用通量向量分裂方法计算得到各基础网格单元对应的通量值以及所述通量值对应的正特征值和负特征值。
根据本申请提供的一种用于动导数数值模拟的滑移动网格生成装置,所述基于所述交界面上各基础网格单元的实时位置,确定各基础网格单元对应的流场参数平均值,具体包括:
基于所述交界面上各基础网格单元的实时位置,利用单元中心的有限体积方法求解积分形式的流体力学方程组以得到各基础网格单元对应的流场参数平均值。
本申请实施例还提供了一种电子设备,所述电子设备包括存储器和处理器,所述存储器存储有计算机程序,所述处理器运行所述计算机程序时,执行上述的用于动导数数值模拟的滑移动网格生成方法中的任一实现方式中的步骤。
本申请实施例还提供了一种可读存储介质,所述可读取存储介质中存储有计算机程序,所述计算机程序在处理器上运行时,执行上述的用于动导数数值模拟的滑移动网格生成方法中的任一实现方式中的步骤。
本申请提供的用于动导数数值模拟的滑移动网格生成方法和装置,所述方法包括:基于预设数值模拟要求将目标计算区域划分为动网格区域和静网格区域,并分别生成所述动网格区域和所述静网格区域的基础网格单元;所述动网格区域为以待研究物体的质心为球心的球形区域,所述动网格区域的外边界为所述静网格区域的内边界;在所述动网格区域运动的情况下,确定所述动网格区域与所述静网格区域交界面上各基础网格单元的实时位置,并基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数;基于所述交界面的流场参数确定所述交界面对应的通量值,通过引入球形核动网格,不但满足了通用的俯仰、偏航、滚转动导数数值模拟需求,同时降低了边界定位的计算复杂性,通过球形动、静边界滑移面守恒型通量重构,实现了动导数高精度数值模拟。
附图说明
为了更清楚地说明本申请或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本申请提供的用于动导数数值模拟的滑移动网格生成方法的流程示意图;
图2是本申请提供的动网格区域和静网格区域的位置关系示意图;
图3是本申请提供的动网格区域和静网格区域的立体结构示意图;
图4是本申请提供的用于动导数数值模拟的滑移动网格生成装置的结构示意图;
图5是本申请提供的电子设备的结构示意图。
图标:201-静网格区域;202-动网格区域;203-动、静网格交界面;204-物面边界;205-计算区域外边界;401-处理器;402-通信接口;403-存储器;404-通信总线。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,下面将结合本申请中的附图,对本申请中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
图1是本申请提供的用于动导数数值模拟的滑移动网格生成方法的流程示意图,如图1所示,该方法包括:
步骤101,基于预设数值模拟要求将目标计算区域划分为动网格区域202和静网格区域201,并分别生成所述动网格区域202和所述静网格区域201的基础网格单元;所述动网格区域202为以待研究物体的质心为球心的球形区域,所述动网格区域202的外边界为所述静网格区域201的内边界。
具体的,所述预设数值模拟要求为相关工作人员根据动导数数值模拟需求预先确定并设置的,基于此,本申请实施例即可基于预设数值模拟要求将目标计算区域划分为动网格区域202和静网格区域201,并分别生成所述动网格区域202和所述静网格区域201的基础网格单元。所述基础网格单元的生成步骤包括:(1)、定义基本几何元素,所述基本几何元素包括点、线、面、块;(2)、定义各基本几何元素的对应关系;(3)、基于各基本几何元素的对应关系构造块网格;(4)、分析块网格的完备性和可生成性;(5)、在各块网络的完备性和可生成性均符合要求的情况下生成对应的基础网格单元并记录基础网格单元信息,例如坐标。基于上述方式,最多能够采用六面体作为基础网格单元,保证了基础网格单元的生成效率,且可适应任意复杂外形的目标计算区域。
图2是本申请提供的动网格区域和静网格区域的位置关系示意图,图3是本申请提供的动网格区域和静网格区域的立体结构示意图,如图2-3所示,静网格区域201和动网格区域202位于计算区域外边界205之内;所述动网格区域202为以待研究物体的质心为球心的球形区域,204为物面边界,即待研究物体的外表面。所述动网格区域202的外边界为所述静网格区域201的内边界(即动、静网格交界面203)。由于变形动网格和超限插值动网格的节点和单元随时间变化会导致网格计算复杂性和数值模拟过程的复杂度增加,重叠网格的动、静网格区域重叠需要随时间判别重叠区域各自对方区域的位置,这也导致计算复杂度增加。基于此,本申请实施例利用球形核动网格的对称性,能够最大限度降低计算复杂性。同时,重叠网格边界定位非常复杂,影响计算效率,而且程序实现难度大。滑移动网格的动网格外边界与静网格内边界是同一边界面,此时边界面上的离散点在动网格运动中,动、静边界离散点相对位置发生变化,需要定位,本申请实施例利用球形核动网格的对称性,简单坐标变换即可实现动、静网格边界定位。进一步的,本申请实施例采用滑移动网格的网格单元不变性,避开了几何守恒律问题,进一步降低了计算复杂性。具体的,动网格区域202与静网格区域201交界面作为静网格区域201的边界时,被赋予静态属性,作为动网格区域202的边界时,被赋予运动属性。当动网格区域202整体绕质心俯仰、侧滑、滚转运动时,动网格区域202的基础网格单元运动中保持几何不变,此时避开了几何守恒性问题,被赋予动、静属性的边界面沿球面发生滑移,由于球面的完全对称性,动、静网格区域的交界面保持几何不变,交界面上的离散点、离散面发生滑移。基于此,能够最大限度提高计算效率。
同时,动导数数值模拟的特点是研究刚体绕质心俯仰、侧滑、滚转运动,不考虑物体结构的柔性变形。本申请实施例的动网格区域202和静网格区域201能够适用于绕质心俯仰、侧滑、滚转运动的动导数数值模拟,提高了通用性,而根据运动所涉及的三个自由度,完全对称的球形区域边界也能满足边界面的运动不变性。
步骤102,在所述动网格区域202运动的情况下,确定所述动网格区域202与所述静网格区域201交界面上各基础网格单元的实时位置,并基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数。
具体的,基于前述内容可知,动网格区域202发生运动后,动、静网格区域的交界面的球面会滑移移动,保证了计算区域的几何连续,而位于交界面上的分别为动、静属性的离散点和离散面发生相对位置滑移,通过已知的运动方程即可快速确定所述动网格区域202与所述静网格区域201交界面上各基础网格单元的实时位置,大大简化了动、静网格边界彼此定位的计算复杂性。传统的重叠网格不考虑边界通量的物理意义,对边界上的变量采用数值逼近的插值方法,并由此构造边界通量,这势必会导致非物理解。针对该问题,本申请实施例在确定所述动网格区域202与所述静网格区域201交界面上各基础网格单元的实时位置之后,进一步基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数。所述基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数,具体包括:
基于所述交界面上各基础网格单元的实时位置,确定各基础网格单元对应的流场参数平均值;
对各基础网格单元对应的流场参数平均值进行重构,以得到各基础网格单元对应的流场参数节点值。
更具体的,所述基于所述交界面上各基础网格单元的实时位置,确定各基础网格单元对应的流场参数平均值,具体包括:
基于所述交界面上各基础网格单元的实时位置,利用单元中心的有限体积方法求解积分形式的流体力学方程组以得到各基础网格单元对应的流场参数平均值。
可以理解的是,对各基础网格单元对应的流场参数平均值进行重构可以采用TVD型重构、通量限制器型重构、ENO重构等重构方式,本申请实施例对此不作具体限定。基于此,能够高效确定各基础网格单元对应的流场参数节点值。
步骤103,基于所述交界面的流场参数确定所述交界面对应的通量值。
具体的,所述基于所述交界面的流场参数确定所述交界面对应的通量值,具体包括:
基于所述各基础网格单元对应的流场参数节点值,利用通量向量分裂方法计算得到各基础网格单元对应的通量值以及所述通量值对应的正特征值和负特征值。所述通量向量分裂方法为本领域常用的差分构造方法,其原理及效果在此不作过多赘述。基于上述方法,本申请实施例能够根据流体微元边界上流动变量Riemann间断原理重构流体微元边界面的通量,以得到守恒型通量,提高计算精度。
本申请实施例提供的方法,基于预设数值模拟要求将目标计算区域划分为动网格区域202和静网格区域201,并分别生成所述动网格区域202和所述静网格区域201的基础网格单元;所述动网格区域202为以待研究物体的质心为球心的球形区域,所述动网格区域202的外边界为所述静网格区域201的内边界;在所述动网格区域202运动的情况下,确定所述动网格区域202与所述静网格区域201交界面上各基础网格单元的实时位置,并基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数;基于所述交界面的流场参数确定所述交界面对应的通量值,通过引入球形核动网格,不但满足了通用的俯仰、偏航、滚转动导数数值模拟需求,同时降低了边界定位的计算复杂性,通过球形动、静边界滑移面守恒型通量重构,实现了动导数高精度数值模拟。
下面对本申请提供的用于动导数数值模拟的滑移动网格生成装置进行描述,下文描述的用于动导数数值模拟的滑移动网格生成装置与上文描述的用于动导数数值模拟的滑移动网格生成方法可相互对应参照。
基于上述任一实施例,图4是本申请提供的用于动导数数值模拟的滑移动网格生成装置的结构示意图,如图4所示,该装置包括:
基础网格生成模块301,用于基于预设数值模拟要求将目标计算区域划分为动网格区域202和静网格区域201,并分别生成所述动网格区域202和所述静网格区域201的基础网格单元;所述动网格区域202为以待研究物体的质心为球心的球形区域,所述动网格区域202的外边界为所述静网格区域201的内边界;
流场参数重构模块302,用于在所述动网格区域202运动的情况下,确定所述动网格区域202与所述静网格区域201交界面上各基础网格单元的实时位置,并基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数;
通量确定模块303,用于基于所述交界面的流场参数确定所述交界面对应的通量值。
本申请实施例提供的装置,基础网格生成模块301基于预设数值模拟要求将目标计算区域划分为动网格区域202和静网格区域201,并分别生成所述动网格区域202和所述静网格区域201的基础网格单元;所述动网格区域202为以待研究物体的质心为球心的球形区域,所述动网格区域202的外边界为所述静网格区域201的内边界;流场参数重构模块302在所述动网格区域202运动的情况下,确定所述动网格区域202与所述静网格区域201交界面上各基础网格单元的实时位置,并基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数;通量确定模块303基于所述交界面的流场参数确定所述交界面对应的通量值,通过引入球形核动网格,不但满足了通用的俯仰、偏航、滚转动导数数值模拟需求,同时降低了边界定位的计算复杂性,通过球形动、静边界滑移面守恒型通量重构,实现了动导数高精度数值模拟。
基于上述实施例,所述基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数,具体包括:
基于所述交界面上各基础网格单元的实时位置,确定各基础网格单元对应的流场参数平均值;
对各基础网格单元对应的流场参数平均值进行重构,以得到各基础网格单元对应的流场参数节点值。
基于上述任一实施例,所述基于所述交界面的流场参数确定所述交界面对应的通量值,具体包括:
基于所述各基础网格单元对应的流场参数节点值,利用通量向量分裂方法计算得到各基础网格单元对应的通量值以及所述通量值对应的正特征值和负特征值。
基于上述任一实施例,所述基于所述交界面上各基础网格单元的实时位置,确定各基础网格单元对应的流场参数平均值,具体包括:
基于所述交界面上各基础网格单元的实时位置,利用单元中心的有限体积方法求解积分形式的流体力学方程组以得到各基础网格单元对应的流场参数平均值。
图5示例了一种电子设备的实体结构示意图,如图5所示,该电子设备可以包括:处理器401、通信接口402、存储器403和通信总线404,其中,处理器401,通信接口402,存储器403通过通信总线404完成相互间的通信。处理器401可以调用存储器403中的逻辑指令,以执行上述各方法所提供的用于动导数数值模拟的滑移动网格生成方法,所述方法包括:基于预设数值模拟要求将目标计算区域划分为动网格区域202和静网格区域201,并分别生成所述动网格区域202和所述静网格区域201的基础网格单元;所述动网格区域202为以待研究物体的质心为球心的球形区域,所述动网格区域202的外边界为所述静网格区域201的内边界;在所述动网格区域202运动的情况下,确定所述动网格区域202与所述静网格区域201交界面上各基础网格单元的实时位置,并基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数;基于所述交界面的流场参数确定所述交界面对应的通量值。
此外,上述的存储器403中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面,本申请还提供一种计算机程序产品,所述计算机程序产品包括计算机程序,计算机程序可存储在非暂态计算机可读存储介质上,所述计算机程序被处理器执行时,计算机能够执行上述各方法所提供的用于动导数数值模拟的滑移动网格生成方法,所述方法包括:基于预设数值模拟要求将目标计算区域划分为动网格区域202和静网格区域201,并分别生成所述动网格区域202和所述静网格区域201的基础网格单元;所述动网格区域202为以待研究物体的质心为球心的球形区域,所述动网格区域202的外边界为所述静网格区域201的内边界;在所述动网格区域202运动的情况下,确定所述动网格区域202与所述静网格区域201交界面上各基础网格单元的实时位置,并基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数;基于所述交界面的流场参数确定所述交界面对应的通量值。
又一方面,本申请还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现以执行上述各方法提供的用于动导数数值模拟的滑移动网格生成方法,所述方法包括:基于预设数值模拟要求将目标计算区域划分为动网格区域202和静网格区域201,并分别生成所述动网格区域202和所述静网格区域201的基础网格单元;所述动网格区域202为以待研究物体的质心为球心的球形区域,所述动网格区域202的外边界为所述静网格区域201的内边界;在所述动网格区域202运动的情况下,确定所述动网格区域202与所述静网格区域201交界面上各基础网格单元的实时位置,并基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数;基于所述交界面的流场参数确定所述交界面对应的通量值。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。
Claims (10)
1.一种用于动导数数值模拟的滑移动网格生成方法,其特征在于,所述方法包括:
基于预设数值模拟要求将目标计算区域划分为动网格区域(202)和静网格区域(201),并分别生成所述动网格区域(202)和所述静网格区域(201)的基础网格单元;所述动网格区域(202)为以待研究物体的质心为球心的球形区域,所述动网格区域(202)的外边界为所述静网格区域(201)的内边界;
在所述动网格区域(202)运动的情况下,确定所述动网格区域(202)与所述静网格区域(201)交界面上各基础网格单元的实时位置,并基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数;
基于所述交界面的流场参数确定所述交界面对应的通量值。
2.根据权利要求1所述的用于动导数数值模拟的滑移动网格生成方法,其特征在于,所述基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数,具体包括:
基于所述交界面上各基础网格单元的实时位置,确定各基础网格单元对应的流场参数平均值;
对各基础网格单元对应的流场参数平均值进行重构,以得到各基础网格单元对应的流场参数节点值。
3.根据权利要求2所述的用于动导数数值模拟的滑移动网格生成方法,其特征在于,所述基于所述交界面的流场参数确定所述交界面对应的通量值,具体包括:
基于所述各基础网格单元对应的流场参数节点值,利用通量向量分裂方法计算得到各基础网格单元对应的通量值以及所述通量值对应的正特征值和负特征值。
4.根据权利要求3所述的用于动导数数值模拟的滑移动网格生成方法,其特征在于,所述基于所述交界面上各基础网格单元的实时位置,确定各基础网格单元对应的流场参数平均值,具体包括:
基于所述交界面上各基础网格单元的实时位置,利用单元中心的有限体积方法求解积分形式的流体力学方程组以得到各基础网格单元对应的流场参数平均值。
5.一种用于动导数数值模拟的滑移动网格生成装置,其特征在于,所述装置包括:
基础网格生成模块,用于基于预设数值模拟要求将目标计算区域划分为动网格区域(202)和静网格区域(201),并分别生成所述动网格区域(202)和所述静网格区域(201)的基础网格单元;所述动网格区域(202)为以待研究物体的质心为球心的球形区域,所述动网格区域(202)的外边界为所述静网格区域(201)的内边界;
流场参数重构模块,用于在所述动网格区域(202)运动的情况下,确定所述动网格区域(202)与所述静网格区域(201)交界面上各基础网格单元的实时位置,并基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数;
通量确定模块,用于基于所述交界面的流场参数确定所述交界面对应的通量值。
6.根据权利要求5所述的用于动导数数值模拟的滑移动网格生成装置,其特征在于,所述基于所述交界面上各基础网格单元的实时位置重构所述交界面的流场参数,具体包括:
基于所述交界面上各基础网格单元的实时位置,确定各基础网格单元对应的流场参数平均值;
对各基础网格单元对应的流场参数平均值进行重构,以得到各基础网格单元对应的流场参数节点值。
7.根据权利要求6所述的用于动导数数值模拟的滑移动网格生成装置,其特征在于,所述基于所述交界面的流场参数确定所述交界面对应的通量值,具体包括:
基于所述各基础网格单元对应的流场参数节点值,利用通量向量分裂方法计算得到各基础网格单元对应的通量值以及所述通量值对应的正特征值和负特征值。
8.根据权利要求7所述的用于动导数数值模拟的滑移动网格生成装置,其特征在于,所述基于所述交界面上各基础网格单元的实时位置,确定各基础网格单元对应的流场参数平均值,具体包括:
基于所述交界面上各基础网格单元的实时位置,利用单元中心的有限体积方法求解积分形式的流体力学方程组以得到各基础网格单元对应的流场参数平均值。
9.一种电子设备,其特征在于,所述电子设备包括存储器和处理器,所述存储器存储有计算机程序,所述处理器运行所述计算机程序时执行权利要求1至4中任一项所述的用于动导数数值模拟的滑移动网格生成方法。
10.一种可读存储介质,其特征在于,所述可读存储介质中存储有计算机程序,所述计算机程序在处理器上运行时执行权利要求1至4中任一项所述的用于动导数数值模拟的滑移动网格生成方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310160122.6A CN115841063B (zh) | 2023-02-24 | 2023-02-24 | 用于动导数数值模拟的滑移动网格生成方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310160122.6A CN115841063B (zh) | 2023-02-24 | 2023-02-24 | 用于动导数数值模拟的滑移动网格生成方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115841063A true CN115841063A (zh) | 2023-03-24 |
CN115841063B CN115841063B (zh) | 2023-04-28 |
Family
ID=85580146
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310160122.6A Active CN115841063B (zh) | 2023-02-24 | 2023-02-24 | 用于动导数数值模拟的滑移动网格生成方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115841063B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116522827A (zh) * | 2023-07-04 | 2023-08-01 | 北京凌云智擎软件有限公司 | 一种网格单元边界面的流动变量重构方法、设备及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000331896A (ja) * | 1999-05-25 | 2000-11-30 | Nec Corp | メッシュ生成装置およびメッシュ生成方法 |
CN101599104A (zh) * | 2009-07-16 | 2009-12-09 | 北京航空航天大学 | 一种航空涡轮发动机叶片颤振边界的模拟方法 |
EP3467836A2 (en) * | 2017-10-06 | 2019-04-10 | Sumitomo Heavy Industries, Ltd. | Particle simulation method, simulation apparatus, and program |
DE102017130593A1 (de) * | 2017-12-19 | 2019-06-19 | Airbus Defence and Space GmbH | Quantitative Blockierungsanalyse eines Gitters in Vereisungsbedingungen |
CN115270650A (zh) * | 2022-06-17 | 2022-11-01 | 大连海洋大学 | 一种基于cfd模拟的养殖池中鱼群运动可视化数值模拟方法 |
CN115272187A (zh) * | 2022-06-27 | 2022-11-01 | 江苏大学 | 一种基于车载动态视场帧间关联性的田间稻麦倒伏全局评估方法 |
CN115329696A (zh) * | 2022-10-13 | 2022-11-11 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于非贴体网格的守恒型固壁边界数值模拟方法及设备 |
-
2023
- 2023-02-24 CN CN202310160122.6A patent/CN115841063B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000331896A (ja) * | 1999-05-25 | 2000-11-30 | Nec Corp | メッシュ生成装置およびメッシュ生成方法 |
CN101599104A (zh) * | 2009-07-16 | 2009-12-09 | 北京航空航天大学 | 一种航空涡轮发动机叶片颤振边界的模拟方法 |
EP3467836A2 (en) * | 2017-10-06 | 2019-04-10 | Sumitomo Heavy Industries, Ltd. | Particle simulation method, simulation apparatus, and program |
DE102017130593A1 (de) * | 2017-12-19 | 2019-06-19 | Airbus Defence and Space GmbH | Quantitative Blockierungsanalyse eines Gitters in Vereisungsbedingungen |
CN115270650A (zh) * | 2022-06-17 | 2022-11-01 | 大连海洋大学 | 一种基于cfd模拟的养殖池中鱼群运动可视化数值模拟方法 |
CN115272187A (zh) * | 2022-06-27 | 2022-11-01 | 江苏大学 | 一种基于车载动态视场帧间关联性的田间稻麦倒伏全局评估方法 |
CN115329696A (zh) * | 2022-10-13 | 2022-11-11 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于非贴体网格的守恒型固壁边界数值模拟方法及设备 |
Non-Patent Citations (5)
Title |
---|
LI, ZOUSHUANG 等: "Computational fluid dynamics modeling of spatial atomic layer deposition on microgroove substrates", 《INTERNATIONAL JOURNAL OF HEAT & MASS TRANSFER》 * |
周星;许厚谦;: "计算含大位移动边界非定常流动的无网格算法", 《力学与实践》 * |
唐静;郑鸣;邓有奇;李彬;: "网格自适应技术在复杂外形流场模拟中的应用", 《计算力学学报》 * |
廖佳文 等: "基于二重网格的动网格松弛法改进", 《计算力学学报》 * |
郭泽宇;陈作钢;: "基于复杂柔性运动流场数值求解的动网格方案及应用", 《上海交通大学学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116522827A (zh) * | 2023-07-04 | 2023-08-01 | 北京凌云智擎软件有限公司 | 一种网格单元边界面的流动变量重构方法、设备及装置 |
CN116522827B (zh) * | 2023-07-04 | 2023-10-20 | 北京凌云智擎软件有限公司 | 一种气动热环境计算的流动变量重构方法、设备及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN115841063B (zh) | 2023-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108073682B (zh) | 基于参数视图函数查询数据库 | |
Holdstein et al. | Three-dimensional surface reconstruction using meshing growing neural gas (MGNG) | |
CN115841063A (zh) | 用于动导数数值模拟的滑移动网格生成方法和装置 | |
US20150127301A1 (en) | Updating A CAD Model To Reflect Global Or Local Shape Changes | |
CN107886569B (zh) | 一种基于离散李导数的测度可控的曲面参数化方法及系统 | |
CN116310048B (zh) | 一种基于曲率细分计算光线跟踪与nurbs曲面交点的方法 | |
Hansen et al. | Fixed-axis tool positioning with built-in global interference checking for NC path generation | |
CN110008387B (zh) | 流场可视化实现方法、装置及电子设备 | |
Gonçalves et al. | A benchmark study on accuracy-controlled distance calculation between superellipsoid and superovoid contact geometries | |
CN110188395B (zh) | 一种基于线面体的维度增加式计算流体网格生成方法 | |
Mörschel et al. | Generating neuron geometries for detailed three-dimensional simulations using anamorph | |
Ramanathan et al. | Interior Medial Axis Transform computation of 3D objects bound by free-form surfaces | |
Rego et al. | Growing self-reconstruction maps | |
Delingette | Adaptive and deformable models based on simplex meshes | |
CN106649992A (zh) | 舰船与尾迹的网格模型的融合与优化方法 | |
Jiang et al. | Reparameterization of B-spline surface and its application in ship hull modeling | |
Makhlouf et al. | Approach for CAD model Reconstruction from a deformed mesh | |
Jenkins et al. | A unit quaternion based som for anatomical joint constraint modelling | |
Tang et al. | A local and fast interpolation method for mesh deformation | |
CN113269890B (zh) | 边界面网格生成方法、智能终端以及存储装置 | |
CN117217130B (zh) | 一种逐层推进的壁面距离确定方法、装置以及介质 | |
Jenkins et al. | Self-Organizing Maps for Anatomical Joint Constraint Modelling | |
Shih et al. | Geometry and mesh generation for high fidelity computational simulations using non-uniform rational B-splines | |
Gandotra et al. | Representation of model for efficient collision detection in virtual reality environment | |
Lemasson et al. | Accurate thickness computation of a B-Rep model on the GPU |
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 |