A Hierarchy of Probability, Fluid and Generalized Densities for the Eulerian Velocivolumetric Description of Fluid Flow, for New Families of Conservation Laws

Abstract

The Reynolds transport theorem occupies a central place in continuum mechanics, providing a generalized integral conservation equation for the transport of any conserved quantity within a fluid or material volume, which can be connected to its corresponding differential equation. Recently, a more generalized framework was presented for this theorem, enabling parametric transformations between positions on a manifold or in any generalized coordinate space, exploiting the underlying continuous multivariate (Lie) symmetries of a vector or tensor field associated with a conserved quantity. We explore the implications of this framework for fluid flow systems, based on an Eulerian velocivolumetric (position-velocity) description of fluid flow. The analysis invokes a hierarchy of five probability density functions, which by convolution are used to define five fluid densities and generalized densities relevant to this description. We derive 11 formulations of the generalized Reynolds transport theorem for different choices of the coordinate space, parameter space and density, only the first of which is commonly known. These are used to generate a table of integral and differential conservation laws applicable to each formulation, for eight important conserved quantities (fluid mass, species mass, linear momentum, angular momentum, energy, charge, entropy and probability). The findings substantially expand the set of conservation laws for the analysis of fluid flow and dynamical systems.

1. Introduction

Near the end of his distinguished career, Osborne Reynolds presented what is now called the Reynolds transport theorem: a generalized conservation equation for the transport of a conserved quantity within a body of fluid (the domain, fluid volume or material volume) as it moves through a prescribed region of space (the control volume) [1]. This provides a universal formulation for the construction of integral conservation equations for any conserved quantity, and can be used to derive the corresponding differential equations for these quantities, e.g., [2,3,4,5,6]. Such conservation laws—founded on the paradigm of the field or Eulerian description of fluid flow—provide the basis for most theoretical and numerical analyses of flow systems. Extensions of the Reynolds transport theorem have been presented for moving and smoothly-deforming control volumes [3,6,7], domains with fixed or moving discontinuities [7,8], irregular and rough domains [9,10], two-dimensional domains [11,12,13,14], and differentiable manifolds using the formalism of exterior calculus [15,16,17,18,19]. The Reynolds transport theorem is also a special case of the Helmholtz transport theorem for flow through an open and moving surface [20], and of stochastic formulations to incorporate random diffusion and uncertainty [21,22].

Traditionally, the Reynolds transport theorem has been viewed exclusively as a continuous one-parameter (temporal) mapping of the density of a conserved quantity in geometric space, along the pathlines described by a time-dependent velocity vector field. Indeed, the above formulations all conform to this tradition. However, Flanders [15] interpreted the theorem more broadly as a generalization of the Leibniz rule for differentiation of an integral, rather than simply a construct of continuum mechanics. It is therefore far more general and powerful than the traditional interpretation might suggest. Using this insight, a generalized framework for the Reynolds transport theorem has recently been presented [23,24], based on continuous multiparametric mappings of a differential form on a manifold—or of a density within a generalized coordinate space—connected by the maximal integral curves or surfaces described by a vector or tensor field. This extends the traditional interpretation to encompass new transformation theorems in any parameter space, not just in time; these exploit previously unreported multiparametric continuous (Lie) symmetries associated with a conserved quantity in the space considered. These can be used, for example, to connect different positions in a velocity space connected by a velocity gradient tensor field, different positions in a Fourier spectral space connected by a velocity-wavenumber tensor field, or different positions in a velocity and chemical species space connected by velocity and concentration gradient tensor fields, for the analysis of chemical reaction systems and plasmas [23]. The generalized framework also yields new forms of the Liouville equation for the conservation of probability in different spaces, and of the Perron-Frobenius and Koopman operators for the extrapolation of probability densities or observable densities in such systems [23].

The aim of this work is to examine the implications of the generalized Reynolds transport theorem—and consequential integral and differential conservation laws—for an Eulerian velocivolumetric (position-velocity) description of fluid flow systems. The analysis commences in Section 2 with a detailed discussion of the extended Eulerian description, and of the properties of the volumetric and velocimetric domains for several well-known classes of fluid flow systems. This leads in Section 3 to a hierarchy of densities, starting in Section 3.1 with five probability density functions (pdfs), which are defined and in which their commutative relations are examined. These are used in Section 3.2 to define five corresponding fluid or material densities, of which only the volumetric density $\rho$ is commonly used in continuum mechanics; the remaining four have many similarities to other densities (such as the phase space density) used in other branches of physics. The fluid densities are formally defined from the pdfs by mathematical convolutions, for which the definitions and philosophical implications are discussed in Appendices Appendix A and Appendix B. The fluid densities are then used in Section 3.3 to define corresponding generalized densities for any conserved quantity. In Section 4, we present the generalized framework for the Reynolds transport theorem in both exterior calculus and vector calculus formulations. As explained in Section 4 and Appendix C, this can also be used to extract Lie or partial differential equations for an individual fluid element, the former expressed in terms of the Lie derivative of a volume form in the domain. In Section 5, these equations are then used to generate 11 formulations of the Reynolds transport theorem arising from the velocivolumetric description, for different choices of the coordinate space, parameter space and density. Of these, only the first accords with the traditional Reynolds transport theorem [1]. For each formulation, a table of integral and differential conservation equations is presented for the eight conserved quantities commonly considered in fluid mechanics (fluid mass, species mass, linear momentum, angular momentum, energy, charge, entropy and probability). The analyses provide a considerable assortment of new conservation laws for the analysis of fluid flow systems.

In the following sections, the mathematical notation is defined when first used and is also listed in

Table 1. Nomenclature used in this study.

Symbol	Description	SI Units
Mathematical Operators
${}^{\top }$	transpose
·	vector scalar product
×	cross product; multiplication symbol (only if there is a line break in the equation)
∂	partial derivative operator; boundary of domain
${\mathbf{\partial }}_{\mathit{C}}$	vector partial derivative operator with respect to $\mathit{C}$
${\nabla }_{\mathit{C}}$	gradient operator with respect to $\mathit{C}$
${\nabla }_{\mathit{u}}$	velocital gradient operator with respect to $\mathit{u}$	(m s−1)−1
${\nabla }_{\mathit{x}}$	spatial gradient operator with respect to $\mathit{x}$	m−1
${\nabla }_{\mathit{X}}$	gradient operator with respect to $\mathit{X}$
${▾}_{\mathit{u}}={\nabla }_{\mathit{u},t}$	velocitotemporal gradient operator with respect to $[\mathit{u},t]$	[(m s−1)$-1$, (m s$-1$)$-1$, (m s$-1$)$-1$, s$-1$]
${▾}_{\mathit{x}}={\nabla }_{\mathit{x},t}$	spatiotemporal gradient operator with respect to $[\mathit{x},t]$	[m$-1$, m$-1$, m$-1$, s$-1$]
$[·]$	expectation over small volume $\mathcal{U}$
$\langle ·\rangle$	expectation over small velocity domain $\mathcal{V}$
$\left[\right[[·]\left]\right]$	integral over volume $\mathsf{\Omega }$
$\langle \langle \langle ·\rangle \rangle \rangle$	integral over velocity domain $\mathcal{D}$
∧	wedge product
∟	augmentation operator, such that $\mathit{V}∟\mathit{C}$ is the tensor $\mathit{V}$ based on coordinates $\mathit{X}$ augmented by parameter $\mathit{C}$
➀	integral path as labeled in Figure 5
Conventions
Vector derivatives are defined by the $\partial (\to )/\partial (\downarrow )$ convention
The product of two vectors implies a tensor, e.g., $\mathit{u}\mathit{u}:=\mathit{u}{\mathit{u}}^{\top }$
The divergence of a tensor is rotated, e.g., ${\nabla }_{\mathit{x}}·\mathit{G}:={({\nabla }_{\mathit{x}}^{\top }\mathit{G})}^{\top }$
Roman symbols
c	index of chemical species; index of components of $\mathit{C}$
$C_c$	cth component of $\mathit{C}$
$\mathit{C}$	generalized m-dimensional parameter vector
$CV$	control volume = reference frame for fluid motion
d	differential of a function; exterior derivative of a differential form
$\widehat{d}$	extended exterior derivative based on augmented coordinates
$d/dt$	total derivative in time	s$-1$
$dA$	infinitesimal area element in volumetric space	m2
$d\mathit{A}$	directed infinitesimal area element in volumetric space	m2
$dB$	infinitesimal surface element in velocimetric space	(m s$-1$)2
$d\mathit{B}$	directed infinitesimal surface element in velocimetric space	(m s$-1$)2
$dm$	infinitesimal element of fluid mass	kg
$dU=dudvdw$	infinitesimal element of velocimetric space	(m s$-1$)3
$dV=dxdydz$	infinitesimal element of volumetric space	m3
$d{X}_j$	cotangent to jth component of generalized vector $\mathit{X}$
$d^{n-1}\mathit{X}$	generalized directed area element on $\partial \mathsf{\Omega }$
$d^n\mathit{X}$	generalized volume element in $\mathsf{\Omega }$
${\mathit{d}}_{\mathit{x},t}$	vector of spatiotemporal SI units	[m, m, m, s]
${\mathit{d}}_{\mathit{u},t}$	vector of velocitotemporal SI units	[m s$-1$, m s$-1$, m s$-1$, s]
$D/Dt$	substantial or material derivative in time	s$-1$
$\mathcal{D}$	velocimetric domain
$\underline{e}$	local specific total energy	J kg$-1$
$\breve{e}$	velocity-distinct specific energy	J kg$-1$
$\breve{\underline{e}}$	velocity-distinct local specific energy	J kg$-1$
E	total energy	J
$\sum \mathit{F}$	sum of forces	N
$\mathit{g}$	acceleration due to gravity	m s$-2$
$\mathit{G}:={\nabla }_{\mathit{x}}\mathit{u}$	velocity gradient tensor field	(m s$-1$) m$-1$ = s$-1$
$\tilde{\mathit{G}}:={▾}_{\mathit{x}}\mathit{u}$	augmented velocity gradient tensor field	[s$-1$, s$-1$, s$-1$, m s$-2$]
$\mathit{i}$	electrical flux	C m$-2$ s$-1$ = A m$-2$
$i_{\mathit{V}}^{\left(\mathit{C}\right)}$	multivariate interior product with respect to $\mathit{V}$ over parameters $\mathit{C}$
I	net inward electrical current (passive sign convention)	C s$-1$ = A
${\mathit{I}}_m$	identity matrix of size m
j	index of components of generalized coordinates $\mathit{X}$
${\mathit{j}}_c$	molar flux of species c	mol m$-2$ s$-1$
${\mathit{j}}_Q$	heat flux	J m$-2$ s$-1$
${\mathit{j}}_S$	entropy flux	J K$-1$ m$-2$ s$-1$
${\mathcal{L}}_{\mathit{V}}^{\left(\mathit{C}\right)}$	multivariate Lie derivative with respect to $\mathit{V}$ over parameters $\mathit{C}$
${\mathcal{L}}_{\mathit{V}∟\mathit{C}}^{\left(\mathit{C}\right)}$	multivariate Lie derivative with respect to $\mathit{V}∟\mathit{C}$ over parameters $\mathit{C}$
LHS	left-hand side
$\mathsf{m}$	small fluid mass domain
m	dimension of vector parameter $\mathit{C}$
$\mathsf{M}$	orientable differentiable manifold; generalized space
M	total fluid mass	kg
$\dot{M}$	rate of change of total fluid mass	kg s$-1$
${\dot{M}}_c$	rate of change of mass of species c	kgc s$-1$
${\mathcal{M}}_c$	molar mass of species c	kgc mol$-1$
n	dimension of manifold M, dimension of coordinates $\mathit{X}$
$\mathit{n}$	outward unit normal to $\mathsf{\Omega }$ in volumetric space
${\mathit{n}}_B$	outward unit normal to $\mathcal{D}$ in velocimetric space
$p(a,b|c)=p_{a,b|c}$	conditional pdf of a and b subject to c	units of ${\left(ab\right)}^{-1}$
P	pressure	Pa = J m$-3$
Q	total conserved quantity (of any type)	qty
${\dot{Q}}_{in}$	net inward heat flow rate	J s$-1$
r	dimension of submanifold $\mathsf{\Omega }$, dimension of differential form ${\omega }^r$
$\mathit{r}={[r_x,r_y,r_z]}^{\top }$	local Cartesian position coordinates	m
$\underline{\mathfrak{r}}$	local radius of a lever arm	m
$\breve{\mathfrak{r}}$	velocity-distinct radius of a lever arm	m
$\underline{\breve{\mathfrak{r}}}$	velocity-distinct local radius of a lever arm	m
RHS	right-hand side
$\underline{s}$	local specific entropy	J K$-1$ kg$-1$
$\breve{s}$	velocity-distinct specific entropy	J K$-1$ kg$-1$
$\breve{\underline{s}}$	velocity-distinct local specific entropy	J K$-1$ kg$-1$
$\mathit{s}={[s_u,s_v,s_w]}^{\top }$	local Cartesian velocity coordinates	m s$-1$
S	total entropy	J K$-1$
${\dot{S}}^{nf}$	total net inward non-fluid entropy flow rate	J K$-1$ s$-1$
t	time	s
$\sum \mathit{T}$	sum of torques	N m
$\mathit{u}={[u,v,w]}^{\top }$	Cartesian velocity field $:=\frac{\partial \mathit{x}}{\partial t}$	m s$-1$
$\dot{\mathit{u}}={[\dot{u},\dot{v},\dot{w}]}^{\top }$	Cartesian local acceleration field $:=\frac{\partial \mathit{u}}{\partial t}$	m s$-2$
$\mathcal{U}$	small velocity domain
${vol}_{\mathit{X}}^n$	volume of an infinitesimal n-dimensional parallelopiped spanned by the cotangents to $\mathit{X}$
$V_{jc}$	$\left(jc\right)$th component of generalized vector or tensor field $\mathit{V}$
$\mathit{V}$	generalized vector or tensor field
$\mathcal{V}$	small fluid volume
${\dot{W}}_{in}$	net inward work flow rate	J s$-1$
$\mathit{x}={[x,y,z]}^{\top }$	Cartesian position coordinates	m
${\mathit{x}}_0={[x_0,y_0,z_0]}^{\top }$	Cartesian Lagrangian position coordinates	m
$X_j$	jth component of vector $\mathit{X}$
$\mathit{X}$	generalized n-dimensional local or global Cartesian coordinates
$z_c$	charge per mass of species c	C kgc$-1$
$\underline{z}$	local specific charge	C kg$-1$
$\breve{z}$	velocity-distinct specific charge	C kg$-1$
$\breve{\underline{z}}$	velocity-distinct local specific charge	C kg$-1$
Z	total charge	C
Greek symbols
$\alpha =\rho \underline{\alpha }$	generalized volumetric density	qty m$-3$
$\underline{\alpha }$	local generalized specific density	qty kg$-1$
$\beta =Д\breve{\beta }$	generalized velocimetric density	qty (m s$-1$)$-3$
$\breve{\beta }$	velocity-distinct generalized specific density	qty kg$-1$
$\mathsf{\Gamma }:={\nabla }_{\mathit{u}}\mathit{x}$	inverse velocity gradient tensor field	m (m s$-1$)$-1$ = s
$\tilde{\mathsf{\Gamma }}:={▾}_{\mathit{u}}\mathit{x}$	augmented inverse velocity gradient tensor field,	[s, s, s, m s$-1$]
$\mathit{\delta }$	Kronecker delta tensor
$ϵ=\xi \underline{\breve{ϵ}}$	generalized conditional volumetric density	qty m$-3$
$\underline{\breve{ϵ}}$	velocity-distinct local generalized specific density	qty kg$-1$
$\zeta$	velocivolumetric fluid mass density	kg m$-3$ (m s$-1$)$-3$
${\zeta }_c$	velocivolumetric mass density of species c	kgc m$-3$ (m s$-1$)$-3$
$\eta$	conditional velocimetric fluid mass density	kg (m s$-1$)$-3$
${\eta }_c$	conditional velocimetric mass density of species c	kgc (m s$-1$)$-3$
$\theta =\eta \underline{\breve{\theta }}$	generalized conditional velocimetric density	qty (m s$-1$)$-3$
$\underline{\breve{\theta }}$	velocity-distinct local generalized specific density	qty kg$-1$
$\xi$	conditional volumetric fluid mass density	kg m$-3$
${\xi }_c$	conditional volumetric mass density of species c	kgc m$-3$
${\widehat{\dot{\xi }}}_c$	molar rate of production of species c	mol m$-3$ s$-1$
$\rho$	volumetric fluid mass density	kg m$-3$
${\rho }_c$	volumetric mass density of species c	kgc m$-3$
$\dot{\sigma }$	total entropy production	J K$-1$ s$-1$
$\widehat{\dot{\sigma }}$	local entropy production	J K$-1$ m$-3$ s$-1$
$\mathit{\tau }$	stress tensor (positive in compression)	Pa = J m$-3$
${\varphi }^{\mathit{C}}$	multivariate flow generated by $\mathit{V}$
${\widehat{\varphi }}^{\mathit{C}}$	augmented multivariate flow generated by $\mathit{V}∟\mathit{C}$
$\phi =\zeta \underline{\breve{\phi }}$	generalized velocivolumetric density	qty m$-3$ (m s$-1$)$-3$
$\underline{\breve{\phi }}$	velocity-distinct local generalized specific density	qty kg$-1$
${\underline{\chi }}_c$	local specific mass density of species c	kgc kg$-1$
${\breve{\chi }}_c$	velocity-distinct specific mass density of species c	kgc kg$-1$
${\breve{\underline{\chi }}}_c$	velocity-distinct local specific mass density of species c	kgc kg$-1$
$\psi$	generalized density of conserved quantity in generalized space
${\omega }^r,{\omega }^n$	r-form, n-form (respectively) in submanifold $\mathsf{\Omega }$
$\mathsf{\Omega }$	general submanifold or domain; volumetric domain (fluid volume or material volume)
Cyrillic symbols
Д	velocimetric fluid mass density	kg (m s$-1$)$-3$
${Д}_c$	velocimetric mass density of species c	kgc (m s$-1$)$-3$

.

2. The Position-Velocity Description and Domains

In continuum mechanics, fluid flow systems are commonly examined using the Eulerian description, in which each local property of the fluid is specified as a function of position (such as in Cartesian coordinates) $\mathit{x}={[x,y,z]}^{\top }\in \mathsf{\Omega }\subset {\mathbb{R}}^3$ and time $t\in \mathbb{R}$ as the fluid moves past, where $\mathsf{\Omega }$ is a three-dimensional geometric space and ${}^{\top }$ is the transpose. Thus, for example, fluid mechanicians commonly consider the three-dimensional velocity $\mathit{u}(\mathit{x},t)$, the volumetric mass density $\rho (\mathit{x},t)$ and the volumetric mass concentration ${\rho }_c(\mathit{x},t)$ of the cth chemical species within the Eulerian description. We here consider an extended velocivolumetric continuum description of a dynamical system based on Eulerian velocity and position coordinates, in which each local property of a fluid is specified as a function of the instantaneous fluid velocity $\mathit{u}={[u,v,w]}^{\top }\in \mathcal{D}\subset {\mathbb{R}}^3$, position $\mathit{x}\in \mathsf{\Omega }\subset {\mathbb{R}}^3$ and time $t\in \mathbb{R}$ as the fluid moves past, where $\mathcal{D}$ is a three-dimensional velocity space. This treatment—somewhat analogous to the phase space description used in many branches of physics—has the advantage of explicitly incorporating the velocity dependence of physical quantities, significantly extending the breadth of physical quantities that can be considered, and the scope and fidelity of the analyses.

We consider two alternative representations of the paired domains $\mathsf{\Omega }$ and $\mathcal{D}$:

(a)

The geometric representation—the usual physical viewpoint—in which $\mathcal{D}(\mathit{x},t)$ is a function of position and time, and $\mathsf{\Omega }\left(t\right)$ is a function of time. In this perspective, as illustrated in Figure 1a, there exists a map between each position $\mathit{x}\in \mathsf{\Omega }\left(t\right)$ and an entire velocity space $\mathcal{D}(\mathit{x},t)$, consisting of all possible velocities for this position and time.

(b)

The velocimetric representation—an alternative viewpoint—in which $\mathsf{\Omega }(\mathit{u},t)$ is a function of velocity and time, and $\mathcal{D}\left(t\right)$ is a function of time. In this perspective, as illustrated in Figure 1b, there exists a map between each velocity $\mathit{u}\in \mathcal{D}\left(t\right)$ and an entire geometric space $\mathsf{\Omega }(\mathit{u},t)$, consisting of all possible positions for this velocity and time.

In principle, the set of all ordered triples $(\mathit{u},\mathit{x},t)$ for a given flow system can be mapped into either of these representations, hence $\mathcal{D}(\mathit{x},t)\subseteq \mathcal{D}\left(t\right)$ for all $\mathit{x}\in \mathsf{\Omega }\left(t\right)$ and $\mathsf{\Omega }(\mathit{u},t)\subseteq \mathsf{\Omega }\left(t\right)$ for all $\mathit{u}\in \mathcal{D}\left(t\right)$.

Examining the velocity domain $\mathcal{D}\left(t\right)$ or $\mathcal{D}(\mathit{x},t)$, we make two important assertions. First, we expect $\mathcal{D}\left(t\right)$ or $\mathcal{D}(\mathit{x},t)$ to be continuous, since for most flows, it is physically impossible for a local velocity to change from ${\mathit{u}}_1$ to ${\mathit{u}}_2$ without passing through all intermediate velocities ${\mathit{u}}_1<\mathit{u}<{\mathit{u}}_2$, however fleetingly. The main exception to this rule are flows that cross a regime threshold, e.g., from laminar to turbulent flow, or subsonic to supersonic flow, leading to two internally continuous velocity domains ${\mathcal{D}}_1$ and ${\mathcal{D}}_2$, which may be disjoint. Second, the velocity cannot be infinite (positive or negative) for any physically realizable flow, as this would require local velocities of infinite kinetic energy. In consequence, for all flows the domains $\mathcal{D}\left(t\right)$ or $\mathcal{D}(\mathit{x},t)$ should be bounded, and for many flows will also be closed, and thence compact. We note in passing that integration of the velocity over ${\mathbb{R}}^3$ is a useful technique for many calculations, but this invokes an approximation that cannot be manifested physically. In its place, integration over a compact velocity domain $\mathcal{D}\left(t\right)$ or $\mathcal{D}(\mathit{x},t)$ is well-defined, while integration over bounded and open domains can be performed by careful consideration of the limits or, if necessary, by transformation to the Lebesgue integral.

In consequence, each domain in Figure 1 is drawn as compact and simply connected, a useful starting assumption for turbulent flow systems, but with many exceptions. For most flows involving a continuous fluid volume with no velocity discontinuities, the primary domains $\mathsf{\Omega }\left(t\right)$ and $\mathcal{D}\left(t\right)$ and the subsidiary domain $\mathcal{D}(\mathit{x},t)$ should be compact and simply connected, but some subsidiary geometric spaces $\mathsf{\Omega }(\mathit{u},t)$ may consist of disjoint subspaces, each associated with a different location (or set of locations) within the flow. Consider, for example, two-dimensional turbulent Poiseuille flow between parallel plates: each pair of positions $y_{+}$ and $y_{-}$ symmetric about the centerline will have the same mean velocity, hence many instantaneous velocities $\mathit{u}\in \mathcal{D}\left(t\right)$ will map to two disjoint geometric subspaces ${\mathsf{\Omega }}_{+}(\mathit{u},t)$ and ${\mathsf{\Omega }}_{-}(\mathit{u},t)$, respectively, containing $y_{+}$ and $y_{-}$. By breaking the original problem into smaller coupled flows, or by a judicious choice of coordinate system, it should be possible to isolate or unite these subspaces. In flows with different flow regimes, leading to disjoint velocity domains $\mathcal{D}\left(t\right)$ and $\mathcal{D}(\mathit{x},t)$, it should be possible to isolate each subdomain using a dimensionless discriminator (such as a Reynolds, Mach or Froude number). For some flows, for example the turbulent boundary layer, there are long-standing arguments over whether the overall geometric space $\mathsf{\Omega }\left(t\right)$ can be considered compact, due to a lack of boundedness or closedness, but despite this the subsidiary geometric spaces $\mathsf{\Omega }(\mathit{u},t)$ and the overall and subsidiary velocity domains $\mathcal{D}\left(t\right)$ and $\mathcal{D}(\mathit{x},t)$ will very likely be bounded and may also be closed. For homogeneous isotropic turbulence—a highly idealized flow—each subsidiary domain $\mathcal{D}(\mathit{x},t)$ or $\mathsf{\Omega }(\mathit{u},t)$ must be independent, respectively, of position or velocity, and so the two representations will collapse to give two separable domains $\mathsf{\Omega }\left(t\right)$ and $\mathcal{D}\left(t\right)$. In laminar flows, each subsidiary velocity domain $\mathcal{D}(\mathit{x},t)$ can be idealized as a single point (in reality, allowing for fluctuations, a small region), while each subsidiary geometric space $\mathsf{\Omega }(\mathit{u},t)$ will consist of single or multiple disconnected points (or small regions).

