龙游石窟洞室群稳定性分析及最佳旅游路线选择

The large ancient underground rock caverns in Longyou is an important component of grotto cultural. Current task facing the long-term preservation of these unmovable cultural relics is arduous and challenging. The deformation failure of the caverns' surrounding rock is deteriorating. The weathering velocity of these caverns is accelerating. With the strength of caverns' surrounding rock worsening, critical rocks were generated in local regions of the caverns' vault and posing a threat to the security of people passing by. Selection of a maximum-security route and construction a aisle in the caverns might be an efficient way to ensure the security of tourists and reach the target of long-term preservation. The deformation and destruction of the ancient underground caverns is primarily dominated by geological conditions and the special structure of caverns. Based on field investigation, several fundamental conditions for deformation and failure are recognized, and nine deformation and fracture patterns of the Longyou grotto are proposed. In order to judge the stability of caverns’ surrounding rock, the element safety coefficient method is presented. An explicit explanation for the meaning of the method is deduced using Mohr-Coulomb strength criterion. Numerical analyses are carried out in the dissertation through FLAC3D code. Through numerical analysis, the stress distribution regularities of the caverns’ roofs, piles and public side wall are analysed, and the stability properties of caverns’ surrounding rock are also assessed. At the same time, the element safety coefficient method is introduced to contrast the stability degree of different regions in caverns. The above analyses are bases for choosing the optimal tourism routes in the caverns of Longyou grotto. The impact of surface load on the stability of shallow buried cavities in Longyou grotto is evaluated, the results show that building load has significant influence on the stability of the No.1 cavern’s roof, pile and public side wall between the No.1 cavern and the No.2 cavern, pedestrian load has less impact on the stability of surrounding rock than building load. The principles for choosing the optimal tourism routes in the caverns are discussed. With these principles, the dissertation makes a systematic research on the geological analytic method, numerical analytic method and meeting tourism requirements method, which are used in selecting the optimal tourism routes in the caverns. In order to achieve the best effect in the process of tourism routes selection, the above three method are integrated through Theory of Engineering Geomechanics Meta-system(EGMS). According to field investigations, numerical analyses, tourism requirements and expert experiences, the optimal tourism routes through No.1 to No.5 cavern are determined preliminarily. The obtained results from the research work are useful for the security aisle's construction, they also have reference value to other projects in practice.

Thermocapillary Convection in Floating Zones

This paper provides an overview of ongoing studies in the area of thermocapillary convection driven by a surface tension gradient parallel to the free surface in a floating zone. Here, research interests are focused around the onset of oscillatory thermocapillary convection, also known as the transition from quasisteady convection to oscillatory convection. The onset of oscillation depends on a set of critical parameters, and the margin relationship can be represented by a complex function of the critical parameters. The experimental results indicate that the velocity deviation of an oscillatory flow has the same order of magnitude as that of an average flow, and the deviations of other quantities, such as temperature and free surface radii fluctuations, are much smaller when compared with their normal counterparts. Therefore, the onset of oscillation should be a result of the dynamic process in a fluid, and the problem is a strongly nonlinear one. In the past few decades, several theoretical models have been introduced to tackle the problem using analytical methods, linear instability analysis methods, energy instability methods, and unsteady 3D numerical methods. The last of the above mentioned methods is known to be the most suitable for a thorough analysis of strong nonlinear processes, which generally leads to a better comparison with the experimental results. The transition from oscillatory thermocapillary convection to turbulence falls under the studies of chaotic behavior in a new system, which opens a fascinating new frontier in nonlinear science, a hot research area drawing many recent works. This paper reviews theoretical models and analysis, and also experimental research, on thermocapillary connection in floating zones. It cites 93 references.

On dispersion-controlled principles for non-oscillatory shock-capturing schemes

The role of dispersions in the numerical solutions of hydrodynamic equation systems has been realized for long time. It is only during the last two decades that extensive studies on the dispersion-controlled dissipative (DCD) schemes were reported. The studies have demonstrated that this kind of the schemes is distinct from conventional dissipation-based schemes in which the dispersion term of the modified equation is not considered in scheme construction to avoid nonphysical oscillation occurring in shock wave simulations. The principle of the dispersion controlled aims at removing nonphysical oscillations by making use of dispersion characteristics instead of adding artificial viscosity to dissipate the oscillation as the conventional schemes do. Research progresses on the dispersion controlled principles are reviewed in this paper, including the exploration of the role of dispersions in numerical simulations, the development of the dispersion-controlled principles, efforts devoted to high-order dispersion-controlled dissipative schemes, the extension to both the finite volume and the finite element methods, scheme verification and solution validation, and comments on several aspects of the schemes from author's viewpoint.

Non-linear theory of a positive-column in a magnetic-field

A nonlinear theory of an intermediate pressure discharge column in a magnetic field is presented. Motion of the neutral gas is considered. The continuity and momentum transfer equations for charged particles and neutral particles are solved by numerical methods. The main result obtained is that the rotating velocities of ionic gas and neutral gas are approximately equal. Bohm's criterion and potential inversion in the presence of neutral gas motion are also discussed.

3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.

Occurrence of spanning of a submarine pipeline with initial embedment

