Simulation of Vibration Caused by Blasting for a Layered Slope

Recent studies showed that vibration caused by blasting mainly reflects the property of geological structure itself neighboring the blasting center. Different vibration signals can be collected for different geological structures under blasting. Hence, vibration signal can be used to identify geological structure, especial for a slope with a weak layer. As the geological structure for a practical slope is usually complicated, the simulation of vibration caused by blasting should be carried out first. Generally, the material in a certain zone near the blasting center will undergo damage, so the physical model to simulate this region is the most concerned. In this paper, the damaged zone near blasting center is neglected, and the blasting load can be considered being applied on the interface between the damaged zone and undamaged zone. Regarding the relations between the weight of explosive, the size of damaged zone, and the dynamic loading to propagate out away, the vibration caused by blasting for a practical layered slope is simulated. Compared with the measured signal in site, it can be seen that the simulating result is in well agreement with that of practical testing. The results also indicate that the farther the testing point apart from the blasting center, the more accurate the simulation is.

Simulation of Natural Gas Production in Hydrate Reservoirs

This paper simulates a one-dimensional physical model of natural gas production from hydrate dissociation in a reservoir by depressurization. According to the principles of solid hydrate decomposition in stratum and flow of natural gas in porous medium, the pressure governing equations for both gas zone and hydrate zone are set up based on the physical production model. Using the approximation reported by N. N. Verigin et al. (1980), the nonlinear governing equations are simplified and the self-similar solutions are obtained. Through calculation, for different reservoir parameters, the distribution characters of pressure are analyzed. The decline trend of natural gas production rate with time is also studied. The simulation results show that production of natural gas from a hydrate reservoir is very sensitive to several reservoir parameters, such as wellbore pressure and stratum porosity and permeability.

Simulation of Vortex-dominant and Turbulent Flows in Natural Environment

In the present talk, the simulation of vortex dominant and turbulent flows are primarily addressed. To cope with complicated circumstances in environmental flows we illustrate the strategy of combining simplified physical model and suitable algorithm by a few examples.

Numerical Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

Computation of turbulence-generated noise by large-eddy simulation

The hybrid method of large eddy simulation (LES) and the Lighthill analogy is being developed to compute the sound radiated from turbulent flows. The results obtained from the hybrid method are often contaminated by the absence of small scales in LES, since the energy level of sound is much smaller than that of turbulent flows. Previous researches investigate the effects of subgrid sacle (SGS) eddies on the frequency spectra of sound radiated by isotropic turbulence and suggest a SGS noise model to represent the SGS contributions to the frequency spectra. Their investigations are conducted in physical space and are unavoidably influenced by boundary conditions. In this paper, we propose to perform such calculations in Fourier space so that the effects of boundary conditions can be correctly treated. Posteriori tests are carried out to investigate the SGS contribution to the sound. The results obtained recover the -7/2 law within certain wave-number ranges, but under-estimate the amplitudes of the frequency spectra. The reason for the underestimation is also discussed.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been devel- oped to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Monte Carlo Simulation Of Thermal Fluctuations Below The Onset Of Rayleigh-Benard Convection

The density fluctuations below the onset of convection in the Rayleigh-Benard problem are studied with the direct simulation Monte Carlo method. The particle simulation results clearly show the connection between the static correlation functions of fluctuations below the critical Rayleigh number and the flow patterns above the onset of convection for small Knudsen number flows (Kn=0.01 and Kn=0.005). Furthermore, the physical nature for no convection in the Rayleigh-Benard problem under large Knudsen number conditions (Kn>0.028) is explained based on the dynamics of fluctuations.

An Improved CE/SE Scheme for Numerical Simulation of Gaseous and Two-Phase Detonations

A new structure of solution elements and conservation elements based on rectangular mesh was pro- posed and an improved space-time conservation element and solution element (CE/SE) scheme with sec- ond-order accuracy was constructed. Furthermore, the application of improved CE/SE scheme was extended to detonation simulation. Three models were used for chemical reaction in gaseous detonation. And a two-fluid model was used for two-phase (gas–droplet) detonation. Shock reflections were simu- lated by the improved CE/SE scheme and the numerical results were compared with those obtained by other different numerical schemes. Gaseous and gas–droplet planar detonations were simulated and the numerical results were carefully compared with the experimental data and theoretical results based on C–J theory. Mach reflection of a cellular detonation was also simulated, and the numerical cellular pat- terns were compared with experimental ones. Comparisons show that the improved CE/SE scheme is clear in physical concept, easy to be implemented and high accurate for above-mentioned problems.

Broad-bandwidth semi-noncollinear optical parametric amplification in periodically poled LiNbO3 based on tilted quasi-phase-matched gratings

On the basis of noncollinear optical parametric amplification in periodically poled lithium niobate (PPLN) which is realized by quasi-phase matching (QPM) technology, we consider the possibility of semi-noncollinear phase matching between collinear and noncollinear geometries by tilting a PPLN-crystal's parallel grating at a sure angle. Numerical simulation with proper parameters shows that we can achieve a broader optical parametric amplification (OPA) bandwidth than that of noncollinear geometry. About 121 nm at a signal wavelength of 800 and 70 nm at a signal wavelength of 1064 nm under optimal conditions are obtained when the crystal length is 9 mm.

Acceleration of ultra-thin electron layer.Analytical treatment compared with 1D-PIC simulation