3. A Hierarchy of Densities

3.1. Probability Density Functions

We can now define the primary probability density functions (pdfs) that underlie continuum systems, and will subsequently be used to define physical densities in these systems. Writing the nonnegative real line as ${\mathbb{R}}_0^{+}$, the velocivolumetric description gives rise to the following five pdfs:

(a)

A volumetric pdf $p\left(\mathit{x}\right|t):\mathsf{\Omega }\times \mathbb{R}\to {\mathbb{R}}_0^{+}$ [SI units: m⁻³];

(b)

A velocimetric pdf $p\left(\mathit{u}\right|t):\mathcal{D}\times \mathbb{R}\to {\mathbb{R}}_0^{+}$ [(m s⁻¹)⁻³];

(c)

A velocivolumetric pdf $p(\mathit{u},\mathit{x}|t):\mathcal{D}\times \mathsf{\Omega }\times \mathbb{R}\to {\mathbb{R}}_0^{+}$ [m⁻³ (m s⁻¹)⁻³];

(d)

A conditional velocimetric (ensemble) pdf $p\left(\mathit{u}\right|\mathit{x},t):\mathcal{D}\times \mathsf{\Omega }\times \mathbb{R}\to {\mathbb{R}}_0^{+}$ [(m s⁻¹)⁻³]; and

(e)

A conditional volumetric pdf $p\left(\mathit{x}\right|\mathit{u},t):\mathsf{\Omega }\times \mathcal{D}\times \mathbb{R}\to {\mathbb{R}}_0^{+}$ [m⁻³];

where the solidus | is the conditional probability symbol, with the conditions listed to the right (these follow the standard pdf notation with a common symbol p, leading to a mixed signature and functional notation). Generally, the pdf at each point forms part of a probability density field defined throughout its domain.

Using the notation $dV=dxdydz=d^3\mathit{x}$ for an infinitesimal volume element and $dU=dudvdw=d^3\mathit{u}$ for an infinitesimal three-dimensional velocity element, the five pdfs will by definition satisfy the following nine relations:

$$\begin{array}{cc}\hfill 1& ={\int \int \int }_{\mathsf{\Omega }\left(t\right)}p\left(\mathit{x}\right|t)dV\hfill \end{array}$$

(1)

$$\begin{array}{cc}\hfill 1& ={\int \int \int }_{\mathcal{D}\left(t\right)}p\left(\mathit{u}\right|t)dU\hfill \end{array}$$

(2)

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill 1& ={\int \int \int }_{\mathsf{\Omega }\left(t\right)}{\int \int \int }_{\mathcal{D}(\mathit{x},t)}p(\mathit{u},\mathit{x}|t)dUdV={\int \int \int }_{\mathcal{D}\left(t\right)}{\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}p(\mathit{u},\mathit{x}|t)dVdU\hfill \end{array}\end{array}$$

(3)

$$\begin{array}{cc}\hfill p\left(\mathit{x}\right|t)& ={\int \int \int }_{\mathcal{D}(\mathit{x},t)}p(\mathit{u},\mathit{x}|t)dU\hfill \end{array}$$

(4)

$$\begin{array}{cc}\hfill p\left(\mathit{u}\right|t)& ={\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}p(\mathit{u},\mathit{x}|t)dV\hfill \end{array}$$

(5)

$$\begin{array}{cc}\hfill p\left(\mathit{u}\right|\mathit{x},t)& =\frac{p(\mathit{u},\mathit{x}|t)}{p\left(\mathit{x}\right|t)}=\frac{p(\mathit{u},\mathit{x}|t)}{{\int \int \int }_{\mathcal{D}(\mathit{x},t)}p(\mathit{u},\mathit{x}|t)dU}\hfill \end{array}$$

(6)

$$\begin{array}{cc}\hfill p\left(\mathit{x}\right|\mathit{u},t)& =\frac{p(\mathit{u},\mathit{x}|t)}{p\left(\mathit{u}\right|t)}=\frac{p(\mathit{u},\mathit{x}|t)}{{\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}p(\mathit{u},\mathit{x}|t)dV}\hfill \end{array}$$

(7)

$$\begin{array}{cc}\hfill 1& ={\int \int \int }_{\mathcal{D}(\mathit{x},t)}p\left(\mathit{u}\right|\mathit{x},t)dU\hfill \end{array}$$

(8)

$$\begin{array}{cc}\hfill 1& ={\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}p\left(\mathit{x}\right|\mathit{u},t)dV\hfill \end{array}$$

(9)

The connections between the five pdfs, and the roles of the different domain representations, are illustrated in the relational diagram in Figure 2. As evident, the geometric representation is adopted for the integration paths on the left-hand side of Figure 2, while the velocimetric representation is required on the right-hand side.

Designating the infinitesimal velocity element, position element and time interval, respectively, by $[\mathit{u},\mathit{u}+d\mathit{u}]$, $[\mathit{x},\mathit{x}+d\mathit{x}]$ and $[t,t+dt]$, we examine each pdf in turn:

(a)

The velocivolumetric pdf $p(\mathit{u},\mathit{x}|t)$ is the most fundamental of the pdfs, giving rise to $p\left(\mathit{x}\right|t)$ or $p\left(\mathit{u}\right|t)$ by the marginalization operations in Equations (4) and (5), and $p\left(\mathit{u}\right|\mathit{x},t)$ and $p\left(\mathit{x}\right|\mathit{u},t)$ by the conditioning operations in Equations (6) and (7). Physically—albeit imprecisely [25,26]—we can interpret $p(\mathit{u},\mathit{x}|t)dUdV$ as the joint probability of an infinitesimal fluid element having a velocity of $[\mathit{u},\mathit{u}+d\mathit{u}]$ and position in $[\mathit{x},\mathit{x}+d\mathit{x}]$, during the time interval $[t,t+dt]$.

(b)

The volumetric pdf $p\left(\mathit{x}\right|t)$ can be recognized as the common probabilistic descriptor for fluid flow systems, forming the basis of the fluid mechanics formulations of the Liouville and Fokker–Planck Equations [27,28,29,30,31], and allied to the volumetric density $\rho (\mathit{x},t)$. Physically, $p\left(\mathit{x}\right|t)dV$ can be interpreted as the probability that a fluid element is situated at the position $[\mathit{x},\mathit{x}+d\mathit{x}]$ in the time interval $[t,t+dt]$, regardless of velocity.

(c)

The velocimetric pdf $p\left(\mathit{u}\right|t)$ is rather strange. Physically, $p\left(\mathit{u}\right|t)dU$ can be interpreted as the probability of fluid elements within the control volume having a velocity of $[\mathit{u},\mathit{u}+d\mathit{u}]$ in the time interval $[t,t+dt]$, regardless of position.

(d)

To understand the conditional velocimetric pdf $p\left(\mathit{u}\right|\mathit{x},t)$, we interpret $p\left(\mathit{u}\right|\mathit{x},t)dU$ as the probability that a fluid element has a velocity of $[\mathit{u},\mathit{u}+d\mathit{u}]$, at the position $[\mathit{x},\mathit{x}+d\mathit{x}]$ and time $[t,t+dt]$. We therefore recognize $p\left(\mathit{u}\right|\mathit{x},t)$—typically but incorrectly written as $p\left(\mathit{u}\right)$—as the basis of the ensemble mean commonly used in the Reynolds-averaged Navier–Stokes (RANS) equations, and of the single-position correlation functions of turbulent fluid mechanics [32,33,34,35,36].

(e)

To understand the conditional volumetric pdf $p\left(\mathit{x}\right|\mathit{u},t)$, we interpret $p\left(\mathit{x}\right|\mathit{u},t)dV$ as the probability that a fluid element has a position in $[\mathit{x},\mathit{x}+d\mathit{x}]$, for a velocity of $[\mathit{u},\mathit{u}+d\mathit{u}]$ and time $[t,t+dt]$.

The infinitesimal intervals are necessitated by the fact that the probability itself vanishes at each point: for example, for the random variable X with values x, by definition $\mathrm{Prob}(X=x)$ $={lim}_{dx\to 0}\mathrm{Prob}(x\le X<x+dx)={lim}_{dx\to 0}{\int }_x^{x+dx}p\left(x\right)dx={\int }_x^{x}p\left(x\right)dx=0$, e.g., [25,37].

In addition to the pdfs in Figure 2, it is possible to consider joint pdfs with respect to time, including $p(\mathit{x},t)$, $p(\mathit{u},t)$ and $p(\mathit{u},\mathit{x},t)$. These require normalization over a time interval in addition to their volume and/or velocity domain(s). Such pdfs are closely associated with path-based formulations for the description of entire histories of events, e.g., [38,39,40], and are not considered further here.

3.2. Fluid or Material Densities

The above five pdfs can be used to define corresponding fluid or material mass densities, four of which are not commonly used for the analysis of continuum systems (these revert to standard signature and functional notation):

(a)

A volumetric fluid density $\rho :\mathsf{\Omega }\times \mathbb{R}\to {\mathbb{R}}_0^{+}$, $(\mathit{x},t)\mapsto \rho (\mathit{x},t)$ [kg m⁻³];

(b)

A velocimetric fluid density Д: $\mathcal{D}\times \mathbb{R}\to {\mathbb{R}}_0^{+}$, $(\mathit{u},t)\mapsto (\mathit{u},t)$ [kg (m s⁻¹)⁻³];

(c)

A velocivolumetric fluid density $\zeta :\mathcal{D}\times \mathsf{\Omega }\times \mathbb{R}\to {\mathbb{R}}_0^{+}$, $(\mathit{u},\mathit{x},t)\mapsto \zeta (\mathit{u},\mathit{x},t)$ [kg m⁻³ (m s⁻¹)⁻³];

(d)

A conditional velocimetric (ensemble) fluid density $\eta :\mathcal{D}\times \mathsf{\Omega }\times \mathbb{R}\to {\mathbb{R}}_0^{+}$, $(\mathit{u},\mathit{x},t)\mapsto \eta (\mathit{u},\mathit{x},t)$ [kg (m s⁻¹)⁻³]; and

(e)

A conditional volumetric fluid density $\xi :\mathcal{D}\times \mathsf{\Omega }\times \mathbb{R}\to {\mathbb{R}}_0^{+}$, $(\mathit{u},\mathit{x},t)\mapsto \xi (\mathit{u},\mathit{x},t)$ [kg m⁻³];

As with the pdfs, each density forms part of a density field, defined throughout its domain. The symbols for the last four densities, including the Cyrillic “de” character (from the transliteration of “density”), are chosen to not conflict with the most common notation of fluid mechanics.

For analysis, consider a material volume containing the fluid or material mass M, which in the absence of sources or sinks of fluid will be constant in time (n.b., systems with sources or sinks of fluid mass will require a more complicated treatment). First, we require the five fluid densities to satisfy the following nine relations:

$$\begin{array}{cc}\hfill M& ={\int \int \int }_{\mathsf{\Omega }\left(t\right)}\rho (\mathit{x},t)dV\hfill \end{array}$$

(10)

$$\begin{array}{cc}\hfill M& ={\int \int \int }_{\mathcal{D}\left(t\right)}Д(\mathit{u},t)dU\hfill \end{array}$$

(11)

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill M& ={\int \int \int }_{\mathsf{\Omega }\left(t\right)}{\int \int \int }_{\mathcal{D}(\mathit{x},t)}\zeta (\mathit{u},\mathit{x},t)dUdV={\int \int \int }_{\mathcal{D}\left(t\right)}{\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}\zeta (\mathit{u},\mathit{x},t)dVdU\hfill \end{array}\end{array}$$

(12)

$$\begin{array}{cc}\hfill \rho (\mathit{x},t)& ={\int \int \int }_{\mathcal{D}(\mathit{x},t)}\zeta (\mathit{u},\mathit{x},t)dU\hfill \end{array}$$

(13)

$$\begin{array}{cc}\hfill Д(\mathit{u},t)& ={\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}\zeta (\mathit{u},\mathit{x},t)dV\hfill \end{array}$$

(14)

$$\begin{array}{cc}\hfill \eta (\mathit{u},\mathit{x},t)& =\frac{\zeta (\mathit{u},\mathit{x},t)M}{\rho (\mathit{x},t)}=\frac{\zeta (\mathit{u},\mathit{x},t)}{\left({\displaystyle \frac{{\int \int \int }_{\mathcal{D}(\mathit{x},t)}\zeta (\mathit{u},\mathit{x},t)dU}{M}}\right)}=\frac{\zeta (\mathit{u},\mathit{x},t)}{\left({\displaystyle \frac{{\int \int \int }_{\mathcal{D}(\mathit{x},t)}\zeta (\mathit{u},\mathit{x},t)dU}{{\int \int \int }_{\mathsf{\Omega }\left(t\right)}{\int \int \int }_{\mathcal{D}(\mathit{x},t)}\zeta (\mathit{u},\mathit{x},t)dUdV}}\right)}\hfill \end{array}$$

(15)

$$\begin{array}{cc}\hfill \xi (\mathit{u},\mathit{x},t)& =\frac{\zeta (\mathit{u},\mathit{x},t)M}{Д(\mathit{u},t)}=\frac{\zeta (\mathit{u},\mathit{x},t)}{\left({\displaystyle \frac{{\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}\zeta (\mathit{u},\mathit{x},t)dV}{M}}\right)}=\frac{\zeta (\mathit{u},\mathit{x},t)}{\left({\displaystyle \frac{{\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}\zeta (\mathit{u},\mathit{x},t)dV}{{\int \int \int }_{\mathcal{D}\left(t\right)}{\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}\zeta (\mathit{u},\mathit{x},t)dVdU}}\right)}\hfill \end{array}$$

(16)

$$\begin{array}{cc}\hfill M& ={\int \int \int }_{\mathcal{D}(\mathit{x},t)}\eta (\mathit{u},\mathit{x},t)dU\hfill \end{array}$$

(17)

$$\begin{array}{cc}\hfill M& ={\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}\xi (\mathit{u},\mathit{x},t)dV\hfill \end{array}$$

(18)

Clearly these are analogs of Equations (1)–(9), modified only by the introduction of M into the definitions of $\eta$ and $\xi$. Note that the time dependence is lost in the integrations in Equations (10)–(12) and Equations (17) and (18); furthermore Equations 17) and (18) also lose their dependence, respectively, on $\mathit{x}$ or $\mathit{u}$. The connections between the different fluid densities are illustrated in the relational diagram in Figure 3.

Second, following a rich line of research, e.g., [41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68], the fluid or material density $\rho$ at each position $\mathit{x}$ and time t can be rigorously defined by integration over a small fluid volume $\mathcal{V}$—or equivalently a small fluid mass $\mathsf{m}$—for which there are two interpretations. In the common viewpoint, $\mathcal{V}$ must be sufficiently large to enable the fluid to be considered a continuum, e.g., [44,45,46,47,56,57]. Thus, in single phase systems, it must provide a “microscopic” scale large enough to average out the molecular phenomena, while for multiphase systems it may need to be larger than the dominant “macroscopic” scales [45,48,58]. However, in the contrary “relativist” viewpoint, $\mathcal{V}$ is not considered a property of the continuum, but is simply a characteristic of the measurement scale [53,54,59,61]. The analyses here are agnostic, encompassing both viewpoints; we also adopt a small velocimetric domain $\mathcal{U}$ for velocity averaging. For each small integral, we adopt the local Cartesian position coordinates $\mathit{r}={[r_x,r_y,r_z]}^{\top }\in \mathsf{\Omega }\subset {\mathbb{R}}^3$ and velocity coordinates $\mathit{s}={[s_u,s_v,s_w]}^{\top }\in \mathcal{D}\subset {\mathbb{R}}^3$ aligned with their corresponding global coordinates, and with their origin, respectively, at $\mathit{x}$ or $\mathit{u}$. The small domains $\mathcal{V}$ and $\mathcal{U}$ are functions, respectively, of $\mathit{x}$ and $\mathit{u}$, as well as time, and must conform to the two domain representations introduced earlier, bringing additional dependencies on $\mathit{r}$ and/or $\mathit{s}$ into the definition of $\zeta$. The five fluid or material densities can then be defined from their underlying pdfs by the following convolutions:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \left[\rho \right](\mathit{x},t)& ={\int }_{\mathsf{m}(\mathit{x},t)}p\left(\mathit{r}\right|t)dm(\mathit{x}+\mathit{r},t)={\int \int \int }_{\mathcal{V}(\mathit{x},t)}p\left(\mathit{r}\right|t)\rho (\mathit{x}+\mathit{r},t)dV(\mathit{r},t)\hfill \end{array}\end{array}$$

(19)

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \langle Д\rangle (\mathit{u},t)& ={\int }_{\mathsf{m}(\mathit{u},t)}p\left(\mathit{s}\right|t)dm(\mathit{u}+\mathit{s},t)={\int \int \int }_{\mathcal{U}(\mathit{u},t)}p\left(\mathit{s}\right|t)(\mathit{u}+\mathit{s},t)dU(\mathit{s},t)\hfill \end{array}\end{array}$$

(20)

$$\begin{array}{cc}\hfill [\langle \zeta \rangle ](u,x,t)& ={\int }_{\mathrm{m}(u,x,t)}p(s,r\mid t)dm(\mathit{u}+\mathit{s},\mathit{x}+\mathit{r},t)\hfill \\ \hfill & =\underset{\mathcal{V}(x,t)}{\int \int \int }\underset{\mathcal{U}(u,x+\mathit{r},t)}{\int \int \int }p(s,\mathit{r}\mid t)\zeta (\mathit{u}+\mathit{s},\mathit{x}+\mathit{r},t)dU(\mathit{s},\mathit{x}+\mathit{r},t)dV(s,\mathit{r},t)\hfill \\ \hfill & =\underset{\mathcal{U}(u,t)}{\int \int \int }\underset{\mathcal{V}(\mathit{u}+\mathit{s},x,t)}{\int \int \int }p(\mathit{s},\mathit{r}\mid t)\zeta (\mathit{u}+\mathit{s},\mathit{x}+\mathit{r},t)dV(\mathit{u}+\mathit{s},\mathit{r},t)dU(\mathit{s},\mathit{r},t)\hfill \end{array}$$

(21)

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \langle \eta \rangle (\mathit{u},\mathit{x},t)& ={\int }_{\mathsf{m}(\mathit{u},\mathit{x},t)}p\left(\mathit{s}\right|\mathit{x},t)dm(\mathit{u}+\mathit{s},\mathit{x},t)={\int \int \int }_{\mathcal{U}(\mathit{u},\mathit{x},t)}p\left(\mathit{s}\right|\mathit{x},t)\eta (\mathit{u}+\mathit{s},\mathit{x},t)dU(\mathit{s},\mathit{x},t)\hfill \end{array}\end{array}$$

(22)

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \left[\xi \right](\mathit{u},\mathit{x},t)& ={\int }_{\mathsf{m}(\mathit{u},\mathit{x},t)}p\left(\mathit{r}\right|\mathit{u},t)dm(\mathit{u},\mathit{x}+\mathit{r},t)={\int \int \int }_{\mathcal{V}(\mathit{u},\mathit{x},t)}p\left(\mathit{r}\right|\mathit{u},t)\xi (\mathit{u},\mathit{x}+\mathit{r},t)dV(\mathit{u},\mathit{r},t)\hfill \end{array}\end{array}$$

(23)

where $[·]$ is a local volumetric expectation, $\langle ·\rangle$ is a local velocimetric expectation, and $dm$ is an infinitesimal element of fluid mass. A more detailed discussion of these definitions is given in Appendix A, while the philosophical implications of the use of pdfs to define material densities are explored in Appendix B. Note that if each pdf is assumed uniformly distributed over its domain, each expected fluid density reduces to the product of its underlying pdf and the fluid mass, as would be obtained from dimensional considerations.

In this study, we adopt the mass integrals in Equations (19)–(23) as the primary definitions of the fluid densities, since the volumetric and velocimetric integrals require knowledge of point density terms that need to be defined. For convenience the expectation notations used in Equations (19)–(23) are now dropped.

How should we physically interpret the different fluid or material densities? Several schematic diagrams to aid their interpretation are given in Figure 4. From their functional forms and underlying pdfs, we see that $\xi$, $\eta$, $\zeta$ and $\rho$ are local densities, i.e., they apply to each infinitesimal volume element within the geometric space. In contrast, Д is a non-local density, applying over the entire geometric space $\mathsf{\Omega }$. Furthermore, $\xi$, $\eta$, $\zeta$ and Д also apply to infinitesimal velocity elements within the velocity domain. Considering each density in turn:

(a)

As shown in Figure 4a, the velocivolumetric density $\zeta$ represents the fluid mass per unit of velocimetric and geometric space carried by an infinitesimal fluid element of velocity $[\mathit{u},\mathit{u}+d\mathit{u}]$ through the infinitesimal control volume element at $[\mathit{x},\mathit{x}+d\mathit{x}]$, during the infinitesimal time interval $[t,t+dt]$. In consequence, $\zeta$ is both a velocity spectral density and a local volumetric density, accounting for the distribution of fluid mass with both velocity and position. As evident in Figure 3, $\zeta$ is central to the current formulation, giving the other fluid densities by marginalization or conditioning.

(b)

As shown in Figure 4b, the well-known volumetric fluid density $\rho$ represents the fluid mass per unit volume carried by the fluid through the infinitesimal control volume element at $[\mathit{x},\mathit{x}+d\mathit{x}]$, during the time interval $[t,t+dt]$. From Equation (13), $\rho$ is obtained by integration (marginalization) of $\zeta$ over the subsidiary velocity domain $\mathcal{D}(\mathit{x},t)$, consisting of all realizable velocities for this position and time. In well-behaved systems, $\mathcal{D}(\mathit{x},t)$ should consist of an infinitesimal trajectory (or trajectory bundle) in velocity space, from which it may be possible to calculate $\rho$ by line integration.

(c)

In contrast, as shown in Figure 4c, the velocimetric density Д represents the fluid mass per unit of velocimetric space transported by fluid elements of velocity $[\mathit{u},\mathit{u}+d\mathit{u}]$ throughout the control volume, during the time interval $[t,t+dt]$. This is a very strange, aggregated density field, representing the distribution of fluid mass across the velocity spectrum rather than with position, but nonetheless both it and its underlying pdf $p\left(\mathit{u}\right|t)$ are well-defined. From Equation (14), Д is obtained by integration (marginalization) of $\zeta$ over the subsidiary geometric space $\mathsf{\Omega }(\mathit{u},t)$, consisting of all realizable positions for this velocity and time. As discussed in Section 2 and illustrated in Figure 4c, in many flow systems $\mathsf{\Omega }(\mathit{u},t)$ will consist of several disjoint but continuous domains, which depending on the flow system may be bounded and may also be closed.

(d)

The conditional ensemble density $\eta$ (not illustrated) represents the fluid mass per unit velocimetric space carried by a fluid element of velocity $[\mathit{u},\mathit{u}+d\mathit{u}]$, at the position $[\mathit{x},\mathit{x}+d\mathit{x}]$ during the time interval $[t,t+dt]$. From Equation (15), $\eta$ is obtained by the ratio of $\zeta M$ and $\rho$, which can be interpreted as a conditioning operation over position. This removes the volume from the dimensions of $\eta$, giving the units of fluid mass per unit velocity space.

(e)

The conditional density $\xi$ (not illustrated) represents the fluid mass per unit volume carried by a fluid element in the position $[\mathit{x},\mathit{x}+d\mathit{x}]$, of velocity $[\mathit{u},\mathit{u}+d\mathit{u}]$ during the time interval $[t,t+dt]$. From Equation (16), $\xi$ is obtained by the ratio of $\zeta M$ and Д, which can be interpreted as a conditioning operation over velocity. This removes the velocity volume from the dimensions of $\xi$, giving the units of fluid mass per unit volume.

