Skip to main content
Engineering LibreTexts

7.10: The Concentration Distribution

  • Page ID
    31789
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \) \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)\(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\) \(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\)\(\newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    7.10.1 The Advection Diffusion Equation

    In chapter 5 sediment transport in open channel flow and pipe flow has been discussed. The advection diffusion equation is derived and solved for low concentrations and for high concentrations including the effect of upwards liquid flow. Wasp et al. (1977) and Doron et al. (1987) use the solution for low concentrations, while Karabelas (1977) and Kaushal & Tomita (2002B) use the Hunt (1954) approach with upwards liquid flow. Hindered settling is not yet included in the basic solutions, but added by replacing the terminal settling velocity by the hindered terminal settling velocity. For the diffusivity and the relation between the sediment diffusivity and the turbulent eddy momentum diffusivity different approaches are possible.

    If we assume the diffusivity is a constant, the advection diffusion equation can be solved. Giving the differential equation in the equilibrium situation with hindered settling, but without the upwards liquid velocity, since this is assumed to be part of the hindered settling (with r the vertical distance from the bottom of the pipe):

    \[\ \mathrm{C_{v s}(r) \cdot v_{t h}+\varepsilon_{s} \cdot \frac{d C_{v s}(r)}{d r}=C_{v s}(r) \cdot v_{t h}+\beta_{s m} \cdot \varepsilon_{m} \cdot \frac{d C_{v s}(r)}{d r}=0}\]

    The coordinate r ranges from 0 to Dp, the pipe diameter. Now the variables have to be separated according to:

    \[\ \frac{\mathrm{d} \mathrm{C}_{\mathrm{v} \mathrm{s}} \mathrm{r} \mathrm{)}}{\mathrm{C}_{\mathrm{v s}}(\mathrm{r})} \cdot=-\frac{\mathrm{v}_{\mathrm{t h}}}{\beta_{\mathrm{s m}} \cdot \varepsilon_{\mathrm{m}}} \cdot \mathrm{d r} \quad \Rightarrow \quad \ln \left(\mathrm{C}_{\mathrm{v s}}(\mathrm{r})\right)=-\frac{\mathrm{v}_{\mathrm{t h}}}{\beta_{\mathrm{s m}} \cdot \varepsilon_{\mathrm{m}}} \cdot \mathrm{r}+\mathrm{C}\]

    With Cvs(0)=CvB, the bottom concentration, the integration constant can be determined, giving:

    \[\ \mathrm{C}_{\mathrm{v s}}(\mathrm{r})=\mathrm{C}_{\mathrm{v B}} \cdot \mathrm{e}^{-\frac{\mathrm{v}_{\mathrm{t h}}}{\beta_{\mathrm{sm}} \cdot \varepsilon_{\mathrm{m}}} \cdot \mathrm{r}}\]

    This basic solution is still equal to the solution for open channel flow. Although this is just an indicative equation for open channel flow, Doron et al. (1987) and Doron & Barnea (1993) used it in their 2 and 3 layer models. The difference between pipe flow and open channel flow is in the determination of the diffusivity. Assuming the Law of the Wall, one can also determine the average diffusivity by integration (Lane & Kalinske (1941)):

    \[\ \varepsilon_{\mathrm{s}}=\beta_{\mathrm{sm}} \cdot \kappa \cdot \mathrm{u_{*}} \cdot \mathrm{r} \cdot\left(\frac{\mathrm{R}-\mathrm{r}}{\mathrm{R}}\right)=\beta_{\mathrm{sm}} \cdot \kappa \cdot \mathrm{u}_{*} \cdot \mathrm{R \cdot \frac{r}{R} }\cdot\left(1-\frac{\mathrm{r}}{\mathrm{R}}\right)=\beta_{\mathrm{sm}} \cdot \kappa \cdot \mathrm{u}_{*} \cdot \mathrm{R} \cdot \tilde{\mathrm{r}} \cdot(1-\tilde{\mathrm{r}})\]

    Integration over the cross section of the pipe gives:

    \[\ \begin{array}{left}\bar{\varepsilon}_{\mathrm{s}}&=\frac{1}{\pi \cdot \mathrm{R}^{2}} \cdot \int_{0}^{2 \cdot \pi} \int_{\mathrm{r=0}}^{\mathrm{r=R}} \varepsilon_{\mathrm{s}} \cdot \mathrm{dr} \cdot \mathrm{r} \cdot \mathrm{d} \phi=\beta_{\mathrm{sm}} \cdot \kappa \cdot \mathrm{u}_{*} \cdot \frac{\mathrm{R}^{3}}{\pi \cdot \mathrm{R}^{2}} \cdot \int_{\mathrm{0}}^{2 \cdot \pi}{ \int_{\tilde{\mathrm{r}}=0}^{\tilde{\mathrm{r}}=1}} \tilde{\mathrm{r}}^{2} \cdot(1-\tilde{\mathrm{r}}) \cdot \mathrm{d} \tilde{\mathrm{r}} \cdot \mathrm{d} \phi\\
    &=\beta_{\mathrm{sm}} \cdot \kappa \cdot \mathrm{u}_{*} \cdot \mathrm{R} \cdot \frac{1}{\pi} \cdot \int_{\mathrm{0}}^{2 \cdot \pi }\int_{\tilde{\mathrm{r}}=0}^{ \tilde{\mathrm{r}}=1} {\mathrm{\tilde { r }}^{2}} \cdot(1-\tilde{\mathrm{r}}) \cdot \mathrm{d} \tilde{\mathrm{r}} \cdot \mathrm{d} \phi=\beta_{\mathrm{s m}} \cdot \mathrm{\kappa} \cdot \mathrm{u}_{*} \cdot \mathrm{R} \cdot \frac{\mathrm{1}}{\pi} \cdot \int_{\mathrm{0}}^{2 \cdot \pi} \int_{\tilde{\mathrm{r}}=\mathrm{0}}^{\tilde{\mathrm{r}}=\mathrm{1}}\left(\tilde{\mathrm{r}}^{2}-\tilde{\mathrm{r}}^{3}\right) \cdot \mathrm{d} \tilde{\mathrm{r}} \cdot \mathrm{d} \phi\\
    &=\beta_{\mathrm{sm}} \cdot \kappa \cdot \mathrm{u}_{*} \cdot \mathrm{R} \cdot \frac{2 \cdot \pi}{\pi} \cdot\left(\frac{1}{3} \cdot \tilde{\mathrm{r}}^{3}-\frac{1}{4} \cdot \tilde{\mathrm{r}}^{4}\right)_{0}^{1}=\frac{\beta_{\mathrm{sm}} \cdot \mathrm{\kappa} \cdot \mathrm{u}_{*} \cdot \mathrm{R}}{\mathrm{6}}=\frac{\beta_{\mathrm{sm}} \cdot \mathrm{\kappa} \cdot \mathrm{u}_{*} \cdot \mathrm{D}_{\mathrm{p}}}{\mathrm{1} \mathrm{2}}\end{array}\]

    This gives for the concentration as a function of the vertical distance from the pipe bottom:

    \[\ \mathrm{C}_{\mathrm{v s}}(\mathrm{r})=\mathrm{C}_{\mathrm{v B}} \cdot \mathrm{e}^{-\mathrm{1 2} \cdot \frac{\mathrm{v}_{\mathrm{t h}}}{\beta_{\mathrm{sm}} \cdot \mathrm{k} \cdot \mathrm{u}_{*}} \cdot \frac{\mathrm{r}}{\mathrm{D_p}}}\]

    Based on integration, assuming open channel flow and settling in the vehicle, the bottom concentration CvB can be found, at the LDV CvB=Cvb:

    \[\ \mathrm{C}_{\mathrm{vB}}=\mathrm{C}_{\mathrm{vs}} \cdot \frac{\left(\frac{\mathrm{1 2} \cdot \mathrm{v}_{\mathrm{tv}}}{\beta_{\mathrm{sm}} \cdot \mathrm{k} \cdot \mathrm{u}_{*}}\right)}{\left(1-\mathrm{e}^{-12 \cdot \frac{\mathrm{v}_{\mathrm{tv}}}{\beta_{\mathrm{sm}} \cdot \kappa \cdot \mathrm{u}_{*}}}\right)}\]

    7.10.2 The Diffusivity Based on the LDV

    In order to have a consistent model, the bottom concentration should be equal to the bed concentration at the LDV, since this is the definition of the LDV. Now most concentration profile equations are not related to the LDV, but make use of a modified diffusivity to match experiments. Here an attempt is made to make the concentration profile match the LDV. At the Limit Deposit Velocity the bottom concentration CvB equals the bed concentration Cvb giving:

    \[\ \mathrm{C}_{\mathrm{v b}}=\mathrm{C _ { v s }} \cdot \frac{\left(\frac{\mathrm{1 2} \cdot \mathrm{v}_{\mathrm{t} \mathrm{v}, \mathrm{l d v}}}{\boldsymbol{\beta}_{\mathrm{s m}, \mathrm{ld} \mathrm{v}} \cdot \mathrm{\kappa} \cdot \mathrm{u}_{*, \mathrm{l d v}}}\right)}{\left(\mathrm{1 - e}^{-\mathrm{1 2} \cdot \frac{\mathrm{v}_{\mathrm{t} \mathrm{v}, \mathrm{ld v}}}{\beta_{\mathrm{sm}} \cdot \kappa \cdot \mathrm{u}_{*, \mathrm{ldv}}}}\right)}\]

    Neglecting the denominator at low concentrations, since it’s close to unity (say the denominator equals a factor αsm), the diffusivity can be derived.

    \[\ \beta_{\mathrm{sm}, \mathrm{ldv}}=12 \cdot \frac{\mathrm{C}_{\mathrm{vs}}}{\mathrm{C}_{\mathrm{vb}}} \cdot \frac{\mathrm{v}_{\mathrm{tv}, \mathrm{ldv}}}{\alpha_{\mathrm{sm}} \cdot \kappa \cdot \mathrm{u}_{*, \mathrm{ldv}}}=12 \cdot \mathrm{C}_{\mathrm{vr}} \cdot \frac{\mathrm{v}_{\mathrm{tv}, \mathrm{ldv}}}{\alpha_{\mathrm{sm}} \cdot \kappa \cdot \mathrm{u}_{*, \mathrm{ldv}}}\]

    This gives for the concentration distribution in the pipe:

    \[\ \mathrm{C}_{\mathrm{v} \mathrm{s}}(\mathrm{r})=\mathrm{C}_{\mathrm{v} \mathrm{B}} \cdot \mathrm{e}^{-\mathrm{1 2} \cdot \frac{\mathrm{v}_{\mathrm{tv}}}{\left(\mathrm{1 2} \cdot \mathrm{C}_{\mathrm{v r}} \cdot \frac{\mathrm{v}_{\mathrm{t v}, \mathrm{ldv}}}{\alpha_{\mathrm{sm}} \cdot \mathrm{\kappa} \cdot \mathrm{u}_{*, \mathrm{dv}}}\right) \cdot \mathrm{\kappa} \cdot \mathrm{u}_{*}} \cdot \frac{\mathrm{r}}{\mathrm{D}_{\mathrm{p}}}}{=\mathrm{C}_{\mathrm{v B}} \cdot \mathrm{e}}^{-\frac{\alpha_{\mathrm{sm}}}{\mathrm{C}_{\mathrm{vr}}}\cdot\frac{\mathrm{u_{*,ldv}}}{\mathrm{u_*}}\cdot \frac{\mathrm{v_{tv}}}{\mathrm{v_{tv,ldv}}}\cdot \frac{\mathrm{r}}{\mathrm{D_p}}}\]

    The bottom concentration is now, using the fact that the denominator equals the factor αsm:

    \[\ \mathrm{C}_{\mathrm{v B}}=\mathrm{C}_{\mathrm{v s}}\cdot \frac{\mathrm{\left( \frac{\alpha_{sm}}{C_{vr}}\cdot \frac{u_{*,ldv}}{u_*}\cdot \frac{v_{tv}}{v_{tv,ldv}} \right)}}{\mathrm{\left(1-e^{-\frac{\alpha_{sm}}{C_{vr}}\cdot \frac{u_{*,ldv}}{u_*}\cdot \frac{v_{tv}}{v_{tv,ldv}}} \right)}}=\mathrm{C_{vb} \cdot \frac{\left( \alpha_{sm}\cdot\frac{u_{*,ldv}}{u_*}\cdot \frac{v_{tv}}{v_{tv,ldv}} \right)}{\left(1-e^ {-\frac{\alpha_{sm}}{C_{vr}}\cdot \frac{u_{*,ldv}}{u_*}\cdot\frac{v_{tv}}{v_{tv,ldv}}} \right)}=C_{vb}\cdot\frac{u_{*,ldv}}{u_*}\cdot \frac{v_{tv}}{v_{tv,ldv}}}\]

    The correction factor has to be determined at the LDV, giving an implicit equation with only the relative volumetric concentration as the parameter:

    \[\ \alpha_{\mathrm{sm}}=\left(1-\mathrm{e}^{-\frac{\alpha_{\mathrm{sm}}}{\mathrm{C}_{\mathrm{vr}}}}\right)\]

    The correction factor appears to depend only on the relative concentration Cvr according to:

    \[\ \alpha_{\mathrm{sm}}=0.9847+0.304 \cdot \mathrm{C}_{\mathrm{vr}}-1.196 \cdot \mathrm{C}_{\mathrm{vr}}^{2}-0.5564 \cdot \mathrm{C}_{\mathrm{vr}}^{3}+0.47 \cdot \mathrm{C}_{\mathrm{vr}}^{4}\]

    At low relative concentrations, Cvr<0.3, this factor is about 1. Based on the diffusivity derived, the portion of the solids in the vehicle according to the Wasp criterion can be determined by:

    \[\ \mathrm{ \frac{C_{\mathrm{v s}, \mathrm{v}}}{\mathrm{C}_{\mathrm{v s}}}=\mathrm{e}^{-(0.92-\mathrm{0 . 5}) \cdot \frac{\alpha_{\mathrm{sm}}}{\mathrm{C}_{\mathrm{v r}}} \cdot \frac{\mathrm{u}_{*, \mathrm{ldv}}}{\mathrm{u}_{*}} \cdot \frac{\mathrm{v}_{\mathrm{tv}}}{\mathrm{v}_{\mathrm{tv}, \mathrm{ldv}}}}}\]

    7.10.3 Simplification of the Equations

    Now in reality, the concentration distribution depends on more than just the advection diffusion processes. At line speeds below the LDV there is a bed with possibly sheet flow and also collisions play a dominant role in the heterogeneous regime. Above the LDV the Darcy Weisbach friction factor is not changing much at higher line speeds. The power of the line speed is about -0.15 for larger pipe diameters. The power of the pipe diameter is about -0.175.

    \[\ \lambda_{\mathrm{l}}=0.01216 \cdot\left(\mathrm{v_{1 s}}\right)^{-0.15} \cdot\left(\mathrm{D_{p}}\right)^{-0.175}\]

    Substituting this for the friction velocities gives:

    \[\ \begin{array}{left} \mathrm{C}_{\mathrm{v s}}(\mathrm{r})&=\mathrm{C}_{\mathrm{v B}} \cdot \mathrm{e}^{-\frac{\alpha_{\mathrm{sm}}}{\mathrm{C}_{\mathrm{v r}}} \cdot \frac{\sqrt{\lambda_{\mathrm{l}, \mathrm{ldv}} / \mathrm{8}} \cdot \mathrm{v}_{\mathrm{ls}, \mathrm{ldv}}}{\sqrt{\lambda_{\mathrm{l}} / \mathrm{8}}{ \cdot \mathrm{v}_{\mathrm{ls}}}} \cdot \frac{\mathrm{v}_{\mathrm{tv}}}{\mathrm{v}_{\mathrm{tv}, \mathrm{ldv}}} \cdot \frac{\mathrm{r}}{\mathrm{D}_{\mathrm{p}}}}\\
    &=\mathrm{C_{vB}{ \cdot e^{-{\frac{\alpha_{sm}}{C_{vr}}\cdot \frac{\sqrt{0.01216 \cdot (v_{ls,ldv})^{-0.15}\cdot(D_p)^{-0.175}}\cdot v_{ls,ldv}}{\sqrt{0.01216 \cdot (v_{ls})^{-0.15}\cdot(D_p)^{-0.175}}\cdot v_{ls}}\cdot\frac{v_{tv}}{v_{tv,ldv}}\cdot\frac{r}{D_p}}}}}\\
    &=\mathrm{C_{vB}\cdot e^{-\frac{\alpha_{sm}}{C_{vr}}\cdot \left(\frac{v_{ls,ldv}}{v_{ls}} \right)^{0.925}\cdot \frac{v_{tv}}{v_{tv,ldv}}\cdot \frac{r}{D_p}}}\end{array}\]

    The bottom concentration is now for line speeds above the LDV:

    \[\ \mathrm{C}_{\mathrm{v B}}=\mathrm{C}_{\mathrm{v b}} \cdot\left(\frac{\mathrm{v}_{\mathrm{l} \mathrm{s}, \mathrm{ld} \mathrm{v}}}{\mathrm{v}_{\mathrm{l s}}}\right)^{0.925} \cdot \frac{\mathrm{v}_{\mathrm{t v}}}{\mathrm{v}_{\mathrm{t} \mathrm{v}, \mathrm{ld v}}} \quad\text{ or }\quad \mathrm{C}_{\mathrm{v B}}=\mathrm{C}_{\mathrm{v b}} \cdot\left(\frac{\mathrm{v}_{\mathrm{l} \mathrm{s}, \mathrm{ld} \mathrm{v}}}{\mathrm{v}_{\mathrm{l} \mathrm{s}}}\right)^{0.925} \cdot \frac{\mathrm{v}_{\mathrm{t h v}}}{\mathrm{v}_{\mathrm{t h v}, \mathrm{ld} \mathrm{v}}}\]

    For line speeds below the LDV there is a bed, so the concentration profile has to be adjusted for the presence of the bed. The portion of the solids in the vehicle according to the Wasp criterion can be determined by:

    \[\ \frac{\mathrm{C}_{\mathrm{vs}, \mathrm{v}}}{\mathrm{C}_{\mathrm{vs}}}=\mathrm{e}^{-(0.92-0.5) \cdot \frac{\alpha_{\mathrm{sm}}}{\mathrm{C}_{\mathrm{vr}}} \cdot\left(\frac{\mathrm{v}_{\mathrm{ls}, \mathrm{ldv}}}{\mathrm{v}_{\mathrm{ls}}}\right)^{0.925} \cdot \frac{\mathrm{v}_{\mathrm{tv}}}{\mathrm{v}_{\mathrm{tv}, \mathrm{ldv}}}} \quad\text{ or }\quad\frac{\mathrm{C}_{\mathrm{vs}, \mathrm{v}}}{\mathrm{C}_{\mathrm{vs}}}=\mathrm{e}^{-(0.92-0.5) \cdot \frac{\alpha_{\mathrm{sm}}}{\mathrm{C}_{\mathrm{vr}}} \cdot\left(\frac{\mathrm{v}_{\mathrm{ls}, \mathrm{ldv}}}{\mathrm{v}_{\mathrm{ls}}}\right)^{0.925} \cdot \frac{\mathrm{v}_{\mathrm{thv}}}{\mathrm{v}_{\mathrm{thv}, \mathrm{ldv}}}}\]

    The above equations give the solution for both hindered settling and non-hindered settling. If the hindered settling is based on the total volumetric concentration, it has no influence, since it gives a constant reduction factor. However if the hindered settling is based on the local volumetric concentration at a certain location in the pipe, it will influence the outcome.

    In chapter 6, the Wasp model, the behavior of the concentration profile based on the friction velocities is shown, compared with the models of Abulnaga (2002), Lahiri (2009) and Kaushal & Tomita (2002B). All models are based on a vertical concentration profile, without correction for the circular shape of a pipe, similar to the original Wasp et al. (1977) model for uniform sands. The advantage of this approach is that the concentration profiles can be determined analytically. The full Kaushal & Tomita (2002B) approach requires a numerical solution.

    Figure 7.10-1, Figure 7.10-3, Figure 7.10-5 and Figure 7.10-7 show the relative excess hydraulic gradient of the Wasp model for 6 different implementations as described in chapter 6. The DHLLDV implementation is based on the simplified equations, including hindered settling. Figure 7.10-2, Figure 7.10-4, Figure 7.10-6 and Figure 7.10-8 show the fractions in suspension of the 6 implementations. The calculations are carried out for uniform sands with particle diameters of d=0.1, 0.2, 0.5 and 1.0 mm in a 0.1524 m diameter pipe. For the 0.1 mm particle, the DHLLDV Framework gives a smaller suspended fraction compared to the other implementations.

    Figure 7.10-1: The Wasp model for a d=0.1 mm sand particle, modified diffusivity and hindered settling.

    Screen Shot 2020-07-30 at 8.56.19 PM.png

    Figure 7.10-2: The suspended fraction for a d=0.1 mm sand particle, modified diffusivity and hindered settling.

    Screen Shot 2020-07-30 at 8.57.10 PM.png

    Figure 7.10-3: The Wasp model for a d=0.2 mm sand particle, modified diffusivity and hindered settling.

    Screen Shot 2020-07-30 at 8.58.04 PM.png

    Figure 7.10-4: The suspended fraction for a d=0.2 mm sand particle, modified diffusivity and hindered settling.

    Screen Shot 2020-07-30 at 8.58.53 PM.png

    For the 0.2 mm particle, the DHLLDV implementation gives a larger suspended fraction compared to the other implementations above a line speed of about 4.5 m/sec. In terms of the Erhg value, the differences are small. In terms of the hydraulic gradient, the differences are very small.

    Figure 7.10-5: The Wasp model for a d=0.5 mm sand particle, modified diffusivity and hindered settling.

    Screen Shot 2020-07-30 at 9.00.20 PM.png

    Figure 7.10-6: The suspended fraction for a d=0.5 mm sand particle, modified diffusivity and hindered settling.

    Screen Shot 2020-07-30 at 9.02.00 PM.png

    For the 0.5 mm particle, the DHLLDV implementation gives a larger suspended fraction compared to the other implementations. In terms of the Erhg value, the differences are small. In terms of the hydraulic gradient, the differences are very small.

    Figure 7.10-7: The Wasp model for a d=1.0 mm sand particle, modified diffusivity and hindered settling.

    Screen Shot 2020-07-30 at 9.03.14 PM.png

    Figure 7.10-8: The suspended fraction for a d=1.0 mm sand particle, modified diffusivity and hindered settling.

    Screen Shot 2020-07-30 at 9.06.01 PM.png

    For the 1.0 mm particle, the DHLLDV implementation gives a much larger suspended fraction compared to the other implementations. In terms of the Erhg value, the differences are not small anymore. In terms of the hydraulic gradient, the differences are still small.

    7.10.4 Numerical Implementation

    The concentration in the pipe can be described according to:

    \[\ \mathrm{C}_{\mathrm{v s}}(\mathrm{r})=\mathrm{C}_{\mathrm{v B}} \cdot \mathrm{e}^{-\frac{\alpha_{\mathrm{sm}}}{\mathrm{C}_{\mathrm{v r}}} \cdot\left(\frac{\mathrm{v}_{\mathrm{l s}, \mathrm{ld} \mathrm{v}}}{\mathrm{v}_{\mathrm{l s}}}\right)^{0.925}} \cdot \frac{\mathrm{v}_{\mathrm{tv}}}{\mathrm{v}_{\mathrm{t} \mathrm{v}, \mathrm{ldv}}} \cdot \frac{\mathrm{r}}{\mathrm{D}_{\mathrm{p}}} \quad \text{or}\quad \mathrm{C}_{\mathrm{v s}}(\mathrm{r})=\mathrm{C}_{\mathrm{v B}} \cdot \mathrm{e}^{-\frac{\alpha_{\mathrm{sm}}}{\mathrm{C}_{\mathrm{vr}}} \cdot\left(\frac{\mathrm{v}_{\mathrm{ls}, \mathrm{ldv}}}{\mathrm{v}_{\mathrm{ls}}}\right)^{0.925} \cdot \frac{\mathrm{v}_{\mathrm{thv}}}{\mathrm{v}_{\mathrm{thv}, \mathrm{ldv}}} \cdot \frac{\mathrm{r}}{\mathrm{D}_{\mathrm{p}}}}\]

    The bottom concentration is now for line speeds above the LDV:

    \[\ \mathrm{C}_{\mathrm{v B}}=\mathrm{C}_{\mathrm{v b}} \cdot\left(\frac{\mathrm{v}_{\mathrm{l} \mathrm{s}, \mathrm{l d v}}}{\mathrm{v}_{\mathrm{l s}}}\right)^{0.925} \cdot \frac{\mathrm{v}_{\mathrm{t v}}}{\mathrm{v}_{\mathrm{t} \mathrm{v}, \mathrm{ld v}}} \quad\text{ or }\quad \mathrm{C}_{\mathrm{v B}}=\mathrm{C}_{\mathrm{v b}} \cdot\left(\frac{\mathrm{v}_{\mathrm{l} \mathrm{s}, \mathrm{l d v}}}{\mathrm{v}_{\mathrm{l s}}}\right)^{0.925} \cdot \frac{\mathrm{v}_{\mathrm{t h v}}}{\mathrm{v}_{\mathrm{t h v}, \mathrm{ld v}}}\]

    Now assuming that the terminal (hindered) settling velocity in the suspension lightly depends on the line speed, these equations can be written as:

    The concentration in the pipe can be described according to:

    \[\ \mathrm{C_{vs}(r)=C_{vB}\cdot e^{-\frac{\alpha_{sm}}{C_{vr}}\cdot \frac{u_{*,ldv}}{u_*}\cdot \frac{v_{tv}}{v_{tv,ldv}}\cdot\frac{r}{D_p}}\approx C_{vB}\cdot e^{-\frac{\alpha_{sm}}{C_{vr}}\cdot\left(\frac{v_{ls,ldv}}{v_{ls}} \right)^{1.15}\cdot\frac{r}{D_p}}} \]

    The bottom concentration is now for line speeds above the LDV:

    \[\ \mathrm{C_{vB}=C_{vb}\cdot \frac{\left(\alpha_{sm}\cdot \frac{u_{*,ldv}}{u_*}\cdot \frac{v_{tv}}{v_{tv,ldv}} \right)}{1-e^{-\frac{\alpha_{sm}}{C_{vr}}\cdot \frac{u_{*,ldv}}{u_*}\cdot \frac{v_{tv}}{v_{tv,ldv}}}}=C_{vb}\cdot \frac{\left(\alpha_{sm} \cdot \left(\frac{v_{ls,ldv}}{v_{ls}} \right)^{1.15} \right)}{\left(1-e^{-\frac{\alpha_{sm}}{C_{vr}}\cdot\left(\frac{v_{ls,ldv}}{v_{ls}} \right)^{1.15}} \right)}} \]

    These equations describe the concentration distribution well for 2D channel flow above the LDV. However for a circular pipe and below the LDV some adjustments have to be made. When the concentration found is integrated over the circular cross section of the pipe, the cross sectional averaged concentration has to be equal to the average concentration that is input to the calculations, which might be the case for a symmetrical concentration distribution, but certainly not for an asymmetrical concentration distribution.

    The concentration in the pipe can be described according to:

    \[\ \mathrm{C}_{\mathrm{v s}}(\mathrm{r})=\mathrm{C}_{\mathrm{v B}} \cdot \mathrm{e}^{-\frac{\alpha_{\mathrm{sm}}}{\mathrm{C}_{\mathrm{v r}}} \cdot\left(\frac{\mathrm{v}_{\mathrm{l s}, \mathrm{ldv}}}{\mathrm{v}_{\mathrm{l s}}}\right)^{1.15} \cdot \frac{\mathrm{r}}{\mathrm{D}_{\mathrm{p}}}}\]

    Now in the case of a circular pipe the vertical coordinate r/Dp should be replaced by the fraction of the cross section f. This gives a much better match with the cross sectional averaged concentration in case there is no bed, so above the LDV. This fraction can be determined by the angle β matching a certain vertical coordinate, similar to the angle β for the stationary and sliding bed.

    \[\ \mathrm{\beta=\operatorname{acos}\left(\frac{0.5-\frac{r}{D_{p}}}{0.5}\right)}\]

    The fraction f is now:

    \[\ \mathrm{f}=\frac{\beta-\sin (\beta) \cdot \cos (\beta)}{\pi}\]

    The concentration at f is now:

    \[\ \mathrm{C}_{\mathrm{v s}}(\mathrm{f})=\mathrm{C}_{\mathrm{v B}} \cdot \mathrm{e}^{-\frac{\alpha_{\mathrm{sm}}}{\mathrm{C}_{\mathrm{vr}}} \cdot\left(\frac{\mathrm{v}_{\mathrm{ls}, \mathrm{ld} \mathrm{v}}}{\mathrm{v}_{\mathrm{ls}}}\right)^{1.15}\cdot \mathrm{f}}\]

    The bottom concentration is now for line speeds above the LDV:

    \[\ \mathrm{C_{vB}=C_{vb}\cdot \frac{\left(\alpha_{sm}\cdot \left(\frac{v_{ls,ldv}}{v_{ls}} \right)^{1.15} \right)}{\left(1-e^{-\frac{\alpha_{sm}}{C_{vr}}\cdot \left(\frac{v_{ls,ldv}}{v_{ls}} \right)^{1.15}}\right)}} \]

    Because the Limit Deposit Velocity is based on the occurrence of some bed at the bottom of the pipe, this bed does not need to have the maximum bed density. A bed may start to occur with a bottom concentration of about 50%, while the maximum bed concentration will be in the range of 60%-65%. In order to find a bottom concentration of about 50% at the LDV assuming a maximum bed concentration of 60%, an additional velocity ratio rLDV is introduced giving:

    \[\ \mathrm{C_{vB}=C_{vb}\cdot \frac{\left( \alpha_{sm}\cdot \left(r_{LDV}{ \cdot \frac{v_{ls,ldv}}{v_{ls}}} \right)^{1.15} \right)}{\left(1-e^{-\frac{\alpha_{sm}}{C_{vr}}\cdot \left(r_{LDV}{ \cdot\frac{v_{ls,ldv}}{v_{ls}}} \right)^{1.15}} \right)} \quad\text{and}\quad C_{vs}(f)=C_{vB}\cdot e^{-\frac{\alpha_{sm}}{C_{vr}}\cdot \left(r_{LDV}{ \cdot\frac{v_{ls,ldv}}{v_{ls}}} \right)^{1.15}\cdot f}} \]

    The additional velocity ratio rLDV can be estimated by, considering the maximum LDV occurs at a concentration of 17.5%:

    \[\ \begin{array}{left}\text{SF = Shape Factor} \quad \text{SF=0.77 for sand}\quad \text{SF=1.0 for spheres}\\
    \mathrm{C_{\text {vrMax }}=\frac{0.175}{C_{v b}}}\\
    \alpha_{\beta}=1.8-56 \cdot \mathrm{v}_{\mathrm{t}} \quad\text{ with: }\quad \alpha_{\beta} \geq 1.1\\
    \text{If } \mathrm{C}_{\mathrm{v r}}<\mathrm{C}_{\mathrm{v r} \mathrm{M a x}} \text{ then}\\
    \mathrm{r_{L D V}=0.6 \cdot \frac{e^{(\beta / 2.34)^{\alpha_{\beta}}}}{e} \cdot\left(\frac{0.0005}{d}\right)^{S F^{6}} \cdot\left(\frac{C_{v r M a x}}{C_{v r}}\right)^{1 / 3} \quad\text{ with: }\quad r_{L D V} \geq 1.2 \cdot\left(\frac{C_{v r M a x}}{C_{v r}}\right)^{1 / 3}}\\
    \text{If } \mathrm{C}_{\mathrm{v r}} \geq \mathrm{C}_{\mathrm{v r M a x}}\text{ then }\\
    \mathrm{r_{L D V}=0.6 \cdot \frac{e^{(\beta / 2.34)^{\alpha_{\beta}}}}{e} \cdot\left(\frac{0.0005}{d}\right)^{S F^{6}} \cdot\left(\frac{C_{v r}}{C_{v r M a x}}\right)^{1 / 6} \quad\text{ with: }\quad r_{L D V} \geq 1.2 \cdot\left(\frac{C_{v r}}{C_{v r M a x}}\right)^{1 / 6}}\end{array}\]

    7.10.5 Examples Numerical Implementation

    In Figure 7.10-9, Figure 7.10-10 and Figure 7.10-11 the bottom concentration is limited to the bed concentration. If equation (7.10-11) gives a higher concentration, the top of the bed is determined at the level where the concentration equals the bed concentration and above the bed a new concentration profile is determined. This is repeated until a concentration profile is found matching the input spatial concentration.

    Figure 7.10-9: The concentration profiles at vls=0.25·vls,ldv.

    Screen Shot 2020-07-31 at 9.22.39 AM.png

    Figure 7.10-10: The concentration profiles at vls=0.50·vls,ldv.

    Screen Shot 2020-07-31 at 9.23.46 AM.png

    The profiles are corrected numerically for the circular shape of a pipe, resulting in the correct cross section averaged spatial concentration. In other words, the cross sectional averaged relative concentration is always equal to the input relative concentration, which it should be. This way the turbulent diffusion equation is not used to determine the concentration in the stationary or sliding bed, which would be inappropriate, since the bed behavior is soil mechanics and not fluid mechanics.

    Figure 7.10-11: The concentration profiles at vls=0.75·vls,ldv.

    Screen Shot 2020-07-31 at 9.25.33 AM.png

    Figure 7.10-12: The concentration profiles at vls=1.00·vls,ldv.

    Screen Shot 2020-07-31 at 9.26.39 AM.png

    Figure 7.10-9, Figure 7.10-10, Figure 7.10-11, Figure 7.10-12, Figure 7.10-13, Figure 7.10-14, Figure 7.10-15 and Figure 7.10-16 show the concentration profiles for relative concentrations, Cvr=Cvs/Cvb, ranging from 0.05 to 0.95 and line speeds of 0.25, 0.5, 0.75, 1.0, 1.5, 2.0, 3.0 and 4.0 times the LDV. When the line speed equals the LDV, the concentrations at the bottom of the pipe are equal to the bed concentration. Since the model is using the relative concentration and the relative line speed, the graphs are applicable for each type of material, particle size or solids density.

    Figure 7.10-13: The concentration profiles at vls=1.50·vls,ldv.

    Screen Shot 2020-07-31 at 9.28.56 AM.png

    Figure 7.10-14: The concentration profiles at vls=2.00·vls,ldv.

    Screen Shot 2020-07-31 at 9.30.49 AM.png

    Figure 7.10-15: The concentration profiles at vls=3.00·vls,ldv.

    Screen Shot 2020-07-31 at 9.31.59 AM.png

    Figure 7.10-16: The concentration profiles at vls=4.00·vls,ldv.

    Screen Shot 2020-07-31 at 9.33.32 AM.png

    The curves are compared with data from Kaushal et al. (2005) with d=0.44 mm and Dp=0.0549 m. The maximum LDV for this sand and pipe diameter is about 2.5 m/s. The data points are at 3 m/s, so a bit above the LDV depending on the concentration, resulting in smaller concentrations at the bottom of the pipe.

    Figure 7.10-17: Data from Kaushal et al. (2005).

    Screen Shot 2020-07-31 at 9.36.04 AM.png

    The data and the curves match pretty well, although at the top of the pipe the measured concentrations are a bit smaller. This is probably the effect of the local value of the hindered terminal settling velocity, which is higher at lower concentration, resulting in a higher downwards flux of the particles.

    Figure 7.10-18: Data from Kaushal et al. (2005) with local hindered settling.

    Screen Shot 2020-07-31 at 9.37.19 AM.png

    7.10.6 Hindered Settling Numerical Implementation

    The concentration profiles found match well except for the influence of local hindered settling, which will decrease the concentration at the top of the pipe and increase the concentration at the bottom of the pipe.

    The concentration in the pipe without the local hindered settling effect can be described according to:

    \[\ \mathrm{C}_{\mathrm{v B}}=\mathrm{C}_{\mathrm{v b}} \cdot \frac{\left(\alpha_{\mathrm{s m}} \cdot\left(\mathrm{r}_{\mathrm{L D V}} \cdot \frac{\mathrm{v}_{\mathrm{l s}, \mathrm{ld} \mathrm{v}}}{\mathrm{v}_{\mathrm{l s}}}\right)^{\mathrm{1 . 1 5}}\right)}{\left(\mathrm{1 - e}^{-\frac{\alpha_{\mathrm{sm}}}{\mathrm{C}_{\mathrm{v r}}} \cdot\left(\mathrm{r}_{\mathrm{L D V}} \cdot \frac{\mathrm{v}_{\mathrm{l s}, \mathrm{ld v}}}{\mathrm{v}_{\mathrm{l s}}}\right)^{1.15}}\right)} \quad\text{and}\quad \mathrm{C_{vs}(f)=C_{vB}\cdot e^{-\frac{\alpha_{sm}}{C_{vr}}\cdot \left(r_{LDV} \cdot \frac{v_{ls,ldv}}{v_{ls}} \right)^{1.15}\cdot f}}\]

    Hindered settling is not explicitly present in this equation. However, the Limit Deposit Velocity depends strongly on the concentration and on hindered settling, so implicitly hindered settling is present. The Limit Deposit Velocity increases with the concentration up to a concentration of 17.5%. For higher concentrations the Limit Deposit Velocity decreases with increasing concentration. This of course influences the concentration distribution.

    In the derivation of the concentration distribution it is assumed that the terminal settling velocity or the hindered terminal settling velocity is a constant, determined by the cross sectional average concentration. In reality the concentration is not a constant but decreases from top to bottom. So, the hindered terminal settling velocity at the top of the pipe will be higher than at the bottom, since at the top the concentration is lower. This could be implemented in the basic advection diffusion equation and solved numerically by iteration. However, one can also add a correction to adjust for this hindered settling effect. Such a correction should take the cross sectional averaged concentration into account, since the effect of hindered settling is larger at higher concentrations. After trial and error, the following method is found to take the effect of hindered settling into account:

    \[\ \begin{array}{left} \mathrm{C}_{\mathrm{v s}}(\mathrm{f}) \cdot \mathrm{v}_{\mathrm{t h}}(\mathrm{f})+\beta_{\mathrm{s m}} \cdot \varepsilon_{\mathrm{m}} \cdot \frac{\mathrm{d} \mathrm{C}_{\mathrm{v s}}(\mathrm{f})}{\mathrm{d f}}=\mathrm{0} \Rightarrow\\
    \frac{\mathrm{d} \mathrm{C}_{\mathrm{v s}}(\mathrm{f})}{\mathrm{d f}}=-\frac{\mathrm{C}_{\mathrm{v s} }(\mathrm{f}) \cdot \mathrm{v}_{\mathrm{t h}}(\mathrm{f})}{\beta_{\mathrm{s m}} \cdot \varepsilon_{\mathrm{m}}}=-\frac{\mathrm{C}_{\mathrm{v s}}(\mathrm{f}) \cdot \mathrm{v}_{\mathrm{t}} \cdot\left(\mathrm{1}-\mathrm{C}_{\mathrm{v s}}(\mathrm{f})\right)^{\beta}}{\beta_{\mathrm{sm}} \cdot \varepsilon_{\mathrm{m}}}\end{array}\]

    This is the advection diffusion equation with location dependent hindered settling. To determine the location dependent hindered settling, the local concentration must be known. In the zero step (index 0) the concentration profile is determined without local hindered settling based on equation (7.10-30). The location dependent concentration is already part of the solution and should not be considered in the first iteration (correction) step. Now at each level in the pipe the corrected concentration gradient can be determined according to:

    \[\ \left(\frac{\mathrm{d} \mathrm{C}_{\mathrm{vs}, 1}(\mathrm{f})}{\mathrm{df}}\right)=\left(\frac{\mathrm{d} \mathrm{C}_{\mathrm{vs}, 0}(\mathrm{f})}{\mathrm{df}}\right) \cdot\left(\frac{1-\mathrm{C}_{\mathrm{vs}, 0}(\mathrm{f})}{1-\mathrm{C}_{\mathrm{vs}}}\right)^{\beta}\]

    The left-hand side is the corrected concentration gradient, the first term on the right-hand side the concentration gradient determined with equation (7.10-30). The second term on the right-hand side gives the correction factor according to the Richardson & Zaki (1954) equation, but with a relative concentration dependent power. The power β is the Richardson & Zaki (1954) power with a value of 4.7 for very small particles and 2.34 for large particles. After determining the corrected concentration gradient at each level in the pipe, the new concentration profile is found by integrating the concentration gradient from bottom to top. It is assumed that the bottom concentration is unchanged. This process can be repeated several times until there is no significant change in the concentration profile. Here 12 iteration steps are used. For the next iteration steps (starting at index 2) also the location dependent concentration ratio is added, because it influences the concentration gradient, giving:

    \[\ \left(\frac{\mathrm{d} \mathrm{C}_{\mathrm{vs}, \mathrm{i}}(\mathrm{f})}{\mathrm{df}}\right)=\left(\frac{\mathrm{d} \mathrm{C}_{\mathrm{vs}, \mathrm{i}-1}(\mathrm{f})}{\mathrm{df}}\right) \cdot\left(\frac{\mathrm{C}_{\mathrm{vs}, \mathrm{i}-1}(\mathrm{f})}{\mathrm{C}_{\mathrm{vs}, \mathrm{i}-2}(\mathrm{f})}\right) \cdot\left(\frac{1-\mathrm{C}_{\mathrm{vs}, \mathrm{i}-1}(\mathrm{f})}{1-\mathrm{C}_{\mathrm{vs}, \mathrm{i}-2}(\mathrm{f})}\right)^{\beta}\]

    Although this method gives good results, there are still issues. One of the main issues is that the Richardson & Zaki (1954) hindered settling equation is based on the spatial volumetric concentration Cvs and not on the relative spatial volumetric concentration Cvr=Cvs/Cvb.

    \[\ \mathrm{ \frac{\mathrm{v}_{t h}}{\mathrm{v}_{t}}=\left(1-\mathrm{C}_{\mathrm{v s}}\right)^{\beta}}\]

    So even when the spatial volumetric concentration reaches a concentration where a bed with maximum porosity occurs, for sand at about Cvs=50%, still a hindered settling velocity is determined, while in reality this hindered settling velocity will be close to zero. Normal sands will have a porosity of about 40%, so Cvb=60%. A fixed bed may have a porosity of 40%, but a sliding bed will have a higher porosity in between 40% and 50%. The porosities mentioned here depend on the type of sand but are mentioned to give a feeling of the order of magnitude. An equation that may work better is:

    \[\ \frac{\mathrm{v}_{\mathrm{t h}}}{\mathrm{v}_{\mathrm{t}}}=\left(1-\mathrm{C}_{\mathrm{v r}}\right)^{\alpha \cdot \frac{\beta}{2.34}}\]

    For small concentrations this equation gives the same result as the original equation, but for concentrations approaching the bed concentration, this equation approaches a zero-settling velocity. This would describe the bed behavior much better. So, for small concentrations this equation describes hindered settling, while for large relative concentrations approaching 1, the behavior is closer to consolidation behavior. The power β in this equation is equal to the original power β.

    To determine the location dependent hindered settling, the local concentration must be known. In the zero step (index 0) the concentration profile is determined without local hindered settling based on equation (7.10-35). The location dependent concentration is already part of the solution and should not be considered in the first iteration (correction) step. Now at each level in the pipe the corrected concentration gradient can be determined according to:

    \[\ \left(\frac{\mathrm{d} \mathrm{C}_{\mathrm{vs}, 1}(\mathrm{f})}{\mathrm{df}}\right)=\left(\frac{\mathrm{d} \mathrm{C}_{\mathrm{vs}, 0}(\mathrm{f})}{\mathrm{df}}\right) \cdot\left(\frac{\left(1-\mathrm{C}_{\mathrm{vr}, 0}\right)}{\left(1-\mathrm{C}_{\mathrm{vr}}\right)}\right)^{\alpha \cdot \frac{\beta}{2.34}}\]

    The left-hand side is the corrected concentration gradient, the first term on the right-hand side the concentration gradient determined with equation (7.10-35). The second term on the right-hand side gives the correction factor according to equation (7.10-35), but with a relative concentration dependent power. The power β is the Richardson & Zaki (1954) power with a value of 4.7 for very small particles and 2.34 for large particles. After determining the corrected concentration gradient at each level in the pipe, the new concentration profile is found by integrating the concentration gradient from bottom to top. It is assumed that the bottom concentration is unchanged. This process can be repeated several times until there is no significant change in the concentration profile. Here 12 iteration steps are used. For the next iteration steps (starting at index 2) also the location dependent relative concentration ratio is added, because it influences the concentration gradient, giving:

    \[\ \left(\frac{\mathrm{d} \mathrm{C}_{\mathrm{vs}, \mathrm{i}}(\mathrm{f})}{\mathrm{df}}\right)=\left(\frac{\mathrm{d} \mathrm{C}_{\mathrm{vs}, \mathrm{i}-1}(\mathrm{f})}{\mathrm{df}}\right) \cdot\left(\frac{\mathrm{C}_{\mathrm{vr}, \mathrm{i}-1}(\mathrm{f})}{\mathrm{C}_{\mathrm{vr}, \mathrm{i}-2}(\mathrm{f})}\right) \cdot\left(\frac{\left(1-\mathrm{C}_{\mathrm{vr}, \mathrm{i}-1}\right)}{\left(1-\mathrm{C}_{\mathrm{vr}, \mathrm{i}-2}\right)}\right)^{\alpha \cdot \frac{\beta}{2.34}}\]

    Using this new equation gives significant different concentration profiles. The concentration profile of a sliding bed with sheet flow is simulated well with this equation. The power α is determined with the following equations:

    \[\ \begin{array}{left}\mathrm{S F}=\text{Shape Factor }\quad \mathrm{S F}=\mathrm{0 . 7 7} \text{ for sand }\quad \mathrm{S F}=\mathrm{1 . 0}\text{ for spheres }\\
    \mathrm{C}_{\text {vrMax }}=\frac{0.175}{\mathrm{C_{v b}}}\\
    \alpha=0.275 \cdot\left(\frac{\mathrm{SF}}{0.77}\right)^{1.5} \cdot\left(\frac{\mathrm{C}_{\mathrm{vr}}}{\mathrm{C}_{\mathrm{vrMax}}}\right)^{3} \cdot\left(\frac{\mathrm{v}_{\mathrm{Is}, \mathrm{LDV}}}{\mathrm{v}_{\mathrm{ls}}}\right)^{0.15}\text{ for }\mathrm{C}_{\mathrm{vr}}<\mathrm{C}_{\mathrm{vrMax}}\\
    \alpha=0.275 \cdot\left(\frac{\mathrm{SF}}{0.77}\right)^{1.5} \cdot\left(\frac{\mathrm{C}_{\mathrm{vr}}}{\mathrm{C}_{\mathrm{vrMax}}}\right)^{2 / 3} \cdot\left(\frac{\mathrm{v}_{\mathrm{Is}, \mathrm{LDV}}}{\mathrm{v}_{\mathrm{ls}}}\right)^{0.15}\text{ for }\mathrm{C}_{\mathrm{vr}} \geq \mathrm{C}_{\mathrm{vrMax}}\end{array}\]
    Figure 7.10-19: Concentration profiles with and without local hindered setting.

    Screen Shot 2020-07-31 at 11.45.07 AM.png

    Figure 7.10-19 shows a concentration profile for a line speed of 1 times the LDV and a relative concentration of 50% of the bed concentration. The red line shows the concentration profile as determined with equation (7.10-35). This results in the 12 iteration steps, where above the bottom the concentration is increased, while at the top the concentration has decreased. It is clear from this figure that 12 iteration steps give enough convergence.

    Figure 7.10-20, Figure 7.10-21, Figure 7.10-22, Figure 7.10-23, Figure 7.10-24, Figure 7.10-25, Figure 7.10-26 and Figure 7.10-27 show the concentration profiles with the effect of local hindered settling after applying the local hindered settling influence. Of course, this is a pragmatic solution, but it does give a closer fit to the data in Figure 7.10-18 compared to Figure 7.10-17. At high concentrations the concentration profile at the top of the pipe may not be accurate, especially at small line speeds. This is caused by the high concentrations at the bottom of the pipe. The standard graphs are determined with an average of β=2.7 for medium sized particles.

    To use these graphs, one should first determine the LDV at different concentrations for a specific particle diameter. Based on the LDV, the line speed to LDV ratio can be determined. Once this ratio is known, the graph the closest to the ratio should be chosen.

    7.10.7 Examples Hindered Settling Numerical Implementation

    Figure 7.10-20: The concentration profiles at vls=0.25·vls,ldv with hindered settling.

    Screen Shot 2020-07-31 at 11.47.39 AM.png

    Figure 7.10-21: The concentration profiles at vls=0.50·vls,ldv with hindered settling.

    Screen Shot 2020-07-31 at 11.48.28 AM.png

    The graphs are determined with a maximum bed concentration of 60%, which is reasonable for dense sand. For line speeds below the LDV sheet flow is observed by many researchers. Sheet flow is a layer of fast moving particles on top of a stationary or sliding bed. The particles closest to the bed move with the speed of the bed, while the particles at the top of the sheet flow layer move with a speed related to the speed of the liquid above the bed. The concentration is assumed to decrease almost linear, starting at the bed with the bed concentration, to almost zero at the top of the sheet flow layer. The start of the sheet flow layer is not at the maximum bed concentration, but at a concentration where the particles still form a solid bed, which would be at about 50% bed concentration, a relative concentration of about 0.83 or 83%. The graphs show an almost linear concentration profile from a relative concentration of 83% to a relative concentration of 10%, except for the highest overall relative concentrations of 90% and 95%. The steepness of the almost linear part of the concentration profiles increases with increasing overall (cross sectional averaged) relative concentration and decreases with decreasing line speed.

    Figure 7.10-22: The concentration profiles at vls=0.75·vls,ldv with hindered settling.

    Screen Shot 2020-07-31 at 11.50.35 AM.png

    Figure 7.10-23: The concentration profiles at vls=1.00·vls,ldv with hindered settling.

    Screen Shot 2020-07-31 at 11.51.41 AM.png

    Figure 7.10-24: The concentration profiles at vls=1.50·vls,ldv with hindered settling.

    Screen Shot 2020-07-31 at 11.52.35 AM.png

    Figure 7.10-25: The concentration profiles at vls=2.00·vls,ldv with hindered settling.

    Screen Shot 2020-07-31 at 11.53.39 AM.png

    The sheet flow behavior is simulated well with the new hindered settling equation and the method used. Above the LDV the bottom concentration is reduced depending on the line speed to LDV ratio. The concentration profiles become steeper with increasing line speed. In the range of the line speed to LDV ratio of 1-1.5 the concentration profiles change significantly. The sheet flow behavior completely disappears. Also the stationary or sliding bed completely disappears. There is not much difference between the use of the original Richardson & Zaki (1954) equation and equation (7.10-35), since they give the same hindered settling velocity at lower concentrations.

    Since the graphs are constructed for a hindered settling power of 3, one should only apply the graphs for medium sands. Fine sands will have a higher power, gravel a lower power. Richardson & Zaki (1954) derived their equation for spheres, so non spherical particles will have a different power. Usually non spherical particles have a higher power than spherical particles. So although the graphs are constructed for medium sands, they do give a good impression of the development of the concentration profiles as a function of the line speed to LDV ratio.

    Figure 7.10-26: The concentration profiles at vls=3.00·vls,ldv with hindered settling.

    Screen Shot 2020-07-31 at 11.56.36 AM.png

    Figure 7.10-27: The concentration profiles at vls=4.00·vls,ldv with hindered settling.

    Screen Shot 2020-07-31 at 11.57.55 AM.png

    7.10.8 Experiments

    Figure 7.10-28: Gillies (1993) data, d=0.29 mm, Dp=0.263 m, vls=4.0-4.7 m/s, β=3.45, LDV=3.43-4.04 m/s.

    Screen Shot 2020-07-31 at 12.00.58 PM.png

    Figure 7.10-29: Gillies (1993) data, d=0.38 mm, Dp=0.263 m, vls=4.1-4.8 m/s, β=3.21, LDV=3.55-4.22 m/s.

    Screen Shot 2020-07-31 at 12.02.52 PM.png

    Figure 7.10-28, Figure 7.10-29, Figure 7.10-30 and Figure 7.10-31 are based on experiments carried out by Gillies (1993). The experiments were carried out at two line speeds for each particle diameter. The bed concentrations are estimated at 0.55-0.60. In general the concentration profiles are similar for the different particle diameters, however the d=0.55 mm particle shows a very low concentration at the top of the pipe, which cannot be predicted.

    Figure 7.10-30: Gillies (1993) data, d=0.55 mm, Dp=0.263 m, vls=3.9-4.4 m/s, β=2.95, LDV=3.89-4.27 m/s.

    Screen Shot 2020-07-31 at 12.05.25 PM.png

    Figure 7.10-31: Gillies (1993) data, d=2.40 mm, Dp=0.263 m, vls=3.3-4.0 m/s, β=2.47, LDV=3.59-3.61 m/s.

    Screen Shot 2020-07-31 at 12.06.17 PM.png

    It is remarkable that the LDV approach, relating the concentration profile to the LDV instead of some diffusivity, gives such good results for the different particle diameters. The only difference between the particle diameters is the Richardson & Zaki power β. The downside of this approach is, that each concentration has a different LDV, so a different line speed to LDV ratio has to be used for each concentration. Here constant ratios are used for convenience.

    7.10.9 Concentration Distribution Sliding Flow Regime

    In the case of Sliding Flow, the bottom concentration decreases with increasing line speed and with decreasing spatial concentration. The bottom concentration can be determined with the following equation, where the bottom concentration can never be larger than the maximum bed concentration Cvb and never smaller than the spatial concentration Cvs. The Limit of Stationary Deposit Velocity (LSDV) has to be determined at a concentration of 17.5%, because the equation is calibrated for Cvs=0.175.

    \[\ \mathrm{C}_{\mathrm{v} \mathrm{B}}=\mathrm{3.1} \cdot \mathrm{C}_{\mathrm{v b}} \cdot\left(\frac{\mathrm{v}_{\mathrm{l} \mathrm{s}, \mathrm{l} \mathrm{s} \mathrm{d v}}}{\mathrm{v}_{\mathrm{l} \mathrm{s}}}\right)^{0.4} \cdot\left(\frac{\mathrm{C}_{\mathrm{v s}}}{\mathrm{C}_{\mathrm{v} \mathrm{b}}}\right)^{0.5} \cdot \mathrm{v}_{\mathrm{t}}^{1 / 6} \cdot \mathrm{e}^{\mathrm{D}_{\mathrm{p}}} \quad \mathrm{w} \mathrm{i t h}: \mathrm{1 . 1} \cdot \mathrm{C}_{\mathrm{v} \mathrm{s}} \leq \mathrm{C}_{\mathrm{v} \mathrm{B}} \leq \mathrm{C}_{\mathrm{v} \mathrm{b}}\]

    To determine the concentration distribution, the procedure outlined in chapter 7.10 should be followed with a line speed to Limit Deposit Velocity ratio of 1 and a bottom/bed concentration as determined with the above equation. Physically this means that a sliding bed will transit to sliding flow by increasing the porosity between the particles with increasing line speed. The Limit Deposit Velocity has no physical meaning in the Sliding Flow Regime.

    The concentration distribution is now:

    \[\ \begin{array}{left} \mathrm{C}_{\mathrm{v s}}(\mathrm{f})=\mathrm{C}_{\mathrm{v} \mathrm{B}} \cdot \mathrm{e}^{-\frac{\alpha_{\mathrm{s m}}}{\mathrm{C}_{\mathrm{v r}}} \cdot\mathrm{f}} \quad\text{ with: }\mathrm{C}_{\mathrm{v r}}=\frac{\mathrm{C}_{\mathrm{v} \mathrm{s}}}{\mathrm{C}_{\mathrm{v} \mathrm{B}}}\\
    \alpha_{\mathrm{sm}}=\mathrm{0.9 8 4 7}+\mathrm{0 .3 0 4} \cdot \mathrm{C}_{\mathrm{v r}}-\mathrm{1 .1 9 6} \cdot \mathrm{C}_{\mathrm{v r}}^{2}-\mathrm{0 .5 5 6 4} \cdot \mathrm{C}_{\mathrm{v r}}^{3}+\mathrm{0 .4 7} \cdot \mathrm{C}_{\mathrm{v r}}^{4}\end{array}\]

    Figure 7.10-32 shows an example of the concentration distribution in the Sliding Flow Regime.

    Figure 7.10-32: Gillies (1993) data, d=2.40 mm, Dp=0.0532 m, vls=1.8-3.1 m/s.

    Screen Shot 2020-07-31 at 12.16.59 PM.png

    7.10.10 Conclusions & Discussion

    The concentration profiles are based on an exponential power, containing the hindered terminal settling velocity. This hindered settling velocity is implemented as a constant over the vertical coordinate. Since in reality and also following from the above graphs, the concentration decreases with an increasing vertical coordinate, one may expect that the hindered settling velocity increases with the vertical coordinate. So in the top part the hindered settling velocity is higher than in the bottom part. This will result in slightly lower concentrations at the top and higher concentrations at the bottom, compared with the graphs without the local hindered settling velocity. The graphs with local hindered settling velocity have a correction for this effect. Also the repulsive effects near the pipe wall will influence the concentration profiles, which is not taken into account here.

    It is also assumed that the concentration profile is 1 dimensional with constant concentration in the transverse direction of the pipe. This is also questionable, since the repulsive effect near the pipe wall is effective at the whole pipe circumference and the local velocities are also influenced by the pipe wall. Recent CFD simulations have shown both effects (Ofei & Ismail (August 2016)).

    Since the model assumes a constant concentration in the transverse direction and is not corrected for the repulsive effect at the pipe wall, experimental data will never match exactly with the model. since usually the experimental data only show the measured concentrations on the vertical axis.

    Figure 7.10-28, Figure 7.10-29, Figure 7.10-30 and Figure 7.10-31 show the calculated curves, including the local hindered settling effect with data of Gillies (1993) in a Dp=0.263 m pipe. The experiments were carried out with particles of d=0.29 mm, d=0.38 mm, d=0.55 mm and d=2.4 mm. The line speeds varied from line speeds just below the LDV to line speeds just above the LDV. Since the LDV differs based on the particle diameter and the volumetric concentrations, these are given in the caption of each figure. In general the correlation is good, especially in the center of the pipe. At the top and bottom of the pipe the model may deviate from the experimental data. Figure 6.22-23, Figure 6.22-24, Figure 6.22-25, Figure 6.22-26, Figure 6.22-27, Figure 6.22-28, Figure 6.22-29 and Figure 6.22-30 show the data of Gillies (1993) in a Dp=0.0532 m pipe. These experiments were carried out with particles of d=0.18 mm, d=0.29 mm, d=0.55 mm and d=2.4 mm. The line speeds varied from 1.8 m/s to 3.1 m/s. The line speed of 1.8 m/s is close to the LDV depending on the particle diameter and the concentration.

    The current model gives a good estimate of the concentration distribution and is corrected for the cross sectional averaged concentration and the local hindered settling velocity, which is not the case for most models from literature.

    The choice of using the fraction of the cross sectional area of the pipe instead of the vertical coordinate seems to be a good choice. Even with the constant hindered settling velocity the concentration profiles match the experimental data much better. Using the adjustment of local hindered settling gives an even better correlation with the experimental data. The choice of an LDV based diffusivity gives good results. Of course this depends on the correctness of the determination of the LDV, but based on the comparison with experimental data, the model is close.

    It should be mentioned that the concentration profiles are determined with 100 steps, so a stepsize of 0.01·r/Dp. It should also be mentioned that the LDV values used are an average based on the DHLLDV Framework, while the LDV is concentration dependent. The DHLLDV Framework is a bit conservative giving high values for the LDV.

    To validate the model, two issues have to be addressed. The first issue is the value of the LDV. If the value of the LDV is not correct, the wrong line speed to LDV ratio will be used. The second issue is the shape of the concentration profile. Using the correct LDV, the shape of the concentration profile should match the experimental data. In general the shapes of the concentration profiles match well given the LDV values used. The LDV values from chapter 7.8 are a bit conservative, giving slightly to high values in order to have safe operations. Especially for the smallest particles of d=0.18 mm in a Dp=0.0532 m pipe, and LDV of 1.2 m/s is applied, while the theoretical LDV varies from 1.2 m/s to 1.8 m/s.

    Figure 7.10-33: The concentration profiles as a function of the line speed for Cvr=0.4.

    Screen Shot 2020-07-31 at 12.20.11 PM.png

    Figure 7.10-34: The concentration profiles as a function of the line speed for Cvr=0.6.

    Screen Shot 2020-07-31 at 12.21.11 PM.png

    Figure 7.10-33 and Figure 7.10-34 show the concentration profiles as a function of the line speed. A new term is introduced, the Maximum Limit Deposit Velocity. The reason for this is, that the observed LDV is based on the line speed where the bed completely disappeared. This will occur for a bed concentration close to 50%. In reality however, bed concentrations up to 60% are possible at lower line speeds. The line speed where this maximum bed concentration (60%) occurs is named the MLDV and the line speed where the 50% bed concentration occurs is named the LDV (the observed LDV). The LDV always occurs at a line speed higher than the MLDV. For a relative cross sectional averaged concentration of 0.4 (Figure 7.10-33) the LDV will occur at about 1.4 times the MLDV. For a relative cross sectional averaged concentration of 0.6 (Figure 7.10-34) the LDV will occur at about 1.9 times the MLDV. The higher the relative cross sectional averaged concentration, the larger this factor. The minimum value of this factor is about 1.25 for very low concentrations.

    Relating the concentration profile to the LDV is possible, but should be corrected for the fact that there is not just one bed concentration. A bed will start with the lowest possible bed concentration of about 50% or even a few percent lower and increases with decreasing line speed up to about 60% or even a few percent higher, depending on the type of sand, the PSD and the angularity.

    7.10.11 Nomenclature Concentration Distribution

    vehicle

    Pseudo liquid formed by carrier liquid and suspended particles

    -

    C

    Integration constant

    -

    Cvs(r)

    Volumetric concentration as a function of vertical coordinate in pipe

    -

    Cvs,v

    Volumetric concentration in vehicle

    -

    Cvb

    Volumetric concentration bed

    -

    CvB

    Volumetric concentration bottom of channel or pipe

    -

    Cvr

    Relative volumetric concentration

    -

    d

    Particle diameter

    m

    Dp

    Pipe diameter

    m

    g

    Gravitational constant (9.81 m/s2)

    m/s2

    r

    Vertical coordinate in pipe

    m

    \(\ \tilde{\mathrm{r}}\) Dimensionless vertical coordinate in pipe -

    R

    Radius of pipe

    m

    Rsd

    Relative submerged density (about 1.65 for sand)

    -

    u*

    Friction velocity

    m/s

    u*,ldv

    Friction velocity at LDV

    m/s

    vls

    Line speed

    m/s

    vls,ldv

    Limit Deposit Velocity

    m/s

    vt

    Terminal settling Velocity

    m/s

    vtv

    Terminal settling velocity in vehicle

    m/s

    vtv,ldv

    Terminal settling velocity in vehicle at LDV

    m/s

    vth

    Terminal settling velocity hindered

    m/s

    vth,ldv

    Terminal settling velocity hindered at LDV

    m/s

    vthv

    Terminal settling velocity hindered in vehicle

    m/s

    vthv,ldv

    Terminal settling velocity hindered in vehicle at LDV

    m/s

    αsm

    Correction factor

    -

    αc1

    Correction factor

    -

    αc2

    Correction factor

    -

    β

    Richardson & Zaki power hindered settling

    -

    βsm

    Relation sediment mass diffusivity to eddy momentum diffusivity.

    -

    βsm,ldv

    Relation sediment mass diffusivity to eddy momentum diffusivity at LDV

    -

    εm

    Eddy momentum diffusivity

    m/s

    εs

    Sediment diffusivity

    m/s

    \(\ \overline{\boldsymbol{\varepsilon}}_{\mathrm{S}}\)

    Average sediment diffusivity

    m/s

    λ1

    Darcy Weisbach friction factor

    -

    λ1,ldv

    Darcy Weisbach friction factor at LDV

    -

    κ

    von Karman constant

    -

    φ

    Pipe angle

    o


    This page titled 7.10: The Concentration Distribution is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Sape A. Miedema (TU Delft Open Textbooks) via source content that was edited to the style and standards of the LibreTexts platform; a detailed edit history is available upon request.