Thermocapillary Migration Of Nondeformable Drops

In this paper, we present a numerical study on the thermocapillary migration of drops. The Navier-Stokes equations coupled with the energy conservation equation are solved by the finite-difference front-tracking scheme. The axisymmetric model is adopted in Our simulations, and the drops are assumed to be perfectly spherical and nondeformable. The benchmark simulation starts from the classical initial condition with a uniform temperature gradient. The detailed discussions and physical explanations of migration phenomena are presented for the different values of (1) the Marangoni numbers and Reynolds numbers of continuous phases and drops and (2) the ratios of drop densities and specific heats to those of continuous phases. It is found that fairly large Marangoni numbers may lead to fluctuations in drop velocities at the beginning part of simulations. Finally, we also discuss the influence of initial conditions on the thermocapillary migrations. (C) 2008 American Institute of Physics.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

Singularity of dynamic stress fields around an interface crack between viscoelastic bodies

The singular nature of the dynamic stress fields around an interface crack located between two dissimilar isotropic linearly viscoelastic bodies is studied. A harmonic load is imposed on the surfaces of the interface crack. The dynamic stress fields around the crack are obtained by solving a set of simultaneous singular integral equations in terms of the normal and tangent crack dislocation densities. The singularity of the dynamic stress fields near the crack tips is embodied in the fundamental solutions of the singular integral equations. The investigation of the fundamental solutions indicates that the singularity and oscillation indices of the stress fields are both dependent upon the material constants and the frequency of the harmonic load. This observation is different from the well-known -1/2 oscillating singularity for elastic bi-materials. The explanation for the differences between viscoelastic and elastic bi-materials can be given by the additional viscosity mismatch in the case of viscoelastic bi-materials. As an example, the standard linear solid model of a viscoelastic material is used. The effects of the frequency and the material constants (short-term modulus, long-term modulus and relaxation time) on the singularity and the oscillation indices are studied numerically.

Impact of cathode evaporation on a free-burning arc

A free-burning, high-intensity argon arc at atmospheric pressure was modelled during the evaporation of copper from the cathode. The effect of cathode evaporation on the temperature, mass flow, current flow and Cu concentration was studied for the entire plasma region. The copper evaporates from the tip of the cathode with an evaporation rate of 1 mg s-1. The copper vapour in the cathode region has a velocity of 210 m s-1 with a mass concentration of above 90% within 0.5 mm from the arc axis. The vapour passes from the cathode toward the anode with a slight diffusion in the argon plasma. Higher temperatures and current densities were calculated in the core of the arc caused by the cathode evaporation.

Phase relaxation model for 2-phase flows

A phase relaxation model (PRM) for 2-phase flows is presented in this paper on the basis of three principal assumptions. The basic equations for PRM arc derived from the Boltzmann equations for gas-partlcle mixture, The general characteristics and solving process of the PRM's basic equations are also presented and discussed. Many terms in the PRM's basic equations contain a factor ε= ρgρp/ρg+ρp2 which is an intrinsic small parameter for 2-phase mixture, with ρg and ρp being respectively the densities of gas and particle phases.This makes it possible to simplify the computation of the PRM's basic equations. The model is applied to for example, studying file steady propagation of shock waves in gas-particle mixture. The analysis shows that with an increase of shock wave strength the relaxation process behind a gasdynamics shock front becomes a kind of dynamics relaxation instead of the standard exponential relaxation process. A method of determining experimentally the velocity and tem...更多perature relaxation rates (or times) of gas-particle flows is suggested and analyzed.

Measurement of Electron Density Profile Behind Strong Shock Waves with a Langmuir Probe

At the shock velocity range of 7~9km/s, the variations of electron density behind strong normal shock waves are measured in a low-density shock tube by using the Langmuir electrostatic probe technique. The electron temperature, calculated based on Park’s three-temperature model, is used in interpreting the probe current data. The peak electron densities determined in the present experiment are shown to be in a good agreement with those predicted by Lin’s calculation. The experimentally obtained ratios of the characteristic ionization distance to the mean free path of freestream ahead of the shock wave are found to be in a good agreement with the existing experiments and Park’s calculation.

Mechanics of metallic foams: A fractal approach

A fractal approach was proposed to investigate the meso structures and size effect of metallic foams: For a series At foams of different relative densities, the information dimension method was applied to measure meso structures. The generalized sierpinski carpet was introduced to map the meso structures of the foam according to specific dimension. The results show that the fractal-based model can not only reveal the variation of yield strength with specimen size, but also bridge the meso structures and mechanical proper-ties of Al foams directly. Key words: metallic foams; fractal; size effect; meso structures

Characteristics of flow fields induced by interfacial waves in two-layer fluid

When designing deep ocean structures, it is necessary to estimate the effects of internal waves on the platform and auxiliary parts such as tension leg, riser and mooring lines. Up to now, only a few studies are concerned with the internal wave velocity fields. By using the most representative two-layer model, we have analyzed the behavior of velocity field induced by interfacial wave in the present paper. We find that there may exist velocity shear of fluid particles in the upper and lower layers so that any structures in the ocean are subjected to shear force nearby the interface. In the meantime, the magnitude of velocity for long internal wave appears spatially uniform in the respective layer although they still decay exponentially. Finally, the temporal variation for Stokes and solitary waves are shown to be of periodical and pulse type.

