Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

When can one assume Incompressible Flow?

  1. Aug 27, 2006 #1


    User Avatar
    Science Advisor
    Gold Member

    Incompressible Flow: Assumptions for its validity

    After a recent hot discussion brought to the board, I think it would be good to clear up this question.

    Firstly, it does not make sense to talk about an incompressible fluid. There are no incompressible fluids in the Nature, we can only model a incompressible flow under the assumptions I clarify here.

    1.1. The Integral Continuity Equation.

    An incompressible flow field is a field such that the divergence of the velocity is zero [tex]\nabla\cdot \overline{u}=0[/tex]. In other words, the velocity field is solenoidal. For a fluid material volume, that is a volume travelling with the local fluid velocity [tex]\overbar{u}[/tex] and getting deformed as the fluid does flowing, the conservation of mass must hold:

    [tex]\frac{d}{dt}\oint_{V_f(t)}\rho dV=0[/tex], that is, there cannot be any mass variation inside the material volume.

    For doing that time derivative one cannot permute the integral sign and the derivative because the domain of integration depends on time. In order to accomplish the integration one must use the Leibnitz rule, or in this context the so-called Reynolds Transport Teorem (RTT):

    [tex]\frac{d}{dt}\oint_{V_f(t)}\rho dV=\oint_{V_f(t)}\frac{\partial \rho}{\partial t} dV+\oint_{S_f(t)}\overline{u}\cdot \overline{dS}=0[/tex] (1)

    Usually the formulation in terms of the fluid volume [tex]V_f[/tex] is useless, because in real applications is hard to know how an initial fluid volume is distorted. On the other hand, it is more feasible do the formulation in terms of a control volume [tex]V_c(t)[/tex] located randomly (usually one chooses the easiest location) and which boundaries translate with a velocity [tex]\overline{u_c}[/tex]. Let us accomplish the same time derivation than before using the RTT:

    [tex]\frac{d}{dt}\oint_{V_c(t)}\rho dV=\oint_{V_c(t)}\frac{\partial \rho}{\partial t} dV+\oint_{S_c(t)}\overline{u_c}\cdot \overline{dS}[/tex] (2)

    Assume both control and fluid volume coincide in space at some instant [tex]t[/tex]. Then, substracting (2) from (1), one obtains the Continuty Equation:

    [tex]\mathbf{\frac{d}{dt}\oint_{V_c(t)}\rho dV+\oint_{S_c(t)}(\overline{u}-\overline{u_c})\cdot \overline{dS}=0}[/tex] (3)

    Note that for constant density and fixed volume control, the equation reduces to the conservation of the volume flux [tex]\oint_{S_c}\overline{u}\cdot \overline{dS}=0[/tex].

    1.2 The Differential Continuity Equation

    Assume a fixed control volume. Applying the Gauss Theorem to the surface integral of the expression (3), one obtains:

    [tex]\oint_{V}\left\{\frac{\partial\rho}{\partial t}+\nabla \cdot (\rho\overline{u})\right\}dV=0[/tex] (4)

    which is true if:

    [tex]\mathbf{\frac{\partial\rho}{\partial t}+\nabla \cdot (\rho\overline{u})=0}[/tex] (5)

    which corresponds to the differential version of the Continuity Equation. In terms of the Material Derivative, that is, the time derivative viewed from a reference frame attached to the fluid volume:

    [tex]\frac{D}{Dt}=\frac{\partial}{\partial t}+\overline{u}\cdot\nabla[/tex] (6)

    equation (5) can be expressed as:

    [tex]\frac{D\rho}{Dt}+\rho \nabla\cdot \overline{u}=0[/tex] (7)

    and rearranging it in a more suitable form:

    [tex]\mathbf{\nabla\cdot \overline{u}=-\frac{1}{\rho}\frac{D\rho}{Dt}}[/tex] (8)

    Equation (8) is a general expression with only one previous assumption made. I am assuming that the flow is continuum, that is, the spatial derivatives are well defined, what is not the case in highly rarified flows. The measure of rarification is the Knudsen Number [tex]Kn=\frac{\lambda}{L}[/tex], where the mean free path [tex]\lambda\sim 10^{-8}m[/tex] in air at standard conditions. Therefore, the Knudsen Number is low enough in practical applications for considering continuum flow, except those of reentry of space vehicles and hypersonic flight. Also, equation (8) is valid in Reactive Flow. The total production of mass in a fluid volume is still zero under chemical reaction.


    Let us introduce the definition of the diferential of the density. According to the Thermodynamics the differential of [tex]\rho[/tex] can be written as a function of two other thermodynamic variables, let say pressure [tex]P[/tex] and entropy [tex]s[/tex] per unit mass. Then, at first order:

    [tex]d\rho=\frac{\partial \rho}{\partial P}\big)_s dP+\frac{\partial \rho}{\partial s}\big)_P ds[/tex] (9)

    where the subscript means holding constant that variable. Both coefficients [tex]\frac{\partial \rho}{\partial P}\big)_s[/tex] and [tex]\frac{\partial \rho}{\partial s}\big)_P [/tex] are called Thermodynamic Coefficients. Turns out that they correspond to:

    [tex]\frac{\partial \rho}{\partial P}\big)_s =\frac{1}{c^2}[/tex] where [tex]c[/tex] is the speed of sound, and

    [tex]\frac{\partial \rho}{\partial s}\big)_P=\frac{T\rho\beta}{c_p}[/tex], where [tex]\beta[/tex] is the coefficient of volumetric expansion.

    Taking the material derivative of (9) and using the above expressions for the Thermodynamic Coefficients:

    [tex] \frac{D\rho}{Dt}=\frac{1}{c^2}\frac{DP}{Dt}+\frac{T^2\beta}{c_p}\frac{Ds}{Dt}[/tex] (10)

    And from the definition (6) of Material Derivative:

    [tex]\frac{DP}{Dt}=\frac{\partial P}{\partial t}+\overline{u}\cdot\nabla P[/tex] (11)

    Also, the Differential Entropy Conservation equation yields:

    [tex]\frac{Ds}{Dt}=\frac{\phi_v}{\rho T}+K\frac{\nabla^2 T}{\rho T}[/tex] (12)

    (see references). Equation (12) is basically stating that the change of entropy in a fluid volume is caused only by the viscous dissipation [tex]\phi_v[/tex] and the divergence of the heat flux.

    Inserting (11) and (12) in (10), and (10) in (8), one obtains the transcendental equation:

    [tex]\mathbf{\nabla\cdot \overline{u}=-\frac{1}{\rho}\left \{\frac{1}{c^2}\left(\frac{\partial P}{\partial t}+\overline{u}\cdot\nabla P\right)+\frac{T\beta}{c_p}\left(\frac{\phi_v}{T}+K\frac{\nabla^2 T}{T}\right)\right\}}[/tex] (13).

    In the next post I will put an example for working out under what conditions equation (13) can be approximated as [tex]\nabla\cdot \overline{u}\sim 0[/tex] with an error asymptotically small.
    Last edited: Aug 27, 2006
  2. jcsd
  3. Aug 27, 2006 #2


    User Avatar
    Science Advisor
    Gold Member

    Let us remind the equation (13):
    [tex]\mathbf{\nabla\cdot \overline{u}=-\frac{1}{\rho}\left \{\frac{1}{c^2}\left(\frac{\partial P}{\partial t}+\overline{u}\cdot\nabla P\right)+\frac{T^2\beta}{c_p}\left(\frac{\phi_v}{T }+K\frac{\nabla^2 T}{T}\right)\right\}}[/tex]

    Let us imagine a two dimensional channel of length [tex]L[/tex], in vertical position under the action of gravity [tex]g[/tex], and width [tex]D[/tex]. Inside the channel the flow has a characteristic velocity [tex]U[/tex]. The temperature of the left wall is [tex]T_o[/tex] whereas the temperature of the right wall is [tex]T_1[/tex]. In addition to that, the channel is fed by a pump which pulses fluid at a frequency [tex]\omega[/tex]. The fluid has a density [tex]\rho[/tex], dynamic viscosity [tex]\mu[/tex], thermal conductivity [tex]K[/tex], speed of sound [tex]c[/tex] and heat pressure constant [tex]c_p[/tex].

    Do not pay to much attention to the geometrical details. What matters here is the intensity of the nondimensional parameters. All the vital information has been given in the above statement.

    I am going to look at the order of magnitude expected of each term of the right hand side of the equation (13), relative to the order of magnitude of the left hand side:

    [tex]\nabla\cdot\overline{u}\sim O\left(\frac{U}{L}\right)[/tex] (14)

    where [tex]O[/tex] means of the order of.

    2.1. Pressure Unsteadiness Compressibility.

    Look at the first term on the right hand side of (13). The term [tex]\partial P/\partial t[/tex] represents the unsteadiness of the pressure field. This term is usually enabled in Turbomachinery. Turbomachines generate pulses of pressure that translate across the fluid generating pulses of density. One should expect that depending on how large are the pulses of pressure the assumption of incompressibility of the velocity field loses accuracy. What matters here is the value of the pressure generated and the frequency of oscillation:

    [tex]\frac{\frac{1}{\rho c^2}\frac{\partial P}{\partial t}}{\nabla\cdot\overline{u}}\sim \frac{O\left(\frac{\omega \Delta P}{\rho c^2}\right)}{O\left(\frac{U}{L}\right)}=St\cdot M^2[/tex] (15)

    where [tex]St=\omega L/U[/tex] is the Strouhal Number, and [tex]M=U^2/c^2[/tex] is the Mach Number. The first condition for incompressibilty is then:

    [tex]\mathbf{St\cdot M^2<<1}[/tex] (16)

    That is, for the Strouhal Number being small the characteristic time of pulsation [tex]t_o\sim 1/\omega[/tex] must be much longer than the characteristic time than a fluid particle resides in the pipe [tex]t_r\sim L/U[/tex]. In addition to that, low velocities compared with the speed of sound also make the Mach Number small collaborating to the accuracy of the Incompressible flow assumption.

    2.2 Compressibility of Large Scale Flows.

    Look at the second term of the right hand side of the equation (13). The term [tex]\overline{u}\cdot \nabla P[/tex] represents the transport of the pressure by the inertia of the fluid. This transport is important in Ocean and Atmospheric Dynamics, where large scales are involved:

    [tex]\frac{\frac{\overline{u}\cdot \nabla P}{\rho c^2}}{\nabla\cdot\overline{u}}\sim \frac{O\left(\frac{g U}{c^2}\right)}{O(U/L)}\sim \frac{gL}{c^2}=\frac{M^2}{Fr^2}[/tex] (17)

    where [tex]Fr=U^2/\sqrt{gL}[/tex] is the Froude Number. Thus, for the accuracy of the incompressibility assumption, [tex] M^2/Fr^2<<1[/tex], or what is the same:

    [tex]\mathbf{\frac{gL}{c^2}<<1}[/tex] (18)

    Note that the larger the system, the lesser the accuracy of the incompressibility assumption. The length [tex]c^2/g[/tex] is called the scale heigth, and it is the characteristic scale in the atmosphere over which there are appreciable variations of pressure. Thus for scales larger than the scale heigth the gradient of pressure causes hydrostatically large gradients of density. In order of (18) to be of order unity, [tex]L\sim c^2/g[/tex]. Thus, the incompressibility assumption begins to fail on systems of scale of 100 km in the Ocean and 1 km in the atmosphere. The incompressibility assumption in the Ocean works better than in the Atmosphere.

    2.3. Compressibility caused by the Coefficient of Volumetric Expansion.

    Let us take the third term of the right hand side of the equation (13). The viscous dissipation function [tex]\phi_v[/tex] is proportional to the contracted product of the strain tensor with itself.

    [tex]\frac{\frac{\phi_v T\beta}{c_p\rho}}{\nabla\cdot\overline{u}}\sim \frac{O\left(\frac{\mu U^2}{L^2}\frac{\beta}{c_p\rho}\right)}{O(U/L)}\sim Re^{-1}\frac{\beta U^2}{c_p}[/tex] (19)

    where [tex]Re=\rho UL/\mu[/tex] is the Reynolds Number. Therefore, for the accuracy of the incompressibility assumption:

    [tex] \mathbf{Re^{-1}\frac{\beta U^2}{c_p}<<1}[/tex] (20)

    which is accomplished for low volumetric expansion coefficients.

    2.4 Compressibility caused by Temperature Gradient.

    The last term on the right hand side of equation (13) represents an irreversible source of entropy due to heat transfer.

    [tex]\frac{T\beta K\frac{\nabla^2 T}{c_p T}}{\nabla\cdot\overline{u}}\sim \frac{O\left(\frac{\beta K \Delta T}{c_p \rho L^2}\right)}{O(U/L)}\sim Re^{-1}Pr^{-1}\beta \Delta T[/tex] (21)

    Where [tex]Pr=\mu c_p/K[/tex] is the Prandtl Number. For air [tex]Pr\sim 0.7[/tex] whereas for water [tex]\Pr\sim 7[/tex], which means that air is a better heat conductor than water. For the accuracy of the Incompressibility assumption:

    [tex]\mathbf{Re^{-1}Pr^{-1}\beta \Delta T<<1}[/tex] (22)

    That is, the flow must be fast, the substance should not transfer heat very well for having a high Prandtl Number, and the temperature differences must not be very large. The critical temperature difference should be

    [tex] \Delta T\sim Re\cdot Pr /\beta[/tex] (23)


    Assume for the problem stated above that the substance is water, with [tex]\rho=1000kg/m^3[/tex], [tex]\mu\sim 10^{-3}Kg/ms[/tex], and
    [tex]\beta\sim 10^{-8} K^{-1}[/tex], [tex]c_p\sim 1000J/KgK[/tex], [tex]c\sim 10^3 m/s[/tex] and an usual flow velocity of [tex]U\sim 10 m/s[/tex]. And leave the length of the channel, the difference of temperature and the frequency of pulsation free. The criterions for incompressible flow are (23), (20), (18) and (16):

    I) [tex] St \cdot M^2=\frac{\omega L}{10}\frac{10^2}{1000^2}=\omega L 10^3<<1[/tex]

    II) [tex] \frac{gL}{c^2}=\frac{10L}{1000^2}=\frac{L}{100000}<<1[/tex]

    This gives a critical length [tex]L=100km[/tex] where the incompressibility assumption begins to fail. Thus, the critical frequency for that length is [tex]\omega\sim 10^{-8}Hz[/tex]. Frequencies larger would cause I) to be of order 1 and the incompressibility assumption fails. For ordinary lengths like [tex]L\sim 1m[/tex], the critical frequency is [tex]\omega \sim10KHz[/tex]. At that frequency the time required for a pressure wave to travel from one end to the other is of the same order than the time of residence of a fluid particle inside the channel. Let us go on:

    III) [tex] Re^{-1}\frac{\beta U^2}{c_p}\sim 10^{-11}<<1[/tex] for ordinary lengths (1m).

    IV) [tex]Re^{-1}Pr^{-1}\beta \Delta T\sim \Delta T \cdot 10^{-14}[/tex] for ordinary lengths. For being of order unity one needs a [tex]\Delta T[/tex] enormous, what is physically inviable in water.


    According to the above analysis, the right hand side of the equation (13) becomes vanishingly small compared with the left hand side under the analyzed conditions.

    Let's call
    [tex]\Lambda=St \cdot M^2[/tex]


    [tex]\Gamma=Re^{-1}\frac{\beta U^2}{c_p}[/tex]

    and [tex]\Theta=Re^{-1}Pr^{-1}\beta \Delta T[/tex]

    being [tex]\Lambda, \Gamma, \Pi, \Theta<<1[/tex]. Let us make an asymptotic expansion of the velocity field in terms of these small parameters:

    [tex]\frac{\overline{u}}{U}=\overline{u}_o+(\Lambda+\Gamma+\Pi+\Theta)\overline{u}_1+ O(\Lambda^2,\Gamma^2,\Pi^2,\Theta^2)[/tex]

    and substitute this expansion into (13) letting all these four nondimensional parameters tend to zero, then at first order (i.e. with errors of those non dimensional parameters), the Continuity Equation reads:

    [tex]\mathbf{\nabla \cdot \overline{u}_o=0+ O(\Lambda+\Gamma+\Pi+\Theta)}[/tex]

    being the second term on the right vanishingly small under the assumptions made above.

    This is what is called the Boussinesq Approximation.

    Finally, to end with this tutorial, you have to read my signature:

    However, in Science, understanding is achieved through simplification. LiƱan A. and Williams F.A. in "Fundamentals Aspects of Combustion".
    Last edited: Aug 27, 2006
  4. Aug 27, 2006 #3


    User Avatar
    Staff Emeritus
    Science Advisor

    Nice job :approve:
  5. Aug 27, 2006 #4
    Wow...I just assumed Ma < 0.3 usually, but I suppose undergrads can get away with it. :biggrin: Very impressive.
  6. Aug 30, 2006 #5
    I have never taken a fluid dynamics class, but this tutorial is magnificent. Clear, concise and with order of magnitude arguments even someone without a fluids background was able to understand your arguments.

    Hats off to you good sir!
  7. Sep 5, 2006 #6
    wow-much efforts-
    good work:approve:
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook