Simulations of impinging droplets with surfactant-dependent dynamic contact angle

An arbitrary Lagrangian-Eulerian (ALE) finite element scheme for computations of soluble surfactant droplet impingement on a horizontal surface is presented. The numerical scheme solves the time-dependent Navier-Stokes equations for the fluid flow, scalar convection-diffusion equation for the surfactant transport in the bulk phase, and simultaneously, surface evolution equations for the surfactants on the free surface and on the liquid-solid interface. The effects of surfactants on the flow dynamics are included into the model through the surface tension and surfactant-dependent dynamic contact angle. In particular, the dynamic contact angle (theta(d)) of the droplet is defined as a function of the surfactant concentration at the contact line and the equilibrium contact angle (theta(0)(e)) of the clean surface using the nonlinear equation of state for surface tension. Further, the surface forces are included into the model as surface divergence of the surface stress tensor that allows to incorporate the Marangoni effects without calculating the surface gradient of the surfactant concentration on the free surface. In addition to a mesh convergence study and validation of the numerical results with experiments, the effects of adsorption and desorption surfactant coefficients on the flow dynamics in wetting, partially wetting and non-wetting droplets are studied in detail. It is observed that the effects of surfactants are more in wetting droplets than in the non-wetting droplets. Further, the presence of surfactants at the contact line reduces the equilibrium contact angle further when theta(0)(e) is less than 90 degrees, and increases it further when theta(0)(e) is greater than 90 degrees. Nevertheless, the presence of surfactants has no effect on the contact angle when theta(0)(e) = 90 degrees. The numerical study clearly demonstrates that the surfactant-dependent contact angle has to be considered, in addition to the Marangoni effect, in order to study the flow dynamics and the equilibrium states of surfactant droplet impingement accurately. The proposed numerical scheme guarantees the conservation of fluid mass and of the surfactant mass accurately. (C) 2015 Elsevier Inc. All rights reserved.

对流扩散方程的变步长摄动有限差分格式

摄动有限差分(PFD)方法是构造高精度差分格式的一种新方法。变步长摄动有限差分方法是等步长摄动有限差分方法的发展和推广。对需要局部加密网格的计算问题，变步长PFD格式不需要对自变量进行数学变换，且和等步长PFD格式一样，具有如下的共同特点：从变步长一阶迎风格式出发，通过把非微商项(对流系数和源项)作变步长摄动展开，展开幂级数系数通过消去摄动格式修正微分方程的截断误差项求出，由此获得高精度变步长PFD格式。该格式在一、二和三维情况下分别仅使用三、五和七个基点，且具有迎风性。文中利用变步长PFD格式对对流扩散反应模型方程，变系数方程及Burgers方程等进行了数值模拟，并与一阶迎风和二阶中心格式及其问题的精确解作了比较。数值试验表明，与一阶迎风和二阶中心格式相比，变步长PFD格式具有精度高，稳定性与收敛性好的特点。变步长PFD格式与等步长PFD格式相比，变步长PFD解在薄边界层型区域的分辨率得到了明显的提高。

Direct Numerical Simulation of Passive Scalar in Decaying Compressible Turbulence

The passive scalars in the decaying compressible turbulence with the initial Reynolds number (defined by Taylor scale and RMS velocity) Re=72, the initial turbulent Mach numbers (defined by RMS velocity and mean sound speed) Mt=0.2-0.9, and the Schmidt numbers of passive scalar Sc=2-10 are numerically simulated by using a 7th order upwind difference scheme and 8th order group velocity control scheme. The computed results are validated with different numerical methods and different mesh sizes. The Batchelor scaling with k(-1) range is found in scalar spectra. The passive scalar spectra decay faster with the increasing turbulent Mach number. The extended self-similarity (ESS) is found in the passive scalar of compressible turbulence.

Mechanism on Progressive Failure of a Faulted Rock Slope Due to Slip-Weakening

Slip-weakening is one of the characteristics of geological materials under certain loadings. Non-uniform rock structure may exist in the vicinity of the slip surface for a rock slope. Some portion of the slip surface may be penetrated but the other not. For the latter case, the crack or the fault surface will undergo shear deformation before it becomes a successive surface under a certain loading. As the slipped portion advances,slip-weakening occurs over a distance behind the crack tip. In the weakening zone, the shear strength will decrease from its peak value to residual friction level. The stress will redistribute along the surface of crack and in the weakening zone. Thus the changed local stress concentration leads the crack to extend and the ratio of penetration of the slip surface to increase. From the view of large-scale for the whole slip surface, the shear strength will decrease due to the damage of interior rock structure, and the faulted rock behaves as a softening material. Such a kind of mechanism performs in a large number of practical landslides in the zones experienced strong earthquakes. It should be noted that the mechanism mentioned above is different from that of the breakage of structural clay,in which the geological material is regarded as a medium containing structural lumps and structural bands. In this paper, the softening behavior of a faulted rock should be regarded as a comprehensive result of the whole complicated process including slip-weakening, redistribution of stress, extension of crack tip, and the penetration of the slip surface. This process is accompanied by progressive failure and abrupt structural damage. The size of slip-weakening zone is related to the undergoing strain. Once the relative slide is initiated (local or integrated), the effect of slip-weakening will behave in a certain length behind the crack tip until the formation of the whole slip surface.

