混合润湿性柱状纳米结构对铜板上纳米氩膜沸腾传热的影响

混合润湿性对固/液相互作用有显著影响,因此对提高相变过程中的传热速率有积极作用.采用分子动力学模拟方法研究了柱状纳米结构表面混合润湿性对池沸腾传热的影响.分析了混合润湿性和纳米结构柱高对液体起始沸腾时间和温度的影响及其机理.结果表明,疏水比例和柱高会影响爆沸的起始温度和时间.与纯亲水壁相比,增加疏水比改变了固液界面性质,可以降低沸腾温度,更容易突破势能壁垒,使液体起始沸腾时间提前,并且随着疏水比的增加,不同柱高下的沸腾温度降低;当疏水比相同时,增加柱高扩大了混合润湿性的影响,也能降低沸腾起始温度并使液体起始沸腾时间提前.这为设计微纳粗糙结构和混合润湿表面以强化沸腾传热提供了思路.     

蒸汽冷凝过程的分子团聚模型及其对不凝性气体影响传热性能的解释

根据相变过程的微观物理机理和热力学特性,提出了冷凝传热过程中,近壁面蒸汽分子经由团聚阶段进而冷凝成宏观液滴的物理模型.并将团聚体分布与滴状冷凝传热性能相联系,从而研究不凝性气体对滴状冷凝传热过程的影响.在改进的Dillmann和Meier(DM)模型基础上,将分子团聚过程中的临界过饱和度与冷凝过程中的过冷度相联系,以及将团聚体的能量特性与液固界面物理化学特性相联系,将团聚模型与考虑固液界面效应的滴状冷凝传热模型相联系,建立了近壁面条件影响的分子团聚模型.利用模型计算了近壁面蒸汽中团簇体尺寸和分布,以及不凝性气体存在导致的蒸汽冷凝团聚体分布的变化,并结合滴状冷凝传热模型,定量解释了少量不凝性气体的存在,极大影响了冷凝传热性能的现象.模型计算结果与实验结果及文献中含不凝气的蒸汽冷凝传热实验数据进行了比较,两者符合较好,验证了所提出模型的合理性.     

喷雾冷却相变传热的建模与模拟计算

喷雾冷却是一种高效的热控技术,为了探索形成完善的喷雾冷却技术设计流程,文章开展了喷雾冷却传热过程的建模研究.针对喷雾冷却传热过程的模拟计算,基于喷雾冷却相变传热的4个传热机制:液膜对流传热、池沸腾传热、二次气泡沸腾传热、二次气泡高频化机制,利用Monte Carlo方法描述了不同粒径与速率分布的液滴撞击液膜并生成二次气...     

Fe合金FCC-BCC原子尺度台阶型马氏体相界面迁移行为的分子动力学模拟研究

两相界面的原子尺度结构对相界面迁移行为具有重要影响.高分辨透射电子显微分析表明钢中马氏体相界面具有高度为若干原子层间距的台阶结构,然而目前Fe合金马氏体相变的模拟研究工作中绝大多数使用非台阶型相界面结构作为模拟初始模型.本文基于拓扑模型和相变位错理论构建了Fe合金FCC/BCC台阶型相界面初始模型,采用分子动力学模拟方法研究了Fe合金马氏体相界面的迁移行为.研究结果表明,当两相界面具有约束共格匹配关系及台阶结构时,体系发生FCC→BCC马氏体相变并呈现典型的非扩散切变特征;相变过程中FCC/BCC宏观尺度相界面沿其法线方向以(4.4±0.3)×10~2 m/s的速度迁移,且相界面在迁移过程中始终保持稳定的台阶结构和相对平直的宏观界面形貌特征;相变位错的滑移速度高达(2.8±0.2)×10~3 m/s,相变位错阵列沿台阶面的协同侧向滑移不仅是马氏体台阶结构宏观相界面迁移的微观机制,也是马氏体相变宏观形状应变的主要来源;采用分子动力学模拟方法获得的Fe合金马氏体相变晶体学特征参量与拓扑模型的解析解数值非常接近,相变产生的整体宏观形状应变由平行于相界面的剪切应变和垂直于相界面的法向应变两部分组成.     

水平管外开孔铜泡沫R134a池沸腾换热实验研究

本文对6根不同几何参数开孔铜金属泡沫覆盖层铜管外R134a的池沸腾换热性能进行了试验研究。实验结果表明铜金属泡沫对R134a池沸腾相变传热能起到很好的强化作用,在低热流密度范围内其沸腾传热系数约为光管的2.6～4.4倍。     

液滴撞击热壁面的相变行为研究

