Magnetic Reconnection, Turbulence, and Collisionless Shock

A short summary of recent progress in measuring and understanding turbulence during magnetic reconnection in laboratory plasmas is given. Magnetic reconnection is considered as a primary process to dissipate magnetic energy in laboratory and astrophysical plasmas. A central question concerns why the observed reconnection rates are much faster than predictions made by classical theories, such as the Sweet–Parker model based on MHD with classical Spitzer resistivity. Often, the local resistivity is conjectured to be enhanced by turbulence to accelerate reconnection rates either in the context of the Sweet–Parker model or by facilitating setup of the Pestchek model. Measurements at a dedicated laboratory experiment, called MRX or Magnetic Reconnection Experiment, have indicated existence of strong electromagnetic turbulence in current sheets undergoing fast reconnection. The origin of the turbulence has been identified as right-hand polarized whistler waves, propagating obliquely to the reconnecting field, with a phase velocity comparable to the relative drift velocity. These waves are consistent with an obliquely propagating electromagnetic lower-hybrid drift instability driven by drift speeds large compared to the Alfven speed in high-beta plasmas. Interestingly, this instability may explain electromagnetic turbulence also observed in collisionless shocks, which are common in energetic astrophysical phenomena.     

THE SIMULATED STORAGE CAPACITY OF FLOOD AND WATERLOGGING IN THE TYPICAL AGRICULTURAL REGION IN JIANGHAN PLAIN

１ＩＮＴＲＯＤＵＣＴＩＯＮＪｉａｎｇｈａｎＰｌａｉｎｉｓｏｎｅｏｆｔｈｅｉｍｐｏｒｔａｎｔｂａｓｅｓｏｆＣｈｉｎａ′ｓａｇｒｉｃｕｌｔｕｒａｌｐｒｏｄｕｃｔｉｖｉｔｙｌｏｃａｔｅｄｉｎｔｈｅｍｉｄｄｌｅＣｈａｎｇｊｉａｎｇＲｉｖｅｒｖａｌｌｅｙｏｆｃｅｎｔｒａｌＣｈｉｎａ．ＴｈｅＪｉａｎｇｈａｎｐｌａｉｎａｐｐｒｏｘｉｍａｔｅｌｙ６．６×１０４ｋｍ２ｉｎａｒｅａａｎｄｏｖｅｒ３０ｍｉｌｌｉｏｎｉｎｐｏｐｕｌａｔｉｏｎｉｓｔｈｅｒｅｎｏｗｎｅｄ“ｃｏｕｎｔｒｙｏｆｆｉｓｈａｎｄｒｉｃｅ＂ｉｎ…     

重力插值方法研究

为了根据离散观测数据构制连续空间重力变化图像，分析和讨论了３种数值插值方法，计算结果表明多面函数方法插值精度最高。由于逐步回归分析筛选核函数中心点的计算繁琐，文中提出根据分形理论和Ｓｈａｎｎｏｎ取样定理来确定核函数中心点。对滇西试验场进行模拟试算，插值精度可达到４～５（１０－８ｍｓ－２）。     

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.     

Optimization of grid-drilling using computer simulation

A computer simulation method has been developed to find efficient drilling grids for mineral deposits. A well-known ore deposit is used as a model to develop an efficient pattern for undiscovered ore bodies in the same area or in other prospects where similar geometry is suspected. The model for this study is the Austinville, Virginia deposit, a Mississippi Valley-type deposit composed of 17 ore bodies totaling 34 million short tons (30 million metric tons). The method employs a computer program that simulates drilling the model deposit with different patterns, including various levels of follow-up drilling. Follow-up holes are drilled in fences at one half the original spacing around holes in the grid that show ore-grade mineralization. Each pattern is drilled 100 times from random starting locations to provide a range of outcomes of drilling, including the best, worst, and most likely. For this study, patterns of 100 drill holes were composed of 10 fences spaced 1000–5000 feet (305–1524 m) apart, each with 10 holes spaced 200–1000 feet (61–305 m) apart. In all, 25 grids were used with zero to three levels of follow-up drilling. The 600/2000 grid, with drill holes spaced 600 feet (183 m) apart in fences spaced 2000 feet (610 m) apart, was compared with the 200/5000 grid because they represented contrasting outcomes. The 600/2000 grid penetrated many ore bodies consistently but with few multiple hits to individual ore bodies; whereas the 200/5000 grid inconsistently penetrated few ore bodies with many multiple hits. The 600/2000 grid was more efficient than the 200/5000 grid at hitting large ore bodies of 1,000,000 short tons or greater (900,000 metric tons or greater) and was made more effective by adding one cycle of follow-up drilling. The 600/2000 grid had a 97% chance of hitting one or more large ore bodies with at least one drill hole per ore body, and the 200/5000 grid had a 64% chance. Once hit, there was an 82% chance that the largest ore body would be penetrated by three or more holes when using the 600/2000 grid and an 88% chance using the 200/5000 grid.     

A Critical review of theoretical models of negatively polarized light scattered by atmosphereless solar system bodies

About a dozen physical mechanisms and models aspire to explain the negative polarization of light scattered by atmosphereless celestial bodies. This is too large a number for the reliable interpretation of observational data. Through a comparative analysis of the models, our main goal is to answer the question: Does any one model have an advantage over the others? Our analysis is based on new laboratory polarimetric and photometric data as well as on theoretical results. We show that the widely used models due to Hopfield and Wolff cannot realistically explain the phase-angle dependence of the degree of polarization observed at small phase angles. The so-called interference or coherent backscattering mechanism is the most promising model. Models based on that mechanism use well-defined physical parameters to explain both negative polarization and the opposition effect. They are supported by laboratory experiments, particularly those showing enhancement of negative polarization with decreasing particle size down to the wavelength of light. According to the interference mechanism, pronounced negative branches of polarization, like those of C-class asteroids, may indicate a high degree of optical inhomogeneity of light-scattering surfaces at small scales. The mechanism also seems appropriate for treating the negative polarization and opposition effects of cometary dust comae, planetary rings, and the zodiacal light.     

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

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