Numerical study of the coupled hydro-mechanical effects in dynamic compaction of saturated granular soils

Dynamic compaction is a widely used method for improvement of loose granular deposits. Its applicability in saturated layers generally considered to be less effective because of the fact that part of the applied energy is absorbed by pore water. Up to now the majority of numerical simulations have focused on the analysis of dynamic compaction in dry/moist soils. In this paper, a fully coupled hydro-mechanical finite element code has been developed and employed to evaluate the dynamic compaction effects on saturated granular soils. After verification of the results by comparing the numerical results with those measured in a real field case of DC treatment in a highway, some sensitivity analyses have been performed to evaluate the effect of water phase on the dimensions of the zone of improvement in the soil beneath the tamper. The results indicate that in the DC process the soil demonstrate two different behaviors. At the very early stage after impact, the soil behaves in an undrained manner and high oscillation of pore pressure occurs. After this phase, consolidation begins during which the pore-water-flow out of the soil mass takes place. The numerical analysis reveals that most of the DC improvement occurs during the undrained phase. The main mechanism responsible for the densification of soil during the undrained phase seems to be the compressibility of pore water. The simulation results indicate that the improvement zone diminishes when the degree of saturation increases.     

Study of pore pressure variation during liquefaction using two constitutive models for sand

Numerical analyses of liquefiable sand are presented in this paper. Liquefaction phenomenon is an undrained response of saturated sandy soils when they are subjected to static or dynamic loads. A fully coupled dynamic computer code is developed to predict the liquefaction potential of a saturated sandy layer. Coupled dynamic field equations of extended Biot's theory with u–P formulation are used to determine the responses of pore fluid and soil skeleton. Generalized Newmark method is employed for integration in time. The soil behavior is modelled by two constitutive models; a critical state two-surface plasticity model, and a densification model. A class ‘B’ analysis of a centrifuge experiment is performed to simulate the dynamic response of level ground sites. The results of the numerical analyses demonstrate the capability of the critical sate two-surface plasticity model in producing pore pressures that are consistent with observations of the behavior of liquefiable sand in the centrifuge test.     

On the impact of atmospheric vs oceanic resolutions on the representation of the sea surface temperature in the South Eastern Tropical Atlantic

Climate Dynamics - Despite the efforts of the modelling community to improve the representation of the sea surface temperature (SST) over the South Eastern Tropical Atlantic, warm biases still...     

Small-Strain Shear Modulus of Cement-Admixed Kaolinite

An experimental study is conducted to measure small-strain shear modulus of clay-cement mixture using bender element apparatus setup in a triaxial cell. Bender element tests were conducted on cement-treated soils and the results were analyzed to study the variation of shear modulus properties of soil specimens at different cement contents, confining pressures, curing times, and compaction moisture contents. Based on the obtained results increasing the cement ratio has a significant effect on the small-strain shear modulus of the treated soils, and this effect signifies with increasing the moisture content and curing time. Rates of shear modulus enhancements due to cement content, curing time, and compaction moisture content are quantified and presented. In this study, a clay–cement–water ratio formulation is proposed that enables one to calculate cement and water contents required to obtain specific small-strain shear modulus.     

A major change in magma sources in late Mesozoic active margin of the circum-Sea of Japan domain: Geochemical constraints from late Paleozoic to Paleogene mafic dykes in the Sergeevka belt,southern Primorye,Russia

More than 30 mafic dykes crop out in the Sergeevka belt in the coastal South Primorye, Far East Russia, of which geologic settings have been unclear for years. This study conducted major- and trace elements characterization, Sr–Nd isotope analyses, and Ar–Ar amphibole and U–Pb zircon datings for these rocks in order to identify their origin. The results demonstrated that all dykes are characterized by high Ba/Yb and low Nb/Y, Zr/Y, and Th/Yb ratios, which suggest their origin from arc melts derived from thin wedge mantle and shallow-dipping slab. These dykes are clearly separated into two distinct age/geochemistry suites; that is, the Paleogene and Early Cretaceous one with dolerites/basalts and adakitic rocks, and the Permian–Triassic one with high-Mg and high-Al gabbro-dolerite varieties. Their geochemistry suggests that the older suite was sourced from a primitive depleted MORB mantle (DMM)-type mantle, whereas the younger suite from an enriched mantle II (EM2)-type mantle domain. The transition in source type from DMM to EM2 occurred during the Jurassic-earliest Cretaceous time, probably by a strong influence of a mantle plume onto the long-continuing subduction-related magmatism. The plume influence reached the maximum when the unique meimechite-picrite complex formed in the region.     

Lidar and Sun photometer observations of atmospheric boundary-layer characteristics over an urban area in a mountain valley

