The Central Idea

A physical field model assigns measurable quantities to positions in space and sometimes to time.

A scalar field assigns one number:

\[ f(x,y,z,t). \]

A vector field assigns a vector:

\[ \mathbf{F}(x,y,z,t). \]

Vector calculus tells us how to extract physically meaningful information from those fields.

The Field-Modeling Toolbox

The main tools of this section can be organized by the question they answer.

Gradient:

\[ \boxed{ \nabla f } \]

asks: In what direction does a scalar field increase fastest?

Divergence:

\[ \boxed{ \nabla\cdot\mathbf{F} } \]

asks: Does the vector field behave locally like a source or sink?

Curl:

\[ \boxed{ \nabla\times\mathbf{F} } \]

asks: Does the field have local rotational tendency?

Integral Questions

Line integrals answer accumulation questions along paths:

\[ \int_C \mathbf{F}\cdot d\mathbf{r}. \]

Surface integrals answer accumulation questions across surfaces:

\[ \iint_S f\,dS. \]

Flux measures vector-field crossing through an oriented surface:

\[ \iint_S \mathbf{F}\cdot\mathbf{n}\,dS. \]

Local Laws and Global Laws

Physical models are often written in two equivalent styles.

A local law uses derivatives such as divergence or curl.

A global law uses integrals over curves, surfaces, or volumes.

Green's Theorem, Stokes' Theorem, and the Divergence Theorem connect these two descriptions.

Scalar Fields in Physical Models

Common scalar fields include:

Vector Fields in Physical Models

Common vector fields include:

Potential Fields

Many physical vector fields arise from scalar potentials.

Mathematically, a gradient field has the form

\[ \mathbf{F} = \nabla f. \]

In physical applications, a minus sign is also common:

\[ \boxed{ \mathbf{F} = -\nabla U }. \]

The sign convention depends on what the scalar potential represents.

Mechanical Force Fields

A position-dependent force can be modeled by

\[ \mathbf{F}(x,y,z). \]

If a particle follows a path \(C\), the work done by the force field is

\[ \boxed{ W = \int_C \mathbf{F}\cdot d\mathbf{r} }. \]

Work as Directional Accumulation

The dot product

\[ \mathbf{F}\cdot d\mathbf{r} \]

counts only the component of force aligned with the particle's displacement.

Force perpendicular to motion contributes no instantaneous work.

Conservative Mechanical Forces

A conservative mechanical force is commonly represented as

\[ \boxed{ \mathbf{F} = -\nabla U }. \]

Here \(U\) is potential energy.

The work from \(A\) to \(B\) is

\[ \boxed{ W = U(A)-U(B) }. \]

Why Potential Energy Is Useful

If the force is conservative, the work depends only on the endpoints.

We do not need to know the complete route taken.

This converts a path-history problem into a state-comparison problem.

Example: Quadratic Potential Energy

Suppose

\[ U(x,y) = \frac{k}{2} (x^2+y^2). \]

Then

\[ \nabla U = \langle kx,ky\rangle. \]

Therefore,

\[ \boxed{ \mathbf{F} = \langle-kx,-ky\rangle }. \]

The force points toward the origin.

Fluid Velocity Fields

A fluid velocity field is

\[ \boxed{ \mathbf{v}(x,y,z,t) }. \]

At each point and time, the vector tells us the local direction and speed of fluid motion.

Streamlines

At a fixed time, a streamline is a curve whose tangent direction aligns with the velocity field.

If a curve is parameterized by

\[ \mathbf{r}(s), \]

then its tangent direction satisfies

