On the theory of laser heating of water drops

It is shown that pressure caused by the law of momentum conservation and a nonuniform distribution of energy sources over a drop volume should be taken into account when describing the evaporation of large water drops subjected to high-intensity optical radiation.     

Evaporation of Charged Drops

A model of evaporation of a multiply charged liquid drop is developed. The model self-consistently takes into account the main factors influencing the charged drop evaporation, including effects of the drop surface curvature and charge on the saturated vapor pressure, repeated fragmentation of drops during evaporation, and the capability of drops having a unit charge and a certain stabilization radius not to evaporate even in an unsaturated vapor medium. Analytical dependences are derived that can be used to calculate an integral lifetime of a charged drop with allowance for its fragmentation into smaller drops. Our estimates demonstrate that the evaporation time of charged drops is much smaller than that of uncharged drops.     

Influence of the accommodation coefficient on the water drop evaporation process in a radiation field

The process of water drop evaporation in a field of intense laser radiation is examined on the basis of the hydro-thermodynamics equations under the assumption of quasistationarity in the conditions in the surrounding medium and in the radius. The influence of the accommodation coefficient on heating of the drops on the magnitude of the surface jump in water vapor density and on the position of the upper bound of the convective evaporation region is analyzed. It is shown that the surface jump in the vapor does not alter the rate of convective evaporation of the drop in practice, but can result in substantially different time dependences of the radius of a diffusely evaporating drop as compared with those found without it. The solution obtained is compared with the solution of the problem of evaporation in a Stefan approximation.     

热质的运动和传递-热子气的守恒方程和傅立叶定律

本文在热质和热子气概念的基础上建立了热子气的质量、动量、能量守恒方程.基于傅立叶导热定律求得了热子气粘性系数和粘性力的近似表达式.分析了傅立叶定律的物理意义:傅立叶定律是在忽略惯性力的条件下对热子气动量方程的近似.在极低温度或极高热流密度条件下傅立叶定律不再适用.     

Dynamics of magnetic vortices

The canonical conservation laws of linear and angular momentum in the ferromagnetic continuum have been known to be plagued by certain ambiguities which are resolved in this paper by constructing conservation laws as suitable moments of a topological density. The resulting canonical structure is then shown to be analogous to that encountered in the familiar Hall effect and explains the unusual features of the dynamics of magnetic vortices without resorting to a detailed solution of the underlying nonlinear equations. Thus, in the absence of external magnetic fields, a magnetic vortex is shown to be spontaneously pinned around a fixed guiding center. The guiding center would drift in a direction perpendicular to an applied magnetic field gradient, provided that dissipation can be neglected, with a Hall velocity that is calculated explicitly in terms of the initial configuration of the vortex. In the presence of dissipation, the vortex undergoes skew deflection at an angle δ ≠ 90° with respect to the applied field gradient. The angle δ is related to the winding number of the vortex according to the well-known golden rule of bubble dynamics.     

Large Eddy Simulation composition equations for single-phase and two-phase fully multicomponent flows

The Large Eddy Simulation (LES) equations for multicomponent (MC) fuel single-phase (SP) flow and two-phase (TP) flow with phase change are derived from the Direct Numerical Simulation (DNS) equations by filtering the DNS equations using a top-hat filter. Additional to the equations solved for single-component (SC) fuels, composition equations enter the formulation. The species composition is represented through a Probability Distribution Function (PDF), and DNS equations for the PDF moments are solved to find the composition. The TP filtered equations contain three categories of subgrid-scale (SGS) terms: (1) SGS–flux terms, (2) filtered source terms (FSTs) and (3) terms representing the ‘LES assumptions’. For SP flows no FSTs exist. The SGS terms in the LES equations must be either shown negligible or modeled. It is shown that for the composition equations, two equivalent forms of the DNS equations lead to two non-equivalent forms of the LES equations. Criteria are proposed to select the form best suited for LES. These criteria are used in conjunction with evaluations based on a DNS database portraying mixing and phase change, and lead to choosing one of the LES forms which satisfies all criteria. It is shown that the LES assumptions lead to additional SGS terms which require modeling. Further considerations are made for reactive flows.     

A spray flamelet/progress variable approach combined with a transported joint PDF model for turbulent spray flames

