Proofs Chapter 7: The (A·∇) Operator

Advection, Material Derivative, and Lamb's Identity

7. The (A·∇) Operator

This final chapter of the vector-calculus identity series treats the advection operator $(\mathbf{A} \cdot \nabla)$. The operator lies at the heart of the material derivative in continuum mechanics and appears directly as the nonlinear term of the Navier–Stokes equations.

Roadmap of this chapter

The proofs are organized as follows. Equation numbers follow the global numbering of the identity reference.

  • 7.1 Definition (58–59): the $(\mathbf{A} \cdot \nabla)$ operator and its component form
  • 7.2 Basic formulas (60–63): action on scalar fields, relation with the curl, gradient of a dot product
  • 7.3 Material derivative (64–65): time derivative along a fluid particle
  • 7.4 Lamb's identity (66): decomposition of $(\mathbf{v} \cdot \nabla)\mathbf{v}$ into kinetic-energy gradient and vortex term

7.1 Definition

(58) Definition of the (A·∇) operator

Definition: $\displaystyle (\mathbf{A} \cdot \nabla)\mathbf{B} = \left(A_x \dfrac{\partial}{\partial x} + A_y \dfrac{\partial}{\partial y} + A_z \dfrac{\partial}{\partial z}\right)\mathbf{B}$
Conditions: $\mathbf{A}, \mathbf{B}$ are differentiable vector fields
Discussion

The symbol $(\mathbf{A} \cdot \nabla)$ should be read as the inner product $\mathbf{A} \cdot \nabla$ regarded as a single operator. Formally substituting the components of $\nabla = (\partial/\partial x, \partial/\partial y, \partial/\partial z)$ and $\mathbf{A} = (A_x, A_y, A_z)$ gives

\begin{equation}\mathbf{A} \cdot \nabla = A_x \dfrac{\partial}{\partial x} + A_y \dfrac{\partial}{\partial y} + A_z \dfrac{\partial}{\partial z} \label{eq:7-1-1}\end{equation}

This is a scalar-valued first-order linear differential operator that can act on either a scalar field or a vector field.

Well-definedness. The right-hand side of $\eqref{eq:7-1-1}$ is a finite sum of linear differential operators, so its value is uniquely determined whenever the target $\mathbf{B}$ is differentiable. We define $(\mathbf{A} \cdot \nabla)\mathbf{B}$ as the vector field obtained by applying the scalar operator $\mathbf{A} \cdot \nabla$ component-by-component to $\mathbf{B}$.

Relation to the directional derivative. Writing $\hat{\mathbf{n}} = \mathbf{A}/|\mathbf{A}|$ for the unit vector along $\mathbf{A}$, we have

\begin{equation}(\mathbf{A} \cdot \nabla)\mathbf{B} = |\mathbf{A}| \, (\hat{\mathbf{n}} \cdot \nabla)\mathbf{B} \label{eq:7-1-2}\end{equation}

The quantity $(\hat{\mathbf{n}} \cdot \nabla)\mathbf{B}$ is the directional derivative of $\mathbf{B}$ along $\hat{\mathbf{n}}$, so $(\mathbf{A} \cdot \nabla)\mathbf{B}$ equals that directional derivative scaled by the magnitude $|\mathbf{A}|$.

Terminology: in general $(\mathbf{A} \cdot \nabla)$ is called the directional derivative operator. The name advection operator is normally reserved for the fluid-mechanical case in which $\mathbf{A}$ is the velocity field, i.e. for $(\mathbf{v} \cdot \nabla)$; one does not call $(\mathbf{E} \cdot \nabla)$ an advection operator. Note also that $(\mathbf{A} \cdot \nabla)$ is the name of the operator, whereas $\mathbf{A} \cdot \nabla f$ is the result of applying it to $f$.
Remark: The order of $\mathbf{A}$ and $\nabla$ matters. The combination $(\nabla \cdot \mathbf{A})$ is a scalar field (the divergence), whereas $(\mathbf{A} \cdot \nabla)$ is a differential operator; they are objects of completely different types.

(59) Component form of (A·∇)B

