Micromechanical Aspects of Liquefaction-Induced Lateral Spreading

This paper reports the results of model-based simulations of 1-g shake table tests of level and sloping saturated granular soils subject to seismic excitations. The simulations utilize a transient fully coupled continuum-fluid discrete-particle model of water-saturated soils. The fluid (water) phase is idealized at a mesoscale using an averaged form of Navier-Stokes equations. The solid particles are modeled at the microscale as an assemblage of discrete spheres using the discrete element method (DEM). The interphase momentum transfer is accounted for using an established relationship. The employed model reproduced a number of response patterns observed in the 1-g experiments. In addition, the simulation results provided valuable information on the mechanics of liquefaction initiation and subsequent occurrence of lateral spreading in sloping ground. Specifically, the simulations captured sliding block failure instances at different depth locations. The DEM simulation also quantified the impact of void redistribution during shaking on the developed water pressure and lateral spreading. Near the surface, the particles dilated and produced an increase in volume, while the particles at deeper depth locations experienced a decrease in volume during shaking.     

Site Investigations for Involvement of Water Films in Lateral Flow in Liquefied Ground

Previous research indicates that if layered sand deposits are liquefied during earthquakes, water films are likely to develop beneath less permeable sublayers and lead to the destabilization of sloping ground. In Niigata City, large lateral flow displacements were reported in almost flat areas during the 1964 Niigata earthquake. The involvement of water films in lateral flow failure during the earthquake is examined in this research based on site investigation data. Soil profiles in the investigated areas estimated from many borehole logs indicate that continuous or partially continuous sublayers of fine soil that cap liquefiable loose sand exist. Elevation contours of 0.1 m increments are drawn based on an in situ leveling survey and local maps. The ground slopes obtained are found to be closely related to flow displacements evaluated in previous research, indicating that a gentle slope of less than 1% results in displacement of several meters. This strongly suggests that water films with literally no shear resistance formed beneath fine soil sublayers were highly responsible for the large lateral flow displacements in these areas during the Niigata earthquake.     

Numerical Modeling of Abutment Scour with the Focus on the Incipient Motion on Sloping Beds

A three-dimensional computational fluid dynamics model is applied to predict local scour around an abutment in a rectangular laboratory flume. When modeling local scour, steep bed slopes up to the angle of repose occur. To predict the depth and the shape of the local scour correctly, the reduction of the critical shear stress due to the sloping bed must be taken into account. The focus of this study is to investigate different formulas for the threshold of noncohesive sediment motion on sloping beds. Some formulas only take the transversal angle (perpendicular to the flow direction) into account, but others also consider the longitudinal angle (streamwise direction). The numerical model solves the transient Reynolds-averaged Navier-Stokes equations in all three dimensions to compute the water flow. Sediment continuity in combination with an empirical formula is used to capture the bed load transport and the resulting bed changes. When the sloping bed exceeds the angle of repose, the bed slope is corrected with a sand-slide algorithm. The results from the numerical simulations are compared with data from physical experiments. The reduction of the bed shear stress on the sloping bed improves the results of the numerical simulation distinctly. The best results are obtained with the formulas that use both the transversal and the longitudinal angle for the reduction of the critical bed shear stress.     

Deformation of Pores in Viscoplastic Soil Material

Changes in soil pore volume and shape in response to internal and external mechanical stresses alter key soil hydrologic and transport properties. The extent of these changes is dependent on details of pore shape and size evolution. We present a model for quantifying rates of deformation and shape evolution of idealized spheroidal pores as functions of macroscopic stresses and soil rheological properties. Previous solutions for shrinkage of spherical pores embedded in a viscoplastic matrix under isotropic stress were extended to spheroidal pore shapes and biaxial stresses using Eshelby’s classical theory. Bulk soil behavior was obtained from upscaling of detailed single pore deformation. Results show that pore closure rates increase with decreasing initial aspect ratio (i.e., oblate pores close faster than spherical pores), and with higher deviatoric stress. Incomplete pore closure is attributed to soil hardening due to pore shape accommodation under biaxial stresses. The model provides a means for approximating pore deformation as input to predictive models for soil hydraulic properties.     

Validation of a Three-Dimensional Numerical Code in the Simulation of Pseudo-Natural Meandering Flows

