Gradient, Divergence, and Curl

In vector calculus, we often deal with scalar fields and vector fields. A scalar field assigns a scalar value to each point in space, like temperature or pressure. A vector field assigns a vector to each point in space, like velocity or force. Gradient, divergence, and curl are fundamental operations that describe how these fields change in space. They are crucial for understanding phenomena in physics and engineering, such as fluid dynamics, electromagnetism, and heat transfer.

1. Gradient

The gradient is an operation that applies to a scalar field and produces a vector field. It tells us the direction and magnitude of the steepest ascent of the scalar field at any given point.

1.1 Definition of Gradient

Let $f(x, y, z)$ be a scalar field defined in a region of space. The gradient of $f$, denoted by $\nabla f$ or grad $f$, is a vector field defined as:

$\nabla f = \frac{\partial f}{\partial x} \mathbf{i} + \frac{\partial f}{\partial y} \mathbf{j} + \frac{\partial f}{\partial z} \mathbf{k}$

where $\mathbf{i}$, $\mathbf{j}$, and $\mathbf{k}$ are the unit vectors along the x, y, and z axes, respectively. The symbol $\nabla$ (nabla or del) is a vector differential operator. In Cartesian coordinates, it is defined as $\nabla = \frac{\partial}{\partial x} \mathbf{i} + \frac{\partial}{\partial y} \mathbf{j} + \frac{\partial}{\partial z} \mathbf{k}$.

1.2 Properties of Gradient

• The gradient vector $\nabla f$ at a point $P$ points in the direction of the greatest rate of increase of $f$ at $P$.
• The magnitude $|\nabla f|$ at $P$ is the maximum rate of increase of $f$ at $P$.
• The gradient is perpendicular to the level surfaces of the scalar field. A level surface is a surface where $f(x, y, z) = c$ for some constant $c$.
• The gradient satisfies certain linearity properties:
  • $\nabla (f+g) = \nabla f + \nabla g$
  • $\nabla (cf) = c \nabla f$, where $c$ is a constant.
  • $\nabla (fg) = f \nabla g + g \nabla f$ (Product rule)

1.3 Example of Gradient

Consider the scalar field representing the temperature in a room: $T(x, y, z) = x^2 + y^2 + z^2$. We want to find the direction and rate of the steepest temperature increase at the point $(1, 2, 3)$.

First, we find the gradient of $T$:

$\nabla T = \frac{\partial T}{\partial x} \mathbf{i} + \frac{\partial T}{\partial y} \mathbf{j} + \frac{\partial T}{\partial z} \mathbf{k}$

$\frac{\partial T}{\partial x} = 2x$

$\frac{\partial T}{\partial y} = 2y$

$\frac{\partial T}{\partial z} = 2z$

So, $\nabla T = 2x \mathbf{i} + 2y \mathbf{j} + 2z \mathbf{k}$.

Now, we evaluate the gradient at the point $(1, 2, 3)$:

$\nabla T(1, 2, 3) = 2(1) \mathbf{i} + 2(2) \mathbf{j} + 2(3) \mathbf{k} = 2\mathbf{i} + 4\mathbf{j} + 6\mathbf{k}$.

This vector indicates that the temperature increases most rapidly in the direction $(2, 4, 6)$ at the point $(1, 2, 3)$. The rate of this increase is the magnitude of the gradient:

$|\nabla T(1, 2, 3)| = \sqrt{2^2 + 4^2 + 6^2} = \sqrt{4 + 16 + 36} = \sqrt{56} = 2\sqrt{14}$.

The level surfaces for this temperature field are spheres centered at the origin ($x^2 + y^2 + z^2 = c$). The gradient vector at any point on a sphere is radial, pointing outward, which is perpendicular to the sphere's surface.

2. Divergence

Divergence is an operation that applies to a vector field and produces a scalar field. It measures the extent to which a vector field flows outward from a point. It can be thought of as the "source strength" or "sink strength" of the field at that point.

2.1 Definition of Divergence

Let $\mathbf{F}(x, y, z) = F_1(x, y, z) \mathbf{i} + F_2(x, y, z) \mathbf{j} + F_3(x, y, z) \mathbf{k}$ be a vector field. The divergence of $\mathbf{F}$, denoted by $\nabla \cdot \mathbf{F}$ or div $\mathbf{F}$, is a scalar field defined as:

