Matly

Càlcul vectorial: integrals de línia i de superfície, Green, Stokes i Gauss

Vols practicar? Aquest tema té 110 exercicis interactius de 3 tipus, de deu en deu, amb correcció immediata i explicació pas a pas.

El càlcul vectorial és on el càlcul de diverses variables es converteix en física. Treball, flux, circulació i caudal deixen de ser paraules i passen a ser integrals concretes, i els tres grans teoremes (Green, Stokes i Gauss) diuen tots la mateixa cosa amb vestits diferents: el que passa a la frontera queda determinat pel que passa a dins.

Camps i operadors

Un camp escalar assigna un número a cada punt; un camp vectorial, un vector. Els tres operadors bàsics, en coordenades cartesianes, són

∇f=(fx,fy,fz),div⁡F=Px+Qy+Rz,\nabla f = (f_x, f_y, f_z), \qquad \operatorname{div}F = P_x + Q_y + R_z, rot⁡F=(Ry−Qz, Pz−Rx, Qx−Py).\operatorname{rot}F = (R_y - Q_z,\ P_z - R_x,\ Q_x - P_y).

El gradient apunta en la direcció de màxim creixement i dona la derivada direccional com a producte escalar amb el vector unitari. La divergència mesura la densitat local de flux sortint: on és positiva, hi ha una font. El rotacional mesura la circulació local, és a dir com de retorçat està el camp.

Dues identitats que surten de la igualtat de derivades creuades per a camps de classe C2C^{2}, i que es fan servir contínuament: el rotacional d'un gradient és zero i la divergència d'un rotacional és zero.

Integrals de línia: dues, i no es confonen

Sobre una corba parametritzada hi ha dues integrals diferents:

∫Cf ds=∫abf(r(t)) ∥r′(t)∥ dt,∫CF⋅dr=∫abF(r(t))⋅r′(t) dt.\int_C f\,ds = \int_a^b f(r(t))\,\lVert r'(t)\rVert\,dt, \qquad \int_C F\cdot dr = \int_a^b F(r(t))\cdot r'(t)\,dt.

