Movatterモバイル変換


[0]ホーム

URL:


Jump to content
WikipediaThe Free Encyclopedia
Search

Incompressible flow

From Wikipedia, the free encyclopedia
Fluid flow in which density remains constant

icon
This articleneeds additional citations forverification. Please helpimprove this article byadding citations to reliable sources. Unsourced material may be challenged and removed.
Find sources: "Incompressible flow" – news ·newspapers ·books ·scholar ·JSTOR
(December 2019) (Learn how and when to remove this message)

Influid mechanics, or more generallycontinuum mechanics,incompressible flow is aflow in which the materialdensity does not vary over time. Equivalently, thedivergence of an incompressible flow velocity is zero. Under certain conditions, the flow ofcompressible fluids can be modelled as incompressible flow to a good approximation.

Derivation

[edit]

The fundamental requirement for incompressible flow is that the density,ρ{\displaystyle \rho }, is constant within a small element volume,dV, which moves at the flow velocityu. Mathematically, this constraint implies that thematerial derivative (discussed below) of the density must vanish to ensure incompressible flow. Before introducing this constraint, we must apply theconservation of mass to generate the necessary relations. The mass is calculated by avolume integral of the density,ρ{\displaystyle \rho }:

m=VρdV.{\displaystyle {m}={\iiint \limits _{V}\!\rho \,\mathrm {d} V}.}

The conservation of mass requires that thetime derivative of the mass inside acontrol volume be equal to themass flux,J, across its boundaries. Mathematically, we can represent this constraint in terms of asurface integral:

mt=SJdS{\displaystyle {\partial m \over \partial t}=-\oint _{S}\mathbf {J} \cdot d\mathbf {S} }

The negative sign in the above expression ensures that outward flow results in a decrease in the mass with respect to time, using the convention that the surface area vector points outward. Now, using thedivergence theorem we can derive the relationship between the flux and the partial time derivative of the density:

VρtdV=V(J)dV,{\displaystyle {\iiint \limits _{V}{\partial \rho \over \partial t}\,\mathrm {d} V}={-\iiint \limits _{V}\left(\nabla \cdot \mathbf {J} \right)\,\mathrm {d} V},}

therefore:

ρt=J.{\displaystyle {\partial \rho \over \partial t}=-\nabla \cdot \mathbf {J} .}

Thepartial derivative of the density with respect to time need not vanish to ensure incompressibleflow. When we speak of the partial derivative of the density with respect to time, we refer to this rate of change within a control volume offixed position. By letting the partial time derivative of the density be non-zero, we are not restricting ourselves to incompressiblefluids, because the density can change as observed from a fixed position as fluid flows through the control volume. This approach maintains generality, and not requiring that the partial time derivative of the density vanish illustrates that incompressible fluids can still undergo compressible flow. What interests us is the change in density of a control volume that moves along with the flow velocity,u. The flux is related to the flow velocity through the following function:

J=ρu.{\displaystyle {\mathbf {J} }={\rho \mathbf {u} }.}

So that the conservation of mass implies that:

ρt+(ρu)=ρt+ρu+ρ(u)=0.{\displaystyle {\partial \rho \over \partial t}+{\nabla \cdot \left(\rho \mathbf {u} \right)}={\partial \rho \over \partial t}+{\nabla \rho \cdot \mathbf {u} }+{\rho \left(\nabla \cdot \mathbf {u} \right)}=0.}

The previous relation (where we have used the appropriateproduct rule) is known as thecontinuity equation. Now, we need the following relation about thetotal derivative of the density (where we apply thechain rule):

dρdt=ρt+ρxdxdt+ρydydt+ρzdzdt.{\displaystyle {\mathrm {d} \rho \over \mathrm {d} t}={\partial \rho \over \partial t}+{\partial \rho \over \partial x}{\mathrm {d} x \over \mathrm {d} t}+{\partial \rho \over \partial y}{\mathrm {d} y \over \mathrm {d} t}+{\partial \rho \over \partial z}{\mathrm {d} z \over \mathrm {d} t}.}

So if we choose a control volume that is moving at the same rate as the fluid (i.e. (dx/dtdy/dtdz/dt) = u), then this expression simplifies to thematerial derivative:

DρDt=ρt+ρu.{\displaystyle {D\rho \over Dt}={\partial \rho \over \partial t}+{\nabla \rho \cdot \mathbf {u} }.}

And so using the continuity equation derived above, we see that:

DρDt=ρ(u).{\displaystyle {D\rho \over Dt}={-\rho \left(\nabla \cdot \mathbf {u} \right)}.}

A change in the density over time would imply that the fluid had either compressed or expanded (or that the mass contained in our constant volume,dV, had changed), which we have prohibited. We must then require that the material derivative of the density vanishes, and equivalently (for non-zero density) so must the divergence of the flow velocity:

u=0.{\displaystyle {\nabla \cdot \mathbf {u} }=0.}

And so beginning with the conservation of mass and the constraint that the density within a moving volume of fluid remains constant, it has been shown that an equivalent condition required for incompressible flow is that the divergence of the flow velocity vanishes.

Relation to compressibility

[edit]

In some fields, a measure of the incompressibility of a flow is the change in density as a result of the pressure variations. This is best expressed in terms of thecompressibility

β=1ρdρdp.{\displaystyle \beta ={\frac {1}{\rho }}{\frac {\mathrm {d} \rho }{\mathrm {d} p}}.}

