Talk:PlanetPhysics/Continuity Equation

Original TeX Content from PlanetPhysics Archive

%%% This file is part of PlanetPhysics snapshot of 2011-09-01 %%% Primary Title: continuity equation %%% Primary Category Code: 47.10.+g %%% Filename: ContinuityEquation.tex %%% Version: 1 %%% Owner: mdo %%% Author(s): mdo %%% PlanetPhysics is released under the GNU Free Documentation License. %%% You should have received a file called fdl.txt along with this file. %%% If not, please write to gnu@gnu.org. \documentclass[12pt]{article} \pagestyle{empty} \setlength{\paperwidth}{8.5in} \setlength{\paperheight}{11in}

\setlength{\topmargin}{0.00in} \setlength{\headsep}{0.00in} \setlength{\headheight}{0.00in} \setlength{\evensidemargin}{0.00in} \setlength{\oddsidemargin}{0.00in} \setlength{\textwidth}{6.5in} \setlength{\textheight}{9.00in} \setlength{\voffset}{0.00in} \setlength{\hoffset}{0.00in} \setlength{\marginparwidth}{0.00in} \setlength{\marginparsep}{0.00in} \setlength{\parindent}{0.00in} \setlength{\parskip}{0.15in}

\usepackage{html}

% this is the default PlanetPhysics preamble. as your knowledge % of TeX increases, you will probably want to edit this, but % it should be fine as is for beginners.

% almost certainly you want these \usepackage{amssymb} \usepackage{amsmath} \usepackage{amsfonts}

% used for TeXing text within eps files %\usepackage{psfrag} % need this for including graphics (\includegraphics) %\usepackage{graphicx} % for neatly defining theorems and propositions %\usepackage{amsthm} % making logically defined graphics %\usepackage{xypic}

% there are many more packages, add them here as you need them

% define commands here

\begin{document}

