27.04.2013 Views

Cap´ıtol 5. Representacions de camps per funcions holomorfes.

Cap´ıtol 5. Representacions de camps per funcions holomorfes.

Cap´ıtol 5. Representacions de camps per funcions holomorfes.

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Capítol <strong>5.</strong> <strong>Representacions</strong> <strong>de</strong> <strong>camps</strong> <strong>per</strong> <strong>funcions</strong> <strong>holomorfes</strong>.<br />

J. Solà-Morales<br />

1.- Mecànica <strong>de</strong> fluids.<br />

El moviment o flux d’un fluid que ocupa una certa regió <strong>de</strong> l’espai queda <strong>de</strong>terminat si coneixem<br />

les tres components <strong>de</strong>l seu camp <strong>de</strong> velocitats V (x, y, z, t) = (u(x, y, z, t), v(x, y, z, t), w(x, y, z, t)),<br />

que dona en el temps t la velocitat V <strong>de</strong> la partícula que en aquell mateix instant ocupa la posició<br />

(x, y, z). Una versió més simple, que pot ser raonable en certs casos, és suposar que es tracta d’un<br />

flux pla o pla-parallel, és a dir que w ≡ 0 i que u i v no <strong>de</strong>penen <strong>de</strong> la coor<strong>de</strong>nada z. Una altra<br />

simplificació és suposar règim estacionari, és a dir que u i v no <strong>de</strong>penen <strong>de</strong> t. Amb aquestes dues<br />

suposicions, doncs, u ≡ u(x, y) i v ≡ v(x, y).<br />

Una altra hipòtesi que fem és que es tracta d’un fluid incompressible, que es veu reflexada en que<br />

el camp <strong>de</strong> velocitats té divergència zero:<br />

ux + vy = 0 (1)<br />

Aquesta equació pot integrar-se: en dimensió dos és equivalent a l’existència d’una funció ψ(x, y),<br />

anomenada funció <strong>de</strong> corrent tal que<br />

ψx = −v, ψy = u.<br />

Aquesta funció pot ser no unívoca, en dominis no simplement connexos. És fàcil <strong>de</strong> veure directament<br />

que ψ és una integral primera <strong>de</strong>l camp <strong>de</strong> velocitats V , <strong>per</strong>ò el significat físic <strong>de</strong> ψ ve donat pel<br />

següent<br />

Teorema (<strong>de</strong>l cabal): Sigui ψ(x, y) la funció <strong>de</strong> corrent, i Γ una corba que uneixi ψ = α amb<br />

ψ = β. El cabal que travessa Γ <strong>per</strong> unitat <strong>de</strong> temps i en sentit positiu val β − α.<br />

Dem.: (vegeu Figura 1) Suposem Γ = Γ(s) = (Γx(s), Γy(s)) parametritzada <strong>per</strong> la longitud <strong>de</strong><br />

l’arc. Tenim que n(s) = (Γ ′ y(s), −Γ ′ x(s)) ⊤ . I el flux net o cabal val<br />

L<br />

0<br />

(u, v) · nds =<br />

L<br />

0<br />

<br />

ψyΓ ′ <br />

y(s) + (−ψx)(−Γx(s)) ds = ψ(Γ(L)) − ψ(Γ(0)) = β − α. <br />

Encara una altra suposició serà que el fluid és no viscós. Matemàticament, això vol dir que es<br />

satisfà el sistema d’equacions d’Euler, que en el nostre cas (pla i estacionari) és:<br />

<br />

uux + vuy = −px/ρ,<br />

(2)<br />

uvx + vvy = −py/ρ,<br />

on p = p(x, y) és una nova incògnita, la pressió, i ρ és la <strong>de</strong>nsitat, que suposem constant.<br />

No discutirem amb <strong>de</strong>tall el significat <strong>de</strong> la suposició <strong>de</strong>l caràcter no viscós. Po<strong>de</strong>n veure’s <strong>de</strong>talls a<br />

[1]. Però direm que aquesta suposició fonamenta l’interès especial <strong>de</strong>ls fluxos irrotacionals (∇ × V ≡<br />

0), és a dir d’aquells pels que<br />

uy − vx = 0,<br />

o, en termes <strong>de</strong> ψ: ∆ψ = 0. Per aquests fluxos, també anomenats potencials, existeix, al menys<br />

en simplement connexos, una funció φ(x, y) anomenada potencial <strong>de</strong> velocitats tal que u = φx i<br />

v = φy, o sigui V = ∇φ. Notem que també a causa <strong>de</strong> la incompressibilitat ∆φ = 0.<br />

1


Figura 1: Corbes ψ = β (dalt), ψ = α (baix) i corba Γ (transversal). Vectors tangent i normal<br />

(unitaris) a Γ.<br />

Els fluxos irrotacionals compleixen les equacions d’Euler amb la pressió donada <strong>per</strong> la fórmula <strong>de</strong><br />

Bernoulli:<br />

p = p0 − ρ<br />

2 (u2 + v 2 ), (3)<br />

on p0 és una constant arbitrària.<br />

