Direct Numerical Simulation Of Three-Dimensional Richtmyer-Meshkov Instability

Direct numerical simulation (DNS) is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interfacial instability is known as the Richtmyer-Meshkov (R-M) instability. The compressible Navier-Stoke equations are discretized with group velocity control (GVC) modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interfacial instability on produced vorticity and mixing is discussed.

Computer simulation of influence of sedimentation on rapid coagulation

A computer simulation was performed to explore the features and effects of sedimentation on rapid coagulation. To estimate the accumulated influence of gravity on coagulation for dispersions, a sedimentation influence ratio is defined. Some factors possibly related to the influence of sedimentation were considered in the simulation and analysed by comparing the size distribution of aggregates, the change in collision number, and coagulation rates at different gravity levels (0 g, 1 g and more with g being the gravitational constant).

Surface-induced melting of metal nanoclusters

We investigate the size effect on melting of metal nanoclusters by molecular dynamics simulation and thermo dynamic theory based on Kofman's melt model. By the minimization of the free energy of metal nanoclusters with respect to the thickness of the surface liquid layer, it has been found that the nanoclusters of the same metal have the same premelting temperature T-pre = T-0 - T-0(gamma(su) - gamma(lv) - gamma(sl))/(rhoLxi) (T-0 is the melting point of bulk metal, gamma(sv) the solid-vapour interfacial free energy, gamma(sl) the liquid-vapour interfacial free energy, gamma(sl),l the solid-liquid interfacial free energy, p the density of metal, L the latent heat of bulk metal, and xi the characteristic length of surface-interface interaction) to be independent of the size of nanoclusters, so that the characteristic length of a metal can be obtained easily by T-pre, which can be obtained by experiments or molecular dynamics (MD) simulations. The premelting temperature T-pre of Cu is obtained by AID simulations, then xi is obtained. The melting point T-cm is further predicted by free energy analysis and is in good agreement with the result of our MD simulations. We also predict the maximum premelting-liquid width of Cu nanoclusters with various sizes and the critical size, below which there is no premelting.

Examination of the LBM in simulation of microchannel flow in transitional regime

The gas flows in micro-electro-mechanical systems possess relatively large Knudsen number and usually belong to the slip flow and transitional flow regimes. Recently the lattice Boltzmann method (LBM) was proposed by Nie et al. in Journal of Statistical Physics, vol. 107, pp. 279-289, in 2002 to simulate the microchannel and microcavity flows in the transitional flow regime. The present article intends to test the feasibility of doing so. The results of using the lattice Boltzmann method and the direct simulation Monte Carlo method show good agreement between them for small Kn (Kn = 0.0194), poor agreement for Kn = 0.194, and large deviation for Kn = 0.388 in simulating microchannel flows. This suggests that the present version of the lattice Boltzmann method is not feasible to simulate the transitional channel flow.

On the computation of space-time correlations by large-eddy simulation

The effect of subgrid-scale (SGS) modeling on velocity (space-) time correlations is investigated in decaying isotropic turbulence. The performance of several SGS models is evaluated, which shows superiority of the dynamic Smagorinsky model used in conjunction with the multiscale large-eddy simulation (LES) procedure. Compared to the results of direct numerical simulation, LES is shown to underpredict the (un-normalized) correlation magnitude and slightly overpredict the decorrelation time scales. This can lead to inaccurate solutions in applications such as aeroacoustics. The underprediction of correlation functions is particularly severe for higher wavenumber modes which are swept by the most energetic modes. The classic sweeping hypothesis for stationary turbulence is generalized for decaying turbulence and used to analyze the observed discrepancies. Based on this analysis, the time correlations are determined by the wavenumber energy spectra and the sweeping velocity, which is the square root of the total energy. Hence, an accurate prediction of the instantaneous energy spectra is most critical to the accurate computation of time correlations. (C) 2004 American Institute of Physics.

Numerical simulation of Richtmyer-Meshkov instability

