Material derivative

The material derivative^[1][2] is a derivative taken along a path moving with velocity v, and is often used in fluid mechanics and classical mechanics. It describes the time rate of change of some quantity (such as heat or momentum) by following it, while moving with a space- and time-dependent velocity field.

For example, in fluid dynamics, take the case that the velocity field under consideration is the flow velocity itself, and the quantity of interest is the temperature of the fluid. Then the material derivative describes the temperature evolution of a certain fluid parcel in time, as it is being moved along its pathline (trajectory) while following the fluid flow.

The material derivative can serve as a link between Eulerian and Lagrangian descriptions of fluid motion.^[3]

Names

There are many other names for this operator, including:

• convective derivative^[4]
• advective derivative
• substantive derivative^[5]
• substantial derivative^[1]
• Lagrangian derivative^[6]
• Stokes derivative^[5]
• particle derivative
• hydrodynamic derivative^[1]
• derivative following the motion^[1]
• total derivative^[1]

Definition

The material derivatives of a scalar field φ( x, t ) and a vector field u( x, t ) are defined respectively as:

$\frac{D\varphi}{Dt} = \frac{\partial \varphi}{\partial t} + \mathbf{v}\cdot\nabla \varphi,$

$\frac{D\mathbf{u}}{Dt} = \frac{\partial \mathbf{u}}{\partial t} + \mathbf{v}\cdot\nabla \mathbf{u},$

where the distinction is that $\nabla \varphi$ is the gradient of a scalar, while $\nabla \mathbf{u}$ is the covariant derivative of a vector. In case of the material derivative of a vector field, the term v•∇u can both be interpreted as v•(∇u) involving the tensor derivative of u, or as (v•∇)u, leading to the same result.^[7]

Confusingly, the term convective derivative is both used for the whole material derivative Dφ/Dt or Du/Dt, and for only the spatial rate of change part, v•∇φ or v•∇u respectively.^[2] For that case, the convective derivative only equals D/Dt for time independent flows.

These derivatives are physical in nature and describe the transport of a scalar or vector quantity in a velocity field v( x, t ). The effect of the time independent terms in the definitions are for the scalar and vector case respectively known as advection and convection.

Development

Consider a scalar quantity φ = φ( x, t ), where t is understood as time and x as position. This may be some physical variable such as temperature or chemical concentration. The physical quantity exists in a fluid, whose velocity is represented by the vector field v( x, t ).

The (total) derivative with respect to time of φ is expanded through the multivariate chain rule:

$\frac{d}{d t}(\varphi(\mathbf x, t)) = \frac{\partial \varphi}{\partial t} + \nabla \varphi \cdot \frac{d \mathbf x}{d t}$

It is apparent that this derivative is dependent on the vector

$\frac{d \mathbf x}{d t} = \left(\frac{d x}{d t}, \frac{d y}{d t}, \frac{d z}{d t}\right)$

which describes a chosen path x(t) in space. For example, if $d \mathbf x/d t = 0$ is chosen, the time derivative becomes equal to the partial derivative, which agrees with the definition of a partial derivative: a derivative taken with respect to some variable (time in this case) holding other variables constant (space in this case). This makes sense because if $d \mathbf x/d t = 0$, then the derivative is taken at some constant position. This static position derivative is called the Eulerian derivative.

An example of this case is a swimmer standing still and sensing temperature change in a lake early in the morning: the water gradually becomes warmer due to heating from the sun.

If, instead, the path x(t) is not a standstill, the (total) time derivative of φ may change due to the path. For example, imagine the swimmer is in a motionless pool of water, indoors and unaffected by the sun. One end happens to be a constant hot temperature and the other end a constant cold temperature. By swimming from one end to the other the swimmer senses a change of temperature with respect to time, even though the temperature at any given (static) point is a constant. This is because the derivative is taken at the swimmer's changing location. A temperature sensor attached to the swimmer would show temperature varying in time, even though the pool is held at a steady temperature distribution.

