高温辐射环境中液滴的沸腾效应

建立了液滴在高温对流和辐射环境中的受热和蒸发模型,结合液滴均质沸腾模型,编制了计算程序。以正十二烷液滴为例,考虑液滴的膨胀效应以及液滴与周围气流的热物性变化,数值模拟了高温辐射与对流加热下的液滴升温和蒸发过程。分析了不同对流和高温辐射条件下,液滴内部是否能够发生沸腾。研究表明,液滴在高温辐射和对流加热下,蒸发伴随热膨胀;高温热辐射加热可导致液滴内部温度高于表面温度,升温到一定程度后可达到液滴内部沸腾状态;影响液滴沸腾的因素有液滴半径、辐射温度、环境气流温度等;同时,随着液滴蒸发,高温环境中液滴的沸腾过热度逐渐增大。     

单液滴蒸发影响因素实验研究

设计了液滴蒸发实验装置,利用高速摄像系统记录高温气流中液滴蒸发过程,并根据实验的数据总结出液滴蒸发规律。在蒸发过程中,液滴蒸发首先经历一个非稳态初始加热阶段,然后液滴蒸发进入一个相对稳定的蒸发阶段,其直径变化基本遵循D2定律。气流温度越高,气流速度越大,液滴蒸发时间越短,液滴蒸发速度越快,但是气流速度对蒸发影响很小。随着液滴初始直径的增大,内部环流加强,液滴后部的尾涡加大,传质阻力减小,液滴的蒸发常数变大。     

盐水液滴降压蒸发析盐过程传热传质特性

针对单个盐水(NaCl溶液)液滴在降压环境下蒸发析盐的传热传质过程建立了数学模型。模型考虑了多孔盐壳在液滴表面的形成过程,降压过程引起的气流运动,液核通过多孔介质的传质扩散,以及液滴表面的蒸发换热和对流换热。将实验数据与计算结果对比,验证了模型的有效性。通过模型计算获得了液滴表面温度及液滴质量随时间的变化。结果表明盐水液滴在降压环境下蒸发析盐过程的温度变化分为4个阶段:温度骤降阶段、温度回升阶段、平衡温度阶段和温度上升阶段。平衡温度阶段,盐壳界面运动较慢,随蒸发进行,液核尺寸逐渐减小,盐壳界面运动速度加快。理论分析了环境压力对盐水液滴蒸发析盐过程的影响,环境压力越低,平衡温度越低,盐分完全析出时间越短。     

多物理场耦合模拟微波蒸馏反应器：升温和沸腾过程

使用COMSOL Multiphysics软件建立了耦合电磁场、流体流动、传热以及物质传输的多物理场模型用于模拟蒸馏型反应器的微波能量利用过程,探究了蒸馏反应器中水负载在微波能辐射作用下从升温至沸腾过程,阐明了在升温阶段,样品温度呈上下层分布,上层温度较高,最大温差达20 K,自然对流的产生改善了温度分布的不均匀性;在沸腾阶段,由于下层温度较低,沸腾现象有延迟,气泡的产生消除了部分过热,其中表面蒸发量更大,最大时约为内部蒸发量的3倍,与此同时湍流现象明显改善了温度均匀性。探究了馈入功率对全沸腾状态的影响,揭示了全沸腾状态的最终温度取决于馈入功率和蒸发损耗功率的相对大小。研究结果可为微波辅助分离、反应等化工过程及装备设计提供理论基础与借鉴。     

尿素水溶液液滴的蒸发特性

在石英管式炉上通过挂滴法来观察单个尿素水溶液（urea-aqueous-solution,UAS）液滴的具体蒸发过程,比较了不同环境温度以及不同初始直径大小下液滴的蒸发特性。结果表明,尿素溶液液滴在100～1300 ℃的温度范围内呈现出了不同的蒸发行为。在较高的温度下,液滴的蒸发行为较为复杂,如气泡的产生、液滴的变形以及发生微爆的现象;但是,随着环境温度的降低,这些现象就变得非常微弱甚至消失。同时,还定量分析了稳态蒸发常数与温度、液滴初始直径之间的变化关系,发现在初始直径为2.5 mm、温度在100～600 ℃之间变化的情况下,稳态蒸发常数从0.02075 mm2/s增加到了0.23953 mm2/s,增大了10倍左右。此外,还对气流流速为0.25～1.25 m/s范围内的液滴蒸发特性作了实验研究。当液滴周围有强迫气流存在时,液滴与气体间的换热方式由导热转变为对流换热,从而增强了液滴表面的传热传质能力,促进了液滴的蒸发。     

