skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: SedFoam-2.0: a 3-D two-phase flow numerical model for sediment transport

Abstract. In this paper, a three-dimensional two-phase flow solver, SedFoam-2.0, is presented for sediment transport applications. The solver is extended from twoPhaseEulerFoam available in the 2.1.0 release of the open-source CFD (computational fluid dynamics) toolbox OpenFOAM. In this approach the sediment phase is modeled as a continuum, and constitutive laws have to be prescribed for the sediment stresses. In the proposed solver, two different intergranular stress models are implemented: the kinetic theory of granular flows and the dense granular flow rheology μ(I). For the fluid stress, laminar or turbulent flow regimes can be simulated and three different turbulence models are available for sediment transport: a simple mixing length model (one-dimensional configuration only), a k − ε, and a k − ω model. The numerical implementation is demonstrated on four test cases: sedimentation of suspended particles, laminar bed load, sheet flow, and scour at an apron. These test cases illustrate the capabilities of SedFoam-2.0 to deal with complex turbulent sediment transport problems with different combinations of intergranular stress and turbulence models.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Geoscientific Model Development
Page Range / eLocation ID:
4367 to 4392
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

    Entrainment and suspension of sediment particles with the size distribution similar to a range of natural sands were simulated with a focus on the vertical size sorting and transport dynamics in response to different wave conditions. The simulations were performed using a two‐phase Eulerian‐Lagrangian model by combining the LIGGGHTS discrete element method solver for sediment and SedFoam solver for the fluid phase. The model was first validated for a range of sand grain sizes from 0.21 to 0.97 mm having well‐sorted and mixed (bimodal) size distributions using laboratory oscillatory flow data. Three sediment bed configurations were studied under a wide range of velocity‐skewed waves with different wave intensity and skewness. It was found that the bimodal distribution having only 30% of coarse fraction and 70% of medium fraction responds similar to a well‐sorted coarse sand configuration. Sediment fluxes of the bimodal distribution were slightly higher than those of well‐sorted coarse sand because of the pronounced inverse grading in the bimodal distribution. Furthermore, for the bimodal distribution the medium fraction acted as a relatively smooth foundation underneath the coarse fraction which facilitated the mobilization of the coarser particles. Under high energy wave conditions, the smoothing feature was exacerbated and further caused the formation of plug flow where a thick layer of intense sediment flux was observed. Model results also showed that under high skewness waves, phase‐lag effect occurred in well‐sorted medium sand which caused lower net onshore sediment transport rates but the effect was significantly reduced for mixed sediments.

    more » « less
  3. One of the key factors in simulating realistic wall-bounded flows at high Reynolds numbers is the selection of an appropriate turbulence model for the steady Reynolds Averaged Navier–Stokes equations (RANS) equations. In this investigation, the performance of several turbulence models was explored for the simulation of steady, compressible, turbulent flow on complex geometries (concave and convex surface curvatures) and unstructured grids. The turbulence models considered were the Spalart–Allmaras model, the Wilcox k- ω model and the Menter shear stress transport (SST) model. The FLITE3D flow solver was employed, which utilizes a stabilized finite volume method with discontinuity capturing. A numerical benchmarking of the different models was performed for classical Computational Fluid Dynamic (CFD) cases, such as supersonic flow over an isothermal flat plate, transonic flow over the RAE2822 airfoil, the ONERA M6 wing and a generic F15 aircraft configuration. Validation was performed by means of available experimental data from the literature as well as high spatial/temporal resolution Direct Numerical Simulation (DNS). For attached or mildly separated flows, the performance of all turbulence models was consistent. However, the contrary was observed in separated flows with recirculation zones. Particularly, the Menter SST model showed the best compromise between accurately describing the physics of the flow and numerical stability. 
    more » « less
  4. null (Ed.)
    Abstract Aerosol jet printing (AJP) is a direct-write additive manufacturing technique, which has emerged as a high-resolution method for the fabrication of a broad spectrum of electronic devices. Despite the advantages and critical applications of AJP in the printed-electronics industry, AJP process is intrinsically unstable, complex, and prone to unexpected gradual drifts, which adversely affect the morphology and consequently the functional performance of a printed electronic device. Therefore, in situ process monitoring and control in AJP is an inevitable need. In this respect, in addition to experimental characterization of the AJP process, physical models would be required to explain the underlying aerodynamic phenomena in AJP. The goal of this research work is to establish a physics-based computational platform for prediction of aerosol flow regimes and ultimately, physics-driven control of the AJP process. In pursuit of this goal, the objective is to forward a three-dimensional (3D) compressible, turbulent, multiphase computational fluid dynamics (CFD) model to investigate the aerodynamics behind: (i) aerosol generation, (ii) aerosol transport, and (iii) aerosol deposition on a moving free surface in the AJP process. The complex geometries of the deposition head as well as the pneumatic atomizer were modeled in the ansys-fluent environment, based on patented designs in addition to accurate measurements, obtained from 3D X-ray micro-computed tomography (μ-CT) imaging. The entire volume of the constructed geometries was subsequently meshed using a mixture of smooth and soft quadrilateral elements, with consideration of layers of inflation to obtain an accurate solution near the walls. A combined approach, based on the density-based and pressure-based Navier–Stokes formation, was adopted to obtain steady-state solutions and to bring the conservation imbalances below a specified linearization tolerance (i.e., 10−6). Turbulence was modeled using the realizable k-ε viscous model with scalable wall functions. A coupled two-phase flow model was, in addition, set up to track a large number of injected particles. The boundary conditions of the CFD model were defined based on experimental sensor data, recorded from the AJP control system. The accuracy of the model was validated using a factorial experiment, composed of AJ-deposition of a silver nanoparticle ink on a polyimide substrate. The outcomes of this study pave the way for the implementation of physics-driven in situ monitoring and control of AJP. 
    more » « less
  5. The effect of slip surfaces on the laminar–turbulent separatrix of plane Poiseuille flow is studied by direct numerical simulation. In laminar flows, the inclusion of the slip surfaces results in a drag reduction of over 10 %, which is in good agreement with previous studies and the theory of laminar slip flows. Turbulence lifetimes, the likelihood that turbulence is sustained, is investigated for transitional flows with various slip lengths. We show that slip surfaces decrease the likelihood of sustained turbulence compared to the no-slip case, and the likelihood is further decreased as slip length is increased. A more deterministic analysis of the effects of slip surfaces on a transition to turbulence is performed by using nonlinear travelling-wave solutions to the Navier–Stokes equations, also known as exact coherent solutions. Two solution families, dubbed P3 and P4, are used since their lower-branch solutions are embedded on the boundary of the basin of attraction of laminar and turbulent flows (Park & Graham, J. Fluid Mech. , vol. 782, 2015, pp. 430–454). Additionally, they exhibit distinct flow structures – the P3 and P4 are denoted as core mode and critical-layer mode, respectively. Distinct effects of slip surfaces on the solutions are observed by the skin-friction evolution, linear growth rate and phase-space projection of transitional trajectories. The slip surface appears to modify the transition dynamics very little for the core mode, but quite considerably for the critical-layer mode. Most importantly, the slip surface promotes different transition dynamics – an early and bypass-like transition for the core mode and a delayed and H- or K-type-like transition for the critical-layer mode. We explain these distinct transition dynamics based on spatio-temporal and quadrant analyses. It is found that slip surfaces promote the prevalence of strong wall-toward motions (sweep-like events) near vortex cores close to the channel centre, inducing an early transition, while long sustained ejection events are present in the region of the $\unicode[STIX]{x1D6EC}$ -shaped vortex cores close to the critical layer, resulting in a delayed transition. This should motivate flow control strategies to fully exploit these distinct transition dynamics for transition to turbulence. 
    more » « less