The distinctions between $\xi$, $\eta$ and $\zeta$ are therefore quite subtle, but since they arise from separate underlying pdfs, each density is mathematically well-defined. We further see that Figure 3 is an analog of Figure 2, exhibiting the same vertical symmetry, with the geometric representation evident on the left-hand side (integration paths in Equations (10), (13) and (17)), and the velocimetric representation evident on the right-hand side (integration paths in Equations (11), (14) and (18)).

3.3. Generalized Densities

Based on the above definitions, we can now construct five generalized densities of any conserved quantity carried by a fluid flow (strictly, density fields of an extensive variable):

(a)

Volumetric densities $\alpha :\mathsf{\Omega }\times \mathbb{R}\to {\mathbb{R}}_0^{+}$, $(\mathit{x},t)\mapsto \alpha (\mathit{x},t)$ [qty m⁻³];

(b)

Velocimetric densities $\beta :\mathcal{D}\times \mathbb{R}\to {\mathbb{R}}_0^{+}$, $(\mathit{u},t)\mapsto \beta (\mathit{u},t)$ [qty (m s⁻¹)⁻³];

(c)

Velocivolumetric densities $\phi :\mathcal{D}\times \mathsf{\Omega }\times \mathbb{R}\to {\mathbb{R}}_0^{+}$, $(\mathit{u},\mathit{x},t)\mapsto \phi (\mathit{u},\mathit{x},t)$ [qty m⁻³ (m s⁻¹)⁻³];

(d)

Conditional velocimetric (ensemble) densities $\theta :\mathcal{D}\times \mathsf{\Omega }\times \mathbb{R}\to {\mathbb{R}}_0^{+}$, $(\mathit{u},\mathit{x},t)\mapsto \theta (\mathit{u},\mathit{x},t)$ [qty (m s⁻¹)⁻³]; and

(e)

Conditional volumetric densities $ϵ:\mathcal{D}\times \mathsf{\Omega }\times \mathbb{R}\to {\mathbb{R}}_0^{+}$, $(\mathit{u},\mathit{x},t)\mapsto ϵ(\mathit{u},\mathit{x},t)$ [qty m⁻³];

where “qty” denotes the units of the conserved quantity. These can be defined by the following relations:

$$\begin{array}{cc}\hfill \alpha (\mathit{x},t)& =\rho (\mathit{x},t)\underline{\alpha }(\mathit{x},t)\hfill \end{array}$$

(24)

$$\begin{array}{cc}\hfill \beta (\mathit{u},t)& =Д(\mathit{u},t)\breve{\beta }(\mathit{u},t)\hfill \end{array}$$

(25)

$$\begin{array}{cc}\hfill \phi (\mathit{u},\mathit{x},t)& =\zeta (\mathit{u},\mathit{x},t)\underline{\breve{\phi }}(\mathit{u},\mathit{x},t)\hfill \end{array}$$

(26)

$$\begin{array}{cc}\hfill \theta (\mathit{u},\mathit{x},t)& =\eta (\mathit{u},\mathit{x},t)\underline{\breve{\theta }}(\mathit{u},\mathit{x},t)\hfill \end{array}$$

(27)

$$\begin{array}{cc}\hfill ϵ(\mathit{u},\mathit{x},t)& =\xi (\mathit{u},\mathit{x},t)\underline{\breve{ϵ}}(\mathit{u},\mathit{x},t)\hfill \end{array}$$

(28)

where $\underline{\alpha }$, $\breve{\beta }$, $\underline{\breve{\phi }}$, $\underline{\breve{\theta }}$ and $\underline{\breve{ϵ}}$ [qty kg⁻¹] are specific quantities, representing the quantity carried per unit fluid or material mass. For precision, these are labeled by an underline or a breve accent to designate their functional dependencies; e.g., for the specific energy, $\underline{e}$ indicates a local density (a function of $\mathit{x}$), $\breve{e}$ indicates a velocity-distinct density (a function of $\mathit{u}$), and $\underline{\breve{e}}$ indicates dependence on both $\mathit{u}$ and $\mathit{x}$ (hence $\underline{e}={\int \int \int }_{\mathcal{D}(\mathit{x},t)}\underline{\breve{e}}dU$ and $\breve{e}={\int \int \int }_{\mathsf{\Omega }\left(t\right)}\underline{\breve{e}}dV$). However, the specific momentum density (the local velocity) $\mathit{u}$ does not require these designations—being already velocity-dependent—provided that care is taken over its dependence on position.

We also define $Q\left(t\right)$ to be the total conserved quantity (of any type) in the volumetric domain, which due to sources or sinks of the conserved quantity will in general be a function of time t. The generalized densities will then satisfy the following nine integral relations:

$$\begin{array}{cc}\hfill Q\left(t\right)& ={\int \int \int }_{\mathsf{\Omega }\left(t\right)}\alpha (\mathit{x},t)dV\hfill \end{array}$$

(29)

$$\begin{array}{cc}\hfill Q\left(t\right)& ={\int \int \int }_{\mathcal{D}\left(t\right)}\beta (\mathit{u},t)dU\hfill \end{array}$$

(30)

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill Q\left(t\right)& ={\int \int \int }_{\mathsf{\Omega }\left(t\right)}{\int \int \int }_{\mathcal{D}(\mathit{x},t)}\phi (\mathit{u},\mathit{x},t)dUdV={\int \int \int }_{\mathcal{D}\left(t\right)}{\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}\phi (\mathit{u},\mathit{x},t)dVdU\hfill \end{array}\end{array}$$

(31)

$$\begin{array}{cc}\hfill \alpha (\mathit{x},t)& ={\int \int \int }_{\mathcal{D}(\mathit{x},t)}\phi (\mathit{u},\mathit{x},t)dU\hfill \end{array}$$

(32)

$$\begin{array}{cc}\hfill \beta (\mathit{u},t)& ={\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}\phi (\mathit{u},\mathit{x},t)dV\hfill \end{array}$$

(33)

$$\begin{array}{cc}\hfill \theta (\mathit{u},\mathit{x},t)& =\frac{\phi (\mathit{u},\mathit{x},t)Q\left(t\right)}{\alpha (\mathit{x},t)}={\displaystyle \frac{\phi (\mathit{u},\mathit{x},t)}{\left({\displaystyle \frac{{\int \int \int }_{\mathcal{D}(\mathit{x},t)}\phi (\mathit{u},\mathit{x},t)dU}{Q\left(t\right)}}\right)}}={\displaystyle \frac{\phi (\mathit{u},\mathit{x},t)}{\left({\displaystyle \frac{{\int \int \int }_{\mathcal{D}(\mathit{x},t)}\phi (\mathit{u},\mathit{x},t)dU}{{\int \int \int }_{\mathsf{\Omega }\left(t\right)}{\int \int \int }_{\mathcal{D}(\mathit{x},t)}\phi (\mathit{u},\mathit{x},t)dUdV}}\right)}}\hfill \end{array}$$

(34)

$$\begin{array}{cc}\hfill ϵ(\mathit{u},\mathit{x},t)& =\frac{\phi (\mathit{u},\mathit{x},t)Q\left(t\right)}{\beta (\mathit{u},t)}={\displaystyle \frac{\phi (\mathit{u},\mathit{x},t)}{\left({\displaystyle \frac{{\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}\phi (\mathit{u},\mathit{x},t)dV}{Q\left(t\right)}}\right)}}={\displaystyle \frac{\phi (\mathit{u},\mathit{x},t)}{\left({\displaystyle \frac{{\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}\phi (\mathit{u},\mathit{x},t)dV}{{\int \int \int }_{\mathcal{D}\left(t\right)}{\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}\phi (\mathit{u},\mathit{x},t)dVdU}}\right)}}\hfill \end{array}$$

(35)

$$\begin{array}{cc}\hfill Q\left(t\right)& ={\int \int \int }_{\mathcal{D}(\mathit{x},t)}\theta (\mathit{u},\mathit{x},t)dU\hfill \end{array}$$

(36)

$$\begin{array}{cc}\hfill Q\left(t\right)& ={\int \int \int }_{\mathsf{\Omega }(\mathit{u},t)}ϵ(\mathit{u},\mathit{x},t)dV\hfill \end{array}$$

(37)

Clearly these are analogs of Equations (10)–(18), containing $Q\left(t\right)$ rather than M. The connections between generalized densities are shown in the relational diagram in Figure 5 (compare Figure 2 and Figure 3). By synthesis of Equations (19)–(23) and (24)–(28), we can also extract the relations between each generalized density and its underlying pdf (see discussion in Appendix A).

Finally, we note that all formulations in this section are functions of time, and in general of both position and velocity. These can be simplified to give quantities and domains that are not functions of time, position or velocity, leading to a considerable assortment of reduced mathematical formulations.

4. Generalized Formulations of Conservation Equations

4.1. Exterior Calculus Formulations

We now examine generalized forms of the Reynolds transport theorem, which can be interpreted more broadly as transformation theorems, i.e., they provide a continuous mapping within a domain, described by the maximal integral curves of a vector or tensor field defined with respect to a parameter space. Traditionally, this is used to define a mapping between positions connected by a velocity field, parameterized by time, to give the usual Reynolds transport theorem [1]. However, as has been shown [23,24], this is not the only possible formulation.

For maximum generality we adopt an exterior calculus formulation, e.g., [15,16,17,18,69,70,71,72,73,74,75,76]. Consider an r-dimensional oriented compact submanifold ${\mathsf{\Omega }}^r$ within an n-dimensional orientable differentiable manifold ${\mathsf{M}}^n$, described using a patchwork of local coordinate systems. Let $\mathit{V}$ be a vector or tensor field in ${\mathsf{M}}^n$, a function of the local coordinates and parameterized by the m-dimensional parameter vector $\mathit{C}$, but not a function of $\mathit{C}$. The components $C_c$ of $\mathit{C}$ are assumed orthogonal. The field trajectories (tangent bundles) of $\mathit{V}$ define the continuous multivariate map (“flow”) ${\varphi }^{\mathit{C}}:{\mathsf{M}}^n\to {\mathsf{M}}^n$ such that $\mathit{V}={(\partial {\varphi }^{\mathit{C}}/\partial \mathit{C})}^{\top }$. This is linear and invertible, and can be used to map the entire submanifold. If ${\omega }^r$ is an r-form representing a conserved quantity, its integral over ${\mathsf{\Omega }}^r$ can be proved to satisfy [23]:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill d\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }{\omega }^r& =\left[\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }{\mathcal{L}}_{\mathit{V}}^{\left(\mathit{C}\right)}{\omega }^r\right]·d\mathit{C}\hfill \end{array}\begin{array}{c}\hfill =\left[\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }{i}_{\mathit{V}}^{\left(\mathit{C}\right)}d{\omega }^r+\underset{\partial \mathsf{\Omega }\left(\mathit{C}\right)}{\oint }{i}_{\mathit{V}}^{\left(\mathit{C}\right)}{\omega }^r\right]·d\mathit{C}=\left[\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }{i}_{\mathit{V}}^{\left(\mathit{C}\right)}d{\omega }^r+d\left(i_{\mathit{V}}^{\left(\mathit{C}\right)}{\omega }^r\right)\right]·d\mathit{C}\end{array}\end{array}$$

(38)

where d is the exterior derivative, $\partial \mathsf{\Omega }$ is the submanifold boundary, “·” is the dot product, ${\mathcal{L}}_{\mathit{V}}^{\left(\mathit{C}\right)}$ is a multivariate Lie derivative with respect to $\mathit{V}$ over parameters $\mathit{C}$, and $i_{\mathit{V}}^{\left(\mathit{C}\right)}$ is a multivariate interior product with respect to $\mathit{V}$ over parameters $\mathit{C}$. The multivariate operators provide vector extensions of their usual one-parameter definitions in exterior calculus [23], while the last step in Equation (38) invokes Stokes’ theorem, imposing a regularity condition on ${\omega }^r$ in the submanifold $\mathsf{\Omega }$. For $\mathit{C}=t$, Equation (38) reduces to the one-parameter exterior calculus formulation of the Reynolds transport theorem [18].

If the vector or tensor field $\mathit{V}$ is also a function of $\mathit{C}$, the problem can be analyzed by augmenting the manifold with the parameter space, to define the flow ${\widehat{\varphi }}^C:{\mathsf{M}}^n\times {\mathbb{R}}^m\to {\mathsf{M}}^n\times {\mathbb{R}}^m$ based on the augmented $(n+m)\times m$ tensor field $\mathit{V}∟\mathit{C}={(\partial {\widehat{\varphi }}^{\mathit{C}}/\partial \mathit{C})}^{\top }$, where ∟ denotes this composition [23]. This is given by $\mathit{V}∟\mathit{C}=\left[\begin{array}{c}\mathit{V}\\ {\mathit{I}}_m\end{array}\right]$, where ${\mathit{I}}_m$ is the identity matrix of size m. Applying Equation (38) based on $\mathit{V}∟\mathit{C}$, this simplifies to the extended theorem [23]:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \widehat{d}\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }{\omega }^r& =\left[\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }{\mathcal{L}}_{\mathit{V}∟\mathit{C}}^{\left(\mathit{C}\right)}{\omega }^r\right]·d\mathit{C}=\left[\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }{\mathbf{\partial }}_{\mathit{C}}{\omega }^r+i_{\mathit{V}}^{\left(\mathit{C}\right)}d{\omega }^r+d\left(i_{\mathit{V}}^{\left(\mathit{C}\right)}{\omega }^r\right)\right]·d\mathit{C}\hfill \end{array}\end{array}$$

(39)

where $\widehat{d}$ is the extended exterior derivative based on the augmented coordinates, and ${\mathbf{\partial }}_{\mathit{C}}$ is a vector partial derivative operator with respect to the components of $\mathit{C}$. For Cartesian parameters $\mathit{C}$, ${\mathbf{\partial }}_{\mathit{C}}={\nabla }_{\mathit{C}}$, while for $\mathit{X}=\mathit{x}$, $\mathit{C}=t$ and $\mathit{V}=\mathit{u}$, Equation (39) reduces to the known exterior calculus formulation for a time-varying velocity field $\mathit{u}\left(t\right)$ [15,18,22].

We can now extract the Lie differential equations applicable to each differential form in $\mathsf{\Omega }$, respectively, from Equations (38) and (39):

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill {\mathcal{L}}_{\mathit{V}}^{\left(\mathit{C}\right)}{\omega }^r& =i_{\mathit{V}}^{\left(\mathit{C}\right)}d{\omega }^r+d\left(i_{\mathit{V}}^{\left(\mathit{C}\right)}{\omega }^r\right)\hfill \end{array}\end{array}$$

(40)

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill {\mathcal{L}}_{\mathit{V}∟\mathit{C}}^{\left(\mathit{C}\right)}{\omega }^r& ={\mathbf{\partial }}_{\mathit{C}}{\omega }^r+{\mathcal{L}}_{\mathit{V}}^{\left(\mathit{C}\right)}{\omega }^r={\mathbf{\partial }}_{\mathit{C}}{\omega }^r+i_{\mathit{V}}^{\left(\mathit{C}\right)}d{\omega }^r+d\left(i_{\mathit{V}}^{\left(\mathit{C}\right)}{\omega }^r\right)\hfill \end{array}\end{array}$$

(41)

The former is the multivariate extension of Cartan’s relation of exterior calculus e.g., [17,18,72], while the latter provides an extended form for a vector field $\mathit{V}\left(\mathit{C}\right)$ that is a function of $\mathit{C}$.

4.2. Vector Calculus Formulations

Equations (38) and (39) provide very general equations applicable to submanifolds of any dimension in a manifold. For a system with global coordinates, these can be simplified to give a generalized parametric Reynolds transport theorem. Consider an n-dimensional compact domain $\mathsf{\Omega }$ within an n-dimensional space $\mathsf{M}$, described by the global Cartesian coordinates $\mathit{X}$. Let $\mathit{V}={(\partial \mathit{X}/\partial \mathit{C})}^{\top }={\left({\nabla }_{\mathit{C}}\mathit{X}\right)}^{\top }$ be a vector or tensor field in $\mathsf{M}$, using the $\partial (\to )/\partial (\downarrow )$ vector derivative convention (where → and ↓ refer to row and column vectors, respectively), in general with $\mathit{V}$ a function of $\mathit{C}$. Let ${\omega }^n$ be an n-dimensional volume form defined by:

$$\begin{array}{c}\hfill \omega =\psi d{X}_1\wedge \dots \wedge d{X}_n=\psi {vol}_{\mathit{X}}^n\end{array}$$

(42)

where $\psi (\mathit{X},\mathit{C})$ is the density of a conserved quantity, ∧ is the wedge product and ${vol}_{\mathit{X}}^n$ is the volume of an infinitesimal n-dimensional parallelopiped spanned by the cotangents to $\mathit{X}$. We assume $\psi$ is continuous and continuously differentiable with respect to $\mathit{X}$ and $\mathit{C}$ throughout $\mathsf{\Omega }$. It can be shown that Equation (39) then reduces to [23]:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill d\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }\psi d^n\mathit{X}& =\left[{\nabla }_{\mathit{C}}\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }\psi d^n\mathit{X}\right]·d\mathit{C}\hfill \end{array}=\left[\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }{\nabla }_{\mathit{C}}\psi d^n\mathit{X}+\underset{\partial \mathsf{\Omega }\left(\mathit{C}\right)}{\oint }\psi \mathit{V}·d^{n-1}\mathit{X}\right]·d\mathit{C}=\left[\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }\left[{\nabla }_{\mathit{C}}\psi +{\nabla }_{\mathit{X}}·\left(\psi \mathit{V}\right)\right]d^n\mathit{x}\right]·d\mathit{C}\end{array}$$

(43)

where d is now the differential operator, ${\nabla }_{\mathit{X}}$ is the n-dimensional nabla operator with respect to $\mathit{X}$, ${\nabla }_{\mathit{C}}$ is the m-dimensional nabla operator with respect to $\mathit{C}$, $d^n\mathit{X}$ is a volume element in $\mathsf{\Omega }$, and $d^{n-1}\mathit{X}$ is a directed area element with outward unit normal on the boundary $\partial \mathsf{\Omega }$. For consistency with the vector derivative convention used, the divergence in Equation (43) is defined by ${\nabla }_{\mathit{X}}·\left(\psi \mathit{V}\right)={\left[{\nabla }_{\mathit{X}}^{\top }\left(\psi \mathit{V}\right)\right]}^{\top }$.

We emphasize that Equation (43) applies to a compact domain $\mathsf{\Omega }$ with smoothly varying densities, with coordinates $\mathit{X}$ measured with respect to a fixed frame of reference. Further extensions can also be derived for moving and smoothly deforming frames of reference (control domains) [3,6,23], domains with jump discontinuities [7,8], irregular and fragmenting domains [9,10] and stochastic flows [21,22].

We now combine the exterior calculus (Equation (39)) and vector calculus (Equation (43)) formulations, to provide a generalized Lie differential equation applicable to each volume form in $\mathsf{\Omega }$. Using the notation in the second part of Equation (42), we rewrite the first part of Equation (39):

$$\begin{array}{c}\hfill \widehat{d}\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }\psi {vol}_{\mathit{X}}^n=\left[\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }{\mathcal{L}}_{\mathit{V}∟\mathit{C}}^{\left(\mathit{C}\right)}\left(\psi {vol}_{\mathit{X}}^n\right)\right]·d\mathit{C}\end{array}$$

(44)

Since the exterior derivative of the integral is equivalent to its differential, and integration of $\psi$ with respect to ${vol}_{\mathit{X}}^n$ is identical to integration over $d^n\mathit{X}$, we also have:

$$\begin{array}{c}\hfill \widehat{d}\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }\psi {vol}_{\mathit{X}}^n=d\underset{\mathsf{\Omega }\left(\mathit{C}\right)}{\int }\psi d^n\mathit{X}\end{array}$$

(45)

Combining Equations (43)–(45) and equating integrands gives, for each volume form in $\mathsf{\Omega }$:

$$\begin{array}{c}\hfill {\mathcal{L}}_{\mathit{V}∟\mathit{C}}^{\left(\mathit{C}\right)}\left(\psi {vol}_{\mathit{X}}^n\right)=\left[{\nabla }_{\mathit{C}}\psi +{\nabla }_{\mathit{X}}·\left(\psi \mathit{V}\right)\right]d^n\mathit{X}\end{array}$$

(46)

This invokes the fundamental lemma of the calculus of variations, thus imposing a regularity assumption on $\psi$ and $\mathit{V}$ within $\mathsf{\Omega }$. The generalized Reynolds transport theorem in Equation (43) thus yields a Lie differential equation in the form of Equation (46), defined using the augmented field $\mathit{V}∟\mathit{C}$. In some flow systems, this can be simplified to give a partial differential equation applicable to each local element $d^n\mathit{X}$.

5. Example Flow Systems

In the following sections we explore different choices of the coordinates $\mathit{X}$, parameters $\mathit{C}$, vector or tensor field $\mathit{V}$ and generalized density $\psi$ in Equations (43) and (46), based on the velocivolumetric description and the hierarchy of densities developed in Section 2 and Section 3. This yields 11 different case study flow systems, for the six integration paths labeled in red in Figure 5, with some choices examined for both time-independent and time-dependent systems. To establish a consistent nomenclature, each system is named using its $\mathit{X}-\mathit{C}$ coordinates and its selected density. The domains $\mathsf{\Omega }$ and $\mathcal{D}$ are assumed compact, and the components of $\mathit{X}$ and $\mathit{C}$ are assumed orthogonal. In all cases we consider $\mathit{V}=\mathit{V}\left(\mathit{C}\right)$, but report only the intrinsic equations, with $\mathit{V}$ measured with respect to a fixed frame of reference. We also consider only smoothly-varying fields within a compact and simply-connected domain.

5.1. Volumetric-Temporal Formulation (Density $\alpha$)

We first consider the well-known volumetric-temporal formulation of the Reynolds transport theorem [1], based on the geometric space $\mathsf{\Omega }$ with Cartesian coordinates $\mathit{X}=\mathit{x}$, time parameter $\mathit{C}=t$, intrinsic velocity vector field $\mathit{V}=\mathit{u}(\mathit{x},t):=\partial \mathit{x}/\partial t$ and generalized density $\psi =\alpha (\mathit{x},t)$, defined from the fluid density $\rho (\mathit{x},t)$. This follows integration path ➀ in Figure 5. From Equation (43) and the definition of $\alpha$ in Equation (29), we obtain:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \frac{dQ}{dt}=\frac{DQ}{Dt}& =\frac{d}{dt}\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\alpha dV=\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\frac{\partial \alpha }{\partial t}dV+{\displaystyle \underset{\partial \mathsf{\Omega }\left(t\right)}{∯}}\alpha \mathit{u}·\mathit{n}dA=\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\left[\frac{\partial \alpha }{\partial t}+{\nabla }_{\mathit{x}}·\left(\alpha \mathit{u}\right)\right]dV\hfill \end{array}\end{array}$$

(47)

using $dV=d^3\mathit{x}$ and $\mathit{n}dA=d\mathit{A}=d^2\mathit{x}$, where $\mathit{n}$ is the outward unit normal, $dA$ is an area element and $d\mathit{A}$ is a directed area element. Equation (47) is commonly written in terms of the substantial derivative $D/Dt=\partial /\partial t+\mathit{u}·{\nabla }_{\mathit{x}}$, expressing the transport of the conserved quantity with the fluid volume.

To extract the differential equation, we require a local form of the left-hand side of Equation (47). Using the continuity equation:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill 0& =\frac{\partial \rho }{\partial t}+{\nabla }_{\mathit{x}}·\left(\rho \mathit{u}\right)\hfill \end{array}\end{array}$$

(48)

a simple manipulation of Equation (47) using the local specific density $\underline{\alpha }$ in Equation (24) gives:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \rho \frac{D\underline{\alpha }}{Dt}& =\frac{\partial \left(\rho \underline{\alpha }\right)}{\partial t}+{\nabla }_{\mathit{x}}·\left(\rho \underline{\alpha }\mathit{u}\right)\hfill \end{array}\end{array}$$

(49)

Further details of this derivation are given in [77] and Appendix C.

The integral and differential Equations (47) and (49) provide generalized forms of the standard conservation laws of fluid mechanics. In these equations, the left-hand terms are generally used as placeholders for any source-sink terms or driving forces for the conserved quantity to enter or leave the fluid volume or differential fluid element. Such equations for the seven common conserved quantities (fluid mass, chemical species mass, linear momentum, angular momentum, energy, charge and entropy) are listed in