La primera és escalar (longitud, massa d'un filferro) i no canvia si inverteixes l'orientació, perquè la norma és sempre positiva. La segona és treball o circulació i canvia de signe en invertir l'orientació. Confondre el diferencial d'arc amb el diferencial vectorial és l'error més freqüent del tema.

Camps conservatius

Teorema fonamental de les integrals de línia

Si el camp és el gradient d'un potencial de classe C1C^{1}, aleshores

∫CF⋅dr=ϕ(r(b))−ϕ(r(a)).\int_C F\cdot dr = \phi(r(b)) - \phi(r(a)).

La prova és la regla de la cadena aplicada a la composició del potencial amb la parametrització, més el teorema fonamental del càlcul d'una variable.

En un obert connex per camins, són equivalents: que la integral sigui independent del camí, que existeixi potencial, i que la circulació sobre tota corba tancada sigui nul·la. La condició pràctica de rotacional nul és suficient només si el domini és simplement connex, és a dir sense forats.

Un camp amb rotacional zero que no és conservatiu

Al pla menys l'origen, considera

F(x,y)=(−yx2+y2, xx2+y2).F(x,y) = \left(\frac{-y}{x^{2}+y^{2}},\ \frac{x}{x^{2}+y^{2}}\right).

Un càlcul directe dona Qx−Py=0Q_x - P_y = 0 a tot el domini. I tanmateix, sobre la circumferència unitat parametritzada de la manera habitual, el producte escalar val exactament 11 a cada instant, de manera que la circulació és 2π2\pi, no zero.

On és el truc? El domini té un forat a l'origen i no és simplement connex. La circumferència no es pot contraure a un punt sense travessar-lo. Per tant no existeix cap potencial global, encara que en local sí que n'hi hagi (l'angle polar, que no es pot definir de manera contínua en tota la volta).

La moral: comprovar que el rotacional és nul no demostra que un camp sigui conservatiu si no has mirat abans la forma del domini.

Per trobar un potencial, integres una component respecte de la seva variable i hi afegeixes una funció desconeguda de les altres; després derives i imposes les components que queden. Per a F=(2xy, x2+2y)F = (2xy,\ x^{2}+2y): integrant la primera respecte de xx surt ϕ=x2y+h(y)\phi = x^{2}y + h(y), i imposant la segona queda h′(y)=2yh'(y) = 2y, és a dir h=y2+Ch = y^{2} + C.

Superfícies, àrea i flux

Amb una parametrització regular d'una superfície,

∬Sf dS=∬Df ∥ru×rv∥ du dv,∬SF⋅n dS=∬DF⋅(ru×rv) du dv.\iint_S f\,dS = \iint_D f\,\lVert r_u\times r_v\rVert\,du\,dv, \qquad \iint_S F\cdot n\,dS = \iint_D F\cdot (r_u\times r_v)\,du\,dv.

No multipliquis dues vegades l'element d'àrea

A la fórmula del flux no hi torna a aparèixer la norma: el producte vectorial ja incorpora l'element d'àrea orientada. Posar-hi la norma a més és un error que multiplica el resultat per un factor variable.

I sempre cal comprovar que el producte vectorial té l'orientació que demana l'enunciat; si apunta al revés, el flux canvia de signe.

Si la superfície és el gràfic d'una funció amb orientació cap amunt, l'element vectorial és (−gx,−gy,1) dx dy(-g_x, -g_y, 1)\,dx\,dy, que és una de les fórmules que més s'utilitzen.

Els tres teoremes

Green relaciona una circulació plana amb una integral doble:

∮∂DP dx+Q dy=∬D(Qx−Py) dx dy,\oint_{\partial D} P\,dx + Q\,dy = \iint_D (Q_x - P_y)\,dx\,dy,

amb l'orientació positiva, que deixa la regió a l'esquerra: frontera exterior en sentit antihorari i forats en sentit horari. Es dedueix primer en rectangles aplicant el teorema fonamental a cada variable, i després es descompon el domini; les vores interiors es cancel·len entre elles, i per això l'orientació és essencial.

Stokes és la versió en l'espai: la circulació sobre la frontera d'una superfície orientable és el flux del rotacional a través d'ella, amb l'orientació induïda per la regla de la mà dreta. Una conseqüència pràctica molt útil: pots canviar la superfície per una altra amb la mateixa frontera orientada, si el camp és regular a la regió que hi ha entremig.

Gauss o de la divergència: el flux exterior a través d'una superfície tancada és la integral triple de la divergència. La condició de tancada no és negociable; si a la teva superfície li falta una tapa, l'has d'afegir, aplicar el teorema, i després restar el flux de la tapa.

El mateix flux, calculat de dues maneres

Calcula el flux exterior del camp de posició F=(x,y,z)F = (x,y,z) a través de l'esfera de radi aa.

Amb Gauss. La divergència val 33 a tot arreu, de manera que

∬∂VF⋅n dS=∭V3 dV=3⋅4πa33=4πa3.\iint_{\partial V}F\cdot n\,dS = \iiint_V 3\,dV = 3\cdot\frac{4\pi a^{3}}{3} = 4\pi a^{3}.

Directament. A l'esfera, el camp de posició és radial i la normal exterior també, amb la qual cosa el producte escalar val exactament aa a tot punt. Multiplicant per l'àrea de l'esfera:

a⋅4πa2=4πa3.a\cdot 4\pi a^{2} = 4\pi a^{3}.

Coincideixen. Aquest tipus de doble càlcul és la millor manera de guanyar confiança amb els teoremes: quan un problema es pot fer de dues maneres, fes-lo de les dues fins que t'hi acostumis.

Stokes estalviant una parametrització

Calcula la circulació de F=(−y/2, x/2, 0)F = (-y/2,\ x/2,\ 0) sobre la circumferència unitat del pla z=0z=0, vista en sentit antihorari des de dalt.

Amb Stokes. El rotacional val (0,0,1)(0,0,1). Prenent com a superfície el disc pla, la normal és (0,0,1)(0,0,1) i el flux és simplement l'àrea del disc:

∬S(rot⁡F)⋅n dS=π.\iint_S (\operatorname{rot}F)\cdot n\,dS = \pi.

Directament. Parametritzant la circumferència, el producte escalar surt constant i igual a 1/21/2, i integrant entre 00 i 2π2\pi dona també π\pi.

Aquest camp concret és el que dona l'àrea encerclada per qualsevol corba tancada simple, i és la base del planímetre.

Exercici resolt: flux pel cub

Calcula el flux exterior de F=(x2,y2,z2)F = (x^{2}, y^{2}, z^{2}) pel cub unitat.

La superfície és tancada, o sigui que Gauss s'hi aplica directament. La divergència val

div⁡F=2x+2y+2z.\operatorname{div}F = 2x + 2y + 2z.

Integrant sobre el cub, cada variable aporta el mateix per simetria:

∭2x dV=2⋅12=1,\iiint 2x\,dV = 2\cdot\frac12 = 1,

i igual per a les altres dues. El flux total és 33.

Compara-ho amb el càlcul directe: caldria parametritzar sis cares, vigilar sis orientacions i fer sis integrals dobles. Gauss ho redueix a una integral triple d'un polinomi. Quan la superfície és tancada, gairebé sempre val la pena.

Exercici resolt: treball d'un camp conservatiu

Calcula el treball de F=(2xy, x2)F = (2xy,\ x^{2}) per anar de l'origen al punt (1,2)(1,2).

Primer, és conservatiu? Qx=2xQ_x = 2x i Py=2xP_y = 2x: coincideixen, i el domini és tot el pla, que és simplement connex. Sí que ho és.

Busquem el potencial. Integrant la primera component respecte de xx: ϕ=x2y+h(y)\phi = x^{2}y + h(y). Derivant respecte de yy i igualant a la segona component: x2+h′(y)=x2x^{2} + h'(y) = x^{2}, o sigui hh constant. Prenem ϕ=x2y\phi = x^{2}y.

El treball és la diferència de potencial entre els extrems:

ϕ(1,2)−ϕ(0,0)=2−0=2,\phi(1,2) - \phi(0,0) = 2 - 0 = 2,

i no depèn del camí. Si l'enunciat t'hagués donat una corba complicada, no hauries hagut ni de mirar-la: comprovar que el camp és conservatiu val més que qualsevol parametrització.

Ara prova-ho tu. 110 exercicis d'aquest tema, de 3 tipus diferents, per resoldre en línia, de deu en deu: escrius la resposta, te la corregeix a l'instant i, si falles, t'explica com es fa pas a pas. Quan acabis una tanda, en tens deu més, amb dades noves.

Fer els exercicisSimulacre d'examen