Hydrodynamic binary coalescence of droplets under air flow in a hydrophobic microchannel

Based on the volume of fluid(VOF) method, we conduct a numerical simulation to study the hydrodynamic binary coalescence of droplets under air flow in a hydrophobic rectangular microchannel. Two distinct regimes, coalescence followed by sliding motion and that followed by detaching motion, are identified and discussed. Additionally, the detailed hydrodynamic information behind the binary coalescence is provided, based on which a dynamic mechanical analysis is conducted to reveal the hydrodynamic mechanisms underlying these two regimes. The simulation results indicate that the sliding motion of droplets is driven by the drag force and restrained by the adhesion force induced by the interfacial tension along the main flow direction. The detachment(i.e., upward motion) of the droplet is driven by the lift force associated with an aerodynamic lifting pressure difference imposed on the coalescent droplet, and also restrained by the adhesion force perpendicular to the main flow direction. Especially, the lift force is mainly induced by an aerodynamic lifting pressure difference imposed on the coalescent droplet. Two typical regimes can be quantitatively recognized by a regime diagram depending on Re and We. The higher Re and We respectively lead to relatively larger lift forces and smaller adhesion forces acting on the droplet, both of which are helpful to detachment of the coalesced droplet.     

Capillary filling for multicomponent fluid using the pseudo-potential Lattice Boltzmann method

We present a systematic study of capillary filling for a binary fluid by using mesoscopic a lattice Boltzmann model describing a diffusive interface moving at a given contact angle with respect to the walls. We compare the numerical results at changing the ratio the typical size of the capillary, H, and the wettability of walls. Numerical results yield quantitative agreement with the theoretical Washburn law, provided that the channel height is sufficiently larger than the interface width and variations of the dynamic contact angle with the capillary number are taken into account.     

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.     

Molecular dynamics study of phase separation kinetics in thin films

We use molecular dynamics to simulate experiments where a symmetric binary fluid mixture (AB), confined between walls that preferentially attract one component (A), is quenched from the one-phase region into the miscibility gap. Surface enrichment occurs during the early stages, yielding a B-rich mixture in the film center with well-defined A-rich droplets. The droplet size grows with time as l(t) proportional t(2/3) after a transient regime. The present atomistic model is also compared to mesoscopic coarse-grained models for this problem.     

偏晶合金液-液相分离的格子玻尔兹曼方法模拟

本文建立了一个模拟在弥散相液滴的扩散长大、碰撞凝并和Ostwald熟化等因素的作用下偏晶合金液-液相分离过程的二维格子玻尔兹曼方法 (lattice Boltzmann method, LBM) 模型.该模型结合了Shan-Chen的两相流模型和Qin的介观粒子相互作用势模型的优点,并在LB演化方程中引入了反映相变的源项.应用该模型模拟研究了偏晶合金液-液相分离过程中单液滴的生长、两液滴的合并和多液滴的生长规律.结果表明在两液相区中第二相单个液滴的生长是一个通过扩散从非平衡态到平衡态过渡的过程.两液滴合并 关键词： 偏晶合金 液-液相分离 格子玻尔兹曼方法     

The theory of droplet diffusiophoresis for concentrated solutions

The theoretical diffusiophoresis velocity is obtained for a droplet of a concentrated solution suspended in a binary gaseous mixture. The droplet is characterized by a high thermal conductivity. The droplet radius is assumed to be much greater than the mean free path for gaseous-mixture molecules. One of the gaseous-mixture molecular components is the vapor of the droplet solvent. According to the formula obtained in this study, the droplet is driven toward lower concentration of the volatile gaseous-mixture component by diffusive slip and in the opposite direction by phase transition. An increase in the relative mass concentration of the volatile solvent in the droplet enhances effects associated with the dependence of surface tension on the volatile-component concentration and the reactive transport due to the surface nonuniformity of phase transition. As the relative mass concentration of the volatile solvent in the droplet approaches unity, the effect of diffusive slip tends to vanish.     

The raspberry model for protein-like particles: Ellipsoids and confinement in cylindrical pores

