33 0 331KB
TD 17 – Électromagnétisme
BLAISE PASCAL PT 2018-2019
Correction
Électromagnétisme en régime variable Exercices
Exercice 1 : Manipulation des équations de Maxwell 1
Par définition,
#” ∂ Ex div E = =0 ∂x ce qui est conforme à l’équation de Maxwell-Gauss car ρ = 0 par hypothèse. 2
De même, #” ∂ By div B = =0 ∂y
ce qui est là aussi conforme à l’équation de Maxwell-Thomson. 3
Compte tenu de la forme des champs, # ” #” rot E = f 0 (z) e−t/τ #” ey
#” 1 ∂B = g(z) e−t/τ #” ey − ∂t τ
et
d’où on déduit de l’équation de Maxwell-Faraday f 0 (z) =
4
1 g(z) τ
Compte tenu de la forme des champs, # ” #” rot B = −g 0 (z) e−t/τ #” ex
#” ∂E 1 = − f (z) e−t/τ #” ex . ∂t τ
et
d’où on déduit de l’équation de Maxwell-Ampère g 0 (z) =
5
ε0 µ0 f (z) e−t/τ . τ
En dérivant par rapport à z l’équation obtenue par Maxwell-Faraday, on obtient f 00 (z) =
1 0 g (z) τ
et en combinant avec la question précédente on obtient f 00 (z) −
ε0 µ0 f (z) = 0 . τ2
Il s’agit d’une équation différentielle du second ordre à coefficients constants portant sur la seule variable z. Son polynôme caractéristique admet pour racines r √ µ0 ε0 µ0 ε0 r± = ± = ± 2 τ τ
1/8
Étienne Thibierge, 17 mars 2019, www.etienne-thibierge.fr
Correction TD 17 : Électromagnétisme en régime variable
Blaise Pascal, PT 2018-2019
Dans un cas aussi simple, pas besoin de passer par le discriminant ! La solution particulière est nulle et la solution homogène est √ √ µ0 ε0 µ0 ε0 f (z) = α exp − z + β exp z τ τ Les conditions données par l’énoncé imposent que β = 0 et que f (z = 0) = E0 = α. On en déduit : 1 √ #” ex . E = E0 exp − (t + µ0 ε0 z) #” τ 6
D’après ce qui précède, 1 √ √ g(z) = τ f 0 (z) = −E0 µ0 ε0 exp − (t + µ0 ε0 z) τ
d’où on déduit 1 √ √ #” B = −E0 µ0 ε0 exp − (t + µ0 ε0 z) #” ey . τ
Exercice 2 : Chauffage par induction 1 Voir la figure 1. En orientant l’axe z par rapport au courant i (règle de la main droite), le champ créé par le solénoïde s’écrit #” B = µ0 ni #” ez avec n le nombre de spires par unité de longueur. L’expression est identique dans l’ARQS et en magnétostatique. i(t) z
Figure 1 – Dispositif de soudure inductive. 2 L’équation de Maxwell-Faraday indique qu’un champ magnétique dépendant du temps est source de champ électrique. Au sein du métal, ce champ électrique va faire apparaître une densité de courant #” et céder de l’énergie au métal par effet Joule. Une autre justification moins orientée « équations » consisterait à dire que le métal est un conducteur placé dans un champ magnétique variable, et est donc le siège d’un phénomène d’induction. Ce sont les courants induits (les courants de Foucault) qui vont chauffer le matériau par effet Joule. 3
D’après l’équation de Maxwell-Faraday, #” ∂B di # ” #” rot E = − = −µ0 n #” e z = +µ0 nI0 ω sin(ωt) #” ez ∂t dt
et en utilisant le formulaire d’où on déduit
# ” #” 1 ∂ rEθ #” ez rot E = r ∂r ∂ rEθ = µ0 nI0 ω sin(ωt) . ∂r
En intégrant par rapport à r, rEθ = µ0 nI0 ω sin(ωt)
r2 + f (t) . 2
Comme on intègre une dérivée partielle pour laquelle certaines variables (ici t) sont fixées, alors la constante d’intégration est une « constante partielle », c’est-à-dire une fonction des variables fixées. On montre que la fonction f est nulle en se plaçant en r = 0 : comme il n’y a « rien de spécial » en ce point, la composante Eθ y prend une valeur finie. Finalement, #” 1 E = µ0 nrI0 ω sin(ωt) #” eθ . 2
2/8
Étienne Thibierge, 17 mars 2019, www.etienne-thibierge.fr
Correction TD 17 : Électromagnétisme en régime variable
Blaise Pascal, PT 2018-2019
L’analyse des invariances est identique à celle menée en électrostatique. En revanche, celle des symétries est beaucoup plus subtile à cause du champ magnétique variable, et hors de portée en PT : comme #” # ” #” c’est rot E qui intervient dans l’équation de Maxwell-Faraday, alors un plan de symétrie de ∂ B/∂t est #” plan d’antisymétrie de E ... ce qui est la propriété inverse par rapport à la charge électrique, car ρ #” apparaît dans l’équation de Maxwell-Gauss qui porte sur div E. 4
D’après la loi d’Ohm locale, #” 1 #” j = γ E = µ0 γnrI0 ω sin(ωt) #” eθ . 2
On en déduit la puissance volumique dissipée par effet Joule, #” #” 1 P = j · E = γµ02 n2 r2 I02 ω 2 sin2 (ωt) . 4 Elle est maximale lorsque r est maximal, c’est donc sur les bords du cylindre que le métal chauffe le plus et qu’il fond en premier.
Exercice 3 : Émission radioactive 1 L’amas émet des charges positifs, il est donc logique que sa charge soit négative et décroissante. L’allure exponentielle est typique de la radioactivité. À l’instant initial, on a q(t = 0) = 0 : aucun atome ne s’est désintégré donc l’amas est neutre. 2 L’émission étant isotrope, le vecteur densité de courant est forcément à symétrie sphérique : #” = jr (r) #” e r . Ainsi, #” tout plan contenant (OM ) est plan de symétrie de la distribution de courant, donc B(M, t) doit être orthogonal à chacun de ces plans et est donc forcément nul : #” #” B(M, t) = 0 . L’équation de Maxwell-Ampère indique qu’un champ électrique variable est également source de champ #” magnétique : dans un cas plus général, il faut donc également étudier les plans de symétrie de ∂ E/∂t. Cependant, d’après le théorème de superposition, les propriétés de symétrie doivent être vérifiées pour #” #” chaque « famille » de sources. Les symétries des courants imposent d’avoir B = 0 , le champ électrique ne pourra donc rien y changer. • Symétries : la distribution de charge possède les mêmes symétries que la distribution de courant, mais #” comme E doit être inclus dans un plan de symétrie il n’est pas nul mais radial, 3
#” E(M, t) = Er (M, t) #” er . L’équation de Maxwell-Faraday indique qu’un champ magnétique variable est également source de champ électrique : dans un cas plus général, il faut donc également étudier les plans de symétrie #” #” #” de ∂ B/∂t ... mais ce n’est pas la peine ici car B = 0 , tout se passe exactement comme en statique. #” • Invariances : à tout instant la distribution de charge est invariante par rotation autour de l’amas, donc E ne dépend que de r. • Théorème de Gauss, qui s’applique aussi en régime variable : on raisonne sur une sphère de rayon r centré sur l’amas. ‹ #” # ” . calcul du flux : E · dS = 4πr2 Er . . calcul de la charge intérieure : il faut un temps r/v0 pour qu’une particule α sorte de la sphère de rayon r, donc la charge contenue dans la sphère à un instant t est la charge de l’amas à l’instant t − r/v0 . r Qint = q t − v0 • Conclusion : #” E(r, t) =
Q0 1 r #” exp − t − er . 4πε0 r2 τ v0 3/8
Étienne Thibierge, 17 mars 2019, www.etienne-thibierge.fr
Correction TD 17 : Électromagnétisme en régime variable 4
Blaise Pascal, PT 2018-2019
D’après l’équation de Maxwell-Gauss, #” ρ(r, t) = ε0 div E t r 1 Q0 d exp − + r2 4π dr τ v0 τ 1 Q0 t r × = exp − + 4πr2 v0 τ τ v0 τ Q0 t r ρ(r, t) = exp − + 4π v0 τ r2 τ v0 τ =
D’après l’équation de Mawxell-Ampère, #” 1 # ” #” ∂E #” (r, t) = rot B − ε0 µ0 ∂t Q0 ∂ t r #” #” = 0− exp − + er 4πr2 ∂t τ v0 τ 1 Q0 t r #” × − exp − + =− er 2 4πr τ τ v0 τ Q0 t r #” #” (r, t) = exp − + er . 4π τ r2 τ v0 τ #” On retrouve j = ρv0 #” e r ... on aurait pu s’y attendre !
Exercice 4 : Modélisation électromagnétique de la charge d’un condensateur Le condensateur étudié est schématisé figure 2. z
#” Π
e
#” ⊗B
#” E
#” Π a
Figure 2 – Schéma des notations. 1
D’après la loi des mailles, U0 = Ri + u
et avec la loi de comportement du condensateur on en déduit RC
du + u = U0 dt
ce qui se met sous forme canonique du 1 U0 + u= dt τ τ
avec
τ = RC .
Cette équation admet U0 comme solution particulière, et les solutions homogènes sont de la forme A e−t/τ , ainsi u(t) = A + U0 e−t/τ Comme le condensateur est initialement déchargé, alors u(t = 0) |{z} = 0 |{z} = A + U0 CI
sol
4/8
Étienne Thibierge, 17 mars 2019, www.etienne-thibierge.fr
Correction TD 17 : Électromagnétisme en régime variable
Blaise Pascal, PT 2018-2019
d’où on a finalement u(t) = U0 1 − e−t/τ . 2
Par définitions, la charge stockée sur le condensateur vaut q = σπa2 = Cu = ε0
et ainsi
πa2 u e
d’où
σ=
ε0 u e
#” u E = #” ez . e
3 Les sources de champ magnétique sont les courants et les variations de champ électrique. Il n’y a pas de courant entre les armatures, la seule source de champ magnétique est donc le champ électrique, ou plus exactement ses variations. D’après la question précédente, #” ∂E 1 du #” = ez . ∂t e dt On se place en un point M quelconque. Le plan passant par M et contenant l’axe (Oz) (plan M, #” e r , #” e z ) est plan de #” #” symétrie de ∂ E/∂t, B est donc orthogonal à ce plan, #” B = Bθ (M ) #” eθ . #” C’est l’équation de Maxwell-Ampère qui montre que les symétries de ∂ E/∂t doivent être traitées comme #” celles de pour prévoir le champ magnétique, #” ∂E # ” #” #” rot B = µ0 + ε0 µ0 . ∂t 4
L’équation de Maxwell-Thomson donne
1 ∂Bθ =0 r ∂θ #” c’est-à-dire que B est indépendant de θ. De plus, l’équation de Maxwell-Ampère donne −
∂ Bθ #” 1 ∂ (rBθ ) #” 1 du #” er + e z = ε0 µ0 ez . ∂z r ∂r e dt
On en déduit d’une part que ∂ Bθ =0 ∂z #” c’est-à-dire que B est indépendant de z, et d’autre part que 1 ∂ (rBθ ) 1 du = ε0 µ0 . r ∂r e dt Cette relation donne
∂ (rBθ ) 1 du = ε0 µ0 r ∂r e dt
ce qui s’intègre en
1 du r2 + f (t) . e dt 2 #” Le champ B ne pouvant diverger en r = 0 où ses sources n’ont pas de singularité, la fonction f (t) est forcément nulle. On en déduit donc 1 du #” #” 1 r eθ . B = ε0 µ0 2 e dt rBθ = ε0 µ0
Comme on intègre une dérivée partielle, la « constante » d’intégration est en fait une fonction qui dépend des variables laissées constantes dans la dérivée partielle : ici, f doit a priori dépendre de z, θ #” et t ... mais elle ne dépend que de t car il a été montré par ailleurs que B ne dépend ni de z ni de θ. 5
Par définition,
#” #” u #” 1 1 du #” #” E ∧ B Π= = ez ∧ ε0 r eθ µ0 e 2 e dt 5/8
Étienne Thibierge, 17 mars 2019, www.etienne-thibierge.fr
Correction TD 17 : Électromagnétisme en régime variable
Blaise Pascal, PT 2018-2019
ce qui se simplifie en
1 1 du #” Π = − ε0 2 u r #” er . 2 e dt Le flux du vecteur de Poynting correspond à une puissance sortant, la puissance entrant dans le condensateur s’écrit donc (on intègre sur la surface en pointillés verts de la figure 2) ‹ #” # ” P=− Π · dS ¨condensateur ¨ ¨ =− Πr #” e r · dS #” ez − Πr #” e r · (−dS #” ez ) − Πr #” e r · dS #” er sup
sup
lat
1 1 du = 0 + 0 + ε0 2 u a × 2πa e 2 e dt 1 du P = ε0 u πa2 . e dt #” Cette puissance vient de l’extérieur du condensateur, et la direction de Π montre qu’elle entre par les parois latérales. Un bilan analogue appliqué au générateur montrerait qu’il cède de la puissance. On peut en fait montrer que le générateur cède une puissance 2P, la moitié étant stockée par le condensateur et l’autre moitié dissipée par la résistance. 6
On peut réécrire la puissance reçue sous la forme P=
du ε0 πa2 du u = Cu e dt dt
P=
soit
d dt
1 2 Cu 2
.
On reconnaît l’expression de l’énergie stockée dans le condensateur. 7
L’énergie reçue s’écrit
ˆ E=
∞
P dt =
0
1 2 Cu 2
∞ 0
Avec l’expression de u déterminée à la première question (ou avec la condition initiale et la solution particulière), on obtient 1 E = CU02 . 2
Annales de concours Exercice 5 : Bilan de puissance d’un conducteur
[oral banque PT]
1 Il s’agit de la conductivité électrique, notée γ. Celle du cuivre, l’un des meilleurs conducteurs, est de 6 · 107 S · m−1 ; à tire de comparaison celle de l’eau de mer n’est que de l’ordre de 5 S · m−1 . 2
D’après la loi d’Ohm locale,
#” #” j = γE ,
d’où on déduit l’intensité par calcul du flux au travers d’une section transverse, c’est-à-dire un cercle, ¨ #” # ” I= j · dS soit I = πR2 γE . Pour le calcul du champ magnétique, plaçons nous en un point M repéré par ses coordonnées cylindriques, voir figure 3. • Invariances et symétries : . La distribution de courant est invariante par translation le long de l’axe (Oz) et par rotation autour de cet axe, le champ magnétique ne dépend donc que du rayon r. . Le plan passant par M est contenant l’axe Oz est plan de symétrie de la distribution de courant, le champ magnétique lui est donc orthogonal : on en déduit qu’il est porté par #” uθ . . Conclusion : #” B = Bθ (r) #” uθ . 6/8
Étienne Thibierge, 17 mars 2019, www.etienne-thibierge.fr
Correction TD 17 : Électromagnétisme en régime variable
Blaise Pascal, PT 2018-2019
#” B M C #” E z
Figure 3 – Bilan de puissance électromagnétique d’un cylindre. • Théorème d’Ampère : Raisonnons sur un contour C circulaire de rayon r centré sur l’axe (Oz) et passant par M . Comme le champ est de même norme sur ce contour, alors ˛ #” #” B · d` = 2πr Bθ (r) C
Le courant enlacé dépend de la position de M : . si r > R alors Ienl = I = πR2 γE ; . si r < R alors Ienl = πr2 j = πr2 γE Finalement, µ0 πR2 γE #” uθ #” 2πr B= 2 µ0 πr γE #” uθ 2πr d’où en simplifiant 2 µ0 R γE #” uθ #” 2r B= µ0 r γE #” uθ 2 3
r≥R r≤R
r≥R r≤R
La puissance volumique dissipée par effet Joule s’écrit #” #” pJ,vol = j · E = γE 2 .
La puissance dissipée dans le tronçon vaut donc PJ = γE 2 `πR2 . 4
Le vecteur de Poynting s’écrit
#” #” #” E ∧ B Π= . µ0 #” #” #” Comme E est porté par #” uz et B par #” uθ alors Π est porté par #” ur . 1 µ0 r γE #” rγE 2 #” #” Π= ×E× ( uz ∧ #” uθ ) = − ur . µ0 2 2 Enfin, la puissance rayonnée au travers des parois du tronçon de cylindre est égale au flux sortant du vecteur de Poynting en r = R. Sur les surfaces verticales du cylindre (colorées de la figure 3), son flux est nul. Ainsi, ‹ ¨ #” # ” #” Pray = Π · dS = Π · dS #” ur , S
d’où Pray = −
Slat
RγE 2 × 2πR` 2
soit Pray = −γE 2 πR2 ` .
7/8
Étienne Thibierge, 17 mars 2019, www.etienne-thibierge.fr
Correction TD 17 : Électromagnétisme en régime variable
Blaise Pascal, PT 2018-2019
La puissance rayonnée est négative, ce qui signifie que l’énergie ne sort pas mais entre dans le cylindre au travers des parois. 5
On constate l’égalité des deux termes, PJ + Pray = 0 .
Cela signifie que toute la puissance dissipée par effet Joule dans le cylindre y entre par rayonnement, ce qui est logique puisqu’il n’y a pas d’autre source de puissance dans le tronçon de cylindre étudié. Ce résultat est par ailleurs cohérent avec l’équation intégrale de Poynting, ‹ dEem #” # ” =− Π · dS − PJ = 0 , dt car en régime stationnaire l’énergie électromagnétique Eem contenue dans le système est constante.
8/8
Étienne Thibierge, 17 mars 2019, www.etienne-thibierge.fr