A stabilized Smoothed Particle Hydrodynamics,Taylor–Galerkin algorithm for soil dynamics problems

Modelling of failure under dynamic conditions in geomaterials with finite elements presents a series of complex problems, among which we can mention those of (i) volumetric locking, which results on higher failure loads, (ii) influence of mesh alignment, resulting to unrealistic failure surfaces, (iii) diffusion of the shear band over some element widths, (iv) nonoptimal propagation properties (numerical diffusion and dispersion), (v) fulfilling Babuska–Brezzi conditions when using the same order of interpolation for displacement and pressures in coupled problems and (vi) large deformation analysis. This paper is based on previous work done by the authors, who developed a mixed approximation based on (i) casting the dynamic problem in the form of a system of first order PDEs and (ii) using stresses and velocities as nodal variables. The equations were discretized following a Taylor–Galerkin algorithm, first in time using a Taylor expansion and then in space using Galerkin method. The model was limited to small deformations. The purpose of this paper is to show how Taylor–Galerkin method can be extended to meshless formulations, such as the SPH method. The algorithm consists of (i) discretizing in time using a Taylor series expansion complemented with integration of source terms using a Runge–Kutta scheme and then (ii) discretizing in space using the SPH method. It is shown how the proposed method keeps the advantages of the Taylor–Galerkin method in Finite Elements (good propagation properties and capturing of shear bands) and avoids the tensile instability. A set of test problems ranging from elastic propagation of a wave in a bar to failure of a slope on a cohesive softening material are used to assess the performance of the method. Copyright © 2011 John Wiley & Sons, Ltd.     

Resolving cosmic structure formation with the Millennium-II Simulation

We present the Millennium-II Simulation (MS-II), a very large N -body simulation of dark matter evolution in the concordance Λ cold dark matter (ΛCDM) cosmology. The MS-II assumes the same cosmological parameters and uses the same particle number and output data structure as the original Millennium Simulation (MS), but was carried out in a periodic cube one-fifth the size  (100  h ⁻¹ Mpc)  with five times better spatial resolution (a Plummer equivalent softening of  1.0  h ⁻¹ kpc  ) and with 125 times better mass resolution (a particle mass of  6.9 × 10⁶  h ⁻¹ M_⊙  ). By comparing results at MS and MS-II resolution, we demonstrate excellent convergence in dark matter statistics such as the halo mass function, the subhalo abundance distribution, the mass dependence of halo formation times, the linear and non-linear autocorrelations and power spectra, and halo assembly bias. Together, the two simulations provide precise results for such statistics over an unprecedented range of scales, from haloes similar to those hosting Local Group dwarf spheroidal galaxies to haloes corresponding to the richest galaxy clusters. The 'Milky Way' haloes of the Aquarius Project were selected from a lower resolution version of the MS-II and were then resimulated at much higher resolution. As a result, they are present in the MS-II along with thousands of other similar mass haloes. A comparison of their assembly histories in the MS-II and in resimulations of 1000 times better resolution shows detailed agreement over a factor of 100 in mass growth. We publicly release halo catalogues and assembly trees for the MS-II in the same format within the same archive as those already released for the MS.     

Effects of seed particle size and composition on heterogeneous nucleation of n-nonane

Heterogeneous nucleation of supersaturated n-nonane vapour on seed particles of different size and composition has been investigated using a fast expansion chamber. Monodisperse seed particle sizes were ranging from about 4 nm up to about 24 nm in diameter. By using different types of particle generators WO_x, Ag and (NH₄)₂SO₄ particles were generated. For direct comparison between different particle compositions overlapping sizes have been generated for WO_x and Ag at about 7 nm particle diameter as well as for Ag and (NH₄)₂SO₄ at about 15 nm. Nucleation temperature was kept constant at about 278 K. Experimental data were compared to Kelvin equation and Fletcher theory including the effect of line tension. It was found that heterogeneous nucleation of n-nonane seems to be independent of seed particle composition and starts well below the Kelvin curve. Good agreement was achieved with Fletcher theory including the effect of line tension.     

The importance of three centuries of land-use change for the global and regional terrestrial carbon cycle

Large amounts of carbon (C) have been released into the atmosphere over the past centuries. Less than half of this C stays in the atmosphere. The remainder is taken up by the oceans and terrestrial ecosystems. Where does the C come from and where and when does this uptake occur? We address these questions by providing new estimates of regional land-use emissions and natural carbon fluxes for the 1700–2000 period, simultaneously considering multiple anthropogenic (e.g. land and energy demand) and biochemical factors in a geographically explicit manner. The observed historical atmospheric CO₂ concentration profile for the 1700 to 2000 period has been reproduced well. The terrestrial natural biosphere has been a major carbon sink, due to changes in climate, atmospheric CO₂, nitrogen and management. Due to land-use change large amounts of carbon have been emitted into the atmosphere. The net effect was an emission of 35 Pg C into the atmosphere for the 1700 to 2000 period. If land use had remained constant at its distribution in 1700, then the terrestrial C uptake would have increased by 142 Pg C. This overall difference of including or excluding land-use changes (i.e. 177 Pg C) comes to more than half of the historical fossil-fuel related emissions of 308 Pg C. Historically, global land-use emissions were predominantly caused by the expansion of cropland and pasture, while wood harvesting (for timber and fuel wood) only played a minor role. These findings are robust even when changing some of the important drivers like the extent of historical land-use changes. Under varying assumptions, land-use emissions over the past three centuries could have increased up to 20%, but remained significantly lower than from other sources. Combining the regional land-use and natural C fluxes, North America and Europe were net C sources before 1900, but turned into sinks during the twentieth century. Nowadays, these fluxes are a magnitude smaller than energy- and industry-related emissions. Tropical regions were C neutral prior to 1950, but then accelerated deforestation turned these regions into major C sources. The energy- and industry-related emissions are currently increasing in many tropical regions, but are still less than the land-use emissions. Based on the presented relevance of the land-use and natural fluxes for the historical C cycle and the significance of fossil-fuel emissions nowadays, there is a need for an integrated approach for energy, nature and land use in evaluating possible climate change mitigation policies.     

