CN111125922A - 一种冷空间内下落液滴冻结过程的计算方法 - Google Patents
一种冷空间内下落液滴冻结过程的计算方法 Download PDFInfo
- Publication number
- CN111125922A CN111125922A CN201911384259.XA CN201911384259A CN111125922A CN 111125922 A CN111125922 A CN 111125922A CN 201911384259 A CN201911384259 A CN 201911384259A CN 111125922 A CN111125922 A CN 111125922A
- Authority
- CN
- China
- Prior art keywords
- temperature
- distribution function
- calculating
- freezing
- liquid drops
- 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
- 239000007788 liquid Substances 0.000 title claims abstract description 100
- 238000000034 method Methods 0.000 title claims abstract description 74
- 230000008014 freezing Effects 0.000 title claims abstract description 68
- 238000007710 freezing Methods 0.000 title claims abstract description 68
- 230000008569 process Effects 0.000 title claims abstract description 51
- 239000007790 solid phase Substances 0.000 claims abstract description 28
- 230000005855 radiation Effects 0.000 claims abstract description 18
- 238000005315 distribution function Methods 0.000 claims abstract description 16
- 239000012530 fluid Substances 0.000 claims abstract description 9
- 230000005012 migration Effects 0.000 claims abstract description 4
- 238000013508 migration Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 18
- 239000002245 particle Substances 0.000 claims description 10
- 239000007787 solid Substances 0.000 claims description 6
- 230000003993 interaction Effects 0.000 claims description 4
- 238000012821 model calculation Methods 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 230000009878 intermolecular interaction Effects 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 239000012071 phase Substances 0.000 abstract description 25
- 230000008859 change Effects 0.000 abstract description 23
- 238000009826 distribution Methods 0.000 abstract description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 9
- 238000004088 simulation Methods 0.000 description 8
- 230000007423 decrease Effects 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 238000001816 cooling Methods 0.000 description 4
- 239000007791 liquid phase Substances 0.000 description 4
- 230000005501 phase interface Effects 0.000 description 4
- 230000007704 transition Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 206010060904 Freezing phenomenon Diseases 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 238000004378 air conditioning Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002425 crystallisation Methods 0.000 description 1
- 230000008025 crystallization Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005057 refrigeration Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000007921 spray Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种冷空间内下落液滴冻结过程的计算方法,在焓法格子Boltzmann相变模型的基础上进行修正,构建了冷空间内下落液滴冻结过程计算方法。该方法使用密度分布函数和温度分布函数分别模拟流体的速度场和温度场,通过求解固相体积分数追踪相变过程中各相的体积分布情况,使用固相体积分数修正焓值以及流场演化方程迁移步,并考虑相变潜热和温差相对较大(大40K)时辐射换热对相变过程的影响。
Description
技术领域
本发明属于液滴冻结过程计算领域,尤其是涉及一种冷空间内下落液滴冻结过程的计算方法。
背景技术
下落液滴冻结现象普遍存在于自然界、制冷、空气调节、航空航天等领域。因此,掌握具体冻结过程对于工程应用来说具有十分重要的研究意义。以喷雾结晶和人工降雪为研究背景,从传热传质以及相变理论出发分析液滴冻结过程。当液滴温度被冷却至一定温度以下时,液滴会发生相变,即由液态转变为固态,该过程称为液滴冻结过程。液滴由液态转变为固态时的转变温度称为液滴冻结温度。在湿空气近冷壁面结霜机理中,液滴冻结过程是整个结霜过程的开始,具体可以概括为:液滴在冷环境中被冷却到冻结温度以下,液滴表面先形成冰壳,液滴由外向内逐渐发生冻结现象,最终整个液滴被冻结成冰粒。当液滴在冷空间内下落时,液滴被周围环境冷却发生冻结现象,液滴表面先形成冰壳,假定该冰壳为刚性,冰壳的存在将阻碍流体流动,但是换热过程依然继续。当液滴与周围环境温差相对较大(大于40K)时,辐射换热对相变过程有一定影响。因此,运动液滴冻结过程将同时受热传导、热对流以及热辐射影响。此外,液滴内部固液相界面处也将释放相变潜热。然而,由于液滴内部固液相界面不规则变化、固液两相区结构复杂以及温度场影响,液滴冻结过程的试验研究是一项极具挑战性和高成本的任务。虽然试验方法可以得到液滴冻结过程的温度分布,但是试验设备造成的误差或干扰对实验结果有很大影响,例如,当利用热电偶测量体积较小的液滴表面或内部温度时,热电偶的热阻会影响液滴的换热过程,因此热电偶对实验结果的准确性将产生一定影响。此外,观测用来判断液滴是否完全冻结的固液相界面位置也十分困难。因此,需要借助其他手段进行辅助研究。目前用来计算非等温冻结过程的方法主要是Jiaung提出的焓法格子Boltzmann模型。对于运动液滴冻结过程而言,该过程是一个涉及流体流动的非等温相变过程,目前用来解决非等温相变过程的数学模型大多未将冻结过程对流体流动的影响考虑在内。此外,当液滴与周围环境间的温差相对较大(大于40K)时,辐射换热对冻结过程也会产生一定影响,但大部分计算相变过程的格子Boltzmann模型都未考虑辐射换热。在目前的研究成果中,关于冷空间内下落液滴冻结过程中温度场、速度场以及液滴冻结状态的研究也十分有限。因此需要建立新的数学模型来解决这一问题。
发明内容
为了更好地解决冷空间内悬浮液滴以及下落液滴冻结问题,得到不同时刻液滴内温度场、液滴冻结状态、冻结参数及影响因素,本发明给出了一种冷空间内下落液滴冻结过程的计算方法。
针对上述问题,本发明给出了一种冷空间内下落液滴冻结过程的计算方法,在焓法格子Boltzmann相变模型的基础上进行修正,构建了冷空间内下落液滴冻结过程计算方法。该方法使用密度分布函数和温度分布函数分别模拟流体的速度场和温度场,通过求解固相体积分数追踪相变过程中各相的体积分布情况,使用固相体积分数修正焓值以及流场演化方程迁移步,并考虑相变潜热和温差相对较大(大于40K)时辐射换热对相变过程的影响。
为解决上述技术问题,本发明创造采用的技术方案是:
一种冷空间内下落液滴冻结过程的计算方法,主要包括如下步骤:
(步骤一)根据物理模型计算所需的计算域和网格尺寸,进行网格划分,如图1所示;
(步骤二)对模型计算域的密度、温度、速度以及固相体积分数进行初始化;
(步骤三)输出计算域的流场、温度场、辐射强度以及辐射热流量;
模型的流场演化方程为:
其中,f为粒子分布函数,τ为松弛时间。
其中,cs为格子声速,ωi为权系数,ρ为密度,ei为离散速度。平衡态速度ueq为:
分子间相互作用力为:
F(x)=Ff(x)+Fg(x)
其中,流体间相互作用力为:
式中,ψ(ρ)=ρ*[1-exp(-ρ/ρ*)]为有效密度。G1和G2分别为相邻和次相邻位置的流体间相互作用系数。
重力为:
Fg(x)=ρ(x)g
演化方程迁移步为:
fi(x,t+δt)=αsfi(x,t)+(1-αs)fi(x-eiδt,t)
其中,fi为碰撞后的粒子分布函数。
温度场演化方程为:
计算辐射强度的演化方程为:
(步骤四)输出计算域内每个节点的密度、速度以及温度;
密度、速度以及温度的计算式为:
(步骤五)输出计算域内每个节点的固相体积分数以及焓值;
固相体积分数αs表达式为:
焓值表达式为:
H=αscsolidT+(1-αs)cgasT+αsL
(步骤六)返回(步骤三)循环计算,直到液滴内部固相体积分数全为1,程序结束。
本发明具有以下优点和有益效果:
本发明所述的一种冷空间内下落液滴冻结过程的计算方法,可以解决冷空间内悬浮液滴以及下落液滴冻结问题,得到不同时刻液滴内温度场、液滴冻结状态、冻结参数及影响因素,使得液滴冻结过程的研究更加完善,并得到如下结论:液滴下部空气冷却效果较强,液滴下表面优先发生冻结现象;垂直初速度对液滴形变有一定影响,随着下落速度的增大,形成的冰粒形状从圆形变为椭圆形;液滴冻结时间随液滴初始温度降低而缩短,随下落速度增大而缩短;上述结论可以用于指导实际生产,即提高冷空间内下落液滴冻结效率的最有效方法是改变其初始温度。
附图说明
图1体现了D2Q9模型的网格划分和速度矢量;
其中,图1(a)为D2Q9模型的网格划分图;图1(b)为D2Q9模型的速度矢量图。
图2为焓法格子Boltzmann相变模型与下落液滴冻结过程计算方法模拟得到的液滴冻结状态的对比图;
其中,图2(a)为焓法格子Boltzmann相变模型;图2(b)为下落液滴冻结过程计算方法。
图3为冷空间内下落液滴冻结过程的物理模型。
图4给出了不同网格体系中下落液滴内A点的固相体积分数。
图5体现了不同时刻冷空间的温度场和速度场;
其中,图5(a)t=0s;图5(b)t=0.6s;图5(c)t=0.75s;图5(d)t=0.9s;图5(e)t=1.2s;
图5(f)t=1.5s;图5(g)t=1.8s;图5(h)t=2.22s;图5(i)为温度云图图例。
图6体现了不同时刻冷空间内下落液滴的固相体积分数;
图6(a)t=0s;图6(b)t=0.6s;图6(c)t=0.75s;图6(d)t=0.9s;图6(e)t=1.2s;图6(f)t=1.5s;图6(g)t=1.8s;图6(h)t=2.22s;图6(i)为固相体积分数云图图例。
图7为t=0.75s和t=1.2s时下落液滴内O点到B点的固相体积分数对比图。
图8体现了相同时间间隔下落液滴内O点和B点的固相体积分数增量对比图。
图9为t=0.9s时不同垂直初速度的下落液滴冻结形状对比图;
其中,图9(a)v0=0.001m/s;图9(b)v0=0.005m/s;图9(c)v0=0.01m/s;图9(d)v0=0.015m/s;图9(e)t=0.02m/s。
图10为不同垂直初速度的下落液滴长短轴比对比图。
图11为不同液滴直径时液滴冻结时间与液滴初始温度的关系对比图。
具体实施方式
下面通过实施例对本发明作进一步详细说明。
实施例1
本实施例分别采用焓法格子Boltzmann模型和下落液滴冻结过程计算方法模拟冷空间内下落液滴冻结过程。选取计算域为0.03m×0.07m,计算网格数为600×1400,四周采用周期性边界。冷空间内有一初始温度为Tl,in=273K的圆形液滴,液滴初始位置为(0.015m,0.06m),液滴初始直径为D=0.003m,液滴内O点为中心点,B点为边界点,A点为O点和B点的中点,液滴垂直初速度为v0=0.001m/s,水密度为ρl=1000kg/m3,液滴初始固相体积分数为αs=0。冷空间内空气初始温度为Tg,in=223K,空气密度为ρg=1.584kg/m3。选取水的冻结温度为Tf=268K,水比热容为cliquid=4.2×103J/(kg·K),冰比热容为csolid=2.01×103J/(kg·K),水导热率为k=0.55W/(m·K),空气Prandtl数为Pr=0.71,Stefan数为Ste=5.625。如图2所示,由对比结果可知,使用焓法格子Boltzmann相变模型模拟得到的下落液滴在表面形成冰壳后依然发生形变,这与实际下落液滴冻结过程不符,进一步说明了下落液滴冻结过程计算方法的正确性以及优势。
实施例2
利用本发明的冷空间内下落液滴冻结过程计算方法模拟冷空间内下落液滴冻结过程。图3给出了相应的物理模型,选取计算域为0.03m×0.07m,计算网格数为600×1400,四周采用周期性边界。冷空间内有一初始温度为Tl,in=273K的圆形液滴,液滴初始位置为(0.015m,0.06m),液滴初始直径为D=0.003m,液滴内O点为中心点,B点为边界点,A点为O点和B点的中点,液滴垂直初速度为v0=0.001m/s,水密度为ρl=1000kg/m3,液滴初始固相体积分数为αs=0。冷空间内空气初始温度为Tg,in=223K,空气密度为ρg=1.584kg/m3。选取水的冻结温度为Tf=268K,水比热容为cliquid=4.2×103J/(kg·K),冰比热容为csolid=2.01×103J/(kg·K),水导热率为k=0.55W/(m·K),空气Prandtl数为Pr=0.71,Stefan数为Ste=5.625。
为了确保网格密度对计算精度没有影响,使用4套均匀网格分别进行数值模拟计算(网格1:150×350,网格2:300×700,网格3:600×1400,网格4:1200×2800)。图4给出了不同网格体系中下落液滴内A点的固相体积分数。模拟结果显示,网格3和网格4得到的模拟结果之间相对误差小于0.1%。考虑计算成本,采用网格3进行数值模拟计算。
图5和图6分别给出了不同时刻冷空间的温度场和速度场以及下落液滴的固相体积分数。冷空间内下落液滴冻结过程可以概括如下:液滴在下落过程中被周围环境冷却,液滴温度随时间不断降低,当液滴温度低于水的冻结温度时,液滴发生冻结现象。对于整个液滴来说,冻结现象首先发生在液滴表面,随着时间的推移,液滴内固液相界面位置逐渐向液滴中心移动,直到液滴完全冻结成冰粒。由模拟结果可知,液滴内部固相体积分数在下落过程后期呈不对称分布,液滴下部先于上部发生冻结现象,主要因为液滴下部空气冷却效果强于上部。
图7给出了t=0.75s和t=1.2s时下落液滴内O点(液滴中心)到B点(液滴边缘)的固相体积分数。模拟结果显示,同一时刻液滴表面的固相体积分数高于液滴中心,主要因为热量的传导由外向内进行,液滴表面先发生冻结现象。图8给出了相同时间间隔下落液滴内O点和B点的固相体积分数增量。初始时刻,液滴未开始冻结,液滴内O点和B点的固相体积分数增量为零。一段时间后,液滴开始冻结,此时液滴内O点和B点的固相体积分数增量先增大后减小。在液滴冻结过程后期,液滴冻结趋势逐渐减缓,主要因为相变驱动力在液滴冻结过程后期逐渐减小。
当下落液滴发生冻结现象时,液滴表面形成冰壳,冰壳的存在将阻碍下落液滴发生形变,随后液滴继续下落但不再发生形变。图9显示了t=0.9s时不同垂直初速度的下落液滴冻结形状。模拟结果显示,随着液滴垂直初速度的增大,液滴形状从圆形变为椭圆形。图10给出了不同垂直初速度(从v0=0.001m/s到v0=0.02m/s)的下落液滴长短轴比。由模拟结果可知,液滴在垂直初速度较大时形变比较明显,因为垂直初速度越大,液滴下落速度越大,液滴下落过程中受到的阻力越大,液滴形变越明显。
冻结时间是决定冻结效率的重要因素之一,通过求解固相体积分数来预测液滴冻结时间。图11给出了不同液滴直径时液滴冻结时间与液滴初始温度的关系。液滴冻结时间随液滴初始温度降低而缩短,主要原因为相变驱动力随液滴初始温度降低而增大,较高的相变驱动力将促进相变现象发生,因此缩短了液滴冻结时间,提高了液滴冻结效率。
通过模拟冷空间内不同下落方式液滴冻结过程,获得了不同时刻冷空间的温度场、速度场以及下落液滴冻结状态及其相应的冻结参数,考察了下落方式和液滴初始温度对液滴冻结现象的影响,得到如下结论:液滴下部空气冷却效果较强,液滴下表面优先发生冻结现象;垂直初速度对液滴形变有一定影响,随着下落速度的增大(从v0=0.001m/s到v0=0.02m/s),形成的冰粒形状从圆形变为椭圆形;液滴冻结时间随液滴初始温度降低而缩短,随下落速度增大而缩短。
上述结论可以用于指导实际生产,即可以通过延长冷却时间来降低液滴初始温度(从269K到278K)。此外,同一液滴初始温度时,减小液滴直径也可以缩短液滴冻结时间(以275K为例,将直径由3mm减小到1mm时,可以缩短2秒左右的冻结时间),主要因为减小液滴直径将导致液滴单位重量表面积增大以及液滴终态速度降低,进而缩短液滴冻结时间。然而,在实际应用中(如人工造雪),如果液滴直径太小,液滴会因终态速度过小而被风吹走。因此,提高冷空间内下落液滴冻结效率的最有效方法是改变其初始温度。
以上对本发明创造的一个实施例进行了详细说明,但所述内容仅为本发明创造的较佳实施例,不能被认为用于限定本发明创造的实施范围。凡依本发明创造申请范围所作的均等变化与改进等,均应仍归属于本发明创造的专利涵盖范围之内。
Claims (1)
1.一种冷空间内下落液滴冻结过程的计算方法,其特征在于,主要包括如下步骤:
(步骤一)根据物理模型计算所需的计算域和网格尺寸,并进行网格划分;
(步骤二)对模型计算域的密度、温度、速度以及固相体积分数进行初始化;
(步骤三)输出计算域的流场、温度场、辐射强度以及辐射热流量;
模型的流场演化方程为:
其中,f为粒子分布函数,τ为松弛时间;
平衡态分布函数fi (eq)(x,t)为:
其中,cs为格子声速,ωi为权系数,ρ为密度,ei为离散速度;平衡态速度ueq为:
分子间相互作用力为:
F(x)=Ff(x)+Fg(x)
其中,流体间相互作用力为:
式中,ψ(ρ)=ρ*[1-exp(-ρ/ρ*)]为有效密度;G1和G2分别为相邻和次相邻位置的流体间相互作用系数;
重力为:
Fg(x)=ρ(x)g
演化方程迁移步为:
fi(x,t+δt)=αsfi(x,t)+(1-αs)fi(x-eiδt,t)
其中,fi为碰撞后的粒子分布函数;
温度场演化方程为:
计算辐射强度的演化方程为:
(步骤四)输出计算域内每个节点的密度、速度以及温度;
密度、速度以及温度的计算式为:
(步骤五)输出计算域内每个节点的固相体积分数以及焓值;
固相体积分数αs表达式为:
焓值表达式为:
H=αscsolidT+(1-αs)cgasT+αsL
(步骤六)返回(步骤三)循环计算,直到液滴内部固相体积分数全为1,程序结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911384259.XA CN111125922B (zh) | 2019-12-28 | 2019-12-28 | 一种冷空间内下落液滴冻结过程的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911384259.XA CN111125922B (zh) | 2019-12-28 | 2019-12-28 | 一种冷空间内下落液滴冻结过程的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111125922A true CN111125922A (zh) | 2020-05-08 |
CN111125922B CN111125922B (zh) | 2024-02-20 |
Family
ID=70505721
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911384259.XA Active CN111125922B (zh) | 2019-12-28 | 2019-12-28 | 一种冷空间内下落液滴冻结过程的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111125922B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140257771A1 (en) * | 2011-11-30 | 2014-09-11 | Ming Lu | Numerical simulation method for aircrasft flight-icing |
DE102015217177A1 (de) * | 2015-09-09 | 2017-03-09 | Bayerische Motoren Werke Aktiengesellschaft | Transiente TBS-Berechnung des Gesamtfahrzeugs, Standardlastfälle und Nachheizen |
CN109766587A (zh) * | 2018-12-17 | 2019-05-17 | 西南科技大学 | 一种大密度比气液相变的混合lbm数值模型及应用 |
CN110390176A (zh) * | 2019-07-31 | 2019-10-29 | 西南交通大学 | 一种无砟轨道冻结与损伤行为计算方法 |
-
2019
- 2019-12-28 CN CN201911384259.XA patent/CN111125922B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140257771A1 (en) * | 2011-11-30 | 2014-09-11 | Ming Lu | Numerical simulation method for aircrasft flight-icing |
DE102015217177A1 (de) * | 2015-09-09 | 2017-03-09 | Bayerische Motoren Werke Aktiengesellschaft | Transiente TBS-Berechnung des Gesamtfahrzeugs, Standardlastfälle und Nachheizen |
CN109766587A (zh) * | 2018-12-17 | 2019-05-17 | 西南科技大学 | 一种大密度比气液相变的混合lbm数值模型及应用 |
CN110390176A (zh) * | 2019-07-31 | 2019-10-29 | 西南交通大学 | 一种无砟轨道冻结与损伤行为计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111125922B (zh) | 2024-02-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108304601B (zh) | 一种高超声速飞行器边界层转捩的判断方法 | |
Nahor et al. | CFD model of the airflow, heat and mass transfer in cool stores | |
CN108304600B (zh) | 一种高超声速飞行器转捩位置预测方法 | |
Hindmarsh et al. | Experimental and numerical analysis of the temperature transition of a suspended freezing water droplet | |
Zhang et al. | Development of a Para-AMR algorithm for simulating dendrite growth under convection using a phase-field–lattice Boltzmann method | |
Sun et al. | Dynamic study of frost formation on cryogenic surface | |
CA2430823C (en) | Morphogenetic modelling of in-flight icing | |
Norde et al. | Eulerian method for ice crystal icing | |
CN111062121A (zh) | 基于高度函数-格子玻尔兹曼法的粉末熔化数值模拟方法 | |
CN105426577B (zh) | 一种大型原油浮顶储罐非稳态传热过程数值模拟方法 | |
Gong et al. | Mesoscopic investigation of frost crystal nucleation on cold surface based on the lattice-Boltzmann method | |
CN111125922A (zh) | 一种冷空间内下落液滴冻结过程的计算方法 | |
Liu et al. | Spreading and freezing of supercooled water droplets impacting an ice surface | |
Gong et al. | A lattice Boltzmann study of frost growth on a cold surface | |
CN111159890B (zh) | 一种用于抑制预冷器结霜的模拟计算方法 | |
Bodaghkhani et al. | Water droplet freezing on cold surfaces with distinct wetabilities | |
CN111144012B (zh) | 一种用于冷空间内冰粒沉积过程的计算方法 | |
Lian et al. | Experimental and numerical investigation of icing process of a liquid droplet | |
Zhao et al. | Analysis of freezing process about falling droplet using the lattice Boltzmann method | |
Berdnikov et al. | Numerical simulation of crystal growth processes by means of horizontal unidirectional crystallization from melts with different Prandtl numbers | |
Szilder et al. | Three‐dimensional modelling of ice accretion density | |
Cui et al. | An Improved Enthalpy-Based Lattice Boltzmann 3D Model with Added Probability Function for the Growth Characteristics of Frost Layer | |
Abbadi et al. | Preliminary analysis of ice accretion prediction on wind turbine blades | |
CN111159942A (zh) | 一种基于定常模拟的带翼飞行器滚转阻尼力矩计算方法 | |
Naterer | Dispersed multiphase flow with air-driven runback of a liquid layer at a moving boundary |
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 |