Green's Theorem

Цели урока

• State and apply Green's theorem to convert boundary line integrals to double integrals
• Use the area formula to compute areas of planar regions from their boundary parametrizations
• Understand why zero curl implies conservativeness on simply connected domains
• See Green's theorem as the 2D case of the generalized Stokes theorem

Предварительные знания

• Line integrals
• Partial derivatives and curl
• Parametric curves

• Line Integrals

How can the area of an arbitrary region be computed from only its boundary - with no integration over the interior whatsoever?

• PostGIS and GIS systems: polygon areas computed from GPS boundary coordinates via the shoelace formula - Green's theorem in O(n)
• COMSOL, Abaqus, FEniCS: FEM converts differential equations to weak form via Green's theorem before any matrix is assembled
• Computational fluid dynamics: finite volume method uses discrete Green's theorem to convert volume sources into face fluxes
• Computer graphics: signed area of a triangle mesh face is computed via the shoelace formula, enabling fast inside-outside testing

Green, the self-taught miller

George Green (1793-1841) received roughly one year of formal schooling as a child. He was a miller in Nottingham who taught himself mathematics from a library subscription. His 1828 essay 'An Essay on the Application of Mathematical Analysis to the Theories of Electricity and Magnetism' contained the theorem now bearing his name, printed in 51 copies. Green died in 1841. Four years later, Lord Kelvin found a copy and sent it to Jacobi, Liouville, Crelle, and Sturm. All were astonished - Green had anticipated work that would take Europe another decade. One essay, 51 copies, foundation of a 2.4 billion USD industry.

Statement and Geometric Meaning

1828. George Green - a self-taught miller's son from Nottingham - publishes an essay in a print run of 51 copies. Lord Kelvin finds a copy 17 years after Green's death and mails it to every major mathematician in Europe. Today Ansys sells FEM simulation licenses generating 2.4 billion USD in annual revenue. Every solver in that portfolio uses Green's theorem internally to exchange boundary integrals for area integrals.

Three equivalent area formulas from Green's theorem: S = oint x dy = -oint y dx = (1/2) oint (x dy - y dx). All give the same result. The shoelace algorithm uses the last form for polygon vertices.

Curl, Conservativeness, and Green's Theorem

GIS systems like PostGIS compute areas of geographic polygons from GPS coordinates. The underlying algorithm is the shoelace formula - a direct discrete version of Green's theorem. One billion polygon area computations per day, all reducing to the observation that dQ/dx - dP/dy = 1 when Q=x, P=0.

Green's theorem requires P and Q to have continuous partial derivatives on the closed region D. If the field has a singularity inside D (like F = (-y,x)/(x^2+y^2) at the origin), the theorem does not apply directly - the singular point must be excluded.

Shoelace algorithm for polygon area: S = (1/2) |sum_{i=0}^{n-1} (x_i * y_{i+1} - x_{i+1} * y_i)|. This is the discrete version of (1/2) oint (x dy - y dx). Runs in O(n) - directly derived from Green's theorem.

Proof Sketch and Connection to Stokes

The proof of Green's theorem for a rectangle is three lines. For a general region: decompose into elementary pieces, apply the rectangle case, note that shared boundaries cancel. This exact same decomposition-and-cancellation idea drives FEM: each element contributes boundary terms, interior ones cancel, and only the global boundary remains.

Cauchy's theorem in complex analysis (oint_C f(z) dz = 0 for analytic f) is a direct consequence of Green's theorem. The Cauchy-Riemann equations du/dx = dv/dy and du/dy = -dv/dx are exactly the zero-curl conditions for the two real-valued fields formed from the real and imaginary parts of f.

Numerical Applications and FEM

Finite Element Method (FEM) is built entirely on Green's theorem. The domain is split into small triangles. On each triangle, the differential equation is multiplied by a test function and integrated. Green's theorem converts volume derivatives into boundary terms. Shared interior edges cancel. Only the global boundary remains. This is the weak formulation - the foundation of every FEM solver from Abaqus to FEniCS.

COMSOL Multiphysics, the industry standard for FEM simulation of heat transfer, fluid flow, and electromagnetics, converts every PDE into a weak form via Green's theorem before discretizing. One research license: around 10,000 USD. Every differential equation transformed by Green's theorem before a single matrix entry is computed.

Green's first identity iint |grad u|^2 dA = oint u * du/dn ds - iint u * Delta u dA is the basis for the Dirichlet principle: the harmonic function (Delta u = 0) minimizes the 'Dirichlet energy' iint |grad u|^2 dA among all functions with the same boundary values.

Connections to other topics

Green's theorem sits at the intersection of real analysis, vector calculus, complex analysis, and numerical methods.

• Complex analysis — Related topic
• Finite Element Method — Related topic
• Computational geometry — Related topic
• Fluid mechanics — Related topic

Итоги

• Green's theorem: oint_C P dx + Q dy = iint_D (dQ/dx - dP/dy) dA - boundary integral = double integral of curl
• Area formulas: S = oint x dy = -oint y dx = (1/2) oint (x dy - y dx) - all from Green with curl = 1
• Conservativeness: oint = 0 iff curl = 0 on simply connected D - Green provides the rigorous proof
• Cauchy's theorem in complex analysis is Green's theorem + Cauchy-Riemann equations
• FEM weak formulation is Green's theorem applied to convert second-order PDE to first-order variational form

Вопросы для размышления

• Green's theorem says boundary data determines interior behavior. What does this suggest about the nature of harmonic functions (which satisfy Delta u = 0)?
• The shoelace formula computes polygon area from vertices alone. Could one generalize this to compute volume of a 3D polyhedron from face data - and what theorem would that be?
• Cauchy's theorem says analytic functions integrate to zero around any loop. What does this say about the topology of the complex plane and the singularities of analytic functions?

Связанные уроки

• calc-24-line-integrals — Line integrals along curves are the left side of Green's theorem
• calc-26-divergence-theorem — Divergence theorem is the 3D analogue
• calc-27-diff-forms — Differential forms state Green's theorem as d omega = boundary integral
• la-06-transformations