乙醇溶液液滴降压蒸发过程传热传质特性

针对单个乙醇溶液液滴在降压环境下蒸发的传热传质过程建立了数学模型。模型基于液相的能量守恒和 传质扩散理论,利用经典拓展模型计算液滴的质量蒸发率,并引入活度系数考虑液滴表面的蒸气分压。采用液 滴悬挂法进行实验,分别记录了乙醇溶液液滴和乙酸溶液液滴在降压蒸发过程中的液滴内温度变化。将实验数 据与计算结果对比,验证了模型的有效性。通过模型计算获得了液滴内部温度分布以及浓度分布随时间的变化。 结果表明:快速降压阶段空气流动较快,加之乙醇工质易挥发,液滴表面温度下降迅速,液滴内部温差和乙醇 浓度梯度较大;压力稳定后,空气流速为零,液滴内部温差和乙醇浓度梯度逐渐减小。由于液滴内部的热扩散 速率大于传质扩散系数,内部温度随时间的变化比浓度随时间的变化更快。     

局部微型加热下液滴表面温度分布特性

液滴表面温度分布直接影响液滴内部微流状况,而目前文献对液滴表面温度研究主要基于平板全局加热模式。采用MEMS集成工艺和红外热像分析手段,对基于中心局部微型加热下的液滴表面温度分布特性进行了实验研究。研究发现：局部加热液滴表面温度分布与平板全局加热液滴表面温度分布不同,呈现顶端温度高、边缘温度低的凸形温度分布规律;随着加热功率增加,液滴表面温度和温度梯度都会随之增加,而当加热功率增加到一定值后,液滴表面温度增幅趋于一致,表面温度梯度趋于稳定分布状态。同时对液滴局部沸腾时气泡破裂前后液滴表面温度分布进行了研究。研究结果有助于理解和控制液滴微流分布。     

低压环境下海水液滴蒸发模型研究

针对低压环境下的多种盐分组成的海水液滴(NaCl、MgCl_2、CaCl_2、KCl的质量分数分别为20.07%、1.29%、0.412%、0.399%)的蒸发过程进行了模型研究,模型将能量守恒方程和传质扩散方程进行耦合。通过模型计算,分析了液滴内组分含量和温度分布随时间的变化,获得了环境压力和环境温度对液滴温度变化的影响。结果表明,在蒸发过程中,液滴内部的各种盐含量升高过程不是成单一的线性增加的,在液滴的蒸发和内部盐含量扩散的共同作用下,盐含量表现出在先快速上升之后,逐渐下降继而又有所回升的趋势。液滴温度下降至最低温度后,也逐渐回升。并且环境压力越低,液滴温度下降越快,最低温度越低,温度回升越快。环境温度越高,液滴最低温度越高,温度回升越快。     

电场作用下基面液滴蒸发与内部流动数值模拟

为探究电场强化基面液滴蒸发的原理,本文采用有限元方法,对外加电场作用下的固体基面上液滴的蒸发过程进行了数值模拟,对比了不同电导率液滴的蒸发过程,分析了电场、液滴蒸发速率和内部流动的影响及其成因,以及液滴在电场作用下的内部流动与液滴传热传质的关系,结果表明,电场力的作用能够显著强化液滴内部的流动,对液滴的传热传质具有促进作用。此外,本文分析了温度对电场下基面液滴蒸发及内部流动的影响,发现温度对电场、液滴内部流动及蒸发的强化作用也有着较为明显的影响：对于电导率较低的纯水液滴,当电场强度低于和高于临界值6kV/cm时,温度对电场强化液滴内部流动和蒸发的影响有所不同;对于电导率较高的盐酸液滴,温度对电场强化液滴内部流动和蒸发的影响随电场强度升高均较大。本文为发展高效静电喷雾冷却技术提供了研究基础。     

平面上固定接触角蒸发液滴内Marangoni对流失稳现象