Identity: $\displaystyle \bigl[(\mathbf{A} \cdot \nabla)\mathbf{B}\bigr]_i = A_x \dfrac{\partial B_i}{\partial x} + A_y \dfrac{\partial B_i}{\partial y} + A_z \dfrac{\partial B_i}{\partial z}$
Conditions: $\mathbf{A}, \mathbf{B}$ differentiable vector fields, $i \in \{x, y, z\}$
Proof

Prerequisites: definition (58); the components of a vector field are differentiated independently.

From definition (58),

\begin{equation}(\mathbf{A} \cdot \nabla)\mathbf{B} = \left(A_x \dfrac{\partial}{\partial x} + A_y \dfrac{\partial}{\partial y} + A_z \dfrac{\partial}{\partial z}\right)\mathbf{B} \label{eq:7-2-1}\end{equation}

The right-hand side acts on the vector field $\mathbf{B} = (B_x, B_y, B_z)$, and its $i$-th component is therefore

\begin{equation}\bigl[(\mathbf{A} \cdot \nabla)\mathbf{B}\bigr]_i = \left(A_x \dfrac{\partial}{\partial x} + A_y \dfrac{\partial}{\partial y} + A_z \dfrac{\partial}{\partial z}\right) B_i \label{eq:7-2-2}\end{equation}

Linear differential operators distribute over their argument, so $\eqref{eq:7-2-2}$ expands to

\begin{equation}\bigl[(\mathbf{A} \cdot \nabla)\mathbf{B}\bigr]_i = A_x \dfrac{\partial B_i}{\partial x} + A_y \dfrac{\partial B_i}{\partial y} + A_z \dfrac{\partial B_i}{\partial z} \label{eq:7-2-3}\end{equation}

which is the desired identity. With the Einstein summation convention this can be written compactly as $\sum_j A_j \,\partial B_i/\partial x_j$.

Remark: The component form reduces $(\mathbf{A} \cdot \nabla)\mathbf{B}$ to three independent scalar directional derivatives. This is why the components of the Navier–Stokes equations can be discretized independently in numerical implementations.

7.2 Basic formulas

(60) Action on a scalar field

Identity: $(\mathbf{A} \cdot \nabla)f = \mathbf{A} \cdot \nabla f$
Conditions: $f$ is a differentiable scalar field, $\mathbf{A}$ a vector field
Proof

Prerequisites: definition (58); the gradient $\nabla f = (\partial f/\partial x, \partial f/\partial y, \partial f/\partial z)$; the definition of the inner product.

Expand the left-hand side using definition (58).

\begin{equation}(\mathbf{A} \cdot \nabla)f = \left(A_x \dfrac{\partial}{\partial x} + A_y \dfrac{\partial}{\partial y} + A_z \dfrac{\partial}{\partial z}\right) f = A_x \dfrac{\partial f}{\partial x} + A_y \dfrac{\partial f}{\partial y} + A_z \dfrac{\partial f}{\partial z} \label{eq:7-3-1}\end{equation}

The right-hand expression is exactly the dot product of $\mathbf{A} = (A_x, A_y, A_z)$ with $\nabla f = (\partial f/\partial x, \partial f/\partial y, \partial f/\partial z)$:

\begin{equation}A_x \dfrac{\partial f}{\partial x} + A_y \dfrac{\partial f}{\partial y} + A_z \dfrac{\partial f}{\partial z} = \mathbf{A} \cdot \nabla f \label{eq:7-3-2}\end{equation}

From $\eqref{eq:7-3-1}$ and $\eqref{eq:7-3-2}$ the identity $(\mathbf{A} \cdot \nabla)f = \mathbf{A} \cdot \nabla f$ follows.

Remark: For scalar fields the operator $(\mathbf{A} \cdot \nabla)$ and the expression $\mathbf{A} \cdot \nabla$ produce the same result. For vector fields, $(\mathbf{A} \cdot \nabla)\mathbf{B}$ and $\mathbf{A} \cdot (\nabla \mathbf{B})$ also denote the same components $\sum_j A_j \,\partial B_i/\partial x_j$, but should not be confused with combinations involving $\nabla\cdot$, e.g. $\nabla(\mathbf{A} \cdot \mathbf{B})$ or $\nabla \cdot (\mathbf{A}\mathbf{B})$.