In this paper, we apply an analytical model [V.V. Kulagin et al., Phys. Plasmas 14, 113101 (2007)] to describe the acceleration of an ultra-thin electron layer by a schematic single-cycle laser pulse and compare with one-dimensional particle-in-cell (1D-PIC) simulations. This is in the context of creating a relativistic mirror for coherent backscattering and supplements two related papers in this EPJD volume. The model is shown to reproduce the 1D-PIC results almost quantitatively for the short time of a few laser periods sufficient for the backscattering of ultra-short probe pulses.

Simulation of growth process of thin film on non-planar substrate

In (2 + 1) dimension, growth process of thin film on non-planar substrate in Kuramoto-Sivashinsky model is studied with numerical simulation approach. 15 x 15 semi-ellipsoids arranged orderly on the surface of substrate are used to represent initial rough surface. The results show that at the initial stage of growth process, the surface morphology of thin film appears to be grid-structure, and the interface width constantly decreases with the growth time, then reaches minimum. However, the grid-structure becomes ambiguous, and granules of different sizes distribute evenly on the surface of thin film with the increase of growth time. Thereafter, the average size of granules and the interface width gradually increase, and the surface morphology of thin film presents fractal properties. The numerical results of height-height correlation functions of thin film verify the surface morphology of thin film to be fractal for a longer growth time. By fitting of the height-height correlation functions of thin film with different growth times, the growth process is described quantitatively. (c) 2004 Elsevier B.V. All rights reserved.

Characterize dynamic conformational space of human CCR5 extracellular domain by molecular modeling and molecular dynamics simulation

The chemokine receptor CCR5 is the receptor for several chemokines and major coreceptor for R5 human immunodeficiency virus type-1 strains entry into cell. Three-dimensional models of CCR5 were built by using homology modeling approach and 1 ns molecular dynamics (MD) simulation, because studies of site-directed mutagenesis and chimeric receptors have indicated that the N-terminus (Nt) and extracellular loops (ECLs) of CCR5 are important for ligands binding and viral fusion and entry, special attention was focused on disulfide bond function, conformational flexibility, hydrogen bonding, electrostatic interactions, and solvent-accessible surface area of Nt and ECLs of this protein part. We found that the extracellular segments of CCR5 formed a well-packet globular domain with complex interactions occurred between them in a majority of time of MID simulation, but Nt region could protrude from this domain sometimes. The disulfide bond Cys20-Cys269 is essential in controlling specific orientation of Nt region and maintaining conformational integrity of extracellular domain. RMS comparison analysis between conformers revealed the ECL1 of CCR5 stays relative rigid, whereas the ECL2 and Nt are rather flexible. Solvent-accessible surface area calculations indicated that the charged residues within Nt and ECL2 are often exposed to solvent. Integrating these results with available experimental data, a two-step gp120-CCR5 binding mechanism was proposed. The dynamic interaction of CCR5 extracellular domain with gp120 was emphasized. (C) 2004 Elsevier B.V. All rights reserved.

Exploring consensus mRNA secondary (folding) structure units by stochastic sampling and folding simulation

It is extremely difficult to explore mRNA folding structure by biological experiments. In this report, we use stochastic sampling and folding simulation to test the existence of the stable secondary structural units of-mRNA, look for the folding units, and explore the probabilistic stabilization of the units. Using this method, We made simulations for all possible local optimum secondary structures of a single strand mRNA within a certain range, and searched for the common parts of the secondary structures. The consensus secondary structure units (CSSUs) extracted from the above method are mainly hairpins, with a few single strands. These CSSUs suggest that the mRNA folding units could be relatively stable and could perform specific biological function. The significance of these observations for the mRNA folding problem in general is also discussed. (c) 2004 Elsevier B.V. All rights reserved.

Mode Simulation for Midinfrared Microsquare Resonators With Sloped Sidewalls and Confined Metals

EQUILATERAL-TRIANGLE; MU-M; LASERS; MICROLASERS; MICRODISK Abstract: Mode characteristics for midinfrared microsquare resonators with sloped sidewalls and confined metal layers are investigated by finite-difference time-domain (FDTD) techniques. For a microsquare with a side length of 10 mu m, the mode quality (Q)-factors of 8329, 4772, and 2053 are obtained for TM5,7 mode at wavelength 7.1 mu m by three-dimensional FDTD simulations, as the tilting angles of the side walls are 90 degrees, 88 degrees, and 86 degrees, respectively. Furthermore, microsquare resonators laterally surrounded by SiO2 and metal layers are investigated by the two-dimensional FDTD technique for the metal layers of Au, Ti-Au, Ag-Au, and Ti-Ag-Au, respectively.

Investigation of gain recovery for InAs/GaAs quantum dot semiconductor optical amplifiers by rate equation simulation

The gain recoveries in quantum dot semiconductor optical amplifiers (QD SOAs) are numerically studied by rate equation simulation. Similar to the optical pump-probe experiment, the injection of double 150 fs optical pulses is used to simulate the gain recovery of a weak continuous signal under different injection levels, inhomogeneous broadenings, detuning wavelengths, and pulse signal energies for the QD SOAs. The obtained gain recoveries are then fitted by a response function with multiple exponential terms to determine the response times. The gain recovery can be described by three exponential terms with the time constants, which can be explained as carrier relaxation from the excited state to the ground state, carrier captured by the excited state from the wetting layer, and the supply of the wetting layer carriers. The fitted lifetimes decrease with the increase of the injection currents under gain unsaturation, slightly decrease with the decrease of inhomogeneous broadening of QDs, and increase with the increase of detuning wavelength between continuous signal and pulse signal and the increase of the pulse energy.