CN114841095A - 一种不可压缩流动的扰动域推进方法 - Google Patents

一种不可压缩流动的扰动域推进方法 Download PDF

Info

Publication number
CN114841095A
CN114841095A CN202210781128.0A CN202210781128A CN114841095A CN 114841095 A CN114841095 A CN 114841095A CN 202210781128 A CN202210781128 A CN 202210781128A CN 114841095 A CN114841095 A CN 114841095A
Authority
CN
China
Prior art keywords
grid
domain
pressure
disturbance
flow field
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
Application number
CN202210781128.0A
Other languages
English (en)
Other versions
CN114841095B (zh
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.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN202210781128.0A priority Critical patent/CN114841095B/zh
Publication of CN114841095A publication Critical patent/CN114841095A/zh
Application granted granted Critical
Publication of CN114841095B publication Critical patent/CN114841095B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12Simultaneous equations, e.g. systems of linear equations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Data Mining & Analysis (AREA)
  • Fluid Mechanics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及计算流体力学技术领域,提供了一种不可压缩流动的扰动域推进方法,该方法包括:基于来流速度,建立初始不可压缩流动的计算域和初始流场;基于初始流场,通过SIMPLE方法,进行第一步迭代,获得更新后的速度流场和压力流场;将第一步迭代计算出的当地网格残差作为数值扰动,通过网格单元的位置和距离物体壁面的距离判断,进行扰动域的推进;在扰动域推进后进行下一步迭代循环,获得收敛的解。本发明通过建立随收敛进度变化的动态扰动域,减少了单次迭代过程中的无效计算,避免了大量计算资源的浪费,提高了不可压缩流动数值模拟求解的计算效率。

Description

一种不可压缩流动的扰动域推进方法
技术领域
本发明涉及计算流体力学技术领域,尤其涉及一种不可压缩流动的扰动域推进方法。
背景技术
计算流体力学是计算科学的一个分支,专注于求解流体力学控制方程。科学家和工程师经常通过计算流体力学方法得到流体运动的定量信息,比如在特定条件下流场内不同位置的速度、压力、密度或者温度等参数的值。当速度比较低时,上述流体力学控制方程为不可压缩的纳维-斯托克斯方程,需要求解的物理参数主要包括速度和压力,而求解该方程的主流算法是以SIMPLE(Semi-Implicit Method for Pressure-Linked Equations)系列算法为代表的压力与速度耦合方法。
尽管计算机硬件性能发展迅速,但上述对于航空航天等领域真实飞行器在飞行包线内的性能计算方法,仍需要耗费大量的计算资源和超长周期,难以满足实际工程应用的需求。为了提高数值计算的效率,研究者开展了代数多重网格、并行计算等多方面的研究,在一定程度上提高了计算效率。但随着计算流体力学应用范围越来越广,应用层次逐渐加深,在工程实践中对计算流体力学方法的效率仍有着更高的要求。
扰动域推进方法提出了一种加速求解的新思路,针对可压缩流动的数值模拟求解,传统方法预设固定计算域存在大量无效计算,而扰动域推进方法着眼于计算域在时间推进求解过程中的自适应变化,以期减少无效计算量,提高单个迭代步的收敛速率。然而,这一方法针对的是可压缩流动,不适用于不可压缩纳维-斯托克斯方程的求解。其原因在于对于不可压缩流场,声速趋于无穷大,流场内任一点的扰动可以传遍整个流场,原扰动域推进方法根据声速和当地速度作为扰动传播的判断依据对不可压缩流动而言是失效的。因此,需要针对不可压缩流动的方程性质和求解方法进行深入研究,提出不可压缩流动的扰动域推进方法。
发明内容
有鉴于此,本发明提供了一种不可压缩流动的扰动域推进方法,以解决现有技术中无法提高单次迭代过程中的无效计算以及不可压缩流动数值模拟求解的计算效率的问题。
本发明提供了一种不可压缩流动的扰动域推进方法,包括:
S1基于来流速度,建立初始不可压缩流动的扰动域和初始流场;
S2基于所述初始流场,通过SIMPLE方法,进行第一步迭代,获得更新后的速度流场和压力流场;
S3将由更新后的速度流场和压力流场计算出的当地网格残差作为数值扰动,通过网格单元的位置和距离物体壁面的距离判断,进行扰动域推进;
S4在所述扰动域推进后进行下一步迭代循环,获得收敛的解。
进一步地,所述S1中所述不可压缩流动,包括二维不可压缩流动和三维不可压缩流动。
进一步地,所述不可压缩流动为二维不可压缩流动时,所述S1包括:
根据来流速度,建立20倍飞行器长度的初始网格,并将全部所述初始网格作为初始扰动域;
将不可压缩纳维-斯托克斯方程组在每一网格上解耦分解和离散表示为动量离散方程和压力修正方程;
通过假定一个速度分布和一个压力场,计算出所述动量离散方程中的系数及常数项。
进一步地,所述不可压缩纳维-斯托克斯方程组、所述动量离散方程和所述压力修正方程,表达式分别如下:
Figure 4067DEST_PATH_IMAGE001
Figure 444276DEST_PATH_IMAGE002
Figure 45022DEST_PATH_IMAGE003
其中,xy为二维的两个方向,u、vxy两个方向上的速度,t为时间,P为压力,ρ为密度,μ为动力粘性系数,a、b、a'、b'均为不可压缩纳维-斯托克斯方程组在当前时刻离散得到的方程系数,其中,a包括a w 、a N 、a E 、a S 、a P b包括b P a'包括a' P 、a' w 、a' N 、a' E 、a' S b'包括b' P ;ϕ表示由速度uv构成的矩阵,包括ϕ P 、ϕ W 、ϕ E 、ϕ S 、ϕ N p'为修正压力,包括p' P 、p' E p' W 、p' N 、p' S ;下标P指当前网格,W、E、S、N指与当前网格相邻的左、右、下、上四个网格。
进一步地,所述S2,包括:
依次在二维方向进行迭代求解动量离散方程,获得二维方向上的新的速度分布;
求解所述压力修正方程,获得修正压力;
通过所述修正压力,获得修正速度;
利用所述修正压力和修正速度,获得更新后的所述压力流场和所述速度流场。
进一步地,所述S3,包括:
根据计算条件和计算需求,设置全局收敛参数;
将所述更新后的压力流场和速度流场,通过不可压缩纳维-斯托克斯方程组,计算求得所述当地网格的残差;
判断所述当地网格的残差对当地网格单元格心与物体壁面的距离的影响是否小于预先设定的删除阈值,若不小于,则不作为;若小于,则进一步判断所述当地网格单元格心与物体壁面的距离是否位于所述删除阈值的扰动域的最外围,若是,删除所有处于所述扰动域中的与所述网格单元紧邻的网格单元;若否,不删除;
循环判断整个扰动域,完成所述扰动域的推进。
进一步地,所述当地网格的残差对当地网格单元格心与物体壁面的距离的影响与所述预先设定的删除阈值的大小关系判断表达式如下:
Figure 305102DEST_PATH_IMAGE004
其中,ɛ loc 为当地网格的残差,ε d 为删除阈值,d loc 为网格单元格心距物体壁面的距离,l为所计算问题的特征尺度。
进一步地,与所述网格单元紧邻的网格单元距离壁面中心的距离,满足小于等于所述网格单元距离壁面中心的距离,表达式如下:
Figure 665807DEST_PATH_IMAGE005
其中,d adj 为紧邻的网格单元距离壁面中心的距离。
进一步地,所述S4,包括:
统计所有所述当地网格残差,并取所述当地网格残差的最大值;
判断所述当地网格残差的最大值是否小于所述全局收敛参数,若是,则输出结果完成计算;若否,跳转至S2。
进一步地,所述收敛的解,包括:所述更新后的速度流场和压力流场。
本发明与现有技术相比存在的有益效果是:
本发明通过考虑不可压缩流动扰动域中残差数值扰动与物面距离成反比的特点,通过建立随收敛进度变化的动态扰动域,有效减少单次迭代过程中的无效计算,有效提高不可压缩流动数值模拟求解的计算效率。
附图说明
为了更清楚地说明本发明中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1是本发明提供的一种不可压缩流动的扰动域推进方法的流程图;
图2是本发明提供的建立初始不可压缩流动的扰动域和初始流场的流程图;
图3是本发明提供的获得更新后的速度流场和压力流场的流程图;
图4是本发明提供的进行扰动域的推进的流程图。
具体实施方式
以下描述中,为了说明而不是为了限定,提出了诸如特定系统结构、技术之类的具体细节,以便透彻理解本发明实施例。然而,本领域的技术人员应当清楚,在没有这些具体细节的其它实施例中也可以实现本发明。在其它情况中,省略对众所周知的系统、装置、电路以及方法的详细说明,以免不必要的细节妨碍本发明的描述。
下面将结合附图详细说明根据本发明的一种不可压缩流动的扰动域推进方法。
图1是本发明提供的一种不可压缩流动的扰动域推进方法的流程图。
如图1所示,该扰动域推进方法包括:
S1,基于来流速度,建立初始不可压缩流动的计算域和初始流场。
S1中不可压缩流动,包括二维不可压缩流动和三维不可压缩流动。
图2是本发明提供的建立初始不可压缩流动的扰动域和初始流场的流程图。
如图2所示,不可压缩流动为二维不可压缩流动时,S1包括:
S11,根据来流速度,建立20倍飞行器长度的初始网格,并将全部初始网格作为初始扰动域。
参考具体计算算例的物理参数,建立20倍特征尺度的初始网格,并将全部网格作为初始扰动域,其中,物理参数指流场初始值来流速度。
S12,将不可压缩纳维-斯托克斯方程组在每一网格上解耦分解和离散表示为动量离散方程和压力修正方程;
S13,通过假定一个速度分布和一个压力场,计算出动量离散方程中的系数及常数项。
假定一个速度分布,分别记为u 0 ,v 0 ,假定一个压力流场P *;以计算动量离散方程中的系数及常数项。
不可压缩纳维-斯托克斯方程组、动量离散方程和压力修正方程,表达式分别如下:
Figure 694943DEST_PATH_IMAGE006
(1)
Figure 997748DEST_PATH_IMAGE007
(2)
Figure 213966DEST_PATH_IMAGE008
(3)
其中,xy为二维的两个方向,u、vxy两个方向上的速度,t为时间,P为压力,ρ为密度,μ为动力粘性系数,a、b、a'、b'均为不可压缩纳维-斯托克斯方程组在当前时刻离散得到的方程系数,其中,a包括a w 、a N 、a E 、a S 、a P b包括b P a'包括a' P 、a' w 、a' N 、a' E 、a' S b'包括b' P ;ϕ表示由速度uv构成的矩阵,包括ϕ P 、ϕ W 、ϕ E 、ϕ S 、ϕ N p'为修正压力,包括p' P 、p' E p' W 、p' N 、p' S ;下标P指当前网格,W、E、S、N指与当前网格相邻的左、右、下、上四个网格。
S2, 基于初始流场,通过SIMPLE方法,进行第一步迭代,获得更新后的速度流场和压力流场。
图3是本发明提供的获得更新后的速度流场和压力流场的流程图。
如图3所示,S2包括:
S21,依次在二维方向进行迭代求解动量离散方程,获得二维方向上的新的速度分布。
依次在xy方向上进行迭代求解动量离散方程式(2),获得ϕ=[u*,v*]。
其中,u*v*xy方向上的新的速度。
S22,求解压力修正方程,获得修正压力。
求解压力修正方程式(3),得修正压力
Figure 627630DEST_PATH_IMAGE009
S23,通过修正压力
Figure 760540DEST_PATH_IMAGE010
,获得修正速度。
按照式(4),由修正压力
Figure 234246DEST_PATH_IMAGE010
,计算得到修正速度
Figure 937760DEST_PATH_IMAGE011
Figure 889536DEST_PATH_IMAGE012
(4)
其中,u P v P 分别为当前网格在x和y方向上的速度修正,
Figure 643996DEST_PATH_IMAGE013
为动量离散方程离散在当前网格所得到的系数,
Figure 23025DEST_PATH_IMAGE014
分别代指左、右、下、上四个网格上的压力修正值。
S24,利用修正压力和修正速度,获得更新后的压力流场和速度流场。
利用修正压力和修正速度,可以得到改进后的压力场和速度场,即
Figure 479414DEST_PATH_IMAGE015
Figure 234880DEST_PATH_IMAGE016
,其中,p*表示新的压力。
S3,将由更新后的速度流场和压力流场计算出的当地网格残差作为数值扰动,通过网格单元的位置和距离物体壁面的距离判断,进行扰动域的推进。
图4是本发明提供的进行扰动域的推进的流程图。
如图4所示,S3包括:
S31,根据计算条件和计算需求,设置全局收敛参数。
在本发明中全局收敛参数设置为10-8
S32,将更新后的压力流场和速度流场,通过不可压缩纳维-斯托克斯方程组,计算求得当地网格的残差。
将S2计算所得更新后的压力场和速度场带入不可压缩纳维-斯托克斯方程组(1)中,计算得到不满足该方程的部分,即为当地网格的残差ɛ loc
S33,判断当地网格的残差对当地网格单元格心与物体壁面的距离的影响是否小于预先设定的删除阈值,若不小于,则不作为;若小于,则进一步判断当地网格单元格心与物体壁面的距离是否位于删除阈值的扰动域的最外围,若是,删除所有处于扰动域中的与网格单元紧邻的网格单元;若否,不删除。
综合考虑当地网格与物体壁面的距离,判断当地网格的残差对与当地网格距离物体壁面距离的影响是否小于可以删除的阈值,若小于,则判断当地网格单元格心距物体壁面的距离是否位于删除阈值的扰动域的最外围,若是,删除所有处于扰动域中的与网格单元紧邻的网格单元;若否,不删除所有处于扰动域中的与网格单元紧邻的网格单元;若不小于,则无需判断当地网格单元格心距物体壁面的距离是否位于删除阈值的扰动域的最外围。
当地网格的残差对当地网格单元格心与物体壁面的距离的影响与预先设定的删除阈值的大小关系判断表达式如下:
Figure 827536DEST_PATH_IMAGE017
其中,ɛ loc 为当地网格的残差,ε d 为删除阈值,d loc 为网格单元格心距物体壁面的距离,l为所计算问题的特征尺度。
S34,循环判断整个扰动域,完成扰动域的推进。
与网格单元紧邻的网格单元距离壁面中心的距离,满足小于等于与网格单元距离壁面中心的距离,表达式如下:
Figure 892313DEST_PATH_IMAGE018
其中,d adj 为紧邻的网格单元距离壁面中心的距离。
S4,在扰动域推进后进行下一步迭代循环,获得收敛的解。
S4,包括:
S41,统计所有当地网格残差,并取当地网格残差的最大值。
统计S3中得到的当地网格残差ɛ loc ,取所有网格残差中的最大值;
S42,判断当地网格残差的最大值是否小于全局收敛参数,若是,则输出结果完成计算;若否,跳转至S2。
本发明通过建立随收敛进度变化的动态扰动域,改善了传统方法对不可压缩流动的数值模拟求解全都是基于固定区域的计算域的固有模式;减少了单次迭代过程中的无效计算,避免了大量计算资源的浪费,提高了不可压缩流动数值模拟求解的计算效率。
上述所有可选技术方案,可以采用任意结合形成本申请的可选实施例,在此不再一一赘述。
应理解,上述实施例中各步骤的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本发明实施例的实施过程构成任何限定。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围,均应包含在本发明的保护范围之内。

Claims (10)

1.一种不可压缩流动的扰动域推进方法,其特征在于,包括:
S1基于来流速度,建立初始不可压缩流动的扰动域和初始流场;
S2基于所述初始流场,通过SIMPLE方法,进行第一步迭代,获得更新后的速度流场和压力流场;
S3将由更新后的速度流场和压力流场计算出的当地网格残差作为数值扰动,通过网格单元的位置和距离物体壁面的距离判断,进行扰动域的推进;
S4在所述扰动域推进后进行下一步迭代循环,获得收敛的解。
2.根据权利要求1所述的扰动域推进方法,其特征在于,所述S1中所述不可压缩流动,包括二维不可压缩流动和三维不可压缩流动。
3.根据权利要求2所述的扰动域推进方法,其特征在于,所述不可压缩流动为二维不可压缩流动时,所述S1包括:
根据来流速度,建立20倍飞行器长度的初始网格,并将全部所述初始网格作为初始扰动域;
将不可压缩纳维-斯托克斯方程组在每一网格上解耦分解和离散表示为动量离散方程和压力修正方程;
通过假定一个速度分布和一个压力场,计算出所述动量离散方程中的系数及常数项。
4.根据权利要求3所述的扰动域推进方法,其特征在于,所述不可压缩纳维-斯托克斯方程组、所述动量离散方程和所述压力修正方程,表达式分别如下:
Figure 50220DEST_PATH_IMAGE001
Figure 334702DEST_PATH_IMAGE002
Figure 389246DEST_PATH_IMAGE003
其中,xy为二维的两个方向,u、vxy两个方向上的速度,t为时间,P为压力,ρ为密度,μ为动力粘性系数,a、b、a'、b'均为不可压缩纳维-斯托克斯方程组在当前时刻离散得到的方程系数,其中,a包括a w 、a N 、a E 、a S 、a P b包括b P a'包括a' P 、a' w 、a' N 、a' E 、a' S b'包括b' P ϕ表示由速度uv构成的矩阵,包括ϕ P 、ϕ W 、ϕ E 、ϕ S 、ϕ N p'为修正压力,包括p' P 、p' E 、p' W p' N 、p' S ;下标P指当前网格,W、E、S、N指与当前网格相邻的左、右、下、上四个网格。
5.根据权利要求4所述的扰动域推进方法,其特征在于,所述S2,包括:
依次在二维方向进行迭代求解动量离散方程,获得二维方向上的新的速度分布;
求解所述压力修正方程,获得修正压力;
通过所述修正压力,获得修正速度;
利用所述修正压力和修正速度,获得更新后的所述压力流场和所述速度流场。
6.根据权利要求5所述的扰动域推进方法,其特征在于,所述S3,包括:
根据计算条件和计算需求,设置全局收敛参数;
将所述更新后的压力流场和速度流场,通过不可压缩纳维-斯托克斯方程组,计算求得所述当地网格的残差;
判断所述当地网格的残差对当地网格单元格心与物体壁面的距离的影响是否小于预先设定的删除阈值,若不小于,则不作为;若小于,则进一步判断所述当地网格单元格心与物体壁面的距离是否位于所述删除阈值的扰动域的最外围,若是,删除所有处于所述扰动域中的与所述网格单元紧邻的网格单元;若否,不删除;
循环判断整个扰动域,完成所述扰动域的推进。
7.根据权利要求6所述的扰动域推进方法,其特征在于,所述当地网格的残差对当地网格单元格心与物体壁面的距离的影响与所述预先设定的删除阈值的大小关系判断表达式如下:
Figure 350249DEST_PATH_IMAGE004
其中,ɛ loc 为当地网格的残差,ε d 为删除阈值,d loc 为网格单元格心距物体壁面的距离,l为所计算问题的特征尺度。
8.根据权利要求7所述的扰动域推进方法,其特征在于,与所述网格单元紧邻的网格单元距离壁面中心的距离,满足小于等于所述网格单元距离壁面中心的距离,表达式如下:
Figure 123033DEST_PATH_IMAGE005
其中,d adj 为紧邻的网格单元距离壁面中心的距离。
9.根据权利要求8所述的扰动域推进方法,其特征在于,所述S4,包括:
统计所有所述当地网格残差,并取所述当地网格残差的最大值;
判断所述当地网格残差的最大值是否小于所述全局收敛参数,若是,则输出结果完成计算;若否,跳转至S2。
10.根据权利要求1所述的扰动域推进方法,其特征在于,所述收敛的解,包括:所述更新后的速度流场和压力流场。
CN202210781128.0A 2022-07-05 2022-07-05 一种不可压缩流动的扰动域推进方法 Active CN114841095B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210781128.0A CN114841095B (zh) 2022-07-05 2022-07-05 一种不可压缩流动的扰动域推进方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210781128.0A CN114841095B (zh) 2022-07-05 2022-07-05 一种不可压缩流动的扰动域推进方法

Publications (2)

Publication Number Publication Date
CN114841095A true CN114841095A (zh) 2022-08-02
CN114841095B CN114841095B (zh) 2022-09-09

Family

ID=82574425

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210781128.0A Active CN114841095B (zh) 2022-07-05 2022-07-05 一种不可压缩流动的扰动域推进方法

Country Status (1)

Country Link
CN (1) CN114841095B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114996634A (zh) * 2022-08-04 2022-09-02 中国空气动力研究与发展中心计算空气动力研究所 一种压力参考点的隐式处理方法、设备及介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108416075A (zh) * 2017-08-29 2018-08-17 沈阳工业大学 基于cfd技术的风力机气动计算模型建模方法
CN108563843A (zh) * 2018-03-26 2018-09-21 北京航空航天大学 定常可压缩流动的扰动区域更新方法
CN110516367A (zh) * 2019-08-28 2019-11-29 西安逆流环保科技有限公司 一种计算较小尺度内人工气流扰动逆温层范围和强度的方法
CN111651831A (zh) * 2020-04-13 2020-09-11 北京航空航天大学 飞行器定常粘性可压缩绕流的分区扰动域更新计算方法
CN113609598A (zh) * 2021-10-09 2021-11-05 北京航空航天大学 飞行器气动特性模拟的rans/les扰动域更新方法
CN113850008A (zh) * 2021-12-02 2021-12-28 北京航空航天大学 飞行器气动特性预测的自适应网格扰动域更新加速方法
CN113962030A (zh) * 2021-12-20 2022-01-21 北京航空航天大学 飞行器多体分离模拟的重叠网格扰动域更新方法
CN114218878A (zh) * 2022-02-17 2022-03-22 北京航空航天大学 一种飞行器机动过程模拟的动网格扰动域更新方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108416075A (zh) * 2017-08-29 2018-08-17 沈阳工业大学 基于cfd技术的风力机气动计算模型建模方法
CN108563843A (zh) * 2018-03-26 2018-09-21 北京航空航天大学 定常可压缩流动的扰动区域更新方法
CN110516367A (zh) * 2019-08-28 2019-11-29 西安逆流环保科技有限公司 一种计算较小尺度内人工气流扰动逆温层范围和强度的方法
CN111651831A (zh) * 2020-04-13 2020-09-11 北京航空航天大学 飞行器定常粘性可压缩绕流的分区扰动域更新计算方法
CN113609598A (zh) * 2021-10-09 2021-11-05 北京航空航天大学 飞行器气动特性模拟的rans/les扰动域更新方法
CN113850008A (zh) * 2021-12-02 2021-12-28 北京航空航天大学 飞行器气动特性预测的自适应网格扰动域更新加速方法
CN113962030A (zh) * 2021-12-20 2022-01-21 北京航空航天大学 飞行器多体分离模拟的重叠网格扰动域更新方法
CN114218878A (zh) * 2022-02-17 2022-03-22 北京航空航天大学 一种飞行器机动过程模拟的动网格扰动域更新方法

Non-Patent Citations (10)

* Cited by examiner, † Cited by third party
Title
CHRISTIAN H.WHITING等: "A stabilized finite element method for the incompressible Nativer-Stokes equations using a hierarchical basis", 《INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS》 *
FB PIDDUCK: "On the Propagation of a Disturbance a Fluid under Gravity", 《ISSN》 *
KARIM DEBBAGH等: "A two-dimensional disturbed flows over a flat plate theoretical and numerical approach", 《IUTAM BOOKSERIES》 *
MAIINSKI FAKULTET等: "A COLLOCATED FINITE VOLUME METHOD FOR PREDICTING FLOWS AT ALL SPEEDS", 《INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS》 *
NACER E. EL KADRI E.等: "A finite element method for compressible and incompressible fows", 《SPRINGER NATURE SWITZERLAND AG 2020》 *
SHANTANU SHAHANE等: "Consistency and Convergence of a High Order Accurate Meshless Method for Solution of Incompressible Fluid Flows", 《ARXIV》 *
SHUYAO HU等: "Disturbance region update method for steady compressible flows", 《COMPUTER PHYSICS COMMUNICATIONS》 *
T. R. AKYLAS,C. C. MEI: "I-campus pro ject School-wide Program on Fluid Mechanics Modules on Waves in fluids", 《SCHOOL-WIDE MODULAR PROGRAM ON FLUID MECHANICS》 *
赵星宇: "临近空间高速飞行器抛罩过程的动网格方法对比", 《工程科技II辑》 *
邱周华: "谱元法求解不可压缩流体流动与传热问题及柱坐标系下线性稳定性分析", 《中国博士论文全文数据库》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114996634A (zh) * 2022-08-04 2022-09-02 中国空气动力研究与发展中心计算空气动力研究所 一种压力参考点的隐式处理方法、设备及介质

Also Published As

Publication number Publication date
CN114841095B (zh) 2022-09-09

Similar Documents

Publication Publication Date Title
Perić Analysis of pressure-velocity coupling on nonorthogonal grids
CN114841095B (zh) 一种不可压缩流动的扰动域推进方法
Chen et al. Improved particle swarm optimization-based form-finding method for suspension bridge installation analysis
CN111651831B (zh) 飞行器定常粘性可压缩绕流的分区扰动域更新计算方法
Groth et al. Parallel implicit adaptive mesh refinement scheme for body-fitted multi-block mesh
CN110826188B (zh) 一种基于gpu加速的天然气管网水力参数仿真方法
Andreianov et al. Qualitative behaviour and numerical approximation of solutions to conservation laws with non-local point constraints on the flux and modeling of crowd dynamics at the bottlenecks
CN111859645B (zh) 冲击波求解的改进musl格式物质点法
CN107562993B (zh) 一种提高低频隔声性能的平板结构优化设计方法
Kanazaki et al. Multi-objective aerodynamic optimization of elements' setting for high-lift airfoil using kriging model
CN108459993A (zh) 基于快速追峰采样的复杂高维系统优化方法
Danvin et al. Laminar to turbulent transition prediction in hypersonic flows with neural networks committee
CN101719194A (zh) 人工基因调控网络的模拟方法
Dai et al. Iterative implementation of an implicit–explicit hybrid scheme for hydrodynamics
Ren et al. Numerical stability analysis of shock-capturing methods for strong shocks I: second-order MUSCL schemes
Haga et al. Efficient solution techniques for high-order methods on 3-D anisotropic hybrid meshes
CN111259589A (zh) 一种考虑破损-安全的连续体频率约束拓扑优化设计方法
CN109117502B (zh) 一种降低平板结构瞬态声辐射的优化设计方法
CN114936522B (zh) 一种天然气系统动态建模方法
CN108170913A (zh) 一种基于图形法和二分法的斜拉桥拉索参数计算方法
CN111191358B (zh) 一种吸气式超声速导弹轨迹优化设计方法
Fournier et al. Multiplicative and additive parallel multigrid algorithms for the acceleration of compressible flow computations on unstructured meshes
David et al. Numerical solution of Euler equations using rotated-hybrid Riemann solver
CN116595745A (zh) 基于粒子水平集的强冲击物质界面高精度捕捉方法
Yamazaki Efficient robust design optimization by variable fidelity kriging model

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