Explicit solutions for the instantaneous undrained contraction of hollow cylinders and spheres in porous elastoplastic medium

In this article we present closed‐form solutions for the undrained variations in stress, pore pressure, deformation and displacement inside hollow cylinders and hollow spheres subjected to uniform mechanical pressure instantaneously applied to their external and internal boundary surfaces. The material is assumed to be a saturated porous medium obeying a Mohr–Coulomb model failure criterion, exhibiting dilatant plastic deformation according to a non‐associated flow rule which accounts for isotropically strain hardening or softening. The instantaneous response of a porous medium submitted to an instantaneous loading is undrained, i.e. without any fluid mass exchange. The short‐term equilibrium problem to be solved is now formally identical to a problem of elastoplasticity where the constitutive equations involve the undrained elastic moduli and particular equivalent plastic parameters. The response of the model is presented (i) for extension and compression undrained triaxial tests, and (ii) for unloading problems of hollow cylinders and spheres through the use of appropriately developed closed‐form solutions. Numerical results are presented for a plastic clay stone with strain hardening and an argilite with strain softening. The effects of plastic dilation, of the strain softening law and also of geometry of the cavity on the behaviour of the porous medium have been underlined. Analytical solutions provide valuable benchmarks enabling various numerical methods in undrained conditions with a finite boundary to be verified. Copyright © 2002 John Wiley & Sons, Ltd.     

Kinematic wave solutions for pollutant transport over an infiltrating plane with finite‐period mixing and mixing zone

Kinematic wave solutions are derived for transport of a conservative non‐point‐source pollutant during a rainfall‐runoff event over an infiltrating plane for two cases: (i) finite‐period mixing and (ii) soil‐mixing zone. Rainfall is assumed to be steady, uniform and finite in duration, and it is assumed to have zero concentration of pollutants. Infiltration is assumed constant in time and space. Prior to the start of rainfall, the pollutant is distributed uniformly over the plane. In the first case, when rainfall occurs, the mixing of pollutant in the runoff water occurs in a finite period of time. In the second case, the chemical concentration is assumed to be a linearly decreasing function of rainfall intensity and overland flow. The solute concentration and discharge are found to depend on the flow characteristics as well as the solute concentration parameters. The characteristics of solute concentration and discharge graphs seem to be similar to those reported in the literature and observed in laboratory experiments. Copyright © 2002 John Wiley & Sons, Ltd.     

The effect of a dense atmosphere on the tidally induced potential of Titan

The effect of the dense atmosphere of Titan on the tidal variations of the external gravitational potential of degree two is quantified. The atmospheric tides perturb the external gravitational potential of Titan in two ways. First, the atmosphere itself contributes directly to the external gravitational potential with a period of 15.945 days. Second, the variable loading of the atmosphere induces mass redistribution within Titan, which also changes the external gravitational potential. It is shown that the relative atmospheric contributions to the tides are most likely less than 2% and vanish almost completely for the most plausible models with a subsurface ocean. This suggest that atmospheric tidal perturbations will contribute only negligibly to Cassini measurements of Titan's gravitational field so that the tidal Love numbers derived from these observations can be directly interpreted in terms of the satellite's interior.     

On Motions Near the Lagrange Equilibrium Point L 4 in the Case of Routh's Critical Mass Ratio

We deal with the planar restricted circular problem of three bodies. We study trajectories in a small neighborhood of the Lagrange equilibrium point L ₄ when mass ratio is close to Routh's value. In particular, we show that the case of proper degeneracy takes place and for most initial conditions trajectories are conditionally-periodic. We obtain an approximate representation of families of periodic solution emanating from the equilibrium point L ₄. We also show that in the case of instability of L ₄ the trajectories starting in a vicinity of L ₄ remain in a finite domain forever. We give an upper bound of this domain. To carry out our investigation, we analyze the dynamics of a general Hamiltonian system with two degrees of freedom in the case of 1 : 1 resonance in detail.This revised version was published online in October 2005 with corrections to the Cover Date.     

An analytical assessment on the impact of covers on the onset of air convection in mine wastes