(61) Cross product with the curl

Identity: $\mathbf{A} \times (\nabla \times \mathbf{B}) = (\nabla \mathbf{B}) \cdot \mathbf{A} - (\mathbf{A} \cdot \nabla)\mathbf{B}$
Conditions: $\mathbf{A}, \mathbf{B}$ differentiable vector fields. The Jacobian $(\nabla \mathbf{B})$ has components $(\nabla \mathbf{B})_{ij} = \partial B_j/\partial x_i$, and we adopt the convention $[(\nabla \mathbf{B}) \cdot \mathbf{A}]_i = \displaystyle\sum_j A_j \, \partial B_j/\partial x_i$.
Proof

Prerequisites: the Levi-Civita symbol $\varepsilon_{ijk}$; the $\varepsilon$-$\delta$ identity $\displaystyle\sum_k \varepsilon_{ijk}\varepsilon_{klm} = \delta_{il}\delta_{jm} - \delta_{im}\delta_{jl}$; the component formulas $[\mathbf{a} \times \mathbf{b}]_i = \varepsilon_{ijk} a_j b_k$ and $[\nabla \times \mathbf{B}]_k = \varepsilon_{klm} \,\partial B_m/\partial x_l$. Einstein summation is in force.

(i) Component expansion of the LHS. Expand the cross product in two stages.

\begin{equation}[\mathbf{A} \times (\nabla \times \mathbf{B})]_i = \varepsilon_{ijk} A_j \, [\nabla \times \mathbf{B}]_k = \varepsilon_{ijk} \varepsilon_{klm} A_j \dfrac{\partial B_m}{\partial x_l} \label{eq:7-4-1}\end{equation}

(ii) Apply the $\varepsilon$-$\delta$ identity. Substituting $\varepsilon_{ijk} \varepsilon_{klm} = \delta_{il}\delta_{jm} - \delta_{im}\delta_{jl}$ into $\eqref{eq:7-4-1}$ yields

\begin{equation}[\mathbf{A} \times (\nabla \times \mathbf{B})]_i = (\delta_{il}\delta_{jm} - \delta_{im}\delta_{jl}) A_j \dfrac{\partial B_m}{\partial x_l} \label{eq:7-4-2}\end{equation}

(iii) Contract the Kronecker deltas. $\delta_{il}$ forces $l = i$ and $\delta_{jm}$ forces $m = j$, so the first piece becomes

\begin{equation}\delta_{il}\delta_{jm} A_j \dfrac{\partial B_m}{\partial x_l} = A_j \dfrac{\partial B_j}{\partial x_i} \label{eq:7-4-3}\end{equation}

Similarly $\delta_{im}\delta_{jl}$ enforces $m = i$ and $l = j$, giving

\begin{equation}\delta_{im}\delta_{jl} A_j \dfrac{\partial B_m}{\partial x_l} = A_j \dfrac{\partial B_i}{\partial x_j} \label{eq:7-4-4}\end{equation}

(iv) Combine. From $\eqref{eq:7-4-2}$, $\eqref{eq:7-4-3}$, $\eqref{eq:7-4-4}$:

\begin{equation}[\mathbf{A} \times (\nabla \times \mathbf{B})]_i = A_j \dfrac{\partial B_j}{\partial x_i} - A_j \dfrac{\partial B_i}{\partial x_j} \label{eq:7-4-5}\end{equation}

(v) Identify each piece. The first sum $\displaystyle\sum_j A_j \,\partial B_j/\partial x_i$ is exactly $[(\nabla \mathbf{B}) \cdot \mathbf{A}]_i$ by the convention stated in the conditions. The second sum $\displaystyle\sum_j A_j \,\partial B_i/\partial x_j$ equals $[(\mathbf{A} \cdot \nabla)\mathbf{B}]_i$ by (59). Hence