Three-Dimensional Propagation of the Transmitted Shock Wave in a Square Cross-Sectional Chamber

An investigation into the three-dimensional propagation of the transmitted shock wave in a square cross-section chamber was described in this paper, and the work was carried out numerically by solving the Euler equations with a dispersion-controlled scheme. Computational images were constructed from the density distribution of the transmitted shock wave discharging from the open end of the square shock tube and compared directly with holographic interferograms available for CFD validation. Two cases of the transmitted shock wave propagating at different Mach numbers in the same geometry were simulated. A special shock reflection system near the corner of the square cross-section chamber was observed, consisting of four shock waves: the transmitted shock wave, two reflection shock waves and a Mach stem. A contact surface may appear in the four-shock system when the transmitted shock wave becomes stronger. Both the secondary shock wave and the primary vortex loop are three-dimensional in the present case due to the non-uniform flow expansion behind the transmitted shock.

Thermocapillary Effect in a Thermal Rheological Jet

The process of die swell in polymer jets is an important feature within polymer processing and can be explained through a study of its rheological effects. The existence of a thermocapillary effect, driven by the gradient of its surface tension, should be considered when examining a thermal jet that has a non-uniform temperature distribution on its free surface, as in various polymer processings. Both the rheological effect and thermocapillary effect on die swell can be studied numerically through a finite element method as used on a two-dimensional and unsteady model, in which a Coleman-Noll second-order fluid model is employed. The results show that the expanding angle depends on both the rheological property of the fluid and the pressure at the vessel exit. Although both the thermocapillary and the rheological effects contribute to the cross-section expansion of the fluid jet, the latter is more important in determining the expansion.

Optimized Group Velocity Control Scheme and DNS of Decaying Compressible Turbulence of Relative High Turbulent Mach Number

To overcome the difficulty in the DNS of compressible turbulence at high turbulent Mach number, a new difference scheme called GVC8 is developed. We have succeeded in the direct numerical simulation of decaying compressible turbulence up to turbulent Mach number 0.95. The statistical quantities thus obtained at lower turbulent Mach number agree well with those from previous authors with the same initial conditions, but they are limited to simulate at lower turbulent Mach numbers due to the so-called start-up problem. The energy spectrum and coherent structure of compressible turbulent flow are analysed. The scaling law of compressible turbulence is studied. The computed results indicate that the extended self-similarity holds in decaying compressible turbulence despite the occurrence of shocklets, and compressibility has little effects on relative scaling exponents when turbulent Mach number is not very high.

Plastic constitutive behavior of short-fiber/particle reinforced composites

A cylindrical cell model based on continuum theory for plastic constitutive behavior of short-fiber/particle reinforced composites is proposed. The composite is idealized as uniformly distributed periodic arrays of aligned cells, and each cell consists of a cylindrical inclusion surrounded by a plastically deforming matrix. In the analysis, the non-uniform deformation field of the cell is decomposed into the sum of the first order approximate field and the trial additional deformation field. The precise deformation field are determined based on the minimum strain energy principle. Systematic calculation results are presented for the influence of reinforcement volume fraction and shape on the overall mechanical behavior of the composites. The results are in good agreement with the existing finite element analyses and the experimental results. This paper attempts to stimulate the work to get the analytical constitutive relation of short-fiber/particle reinforced composites.

一个基于非均匀网格上的二阶基本无振荡差分格式

在基本无振荡格式的构造中，将通常的对流通量f的逼近方式推广到对通量导数的逼近，这一构造方法可以有效地应用到非均匀或非结构网格。直接基于非均匀网格上，构造了一个二阶的基本无振荡(ENO)差分格式。该格式具有形式简单，对网格的划分灵活，与传统格式相比不增加计算量等优点，儿个数值算例证明了格式的有效性。

Investigation on a simulation model of floating half zone convection .2. Experiment

A simulation model of floating half zone with non-uniform temperature distribution at the upper rod and uniform temperature distribution at lower rod was discussed by numerical investigation in a previous paper. In the present paper, the experimental investigation of the simulation model is given generally. The results of the present model show that the temperature profile is quite different and the critical applied temperature difference is lower than the one of usual model with same geometrical parameters in most cases. The features of critical Marangoni number depending on the liquid bridge volume are also different from the ones of usual model.

Triple-region structure for turbulent-flow in a square duct - a finite-element approach

