Uncoupling of coupled flows in soil—a finite element method

Coupled flow of water, chemicals, heat and electrical potential in soil are of significance in a variety of circumstances. The problem is characterized by the coupling between different flows, i.e. a flow of one type driven by gradients of other types, and by the dual nature of certain flows, i.e. combined convection and conduction. Effective numerical solutions to the problem are challenged due to the coupling and the dual nature. In this paper, we first present a general expression that can be used to represent various types of coupled flows in soil. A finite element method is then proposed to solve the generalized coupled flows of convection-conduction pattern. The unknown vector is first decomposed into two parts, a convective part forming a hyperbolic system and a conductive part forming a parabolic system. At each time step, the hyperbolic system is solved analytically to give an initial solution. To solve the multi-dimensional hyperbolic system, we assume that a common eigenspace exists for the coefficient matrices, so that the system can be uncoupled by transforming the unknown vector to the common eigenspace. The uncoupled system is solved by the method of characteristics. Using the solution of the hyperbolic system as the initial condition, we then solve the parabolic system by a Galerkin finite element method for space discretization and a finite difference scheme for time stepping. The proposed technique can be used for solving multi-dimensional, transient, coupled or simultaneous flows of convection-conduction type. Application to a flow example shows that the technique indeed exhibits optimality in convergence and in stability.     

Simulation of individual particle movement in a gravel streambed

A two‐dimensional simulation model of travel distances of individual particles in a gravel‐bed river is presented. The model is based on a number of rules, which include particle size, entrainment, trajectory, distance of movement and entrapment. Particle interactions are controlled by resistance fields defined about each obstacle and critical elevation defined in the model. Resistance fields, particle dropping and critical elevation rules control particle interactions. The interaction rules cause the particles to develop pebble clusters, stone cells and transverse structures (transverse ribs). The simulated travel distances of individual particles are consistent with reported field results. Individual particle travel distances were simulated using two different models; one without interactions between the individual particles and the stationary bed and one with interactions. The case without interactions demonstrates the random nature of sediment transport, and narrow ranges of travel distances. Wider ranges of travel distances, similar to those for natural situations, were obtained for the cases with interactions. The more intense the interaction between the mobile stones and the stationary ones, the wider the range of distances of travel for a given particle size. Modelling the mean travel distance yielded a result similar to that published previously, which was based on empirical data. Well developed bed‐surface structures were obtained for relatively poorly sorted sediment with intense interactions between particles. Transverse structures developed when relatively large particles were allowed to move. Copyright © 2002 John Wiley & Sons, Ltd.     

Nonlinear interaction between acoustic gravity waves

The resonant interaction between three acoustic gravity waves is considered. We improve on the results of previous authors and write the new coupling coefficients in a symmetric form. Particular attention is paid to the low-frequency limit.     

Analyses of trace elements on quartz surfaces in sulfidic mine tailings from Kristineberg (Sweden) a few years after remediation

Laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) has been applied to determine the elemental composition of the surface layer, as well as of the first interior layer, of quartz grains from the mine tailings from Kristineberg (northern Sweden) in order to determine concentration gradients between these two layers. The quartz grains were collected from the oxidized and unoxidized zones within the tailings. The aim of this study is to assess the role of quartz surfaces as sites for the attenuation of solutes from the mine-tailings porewater. Concentrations of Cu, Ag, Sb, Pb and Bi are highest near the surface of each grain and decrease towards the interior. The surface concentration of Cu, Zn and Pb is more pronounced within the unoxidized than within the oxidized zone of the tailings. Cu exhibits a distinct concentration peak at the surface of the quartz grains below the pre-remediation oxidation front. For Zn and Ce the trend of high surface concentration is less pronounced than for Cu or Pb. Silver, Bi and As are preferably adsorbed within the uppermost layers of the oxidized zone where the pH is as high as 6.2. The conversion of intensity signals of the elements to concentration values in ppm was done by using external standards (NIST silicate glass).     

The use of press data in the development of a database for rainfall-induced landslides in Tegucigalpa,Honduras, 1980–2005