A spray flamelet/progress variable approach is developed for use in spray combustion with partly pre-vaporised liquid fuel, where a laminar spray flamelet library accounts for evaporation within the laminar flame structures. For this purpose, the standard spray flamelet formulation for pure evaporating liquid fuel and oxidiser is extended by a chemical reaction progress variable in both the turbulent spray flame model and the laminar spray flame structures, in order to account for the effect of pre-vaporised liquid fuel for instance through use of a pilot flame. This new approach is combined with a transported joint probability density function (PDF) method for the simulation of a turbulent piloted ethanol/air spray flame, and the extension requires the formulation of a joint three-variate PDF depending on the gas phase mixture fraction, the chemical reaction progress variable, and gas enthalpy. The molecular mixing is modelled with the extended interaction-by-exchange-with-the-mean (IEM) model, where source terms account for spray evaporation and heat exchange due to evaporation as well as the chemical reaction rate for the chemical reaction progress variable. This is the first formulation using a spray flamelet model considering both evaporation and partly pre-vaporised liquid fuel within the laminar spray flamelets. Results with this new formulation show good agreement with the experimental data provided by A.R. Masri, Sydney, Australia. The analysis of the Lagrangian statistics of the gas temperature and the OH mass fraction indicates that partially premixed combustion prevails near the nozzle exit of the spray, whereas further downstream, the non-premixed flame is promoted towards the inner rich-side of the spray jet since the pilot flame heats up the premixed inner spray zone. In summary, the simulation with the new formulation considering the reaction progress variable shows good performance, greatly improving the standard formulation, and it provides new insight into the local structure of this complex spray flame.     

三维煤粉燃烧全双流体模型的数值模拟

1引言目前煤粉燃烧综合模型最为广泛的是以轨道模型为基础，美国杨伯翰大学的先进燃烧工程中心（CERC）自1980年起研制二维煤粉燃烧程序PCGC－2，从1990年起研制三维煤粉燃烧的PCGC－3程序[1]。与轨道模型发展的同时，美国Rabcock＆Wilcox公司Fiveland[2]等人研制了FURMO程序，用无滑移模型对560MW侧墙喷燃煤粉炉进行了三维全模拟。其特点是首次用全欧拉的处理方法计算三维煤粉燃烧过程，其不足之处是不考虑气粒两相间的速度滑移和温度滑移．总的看来，用轨道模型模拟煤粉燃烧，易于考虑颗粒反应经历，也可给出两相之间的速度及温…     

声子气的状态方程和声子气运动的守恒方程

根据爱因斯坦的狭义相对论中质能的等效关系，把固体(本文指非导体)晶格(原子)的质量分为晶格(原子)的静质量和晶格热振动能量的等效质量两个部分，后者就是固体中声子气的等效质量.晶格(原子)热振动的能量则分为晶格(原子)静质量具有的热能以及声子气质量具有的热能.基于固体的状态方程，导得了晶格静质量热振动的状态方程和声子气的状态方程.声子气在固体介质中的宏观运动就是热量在固体中的传递过程.建立了声子气运动的守恒方程组，分析表明，忽略惯性力时声子气的动量守恒方程就退化为傅里叶导热定律，阐明了傅里叶导热定律的物理本质是声子气驱动力与阻力的平衡方程.当热流密度很大惯性力不能忽略时，傅里叶导热定律不再适用. 关键词： 非傅里叶导热 声子气 声子气质量 状态方程 守恒方程     

燃料抛撒成雾及其燃烧爆炸的光滑离散颗粒流体动力学方法数值模拟研究

燃料在炸药爆炸驱动下形成燃料空气爆炸云团, 进而引燃爆炸, 对目标造成毁伤. 本文在前期提出的光滑离散颗粒流体动力学方法(SDPH)的基础上, 引入描述炸药由爆轰到膨胀整个过程的Jones-Wilkins-Lee状态方程及描述气体快速燃烧过程的EBU-Arrhenius燃烧模型, 建立了求解战斗部起爆、燃料抛撒和燃料二次引燃爆炸问题的新型SDPH方法. 设计了圆环形燃料颗粒在炸药爆炸驱动下运动抛撒的算例进行数值验证, 结果与理论相符; 对燃料空气炸药(FAE)云雾的形成和发展过程进行了数值模拟, 分析了云雾的形态, 并与实验结果进行对比, 符合较好, 同时分析了不同起爆方式对云雾团成型的影响; 最后, 在云雾团成型的基础上, 引入蒸发燃烧模型对FAE的燃烧爆炸过程进行了模拟研究. 结果表明, 本文建立的数学模型和计算方法可以较好的模拟燃料空气炸药抛撒成雾及云雾燃烧爆炸过程, 为该类武器装备的设计研究提供了较好的数值方法.     