Table 2. Conservation Laws for the (Well-Known) Volumetric-Temporal Formulation (based on the Volumetric Fluid Density $\rho (\mathit{x},t)$) after [2,3,4,6,77,78].

Conserved Quantity	Density	Integral Equation			Differential Equation
	$\alpha (\mathit{x},t)$$=\rho \underline{\alpha }$	LHS = ${\displaystyle \frac{dQ}{dt}}$ = ${\displaystyle \frac{DQ}{Dt}}$	= RHS	SI Units	LHS	= RHS	SI Units
Fluid mass	$\rho$	$0=\dot{M}$$={\displaystyle \frac{d\left[\right[\left[\rho \right]\left]\right]}{dt}}$	$=\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial \rho }{\partial t}}+{\nabla }_{\mathit{x}}·\left(\rho \mathit{u}\right)\right]dV$	[kg s$-1$]	0	$={\displaystyle \frac{\partial \rho }{\partial t}}+{\nabla }_{\mathit{x}}·\left(\rho \mathit{u}\right)$	[kg s$-1$ m$-3$]
Species mass	$\rho {\underline{\chi }}_c$$={\rho }_c$	$\left({\dot{M}}_c\right)={\displaystyle \frac{d\left[\left[\left[\rho {\underline{\chi }}_c\right]\right]\right]}{dt}}$	$=\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial \rho {\underline{\chi }}_c}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\rho {\underline{\chi }}_c\mathit{u}\right)\right]dV$	[kgc s$-1$]	$\rho {\displaystyle \frac{D{\underline{\chi }}_c}{Dt}}={\mathcal{M}}_c({\widehat{\dot{\xi }}}_c-{\nabla }_{\mathit{x}}·{\mathit{j}}_c)$	$={\displaystyle \frac{\partial \rho {\underline{\chi }}_c}{\partial t}}$$+{\nabla }_{\mathit{x}}·\left(\rho {\underline{\chi }}_c\mathit{u}\right)$	[kgc s$-1$ m$-3$]
Linear momentum	$\rho \mathit{u}$	$\sum \mathit{F}={\displaystyle \frac{d\left[\right[\left[\rho \mathit{u}\right]\left]\right]}{dt}}$	$=\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial \rho \mathit{u}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\rho \mathit{u}\mathit{u}\right)\right]dV$	[(kg m s$-1$) s$-1$ = N]	$\rho {\displaystyle \frac{D\mathit{u}}{Dt}}=-{\nabla }_{\mathit{x}}P-{\nabla }_{\mathit{x}}·\mathit{\tau }+\rho \mathit{g}$	$={\displaystyle \frac{\partial \rho \mathit{u}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\rho \mathit{u}\mathit{u}\right)$	[N m$-3$]
Angular momentum ${}^a$	$\rho (\underline{\mathfrak{r}}\times \mathit{u})$	$\sum \mathit{T}={\displaystyle \frac{d\left[\left[\left[\rho (\underline{\mathfrak{r}}\times \mathit{u})\right]\right]\right]}{dt}}$	$=\left\{\begin{array}{c}\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }[{\displaystyle \frac{\partial \rho (\underline{\mathfrak{r}}\times \mathit{u})}{\partial t}}\hfill \\ +{\nabla }_{\mathit{x}}·\left(\rho (\underline{\mathfrak{r}}\times \mathit{u})\mathit{u}\right)]dV\hfill \end{array}\right.$	[(kg m${}^2$ s$-1$) s$-1$ = N m]	$\left.\begin{array}{c}\hfill \rho {\displaystyle \frac{D(\underline{\mathfrak{r}}\times \mathit{u})}{Dt}}=-{\nabla }_{\mathit{x}}·(\underline{\mathfrak{r}}\times P\delta )\\ \hfill -{\nabla }_{\mathit{x}}·(\underline{\mathfrak{r}}\times \mathit{\tau })+(\underline{\mathfrak{r}}\times \rho \mathit{g})\end{array}\right\}$	$=\left\{\begin{array}{c}{\displaystyle \frac{\partial \rho (\underline{\mathfrak{r}}\times \mathit{u})}{\partial t}}\hfill \\ +{\nabla }_{\mathit{x}}·\left(\rho (\underline{\mathfrak{r}}\times \mathit{u})\mathit{u}\right)\hfill \end{array}\right.$	[N m$-2$ = (N m) m$-3$]
Energy	$\rho \underline{e}$	$\left.\begin{array}{c}\hfill \frac{DE}{Dt}=\frac{d\left[\left[\left[\rho \underline{e}\right]\right]\right]}{dt}\\ \hfill =({\dot{Q}}_{in}+{\dot{W}}_{in})\end{array}\right\}$	$=\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial \rho \underline{e}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\rho \underline{e}\mathit{u}\right)\right]dV$	[J s$-1$ = W]	$\left.\begin{array}{c}\hfill \rho {\displaystyle \frac{D\underline{e}}{Dt}}=-{\nabla }_{\mathit{x}}·{\mathit{j}}_Q-{\nabla }_{\mathit{x}}·\left(P\mathit{u}\right)\\ \hfill -{\nabla }_{\mathit{x}}·(\mathit{\tau }·\mathit{u})+\rho \mathit{g}·\mathit{u}\end{array}\right\}$	$={\displaystyle \frac{\partial \rho \underline{e}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\rho \underline{e}\mathit{u}\right)$	[J s$-1$ m$-3$ = W m$-3$]
Charge (in solution)	$\rho \underline{z}$	$\left.\begin{array}{c}\hfill \frac{DZ}{Dt}=\frac{d\left[\left[\left[\rho \underline{z}\right]\right]\right]}{dt}\\ \hfill =I+\left({\sum }_c{z}_c{\dot{M}}_c\right)\end{array}\right\}$	$=\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial \rho \underline{z}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\rho \underline{z}\mathit{u}\right)\right]dV$	[C s$-1$ = A]	$\left.\begin{array}{c}\hfill \rho {\displaystyle \frac{D\underline{z}}{Dt}}=-{\nabla }_{\mathit{x}}·\mathit{i}\\ \hfill +{\sum }_c{\mathcal{M}}_c{z}_c({\widehat{\dot{\xi }}}_c-{\nabla }_{\mathit{x}}·{\mathit{j}}_c)\end{array}\right\}$	$={\displaystyle \frac{\partial \rho \underline{z}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\rho \underline{z}\mathit{u}\right)$	[C s$-1$ m$-3$ = A m$-3$]
Entropy	$\rho \underline{s}$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{DS}{Dt}}={\displaystyle \frac{d\left[\left[\left[\rho \underline{s}\right]\right]\right]}{dt}}\\ \hfill =\dot{\sigma }+{\dot{S}}^{nf}\end{array}\right\}$	$=\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial \rho \underline{s}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\rho \underline{s}\mathit{u}\right)\right]dV$	[J K$-1$ s$-1$]	$\rho {\displaystyle \frac{D\underline{s}}{Dt}}=$$\widehat{\dot{\sigma }}-{\nabla }_{\mathit{x}}·{\mathit{j}}_S$	$={\displaystyle \frac{\partial \rho \underline{s}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\rho \underline{s}\mathit{u}\right)$	[J K$-1$ s$-1$ m$-3$]
Probability	$p\left(\mathit{x}\right|t)$$=p_{\mathit{x}|t}$	$0={\displaystyle \frac{d\left[\left[\left[p_{\mathit{x}|t}\right]\right]\right]}{dt}}$	$=\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial p_{\mathit{x}|t}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{x}|t}\mathit{u}\right)\right]dV$	[s$-1$]	0	$={\displaystyle \frac{\partial p_{\mathit{x}|t}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{x}|t}\mathit{u}\right)$	[s$-1$ m$-3$]

^a It can be shown that the differential equation for angular momentum reduces to τ = τ^⊤ [3,77].

. All symbols used are listed in Table 1 (note that some minor overlaps of symbols could not be avoided). The equations given in Table 2 contain typical source-sink terms for the left-hand side of each equation. More comprehensive versions for different representations or coupled phenomena can also be derived, e.g., [2,3,4,6,77,78]. To enable dimensional comparisons, the SI units of each integral and differential equation are also included in Table 2.

To the seven conservation laws, we can also add an eighth, by assigning $\alpha (\mathit{x},t)$ to its underlying pdf $p\left(\mathit{x}\right|t)$ (strictly, following the probabilistic averaging method of Appendix A, this is achieved by assigning $\alpha (\mathit{x}+\mathit{r},t)=p(\mathit{x}+\mathit{r}|t)$ in Equations (A3) and (A4) or $\underline{\alpha }(\mathit{x}+\mathit{r},t)=1$ [kg⁻¹] in Equation (A5), producing a probabilistic convolution). By normalization ${\int \int \int }_{\mathsf{\Omega }\left(t\right)}p\left(\mathit{x}\right|t)dV=1$, the left-hand term of the Reynolds transport theorem in Equation (47) then vanishes. In consequence, for a compactly supported continuous and continuously differentiable pdf $p\left(\mathit{x}\right|t)$ [79] we obtain the differential Equation [23]:

$$\begin{array}{c}\hfill 0=\frac{\partial p\left(\mathit{x}\right|t)}{\partial t}+{\nabla }_{\mathit{x}}·\left(p\left(\mathit{x}\right|t)\mathit{u}\right)\end{array}$$

(50)

This is the Liouville equation of fluid mechanics [27,29,30,31]. This and its integral form are included in Table 2. As will be shown, other Liouville equations based on different pdfs can be derived for other representations.

5.2. Velocimetric-Temporal Formulation (Density $\beta$)

Now consider a velocimetric-temporal formulation of the Reynolds transport theorem, based on the Eulerian velocity space $\mathcal{D}$ with Cartesian velocity coordinates $\mathit{X}=\mathit{u}$, time parameter $\mathit{C}=t$, local acceleration vector field $\mathit{V}=\dot{\mathit{u}}(\mathit{u},t):=\partial \mathit{u}/\partial t$ and generalized density $\psi =\beta (\mathit{u},t)$, defined from the velocimetric fluid density $Д(\mathit{u},t)$. This follows integration path ➁ in Figure 5. Recall that this requires $\dot{\mathit{u}}$, $\beta$ and Д to be defined in the velocimetric representation, for fluid elements of velocity $\mathit{u}$ aggregated over all positions. From Equation (43) and the definition of $\beta$ in Equation (30) [24]:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \frac{dQ}{dt}& =\frac{d}{dt}\underset{\mathcal{D}\left(t\right)}{\int \int \int }\beta dU=\underset{\mathcal{D}\left(t\right)}{\int \int \int }\frac{\partial \beta }{\partial t}dU{+}_{\partial \mathcal{D}\left(t\right)}\beta \dot{\mathit{u}}·{\mathit{n}}_{B}dB=\underset{\mathcal{D}\left(t\right)}{\int \int \int }\left[\frac{\partial \beta }{\partial t}+{\nabla }_{\mathit{u}}·\left(\beta \dot{\mathit{u}}\right)\right]dU\hfill \end{array}\end{array}$$

(51)

using $dU=d^3\mathit{u}$ and ${\mathit{n}}_{B}dB=d\mathit{B}=d^2\mathit{u}$, where ${\mathit{n}}_B$ is the outward unit normal, $dB$ is a velocimetric boundary element and $d\mathit{B}$ is a directed velocimetric boundary element.

Using the previous manipulation (Appendix C), we recover a velocimetric analog of the continuity equation at each velocity:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill 0& =\frac{\partial }{\partial t}+{\nabla }_{\mathit{u}}·\left(\dot{\mathit{u}}\right)\hfill \end{array}\end{array}$$

(52)

and in general the differential equation based on $\breve{\beta }$:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \frac{Дd\breve{\beta }}{dt}& =\frac{\partial \left(\breve{\mathrm{Д\beta }}\right)}{\partial t}+{\nabla }_{\mathit{u}}·\left(\breve{\mathrm{Д\beta }}\dot{\mathit{u}}\right)\hfill \end{array}\end{array}$$

(53)

The conservation laws derived from Equations (51) and (53) are listed with their SI units in

Table 3. Conservation Laws for the Velocimetric-Temporal Formulation (based on the Velocimetric Fluid Density $Д(\mathit{u},t)$).

Conserved Quantity	Density	Integral Equation			Differential Equation
	$\beta (\mathit{u},t)=Д\breve{\beta }$	LHS = ${\displaystyle \frac{dQ}{dt}}$ = ${\displaystyle \frac{DQ}{Dt}}$	= RHS	SI Units	LHS	= RHS	SI Units
Fluid mass	Д	$0=\dot{M}$$={\displaystyle \frac{d\langle \langle \langle Д\rangle \rangle \rangle }{dt}}$	$=\underset{\mathcal{D}\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial }{\partial t}}+{\nabla }_{\mathit{u}}·(Д\dot{\mathit{u}})\right]dU$	[kg s$-1$]	0	$={\displaystyle \frac{\partial Д}{\partial t}}+{\nabla }_{\mathit{u}}·(Д\dot{\mathit{u}})$	[kg s$-1$ (m s$-1$)$-3$]
Species mass	$Д{\breve{\chi }}_c{Д}_c$	$\left({\dot{M}}_c\right)={\displaystyle \frac{d\langle \langle \langle Д{\breve{\chi }}_c\rangle \rangle \rangle }{dt}}$	$=\underset{\mathcal{D}\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial Д{\breve{\chi }}_c}{\partial t}}+{\nabla }_{\mathit{u}}·(Д{\breve{\chi }}_c\dot{\mathit{u}})\right]dU$	[kgc s$-1$]	$Д{\displaystyle \frac{d{\breve{\chi }}_c}{dt}}$	$={\displaystyle \frac{\partial Д{\breve{\chi }}_c}{\partial t}}+{\nabla }_{\mathit{u}}·(Д{\breve{\chi }}_c\dot{\mathit{u}})$	[kgc s$-1$ (m s$-1$)$-3$]
Linear momentum	$Д\mathit{u}$	$\sum \mathit{F}={\displaystyle \frac{d\langle \langle \langle Д\mathit{u}\rangle \rangle \rangle }{dt}}$	$=\underset{\mathcal{D}\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial Д\mathit{u}}{\partial t}}+{\nabla }_{\mathit{u}}·(Д\mathit{u}\dot{\mathit{u}})\right]dU$	[(kg m s$-1$) s$-1$ = N]	$Д{\displaystyle \frac{d\mathit{u}}{dt}}$	$={\displaystyle \frac{\partial Д\mathit{u}}{\partial t}}+{\nabla }_{\mathit{u}}·(Д\mathit{u}\dot{\mathit{u}})$	[N (m s$-1$)$-3$]
Angular momentum	$Д(\breve{\mathfrak{r}}\times \mathit{u})$	$\sum \mathit{T}={\displaystyle \frac{d\langle \langle \langle Д(\breve{\mathfrak{r}}\times \mathit{u})\rangle \rangle \rangle }{dt}}$	$=\underset{\mathcal{D}\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial Д(\breve{\mathfrak{r}}\times \mathit{u})}{\partial t}}+{\nabla }_{\mathit{u}}·\left((Д\breve{\mathfrak{r}}\times \mathit{u})\dot{\mathit{u}}\right)\right]dU$	[(kg m${}^2$ s$-1$) s$-1$ = N m]	$Д{\displaystyle \frac{d(\breve{\mathfrak{r}}\times \mathit{u})}{dt}}$	$=\left\{\begin{array}{c}{\displaystyle \frac{\partial Д(\breve{\mathfrak{r}}\times \mathit{u})}{\partial t}}\hfill \\ +{\nabla }_{\mathit{u}}·\left(Д(\breve{\mathfrak{r}}\times \mathit{u})\dot{\mathit{u}}\right)\hfill \end{array}\right.$	[(N m) (m s$-1$)$-3$]
Energy	$Д\breve{e}$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{DE}{Dt}}={\displaystyle \frac{d\langle \langle \langle Д\breve{e}\rangle \rangle \rangle }{dt}}\\ \hfill =({\dot{Q}}_{in}+{\dot{W}}_{in})\end{array}\right\}$	$=\underset{\mathcal{D}\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial Д\breve{e}}{\partial t}}+{\nabla }_{\mathit{u}}·(Д\breve{e}\dot{\mathit{u}})\right]dU$	[J s$-1$ = W]	$Д{\displaystyle \frac{d\breve{e}}{dt}}$	$={\displaystyle \frac{\partial Д\breve{e}}{\partial t}}+{\nabla }_{\mathit{u}}·(Д\breve{e}\dot{\mathit{u}})$	[J s$-1$ (m s$-1$)$-3$ = W (m s$-1$)$-3$]
Charge (in solution)	$Д\breve{z}$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{DZ}{Dt}}={\displaystyle \frac{d\langle \langle \langle Д\breve{z}\rangle \rangle \rangle }{dt}}\\ \hfill =I+\left({\sum }_c{z}_c{\dot{M}}_c\right)\end{array}\right\}$	$=\underset{\mathcal{D}\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial Д\breve{z}}{\partial t}}+{\nabla }_{\mathit{u}}·(Д\breve{z}\dot{\mathit{u}})\right]dU$	[C s$-1$ = A]	$Д{\displaystyle \frac{d\breve{z}}{dt}}$	$={\displaystyle \frac{\partial Д\breve{z}}{\partial t}}+{\nabla }_{\mathit{u}}·(Д\breve{z}\dot{\mathit{u}})$	[C s$-1$ (m s$-1$)$-3$ = A (m s$-1$)$-3$]
Entropy	$Д\breve{s}$	${\displaystyle \frac{DS}{Dt}}={\displaystyle \frac{d\langle \langle \langle Д\breve{s}\rangle \rangle \rangle }{dt}}=\dot{\sigma }$$+{\dot{S}}^{nf}$	$=\underset{\mathcal{D}\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial Д\breve{s}}{\partial t}}+{\nabla }_{\mathit{u}}·(Д\breve{s}\dot{\mathit{u}})\right]dU$	[J K$-1$ s$-1$]	$Д{\displaystyle \frac{d\breve{s}}{dt}}$	$={\displaystyle \frac{\partial Д\breve{s}}{\partial t}}+{\nabla }_{\mathit{u}}·(Д\breve{s}\dot{\mathit{u}})$	[J K$-1$ s$-1$ (m s$-1$)$-3$]
Probability	$p\left(\mathit{u}\right|t)$$=p_{\mathit{u}|t}$	$0={\displaystyle \frac{d\langle \langle \langle p_{\mathit{u}|t}\rangle \rangle \rangle }{dt}}$	$=\underset{\mathcal{D}\left(t\right)}{\int \int \int }\left[{\displaystyle \frac{\partial p_{\mathit{u}|t}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u}|t}\dot{\mathit{u}}\right)\right]dU$	[s$-1$]	0	$={\displaystyle \frac{\partial p_{\mathit{u}|t}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u}|t}\dot{\mathit{u}}\right)$	[s$-1$ (m s$-1$)$-3$]

. Note that from Equations (47) and (51) (see also Figure 5), each integral on path ➁ is equal to the same rate of change as its volumetric counterpart on path ➀ (see Table 2). In contrast, the differential equations are localized by velocity rather than position, and so the source-sink terms are not immediately identifiable, but nonetheless can be written in the form of Equation (53), providing convenient placeholders for all source-sink terms. The corresponding integral and differential Liouville equations, obtained by assigning $\beta (\mathit{u},t)=p\left(\mathit{u}\right|t)$, are also listed in Table 3.

5.3. Velocivolumetric-Temporal Formulation (Density $\phi$)

Now consider a velocivolumetric-temporal formulation of the Reynolds transport theorem, based on the Eulerian velocity-position space $\mathcal{D}\times \mathsf{\Omega }$ with six-dimensional Cartesian coordinates $\mathit{X}=\left[\begin{array}{c}\mathit{u}\\ \mathit{x}\end{array}\right]$, time parameter $\mathit{C}=t$, composite vector field $\mathit{V}=\left[\begin{array}{c}\dot{\mathit{u}}\\ \mathit{u}\end{array}\right](\mathit{u},\mathit{x},t):=\partial \left[\begin{array}{c}\mathit{u}\\ \mathit{x}\end{array}\right]/\partial t$ and generalized density $\psi =\phi (\mathit{u},\mathit{x},t)$, defined from the velocivolumetric density $\zeta (\mathit{u},\mathit{x},t)$. This follows the double integration path ➂–➀ or ➃–➁ in Figure 5. For path ➂–➀, from Equation (43) and the definition of $\phi$ in Equation (31):

$$\begin{array}{cc}\hfill \frac{dQ}{dt}& =\frac{d}{dt}\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\underset{\mathcal{D}(x,t)}{\int \int \int }\phi dUdV\hfill \\ \hfill & =\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\underset{\mathcal{D}(x,t)}{\int \int \int }\frac{\partial \phi }{\partial t}dUdV+{\displaystyle \underset{\partial \mathsf{\Omega }\left(t\right)}{∯}}{\displaystyle \underset{\partial \mathcal{D}(\mathit{x},t)}{∯}}\phi \left[\begin{array}{c}\dot{u}\\ u\end{array}\right]·\left[\begin{array}{c}{n}_B\\ n\end{array}\right]dBdA\hfill \\ \hfill & =\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\underset{\mathcal{D}(x,t)}{\int \int \int }\left[\frac{\partial \phi }{\partial t}+{\nabla }_{u,x}·\left(\phi \left[\begin{array}{c}\dot{u}\\ u\end{array}\right]\right)\right]dUdV\hfill \\ \hfill & =\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\underset{\mathcal{D}(x,t)}{\int \int \int }\left[\frac{\partial \phi }{\partial t}+{\nabla }_u·\left(\phi \dot{u}\right)+{\nabla }_x·\left(\phi u\right)\right]dUdV\hfill \end{array}$$

(54)