Much like the \htmladdnormallink{mechanics}{http://planetphysics.us/encyclopedia/Mechanics.html} of \htmladdnormallink{solid}{http://planetphysics.us/encyclopedia/CoIntersections.html} bodies, the \htmladdnormallink{concept}{http://planetphysics.us/encyclopedia/PreciseIdea.html} of \htmladdnormallink{mass}{http://planetphysics.us/encyclopedia/CosmologicalConstant.html} conservation plays an important role in the study of fluids. In conceptual terms, mass conservation for a fluid states that for a given \htmladdnormallink{volume}{http://planetphysics.us/encyclopedia/Volume.html} in fixed space, the rate of change of fluid mass contained within that volume is equivalent to the difference between the rate of change of mass into and out of that volume. A typical way to conceptualise this is to consider the flow of fluid in a pipe. At one end of the pipe, water is flowing in at some constant rate, and at the other end, water pours out at the same rate. What does this mean if the pipe is of variable width? It means that water must flow at a greater rate through the narrow \htmladdnormallink{sections}{http://planetphysics.us/encyclopedia/IsomorphicObjectsUnderAnIsomorphism.html} of pipe than through the broader sections. For a given unit of \htmladdnormallink{velocity}{http://planetphysics.us/encyclopedia/Velocity.html}, a greater amount of mass will flow through the wide section than through the narrow section. Thus, to conserve mass, the velocity of the water must increase as the radius of the pipe decreases.

The mass of a fluid can be written mathematically in terms of the fluid density and volume. Suppose that the containing volume is a sphere of radius $r$, and the fluid is of constant density. Then, the mass of fluid, $M$, is \textit{density of fluid}$\times$\textit{volume of fluid}. In the case presented in this example, $M=\rho\times(4\pi r^{3}/3)$, where $\rho$ is used to denote density. In general, we cannot assume constant density, nor can we assume ideal geometric forms for the containing volume. Thus, we write mass as \begin{equation} M = \int_{V}\rho dV, \end{equation} with $V$ the volume. Mass conservation then takes the form \begin{equation} \frac{\partial}{\partial t}\int_{V}\rho dV = -\int_{S}\rho\mathbf{u}\cdot d\mathbf{S}. \end{equation} The left hand side is simply a mathematical \htmladdnormallink{representation}{http://planetphysics.us/encyclopedia/CategoricalGroupRepresentation.html} of the statement \textit{the rate of change of mass in volume $V$}. The term $\rho\mathbf{u}$ is the mass \htmladdnormallink{flux}{http://planetphysics.us/encyclopedia/AbsoluteMagnitude.html} across some area on the surface of the volume, $V$, where $\mathbf{u}$ is the velocity \htmladdnormallink{vector}{http://planetphysics.us/encyclopedia/Vectors.html}. Thus, the surface integral of the right hand side gives the total mass exiting the volume $V$ in a unit time. To bring this equation into its common state, we must apply the \htmladdnormallink{divergence}{http://planetphysics.us/encyclopedia/DivergenceOfAVectorField.html} \htmladdnormallink{theorem}{http://planetphysics.us/encyclopedia/Formula.html}. This allows us to convert the flux through a surface term on the right hand side into a divergence within a volume. In particular, \begin{equation} \int_{S}\rho\mathbf{u}\cdot d\mathbf{S} = \int_{V}\nabla\cdot(\rho\mathbf{u})dV. \end{equation} Substituting this result into the original statement of mass conservation, and moving all terms to the left hand side, we have \begin{equation} \frac{\partial}{\partial t}\int_{V}\rho dV + \int_{V}\nabla\cdot(\rho\mathbf{u})dV = \int_{V}\left( \frac{\partial\rho}{\partial t} + \nabla\cdot(\rho\mathbf{u})\right) dV = 0. \end{equation} The time derivative was taken inside the integral because the setting here is an Eularian one \emph{i.e.} the volume is fixed in space. This equation abstracts directly to the familiar form of mass conservation for a fluid. Since we initially placed no restrictions or constraints on the form the volume $V$ takes, it must be true that the integral holds over all volumes $V$. Thus, it follows that \begin{equation}\label{continuity_eularian} \frac{\partial\rho}{\partial t} + \nabla\cdot(\rho\mathbf{u}) = 0. \end{equation} This is the expression of mass conservation in the Eularian setting. The equation can be equivalently formulated in the \htmladdnormallink{Lagrangian}{http://planetphysics.us/encyclopedia/LagrangesEquations.html} setting, and is written as \begin{equation}\label{continuity_lagrangian} \frac{D\rho}{Dt} + \rho\nabla\cdot\mathbf{u} = 0. \end{equation}

These expressions for the conservation of mass are typically referred to as continuity equations. The continuity equation finds application in nearly all areas of fluid mechanics. For instance, in physical oceanography, upwelling is an area of great interest. Upwelling is the upward \htmladdnormallink{motion}{http://planetphysics.us/encyclopedia/CosmologicalConstant.html} of deep water through the water column. Upwelling will, for instance, bring cooler bottom water to the surface, or circulate nutrient rich and nutrient poor waters. One mechanism through which upwelling can occur is mass conservation. Suppose that the horizontal flow of water converges at some point along the ocean floor. At that point, the water must move upward, in order to conserve mass, thus giving a region of upwelling.

As a final point, it is important that we remark upon the form of the continuity equation for an incompressible fluid. A fluid is considered incompressible when changes in pressure do not result in changes in density. For many practical applications, water is considered incompressible, for example. By enforcing the incompressibility restriction, it follows that the rate of change of density following the motion (the Lagrangian setting) must be zero \emph{i.e} the density following the motion of the fluid parcel is constant. Thus, Equation \ref{continuity_lagrangian} reduces to \begin{equation} \nabla\cdot\mathbf{u} = 0. \end{equation}

This states that for mass to be conserved in an incompressible fluid, the velocity \htmladdnormallink{field}{http://planetphysics.us/encyclopedia/CosmologicalConstant.html} must be divergence free.

\end{document}