Reconstructing The Vertical Distribution Of The Aeolian Saltation Mass Flux Based On The Probability Distribution Of Lift-Off Velocity

The probability distribution of lift-off velocity of the saltating grains is a bridge to linking microscopic and macroscopic research of aeolian sand transport. The lift-off parameters of saltating grains (i.e., the horizontal and vertical lift-off velocities, resultant lift-off velocity, and lift-off angle) in a wind tunnel are measured by using a Phase Doppler Particle Analyzer (PDPA). The experimental results show that the probability distribution of horizontal lift-off velocity of saltating particles on a bed surface is a normal function, and that of vertical lift-off velocity is an exponential function. The probability distribution of resultant lift-off velocity of saltating grains can be expressed as a log-normal function, and that of lift-off angle complies with an exponential function. A numerical model for the vertical distribution of aeolian mass flux based on the probability distribution of lift-off velocity is established. The simulation gives a sand mass flux distribution which is consistent with the field data of Namikas (Namikas, S.L., 2003. Field measurement and numerical modelling of acolian mass flux distributions on a sandy beach, Sedimentology 50, 303-326). Therefore, these findings are helpful to further understand the probability characteristics of lift-off grains in aeolian sand transport. (c) 2007 Elsevier B.V. All rights reserved.

Thermal fatigue on pistons induced by shaped high power laser. Part II: Design of spatial intensity distribution via numerical simulation

In the laser induced thermal fatigue simulation test on pistons, the high power laser was transformed from the incident Gaussian beam into a concentric multi-circular pattern with specific intensity ratio. The spatial intensity distribution of the shaped beam, which determines the temperature field in the piston, must be designed before a diffractive optical element (DOE) can be manufactured. In this paper, a reverse method based on finite element model (FEM) was proposed to design the intensity distribution in order to simulate the thermal loadings on pistons. Temperature fields were obtained by solving a transient three-dimensional heat conduction equation with convective boundary conditions at the surfaces of the piston workpiece. The numerical model then was validated by approaching the computational results to the experimental data. During the process, some important parameters including laser absorptivity, convective heat transfer coefficient, thermal conductivity and Biot number were also validated. Then, optimization procedure was processed to find favorable spatial intensity distribution for the shaped beam, with the aid of the validated FEM. The analysis shows that the reverse method incorporated with numerical simulation can reduce design cycle and design expense efficiently. This method can serve as a kind of virtual experimental vehicle as well, which makes the thermal fatigue simulation test more controllable and predictable. (C) 2007 Elsevier Ltd. All rights reserved.

Interface Layer Effect On The Stress Distribution Of A Wafer-Bonded Bilayer Structure

The interface layer plays an important role in stress transfer in composite structures. However, many interface layer properties such as the modulus, thickness, and uniformity are difficult to determine. The model developed in this article links the influence of the interface layer on the normal stress distribution along the layer thickness with the layer surface morphology before bonding. By doing so, a new method of determining the interfacial parameter(s) is suggested. The effects of the layer thickness and the surface roughness before bonding on the normal stress distribution and its depth profile are also discussed. For ideal interface case with no interfacial shear stress, the normal stress distribution pattern can only be monotonically decreased from the interface. Due to the presence of interfacial shear stress, the normal stress distribution is much more complex, and varies dramatically with changes in the properties of the interface layer, or the dimensions of the bonding layers. The consequence of this dramatic stress field change, such as the shift of the maximum stress from the interface is also addressed. The size-dependent stress distribution in the thickness direction due to the interface layer effect is presented. When the interfacial shear stress is reduced to zero, the model presented in this article is also demonstrated to have the same normal stress distribution as obtained by the previous model, which does not consider the interface layer effect.

Regularity of strain distribution in short-fiber/whisker reinforced composites

Based on studies on the strain distribution in short-fiber/whisker reinforced metal matrix composites, a deformation characteristic parameter, lambda is defined as a ratio of root-mean-square strain of the reinforcers identically oriented to the macro-linear strain along the same direction. Quantitative relation between lambda and microstructure parameters of composites is obtained. By using lambda, the stiffness moduli of composites with arbitrary reinforcer orientation density function and under arbitrary loading condition are derived. The upper-bound and lower-bound of the present prediction are the same as those from the equal-strain theory and equal-stress theory, respectively. The present theory provides a physical explanation and theoretical base for the present commonly-used empirical formulae. Compared with the microscopic mechanical theories, the present theory is competent for stiffness modulus prediction of practical engineering composites in accuracy and simplicity.