\begin{equation}[\mathbf{A} \times (\nabla \times \mathbf{B})]_i = [(\nabla \mathbf{B}) \cdot \mathbf{A}]_i - [(\mathbf{A} \cdot \nabla)\mathbf{B}]_i \label{eq:7-4-6}\end{equation}

Since $i$ was arbitrary, we obtain the vector identity $\mathbf{A} \times (\nabla \times \mathbf{B}) = (\nabla \mathbf{B}) \cdot \mathbf{A} - (\mathbf{A} \cdot \nabla)\mathbf{B}$, as required.

Remark: The exact meaning of $(\nabla \mathbf{B}) \cdot \mathbf{A}$ depends on convention. We use $[(\nabla \mathbf{B}) \cdot \mathbf{A}]_i = \displaystyle\sum_j A_j \,\partial B_j/\partial x_i$, which is the right-multiplication of the Jacobian by $\mathbf{A}$. The matrix used on this page is $$\nabla\mathbf{B} = \begin{pmatrix} \partial_x B_x & \partial_x B_y & \partial_x B_z \\ \partial_y B_x & \partial_y B_y & \partial_y B_z \\ \partial_z B_x & \partial_z B_y & \partial_z B_z \end{pmatrix}$$ whose $(i,j)$ entry is $\partial B_j/\partial x_i$ (some texts use the transpose, so always check when comparing with other sources). Thus $(\nabla\mathbf{B}) \cdot \mathbf{A}$ is simply the linear map defined by this matrix applied to $\mathbf{A}$.

(62) Component form of A×(∇×B)

Identity: $\displaystyle \bigl[\mathbf{A} \times (\nabla \times \mathbf{B})\bigr]_i = \displaystyle\sum_j A_j \dfrac{\partial B_j}{\partial x_i} - \displaystyle\sum_j A_j \dfrac{\partial B_i}{\partial x_j}$
Conditions: $\mathbf{A}, \mathbf{B}$ differentiable vector fields
Proof

Prerequisites: equation $\eqref{eq:7-4-5}$ from the proof of (61).

Steps (i)–(iv) of the proof of (61) produced

\begin{equation}[\mathbf{A} \times (\nabla \times \mathbf{B})]_i = \displaystyle\sum_j A_j \dfrac{\partial B_j}{\partial x_i} - \displaystyle\sum_j A_j \dfrac{\partial B_i}{\partial x_j} \label{eq:7-5-1}\end{equation}

which is precisely identity (62). In short, (62) is the component form of the vector identity (61).

Remark: Identity (62) is not independent of (61); it is its component restatement. For explicit calculations the component form is often the more convenient one.

(63) Alternate form of the gradient of a dot product

Identity: $(\mathbf{A} \cdot \nabla)\mathbf{B} + (\mathbf{B} \cdot \nabla)\mathbf{A} = \nabla(\mathbf{A} \cdot \mathbf{B}) - \mathbf{A} \times (\nabla \times \mathbf{B}) - \mathbf{B} \times (\nabla \times \mathbf{A})$
Conditions: $\mathbf{A}, \mathbf{B}$ differentiable vector fields
Proof

Prerequisites: the gradient-of-dot-product formula (16): $\nabla(\mathbf{A} \cdot \mathbf{B}) = (\mathbf{B} \cdot \nabla)\mathbf{A} + (\mathbf{A} \cdot \nabla)\mathbf{B} + \mathbf{B} \times (\nabla \times \mathbf{A}) + \mathbf{A} \times (\nabla \times \mathbf{B})$, which we take as already established.

Start from (16):

\begin{equation}\nabla(\mathbf{A} \cdot \mathbf{B}) = (\mathbf{B} \cdot \nabla)\mathbf{A} + (\mathbf{A} \cdot \nabla)\mathbf{B} + \mathbf{B} \times (\nabla \times \mathbf{A}) + \mathbf{A} \times (\nabla \times \mathbf{B}) \label{eq:7-6-1}\end{equation}

Subtracting $\mathbf{A} \times (\nabla \times \mathbf{B}) + \mathbf{B} \times (\nabla \times \mathbf{A})$ from both sides gives

