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.     

A new treatment of the albedo radiation pressure in the case of a uniform albedo and of a spherical satellite

The purpose of this paper is to present a model for the radiation pressure acceleration of a spherical satellite, due to the radiation reflected by a planet with a uniform albedo. A particular choice of variables allows one to reduce the surface integrals over the lit portion of the planet visible to the satellite to one-dimensional integrals. Exact analytical expressions are found for the integrals corresponding to the case where the spacecraft does not "see" the terminator. The other integrals can be computed either numerically, or analytically in an approximate form. The results are compared with those of Lochry (1966). The model is applied to Magellan, a spacecraft orbiting Venus.     

云南大地热流及地热地质问题

本文报道了云南20多个大地热流测试数据。据热流分布特点并参考前人有关资料,结合地球物理测深资料、岩石热物性与放射性生热率测试数据,对各区地壳深度范围的温度、热流配分组成作了粗略分析,划分出各有特征的若干地质-热流区。云南普遍地热偏高(滇东南小范围除外)的主要原因是新构造运动强烈;怒江以西近代岩浆活动亦强烈,具有典型的板块交汇带的基本地热特点     

Reservoir description using dynamic parameterisation selection with a combined stochastic and gradient search

There is a correspondence between flow in a reservoir and large scale permeability trends. This correspondence can be derived by constraining reservoir models using observed production data. One of the challenges in deriving the permeability distribution of a field using production data involves determination of the scale of resolution of the permeability. The Adaptive Multiscale Estimation (AME) seeks to overcome the problems related to choosing the resolution of the permeability field by a dynamic parameterisation selection. The standard AME uses a gradient algorithm in solving several optimisation problems with increasing permeability resolution. This paper presents a hybrid algorithm which combines a gradient search and a stochastic algorithm to improve the robustness of the dynamic parameterisation selection. At low dimension, we use the stochastic algorithm to generate several optimised models. We use information from all these produced models to find new optimal refinements, and start out new optimisations with several unequally suggested parameterisations. At higher dimensions we change to a gradient-type optimiser, where the initial solution is chosen from the ensemble of models suggested by the stochastic algorithm. The selection is based on a predefined criterion. We demonstrate the robustness of the hybrid algorithm on sample synthetic cases, which most of them were considered insolvable using the standard AME algorithm.     

Rock fracturing by explosive energy: review of state-of-the-art

A study of the dynamic rock fracture initiation and propagation due to explosive energy is presented through a detailed state-of-the-art review. Explosive energy dissipation in crushing and fracturing is examined and the various means to enhance the explosive energy utilization for dynamic rock fracturing are reviewed. The study highlights the need for a better understanding of the dynamic fracturing process particularly in the presence of in situ stresses in the rock mass.     

Modeling the cyclic swelling pressure of mudrock using artificial neural networks

The stochastic nature of the cyclic swelling behavior of mudrock and its dependence on a large number of interdependent parameters was modeled using Time Delay Neural Networks (TDNNs). This method has facilitated predicting cyclic swelling pressure with an acceptable level of accuracy where developing a general mathematical model is almost impossible. A number of total pressure cells between shotcrete and concrete walls of the powerhouse cavern at Masjed–Soleiman Hydroelectric Powerhouse Project, South of Iran, where mudrock outcrops, confirmed a cyclic swelling pressure on the lining since 1999. In several locations, small cracks are generated which has raised doubts about long term stability of the powerhouse structure. This necessitated a study for predicting future swelling pressure. Considering the complexity of the interdependent parameters in this problem, TDNNs proved to be a powerful tool. The results of this modeling are presented in this paper.     

The distribution of radiogenic heat production as a function of depth in the Sierra Nevada Batholith, California

Geochemical analyses and geobarometric determinations have been combined to create a depth vs. radiogenic heat production database for the Sierra Nevada batholith, California. This database shows that mean heat production values first increase, then decrease, with increasing depth. Heat production is 2 μW/m³ within the 3-km-thick volcanic pile at the top of the batholith, below which it increases to an average value of 3.5 μW/m³ at 5.5 km depth, then decreases to 0.5–1 μW/m³ at 15 km depth and remains at these values through the entire crust below 15 km. Below the crust, from depths of 40–125 km, the batholith's root and mantle wedge that coevolved beneath the batholith appears to have an average radiogenic heat production rate of 0.14 μW/m³. This is higher than the rates from most published xenolith studies, but reasonable given the presence of crustal components in the arc root assemblages. The pattern of radiogenic heat production interpreted from the depth vs. heat production database is not consistent with the downward-decreasing exponential distribution predicted from modeling of surface heat flow data. The interpreted distribution predicts a reasonable range of geothermal gradients and shows that essentially all of the present day surface heat flow from the Sierra Nevada could be generated within the 35 km thick crust. This requires a very low heat flux from the mantle, which is consistent with a model of cessation of Sierran magmatism during Laramide flat-slab subduction, followed by conductive cooling of the upper mantle for 70 m.y. The heat production variation with depth is principally due to large variations in uranium and thorium concentration; potassium is less variable in concentration within the Sierran crust, and produces relatively little of the heat in high heat production rocks. Because silica content is relatively constant through the upper 30 km of the Sierran batholith, while U, Th, and K concentrations are highly variable, radiogenic heat production does not vary directly with silica content.     

New Approach for Estimation of Static and Seismic Active Earth Pressure

To estimate static and seismic active earth pressure (P_ad) on a rigid retaining wall, numerical analyses using different step sizes have been carried out in this paper, based on the modified Culmann line method by considering Coulomb’s planar rupture surface. Equivalent pseudo-static seismic forces are considered in the analysis. A new concept of modified unit weight by considering ground surcharge is introduced under static and seismic conditions. By numerical analysis, area of soil (A) has been estimated to obtain the ratio of A/A₀ where A₀ is θh², θ is the angle between retaining structure and ground surface and h is the vertical height of the wall. This ratio remains constant for a particular type of soil and has been used to estimate the maximum active earth pressure using force diagram. Results are provided in tabular form for easy calculation of the coefficient of static and seismic active earth pressure. Present results by considering the new technique, compares well with the results obtained by earlier researchers.     

Effect of a submerged bay-mouth breakwater on tsunami behavior analyzed by 2D/3D hybrid model simulation

The hybrid numerical model had been developed to simulate a complicated 3D flow around structures generated by tsunami. In the model, the conventional 2D model is adopted for the wide region far from structures and the 3D non-hydrostatic pressure model is used in the limited region adjacent to structures. The applicability of the model is shown by comparisons of the numerical results with the experimental results of the laboratory model tests and the numerical analysis results of the conventional whole 2D simulation. In addition, the effect of a submerged structure at the opening of a breakwater is discussed from the numerical simulations by the hybrid model. The submerged structure improves the stability of the rubble mound and reduces the tsunami inflow into the bay, while it increases the water surface velocity around the opening of the breakwater. The increase of surface velocity causes the increases of impulsive forces by collision with drafts and so on.     

Real Time Analysis and Forecasting of Strata Caving Behaviour during Longwall Operations

Summary. Discontinuous manual observations and irregular caving characteristics of roof rocks often lead to improper decisions resulting in accidents and production loss. Hence, systematic monitoring of the hanging roof behind the chock shields is necessary for safe and productive mining operations. A real-time application was successfully implemented in an Indian mine for forecasting of hanging roof behaviour to enhance safety and productivity. This paper reports the functioning of real-time TWAP (time weighted average pressure) analysis in the forecasting of hanging roof behaviour in real time.