CN105117558B - 一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法 - Google Patents
一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法 Download PDFInfo
- Publication number
- CN105117558B CN105117558B CN201510580647.0A CN201510580647A CN105117558B CN 105117558 B CN105117558 B CN 105117558B CN 201510580647 A CN201510580647 A CN 201510580647A CN 105117558 B CN105117558 B CN 105117558B
- Authority
- CN
- China
- Prior art keywords
- liquid film
- liquid
- node
- spray
- thickness
- 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.)
- Expired - Fee Related
Links
- 239000007788 liquid Substances 0.000 title claims abstract description 129
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000001816 cooling Methods 0.000 title claims abstract description 17
- 238000005507 spraying Methods 0.000 claims abstract description 20
- 230000001133 acceleration Effects 0.000 claims abstract description 14
- 238000004088 simulation Methods 0.000 claims abstract description 10
- 239000007921 spray Substances 0.000 claims description 43
- 238000009825 accumulation Methods 0.000 claims description 13
- 230000002411 adverse Effects 0.000 claims description 10
- 230000005484 gravity Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims 1
- 238000012856 packing Methods 0.000 abstract 1
- 238000011160 research Methods 0.000 description 8
- 238000002474 experimental method Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 230000003068 static effect Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 239000012530 fluid Substances 0.000 description 2
- 230000003116 impacting effect Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000005486 microgravity Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
Abstract
本发明提供了一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法,其特征在于包括:过载加速度的模拟方法;喷雾状态模拟方法;不利高过载导致液膜堆积现象发生的判断方法;不利高过载导致液膜破碎现象发生的判断方法。
Description
技术领域
本发明涉及一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法。
背景技术
喷雾冷却作为一种高效的换热方式,在航空航天领域的大功率电子系统中有十分广阔的应用前景。电子设备在机载环境下将面临不同方向上的过载影响。以飞机为例,机动飞行过程中所受载荷比水平直线稳定飞行时大很多,而且飞机机动性要求越高,所承受的过载也越大。飞机过载系数是指飞机所受除重力之外的外力总和与飞机重量之比,战斗机在俯冲时的过载系数可以达到9~10。机动飞行动作不同则所受的过载力大小和方向也不相同。喷雾液滴撞击热源面后将在壁面上形成一层薄液膜,该液膜对喷雾冷却的换热性能有重要作用。过载力的存在将会引发表面液膜畸变,改变液膜在热源面上的分布和流动形态,进而影响喷雾冷却换热效果。对于过载系数较大的高过载环境,这一作用效果将更为严重,且有可能对换热造成不利影响。因此,针对喷雾冷却应用于航空航天领域的可行性,有必要研究不利高过载对喷雾冷却表面液膜畸变的影响。
目前研究重力因素对喷雾冷却性能的影响多以实验方法为主,常用的方法有两种。一是利用飞机的抛物飞行,结合调整热源面的朝向,为喷雾冷却创造超重和失重环境。二是利用落塔装置实现微重力环境的模拟。两种方法可以获得的实验时间很短,如飞机的抛物飞行实验约为20s,落塔实验只有几秒。相关实验可以直观地反映重力因素对喷雾冷却换热的影响,但由于实验时间很短且实验条件有差异,实验现象观察不充分,结论存在不一致性,不利于非常重力环境下喷雾冷却换热机理的研究。
另一方面,针对喷雾冷却换热机理已开展大量的数值模拟研究,多以单个或有限个液滴撞击壁面或表面液膜的流动过程为研究对象,观察不同初始条件对结果的影响。喷雾冷却是大量液滴源源不断撞击壁面后形成液膜的流动过程,有限个液滴撞击并不能很好反映出热源表面液膜的形成和流动形态。目前较为常用的计算模型为VOF多相流动模型。VOF模型多被用于计算常重力环境下的液滴撞击热源面问题,几乎没有被用于高过载环境。
VOF模型是一种多相流模型,通过求解单一的动量方程和追踪控制体内流体的体积分数来模拟两种或多种互不相融的流体的流动。当在模型中引入新的一相,则会加入该相的体积分数作为新的变量,所有相的体积分数之和为1。若各相的体积分数已知,则所有变量和属性由各相共享且以均值来体现。
质量守恒方程:
体积分数方程:
物质属性方程:ρ=∑αiρi
动量方程:
其中,为速度,α为体积分数,t为时间,S为质量源项,ρ为密度,m为质量传递,φ为混合属性,P为静态压力,μ为动力粘度,为动力源项。
发明内容
根据本发明的一个方面,提供了一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法,其特征在于包括:提出过载加速度的模拟方法;建立喷雾状态模拟方法;建立不利高过载导致液膜堆积现象发生的判断方法;建立不利高过载导致液膜破碎现象发生的判断方法。
根据本发明的另一个方面,提供了一种高过载环境下液膜畸变的确定方法,其特征在于:
提出过载加速度的模拟方法,在VOF多相流动模型中引入过载加速度,建立过载环境下的VOF计算模型。
建立喷雾状态模拟方法,将喷雾模型简化为多列液滴撞击热源面模型,以三个液滴参数控制工质流量和喷雾疏密程度;
建立不利过载导致液膜堆积现象发生的判断方法,通过对热源壁面节点对应的液膜厚度进行检验,判断液膜堆积现象发生的位置;
建立不利过载导致液膜破碎现象发生的判断方法,通过对热源壁面节点对应的液膜厚度进行检验,判断液膜断裂现象发生的位置。
说明书附图
图1显示了根据本发明的一个实施例的方法的组成步骤。
图2为根据本发明的一个实施例的模型实例。
图3为根据本发明的一个实施例的喷雾状态简化图。
图4为根据本发明的一个实施例的不利过载导致液膜堆积现象发生的判断流程。
具体实施方式
本发明提供了一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法,解决实验条件限制导致的实验时间过短及现象观察不充分等问题,对常重力条件下的VOF模型加以修正,建立过载条件下的计算模型,并以液膜堆积和破碎的判断方法为依据,确定不利过载引发液膜畸变现象发生的位置,研究过载影响。
根据本发明的一个方面,建立了过载条件下的数值计算模型;对喷雾状态进行模拟,其中对喷雾模型进行了简化,以三个液滴参数控制工质流量和喷雾疏密程度;对液膜畸变现象发生进行判断,确定在一定过载条件和喷雾状态下,液膜堆积和破裂现象发生的位置。
图1显示了根据本发明的一个实施例的一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法,包括:
(1)建立过载环境下的算法模型
以常重力环境下的VOF模型为基础,将参考系建立在水平加速运动的热源面上,引入平移惯性力,建立过载环境下的算法模型。确定计算区域,进行时间及空间上的网格划分,设置边界条件。(图1的步骤101)
(2)对喷雾状态进行模拟
将喷雾模型简化为多列液滴撞击壁面模型。以液滴直径、液滴速度和液滴间距为描述喷雾状态的三个重要参数,控制工质流量和喷雾疏密程度。(图1的步骤102)
(3)判断液膜畸变
根据计算的液膜厚度,判断液膜畸变,确定由不利高过载引发的液膜堆积和破裂现象发生的位置。(图1的步骤103)
下面结合二维实例介绍本发明研究液膜畸变方法的具体实施步骤。如图2所示,计算区域宽度为10mm,高度为5mm,空间步长Δx和Δy均为0.1mm,网格数量为5×103个。左右两端为压力出口,顶端为压力入口,底部为热源面壁面。边界大气压力为0.1MPa,温度为27℃。空气为主相,液态水为第二相,水的表面张力系数为0.073N/m,水与壁面的静态接触角为90°。
(1)过载环境下算法模型的建立
在地面惯性系中热源面以加速度a0运动,将坐标系建立在热源壁面上,则确定了一个非惯性系。在此非惯性系中引入惯性加速度a,其大小与a0相等,方向与a0相反。据此,液滴撞击以加速度a0运动的热源面问题,转换为了处于惯性力场a中的液滴撞击固定热源面问题。
过载环境下VOF模型中的动量方程变为如下形式
其中ρ为密度;
μ为动力粘度;
为速度;
t为时间;
P为静态压力;
为重力加速度;
为引入的惯性加速度;
为表面张力产生的动量源相。
(2)模拟喷雾状态
喷雾状态直接影响液膜在表面上的流动形态,与多个参数有关。根据本发明的一个实施例,将喷雾简化为多列液滴撞击热源面问题,对于某一喷雾状态,假设所有液滴直径、液滴间距和液滴速度固定,且液滴撞击热源面时的速度垂直于热源面,如图3所示。其中,喷雾过程中液滴源源不断以一定初始条件撞击热源面。把液滴直径、液滴速度和液滴间距作为描述喷雾状态的三个重要参数,通过三个参数控制热源面上的工质流量和喷雾疏密程度。把单列液滴对应的微小热源面作为研究对象,则可计算该微小热源面上的液滴数流量、体积流量和喷雾疏密系数。
液滴的数流量、体积流量和喷雾疏密系数由下式确定
其中,
d为液滴直径;
v为液滴速度;
l为液滴间距;
N为液滴数流量;
Q为液滴体积流量;
ξρ为喷雾疏密系数。
通过液滴直径、液滴速度和液滴间距三个参数的设置,可计算液滴数流量、体积流量和喷雾疏密系数这三个直接反映喷雾状态的复合参数,进而模拟多种喷雾状态,为不同喷雾状态在高过载环境下对液膜形态的影响研究做铺垫。如在本实例中可取d=0.5mm,v=7m/s,l=1.2mm,则N=5.83×103/s,Q=0.38×10-6m3/s,ξρ=2.4。
(3)液膜畸变判断
以建立的过载环境下的算法模型为基础,结合液滴动态模拟设定不同的喷雾形态,对液滴撞击热源面的两相流动过程进行模拟。不利高过载环境下的液膜畸变主要有液膜堆积和液膜破裂,本发明针对两现象提出了对应的判断方法。
1.不利过载导致液膜堆积现象发生的判断(见图4)
(a)根据建立的过载环境下的计算模型,计算得到k时刻各节点的体积分数。根据体积分数获得该时刻壁面节点所对应的液膜厚度
(b)根据k时刻壁面节点的液膜厚度值计算该时刻液膜的平均厚度:
其中k代表时刻,j代表壁面节点,n为节点个数,为j节点在k时刻的液膜厚度,为k时刻液膜平均厚度。
(c)比较k时刻壁面节点的液膜厚度是否满足下列条件:
且
若满足,则记录节点位置j,继续比较下一个节点j+1,直至搜寻完所有壁面节点n。
(d)重复步骤(a)(b),若某节点在某时间段(如r=20s)内一直满足步骤(c)中条件,则判定该点发生液膜堆积现象。
2.不利过载导致液膜断裂现象发生的判断
(a)根据建立的过载环境下的VOF模型,计算得到k时刻各节点的体积分数。根据体积分数获得该时刻壁面节点所对应的液膜厚度
(b)比较k时刻壁面节点的液膜厚度是否满足下列条件:
其中为j节点在k时刻的液膜厚度。
若满足,则记录节点位置j,判定该节点位置发生液膜断裂,继续比较下一节点j+1,直至搜寻完所有壁面节点。
根据以上方法可以找到热源面发生液膜畸变的位置和时刻。
本发明与现有的研究方法相比具有以下优点:
(1)建立了过载加速度的模拟方法,将液滴撞击加速运动的热源面问题,转变为处于惯性力场中的液滴撞击固定热源面问题。对常重力环境下的VOF模型加以修正,建立了过载环境下的计算模型,为不利高过载环境下的液膜畸变研究做好铺垫。
(2)提出喷雾状态模拟方法,将喷雾模型简化为多列液滴源源不断撞击热源面模型。以液滴直径、速度和间距作为描述喷雾状态的三个重要参数,通过三个参数控制喷雾流量和喷雾疏密程度。此方法将克服单个或有限个液滴撞击壁面研究的局限性,在更加长的时间段内观察和研究液膜的流动形态。
(3)本发明提出了不利高过载导致液膜畸变的判断方法,包括液膜堆积和破碎判断,可以更好地估计液膜畸变发生的位置,方便分析过载因素对液膜畸变的影响。
Claims (5)
1.不利高过载诱发的喷雾冷却表面液膜畸变的建模方法,其特征在于包括:
A)建立过载环境下的算法模型,
B)对喷雾状态进行模拟,
C)判断液膜畸变,包括:
根据计算的液膜厚度,判断液膜畸变,确定由不利高过载引发的液膜堆积现象和液膜破裂现象的发生和位置,
其中,
所述液膜堆积现象的发生和位置的确定包括:
计算获得某一时刻热源面各节点位置对应的液膜厚度,以其平均值作为该时刻热源面上液膜的平均厚度,
若某壁面节点在一定时间段内所对应的液膜厚度均高于该时间段内同一时刻热源面上液膜的平均厚度,且液膜厚度持续增加,则认为不利高过载导致该位置发生液膜堆积现象,包括:
i)根据k时刻各壁面节点的液膜厚度值计算该时刻液膜的平均厚度:其中k为时刻,j为节点位置,n为节点个数,为j节点在k时刻的液膜厚度,为k时刻液膜平均厚度,
ii)比较k时刻各壁面节点的液膜厚度是否满足下列条件:
且
若满足,则记录节点位置j,继续比较下一个节点j+1,直至搜寻完所有节点,
iii)k=k+1,重复步骤i)和ii),若某壁面节点在一个预定时间段内一直满足步骤ii)中条件,判定该点发生液膜堆积现象,
所述液膜破裂现象的发生和位置的确定包括:
确定某一时刻热源面各节点位置对应的液膜厚度,若k时刻节点j对应的液膜厚度等于或接近于零,则判定不利高 过载导致该位置在该时刻发生液膜破裂现象,
若其中为j节点在k时刻的液膜厚度,则节点位置j在k时候发生液膜破裂现象。
2.根据权利要求1的不利高过载诱发的喷雾冷却表面液膜畸变的建模方法,其特征在于所述步骤A)包括:
以常重力环境下的VOF模型为基础,将参考系建立在水平加速运动的热源面上,引入惯性力从而建立起非惯性系,建立过载环境下的算法模型;
确定计算区域,进行时间及空间上的网格划分,设置边界条件。
3.根据权利要求1的不利高过载诱发的喷雾冷却表面液膜畸变的建模方法,其特征在于所述步骤B)包括:
将喷雾模型简化为多列液滴撞击壁面模型;
以液滴直径、液滴速度和液滴间距为描述喷雾状态的三个重要参数,控制工质流量和喷雾疏密程度。
4.根据权利要求2的不利高过载诱发的喷雾冷却表面液膜畸变的建模方法,其特征在于:
在所述非惯性系中,引入的所述惯性力的加速度大小和热源面运动的加速度大小相等,方向相反,
F+f=ma
f=-ma0
其中,
F为质点所受真实外力;
f为引入的平移惯性力;
m为质点质量;
a为非惯性系中质点运动加速度;
a0为惯性系中热源面运动加速度。
5.根据权利要求3的不利高过载诱发的喷雾冷却表面液膜畸变的建模方法,其特征在于:
对于某一喷雾状态,液滴直径、速度和间距固定,且液滴撞击速度垂直于热源面,
通过液滴直径、液滴速度和液滴间距作为描述喷雾状态三个参数控制工质流量和喷雾疏密程度:
其中,
d为液滴直径;
v为液滴速度;
l为液滴间距;
N为液滴数流量;
Q为液滴体积流量;
ξρ为喷雾疏密系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510580647.0A CN105117558B (zh) | 2015-09-13 | 2015-09-13 | 一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510580647.0A CN105117558B (zh) | 2015-09-13 | 2015-09-13 | 一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105117558A CN105117558A (zh) | 2015-12-02 |
CN105117558B true CN105117558B (zh) | 2018-05-04 |
Family
ID=54665546
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510580647.0A Expired - Fee Related CN105117558B (zh) | 2015-09-13 | 2015-09-13 | 一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105117558B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020045A (zh) * | 2014-05-20 | 2014-09-03 | 北京航空航天大学 | 一种测定飞机/直升机动部件高周疲劳特性的步进加速试验方法 |
CN104354867A (zh) * | 2014-09-26 | 2015-02-18 | 空气动力学国家重点实验室 | 一种过冷大水滴结冰探测器的设计方法及探测器 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8466649B2 (en) * | 2010-05-19 | 2013-06-18 | The Invention Science Fund I Llc | Heat removal from motor components |
-
2015
- 2015-09-13 CN CN201510580647.0A patent/CN105117558B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020045A (zh) * | 2014-05-20 | 2014-09-03 | 北京航空航天大学 | 一种测定飞机/直升机动部件高周疲劳特性的步进加速试验方法 |
CN104354867A (zh) * | 2014-09-26 | 2015-02-18 | 空气动力学国家重点实验室 | 一种过冷大水滴结冰探测器的设计方法及探测器 |
Non-Patent Citations (1)
Title |
---|
多喷嘴喷雾冷却实验研究与数值模拟;侯燕;《中国博士学位论文全文数据库 工程科技II辑》;20141015;摘要、第1-98页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105117558A (zh) | 2015-12-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shinan et al. | Numerical investigation on impingement dynamics and freezing performance of micrometer-sized water droplet on dry flat surface in supercooled environment | |
Mostaghimi et al. | Modeling thermal spray coating processes: a powerful tool in design and optimization | |
Wu et al. | Time‐dependent turbulent cavitating flow computations with interfacial transport and filter‐based models | |
Wang et al. | Asymmetric heat transfer characteristics of a double droplet impact on a moving liquid film | |
Wang | Numerical study of a droplet impact on cylindrical objects: Towards the anti-icing property of power transmission lines | |
Raman et al. | Dynamics of simultaneously impinging drops on a dry surface: Role of impact velocity and air inertia | |
Li et al. | Incompressible multiphase flow and encapsulation simulations using the moment‐of‐fluid method | |
Lane et al. | Lagrangian trajectory modeling of lunar dust particles | |
Hu et al. | 3D simulation of micro droplet impact on the structured superhydrophobic surface | |
KR20140115835A (ko) | 공기 흡입구의 표면 결빙 시뮬레이션 방법 | |
Balakrishnan et al. | Fluid density effects in supersonic jet-induced cratering in a granular bed on a planetary body having an atmosphere in the continuum regime | |
Vahab et al. | An adaptive coupled level set and moment-of-fluid method for simulating droplet impact and solidification on solid surfaces with application to aircraft icing | |
Wang et al. | Numerical investigation of Weber number and gravity effects on fluid flow and heat transfer of successive droplets impacting liquid film | |
CN105117558B (zh) | 一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法 | |
Lian et al. | A modified spongy icing model considering the effect of droplets retention on the ice accretion process | |
Quero et al. | Analysis of super-cooled water droplet impact on a thin water layer and ice growth | |
Szilder et al. | Progress towards a 3D numerical simulation of ice accretion on a swept wing using the morphogenetic approach | |
GB2572352A (en) | Methods and apparatus for simulating liquid collection on aerodynamic components | |
Lian et al. | Experimental and numerical investigation of icing process of a liquid droplet | |
Yue et al. | Study of droplets distribution on canopy of ringsail parachute in light rain | |
Dastan et al. | Microfiber motion and web formation in a microchannel heat sink: A numerical approach | |
Szilder et al. | Numerical simulation of complex ice shapes on swept wings | |
Shimura et al. | Comparative Investigation on Effect of Droplet Deformation Models on SLD Icing | |
Blake et al. | Effects of surface characteristics and droplet diameter on the freezing of supercooled water droplets impacting a cooled substrate | |
Mawarsih et al. | Numerical simulation of single droplet phenomenon using method finite difference and front-tracking |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180504 Termination date: 20190913 |