També direm que el caràcter no viscós és el fonament <strong>de</strong>l fet que la força exercida pel flux sobre<br />

un obstacle S vingui donada <strong>per</strong> <br />

−pndℓ (4)<br />

∂S<br />

En aquest cas pla en realitat estem parlant <strong>de</strong> força <strong>per</strong> unitat <strong>de</strong> longitud <strong>de</strong> l’obstacle (cil·lindric,<br />

en la direcció <strong>de</strong> z).<br />

Teorema:<br />

i) En les hipòtesis anteriors (flux pla, incompressible, no viscós i irrotacional) u − iv és una funció<br />

holomorfa (la velocitat conjugada) i u−iv = d/dz [Ω(z)], on Ω(z) = φ+iψ (que s’anomena potencial<br />

complex).<br />

ii) Si Ω(z) és una funció holomorfa, aleshores <strong>de</strong>finint (u, v) <strong>per</strong> u − iv = Ω ′ (z) es té que (u, v) és un<br />

flux que satisfà aquestes hipòtesis.<br />

Dem: i) són les equacions <strong>de</strong> Cauchy-Riemann i les <strong>de</strong>finicions <strong>de</strong> φ i ψ. ii) és també evi<strong>de</strong>nt <strong>per</strong>ò<br />

cal comprovar a més a més que, com ja hem dit, (u, v) compleix les equacions d’Euler, amb la pressió<br />

p donada <strong>per</strong> (3). <br />

En mecànica <strong>de</strong> fluids una quantitat important que serveix <strong>per</strong> a mesurar en cert sentit la rotació<br />

<strong>de</strong>l moviment és la circulació <strong>de</strong>l camp <strong>de</strong> velocitats V al llarg d’una corba tancada γ:<br />

circ (γ) =<br />

S<br />

0<br />

V · γ ′ (s)ds.<br />

Pel conegut teorema <strong>de</strong> Stokes, si V = (u, v) ⊤ , aleshores<br />

<br />

circ (∂W +) = (uy − vx) dV<br />

W<br />

i <strong>per</strong> tant en els fluxos irrotacionals aquesta integral valdrà zero. Tot i això tindrà molt d’interès<br />

el calcular la circulació al llarg <strong>de</strong> corbes tanca<strong>de</strong>s <strong>de</strong> fluxos plans incompressibles i irrotacionals<br />

2


que tinguin alguna singularitat aïllada a l’interior <strong>de</strong>l recite d’integració W . A aquests efectes convé<br />

relacionar la circulació amb una integral <strong>de</strong> línia complexa, que és la següent:<br />

<br />

circ (γ) = Re (u − iv)dz. (5)<br />

Potencials complexos <strong>de</strong> fluxos senzills.<br />

1) Flux homogeni. Es diu flux homogeni a aquell en el que la velocitat és la mateixa a tots els punts<br />

<strong>de</strong>l pla. Si el vector velocitat és α + iβ aleshores el potencial complex és Ω(z) = (α − iβ)z. Tenim<br />

ψ = −βx + αy.<br />

2) Vòrtex potencial. És el flux en que tota partícula roda al voltant d’un punt z0 amb velocitat<br />

angular que <strong>de</strong>pèn només <strong>de</strong>l radi <strong>de</strong> gir. Pot <strong>de</strong>mostrar-se que necessàriament és <strong>de</strong> la forma Ω(z) =<br />

ik log(z − z0), k ∈ R amb el que ψ = k log |z − z0|. Usant (5) veiem que la circulació d’un vòrtex<br />

potencial al llarg d’una corba γ que encercli una vegada al punt z = z0 és −2πk. Aquesta quantitat<br />

s’anomena intensitat <strong>de</strong>l vòrtex, i s’acostuma a escriure Γ. Per això el potencial complex s’escriu<br />

Si Γ > 0 el vòrtex gira en sentit positiu.<br />

Ω(z) = Γ<br />

log(z − z0).<br />

2πi<br />

3) Font o xuclador. És el flux que es produeix quan surt (o entra) <strong>de</strong>l punt z = z0 <strong>de</strong> forma <strong>per</strong>manent<br />

una quantitat <strong>de</strong> fluid Q (unitats <strong>de</strong> volum <strong>per</strong> unitat <strong>de</strong> temps, o sigui cabal). Naturalment, si Q > 0<br />

és una font, i si Q < 0 és un xuclador. Q rep el nom <strong>de</strong> potència <strong>de</strong> la font. És<br />

Ω(z) = Q<br />

log(z − z0).<br />

2π<br />

Naturalment també pot consi<strong>de</strong>rar-se la font-vòrtex<br />

Ω(z) =<br />

Γ + iQ<br />

2πi<br />

γ<br />

log(z − z0).<br />

4) Su<strong>per</strong>posició <strong>de</strong> fluxos: el dipol. Donat a ∈ C, a = 0, po<strong>de</strong>m consi<strong>de</strong>rar el flux donat <strong>per</strong> la<br />

su<strong>per</strong>posició d’una font a z = z0 − ɛa i un xuclador d’igual potència a z = z0 + ɛa:<br />

Ω(z) = Q<br />

