A Block-Based Theoretical Model Suited to Gravitational Sliding

An original theoretical model has been devised to simulate mass flow over hill slopes due to gravitational sliding. The sliding mass is discretized into a sequence of contiguous blocks which are subjected to gravitational forces, to bottom friction and to surface resistance stresses that are generally negligible for subaerial flows, but are relevant for submarine slides. The blocks interact with each other while sliding down the hill flanks because of internal forces that dissipate mechanical energy and produce a momentum exchange between the individual blocks, yet conserving the total momentum of the mass. Internal forces are expressed in terms of interaction coefficients depending on the instantaneous distance between the block centers of mass, which is a measure of the deformation experienced by the blocks: the functional dependence includes three parameters, namely the interaction intensity ¯, the deformability parameter and the shape parameter , by means of which a wide range of interaction types can be fully accounted for. The time integration is performed numerically by solving the equations for the block velocities and positions at any time ti by means of the block accelerations at the previous time t_i-1, and by subsequently updating the block accelerations, which allows to proceed iteratively to the following times. The model has been tested against laboratory results available from literature and by means of several numerical experiments involving a simplified geometry both for the sliding body and the basal surface, with the purpose of clarifying the influence of the model parameters on the slide dynamics. The model improves the performance of the existing kinematic models for slides, moreover preserving an equivalent numerical simplicity. Future applications and possible improvements of this model are suggested.     

含煤层地质环境下地震波场的数值模拟

本文对含低速煤层地质环境下弹性波场多波多分量地震资料进行了二维数值模拟研究,对人工边界反射进行了有效处理,频散效应得到了有效的压制,对几种不同激发与观测排列方式下的弹性波资料进行了模型计算与分析。     

基于GIS城市洪水淹没模拟分析

城市洪水淹没范围动态模拟分析是城市防洪规划与防治决策的一个重要基础工作。本文主要研究洪水从翻堤口出发在地理空间蔓延、扩散、动态行进及确定淹没范围的数字模拟方法。为此,根据数学形态学及测地圆概念,研究设计了洪水扩散范围的“膨胀”模拟算法和淹没范围搜索算法,用于查询淹没通块中从翻堤处到任一点之间的淹没路径和ｔ时刻洪水扩散范围。     

Near-infrared imaging in H2 of molecular (CO) outflows from young stars

We have imaged several known molecular (CO) outflows in H₂ v=1-0 S(1) and wide-band K in order to identify the molecular shocks associated with the acceleration of ambient gas by outflows from young stars. We detected H₂ line emission in all the flows we observed: L 1157, VLA 1623, NGC 6334I, NGC 2264G, L 1641N and Haro 4-255. A comparison of the H₂ data with CO outflow maps strongly suggests that prompt entrainment near the head of a collimated jet probably is the dominant mechanism for producing the CO outflows in these sources.     

Dynamical evolution of comets and the problem of cometary fading

Possibilities to explain the observed 1/a-distribution are discussed in the light of improved understanding of the dynamical evolution of long-period comets. It appears that the fading problem applies both to single-injection and continuous-injection models. Although uncertainties due to nongravitational effects do not allow detailed results to be drawn from the observed 1/a-distribution at small perihelion distance q, that for q 1.5 AU shows that a constant fading probability cannot explain all the features of the observed distribution. Assuming that comets can reappear following a period of fading, values for the assumed constant fading and renewal probabilities, and the total cometary flux have been estimated for q > 1.5 AU.     

Dynamics of Mars-orbiting dust: Effects of light pressure and planetary oblateness

This paper deals with dynamics of impact ejecta from Phobos and Deimos initially on near-circular equatorial orbits around Mars. For particles emitted in a wide size regime of 1 micron and greater, and taking into account the typical particle lifetimes to be less than 100 years, the motion is governed by two perturbing forces: solar radiation pressure and influence of Mars' oblateness. The equations of motion of particles in Lagrangian non-singular elements are deduced and solved, both analytically and numerically, for different-sized ejecta. We state that the coupled effect of both forces above is essential so that on no account can the oblateness of Mars are be neglected. The dynamics of grains prove to be quite different for the ejecta of Phobos and Deimos. For Deimos, the qualitative results are relatively simple and imply oscillations of eccentricity and long-term variations of orbital inclination, with amplitudes and periods both depending on grain size. For Phobos, the dynamics are shown to be much more complicated, and we discuss it in detail. We have found an intriguous peculiar behavior of debris near 300 µm in size. Another finding is that almost all the Phobos ejecta with radii less than 30 µm (against the values of 5 to 20 µm adopted earlier by many authors) should be rapidly lost by collisions with martian surface. The results of the paper may be the base for constructing an improved model of dust belts that presumably exist around Mars.     