The study of protein mass transport via atomistic simulation requires time and length scales beyond the computational capabilities of modern computer systems. The raspberry model for colloidal particles in combination with the mesoscopic hydrodynamic method of lattice Boltzmann facilitates coarse-grained simulations that are on the order of microseconds and hundreds of nanometers for the study of diffusive transport of protein-like colloid particles. The raspberry model reproduces linearity in resistance to motion versus particle size and correct enhanced drag within cylindrical pores at off-center coordinates for spherical particles. Owing to the high aspect ratio of many proteins, ellipsoidal raspberry colloid particles were constructed and reproduced the geometric resistance factors of Perrin and of Happel and Brenner in the laboratory-frame and in the moving body-frame. Accurate body-frame rotations during diffusive motion have been captured for the first time using projections of displacements. The spatial discretization of the fluid leads to a renormalization of the hydrodynamic radius, however, the data describes a self-consistent hydrodynamic frame within this renormalized system.     

旋涡的场特征与物质性──离散涡方法基础

旋涡的场特征与物质性──离散涡方法基础吴文权（华东工业大学，上海，２０００９３）关键词旋涡动力学，旋涡场性，旋涡物质性众人关注的不稳定、非定常、分离流动无不与旋涡有关。因而对旋涡的研究吸引了人们的注意。尤其是近来非线性科学的发展，提供了新的思想、方法...     

Capillary filling with pseudo-potential binary Lattice-Boltzmann model

We present a systematic study of capillary filling for a binary fluid by using a mesoscopic lattice Boltzmann model for immiscible fluids describing a diffusive interface moving at a given contact angle with respect to the walls. The phenomenological way to impose a given contact angle is analysed. Particular attention is given to the case of complete wetting, that is contact angle equal to zero. Numerical results yield quantitative agreement with the theoretical Washburn’s law, provided that the correct ratio of the dynamic viscosities between the two fluids is used. Finally, the presence of precursor films is experienced and it is shown that these films advance in time with a square-root law but with a different prefactor with respect to the bulk interface.     

基于范德瓦尔斯表面张力模式液滴撞击疏水壁面过程的研究

液滴撞击疏水壁面过程的研究在介观流体力学和微流体作用材料科学的研究中具有重要的理论意义和工程价值. 论文在SPH方法中引入范德瓦尔斯状态方程处理液滴表面张力, 考虑流体粒子之间远程吸引, 近程排斥的内部作用力, 提出了流体粒子与疏水壁面粒子间势能函数与表面张力相结合的作用模式. 通过模拟真空条件下两个静止的等体积液滴相互融合的过程, 验证了计算模式在模拟液滴的表面张力中的有效性. 采用该模式模拟的液滴撞击疏水壁面过程, 不仅能够有效地模拟液滴撞击壁面后的变形过程, 而且清晰地模拟出液滴的回弹、腾空以及二次撞壁现象的完整过程. 模拟结果与液滴撞击疏水壁面的实验结果以及VOF模拟结果符合较好, 表明本文所提出的表面张力和疏水壁面作用力处理模式对模拟液滴撞壁过程具有实际应用价值.     

微管道内两相流数值算法及在电浸润液滴控制中的应用

叙述了一种模拟电介质电润湿(electrowetting on dielectric，EWOD)下的微液滴的运动的数值方法. 采用二阶投影法求解N-S方程和level set 函数，并利用零level set函数俘获液滴运动界面，在液体与固体接触的边界上，通过引入动态接触角表征电介质表面润湿性随电势的改变. 数值计算基于MAC网格，模拟了2维微管道内与固体壁面接触的变润湿性的两种液体的分界面形状、平板上的微液滴在不同电势作用下处于不同湿润性的形态，以及微管道内改变接触角液滴的运动变形过程等算例. 关键词： 电浸润 接触角 level set函数 投影法     

Formation and mobility of droplets on composite layered substrates

A mesoscale fluid film placed on a solid support may break up and form droplets. In addition, droplets may exhibit spontaneous translation by modifying the wetting properties of the substrate, resulting in asymmetry in the contact angles. We examine mechanisms for droplet formation and motion on uniform and terraced landscapes, i.e., composite substrates. The fluid film stability, droplet formation and velocity are studied theoretically in the isothermal case using a lubrication approach in one spatial dimension. The droplet properties are found to involve contributions from both the terraced layer thickness and molecular interactions via the disjoining potential.     

Thermodynamically stable pickering emulsions