2π log(z + ɛa − z0) − Q<br />

2π log(z − ɛa − z0) = Q z + ɛa − z0<br />

log<br />

2π z − ɛa − z0<br />

(Q > 0). Si fem ɛ → 0 la font i el xuclador s’acosten. A fi que no <strong>de</strong>saparegui el potencial Ω cal fer<br />

tendir també Q → ∞ <strong>de</strong> manera que Q(ɛ)ɛa/π → µ (finit). Aleshores es comprova que el flux límit<br />

tindria<br />

Ω(z) = µ<br />

.<br />

z − z0<br />

Direm que µ és el moment <strong>de</strong>l dipol. Notem que si z0 = 0 i µ = 1 aleshores ψ = −y/(x 2 + y 2 ).<br />

5) Fluxos <strong>de</strong>finits <strong>per</strong> <strong>funcions</strong> meromorfes. Si z = z0 és un pol <strong>de</strong> la funció u − iv, en un entorn <strong>de</strong><br />

z = z0 podrem escriure (<strong>de</strong>senvolupant en sèrie <strong>de</strong> Laurent i integrant)<br />

Ω(z) = C−n<br />

µ<br />

+ · · · + +<br />

(z − z0) n z − z0<br />

Γ + iQ<br />

log(z − z0) + C0 + c1(z − z0) + . . .<br />

2πi<br />

Direm que el terme Γ+iQ<br />

2πi log(z − z0) <strong>de</strong>termina una font-vòrtex a z = z0. Direm també que<br />

µ<br />

z−z0<br />

és un<br />

dipol <strong>de</strong> moment µ i que els termes C−k<br />

(z−z0) k <strong>de</strong>terminen multipols d’ordre 2k i intensitat C−k a z = z0.<br />

3


Igualment, <strong>de</strong>senvolupant Ω(z) prop <strong>de</strong> z = ∞ tindrem<br />

Ω(z) = C−nz n + · · · + µz +<br />

Γ + iQ<br />

2πi<br />

log z + C1<br />

z<br />

i direm que a z = ∞ tenim una font <strong>de</strong> potència −Q, un vòrtex d’intensitat −Γ, un dipol <strong>de</strong> moment<br />

µ (les línies <strong>de</strong> corrent tenen la direcció <strong>de</strong> µ) i multipols d’ordre 2k i intensitat C−k. Noti’s que si<br />

C−n = · · · = C−2 = 0 aleshores el flux té una velocitat a l’infinit igual a µ.<br />

Exercicis ([3]).<br />

En els exercicis següents es dóna el potencial complex i es <strong>de</strong>mana: línies equipotencials (φ =Ct.),<br />

línies <strong>de</strong> corrent i representació gràfica, camp <strong>de</strong> velocitats V , punts singulars (potencies, intensitats,moments...),<br />

punts crítics ( V = 0) i comportament a z = ∞<br />

1.- Ω(z) = z n (en particular, n = 2, 3).<br />

2.- Ω(z) =<br />

Γ + iQ<br />

2πi<br />

log<br />

z − a<br />

z − b .<br />

3.- i) Ω(z) = z + R2<br />

R2<br />

; ii) Ω(z) = z −<br />

z z .<br />

4.- Ω(z) = 1<br />

.<br />

z2 <strong>5.</strong>- Ω(z) = log(z 2 − a 2 ) (a > 0).<br />

6.- Ω(z) = log z2 − a2 z2 (a > 0).<br />

+ a2 Exercici 7: Pot sortir una línia <strong>de</strong> corrent d’un punt que sigui: i) un vòrtex? ii) un dipol? iii) un<br />

vòrtex i un dipol a la vegada?<br />

Exercici 8: L’òval <strong>de</strong> Rankine és una regió d’aigües aïlla<strong>de</strong>s que es produeix quan hi ha una font<br />

a z = −1 un xuclador d’igual potència a z = 1 i la velocitat a z = ∞ és (1, 0). Mireu <strong>de</strong> dibuixar-lo...<br />

(numèricament).<br />

Mèto<strong>de</strong> <strong>de</strong> les imatges o <strong>de</strong> reflexió.<br />

Teorema: Suposem que Ω(z) és el potencial complex d’un flux a tot C <strong>de</strong>finit <strong>per</strong> les seves singularitats<br />

(fonts, vòrtexs, dipols...), i que aquestes són aïlla<strong>de</strong>s, i amb velocitat zero a l’infinit. Suposem<br />

que totes les sigularitats es troben al semiplà Im(z) > 0. Aleshores el flux produït <strong>per</strong> les mateixes<br />

singularitats a Im(z) > 0 <strong>per</strong>ò pel que Im(z) = 0 és una paret sòlida és Ω(z) + Ω(z).<br />

Dem.: S’ha <strong>de</strong> comprovar primer que Ω(z) és holomorfa (amb singularitats) i <strong>de</strong>sprés que Im(z) = 0<br />

és una línia <strong>de</strong> corrent <strong>de</strong> Ω(z) + Ω(z). <br />

Teorema: (<strong>de</strong> Milne-Thompson)Suposem que Ω(z) és el potencial complex d’un flux a tot C <strong>de</strong>finit<br />