The material derivative finally is obtained when the path x(t) is chosen to have a velocity equal to the fluid velocity:

$\frac{d \mathbf x}{d t} = \mathbf v$

That is, the path follows the fluid current described by the fluid's velocity field v. So, the material derivative of the scalar φ is:

$\frac{D \varphi}{D t} = \frac{\partial \varphi}{\partial t} + \nabla \varphi \cdot \mathbf v$

An example of this case is a lightweight, neutrally buoyant particle swept around in a flowing river undergoing temperature changes, maybe due to one portion of the river being sunny and the other in a shadow. The water as a whole may be heating as the day progresses. The changes due to the particle's motion (itself caused by fluid motion) is called advection (or convection if a vector is being transported).

The definition above relied on the physical nature of fluid current; however no laws of physics were invoked (for example, it hasn't been shown that a lightweight particle in a river will follow the velocity of the water). It turns out, however, that many physical concepts can be written concisely with the material derivative. The general case of advection, however, relies on conservation of mass in the fluid stream; the situation becomes slightly different if advection happens in a non-conservative medium.

Only a path was considered for the scalar above. For a vector, the gradient becomes a tensor derivative; for tensor fields we may want to take into account not only translation of the coordinate system due to the fluid movement but also its rotation and stretching. This is achieved by the upper convected time derivative.

Orthogonal coordinates

It may be shown that, in orthogonal coordinates, the j^th component of convection is given by^[8]:

$[\mathbf{v}\cdot\nabla \mathbf{u}]_j = \sum_i \frac{v_i}{h_i} \frac{\partial u_j}{\partial q^i} + \frac{u_i}{h_i h_j}\left(v_j \frac{\partial h_j}{\partial q^i} - v_i \frac{\partial h_i}{\partial q^j}\right),$

where the h_i's are related to the metric tensors by

$h_i=\sqrt{g_{ii}}.$

See also

• Navier–Stokes equations
• Euler equations
• Derivative (generalizations)
• Lie derivative

References

1. ^ ^{a b c d e} Bird, R.B., Stewart, W.E. and Lightfoot, E.N. (2007). Transport Phenomena (Revised Second Edition ed.). John Wiley & Sons. ISBN 978-0-470-11539-8.  p. 83.
2. ^ ^{a b} Batchelor, G.K. (1967). An Introduction to Fluid Dynamics. Cambridge University Press. ISBN 0521663962.  p. 72–73.
3. ^ Trenberth, K. E. (1993). Climate System Modeling. Cambridge University Press. ISBN 0521432316.  p. 99
4. ^ Ockendon, H.; Ockendon, J.R. (2004). Waves and Compressible Flow. Springer. ISBN 038740399X.  p. 6.
5. ^ ^{a b} Granger, R.A. (1995). Fluid Mechanics. Courier Dover Publications. ISBN 0486683567.  p. 30.
6. ^ Mellor, G.L. (1996). Introduction to Physical Oceanography. Springer. ISBN 1563962101.  p. 19.
7. ^ Emanuel, G. (2001). Analytical fluid dynamics (second ed.). CRC Press. ISBN 0849391148.  pp. 6–7.
8. ^ Eric W. Weisstein. "Convective Operator". MathWorld. http://mathworld.wolfram.com/ConvectiveOperator.html. Retrieved 2008-22-7. 

Further reading

• Sussman, G.J.; Wisdom, J.; Mayer, M.E.. "1.6 How to Find Lagrangians". Structure and Interpretation of Classical Mechanics. MIT Press. http://mitpress.mit.edu/SICM/book-Z-H-13.html#%_sec_Temp_122. 
• Cohen, Ira M.; Kundu, Pijush K. Fluid Mechanics (4 ed.). Academic Press. ISBN 978-0-12-373735-9. 
• Lai, Michael; Krempl, Erhard; Ruben, David. Introduction to Continuum Mechanics (4 ed.). Elsevier. ISBN 978-0-7506-8560-3.