The alternative path ➃–➁ can also be written using the second part of Equation (31). For separable integrals, further simplification is possible. From Equation (54) we can extract the continuity equation and differential equation based on $\underline{\breve{\phi }}$ (Appendix C), respectively:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill 0& =\frac{\partial \zeta }{\partial t}+{\nabla }_{\mathit{u}}·\left(\zeta \dot{\mathit{u}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \mathit{u}\right)\hfill \end{array}\end{array}$$

(55)

$$\begin{array}{cc}\hfill \zeta \frac{d\underline{\breve{\phi }}}{dt}& =\frac{\partial \left(\zeta \underline{\breve{\phi }}\right)}{\partial t}+{\nabla }_{\mathit{u}}·\left(\zeta \underline{\breve{\phi }}\dot{\mathit{u}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \underline{\breve{\phi }}\mathit{u}\right)\hfill \end{array}$$

(56)

The conservation laws derived from Equations (54) and (56) are listed with their SI units in

Table 4. Conservation Laws for the Velocivolumetric-Temporal Formulation (based on the Velocivolumetric Fluid Density $\zeta (\mathit{u},\mathit{x},t)$).

Conserved Quantity	Density	Integral Equation			Differential Equation
	$\phi (\mathit{u},\mathit{x},t)$$=\zeta \underline{\breve{\phi }}$	LHS = ${\displaystyle \frac{dQ}{dt}}$ = ${\displaystyle \frac{DQ}{Dt}}$	= RHS	SI Units	LHS	= RHS	SI Units
Fluid mass	$\zeta$	$0=\dot{M}$$={\displaystyle \frac{d\left[\right[[\langle \langle \langle \zeta \rangle \rangle \rangle ]\left]\right]}{dt}}$	$=\left\{\begin{array}{c}\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }[{\displaystyle \frac{\partial \zeta }{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \dot{\mathit{u}}\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta \mathit{u}\right)]dUdV\hfill \end{array}\right.$	[kg s$-1$]	0	$={\displaystyle \frac{\partial \zeta }{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \dot{\mathit{u}}\right)$$+{\nabla }_{\mathit{x}}·\left(\zeta \mathit{u}\right)$	[kg s$-1$ m$-3$ (m s$-1$)$-3$]
Species mass	$\zeta {\underline{\breve{\chi }}}_c$ = ${\zeta }_c$	$\left({\dot{M}}_c\right)={\displaystyle \frac{d\left[\left[\left[\langle \langle \langle \zeta {\underline{\breve{\chi }}}_c\rangle \rangle \rangle \right]\right]\right]}{dt}}$	$=\left\{\begin{array}{c}\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }[{\displaystyle \frac{\partial \zeta {\underline{\breve{\chi }}}_c}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta {\underline{\breve{\chi }}}_c\dot{\mathit{u}}\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta {\underline{\breve{\chi }}}_c\mathit{u}\right)]dUdV\hfill \end{array}\right.$	[kgc s$-1$]	$\zeta {\displaystyle \frac{d{\underline{\breve{\chi }}}_c}{dt}}$	$=\left\{\begin{array}{c}{\displaystyle \frac{\partial \zeta {\underline{\breve{\chi }}}_c}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta {\underline{\breve{\chi }}}_c\dot{\mathit{u}}\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta {\underline{\breve{\chi }}}_c\mathit{u}\right)\hfill \end{array}\right.$	[kgc s$-1$ m$-3$ (m s$-1$)$-3$]
Linear momentum	$\zeta \mathit{u}$	$\sum \mathit{F}={\displaystyle \frac{d\left[\right[[\langle \langle \langle \zeta \mathit{u}\rangle \rangle \rangle ]\left]\right]}{dt}}$	$=\left\{\begin{array}{c}\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }[{\displaystyle \frac{\partial \zeta \mathit{u}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \mathit{u}\dot{\mathit{u}}\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta \mathit{u}\mathit{u}\right)]dUdV\hfill \end{array}\right.$	[(kg m s$-1$) s$-1$ = N]	$\zeta {\displaystyle \frac{d\mathit{u}}{dt}}$	$=\left\{\begin{array}{c}{\displaystyle \frac{\partial \zeta \mathit{u}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \mathit{u}\dot{\mathit{u}}\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta \mathit{u}\mathit{u}\right)\hfill \end{array}\right.$	[N m$-3$ (m s$-1$)$-3$]
Angular momentum	$\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})$	$\left.\begin{array}{c}\hfill \sum \mathit{T}\\ \hfill ={\displaystyle \frac{d\left[\left[\left[\langle \langle \langle \zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\rangle \rangle \rangle \right]\right]\right]}{dt}}\end{array}\right\}$	$=\left\{\begin{array}{c}\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }[{\displaystyle \frac{\partial \zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\dot{\mathit{u}}\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\mathit{u}\right)]dUdV\hfill \end{array}\right.$	[(kg m${}^2$ s$-1$) s$-1$ = N m]	$\zeta {\displaystyle \frac{d(\underline{\breve{\mathfrak{r}}}\times \mathit{u})}{dt}}$	$=\left\{\begin{array}{c}{\displaystyle \frac{\partial \zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\dot{\mathit{u}}\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\mathit{u}\right)\hfill \end{array}\right.$	[(N m) m$-3$ (m s$-1$)$-3$]
Energy	$\zeta \underline{\breve{e}}$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{DE}{Dt}}={\displaystyle \frac{d\left[\left[\left[\langle \langle \langle \zeta \underline{\breve{e}}\rangle \rangle \rangle \right]\right]\right]}{dt}}\\ \hfill =({\dot{Q}}_{in}+{\dot{W}}_{in})\end{array}\right\}$	$=\left\{\begin{array}{c}\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }[{\displaystyle \frac{\partial \zeta \underline{\breve{e}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \underline{\breve{e}}\dot{\mathit{u}}\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta \underline{\breve{e}}\mathit{u}\right)]dUdV\hfill \end{array}\right.$	[J s$-1$ = W]	$\zeta {\displaystyle \frac{d\underline{\breve{e}}}{dt}}$	$={\displaystyle \frac{\partial \zeta \underline{\breve{e}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \underline{\breve{e}}\dot{\mathit{u}}\right)$$+{\nabla }_{\mathit{x}}·\left(\zeta \underline{\breve{e}}\mathit{u}\right)$	[J s$-1$ m$-3$ (m s$-1$)$-3$]
Charge (in solution)	$\zeta \underline{\breve{z}}$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{DZ}{Dt}}={\displaystyle \frac{d\left[\left[\left[\langle \langle \langle \zeta \underline{\breve{z}}\rangle \rangle \rangle \right]\right]\right]}{dt}}\\ \hfill =I+\left({\sum }_c{z}_c{\dot{M}}_c\right)\end{array}\right\}$	$=\left\{\begin{array}{c}\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }[{\displaystyle \frac{\partial \zeta \underline{\breve{z}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \underline{\breve{z}}\dot{\mathit{u}}\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta \underline{\breve{z}}\mathit{u}\right)]dUdV\hfill \end{array}\right.$	[C s$-1$ = A]	$\zeta {\displaystyle \frac{d\underline{\breve{z}}}{dt}}$	$={\displaystyle \frac{\partial \zeta \underline{\breve{z}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \underline{\breve{z}}\dot{\mathit{u}}\right)$$+{\nabla }_{\mathit{x}}·\left(\zeta \underline{\breve{z}}\mathit{u}\right)$	[C s$-1$ m$-3$ (m s$-1$)$-3$]
Entropy	$\zeta \underline{\breve{s}}$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{DS}{Dt}}={\displaystyle \frac{d\left[\left[\left[\langle \langle \langle \zeta \underline{\breve{s}}\rangle \rangle \rangle \right]\right]\right]}{dt}}\\ \hfill =\dot{\sigma }+{\dot{S}}^{nf}\end{array}\right\}$	$=\left\{\begin{array}{c}\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }[{\displaystyle \frac{\partial \zeta \underline{\breve{s}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \underline{\breve{s}}\dot{\mathit{u}}\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta \underline{\breve{s}}\mathit{u}\right)]dUdV\hfill \end{array}\right.$	[J K$-1$ s$-1$]	$\zeta {\displaystyle \frac{d\underline{\breve{s}}}{dt}}$	$=\left\{\begin{array}{c}{\displaystyle \frac{\partial \zeta \underline{\breve{s}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \underline{\breve{s}}\dot{\mathit{u}}\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta \underline{\breve{s}}\mathit{u}\right)\hfill \end{array}\right.$	[J K$-1$ s$-1$ m$-3$ (m s$-1$)$-3$]
Probability	$p(\mathit{u},\mathit{x}|t)$$=p_{\mathit{u},\mathit{x}|t}$	$0={\displaystyle \frac{d\left[\left[\left[\langle \langle \langle p_{\mathit{u},\mathit{x}|t}\rangle \rangle \rangle \right]\right]\right]}{dt}}$	$=\left\{\begin{array}{c}\underset{\mathsf{\Omega }\left(t\right)}{\int \int \int }\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }[{\displaystyle \frac{\partial p_{\mathit{u},\mathit{x}|t}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u},\mathit{x}|t}\dot{\mathit{u}}\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(p_{\mathit{u},\mathit{x}|t}\mathit{u}\right)]dUdV\hfill \end{array}\right.$	[s$-1$]	0	$=\left\{\begin{array}{c}{\displaystyle \frac{\partial p_{\mathit{u},\mathit{x}|t}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u},\mathit{x}|t}\dot{\mathit{u}}\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(p_{\mathit{u},\mathit{x}|t}\mathit{u}\right)\hfill \end{array}\right.$	[s$-1$ m$-3$ (m s$-1$)$-3$]

. Again the integrals equate to $dQ/dt$, hence to the same rates of change as the volumetric form (Table 2). The left-hand sides of the differential equations are written in the form of Equation (56). The corresponding Liouville equations, obtained from $\phi (\mathit{u},\mathit{x},t)=p(\mathit{u},\mathit{x}|t)$, are also listed in Table 4.

5.4. Velocimetric-Temporal Formulation (Density $\phi$)

Now consider a different velocimetric-temporal formulation, defined as in Section 5.2 but using the generalized density $\psi =\phi (\mathit{u},\mathit{x},t)$ based on $\zeta (\mathit{u},\mathit{x},t)$. This follows the partial integration path ➂ in Figure 5. From relation Equation (32 between $\phi$ and $\alpha$:

$$\begin{array}{c}\hfill \begin{array}{c}\hfill d\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\phi dU=d\alpha ={\nabla }_{\mathit{x}}\alpha ·d\mathit{x}+\frac{\partial \alpha }{\partial t}dt\end{array}\end{array}$$

(57)

Examining the time derivative term using Equation (43):

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \frac{\partial \alpha }{\partial t}& =\frac{\partial }{\partial t}\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\phi dU=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\frac{\partial \phi }{\partial t}dU{+}_{\partial \mathcal{D}(\mathit{x},t)}\phi \dot{\mathit{u}}·{\mathit{n}}_{B}dB=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[\frac{\partial \phi }{\partial t}+{\nabla }_{\mathit{u}}·\left(\phi \dot{\mathit{u}}\right)\right]dU\hfill \end{array}\end{array}$$

(58)

For the conserved quantities considered, it is not straightforward to extract a differential equation from Equation (58). Instead, it is necessary to write a Lie differential Equation (46) in terms of the Lie derivative ${\mathcal{L}}_{\dot{\mathit{u}}∟t}^{\left(t\right)}$ with respect to the augmented local acceleration field $\dot{\mathit{u}}∟t:=\left[\begin{array}{c}\dot{\mathit{u}}\\ 1\end{array}\right]$. The left-hand term can then be treated as a placeholder for sources, sinks or drivers of the conserved quantity in the volume form ${vol}_{\mathit{u}}^3$. The conservation laws derived from Equation (58) are listed in

Table 5. Conservation Laws for the Velocimetric-Temporal Formulation (based on the Velocivolumetric Fluid Density $\zeta (\mathit{u},\mathit{x},t)$).

Conserved Quantity	Density	Integral Equation			Differential Equation
	$\phi (\mathit{u},\mathit{x},t)$$=\zeta \underline{\breve{\phi }}$	LHS = ${\displaystyle \frac{\partial \alpha }{\partial t}}={\displaystyle \frac{\partial \left(\rho \underline{\alpha }\right)}{\partial t}}$	= RHS	SI Units	LHS	= RHS	SI Units
Fluid mass	$\zeta$	${\displaystyle \frac{\partial \rho }{\partial t}}$$={\displaystyle \frac{\partial \langle \langle \langle \zeta \rangle \rangle \rangle }{\partial t}}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \zeta }{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \dot{\mathit{u}}\right)\right]dU$	[kg m$-3$ s$-1$]	${\mathcal{L}}_{\dot{\mathit{u}}∟t}^{\left(t\right)}\left(\zeta {vol}_{\mathit{u}}^3\right)$	$=\left[{\displaystyle \frac{\partial \zeta }{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \dot{\mathit{u}}\right)\right]dU$	[kg m$-3$ s$-1$]
Species mass	$\zeta {\underline{\breve{\chi }}}_c$ = ${\zeta }_c$	${\displaystyle \frac{\partial \left(\rho {\underline{\chi }}_c\right)}{\partial t}}={\displaystyle \frac{\partial \langle \langle \langle \zeta {\underline{\breve{\chi }}}_c\rangle \rangle \rangle }{\partial t}}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \zeta {\underline{\breve{\chi }}}_c}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta {\underline{\breve{\chi }}}_c\dot{\mathit{u}}\right)\right]dU$	[kgc m$-3$ s$-1$]	${\mathcal{L}}_{\dot{\mathit{u}}∟t}^{\left(t\right)}\left(\zeta {\underline{\breve{\chi }}}_c{vol}_{\mathit{u}}^3\right)$	$=\left[{\displaystyle \frac{\partial \zeta {\underline{\breve{\chi }}}_c}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta {\underline{\breve{\chi }}}_c\dot{\mathit{u}}\right)\right]dU$	[kgc m$-3$ s$-1$]
Linear momentum	$\zeta \mathit{u}$	${\displaystyle \frac{\partial \left(\rho \mathit{u}\right)}{\partial t}}={\displaystyle \frac{\partial \langle \langle \langle \zeta \mathit{u}\rangle \rangle \rangle }{\partial t}}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \zeta \mathit{u}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \mathit{u}\dot{\mathit{u}}\right)\right]dU$	[(kg m s$-1$) m$-3$ s$-1$ = N m$-3$]	${\mathcal{L}}_{\dot{\mathit{u}}∟t}^{\left(t\right)}\left(\zeta \mathit{u}{vol}_{\mathit{u}}^3\right)$	$=[{\displaystyle \frac{\partial \zeta \mathit{u}}{\partial t}}$$+{\nabla }_{\mathit{u}}·\left(\zeta \mathit{u}\dot{\mathit{u}}\right)]dU$	[N m$-3$]
Angular momentum	$\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{\partial (\rho (\underline{\mathfrak{r}}\times \mathit{u}))}{\partial t}}\\ \hfill ={\displaystyle \frac{\partial \langle \langle \langle \zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\rangle \rangle \rangle }{\partial t}}\end{array}\right\}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\dot{\mathit{u}}\right)\right]dU$	[(kg m${}^2$ s$-1$) m$-3$ s$-1$ = (N m) m$-3$]	$\left.\begin{array}{c}\hfill {\mathcal{L}}_{\dot{\mathit{u}}∟t}^{\left(t\right)}\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\right)\\ \hfill {vol}_{\mathit{u}}^3)\end{array}\right\}$	$=\left[{\displaystyle \frac{\partial \zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\dot{\mathit{u}}\right)\right]dU$	[(N m) m$-3$]
Energy	$\zeta \underline{\breve{e}}$	${\displaystyle \frac{\partial \left(\rho \underline{e}\right)}{\partial t}}={\displaystyle \frac{\partial \langle \langle \langle \zeta \underline{\breve{e}}\rangle \rangle \rangle }{\partial t}}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \zeta \underline{\breve{e}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \underline{\breve{e}}\dot{\mathit{u}}\right)\right]dU$	[J m$-3$ s$-1$ = W m$-3$]	${\mathcal{L}}_{\dot{\mathit{u}}∟t}^{\left(t\right)}\left(\zeta \underline{\breve{e}}{vol}_{\mathit{u}}^3\right)$	$=\left[{\displaystyle \frac{\partial \zeta \underline{\breve{e}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \underline{\breve{e}}\dot{\mathit{u}}\right)\right]dU$	[J m$-3$ s$-1$ = W m$-3$]
Charge (in solution)	$\zeta \underline{\breve{z}}$	${\displaystyle \frac{\partial \left(\rho \underline{z}\right)}{\partial t}}={\displaystyle \frac{\partial \langle \langle \langle \zeta \underline{\breve{z}}\rangle \rangle \rangle }{\partial t}}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \zeta \underline{\breve{z}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \underline{\breve{z}}\dot{\mathit{u}}\right)\right]dU$	[C m$-3$ s$-1$ = A m$-3$]	${\mathcal{L}}_{\dot{\mathit{u}}∟t}^{\left(t\right)}\left(\zeta \underline{\breve{z}}{vol}_{\mathit{u}}^3\right)$	$=\left[{\displaystyle \frac{\partial \zeta \underline{\breve{z}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \underline{\breve{z}}\dot{\mathit{u}}\right)\right]dU$	[C m$-3$ s$-1$ = A m$-3$]
Entropy	$\zeta \underline{\breve{s}}$	${\displaystyle \frac{\partial \left(\rho \underline{s}\right)}{\partial t}}={\displaystyle \frac{\partial \langle \langle \langle \zeta \underline{\breve{s}}\rangle \rangle \rangle }{\partial t}}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \zeta \underline{\breve{s}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \underline{\breve{s}}\dot{\mathit{u}}\right)\right]dU$	[J K$-1$ m$-3$ s$-1$]	${\mathcal{L}}_{\dot{\mathit{u}}∟t}^{\left(t\right)}\left(\zeta \underline{\breve{s}}{vol}_{\mathit{u}}^3\right)$	$=\left[{\displaystyle \frac{\partial \zeta \underline{\breve{s}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\zeta \underline{\breve{s}}\dot{\mathit{u}}\right)\right]dU$	[J K$-1$ m$-3$ s$-1$]
Probability	$p(\mathit{u},\mathit{x}|t)$$=p_{\mathit{u},\mathit{x}|t}$	${\displaystyle \frac{\partial p_{\mathit{x}|t}}{\partial t}}={\displaystyle \frac{\partial \langle \langle \langle p_{\mathit{u},\mathit{x}|t}\rangle \rangle \rangle }{\partial t}}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial p_{\mathit{u},\mathit{x}|t}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u},\mathit{x}|t}\dot{\mathit{u}}\right)\right]dU$	[m$-3$ s$-1$]	${\mathcal{L}}_{\dot{\mathit{u}}∟t}^{\left(t\right)}\left(p_{\mathit{u},\mathit{x}|t}{vol}_{\mathit{u}}^3\right)$	$=\left[{\displaystyle \frac{\partial p_{\mathit{u},\mathit{x}|t}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u},\mathit{x}|t}\dot{\mathit{u}}\right)\right]dU$	[m$-3$ s$-1$]
	$p\left(\mathit{u}\right|\mathit{x},t)$$=p_{\mathit{u}|\mathit{x},t}$	$0={\displaystyle \frac{\partial \langle \langle \langle p_{\mathit{u}|\mathit{x},t}\rangle \rangle \rangle }{\partial t}}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial p_{\mathit{u}|\mathit{x},t}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u}|\mathit{x},t}\dot{\mathit{u}}\right)\right]dU$	[s$-1$]	0	$={\displaystyle \frac{\partial p_{\mathit{u}|\mathit{x},t}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u}|\mathit{x},t}\dot{\mathit{u}}\right)$	[(m s$-1$)$-3$ s$-1$]

.

For this formulation, it is possible to derive two sets of temporal Liouville equations, by assigning the density $\phi (\mathit{u},\mathit{x},t)$, respectively, to $p(\mathit{u},\mathit{x}|t)$ or $p\left(\mathit{u}\right|\mathit{x},t)$. The former gives a Lie differential equation, while the latter reduces by normalization in Equation (8 to a partial differential equation. Both sets are listed in Table 5.

5.5. Volumetric-Temporal Formulation (Density $\phi$)

We here consider a different volumetric-temporal formulation, defined as in Section 5.1 but using the generalized density $\psi =\phi (\mathit{u},\mathit{x},t)$ based on $\zeta (\mathit{u},\mathit{x},t)$. This follows the partial integration path ➃ in Figure 5. From Equation (33) between $\phi$ and $\beta$:

$$\begin{array}{c}\hfill \begin{array}{c}\hfill d\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\phi dV=d\beta ={\nabla }_{\mathit{u}}\beta ·d\mathit{u}+\frac{\partial \beta }{\partial t}dt\end{array}\end{array}$$

(59)

Examining the time derivative term using Equation (43):

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \frac{\partial \beta }{\partial t}& =\frac{\partial }{\partial t}\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\phi dV=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\frac{\partial \phi }{\partial t}dV{+}_{{\displaystyle \underset{\partial \mathsf{\Omega }\left(t\right)}{∯}}}\phi \mathit{u}·\mathit{n}dA=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[\frac{\partial \phi }{\partial t}+{\nabla }_{\mathit{x}}·\left(\phi \mathit{u}\right)\right]dV\hfill \end{array}\end{array}$$

(60)

We again extract Lie differential Equation (46), here written in terms of the Lie derivative ${\mathcal{L}}_{\mathit{u}∟t}^{\left(t\right)}$ with respect to the augmented velocity field $\mathit{u}∟t:=\left[\begin{array}{c}\mathit{u}\\ 1\end{array}\right]$. The conservation laws derived from Equation (60) are listed in

Table 6. Conservation Laws for the Volumetric-Temporal Formulation (based on the Velocivolumetric Fluid Density $\zeta (\mathit{u},\mathit{x},t)$).

Conserved Quantity	Density	Integral Equation			Differential Equation
	$\phi (\mathit{u},\mathit{x},t)=\zeta \underline{\breve{\phi }}$	LHS = ${\displaystyle \frac{\partial \beta }{\partial t}}={\displaystyle \frac{\partial (Д\breve{\beta })}{\partial t}}$	= RHS	SI Units	LHS	= RHS	SI Units
Fluid mass	$\zeta$	${\displaystyle \frac{\partial Д}{\partial t}}$$={\displaystyle \frac{\partial \left[\right[\left[\zeta \right]\left]\right]}{\partial t}}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \zeta }{\partial t}}+{\nabla }_{\mathit{x}}·\left(\zeta \mathit{u}\right)\right]dV$	[kg (m s$-1$)$-3$ s$-1$]	${\mathcal{L}}_{\mathit{u}∟t}^{\left(t\right)}\left(\zeta {vol}_{\mathit{x}}^3\right)$	$=\left[{\displaystyle \frac{\partial \zeta }{\partial t}}+{\nabla }_{\mathit{x}}·\left(\zeta \mathit{u}\right)\right]dV$	[kg (m s$-1$)$-3$ s$-1$]
Species mass	$\zeta {\underline{\breve{\chi }}}_c$ = ${\zeta }_c$	${\displaystyle \frac{\partial (Д{\breve{\chi }}_c)}{\partial t}}={\displaystyle \frac{\partial \left[\left[\left[\zeta {\underline{\breve{\chi }}}_c\right]\right]\right]}{\partial t}}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \zeta {\underline{\breve{\chi }}}_c}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\zeta {\underline{\breve{\chi }}}_c\mathit{u}\right)\right]dV$	[kgc (m s$-1$)$-3$ s$-1$]	${\mathcal{L}}_{\mathit{u}∟t}^{\left(t\right)}\left(\zeta {\underline{\breve{\chi }}}_c{vol}_{\mathit{x}}^3\right)$	$=\left[{\displaystyle \frac{\partial \zeta {\underline{\breve{\chi }}}_c}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\zeta {\underline{\breve{\chi }}}_c\mathit{u}\right)\right]dV$	[kgc (m s$-1$)$-3$ s$-1$]
Linear momentum	$\zeta \mathit{u}$	${\displaystyle \frac{\partial (Д\mathit{u})}{\partial t}}={\displaystyle \frac{\partial \left[\right[\left[\zeta \mathit{u}\right]\left]\right]}{\partial t}}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \zeta \mathit{u}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\zeta \mathit{u}\mathit{u}\right)\right]dV$	[(kg m s$-1$) (m s$-1$)$-3$ s$-1$ = N (m s$-1$)$-3$]	${\mathcal{L}}_{\mathit{u}∟t}^{\left(t\right)}\left(\zeta \mathit{u}{vol}_{\mathit{x}}^3\right)$	$=[{\displaystyle \frac{\partial \zeta \mathit{u}}{\partial t}}$$+{\nabla }_{\mathit{x}}·\left(\zeta \mathit{u}\mathit{u}\right)]dV$	[N (m s$-1$)$-3$]
Angular momentum	$\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{\partial \left(Д(\breve{\mathfrak{r}}\times \mathit{u})\right)}{\partial t}}\\ \hfill ={\displaystyle \frac{\partial \left[\left[\left[\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\right]\right]\right]}{\partial t}}\end{array}\right\}$	$=\left\{\begin{array}{c}\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }[{\displaystyle \frac{\partial \zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})}{\partial t}}\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\mathit{u}\right)]dV\hfill \end{array}\right.$	[(kg m${}^2$ s$-1$) (m s$-1$)$-3$ s$-1$ = (N m) (m s$-1$)$-3$]	$\left.\begin{array}{c}\hfill {\mathcal{L}}_{\mathit{u}∟t}^{\left(t\right)}\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\right)\\ \hfill {vol}_{\mathit{x}}^3)\end{array}\right\}$	$=\left\{\begin{array}{c}[{\displaystyle \frac{\partial \zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})}{\partial t}}\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\mathit{u}\right)]dV\hfill \end{array}\right.$	[(N m) (m s$-1$)$-3$]
Energy	$\zeta \underline{\breve{e}}$	${\displaystyle \frac{\partial (Д\breve{e})}{\partial t}}={\displaystyle \frac{\partial \left[\left[\left[\zeta \underline{\breve{e}}\right]\right]\right]}{\partial t}}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \zeta \underline{\breve{e}}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\zeta \underline{\breve{e}}\mathit{u}\right)\right]dV$	[J (m s$-1$)$-3$ s$-1$ = W (m s$-1$)$-3$]	${\mathcal{L}}_{\mathit{u}∟t}^{\left(t\right)}\left(\zeta \underline{\breve{e}}{vol}_{\mathit{x}}^3\right)$	$=\left[{\displaystyle \frac{\partial \zeta \underline{\breve{e}}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\zeta \underline{\breve{e}}\mathit{u}\right)\right]dV$	[J (m s$-1$)$-3$ s$-1$ = W (m s$-1$)$-3$]
Charge (in solution)	$\zeta \underline{\breve{z}}$	${\displaystyle \frac{\partial (Д\breve{z})}{\partial t}}={\displaystyle \frac{\partial \left[\left[\left[\zeta \underline{\breve{z}}\right]\right]\right]}{\partial t}}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \zeta \underline{\breve{z}}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\zeta \underline{\breve{z}}\mathit{u}\right)\right]dV$	[C (m s$-1$)$-3$ s$-1$ = A (m s$-1$)$-3$]	${\mathcal{L}}_{\mathit{u}∟t}^{\left(t\right)}\left(\zeta \underline{\breve{z}}{vol}_{\mathit{x}}^3\right)$	$=\left[{\displaystyle \frac{\partial \zeta \underline{\breve{z}}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\zeta \underline{\breve{z}}\mathit{u}\right)\right]dV$	[C (m s$-1$)$-3$ s$-1$ = A (m s$-1$)$-3$]
Entropy	$\zeta \underline{\breve{s}}$	${\displaystyle \frac{\partial (Д\breve{s})}{\partial t}}={\displaystyle \frac{\partial \left[\left[\left[\zeta \underline{\breve{s}}\right]\right]\right]}{\partial t}}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \zeta \underline{\breve{s}}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\zeta \underline{\breve{s}}\mathit{u}\right)\right]dV$	[J K$-1$ (m s$-1$)$-3$ s$-1$]	${\mathcal{L}}_{\mathit{u}∟t}^{\left(t\right)}\left(\zeta \underline{\breve{s}}{vol}_{\mathit{x}}^3\right)$	$=\left[{\displaystyle \frac{\partial \zeta \underline{\breve{s}}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\zeta \underline{\breve{s}}\mathit{u}\right)\right]dV$	[J K$-1$ (m s$-1$)$-3$ s$-1$]
Probability	$p(\mathit{u},\mathit{x}|t)$$=p_{\mathit{u},\mathit{x}|t}$	${\displaystyle \frac{\partial p_{\mathit{u}|t}}{\partial t}}={\displaystyle \frac{\partial \left[\left[\left[p_{\mathit{u},\mathit{x}|t}\right]\right]\right]}{\partial t}}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial p_{\mathit{u},\mathit{x}|t}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{u},\mathit{x}|t}\mathit{u}\right)\right]dV$	[(m s$-1$)$-3$ s$-1$]	${\mathcal{L}}_{\mathit{u}∟t}^{\left(t\right)}\left(p_{\mathit{u},\mathit{x}|t}{vol}_{\mathit{x}}^3\right)$	$=\left[{\displaystyle \frac{\partial p_{\mathit{u},\mathit{x}|t}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{u},\mathit{x}|t}\mathit{u}\right)\right]dV$	[(m s$-1$)$-3$ s$-1$]
	$p\left(\mathit{x}\right|\mathit{u},t)$$=p_{\mathit{x}|\mathit{u},t}$	$0={\displaystyle \frac{\partial \left[\left[\left[p_{\mathit{x}|\mathit{u},t}\right]\right]\right]}{\partial t}}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial p_{\mathit{x}|\mathit{u},t}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{x}|\mathit{u},t}\mathit{u}\right)\right]dV$	[s$-1$]	0	$={\displaystyle \frac{\partial p_{\mathit{x}|\mathit{u},t}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{x}|\mathit{u},t}\mathit{u}\right)$	[m$-3$ s$-1$]