$\nabla \cdot \mathbf{F} = \frac{\partial F_1}{\partial x} + \frac{\partial F_2}{\partial y} + \frac{\partial F_3}{\partial z}$

The dot in $\nabla \cdot \mathbf{F}$ is suggestive of the dot product between the vector operator $\nabla$ and the vector field $\mathbf{F}$.

2.2 Interpretation of Divergence

• If $\nabla \cdot \mathbf{F} > 0$ at a point, the point acts as a source for the field, meaning there is a net outward flow from that point.
• If $\nabla \cdot \mathbf{F} < 0$ at a point, the point acts as a sink for the field, meaning there is a net inward flow towards that point.
• If $\nabla \cdot \mathbf{F} = 0$ at a point, the field is said to be solenoidal or divergence-free at that point. This means that the net flow into or out of any small volume around the point is zero. For example, the magnetic field in Maxwell's equations is divergence-free ($\nabla \cdot \mathbf{B} = 0$).

2.3 Properties of Divergence

• Divergence is a linear operator:
  • $\nabla \cdot (\mathbf{F} + \mathbf{G}) = \nabla \cdot \mathbf{F} + \nabla \cdot \mathbf{G}$
  • $\nabla \cdot (c\mathbf{F}) = c (\nabla \cdot \mathbf{F})$, where $c$ is a constant.
• Product rule for divergence:
  • $\nabla \cdot (f\mathbf{F}) = (\nabla f) \cdot \mathbf{F} + f (\nabla \cdot \mathbf{F})$, where $f$ is a scalar field.

2.4 Example of Divergence

Consider the vector field representing fluid velocity: $\mathbf{v}(x, y, z) = x\mathbf{i} + y\mathbf{j} + z\mathbf{k}$. We want to find the divergence of this field.

Here, $F_1 = x$, $F_2 = y$, and $F_3 = z$.

$\nabla \cdot \mathbf{v} = \frac{\partial (x)}{\partial x} + \frac{\partial (y)}{\partial y} + \frac{\partial (z)}{\partial z}$

$\frac{\partial (x)}{\partial x} = 1$

$\frac{\partial (y)}{\partial y} = 1$

$\frac{\partial (z)}{\partial z} = 1$

So, $\nabla \cdot \mathbf{v} = 1 + 1 + 1 = 3$.

Since the divergence is a positive constant (3), this vector field represents a fluid that is expanding uniformly in all directions from every point. Every point acts as a source. This is a pure source field.

Now consider another vector field: $\mathbf{F}(x, y, z) = x^2\mathbf{i} + y^2\mathbf{j} + z^2\mathbf{k}$.

Here, $F_1 = x^2$, $F_2 = y^2$, and $F_3 = z^2$.

$\nabla \cdot \mathbf{F} = \frac{\partial (x^2)}{\partial x} + \frac{\partial (y^2)}{\partial y} + \frac{\partial (z^2)}{\partial z}$

$\frac{\partial (x^2)}{\partial x} = 2x$

$\frac{\partial (y^2)}{\partial y} = 2y$

$\frac{\partial (z^2)}{\partial z} = 2z$

So, $\nabla \cdot \mathbf{F} = 2x + 2y + 2z$.

The divergence is not constant. At points where $x+y+z > 0$, the divergence is positive, indicating a source. At points where $x+y+z < 0$, the divergence is negative, indicating a sink.

3. Curl

Curl is an operation that applies to a vector field and produces another vector field. It measures the tendency of the vector field to rotate or circulate around a point.

3.1 Definition of Curl

Let $\mathbf{F}(x, y, z) = F_1(x, y, z) \mathbf{i} + F_2(x, y, z) \mathbf{j} + F_3(x, y, z) \mathbf{k}$ be a vector field. The curl of $\mathbf{F}$, denoted by $\nabla \times \mathbf{F}$ or curl $\mathbf{F}$, is a vector field defined as:

$\nabla \times \mathbf{F} = \left( \frac{\partial F_3}{\partial y} - \frac{\partial F_2}{\partial z} \right) \mathbf{i} + \left( \frac{\partial F_1}{\partial z} - \frac{\partial F_3}{\partial x} \right) \mathbf{j} + \left( \frac{\partial F_2}{\partial x} - \frac{\partial F_1}{\partial y} \right) \mathbf{k}$