Fracture statistics of brittle materials: Weibull or normal distribution

The fit of fracture strength data of brittle materials (Si3N4, SiC, and ZnO) to the Weibull and normal distributions is compared in terms of the Akaike information criterion. For Si3N4, the Weibull distribution fits the data better than the normal distribution, but for ZnO the result is just the opposite. In the case of SiC, the difference is not large enough to make a clear distinction between the two distributions. There is not sufficient evidence to show that the Weibull distribution is always preferred to other distributions, and the uncritical use of the Weibull distribution for strength data is questioned.

First-passage time of Duffing oscillator under combined harmonic and white-noise excitations

The first-passage time of Duffing oscillator under combined harmonic and white-noise excitations is studied. The equation of motion of the system is first reduced to a set of averaged Ito stochastic differential equations by using the stochastic averaging method. Then, a backward Kolmogorov equation governing the conditional reliability function and a set of generalized Pontryagin equations governing the conditional moments of first-passage time are established. Finally, the conditional reliability function, and the conditional probability density and moments of first-passage time are obtained by solving the backward Kolmogorov equation and generalized Pontryagin equations with suitable initial and boundary conditions. Numerical results for two resonant cases with several sets of parameter values are obtained and the analytical results are verified by using those from digital simulation.

Optimal bounded control of first-passage failure of strongly non-linear oscillators under combined harmonic and white-noise excitations

A procedure for designing the optimal bounded control of strongly non-linear oscillators under combined harmonic and white-noise excitations for minimizing their first-passage failure is proposed. First, a stochastic averaging method for strongly non-linear oscillators under combined harmonic and white-noise excitations using generalized harmonic functions is introduced. Then, the dynamical programming equations and their boundary and final time conditions for the control problems of maximizing reliability and of maximizing mean first-passage time are formulated from the averaged Ito equations by using the dynamical programming principle. The optimal control law is derived from the dynamical programming equations and control constraint. Finally, the conditional reliability function, the conditional probability density and mean of the first-passage time of the optimally controlled system are obtained from solving the backward Kolmogorov equation and Pontryagin equation. An example is given to illustrate the proposed procedure and the results obtained are verified by using those from digital simulation. (C) 2003 Elsevier Ltd. All rights reserved.

NUMERICAL ANALYSIS OF THE FLOWFIELD IN A SUPERSONIC COIL WITH AN INTERLEAVED JET CONFIGURATION AND ITS EFFECT ON THE GAIN DISTRIBUTION

In a supersonic chemical oxygen-iodine laser (COIL) operating without primary buffer gas, the features of flowfield have significant effects on the Laser efficiency and beam quality. In this paper three-dimensional, multi-species, chemically reactive CFD technology was used to study the flowfield in mixing nozzle implemented with a supersonic interleaving jet configuration. The features of the flowfield as well as its effect on the spatial distribution of small signal gain were analyzed.

Numerical Study on Flow Field and Temperature Distribution in Growth Process of 200 mm Czochralski Silicon Crystals

The melt flow and temperature distribution in a 200 mm silicon Czochralski furnace with a cusp magnetic field was modeled and simulated by using a finite-volume based FLUTRAPP ( Fluid Flow and Transport Phenomena Program) code. The melt flow in the crucible was focused, which is a result of the competition of buoyancy, the centrifugal forces caused by the rotations of the crucible and crystal, the thermocapillary force on the free surfaces and the Lorentz force induced by the cusp magnetic field. The zonal method for radiative heat transfer was used in the growth chamber, which was confined by the crystal surface, melt surface, crucible, heat shield, and pull chamber. It was found that the cusp magnetic field could strength the dominant counter-rotating swirling flow cell in the crucible and reduce the flow oscillation and the pulling-rate fluctuation. The fluctuation of dopant and oxygen concentration in the growing crystal could thus be smoothed.

Experimental investigation on bubble coalescence under nonuniform temperature distribution in reduced gravity

Results on bubble coalescences from the space experiment of thermocapillary bubble migration conducted on board the Chinese 22nd recoverable satellite are presented in this paper. Some coalescences of large spherical bubbles under microgravity are observed through bubbles staying at the upper side of the test cell. The data of bubble coalescence time are recorded and compared with theoretical predictions, which is based on a theory to describe the tendency of coalescence connected to chemical potential difference. It is implied that the theory is applicable for the experimental data of bubble coalescence. Moreover, the angle between the line of two bubble centers and temperature gradient falled mostly in the range 20 degrees-40 degrees. (C) 2007 Elsevier Inc. All rights reserved.