\begin{equation}\nabla(\mathbf{A} \cdot \mathbf{B}) - \mathbf{A} \times (\nabla \times \mathbf{B}) - \mathbf{B} \times (\nabla \times \mathbf{A}) = (\mathbf{A} \cdot \nabla)\mathbf{B} + (\mathbf{B} \cdot \nabla)\mathbf{A} \label{eq:7-6-2}\end{equation}

which is exactly (63).

Remark: Identity (63) is just (16) solved for $(\mathbf{A} \cdot \nabla)\mathbf{B} + (\mathbf{B} \cdot \nabla)\mathbf{A}$. It expresses advection in terms of a gradient and two vortex terms. Setting $\mathbf{A} = \mathbf{B}$ produces Lamb's identity (66).

7.3 Material derivative

In continuum mechanics, the rate of change observed while travelling with a fluid particle is called the material derivative (or the Lagrangian derivative) and is denoted $D/Dt$. It differs from the partial time derivative $\partial/\partial t$ measured at a fixed point in space.

(64) Material derivative of a scalar field

Identity: $\dfrac{Df}{Dt} = \dfrac{\partial f}{\partial t} + (\mathbf{v} \cdot \nabla)f$
Conditions: $f(\mathbf{x}, t)$ is a differentiable scalar field; $\mathbf{v}(\mathbf{x}, t)$ is the velocity field of the fluid
Proof

Prerequisites: a fluid-particle trajectory $\mathbf{x}(t) = (x(t), y(t), z(t))$ satisfies $d\mathbf{x}/dt = \mathbf{v}(\mathbf{x}(t), t)$; the multivariate chain rule; identity (60).

(i) Definition of the material derivative. The material derivative $Df/Dt$ is defined as the total derivative of the composite function $f(\mathbf{x}(t), t)$ along a fluid trajectory $\mathbf{x}(t)$:

\begin{equation}\dfrac{Df}{Dt} = \dfrac{d}{dt} f(\mathbf{x}(t), t) \label{eq:7-7-1}\end{equation}

(ii) Multivariate chain rule. Since $f$ is a function of the four variables $x, y, z, t$ and $x, y, z$ themselves depend on $t$, the chain rule gives

\begin{equation}\dfrac{d}{dt} f(\mathbf{x}(t), t) = \dfrac{\partial f}{\partial t} + \dfrac{\partial f}{\partial x} \dfrac{dx}{dt} + \dfrac{\partial f}{\partial y} \dfrac{dy}{dt} + \dfrac{\partial f}{\partial z} \dfrac{dz}{dt} \label{eq:7-7-2}\end{equation}

(iii) Substitute the velocity. By definition of the trajectory, $dx/dt = v_x$, $dy/dt = v_y$, $dz/dt = v_z$, so

\begin{equation}\dfrac{d}{dt} f(\mathbf{x}(t), t) = \dfrac{\partial f}{\partial t} + v_x \dfrac{\partial f}{\partial x} + v_y \dfrac{\partial f}{\partial y} + v_z \dfrac{\partial f}{\partial z} \label{eq:7-7-3}\end{equation}

(iv) Identify the advective term. The trailing terms of $\eqref{eq:7-7-3}$ form the dot product $\mathbf{v} \cdot \nabla f$, which by (60) equals $(\mathbf{v} \cdot \nabla)f$. Therefore

\begin{equation}\dfrac{Df}{Dt} = \dfrac{\partial f}{\partial t} + (\mathbf{v} \cdot \nabla)f \label{eq:7-7-4}\end{equation}

as claimed.

Remark: $\partial f/\partial t$ measures the time variation at a fixed Eulerian point, while $Df/Dt$ measures it along a moving fluid particle (Lagrangian description). The difference $(\mathbf{v} \cdot \nabla)f$ is the advective term and represents the change of $f$ caused by the motion of the particle.

(65) Material derivative of a vector field

Identity: $\dfrac{D\mathbf{A}}{Dt} = \dfrac{\partial \mathbf{A}}{\partial t} + (\mathbf{v} \cdot \nabla)\mathbf{A}$
Conditions: $\mathbf{A}(\mathbf{x}, t)$ is a differentiable vector field; $\mathbf{v}(\mathbf{x}, t)$ is the velocity field
Proof