本文基于格子Boltzmann方法建立了液滴撞击热壁面的数学模型,研究了撞击过程液滴相变传热机制,重点分析了壁面温度与We数的影响规律.结果表明,随着壁面温度升高,液滴撞击壁面出现接触沸腾、膜态沸腾相变模式.接触沸腾过程中液滴内部与壁面接触的三相线附近出现最大局部热流.当壁面温度超过Leidenfrost点后,撞击液滴因界面蒸发产生的气膜托举作用而与壁面无接触地铺展、回缩直至完全反弹形成膜态沸腾,这一阶段液滴所能达到的最大铺展因子与液滴We数呈幂次关系(指数为0.313).     

固-液相变糊状区的格子Boltzmann研究

为深入探讨固-液相变材料在相界面内(糊状区)的流动和传热对固-液相变过程的影响,本文基于"焓-多孔介质"模型,提出研究糊状区的"多相流-多孔介质"两区域复合模型,将糊状区的高含液率视为多相流区域,而低含液率区域视为多孔介质,采用格子Boltzmann方法对固-液相变糊状区的流动和传热过程进行研究.重点讨论了相变温度半径,不同高、低含液率区的分界点γ_(tr)对糊状区的发展、迁移以及对相变过程的影响.研究结果表明糊状区厚度随相变半径的增大而增加,进而影响相变过程中流动和换热;对比复合模型的不同高、低含液率分界点γ_(tr)时和单纯多孔介质模型时的糊状区的边界、流动和换热情况,表明将糊状区高液相率区采用多相流模型进行分析,可以正确反映出高含液率区域内固相微粒随液相流动的特性及其对相变过程的影响。     

光纤探针测量及多尺度熵鉴别超临界类沸腾传热模式

亚临界沸腾包括界面蒸发和气泡动力学诱导的传热,但超临界传热是否存在类界面蒸发和类气泡传热以及两者间的转换缺少直接的实验证据.本文进行了超临界CO₂液池传热的实验研究,压力和液池温度分别为8—10 MPa和15℃.作为加热元件和感温元件,22 mm长和70μm直径的镍铬丝水平放置在液池中,光纤探针垂直放置,其顶端高于镍铬丝200μm.发现随热流密度或壁面过热度的持续增大,依次发生自然对流、类界面蒸发、类蒸发-沸腾转换、类沸腾4种传热模式.本文重点关注类界面蒸发和类沸腾传热以及两者间的转换.在类界面蒸发模式下,传热系数随壁面过热度增大略有下降,光纤输出小幅/高频信号,不存在主频,多尺度熵大,表征随机信号波动.在类蒸发-沸腾转换模式下,光纤输出大幅/低频周期信号,存在明显主频,多尺度熵小,代表有序的周期性脉动传热.在以类气泡为特征的类沸腾模式下,光纤信号波动幅度介于类蒸发和转换模式之间,主频不明显,多尺度熵也介于类蒸发和转换模式之间.研究获得了超临界类沸腾直接的实验证据,加深了对超临界传热机理的理解,为后续理论研究和工程应用提供了基础.     

相变微胶囊复合板相变传热数值研究

本文以相变微胶囊复合板隔热组件为研究对象,基于焓法方程建立了蓄热/隔热型耦合相变传热数值模型,预测了板内的温度分布、相界面移动等响应特征,系统地分析了Stefan数及相变温度区间对传热的影响。结果表明:相变微胶囊板内固相界面的移动速度大于液相界面,固相界面消失后,液相界面移动明显加速;固、液相界面移动速度随Stefan数的增大而增快,同时板内整体温度水平升高;固相界面移动速度随着相变温度区间的增大而增快,但相变温度区间对相变材料完全熔化所需时间影响不大。     

纳米多晶体的热力学函数及其在相变热力学中的应用

以往关于纳米材料热力学的研究，绝大多数以界面的热力学函数表征整体纳米材料的热力学性质，这种近似处理，对于尺寸超过几十纳米的较粗纳米材料，在相变热力学中对特征转变温度和临界尺寸等重要参量的预测，将导致很大误差. 应用“界面膨胀模型”和普适状态方程，研究了纳米晶界的热力学特性，进一步发展了纳米晶整体材料热力学函数的计算模型，给出了单相纳米多晶体的焓、熵和吉布斯自由能随界面过剩体积、温度，以及晶粒尺寸发生变化的明确表达式. 以Co纳米晶为例，分析了界面与整体纳米多晶体热力学函数的差异，确定了相变温度与晶粒尺寸的依赖关系，以及一定温度下可能发生相变的临界尺寸. 关键词： 纳米多晶体 热力学函数 相变热力学     