The densest prestellar condensations: H2CO Studies ofρ Oph B1

We present new images of the well-known molecular outflow and Herbig-Haro complex L 1551-IRS 5. Deep, high-resolution images of the central region of the flow in [SII] 6716,6731 and H (6565 Å) are complemented by a mosaic of much of the CO outflow in H₂ v=1-0 S(1). While the optical data trace the intermediate-to-high excitation shocks in the flow (v _shock > 30 – 50 km s^–1), the near-IR data reveal the lower-excitation, molecular shocks (v _shock 10–50 km s^–1). In particular, the H₂ data highlight the regions where the flow impacts and shocks ambient molecular gas.     

Phase transitions in the ISM a source of dissipative behaviour

By considering a simple fluid model, we investigate the role of phase transitions in the ISM on the galaxy- scale gas dynamics. Cooling and heating timescales in the ISM are typically shorter than typical galactic rotation timescales, so the individual phases in the ISM can be assumed to be in temperature equilibrium with the radiation field. Using this approximation we can construct an equation of state which depends upon the average density and mass fractions in the individual phases. Previous studies suggest that there is an equilibrium phase fraction as a function of pressure. We incorporate evolution towards this equilibrium state as a relaxation term with a time to obtain equilibrium . We derive a condition in terms of a critical Mach number when one dimensional shocks should be continuous. For small values of the relaxation time we show that the relaxation term acts like a viscosity. We show with one dimensional simulations that increasing causes shocks to become smoother. This work suggests that phase changes can strongly effect the gas dynamics of the ISM across spiral arms and bars.     

Large eddy simulation of hydrocyclone—prediction of air-core diameter and shape

Hydrocyclones are widely used in the mining and chemical industries. An attempt has been made in this study, to develop a CFD (computational fluid dynamics) model, which is capable of predicting the flow patterns inside the hydrocyclone, including accurate prediction of flow split as well as the size of the air-core. The flow velocities and air-core diameters are predicted by DRSM (differential Reynolds stress model) and LES (large eddy simulations) models were compared to experimental results. The predicted water splits and air-core diameter with LES and RSM turbulence models along with VOF (volume of fluid) model for the air phase, through the outlets for various inlet pressures were also analyzed. The LES turbulence model led to an improved turbulence field prediction and thereby to more accurate prediction of pressure and velocity fields. This improvement was distinctive for the axial profile of pressure, indicating that air-core development is principally a transport effect rather than a pressure effect.     

The effects of carbon dioxide leakage on fractures, and water quality of potable aquifers during geological sequestration of CO2

Since industrial revolution, the ＂greenhouse effect＂ is one of the most important global environmental issues. Of all the greenhouse gases, CO2 is responsible for about 64% of the enhanced ＂greenhouse effect＂, making it the target for mitigation, so reducing anthropogenic discharge of carbon dioxide attracts more and more attention. Geological sequestration of CO2 in deep saline aquifers is one of the most promising options. But because unknown fractures and faults may exist in the caprock layers which can prevent the leakage of CO2, CO2 will leak upward into upper potable aquifers, and lead to adverse impacts on the shallow potable aquifers. In order to assess the potential effect of CO2 leakage from underground storage reservoirs on fractures and water quality of potable aquifers, this study used the non-isothermal reactive geochemical transport code TOUGHREACT developed by Xu et al to establish a simplified 2-D model of CO2 underground sequestration system, which includes deep saline aquifers, caprock layers, and shallow potable aquifers, and study and analyze the changes of mineral and aqueous components. The simulation results indicated that the minerals of deep saline aquifers and fractures should be mainly composed of aluminosilicate and silicate minerals, which not only enhance the mass of CO2 sequestrated by mineral trapping, but also decrease the porosity and permeability of caprock layers and fractures to prevent and reduce CO2 leakage. The results from deep saline aquifers showed that the mass of carbon dioxide trapped by minerals and solution phases is limited, the rest remained as a supercritical phase, and so once the caprock aquifers have some unknown fractures, the free carbon dioxide phase may leak from CO2 geologic sequestration reservoirs by buoyancy.