We present measurements of the vertical aerosol structure and the aerosol optical depth in the lower troposphere performed above the city of Sofia (an urban area situated in a mountain valley), western Bulgaria by means of a ground-based aerosol lidar operating continuously for a number of years. The lidar measurements were accompanied by measurements of the aerosol optical depth (AOD) in the visible and near infrared regions of the spectrum performed in October 2004 using Microtops II radiometers. The maximum values of the AOD were found to occur 1–2 h before the complete development of the atmospheric boundary layer, i.e. during the residual layer destruction, which confirms our hypothesis concerning the slope circulation effect on the processes taking place in the atmospheric boundary layer. The AOD values obtained by the lidar are lower than those taken by the sun photometer. Further, the AOD exhibits two different types of behaviour. In the case of a ‘clear atmosphere’ (i.e. in the absence of volcanic eruptions and/or dust transport from the Sahara) most of the aerosol accumulated within the atmospheric boundary layer over the urban area considered. The combined use of the two instruments allows the comparison between the optical characteristics of the atmospheric aerosol (e.g. aerosol extinction coefficient, etc.) obtained by the lidar and through an independent method (sun photometer).     

Seasonal, interannual and long-term variability of precipitation and snow depth in the region of the Barents and Kara seas

Observation data of temperature, precipitation and snow depth have been compiled and generalized climatologically for a network of 38 stations in and around the Barents and Kara seas, for the period 1951–1992. The monthly precipitation totals were corrected for measuring errors, and the correction method is described in detail. The corrected precipitation values show that the annual precipitation in the region ranges from more than 500 mm along the coast of the Kola Peninsula to less than 200 mm in parts of the north-eastern Kara Sea. The solid fraction of the annual precipitation ranges from 70% in northern parts to 35% in southern parts. For the period 1951–1992 the analysis indicates decreasing trends in annual values of temperature, precipitation and snow depths in the north-eastern parts of the region.     

Element partitioning between immiscible borosilicate liquids: A high-temperature centrifuge study

The main objectives of this study were to investigate conditions for stable and metastable liquid immiscibility in dry borosilicate synthetic systems and to evaluate effects of temperature and bulk melt composition on two-liquid element partitioning and boron speciation. To distinguish between the stable immiscibility and quench heterogeneity, we used high-temperature centrifuge phase separation. For the case of stable liquid immiscibility, silica-rich (LS) and borate-rich (LB) conjugate liquids formed two distinct layers separated by a sharp meniscus. The liquids were quenched into glasses, which were analysed by electron microprobe. Some of the glasses were also studied by Raman spectroscopy. We used several synthetic mixtures along the danburite-anorthite (CaB₂Si₂O₈-CaAl₂Si₂O₈) and danburite-reedmergnerite (CaB₂Si₂O₈-NaBSi₃O₈) joins. In addition, we studied four complex, six-component, Mg-bearing compositions with variable Na₂O and Al₂O₃ contents. The experiments show that the width of the LS-LB miscibility gap decreases more rapidly with the B-Al substitution (in the danburite-anorthite join) than with the Ca-Na substitution, implying that interactions between network-forming elements have a greater effect on borate-silicate unmixing than the nature of network-modifying cations. Ca and Mg partition strongly to the depolymerised borate-rich liquid with LB-LS partition coefficients of ∼40 and higher. On the other hand, two-liquid partition coefficients of Na and Al in most cases are close to 1 and show complex variations with temperature and bulk melt composition. Raman spectra of LB glasses quenched at different temperatures suggest that the proportion of trigonal boron in bulk boron content decreases with decreasing temperature. The change in boron speciation appears to affect Al and Na two-liquid partitioning in such a way that at low temperatures, the latter element becomes more compatible with LS.     

Multiphase flow modeling with density functional method

This paper is a review of applications of density functional theory (DFT) in compositional hydrodynamics. The basic idea is representation of the entropy or the Helmholtz energy of the mixture as the functional depending on the molar densities of chemical components (density functional). The hydrodynamics is governed by local conservation laws of chemical components, momentum, and energy, while constitutive relations and boundary conditions are introduced in accordance with the explicit form of the density functional. The general ideas and the history of the DFT in compositional hydrodynamics are discussed. Then the DFT for multiphase multicomponent mixtures is presented including the exposition of the first principles, governing equations and constitutive relations, and explicit expressions of density functional depending on physical situation. The DFT-based numerical simulator is described, and several multiphase simulation results are presented to illustrate the scope and effectiveness of DFT: sessile drop with and without surfactant, droplet breakup in shear flow, and three-phase hydrodynamics with mobile solid phase. Also, two practical scenarios with multiphase simulations in micro-CT porous rock models are presented: two-phase immiscible water-oil flow and three-phase water-gas-condensate flow with phase transitions. All numerical results are obtained by essentially the same code; both the number of chemical components and the Helmholtz energy have been set up in accordance with physical situation.