Problem 35 - Thomas Calculus

Question

Conservation of mass Let \(\mathbf{v}(t, x, y, z)\) be a continuously differentiable vector field over the region \(D\) in space and let \(p(t, x)\) \(y, z )\) be a continuously differentiable scalar function. The variable \(t\) represents the time domain. The Law of Conservation of Mass asserts that $$\frac{d}{d t} \iiint_{D} p(t, x, y, z) d V=-\iint_{S} p \mathbf{v} \cdot \mathbf{n} d \sigma$$ where \(S\) is the surface enclosing \(D\) a. Give a physical interpretation of the conservation of mass law if \(\mathbf{v}\) is a velocity flow field and \(p\) represents the density of the fluid at point \((x, y, z)\) at time \(t\) b. Use the Divergence Theorem and Leibniz's Rule, $$\frac{d}{d t} \iiint_{D} p(t, x, y, z) d V=\iiint_{D} \frac{\partial p}{\partial t} d V$$ to show that the Law of Conservation of Mass is equivalent to the continuity equation, $$\nabla \cdot p \mathbf{v}+\frac{\partial p}{\partial t}=0$$ (In the first term \(\nabla \cdot p \mathbf{v},\) the variable \(t\) is held fixed, and in the second term \(\partial p / \partial t,\) it is assumed that the point \((x, y, z)\) in \(D\) is held fixed.)

Step-by-Step Solution

Answer

The conservation of mass states mass change in a region equals net flow through its boundary, leading to the continuity equation \(\nabla \cdot (p \mathbf{v}) + \frac{\partial p}{\partial t} = 0\).

Step 1: Physical Interpretation

The conservation of mass law states that the rate of change of mass within a volume \(D\) is equal to the net flux of mass across the boundary \(S\) of that volume. Physically, if \(\mathbf{v}\) represents the flow velocity of the fluid, and \(p\) represents the density at a point \((x, y, z)\) and time \(t\), this law indicates that any increase or decrease in mass inside the region \(D\) is accounted for by mass entering or leaving through the surface \(S\).

Step 2: Applying the Divergence Theorem

The Divergence Theorem allows the transformation from a surface integral to a volume integral: \[ \iint_{S} p \mathbf{v} \cdot \mathbf{n} \, d \sigma = \iiint_{D} abla \cdot (p \mathbf{v}) \, dV \]This means that the total flow of density through the surface \(S\) is the integral of the divergence over the volume \(D\).

Step 3: Using Leibniz's Rule

By using Leibniz's Rule for differentiating under the integral sign, we interpret the change in integral of \(p\) over \(D\) as:\[ \frac{d}{d t} \iiint_{D} p(t, x, y, z) \, dV = \iiint_{D} \frac{\partial p}{\partial t} \, dV \]This states that the rate of change of mass in \(D\) is the integral of the change in density over time within the volume.

Step 4: Equating Divergence and Time Derivative

Now, setting the results from Steps 2 and 3 in the conservation equation:\[ \iiint_{D} \frac{\partial p}{\partial t} \, dV = -\iiint_{D} abla \cdot (p \mathbf{v}) \, dV \]We can combine them into a single integral, considering that it holds true for all configurations of \(D\):\[ \iiint_{D} \left( \frac{\partial p}{\partial t} + abla \cdot (p \mathbf{v}) \right) \, dV = 0 \]

Step 5: Concluding with the Continuity Equation

From the integral being zero over any arbitrary region \(D\), it follows directly that:\[ \frac{\partial p}{\partial t} + abla \cdot (p \mathbf{v}) = 0 \]This is called the continuity equation, which expresses the local conservation of mass. The terms indicate how the density \(p\) changes with time due to divergence in the velocity field \(\mathbf{v}\).

Key Concepts

Continuity Equation

The continuity equation is an essential mathematical representation that stems from the principle of conservation of mass. When we study a flowing fluid, such as water or air, we need to ascertain that no mass is lost or gained in a region unless it is flowing in or out. The continuity equation states that the change in mass within a control volume is balanced by the mass flowing across its boundaries.
This can be mathematically expressed as:

• \( \frac{\partial p}{\partial t} + abla \cdot (p \mathbf{v}) = 0 \)

Here, \(p\) denotes the density, and \(\mathbf{v}\) represents the fluid's velocity. This equation ensures that as time progresses, the relationship between density changes and the flow of material adheres to mass conservation.
If we simplify, this equation determines when and where fluid will compress or expand in response to changes in velocity or density, giving engineers and scientists critical insights into fluid behavior.

Divergence Theorem

The Divergence Theorem acts as a bridge between the flux of a vector field across a surface and the behavior of the field inside a volume. When dealing with the conservation of mass, this theorem enables the transformation of a surface integral—describing the mass flux through a closed surface—into a volume integral which accounts for the divergence of the field within the volume.
Mathematically, the theorem is expressed as:

• \( \iint_{S} \mathbf{F} \cdot \mathbf{n} \, d \sigma = \iiint_{D} abla \cdot \mathbf{F} \, dV \)

In this equation, \(S\) stands for the surface enclosing the volume \(D\), \(\mathbf{F}\) represents a vector field (such as \(p \mathbf{v}\) in our context), and \(\mathbf{n}\) is the outward normal vector.
This powerful theorem simplifies calculations by converting surface-based properties to volume-based ones, helping to analyze changes and predict outcomes effectively.

Leibniz's Rule

Leibniz's Rule is vital for differentiating an integral when its limits are constants and its integrand is a function of time. This is particularly crucial in dynamic settings like fluid mechanics, where the properties change over time.
The Leibniz's Rule for differentiating under the integral is given by:

• \( \frac{d}{dt} \iiint_{D} p(t, x, y, z) \, dV = \iiint_{D} \frac{\partial p}{\partial t} \, dV \)

This formula is beneficial when you want to focus on how a quantity distributed over a region changes as time goes by, without altering the spatial boundaries.
Applying Leibniz's Rule simplifies the understanding of how the internal properties of a mass change temporally, distinguishing time factors from spatial variables efficiently.

Density

Density is a scalar quantity that reflects the mass per unit volume at a given point in space and time. In the conservation of mass context, density changes signal how mass is distributed across a flow field.
For instance, in the equation:

• \( \frac{\partial p}{\partial t} \)

\(p\) symbolizes the fluid's density, and the term tells us how density varies with time at a consistent spatial point.
A higher density often indicates more mass in a smaller area, potentially implying higher pressure or reduced velocity depending on constraints.
Understanding density's role is critical, as it affects how fluids behave under various forces and environmental conditions.

Vector Field

A vector field assigns a vector to every point in space, representing both magnitude and direction. In fluid dynamics, such as in the study of mass conservation, these vectors often illustrate flow velocity at each point.
Considering the vector field in our density-flow relationship:

• \( \mathbf{v} = \text{{velocity of the fluid}} \)

The product \(p \mathbf{v} \) is then a vector field that reflects how the mass (density times velocity) flows through each space point.
Vector fields help visualize complex fluid behaviors and simplify the understanding of how fluids interact with their environment, making them indispensable tools in modeling physical systems.