skip to main content


Title: Three-phase flow simulation of local scour around a submerged horizontal cylinder
Wave-induced scour plays a key role in the stability analysis of coastal structures, submarine pipelines or cables. There is a rich literature in current-induced scour, but more research is needed to understand the characteristics of wave-induced scour and the mechanisms that are important to the scour process. Sediment transport and flow-induced scour are three-phase (air-water-sediment) flow problems in nature and multi-phase flow simulation is a useful tools that can provide information difficult to obtain from physical tests. Most existing numerical models developed for simulating local scours are based on one-way coupling, which neglects effects of sediment phase on hydrodynamics of the flow. The present study uses a three-phase (air, water and sediment) flow model, which allows for a two-way coupling, to simulate wave-induced local scour problems. The three-phase flow model captures the air-water interface using a modified VOF method, and uses an improved rheology for the sediment phase for better results. The model is validated and verified using one set of existing experiment results for local scour around a submerged horizontal pipe. The detailed flow fields of both the sediment phase and the water phase around the scour are analyzed to understand the scour process. All three-phase flow simulations flow simulations on XSEDE’s Stampede2 supercomputers. The applicability of the model to other local scour problems is also discussed.  more » « less
Award ID(s):
1706938
NSF-PAR ID:
10299436
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Proceedings of the Thirtieth (2020) International Ocean and Polar Engineering Conference
Page Range / eLocation ID:
1612-1617
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. A three-dimensional Eulerian two-phase flow solver, SedFoam, has been developed for various sediment transport applications. The solver has demonstrated success in modeling sheet flow and bedforms driven by oscillatory flows using a Reynolds-averaged Navier–Stokes (RANS) formulation. However, the accuracy of the RANS formulation for more complex flows, such as scour around structures, requires further evaluation. SedFoam has recently been enhanced to incorporate two-phase large-eddy simulation (LES) capability. In this study, RANS and LES approaches are tested via a three-dimensional case of wave-induced local scour around a single vertical circular pile. Two laboratory experiments, one with an erodible bed and the other with a rigid bed, were chosen for simulation, with both experiments having a Keulegan-Carpenter (KC) number of 10. The k-ω turbulence closure was selected for the RANS simulation, and the dynamic Lagrangian subgrid closure was chosen for the LES simulation. Numerical results reveal that both RANS and LES simulations can resolve lee-wake vortices, although the vortices are significantly weaker in the RANS simulation. In comparison with the LES results, the RANS approach fails to predict horseshoe vortex with sufficient intensity, leading to an underestimation of scour hole depth development. Although the scour depths develop at a very similar rate in the early stage, the scour depth predicted by the RANS simulation quickly reaches equilibrium, while the LES simulation follows the measured trend. These findings indicate that a turbulence-resolving methodology, i.e. LES, is necessary for accurate scour simulations. 
    more » « less
  2. Abstract

    For most of the year, a dry‐bed desert wash is void of water flow. Intensive rain events, however, could trigger significant flash floods that bring about highly complicated hydrodynamics and morphodynamics processes within a desert stream. We present a fully coupled three‐phase flow model of air, water, and sediment to simulate numerically the propagation of a flash flood in a field‐scale fluvial desert stream, the so‐called Tex Wash located in the Mojave Desert, California, United States. The turbulent flow of the flash flood is computed using the three‐dimensional unsteady Reynolds‐averaged Navier–Stokes equations closed with the shear stress transportkωmodel. The free surface of the flash flood at the interface of air and water phases is computed with the level‐set method, which enables instantaneous tracking of the water surface as the flash flood propagates over the dry bed of the desert stream. The evolution of the desert fluvial stream's morphology, due to the action of the propagating flash flood on the mobile bed, is calculated using a Eulerian morphodynamics model based on the curvilinear immersed boundary method. The capabilities of the proposed numerical framework are demonstrated by applying it to simulate a flash flood event in a 0.65‐km‐long reach of the Tex Wash, the intricate channel morphology of which is obtained using light imaging detection and ranging technology. The simulated region of the stream includes a number of bridge foundations. The simulation results of the model for the flash flood event revealed the formation of a highly complex flow field and scour patterns within the stream. Moreover, our simulation results showed that most scour processes take place during the steady phase of the flash flood, that is, after the flash flood fills the stream. The transient phase of the flash flood is rather short and contributes to a very limited amount of erosion within the desert stream.

     
    more » « less
  3. null (Ed.)
    Porous hydraulic structures, such as Large Woody Debris (LWD) and Engineered Log Jams (ELJs), play a very important role in erosion control and habitation conservation in rivers. Previous experimental research has shed some light on the flow and sediment dynamics through and around porous structures. It was found that the scour process is strongly dependent on porosity. Computational models have great value in revealing more details of the processes which are difficult to capture in laboratory experiments. For example, previous computational modeling work has shown that the level of resolution of the complex hydraulic structures in computer models has great effect on the simulated flow dynamics. The less computationally expensive porosity model, instead of resolving all geometric details, can capture the bulk behavior for the flow field, especially in the far field. In the near field where sediment transport is most intensive, the flow result is inaccurate. The way in which this error is translated to the sediment transport results is unknown. This work aims to answer this question. More specifically, the suitability and limitations of using a porosity model in simulating scour around porous hydraulic structures are investigated. To capture the evolution of the sediment bed, an immersed boundary method is implemented. The computational results are compared against flume experiments to evaluate the performance of the porosity model. 
    more » « less
  4. null (Ed.)
    Abstract

    Shock waves from underwater and air explosions are significant threats to surface and underwater vehicles and structures. Recent studies on the mechanical and thermal properties of various phase-separated elastomers indicate the possibility of applying these materials as a coating to mitigate shock-induced structural failures. To demonstrate this approach and investigate its efficacy, this paper presents a fluid-structure coupled computational model capable of predicting the dynamic response of air-backed bilayer (i.e. elastomer coating – metal substrate) structures submerged in water to hydrostatic and underwater explosion loads. The model couples a three-dimensional multiphase finite volume computational fluid dynamics model with a nonlinear finite element computational solid dynamics model using the FIVER (FInite Volume method with Exact multi-material Riemann solvers) method. The kinematic boundary condition at the fluid-structure interface is enforced using an embedded boundary method that is capable of handling large structural deformation and topological changes. The dynamic interface condition is enforced by formulating and solving local, one-dimensional fluid-solid Riemann problems, which is well-suited for transferring shock and impulsive loads. The capability of this computational model is demonstrated through a numerical investigation of hydrostatic and shock-induced collapse of aluminum tubes with polyurea coating on its inner surface. The thickness of the structure is resolved explicitly by the finite element mesh. The nonlinear material behavior of polyurea is accounted for using a hyper-viscoelastic constitutive model featuring a modified Mooney-Rivlin equation and a stress relaxation function in the form of prony series. Three numerical experiments are conducted to simulate and compare the collapse of the structure in different loading conditions, including a constant pressure, a fluid environment initially in hydrostatic equilibrium, and a two-phase fluid flow created by a near-field underwater explosion.

     
    more » « less
  5. The response of the Greenland Ice Sheet (GrIS) to a warmer climate is uncertain on long time scales. Climate models, such as those participating in the Coupled Model Intercomparison Project phase 6 (CMIP6), are used to assess this uncertainty. The Community Earth System Model version 2.1 (CESM2) is a CMIP6 model capable of running climate simulations with either one‐way coupling (fixed ice sheet geometry) or two‐way coupling (dynamic geometry) to the GrIS. The model features prognostic snow albedo, online downscaling using elevation classes, and a firn pack to refreeze percolating melt water. Here we evaluate the representation of the GrIS surface energy balance and surface mass balance in CESM2 at 1° resolution with fixed GrIS geometry. CESM2 agrees closely with ERA‐Interim reanalysis data for key controls on GrIS SMB: surface pressure, sea ice extent, 500 hPa geopotential height, wind speed, and 700 hPa air temperature. Cloudsat‐CALIPSO data show that supercooled liquid‐containing clouds are adequately represented, whereas comparisons to Moderate Resolution Imaging Spectroradiometer and CM SAF Cloud, Albedo, and Surface Radiation data set from Advanced Very High Resolution Radiometer data second edition data suggest that CESM2 underestimates surface albedo. The seasonal cycle and spatial patterns of surface energy balance and surface mass balance components in CESM2 agree well with regional climate model RACMO2.3p2, with GrIS‐integrated melt, refreezing, and runoff bracketed by RACMO2 counterparts at 11 and 1 km. Time series of melt, runoff, and SMB show a break point around 1990, similar to RACMO2. These results suggest that GrIS SMB is realistic in CESM2, which adds confidence to coupled ice sheet‐climate experiments that aim to assess the GrIS contribution to future sea level rise.

     
    more » « less