The diagnostics of density distribution for inhomogeneous dense DT plasmas using fast protons

The density distribution of inhomogeneous dense deuterium-tritium plasmas in laser fusion is revealed by the energy loss of fast protons going through the plasma. In our simulation of a plasma density diagnostics, the fast protons used for the diagnostics may be generated in the laser-plasma interaction. Dividing a two-dimensional area into grids and knowing the initial and final energies of the protons, we can obtain a large linear and ill-posed equation set. for the densities of all grids, which is solved with the Tikhonov regularization method. We find that the accuracy of the set plan with four proton sources is better than those of the set plans with less than four proton sources. Also we have done the density reconstruction especially. for four proton sources with and without assuming circularly symmetrical density distribution, and find that the accuracy is better for the reconstruction assuming circular symmetry. The error is about 9% when no noise is added to the final energy for the reconstruction of four proton sources assuming circular symmetry. The accuracies for different random noises to final proton energies with four proton sources are also calculated.

Neutron source from thin foil confined by two laser pulses

Neutron production from a thin deuterium-tritium (D-T) foil irradiated by two intense femtosecond laser pulses from opposite sides with zero phase difference is studied analytically and numerically. For the interaction of a laser pulse of amplitude a = 7, focal area 100 mu m(2) and areal density 4.4 x 10(18) cm(-2) with a D-T plasma foil, about 1.17 x 10(21) neutron s(-1) can be obtained, much more than from other methods. The profiles of the ion and electron densities are also calculated.

The propagation and deposition of fast muons in dense DT mixture

The propagation of the fast muon population mainly due to collisional effect in a dense deuterium-tritium (DT for short) mixture is investigated and analysed within the framework of the relativistic Fokker-Planck equation. Without the approximation that the muons propagate straightly in the DT mixture, the muon penetration length, the straggling length, and the mean transverse dispersion radius are calculated for different initial energies, and especially for different densities of the densely compressed DT mixture in our suggested muon-driven fast ignition (FI). Unlike laser-driven FI requiring super-high temperature, muons can catalyze DT fusion at lower temperatures and may generate an ignition sparkle before the self-heating fusion follows. Our calculation is important for the feasibility and the experimental study of muon-driven FI.

Reflection and transmission of Gaussian beam from a uniaxial crystal slab

We investigate the characteristics of Gaussian beams reflected and transmitted from a uniaxial crystal slab with an arbitrary orientation of its optical axis. The formulas of the total electric and magnetic fields inside and outside the slab are derived by use of Maxwell's equations and by matching the boundary conditions at the interfaces. Numerical simulations are presented and the field values as well as the power densities are computed. Negative refractions are demonstrated when the beam is transmitted through a uniaxial crystal slab. Beam splitting of the reflected beam is observed and is explained by the resonant transmission for plane waves. Dependences of the lateral shift on the incident angle and beam width are discussed. Negative and positive lateral shifts are observed due to the spatial anisotropic properties.

Prescription for optimizing holograms in LiNbO3 : Fe : Mn

The photorefractive holographic dynamics of grating formation in photochromic doubly doped LiNbO3:Fe:Mn crystal is studied numerically and analytically in terms of the two-center model of Kukhtarev Et al. [Ferroelectrics 22, 949 (1979)]. The relations among the recorded and fixed space-charge fields and the doping densities, the oxidation-reduction states of the fields, and the intensities of UV-sensitizing and red recording beams are studied. Important conditions and effects are feued, and an optimal prescription for material doping and oxidation-reduction processing is suggested in which the crystal can be strongly oxidized and the Mn-doping density is smaller than the Fe-doping density. (C) 2000 Optical Society of America. OCIS codes: 050.7330, 190.5330, 090.2900.

Optimal recording wavelength for nonvolatile photorefractive holograms in Ce : Cu : LiNbO

Three wavelengths of red, green and blue of recording beams are systemically tested for the UV-assistant recording and optical fixing of holograms in a strongly oxidized Ce:Cu:LiNbO3 crystal. Three different photorefractive phenomena are observed. It is shown that the green beams will optimally generate a critical strong nonvolatile hologram with quick sensitivity and the optimal switching technique could be jointly used to obtain a nearly 100% high diffraction. Theoretical verification is given, and a prescription on the doping densities and on the oxidation/reduction states of the material to match a defined recording wavelength for high diffraction is suggested.

脉冲激光晶化非晶硅薄膜的有限差分模拟

根据热传导原理，建立了脉冲激光晶化非晶硅薄膜的理论模型。运用有限差分方法研究了不同激光波长、能量密度等因素对薄膜温度变化及相变过程的影响。计算了不同波长激光器对厚度500nm非晶硅晶化的阈值能量密度。结果发现，准分子晶化的阈值能量密度最低，但是在同样的能量密度下，熔融深度却不及使用更长波长的激光器。计算并分析了升高衬底温度对结晶速度和晶粒尺寸的影响，模拟结果较好地验证了实验结论和规律。