CN112861263B - 一种适用于可压缩两相流的计算模拟方法 - Google Patents
一种适用于可压缩两相流的计算模拟方法 Download PDFInfo
- Publication number
- CN112861263B CN112861263B CN202110198896.9A CN202110198896A CN112861263B CN 112861263 B CN112861263 B CN 112861263B CN 202110198896 A CN202110198896 A CN 202110198896A CN 112861263 B CN112861263 B CN 112861263B
- Authority
- CN
- China
- Prior art keywords
- equation
- flux
- phase
- mixing
- flow
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000004364 calculation method Methods 0.000 title claims abstract description 29
- 230000005514 two-phase flow Effects 0.000 title claims abstract description 21
- 238000004088 simulation Methods 0.000 title abstract description 21
- 230000004907 flux Effects 0.000 claims abstract description 29
- 238000005094 computer simulation Methods 0.000 claims abstract description 5
- 239000012071 phase Substances 0.000 claims description 37
- 238000002156 mixing Methods 0.000 claims description 25
- 239000012530 fluid Substances 0.000 claims description 18
- 239000007788 liquid Substances 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 4
- 230000014509 gene expression Effects 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 2
- 239000007791 liquid phase Substances 0.000 claims description 2
- 230000035939 shock Effects 0.000 abstract description 14
- 238000011160 research Methods 0.000 abstract description 7
- 239000000126 substance Substances 0.000 abstract description 2
- 238000004880 explosion Methods 0.000 description 10
- XSNQEMWVLMRPFR-UHFFFAOYSA-N silver nitride Chemical compound [N-3].[Ag+].[Ag+].[Ag+] XSNQEMWVLMRPFR-UHFFFAOYSA-N 0.000 description 8
- 239000000243 solution Substances 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000000889 atomisation Methods 0.000 description 3
- 238000002485 combustion reaction Methods 0.000 description 3
- 239000002360 explosive Substances 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 230000003628 erosive effect Effects 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 1
- 238000012271 agricultural production Methods 0.000 description 1
- 239000003570 air Substances 0.000 description 1
- 238000012824 chemical production Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000012377 drug delivery Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005293 physical law Methods 0.000 description 1
- 238000009834 vaporization Methods 0.000 description 1
- 230000008016 vaporization Effects 0.000 description 1
Classifications
-
- 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
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Algebra (AREA)
- Mathematical Physics (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出了一种适用于可压缩两相流的计算模拟方法,该方法采用五方程模型描述可压缩两相界面的流动,将五方程与刚性气体状态方程耦合使得方程组封闭。本发明采用有限体积法离散化方程组,采用三阶TVD格式对原始变量进行空间重构,采用HLLC近似黎曼求解器求解单元边界的数值通量,采用三阶显式的龙格‑库塔格式进行时间推进。本发明的计算模拟方法可以准确地捕捉到激波间断、物质接触间断和膨胀波,为两相流的数值模拟研究提供了可靠的数值模拟平台。
Description
技术领域
本发明属于计算流体力学技术领域,具体涉及一种适用于可压缩两相流的计算模拟方法。
背景技术
可压缩两相流是一个活跃且具有挑战性的研究领域,在理论研究以及工程应用中都具有重要的意义,并且广泛存在于航空航天、化工生产、生物医学、农业生产等领域中。
在航空航天等领域,激波与液滴的相互作用是典型的两相流问题。当飞行器高速地穿过云层或者雨区时,飞机表面附近产生的激波与空气中悬浮的雨滴发生相互作用,雨滴在破碎过程中会使飞行器表面受到损伤,即雨滴侵蚀现象。深入研究激波与液滴相互作用的物理机制,可有效地减小雨滴侵蚀对飞行器带来的影响。此外,超燃冲压发动机燃烧室在超声速条件下的雾化研究是保证其在启动工况下可靠运行的关键。当气流以超声速流入燃烧室时,液体燃料喷射需要在极短的时间内完成雾化、蒸发以及复杂的燃烧过程,这涉及到激波与液体燃料的复杂相互作用。因此,研究并获得激波与液滴的物理规律对指导高速飞行器、超燃冲压发动机以及雾化装置的设计有着重要意义。
在叶轮机械与舰船结构等领域,空化是典型的两相流流动过程。空化是流体中由于压力极低而产生的气泡快速演化现象。气泡的突然产生以及突然溃灭会造成压力突变,进而导致严重的机械损伤,降低机械的使用寿命。深入研究空化现象的机理,对减缓机械磨损有着重要的指导意义。
在生物医疗等领域,可压缩两相流广泛存在于激波碎石、药物输运等方面。激波碎石是利用冲击波对体内的结石进行无创粉碎,所产生的压力脉冲与结石附近的气泡相互作用,使得结石碎裂,这是否成功取决于激波与气泡相互作用后气泡的动力学行为。因此,可压缩两相流的数值模拟具有重要的研究意义。
相比于单相流动,可压缩两相流动存在着复杂的波系结构,涉及到强激波和大密度比两相界面的时空演化过程,对数值模拟研究提出了巨大的挑战。因此,针对可压缩两相流界面问题的数值模拟方法,发展可靠的数值方法尤为重要。
发明内容
本发明要解决的问题是提供一种适用于可压缩两相流的计算模拟方法,该方法能够高效高精度地模拟两相流的流场演化过程。
为了实现上述目的,本发明提出的适用于可压缩两相流的计算模拟方法,具体步骤如下:
步骤1:根据需要计算模拟的可压缩两相流确定计算域,根据计算网格点数将计算域离散为均匀的网格点,并输入网格点信息,包括所有网格点上的速度、密度和压力;
步骤2:采用有限体积法对控制方程组进行离散,获得方程组的半离散形式;
步骤3:采用三阶TVD格式对离散后的控制方程组进行空间重构,获得单元边界两侧的原始变量;
步骤4:采用HLLC黎曼求解器求解单元边界处的通量;
步骤5:将刚性气体状态方程与流动求解器相结合以模拟气液两相;
步骤6:模拟步骤5中气液两相的混合规律;
步骤7:根据界面通量求解残差,并将半离散有限体积格式转化为时空全离散有限体积格式,全流场进行时间推进求解,得到最终的流场解。
进一步的,步骤2中,采用五方程模型描述两相流体的流动,所述五方程模型包括两相流体各自的连续方程、统一的动量方程、统一的能量方程和体积分数输运方程。
进一步的,步骤2中,不考虑粘性效应和表面张力的控制方程写为以下形式:
其中,U为控制方程的守恒变量,F为x方向的通量,G为y方向的通量,S为源项,定义如下:
ρ1和ρ2分别表示两相的密度,α1和α2分别表示两相的体积分数,下标1和2分别表示两相的物理量。ρ,u,v,p分别表示两相流体的混合密度,混合速度、混合压力,ρu和ρv表示两相流体的总动量,E表示两相流体的总动能。
进一步的,步骤2中,考虑计算的网格单元为[xi-Δx/2,xi+Δx/2],不考虑体积输运方程的控制方程组的半离散形式为:
其中,表示守恒变量的单元平均值,Fi-1/2和Fi+1/2分别表示单元左边界i-1/2和右边界i+1/2的数值通量;对于体积分数输运方程进行链式求导
得到方程的半离散形式为:
其中为体积分数的单元平均值;公式中右端第一项为体积分数在单元边界处的数值通量,ui+1/2和ui-1/2为单元边界的速度;Δx为网格单元尺寸。
进一步的,步骤3中,满足TVD性质的三阶重构方法形式为
其中,χL、χR和的计算公式如下:
进一步的,步骤4中,通过空间重构得到单元边界处的守恒变量Ul、Ur和数值通量Fl、Fr,再根据波系近似得到HLLC黎曼求解器求解的单元边界数值通量FHLLC为:
Fl *和分别为中间通量,S*为中间波速:
Ul *和Ur *分别是中间状态向量,计算方式为:
Sl和Sr分别为左波速和右波速,根据公式
Sl=min(Sl *,0)Sr=max(Sr *,0)
进行计算:
而方程离散速度项uHLLC根据公式
计算。
进一步的,步骤5中由刚性气体状态方程可知
p=(γ-1)ρe-γπ∞
而刚性气体状态方程的声速公式为:
得到总能E为
其中γ是比热比,π∞为参考压力。
进一步的,步骤6中两相流体混合密度、混合体积分数和混合内能的表达式为:
ρ=ρ1α1+ρ2α2
α1+α2=1
ρe=ρ1α1e1+ρ2α2e2
刚性气体状态方程中的两个状态参数Γ和Π的混合规则如下:
其中,π∞,1和π∞,2分别表示两相的参考压力。
有益效果
本发明提出了一种适用于可压缩两相流的计算模拟方法,该方法采用五方程模型描述可压缩两相界面的流动,将五方程与刚性气体状态方程耦合使得方程组封闭。本发明采用有限体积法离散化方程组,采用三阶TVD格式对原始变量进行空间重构,采用HLLC近似黎曼求解器求解单元边界的数值通量,采用三阶显式的龙格-库塔格式进行时间推进。本发明的计算模拟方法可以准确地捕捉到激波间断、物质接触间断和膨胀波,为两相流的数值模拟研究提供了可靠的数值模拟平台。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1为水下三氮化银***问题的示意图。
图2为t=14μs时水下三氮化银***问题的实验结果(a)和计算模拟结果(b)。
图3为t=28μs时水下三氮化银***问题的实验结果(a)和计算模拟结果(b)。
图4为t=42μs时水下三氮化银***问题的实验结果(a)和计算模拟结果(b)。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
水下***是典型的可压缩两相流问题。例如本实施例中给出的水下三氮化银***问题的数值模拟计算区域如图1所示,单位为m。计算区域为{(x,y)|0≤x≤0.35,0≤y≤0.35}(单位:m),分隔空气与水的水平面位于y=0.175m,***核的直径为0.01m,中心位置为(0.175m,0.125m)。计算采用的网格分辨率为3000×1000,CFL数取值为0.1。该算例具体的初始条件如表1所示,其中,水和***核的参考压力π∞均为6×108,空气的参考压力π∞为0;水和***核的比热γ取值为4.4,空气的比热比γ取值为1.4。计算区域的下边界采用反射边界条件,左右边界和上边界采用远场边界条件。
表1水下三氮化银***问题的计算参数
应用本发明的计算方法,对水下三氮化银***问题进行数值模拟,具体过程如下:
步骤1、网格信息输入以及初始化。
首先确定计算域{(x,y)|0≤x≤0.35,0≤y≤0.35},令网格单元数为1000×1000,对计算域进行网格划分,将所有网格点上的速度、密度和压力信息输入计算方法程序中并完成初始化。
步骤2:采用有限体积法对控制方程组进行离散,获得方程组的半离散形式。
本发明采用五方程模型描述两相流体的流动,即假设两相具有相同的速度以及压力,采用各自的连续方程、统一的动量方程、统一的能量方程和体积分数输运方程。以二维为例,不考虑粘性效应和表面张力的控制方程可写为以下形式:
其中,U为控制方程的守恒变量,F为x方向的通量,G为y方向的通量,S为源项,上述4个量的具体定义如下:
上式中,ρ1和ρ2分别表示两相的密度,α1和α2分别表示两相的体积分数,下标1和2分别表示两相的物理量。ρ,u,v,p分别表示两相流体的混合密度,混合速度、混合压力,ρu和ρv表示两相流体的总动量,E表示两相流体的总动能。
以一维为例,考虑计算的网格单元为[xi-Δx/2,xi+Δx/2],Δx为网格单元尺寸,则不考虑体积输运方程时,控制方程组的半离散形式如下:
其中,表示守恒变量的单元平均值,Fi-1/2和Fi+1/2分别表示单元左边界i-1/2和右边界i+1/2的数值通量。对于体积分数输运方程,为使得其求解与欧拉方程组保持一致,对体积分数输运方程进行链式求导如下:
则上述方程的半离散形式可写为:
其中,为体积分数的单元平均值。上式中右端第一项为体积分数在单元边界处的数值通量,ui+1/2和ui-1/2为单元边界的速度。
步骤3:采用三阶TVD格式对离散后的控制方程组进行空间重构,获得单元边界两侧的原始变量。
满足TVD性质的三阶重构方法形式如下:
其中,χL、χR和ri L、ri R的计算公式如下:
步骤4:采用HLLC黎曼求解器求解单元边界处的通量
通过使用HLLC黎曼求解器,可将网格单元附近的流动演变划分为三个波系,左右两测的波系分别为激波或膨胀波,中间的波系为接触间断。本步骤通过空间重构得到单元边界处的守恒变量Ul、Ur和数值通量Fl、Fr后,再根据波系近似可以得到HLLC黎曼求解器求解的单元边界数值通量FHLLC为:
上式中,Fl *和分别为中间通量,S*为中间波速,计算公式如下:
Fk *=Fk+Sk(Uk *-Uk),k=l,r
上式中,Ul *和Ur *分别是中间状态向量,计算方式如下:
其中,Sl和Sr分别为左波速和右波速,可以采用以下公式进行计算:
Sl=min(Sl *,0) Sr=max(Sr *,0)
方程离散速度项uHLLC的计算公式如下:
步骤5:将刚性气体状态方程与流动求解器相结合以模拟气液两相
刚性气体状态方程将流体的密度、压力以及内能紧密联系起来,当前被广泛应用于可压缩多组分流动模拟。本发明将刚性气体状态方程与流动求解器相结合来模拟气液两相。由刚性气体状态方程可知:
p=(γ-1)ρe-γπ∞
刚性气体状态方程的声速计算公式如下:
因此总能E可以表示为:
γ和π∞均是根据经验确定的拟合参数。其中,γ是比热比,π∞为参考压力,是与分子间引力相关的参数。考虑到所模拟的流体是不可混溶的,每一种流体均分别遵循刚性气体状态方程。对于空气,由于其参考压力为0,刚性气体状态方程即为理想气体方程。
步骤6:模拟两相的混合规律
在扩散界面区域内,流体混合物需满足一定的混合规则。对于空气、水和***核,混合密度、混合体积分数和混合内能的表达式定义如下:
ρ=ρ1α1+ρ2α2
α1+α2=1
ρe=ρ1α1e1+ρ2α2e2
刚性气体状态方程中的两个状态参数Γ和Π的混合规则如下:
其中,π∞,1和π∞,2分别表示两相的参考压力,其为与分子之间引力相关的参数,本实施例中均为零。
步骤7:根据界面通量求解残差,并将半离散有限体积格式转化为时空全离散有限体积格式,全流场进行时间推进求解,得到最终的流场解。
图2-图4分别给出了在t=14μs、t=28μs和t=42μs时刻水下三氮化银***问题的实验结果和计算模拟结果。可以看出,本发明的数值模拟方法可准确地捕捉到激波、膨胀波和两相界面等复杂的流动结构。将实验结果和数值模拟结果进行对比,可以看出数值模拟结果与实验情况吻合良好,可以很好地反应真实流动,验证了本发明所采用的数值模拟方法的有效性。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明未详细描述内容为本领域技术人员公知技术。
Claims (1)
1.一种适用于可压缩两相流的计算模拟方法,其特征在于:包括以下步骤:
步骤1:根据需要计算模拟的可压缩两相流确定计算域,根据计算网格点数将计算域离散为均匀的网格点,并输入网格点信息,包括所有网格点上的速度、密度和压力;
步骤2:采用有限体积法对控制方程组进行离散,获得方程组的半离散形式:
采用五方程模型描述两相流体的流动,所述五方程模型包括两相流体各自的连续方程、统一的动量方程、统一的能量方程和体积分数输运方程;
针对二维情况,不考虑粘性效应和表面张力的控制方程写为以下形式:
其中,U为控制方程的守恒变量,F为x方向的通量,G为y方向的通量,S为源项,定义如下:
ρ1和ρ2分别表示两相的密度,α1和α2分别表示两相的体积分数,下标1和2分别表示两相的物理量;ρ,u,v,p分别表示两相流体的混合密度,混合速度、混合压力,ρu和ρv表示两相流体的总动量,E表示两相流体的总动能;
针对一维情况,考虑计算的网格单元为[xi-Δx/2,xi+Δx/2],不考虑体积输运方程的控制方程组的半离散形式为:
其中,表示守恒变量的单元平均值,Fi-1/2和Fi+1/2分别表示单元左边界i-1/2和右边界i+1/2的数值通量;对于体积分数输运方程进行链式求导
得到方程的半离散形式为:
其中为体积分数的单元平均值;公式中右端第一项为体积分数在单元边界处的数值通量,ui+1/2和ui-1/2为单元边界的速度;Δx为网格单元尺寸;
步骤3:采用三阶TVD格式对离散后的控制方程组进行空间重构,获得单元边界两侧的原始变量;
满足TVD性质的三阶重构方法形式为:
其中,χL、χR和ri L、ri R的计算公式如下:
步骤4:采用HLLC黎曼求解器求解单元边界处的通量:
通过空间重构得到单元边界处的守恒变量Ul、Ur和数值通量Fl、Fr,再根据波系近似得到HLLC黎曼求解器求解的单元边界数值通量FHLLC为:
Fl *和Fr *分别为中间通量,S*为中间波速:
Fk *=Fk+Sk(Uk *-Uk),k=l,r
Ul *和Ur *分别是中间状态向量,计算方式为:
Sl和Sr分别为左波速和右波速,根据公式
Sl=min(Sl *,0) Sr=max(Sr *,0)
进行计算:
而方程离散速度项uHLLC根据公式
计算;
步骤5:将刚性气体状态方程与流动求解器相结合以模拟气液两相:
由刚性气体状态方程得到:
p=(γ-1)ρe-γπ∞
而刚性气体状态方程的声速公式为:
得到总能E为
其中γ是比热比,π∞为参考压力;
步骤6:模拟步骤5中气液两相的混合规律:
两相流体混合密度、混合体积分数和混合内能的表达式为:
ρ=ρ1α1+ρ2α2
α1+α2=1
ρe=ρ1α1e1+ρ2α2e2
刚性气体状态方程中的两个状态参数Γ和Π的混合规则如下:
其中,π∞,1和π∞,2分别表示两相的参考压力;
步骤7:根据界面通量求解残差,并将半离散有限体积格式转化为时空全离散有限体积格式,全流场进行时间推进求解,得到最终的流场解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110198896.9A CN112861263B (zh) | 2021-02-22 | 2021-02-22 | 一种适用于可压缩两相流的计算模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110198896.9A CN112861263B (zh) | 2021-02-22 | 2021-02-22 | 一种适用于可压缩两相流的计算模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112861263A CN112861263A (zh) | 2021-05-28 |
CN112861263B true CN112861263B (zh) | 2024-02-13 |
Family
ID=75988554
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110198896.9A Active CN112861263B (zh) | 2021-02-22 | 2021-02-22 | 一种适用于可压缩两相流的计算模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112861263B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113688580B (zh) * | 2021-07-21 | 2024-03-08 | 上海交通大学 | 气液两相流相界面密度计算方法、装置、设备和存储介质 |
CN114169184B (zh) * | 2021-10-19 | 2024-06-21 | 北京理工大学 | 一种高精度自适应有限体积-有限差分耦合数值模拟方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012031398A1 (en) * | 2010-09-09 | 2012-03-15 | Tianjin Aerocode Engineering Application Software Development Inc. | Numerical method for simulating subsonic flows based on euler equations in lagrangian formulation |
NO20121424A1 (no) * | 2012-11-27 | 2014-05-28 | Sinvent As | Fremgangsmåte for simulering av flerfasefase fluidstrømninger i rørledninger |
NO20140079A1 (no) * | 2014-01-24 | 2015-07-27 | Ledaflow Technologies Da | Fremgangsmåte for bestemmelse av strømningsparametere for en flerfasestrømning i en rørledningsseksjon |
CN105468844A (zh) * | 2015-11-23 | 2016-04-06 | 河海大学 | 管道内水-气耦合瞬变流的模拟方法 |
CN105653860A (zh) * | 2015-12-30 | 2016-06-08 | 南京航空航天大学 | 可压缩气体与不可压缩液体多介质界面追踪数值方法 |
CN109214082A (zh) * | 2018-09-03 | 2019-01-15 | 北京理工大学 | 一种近场水下***冲击波载荷的高精度数值模拟方法 |
CN109918787A (zh) * | 2019-03-08 | 2019-06-21 | 河海大学 | 基于有限体积法的输水管道内水气两相均质流的模拟方法 |
CN111428401A (zh) * | 2020-03-18 | 2020-07-17 | 长江水利委员会长江科学院 | 一种堰塞湖溃决过程模拟方法 |
-
2021
- 2021-02-22 CN CN202110198896.9A patent/CN112861263B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012031398A1 (en) * | 2010-09-09 | 2012-03-15 | Tianjin Aerocode Engineering Application Software Development Inc. | Numerical method for simulating subsonic flows based on euler equations in lagrangian formulation |
NO20121424A1 (no) * | 2012-11-27 | 2014-05-28 | Sinvent As | Fremgangsmåte for simulering av flerfasefase fluidstrømninger i rørledninger |
NO20140079A1 (no) * | 2014-01-24 | 2015-07-27 | Ledaflow Technologies Da | Fremgangsmåte for bestemmelse av strømningsparametere for en flerfasestrømning i en rørledningsseksjon |
CN105468844A (zh) * | 2015-11-23 | 2016-04-06 | 河海大学 | 管道内水-气耦合瞬变流的模拟方法 |
CN105653860A (zh) * | 2015-12-30 | 2016-06-08 | 南京航空航天大学 | 可压缩气体与不可压缩液体多介质界面追踪数值方法 |
CN109214082A (zh) * | 2018-09-03 | 2019-01-15 | 北京理工大学 | 一种近场水下***冲击波载荷的高精度数值模拟方法 |
CN109918787A (zh) * | 2019-03-08 | 2019-06-21 | 河海大学 | 基于有限体积法的输水管道内水气两相均质流的模拟方法 |
CN111428401A (zh) * | 2020-03-18 | 2020-07-17 | 长江水利委员会长江科学院 | 一种堰塞湖溃决过程模拟方法 |
Non-Patent Citations (4)
Title |
---|
R.Abgrall等.Stochastic Discrete Equation Method (sDEM) for two-phase flows.Journal of Computational Physics.2015,第2015卷(第299期),281-306. * |
基于HLLC格式的理想MHD方程数值模拟;张向洪;;航空动力学报(05);全文 * |
基于近似黎曼求解器的三维浅水方程组求解方法;卢新华;;人民长江(20);全文 * |
复杂地形上非恒定浅水二维流动数值模拟;毕胜;周建中;张华杰;宋利祥;陈生水;;水动力学研究与进展A辑(01);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112861263A (zh) | 2021-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112861263B (zh) | 一种适用于可压缩两相流的计算模拟方法 | |
Bardina et al. | Improved subgrid-scale models for large-eddy simulation | |
Deere et al. | A computational study of a dual throat fluidic thrust vectoring nozzle concept | |
CN110222467B (zh) | 一种水下发动机喷管摆动高速射流特性的分析方法 | |
CN109359325B (zh) | 关于多喷管火箭流场及对流/辐射耦合换热的仿真方法 | |
CN108416127A (zh) | 潜艇高压吹除管路***多物理场耦合仿真方法及*** | |
Chahine | Numerical simulation of bubble flow interactions | |
CN105302997A (zh) | 一种基于三维cfd的液柱分离-弥合水锤的模拟方法 | |
Drikakis et al. | Azure: An advanced cfd software suite based on high-resolution and high-order methods | |
Nezamirad et al. | Utilization of Schnerr-Sauer cavitation model for simulation of cavitation inception and super cavitation | |
Nezamirad et al. | Effect of Needle Height on Formation of Cavitation in a Six-Hole Diesel Injector Nozzle | |
Naitoh et al. | Simultaneous attainment of light-weight, high-efficiency, and low noise: by the supermultijet-twister engine working from startup to hypersonic scram mode | |
Zhu et al. | Numerical investigation about the unsteady behavior of a free submerged cavitation jet using the SBES approach | |
Yang et al. | Characteristics and load reduction method of the deep-sea implosion of the ellipsoidal and egg-shaped pressure hulls | |
Zmijanovic et al. | Thrust vectoring effects of a transverse gas injection into a supersonic cross flow of an axisymmetric convergent-divergent nozzle | |
Wang et al. | Investigation of self-excited oscillation chamber cavitation effect with special emphasis on wall shape | |
Zhang et al. | Research on the thrust of a high-pressure water jet propulsion system | |
John et al. | Energy Addition–Based Virtual Cowl for Performance Enhancement of Scramjet Intake | |
Ogawa et al. | Two-dimensional computational fluid dynamics analysis of internal flow in a simplified solid rocket motor combustion chamber | |
Zhao et al. | Investigations on the effects of obstacles on the surfaces of blades of the centrifugal pump to suppress cavitation development | |
Sidorenko et al. | Two-dimensional gas dynamics modeling of the relaxation of particles behind the transmitted shock wave | |
He et al. | A refined numerical investigation of a large equivalent shallow-depth underwater explosion | |
Xutao et al. | Numerical Study on Flow-Induced Vibration of Ejector Structure. | |
Boyd et al. | Shock wave interaction with a transcritical fuel droplet | |
Gopalakrishnan et al. | Numerical simulation of flow in fluidic valves in rotating detonation engines |
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 |