8.14: The Concentration Distribution
- Page ID
- 32295
\( \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}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)The fraction can be determined by the angle β matching a certain vertical coordinate, similar to the angle β for the stationary and sliding bed.
\[\ \beta=\operatorname{acos}\left(\frac{0.5-\frac{\mathrm{r}}{\mathrm{D_{p}}}}{0.5}\right) \Leftrightarrow \frac{\mathrm{r}}{\mathrm{D_{p}}}=\frac{1-\cos (\beta)}{2}\]
The fraction f is now:
\[\ \mathrm{f}=\frac{\beta-\sin (\beta) \cdot \cos (\beta)}{\pi}\]
The concentration distribution is (without or with hindered settling):
\[\ \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{l} \mathrm{d v}}}{\mathrm{v}_{\mathrm{ls}}}\right)^{0.925} \cdot \frac{\mathrm{v}_{\mathrm{tv}}}{\mathrm{v}_{\mathrm{tv}, \mathrm{ld v}}} \cdot \mathrm{f}}{\quad \text{or}\quad} \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 v}}}{\mathrm{v}_{\mathrm{ls}}}\right)^{0.925} \cdot \frac{\mathrm{v}_{\mathrm{thv}}}{\mathrm{v}_{\mathrm{th v,ldv}}}\cdot\mathrm{f}}\]
The correction factor appears to depend only on the relative concentration Cvr accordoing 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}\]
The bottom concentration is now (without or with hindered settling):
\[\ \mathrm{C_{vB}=C_{vb}\cdot \frac{\left(\alpha_{sm}\cdot \left(\frac{v_{ls,ldv}}{v_{ls}} \right)^{0.925}\cdot \frac{v_{tv}}{v_{tv,ldv}} \right)}{\left(1-e^{\frac{\alpha_{sm}}{C_{vr}}\cdot \left(\frac{v_{ls,ldv}}{v_{ls}} \right)^{0.925}\cdot \frac{v_{tv}}{v_{tv,ldv}}}\right)}\text{ or }C_{vB}=C_{vb}\cdot\frac{\left(\alpha_{sm}\cdot\left(\frac{v_{ls,ldv}}{v_{ls}} \right)^{0.925}\cdot \frac{v_{thv}}{v_{thv,ldv}} \right)}{\left(1-e^{-\frac{\alpha_{sm}}{C_{vr}}\cdot \left(\frac{v_{ls,ldv}}{v_{ls}} \right)^{0.925}\cdot \frac{v_{thv}}{v_{thv,ldv}}} \right)}} \]
If the bottom concentration is higher than the bed concentration, the concentration profile has to be adjusted. If it is assumed that the settling velocity in the suspension hardly changes as a function of the line speed, the equation for the concentration distribution becomes:
\[\ \mathrm{C}_{\mathrm{v s}, \mathrm{0}}(\mathrm{f})=\mathrm{C}_{\mathrm{v} \mathrm{B}} \cdot \mathrm{e}^{-\frac{\mathrm{\alpha}_{\mathrm{s m}}}{\mathrm{C}_{\mathrm{v r}}} \cdot\left(\frac{\mathrm{v}_{\mathrm{ls}, \mathrm{ld v}}}{\mathrm{v}_{\mathrm{l s}}}\right)^{1.15}\cdot \mathrm{f}}\]
The bottom concentration is now:
\[\ \mathrm{C}_{\mathrm{v B}}=\mathrm{C}_{\mathrm{v b}} \cdot \frac{\left(\alpha_{\mathrm{s m}} \cdot\left(\frac{\mathrm{v}_{\mathrm{l s}, \mathrm{ld} \mathrm{v}}}{\mathrm{v}_{\mathrm{l s}}}\right)^{1.15}\right)}{\left(1-\mathrm{e}^{-\frac{\alpha_{\mathrm{sm}}}{\mathrm{C}_{\mathrm{vr}}} \cdot\left(\frac{\mathrm{v}_{\mathrm{ls}, \mathrm{ldv}}}{\mathrm{v}_{\mathrm{ls}}}\right)^{1.15}}\right)}\]
At each level in the pipe the corrected concentration gradient can be determined for the first iteration step according to:
\[\ \left(\frac{\mathrm{d} \mathrm{C}_{\mathrm{v s}, \mathrm{1}}(\mathrm{f})}{\mathrm{d r}}\right)=\left(\frac{\mathrm{d} \mathrm{C}_{\mathrm{v s}, \mathrm{0}}(\mathrm{f})}{\mathrm{d r}}\right) \cdot\left(\frac{\left(1-\mathrm{C}_{\mathrm{v r}, 0}\right)}{\left(1-\mathrm{C}_{\mathrm{v r}}\right)}\right)^{\alpha \cdot \frac{\beta}{2.34}}\]
For the following iteration steps the concentration gradient has to be adjusted according to:
\[\ \left(\frac{\left.\mathrm{d} \mathrm{C}_{\mathrm{v} \mathrm{s}, \mathrm{i}} \mathrm{(} \mathrm{f}\right)}{\mathrm{d r}}\right)=\left(\frac{\left.\mathrm{d} \mathrm{C}_{\mathrm{v s}, \mathrm{i}-\mathrm{1}} \mathrm{( f}\right)}{\mathrm{d} \mathrm{r}}\right) \cdot\left(\frac{\mathrm{C}_{\mathrm{v r}, \mathrm{i}-\mathrm{1}}(\mathrm{f})}{\mathrm{C}_{\mathrm{v r}, \mathrm{i}-\mathrm{2}}(\mathrm{f})}\right) \cdot\left(\frac{\left(1-\mathrm{C}_{\mathrm{v} \mathrm{r}, \mathrm{i}-1}\right)}{\left(1-\mathrm{C}_{\mathrm{v r}, \mathrm{i}-2}\right)}\right)^{\alpha \cdot \frac{\beta}{2.34}}\]
Integrating the concentration profile again, starting at the bottom with either the bottom concentration (above the LDV) or the bed concentration (below the LDV), gives a concentration profile adjusted for local hindered settling. The power α is determined with the following equations:
\[\ \begin{array}{left}\text{SF = Shape Factor} \quad \text{SF=0.77 for sand} \quad \text{SF=1.0 for spheres}\\
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{ls}, \mathrm{LDV}}}{\mathrm{v}_{\mathrm{ls}}}\right)^{0.15} \quad \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{ls}, \mathrm{LDV}}}{\mathrm{v}_{\mathrm{ls}}}\right)^{0.15} \quad \mathrm{C}_{\mathrm{vr}} \geq \mathrm{C}_{\mathrm{vrMax}}\end{array}\]
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, 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_{s,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:
\[\ \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}}}\\
\mathrm{\alpha_{\beta}=1.8-56 \cdot v_{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} \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}}{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}\]
The factor 1.2 is based on the ratio 60% to 50%, the maximum bed concentration to the minimum bed concentration. For most of the experimental data analysed, a maximum bed concentration of 60% gives very good results. The LDV has a maximum at a concentration of 17.5%.