Hydrodynamics and Platoon Formation for a Totally Asymmetric Exclusion Model with Particlewise Disorder

We consider a one-dimensional totally asymmetric exclusion model with quenched random jump rates associated with the particles, and an equivalent interface growth process on the square lattice. We obtain rigorous limit theorems for the shape of the interface, the motion of a tagged particle, and the macroscopic density profile on the hydrodynamic scale. The theorems are valid under almost every realization of the disordered rates. Under suitable conditions on the distribution of jump rates the model displays a disorder-dominated low-density phase where spatial inhomogeneities develop below the hydrodynamic resolution. The macroscopic signature of the phase transition is a density discontinuity at the front of the rarefaction wave moving out of an initial step-function profile. Numerical simulations of the density fluctuations ahead of the front suggest slow convergence to the predictions of a deterministic particle model on the real line, which contains only random velocities but no temporal noise.     

Stress-driven phase transformation and the roughening of solid-solid interfaces

The application of stress to multiphase solid-liquid systems often results in morphological instabilities. Here we propose a solid-solid phase transformation model for roughening instability in the interface between two porous materials with different porosities under normal compression stresses. This instability is triggered by a finite jump in the free energy density across the interface, and it leads to the formation of fingerlike structures aligned with the principal direction of compaction. The model is proposed as an explanation for the roughening of stylolites-irregular interfaces associated with the compaction of sedimentary rocks that fluctuate about a plane perpendicular to the principal direction of compaction.     

Flame balls with thermally sensitive intermediate kinetics

Spherical flame balls are studied using a model for the chemical kinetics which involves a non-exothermic autocatalytic reaction, describing the chain-branching generation of a chemical radical and an exothermic completion reaction, the rate of which does not depend on temperature. When the chain-branching reaction has a large activation temperature, an asymptotic structure emerges in which the branching reaction generates radicals and consumes fuel at a thin flame interface, although heat is produced and radicals are consumed on a more distributed scale. Another model, based more simply, but less realistically, on the generation of radicals by decomposition of the fuel, provides exactly the same leading order matching conditions. These can be expressed in terms of jump conditions across a reaction sheet that are linear in the dependent variables and their normal gradients. Using these jump conditions, a reactive–diffusive model with linear heat loss then leads to analytical solutions that are multivalued for small enough levels of heat loss, having either a larger or a smaller radius of the interface where fuel is consumed. The same properties are found, numerically, to persist as the activation temperature of the branching reaction is reduced to values that seem to be typical for hydrocarbon chemistry. Part of the solution branch with larger radius is shown to become stable for low enough values of the Lewis number of the fuel.     

Embedding sharp interfaces within the lattice Boltzmann method for fluids with arbitrary density ratios

Current lattice Boltzmann methods for simulating two fluids create a diffuse interface between the fluids. In this work, we develop a novel technique for embedding sharp interfaces between fluids with unbounded density ratios for the LB method. Distribution functions streamed across an interface are transformed so that the receiving node is passed information corresponding to its fluid phase. Two different methods are employed to determine the transformation. The first uses analytical distribution functions from steady Poiseuille flow to determine the jump in moments of the distribution functions across the interface. The second uses approximate expansions of distribution functions to determine jumps in distribution functions. The accuracy and stability of the methods are examined in simulations of Poiseuille-Couette flows with an interface parallel to the walls. Both methods show linear convergence to the analytical solution.     

Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves

This paper is devoted to time domain numerical solutions of two-dimensional (2D) material interface problems governed by the transverse magnetic (TM) and transverse electric (TE) Maxwell's equations with discontinuous electromagnetic solutions. Due to the discontinuity in wave solutions across the interface, the usual numerical methods will converge slowly or even fail to converge. This calls for the development of advanced interface treatments for popular Maxwell solvers. We will investigate such interface treatments by considering two typical Maxwell solvers – one based on collocation formulation and the other based on Galerkin formulation. To restore the accuracy reduction of the collocation finite-difference time-domain (FDTD) algorithm near an interface, the physical jump conditions relating discontinuous wave solutions on both sides of the interface must be rigorously enforced. For this purpose, a novel matched interface and boundary (MIB) scheme is proposed in this work, in which new jump conditions are derived so that the discontinuous and staggered features of electric and magnetic field components can be accommodated. The resulting MIB time-domain (MIBTD) scheme satisfies the jump conditions locally and suppresses the staircase approximation errors completely over the Yee lattices. In the discontinuous Galerkin time-domain (DGTD) algorithm – a popular Galerkin Maxwell solver, a proper numerical flux can be designed to accurately capture the jumps in the electromagnetic waves across the interface and automatically preserves the discontinuity in the explicit time integration. The DGTD solution to Maxwell interface problems is explored in this work, by considering a nodal based high order discontinuous Galerkin method. In benchmark TM and TE tests with analytical solutions, both MIBTD and DGTD schemes achieve the second order of accuracy in solving circular interfaces. In comparison, the numerical convergence of the MIBTD method is slightly more uniform, while the DGTD method is more flexible and robust.     