This can be remembered using the determinant form:

$$ \nabla \times \mathbf{F} = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\ F_1 & F_2 & F_3 \end{vmatrix} $$

3.2 Interpretation of Curl

• If $\nabla \times \mathbf{F} \neq \mathbf{0}$ at a point, the field has a rotational component at that point. The direction of the curl vector indicates the axis of rotation (using the right-hand rule), and its magnitude indicates the strength of the rotation.
• If $\nabla \times \mathbf{F} = \mathbf{0}$ at a point, the field is said to be irrotational at that point. This means there is no tendency for circulation. For example, the electric field in electrostatics is irrotational ($\nabla \times \mathbf{E} = \mathbf{0}$).

3.3 Properties of Curl

• Curl is a linear operator:
  • $\nabla \times (\mathbf{F} + \mathbf{G}) = \nabla \times \mathbf{F} + \nabla \times \mathbf{G}$
  • $\nabla \times (c\mathbf{F}) = c (\nabla \times \mathbf{F})$, where $c$ is a constant.
• Curl of a gradient is always zero: $\nabla \times (\nabla f) = \mathbf{0}$. This is because any scalar field is conservative and thus irrotational.
• Divergence of a curl is always zero: $\nabla \cdot (\nabla \times \mathbf{F}) = 0$. This means that any vector field that can be expressed as the curl of another vector field must be divergence-free (solenoidal).
• Product rule for curl:
  • $\nabla \times (f\mathbf{F}) = (\nabla f) \times \mathbf{F} + f (\nabla \times \mathbf{F})$, where $f$ is a scalar field.

3.4 Example of Curl

Consider the vector field $\mathbf{F}(x, y, z) = -y\mathbf{i} + x\mathbf{j}$. This field represents a rotation around the z-axis. Let's calculate its curl.

Here, $F_1 = -y$, $F_2 = x$, $F_3 = 0$.

$\frac{\partial F_3}{\partial y} = \frac{\partial (0)}{\partial y} = 0$

$\frac{\partial F_2}{\partial z} = \frac{\partial (x)}{\partial z} = 0$

$\frac{\partial F_1}{\partial z} = \frac{\partial (-y)}{\partial z} = 0$

$\frac{\partial F_3}{\partial x} = \frac{\partial (0)}{\partial x} = 0$

$\frac{\partial F_2}{\partial x} = \frac{\partial (x)}{\partial x} = 1$

$\frac{\partial F_1}{\partial y} = \frac{\partial (-y)}{\partial y} = -1$

Now, substitute these into the curl formula:

$\nabla \times \mathbf{F} = (0 - 0)\mathbf{i} + (0 - 0)\mathbf{j} + (1 - (-1))\mathbf{k}$

$\nabla \times \mathbf{F} = 0\mathbf{i} + 0\mathbf{j} + 2\mathbf{k} = 2\mathbf{k}$.

The curl is $2\mathbf{k}$. This vector points along the positive z-axis. The magnitude of 2 indicates the strength of the rotation. This result confirms that the field $\mathbf{F} = -y\mathbf{i} + x\mathbf{j}$ is indeed rotating. The rotation is counter-clockwise in the xy-plane, consistent with the right-hand rule applied to the $\mathbf{k}$ direction.

Consider the vector field $\mathbf{G}(x, y, z) = xy\mathbf{i} + yz\mathbf{j} + zx\mathbf{k}$.

$F_1 = xy$, $F_2 = yz$, $F_3 = zx$.

$\frac{\partial F_3}{\partial y} = \frac{\partial (zx)}{\partial y} = 0$

$\frac{\partial F_2}{\partial z} = \frac{\partial (yz)}{\partial z} = y$

$\frac{\partial F_1}{\partial z} = \frac{\partial (xy)}{\partial z} = 0$

$\frac{\partial F_3}{\partial x} = \frac{\partial (zx)}{\partial x} = z$

$\frac{\partial F_2}{\partial x} = \frac{\partial (yz)}{\partial x} = 0$

$\frac{\partial F_1}{\partial y} = \frac{\partial (xy)}{\partial y} = x$