<strong>per</strong> les seves singularitats (fonts, vòrtexs, dipols...), i que aquestes són aïlla<strong>de</strong>s, i amb velocitat donada<br />

a l’infinit. Suposem que totes les sigularitats es troben a |z| > R. Aleshores el flux produït <strong>per</strong> la<br />

introducció en el corrent <strong>de</strong>l disc sòlid |z| ≤ R (és a dir, el flux que <strong>per</strong> |z| > R té les mateixes<br />

4<br />

+ . . .


singularitats que Ω i la mateixa velocitat a z = ∞ <strong>per</strong>ò pel que |z| = R és una línia <strong>de</strong> corrent) és<br />

Ω(z) + Ω(R 2 /z).<br />

Dem.: Semblant a la d’abans. <br />

Observi‘s que en aquests dos teoremes es dóna un procediment <strong>per</strong> a fer que una <strong>de</strong>terminada<br />

corba sigui una línia <strong>de</strong> corrent introduint singularitats fictícies a l’altra banda <strong>de</strong> la corba. Aquestes<br />

son les singularitats anomena<strong>de</strong>s d’imatge.<br />

Exemple: (molt important) El potencial complex <strong>de</strong>l flux al voltant d’un cil·lindre <strong>de</strong> radi R<br />

centrat a z = 0 i amb velocitat (a, 0) a l’infinit és Ω(z) = a(z + R 2 /z). Vegi’s Figura 2.<br />

Figura 2: Línies <strong>de</strong> corrent <strong>de</strong>l flux al voltant <strong>de</strong>l cil·lindre (sense circulació).<br />

Una variació molt important d’aquest problema és la que correspon a que el flux tingui una<br />

circulació donada Γ al llarg d’una corba que encercli el cil·lindre. En aquest cas és<br />

Ω(z) = a(z + R 2 /z) + Γ<br />

log z. (6)<br />

2πi<br />

Es veu fàcilment que en aquest flux les partícules situa<strong>de</strong>s a |z| = R es mouen amb una velocitat<br />

angular promig igual a Γ/(2πR 2 ), i girant en sentit positiu si Γ > 0. Per aquesta raó es diu que (6)<br />

representa el flux al voltant d’un ci·lindre que gira amb velocitat angular Γ/(2πR 2 ). A la Figura 3 hi<br />

ha un croquis <strong>de</strong> les línies <strong>de</strong> corrent <strong>de</strong> (6) <strong>per</strong> dos valors diferents <strong>de</strong> Γ.<br />

Figura 3: Dos fluxos amb circulacions negatives. Més negativa a la dreta que a l’esquerra.<br />

5


Exercicis<br />

Usant el mèto<strong>de</strong> <strong>de</strong> les imatges trobeu els potencials complexos <strong>de</strong>ls fluxos següents, en els que es<br />

suposa que les parets <strong>de</strong>ls recintes són sòlids impenetrables.<br />

9.- Una font en el recinte Re(z) > 0. Un nombre finit <strong>de</strong> fonts o vòrtexs en el mateix recinte. Pot<br />

afegir-se una font en un punt <strong>de</strong> la frontera? Amb quin cabal? I un vòrtex?<br />

10.- Nombre finit <strong>de</strong> fonts o vòrtexs en el sector angular {0 < Arg(z) < π/4}.<br />

11.- El mateix en el recinte {0 < Arg(z) < π/4, |z| > R}<br />

12.- Una font en el recinte {|z| > R} amb velocitat donada a l’infinit.<br />

13.- Una font i un xuclador amb els mateixos cabals (<strong>per</strong>ò signes diferents) en el recinte {|z| < R}.<br />

Què passa si tenen cabals diferents?<br />

Mèto<strong>de</strong> <strong>de</strong> transformació conforme.<br />

El mèto<strong>de</strong> <strong>de</strong> la transformació conforme es basa en la observació següent. Suposem que W i W ′ són<br />

dos recintes <strong>de</strong>l pla complex i F : W → W ′ és una equivalència conforme, que s’estén diferenciablement<br />

fins a les fronteres respectives. Suposem que Ω ′ : W ′ → C és el potencial d’un flux a W ′ tangent a la<br />

seva frontera. Aleshores si <strong>de</strong>finim Ω ′ = Ω ◦ F resulta que Ω és el potencial complex d’un flux a W<br />

que també és tangent a la frontera <strong>de</strong> W .<br />

És fàcil <strong>de</strong> veure que <strong>per</strong> a les velocitats a l’infinit respectives V ′ i V es té que V ′ = V F ′ (∞). I<br />

també que si Ω ′ té un vòrtex o una font en el punt z = a ′ aleshores Ω té un vòrtex o una font, <strong>de</strong> les<br />

mateixes intensitat o potència, en el punt F −1 (a ′ ).<br />

Exercicis ([3]).<br />

14.- Trobeu el potencial complex <strong>de</strong>l flux al voltant d’un cil·lindre <strong>de</strong> secció el·líptica, amb velocitat<br />