Validation of a three-dimensional finite volume code solving the Navier–Stokes equations with the standard k-ε turbulence model is conducted using a high quality and high spatial resolution data set. The data set was collected from a large-scale meandering channel with a self-formed fixed bed, and comprises detailed bed profiling and laser Doppler anemometer velocity measurements. Comparisons of the computed primary and secondary velocities are made with those observed and it is found that the lateral momentum transfer is generally under predicted. At the apices this results in the predicted position of the primary velocity maximum having a bias towards the channel center, compared to the position where it has been measured. Using a simplified two zone roughness distribution whereby a separate roughness height was prescribed for the channel center and channel sides relative to a single distributed roughness height, generally led to a slightly improved longitudinal velocity distribution; the higher velocities were located nearer to the outside of the bend. Improving both the free surface calculation and scheme for discretization of the convection terms led to no appreciable difference in the computed velocity distributions. A more detailed study involving turbulence measurements and bed form height distribution should discriminate whether using distributed roughness height is a precursor to using an anisotropic turbulence representation for the accurate prediction of three-dimensional river flows.     

三层桨搅拌槽内三维流场的数值模拟

采用计算流体动力学(CFD)的方法,对稀土萃取过程中上两层为平直叶、底层为涡轮桨叶的三层组合桨搅拌槽内三维流场进行了研究。利用标准的k-epsilon双方程模型对无机相(水)和有机相(P507)的混合液在搅拌槽中产生的流场进行数值计算,得到这种搅拌桨以恒定转速300r/min在搅拌槽内转动时产生的速度场和压力场,以及速度分布云图、速度矢量图以及压力云图,为搅拌桨的设计与改进提供理论基础。     

Mechanics of Lateral Spreading Observed in a Full-Scale Shake Test

This paper examines in detail the mechanics of lateral spreading observed in a full-scale test of a sloping saturated fine sand deposit, representative of liquefiable, young alluvial and hydraulic fill sands in the field. The test was conducted using a 6-m tall inclined laminar box shaken at the base. At the end of shaking, nearly the whole deposit was liquefied, and the ground surface displacement had reached 32 cm. The presented analysis of lateral spreading mechanics utilizes a unique set of lateral displacement results, DH, from three independent techniques. One of these techniques—motion tracking analysis of the experiment video recording—is especially useful as it produced DH time histories for all laminar box rings and a complete picture of the lateral spreading initiation with an unprecedented degree of resolution in time and space. A systematic study of the data identifies the progressive stages of initiation and accumulation of lateral spreading, lateral spread contribution of various depth ranges and sliding zones, their relation to the simultaneous pore pressure buildup, and the soil shear strength response during sliding.     

Installation of Torpedo Anchors: Numerical Modeling

Torpedo anchors are used as foundations for mooring deep-water offshore facilities, including risers and floating structures. They are cone-tipped cylindrical steel pipes ballasted with concrete and scrap metal and penetrate the seabed by the kinetic energy they acquire during free fall through the water. A mooring line is usually connected at the top of the anchor. The design of such anchors involves estimation of the embedment depth as well as short-term and long-term pullout capacities. This paper describes the development of a computational procedure that leads to prediction of torpedo-anchor embedment depth. The procedure relies on a computational fluid dynamics (CFD) model for evaluation of the resisting forces on the anchor. In the model, the soil is represented as a viscous fluid and the procedure is applied to axially symmetric penetration of the seabed. The CFD approach provides estimates of not only the embedment depth but the pressure and shear distributions on the soil-anchor interface and in the soil.     

Three-Dimensional Numerical Simulation and Analysis of Flows around a Submerged Weir in a Channel Bendway

To improve navigation conditions for barges passing through river channels, many submerged weirs (SWs) have been installed along the bendways of many waterways by the U.S. Army Corps of Engineers. This paper presents results from three-dimensional numerical simulations that were conducted to study the helical secondary current (HSC) and the near-field flow distribution around one SW. The simulated flow fields around a SW in a scale physical model were validated using experimental data. The three-dimensional flow fields around a SW, the influence of the SW on general HSC, and the implication of effectiveness of submerged weirs to realign the flow field and improve navigability in bendways were analyzed. The numerical simulations indicated that the SW significantly altered the general HSC. Its presence induced a skewed pressure difference cross its top and a triangular-shaped recirculation to the downstream side. The over-top flow tends to realign toward the inner bank and therefore improves conditions for navigation.     