$\nabla \times \mathbf{G} = (0 - y)\mathbf{i} + (0 - z)\mathbf{j} + (0 - x)\mathbf{k} = -y\mathbf{i} - z\mathbf{j} - x\mathbf{k}$.

This field has a non-zero curl, indicating it is rotational.

4. Relationships between Gradient, Divergence, and Curl

These three operators are interconnected, and their relationships are fundamental in vector calculus and its applications.

4.1 Key Identities

Let $f$ be a scalar field and $\mathbf{F}$, $\mathbf{G}$ be vector fields.

1. Divergence of a gradient: $\nabla \cdot (\nabla f)$ is the Laplacian of $f$, denoted $\nabla^2 f$. $$ \nabla^2 f = \frac{\partial^2 f}{\partial x^2} + \frac{\partial^2 f}{\partial y^2} + \frac{\partial^2 f}{\partial z^2} $$ The Laplacian appears in the heat equation, wave equation, and Poisson's equation.
2. Curl of a gradient: $\nabla \times (\nabla f) = \mathbf{0}$. This implies that any irrotational vector field (a field whose curl is zero) can be expressed as the gradient of a scalar potential.
3. Divergence of a curl: $\nabla \cdot (\nabla \times \mathbf{F}) = 0$. This implies that any solenoidal vector field (a field whose divergence is zero) can be expressed as the curl of a vector potential.
4. Curl of a curl: $\nabla \times (\nabla \times \mathbf{F}) = \nabla (\nabla \cdot \mathbf{F}) - \nabla^2 \mathbf{F}$. The term $\nabla^2 \mathbf{F}$ is the vector Laplacian, defined as $\nabla^2 \mathbf{F} = \nabla^2 F_1 \mathbf{i} + \nabla^2 F_2 \mathbf{j} + \nabla^2 F_3 \mathbf{k}$. This identity is crucial in deriving the wave equation for electromagnetic fields.

4.2 Examples of Identities

Identity 1: Curl of a gradient Let $f(x, y, z) = x^2y + y^2z + z^2x$. First, find $\nabla f$: $\frac{\partial f}{\partial x} = 2xy + z^2$ $\frac{\partial f}{\partial y} = x^2 + 2yz$ $\frac{\partial f}{\partial z} = y^2 + 2zx$ So, $\nabla f = (2xy + z^2)\mathbf{i} + (x^2 + 2yz)\mathbf{j} + (y^2 + 2zx)\mathbf{k}$. Now, find the curl of $\nabla f$. Let $F_1 = 2xy + z^2$, $F_2 = x^2 + 2yz$, $F_3 = y^2 + 2zx$. $\nabla \times (\nabla f) = \left( \frac{\partial F_3}{\partial y} - \frac{\partial F_2}{\partial z} \right) \mathbf{i} + \left( \frac{\partial F_1}{\partial z} - \frac{\partial F_3}{\partial x} \right) \mathbf{j} + \left( \frac{\partial F_2}{\partial x} - \frac{\partial F_1}{\partial y} \right) \mathbf{k}$ $\frac{\partial F_3}{\partial y} = \frac{\partial (y^2 + 2zx)}{\partial y} = 2y$ $\frac{\partial F_2}{\partial z} = \frac{\partial (x^2 + 2yz)}{\partial z} = 2y$ $\frac{\partial F_1}{\partial z} = \frac{\partial (2xy + z^2)}{\partial z} = 2z$ $\frac{\partial F_3}{\partial x} = \frac{\partial (y^2 + 2zx)}{\partial x} = 2z$ $\frac{\partial F_2}{\partial x} = \frac{\partial (x^2 + 2yz)}{\partial x} = 2x$ $\frac{\partial F_1}{\partial y} = \frac{\partial (2xy + z^2)}{\partial y} = 2x$ Substituting these values: $\nabla \times (\nabla f) = (2y - 2y)\mathbf{i} + (2z - 2z)\mathbf{j} + (2x - 2x)\mathbf{k} = 0\mathbf{i} + 0\mathbf{j} + 0\mathbf{k} = \mathbf{0}$. This confirms the identity $\nabla \times (\nabla f) = \mathbf{0}$.