实验观察了水平加热基板上1cSt硅油液滴在固定接触角蒸发模式下的Marangoni对流失稳模式的演化规律,分析了接触角和基板温度对Marangoni对流不稳定性的影响。结果表明,随着基板温度的升高液滴内依次呈现热毛细对流、稳定的Bénard-Marangoni（BM）对流和不规则振荡的BM对流。对于稳定的BM对流,涡胞数随润湿半径的减小逐渐减少;当涡胞数少于5时,涡胞变为圆形;随着接触角的增大,由热毛细对流转捩为稳定的BM对流时的临界Marangoni数（Ma_c）增大;蒸发过程中,液滴内无量纲涡胞数随无量纲润湿半径的增大而线性增大,与接触角无关。     

水平及竖直基底上微小固着液滴的蒸发特性分析

为揭示非水平表面上微小蒸发液滴的传热传质特性,本文在准稳态模型的假设下构造三维液滴模型,综合考虑了蒸气扩散、蒸发冷却以及气相域中的自然对流这3种传输机理,对水平以及竖直基底上液滴的蒸发过程进行数值研究。通过分析气液界面上温度分布、蒸发通量分布及总蒸发率的变化,重点探究了基底过热度以及重力的改变对液滴蒸发特性的影响。结果表明：与水平基底上温度的对称分布不同,竖直基底上气液界面温度分布表现出明显的非对称性,且非对称性随基底过热度的升高而增强,最低温度点不再位于液滴顶点,而向一侧偏移。此外,水平基底上气液界面局部蒸发通量呈对称分布,各截面分布相似,而竖直基底上局部蒸发通量分布则呈现出显著的非对称性以及各截面异性,非对称性随着基底过热度的升高而增强,这是重力改变后气相域自然对流发生改变的结果。与水平基底相比,竖直基底上蒸发率更高,总蒸发时间更少。最后,基底由水平变为竖直时,液滴内部流场由对称双涡转变为非对称单涡,单涡流速显著大于双涡流速,液滴内流速随基底过热度的上升而增大,单涡环流造成了气液界面温度分布的改变以及最低温度点的偏移。     

Monodisperse monocomponent fuel droplet heating and evaporation

The results of numerical and experimental studies of heating and evaporation of monodisperse acetone, ethanol, 3-pentanone, n-heptane, n-decane and n-dodecane droplets in an ambient air of fixed temperature and atmospheric pressure are reported. The numerical model took into account the finite thermal conductivity of droplets and recirculation inside them based on the effective thermal conductivity model and the analytical solution to the heat conduction equation inside droplets. The effects of interaction between droplets are taken into account based on the experimentally determined corrections to Nusselt and Sherwood numbers. It is pointed out that the interactions between droplets lead to noticeable reduction of their heating in the case of ethanol, 3-pentanone, n-heptane, n-decane and n-dodecane droplets, and reduction of their cooling in the case of acetone. Although the trends of experimentally observed droplet temperatures and radii are the same as predicted by the model taking into account the interaction between droplets, the actual values of the predicted droplet temperatures can differ from the observed ones by up to about 8 K, and the actual values of the predicted droplet radii can differ from the observed ones by up to about 2%. It is concluded that the effective thermal conductivity model, based on the analytical solution to the heat conduction equation inside droplets, can predict the observed average temperature of droplets with possible errors not exceeding several K, and observed droplet radii with possible errors not exceeding 2% in most cases. These results allow us to recommend the implementation of this model into CFD codes and to use it for multidimensional modelling of spray heating and evaporation based on these codes.     

Effects of thermal airflow and mucus-layer interaction on hygroscopic droplet deposition in a simple mouth–throat model

Hygroscopic growth of inhaled aerosols plays an important role in determining particle trajectories and hence local deposition sites. Accurate predictions of airway temperature and humidity as well as droplet–vapor interaction are critical for the calculation of hygroscopic growth. Employing a simple mouth–throat (MT) model as a computer simulation test bed, the effects of interactive heat transfer between air–droplet flow and mucus-tissue-layer have been analyzed. For a steady inhalation flow rate of 15?L/min, air temperature and relative humidity distributions affecting droplet growth, deposition efficiency (DE), and deposition pattern have been compared for different thermal airway-wall conditions. The effects considered include: (i) the latent heat of mucus-layer evaporation and convection heat transfer; (ii) convection heat transfer only; and (iii) mucus-tissue layer with constant temperature. As the most important outcome, the validated modeling results show that thermal airflow and mucus-layer interaction can significantly reduce hygroscopic growth and thereby decrease the DE of multicomponent droplets up to 10%. The modeling framework presented can be readily expanded to other systems.

