Reynolds transport theorem

Reynolds' transport theorem (also known as the Leibniz-Reynolds' transport theorem), or in short Reynolds theorem, is a three-dimensional generalization of the Leibniz integral rule which is also known as differentiation under the integral sign. The theorem is named after Osborne Reynolds (1842–1912). It is used to recast derivatives of integrated quantities and is useful in formulating the basic equations of continuum mechanics.

Consider integrating \(\mathbf{f} = \mathbf{f}(\mathbf{x},t)\) over the time-dependent region \(\Omega(t)\) that has boundary \(\partial \Omega(t)\), then taking the derivative with respect to time: \[ \cfrac{\mathrm{d}}{\mathrm{d}t}\int_{\Omega(t)} \mathbf{f}~\text{dV} ~. \] If we wish to move the derivative under the integral sign there are two issues: the time dependence of \(\mathbf{f}\), and the introduction of and removal of space from \(\Omega\) due to its dynamic boundary. Reynolds' transport theorem provides the necessary framework.

General form

Reynolds' transport theorem, derived in , is: \[ \cfrac{\mathrm{d}}{\mathrm{d}t}\int_{\Omega(t)} \mathbf{f}~\text{dV} = \int_{\Omega(t)} \frac{\partial \mathbf{f}}{\partial t}~\text{dV} + \int_{\partial \Omega(t)} (\mathbf{v}^{b}\cdot\mathbf{n})\mathbf{f}~\text{dA} ~ \] in which \(\mathbf{n}(\mathbf{x},t)\) is the outward-pointing unit-normal, \(\mathbf{x}\) is a point in the region and is the variable of integration, \(\text{dV}\) and \(\text{dA}\) are volume and surface elements at \(\mathbf{x}\), and \(\mathbf{v}^{b}(\mathbf{x},t)\) is the velocity of the area element – so not necessarily the flow velocity. The function \(\mathbf{f}\) may be vector or scalar valued. Note that the integral on the left hand side is a function solely of time, and so the total derivative has been used.

Form for a material element

In continuum mechanics this theorem is often used for material elements, which are parcels of fluids or solids which no material enters or leaves. If \(\Omega(t)\) is a material element then there is a velocity function \(\mathbf{v}=\mathbf{v}(\mathbf{x},t)\) and the boundary elements obey \[\mathbf{v}^{b}\cdot\mathbf{n}=\mathbf{v}\cdot\mathbf{n}.\] This condition may be substituted to obtain \[ \cfrac{\mathrm{d}}{\mathrm{d}t}\left(\int_{\Omega(t)} \mathbf{f}~\text{dV}\right) = \int_{\Omega(t)} \frac{\partial \mathbf{f}}{\partial t}~\text{dV} + \int_{\partial \Omega(t)} (\mathbf{v}\cdot\mathbf{n})\mathbf{f}~\text{dA} ~. \]

Proof for a material element

