Una primera simplificación es que en el teorema de Stokes, todo ocurre en la superficie. Así que, en cierto sentido, se trata de un teorema 2D, no 3D. Si nos limitamos a la parametrización en la superficie, el teorema de Stokes es precisamente el teorema de Green.
Entonces: ¿cómo descubrirías el teorema de Green?
Empecemos con un teorema más intuitivo: el teorema de la divergencia. El teorema de la divergencia dice que si tenemos un flujo constante de agua con un montón de fuentes y sumideros, entonces para cualquier volumen V que elijamos, la cantidad neta de agua creada en el interior V es igual a la cantidad neta de agua que fluye a través de la superficie de V .
Espero que la idea que subyace a este teorema sea suficientemente intuitiva, aunque traducirla al lenguaje formal es otra cuestión.
Fíjate ahora en que seguramente existe una versión 2D del teorema de la divergencia, en la que todo el fluir del agua ocurre en el plano. Ese teorema dice que la fuente neta en un zona es igual al flujo neto que atraviesa la curva límite de la zona.
Resulta que éste es precisamente el teorema de Green. Esto es muy poco claro de la forma en que se suele presentar, donde el teorema de Green dice algo así como "la cantidad de swirly en un área es igual a la integral de línea a lo largo de la frontera", signifique lo que signifique.
Afortunadamente, resulta mucho más claro cuando rotamos cada vector del campo vectorial por 90 grados. Rotando el campo vectorial por 90 ¡grados convierte cada poco de remolino en un poco de divergencia! Y convierte la integral de línea a lo largo de la curva en el flujo a través de la curva, un concepto mucho más intuitivo.
Si nos fijamos en la declaración, queda más claro. El enunciado habitual es que si tienes un campo vectorial (L,M) entonces tenemos
∮Ldx+Mdy=∬
Obsérvese que la rotación por 90 grados da (L,M) \mapsto (-M,L) por lo que obtenemos
\oint Ldy - Mdx = \iint \left(\frac{\partial L}{\partial x} + \frac{\partial M}{\partial y}\right)\,dxdy
Vemos la divergencia 2D \frac{\partial L}{\partial x} + \frac{\partial M}{\partial y} a la derecha, y a la izquierda vemos Ldy - Mdx = (L,M) \cdot (dy,-dx) = (L,M) \cdot \hat{n}, es decir, el flujo a través de la frontera.