Diffusion
    • 22 Sep 2022
    • 2 Minutes to read

    Diffusion


    Article summary

    The unit models the longitudinal dispersion of any transported water quality variable in 1D Water Quality Solver.

    Description

    The transport and dispersion of water quality variables is governed by the mass balance equation. One term of the equation represents the effective longitudinal dispersion. The amount of dispersion is governed by the longitudinal dispersion coefficient D. 1D Water Quality Solver offers three ways of calculating the longitudinal dispersion coefficient. Different methods of calculating the longitudinal dispersion coefficient may be applied to different parts of the river system.

    The longitudinal dispersion may be set to a constant value, D0, irrespective of discharge or shear velocity, which is applied at all nodes throughout the model. This is set by using the keyword CONSTANT. Alternatively the constant value may be allowed to vary within the model area, in which case the keyword VARIABLE is used.

    For estuary networks the longitudinal dispersion coefficient may be considered to vary with the tidal discharge and/or the longitudinal salinity gradient (mainly applicable to large delta systems). In these cases the keyword ESTUARY is used. The longitudinal dispersion coefficient is represented by a three part empirical function (see equation below). The size of the coefficient is governed by three empirical constants D0, D1 and D2. D0 is a constant term which is used to prevent unrealistically low values of dispersion from arising. D1 is a term which is used to represent the amount of dispersion caused by tidal mixing in the estuary. D2 is the term used to determine the effect of salinity gradient on dispersion.

    For river applications the longitudinal dispersion is related to the shear velocity (see below). In these cases the keyword RIVER is used.

    It is possible that user may wish to have varying values of ESTUARY or RIVER dispersion coefficients throughout the model. The user may even wish to mix the dispersion type. In these cases the keyword DYNAMIC should be used. This keyword enables the user to set individual values of D0, D1 and D2 and diffusion type at each model node.

    Equations

    For estuary networks the dispersion coefficient has the form:

    image                    (1)

    Where:
    D0, D1 and D2= Empirical constants (m2/s)
    Q = Discharge (m3/s)
    U0= Peak velocity at the seaward boundary (m/s)
    A0= Representative cross-sectional area at the seaward boundary (m2)
    S0= Representative salinity at the seaward boundary (ppt)
    L0= Representative limit of saline intrusion from the seaward boundary (m2).

    For river applications the dispersion coefficient has the form:

    image                    

    image                     (2)

    Where:
    u* = Shear velocity (m/s)
    f = Colebrook-White friction factor
    u = Mean flow velocity (m/s)
    b = Channel width/breadth (m)
    D1 = Empirical diffusion coefficient


    Was this article helpful?

    Changing your password will log you out immediately. Use the new password to log back in.
    First name must have atleast 2 characters. Numbers and special characters are not allowed.
    Last name must have atleast 1 characters. Numbers and special characters are not allowed.
    Enter a valid email
    Enter a valid password
    Your profile has been successfully updated.
    ESC

    Eddy AI, facilitating knowledge discovery through conversational intelligence