Let \(\Omega_0\) be reference configuration of the region \(\Omega(t)\). Let the motion and the deformation gradient be given by \[ \mathbf{x} = \boldsymbol{\varphi}(\mathbf{X}, t)~; \qquad\implies\qquad \boldsymbol{F}(\mathbf{X},t) = \boldsymbol{\nabla}_{\circ} \boldsymbol{\varphi} ~. \] Let \(J(\mathbf{X},t) = \det[\boldsymbol{F}(\mathbf{X},t)]\). Then, integrals in the current and the reference configurations are related by \[ \int_{\Omega(t)} \mathbf{f}(\mathbf{x},t)~\text{dV} = \int_{\Omega_0} \mathbf{f}[\boldsymbol{\varphi}(\mathbf{X},t),t]~J(\mathbf{X},t)~\text{dV}_0 = \int_{\Omega_0} \hat{\mathbf{f}}(\mathbf{X},t)~J(\mathbf{X},t)~\text{dV}_0 ~. \] That this derivation is for a material element is implicit in the time constancy of the reference configuration: it is constant in material coordinates. The time derivative of an integral over a volume is defined as \[ \cfrac{d}{dt}\left( \int_{\Omega(t)} \mathbf{f}(\mathbf{x},t)~\text{dV}\right) = \lim_{\Delta t \rightarrow 0} \cfrac{1}{\Delta t} \left(\int_{\Omega(t + \Delta t)} \mathbf{f}(\mathbf{x},t+\Delta t)~\text{dV} - \int_{\Omega(t)} \mathbf{f}(\mathbf{x},t)~\text{dV}\right) ~. \] Converting into integrals over the reference configuration, we get \[ \cfrac{d}{dt}\left( \int_{\Omega(t)} \mathbf{f}(\mathbf{x},t)~\text{dV}\right) = \lim_{\Delta t \rightarrow 0} \cfrac{1}{\Delta t} \left(\int_{\Omega_0} \hat{\mathbf{f}}(\mathbf{X},t+\Delta t)~J(\mathbf{X},t+\Delta t)~\text{dV}_0 - \int_{\Omega_0} \hat{\mathbf{f}}(\mathbf{X},t)~J(\mathbf{X},t)~\text{dV}_0\right) ~. \] Since \(\Omega_0\) is independent of time, we have \[ \begin{align} \cfrac{d}{dt}\left( \int_{\Omega(t)} \mathbf{f}(\mathbf{x},t)~\text{dV}\right) & = \int_{\Omega_0} \left[\lim_{\Delta t \rightarrow 0} \cfrac{ \hat{\mathbf{f}}(\mathbf{X},t+\Delta t)~J(\mathbf{X},t+\Delta t) - \hat{\mathbf{f}}(\mathbf{X},t)~J(\mathbf{X},t)}{\Delta t} \right]~\text{dV}_0 \\ & = \int_{\Omega_0} \frac{\partial }{\partial t}[\hat{\mathbf{f}}(\mathbf{X},t)~J(\mathbf{X},t)]~\text{dV}_0 \\ & = \int_{\Omega_0} \left( \frac{\partial }{\partial t}[\hat{\mathbf{f}}(\mathbf{X},t)]~J(\mathbf{X},t)+ \hat{\mathbf{f}}(\mathbf{X},t)~\frac{\partial }{\partial t}[J(\mathbf{X},t)]\right) ~\text{dV}_0 \end{align} \] Now, the time derivative of \(\det\boldsymbol{F}\) is given by \[ \frac{\partial J(\mathbf{X},t)}{\partial t} = \frac{\partial }{\partial t}(\det\boldsymbol{F}) = (\det\boldsymbol{F})(\boldsymbol{\nabla} \cdot \mathbf{v}) = J(\mathbf{X},t)~\boldsymbol{\nabla} \cdot \mathbf{v}(\boldsymbol{\varphi}(\mathbf{X},t),t) = J(\mathbf{X},t)~\boldsymbol{\nabla} \cdot \mathbf{v}(\mathbf{x},t) ~. \] Therefore, \[ \begin{align} \cfrac{d}{dt}\left( \int_{\Omega(t)} \mathbf{f}(\mathbf{x},t)~\text{dV}\right) & = \int_{\Omega_0} \left( \frac{\partial }{\partial t}[\hat{\mathbf{f}}(\mathbf{X},t)]~J(\mathbf{X},t)+ \hat{\mathbf{f}}(\mathbf{X},t)~J(\mathbf{X},t)~\boldsymbol{\nabla} \cdot \mathbf{v}(\mathbf{x},t)\right) ~\text{dV}_0 \\ & = \int_{\Omega_0} \left(\frac{\partial }{\partial t}[\hat{\mathbf{f}}(\mathbf{X},t)]+ \hat{\mathbf{f}}(\mathbf{X},t)~\boldsymbol{\nabla} \cdot \mathbf{v}(\mathbf{x},t)\right)~J(\mathbf{X},t) ~\text{dV}_0 \\ & = \int_{\Omega(t)} \left(\dot{\mathbf{f}}(\mathbf{x},t)+ \mathbf{f}(\mathbf{x},t)~\boldsymbol{\nabla} \cdot \mathbf{v}(\mathbf{x},t)\right)~\text{dV} \end{align} \] where \(\dot{\mathbf{f}}\) is the material time derivative of \(\mathbf{f}\). Now, the material derivative is given by \[ \dot{\mathbf{f}}(\mathbf{x},t) = \frac{\partial \mathbf{f}(\mathbf{x},t)}{\partial t} + [\boldsymbol{\nabla} \mathbf{f}(\mathbf{x},t)]\cdot\mathbf{v}(\mathbf{x},t) ~. \] Therefore, \[ \cfrac{d}{dt}\left( \int_{\Omega(t)} \mathbf{f}(\mathbf{x},t)~\text{dV}\right) = \int_{\Omega(t)} \left( \frac{\partial \mathbf{f}(\mathbf{x},t)}{\partial t} + [\boldsymbol{\nabla} \mathbf{f}(\mathbf{x},t)]\cdot\mathbf{v}(\mathbf{x},t) + \mathbf{f}(\mathbf{x},t)~\boldsymbol{\nabla} \cdot \mathbf{v}(\mathbf{x},t)\right)~\text{dV} \] or, \[ \cfrac{d}{dt}\left( \int_{\Omega(t)} \mathbf{f}~\text{dV}\right) = \int_{\Omega(t)} \left( \frac{\partial \mathbf{f}}{\partial t} + \boldsymbol{\nabla} \mathbf{f}\cdot\mathbf{v} + \mathbf{f}~\boldsymbol{\nabla} \cdot \mathbf{v}\right)~\text{dV} ~. \] Using the identity \[ \boldsymbol{\nabla} \cdot (\mathbf{v}\otimes\mathbf{w}) = \mathbf{v}(\boldsymbol{\nabla} \cdot \mathbf{w}) + \boldsymbol{\nabla}\mathbf{v}\cdot\mathbf{w} \] we then have \[ \cfrac{d}{dt}\left( \int_{\Omega(t)} \mathbf{f}~\text{dV}\right) = \int_{\Omega(t)} \left(\frac{\partial \mathbf{f}}{\partial t} + \boldsymbol{\nabla} \cdot (\mathbf{f}\otimes\mathbf{v})\right)~\text{dV} ~. \] Using the divergence theorem and the identity \((\mathbf{a}\otimes\mathbf{b})\cdot\mathbf{n} = (\mathbf{b}\cdot\mathbf{n})\mathbf{a}\) we have \[ { \cfrac{d}{dt}\left( \int_{\Omega(t)} \mathbf{f}~\text{dV}\right) = \int_{\Omega(t)}\frac{\partial \mathbf{f}}{\partial t}~\text{dV} + \int_{\partial \Omega(t)}(\mathbf{f}\otimes\mathbf{v})\cdot\mathbf{n}~\text{dA} = \int_{\Omega(t)}\frac{\partial \mathbf{f}}{\partial t}~\text{dV} + \int_{\partial \Omega(t)}(\mathbf{v}\cdot\mathbf{n})\mathbf{f}~\text{dA} \qquad \square } \]