. We also list two sets of temporal Liouville equations based on $p(\mathit{u},\mathit{x}|t)$ or $p\left(\mathit{x}\right|\mathit{u},t)$.

5.6. Velocimetric-Spatial (Time-Independent) Formulation (Density $\phi$)

Now consider a time-independent velocimetric-spatial formulation of the Reynolds transport theorem, based on the velocity space $\mathcal{D}$ with Cartesian velocity coordinates $\mathit{X}=\mathit{u}$, position parameter vector $\mathit{C}=\mathit{x}$, velocity gradient tensor field $\mathit{V}={(\mathit{G}(\mathit{u},\mathit{x}))}^{\top }:={\left({\nabla }_{\mathit{x}}\mathit{u}\right)}^{\top }$ and generalized density $\psi =\phi (\mathit{u},\mathit{x})$, defined from the velocivolumetric fluid density $\zeta (\mathit{u},\mathit{x})$. This follows integration path ➂ in Figure 5, but is independent of time, representing a stationary flow system (this formulation can also be applied to statistically stationary flow systems—for example a turbulent flow at steady state—by mapping the time dependence to a velocity dependence, thus $\phi (\mathit{u},\mathit{x},t)\mapsto \phi (\mathit{u},\mathit{x},t(\mathit{u}\left)\right)$). For the spatial gradient term in Equation (57), from Equation (43) and relation (32) between $\phi$ and $\alpha$ [23,24]:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill {\nabla }_{\mathit{x}}\alpha & ={\nabla }_{\mathit{x}}\underset{\mathcal{D}\left(\mathit{x}\right)}{\int \int \int }\phi dU=\underset{\mathcal{D}\left(\mathit{x}\right)}{\int \int \int }{\nabla }_{\mathit{x}}\phi dU{+}_{\partial \mathcal{D}\left(\mathit{x}\right)}\phi {\mathit{G}}^{\top }·{\mathit{n}}_{B}dB=\underset{\mathcal{D}\left(\mathit{x}\right)}{\int \int \int }\left[{\nabla }_{\mathit{x}}\phi +{\nabla }_{\mathit{u}}·\left(\phi {\mathit{G}}^{\top }\right)\right]dU\hfill \end{array}\end{array}$$

(61)

This gives an integral equation for the spatial gradient of the volumetric density $\alpha$. We can also write a Lie differential Equation (46), expressed in terms of the multivariate Lie derivative ${\mathcal{L}}_{{\mathit{G}}^{\top }∟\mathit{x}}^{\left(\mathit{x}\right)}$ with respect to the augmented field $\mathit{G}∟\mathit{x}:={\nabla }_{\mathit{x}}[\mathit{u},\mathit{x}]$ over parameters $\mathit{x}$. The conservation laws derived from Equation (61) are listed in

Table 7. Conservation Laws for the Velocimetric-Spatial (Time-Independent) Formulation (based on the Velocivolumetric Fluid Density $\zeta (\mathit{u},\mathit{x})$).

Conserved Quantity	Density	Integral Equation			Differential Equation
	$\phi (\mathit{u},\mathit{x})$$=\zeta \breve{\underline{\phi }}$	LHS = ${\nabla }_{\mathit{x}}\alpha ={\nabla }_{\mathit{x}}\left(\rho \underline{\alpha }\right)$	= RHS	SI Units	LHS	= RHS	SI Units
Fluid mass	$\zeta$	${\nabla }_{\mathit{x}}\rho$$={\nabla }_{\mathit{x}}\langle \langle \langle \zeta \rangle \rangle \rangle$	$=\underset{\mathcal{D}\left(\mathit{x}\right)}{\int \int \int }\left({\nabla }_{\mathit{x}}\zeta +{\nabla }_{\mathit{u}}·\left(\zeta {\mathit{G}}^{\top }\right)\right)dU$	[kg m−4 = kg m$-3$ m$-1$]	${\mathcal{L}}_{{\mathit{G}}^{\top }∟\mathit{x}}^{\left(\mathit{x}\right)}\left(\zeta {vol}_{\mathit{u}}^3\right)$	$=\left({\nabla }_{\mathit{x}}\zeta +{\nabla }_{\mathit{u}}·\left(\zeta {\mathit{G}}^{\top }\right)\right)dU$	[kg m−4 = kg m$-3$ m$-1$]
Species mass	$\zeta {\breve{\underline{\chi }}}_c$$={\zeta }_c$	$\left.\begin{array}{c}\hfill {\nabla }_{\mathit{x}}\left(\rho {\underline{\chi }}_c\right)={\nabla }_{\mathit{x}}\left({\rho }_c\right)\\ \hfill ={\nabla }_{\mathit{x}}\langle \langle \langle \zeta {\breve{\underline{\chi }}}_c\rangle \rangle \rangle \end{array}\right\}$	$=\underset{\mathcal{D}\left(\mathit{x}\right)}{\int \int \int }\left({\nabla }_{\mathit{x}}\left(\zeta {\breve{\underline{\chi }}}_c\right)+{\nabla }_{\mathit{u}}·\left(\zeta {\breve{\underline{\chi }}}_c{\mathit{G}}^{\top }\right)\right)dU$	[kgc m−4 = kgc m$-3$ m$-1$]	${\mathcal{L}}_{{\mathit{G}}^{\top }∟\mathit{x}}^{\left(\mathit{x}\right)}\left(\zeta {\breve{\underline{\chi }}}_c{vol}_{\mathit{u}}^3\right)$	$=\left({\nabla }_{\mathit{x}}\left(\zeta {\breve{\underline{\chi }}}_c\right)+{\nabla }_{\mathit{u}}·\left(\zeta {\breve{\underline{\chi }}}_c{\mathit{G}}^{\top }\right)\right)dU$	[kgc m−4 = kgc m$-3$ m$-1$]
Linear momentum	$\zeta \mathit{u}$	${\nabla }_{\mathit{x}}\left(\rho \mathit{u}\right)$$={\nabla }_{\mathit{x}}\langle \langle \langle \zeta \mathit{u}\rangle \rangle \rangle$	$=\underset{\mathcal{D}\left(\mathit{x}\right)}{\int \int \int }\left({\nabla }_{\mathit{x}}\left(\zeta \mathit{u}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \mathit{u}{\mathit{G}}^{\top }\right)\right)dU$	[kg m$-3$ s$-1$ = (kg m s$-1$) m$-3$ m$-1$]	${\mathcal{L}}_{{\mathit{G}}^{\top }∟\mathit{x}}^{\left(\mathit{x}\right)}\left(\zeta \mathit{u}{vol}_{\mathit{u}}^3\right)$	$=\left({\nabla }_{\mathit{x}}\left(\zeta \mathit{u}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \mathit{u}{\mathit{G}}^{\top }\right)\right)dU$	[kg m$-3$ s$-1$]
Angular momentum	$\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})$	$\left.\begin{array}{c}\hfill {\nabla }_{\mathit{x}}\left(\rho (\underline{\mathfrak{r}}\times \mathit{u})\right)\\ \hfill ={\nabla }_{\mathit{x}}\langle \langle \langle \zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\rangle \rangle \rangle \end{array}\right\}$	$=\left\{\begin{array}{c}\underset{\mathcal{D}\left(\mathit{x}\right)}{\int \int \int }({\nabla }_{\mathit{x}}\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\right)\hfill \\ +{\nabla }_{\mathit{u}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u}){\mathit{G}}^{\top }\right))dU\hfill \end{array}\right.$	[kg m$-2$ s$-1$ = (kg m2 s$-1$) m$-3$ m$-1$]	$\left.\begin{array}{c}\hfill {\mathcal{L}}_{{\mathit{G}}^{\top }∟\mathit{x}}^{\left(\mathit{x}\right)}(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\\ \hfill {vol}_{\mathit{u}}^3)\end{array}\right\}$	$=\left\{\begin{array}{c}({\nabla }_{\mathit{x}}\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\right)\hfill \\ +{\nabla }_{\mathit{u}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u}){\mathit{G}}^{\top }\right))dU\hfill \end{array}\right.$	[kg m$-2$ s$-1$]
Energy	$\zeta \breve{\underline{e}}$	${\nabla }_{\mathit{x}}\left(\rho \underline{e}\right)$$={\nabla }_{\mathit{x}}\langle \langle \langle \zeta \breve{\underline{e}}\rangle \rangle \rangle$	$=\underset{\mathcal{D}\left(\mathit{x}\right)}{\int \int \int }\left({\nabla }_{\mathit{x}}\left(\zeta \breve{\underline{e}}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \breve{\underline{e}}{\mathit{G}}^{\top }\right)\right)dU$	[J m$-3$ m$-1$]	${\mathcal{L}}_{{\mathit{G}}^{\top }∟\mathit{x}}^{\left(\mathit{x}\right)}\left(\zeta \breve{\underline{e}}{vol}_{\mathit{u}}^3\right)$	$=\left({\nabla }_{\mathit{x}}\left(\zeta \breve{\underline{e}}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \breve{\underline{e}}{\mathit{G}}^{\top }\right)\right)dU$	[J m$-3$ m$-1$]
Charge (in solution)	$\zeta \breve{\underline{z}}$	${\nabla }_{\mathit{x}}\left(\rho \underline{z}\right)$$={\nabla }_{\mathit{x}}\langle \langle \langle \zeta \breve{\underline{z}}\rangle \rangle \rangle$	$=\underset{\mathcal{D}\left(\mathit{x}\right)}{\int \int \int }\left({\nabla }_{\mathit{x}}\left(\zeta \breve{\underline{z}}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \breve{\underline{z}}{\mathit{G}}^{\top }\right)\right)dU$	[C m$-3$ m$-1$]	${\mathcal{L}}_{{\mathit{G}}^{\top }∟\mathit{x}}^{\left(\mathit{x}\right)}\left(\zeta \breve{\underline{z}}{vol}_{\mathit{u}}^3\right)$	$=\left({\nabla }_{\mathit{x}}\left(\zeta \breve{\underline{z}}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \breve{\underline{z}}{\mathit{G}}^{\top }\right)\right)dU$	[C m$-3$ m$-1$]
Entropy	$\zeta \breve{\underline{s}}$	${\nabla }_{\mathit{x}}\left(\rho \underline{s}\right)$$={\nabla }_{\mathit{x}}\langle \langle \langle \zeta \breve{\underline{s}}\rangle \rangle \rangle$	$=\underset{\mathcal{D}\left(\mathit{x}\right)}{\int \int \int }\left({\nabla }_{\mathit{x}}\left(\zeta \breve{\underline{s}}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \breve{\underline{s}}{\mathit{G}}^{\top }\right)\right)dU$	[J K$-1$ m$-3$ m$-1$]	${\mathcal{L}}_{{\mathit{G}}^{\top }∟\mathit{x}}^{\left(\mathit{x}\right)}\left(\zeta \breve{\underline{s}}{vol}_{\mathit{u}}^3\right)$	$=\left({\nabla }_{\mathit{x}}\left(\zeta \breve{\underline{s}}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \breve{\underline{s}}{\mathit{G}}^{\top }\right)\right)dU$	[J K$-1$ m$-3$ m$-1$]
Probability	$p(\mathit{u},\mathit{x})$$=p_{\mathit{u},\mathit{x}}$	${\nabla }_{\mathit{x}}{p}_{\mathit{x}}$$={\nabla }_{\mathit{x}}\langle \langle \langle p_{\mathit{u},\mathit{x}}\rangle \rangle \rangle$	$=\underset{\mathcal{D}\left(\mathit{x}\right)}{\int \int \int }\left({\nabla }_{\mathit{x}}{p}_{\mathit{u},\mathit{x}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u},\mathit{x}}{\mathit{G}}^{\top }\right)\right)dU$	[m$-3$ m$-1$]	${\mathcal{L}}_{{\mathit{G}}^{\top }∟\mathit{x}}^{\left(\mathit{x}\right)}\left(p_{\mathit{u},\mathit{x}}{vol}_{\mathit{u}}^3\right)$	$=\left({\nabla }_{\mathit{x}}{p}_{\mathit{u},\mathit{x}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u},\mathit{x}}{\mathit{G}}^{\top }\right)\right)dU$	[m$-3$ m$-1$]
	$p\left(\mathit{u}\right|\mathit{x})$$=p_{\mathit{u}|\mathit{x}}$	$\mathbf{0}={\nabla }_{\mathit{x}}\langle \langle \langle p_{\mathit{u}|\mathit{x}}\rangle \rangle \rangle$	$=\underset{\mathcal{D}\left(\mathit{x}\right)}{\int \int \int }\left({\nabla }_{\mathit{x}}{p}_{\mathit{u}|\mathit{x}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u}|\mathit{x}}{\mathit{G}}^{\top }\right)\right)dU$	[m$-1$]	$\mathbf{0}$	$={\nabla }_{\mathit{x}}{p}_{\mathit{u}|\mathit{x}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u}|\mathit{x}}{\mathit{G}}^{\top }\right)$	[(m s$-1$)$-3$ m$-1$]

. We also list two sets of spatial Liouville equations based on $p(\mathit{u},\mathit{x})$ or $p\left(\mathit{u}\right|\mathit{x})$ [23].

5.7. Volumetric-Velocital (Time-Independent) Formulation (Density $\phi$)

Now consider a time-independent volumetric-velocital formulation of the Reynolds transport theorem, based on the geometric space $\mathsf{\Omega }$ with Cartesian coordinates $\mathit{X}=\mathit{x}$, velocity parameter vector $\mathit{C}=\mathit{u}$, inverse velocity gradient tensor field $\mathit{V}={\left(\mathsf{\Gamma }(\mathit{u},\mathit{x})\right)}^{\top }:={\left({\nabla }_{\mathit{u}}\mathit{x}\right)}^{\top }$ and generalized density $\psi =\phi (\mathit{u},\mathit{x})$ defined from $\zeta (\mathit{u},\mathit{x})$. (The descriptor velocital, from Latin velocitas and -al for “pertaining to velocity”, follows [80].) This formulation follows integration path ➃ in Figure 5; recall that this adopts the velocimetric representation, which is integrated over the volumetric space for a distinct velocity. We again consider the time-independent case. From Equations (43), (59) and relation (33) between $\phi$ and $\beta$:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill {\nabla }_{\mathit{u}}\beta & ={\nabla }_{\mathit{u}}\underset{\mathsf{\Omega }\left(\mathit{u}\right)}{\int \int \int }\phi dV=\underset{\mathsf{\Omega }\left(\mathit{u}\right)}{\int \int \int }{\nabla }_{\mathit{u}}\phi dV{+}_{\partial \mathsf{\Omega }\left(\mathit{u}\right)}\phi {\mathsf{\Gamma }}^{\top }·\mathit{n}dA=\underset{\mathsf{\Omega }\left(\mathit{u}\right)}{\int \int \int }\left[{\nabla }_{\mathit{u}}\phi +{\nabla }_{\mathit{x}}·\left(\phi {\mathsf{\Gamma }}^{\top }\right)\right]dV\hfill \end{array}\end{array}$$

(62)

This gives the gradient in velocity space of the velocimetric density $\beta$. The conservation laws derived from Equation (62) are listed in

Table 8. Conservation Laws for the Volumetric-Velocital (Time-Independent) Formulation (based on the Velocivolumetric Fluid Density $\zeta (\mathit{u},\mathit{x})$).

Conserved Quantity	Density	Integral Equation			Differential Equation
	$\phi (\mathit{u},\mathit{x})$$=\zeta \breve{\underline{\phi }}$	LHS = ${\nabla }_{\mathit{u}}\beta$ $={\nabla }_{\mathit{u}}(Д\breve{\beta })$	= RHS	SI Units	LHS	= RHS	SI Units
Fluid mass	$\zeta$	${\nabla }_{\mathit{u}}Д$$={\nabla }_{\mathit{u}}\left[\left[\left[\zeta \right]\right]\right]$	$=\underset{\mathsf{\Omega }\left(\mathit{u}\right)}{\int \int \int }\left({\nabla }_{\mathit{u}}\zeta +{\nabla }_{\mathit{x}}·\left(\zeta {\mathsf{\Gamma }}^{\top }\right)\right)dV$	[kg (m s$-1$)−4 = kg (m s$-1$)$-3$ (m s$-1$)$-1$]	${\mathcal{L}}_{{\mathsf{\Gamma }}^{\top }∟\mathit{u}}^{\left(\mathit{u}\right)}\left(\zeta {vol}_{\mathit{x}}^3\right)$	$=\left({\nabla }_{\mathit{u}}\zeta +{\nabla }_{\mathit{x}}·\left(\zeta {\mathsf{\Gamma }}^{\top }\right)\right)dV$	[kg (m s$-1$)−4 = kg (m s$-1$)$-3$ (m s$-1$)$-1$]
Species mass	$\zeta {\breve{\underline{\chi }}}_c$$={\zeta }_c$	$\left.\begin{array}{c}\hfill {\nabla }_{\mathit{u}}(Д{\breve{\chi }}_c)={\nabla }_{\mathit{u}}\left({Д}_c\right)\\ \hfill ={\nabla }_{\mathit{u}}\left[\left[\left[\zeta {\breve{\underline{\chi }}}_c\right]\right]\right]\end{array}\right\}$	$=\underset{\mathsf{\Omega }\left(\mathit{u}\right)}{\int \int \int }\left({\nabla }_{\mathit{u}}\left(\zeta {\breve{\underline{\chi }}}_c\right)+{\nabla }_{\mathit{x}}·\left(\zeta {\breve{\underline{\chi }}}_c{\mathsf{\Gamma }}^{\top }\right)\right)dV$	[kgc (m s$-1$)−4 = kgc (m s$-1$)$-3$ (m s$-1$)$-1$]	${\mathcal{L}}_{{\mathsf{\Gamma }}^{\top }∟\mathit{u}}^{\left(\mathit{u}\right)}\left(\zeta {\breve{\underline{\chi }}}_c{vol}_{\mathit{x}}^3\right)$	$=\left({\nabla }_{\mathit{u}}\left(\zeta {\breve{\underline{\chi }}}_c\right)+{\nabla }_{\mathit{x}}·\left(\zeta {\breve{\underline{\chi }}}_c{\mathsf{\Gamma }}^{\top }\right)\right)dV$	[kgc (m s$-1$)−4 = kgc (m s$-1$)$-3$ (m s$-1$)$-1$]
Linear momentum	$\zeta \mathit{u}$	${\nabla }_{\mathit{u}}(Д\mathit{u})$$={\nabla }_{\mathit{u}}\left[\left[\left[\zeta \mathit{u}\right]\right]\right]$	$=\underset{\mathsf{\Omega }\left(\mathit{u}\right)}{\int \int \int }\left({\nabla }_{\mathit{u}}\left(\zeta \mathit{u}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \mathit{u}{\mathsf{\Gamma }}^{\top }\right)\right)dV$	[kg (m s$-1$)$-3$ = (kg m s$-1$) (m s$-1$)$-3$ (m s$-1$)$-1$]	${\mathcal{L}}_{{\mathsf{\Gamma }}^{\top }∟\mathit{u}}^{\left(\mathit{u}\right)}\left(\zeta \mathit{u}{vol}_{\mathit{x}}^3\right)$	$=\left({\nabla }_{\mathit{u}}\left(\zeta \mathit{u}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \mathit{u}{\mathsf{\Gamma }}^{\top }\right)\right)dV$	[kg (m s$-1$)$-3$]
Angular momentum	$\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})$	$\left.\begin{array}{c}\hfill {\nabla }_{\mathit{u}}\left(Д(\breve{\mathfrak{r}}\times \mathit{u})\right)\\ \hfill ={\nabla }_{\mathit{u}}\left[\left[\left[\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\right]\right]\right]\end{array}\right\}$	$=\left\{\begin{array}{c}\underset{\mathsf{\Omega }\left(\mathit{u}\right)}{\int \int \int }({\nabla }_{\mathit{u}}\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u}){\mathsf{\Gamma }}^{\top }\right))dV\hfill \end{array}\right.$	[kg m (m s$-1$)$-3$ = (kg m2 s$-1$) (m s$-1$)$-3$ (m s$-1$)$-1$]	$\left.\begin{array}{c}\hfill {\mathcal{L}}_{{\mathsf{\Gamma }}^{\top }∟\mathit{u}}^{\left(\mathit{u}\right)}(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\\ \hfill {vol}_{\mathit{x}}^3)\end{array}\right\}$	$=\left\{\begin{array}{c}({\nabla }_{\mathit{u}}\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u}){\mathsf{\Gamma }}^{\top }\right))dV\hfill \end{array}\right.$	[kg m (m s$-1$)$-3$]
Energy	$\zeta \breve{\underline{e}}$	${\nabla }_{\mathit{u}}(Д\breve{e})$$={\nabla }_{\mathit{u}}\left[\left[\left[\zeta \breve{\underline{e}}\right]\right]\right]$	$=\underset{\mathsf{\Omega }\left(\mathit{u}\right)}{\int \int \int }\left({\nabla }_{\mathit{u}}\left(\zeta \breve{\underline{e}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \breve{\underline{e}}{\mathsf{\Gamma }}^{\top }\right)\right)dV$	[J (m s$-1$)$-3$ (m s$-1$)$-1$]	${\mathcal{L}}_{{\mathsf{\Gamma }}^{\top }∟\mathit{u}}^{\left(\mathit{u}\right)}\left(\zeta \breve{\underline{e}}{vol}_{\mathit{x}}^3\right)$	$=\left({\nabla }_{\mathit{u}}\left(\zeta \breve{\underline{e}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \breve{\underline{e}}{\mathsf{\Gamma }}^{\top }\right)\right)dV$	[J (m s$-1$)$-3$ (m s$-1$)$-1$]
Charge (in solution)	$\zeta \breve{\underline{z}}$	${\nabla }_{\mathit{u}}(Д\breve{z})$$={\nabla }_{\mathit{u}}\left[\left[\left[\zeta \breve{\underline{z}}\right]\right]\right]$	$=\underset{\mathsf{\Omega }\left(\mathit{u}\right)}{\int \int \int }\left({\nabla }_{\mathit{u}}\left(\zeta \breve{\underline{z}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \breve{\underline{z}}{\mathsf{\Gamma }}^{\top }\right)\right)dV$	[C (m s$-1$)$-3$ (m s$-1$)$-1$]	${\mathcal{L}}_{{\mathsf{\Gamma }}^{\top }∟\mathit{u}}^{\left(\mathit{u}\right)}\left(\zeta \breve{\underline{z}}{vol}_{\mathit{x}}^3\right)$	$=\left({\nabla }_{\mathit{u}}\left(\zeta \breve{\underline{z}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \breve{\underline{z}}{\mathsf{\Gamma }}^{\top }\right)\right)dV$	[C (m s$-1$)$-3$ (m s$-1$)$-1$]
Entropy	$\zeta \breve{\underline{s}}$	${\nabla }_{\mathit{u}}(Д\breve{s})$$={\nabla }_{\mathit{u}}\left[\left[\left[\zeta \breve{\underline{s}}\right]\right]\right]$	$=\underset{\mathsf{\Omega }\left(\mathit{u}\right)}{\int \int \int }\left({\nabla }_{\mathit{u}}\left(\zeta \breve{\underline{s}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \breve{\underline{s}}{\mathsf{\Gamma }}^{\top }\right)\right)dV$	[J K$-1$ (m s$-1$)$-3$ (m s$-1$)$-1$]	${\mathcal{L}}_{{\mathsf{\Gamma }}^{\top }∟\mathit{u}}^{\left(\mathit{u}\right)}\left(\zeta \breve{\underline{s}}{vol}_{\mathit{x}}^3\right)$	$=\left({\nabla }_{\mathit{u}}\left(\zeta \breve{\underline{s}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \breve{\underline{s}}{\mathsf{\Gamma }}^{\top }\right)\right)dV$	[J K$-1$ (m s$-1$)$-3$ (m s$-1$)$-1$]
Probability	$p(\mathit{u},\mathit{x})$$=p_{\mathit{u},\mathit{x}}$	${\nabla }_{\mathit{u}}{p}_{\mathit{u}}$$={\nabla }_{\mathit{u}}\left[\left[\left[p_{\mathit{u},\mathit{x}}\right]\right]\right]$	$=\underset{\mathsf{\Omega }\left(\mathit{u}\right)}{\int \int \int }\left({\nabla }_{\mathit{u}}{p}_{\mathit{u},\mathit{x}}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{u},\mathit{x}}{\mathsf{\Gamma }}^{\top }\right)\right)dV$	[(m s$-1$)$-3$ (m s$-1$)$-1$]	${\mathcal{L}}_{{\mathsf{\Gamma }}^{\top }∟\mathit{u}}^{\left(\mathit{u}\right)}\left(p_{\mathit{u},\mathit{x}}{vol}_{\mathit{x}}^3\right)$	$=\left({\nabla }_{\mathit{u}}{p}_{\mathit{u},\mathit{x}}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{u},\mathit{x}}{\mathsf{\Gamma }}^{\top }\right)\right)dV$	[(m s$-1$)$-3$ (m s$-1$)$-1$]
	$p\left(\mathit{x}\right|\mathit{u})$$=p_{\mathit{x}|\mathit{u}}$	$\mathbf{0}={\nabla }_{\mathit{u}}\left[\left[\left[p_{\mathit{x}|\mathit{u}}\right]\right]\right]$	$=\underset{\mathsf{\Omega }\left(\mathit{u}\right)}{\int \int \int }\left({\nabla }_{\mathit{u}}{p}_{\mathit{x}|\mathit{u}}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{x}|\mathit{u}}{\mathsf{\Gamma }}^{\top }\right)\right)dV$	[(m s$-1$)$-1$]	$\mathbf{0}$	$={\nabla }_{\mathit{u}}{p}_{\mathit{x}|\mathit{u}}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{x}|\mathit{u}}{\mathsf{\Gamma }}^{\top }\right)$	[m$-3$ (m s$-1$)$-1$]