If the compressibility is acceptably small, the flow is considered incompressible.

Relation to solenoidal field

[edit]

An incompressible flow is described by asolenoidal flow velocity field. But a solenoidal field, besides having a zerodivergence, also has the additional connotation of having non-zerocurl (i.e., rotational component).

Otherwise, if an incompressible flow also has a curl of zero, so that it is alsoirrotational, then the flow velocity field is actuallyLaplacian.

Difference from material

[edit]

As defined earlier, an incompressible (isochoric) flow is the one in which

u=0.{\displaystyle \nabla \cdot \mathbf {u} =0.\,}

This is equivalent to saying that

DρDt=ρt+uρ=0{\displaystyle {\frac {D\rho }{Dt}}={\frac {\partial \rho }{\partial t}}+\mathbf {u} \cdot \nabla \rho =0}

i.e. thematerial derivative of the density is zero. Thus if one follows a material element, its mass density remains constant. Note that the material derivative consists of two terms. The first termρt{\displaystyle {\tfrac {\partial \rho }{\partial t}}} describes how the density of the material element changes with time. This term is also known as theunsteady term. The second term,uρ{\displaystyle \mathbf {u} \cdot \nabla \rho } describes the changes in the density as the material element moves from one point to another. This is theadvection term (convection term for scalar field). For a flow to be accounted as bearing incompressibility, the accretion sum of these terms should vanish.

On the other hand, ahomogeneous, incompressible material is one that has constant density throughout. For such a material,ρ=constant{\displaystyle \rho ={\text{constant}}}. This implies that,

ρt=0{\displaystyle {\frac {\partial \rho }{\partial t}}=0} and
ρ=0{\displaystyle \nabla \rho =0}independently.

From the continuity equation it follows that

DρDt=ρt+uρ=0  u=0{\displaystyle {\frac {D\rho }{Dt}}={\frac {\partial \rho }{\partial t}}+\mathbf {u} \cdot \nabla \rho =0\ \Rightarrow \ \nabla \cdot \mathbf {u} =0}

Thus homogeneous materials always undergo flow that is incompressible, but the converse is not true. That is, compressible materials might not experience compression in the flow.

Related flow constraints

[edit]

In fluid dynamics, a flow is considered incompressible if the divergence of the flow velocity is zero. However, related formulations can sometimes be used, depending on the flow system being modelled. Some versions are described below:

  1. Incompressible flow:u=0{\displaystyle {\nabla \cdot \mathbf {u} =0}}. This can assume either constant density (strict incompressible) or varying density flow. The varying density set accepts solutions involving small perturbations indensity, pressure and/or temperature fields, and can allow for pressurestratification in the domain.
  2. Anelastic flow:(ρou)=0{\displaystyle {\nabla \cdot \left(\rho _{o}\mathbf {u} \right)=0}}. Principally used in the field ofatmospheric sciences, the anelastic constraint extends incompressible flow validity to stratified density and/or temperature as well as pressure. This allows the thermodynamic variables to relax to an 'atmospheric' base state seen in the lower atmosphere when used in the field of meteorology, for example. This condition can also be used for various astrophysical systems.[1]
  3. Low Mach-number flow, orpseudo-incompressibility:(αu)=β{\displaystyle \nabla \cdot \left(\alpha \mathbf {u} \right)=\beta }. The lowMach-number constraint can be derived from the compressible Euler equations using scale analysis of non-dimensional quantities. The restraint, like the previous in this section, allows for the removal of acoustic waves, but also allows forlarge perturbations in density and/or temperature. The assumption is that the flow remains within a Mach number limit (normally less than 0.3) for any solution using such a constraint to be valid. Again, in accordance with all incompressible flows the pressure deviation must be small in comparison to the pressure base state.[2]

These methods make differing assumptions about the flow, but all take into account the general form of the constraint(αu)=β{\displaystyle \nabla \cdot \left(\alpha \mathbf {u} \right)=\beta } for general flow dependent functionsα{\displaystyle \alpha } andβ{\displaystyle \beta }.

Numerical approximations

[edit]

The stringent nature of incompressible flow equations means that specific mathematical techniques have been devised to solve them. Some of these methods include:

  1. Theprojection method (both approximate and exact)
  2. Artificial compressibility technique (approximate)
  3. Compressibility pre-conditioning

See also

[edit]

References

[edit]
  1. ^Durran, D.R. (1989)."Improving the Anelastic Approximation"(PDF).Journal of the Atmospheric Sciences.46 (11):1453–1461.Bibcode:1989JAtS...46.1453D.doi:10.1175/1520-0469(1989)046<1453:ITAA>2.0.CO;2.ISSN 1520-0469.[dead link]
  2. ^Almgren, A.S.; Bell, J.B.; Rendleman, C.A.; Zingale, M. (2006)."Low Mach Number Modeling of Type Ia Supernovae. I. Hydrodynamics"(PDF).Astrophysical Journal.637 (2):922–936.arXiv:astro-ph/0509892.Bibcode:2006ApJ...637..922A.doi:10.1086/498426. Archived fromthe original(PDF) on 2008-10-31. Retrieved2008-12-04.
Retrieved from "https://en.wikipedia.org/w/index.php?title=Incompressible_flow&oldid=1313025019"
Category:
Hidden categories:

[8]ページ先頭

©2009-2026 Movatter.jp