The divergence theorem states ∬_S F · n dS = ∭_W ∇ · F dV, where S is the closed surface bounding W with outward normal n. Total flux out of region W equals the integral of divergence throughout W. Useful for computing flux through closed surfaces without parametrization.
From your study of the divergence operator, you know that ∇ · F at a point measures the local rate of expansion of the vector field — positive divergence means the field is spreading outward from that point (a source), negative divergence means it's converging (a sink). From surface integrals, you know that ∬_S F · n dS measures the flux — the net amount of the field passing outward through a surface S. The divergence theorem links these two ideas: the total flux through the outer boundary of a region equals the total source strength accumulated throughout the interior.
The physical intuition is easiest to grasp with fluid flow. Imagine a region W of space filled with a fluid whose velocity field is F. The flux integral ∬_S F · n dS counts the net volume of fluid per unit time flowing out through the boundary surface S. If the divergence ∇ · F is positive throughout W, the fluid is expanding — sources inside W are pumping fluid outward, and that fluid exits through S. The divergence theorem says these two perspectives — tallying what exits the boundary versus tallying what's produced inside — must agree. The equation ∬_S F · n dS = ∭_W ∇ · F dV is an exact bookkeeping identity.
The practical power of the theorem is computational flexibility. Surface integrals over complicated closed surfaces can be nightmarish to set up directly — parametrizing a sphere or a cylinder with caps requires careful bookkeeping of orientation and limits. The divergence theorem lets you replace the surface integral with a volume integral, which is often much easier. For example, if F = ⟨x, y, z⟩, then ∇ · F = 3, and the flux through any closed surface bounding a region W is simply 3 · Vol(W) — no surface parametrization needed at all.
The same swap works in the other direction: a volume integral over a complicated region W can sometimes be converted to a surface integral over the simpler boundary ∂W. Choosing the direction of the conversion depends on which integral is easier to evaluate. This flexibility — volume ↔ surface — is the defining feature of all the theorems in vector calculus (Green's, Stokes', and divergence): they all trade one type of integral for another across one dimension of boundary. Understanding the divergence theorem as a higher-dimensional analog of the Fundamental Theorem of Calculus (where ∫_a^b f' dx = f(b) − f(a) trades interior information for boundary information) is the deepest way to understand why these theorems hold.