Conservative Fields, Green's Theorem, the Divergence Theorem & Generalized Stokes Mechanics

Vector calculus is the mathematical language of fluid mechanics, electromagnetism, gravitational fields, computer graphics physics, and general relativity. By connecting differential operators (divergence and curl) to boundary integrals over curves, surfaces, and volumes, vector calculus reveals how localized micro-scale rates of change accumulate into global macro-scale flux and circulation.

1. Vector Fields & Line Integrals#

In multi-dimensional Euclidean space, physical quantities like gravity, electrostatic force, and fluid flow vary continuously across spatial coordinates.

Definition of a Vector Field#

A Vector Field in R^n is an assignment of an n-dimensional vector F(x) to every point x in a subset of space. In 3-dimensional Euclidean space R^3:

F(x, y, z) = [ P(x, y, z), Q(x, y, z), R(x, y, z) ]^T

where P, Q, and R are scalar field component functions representing forces or velocities along the x, y, and z axes.

The Line Integral of a Vector Field#

Let C be a smooth, oriented curve in R^3 parameterized by the position vector r(t) for a <= t <= b:

r(t) = [ x(t), y(t), z(t) ]^T, with velocity vector r'(t) = [ x'(t), y'(t), z'(t) ]^T.

The line integral of vector field F along curve C evaluates the cumulative tangential work performed by the field on a particle traversing the trajectory:

integral_C F · dr = integral_a^b F(r(t)) · r'(t) dt