Changing hillslope and fluvial Holocene sediment dynamics in a Belgian loess catchment

Floodplain deposition is an essential part of the Holocene sediment dynamics of many catchments and a thorough dating control of these floodplain deposits is therefore essential to understand the driving forces of these sediment dynamics. In this paper we date floodplain and colluvial deposition in the Belgian Dijle catchment using accelerator mass spectrometric radiocarbon and optical stimulated luminescence dating. Relative mass accumulation curves for the Holocene were constructed for three colluvial sites and 12 alluvial sites. A database was constructed of all available radiocarbon ages of the catchment and this database was analysed using relative sediment mass accumulation rates and cumulative probability functions of ages and site‐specific sedimentation curves. Cumulative probability functions of ages were split into different depositional environments representing stable phases and phases of accelerated clastic deposition. The results indicate that there is an important variation between the different dated sites. After an initial stable early and middle Holocene phase with mainly peat growth in the floodplains, clastic sedimentation rates increased from 4000 BC on. This first phase was more pronounced and started somewhat earlier for colluvial deposits then for alluvial deposits. The main part of the Holocene deposits, both in colluvial and alluvial valleys, was deposited during the last 1 ka. The sedimentation pattern of the individual dated sites and the catchment‐wide pattern indicate that land use changes are responsible for the main variations in the Holocene sediment dynamics of this catchment, while the field data do not provide indications for a climatological influence on the sediment dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling mean annual sediment yield using a distributed approach

In this paper a spatially distributed model for the calculation of sediment delivery to river channels is presented (SEDEM: SEdiment DElivery Model). The model consists of two components: (1) the calculation of a spatial pattern of mean annual soil erosion rates in the catchment using a RUSLE (Revised Soil Erosion Equation) approach; and (2) the routing of the eroded sediment to the river channel network taking into account the transport capacity of each spatial unit. If the amount of routed sediment exceeds the local transport capacity, sediment deposition occurs. An existing dataset on sediment yield for 24 catchments in central Belgium was used to calibrate the transport capacity parameters of the model. A validation of the model results shows that the sediment yield for small and medium sized catchments (10–5000 ha) can be predicted with an average accuracy of 41 per cent. The predicted sediment yield values with SEDEM are significantly more accurate than the predictions using a lumped regression model. Moreover a spatially distributed approach allows simulation of the effect of different land use scenarios and soil conservation techniques. Copyright © 2001 John Wiley & Sons, Ltd.     

Crustal structure of northeastern Tibet inferred from receiver function analysis

Sponsored by National Science & Technology Committee, the cooperation between China Academy of Geoscience and Institute of Geophysics and Tectonics, University Joseph, France conducted a lithospherical experiment using 40 Minititan 3-component and 13 CEIS 1-component seismometers along the road from Gonghe to Yushu in Qinghai Province during 5 months after June, 1998. The interested area is on the north of Bangong- Nujiang fault, the east of Qaidam basin, the south of Center Qilian fault and the west of Longmenshan fault. And the profile across most tectonic parts of Eastern Tibet such as Southern Qilian, Eastern Kunlun fault, Bayan Har terrane, Jinshajiang suture (Figure 1), which is the first seismological profile across Eastern Tibet (Qinghai-Tibet) and will be beneficial on the comparison with the results of its center parts, especially on the understanding of the effect of the thousands-kilometer-faraway collision between Eurasia Plate and Indian Plate on the uplifting of south and north part of Eastern Kunlun fault, and on the thickening of crust and the feature of deep structure of Qilian mountain on the north of Tibet Plateau.     

The formulation of an effective demulsifier for oil spill emulsions

The formation of stable sea water-in-crude oil emulsions after an oil spill was first reported after the Torrey Canyon wreck in 1967. The problem of handling this stable water (65%)-in-oil (35%) emulsion rather than oil is evidenced in the additional expensive on-site storage requirement as well as the handling difficulty with the viscous semi-solid emulsion. This paper details the results of an R&D program that has established the mechanism for the formation of this emulsion. This insight enabled a specific effective chemical demulsifier to be formulated. A treatment technique was devised that entails injecting the demulsifier into the oil spill emulsion as it is skimmed from the water surface. The emulsion is then rapidly separated into water-free oil and a clean-water phase that can be decanted back into the sea. Full-scale tests were successfully conducted by JBF Scientific Corp. using two types of commercial JBF DIP skimmers in a small concrete lined pond.