The compressible Navier-Stokes equations discretized with a fourth order accurate compact finite difference scheme with group velocity control are used to simulate the Richtmyer-Meshkov (R-M) instability problem produced by cylindrical shock-cylindrical material interface with shock Mach number Ms = 1.2 and density ratio 1:20 (interior density/outer density). Effect of shock refraction, reflection, interaction of the reflected shock with the material interface, and effect of initial perturbation modes on R-M instability are investigated numerically. It is noted that the shock refraction is a main physical mechanism of the initial phase changing of the material surface. The multiple interactions of the reflected shock from the origin with the interface and the R-M instability near the material interface are the reason for formation of the spike-bubble structures. Different viscosities lead to different spike-bubble structure characteristics. The vortex pairing phenomenon is found in the initial double mode simulation. The mode interaction is the main factor of small structures production near the interface.

Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

Numerical simulation of drop Marangoni migration under microgravity

Thermocapillary motion of a drop in a uniform temperature gradient is investigated numerically. The three-dimensional incompressible Navier-Stokes and energy equations are solved by the finite-element method. The front tracking technique is employed to describe the drop interface. To simplify the calculation, the drop shape is assumed to be a sphere. It has been verified that the assumption is reasonable under the microgravity environment. Some calculations have been performed to deal with the thermocapillary motion for the drops of different sizes. It has been verified that the calculated results are in good agreement with available experimental and numerical results. (C) 2003 Elsevier Ltd. All rights reserved.

Molecular dynamics simulation of plastic deformation of nanotwinned copper

The plastic deformation of polycrystalline Cu with ultrathin lamella twins has been studied using molecular dynamics simulations. The results of uniaxial tensile deformation simulation show that the abundance of twin boundaries provides obstacles to dislocation motion, which in consequence leads to a high strain hardening rate in the nanotwinned Cu. We also show that the twin lamellar spacing plays a vital role in controlling the strengthening effects, i.e., the thinner the thickness of the twin lamella, the harder the material. Additionally, twin boundaries can act as dislocation nucleation sites as they gradually lose coherency at large strain. These results indicate that controlled introduction of nanosized twins into metals can be an effective way of improving strength without suppression tensile ductility. (C) 2007 American Institute of Physics.

Atomistic investigation of the effects of temperature and surface roughness on diffusion bonding between Cu and Al

Molecular dynamics (MD) simulations are carried out to analyze the diffusion bonding at Cu/Al interfaces. The results indicate that the thickness of the interfacial layer is temperature-dependent, with higher temperatures yielding larger thicknesses. At temperatures below 750 K, the interface thickness is found to increase in a stepwise manner as a function of time. At temperatures above 750 K, the thickness increases rapidly and smoothly. When surface roughness is present, the bonding process consists of three stages. In the first stage, surfaces deform under stress, resulting in increased contact areas. The second stage involves significant plastic deformation at the interface as temperature increases, resulting in the disappearance of interstices and full contact of the surface pair. The last stage entails the diffusion of atoms under constant temperature. The bonded specimens show tensile strengths reaching 88% of the ideal Cu/Al contact strength. (c) 2007 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Computer simulation on the collision-sticking dynamics of two colloidal particles in an optical trap

Collisions of a particle pair induced by optical tweezers have been employed to study colloidal stability. In order to deepen insights regarding the collision-sticking dynamics of a particle pair in the optical trap that were observed in experimental approaches at the particle level, the authors carry out a Brownian dynamics simulation. In the simulation, various contributing factors, including the Derjaguin-Landau-Verwey-Overbeek interaction of particles, hydrodynamic interactions, optical trapping forces on the two particles, and the Brownian motion, were all taken into account. The simulation reproduces the tendencies of the accumulated sticking probability during the trapping duration for the trapped particle pair described in our previous study and provides an explanation for why the two entangled particles in the trap experience two different statuses. (c) 2007 American Institute of Physics.