Also known as SedFoam
Model type Modular
Model part of larger framework
Incorporated models or components:
Spatial dimensions 3D
Spatial extent Patch-Scale, Grain-Scale
Model domain Hydrology, Terrestrial
One-line model description A four-way coupled two-phase Eulerian model for sediment transport
Extended model description A multi-dimensional numerical model for sediment transport based on the two-phase

flow formulation is developed. With closures of particle stresses and fluid-particle interaction, the model is able to resolve processes in the concentrated region of sediment transport and hence does not require conventional bedload/suspended load assumptions. The numerical model is developed in three spatial dimensions. However, in this version, the model is only validated for Reynolds-averaged two-dimensional vertical (2DV) formulation (with the k − epsilon closure for carrier flow turbulence) for sheet flow in steady and oscillatory flows. This numerical model is developed via the open-source CFD library of solvers, OpenFOAM and the new solver is called twoPhaseEulerSedFoam.


sediment transport, two-phase model, multi-dimensional,

First name Zhen
Last name Cheng
Type of contact Model developer
Institute / Organization Ocean Engineering Lab, University of Delaware
Postal address 1 259 Academy Street
Postal address 2
Town / City Newark
Postal code 19716
State Delaware
Country United States
Email address
Phone +1 (302)4196006

First name Tian-Jian
Last name Hsu
Type of contact Model developer
Institute / Organization Ocean Engineering Lab, University of Delaware
Postal address 1 259 Academy St.
Postal address 2
Town / City Newark
Postal code 19716
State Delaware
Country United States
Email address
Phone +1 (302) 831-4172

Supported platforms Unix, Linux, Mac OS, Windows
Other platform
Programming language C++
Other program language
Code optimized Single Processor, Multiple Processors
Multiple processors implemented Distributed memory
Nr of distributed processors
Nr of shared processors
Start year development 2011
Does model development still take place? Yes
If above answer is no, provide end year model development
Code development status
When did you indicate the 'code development status'?
Model availability As code
Source code availability
(Or provide future intension)
Through CSDMS repository
Source web address
Source csdms web address
Program license type CDDL
Program license type other
Memory requirements
Typical run time one day for typical 50,000 grid points using 1 processor, model scalability has been tested to be very good.

Describe input parameters flow forcing; sediment properties ( density, grain size, etc.); fluid properties; coefficients for carrier fluid turbulence, and parameters for kinetic theory for granular flow; model selection for kinetic theory, such as granular pressure, conductivity, and viscosity model, etc. More details are described in the user maunal.
Input format ASCII
Other input format
Describe output parameters sediment concentration (alpha); carrier fluid turbulence (k, epsilon); granular temperature (Theta); fluid pressure (p); sediment and fluid velocities (Ua, Ub); turbulence modulation factor (tmf); particle pressure ( kinetic part, pa), more details are described in the user maunal.
Output format ASCII
Other output format
Pre-processing software needed? No
Describe pre-processing software
Post-processing software needed? No
Describe post-processing software
Visualization software needed? No
If above answer is yes
Other visualization software

Describe processes represented by the model Sediment transport under steady channel flow, oscillatory flow (sinusoidal and Stokes 2nd order waves)
Describe key physical parameters and equations please see more details in the user manual
Describe length scale and resolution constraints domain needs to be large enough to capture the boundary layer, and the resolution ( on the order of grain size) should be fine enough to capture the sediment-fluid interaction near the bed level.
Describe time scale and resolution constraints The flow should be run to be fully developed, and this time scale depends on the time scale of sediment settling time, and flow periods. the time step should be fine enough not to cause numerical instability, and also capture the varying the flow forcing.
Describe any numerical limitations and issues This model is only tested in rectangle domains, and compared the results with idealized experiments;

Sediment bed state is affected by the initial condition (mainly due to the frictional stess closure in this model). For 2D or 3D runs, it is suggested that first run 1DV to steady or quasi-steady state, and map the 1DV results to 2DV or 3D, in this way, the initial instability of the sediment bed can be avoided.

Describe available calibration data sets Experiment data for steady channel flow can be found in:

Sumer, B. M., Kozakiewicz, A., Fredsoe, J., Deigaard, R., 1996. Velocity and concentration profiles in sheet-flow layer of movable bed. Journal of Hydraulic Engineering, (1996) 549-558.

Experiment for oscillatory flow can be found in: O'Donoghue, T., Wright, S., 2004. Concentrations in oscillatory sheet flow for well sorted and graded sands. Coastal Engineering 50 (2004) 117-138.

Upload calibration data sets if available:
Describe available test data sets
Upload test data sets if available:
Describe ideal data for testing

Do you have current or future plans for collaborating with other researchers?
Is there a manual available? Yes
Upload manual if available: Media:User Manual.pdf
Model website if any
Model forum / discussion board

This part will be filled out by CSDMS staff

OpenMI compliant No but possible
BMI compliant No but possible
WMT component No but possible
PyMT component
Is this a data component
Can be coupled with:
Model info
Nr. of publications: --
Total citations: 0
h-index: --"--" is not a number.
Qrcode TwoPhaseEulerSedFoam.png
Link to this page




Nr. of publications: --
Total citations: 0
h-index: --"--" is not a number.

See more publications of TwoPhaseEulerSedFoam



Input Files

Output Files