= integral_a^b [ P(x(t), y(t), z(t)) · x'(t) + Q(x(t), y(t), z(t)) · y'(t) + R(x(t), y(t), z(t)) · z'(t) ] dt.

Conservative Vector Fields & Potential Functions#

A vector field F is defined to be conservative (or irrotational / gradient field) if there exists a scalar potential function φ: R^n → R such that F is the exact gradient of φ:

F = ∇φ.

(In physical mechanics, force is often defined as F = -∇U, where U is potential energy).

The Fundamental Theorem of Line Integrals (Path Independence)#

Theorem: Let C be any piecewise-smooth curve starting at point A = r(a) and terminating at point B = r(b). If F is a conservative vector field (F = ∇φ), then the line integral depends strictly on the endpoints:

integral_C F · dr = integral_C ∇φ · dr = φ(B) - φ(A) = φ(r(b)) - φ(r(a)).

Analytical Proof:

By the multivariable chain rule:

d/dt [ φ(r(t)) ] = ∇φ(r(t)) · r'(t) = F(r(t)) · r'(t).

Substituting into the line integral definition:

integral_C F · dr = integral_a^b [ d/dt (φ(r(t))) ] dt.

By the single-variable Fundamental Theorem of Calculus:

= φ(r(b)) - φ(r(a)) = φ(B) - φ(A).

Corollaries of Conservative Fields:

1. Complete Path Independence: The work done moving between any two points A and B is identical along any path, whether straight or circuitous.
2. Zero Closed-Loop Circulation: The line integral around ANY closed loop γ (where start A = end B) is identically zero:

oint_γ F · dr = 0.

2. Differential Operators of Vector Analysis: Divergence and Curl#

Vector calculus is governed by two spatial differential operators that measure local expansion and local rotation.

1. Divergence: Volumetric Source Flux Density#

The divergence of a vector field F = [ P, Q, R ]^T, denoted div F or ∇ · F, is a scalar field defined by the formal dot product of the del operator with F:

div F = ∇ · F = ∂P/∂x + ∂Q/∂y + ∂R/∂z.

Physical & Geometric Meaning:

Divergence measures the net rate of outward fluid flux per unit volume at an infinitesimal point (expansion versus contraction):

• ∇ · F > 0: Net Source. Fluid is being generated or expanding outward from that point.
• ∇ · F < 0: Net Sink. Fluid is being absorbed, compressed, or draining inward.
• ∇ · F = 0: Incompressible (Solenoidal) Field. Whatever fluid volume enters an infinitesimal box must exit without accumulation.

2. Curl: Microscopic Circulation & Vorticity Density#

The curl of a vector field F, denoted curl F or ∇ × F, is a vector field defined by the formal cross product of the del operator with F:

curl F = ∇ × F = det([

[ i, j, k ],

[ ∂/∂x, ∂/∂y, ∂/∂z ],

[ P, Q, R ]

])

Evaluating the 3×3 determinant:

curl F = [ (∂R/∂y - ∂Q/∂z), (∂P/∂z - ∂R/∂x), (∂Q/∂x - ∂P/∂y) ]^T.

Physical & Geometric Meaning:

Curl measures the rotational tendency or microscopic vorticity of the field at a point:

• Direction: The direction of curl F defines the axis of maximum local rotation, adhering to the right-hand rule.
• Magnitude: ||curl F|| equals twice the local angular rotational speed. If a tiny paddlewheel were placed in the fluid at point (x, y, z), it would spin if and only if curl F ≠ 0!
• Irrotational Field: A vector field is irrotational if curl F = 0 everywhere.

+-----------------------------------------------------------------------------+
|                     DIVERGENCE VERSUS CURL GEOMETRY                         |
+-----------------------------------------------------------------------------+
|                                                                             |
|         DIVERGENCE (Flux Density)               CURL (Vorticity Density)    |
|                                                                             |
|                    ^                                     ^                  |
|                    |                                    / \                 |
|               <----+---->                           <--*   *-->             |
|                    |                                    \ /                 |
|                    v                                     v                  |
|          Net Outward Expansion                 Local Spinning Rotation      |
|           div(F) = nabla · F                     curl(F) = nabla × F        |
|                                                                             |
+-----------------------------------------------------------------------------+

Fundamental Vector Calculus Null Identities#

Two universal identities govern smooth vector fields:

1. The Divergence of Any Curl is Identically Zero:

∇ · (∇ × F) = 0.

Vorticity produces closed circulating loops; it possesses no net volumetric sources or sinks.

1. The Curl of Any Gradient is Identically Zero:

∇ × (∇φ) = 0.

Conservative gradient fields are purely translational; they possess zero microscopic rotation.

Theorem (Component Test for Conservative Fields):

Let F be a C^1 vector field on a simply connected 3D domain (a domain with no holes).

F is conservative (F = ∇φ) if and only if curl F = 0:

∂R/∂y = ∂Q/∂z, ∂P/∂z = ∂R/∂x, ∂Q/∂x = ∂P/∂y.

3. The Triad of Classical Vector Theorems#

The classical integral theorems demonstrate that integrating a spatial derivative over an interior region equals the evaluation of the field across the enclosing boundary.

1. Green's Theorem (2D Planar Regions)#

Let D be a bounded region in the 2D plane R^2, whose boundary C = ∂D is a piecewise-smooth, simple closed curve oriented counterclockwise (positive orientation: traversing C keeps region D on the left).

Let F = [ P(x, y), Q(x, y) ]^T be a C^1 vector field on D.

Theorem (Green's Theorem):

oint_C (P dx + Q dy) = iint_D (∂Q/∂x - ∂P/∂y) dA.

The microscopic 2D curl (∂Q/∂x - ∂P/∂y) integrated across the entire area of D equals the macroscopic circulation around the enclosing perimeter loop C!

Planar Area Computation via Green's Boundary Integrals:

Setting (∂Q/∂x - ∂P/∂y) = 1 yields boundary formulas for the area of D:

• Choice 1: P = 0, Q = x => Area(D) = oint_C x dy.
• Choice 2: P = -y, Q = 0 => Area(D) = -oint_C y dx.
• Choice 3 (Symmetric Form): Area(D) = (1/2) · oint_C (x dy - y dx).

In computational geometry, this symmetric formula forms the mathematical basis of Gauss's Shoelace Algorithm for computing the exact area of arbitrary irregular polygons from ordered vertex coordinates.

2. Gauss' Divergence Theorem (3D Solid Volumes)#

Let V be a bounded 3-dimensional solid region in R^3, whose boundary S = ∂V is a closed, piecewise-smooth surface oriented with an outward-pointing unit normal vector n.

Let F be a C^1 vector field defined on V.

Theorem (Gauss' Divergence Theorem):

iint_{∂V} (F · n) dS = iiint_V (∇ · F) dV.

The total net outward flux of vector field F passing through the enclosing boundary surface S equals the volume integral of internal source divergence within V!

Physical Invariant (Gauss's Law of Electrostatics):

In electromagnetism, electric flux density D satisfies ∇ · D = ρ_free (charge density).

Integrating over volume V:

iint_{∂V} (D · n) dS = iiint_V (∇ · D) dV = iiint_V ρ_free dV = Q_enclosed.

The net electric field bursting through any enclosing closed surface depends exclusively on the total net electric charge trapped inside, regardless of how the charge is distributed!

3. Classical Stokes' Theorem (Oriented 3D Surfaces)#

Let S be an oriented, piecewise-smooth 3-dimensional surface in R^3, bounded by a closed boundary curve C = ∂S.

The orientation of C is positive relative to the surface unit normal n according to the Right-Hand Rule (curling the fingers of the right hand along C points the thumb in the direction of n).

Let F be a C^1 vector field.

Theorem (Classical Stokes' Theorem):

oint_{∂S} F · dr = iint_S (∇ × F) · n dS.

The total circulation of F around the 1D boundary curve ∂S equals the surface flux of the curl (∇ × F) piercing through the 2D interior surface S!

Surface Independence Property:

Because the boundary curve C = ∂S is fixed, Stokes' Theorem holds for ANY open surface S capped by loop C. If S1 and S2 share the same boundary curve C, then:

iint_S1 (∇ × F) · n1 dS1 = iint_S2 (∇ × F) · n2 dS2 = oint_C F · dr.

The flux of curl through a balloon depends only on the perimeter loop at its mouth!

4. The Generalized Stokes' Theorem in Exterior Calculus & Differential Forms#

In advanced mathematical physics, differential geometry, and general relativity, the Fundamental Theorem of Calculus, Green's Theorem, Gauss' Divergence Theorem, and Classical Stokes' Theorem are revealed to be identical manifestations of a single master equation.

The Master Formulation#

Theorem (Generalized Stokes' Theorem):

Let Ω be an oriented, compact k-dimensional smooth manifold with boundary ∂Ω.

Let ω be a smooth differential (k - 1)-form defined on Ω.

Then:

integral_{∂Ω} ω = integral_Ω dω

where dω is the exterior derivative of differential form ω.

+-----------------------------------------------------------------------------+
|                 THE GRAND UNIFICATION OF INTEGRAL CALCULUS                  |
+-----------------------------------------------------------------------------+
|                                                                             |
|               /                 /                                           |
|              |      omega  =   |     d(omega)                               |
|             /                 /                                             |
|            Boundary d_Omega   Interior Omega                                |
|                                                                             |
|  "Integrating a differential form over the boundary of a region             |
|   is mathematically identical to integrating its exterior derivative        |
|   over the entire interior of the region."                                  |
+-----------------------------------------------------------------------------+

Complete Taxonomy of the Classical Theorems#

The wedge product ∧ enforces antisymmetry (dx ∧ dy = -dy ∧ dx), automatically encoding geometric orientation.

The nilpotency of the exterior derivative, d(dω) = d^2 = 0, identically unifies div(curl F) = 0 and curl(grad φ) = 0!

5. Computational Physics, Fluid Dynamics & Computer Graphics Applications#

Vector calculus is the foundation of physical modeling in modern computing.

1. Incompressible Fluid Simulation: The Navier-Stokes Equations#

In computer graphics (Pixar, Disney, video game physics), simulating smoke, water, and fire uses the incompressible Navier-Stokes equations:

∂u/∂t + (u · ∇) u = -(1 / ρ) ∇p + ν ∇^2 u + g

subject to the continuity constraint: ∇ · u = 0.

• (u · ∇) u: Non-linear advection (velocity transporting itself).
• -(1 / ρ) ∇p: Pressure gradient force pushing fluid from high to low pressure.
• ν ∇^2 u: Viscous diffusion (Laplacian operator ∇^2 = ∇ · ∇).
• ∇ · u = 0: Incompressibility constraint (divergence-free velocity field).

Jos Stam's Stable Fluids Algorithm (Helmholtz-Hodge Decomposition):

By Helmholtz's Theorem, any smooth vector field w can be uniquely decomposed into an irrotational gradient field and a solenoidal (divergence-free) field:

w = u + ∇p, where ∇ · u = 0.

Taking the divergence of both sides:

∇ · w = ∇ · u + ∇ · ∇p = 0 + ∇^2 p.

This yields the Poisson Equation for Pressure:

∇^2 p = ∇ · w.

Solving this linear sparse system on a GPU via Conjugate Gradients and subtracting u = w - ∇p projects the velocity field onto a divergence-free manifold, guaranteeing volume-conserving fluid simulations!

2. Mesh Geometry: Volume and Center of Mass via the Divergence Theorem#

In 3D CAD modeling, 3D printing (slicers), and game engines, calculating the exact volume and center of mass of a complex 3D triangle mesh (e.g., a character model with 500,000 faces) is required.

Voxelizing the interior of the mesh is computationally slow and introduces discretization error.

The Divergence Theorem Solution:

Choose the simple vector field F(x, y, z) = [ x, 0, 0 ]^T.

Its divergence is: ∇ · F = ∂x/∂x = 1.

Applying the Divergence Theorem:

iiint_V (∇ · F) dV = iiint_V 1 dV = Volume(V).

Therefore:

Volume(V) = iint_{∂V} (F · n) dS = sum_{triangles T_i} iint_{T_i} (x · n_x) dS.

The 3D volume integral collapses into a sum over the 2D surface triangles! A 500,000-polygon mesh can have its exact volume computed in milliseconds by summing dot products over surface faces.

3. Electrodynamics: Maxwell's Equations#

James Clerk Maxwell unified electricity and magnetism into four vector calculus equations:

1. Gauss's Law: ∇ · E = ρ / ε0 (Electric charges are divergence sources of E).
2. Gauss's Law for Magnetism: ∇ · B = 0 (No magnetic monopoles exist; B is solenoidal).
3. Faraday's Law of Induction: ∇ × E = -∂B/∂t (A time-varying magnetic field induces electric curl).
4. Ampère-Maxwell Law: ∇ × B = μ0 J + μ0 ε0 (∂E/∂t) (Currents and time-varying electric fields induce magnetic curl).

These four vector relationships predict the propagation of electromagnetic waves (light, Wi-Fi, radio) through vacuum at speed c = 1 / sqrt(μ0 ε0).

6. Common Architectural Pitfalls & Misconceptions#

• Pitfall 1: Applying Path Independence in Domains with Topological Holes.

Consider the vortex field F(x, y) = [ -y / (x^2 + y^2), x / (x^2 + y^2) ]^T.

Direct calculation shows:

∂Q/∂x - ∂P/∂y = [ (y^2 - x^2) / (x^2 + y^2)^2 ] - [ (y^2 - x^2) / (x^2 + y^2)^2 ] = 0.

The curl is identically zero everywhere in R^2 except at the origin (where it is undefined).

Engineers frequently assume F is conservative. However, computing the line integral around the unit circle x = cos(t), y = sin(t) yields:

oint F · dr = integral_0^{2π} [ (-sin(t))(-sin(t)) + (cos(t))(cos(t)) ] dt = integral_0^{2π} 1 dt = 2π ≠ 0!

Because the domain R^2 \ {(0, 0)} is NOT simply connected (it has a puncture hole at the origin), zero curl does NOT guarantee path independence.

• Pitfall 2: Neglecting Boundary Loop Orientation (Right-Hand Rule).

In Green's Theorem and Stokes' Theorem, the boundary curve must be positively oriented. If the loop is traversed clockwise instead of counterclockwise, the line integral flips sign (producing -Area instead of +Area).

• Pitfall 3: Applying the Divergence Theorem to an Open Surface.

Gauss' Divergence Theorem strictly requires a CLOSED surface that completely envelops a 3D solid volume (like an intact balloon). You cannot apply the Divergence Theorem to an open bowl or parachute without adding a "cap" surface and accounting for the flux through the cap.

• Pitfall 4: Neglecting Outward Unit Normal Vectors.

In surface flux integrals iint (F · n) dS, the normal vector n must point strictly OUTWARD from the enclosed volume. An inward-pointing normal flips the sign of the entire integral.

7. Summary & Checkpoint Exercises#

Solidify your comprehension of line integrals, conservative fields, Green's Theorem, Gauss' Divergence Theorem, and Stokes' Theorem with these four checkpoint problems. Full step-by-step mathematical solutions are provided.

Problem 1 (Conservative Field Verification, Potential Reconstruction & Work):

Consider the 3D vector field:

F(x, y, z) = [ 2 x y + z^3, x^2 + 3 y^2, 3 x z^2 + 1 ]^T.

1. Prove that F is a conservative vector field by computing curl F = ∇ × F.
2. Find the scalar potential function φ(x, y, z) such that F = ∇φ.
3. Compute the work done by field F along a smooth spiral path C parameterized by r(t) = [ cos(π t), sin(π t), t ]^T from t = 0 to t = 2.

Solution:

Part 1: Prove curl F = 0.

F = [ P, Q, R ]^T with P = 2 x y + z^3, Q = x^2 + 3 y^2, R = 3 x z^2 + 1.

Compute components of curl F:

• (curl F)_x = ∂R/∂y - ∂Q/∂z = d/dy [ 3 x z^2 + 1 ] - d/dz [ x^2 + 3 y^2 ] = 0 - 0 = 0.
• (curl F)_y = ∂P/∂z - ∂R/∂x = d/dz [ 2 x y + z^3 ] - d/dx [ 3 x z^2 + 1 ] = 3 z^2 - 3 z^2 = 0.
• (curl F)_z = ∂Q/∂x - ∂P/∂y = d/dx [ x^2 + 3 y^2 ] - d/dy [ 2 x y + z^3 ] = 2 x - 2 x = 0.

curl F = [ 0, 0, 0 ]^T.

Since the domain is all of R^3 (which is simply connected) and curl F = 0, field F is strictly conservative.

Part 2: Reconstruct Potential Function φ(x, y, z).

We require ∇φ = [ ∂φ/∂x, ∂φ/∂y, ∂φ/∂z ]^T = [ P, Q, R ]^T.

1. Integrate P with respect to x:

∂φ/∂x = 2 x y + z^3 => φ(x, y, z) = integral (2 x y + z^3) dx = x^2 y + x z^3 + g(y, z).

1. Differentiate with respect to y and match Q:

∂φ/∂y = x^2 + ∂g/∂y = x^2 + 3 y^2.

∂g/∂y = 3 y^2 => g(y, z) = y^3 + h(z).

φ(x, y, z) = x^2 y + x z^3 + y^3 + h(z).

1. Differentiate with respect to z and match R:

∂φ/∂z = 3 x z^2 + h'(z) = 3 x z^2 + 1.

h'(z) = 1 => h(z) = z + C.

Scalar potential function:

φ(x, y, z) = x^2 y + x z^3 + y^3 + z + C.

Part 3: Evaluate Line Integral.

By the Fundamental Theorem of Line Integrals, the path integral depends strictly on the endpoints:

• Start point A = r(0):

r(0) = [ cos(0), sin(0), 0 ]^T = [ 1, 0, 0 ]^T.

• End point B = r(2):

r(2) = [ cos(2π), sin(2π), 2 ]^T = [ 1, 0, 2 ]^T.

Evaluate potential φ at endpoints (setting C = 0):

• φ(A) = φ(1, 0, 0) = (1^2)(0) + (1)(0^3) + (0^3) + 0 = 0.
• φ(B) = φ(1, 0, 2) = (1^2)(0) + (1)(2^3) + (0^3) + 2 = 0 + 8 + 0 + 2 = 10.

Total Work:

integral_C F · dr = φ(B) - φ(A) = 10 - 0 = 10.

Problem 2 (Applying Green's Theorem to Compute Work and Area):

Consider the 2D vector field F(x, y) = [ y^2, 3 x y ]^T.

Compute the line integral oint_C F · dr around the boundary of the triangular region D with vertices (0, 0), (2, 0), and (2, 4), oriented counterclockwise.

Solution:

Step 1: Apply Green's Theorem.

By Green's Theorem:

oint_C (P dx + Q dy) = iint_D (∂Q/∂x - ∂P/∂y) dA

where P = y^2 and Q = 3 x y.

Compute partial derivatives:

• ∂Q/∂x = d/dx [ 3 x y ] = 3 y.
• ∂P/∂y = d/dy [ y^2 ] = 2 y.

Integrand:

∂Q/∂x - ∂P/∂y = 3 y - 2 y = y.

Step 2: Set up the Double Integral over Triangular Domain D.

The region D is bounded by:

• Bottom: y = 0
• Right: x = 2
• Hypotenuse: connects (0, 0) to (2, 4), slope m = (4 - 0)/(2 - 0) = 2 => y = 2 x, or x = y/2.

Express D as Type I region:

x ranges from 0 to 2.

For a given x, y ranges from y = 0 up to y = 2 x.

Step 3: Evaluate the Double Integral.

iint_D y dA = integral_0^2 [ integral_0^{2x} y dy ] dx.

Inner integral with respect to y:

integral_0^{2x} y dy = [ y^2 / 2 ]_0^{2x} = (2x)^2 / 2 - 0 = (4 x^2) / 2 = 2 x^2.

Outer integral with respect to x:

integral_0^2 2 x^2 dx = [ (2/3) x^3 ]_0^2 = (2/3)(2^3) - 0 = (2/3)(8) = 16/3.

Conclusion:

oint_C (y^2 dx + 3 x y dy) = 16/3.

(Calculating this via 3 separate line integrals along each triangle edge yields the identical 16/3).

Problem 3 (Gauss' Divergence Theorem Surface Flux Calculation):

Let S be the closed surface of the cylinder defined by x^2 + y^2 <= 4 and 0 <= z <= 3 (including the circular top cap, circular bottom cap, and cylindrical lateral wall), oriented with outward normal n.

Let vector field F(x, y, z) = [ x^3, y^3, z^2 ]^T.

Compute the total outward flux:

Phi = iint_S (F · n) dS.

Solution:

Step 1: Apply Gauss' Divergence Theorem.

Because S is a closed surface enclosing the cylindrical solid volume V, the Divergence Theorem states:

Phi = iint_S (F · n) dS = iiint_V (∇ · F) dV.

Step 2: Compute Divergence ∇ · F.

∇ · F = ∂/∂x [ x^3 ] + ∂/∂y [ y^3 ] + ∂/∂z [ z^2 ]

= 3 x^2 + 3 y^2 + 2 z = 3(x^2 + y^2) + 2 z.

Step 3: Convert to Cylindrical Coordinates.

In cylindrical coordinates:

x = r cos(θ), y = r sin(θ), z = z.

x^2 + y^2 = r^2.

Volume element: dV = r dr dθ dz.

Integration Bounds for Solid Cylinder V:

• r ranges from 0 to 2 (radius = sqrt(4) = 2).
• θ ranges from 0 to 2π.
• z ranges from 0 to 3.

Step 4: Evaluate the Triple Integral.

iiint_V [ 3 r^2 + 2 z ] r dr dθ dz = integral_0^{2π} dθ · integral_0^3 dz · integral_0^2 [ 3 r^3 + 2 z r ] dr.

1. Angular integral:

integral_0^{2π} dθ = 2π.

1. Radial integral with respect to r:

integral_0^2 [ 3 r^3 + 2 z r ] dr = [ (3/4) r^4 + z r^2 ]_0^2

= (3/4)(16) + z(4) - 0 = 12 + 4 z.

1. Height integral with respect to z:

integral_0^3 [ 12 + 4 z ] dz = [ 12 z + 2 z^2 ]_0^3

= 12(3) + 2(3^2) - 0 = 36 + 18 = 54.

1. Multiply by angular factor 2π:

Phi = 2π · 54 = 108π.

Total outward flux: Phi = 108π ≈ 339.292.

Applying the Divergence Theorem avoided calculating three separate surface integrals (top disc, bottom disc, side wall).

Problem 4 (Classical Stokes' Theorem Surface Curl Flux Verification):

Let C be the circular boundary curve in the plane z = 2 defined by x^2 + y^2 = 9, traversed counterclockwise when viewed from above.

Let vector field F(x, y, z) = [ -y, x, z^2 ]^T.

1. Evaluate the line integral oint_C F · dr directly using curve parameterization.
2. Evaluate the line integral using Classical Stokes' Theorem by choosing the flat circular disc S capped by curve C. Verify that both methods yield identical results.

Solution:

Part 1: Direct Line Integral Evaluation.

Curve C lies in plane z = 2, with radius R = 3.

Parameterization:

r(t) = [ 3 cos(t), 3 sin(t), 2 ]^T for 0 <= t <= 2π.

Velocity vector:

r'(t) = [ -3 sin(t), 3 cos(t), 0 ]^T.

Evaluate F along r(t):

F(r(t)) = [ -y(t), x(t), z(t)^2 ]^T = [ -3 sin(t), 3 cos(t), 4 ]^T.

Compute dot product F · r'(t):

F · r'(t) = (-3 sin(t))(-3 sin(t)) + (3 cos(t))(3 cos(t)) + (4)(0)

= 9 sin^2(t) + 9 cos^2(t) + 0 = 9(sin^2(t) + cos^2(t)) = 9.

Integrate:

oint_C F · dr = integral_0^{2π} 9 dt = 9(2π - 0) = 18π.

Part 2: Evaluation via Classical Stokes' Theorem.

Stokes' Theorem:

oint_C F · dr = iint_S (∇ × F) · n dS.

1. Compute curl F:

curl F = det([

[ i, j, k ],

[ ∂/∂x, ∂/∂y, ∂/∂z ],

[ -y, x, z^2 ]

])

- x-component: ∂(z^2)/∂y - ∂(x)/∂z = 0 - 0 = 0.

- y-component: ∂(-y)/∂z - ∂(z^2)/∂x = 0 - 0 = 0.

- z-component: ∂(x)/∂x - ∂(-y)/∂y = 1 - (-1) = 2.

curl F = [ 0, 0, 2 ]^T.

1. Identify Surface S and Normal Vector n:

Choose the flat circular disc S: x^2 + y^2 <= 9 in plane z = 2.

Because C is oriented counterclockwise viewed from above, the right-hand rule dictates that the surface normal points straight upward:

n = [ 0, 0, 1 ]^T.

1. Compute dot product (∇ × F) · n:

(∇ × F) · n = [ 0, 0, 2 ] · [ 0, 0, 1 ]^T = 2.

1. Evaluate Surface Integral:

iint_S 2 dS = 2 · iint_S 1 dS = 2 · Area(Disc S).

Radius of disc is R = 3:

Area(S) = π · R^2 = π · (3^2) = 9π.

Surface integral value:

iint_S (∇ × F) · n dS = 2 · (9π) = 18π.

Conclusion:

Both direct path parameterization and Stokes' surface curl integration evaluate to exactly 18π. Stokes' theorem is rigorously verified.