Copyright © 2018 American Association for Aerosol Research     

The separation of two different sized particles in an evaporating droplet

The separation of two different sized particles during evaporation of a dilute droplet is examined both computationally and experimentally. A transport model of the evaporating droplet system was solved using the finite element method to determine the fluid velocity, pressure, vapor concentration surrounding the droplet, temperature, and both particle concentrations. Experimentally, 1 μm and 3 μm polystyrene particles were used during the evaporation of a sessile water droplet. It was determined that to accurately model particle deposition, thermal effects need to be considered. The Marangoni currents in evaporating droplets keep particles suspended in the droplet until the end of the evaporation. Previous models of particle deposition during droplet evaporation have rapid accumulation of particles at the contact line. Our experiments and the experiments of others demonstrate that this is not accurate physically. In addition, to model the separation of two different sized particles the consideration of thermal effects is essential. © 2015 American Institute of Chemical Engineers AIChE J, 61: 3547–3556, 2015     

Internal fluid motion and particle transport in externally heated sessile droplets

Experiments and numerical simulations were carried out for an evaporating sessile droplet. The droplet was confined in the narrow gap between two glass plates, making it a “Hele‐Shaw” droplet and particle image velocimetry technique was used. In case of the evaporating droplet with pinned contact line and exposed to ambient condition, two symmetric but counterrotating convection cells were observed. After complete evaporation, the particles deposited on the substrate near the contact line. The direction of the flow was reversed for a droplet placed on uniformly heated substrate, and the final deposition pattern was a large spot at the center with a thin line at the periphery. For asymmetrically heated substrate a single convection cell appeared, and the final deposition was also asymmetric. When the liquid was subjected to localized heating, the contact line no longer remains pinned and a relatively uniform deposition was obtained after complete drying. © 2015 American Institute of Chemical Engineers AIChE J, 62: 1308–1321, 2016     

Influence of Plate Size on the Evaporation Rate of a Heated Droplet

The purpose of this study is to numerically investigate how the width of a plate influences natural convection around a droplet. Droplets evaporating on hot surfaces have many applications including drying of dishes and paint. Evaporation rate and deposition of particles withheld in the fluid are of great importance in both cases. As a first step to investigate how the drying rate and deposition mechanisms can be controlled, this work aims to investigate how the external flow around a water droplet influences the evaporation rate. Natural convection caused by the hot plate on which the droplet rests is considered and the effect of different widths is examined. Results show that an extension of the plate past the droplet will increase the maximum velocity in the domain due to natural convection while the flow close to the surface is decreased due to the no-slip condition and temperature gradient. A decrease of the evaporation rate is therefore observed when the plate is extended past the droplet as compared to the case when the plate and droplet have the same diameter. Simulations furthermore show that the results from the heat and mass transfer analogy only compare well to the results of Fick's law when the droplet and plate have the same width.     

Flow visualization and solute transport in evaporating droplets

We have investigated the velocity field and associated particle transport in an evaporating water droplet using the tool of particle image velocimetry. Experiments were performed where single droplets containing polystyrene particles were exposed to evaporation. Our method applicable to droplets confined between two parallel surfaces differs from the conventional PIV techniques on the 3D droplets and removes many of the limitations associated with mapping of velocity field. To avoid refraction of light at the droplet surface we have studied the motion in a disc‐shaped droplet which was prepared by confining the drop between two nonwetting surfaces and its base is pinned to a wetting surface. Experiments were carried out under the conditions where Marangoni flow creates convection cells and finally leading to deposition of particles toward the pinned edge. The contact angle, height of the droplet, velocity field, and the particle concentration inside the evaporating droplet was measured and its time evolution was recorded. © 2009 American Institute of Chemical Engineers AIChE J, 56: 1674–1683, 2010