. We also list two sets of spatial Liouville equations based on $p(\mathit{u},\mathit{x})$ or $p\left(\mathit{x}\right|\mathit{u})$.

5.8. Velocimetric-Spatiotemporal Formulation (Density $\phi$)

Now consider the complete velocimetric-spatiotemporal formulation of the Reynolds transport theorem, a time-dependent extension of Section 5.6 (path ➂ in Figure 5). This examines the velocity space $\mathcal{D}$ with $\mathit{X}=\mathit{u}$, $\mathit{C}={[\mathit{x},t]}^{\top }$, $\mathit{V}={\left(\tilde{\mathit{G}}(\mathit{u},\mathit{x},t)\right)}^{\top }:={\left({▾}_{\mathit{x}}\mathit{u}\right)}^{\top }$, where ${▾}_{\mathit{x}}={\nabla }_{\mathit{x},t}$ is the spatiotemporal gradient, and density $\psi =\phi (\mathit{u},\mathit{x},t)$ defined from $\zeta (\mathit{u},\mathit{x},t)$. Here $\tilde{\mathit{G}}={[\mathit{G},\dot{\mathit{u}}]}^{\top }$ can be recognized as the velocity gradient tensor field $\mathit{G}$ augmented with the local acceleration vector field $\dot{\mathit{u}}:=\partial \mathit{u}/\partial t$. From Equations (43 and (32 [23]:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill {▾}_{\mathit{x}}\alpha & ={▾}_{\mathit{x}}\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\phi dU=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }{▾}_{\mathit{x}}\phi dU{+}_{\partial \mathcal{D}(\mathit{x},t)}\phi {\tilde{\mathit{G}}}^{\top }·{\mathit{n}}_{B}dB=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{▾}_{\mathit{x}}\phi +{\nabla }_{\mathit{u}}·\left(\phi {\tilde{\mathit{G}}}^{\top }\right)\right]dU\hfill \end{array}\end{array}$$

(63)

This gives the spatiotemporal gradient ${▾}_{\mathit{x}}\alpha ={[{\nabla }_{\mathit{x}}\alpha ,\dot{\alpha }]}^{\top }$ of the volumetric density $\alpha$. The conservation laws derived from Equation (63) are listed in

Table 9. Conservation Laws for the Velocimetric-Spatiotemporal Formulation (based on the Velocivolumetric Fluid Density $\zeta (\mathit{u},\mathit{x},t)$).

Conserved Quantity	Density	Integral Equation			Differential Equation
	$\phi (\mathit{u},\mathit{x},t)$$=\zeta \breve{\underline{\phi }}$	LHS = ${▾}_{\mathit{x}}\alpha$ $={▾}_{\mathit{x}}\left(\rho \underline{\alpha }\right)$	= RHS	SI Units	LHS	= RHS	SI Units
