Computer Simulation Of Random Sphere Packing In An Arbitrarily Shaped Container

Most simulations of random sphere packing concern a cubic or cylindric container with periodic boundary, containers of other shapes are rarely studied. In this paper, a new relaxation algorithm with pre-expanding procedure for random sphere packing in an arbitrarily shaped container is presented. Boundaries of the container are simulated by overlapping spheres which covers the boundary surface of the container. We find 0.4 similar to 0.6 of the overlap rate is a proper value for boundary spheres. The algorithm begins with a random distribution of small internal spheres. Then the expansion and relaxation procedures are performed alternately to increase the packing density. The pre-expanding procedure stops when the packing density of internal spheres reaches a preset value. Following the pre-expanding procedure, the relaxation and shrinking iterations are carried out alternately to reduce the overlaps of internal spheres. The pre-expanding procedure avoids the overflow problem and gives a uniform distribution of initial spheres. Efficiency of the algorithm is increased with the cubic cell background system and double link data structure. Examples show the packing results agree well with both computational and experimental results. Packing density about 0.63 is obtained by the algorithm for random sphere packing in containers of various shapes.

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.

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.

The strength distribution of brittle materials with a high-concentration of cavities

The maximum stress concentration factor in brittle materials with a high concentration of cavities is obtained. The interaction between the nearest cavities, in addition to the far field interactions, is taken into account to evaluate the strength distribution based on the statistical analysis of the nearest distance distribution. Through this investigation, it is found that the interaction between the nearest neighbors is much more important than the far field interactions, and one has to consider it in calculating the strength of brittle materials even if the volume fraction of cavities it contains is small. The other important conclusion is that the maximum stress concentration factor has a wide scattered distribution.

Distribution of Temperature in GlidArc Discharge Plasma

The GlidArc discharge is one of the main generation methods of non-equilibrium plasma near atmospheric pressures. In general, Gliding Arc discharge is driven by gas flow [1] in axial direction or by magnetic field in circumferential direction. [2] In this paper, a GlidArc discharge driven by rotating-gas-flow in circumferential direction is presented. The principle of the plasma generator is analyzed. The distribution of the temperature in axial direction is measured by a digital thermometer for three different gases. The experimental set-up of the GlidArc plasma is shown in Fig.1. It consists of a center electrode, an outside electrode, a power supply and a gas supply. The shortest distance between the electrodes is 2-3 mm. When a power supply with 10000 volts is attached to the electrodes, the arc will be ignited at the shortest distance. The small plasma column is rotated by the rotating gas flow in circumferential direction and then the rotating arc is driven towards the exit of the setup by the gas flow.

Effects of Thermocapillary Convection on the Growth Process in Industrial 8-inch Silicon Czochralski Growth

Czochralski (CZ) crystal growth process is a widely used technique in manufacturing of silicon crystals and other semiconductor materials. The ultimate goal of the IC industry is to have the highest quality substrates, which are free of point defect, impurities and micro defect clusters. The scale up of silicon wafer size from 200 mm to 300 mm requires large crucible size and more heat power. Transport phenomena in crystal growth processes are quite complex due to melt and gas flows that may be oscillatory and/or turbulent, coupled convection and radiation, impurities and dopant distributions, unsteady kinetics of the growth process, melt crystal interface dynamics, free surface and meniscus, stoichiometry in the case of compound materials. A global model has been developed to simulate the temperature distribution and melt flow in an 8-inch system. The present program features the fluid convection, magnetohydrodynamics, and radiation models. A multi-zone method is used to divide the Cz system into different zones, e.g., the melt, the crystal and the hot zone. For calculation of temperature distribution, the whole system inside the stainless chamber is considered. For the convective flow, only the melt is considered. The widely used zonal method divides the surface of the radiation enclosure into a number of zones, which has a uniform distribution of temperature, radiative properties and composition. The integro-differential equations for the radiative heat transfer are solved using the matrix inversion technique. The zonal method for radiative heat transfer is used in the growth chamber, which is confined by crystal surface, melt surface, heat shield, and pull chamber. Free surface and crystal/melt interface are tracked using adaptive grid generation. The competition between the thermocapillary convection induced by non-uniform temperature distributions on the free surface and the forced convection by the rotation of the crystal determines the interface shape, dopant distribution, and striation pattern. The temperature gradients on the free surface are influenced by the effects of the thermocapillary force on the free surface and the rotation of the crystal and the crucible.

Formation mechanisms of uniform arrays of periodic nanoparticles and nanoripples on 6H-SiC crystal surface induced by femtosecond laser ablation

Uniform arrays of periodic nanoparticles with 80-nm period are formed on 6H-SiC crystal irradiated by circularly polarized 400-nm femtosecond laser pulses. In order to understand the formation mechanism, the morphology evolvement as a function of laser pulse energy and number is studied. Periodic nanoripples are also formed on the sample surface irradiated by linearly polarized 400-, 510- and 800-nm femtosecond laser pulses. All these results support well the mechanism that second-harmonic generation plays an important role in the formation of periodic nanostructures.

Distribution of the microcrystallites generated in borate glass irradiated by femtosecond laser pulses

This paper deals with the distribution of generated microcrystallites in borate glass irradiated by 120 fs laser pulses at a central wavelength of 800 nm. Raman spectroscopy is used to investigate the distribution of the high and low temperature phases of barium metaborate crystals generated in the borate glass. In combination with a microexplosion model, bond-breaking induced by laser irradiation is served as the origin of the formation of BBO crystals. Depending on the laser fluence and cooling conditions, the distribution mechanisms have been discussed. (c) 2006 Elsevier B.V. All rights reserved.

Far-field intensity distribution of a beam generated by a resonator with a phase-unifying mirror

The far-field intensity distribution (FFID) of a beam generated by a phase-unifying mirror resonator was investigated based on scalar diffraction theory. Attention was paid to the parameters, such as obscuration ratio and reflectivity of the phase-unifying mirror, that determine the FFID. All analyses were limited to the TEM00 fundamental mode. (c) 2005 Optical Society of America.

Far-field intensity distribution of beam generated by Gaussian mirror resonator

Based on scalar diffraction theory, we investigated far-field intensity distribution (FFID) of beam generated by Gaussian mirror resonator. We found usable analytical expressions of diffracted field with respect to variation of diffraction parameters. Particular attention was paid to the parameters such as mirror spot size and radius of the Gaussian mirror, which determine the FFID. All analyses were limited to TEM00 fundamental mode. (c) 2004 Elsevier B.V. All rights reserved.