Prerequisites: the scalar material derivative (64); the components of a vector field have independent material derivatives; definition (58).

Each component $A_i$ of $\mathbf{A} = (A_x, A_y, A_z)$ is a scalar field on space-time, so identity (64) applies to it:

\begin{equation}\dfrac{DA_i}{Dt} = \dfrac{\partial A_i}{\partial t} + (\mathbf{v} \cdot \nabla) A_i \quad (i = x, y, z) \label{eq:7-8-1}\end{equation}

Stacking these three scalar equations into a single vector equation:

\begin{equation}\dfrac{D\mathbf{A}}{Dt} = \begin{pmatrix} DA_x/Dt \\ DA_y/Dt \\ DA_z/Dt \end{pmatrix} = \begin{pmatrix} \partial A_x/\partial t \\ \partial A_y/\partial t \\ \partial A_z/\partial t \end{pmatrix} + \begin{pmatrix} (\mathbf{v} \cdot \nabla) A_x \\ (\mathbf{v} \cdot \nabla) A_y \\ (\mathbf{v} \cdot \nabla) A_z \end{pmatrix} \label{eq:7-8-2}\end{equation}

The first vector on the right is $\partial \mathbf{A}/\partial t$, and the second one is $(\mathbf{v} \cdot \nabla)\mathbf{A}$ by the definition (58). Thus

\begin{equation}\dfrac{D\mathbf{A}}{Dt} = \dfrac{\partial \mathbf{A}}{\partial t} + (\mathbf{v} \cdot \nabla)\mathbf{A} \label{eq:7-8-3}\end{equation}

as required.

Remark: Setting $\mathbf{A} = \mathbf{v}$ gives $D\mathbf{v}/Dt = \partial \mathbf{v}/\partial t + (\mathbf{v} \cdot \nabla)\mathbf{v}$, which appears on the left-hand side of Newton's second law for a fluid particle, $\rho \, D\mathbf{v}/Dt = \mathbf{f}$. The nonlinearity $(\mathbf{v} \cdot \nabla)\mathbf{v}$ in the Navier–Stokes equations originates here.

7.4 Lamb's identity

(66) Lamb's identity

Identity: $(\mathbf{v} \cdot \nabla)\mathbf{v} = \nabla\!\left(\dfrac{|\mathbf{v}|^2}{2}\right) - \mathbf{v} \times (\nabla \times \mathbf{v})$
Conditions: $\mathbf{v}$ is a differentiable vector field
Proof

Prerequisites: the alternate form of the gradient of a dot product (63); $\mathbf{v} \cdot \mathbf{v} = |\mathbf{v}|^2$.

(i) Set $\mathbf{A} = \mathbf{B} = \mathbf{v}$ in (63):

\begin{equation}(\mathbf{v} \cdot \nabla)\mathbf{v} + (\mathbf{v} \cdot \nabla)\mathbf{v} = \nabla(\mathbf{v} \cdot \mathbf{v}) - \mathbf{v} \times (\nabla \times \mathbf{v}) - \mathbf{v} \times (\nabla \times \mathbf{v}) \label{eq:7-9-1}\end{equation}

(ii) Combine like terms:

\begin{equation}2 (\mathbf{v} \cdot \nabla)\mathbf{v} = \nabla(\mathbf{v} \cdot \mathbf{v}) - 2 \, \mathbf{v} \times (\nabla \times \mathbf{v}) \label{eq:7-9-2}\end{equation}

(iii) Divide both sides by 2:

\begin{equation}(\mathbf{v} \cdot \nabla)\mathbf{v} = \dfrac{1}{2} \nabla(\mathbf{v} \cdot \mathbf{v}) - \mathbf{v} \times (\nabla \times \mathbf{v}) \label{eq:7-9-3}\end{equation}

(iv) Substitute $\mathbf{v} \cdot \mathbf{v} = |\mathbf{v}|^2$. By linearity of the gradient $\dfrac{1}{2} \nabla(\mathbf{v} \cdot \mathbf{v}) = \nabla\!\left(\dfrac{|\mathbf{v}|^2}{2}\right)$, so

