Propagation of seismic waves through liquefied soils

To predict the earthquake response of saturated porous media it is essential to correctly simulate the generation, redistribution, and dissipation of excess pore water pressure during and after earthquake shaking. To this end, a reliable numerical tool requires a dynamic, fully coupled formulation for solid–fluid interaction and a versatile constitutive model. Presented in this paper is a 3D finite element framework that has been developed and utilized for this purpose. The framework employs fully coupled dynamic field equations with a u–p–U formulation for simulation of pore fluid and solid skeleton interaction and a SANISAND constitutive model for response of solid skeleton. After a detailed verification and validation of the formulation and implementation of the developed numerical tool, it is employed in the seismic response of saturated porous media. The study includes examination of the mechanism of propagation of the earthquake-induced shear waves and liquefaction phenomenon in uniform and layered profiles of saturated sand deposits.     

A destructuration theory and its application to SANICLAY model

Many natural clays have an undisturbed shear strength in excess of the remoulded strength. Destructuration modeling provides a means to account for such sensitivity in a constitutive model. This paper extends the SANICLAY model to include destructuration. Two distinct types of destructuration are considered: isotropic and frictional. The former is a concept already presented in relation to other models and in essence constitutes a mechanism of isotropic softening of the yield surface with destructuration. The latter refers to the reduction of the critical stress ratio reflecting the effect of destructuration on the friction angle, and is believed to be a novel proposition. Both the types depend on a measure of destructuration rate expressed in terms of combined plastic volumetric and deviatoric strain rates. The SANICLAY model itself is generalized from its previous form by additional dependence of the yield surface on the third isotropic stress invariant. Such a generalization allows to obtain as particular cases simplified model versions of lower complexity including one with a single surface and associative flow rule, by simply setting accordingly parameters of the generalized version. A detailed calibration procedure of the relatively few model constants is presented, and the performance of three versions of the model, in descending order of complexity, is validated by comparison of simulations to various data for oedometric consolidation followed by triaxial undrained compression and extension tests on two structured clays. Copyright © 2009 John Wiley & Sons, Ltd.     

Gulf Stream evolution through the Straits of Florida: the role of eddies and upwelling near Cuba

Ocean Dynamics - The meandering of the Gulf Stream through the Straits of Florida is associated with eddy activity to the north (along the Florida Keys) and the south (along the Cuban coast). This...     

Mixing efficiency from microstructure measurements in the Sicily Channel

Ocean Dynamics - The dissipation flux coefficient, a measure of the mixing efficiency of a turbulent flow, was computed from microstructure measurements collected with a vertical microstructure...     

Modeling cyclic shearing of sands in the semifluidized state

Liquefaction is associated with the loss of mean effective stress and increase of the pore water pressure in saturated granular materials due to their contractive tendency under cyclic shear loading. The loss of mean effective stress is linked to loss of grain contacts, bringing the granular material to a “semifluidized state” and leading to development and accumulation of large cyclic shear strains. Constitutive modeling of the cyclic stress-strain response in earthquake-induced liquefaction and post-liquefaction is complex and yet very important for stress-deformation and performance-based analysis of sand deposits. A new state internal variable named strain liquefaction factor is introduced that evolves at low mean effective stresses, and its constitutive role is to reduce the plastic shear stiffness and dilatancy while maintaining the same plastic volumetric strain rate in the semifluidized state. This new constitutive ingredient is added to an existing critical state compatible, bounding surface plasticity reference model, that is well established for constitutive modeling of cyclic response of sands in the pre-liquefaction state. The roles of the key components of the proposed formulation are examined in a series of sensitivity analyses. Their combined effects in improving the performance of the reference model are examined by simulating undrained cyclic simple shear tests on Ottawa sand, with focus on reproducing the increasing shear strain amplitude as well as its saturation in the post-liquefaction response.     

A critical state sand plasticity model accounting for fabric evolution

Fabric and its evolution need to be fully considered for effective modeling of the anisotropic behavior of cohesionless granular sand. In this study, a three‐dimensional anisotropic model for granular material is proposed based on the anisotropic critical state theory recently proposed by Li & Dafalias [2012], in which the role of fabric evolution is highlighted. An explicit expression for the yield function is proposed in terms of the invariants and joint invariants of the normalized deviatoric stress ratio tensor and the deviatoric fabric tensor. A void‐based fabric tensor that characterizes the average void size and its orientation of a granular assembly is employed in the model. Upon plastic loading, the material fabric is assumed to evolve continuously with its principal direction tending steadily towards the loading direction. A fabric evolution law is proposed to describe this behavior. With these considerations, a non‐coaxial flow rule is naturally obtained. The model is shown to be capable of characterizing the complex anisotropic behavior of granular materials under monotonic loading conditions and meanwhile retains a relatively simple formulation for numerical implementation. The model predictions of typical behavior of both Toyoura sand and Fraser River sand compare well with experimental data. Copyright © 2013 John Wiley & Sons, Ltd.     

The role of the West Florida Shelf topography on the Loop Current system variability

Ocean Dynamics - The evolution of the Loop Current (LC) system under the interaction with the complex topography of the Gulf of Mexico (GoM) is examined. Focusing on the eastern GoM, we study the...     

On the processes that influence the transport and fate of Mississippi waters under flooding outflow conditions

The Mississippi River (MR) freshwater outflow is a major circulation forcing mechanism for the Northern Gulf of Mexico. We investigate the transport and fate of the brackish waters under flood conditions. The largest outflow in history (45,000 m³/s in 2011) is compared with the second largest outflow in the last 8 years (41,000 m³/s in 2008). Realistically forced simulations reveal the synergistic effect of enhanced discharge, winds, stratification of ambient shelf waters, and offshore circulation over the transport of plume waters. The strongest impact is attributed to the evolution of the Loop Current (LC) and associated frontal cyclonic eddies and anticyclonic rings, which exhibited distinctly different influence during the two study periods. The northward LC intrusion in the summer of 2011 weakened and blocked the buoyancy-driven downstream (westward) transport of brackish waters. The 2011 flood was thus characterized by upstream (eastward) flow and an extensive coverage of the Mississippi–Alabama–Florida shelf. An immediate response between the LC and the brackish offshore eastward spreading is computed during and after this historic event. The absence of a LC northward intrusion during the 2008 flood, in combination with wind effects, promotes downstream advection of MR waters towards the Louisiana–Texas shelf; large amounts of buoyant waters are also retained near the Delta, subject to local offshore advection under the synergistic action of LC-associated counter-rotating eddies.