High-order FDTD methods for transverse electromagnetic systems in dispersive inhomogeneous media

This Letter introduces a novel finite-difference time-domain (FDTD) formulation for solving transverse electromagnetic systems in dispersive media. Based on the auxiliary differential equation approach, the Debye dispersion model is coupled with Maxwell's equations to derive a supplementary ordinary differential equation for describing the regularity changes in electromagnetic fields at the dispersive interface. The resulting time-dependent jump conditions are rigorously enforced in the FDTD discretization by means of the matched interface and boundary scheme. High-order convergences are numerically achieved for the first time in the literature in the FDTD simulations of dispersive inhomogeneous media.     

Hydrodynamics of binary fluid phase segregation

Starting with the Vlasov-Boltzmann equation for a binary fluid mixture, we derive an equation for the velocity field u when the system is segregated into two phases (at low temperatures) with a sharp interface between them. u satisfies the incompressible Navier-Stokes equations together with a jump boundary condition for the pressure across the interface which, in turn, moves with a velocity given by the normal component of u. Numerical simulations of the Vlasov-Boltzmann equations for shear flows parallel and perpendicular to the interface in a phase segregated mixture support this analysis. We expect similar behavior in real fluid mixtures.     

Axisymmetric convection flow of fractional Maxwell fluid past a vertical cylinder with velocity slip and temperature jump

A novel finite volume method is developed to investigate the axisymmetric convection flow and heat transfer of fractional viscoelastic fluid past a vertical cylinder. Fractional cylindrical governing equations are formulated by fractional Maxwell model and generalized Fourier's law. The velocity slip and temperature jump boundary conditions are considered across the fluid-solid interface. Numerical results are validated by exact solutions of special case with source terms. The effects of fractional derivative parameter and boundary condition parameters on flow and heat transfer characteristics are discussed. The viscoelastic fluid performs evident shear thickening property in the fractional Maxwell constitutive relation. Moreover, the boundary condition parameters have remarkable influence on velocity and temperature distributions.     

A matched interface and boundary method for solving multi-flow Navier–Stokes equations with applications to geodynamics

We have developed a second-order numerical method, based on the matched interface and boundary (MIB) approach, to solve the Navier–Stokes equations with discontinuous viscosity and density on non-staggered Cartesian grids. We have derived for the first time the interface conditions for the intermediate velocity field and the pressure potential function that are introduced in the projection method. Differentiation of the velocity components on stencils across the interface is aided by the coupled fictitious velocity values, whose representations are solved by using the coupled velocity interface conditions. These fictitious values and the non-staggered grid allow a convenient and accurate approximation of the pressure and potential jump conditions. A compact finite difference method was adopted to explicitly compute the pressure derivatives at regular nodes to avoid the pressure–velocity decoupling. Numerical experiments verified the desired accuracy of the numerical method. Applications to geophysical problems demonstrated that the sharp pressure jumps on the clast-Newtonian matrix are accurately captured for various shear conditions, moderate viscosity contrasts and a wide range of density contrasts. We showed that large transfer errors will be introduced to the jumps of the pressure and the potential function in case of a large absolute difference of the viscosity across the interface; these errors will cause simulations to become unstable.     

Interface model for non-equilibrium evaporation

A microscopic interface condition for condensing/evaporating interfaces is developed by combining a velocity dependent condensation probability [T. Tsuruta, H. Tanaka, T. Masuoka, Int. J. Heat Mass Transfer 42 (1999) 4107] and Maxwell type interface conditions with accommodation. Using methods from kinetic theory, macroscopic interface conditions for mass and energy transport across the phase boundary are derived. This model only applies to simple substances, where diffusive effects in the bulk phases are not present. The results are compared to classical non-equilibrium thermodynamics. The interface conditions are considered for the limit of small deviation from equilibrium, and the corresponding Onsager coefficients are computed. These results are useful as boundary conditions for non-equilibrium evaporation and condensation problems, as done previously by our group [M. Bond, H. Struchtrup, Phys. Rev. E 70 (2004) 061605].