The Reynolds-averaged Navier-Stokes equations for describing the turbulent flow in a straight square duct are formulated with two different turbulence models. The governing equations are then expanded as a multi-deck structure in a plane perpendicular to the streamwise direction, with each deck characterized by its dominant physical forces as commonly carried out in analytical work using triple-deck expansion. The resulting equations are numerically integrated using higher polynomial (H-P) finite element technique for each cross-sectional plane to be followed by finite difference representation in the streamwise direction until a fully developed state is reached. The computed results using the two different turbulence models show fair agreement with each other, and concur with the vast body of available experimental data. There is also general agreement between our results and the recent numerical works anisotropic k-epsilon turbulence model.

Behavior of ar plasma formed in a high-density plasma source - an ecr reactor

In order to develop the ultra-large scale integration(ULSI), low pressure and high density plasma apparatus are required for etching and deposit of thin films. To understand critical parameters such as the pressure, temperature, electrostatic potential and energy distribution of ions impacting on the wafer, it is necessary to understand how these parameters are influenced by the power input and neutral gas pressure. In the present work, a 2-D hybrid electron fluid-particle ion model has been developed to simulate one of the high density plasma sources-an Electron Cyclotron Resonance (ECR) plasma system with various pressures and power inputs in a non-uniform magnetic field. By means of numerical simulation, the energy distributions of argon ion impacting on the wafer are obtained and the plasma density, electron temperature and plasma electrostatic potential are plotted in 3-D. It is concluded that the plasma density depends mainly on both the power input and neutral gas pressure. However, the plasma potential and electron temperature can hardly be affected by the power input, they seem to be primarily dependent on the neutral gas pressure. The comparison shows that the simulation results are qualitatively in good agreement with the experiment measurements.

Domain decomposition hybrid method for numerical simulation of bluff body flows:theoretical model and application

The discrete vortex method is not capable of precisely predicting the bluff body flow separation and the fine structure of flow field in the vicinity of the body surface. In order to make a theoretical improvement over the method and to reduce the difficulty in finite-difference solution of N-S equations at high Reynolds number, in the present paper, we suggest a new numerical simulation model and a theoretical method for domain decomposition hybrid combination of finite-difference method and vortex method. Specifically, the full flow. field is decomposed into two domains. In the region of O(R) near the body surface (R is the characteristic dimension of body), we use the finite-difference method to solve the N-S equations and in the exterior domain, we take the Lagrange-Euler vortex method. The connection and coupling conditions for flow in the two domains are established. The specific numerical scheme of this theoretical model is given. As a preliminary application, some numerical simulations for flows at Re=100 and Re-1000 about a circular cylinder are made, and compared with the finite-difference solution of N-S equations for full flow field and experimental results, and the stability of the solution against the change of the interface between the two domains is examined. The results show that the method of the present paper has the advantage of finite-difference solution for N-S equations in precisely predicting the fine structure of flow field, as well as the advantage of vortex method in efficiently computing the global characteristics of the separated flow. It saves computer time and reduces the amount of computation, as compared with pure N-S equation solution. The present method can be used for numerical simulation of bluff body flow at high Reynolds number and would exhibit even greater merit in that case.

An asymptotic analysis for one dimensional stress wave propagation in damaged viscoelastic media

A regular perturbation technique is suggested to deal with the problem of one dimensional stress wave propagation in viscoelastic media with damage. Based upon the first order asymptotic solution obtained, the characteristics of wave attenuation are studied. In fact, there exist three different time-dependent phenomena featuring the dynamic response of the materials, the first expressing the characteristics of wave propagation, the second indicating the innate effect of visco-elastic matrix and the third coming from the time dependent damage. The comparision of first order asymptotic solution with the numerical results calculated by a finite difference procedure shows that the perturbation expansion technique may offer a useful approach to the problem concerned.

Compressible laminar boundary-layer flows of a dusty gas over a semi-infinite flat-plate

The compressible laminar boundary-layer flows of a dilute gas-particle mixture over a semi-infinite flat plate are investigated analytically. The governing equations are presented in a general form where more reasonable relations for the two-phase interaction and the gas viscosity are included. The detailed flow structures of the gas and particle phases are given in three distinct regions : the large-slip region near the leading edge, the moderate-slip region and the small-slip region far downstream. The asymptotic solutions for the two limiting regions are obtained by using a seriesexpansion method. The finite-difference solutions along the whole length of the plate are obtained by using implicit four-point and six-point schemes. The results from these two methods are compared and very good agreement is achieved. The characteristic quantities of the boundary layer are calculated and the effects on the flow produced by the particles are discussed. It is found that in the case of laminar boundary-layer flows, the skin friction and wall heat-transfer are higher and the displacement thickness is lower than in the pure-gas case alone. The results indicate that the Stokes-interaction relation is reasonable qualitatively but not correct quantitatively and a relevant non-Stokes relation of the interaction between the two phases should be specified when the particle Reynolds number is higher than unity.