Conduit Model

Closed

The conduit model is used to calculate the transport of suspended sediment and dissolved pollutant, and the erosion and deposition of sediment, in conduits. The transport process and the sediment erosion and deposition process are solved separately within each time step.

As with the hydraulic conduit model, a conduit is represented as a conceptual link of defined length between two nodes in the network. Control structures are treated as links of zero length in which no erosion or deposition takes place.

It is assumed that:

Transport

The equation describing the transport of the suspended sediment and the dissolved pollutant is based on conservation of mass. With the assumptions listed above, this leads to the one-dimensional advection equation as described in, for example, Cunge J A et al (1980).

where

c is the concentration (kg/m3)

u is the flow velocity (m/s) - obtained from the hydraulic simulation

t is time (s)

x is the spatial co-ordinate (m).

The boundary condition at the upstream end required for this equation is generated by the Network Model.

The carrying capacity of the flow is calculated using one of the three erosion/deposition models available in InfoWorks ICM . See Sediment Erosion and Deposition for more details.

Numerical Techniques

Advection

The advection equation is solved in each conduit by the Holly-Preissmann scheme (Holly F.M. & Preissmann A. (1977)). This is a semi-Lagrangian method. It tracks conceptual parcels of pollutant in the flow.

The solution of the advection equation is known to be constant:

along trajectories, X, given by:

This differential equation is solved for the trajectories by the mid-point rule:

with the velocities at the half time level generated by averaging the velocities at the two most recent time levels, and the position of the trajectory at the half time level generated by iteration. This gives a method which is second order in time.

Since the position of the foot of a trajectory at the previous known time level, Xn, does not generally coincide with a computational point, the value at the foot is generated by interpolation. The solution generated by this scheme is unconditionally stable - there is no restriction on the allowable size of the time step.

Cubic interpolation at the foot of a trajectory gives a solution which is third order in space but it can introduce un-physical oscillations to the numerical solution. This is avoided by the use of flux-corrected transport (FCT) (Boris J.P. and Book D.L. (1973)) - the local weighted averaging of two different numerical solutions to impose monotonicity.

The flux-corrected transport solution is not conservative. In the InfoWorks ICMwater quality model, conservation is recovered by the use of sub-optimal weights in the FCT local averaging. See Priestley A. (1993).

1D Diffusion

It is also possible to simulate 1d diffusion, if required. Once the 1D engine has solved advection of the determinants, it can then calculate the diffusion. The equations for 1D diffusion are described in the topic 1D Diffusion.

Sediment Erosion and Deposition

The erosion/deposition model to be used is set in the Water Quality and Sediment Parameters for the network.

The following assumptions and limitations apply to erosion and deposition of sediment:

InfoWorks ICM supports three different models for calculating erosion and deposition in pipes. These models are:

Dependent Sediment Fractions

Multiple sediment fractions can be modelled independently or dependently. If multiple sediment fractions are modelled independently then the above algorithm is used for each sediment fraction. If the sediment fractions are modelled dependently then the following algorithm is used at the end of each water quality timestep.

  1. Calculate the total load in the flow

  2. Calculate single representative d50 and s for the fractions in the flow as a weighted average based on the concentration of the fractions present in the flow

  3. Calculate the flow capacity for the representative fraction using the carrying capacity equation above.
  4. If cT > rsc*n (flow capacity is exceeded) then deposition occurs, and the excess is spread between the fractions in the flow

  5. If cT < rsc*n (there is available flow capacity) then erosion occurs, and the available capacity is spread between the fractions on the bed

Water Quality Model

Network Model

1D Diffusion