We show that under appropriate conditions, mixtures of oil, water, and nanoparticles form thermodynamically stable oil-in-water emulsions with monodisperse droplet diameters in the range of 30-150 nm. This observation challenges current wisdom that so-called Pickering emulsions are at most metastable and points to a new class of mesoscopic equilibrium structures. Thermodynamic stability is demonstrated by the spontaneous evolution of binary droplet mixtures towards one intermediate size distribution. Equilibrium interfacial curvature due to an asymmetric charge distribution induced by adsorbed colloids explains the growth of emulsion droplets upon salt addition. Moreover, the existence of a minimal radius of curvature with a concomitant expulsion of excess oil is in close analogy with microemulsions.     

The quantum Kirchhoff equation and quantum current and energy spectrum of a homogeneous mesoscopic dissipation transmission line

On the basis of quantization of charge, the loop equations of quantum circuits are investigated by using the Heisenberg motion equation for a mesoscopic dissipation transmission line. On the supposition that the system has a symmetry under translation in charge space, the quantum current and the quantum energy spectrum in the mesoscopic transmission line are given by solving their eigenvalue equations. Results show that the quantum current and the quantum energy spectrum are not only related to the parameters of the transmission line, but also dependent on the quantized character of the charge obviously.     

液滴溅落问题的光滑粒子动力学模拟

对传统的光滑粒子动力学方法进行了改进, 改进的光滑粒子动力学方法对传统粒子方法中的核近似式和粒子近似式进行了修正, 采用Riemann 算法求解光滑粒子动力学流体控制方程, 添加了表面张力的计算程序, 考虑了表面张力对液滴溅落的影响. 应用改进的光滑粒子动力学方法对液滴静止状态下冲击液面的飞溅过程进行了数值模拟. 计算结果表明, 改进的光滑粒子动力学方法能够有效地描述液滴溅落液面的动力学特性和自由表面变化特征, 能够得到稳定精度的结果.     

Viscosity and Diffusion Effects at the Boundary Surface of Viscous Fluid and Thermoelastic Diffusive Solid Medium

This paper concentrates on the wave motion at the interface of viscous compressible fluid half-space and homogeneous isotropic, generalized thermoelastic diffusive half-space. The wave solutions in both the fluid and thermoelastic diffusive half-spaces have been investigated; and the complex dispersion equation of leaky Rayleigh wave motion have been derived. The phase velocity and attenuation coefficient of leaky Rayleigh waves have been computed from the complex dispersion equation by using the Muller's method. The amplitudes of displacements, temperature change and concentration have been obtained. The effects of viscosity and diffusion on phase velocity and attenuation coefficient of leaky Rayleigh waves motion for different theories of thermoelastic diffusion have been depicted graphically. The magnitude of heat and mass diffusion flux vectors for different theories of thermoelastic diffusion have also been computed and represented graphically.     

Dynamical Evolution and Leading Order Gravitational Wave Emission of Riemann-S Binaries

An approximate strategy for studying the evolution of binary systems of extended objects is introduced. The stars are assumed to be polytropic ellipsoids. The surfaces of constant density maintain their ellipsoidal shape during the time evolution. The equations of hydrodynamics then reduce to a system of ordinary differential equations for the internal velocities, the principal axes of the stars and the orbital parameters. The equations of motion are given within Lagrangian and Hamiltonian formalism. The special case when both stars are axially symmetric fluid configurations is considered. Leading order gravitational radiation reaction is incorporated, where the quasi-static approximation is applied to the internal degrees of freedom of the stars. The influence of the stellar parameters, in particular the influence of the polytropic index n, on the leading order gravitational waveforms is studied.     

Linear stability analysis on the onset of Soret-driven motion in a nanoparticles suspension

The onset of Soret-driven instability in binary mixture heated from above is analysed using the linear stability theory. The horizontal fluid layer placed between two plates is in a thermally stable state but the Soret diffusion can induce buoyancy-driven convection in the case of a negative Soret coefficient. It is well known that convective motion sets in from both boundaries if the Rayleigh number exceeds its critical value. For the case of highly unstable density stratification the buoyancy-driven motion sets in during the transient diffusion stage. The new stability equations are derived and are solved analytically and numerically. Here the stability limits which are related to the onset time of instabilities and wave number are presented as a function of the Rayleigh number, Lewis number and the separation ratio. The present stability criteria are compared with the existing experimental data.