donada a l’infinit i amb circulació Γ.<br />

1<strong>5.</strong>- Trobeu el potencial complex <strong>de</strong>l flux al voltant <strong>de</strong> la placa {|x| ≤ C, y = 0} amb velocitat a<br />

l’infinit donada.<br />

16.- Trobeu el potencial complex <strong>de</strong>l flux al voltant d’una paràbola, <strong>per</strong> l’interior i <strong>per</strong> l’exterior.<br />

17.- Trobeu el potencial complex <strong>de</strong>l flux inci<strong>de</strong>nt sobre les semirectes {1 < |x| < ∞, y = 0}.<br />

Forces i Paradoxa <strong>de</strong> d’Alembert.<br />

Teorema (1r T. <strong>de</strong> Blasius): La força F exercida <strong>per</strong> un flux <strong>de</strong> potencial complex Ω(z) sobre un<br />

cos submergit <strong>de</strong> frontera γ ve donada <strong>per</strong>:<br />

F = ρ i<br />

<br />

(dΩ/dz)<br />

2<br />

2 dz, (7)<br />

on ρ és la <strong>de</strong>nsitat <strong>de</strong>l fluid.<br />

Dem.: Usant (3) i (4) tenim<br />

γ<br />

<br />

F = − (p0 −<br />

γ<br />

ρ<br />

2 |V |2 )dℓ.<br />

Suposem γ = γ(s) parametritzada <strong>per</strong> la longitud <strong>de</strong> l’arc. Aleshores<br />

−<br />

L<br />

0<br />

p0(−iγ ′ (s))ds = i<br />

6<br />

L<br />

0<br />

p0γ ′ (s)ds = 0


L ρ<br />

F =<br />

0 2 |V |2 (−iγ ′ (s))ds<br />

F = i ρ<br />

2<br />

L<br />

|V |<br />

0<br />

2 γ ′ (s)ds.<br />

Usem ara que V (γ(s)) = γ ′ (s)|V |σ(s) amb σ(s) ∈ {−1, 1}, i tenim<br />

<br />

F = i ρ<br />

2<br />

= i ρ<br />

= i ρ<br />

2<br />

L<br />

|V |γ<br />

0<br />

′ (s)σ(s)|V |γ ′ (s)σ(s)γ ′ (s)ds =<br />

L<br />

(|V |γ<br />

2 0<br />

′ (s)σ(s)) 2 γ ′ (s)ds =<br />

L<br />

(V )<br />

0<br />

2 γ ′ (s)ds = i ρ<br />

<br />

(<br />

2 γ<br />

dΩ<br />

dz )2dz. Teorema (Kutta-Jukowsky): En un flux al voltant d’un obstacle amb velocitat donada a l’infinit,<br />

la força exercida sobre el cos ve donada <strong>per</strong> F = −iρV∞Γ, on Γ és la circulació <strong>de</strong>l flux al llarg d’una<br />

corba que envolti el cos.<br />

Dem.: S’usa el primer teorema <strong>de</strong> Blasius i es <strong>de</strong>forma el contorn d’integració fins a l’infinit. <br />

Corol·lari: (Paradoxa <strong>de</strong> D’Alembert). Un flux sense circulació no exerceix cap força sobre un<br />

obstacle submergit.<br />

Sustentació.<br />

El Teorema <strong>de</strong> Kutta-Jukowsky va inspirar la i<strong>de</strong>a <strong>de</strong> mirar d’interpretar les forces <strong>de</strong> sustentació<br />

(ja que no les d’arrossegament) <strong>de</strong>ls sòlids submergits en corrents, com ara les exerci<strong>de</strong>s a les ales<br />

<strong>de</strong>ls avions, com a conseqüència <strong>de</strong> l’aparició d’una circulació. Es tracta <strong>de</strong> l’anàleg als fluxos donats<br />

<strong>per</strong> (5) al voltant <strong>de</strong>l cil·lindre. Però es manté la paradoxa <strong>de</strong> que no po<strong>de</strong>n explicar-se les forces<br />

d’arrossegament. La causa d’aquesta paradoxa està en el caràcter no viscós <strong>de</strong>l mo<strong>de</strong>l que estem<br />

consi<strong>de</strong>rant.<br />

Llei <strong>de</strong> Coulomb i <strong>camps</strong> elèctrics plans.<br />

2.- Electrostàtica.<br />

La Llei <strong>de</strong> Coulomb diu que sobre dues càrregues elèctriques puntuals q1 i q2 situa<strong>de</strong>s en un medi<br />

tridimensional es produeixen unes forces <strong>de</strong> magnituds q1q2/d 2 , on d és la distància que les separa. Les<br />

forces estan dirigi<strong>de</strong>s <strong>de</strong> l’una a l’altra partícula, i són d’atracció si les càrregues tenen signes diferents,<br />

i <strong>de</strong> repulsió en cas contrari.<br />

Una altra manera <strong>de</strong> dir el mateix és dir que una càrrega q situada en el punt M = (ξ, η, ζ) crea<br />

en tots els punts P = (x, y, z) <strong>de</strong> l’espai el camp elèctric<br />