\[ \boxed{ \mathbf{r}'(s) \parallel \mathbf{v}(\mathbf{r}(s)) }. \]

Fluid Divergence

The divergence

\[ \boxed{ \nabla\cdot\mathbf{v} } \]

measures local volumetric expansion or compression of the velocity field.

Positive divergence corresponds to local expansion.

Negative divergence corresponds to local compression.

Incompressible Flow

An idealized incompressible velocity field satisfies

\[ \boxed{ \nabla\cdot\mathbf{v}=0 }. \]

Under suitable assumptions, no net volume is created or destroyed locally by the flow.

Fluid Curl and Vorticity

The vorticity field is

\[ \boxed{ \boldsymbol{\omega} = \nabla\times\mathbf{v} }. \]

Vorticity measures local rotational character of the flow.

Circulation in a Fluid

Circulation around a closed loop is

\[ \boxed{ \Gamma = \oint_C \mathbf{v}\cdot d\mathbf{r} }. \]

Stokes' Theorem connects this to vorticity flux:

\[ \boxed{ \Gamma = \iint_S \boldsymbol{\omega}\cdot\mathbf{n}\,dS }. \]

Fluid Flux

The volumetric flow through a surface \(S\) is modeled by

\[ \boxed{ \iint_S \mathbf{v}\cdot\mathbf{n}\,dS }. \]

Positive values indicate net flow in the chosen normal direction.

Mass Flux

If \(\rho\) is mass density and \(\mathbf{v}\) is velocity, then

\[ \boxed{ \rho\mathbf{v} } \]

is the mass-flux density.

The mass flow rate through a surface is

\[ \boxed{ \iint_S \rho\mathbf{v}\cdot\mathbf{n}\,dS }. \]

The Continuity Equation

Local conservation of mass is expressed by

\[ \boxed{ \frac{\partial\rho}{\partial t} + \nabla\cdot(\rho\mathbf{v}) = 0 }. \]

This equation says that local accumulation of mass is balanced by net mass flow.

Understanding the Continuity Equation

If

\[ \nabla\cdot(\rho\mathbf{v})>0, \]

more mass is flowing outward locally than inward.

Therefore, unless mass is supplied by another mechanism, local density must decrease.

Integral Conservation of Mass

For a fixed control volume \(V\),

\[ \boxed{ \frac{d}{dt} \iiint_V \rho\,dV = - \oiint_{\partial V} \rho\mathbf{v}\cdot\mathbf{n}\,dS }. \]

The rate of mass accumulation inside equals negative net outward mass flux.

The Divergence Theorem Connects the Two Mass Laws

Apply

\[ \oiint_{\partial V} \rho\mathbf{v}\cdot\mathbf{n}\,dS = \iiint_V \nabla\cdot(\rho\mathbf{v})\,dV. \]

This converts the global conservation law into the local continuity equation under suitable smoothness assumptions.

Gravitational Potential

In a conservative gravitational model, a gravitational field

\[ \mathbf{g} \]

can be represented locally by a potential

\[ \Phi \]

through

\[ \boxed{ \mathbf{g} = -\nabla\Phi }. \]

Radial Gravitational Fields

Outside an idealized spherically symmetric source, a common radial pattern has the form

\[ \boxed{ \mathbf{g} = -\frac{C}{\rho^2}\mathbf{e}_\rho }. \]

The minus sign indicates inward radial direction.

Why Spherical Coordinates Fit Gravitation

A radial field depends mainly on distance from the origin.

In spherical coordinates,

\[ \rho \]

is exactly that distance and

\[ \mathbf{e}_\rho \]

is exactly the radial direction.

Flux of an Inverse-Square Field

For

\[ \mathbf{F} = \frac{C}{\rho^2}\mathbf{e}_\rho, \]

the outward flux through a sphere of radius \(\rho\) is

\[ \boxed{ 4\pi C }. \]

The result does not depend on the sphere radius.

Electric Potential

In electrostatics, the electric field can be expressed in terms of electric potential \(V\):

\[ \boxed{ \mathbf{E} = -\nabla V }. \]

The electric field points in the direction of greatest decrease of potential.

Equipotential Surfaces

An equipotential surface satisfies

\[ V(x,y,z)=c. \]

Since

\[ \mathbf{E} = -\nabla V, \]

the electric field is normal to regular equipotential surfaces.

Electric Work and Potential Difference

For an electrostatic field,

\[ \int_A^B \mathbf{E}\cdot d\mathbf{r} = V(A)-V(B). \]

Thus potential difference replaces a path integral when the field is conservative.

Electric Flux

Electric flux through an oriented surface is

\[ \boxed{ \Phi_E = \iint_S \mathbf{E}\cdot\mathbf{n}\,dS }. \]

For a closed surface, outward electric flux is especially important.

Gauss-Type Structure

A local electric-field law involving divergence can be converted into a global closed-surface flux law using the Divergence Theorem.

The mathematical pattern is

\[ \boxed{ \nabla\cdot\mathbf{E} \quad\longleftrightarrow\quad \oiint_{\partial V} \mathbf{E}\cdot\mathbf{n}\,dS }. \]

Magnetic Fields

A magnetic field is represented by

\[ \boxed{ \mathbf{B}(x,y,z,t) }. \]

Magnetic field models frequently involve both curl and flux.

Magnetic Flux

Magnetic flux through a surface is

\[ \boxed{ \Phi_B = \iint_S \mathbf{B}\cdot\mathbf{n}\,dS }. \]

Divergence-Free Magnetic Field Structure

A standard magnetic-field law has the mathematical form

\[ \boxed{ \nabla\cdot\mathbf{B}=0 }. \]

By the Divergence Theorem, this implies zero net magnetic flux through a suitable closed surface:

\[ \boxed{ \oiint_{\partial V} \mathbf{B}\cdot\mathbf{n}\,dS = 0 }. \]

Curl Laws in Electromagnetism

Electromagnetic field models also contain relations involving curl.

Stokes' Theorem converts a local curl law

\[ \nabla\times\mathbf{F} \]

into a circulation law around the boundary:

\[ \boxed{ \oint_C \mathbf{F}\cdot d\mathbf{r} = \iint_S (\nabla\times\mathbf{F}) \cdot\mathbf{n}\,dS }. \]

Current Density

Electric current distributed through a material can be modeled by a current density vector field

\[ \boxed{ \mathbf{J}(x,y,z,t) }. \]

The current crossing a surface is modeled by

\[ \boxed{ I = \iint_S \mathbf{J}\cdot\mathbf{n}\,dS }. \]

Charge Conservation

If \(\rho_q\) is charge density and \(\mathbf{J}\) is current density, local charge conservation has the form

\[ \boxed{ \frac{\partial\rho_q}{\partial t} + \nabla\cdot\mathbf{J} = 0 }. \]

This has the same mathematical structure as mass conservation.

Heat as a Scalar Field

Temperature is modeled by a scalar field

\[ \boxed{ T(x,y,z,t) }. \]

The gradient

\[ \nabla T \]

points toward greatest local temperature increase.

Heat-Flux Vector

In an isotropic conduction model, heat flux is commonly modeled by Fourier's law:

\[ \boxed{ \mathbf{q} = -k\nabla T }. \]

Here \(k>0\) is thermal conductivity.

Heat flows toward decreasing temperature.

Heat Flow Through a Surface

The rate of heat crossing an oriented surface is modeled by

\[ \boxed{ \iint_S \mathbf{q}\cdot\mathbf{n}\,dS }. \]

Positive sign means heat flows in the chosen normal direction.

Why the Laplacian Appears in Heat Models

Since

\[ \mathbf{q} = -k\nabla T, \]

taking divergence gives

\[ \nabla\cdot\mathbf{q} = -k\nabla^2T \]

when \(k\) is constant.

This is one reason the Laplacian naturally appears in heat equations.

Steady-State Heat

In a simple source-free steady-state conduction model with constant conductivity,

\[ \boxed{ \nabla^2T=0 }. \]

The temperature field is then harmonic.

Diffusion and Concentration Fields

Let

\[ c(x,y,z,t) \]

represent concentration.

A simple diffusion-flux model has the form

\[ \boxed{ \mathbf{J}_c = -D\nabla c }. \]

Here \(D>0\) is a diffusion coefficient.

Why Diffusion Resembles Heat Conduction

Both models take the form

\[ \boxed{ \text{flux} = -(\text{coefficient}) \times \text{gradient} }. \]

The gradient determines the direction of steepest increase, while the minus sign drives transport toward lower values.

General Conservation Law

Many physical conservation laws can be written as

\[ \boxed{ \frac{\partial u}{\partial t} + \nabla\cdot\mathbf{J} = s }. \]

Here:

Reading the General Conservation Law

The term

\[ \frac{\partial u}{\partial t} \]

measures local accumulation.

The term

\[ \nabla\cdot\mathbf{J} \]

measures local net outward transport.

The term

\[ s \]

represents local production or removal not caused by transport across the boundary.

Integral Form of a Conservation Law

Integrating over a fixed volume \(V\),

\[ \frac{d}{dt} \iiint_V u\,dV + \oiint_{\partial V} \mathbf{J}\cdot\mathbf{n}\,dS = \iiint_V s\,dV. \]

This says:

\[ \boxed{ \text{accumulation} + \text{outflow} = \text{production} }. \]

Why Conservation Laws Are So General

The same mathematical architecture appears in models of:

Coordinate Systems in Physical Models

The geometry of the physical system should guide coordinate choice.

Use Cartesian coordinates for rectangular geometry.

Use cylindrical coordinates for axial symmetry.

Use spherical coordinates for radial symmetry.

Example: Pipe Flow

A straight circular pipe naturally suggests cylindrical coordinates.

A simplified axial velocity field might have the form

\[ \boxed{ \mathbf{v} = v_z(r)\mathbf{e}_z }. \]

The velocity depends on distance from the pipe axis but points along the pipe.

Example: Rotational Flow

A purely rotating planar flow may be written

\[ \boxed{ \mathbf{v} = \Omega r\mathbf{e}_\theta }. \]

In Cartesian coordinates, this becomes

\[ \mathbf{v} = \langle -\Omega y, \Omega x, 0 \rangle. \]

Curl of Rigid Rotation

For

\[ \mathbf{v} = \langle -\Omega y, \Omega x, 0 \rangle, \]

we have

\[ \boxed{ \nabla\times\mathbf{v} = \langle0,0,2\Omega\rangle }. \]

The vorticity points along the axis of rotation.

Example: Uniform Expansion

Let

\[ \mathbf{v} = a\langle x,y,z\rangle. \]

Then

\[ \nabla\cdot\mathbf{v} = 3a. \]

If \(a>0\), the flow expands locally.

If \(a<0\), it contracts locally.

Example: Radial Source Flow

Consider

\[ \mathbf{v} = \frac{C}{\rho^2} \mathbf{e}_\rho. \]

Away from the origin,

\[ \nabla\cdot\mathbf{v}=0. \]

Yet a sphere around the origin has nonzero flux.

The origin is a singular point and must be treated separately.

Why Singular Sources Matter

A field can be divergence-free everywhere it is smooth and still have nonzero flux through a closed surface enclosing a singularity.

This does not contradict the Divergence Theorem because the standard theorem requires sufficient smoothness throughout the enclosed volume.

Boundary Conditions

Physical field models are rarely determined by differential equations alone.

Boundary conditions specify what happens on the edge of the modeled region.

Examples include:

Initial Conditions

Time-dependent field models also require an initial state.

For example,

\[ \boxed{ T(x,y,z,0)=T_0(x,y,z) } \]

specifies the initial temperature distribution.

Units Are Part of the Model

Every derivative changes units.

If \(f\) has units \(U\), then

\[ [\nabla f] = \frac{U}{\text{length}}. \]

If \(\mathbf{F}\) has units \(V\), then

\[ [\nabla\cdot\mathbf{F}] = \frac{V}{\text{length}}. \]

Dimensional consistency is a powerful error check.

Integral Units

A line integral has units

\[ [\mathbf{F}] \times \text{length}. \]

A flux integral has units

\[ [\mathbf{F}] \times \text{area}. \]

A volume integral of divergence has the same units as the corresponding closed-surface flux.

Steady Versus Time-Dependent Fields

A steady field does not explicitly depend on time.

For example,

\[ \mathbf{v}(x,y,z). \]

An unsteady field does:

\[ \mathbf{v}(x,y,z,t). \]

Time dependence introduces another layer of change beyond spatial derivatives.

Eulerian Description

A field description such as

\[ \mathbf{v}(x,y,z,t) \]

asks what is happening at each spatial location and time.

This is commonly called an Eulerian viewpoint in continuum models.

Following a Moving Particle

If a particle follows

\[ \mathbf{r}(t), \]

then a scalar field \(f(x,y,z,t)\) experienced by the particle becomes

\[ f(\mathbf{r}(t),t). \]

Differentiating gives

\[ \boxed{ \frac{d}{dt} f(\mathbf{r}(t),t) = \frac{\partial f}{\partial t} + \nabla f\cdot\mathbf{r}'(t) }. \]

Material Derivative Preview

If the particle moves with a velocity field

\[ \mathbf{v}, \]

then

\[ \mathbf{r}'(t)=\mathbf{v}. \]

Therefore,

\[ \boxed{ \frac{Df}{Dt} = \frac{\partial f}{\partial t} + \mathbf{v}\cdot\nabla f }. \]

This is the rate of change experienced while moving with the flow.

Why the Material Derivative Matters

A field can be steady at fixed locations while a moving particle still experiences change because it travels into regions with different field values.

The term

\[ \mathbf{v}\cdot\nabla f \]

captures that advective change.

Field Lines and Physical Interpretation

Field-line diagrams are useful, but they are qualitative.

They can suggest:

Quantitative conclusions still require the field formula and calculus.

Common Error: Treating a Diagram as a Proof

A vector plot can suggest that divergence or curl is positive, negative, or zero.

The actual value is determined by derivatives:

\[ \nabla\cdot\mathbf{F} \]

and

\[ \nabla\times\mathbf{F}. \]

Common Error: Ignoring the Meaning of the Field

The same mathematical vector field can represent velocity, force, heat flux, current density, or another quantity.

The mathematics may be identical while the units and interpretation differ.

Common Error: Mixing Potential Conventions

A mathematical gradient field may be written

\[ \mathbf{F}=\nabla f. \]

Many physical force or transport laws instead use

\[ \mathbf{F}=-\nabla U. \]

Always preserve the sign convention of the model.

Common Error: Confusing Flux with Circulation

Flux measures crossing through a surface:

\[ \iint_S \mathbf{F}\cdot\mathbf{n}\,dS. \]

Circulation measures alignment around a path:

\[ \oint_C \mathbf{F}\cdot d\mathbf{r}. \]

They answer different physical questions.

Common Error: Applying a Theorem Across a Singularity

Green's, Stokes', and Divergence theorems require suitable smoothness on the relevant region.

Singularities must be identified before a global integral theorem is used.

A Reliable Physical-Field Workflow

  1. Identify what the field represents.
  2. Identify whether it is scalar or vector-valued.
  3. Record the units.
  4. Identify the geometry and coordinate system.
  5. Determine the physical question: change, source, rotation, work, circulation, or flux.
  6. Choose gradient, divergence, curl, line integral, or surface integral accordingly.
  7. Check whether an integral theorem can simplify the calculation.
  8. Check domain, singularities, orientation, and boundary conditions.
  9. Interpret the sign and magnitude physically.

A Reliable Conservation-Law Workflow

  1. Identify the stored quantity \(u\).
  2. Identify its flux \(\mathbf{J}\).
  3. Identify any source term \(s\).
  4. Write \(\partial u/\partial t+\nabla\cdot\mathbf{J}=s\).
  5. Check units of all terms.
  6. Integrate over a control volume if a global balance is needed.
  7. Use the Divergence Theorem to convert between local and global forms.

A Reliable Theorem-Selection Workflow

Use the Fundamental Theorem for Line Integrals when:

\[ \mathbf{F} = \nabla f. \]

Use Green's Theorem for closed planar curves.

Use Stokes' Theorem for circulation around the boundary of an oriented surface.

Use the Divergence Theorem for outward flux through a closed surface.

Example: Heat-Flux Model

Suppose

\[ T(x,y,z)=100-x^2-y^2-z^2 \]

and conductivity is a positive constant \(k\).

Then

\[ \nabla T = \langle-2x,-2y,-2z\rangle. \]

Fourier's law gives

\[ \mathbf{q} = -k\nabla T = \boxed{ \langle 2kx, 2ky, 2kz \rangle }. \]

Heat flows outward from the hotter center toward cooler surroundings.

Example: Heat Flux Through a Sphere

On the sphere

\[ x^2+y^2+z^2=R^2, \]

the outward unit normal is

\[ \mathbf{n} = \frac{\langle x,y,z\rangle}{R}. \]

Therefore,

\[ \mathbf{q}\cdot\mathbf{n} = \frac{2k(x^2+y^2+z^2)}{R} = 2kR. \]

The total outward heat flux is

\[ 2kR(4\pi R^2) = \boxed{ 8\pi kR^3 }. \]

Example: Fluid Source and Rotation Together

Consider

\[ \mathbf{v} = \langle ax-\Omega y, \Omega x+ay, az \rangle. \]

Its divergence is

\[ \nabla\cdot\mathbf{v} = a+a+a = \boxed{3a}. \]

Its curl is

\[ \boxed{ \nabla\times\mathbf{v} = \langle0,0,2\Omega\rangle }. \]

The model combines uniform expansion or contraction with rigid rotational tendency.

Example: Conservative Force Test

Let

\[ \mathbf{F}(x,y,z) = \langle 2x, 2y, 2z \rangle. \]

The curl is

\[ \mathbf{0}. \]

A potential is

\[ f=x^2+y^2+z^2. \]

Therefore,

\[ \boxed{ \mathbf{F}=\nabla f }. \]

If interpreted as a mechanical force with potential energy \(U\), the physical convention would instead write

\[ \mathbf{F}=-\nabla U. \]

Practice 1: Temperature Gradient

Let

\[ T(x,y,z)=20+x^2+2y^2+3z^2. \]

Find \(\nabla T\).

Practice 2: Heat-Flux Direction

If

\[ \mathbf{q}=-k\nabla T \]

with \(k>0\), does heat flow toward increasing or decreasing temperature?

Practice 3: Incompressible Flow

Determine whether

\[ \mathbf{v} = \langle -y, x, 0 \rangle \]

is divergence-free.

Practice 4: Vorticity

Find the vorticity of the field in Practice 3.

Practice 5: Mechanical Work

A conservative force satisfies

\[ \mathbf{F}=-\nabla U. \]

If

\[ U(A)=12 \]

and

\[ U(B)=5, \]

find the work done by the force from \(A\) to \(B\).

Practice 6: Mass Flux

If fluid density is constant \(\rho_0\) and velocity is \(\mathbf{v}\), what vector field represents mass flux?

Practice 7: Continuity Equation

Write the local conservation-of-mass equation for density \(\rho\) and velocity \(\mathbf{v}\).

Practice 8: Radial Coordinate Choice

Which coordinate system is usually most natural for an inverse-square field centered at the origin?

Practice 9: Theorem Selection

Which theorem is most directly designed to convert outward flux through a closed surface into an integral over the enclosed volume?

Practice 10: Conservation-Law Interpretation

In

\[ \frac{\partial u}{\partial t} + \nabla\cdot\mathbf{J} = s, \]

what do \(u\), \(\mathbf{J}\), and \(s\) represent?

Practice Answers

Practice 1

\[ \boxed{ \nabla T = \langle 2x, 4y, 6z \rangle }. \]

Practice 2

Because of the negative sign,

\[ \boxed{ \text{heat flows toward decreasing temperature} }. \]

Practice 3

We have

\[ \nabla\cdot\mathbf{v} = 0+0+0 = \boxed{0}. \]

Therefore the field is divergence-free.

Practice 4

\[ \boxed{ \nabla\times\mathbf{v} = \langle0,0,2\rangle }. \]

Practice 5

\[ W = U(A)-U(B) = 12-5 = \boxed{7}. \]

Practice 6

\[ \boxed{ \rho_0\mathbf{v} }. \]

Practice 7

\[ \boxed{ \frac{\partial\rho}{\partial t} + \nabla\cdot(\rho\mathbf{v}) = 0 }. \]

Practice 8

\[ \boxed{ \text{spherical coordinates} }. \]

Practice 9

\[ \boxed{ \text{The Divergence Theorem} }. \]

Practice 10

What to Remember

Physical field models use the same vector-calculus operators repeatedly:

\[ \boxed{ \nabla f } \]

for directional change,

\[ \boxed{ \nabla\cdot\mathbf{F} } \]

for local source-sink behavior,

\[ \boxed{ \nabla\times\mathbf{F} } \]

for local rotational behavior,

\[ \boxed{ \int_C \mathbf{F}\cdot d\mathbf{r} } \]

for work or circulation, and

\[ \boxed{ \iint_S \mathbf{F}\cdot\mathbf{n}\,dS } \]

for flux.

Section Synthesis

Vector calculus can now be viewed as one connected system.

Scalar fields describe distributed scalar quantities.

Vector fields describe distributed directions and magnitudes.

Gradient converts scalar variation into a vector field.

Divergence measures local outward imbalance.

Curl measures local rotational imbalance.

Line and surface integrals accumulate field effects over geometric objects.

Flux measures transport through boundaries.

Potentials identify path-independent structure.

Green's Theorem, Stokes' Theorem, and the Divergence Theorem connect local derivatives to global boundary behavior.

Where We Go Next

This lesson completes the core sequence on Vector Calculus and Fields.

The tools developed here form a foundation for differential equations, continuum mechanics, fluid dynamics, electromagnetism, heat transfer, transport theory, mathematical physics, and advanced applied modeling.