Fluid mass	$\zeta$	${▾}_{\mathit{x}}\rho$$={▾}_{\mathit{x}}\langle \langle \langle \zeta \rangle \rangle \rangle$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left({▾}_{\mathit{x}}\zeta +{\nabla }_{\mathit{u}}·\left(\zeta {\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[kg m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]	${\mathcal{L}}_{{\tilde{\mathit{G}}}^{\top }∟\mathit{x},t}^{(\mathit{x},t)}\left(\zeta {vol}_{\mathit{u}}^3\right)$	$=\left({▾}_{\mathit{x}}\zeta +{\nabla }_{\mathit{u}}·\left(\zeta {\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[kg m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]
Species mass	$\zeta {\breve{\underline{\chi }}}_c$$={\zeta }_c$	$\left.\begin{array}{c}\hfill {▾}_{\mathit{x}}\left(\rho {\underline{\chi }}_c\right)={▾}_{\mathit{x}}\left({\rho }_c\right)\\ \hfill ={▾}_{\mathit{x}}\langle \langle \langle \zeta {\breve{\underline{\chi }}}_c\rangle \rangle \rangle \end{array}\right\}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left({▾}_{\mathit{x}}\left(\zeta {\breve{\underline{\chi }}}_c\right)+{\nabla }_{\mathit{u}}·\left(\zeta {\breve{\underline{\chi }}}_c{\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[kgc m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]	${\mathcal{L}}_{{\tilde{\mathit{G}}}^{\top }∟\mathit{x},t}^{(\mathit{x},t)}\left(\zeta {\breve{\underline{\chi }}}_c{vol}_{\mathit{u}}^3\right)$	$=\left({▾}_{\mathit{x}}\left(\zeta {\breve{\underline{\chi }}}_c\right)+{\nabla }_{\mathit{u}}·\left(\zeta {\breve{\underline{\chi }}}_c{\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[kgc m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]
Linear momentum	$\zeta \mathit{u}$	${▾}_{\mathit{x}}\left(\rho \mathit{u}\right)$$={▾}_{\mathit{x}}\langle \langle \langle \zeta \mathit{u}\rangle \rangle \rangle$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left({▾}_{\mathit{x}}\left(\zeta \mathit{u}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \mathit{u}{\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[(kg m s$-1$) m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]	${\mathcal{L}}_{{\tilde{\mathit{G}}}^{\top }∟\mathit{x},t}^{(\mathit{x},t)}\left(\zeta \mathit{u}{vol}_{\mathit{u}}^3\right)$	$=\left({▾}_{\mathit{x}}\left(\zeta \mathit{u}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \mathit{u}{\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[(kg m s$-1$) m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]
Angular momentum	$\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})$	$\left.\begin{array}{c}\hfill {▾}_{\mathit{x}}\left(\rho (\underline{\mathfrak{r}}\times \mathit{u})\right)\\ \hfill ={▾}_{\mathit{x}}\langle \langle \langle \zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\rangle \rangle \rangle \end{array}\right\}$	$=\left\{\begin{array}{c}\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }({▾}_{\mathit{x}}\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\right)\hfill \\ +{\nabla }_{\mathit{u}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u}){\tilde{\mathit{G}}}^{\top }\right))dU\hfill \end{array}\right.$	[(kg m2 s$-1$) m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]	$\left.\begin{array}{c}\hfill {\mathcal{L}}_{{\tilde{\mathit{G}}}^{\top }∟\mathit{x},t}^{(\mathit{x},t)}(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\\ \hfill {vol}_{\mathit{u}}^3)\end{array}\right\}$	$=\left\{\begin{array}{c}({▾}_{\mathit{x}}\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\right)\hfill \\ +{\nabla }_{\mathit{u}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u}){\tilde{\mathit{G}}}^{\top }\right))dU\hfill \end{array}\right.$	[(kg m2 s$-1$) m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]
Energy	$\zeta \breve{\underline{e}}$	${▾}_{\mathit{x}}\left(\rho \underline{e}\right)$$={▾}_{\mathit{x}}\langle \langle \langle \zeta \breve{\underline{e}}\rangle \rangle \rangle$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left({▾}_{\mathit{x}}\left(\zeta \breve{\underline{e}}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \breve{\underline{e}}{\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[J m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]	${\mathcal{L}}_{{\tilde{\mathit{G}}}^{\top }∟\mathit{x},t}^{(\mathit{x},t)}\left(\zeta \breve{\underline{e}}{vol}_{\mathit{u}}^3\right)$	$=\left({▾}_{\mathit{x}}\left(\zeta \breve{\underline{e}}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \breve{\underline{e}}{\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[J m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]
Charge (in solution)	$\zeta \breve{\underline{z}}$	${▾}_{\mathit{x}}\left(\rho \underline{z}\right)$$={▾}_{\mathit{x}}\langle \langle \langle \zeta \breve{\underline{z}}\rangle \rangle \rangle$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left({▾}_{\mathit{x}}\left(\zeta \breve{\underline{z}}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \breve{\underline{z}}{\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[C m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]	${\mathcal{L}}_{{\tilde{\mathit{G}}}^{\top }∟\mathit{x},t}^{(\mathit{x},t)}\left(\zeta \breve{\underline{z}}{vol}_{\mathit{u}}^3\right)$	$=\left({▾}_{\mathit{x}}\left(\zeta \breve{\underline{z}}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \breve{\underline{z}}{\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[C m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]
Entropy	$\zeta \breve{\underline{s}}$	${▾}_{\mathit{x}}\left(\rho \underline{s}\right)$$={▾}_{\mathit{x}}\langle \langle \langle \zeta \breve{\underline{s}}\rangle \rangle \rangle$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left({▾}_{\mathit{x}}\left(\zeta \breve{\underline{s}}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \breve{\underline{s}}{\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[J K$-1$ m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]	${\mathcal{L}}_{{\tilde{\mathit{G}}}^{\top }∟\mathit{x},t}^{(\mathit{x},t)}\left(\zeta \breve{\underline{s}}{vol}_{\mathit{u}}^3\right)$	$=\left({▾}_{\mathit{x}}\left(\zeta \breve{\underline{s}}\right)+{\nabla }_{\mathit{u}}·\left(\zeta \breve{\underline{s}}{\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[J K$-1$ m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]
Probability	$p(\mathit{u},\mathit{x}|t)$$=p_{\mathit{u},\mathit{x}|t}$	${▾}_{\mathit{x}}{p}_{\mathit{x}|t}$$={▾}_{\mathit{x}}\langle \langle \langle p_{\mathit{u},\mathit{x}|t}\rangle \rangle \rangle$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left({▾}_{\mathit{x}}{p}_{\mathit{u},\mathit{x}|t}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u},\mathit{x}|t}{\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]	$\left.\begin{array}{c}\hfill {\mathcal{L}}_{{\tilde{\mathit{G}}}^{\top }∟\mathit{x},t}^{(\mathit{x},t)}(p_{\mathit{u},\mathit{x}|t}\\ \hfill {vol}_{\mathit{u}}^3)\end{array}\right\}$	$=\left\{\begin{array}{c}({▾}_{\mathit{x}}{p}_{\mathit{u},\mathit{x}|t}\hfill \\ +{\nabla }_{\mathit{u}}·\left(p_{\mathit{u},\mathit{x}|t}{\tilde{\mathit{G}}}^{\top }\right))dU\hfill \end{array}\right.$	[m$-3$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]
	$p\left(\mathit{u}\right|\mathit{x},t)$$=p_{\mathit{u}|\mathit{x},t}$	$\mathbf{0}={▾}_{\mathit{x}}\langle \langle \langle p_{\mathit{u}|\mathit{x},t}\rangle \rangle \rangle$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left({▾}_{\mathit{x}}{p}_{\mathit{u}|\mathit{x},t}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u}|\mathit{x},t}{\tilde{\mathit{G}}}^{\top }\right)\right)dU$	[${\mathit{d}}_{\mathit{x},t}^{-1}$]	$\mathbf{0}$	$={▾}_{\mathit{x}}{p}_{\mathit{u}|\mathit{x},t}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u}|\mathit{x},t}{\tilde{\mathit{G}}}^{\top }\right)$	[(m s$-1$)$-1$ ${\mathit{d}}_{\mathit{x},t}^{-1}$]

. We also present two sets of composite Liouville equations based on $p(\mathit{u},\mathit{x}|t)$ or $p\left(\mathit{u}\right|\mathit{x},t)$.

5.9. Volumetric-Velocitemporal Formulation (Density $\phi$)

Now consider the complete volumetric-velocitemporal formulation of the Reynolds transport theorem, a time-dependent extension of that in Section 5.7 (path ➃ in Figure 5). This examines the geometric space $\mathsf{\Omega }$ with $\mathit{X}=\mathit{x}$, $\mathit{C}={[\mathit{u},t]}^{\top }$, $\mathit{V}={\left(\tilde{\mathsf{\Gamma }}(\mathit{u},\mathit{x},t)\right)}^{\top }:={\left({▾}_{\mathit{u}}\mathit{x}\right)}^{\top }$, where ${▾}_{\mathit{u}}={\nabla }_{\mathit{u},t}$ is a velocitemporal gradient operator, and density $\psi =\phi (\mathit{u},\mathit{x},t)$ defined from $\zeta (\mathit{u},\mathit{x},t)$. We recognize $\tilde{\mathsf{\Gamma }}={[\mathsf{\Gamma },\mathit{u}]}^{\top }$ as the inverse velocity gradient tensor field $\mathsf{\Gamma }$ augmented with the local velocity vector field $\mathit{u}:=\partial \mathit{x}/\partial t$. From Equations (43 and (33):

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill {▾}_{\mathit{u}}\beta & ={▾}_{\mathit{u}}\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\phi dV=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }{▾}_{\mathit{u}}\phi dV{+}_{\partial \mathsf{\Omega }(\mathit{u},t)}\phi {\tilde{\mathsf{\Gamma }}}^{\top }·\mathit{n}dA=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{▾}_{\mathit{u}}\phi +{\nabla }_{\mathit{x}}·\left(\phi {\tilde{\mathsf{\Gamma }}}^{\top }\right)\right]dV\hfill \end{array}\end{array}$$

(64)

This provides the velocitemporal gradient ${▾}_{\mathit{u}}\beta ={[{\nabla }_{\mathit{u}}\beta ,\dot{\beta }]}^{\top }$ of the velocimetric density $\beta$. The conservation laws derived from Equation (64) are listed in

Table 10. Conservation Laws for the Volumetric-Velocitotemporal Formulation (based on the Velocivolumetric Fluid Density $\zeta (\mathit{u},\mathit{x},t)$).

Conserved quantity	Density	Integral Equation			Differential Equation
	$\phi (\mathit{u},\mathit{x},t)$$=\zeta \breve{\underline{\phi }}$	LHS = ${▾}_{\mathit{u}}\beta$ $={▾}_{\mathit{u}}(Д\breve{\beta })$	= RHS	SI Units	LHS	= RHS	SI Units
Fluid mass	$\zeta$	${▾}_{\mathit{u}}Д$$={▾}_{\mathit{u}}\left[\left[\left[\zeta \right]\right]\right]$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left({▾}_{\mathit{u}}\zeta +{\nabla }_{\mathit{x}}·\left(\zeta {\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[kg (m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]	${\mathcal{L}}_{{\tilde{\mathsf{\Gamma }}}^{\top }∟\mathit{u},t}^{(\mathit{u},t)}\left(\zeta {vol}_{\mathit{x}}^3\right)$	$=\left({▾}_{\mathit{u}}\zeta +{\nabla }_{\mathit{x}}·\left(\zeta {\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[kg (m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]
Species mass	$\zeta {\breve{\underline{\chi }}}_c$$={\zeta }_c$	$\left.\begin{array}{c}\hfill {▾}_{\mathit{u}}(Д{\breve{\chi }}_c)={▾}_{\mathit{u}}\left({Д}_c\right)\\ \hfill ={▾}_{\mathit{u}}\left[\left[\left[\zeta {\breve{\underline{\chi }}}_c\right]\right]\right]\end{array}\right\}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left({▾}_{\mathit{u}}\left(\zeta {\breve{\underline{\chi }}}_c\right)+{\nabla }_{\mathit{x}}·\left(\zeta {\breve{\underline{\chi }}}_c{\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[kgc (m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]	${\mathcal{L}}_{{\tilde{\mathsf{\Gamma }}}^{\top }∟\mathit{u},t}^{(\mathit{u},t)}\left(\zeta {\breve{\underline{\chi }}}_c{vol}_{\mathit{x}}^3\right)$	$=\left({▾}_{\mathit{u}}\left(\zeta {\breve{\underline{\chi }}}_c\right)+{\nabla }_{\mathit{x}}·\left(\zeta {\breve{\underline{\chi }}}_c{\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[kgc (m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]
Linear momentum	$\zeta \mathit{u}$	${▾}_{\mathit{u}}(Д\mathit{u})$$={▾}_{\mathit{u}}\left[\left[\left[\zeta \mathit{u}\right]\right]\right]$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left({▾}_{\mathit{u}}\left(\zeta \mathit{u}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \mathit{u}{\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[(kg m s$-1$) (m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]	${\mathcal{L}}_{{\tilde{\mathsf{\Gamma }}}^{\top }∟\mathit{u},t}^{(\mathit{u},t)}\left(\zeta \mathit{u}{vol}_{\mathit{x}}^3\right)$	$=\left({▾}_{\mathit{u}}\left(\zeta \mathit{u}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \mathit{u}{\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[(kg m s$-1$) (m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]
Angular momentum	$\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})$	$\left.\begin{array}{c}\hfill {▾}_{\mathit{u}}\left(Д(\breve{\mathfrak{r}}\times \mathit{u})\right)\\ \hfill ={▾}_{\mathit{u}}\left[\left[\left[\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\right]\right]\right]\end{array}\right\}$	$=\left\{\begin{array}{c}\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }({▾}_{\mathit{u}}\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u}){\tilde{\mathsf{\Gamma }}}^{\top }\right))dV\hfill \end{array}\right.$	[(kg m2 s$-1$) (m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]	$\left.\begin{array}{c}\hfill {\mathcal{L}}_{{\tilde{\mathsf{\Gamma }}}^{\top }∟\mathit{u},t}^{(\mathit{u},t)}(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\\ \hfill {vol}_{\mathit{x}}^3)\end{array}\right\}$	$=\left\{\begin{array}{c}({▾}_{\mathit{u}}\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u})\right)\hfill \\ +{\nabla }_{\mathit{x}}·\left(\zeta (\underline{\breve{\mathfrak{r}}}\times \mathit{u}){\tilde{\mathsf{\Gamma }}}^{\top }\right))dV\hfill \end{array}\right.$	[(kg m2 s$-1$) (m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]
Energy	$\zeta \breve{\underline{e}}$	${▾}_{\mathit{u}}(Д\breve{e})$$={▾}_{\mathit{u}}\left[\left[\left[\zeta \breve{\underline{e}}\right]\right]\right]$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left({▾}_{\mathit{u}}\left(\zeta \breve{\underline{e}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \breve{\underline{e}}{\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[J (m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]	${\mathcal{L}}_{{\tilde{\mathsf{\Gamma }}}^{\top }∟\mathit{u},t}^{(\mathit{u},t)}\left(\zeta \breve{\underline{e}}{vol}_{\mathit{x}}^3\right)$	$=\left({▾}_{\mathit{u}}\left(\zeta \breve{\underline{e}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \breve{\underline{e}}{\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[J (m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]
Charge (in solution)	$\zeta \breve{\underline{z}}$	${▾}_{\mathit{u}}(Д\breve{z})$$={▾}_{\mathit{u}}\left[\left[\left[\zeta \breve{\underline{z}}\right]\right]\right]$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left({▾}_{\mathit{u}}\left(\zeta \breve{\underline{z}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \breve{\underline{z}}{\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[C (m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]	${\mathcal{L}}_{{\tilde{\mathsf{\Gamma }}}^{\top }∟\mathit{u},t}^{(\mathit{u},t)}\left(\zeta \breve{\underline{z}}{vol}_{\mathit{x}}^3\right)$	$=\left({▾}_{\mathit{u}}\left(\zeta \breve{\underline{z}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \breve{\underline{z}}{\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[C (m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]
Entropy	$\zeta \breve{\underline{s}}$	${▾}_{\mathit{u}}(Д\breve{s})$$={▾}_{\mathit{u}}\left[\left[\left[\zeta \breve{\underline{s}}\right]\right]\right]$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left({▾}_{\mathit{u}}\left(\zeta \breve{\underline{s}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \breve{\underline{s}}{\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[J K$-1$ (m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]	${\mathcal{L}}_{{\tilde{\mathsf{\Gamma }}}^{\top }∟\mathit{u},t}^{(\mathit{u},t)}\left(\zeta \breve{\underline{s}}{vol}_{\mathit{x}}^3\right)$	$=({▾}_{\mathit{u}}\left(\zeta \breve{\underline{s}}\right)+{\nabla }_{\mathit{x}}·\left(\zeta \breve{\underline{s}}{\tilde{\mathsf{\Gamma }}}^{\top }\right))dV$	[J K$-1$ (m s$-1$)$-3$${\mathit{d}}_{\mathit{u},t}^{-1}$]
Probability	$p(\mathit{u},\mathit{x}|t)$$=p_{\mathit{u},\mathit{x}|t}$	$\left.\begin{array}{c}\hfill {▾}_{\mathit{u}}{p}_{\mathit{u}|t}\\ \hfill ={▾}_{\mathit{u}}\left[\left[\left[p_{\mathit{u},\mathit{x}|t}\right]\right]\right]\end{array}\right\}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left({▾}_{\mathit{u}}{p}_{\mathit{u},\mathit{x}|t}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{u},\mathit{x}|t}{\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[(m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]	$\left.\begin{array}{c}\hfill {\mathcal{L}}_{{\tilde{\mathsf{\Gamma }}}^{\top }∟\mathit{u},t}^{(\mathit{u},t)}\\ \hfill \left(p_{\mathit{u},\mathit{x}|t}{vol}_{\mathit{x}}^3\right)\end{array}\right\}$	$=\left({▾}_{\mathit{u}}{p}_{\mathit{u},\mathit{x}|t}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{u},\mathit{x}|t}{\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[(m s$-1$)$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]
	$p\left(\mathit{x}\right|\mathit{u},t)$$=p_{\mathit{x}|\mathit{u},t}$	$\mathbf{0}={▾}_{\mathit{u}}\left[\left[\left[p_{\mathit{x}|\mathit{u},t}\right]\right]\right]$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left({▾}_{\mathit{u}}{p}_{\mathit{x}|\mathit{u},t}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{x}|\mathit{u},t}{\tilde{\mathsf{\Gamma }}}^{\top }\right)\right)dV$	[${\mathit{d}}_{\mathit{u},t}^{-1}$]	$\mathbf{0}$	$={▾}_{\mathit{u}}{p}_{\mathit{x}|\mathit{u},t}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{x}|\mathit{u},t}{\tilde{\mathsf{\Gamma }}}^{\top }\right)$	[m$-3$ ${\mathit{d}}_{\mathit{u},t}^{-1}$]

. We also list two sets of composite Liouville equations, here based on $p(\mathit{u},\mathit{x}|t)$ or $p\left(\mathit{x}\right|\mathit{u},t)$.

5.10. Velocimetric-Temporal Formulation (Density $\theta$)

Now consider the alternative velocimetric-temporal formulation on path ➄ in Figure 5, based on the generalized density $\psi =\theta (\mathit{u},\mathit{x},t)$ defined from the fluid density $\eta (\mathit{u},\mathit{x},t)$. From the definition of $\theta$ in Equations (34) and (36):

$$\begin{array}{c}\hfill \begin{array}{c}\hfill d\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\theta dU=dQ\left(t\right)=\frac{dQ}{dt}dt\end{array}\end{array}$$

(65)

Note the peculiar property that integrating $\theta (\mathit{u},\mathit{x},t)$ with respect to $\mathit{u}$ also eliminates $\mathit{x}$, giving a velocimetric-temporal formulation with $\mathit{x}=\mathit{u}$, $\mathit{C}=t$ and $\mathit{V}=\dot{\mathit{u}}$. From Equation (43):

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \frac{dQ}{dt}=\frac{d}{dt}\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\theta dU& =\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[\frac{\partial \theta }{\partial t}+{\nabla }_{\mathit{u}}·\left(\theta \dot{\mathit{u}}\right)\right]dU\hfill \end{array}\end{array}$$

(66)

The left-hand terms are again equivalent to those in Table 2, Table 3 and Table 4. The extracted continuity equation and differential equation based on $\breve{\underline{\theta }}$ (Appendix C) are:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill 0& =\frac{\partial \eta }{\partial t}+{\nabla }_{\mathit{u}}·\left(\eta \dot{\mathit{u}}\right)\hfill \end{array}\end{array}$$

(67)

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \eta \frac{d\breve{\underline{\theta }}}{dt}& =\frac{\partial \left(\eta \breve{\underline{\theta }}\right)}{\partial t}+{\nabla }_{\mathit{u}}·\left(\eta \breve{\underline{\theta }}\dot{\mathit{u}}\right)\hfill \end{array}\end{array}$$

(68)

The conservation laws derived from Equations (66) and (68) are listed in

Table 11. Conservation Laws for the Velocimetric-Temporal Formulation (based on the Velocimetric Fluid Density $\eta (\mathit{u},\mathit{x},t)$).

Conserved Quantity	Density	Integral Equation			Differential Equation
	$\theta (\mathit{u},\mathit{x},t)$$=\eta \breve{\underline{\theta }}$	LHS = ${\displaystyle \frac{dQ}{dt}}$ = ${\displaystyle \frac{DQ}{Dt}}$	= RHS	SI Units	LHS	= RHS	SI Units
Fluid mass	$\eta$	$0=\dot{M}$$={\displaystyle \frac{d\langle \langle \langle \eta \rangle \rangle \rangle }{dt}}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \eta }{\partial t}}+{\nabla }_{\mathit{u}}·\left(\eta \dot{\mathit{u}}\right)\right]dU$	[kg s$-1$]	0	$={\displaystyle \frac{\partial \eta }{\partial t}}+{\nabla }_{\mathit{u}}·\left(\eta \dot{\mathit{u}}\right)$	[kg s$-1$ (m s$-1$)$-3$]
Species mass	$\eta {\breve{\underline{\chi }}}_c$ = ${\eta }_c$	$\left({\dot{M}}_c\right)={\displaystyle \frac{d\langle \langle \langle \eta {\breve{\underline{\chi }}}_c\rangle \rangle \rangle }{dt}}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \eta {\breve{\underline{\chi }}}_c}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\eta {\breve{\underline{\chi }}}_c\dot{\mathit{u}}\right)\right]dU$	[kgc s$-1$]	$\eta {\displaystyle \frac{d{\breve{\underline{\chi }}}_c}{dt}}$	$={\displaystyle \frac{\partial \eta {\breve{\underline{\chi }}}_c}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\eta {\breve{\underline{\chi }}}_c\dot{\mathit{u}}\right)$	[kgc s$-1$ (m s$-1$)$-3$]
Linear momentum	$\eta \mathit{u}$	$\sum \mathit{F}={\displaystyle \frac{d\langle \langle \langle \eta \mathit{u}\rangle \rangle \rangle }{dt}}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \eta \mathit{u}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\eta \mathit{u}\dot{\mathit{u}}\right)\right]dU$	[(kg m s$-1$) s$-1$ = N]	$\eta {\displaystyle \frac{d\mathit{u}}{dt}}$	$={\displaystyle \frac{\partial \eta \mathit{u}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\eta \mathit{u}\dot{\mathit{u}}\right)$	[N (m s$-1$)$-3$]
Angular momentum	$\eta (\breve{\underline{\mathfrak{r}}}\times \mathit{u})$	$\sum \mathit{T}={\displaystyle \frac{d\langle \langle \langle \eta (\breve{\underline{\mathfrak{r}}}\times \mathit{u})\rangle \rangle \rangle }{dt}}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }[{\displaystyle \frac{\partial \eta (\breve{\underline{\mathfrak{r}}}\times \mathit{u})}{\partial t}}$$+{\nabla }_{\mathit{u}}·\left(\eta (\breve{\underline{\mathfrak{r}}}\times \mathit{u})\dot{\mathit{u}}\right)]dU$	[(kg m${}^2$ s$-1$) s$-1$ = N m]	$\eta {\displaystyle \frac{d(\breve{\underline{\mathfrak{r}}}\times \mathit{u})}{dt}}$	$={\displaystyle \frac{\partial \eta (\breve{\underline{\mathfrak{r}}}\times \mathit{u})}{\partial t}}$$+{\nabla }_{\mathit{u}}·\left(\eta (\breve{\underline{\mathfrak{r}}}\times \mathit{u})\dot{\mathit{u}}\right)$	[(N m) (m s$-1$)$-3$]
Energy	$\eta \breve{\underline{e}}$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{DE}{Dt}}={\displaystyle \frac{d\langle \langle \langle \eta \breve{\underline{e}}\rangle \rangle \rangle }{dt}}\\ \hfill =({\dot{Q}}_{in}+{\dot{W}}_{in})\end{array}\right\}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \eta \breve{\underline{e}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\eta \breve{\underline{e}}\dot{\mathit{u}}\right)\right]dU$	[J s$-1$ = W]	$\eta {\displaystyle \frac{d\breve{\underline{e}}}{dt}}$	$={\displaystyle \frac{\partial \eta \breve{\underline{e}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\eta \breve{e}\dot{\mathit{u}}\right)$	[J s$-1$ (m s$-1$)$-3$ = W (m s$-1$)$-3$]
Charge (in solution)	$\eta \breve{\underline{z}}$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{DZ}{Dt}}={\displaystyle \frac{d\langle \langle \langle \eta \breve{\underline{z}}\rangle \rangle \rangle }{dt}}\\ \hfill =I+\left({\sum }_c{z}_c{\dot{M}}_c\right)\end{array}\right\}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \eta \breve{\underline{z}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\eta \breve{\underline{z}}\dot{\mathit{u}}\right)\right]dU$	[C s$-1$ = A]	$\eta {\displaystyle \frac{d\breve{\underline{z}}}{dt}}$	$={\displaystyle \frac{\partial \eta \breve{\underline{z}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\eta \breve{\underline{z}}\dot{\mathit{u}}\right)$	[C s$-1$ (m s$-1$)$-3$ = A (m s$-1$)$-3$]
Entropy	$\eta \breve{\underline{s}}$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{DS}{Dt}}={\displaystyle \frac{d\langle \langle \langle \eta \breve{\underline{s}}\rangle \rangle \rangle }{dt}}\\ \hfill =\dot{\sigma }+{\dot{S}}^{nf}\end{array}\right\}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \eta \breve{\underline{s}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\eta \breve{\underline{s}}\dot{\mathit{u}}\right)\right]dU$	[J K$-1$ s$-1$]	$\eta {\displaystyle \frac{d\breve{\underline{s}}}{dt}}$	$={\displaystyle \frac{\partial \eta \breve{\underline{s}}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(\eta \breve{\underline{s}}\dot{\mathit{u}}\right)$	[J K$-1$ s$-1$ (m s$-1$)$-3$]
Probability	$p\left(\mathit{u}\right|\mathit{x},t)$$=p_{\mathit{u}|\mathit{x},t}$	$0={\displaystyle \frac{d\langle \langle \langle p_{\mathit{u}|\mathit{x},t}\rangle \rangle \rangle }{dt}}$	$=\underset{\mathcal{D}(\mathit{x},t)}{\int \int \int }\left[{\displaystyle \frac{\partial p_{\mathit{u}|\mathit{x},t}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u}|\mathit{x},t}\dot{\mathit{u}}\right)\right]dU$	[s$-1$]	0	$={\displaystyle \frac{\partial p_{\mathit{u}|\mathit{x},t}}{\partial t}}+{\nabla }_{\mathit{u}}·\left(p_{\mathit{u}|\mathit{x},t}\dot{\mathit{u}}\right)$	[(m s$-1$)$-3$ s$-1$]

. Apart from the different fluid density $\eta$ and specific density labels, these are mathematically identical to the conservation laws based on the velocimetric density Д in Table 3. The Liouville equations based on equating $\theta$ to its underlying pdf $p\left(\mathit{u}\right|\mathit{x},t)$ are also listed, giving an alternative route to these equations to that given in Table 5.

5.11. Volumetric-Temporal Formulation (Density $ϵ$)

Now consider the alternative volumetric-temporal formulation on path ⑥ in Figure 5, based on the generalized density $\psi =ϵ(\mathit{u},\mathit{x},t)$ defined from the fluid density $\xi (\mathit{u},\mathit{x},t)$. From the definition of $ϵ$ in Equations (35) and (37):

$$\begin{array}{c}\hfill \begin{array}{c}\hfill d\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }ϵdV=dQ\left(t\right)=\frac{dQ}{dt}dt\end{array}\end{array}$$

(69)

This also has the peculiar property that integration with respect to $\mathit{x}$ eliminates $\mathit{u}$, reducing this to a volumetric-temporal formulation with $\mathit{X}=\mathit{x}$, $\mathit{C}=t$ and $\mathit{V}=\mathit{u}$. From Equation (43):

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \frac{dQ}{dt}=\frac{d}{dt}\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }ϵdV& =\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[\frac{\partial ϵ}{\partial t}+{\nabla }_{\mathit{x}}·\left(ϵ\mathit{u}\right)\right]dV\hfill \end{array}\end{array}$$

(70)

The left-hand terms are again equivalent to those in Table 2, Table 3 and Table 4. The extracted continuity equation and differential equation based on $\breve{\underline{ϵ}}$ (Appendix C) are:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill 0& =\frac{\partial \xi }{\partial t}+{\nabla }_{\mathit{x}}·\left(\xi \mathit{u}\right)\hfill \end{array}\end{array}$$

(71)

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \xi \frac{d\breve{\underline{ϵ}}}{dt}& =\frac{\partial \left(\xi \breve{\underline{ϵ}}\right)}{\partial t}+{\nabla }_{\mathit{x}}·\left(\xi \breve{\underline{ϵ}}\mathit{u}\right)\hfill \end{array}\end{array}$$

(72)

The conservation laws derived from Equations (70) and (72) are listed in

Table 12. Conservation Laws for the Volumetric-Temporal Formulation (based on the Volumetric Fluid Density $\xi (\mathit{u},\mathit{x},t)$).

Conserved quantity	Density	Integral Equation			Differential Equation
	$ϵ(\mathit{u},\mathit{x},t)$$=\xi \breve{\underline{ϵ}}$	LHS = ${\displaystyle \frac{dQ}{dt}}$ = ${\displaystyle \frac{DQ}{Dt}}$	= RHS	SI Units	LHS	= RHS	SI Units
Fluid mass	$\xi$	$0=\dot{M}$$={\displaystyle \frac{d\left[\right[\left[\xi \right]\left]\right]}{dt}}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \xi }{\partial t}}+{\nabla }_{\mathit{x}}·\left(\xi \mathit{u}\right)\right]dV$	[kg s$-1$]	0	$={\displaystyle \frac{\partial \xi }{\partial t}}+{\nabla }_{\mathit{x}}·\left(\xi \mathit{u}\right)$	[kg s$-1$ m$-3$]
Species mass	$\xi {\breve{\underline{\chi }}}_c$$={\xi }_c$	$\left({\dot{M}}_c\right)={\displaystyle \frac{d\left[\left[\left[\xi {\breve{\underline{\chi }}}_c\right]\right]\right]}{dt}}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \xi {\breve{\underline{\chi }}}_c}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\xi {\breve{\underline{\chi }}}_c\mathit{u}\right)\right]dV$	[kgc s$-1$]	$\xi {\displaystyle \frac{D{\breve{\underline{\chi }}}_c}{Dt}}$	$={\displaystyle \frac{\partial \xi {\breve{\underline{\chi }}}_c}{\partial t}}$$+{\nabla }_{\mathit{x}}·\left(\xi {\breve{\underline{\chi }}}_c\mathit{u}\right)$	[kgc s$-1$ m$-3$]
Linear momentum	$\xi \mathit{u}$	$\sum \mathit{F}={\displaystyle \frac{d\left[\right[\left[\xi \mathit{u}\right]\left]\right]}{dt}}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \xi \mathit{u}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\xi \mathit{u}\mathit{u}\right)\right]dV$	[(kg m s$-1$) s$-1$ = N]	$\xi {\displaystyle \frac{D\mathit{u}}{Dt}}$	$={\displaystyle \frac{\partial \xi \mathit{u}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\xi \mathit{u}\mathit{u}\right)$	[N m$-3$]
Angular momentum	$\xi (\breve{\underline{\mathfrak{r}}}\times \mathit{u})$	$\sum \mathit{T}={\displaystyle \frac{d\left[\left[\left[\xi (\breve{\underline{\mathfrak{r}}}\times \mathit{u})\right]\right]\right]}{dt}}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }[{\displaystyle \frac{\partial \xi (\breve{\underline{\mathfrak{r}}}\times \mathit{u})}{\partial t}}$$+{\nabla }_{\mathit{x}}·\left(\xi (\breve{\underline{\mathfrak{r}}}\times \mathit{u})\mathit{u}\right)]dV$	[(kg m${}^2$ s$-1$) s$-1$ = N m]	$\xi {\displaystyle \frac{D(\breve{\underline{\mathfrak{r}}}\times \mathit{u})}{Dt}}$	$={\displaystyle \frac{\partial \xi (\breve{\underline{\mathfrak{r}}}\times \mathit{u})}{\partial t}}$$+{\nabla }_{\mathit{x}}·\left(\xi (\breve{\underline{\mathfrak{r}}}\times \mathit{u})\mathit{u}\right)$	[N m$-2$ = (N m) m$-3$]
Energy	$\xi \breve{\underline{e}}$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{DE}{Dt}}={\displaystyle \frac{d\left[\left[\left[\xi \breve{\underline{e}}\right]\right]\right]}{dt}}\\ \hfill =({\dot{Q}}_{in}+{\dot{W}}_{in})\end{array}\right\}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \xi \breve{\underline{e}}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\xi \breve{\underline{e}}\mathit{u}\right)\right]dV$	[J s$-1$ = W]	$\xi {\displaystyle \frac{D\breve{\underline{e}}}{Dt}}$	$={\displaystyle \frac{\partial \xi \breve{\underline{e}}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\xi \breve{\underline{e}}\mathit{u}\right)$	[J s$-1$ m$-3$ = W m$-3$]
Charge (in solution)	$\xi \breve{\underline{z}}$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{DZ}{Dt}}={\displaystyle \frac{d\left[\left[\left[\xi \breve{\underline{z}}\right]\right]\right]}{dt}}\\ \hfill =I+\left({\sum }_c{z}_c{\dot{M}}_c\right)\end{array}\right\}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \xi \breve{\underline{z}}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\xi \breve{\underline{z}}\mathit{u}\right)\right]dV$	[C s$-1$ = A]	$\xi {\displaystyle \frac{D\breve{\underline{z}}}{Dt}}$	$={\displaystyle \frac{\partial \xi \breve{\underline{z}}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\xi \breve{\underline{z}}\mathit{u}\right)$	[C s$-1$ m$-3$ = A m$-3$]
Entropy	$\xi \breve{\underline{s}}$	$\left.\begin{array}{c}\hfill {\displaystyle \frac{DS}{Dt}}={\displaystyle \frac{d\left[\left[\left[\xi \breve{\underline{s}}\right]\right]\right]}{dt}}\\ \hfill =\dot{\sigma }+{\dot{S}}^{nf}\end{array}\right\}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial \xi \breve{\underline{s}}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\xi \breve{\underline{s}}\mathit{u}\right)\right]dV$	[J K$-1$ s$-1$]	$\xi {\displaystyle \frac{D\breve{\underline{s}}}{Dt}}$	$={\displaystyle \frac{\partial \xi \breve{\underline{s}}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(\xi \breve{\underline{s}}\mathit{u}\right)$	[J K$-1$ s$-1$ m$-3$]
Probability	$p\left(\mathit{x}\right|\mathit{u},t)$$=p_{\mathit{x}|\mathit{u},t}$	$0={\displaystyle \frac{d\left[\left[\left[p_{\mathit{x}|\mathit{u},t}\right]\right]\right]}{dt}}$	$=\underset{\mathsf{\Omega }(\mathit{u},t)}{\int \int \int }\left[{\displaystyle \frac{\partial p_{\mathit{x}|\mathit{u},t}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{x}|\mathit{u},t}\mathit{u}\right)\right]dV$	[s$-1$]	0	$={\displaystyle \frac{\partial p_{\mathit{x}|\mathit{u},t}}{\partial t}}+{\nabla }_{\mathit{x}}·\left(p_{\mathit{x}|\mathit{u},t}\mathit{u}\right)$	[s$-1$ m$-3$]

. Apart from the fluid density $\xi$ and specific density labels, these are mathematically identical to the volumetric conservation laws based on $\rho$ in Table 2. The Liouville equations based on the underlying pdf $p\left(\mathit{x}\right|\mathit{u},t)$ are also listed, providing an alternative route to that in Table 6.

5.12. Discussion

Comparing the above analyses, we see that the first three and last two formulations in Section 5.1, Section 5.2 and Section 5.3 and Section 5.10 and Section 5.11 are connected by their equivalence to the rate of change of the conserved quantity $dQ/dt$, and thence to the source-sink terms of the standard (volumetric) integral equations of fluid mechanics. These are based on complete integrations, respectively, of the densities $\alpha$, $\beta$, $\phi$, $\theta$ and $ϵ$.

The remaining examples involve partial integrations of the velocivolumetric density $\phi$. Of these, the velocimetric-spatiotemporal formulation in Section 5.8 combines the velocimetric-temporal and velocimetric-spatial formulations in Section 5.4 and Section 5.6, connected by Equation (57), based on $\alpha$ and $\phi$. Similarly, the volumetric-velocitemporal formulation in Section 5.9 combines the volumetric-temporal and volumetric-velocital formulations in Section 5.5 and Section 5.7, connected by Equation (59), based on $\beta$ and $\phi$. For these formulations, it is possible to extract a Lie differential equation based on a Lie operator specific to that formulation.

All formulations can be used to derive a number of Liouville equations based on the underlying pdfs defined in Section 3.1 and Figure 2. In general, these take the form of Lie differential equations, but for cases in which the pdf is normalized by the integration, the analysis yields a homogeneous integral equation and associated partial differential equation.

6. Conclusions

This study examines a generalized framework for the Reynolds transport theorem [23,24], which enables continuous multiparametric mappings of a differential form on a manifold—or of a density within a generalized coordinate space—connected by the maximal integral curves or surfaces described by a vector or tensor field. These extend the formulation of [1] to encompass new transformation theorems, which exploit previously unreported multiparametric continuous (Lie) symmetries of a vector or tensor field associated with a conserved quantity. In this study, we explore the implications of the generalized framework for fluid flow systems, using an extended Eulerian velocivolumetric description of fluid flow in place of the standard Eulerian description.

The analysis commences in Section 2 with a detailed discussion of the extended velocity-position description, and of the form and connection between the geometric and velocity domains for different classes of fluid flow systems. In Section 3.1, we then define a hierarchy of five pdfs $\left\{p\right(\mathit{x}\left|t\right),p\left(\mathit{u}\right|t),p(\mathit{u},\mathit{x}\left|t\right),p\left(\mathit{u}\right|\mathit{x},t),p(\mathit{x}|\mathit{u},t)\}$ within this description, the properties of which are then explored. In Section 3.2 and Section 3.3, these are used by convolution to define analogous hierarchies of fluid densities $\{\rho ,Д,\zeta ,\eta ,\xi \}$ and generalized densities $\{\alpha ,\beta ,\phi ,\theta ,ϵ\}$, of which only the first ($\rho$ and $\alpha$) are commonly known. The generalized framework for the Reynolds transport theorem is presented in Section 4, in both exterior calculus and vector calculus formulations. Its connection to underlying partial differential equations and Lie differential equations—the latter containing the Lie derivative of a volume form in the domain—is also examined. In Section 5, the densities and theorems are used to obtain 11 formulations of the Reynolds transport theorem arising from the veloci-volumetric description, for different choices of the coordinate space, parameter space and density. These are reported in the form of 11 tables of integral and differential conservation laws applicable to these systems, for the eight common conserved quantities of interest in fluid mechanics (fluid mass, species mass, linear momentum, angular momentum, energy, charge, entropy and probability). The equations for conservation of probability can be interpreted as analogs of the Liouville equation, applicable to different spaces. The analyses provide a considerable assortment of new conservation laws for the analysis of fluid flow systems.

While every effort has been made to provide a comprehensive treatment, this study exclusively considers compact domains with Cartesian position and velocity coordinates, and no attempt is made to examine orthogonal or non-orthogonal curvilinear coordinate systems or parameter spaces (however, care is taken to distinguish the vector of partial derivatives ${\mathbf{\partial }}_{\mathit{x}}$ from the gradient ${\nabla }_{\mathit{x}}$). Further work is required to extend these analyses to more general domains and coordinate systems. Further detailed study is also required of the partial and Lie differential equations reported in Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9, Table 10, Table 11 and Table 12, to identify the appropriate source-sink terms for each expression. The analyses are also restricted to fixed frames of reference (control volumes or domains), but can readily be extended to moving and deforming frames of reference using relative vector or tensor fields ${\mathit{V}}_{rel}$, e.g., [3,6,23]. The analyses could also be extended to consider domains with jump discontinuities [7,8], irregular and fragmenting domains [9,10], or special or general relativity [81]. The generalized Reynolds theorem framework can also be used to generate conservation laws for other dynamical systems containing conserved quantities [23].