Nonlinear vortex dynamics in open nonequilibrium systems with bulk mass loss and a generation mechanism for tornadoes and typhoons

Based on a general model of nonlinear vortex dynamics in open thermodynamically nonequilibrium systems with bulk or surface mass losses, an analysis is presented of the mechanism of generation of violent atmospheric vortices (tornadoes, typhoons, cyclones) associated with the formation of deep cloud systems by intense condensation of water vapor from moist air cooled below the dew point. Simple particular solutions to the Navier-Stokes equations are found that describe both axisymmetric and nonaxisymmetric incompressible vortex motions involving radial and vertical flows with viscous dissipation vanishing identically everywhere except for a thin shear layer at the boundary of the condensation region. It is shown that the nonlinear convective and local Coriolis forces generated by radial inflow in the presence of a background vorticity due to a global Coriolis force (the Earth’s rotation) accelerate the solid-body rotation in the vortex core either exponentially or in a nonlinear regime of finite-time blow-up. Due to updrafts, such a vortex is characterized by a strong helicity. This mechanism explains a number of observed properties and characteristics of the structure and evolution of tornadoes and typhoons. Upper estimates are found for the kinetic energies of violent atmospheric vortices. It is shown that increase in rotational kinetic energy of atmospheric vortices with constant vortex-core radii is consistent with energy and momentum conservation, because radial inflow continually supplies the required amount of rotational kinetic energy drawn from the ambient atmosphere to an open system.     

Relationship between the diffusiophoresis rate, evaporation coefficient, and size of a large volatile drop

A new theory of diffusiophoresis of large volatile spherical aerosol drops that is an extension of investigations [1–8] is developed. The influence of the radius of the drop, the surface tension coefficient varying over the surface of the drop, the evaporation coefficient α of the liquid, and the flows inside the drop on the diffusiophoresis rate are taken into account. Expressions obtained allow for direct determination of the velocity of large individual aerosol drops in a binary gas mixture nonuniform in component concentration. It is shown that both the magnitude and the direction of the diffusiophoresis velocity depend on α and the size of the drop. It is assumed that the size of the drop varies but remains considerably greater than the mean free path of gas molecules.     

Description of evaporation of a spherical liquid drop by a non-Markovian random process based on integral stochastic equations

Processes of evaporation (condensation) of vapor particles from the surface of a spherical drop and processes of their diffusion into surrounding volume are considered. Special features of evaporation are investigated taking into account vapor particle fluctuations caused by random changes in the temperature, concentration, etc. Statistical characteristics of fluctuations of the corresponding quantities, including the mass flow through the liquid-vapor boundary and concentration on the liquid surface, are presented. The distribution of completely evaporated drop number versus time is presented.     

Shear-thinning-induced chaos in taylor-couette flow

The effect of weak shear thinning on the stability of the Taylor-Couette flow is explored for a Carreau-Bird fluid in the narrow-gap limit. The Galerkin projection method is used to derive a low-order dynamical system from the conservation of mass and momentum equations. In comparison with the Newtonian system, the present equations include additional nonlinear coupling in the velocity components through the viscosity. It is found that the critical Taylor number, corresponding to the loss of stability of the base (Couette) flow, becomes lower as the shear-thinning effect increases. That is, shear thinning tends to precipitate the onset of Taylor vortex flow. Similar to Newtonian fluids, there is an exchange of stability between the Couette and Taylor vortex flows, which coincides with the onset of a supercritical bifurcation. However, unlike the Newtonian model, the Taylor vortex cellular structure loses its stability in turn as the Taylor number reaches a critical value. At this point, a Hopf bifurcation emerges, which exists only for shear-thinning fluids.     

Sound propagation in saturated gas–vapor–droplet suspensions considering the effect of transpiration on droplet evaporation

The sound attenuation and dispersion in saturated gas–vapor–droplet mixtures with evaporation has been investigated theoretically. The theory is based on an extension of the work of Davidson (1975) to accommodate the effects of transpiration on the linear particle relaxation processes of mass, momentum and energy transfer. It is shown that the inclusion of transpiration in the presence of mass transfer improves the agreement between the theory and the experimental data of Cole and Dobbins (1971) for sound attenuation in air–water fogs at low droplet mass concentrations. The results suggest that transpiration has an appreciable effect on both sound absorption and dispersion for both low and high droplet mass concentrations.     