If a mine waste pile is left open, an active chemical reaction of oxidation is often found due to the commonly high content of pyritic materials. The oxidation of pyrites is an exothermic process and the released heat will promote the flow of fresh oxygen from the surrounding atmosphere into the waste dump. As a result, oxidation reaction will accelerate and temperature within the dump can increase to as high as 60°C above the ambient temperature. The oxidation process also releases sulphuric acid and hydrogen ions into ground water to cause water contamination. Low‐permeability covers such as clay liners have been recently proposed to abate the oxidation process in mine wastes. The effectiveness of using low‐permeability materials to cover mine wastes in order to suppress the pyrite oxidation is examined. By conducting the theoretical analysis of the onset of convective air flow within waste rocks, the conditions under which soil gas flow is significant are identified. By comparing the results with previous field measurements and theoretical analysis for the uncovered conditions, it is shown that low‐permeability covers can effectively suppress soil gas flow and slow down the pyrite oxidation process in mine wastes. Copyright © 2001 John Wiley & Sons, Ltd.     

The general solution of the HENON–HEILES problem

The general solution of the Henon–Heiles system is approximated inside a domain of the (x, C) of initial conditions (C is the energy constant). The method applied is that described by Poincaré as ‘the only “crack” permitting penetration into the non-integrable problems’ and involves calculation of a dense set of families of periodic solutions that covers the solution space of the problem. In the case of the Henon–Heiles potential we calculated the families of periodic solutions that re-enter after 1–108 oscillations. The density of the set of such families is defined by a pre-assigned parameter ε (Poincaré parameter), which ascertains that at least one periodic solution is computed and available within a distance ε from any point of the domain (x, C) for which the approximate general solution computed. The approximate general solution presented here corresponds to ε = 0.07. The same solution is further improved by “zooming” into four square sub-domain of (x, C), i.e. by computing sufficient number of families that reduce the density parameter to ε = 0.003. Further zooming to reduce the density parameter, say to ε = 10⁻⁶, or even smaller, although easily performable in both areas occupied by stable as well as unstable solutions, was found unnecessary. The stability of all members of each and all families computed was calculated and presented in this paper for both the large solution domain and for the sub-domains. The correspondence between areas of the approximate general solution occupied by stable periodic solutions and Poincaré sections with well-aligned section points and also correspondence between areas occupied by unstable solutions and Poincaré sections with randomly scattered section points is shown by calculating such sections. All calculations were performed using the Runge-Kutta (R-K) 8th order direct integration method and the large output received, consisting of many thousands of families is saved as “Atlas of the General Solution of the Henon–Heiles Problem,” including their stability and is available at request. It is concluded that approximation of the general solution of this system is straightforward and that the chaotic character of its Poincaré sections imposes no limitations or difficulties.     

Analytical solutions for the equations of motion of a space vehicle during the atmospheric re-entry phase on a 2-D trajectory

A practical and important problem encountered during the atmospheric re-entry phase is to determine analytical solutions for the space vehicle dynamical equations of motion. The author proposes new solutions for the equations of trajectory and flight-path angle of the space vehicle during the re-entry phase in Earth’s atmosphere. Explicit analytical solutions for the aerodynamic equations of motion can be effectively applied to investigate and control the rocket flight characteristics. Setting the initial conditions for the speed, re-entering flight-path angle, altitude, atmosphere density, lift and drag coefficients, the nonlinear differential equations of motion are linearized by a proper choice of the re-entry range angles. After integration, the solutions are expressed with the Exponential Integral, and Generalized Exponential Integral functions. Theoretical frameworks for proposed solutions as well as, several numerical examples, are presented.     

The librations, shape, and icy shell of Europa

Europa, the smallest of the Galilean satellites, has a young icy surface and most likely contains an internal ocean. The primary objective of possible future missions to Europa is the unambiguous detection and characterization of a subsurface ocean. The thickness of the overlying icy shell provides important information on the thermal evolution of the satellite and on the interaction between the ocean and the surface, the latter being fundamental for astrobiology. However, the thickness is not well known, and estimates range from several hundred of meters to some ten of kilometers. Here, we investigate the use of libration (rotation variation) observations to study the interior structure of Europa and in particular its icy shell. A dynamical libration model is developed, which includes gravitational coupling between the icy shell and the heavy solid interior. The amplitude of the main libration signal at 3.55 days (the orbital period) is shown to depend on Europa's shape and structure. Models of the interior structure of Europa are constructed and the equatorial flattening of the internal layers, which are key parameters for the libration, are calculated by assuming that Europa is in hydrostatic equilibrium. Europa's flattened shape is assumed to be due to rotation and permanent tides, and we extend the classical Radau equation for rotationally flattened bodies to include also tidal deformation. We show that the presence of an ocean increases the amplitude of libration by about 10%, depending mainly on the thickness of the icy shell. Therefore, libration observations offer possibility of detection of a subsurface ocean in Europa and estimation of the thickness of its overlying icy shell.