Erroneous sources

This theorem is widely quoted, incorrectly, as being the form specific to material volumes. See the planetmath external link below for an example. Clearly, if the material volume form is applied to regions other than material volumes, errors will ensue.

A special case

If we take \(\Omega\) to be constant with respect to time, then \(\mathbf{v}_b=0\) and the identity reduces to \[ \cfrac{\mathrm{d}}{\mathrm{d}t}\int_{\Omega} f~\text{dV} = \int_{\Omega} \frac{\partial f}{\partial t}~\text{dV} ~, \] as expected. This simplification is not possible if an incorrect form of the Reynolds transport theorem is used.

Interpretation and reduction to one dimension

The theorem is the higher dimensional extension of Differentiation under the integral sign and should reduce to that expression in some cases. Suppose \(f\) is independent of \(y\) & \(z\), and that \(\Omega(t)\) is a unit square in the \(y-z\) plane and has \(x\) limits \(a(t)\) and \(b(t)\). Then Reynolds transport theorem reduces to \[ \cfrac{\mathrm{d}}{\mathrm{d}t}\int_{a(t)}^{b(t)} f~\text{dx} = \int_{a(t)}^{b(t)} \frac{\partial f}{\partial t}~\text{dx} + \frac{\partial b(t)}{\partial t} f(b(t),t) -\frac{\partial a(t)}{\partial t} f(a(t),t) ~, \] which is the expression given on Differentiation under the integral sign, except that there the variables x and t have been swapped.

See also

• Differentiation under the integral sign
• Leibniz integral rule

References

L. G. Leal, 2007, Advanced transport phenomena: fluid mechanics and convective transport processes, Cambridge University Press, p. 912.

O. Reynolds, 1903, Papers on Mechanical and Physical Subjects, Vol. 3, The Sub-Mechanics of the Universe, Cambridge University Press, Cambridge.

External links

• Osborne Reynolds, Collected Papers on Mechanical and Physical Subjects, in three volumes, published circa 1903, now fully and freely available in digital format:

Volume 1

Volume 2

Volume 3

• http://www.catea.org/grade/mecheng/mod6/mod6.html#slide1
• http://planetmath.org/encyclopedia/ReynoldsTransportTheorem.htmlde:Reynolds’scher Transportsatz