E(x, y, z) = q(<br />

x − ξ y − η z − ζ<br />

, , ), (8)<br />

R3 R3 R3 on R = (x − ξ) 2 + (y − η) 2 + (z − ζ) 2 , en el sentit que si situéssim una altra càrrega elèctrica q ′ en<br />

el punt P = (x, y, z) aleshores aquesta es veuria afectada <strong>per</strong> la força q ′ E(x, y, z).<br />

7


Fàcilment es comprova que si diem V = q/R = q/( (x − ξ) 2 + (y − η) 2 + (z − ζ) 2 ), aleshores<br />

E(x, y, z) = −( ∂V ∂V ∂V<br />

, , ) = −∇V, (9)<br />

∂x ∂y ∂z<br />

i <strong>per</strong> aquesta raó la funció V rep el nom <strong>de</strong> potencial elèctric creat <strong>per</strong> la càrrega q.<br />

Si tenim diverses càrregues q1, q2, . . . , qn situa<strong>de</strong>s en els punts M1, M2, . . . , Mn aleshores el camp<br />

elèctric creat seria la suma <strong>de</strong>ls <strong>camps</strong> elèctrics individuals:<br />

E(P ) = −∇(<br />

n<br />

j=1<br />

i es compliria el ben conegut Teorema <strong>de</strong> Gauss<br />

<br />

F · n dS = 4π <br />

si cap <strong>de</strong>ls punts Mj es troba sobre ∂Ω.<br />

∂Ω<br />

dj<br />

) (10)<br />

|Mj − P |<br />

Igualment, si tenim infinites càrregues distribuï<strong>de</strong>s <strong>de</strong> manera continua i uniforme sobre una recta<br />

<strong>per</strong>pendicular al pla (x, y), com ara la recta {(ξ, η, ζ)| − ∞ < ζ < ∞} po<strong>de</strong>m calcular el camp elèctric<br />

com la integral <strong>de</strong>ls <strong>camps</strong> elèctrics individuals. Suposem que q és la càrrega <strong>per</strong> unitat <strong>de</strong> longitud.<br />

mj∈Ω<br />

Per simplicitat, po<strong>de</strong>m suposar ξ = η = 0 i anomenar ρ = x2 + y2 , amb el que<br />

∞ x<br />

E(x, y, z) = − q (<br />

−∞ (ρ2 + ζ2 y<br />

,<br />

) 3/2 (ρ2 + ζ2 z − ζ<br />

,<br />

) 3/2 (ρ2 + ζ2 )<br />

qj<br />

(11)<br />

x y<br />

) dz = −2q( , , 0). (12)<br />

3/2 ρ2 ρ2 (Per calcular les dues primeres integrals pot fer-se en canvi tan α = ζ/ρ; la tercera és zero en el sentit<br />

<strong>de</strong>l valor principal.)<br />

L’expressió (12) ens dona un primer exemple <strong>de</strong> camp electric pla, en el sentit que no <strong>de</strong>pèn <strong>de</strong><br />

la tercera coor<strong>de</strong>nada i que la seva tercera component és zero. Com hem vist, els <strong>camps</strong> elèctrics<br />

plans po<strong>de</strong>n interpretar-se com els creats <strong>per</strong> càrregues distribuï<strong>de</strong>s contínuament i uniformement en<br />

la direcció <strong>per</strong>pendicular al pla. Com que prenent V (x, y) = 2q log (x − ξ) 2 + (y − η) 2 es segueix<br />

tenint que E = −∇V , direm que 2q log |P − M| és el potencial elèctric creat en el punt P = (x, y) <strong>per</strong><br />

una càrrega plana situada en el punt M = (ξ, η).<br />

Com que la funció φ(x, y) = log (x − ξ) 2 + (y − η) 2 és harmònica (φxx + φyy = 0) excepte <strong>per</strong><br />

(x, y) = (ξ, η), acceptem que en general la funció potencial elèctric v(x, y), que pot ser creat <strong>per</strong><br />

diverses càrregues, o <strong>per</strong> distribucions contínues d’aquestes, serà també harmònica, excepte a les<br />

regions ocupa<strong>de</strong>s <strong>per</strong> les càrregues. Si W és una regió sense càrregues, tindrem que el camp elèctric<br />

E = −∇v = −(vx + ivy) és el conjugat d’una funció holomorfa a W . I com que −(vx − ivy) és<br />

holomorfa, localment podrem escriure E = −(vx − ivy) = −iΩ ′ (z). A aquesta funció holomorfa Ω(z)<br />

s’en diu el potencial elèctric complex. La seva part imaginària v(z) és el potencial elèctric real, que<br />

hem consi<strong>de</strong>rat fins ara, i la seva part real u(z), que globalment pot ser multivaluada, rep el nom <strong>de</strong><br />

funció <strong>de</strong> força, <strong>per</strong>què és constant sobre les corbes integrals <strong>de</strong>l camp vectorial E. Les seves corbes<br />

<strong>de</strong> nivell s’anomenen les línies <strong>de</strong> força.<br />

