Knowing whether a vector field is expanding, contracting, or in perfect balance at any given point is a core requirement in fluid dynamics, electromagnetics, and heat transfer — and it's not always obvious from the field equations alone. Use this Divergence Vector Calculator to calculate the divergence of a vector field using component expressions and a choice of Cartesian, cylindrical, or spherical coordinate systems. Getting this right matters in CFD mass conservation checks, Gauss's law applications in antenna design, and thermal diffusion modeling in electronics cooling. This page includes the divergence formula for all 3 coordinate systems, a worked nozzle flow example, full theory, and an FAQ covering practical engineering use cases.
What is divergence of a vector field?
Divergence is a single number — a scalar — that tells you whether a vector field is spreading outward (positive), converging inward (negative), or perfectly balanced (zero) at a specific point in space. It measures how much the field acts as a source or a sink at that location.
Simple Explanation
Think of a vector field like water flowing through a room — arrows show which way the water moves at every point. Divergence tells you whether water is gushing out from a location like a sprinkler (positive), draining into it like a drain (negative), or flowing through without any net gain or loss (zero). No math required to grasp it: it's just a measure of how "sprinkler-like" or "drain-like" any point in the field is.
📐 Browse all 1000+ Interactive Calculators
Vector Field Diagram
Divergence Vector Calculator
How to Use This Calculator
- Select your coordinate system from the Calculation Mode dropdown — Cartesian, cylindrical, spherical, or one of the specialized modes.
- Enter the mathematical expressions for each vector field component in the input fields shown (e.g., Fx, Fy, Fz for Cartesian).
- If using "Divergence at Specific Point," enter the coordinates of the evaluation point as comma-separated values (x, y, z).
- Click Calculate to see your result.
📹 Video Walkthrough — Divergence Vector Interactive Calculator
Divergence Vector Interactive Visualizer
Visualize how vector fields expand, contract, or remain balanced at different points in space. Adjust field components to see real-time divergence calculations and field line behavior.
DIVERGENCE
+3.5
FIELD TYPE
SOURCE
EXPANSION
350%
FIRGELLI Automations — Interactive Engineering Calculators
Divergence Equations
Use the formula below to calculate divergence in your chosen coordinate system.
Cartesian Coordinates
div F = ∇ · F = ∂Fx/∂x + ∂Fy/∂y + ∂Fz/∂z
Where:
- F = vector field with components (Fx, Fy, Fz)
- ∂Fx/∂x = partial derivative of x-component with respect to x
- ∂Fy/∂y = partial derivative of y-component with respect to y
- ∂Fz/∂z = partial derivative of z-component with respect to z
- div F = scalar field representing source/sink density (units depend on F)
Cylindrical Coordinates
div F = (1/r)·∂(rFr)/∂r + (1/r)·∂Fθ/∂θ + ∂Fz/∂z
Where:
- r = radial distance from z-axis (m)
- θ = azimuthal angle (radians)
- z = axial coordinate (m)
- Fr, Fθ, Fz = vector field components in cylindrical basis
Spherical Coordinates
div F = (1/r²)·∂(r²Fr)/∂r + (1/(r·sin θ))·∂(sin θ·Fθ)/∂θ + (1/(r·sin θ))·∂Fφ/∂φ
Where:
- r = radial distance from origin (m)
- θ = polar angle from z-axis (radians, 0 to π)
- φ = azimuthal angle in xy-plane (radians, 0 to 2π)
- Fr, Fθ, Fφ = vector field components in spherical basis
Divergence Theorem (Gauss's Theorem)
∭V (∇ · F) dV = ∬S F · n dS
Where:
- V = enclosed volume (m³)
- S = closed surface bounding volume V (m²)
- n = outward-pointing unit normal vector
- Left side = volume integral of divergence (total source/sink)
- Right side = surface integral of flux (net outward flow)
Simple Example
Take a simple Cartesian vector field: F = (2x, 3y, z). The divergence is ∂(2x)/∂x + ∂(3y)/∂y + ∂(z)/∂z = 2 + 3 + 1 = 6. At every point in space, this field has a divergence of 6 — a uniform source density throughout the volume.
Theory & Engineering Applications
Divergence quantifies the rate at which a vector field expands or contracts at each point in space, providing a scalar measure of field behavior crucial for conservation laws across physics and engineering. Unlike gradient or curl, divergence reduces a vector field to a scalar field, revealing sources (positive divergence) where field lines originate and sinks (negative divergence) where they terminate. This fundamental operation underlies Maxwell's equations in electromagnetism, the continuity equation in fluid dynamics, and heat diffusion analysis in thermodynamics.
Mathematical Foundation and Coordinate System Dependencies
The Cartesian divergence formula ∂Fx/∂x + ∂Fy/∂y + ∂Fz/∂z represents the most intuitive formulation, directly summing the rate of change of each component along its respective axis. However, this simplicity masks critical geometric considerations that emerge in curvilinear coordinates. In cylindrical coordinates, the term (1/r) multiplying ∂(rFr)/∂r accounts for the geometric fact that radial field lines spread as they move outward — the same flux passing through a cylinder of radius r must traverse a larger surface area at radius 2r, reducing the effective field strength. This geometric dilution appears explicitly in the divergence formula and is non-negotiable for accurate calculations in pipe flow, rotating machinery, and cylindrical electromagnetic systems.
Spherical coordinates introduce even more dramatic geometric factors with both (1/r²) and (1/sin θ) terms. The r² factor in ∂(r²Fr)/∂r reflects the rapid expansion of spherical surfaces with radius — area scales as 4πr², so a constant radial flux density requires the field magnitude to decay as 1/r². The sin θ term corrects for the varying "width" of longitude lines at different latitudes, vanishing at the poles where longitude lines converge. These corrections are essential for planetary science (atmospheric circulation), astrophysics (stellar winds), and antenna theory (radiation patterns). A common error in spherical calculations is neglecting these geometric factors, producing divergence values that violate conservation of mass or energy by orders of magnitude.
Physical Interpretation and Conservation Laws
The divergence at a point measures the net flux per unit volume emerging from an infinitesimal region surrounding that point. Positive divergence indicates expansion — field lines spreading apart as if originating from a source. In fluid dynamics, this represents mass generation (injection) or thermal expansion. Negative divergence indicates compression or convergence toward a sink, representing mass removal (suction) or thermal contraction. Zero divergence characterizes solenoidal or incompressible fields where every bit of flux entering a region immediately exits, maintaining constant density or charge conservation.
The divergence theorem connects this local property to global behavior by equating the volume integral of divergence to the surface integral of flux. This is not merely a mathematical curiosity but the foundation of all conservation laws in continuum mechanics. Consider a control volume in a fluid flow: if the total divergence inside is positive, more fluid exits through the surface than enters, requiring either mass generation or decreasing density. For incompressible flow (∇ · v = 0 everywhere), the theorem guarantees that inflow exactly balances outflow — the cornerstone of hydraulic system design.
Engineering Applications Across Disciplines
In computational fluid dynamics (CFD), divergence appears in the continuity equation ∂ρ/∂t + ∇ · (ρv) = 0, which governs mass conservation. Numerical solvers compute divergence at every grid cell, adjusting pressure and velocity fields to drive divergence toward zero (for incompressible flow) or balance it with density changes (for compressible flow). High divergence values in a CFD solution signal numerical instability, unresolved shocks, or physically impossible mass creation, triggering solver refinement or remeshing. Engineers designing aircraft wings analyze divergence of the velocity field to predict shock wave formation and drag — regions of high positive divergence indicate compression shocks that dramatically increase drag above Mach 1.
Electromagnetic systems rely on Maxwell's divergence equations: ∇ · E = ρ/ε₀ (Gauss's law) and ∇ · B = 0 (no magnetic monopoles). The first equation directly relates electric field divergence to charge density, allowing calculation of electric fields from known charge distributions via the divergence theorem. Antenna engineers use this to compute radiation patterns — regions with high electric field divergence correspond to charge accumulation on the antenna structure, determining impedance and efficiency. The zero divergence of magnetic fields has profound consequences: magnetic field lines must form closed loops, never beginning or ending, fundamentally distinguishing magnetism from electrostatics and enabling transformer operation and magnetic confinement fusion.
Heat transfer analysis uses divergence in Fourier's law combined with energy conservation: ρcp∂T/∂t = -∇ · q + Q, where q is heat flux and Q is internal generation. The divergence of heat flux determines whether a region heats or cools. Thermal engineers designing electronics cooling systems calculate ∇ · q to identify hot spots — high positive divergence indicates insufficient heat removal, requiring enhanced cooling. In geothermal energy extraction, negative divergence zones in the subsurface thermal field reveal aquifers drawing heat from surrounding rock, optimal for well placement.
Worked Example: Incompressible Flow Through a Contracting Nozzle
A water nozzle transitions from a circular cross-section of diameter D₁ = 0.050 m at the inlet to D₂ = 0.025 m at the exit over a length L = 0.150 m. The flow velocity at the inlet is V₁ = 2.4 m/s, axially directed. Assuming steady, axisymmetric, incompressible flow, determine:
(a) The exit velocity V₂ using continuity
(b) The divergence of the velocity field at a point halfway through the nozzle
(c) Interpretation of the divergence in terms of flow physics
Solution:
(a) Exit velocity from continuity equation:
For incompressible flow, the continuity equation requires constant volumetric flow rate Q:
Q = A₁V₁ = A₂V₂
Calculate inlet area:
A₁ = π(D₁/2)² = π(0.050/2)² = π(0.025)² = 1.963 × 10⁻³ m²
Calculate exit area:
A₂ = π(D₂/2)² = π(0.025/2)² = π(0.0125)² = 4.909 × 10⁻⁴ m²
Solve for exit velocity:
V₂ = V₁(A₁/A₂) = 2.4 × (1.963 × 10⁻³)/(4.909 × 10⁻⁴) = 2.4 × 4.0 = 9.6 m/s
(b) Divergence at midpoint of nozzle:
At the midpoint (z = L/2 = 0.075 m from inlet), the diameter is:
D(z) = D₁ - (D₁ - D₂)(z/L) = 0.050 - (0.050 - 0.025)(0.075/0.150) = 0.050 - 0.0125 = 0.0375 m
Midpoint area:
Amid = π(0.0375/2)² = 1.104 × 10⁻³ m²
Velocity at midpoint from continuity:
Vmid = V₁(A₁/Amid) = 2.4 × (1.963 × 10⁻³)/(1.104 × 10⁻³) = 4.27 m/s
For axisymmetric flow in cylindrical coordinates with only axial velocity component vz(z) and radial symmetry (no θ or r dependence for the axial component):
∇ · v = (1/r)·∂(rvr)/∂r + ∂vz/∂z
The radial velocity vr exists because the flow must follow the contracting walls. From continuity in integral form and geometry, the axial velocity gradient is:
∂vz/∂z ≈ (V₂ - V₁)/L = (9.6 - 2.4)/0.150 = 48.0 s⁻¹
However, for truly incompressible flow, ∇ · v = 0 everywhere. The positive axial acceleration must be exactly balanced by negative radial velocity contribution (fluid moving radially inward as it accelerates axially). To verify:
The radial component satisfies (1/r)·∂(rvr)/∂r = -∂vz/∂z
Therefore, the divergence is:
∇ · v = 0 (exactly, for incompressible flow)
(c) Interpretation:
The zero divergence confirms incompressibility — no volume expansion or contraction occurs despite the 4× velocity increase from inlet to exit. The acceleration happens purely through geometric constriction, converting pressure head to velocity head as described by Bernoulli's equation. The radial inward velocity component vr (fluid moving toward the centerline as it flows downstream) exactly compensates for the axial stretching represented by ∂vz/∂z.
This balance is fundamental to nozzle design: any deviation from zero divergence would indicate compressibility effects (important above Mach 0.3) or numerical errors in CFD simulations. In practical terms, the zero divergence guarantees that the mass flow rate ρQ is conserved, with the density ρ = 998 kg/m³ for water at 20°C remaining constant throughout the nozzle. The pressure drop from inlet to exit can be calculated from Bernoulli's equation, yielding ΔP = (1/2)ρ(V₂² - V₁²) = (1/2)(998)(9.6² - 2.4²) = 43,100 Pa = 43.1 kPa, representing the conversion of pressure energy to kinetic energy with no energy dissipation to divergence-related compression work.
For further exploration of vector calculus in engineering contexts, visit our comprehensive engineering calculator library, which includes tools for gradient, curl, and line integral calculations across multiple coordinate systems.
Practical Applications
Scenario: HVAC System Design for Commercial Building
Marcus, a mechanical engineer designing the ventilation system for a 12-story office building, needs to ensure proper air circulation in the main lobby — a large atrium space with complex geometry. Using CFD software, he simulates the airflow from multiple ceiling diffusers and computes the divergence field throughout the 800 m³ lobby volume. The simulation reveals a region near the center with divergence of +0.15 s⁻¹, indicating excessive air expansion and explaining temperature complaints from occupants in that zone. By adjusting diffuser angles and flow rates to drive divergence closer to zero everywhere (target below 0.02 s⁻¹), Marcus achieves uniform temperature distribution and meets ASHRAE ventilation standards. This divergence-based optimization reduces energy consumption by 18% compared to the initial design by eliminating unnecessary recirculation zones and ensuring each air change effectively ventilates the space rather than creating stagnant pockets.
Scenario: Antenna Radiation Pattern Analysis
Dr. Chen, an RF engineer developing a 5.8 GHz patch antenna for a wireless communication system, must verify that the antenna impedance matches the 50Ω transmission line to maximize power transfer. She uses electromagnetic simulation software to compute the electric field distribution near the antenna surface and calculates ∇ · E to map the surface charge density via Gauss's law (ρ = ε₀∇ · E). The divergence calculation reveals unexpected charge accumulation (∇ · E = 4.2 × 10⁷ V/m²) at the feed point due to a geometry discontinuity in her initial design, creating a reactive impedance component that causes 35% of the input power to reflect. By smoothing the feed transition to reduce the peak divergence to 1.1 × 10⁷ V/m², Dr. Chen achieves a voltage standing wave ratio (VSWR) of 1.3:1, meaning 96% of power transmits effectively. This divergence-guided optimization is essential because modern communication systems demand strict impedance matching — even small reflections accumulate in multi-antenna arrays and beam-forming systems, degrading signal quality.
Scenario: Groundwater Contamination Modeling
Jennifer, an environmental engineer investigating a potential fuel leak near a municipal water supply, uses hydrogeological modeling to predict contaminant migration through the subsurface aquifer. She models the Darcy velocity field (groundwater flow) and computes its divergence across a 500m × 500m × 30m control volume centered on the suspected leak site. Most of the aquifer shows near-zero divergence (±0.0002 day⁻¹), consistent with incompressible saturated flow, but a localized zone 15 meters below the surface exhibits divergence of -0.012 day⁻¹ — indicating a convergence zone where groundwater is being withdrawn. Cross-referencing property records reveals an unpermitted irrigation well operating in that zone, drawing 2,400 liters per hour. The negative divergence calculation allows Jennifer to quantify the capture zone radius (approximately 38 meters) where contaminants would be drawn toward the well rather than flowing toward the municipal supply intake 220 meters downstream. This analysis demonstrates that the well actually provides hydraulic containment, preventing contaminant migration — a counterintuitive finding that saves the municipality $1.2 million in remediation costs by redirecting cleanup efforts to the actual at-risk pathway identified through divergence mapping.
Frequently Asked Questions
What does positive vs. negative divergence physically mean? +
Why do cylindrical and spherical coordinate formulas look more complicated than Cartesian? +
How is divergence used to verify computational fluid dynamics (CFD) simulations? +
What is the relationship between divergence and the divergence theorem in practical calculations? +
Can a vector field have zero curl but non-zero divergence, or vice versa? +
How do numerical methods approximate divergence in discrete computational grids? +
Free Engineering Calculators
Explore our complete library of free engineering and physics calculators.
Browse All Calculators →🔗 Explore More Free Engineering Calculators
About the Author
Robbie Dickson — Chief Engineer & Founder, FIRGELLI Automations
Robbie Dickson brings over two decades of engineering expertise to FIRGELLI Automations. With a distinguished career at Rolls-Royce, BMW, and Ford, he has deep expertise in mechanical systems, actuator technology, and precision engineering.
Need to implement these calculations?
Explore the precision-engineered motion control solutions used by top engineers.
