CN105183965B - 用于预测雾化过程的大涡模拟方法 - Google Patents

用于预测雾化过程的大涡模拟方法 Download PDF

Info

Publication number
CN105183965B
CN105183965B CN201510534392.4A CN201510534392A CN105183965B CN 105183965 B CN105183965 B CN 105183965B CN 201510534392 A CN201510534392 A CN 201510534392A CN 105183965 B CN105183965 B CN 105183965B
Authority
CN
China
Prior art keywords
liquid
velocity field
equation
liquid velocity
gas
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
Application number
CN201510534392.4A
Other languages
English (en)
Other versions
CN105183965A (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201510534392.4A priority Critical patent/CN105183965B/zh
Publication of CN105183965A publication Critical patent/CN105183965A/zh
Application granted granted Critical
Publication of CN105183965B publication Critical patent/CN105183965B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种用于预测雾化过程的大涡模拟方法,基于真实速度场构造虚拟的液相速度场,并将构造的液相速度场应用于流动控制方程和界面输运方程的求解,模拟气液两相流的实时动态过程,以准确预测液滴破碎过程和液体射流的雾化过程。本发明提高了两相流模拟的计算精度和稳定性,可以准确计算和预测液滴的破碎过程和液柱射流的雾化过程。

Description

用于预测雾化过程的大涡模拟方法
技术领域
本发明涉及流体控制领域,特别地,涉及一种用于预测雾化过程的大涡模拟方法。
背景技术
在发动机燃烧室内,液体燃料的雾化决定了燃料与空气的混合效果,进而影响燃烧性能。雾化过程非常复杂,多种不稳定性(Kelvin-Helmholtz instability,Rayleigh-Taylor instability,Plateau-Rayleigh instability)同时存在,并伴有强烈的湍流,使得理论分析不具可行性。关于雾化已开展了大量的试验研究,但是由于雾化过程形成的液雾遮挡了液柱初始破碎过程,在观察和测量上造成了很大的困难。自从二十世纪七十年代,两相流的数值仿真取得了很大的进步,加深了对雾化机理的认识。
流体计算力学方法分为三种:雷诺平均方法,大涡模拟,直接数值模拟。雷诺平均方法只求解平均速度场,模化湍流运动对流场的影响。大涡模拟求解大尺度涡结构,模化小尺度涡结构对流动的影响。直接数值模拟求解所有尺度的涡结构。雾化过程中,湍流中的大尺度涡可扰动两相流界面,显著地影响液体射流的破碎过程,限制了雷诺平均方法的应用。由于大涡模拟的计算量比直接数值模拟小很多,大涡模拟更适合工程应用。
为了精确求解液柱和液滴的破碎过程,需跟踪液气界面。流行的界面跟踪方法有:Volume of fluid(VOF),Level Set(LS),Coupled LS and VOF(CLSVOF)。其中,VOF方法(流体体积法)在20世纪70年代末由Hirt和Nichols等最先提出,基本思想是在欧拉网格系统上定义一个函数,根据每个网格内所含某种物质的体积量来定义在此网格上的值,然后用体积跟踪的方法求解方程,VOF方法可以精确地保证质量守恒,但是VOF函数的不连续性导致界面的构造很复杂且容易破碎。LS方法(水平集方法)可以很容易地构造界面,但是所得界面包围的液体质量不守恒。CLSVOF方法(水平集复合流体体积方法)可以很好地结合VOF和LS方法的优点,得到了广泛的应用。
由于液气界面两侧密度和流体粘性系数的不连续性,在求解控制方程时,常规的数值离散方法误差大,并造成算法的不稳定,液气密度比越大,算法越不稳定。为了得到收敛的结果,很多已发表的文献在数值仿真中采用较低的液气密度比,但是大多数雾化试验是采用高密度液体(如水、煤油、酒精)在大气环境中进行的,具有较高的液气密度比,现有技术无法将数值仿真结果与实验结果进行比较。因此,现有的两相流大涡模拟存在算法复杂、且外延液体速度不满足连续性方程,离散数值误差大等问题导致无法准确预测液滴破碎过程和液体射流的雾化过程的缺陷。
发明内容
本发明提供了一种用于预测雾化过程的大涡模拟方法,以解决现有两相流模拟方法导致的两相流液滴破碎及液体射流的雾化过程难以准确预测的技术问题。
本发明采用的技术方案如下:
一种用于预测雾化过程的大涡模拟方法,本发明方法基于真实速度场构造虚拟的液相速度场,并将构造的液相速度场应用于流动控制方程和界面输运方程的求解,模拟气液两相流的实时动态过程,以准确预测液滴破碎过程和液体射流的雾化过程。
进一步地,本发明用于预测雾化过程的大涡模拟方法包括:
步骤S10,根据水平集LS函数Φn表示两相流界面,通过真实速度场Un和液相速度场UL n求解两相流控制方程,获得下一时间步对应的真实速度场Un+1
步骤S20,通过外延方法构造n+1时间步的液相速度场UL n+1
步骤S30,通过外延液相速度场去散度化方法使UL n+1满足连续性方程:
步骤S40,利用构造的液相速度场UL n+1,通过水平集复合流体体积CLSVOF方法求解LS函数和VOF函数的输运方程,获得下一时刻的LS函数Φn+1和流体体积VOF函数Fn+1
步骤S50,由气相转变为液相的控制体中,将真实速度场Un+1重置为液相速度场UL n +1,即Un+1=UL n+1
重复以上步骤S10至S50,以模拟两相流的实时动态过程。
进一步地,所述步骤S10中,对所述两相流控制方程进行空间过滤处理。
本发明具有以下有益效果:
本发明用于预测雾化过程的大涡模拟方法,基于真实速度场构造虚拟的液相速度场,并将构造的液相速度场应用于流动控制方程和界面输运方程的求解,模拟气液两相流的实时动态过程,提高了计算精度和稳定性,可以准确计算和预测液滴的破碎过程和液柱射流的雾化过程。
除了上面所描述的目的、特征和优点之外,本发明还有其它的目的、特征和优点。下面将参照图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本发明优选实施例大涡模拟方法的流程示意图;
图2是本发明优选实施例计算变量的分布示意图;
图3是本发明优选实施例液相速度场的构造示意图。
具体实施方式
以下结合附图对本发明的实施例进行详细说明,但是本发明可以由权利要求限定和覆盖的多种不同方式实施。
本发明的优选实施例提供了一种用于预测雾化过程的大涡模拟方法,本发明方法基于真实速度场构造虚拟的液相速度场,并将构造的液相速度场应用于流动控制方程和界面输运方程的求解,模拟气液两相流的实时动态过程,以准确预测液滴破碎过程和液体射流的雾化过程。
参照图1,本实施例控制方法包括:
步骤S10,根据水平集LS函数Φn表示两相流界面,通过真实速度场Un和液相速度场UL n求解两相流控制方程,获得下一时间步对应的真实速度场Un+1
步骤S20,通过外延方法构造n+1时间步的液相速度场UL n+1
步骤S30,通过外延液相速度场去散度化方法使UL n+1满足连续性方程:
步骤S40,利用构造的液相速度场UL n+1,通过水平集复合流体体积CLSVOF方法求解LS函数和VOF函数的输运方程,获得下一时刻的LS函数Φn+1和流体体积VOF函数Fn+1
步骤S50,由气相转变为液相的控制体中,将真实速度场Un+1重置为液相速度场UL n +1,即Un+1=UL n+1
重复以上步骤S10至S50,以模拟两相流的实时动态过程。
作为一种较佳的方式,为了跟踪液气边界,引入两个函数:LS函数Φ和VOF函数F。LS函数Φ是到液气界面的变号距离函数。Φ=0代表液气界面;在液体中Φ>0;在气体中Φ≤0。VOF函数F是每个计算单元中液体体积百分比。
优选地,对所述两相流控制方程(Navier-Stokes方程)进行空间过滤处理。经空间过滤后,连续方程变为:
其中,Ui为速度分量,xi为位置坐标。
通过Smagorinsky涡粘性模型模化亚格子应力,动量方程变为:
其中,P为压强,t为时间,ρ为密度,gi为重力分量,为表面张力。τij分别为粘性应力张量和亚格子应力张量,且计算如下:
ρ=ρG+(ρLG)H(φ) μ=μG+(μLG)H(φ)
μ和μr分别表示动力粘性系数和亚格子粘性系数,过滤宽度Δ取为当地计算单元体积的立方根,Sij为应变张量,CS为Smagorinsky系数。下标G和L分别表示气体和液体,H(φ)为Heaviside函数,表示如下:
表面张力为:
其中σ为表面张力系数,κ为曲率,ni为法向量分量。
VOF函数的控制方程为:
LS函数的控制方程为:
作为一种较佳的方式,为了计算下一时间步的速度场:
首先,通过对流项、扩散项和重力项计算中间速度场:
其次,中间速度场通过压力项修正得到n+1时间步的速度场:
因为n+1时间步的速度场满足连续性方程,通过对上一方程求散度可得如下压力泊松方程(通过此泊松方程可求解n+1时间步的的压力场):
参照图2,图2展示了计算变量的分布示意图,LS函数Φ、压强P、VOF函数F位于计算单元中心,速度以交错方式分布于相应的计算单元表面上。u,v是速度场在x和y方向上的分量。
参照图3,液相速度初始化为动量方程得到的速度场
如果
如果
气体中(φ≤0)的液相速度是通过将沿界面法向方向从液体向气体外延得到的,求解下面的外延方程到稳态:
一阶前向欧拉方法用于时间离散,以液相速度分量uL为例:
伪时间步长Δτ=0.3min(Δxi-1,Δxi,Δyj-1,Δyj,Δyj+1)。一阶迎风格式用于空间离散:
作为较佳的方式,外延的液体速度应满足连续性条件
首先计算单元(i,j)中的速度源项:
修正气体中的液相速度以满足连续性条件:
A=ae|nx|i-1/2,jΔyj+aw|nx|i+1/2,jΔyj+as|ny|i,j-1/2Δxi+an|ny|i,j+1/2Δxi
从以上的描述可以得知,本发明基于真实速度场构造虚拟的液相速度场,并将构造的液相速度场应用于流动控制方程和界面输运方程,模拟气液两相流的实时动态过程,提高了计算精度和稳定性,可以准确计算和预测液滴的破碎过程和液柱射流的雾化过程。且本发明通过液相速度外延算法构造下一时间步的液相速度场,通过外延液相速度无散度化处理,进一步控制了两相流模拟时由于离散数值导致的误差,提高了计算精度和稳定性。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.一种用于预测雾化过程的大涡模拟方法,其特征在于,基于真实速度场构造虚拟的液相速度场,并将构造的液相速度场应用于流动控制方程和界面输运方程的求解,模拟气液两相流的实时动态过程,以准确预测液滴破碎过程和液体射流的雾化过程;
所述基于真实速度场构造虚拟的液相速度场包括以下步骤:
液体中的液相速度初始化为动量方程得到的速度场
如果
如果u和v分别是在x和y方向上的分量;
气体中的液相速度是通过将沿界面法向方向从液体向气体外延得到的,求解下面的外延方程到稳态:
一阶前向欧拉方法用于时间离散
伪时间步长Δτ=0.3min(Δxi-1,Δxi,Δyj-1,Δyj,Δyj+1),
一阶迎风格式用于空间离散:
外延的液体速度满足连续性条件:
首先计算单元(i,j)中的速度源项:
修正气体中的液相速度以满足连续性条件:
A=ae|nx|i-1/2,jΔyj+aw|nx|i+1/2,jΔyj+as|ny|i,j-1/2Δxi+an|ny|i,j+1/2Δxi
其中,函数Φ是到液气界面的变号距离函数,xi为位置坐标,ni为法向量分量。
2.根据权利要求1所述的用于预测雾化过程的大涡模拟方法,其特征在于,所述控制方法包括:
步骤S10,根据水平集LS函数Φn表示两相流界面,通过真实速度场Un和液相速度场ULn求解两相流控制方程,获得下一时间步对应的真实速度场Un+1
步骤S20,通过外延方法构造n+1时间步的液相速度场ULn+1
步骤S30,通过外延液相速度场去散度化方法使ULn+1满足连续性方程:▽·ULn+1=0;
步骤S40,利用构造的液相速度场ULn+1,通过水平集复合流体体积CLSVOF方法求解LS函数和VOF函数的输运方程,获得下一时刻的LS函数Φn+1和流体体积VOF函数Fn+1
步骤S50,由气相转变为液相的控制体中,将真实速度场Un+1重置为液相速度场ULn+1,即Un+1=ULn+1
重复以上步骤S10至S50,以模拟两相流的实时动态过程。
3.根据权利要求2所述的用于预测雾化过程的大涡模拟方法,其特征在于,
所述步骤S10中,对所述两相流控制方程进行空间过滤处理。
CN201510534392.4A 2015-08-27 2015-08-27 用于预测雾化过程的大涡模拟方法 Active CN105183965B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510534392.4A CN105183965B (zh) 2015-08-27 2015-08-27 用于预测雾化过程的大涡模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510534392.4A CN105183965B (zh) 2015-08-27 2015-08-27 用于预测雾化过程的大涡模拟方法

Publications (2)

Publication Number Publication Date
CN105183965A CN105183965A (zh) 2015-12-23
CN105183965B true CN105183965B (zh) 2019-07-12

Family

ID=54906044

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510534392.4A Active CN105183965B (zh) 2015-08-27 2015-08-27 用于预测雾化过程的大涡模拟方法

Country Status (1)

Country Link
CN (1) CN105183965B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110320189B (zh) * 2019-06-27 2020-04-24 中国科学院力学研究所 一种液体燃料雾化过程中两相测量方法及系统
CN110414141B (zh) * 2019-07-30 2022-11-04 辽宁工程技术大学 可压流体跨音速流动过程中的液滴雾化三维数值模拟方法
CN112069689B (zh) * 2020-09-10 2022-02-08 西北工业大学 一种航空发动机燃油雾化特性的仿真方法及系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104239640A (zh) * 2014-09-18 2014-12-24 中国人民解放军国防科学技术大学 不可压流大涡模拟的湍流入口条件生成方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104239640A (zh) * 2014-09-18 2014-12-24 中国人民解放军国防科学技术大学 不可压流大涡模拟的湍流入口条件生成方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A coupled volume-of-fluid and level set (VOSET) method for computing incompressible two-phase flows;D.L.Sun 等;《International Journal of Heat and Mass Transfer》;20100131;第53卷(第4期);第645-655页
基于LES-VOF模型的燃油射流雾化过程模拟;彭天鹏等;《南昌大学学报(工科版)》;20091231;第31卷(第4期);正文第1-4节
追踪不可压缩两相流相界面的CLSVOF方法;宋云超等;《农业机械学报》;20110731;第42卷(第7期);正文第1-5节

Also Published As

Publication number Publication date
CN105183965A (zh) 2015-12-23

Similar Documents

Publication Publication Date Title
Eisenschmidt et al. Direct numerical simulations for multiphase flows: An overview of the multiphase code FS3D
Rosti et al. Numerical simulations of emulsions in shear flows
Ullrich et al. High-order finite-volume methods for the shallow-water equations on the sphere
Zu et al. Phase-field-based lattice Boltzmann model for incompressible binary fluid systems with density and viscosity contrasts
Balabel et al. On the performance of linear and nonlinear k–ε turbulence models in various jet flow applications
Nonomura et al. A simple interface sharpening technique with a hyperbolic tangent function applied to compressible two-fluid modeling
de Langavant et al. Level-set simulations of soluble surfactant driven flows
CN105183965B (zh) 用于预测雾化过程的大涡模拟方法
JP2022022999A5 (zh)
Lee et al. Modeling and simulation of droplet evaporation using a modified Cahn–Hilliard equation
Sverdrup et al. An embedded boundary approach for efficient simulations of viscoplastic fluids in three dimensions
Rhea et al. RANS modelling and LES of a single-phase, impinging plane jet
Rahmati et al. Application of a modified pseudopotential lattice Boltzmann model for simulation of splashing phenomenon
Kaiser et al. A multiresolution local-timestepping scheme for particle–laden multiphase flow simulations using a level-set and point-particle approach
Fontes et al. Numerical simulation of a water droplet splash: Effects of density interpolation schemes
Tang et al. J1. 8 APPLICATION OF CFD SIMULATIONS FOR SHORT-RANGE ATMOSPHERIC DISPERSION OVER OPEN FIELDS AND WITHIN ARRAYS OF BUILDINGS
Brambilla et al. Assessment of dynamic adaptive grids in volume-of-fluid simulations of oblique drop impacts onto liquid films
Lu et al. A front tracking method for the simulation of compressible multimedium flows
Bernard et al. Simulation of particle dynamics for rarefied flows: Backflow in thruster plumes
Ong et al. An immersed boundary projection method for incompressible interface simulations in 3D flows
Subbareddy et al. A synthetic inflow generation method using the attached eddy hypothesis
Owkes et al. Large-eddy simulation study of injector geometry on liquid jet in cross-flow and validation with experiments
Bækkedal Alternative methods of realizing the sea spectrum for time-domain simulations of marine structures in irregular seas
de’Michieli Vitturi et al. An immersed boundary method for compressible multiphase flows: application to the dynamics of pyroclastic density currents
Liu et al. Simulation of incompressible multiphase flows with complex geometry using etching multiblock method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant