8.8: The Transition Heterogeneous Regime - Homogeneous Regime
- Page ID
- 32208
\( \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}\)8.8.1 Introduction
Some general equations that are required:
The Reynolds number in general is:
\[\ \mathrm{R} \mathrm{e}=\frac{\mathrm{v}_{\mathrm{l s}} \cdot \mathrm{D}_{\mathrm{p}}}{v_{\mathrm{l}}}\]
The Darcy Weisbach friction factor can now be determined with:
\[\ \lambda_{\mathrm{l}}=\frac{1.325}{\left(\ln \left(\frac{\varepsilon}{3.7 \cdot \mathrm{D}_{\mathrm{p}}}+\frac{5.75}{\mathrm{Re}^{0.9}}\right)\right)^{2}}=\frac{0.25}{\left(\log _{10}\left(\frac{\varepsilon}{3.7 \cdot \mathrm{D}_{\mathrm{p}}}+\frac{5.75}{\mathrm{Re}^{0.9}}\right)\right)^{2}}\]
The friction velocity is:
\[\ \mathrm{u}_{*}=\sqrt{\frac{\lambda_{\mathrm{l}}}{\mathrm{8}}} \cdot \mathrm{v}_{\mathrm{l} \mathrm{s}}\]
8.8.2 The Lift Ratio
The near wall lift pushes particles away from the pipe wall. The ratio of the lift force to the sum of submerged weight and momentum of a particle is named the Lift Ratio.
\[\ \mathrm{L}_{\mathrm{R}}=\frac{\mathrm{C}_{\mathrm{L}} \cdot \mathrm{u}_{*}^{2}}{\mathrm{d} \cdot\left(\mathrm{R}_{\mathrm{s d}} \cdot \mathrm{g}+\frac{1}{2} \cdot \frac{\rho_{\mathrm{s}}}{\rho_{\mathrm{l}}} \cdot \frac{\mathrm{v}_{\mathrm{t h}}^{2} \cdot \mathrm{u}_{*}}{\boldsymbol{\alpha} \cdot \mathrm{1 1 .6} \cdot v_{\mathrm{l}}}\right)} \cdot\left(\mathrm{1}-\frac{\mathrm{C}_{\mathrm{v} \mathrm{s}}}{\mathrm{C}_{\mathrm{v} \mathrm{b}}}\right)\]
8.8.3 The Heterogeneous Equation
The heterogeneous losses reduce (collapse) if the near wall lift pushes harder than the sum of the submerged weight and the momentum of a particle, resulting in:
For LR<0.7 the heterogeneous equation becomes:
\[\ \mathrm{E}_{\mathrm{rhg}, \mathrm{He}}=\frac{\mathrm{v}_{\mathrm{t}} \cdot\left(1-\frac{\mathrm{C}_{\mathrm{vs}}}{\mathrm{0 . 1 7 5} \cdot(1+\beta)}\right)^{\beta}}{\mathrm{v}_{\mathrm{ls}}}+\mathrm{8 .5}^{2} \cdot\left(\frac{\mathrm{1}}{\lambda_{\mathrm{l}}}\right) \cdot\left(\frac{\mathrm{v}_{\mathrm{t}}}{\sqrt{\mathrm{g} \cdot \mathrm{d}}}\right)^{10 / 3} \cdot\left(\frac{\left(v_{\mathrm{l}} \cdot \mathrm{g}\right)^{1 / 3}}{\mathrm{v}_{\mathrm{ls}}}\right)^{2} \cdot\left(1-\mathrm{L}_{\mathrm{R}}^{2}\right)\]
When the lift ratio has a value close to 1, theoretically there are no more collisions with the wall. However not all particles will have exactly the same kinetic energy, so even when the lift ratio is larger than 1, still some particles will have collisions. Therefore a smoothing function is proposed for lift ratio’s larger than 70% (ζ=0.5), giving:
\[\ \mathrm{E}_{\mathrm{rhg}, \mathrm{He}}=\frac{\mathrm{v}_{\mathrm{t}} \cdot\left(1-\frac{\mathrm{C}_{\mathrm{vs}}}{\mathrm{0 . 1 7 5} \cdot(1+\beta)}\right)^{\beta}}{\mathrm{v}_{\mathrm{ls}}}+\mathrm{8 .5}^{2} \cdot\left(\frac{\mathrm{l}}{\lambda_{\mathrm{l}}}\right) \cdot\left(\frac{\mathrm{v}_{\mathrm{t}}}{\sqrt{\mathrm{g} \cdot \mathrm{d}}}\right)^{10 / 3} \cdot\left(\frac{\left(v_{\mathrm{l}} \cdot \mathrm{g}\right)^{1 / 3}}{\mathrm{v}_{\mathrm{ls}}}\right)^{2} \cdot(1-\zeta) \cdot \frac{\zeta}{\mathrm{L}_{\mathrm{R}}^{2}}\]
8.8.4 The Homogeneous Equation
The mobilization factor of the homogeneous equation is:
\[\ \mathrm{m}=\mathrm{e}^{-\mathrm{0} \cdot \mathrm{1} \cdot \frac{\alpha_{\mathrm{sm}}}{\mathrm{C}_{\mathrm{v r}}} \cdot \frac{\mathrm{u}_{\mathrm{*,ldv}}}{\mathrm{u}_{*}} \cdot \frac{\mathrm{v}_{\mathrm{t h}}}{\mathrm{v}_{\mathrm{t h}, \mathrm{ldv}}}}\]
Basically this shows the concentration gradient at the center of the pipe. The mobilized homogeneous flow Erhg is now:
\[\ \mathrm{E_{rhg,Ho}=m \cdot i_l \cdot\left(1-\left(1-\frac{1+ R_{sd} \cdot C_{vs}-\left(\frac{A_{C_v}}{\kappa}\cdot ln \left(\frac{\rho_m}{\rho_l} \right)\cdot \sqrt{\frac{\lambda_l}{8}}+1 \right)^2}{R_{sd}\cdot C_{vs}\cdot \left(\frac{A_{C_v}}{\kappa}\cdot ln \left(\frac{\rho_m}{\rho_l} \right)\cdot \sqrt{\frac{\lambda_l}{8}}+1 \right)^2} \right) \left(1-\left(\frac{\delta_v}{d}\right) \right)\right)} \]
8.8.5 The Resulting Relative Excess Hydraulic Gradient
The resulting relative excess hydraulic gradient near the transition area can now be determined by:
\[\ \mathrm{E}_{\mathrm{r h g}, \mathrm{H e H o}}=\mathrm{E}_{\mathrm{r h g}, \mathrm{H e}}+\mathrm{E}_{\mathrm{r h g}, \mathrm{H o}}\]
The hydraulic gradient resulting from this relative excess hydraulic gradient is now:
\[\ \mathrm{i}_{\mathrm{m}, \mathrm{H} \mathrm{e} \mathrm{H} \mathrm{o}}=\mathrm{i}_{\mathrm{l}}+\mathrm{E}_{\mathrm{r} \mathrm{h} \mathrm{g}, \mathrm{H} \mathrm{e} \mathrm{H} \mathrm{o}} \cdot \mathrm{R}_{\mathrm{s} \mathrm{d}} \cdot \mathrm{C}_{\mathrm{v} \mathrm{s}}\]
This hydraulic gradient is now valid for the combination of the heterogeneous flow regime and the homogeneous flow regime, but still has to be combined with the stationary bed regime and the sliding bed regime. In the case of d/Dp>0.015 the sliding flow regime will occur and the above equation is not used.
The pressure difference in the transition zone can now be determined with:
\[\ \Delta \mathrm{p}_{\mathrm{m}, \mathrm{H} \mathrm{e} \mathrm{H} \mathrm{o}}=\Delta \mathrm{p}_{\mathrm{l}}+\rho_{\mathrm{l}} \cdot \mathrm{g} \cdot \Delta \mathrm{L} \cdot \mathrm{E}_{\mathrm{r h g}, \mathrm{H} \mathrm{e H o}} \cdot \mathrm{R}_{\mathrm{s d}} \cdot \mathrm{C}_{\mathrm{v s}}\]