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.
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