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.     

武汉区域气象中心数值预报业务并行计算机系统

对武汉区域气象中心并行计算机系统进行了详细地介绍，分析了并行计算机体系结构、网络和存储系统特点；给出了在并行计算机SP上实现数值预报业务并行化的部分结果；对数值预报模式在串、并行编程环境下的结果进行了分析比较。     

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.     

Geo-effectiveness of solar wind extremes

Examples of extreme events of solar wind and their effect on geomagnetic conditions are discussed here. It is found that there are two regimes of high speed solar wind streams with a threshold of ∼ 850 km s^-1. Geomagnetic activity enhancement rate (GAER) is defined as an average increase in Ap value per unit average increase in the peak solar wind velocity (Vp) during the stream. GAER was found to be different in the two regimes of high speed streams with +ve and-ve IMF. GAER is 0.73 and 0.53 for solar wind streams with +ve and -ve IMF respectively for the extremely high speed streams (< 850 km s^-1). This indicates that streams above the threshold speed with +ve IMF are 1.4 times more effective in enhancing geomagnetic activity than those with -ve IMF. However, the high speed streams below the threshold with -ve IMF are 1.1 times more effective in enhancing geomagnetic activity than those with +ve IMF. The violent solar activity period (October–November 2003) of cycle 23 presents a very special case during which many severe and strong effects were seen in the environment of the Earth and other planets; however, the z-component of IMF (Bz) is mostly positive during this period. The most severe geomagnetic storm of this cycle occurred when Bz was positive.     

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.     

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

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