\begin{equation}(\mathbf{v} \cdot \nabla)\mathbf{v} = \nabla\!\left(\dfrac{|\mathbf{v}|^2}{2}\right) - \mathbf{v} \times (\nabla \times \mathbf{v}) \label{eq:7-9-4}\end{equation}

as desired.

Lamb identity: the advective acceleration splits into the kinetic-energy gradient and the Lamb vector streamline ω = ∇×v P v ( v · ∇ ) v ∇( |v|² / 2 ) − v × ω advective acceleration (black) = energy gradient (green) + vorticity term −v × ω (orange)
Figure 1. Geometric meaning of Lamb's identity. The velocity $\mathbf{v}$ (blue) is tangent to the streamline, and the orange arc shows the sense of rotation carried by the vorticity (here $\boldsymbol{\omega}$ points into the page, so the rotation appears clockwise). The advective acceleration $(\mathbf{v}\cdot\nabla)\mathbf{v}$ (black) of a fluid particle $P$ splits into the vector sum of the kinetic-energy gradient $\nabla(|\mathbf{v}|^2/2)$ (green) and the vorticity term $-\mathbf{v}\times\boldsymbol{\omega}$ (orange, with $\boldsymbol{\omega}=\nabla\times\mathbf{v}$). The orange arrow is drawn head-to-tail from the tip of the green one, so the black arrow closes the triangle. Being a cross product, $-\mathbf{v}\times\boldsymbol{\omega}$ is perpendicular to $\mathbf{v}$.
Remark: The advective term $(\mathbf{v} \cdot \nabla)\mathbf{v}$ is split into the kinetic-energy gradient $\nabla(|\mathbf{v}|^2/2)$ and the vortex term (Lamb vector) $-\mathbf{v} \times (\nabla \times \mathbf{v}) = -\mathbf{v} \times \boldsymbol{\omega}$ with vorticity $\boldsymbol{\omega} = \nabla \times \mathbf{v}$. Combined with the Euler equation in steady, inviscid, incompressible flow it yields Bernoulli's theorem $|\mathbf{v}|^2/2 + p/\rho + g h = \text{const}$ along streamlines or vortex lines. Lamb's identity is one of the most-used identities in fluid dynamics, magnetohydrodynamics, and meteorology.

References

  • Batchelor, G. K. (2000). An Introduction to Fluid Dynamics. Cambridge University Press. (Section 3.5: Bernoulli's theorem and the Lamb vector)
  • Lamb, H. (1932). Hydrodynamics (6th ed.). Cambridge University Press. (the original source of Lamb's identity)
  • Arfken, G. B., Weber, H. J., & Harris, F. E. (2013). Mathematical Methods for Physicists (7th ed.). Academic Press. (Chapter 3: vector analysis)
  • Vector calculus identities - Wikipedia
  • Material derivative - Wikipedia

Frequently Asked Questions

What does the (A·∇) operator mean?

$(\mathbf{A}\cdot\nabla)$ is the directional derivative operator with direction vector $\mathbf{A}$. It acts on both scalar fields and vector fields and physically measures the rate of change along the direction $\mathbf{A}$. It is the unnormalised directional derivative: the ordinary directional derivative along the unit vector $\hat{\mathbf{A}}$, scaled by $\lVert\mathbf{A}\rVert$.

What is the material derivative D/Dt?

It is the rate of change observed while moving with a fluid particle. It equals the partial time derivative $\partial/\partial t$ in fixed Eulerian coordinates plus the advective term $(\mathbf{v}\cdot\nabla)$, giving $D/Dt = \partial/\partial t + (\mathbf{v}\cdot\nabla)$.

Where is Lamb's identity used?

It decomposes $(\mathbf{v}\cdot\nabla)\mathbf{v}$ into the kinetic-energy gradient $\nabla(|\mathbf{v}|^2/2)$ and the vortex term $-\mathbf{v}\times(\nabla\times\mathbf{v})$. It is essential for deriving Bernoulli's theorem in inviscid flow and for analysing the momentum term of the Navier–Stokes equations.