En aquests termes, el potencial complex d’una càrrega <strong>de</strong> magnitud 2q situada en el punt z = a és<br />

Ω(z) = 2qi log(1/(z − a). De manera similar a com s’ha fet en el cas hidrodinàmic, es pot consi<strong>de</strong>rar<br />

un dipol <strong>de</strong> moment p, que situat en el punt z = a tindrà com a potencial la funció Ω(z) = ip/(z − a).<br />

Noti’s que p és un nombre complex, i que p és la direcció <strong>de</strong> l’eix <strong>de</strong>l dipol. En general, si z = a és un<br />

pol <strong>de</strong> la funció Ω ′ (z) i tenim el <strong>de</strong>senvolupament<br />

Ω(z) = c−n<br />

ip<br />

1<br />

+ · · · + + 2iq log<br />

(z − a) n z − a z − a + c0 + c1(z − a) + . . .<br />

8


el terme 2iq log(1/(z − a)) representa en el punt z = a una càrrega puntual plana <strong>de</strong> magnitud 2q,<br />

com la consi<strong>de</strong>rada més amunt. El terme ip/(z − a) representa un dipol <strong>de</strong> moment p (p és un nombre<br />

complex i p representa l’eix <strong>de</strong>l dipol. Els termes c−k/(z − a) k representen multipols d’ordre 2k.<br />

Igualment, si en el infinit es té<br />

Ω(z) = cnz n + · · · + ipz + 2iq log z + c0 + c−1<br />

z<br />

el terme 2iq log z representa a l’∞ una càrrega puntual plana <strong>de</strong> magnitud 2q, i el terme ipz un dipol<br />

<strong>de</strong> moment p.<br />

Que la funció u no sigui mai multivaluada fa que en electrostàtica no hi hagi l’anàleg al vòrtex <strong>de</strong><br />

l’hidrodinàmica. Com es diu clàssicament, les línies <strong>de</strong> força sempre neixen i moren a les càrregues.<br />

Exercicis ([3]). En els problemes següents a partir <strong>de</strong>ls potencials complexos donats s’han <strong>de</strong><br />

<strong>de</strong>terminar les <strong>funcions</strong> potencial i <strong>de</strong> força, el caràcter <strong>de</strong> les singularitats (inclòs ∞) i construir<br />

esquemàticament les famílies <strong>de</strong> les línies equipotencials i <strong>de</strong> força.<br />

18.- Ω(z) = 2iq log z−a<br />

z−b .<br />

19.- Ω(z) = 2qi log(z 2 − a 2 ) (a > 0).<br />

20.-Ω(z) = z ± R2<br />

z .<br />

Conductors en <strong>camps</strong> elèctrics. Funció <strong>de</strong> Green.<br />

Quan un cos conductor s’introdueix en el si d’un camp elèctric E, les seves càrregues lliures es<br />

mouen dins <strong>de</strong>l conductor <strong>per</strong> efecte <strong>de</strong>l camp i <strong>de</strong> la seva interacció mútua (en principi, suposem que<br />

el <strong>de</strong>splaçament d’aquestes no afecta a les càrregues llunyanes que han produït el camp E). La situació<br />

estacionària és que totes les càrregues lliures acaben situant-se sobre la su<strong>per</strong>fície <strong>de</strong>l conductor i <strong>de</strong><br />

tal manera que el nou camp elèctric E +E ′ (E ′ és el creat <strong>per</strong> aquestes càrregues lliures <strong>de</strong>l conductor)<br />

té un potencial elèctric que és constant sobre la su<strong>per</strong>fície <strong>de</strong>l conductor. Si el conductor està unit<br />

a terra, aquest valor constant ha <strong>de</strong> ser zero. Per tant, po<strong>de</strong>m dir que el camp elèctric resultant<br />

a l’exterior <strong>de</strong>l conductor és aquell que té les mateixes singularitats que E (les mateixes càrregues<br />

exteriors) <strong>per</strong>ò modificat <strong>de</strong> manera que v = constant a la frontera <strong>de</strong>l conductor.<br />

La funció <strong>de</strong> Green G(z, a) d’un recinte D és una funció molt útil en la resolució d’equacions en<br />

<strong>de</strong>riva<strong>de</strong>s parcials. És una funció real que es <strong>de</strong>fineix en funció <strong>de</strong> la variable z ∈ D i <strong>de</strong>l paràmetre<br />

a ∈ D <strong>per</strong><br />

⎧<br />

1<br />

⎪⎨ ∆(G(z, a) − log ) = 0,<br />

|z − a|<br />

G(z, a) = 0 si z ∈ ∂D<br />

(13)<br />

⎪⎩<br />

G(z, a) fitada prop <strong>de</strong> z = ∞<br />

(la tercera condició és automàtica si D és fitat; si es vol admetre que ∞ ∈ D, llavors s’ha d’excloure el<br />

cas a = ∞ en la tercera condició). Del que s’ha dit anteriorment, si Ω = u + iv és el potencial elèctric<br />

complex <strong>de</strong>gut a una càrrega situada a z = a ∈ D, quan suposem que C \ D és un conductor unit a<br />

terra, tindrem v = 2qG(z, a).<br />

Exercicis ([3]). En els problemes següents es <strong>de</strong>mana trobar els potencials elèctrics complexos a<br />

partir <strong>de</strong> les seves singularitats i en certs recintes. Es suposa que l’exterior <strong>de</strong>l recinte és un conductor.<br />

21.- Una càrrega en el punt a en el recinte Im(z) > 0.<br />

22.- Una càrrega en el punt a en el recinte |z| < R.<br />

23.- Una càrrega en el punt a en el recinte |z| > R.<br />

9<br />

+ . . .


24.- Una càrrega en el punt a = ∞ en el recinte (x/α) 2 + (y/β) 2 > 1.<br />

3.- Problemes Tèrmics.<br />

La funció u(x, y) que ens dona la tem<strong>per</strong>atura en cada punt en un problema pla en règim estacionari<br />

és la incògnita principal amb la que treballem. La conductivitat tèrmica <strong>de</strong>l medi es representa pel<br />

coeficient k > 0, i <strong>de</strong>finim vector flux <strong>de</strong> calor <strong>per</strong> Φ = −k∇u. El vector flux <strong>de</strong> calor té la propietat<br />

que la quantitat <strong>de</strong> calor que travessa <strong>de</strong> dins a fora una corba simple Γ <strong>de</strong>l pla i <strong>per</strong> unitat <strong>de</strong> temps<br />

(i unitat <strong>de</strong> longitud en la tercera coor<strong>de</strong>nada) és<br />

<br />

Γ<br />

<br />

Φ · n |dz| = −k<br />

si és que dins <strong>de</strong> la corba no hi ha fonts <strong>de</strong> calor.<br />

Γ<br />

∂u<br />

∂n<br />

<br />

dℓ =<br />

int.(Γ)<br />

∆u dxdy<br />

Per tant, a les regions on no hi ha fonts <strong>de</strong> calor, u és harmònica, i consi<strong>de</strong>rant la seva funció<br />

harmònica conjugada v(x, y) po<strong>de</strong>m <strong>de</strong>finir el potencial complex <strong>de</strong>l camp tèrmic <strong>per</strong> Ω(z) = u + iv.<br />

Les corbes u =ct. són les isotermes i les v =ct. s’anomenen linies <strong>de</strong> corrent. També es té que<br />

Φ = −kΩ ′ (z).<br />

Si en un entorn <strong>de</strong> z = a tenim<br />

Ω(z) = c−n<br />

c−1 q<br />

+ · · · + +<br />

(z − a) n z − a 2πk log<br />

1<br />

z − a + c0 + c1(z − a) + . . .<br />

el terme q/(2πk) log(1/(z − a)) representa en el punt z = a una font puntual <strong>de</strong> calor <strong>de</strong> potència q, i<br />

el terme c−1/(z − a) un dipol.<br />

Trobeu la distribució <strong>de</strong> tem<strong>per</strong>atures u(x, y) en el recinte |x| < a, y > 0 suposant que la tem<strong>per</strong>atura<br />

és zero a la frontera i que hi ha situada una font <strong>de</strong> potència q en un punt <strong>de</strong> la forma<br />

z0 = ih, h > 0.<br />

Exercicis.<br />

2<strong>5.</strong>- Trobeu la distribució <strong>de</strong> tem<strong>per</strong>atures en el interior <strong>de</strong> l’anell circular r1 < |z| < r2 si es<br />

sap que a l’interior <strong>de</strong> l’anell hi ha una font puntual <strong>de</strong> calor i que en les circumferències frontera la<br />

tem<strong>per</strong>atura pren uns valors constants coneguts u1 i u2.<br />

26.- Problema <strong>de</strong>l cable enterrat. Un cable d’alta tensió <strong>de</strong> radi r està enterrat, i hi ha una distància<br />

d > r <strong>de</strong>l centre <strong>de</strong>l cable a la su<strong>per</strong>fície <strong>de</strong>l terreny. Per efecte Joule, el cable emet q unitats <strong>de</strong> calor<br />

<strong>per</strong> unitat <strong>de</strong> temps i unitat <strong>de</strong> longitud <strong>de</strong>l cable. Suposant coneguda la conductivitat tèrmica k <strong>de</strong>l<br />

terreny, i la tem<strong>per</strong>atura <strong>de</strong> l’atmosfera, trobeu la tem<strong>per</strong>atura a la su<strong>per</strong>fície <strong>de</strong>l cable. (Aquesta és<br />

la tem<strong>per</strong>atura que hauran <strong>de</strong> suportar sense fer-se malbé els materials aïllants que l’envolten.)<br />

Referències<br />

[1] Chorin, A.J. i J. Mars<strong>de</strong>n: A Mathematical Introduction to Fluid Mechanics. Springer.<br />

[2] Spiegel, M.R.: Variable Compleja. McGraw-Hill. Col. Schaum.<br />

[3] Volkovyski, L., G. Lunts i I. Aramanovich: Problemas sobre la teoria <strong>de</strong> funciones <strong>de</strong> variable<br />

compleja. MIR, Moscou.<br />

10

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!