Kinetic boundary layers in gas mixtures: Systems described by nonlinearly coupled kinetic and hydrodynamic equations and applications to droplet condensation and evaporation

We consider a mixture of heavy vapor molecules and a light carrier gas surrounding a liquid droplet. The vapor is described by a variant of the Klein-Kramers equation, a kinetic equation for Brownian particles moving in a spatially inhomogeneous background; the gas is described by the Navier-Stokes equations; the droplet acts as a heat source due to the released heat of condensation. The exchange of momentum and energy between the constituents of the mixture is taken into account by force terms in the kinetic equation and source terms in the Navier-Stokes equations. These are chosen to obtain maximal agreement with the irreversible thermodynamics of a gas mixture. The structure of the kinetic boundary layer around the sphere is then determined from the self-consistent solution of this set of coupled equations with appropriate boundary conditions at the surface of the sphere. For this purpose the kinetic equation is rewritten as a set of coupled moment equations. A complete set of solutions of these moment equations is constructed by numerical integration inward from the region far away from the droplet, where the background inhomogeneities are small. A technique developed in an earlier paper is used to deal with the severe numerical instability of the moment equations. The solutions so obtained for given temperature and pressure profiles in the gas are then combined linearly in such a way that they obey the boundary conditions at the droplet surface; from this solution source terms for the Navier-Stokes equation of the gas are constructed and used to determine improved temperature and pressure profiles for the background gas. For not too large temperature differences between the droplet and the gas at infinity, self-consistency is reached after a few iterations. The method is applied to the condensation of droplets from a supersaturated vapor, where small but significant corrections to an earlier, not fully consistent version of the theory are found, as well as to strong evaporation of droplets under the influence of an external heat source, where corrections of up to 40 % are obtained.     

Mechanism of electrohydrodynamically induced flow in a wire-non-parallel plate electrode type gas pump

An experimental investigation and one-dimensional modeling have been conducted to study the mechanism of net flow direction induced by electrohydrodynamic (EHD) forces in a wire-non-parallel plate electrode type EHD gas pump. The experiments were conducted with various different locations of corona wire electrode for negative and positive applied voltage from 0 to 14 kV at atmospheric pressure and room temperature, where air was used as the working fluid. A one-dimensional cross-sectional averaged model based on mass and momentum conservation as well as Poisson electric field and ion transport equations was also developed. The results show that the net flow direction of electrohydrodynamically induced gas flow in a wire-non-parallel plate electrode system significantly depends on the location of the corona wire electrode relative to the grounded electrode position. The effect of conversion angle of non-parallel plate electrode on the net flow direction and pressure drop also was investigated and discussed in detail.     

Continuity and momentum equations for moist atmospheres

The moist atmosphere with occurring precipitation is considered to be a multiphase fluid composed of dry air, water vapor and hydrometeors. These compositions move with different velocities: they take a macroscopic motion with a reference velocity and a relative motion with a velocity deviated from the reference velocity. The reference velocity can be chosen as the velocities of dry air, a gas mixture and the total air mixture. The budget equations of continuity and momentum are formulated in the three reference-velocity frames. It is shown that the resulting equations are dependent on the chosen reference velocity. The diffusive flux due to compositions moving with velocities deviated from the reference velocity and the internal sources due to the phase transitions of water substances result in additional source terms in continuity and momentum equations. A continuity equation of the total mass is conserved and free of diffusive flux divergence if the reference velocity is referred to the velocity of the total air mixture. However, continuity equations in the dry-air and gasmixture frames are not conserved due to the mass diffusive flux divergence. The diffusive flux introduces additional source terms in the momentum equation. In the dry-air frame, the diffusive flux of water substances and the phase transitions of water substances contribute to the change of the total momentum. The additional sources of total momentum in the frame of a gas mixture are associated with the diffusive flux of hydrometeors, the phase transitions of hydrometeors and the gasmixture diffusive flux. In the frame of total air mixture, the contribution to the total momentum comes from the diffusive flux of all atmospheric compositions instead of the phase transitions. The continuity and momentum equations derived here are more complicated than the traditional model equations. With increasing computing power, it becomes possible to simulate atmospheric processes with these sophisticated equations. It is helpful to the improvement of precipitation forecast.