Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

Mass conservation is usually written as

$$\frac{\partial \rho}{\partial t} + \operatorname{div}(\rho \boldsymbol v) = 0$$

$\rho$ is the density and $\boldsymbol v$ is the fluid velocity. My attempt to rewrite it in notions used in differential geometry:

$$\frac{\partial \rho}{\partial t} + \star\, \text{d} \star (\rho v^\flat) = 0$$

But I have a doubt that it may have another form, probably like in classical mechanics Liouville equation.

Maybe there is a book or something about fluid dynamics on manifolds? (I've asked this last question on Physics.SE, but if somebody knows it here let me know).

share|improve this question
1  
The following (non-English) presentation gives answers to this question on the last two slides. –  Thomas Klimpel Aug 3 '12 at 8:10
add comment

1 Answer 1

up vote 2 down vote accepted

Let $M$ be a smooth manifold with a given volume form $\mu$ on it.
Let $\mathbf v_t$ be the time-dependent vector-field on $M$ describing the fluid velocity, and $\rho_t$ the time-dependent function on $M$ describing the fluid density.

The continuity equation $\partial_t\rho_t+\text{div}(\rho_t\mathbf v_t)$ states the invariance of the volume form $\rho_t\mu$ under $\text{Fl}^{\mathbf v}_{t,s},$ the time-dependent flow of $\mathbf{v},$ i.e., explicitly $$\frac{d}{dt}(\text{Fl}^{\mathbf v}_{t,t_0})^\ast(\rho_t\mu)=0.\tag{*}$$

This interpretation of the continuity equation is based on:

  1. Lie derivative theorem, i.e. $\frac{d}{dt}(\text{Fl}^{\mathbf v}_{t,t_0})^\ast T_t|_p=\left.\left[(\text{Fl}^{\mathbf v}_{t,t_0})^\ast\left(\mathcal{L}_{\mathbf v_t}T_t+\partial_t T_t\right)\right]\right|_p,$ where $T_t$ is an arbitrary time-dependent tensor-field on $M.$
  2. the definition of divergence of a vector-field w.r.t. the volume form $\mu,$ by which $\mathcal{L}_{\mathbf v_t}(\rho_t\mu)=d(\mathbf v_t\rfloor\rho_t\mu)=\text{div}(\rho_t\mathbf v_t)\mu.$

Infact, putting 1. and 2. into (*), we get $(\text{Fl}^{\mathbf v}_{t,t_0})^\ast\left[\left(\text{div}(\rho_t\mathbf v_t)+\partial_t \rho_t\right)\mu\right]=0,$ i.e. $$\text{div}(\rho_t\mathbf v_t)+\partial_t \rho_t=0.$$

Edit Note that the (*) can be rewritten in integral form as the law of conservation of mass: $$\tag{**}\frac{d}{dt}\int_{\text{Fl}_{t,t_0}^{\mathbf v}D}\rho_t\mu=0,\text{ for all compact }D\subseteq M.$$

share|improve this answer
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.