For better understanding the mechanism of the occurrence of pipeline span for a pipeline with initial embedment, physical and numerical methods are adopted in this study. Experimental observations show that there often exist three characteristic phases in the process of the partially embedded pipeline being suspended: (a) local scour around pipe; (b) onset of soil erosion beneath pipe; and (c) complete suspension of pipe. The effects of local scour on the onset of soil erosion beneath the pipe are much less than those of soil seepage failure induced by the pressure drop. Based on the above observations and analyses, the mechanism of the occurrence of pipeline spanning is analyzed numerically in view of soil seepage failure. In the numerical analyses, the current-induced pressure along the soil surface in the vicinity of the pipe (i.e. the pressure drop) is firstly obtained by solving the N-S equations, thereafter the seepage flow in the soil is calculated with the obtained pressure drop as the boundary conditions along the soil surface. Numerical results indicate that the seepage failure (or piping) may occur at the exit of the seepage path when the pressure gradient gets larger than the critical value. The numerical treatment provides a practical tool for evaluating the potentials for the occurrence of pipe span due to the soil seepage failure.

A comparative study of Young's modulus of single-walled carbon nanotube by CPMD, MD, and first principle simulations

Carbon nanotubes (CNTs), due to their exceptional magnetic, electrical and mechanical properties, are promising candidates for several technical applications ranging from nanoelectronic devices to composites. Young's modulus holds the special status in material properties and micro/nano-electromechanical systems (MEMS/NEMS) design. The excellently regular structures of CNTs facilitate accurate simulation of CNTs' behavior by applying a variety of theoretical methods. Here, three representative numerical methods, i.e., Car-Parrinello molecular dynamics (CPMD), density functional theory (DFT) and molecular dynamics (MD), were applied to calculate Young's modulus of single-walled carbon nanotube (SWCNT) with chirality (3,3). The comparative studies showed that the most accurate result is offered by time consuming DFT simulation. MID simulation produced a less accurate result due to neglecting electronic motions. Compared to the two preceding methods the best performance, with a balance between efficiency and precision, was deduced by CPMD.

填土边坡水文作用机理与破坏过程研究—以深圳羊宝地填土滑坡为例

A full understanding of failure mechanism, critical hydrological condition, and process of mobilization and deposition of a landslide is essential for optimal design of stabilization measure and forecasting of landslide hazard. This requires a quantitative study of hydrological response of a slope to rainfall through field monitoring, laboratory test and numerical modelling. At 13:40 on September 18, 2002, a fill slope failed following a period of prolonged rain in Shenzhen, resulting in 5 fatalities and 31 injuries. The failed mass with a volume about 2.5×104m3 traveled about 140m on level ground. Field monitoring, laboratory test, theoretical analysis and numerical modelling were carried out to undestand the hydrological response and failure mechanism of this fill slope. This thesis mainly focuses on the following aspects: (1) The hydrological responses and failure processes of slopes under rainfall infiltration were reviewed. Firstly, the factors influencing on the hydrological responses of slopes were analysed. Secondly, the change of stress state of slope soil and modelling methods of slope failure under rainfall infiltration were reviewed. (2) The characteristics of the Yangbaodi landslide and associated rainfall triggering the failure were presented. The failure was characterized by shallow flowslide, due to an increase of ground water table caused by rainfall infiltration. (3) A fully automated instrumentation was carried out to monitor rainfall, and saturated – unsaturated hydrological response of the fill slope, using a raingauge, piezometers, tensiometers and moisture probes. A conceptual hydrogeological model was presented based on field monitoring and borehole data. Analysis of monitoring data showed that the high pore water pressure in fill slope was caused by upward flow of semiconfined groundwater in the moderately decomposed granite. (4) Laboratory and in-situ testing was performed to study the physical and mechanical properties of fills. Isotropically consolidated undrained compression tests and anisotropically consolidated constant shear stress tests were carried out to understand the failure mechanism of the fill slope. It is indicated that loosely compacted soil is of strain-softening behaviour under undrained conditions, accompanied with a rapid increase in excess pore water pressure. In anisotropically consolidated constant shear stress tests, a very small axial strain was required to induce the failure and the excess pore water pressure increased quickly at failure. This indicated that static liquefaction caused by rise in groundwater table due to rainfall infiltration occurred. (5) The hydraulic conductivity of the highly and moderately decomposed granite was estimated using monitering data of pore water pressure. A saturated – unsaturated flow was modeled to study the hydrological response of the fill slope using rainfall records. It was observed that the lagged failure was due to the geological conditions and the discrepancy of hydraulic conductivity of slope soils. The hydraulic conductivity of moderately decomposed granite is relatively higher than the other materials, resulting in a semiconfied groundwater flow in the moderately decomposed granite, and subsequent upward flow into the upper fill layer. When the ground water table in the fill layer was increased to the critical state, the fill slope failed. (6) Numerical exercises were conducted to replay the failure process of the fill slope, based on field monitoring, laboratory and in-situ testing. It was found that the fill slope was mobilized by a rapid transfer of the concentrated shear stress. The movement of failure mass was characterized by viscosity fluid with a gradual increase in velocity. The failure process, including mobilization and subsequent movement and deposition, was studied using numerical methods.