The capital city of Honduras, Tegucigalpa, suffers from the occurrence of destructive landslides on a yearly basis. During the rainy season, damages to infrastructure as well as injuries, casualties and homeless individuals resulting from landslides are reported in the press. This paper presents the development of a database for rainfall-induced landslides for the period 1980–2005, based on the news reported by two local newspapers. The editions comprehended during the study period have been scrutinized, and articles focusing on landslides, tropical storms, hurricanes, floods and vulnerability of the city have been collected. The interpretation of these archives has allowed the compilation of valuable data of approximately 400 landslides. The analyses of monthly and annual precipitation during the study period show how extreme rainfall events like Hurricane Mitch in October of 1998 have significantly contributed to the initiation of landslides. In addition, the assessment of the slums and neighborhoods affected by landslides during the study period reveals an evident link between the social and physical vulnerability of Tegucigalpa. In order to estimate the reliability of this press-based database, the set of landslides that have been reported as a result of Hurricane Mitch in the press archives has been compared with two inventories based on the interpretation of aerial photographs taken in 1999 and 2001. It is shown that the analysis of the landslide damage left after the hurricane can be enriched with the detailed temporal data provided in the archives and the precise location of these events determined by the aerial photographs. Despite the difficulties faced in the compilation of this database, a good comprehension of the temporal and spatial distribution of landslides in Tegucigalpa has been achieved.     

Optimization routine for identification of model parameters in soil plasticity

The paper presents an optimization routine especially developed for the identification of model parameters in soil plasticity on the basis of different soil tests. Main focus is put on the mathematical aspects and the experience from application of this optimization routine. Mathematically, for the optimization, an objective function and a search strategy are needed. Some alternative expressions for the objective function are formulated. They capture the overall soil behaviour and can be used in a simultaneous optimization against several laboratory tests. Two different search strategies, Rosenbrock's method and the Simplex method, both belonging to the category of direct search methods, are utilized in the routine. Direct search methods have generally proved to be reliable and their relative simplicity make them quite easy to program into workable codes. The Rosenbrock and simplex methods are modified to make the search strategies as efficient and user‐friendly as possible for the type of optimization problem addressed here. Since these search strategies are of a heuristic nature, which makes it difficult (or even impossible) to analyse their performance in a theoretical way, representative optimization examples against both simulated experimental results as well as performed triaxial tests are presented to show the efficiency of the optimization routine. From these examples, it has been concluded that the optimization routine is able to locate a minimum with a good accuracy, fast enough to be a very useful tool for identification of model parameters in soil plasticity. Copyright © 2001 John Wiley & Sons, Ltd.     

Growth Rates of Global Energy Systems and Future Outlooks

The world is interconnected and powered by a number of global energy systems using fossil, nuclear, or renewable energy. This study reviews historical time series of energy production and growth for various energy sources. It compiles a theoretical and empirical foundation for understanding the behaviour underlying global energy systems’ growth. The most extreme growth rates are found in fossil fuels. The presence of scaling behaviour, i.e. proportionality between growth rate and size, is established. The findings are used to investigate the consistency of several long-range scenarios expecting rapid growth for future energy systems. The validity of such projections is questioned, based on past experience. Finally, it is found that even if new energy systems undergo a rapid ‘oil boom’-development—i.e. they mimic the most extreme historical events—their contribution to global energy supply by 2050 will be marginal.     

The aperiodic broad-band X-ray variability of Cygnus X-3

We study the soft X-ray variability of Cygnus X-3. By combining data from the All-Sky Monitor and Proportional Counter Array instruments on the RXTE satellite with EXOSAT /Medium Energy (ME) detector observations, we are able to analyse the power density spectrum (PDS) of the source from 10⁻⁹ to 0.1 Hz, thus covering time-scales from seconds to years. As the data on the longer time-scales are unevenly sampled, we combine traditional power spectral techniques with simulations to analyse the variability in this range. The PDS at higher frequencies  (≳10⁻³ Hz)  are for the first time compared for all states of this source. We find that it is for all states well described by a power law, with index  ∼−2  in the soft states and a tendency for a less steep power law in the hard state. At longer time-scales, we study the effect of the state transitions on the PDS, and find that the variability below  ∼10⁻⁷ Hz  is dominated by the transitions. Furthermore, we find no correlation between the length of a high/soft-state episode and the time since the previous high/soft state. On intermediate time-scales, we find evidence for a break in the PDS at time-scales of the order of the orbital period. This may be interpreted as evidence for the existence of a tidal resonance in the accretion disc around the compact object, and constraining the mass ratio to   M ₂/ M ₁≲ 0.3  .     

Asteroid orbital ranging using Markov‐Chain Monte Carlo

Abstract— We present a novel Markov‐Chain Monte‐Carlo orbital ranging method (MCMC) for poorly observed single‐apparition asteroids with two or more observations. We examine the Bayesian a posteriori probability density of the orbital elements using methods that map a volume of orbits in the orbital‐element phase space. In particular, we use the MCMC method to sample the phase space in an unbiased way. We study the speed of convergence and also the efficiency of the new method for the initial orbit computation problem. We present the results of the MCMC ranging method applied to three objects from different dynamical groups. We conclude that the method is applicable to initial orbit computation for near‐Earth, main‐belt, and transneptunian objects.