Motivation
One of the key challenges in Pipe Flow Dynamics is to predict the pressure distribution along the pipe during the stationary fluid transport.
In many practical cases the stationary pressure distribution can be approximated by Isothermal or Quasi-isothermal homogenous fluid flow model.
Pipeline Flow Pressure Model is addressing this problem with account of the varying pipeline trajectory, gravity effects and fluid friction with pipeline walls.
...
Numerical quadrature solution of Pressure Profile in Homogeneous Steady-State Pipe Flow @model
Outputs
...
Inputs
...
Assumptions
...
Stationary Homogenous flow | Isothermal or conditionsIncompressible fluid \rho(T, p)=\rho_s = \rm const | Isoviscous flow --uriencoded--\displaystyle \frac%7B\partial p%7D%7B\partial t%7D = 0 |
|
|
\mu(T, p) = \mu_s = \rm const--uriencoded--\displaystyle \frac%7B\partial T%7D%7B\partial t%7D =0 \rightarrow T(t,l) = T(l) |
|
|
Homogenous flow | |
...
Pressure profile | Pressure gradient profile |
---|
LaTeX Math Block |
---|
anchor | PPconst |
---|
alignment | left |
---|
|
j_m = j_m[p(l)] \rightarrow p = p(l) |
LaTeX Math Block |
---|
|
\frac{dp}{dl} = {\rm Numerical} \ {\rm Derivative} |
Mass Flux | Mass Flowrate |
LaTeX Math Block |
---|
anchor | MassFlux |
---|
alignment | left |
---|
|
j_m = \sqrt{ 2 \cdot \frac
{
g \, \Delta z + \int_p^{p_0} \frac{dp}{\rho}
}
{
\left( \frac{1}{\rho^2} - \frac{1}{\rho_0^2} \right)
+ \left( \frac{1}{\rho^2} + \frac{1}{\rho_0^2} \right)
\cdot \frac{f \, \cdot \, l}{ 2 \, d}
}
}
|
LaTeX Math Block |
---|
anchor | MassFlowrate |
---|
alignment | left |
---|
|
\dot m = j_m \cdot A =
A \cdot \sqrt{ 2 \cdot \frac
{
g \, \Delta z + \int_p^{p_0} \frac{dp}{\rho}
}
{
\left( \frac{1}{\rho^2} - \frac{1}{\rho_0^2} \right)
+ \left( \frac{1}{\rho^2} + \frac{1}{\rho_0^2} \right)
\cdot \frac{f \, \cdot \, l}{ 2 \, d}
}
} |
Volumetric Flowrate | Intake Fluid velocity |
LaTeX Math Block |
---|
anchor | VolumtericFlowrate |
---|
alignment | left |
---|
|
q_s = \dot m / \rho_s =
\frac{A}{\rho_s} \cdot \sqrt{ 2 \cdot \frac
{
g \, \Delta z + \int_p^{p_0} \frac{dp}{\rho}
}
{
\left( \frac{1}{\rho^2} - \frac{1}{\rho_0^2} \right)
+ \left( \frac{1}{\rho^2} + \frac{1}{\rho_0^2} \right)
\cdot \frac{f \, \cdot \, l}{ 2 \, d}
}
}
|
LaTeX Math Block |
---|
anchor | IntakeFluidVelocity |
---|
alignment | left |
---|
|
u_s = j_m/ \rho_s =q_s / A =
\frac{1}{\rho_s} \cdot \sqrt{ 2 \cdot \frac
{
g \, \Delta z + \int_p^{p_0} \frac{dp}{\rho}
}
{
\left( \frac{1}{\rho^2} - \frac{1}{\rho_0^2} \right)
+ \left( \frac{1}{\rho^2} + \frac{1}{\rho_0^2} \right)
\cdot \frac{f \, \cdot \, l}{ 2 \, d}
}
}
|
where
body | --uriencoded--\displaystyle \frac%7B\partial p%7D%7B\partial \tau_x%7D =\frac%7B\partial p%7D%7B\partial \tau_y%7D =0 \rightarrow p(t, \tau_x,\tau_y,l) = p(l) |
---|
|
| |
Constant inclination | Constant friction along hole |
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle \theta(l) = \theta = %7B\rm const%7D \rightarrow \cos \theta = \frac%7Bdz%7D%7Bdl%7D = %7B\rm const%7D |
---|
|
| |
Equations
...
Pressure profile along the pipe |
---|
LaTeX Math Block |
---|
anchor | PressureProfile |
---|
alignment | left |
---|
| L =\int_{p_0}^{p} \frac{\rho \, dp}{G \, \rho^2 - F} -\frac{d}{f} \cdot \ln \frac{F/\rho^2 - G}{ F/\rho_0^2-G} |
|
where
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle j_m = \frac%7B \dot m %7D%7B A%7D = \rm const |
---|
|
| mass flux |
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle \dot m = \frac%7Bdm %7D%7B dt%7D = \rm const |
---|
|
| mass flowrate |
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle q_0 = \frac%7BdV_0%7D%7Bdt%7D = \frac%7B \dot m %7D%7B \rho_0%7D |
---|
|
| Intake volumetric flowrate |
LaTeX Math Inline |
---|
body | \rho_0 = \rho(T_0, p_0) |
---|
|
| Intake fluid density |
| Intake mass flux |
| mass flowrate |
| Intake Fluid velocity |
LaTeX Math Inline |
---|
body | \Delta z(l) = z(l)-z(0) |
---|
|
| elevation drop along pipe trajectory |
|
_s _s at intake point |
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle %7B\rm Re%7D |
---|
|
|
_s = \frac%7Bu(l) (T,\rho) = \frac%7Bj_m \cdot d%7D%7B\ |
|
|
nu(l)%7D = \frac%7B4 \rho_s q_s%7D%7B\pi d%7D \frac%7B1%7D%7B\mu_s%7DReynolds number at intake point | Reynolds number in Pipe Flow |
| dynamic viscosity as function of fluid temperature and density |
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle d = \sqrt%7B \frac%7B4 A%7D%7B\pi%7D%7D = \rm const |
---|
|
| characteristic linear dimension of the pipe (or exactly a pipe diameter in case of a circular pipe) |
LaTeX Math Inline |
---|
body | G = g \, \cos \theta = \rm const |
---|
|
| gravity acceleration along pipe |
LaTeX Math Inline |
---|
body | --uriencoded--F = j_m%5e2 \cdot f/(2d) = \rm const |
---|
|
|
|
...
Alternative forms
...
...
...
| L = L(\rho) =\int_{\rho_0}^{\rho} \frac{ 1}{G \, \rho^2 - F} \, \frac{ d\rho}{c(\rho)} -\frac{d}{f} \cdot \ln \frac{F/\rho^2 - G}{ F/\rho_0^2-G}
|
|
Approximations
...
LaTeX Math Block |
---|
| \rho(p) = \tilde \rho \cdot \sqrt{ 1- \frac{f}{2d} \frac{j_m^2}{G} ( \rho_0^2 - {\tilde \rho}^2) }
|
|
LaTeX Math Block |
---|
| \int_{\rho_0}^{\tilde \rho} \frac{d \rho}{ \rho^2 \, c(\rho)} = G \cdot L
|
|
where
| no-flow pressure at the pipe end ( LaTeX Math Inline |
---|
body | --uriencoded--%7B\tilde j%7D_m = 0 |
---|
| ) |
...
LaTeX Math Inline |
---|
body | f(l) = f_s = \rm const |
---|
|
...
See also
References
...
...
bgColor | papayawhip |
---|
title | ARAX |
---|
...