Identity 2: Divergence of a curl Let $\mathbf{F}(x, y, z) = yz\mathbf{i} + zx\mathbf{j} + xy\mathbf{k}$. First, find $\nabla \times \mathbf{F}$: $F_1 = yz$, $F_2 = zx$, $F_3 = xy$. $\frac{\partial F_3}{\partial y} = \frac{\partial (xy)}{\partial y} = x$ $\frac{\partial F_2}{\partial z} = \frac{\partial (zx)}{\partial z} = x$ $\frac{\partial F_1}{\partial z} = \frac{\partial (yz)}{\partial z} = y$ $\frac{\partial F_3}{\partial x} = \frac{\partial (xy)}{\partial x} = y$ $\frac{\partial F_2}{\partial x} = \frac{\partial (zx)}{\partial x} = z$ $\frac{\partial F_1}{\partial y} = \frac{\partial (yz)}{\partial y} = z$ $\nabla \times \mathbf{F} = (x-x)\mathbf{i} + (y-y)\mathbf{j} + (z-z)\mathbf{k} = 0\mathbf{i} + 0\mathbf{j} + 0\mathbf{k} = \mathbf{0}$. Now, find the divergence of $\nabla \times \mathbf{F}$. Since $\nabla \times \mathbf{F}$ is the zero vector, its divergence is: $\nabla \cdot (\mathbf{0}) = \frac{\partial (0)}{\partial x} + \frac{\partial (0)}{\partial y} + \frac{\partial (0)}{\partial z} = 0$. This confirms the identity $\nabla \cdot (\nabla \times \mathbf{F}) = 0$.

5. Applications in Physics and Engineering

Gradient, divergence, and curl are indispensable tools in various scientific disciplines.

5.1 Electromagnetism

Maxwell's equations, which govern electromagnetic fields, are expressed using these operators:

• Gauss's Law for Electricity: $\nabla \cdot \mathbf{E} = \frac{\rho}{\epsilon_0}$ (The divergence of the electric field is proportional to the charge density. Charges are sources/sinks of the electric field).
• Gauss's Law for Magnetism: $\nabla \cdot \mathbf{B} = 0$ (The divergence of the magnetic field is zero. There are no magnetic monopoles; magnetic field lines form closed loops).
• Faraday's Law of Induction: $\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}$ (A changing magnetic field induces an electric field, causing circulation).
• Ampère-Maxwell Law: $\nabla \times \mathbf{H} = \mathbf{J} + \frac{\partial \mathbf{D}}{\partial t}$ (Electric currents and changing electric fields produce magnetic fields, causing circulation).

The concept of scalar and vector potentials is also heavily reliant on these operators. For instance, an irrotational field like the electrostatic field can be written as $\mathbf{E} = -\nabla V$, where $V$ is the scalar electric potential. A solenoidal field like the magnetic field can be written as $\mathbf{B} = \nabla \times \mathbf{A}$, where $\mathbf{A}$ is the vector magnetic potential.

5.2 Fluid Dynamics

The behavior of fluids is described using these operators:

• Continuity Equation (Conservation of Mass): $\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0$, where $\rho$ is density and $\mathbf{v}$ is velocity. For an incompressible fluid ($\rho$ = constant), this simplifies to $\nabla \cdot \mathbf{v} = 0$, meaning the velocity field is solenoidal.
• Navier-Stokes Equations: These complex equations describe the motion of viscous fluid substances and involve terms with divergence and curl of the velocity field, as well as pressure gradients (gradients of a scalar field).

5.3 Heat Transfer

The rate of heat flow is often described by a heat flux vector $\mathbf{q}$. Fourier's Law states that $\mathbf{q} = -k \nabla T$, where $k$ is the thermal conductivity and $T$ is the temperature. This shows the gradient operator's role in describing heat flow from hotter to colder regions. The heat equation itself, $\frac{\partial T}{\partial t} = \alpha \nabla^2 T$, involves the Laplacian operator, which is the divergence of the gradient.

5.4 Mechanics

In mechanics, conservative force fields (like gravitational force) are gradients of potential energy functions. For example, $\mathbf{F} = -\nabla U$, where $U$ is the potential energy. This implies that conservative forces are irrotational ($\nabla \times \mathbf{F} = \mathbf{0}$).