Study on numerical simulation to water lowering in an open coal mine digging

In view of the situation of excavation of open coal mine for the underground water disaster,we should carry out simulation studies for the numerical value of the water lowering project and improve the accuracy and the level of the water lowering project.On the basis of the hydrological geological conditions of certain open mine digging,a more reasonable seepage numerical model was built according to MODFLOW.It was simulated in advance that the process of the confined water level descending with the time,and combining with the actual observations to test the correctness of the model.The calculation showed that the results coincided well with the results of actual measurement.Based on this,different water lowering numerical simulations were built for the open coal mine digging.It could be simulated and forecast that the changes of the groundwater level in drainage process within and outside the mine pit,and it was quantitatively assessed that the possible water lowering result of the opencast water drainage process,which provide an important basis for the actual water lowering project and the possible project disposal.     

Upscaling of elastic properties for large scale geomechanical simulations

Large scale geomechanical simulations are being increasingly used to model the compaction of stress dependent reservoirs, predict the long term integrity of under‐ground radioactive waste disposals, and analyse the viability of hot‐dry rock geothermal sites. These large scale simulations require the definition of homogenous mechanical properties for each geomechanical cell whereas the rock properties are expected to vary at a smaller scale. Therefore, this paper proposes a new methodology that makes possible to define the equivalent mechanical properties of the geomechanical cells using the fine scale information given in the geological model. This methodology is implemented on a synthetic reservoir case and two upscaling procedures providing the effective elastic properties of the Hooke's law are tested. The first upscaling procedure is an analytical method for perfectly stratified rock mass, whereas the second procedure computes lower and upper bounds of the equivalent properties with no assumption on the small scale heterogeneity distribution. Both procedures are applied to one geomechanical cell extracted from the reservoir structure. The results show that the analytical and numerical upscaling procedures provide accurate estimations of the effective parameters. Furthermore, a large scale simulation using the homogenized properties of each geomechanical cell calculated with the analytical method demonstrates that the overall behaviour of the reservoir structure is well reproduced for two different loading cases. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

Cosmic structure growth and dark energy

Dark energy has a dramatic effect on the dynamics of the Universe, causing the recently discovered acceleration of the expansion. The dynamics are also central to the behaviour of the growth of large-scale structure, offering the possibility that observations of structure formation provide a sensitive probe of the cosmology and dark energy characteristics. In particular, dark energy with a time-varying equation of state can have an influence on structure formation stretching back well into the matter-dominated epoch. We analyse this impact, first calculating the linear perturbation results, including those for weak gravitational lensing. These dynamical models possess definite observable differences from constant equation of state models. Then we present a large-scale numerical simulation of structure formation, including the largest volume to date involving a time-varying equation of state. We find the halo mass function is well described by the Jenkins et al. mass function formula. We also show how to interpret modifications of the Friedmann equation in terms of a time-variable equation of state. The results presented here provide steps toward realistic computation of the effect of dark energy in cosmological probes involving large-scale structure, such as cluster counts, the Sunyaev–Zel'dovich effect or weak gravitational lensing.     

Numerical Stabilization of Orbital Motion

Mainly, the author focuses on Baumgarte's method and its applications in satellite, asteroid, stellar and planetary problems. In the paper arguments are given for the use of energy relations for stabilization in the elliptical two-body problem. Stabilizing properties of Baumgarte's equations and others are discussed. A simple approach is proposed for stabilizing the equations of almost circular motion. By using Baumgarte's technique, the author derives stabilized equations of perturbed restricted three-body problem. It is shown experimentally that stabilization in the problems mentioned above can raise the accuracy of numerical integration by several orders.     

Coronal Holes and Icosahedral Symmetry

The study of the expansion of the solar wind out of a system of coronal holes is continued. To this end, we consider the numerical integration of partial differential equations for problems with icosahedral symmetry, in general. First, employing Weyl theory, orbifold coordinates are introduced. Second, the icosahedral coordinates are discussed in detail. Third, following an analysis of the properties of these coordinates and the derivation of a few expressions useful for grid construction, various alternatives for the distribution of lattice points required for numerical integration are considered. A comparison of these numerical grids motivates the choice of a specific grid optimized for the numerical integration carried out in the accompanying paper by Kalish et al.(2002). This revised version was published online in July 2006 with corrections to the Cover Date.     

介质横向不均匀对震后地表位移场的影响

利用三维有限元模拟两层地球模型(包括弹性层和粘弹性层)介质的横向不均匀和断层错动的复杂性对震后地表位移场的影响，其数值结果包括弹性解和粘弹性解。结果表明，震后位移场的分布特征不仅与粘弹性层的流变特性有关，还与弹性层介质的剪切模量以及断层的同震错动有关。     

Studying Hydrodynamic Instability Using Shock-Tube Experiments

The hydrodynamic instability, which develops on the contact surface between two fluids, has great importance in astrophysical phenomena such as the inhomogeneous density distribution following a supernova event. In this event acceleration waves pass across a material interface and initiate and enhance unstable conditions in which small perturbations grow dramatically. In the present study, an experimental technique aimed at investigating the above-mentioned hydrodynamic instability is presented. The experimental investigation is based on a shock-tube apparatus by which a shock wave is generated and initiates the instability that develops on the contact surface between two gases. The flexibility of the system enables one to vary the initial shape of the contact surface, the shock-wave Mach number, and the density ratio across the contact surface. Three selected sets of shock-tube experiments are presented in order to demonstrate the system capabilities: (1) large-initial amplitudes with low-Mach-number incident shock waves; (2) small-initial amplitudes with moderate-Mach-number incident shock waves; and (3) shock bubble interaction. In the large-amplitude experiments a reduction of the initial velocity with respect to the linear growth prediction was measured. The results were compared to those predicted by a vorticity-deposition model and to previous experiments with moderate- and high-Mach number incident shock waves that were conducted by others. In this case, a reduction of the initial velocity was noted. However, at late times the growth rate had a 1/t behavior as in the small-amplitude low-Mach number case. In the small-amplitude moderate-Mach number shock experiments a reduction from the impulsive theory was noted at the late stages. The passage of a shock wave through a spherical bubble results in the formation of a vortex ring. Simple dimensional analysis shows that the circulation depends linearly on the speed of sound of the surrounding material and on the initial bubble radius.