Theory of Three-Dimensional Discontinuous Deformation Analysis and Its Application to a Slope Toppling at Amatoribashi, Japan

A three-dimensional (3D) rock slope toppling occurred in a discontinuous rock mass. To simulate the failure process and study the mechanism of this rock failure with contact and large displacement in 3D, a new discrete numerical method has been developed called the 3D discontinuous deformation analysis (DDA). This article first introduces the basic principles and then derives the formulas in detail. Finally, the slope failure simulation is applied as an example to investigate the applicability of this new method to rock slope failure research. The simulation results indicate the advantages of using this new method to study the mechanism of a rock slope failure with 3D behavior.     

Earth Dam on Liquefiable Foundation and Remediation: Numerical Simulation of Centrifuge Experiments

A series of four dynamic centrifuge model tests was performed to investigate the effect of foundation densification on the seismic performance of a zoned earth dam with a saturated sand foundation. In these experiments, thickness of the densified foundation layer was systematically increased, resulting in a comprehensive set of dam-foundation response data. Herein, Class-A and Class-B numerical simulations of these experiments are conducted using a two-phase (solid and fluid) fully coupled finite element code. This code incorporates a plasticity-based soil stress–strain model with the modeling parameters partially calibrated based on earlier studies. The physical and numerical models both indicate reduced deformations and increased crest accelerations with the increase in densified layer thickness. Overall, the differences between the computed and recorded dam displacements are under 50%. At most locations, the computed excess pore pressure and acceleration match the recorded counterparts reasonably well. Based on this study, directions for further improvement of the numerical model are suggested.     

Coherent Structure Dynamics in Turbulent Flows Past In-Stream Structures: Some Insights Gained via Numerical Simulation

Large-scale coherent vortical structures in natural streams and rivers dominate flow and transport processes and impact the stability of stream banks, the diversity and abundance of organisms, and the quality of running waters in aquatic ecosystems. Thus, understanding and being able to model the dynamics of energetic coherent structures in such flows at ecologically relevant scales are crucial prerequisites for developing a science-based ecosystem restoration framework. We review recent progress toward the development of coherent-structure-resolving (CSR) computational fluid dynamics techniques, based on hybrid URANS/LES modeling strategies, for simulating turbulent flows in open-channels with hydraulic structures. CSR simulations of the turbulent horseshoe vortex (THSV) past bed-mounted piers explained the physical mechanism leading to the experimentally documented bimodal velocity fluctuations of the vortex and underscored the importance of the Reynolds number as a key parameter governing the THSV dynamics. Simulations of high Reynolds number flows past surface-piercing, groynelike structures in open channels revealed the complexity of the recirculating region at the upstream face of the groyne, underscored the interaction of the flow in this region with the energetic shear layer shed from the point of separation at the upstream side wall, and demonstrated the importance of flow depth in the vorticity dynamics of such flows. The paper also identifies areas for future work and modeling challenges that need to be addressed for the computational tools to be able to accurately predict flow and transport processes in real-life aquatic environments.     

Direct Numerical Simulation of Turbulent Flow in a Square Duct: Analysis of Secondary Flows

A direct numerical simulation of turbulent flow in a square duct was performed for a Reynolds number based on bulk streamwise velocity and duct height equal to 4,440. The mechanism by which secondary flows are generated in a square duct was investigated. Two counterrotating secondary flows occur around the duct corner. These secondary flows were found to play a key role in momentum transfer between the corner and center of the duct. A conditional quadrant analysis was performed in the local maximum and minimum regions of the wall shear stress in order to characterize the pattern of the mean secondary flows.     

Case Study: Refinement of Hydraulic Operation of a Complex CSO Storage/Treatment Facility by Numerical and Physical Modeling

The performance of a combined sewer overflow (CSO) storage/treatment facility in North Toronto, Ont., Canada, was investigated by conjunctive numerical and physical (hydraulic) modeling. The main objectives of the study were to (1) assess the feasibility of increasing the hydraulic loading of the CSO facility without bypassing; and (2) establish a verified numerical model of the facility for future work. The numerical model [a commercial computational fluid dynamics (CFD), PHOENICS] was validated and verified using results from a hydraulic scale model (1:11.6). The results obtained show that the CFD model can simulate hydraulic conditions in the facility well, as demonstrated by accurate reproduction of the filling rate, water levels at various locations, flow velocities in feed pipes, and overflows from the inflow channel. Numerical simulations identified excessive local head losses and helped select structural changes to reduce such losses. The analysis of the facility showed that with respect to hydraulic operation, the facility is a complex, highly nonlinear hydraulic system. Within the existing constraints, a few structural changes examined by numerical simulation could increase the maximum treatment flow rate in the CSO storage/treatment facility by up to 31%.     

