One-dimensional flow with friction
Flow-station data
The starting point is the governing equations for one-dimensional steady-state flow
Continuity
Momentum
where is the average wall-shear stress
is the tube perimeter, and is the tube length. For circular cross sections
and thus
Energy
Differential Form
In order to remove the integral term in the momentum equation, the governing equations are written in differential form\\
Continuity
\[\rho_1 u_1=\rho_2 u_2=const\Rightarrow\]\\
\begin{equation} \frac{d}{dx}(\rho u)=0 \label{eq:governing:cont} \end{equation}\\
Momentum
\[(\rho_2 u_2^2+p_2-\rho_1 u_1^2+p_1)=-\frac{4}{D}\int_0^L\tau_w dx\Rightarrow\]\\
\begin{equation} \frac{d}{dx}(\rho u^2+p)=-\frac{4}{D}\tau_w \label{eq:governing:mom:c} \end{equation}\\
\[\frac{d}{dx}(\rho u^2+p)=\rho u\frac{du}{dx}+u\frac{d}{dx}(\rho u)+\frac{dp}{dx}=\left\{\frac{d}{dx}(\rho u)=0\right\}=\rho u\frac{du}{dx}+\frac{dp}{dx}\]\\
\begin{equation} \rho u\frac{du}{dx}+\frac{dp}{dx}=-\frac{4}{D}\tau_w \label{eq:governing:mom:d} \end{equation}\\
\noindent The wall shear stress is often approximated using a shear-stress factor, $f$, according to\\
\begin{equation} \tau_w=f\frac{1}{2}\rho u^2 \label{eq:tauw:b} \end{equation}\\
\noindent and thus\\
\begin{equation} \rho u\frac{du}{dx}+\frac{dp}{dx}=-\frac{2}{D}f\rho u^2 \label{eq:governing:mom:d} \end{equation}\\
Energy
\[h_1 + \frac{1}{2}u_1^2=h_2 + \frac{1}{2}u_2^2=const\]\\
\[h_{o_1}=h_{o_2}=const\]\\
\begin{equation} \frac{d}{dx}h_o=0 \label{eq:governing:energy} \end{equation}\\
Summary
\noindent continuity:\\
\[\frac{d}{dx}(\rho u)=0\]\\
\noindent momentum:\\
\[\rho u\frac{du}{dx}+\frac{dp}{dx}=-\frac{2}{D}f\rho u^2\]\\
\noindent energy:\\
\[\frac{d}{dx}h_o=0\]
\noindent From chapter 3.9 we have the following expression for the momentum equation for one-dimensional flow with friction (equation (3.95))\\
\begin{equation*} dp+\rho u du=-\frac{1}{2}\rho u^2 \frac{4 f dx}{D}\ (3.95) \end{equation*}\\
\noindent For cases dealing with calorically perfect gas, (3.95) can be recast completely in terms of Mach number using the following relations\\
\begin{itemize} \item[] speed of sound: $a^2=\gamma p/\rho$ \item[] the definition of Mach number: $M^2=u^2/a^2$ \item[] the ideal gas law for thermally perfect gas: $p=\rho R T$ \item[] the continuity equation: $\rho u=const$ \item[] energy equation: $c_p T+u^2/2=const$ \end{itemize}
Continuity equation
\noindent We start with the continuity equation which for one-dimensional steady flows reads\\
\begin{equation} \rho u=const \label{eqn:cont:a} \end{equation}\\
\noindent Differentiating (\ref{eqn:cont:a}) gives\\
\begin{equation} d(\rho u)=0. \Leftrightarrow \rho du + ud\rho=0. \label{eqn:cont:b} \end{equation}\\
\noindent If $u\neq 0.$ we can divide by $\rho u$ which gives us\\
\begin{equation} \frac{du}{u}+\frac{d\rho}{\rho}=0. \label{eqn:cont:b} \end{equation}\\
\noindent Now, if we divide and multiply the first term in (\ref{eqn:cont:b}) by $2u$ and use the chain rule for derivatives we get\\
\begin{equation} \frac{d(u^2)}{2u^2}+\frac{d\rho}{\rho}=0. \label{eqn:cont} \end{equation}\\
Energy equation
\noindent For an adiabatic one-dimensional flow we have that \\
\begin{equation} c_p T+\frac{u^2}{2}=const \label{eqn:ttot:a} \end{equation}\\
\noindent If we differentiate (\ref{eqn:ttot:a}) we get\\
\begin{equation} c_p dT+\frac{1}{2}d(u^2)=0. \label{eqn:ttot:b} \end{equation}\\
\noindent We replace $c_p$ with $\gamma R/(\gamma-1)$ and multiply and divide the first term with $T$ which gives us\\
\begin{equation} \frac{\gamma RT}{(\gamma-1)}\frac{dT}{T}+\frac{1}{2}d(u^2)=0. \label{eqn:ttot:c} \end{equation}\\
\noindent Now, divide by $\gamma RT/(\gamma-1)$ and multiply and divide the second term by $u^2$ gives\\
\begin{equation} \frac{dT}{T}+\frac{(\gamma-1)}{2}M^2\frac{d(u^2)}{u^2}=0. \label{eqn:ttot} \end{equation}\\
\noindent We want to remove the $dT/T$-term in (\ref{eqn:ttot}). From the definition of Mach number we have that\\
\begin{equation} a^2M^2=u^2 \label{eqn:Mach:a} \end{equation}\\
\noindent which we can rewrite using the expression for speed of sound $(a^2=\gamma RT)$ according to\\
\begin{equation} \gamma RTM^2=u^2 \label{eqn:Mach:b} \end{equation}\\
\noindent Differentiating (\ref{eqn:Mach:b}) gives us\\
\begin{equation} \gamma RM^2 dT+\gamma RT d(M^2)=d(u^2) \label{eqn:Mach:c} \end{equation}\\
\noindent Now, if we divide (\ref{eqn:Mach:c}) by $\gamma RT M^2$ and use $a^2=\gamma RT$ and $a^2M^2=u^2$ we get\\
\begin{equation} \frac{dT}{T}+\frac{d(M^2)}{M^2}=\frac{d(u^2)}{u^2} \label{eqn:Mach} \end{equation}\\
\noindent Equation (\ref{eqn:Mach}) may now be used to replace the $dT/T$-term in equation (\ref{eqn:ttot})\\
\begin{equation} -\frac{d(M^2)}{M^2}+\frac{d(u^2)}{u^2}+\frac{(\gamma-1)}{2}M^2\frac{d(u^2)}{u^2}=0. \label{eqn:ttot:Mach:a} \end{equation}\\
\noindent which can be rewritten according to\\
\begin{equation} \frac{d(u^2)}{u^2}=\left[1+\frac{(\gamma-1)}{2}M^2\right]^{-1}\frac{d(M^2)}{M^2} \label{eqn:ttot:Mach:b} \end{equation}\\
\noindent Using the chain rule for derivatives, the last term may be rewritten according to\\
\[\frac{d(M^2)}{M^2}=2M\frac{dM}{M^2}=2\frac{dM}{M}\]\\
\noindent which gives\\
\begin{equation} \frac{d(u^2)}{u^2}=2\left[1+\frac{(\gamma-1)}{2}M^2\right]^{-1}\frac{dM}{M} \label{eqn:ttot:Mach} \end{equation}\\
\subsection{The ideal gas law} \noindent For a perfect gas the ideal gas law reads\\
\begin{equation} p=\rho R T \label{eqn:gaslaw:a} \end{equation}\\
\noindent Differentiating (\ref{eqn:gaslaw:a}) gives:\\
\begin{equation} dp=\rho R dT+RT d\rho \label{eqn:gaslaw:b} \end{equation}\\
\noindent If $p\neq0.$, we can divide (\ref{eqn:gaslaw:b}) by $p$ which gives\\
\begin{equation} \frac{dp}{p}=\frac{dT}{T}+\frac{d\rho}{\rho} \label{eqn:gaslaw:c} \end{equation}\\
\noindent which can be rearranged according to\\
\begin{equation} \left[\frac{dp}{p}-\frac{d\rho}{\rho}\right]=\frac{dT}{T} \label{eqn:gaslaw:d} \end{equation}\\
\noindent Now, inserting $dT/T$ from equation (\ref{eqn:ttot}) gives\\
\begin{equation} \left[\frac{dp}{p}-\frac{d\rho}{\rho}\right]+\frac{(\gamma-1)}{2}M^2\frac{d(u^2)}{u^2}=0. \label{eqn:gaslaw:b} \end{equation}\\
\noindent The $d\rho/\rho$-term can be replaced using equation (\ref{eqn:cont})\\
\begin{equation} \frac{dp}{p}+\frac{d(u^2)}{2u^2}+\frac{(\gamma-1)}{2}M^2\frac{d(u^2)}{u^2}=0. \label{eqn:gaslaw:c} \end{equation}\\
\noindent Collect terms and rewrite gives\\
\begin{equation} \frac{dp}{p}+\left[\frac{1+(\gamma-1)M^2}{2}\right]\frac{d(u^2)}{u^2}=0. \label{eqn:gaslaw} \end{equation}\\
Momentum equation
By combining the above derived relations and the momentum equation on the form given by (3.95), we can get an expression where the friction force is a function of Mach number only\\
\noindent For convenience equation (3.95) is written again here\\
\[dp+\rho u du=-\frac{1}{2}\rho u^2 \frac{4 f dx}{D}\ (3.95)\]\\
\noindent if $u\neq 0.$, we can divide by $0.5\rho u^2$ which gives\\
\begin{equation} 2\frac{dp}{\rho u^2}+2\frac{\rho u du}{\rho u^2}=-\frac{4 f dx}{D} \label{eqn:mom:a} \end{equation}\\
\noindent using $M^2=u^2/a^2$, $a^2=\gamma p/\rho$ and the chain rule in (\ref{eqn:mom:a}) gives\\
\begin{equation} \frac{2}{\gamma M^2}\frac{dp}{p}+\frac{d(u^2)}{u^2}=-\frac{4 f dx}{D} \label{eqn:mom:b} \end{equation}\\
\noindent From equation (\ref{eqn:gaslaw}) we can get a relation that expresses the pressure derivative term, $dp/p$, in terms of Mach number and $d(u^2)/u^2$. Inserting this in (\ref{eqn:mom:b}) gives\\
\begin{equation} \frac{2}{\gamma M^2}\left\{-\left[\frac{1+(\gamma-1)M^2}{2}\right]\frac{d(u^2)}{u^2}\right\}+\frac{d(u^2)}{u^2}=-\frac{4 f dx}{D} \label{eqn:mom:c} \end{equation}\\
\noindent collecting terms and rearranging gives\\
\begin{equation} \frac{M^2-1}{\gamma M^2}\frac{d(u^2)}{u^2}=\frac{4 f dx}{D} \label{eqn:mom:d} \end{equation}\\
\noindent if we now use equation (\ref{eqn:ttot:Mach}) to get rid of the $d(u^2)/u^2$-term we end up with the following expression\\
\begin{equation*} \frac{4 f dx}{D}=\frac{2}{\gamma M^2}(1-M^2)\left[1+\frac{(\gamma-1)}{2}M^2\right]^{-1}\frac{dM}{M} \end{equation*}
Differential Relations
\noindent In analogy with the heat addition process discussed in the previous section, one-dimensional flow with heat addition is a continuous process. We will derive the differential relations for one-dimensional flow with friction, which will lead to trends for supersonic and supersonic flow with friction.
\begin{figure}[ht!] \begin{center} \includegraphics[]{figures/standalone-figures/Chapter04/pdf/incremental-friction-addition.pdf} \end{center} \caption{Change in flow quantities due to the addition of an infinitesimal pipe segment with the length $dx$} \label{fig:fanno:dx} \end{figure}
\noindent The continuity equation gives
\[d(\rho u)=ud\rho+\rho du \Rightarrow\]
\begin{equation} \dfrac{d\rho}{\rho}=-\dfrac{du}{u} \label{eq:fanno:drho} \end{equation}
\noindent The addition of friction does not affect total temperature and thus the total temperature is constant
\[T_o=T+\dfrac{u^2}{2C_p}=const\]
\noindent differentiating gives
\[dT_o=dT+\dfrac{1}{Cp}udu=0\]
\noindent with $u=M\sqrt{\gamma RT}$, we get
\begin{equation}
\dfrac{dT}{T}=-(\gamma-1)M^2\dfrac{du}{u}
\label{eq:fanno:dT}
\end{equation}\\
\noindent A differential relation for pressure can be obtained from the ideal gas relation
\[p=\rho RT\Rightarrow dp=R(Td\rho+\rho dT)\Rightarrow\]
\begin{equation} \dfrac{dp}{p}=-\left(1+(\gamma-1)M^2\right)\dfrac{du}{u} \label{eq:fanno:dp} \end{equation}\\
\noindent The entropy increase can be obtained from
\[ds=C_v \dfrac{dp}{p}-C_p\dfrac{d\rho}{\rho}\]
\noindent and thus
\begin{equation} ds=-R(1-M^2)\dfrac{du}{u} \label{eq:fanno:ds} \end{equation}\\
\noindent Finally, a relation describing the change in Mach number can be obtained from \\
\[M=\dfrac{u}{\sqrt{\gamma RT}}\Rightarrow dM=M\dfrac{du}{u}-\dfrac{M}{2}\dfrac{dT}{T}\]
\noindent which can be rewritten as
\begin{equation} \dfrac{dM}{M}=\left(1+(\gamma-1)M^2\right)\dfrac{du}{u} \label{eq:fanno:dM} \end{equation}\\
\noindent Eqns.~\ref{eq:fanno:drho} - \ref{eq:fanno:dM} are expressed as functions of $du$ and in order to get a direct relation to the addition of friction caused by the increase in pipe length $dx$, the equations are rewritten so that all variable changes are functions of the entropy increase $ds$.\\
\begin{equation} \dfrac{d\rho}{\rho}=-\dfrac{1}{R(1-M^2)}ds \label{eq:fanno:drho:b} \end{equation}
\begin{equation} \dfrac{dT}{T}=-(\gamma-1)M^2\dfrac{1}{R(1-M^2)}ds \label{eq:fanno:dT:b} \end{equation}
\begin{equation} \dfrac{dp}{p}=-\left(1+(\gamma-1)M^2\right)\dfrac{1}{R(1-M^2)}ds \label{eq:fanno:dp:b} \end{equation}
\begin{equation} \dfrac{dM}{M}=\left(1+\dfrac{(\gamma-1)}{2}M^2\right)\dfrac{1}{R(1-M^2)}ds \label{eq:fanno:dp:b} \end{equation}
\begin{equation} \dfrac{du}{u}=\dfrac{1}{R(1-M^2)}ds \label{eq:fanno:drho:b} \end{equation}\\
\noindent A relation for the change in total pressure can be obtained from
\[ds=C_p\dfrac{dT_o}{T_o}-R\dfrac{dp_o}{p_o}\]\\
\noindent Since total temperature is constant the relation above gives
\begin{equation} \dfrac{dp_o}{p_o}=-\dfrac{ds}{R} \label{eq:fanno:dpo:b} \end{equation}\\
\noindent Using the differential relations above, we can get a good picture of the development of flow variables as friction is continuously added to the flow (see Figure~\ref{fig:fanno:trends}).\\
\begin{figure}[ht!] \begin{subfigure}[b]{0.5\textwidth} \centering \includegraphics[]{figures/standalone-figures/Chapter04/pdf/fanno-subsonic-trends.pdf} \caption{trends (subsonic flow)} \end{subfigure} \begin{subfigure}[b]{0.5\textwidth} \centering \includegraphics[]{figures/standalone-figures/Chapter04/pdf/fanno-supersonic-trends.pdf} \caption{trends (supersonic flow)} \end{subfigure} \caption{Change in flow variables as friction is continuously added to a subsonic flow (left) and supersonic flow (right).} \label{fig:fanno:trends} \end{figure}
Friction Choking
\begin{figure}[ht!] \centering \includegraphics[]{figures/standalone-figures/Chapter04/pdf/fanno-Ts.pdf} \caption{The Fanno flow process illustrated in a $Ts$-diagram. The dashed line shows the critical temperature, $T^\ast$, the blue line is the Fanno flow process, subsonic above $T^\ast$ and supersonic below $T^\ast$. The gray lines are isobars.} \label{fig:friction:Ts} \end{figure}
\noindent Figure~\ref{fig:friction:Ts} shows the Fanno flow process in a $Ts$-diagram. The dashed line represents the sonic temperature, which means that the flow states along the process line above the dashed line are subsonic flow states and the part of the line below the dashed line represents supersonic flow states. In both subsonic and supersonic flow addition of friction leads to a change in temperature in the direction towards the sonic temperature, i.e. the flow approaches sonic conditions ($M=1$). When the length of the pipe through which the fluid flows is equal to the length at which the flow is sonic, the flow is choked (friction choking) and further pipe length cannot be added without a change in the flow conditions. For an initially subsonic flow, a pipe longer than $L^\ast$, the change in flow conditions is analogous to the what happens for addition of heat to a subsonic flow that has reached sonic state discussed in the previous section. The inlet conditions will change such that the massflow is reduced without changing the inlet total conditions such as the pipe length is equal to $L^\ast$ for the new inlet conditions.
\[ \begin{aligned} M_{1'} & = f(L^\ast)\\ T_{1'} & = f(T_o, M_{1'})\\ p_{1'} & = f(p_o, M_{1'})\\ \rho_{1'} & = f(p_{1'}, T_{1'})\\ a_{1'} & = f(T_{1'})\\ u_{1'} & = M_{1'}a_{1'}\\ \end{aligned} \] \\
\begin{figure}[ht!] \begin{subfigure}[b]{0.5\textwidth} \centering \includegraphics[]{figures/standalone-figures/Chapter04/pdf/fanno-Ts-subsonic-choked-mod.pdf} \caption{$Ts$-diagram} \end{subfigure} \begin{subfigure}[b]{0.5\textwidth} \centering \includegraphics[]{figures/standalone-figures/Chapter04/pdf/fanno-Ts-subsonic-choked-mod-close-up.pdf} \caption{$Ts$-diagram (close up)} \end{subfigure} \caption{For friction-choked subsonic flow, further increase of the length of the pipe/duct will lead to an update of the inlet static conditions such that the massflow per unit area is changed and the length corresponding to friction choking $L^\ast$ is increased.} \label{fig:friction:choking:sub} \end{figure}
\begin{figure}[ht!] \centering \includegraphics[]{figures/standalone-figures/Chapter04/pdf/fanno-Ts-supersonic-choked.pdf} \caption{For friction-choked supersonic flow, further increase of the length of the pipe/duct may lead to the generation of a shock inside of the pipe. The location of the shock will be such that $L^\ast$ downstream of the shock equals the remaining length of the pipe at the shock location.} \label{fig:friction:choking:sup} \end{figure}
\noindent For a choked supersonic flow, addition of more friction (increasing the length of the pipe such that $L>L^\ast$) may lead to the generation of a shock inside the pipe. In contrast to the one-dimensional flow with heat addition where a shock does not change $q^\ast$, $L^\ast$ is increased over a shock. The internal shock will be generated in an axial location such that $L^\ast$ downstream of the shock equals the remaining pipe length at the shock location (see Figure~\ref{fig:friction:choking:sup}). As more length is added to the pipe, the shock will move further and further upstream in the pipe until it stands at the pipe entrance. If the pipe is longer than $L^\ast$ after o shock standing at the inlet, the shock will move to the upstream system and the pipe flow will be subsonic and the massflow will be adjusted such that $L=L^\ast$ according to the process described for subsonic choking above.
\noindent From prvevious derivations, we know that $L^\ast$ is a function of mach number according to\\
\[\dfrac{4\bar{f}L^\ast}{D}=\dfrac{1-M^2}{\gamma M^2}+\left(\dfrac{\gamma+1}{2\gamma}\right)\ln\left(\dfrac{(\gamma+1)M^2}{2+(\gamma-1)M^2}\right)\]\\
\noindent by dividing both the numerator and denominator in the fractions by $M^2$ it is easy to see that the choking length (Figure~\ref{fig:friction:factor}) approaches a finite length for great Mach numbers and thus the upper limit for the choking length $L^\ast_1$ is given by\\
\[\left.\dfrac{4\bar{f}L_1^\ast}{D}(M_1)\right|_{M_1\rightarrow \infty}=-\dfrac{1}{\gamma}+\left(\dfrac{\gamma+1}{2\gamma}\right)\ln\left(\dfrac{\gamma+1}{\gamma-1}\right)\]\\
\begin{figure}[ht!] \centering \includegraphics[]{figures/standalone-figures/Chapter04/pdf/fanno-friction-factor.pdf} \caption{$L^\ast$ as a function of Mach number. For high supersonic Mach numbers, the choking length approaches a finite value.} \label{fig:friction:factor} \end{figure}
\noindent From the normal shock relations we know that the downstream Mach number approaches the finite value $\sqrt{(\gamma-1)/2\gamma}$ large Mach numbers and thus the choking length downstream the shock is limited to \\
\[\left.\dfrac{4\bar{f}L_2^\ast}{D}(M_2)\right|_{M_1\rightarrow \infty}=\left(\dfrac{\gamma+1}{\gamma(\gamma-1)}\right)+\left(\dfrac{\gamma+1}{2\gamma}\right)\ln\left(\dfrac{(\gamma+1)(\gamma-1)}{4\gamma+(\gamma-1)^2}\right)\]\\
\noindent From the relations above we get \\
\[\left.\left(\dfrac{4\bar{f}L_2^\ast}{D}(M_2)-\dfrac{4\bar{f}L_1^\ast}{D}(M_1)\right)\right|_{M_1\rightarrow \infty}=\left(\dfrac{2}{\gamma-1}\right)+\left(\dfrac{\gamma+1}{2\gamma}\right)\ln\left[\left(\dfrac{(\gamma+1)(\gamma-1)}{4\gamma+(\gamma-1)^2}\right)\left(\dfrac{\gamma-1}{\gamma+1}\right)\right]\]\\
\noindent Figure~\ref{fig:friction:factor:shock} shows the development of choking length $L_1^\ast$ in a supersonic flow as a function of Mach number in relation to the corresponding choking length $L_2^\ast$ downstream of a normal shock generated at the same Mach number. As can be seen from the figure, a normal shock will always increase the choking length.
\begin{figure}[ht!] \centering \includegraphics[]{figures/standalone-figures/Chapter04/pdf/fanno-friction-factor-shock.pdf} \caption{$L^\ast$ as a function of $M_1$ (Mach number in station 1) for initially supersonic flow. $L_1^\ast$ is the choking length corresponding to $M_1$ and $L_2^\ast$ the choking length downstream of a normal shock at station 1. A normal shock will always increase the choking length.} \label{fig:friction:factor:shock} \end{figure}