Numerical Simulation of Obstructed Round Buoyant Jets in a Static Uniform Ambient

The k-ε turbulence closure model is used to simulate obstructed round buoyant jets in a static uniform ambient, and the results compare well with available experimental data. On the basis of the axial line velocity distribution, three regions in the flow behind the disk are identified: the wake region, the transitional region, and the self-similarity region. The length of the wake region, which varies with flow and geometrical parameters, and the existence of self-similarity are also addressed.     

Numerical Analysis of Peristaltic Transport of a Tangent Hyperbolic Fluid in an Endoscope

This paper presents a tangent hyperbolic fluid in a cylindrical coordinate system. The governing equations are simplified using long-wavelength and low Reynolds number approximations. The solutions of the problem in simplified form are calculated with three methods: (1)?perturbation method, (2)?homotopy analysis method, and (3)?shooting method. Comparison of the three solutions shows very good agreement among them. Pressure rise and frictional force are calculated with the help of numerical integration. Graphical results for pressure rise and frictional forces are presented to show the physical behavior of the Weissenberg number We, amplitude ratio ?, tangent hyperbolic power law index n, and radius ratio ?.     

Lateral Spreading Forces on Bridge Piers and Pile Caps in Laterally Spreading Soil: Effect of Angle of Incidence

In this paper, the kinematic forces which may be applied to bridge piers or pile caps from laterally spreading surficial cohesive soil layers (nonliquefied crusts) through which they pass are considered. Such forces often represent the largest load component acting on a structure and/or foundation during liquefaction-induced lateral spreading. Both circular and square structural inclusions are considered, and particular attention is paid to the orientation of the inclusion to the direction of spreading, here defined as the angle of incidence (θ). Experimental modeling was conducted using a modified direct shearbox to simulate the spreading of kaolin past structural inclusions at various θ. Load-displacement data and particle image velocimetry analysis revealed that the ultimate load for both square and circular cases may be determined using a wedge-based upper-bound plasticity analysis. For circular sections, this ultimate load is independent of θ due to radial symmetry. The ultimate load on square sections was found to depend more significantly on θ and a simple analytical method is presented to account for this. The method suggests that the ultimate loads acting on square bridge piers or pile caps will be a maximum when the spreading soil impinges on the corners of the inclusion, at which time the ultimate load will be 19–26% larger (depending on the soil-structure interface roughness) than for spreading impinging on the edge of the inclusion. Experimental tests suggested a value of 22%. Finally, the tests support previous results suggesting that when the underlying soil is unable to carry redistributed shear stress (i.e., when it is liquefied) load-displacement curves in the crustal layers are less stiff than for typical retaining structures under static conditions. The displacement at soil yield was found to be between 20–30% of the height of the inclusion in the layer, and also depends on θ in the case of square inclusions.     

Rapid Graphical Detection of Weakness Problems in Numerical Simulation Infiltration Models Using a Linearized Form Equation

Recently, Valiantzas proposed a new two-parameter vertical infiltration equation that can be transformed to a linearized-form equation that essentially states that the shape of the cumulative infiltration data, when presented in the form of (i2/t) versus i, is linear. In this paper, the presentation of the numerical data to the Valiantzas linearized-form equation is proposed as an additional criterion to detect easily and rapidly possible errors of the numerical solutions and eventually to choose the best spatial discretization for a simulated infiltration event that is used as setup parameter to the numerical infiltration models. Numerical data and analytical solutions were used to validate the proposed method.     

Analysis of Traffic-Induced Vibrations in a Cable-Stayed Bridge. Part II: Numerical Modeling and Stochastic Simulation

This paper describes the development of a numerical model to simulate the dynamic response of the bridge–vehicle system of Salgueiro Maia cable-stayed bridge, using the results from an extensive experimental investigation to calibrate this model. Further, a set of stochastic Monte Carlo simulations of the bridge–vehicle dynamic response is also presented, with the purpose of evaluating dynamic amplification factors, taking into account the randomness of different factors associated to characteristics of the pavement, of the vehicles and of the traffic flow.