2019 PC Sujets Pour Oral de Physique

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 182

École Polytechnique - ESPCI Concours d’admission 2019

Filière PC

Épreuves orales de Physique

Nous réunissons dans ce document 59 sujets d’oral de physique correspondant à environ


70% des sujets donnés. Ils sont accompagnés de propositions de résolution et d’indications
sur la manière dont l’oral a pu être orienté, permettant aux futurs candidats de se préparer
à cette épreuve. Ces éléments de résolution sont fournis sans garantie d’absence d’erreurs
typographiques.
Induction et ressorts

On considère une barre métallique conductrice de masse m, longueur `, suspendue hori-


zontalement par deux ressorts (raideur k) attachés à ses extrémités. Les ressorts sont des
bobines d’inductance L. La barre reste toujours dans la zone où est appliqué un champ
magnétique constant uniforme, horizontal et orthogonal à la direction de la barre.

Quelle est l’équation du mouvement ?

Quelle est l’évolution de la position de la barre suite à un petit déplacement vers le bas ?
Et celle du courant ?

Faire un bilan énergétique et interpréter ce qu’il se passe.


Eléments de solution

On considère un axe z orienté vers le bas. On note ż = v.

~ soit on utilise la loi de Lenz pour


1) Soit on utilise que le champ électromoteur est ~v ∧ B,
trouver que
e = vB` . (1)
L’équation du circuit électrique est

e = 2Ldi/dt = −vB` . (2)

Ici il peut être utile de se demander le sens du courant en fonction de la vitesse. Le PFD
donne ensuite (on utilise la force de Laplace et le poids)

mv̇ = mg − i`B − 2k(z − z0 ) (3)

où z0 est l’altitude où les ressorts n’exercent pas de force de rappel. En dérivant on trouve
 2 2 
` B
mv̈ + + 2k v = 0 . (4)
2L

On posera donc naturellement


 
`2 B 2 2k
ω2 = + (5)
2Lm m

2) La position d’équilibre est en


mg
zeq = z0 + . (6)
2k

On pose z = zeq + Z et on a Ż = v. C’est un peu plus pénible que d’habitude pour


résoudre avec les condition initiales. Initialement on a Z(t = 0) = Zi , V (t = 0) = 0, mais
on a aussi une condition sur V̇ initialement qui est en utilisant (3) V̇ (t = 0) = −2kZi /m.
Donc on a (en posant ω02 = 2k/m)

ω02 Zi
V =− sin(ω0 t) . (7)
ω
 
ω02
Z = Zi 1 + 2 [cos(ωt) − 1] . (8)
ω
Pour le courant on intègre (2) avec la condition initiale i(t = 0) = 0 et donc

B`Zi k
i = I0 [1 − cos(ωt)] I0 = (9)
2kL + `2 B 2 /2
Les oscillations sont plus petites que dans le cas sans champ magnétique. On va com-
prendre énergétiquement. Mais on peut déjà voir que les forces de Laplace freinent quand
on monte et accélèrent quand on descend car elles sont dirigées vers le bas.

3) Il faut multiplier (3) par v et intégrer pour faire apparaître


 
d 1 2 2 2
mv − mgz + Li + k(z − z0 ) = 0 (10)
dt 2

Par simplicité on discute sans regarder le poids et on pose g = 0 et donc z0 = zeq . On a


donc
d  
Ec + Li2 + kZ 2 = 0 . (11)
dt
On n’a pas de facteur 1/2 dans les deux derniers termes, car on a deux ressorts. Ce qu’il
faut voir c’est que Li2 est maximal quand la barre est au plus haut. C’est pour ça que
l’amplitude des oscillations est plus faible. Une partie de l’énergie cinétique est prise par
cette énergie d’induction quand on monte, mais elle est restituée quand on descend.
Ondes acoustiques sphériques

On modélise une source d’ondes acoustiques par une sphère dont le rayon peut varier.
Le rayon de celle-ci au repos est r0 et lorsque la source est active son rayon oscille et est
donné par r = r0 +  cos(ωt).

— On adopte une notation complexe. Montrer que δP = δP ei(ωt−kr) avec δP = P0 /r


est solution de l’équation d’ondes acoustiques pour la pression. Justifier qu’il s’agit
d’une onde sphérique.

— Quelle est la forme du champ de vitesses ?

— Quelle est l’intensité acoustique ? Quelle est la puissance totale transmise ?

— Quelle est l’efficacité de cette source ?


Eléments de solution

1) Avec l’expression du Laplacien en sphérique sans dépendance angulaire ∆f (r) =


1/r2 ∂r (r2 ∂r f ) on vérifie
∂tt δP − c2 ∆δP = 0 (1)
si ω = kc. On a bien la phase i(ωt − kr) pour la propagation selon la direction radiale et
la décroissance en 1/r caractéristique des ondes sphériques.

2) Avec l’équation d’Euler linéarisée on a

ρ0 ∂t vr = −∂r δP (2)

et donc avec la notation complexe

vr = vr ei(ωt−kr) (3)

on trouve  
1 i
vr = 1− δP . (4)
ρ0 c kr
On voit qu’on a deux termes. Un qui correspond au champ proche et qui est négligeable
quand kr  1 et la solution en champ lointain qui va en 1/r. Mais maintenant on veut
mettre les conditions aux bords de la sphère. Le fluide au niveau de la sphère doit avoir
la vitesse imposée par celle-ci et donc avec r = r0 + eiωt on a vrsphere = iωeiωt . On va
toujours supposer que  est très petit devant la longueur d’onde et donc considérer que
cette vitesse est la vitesse du fluide en r0 et donc on obtient de (4)

kr0
P0 = −ωr0 cρ0 (5)
1 + ikr0
Ici on ne suppose pas kr0  1 (et donc la source pas dans le champ proche) et on donne
bien le cas général.

3) L’intensité acoustique est 1/2Re(δP vr ? ) donc

|P0 |2 ρ0 c (ωr0 )2 ω 2 r02 /c2


I= = (6)
2ρ0 cr2 2 r2 1 + ω 2 r02 /c2

On a à toute distance une décroissance en 1/r2 , et donc une puissance totale à travers une
sphère de rayon r constante. En effet on n’a aucun phénomène dissipatif. La puissance
totale est 4πr2 I.

4) Pour définir une efficacité, il faut avoir en tête la propagation pour une onde plane
progressive. On a en effet que l’intensité acoustique est la densité d’énergie acoustique
multipliée par la vitesse du son. On va donc d’abord calculer la densité d’energie acous-
tique moyenne

1 |δP |2 1 1
u = hρ0 |vr |2 + i = ρ|vr |2 + 2 |δP |2 (7)
2 ρ0 c2 4 4ρc
 r 2 ρ (ω)2  1 + 2(kr)2   (kr )2 
0 0 0
= (8)
r 2 2(kr)2 1 + (kr0 )2

Pour répondre à la question de l’efficacité on n’a jamais besoin de cette dernière expression
générale (qui pourtant est instructive car on voit que pour kr  1, la décroissance en
r n’est pas juste en 1/r2 ). Pour évaluer l’efficacité on regarde juste la densité d’énergie
acoustique au niveau de la sphère donc en r0 .
Au niveau de la sphère on a alors (il faut pas spécialiser le calcul précédent mais
plutôt le faire directement en r0 ce qui est beaucoup plus court)
 
2 ω 2 ρ0 1 + 2k 2 r02
u|r=r0 = . (9)
2 2 + 2k 2 r02

On définit l’efficacité comme


I
α≡ (10)
uc
évalué en r0 . Si c’est unité c’est que tout se passe comme les ondes planes progressives.
Sinon c’est pas le cas. On trouve

2(kr0 )2
α= . (11)
1 + 2(kr0 )2

Si la longueur d’onde est très inférieure à la taille de la source (kr0  1) on a une


efficacité à un et c’est une bonne source. Dans le cas contraire on a une suppression de
l’efficacité en (r0 /λ)2 . Pour des grandes longueurs d’ondes et donc basses fréquences, il
faut des grosses enceintes sinon cela ne fonctionne pas bien.
Tore de fluide

On considère un liquide parfait incompressible en forme de tore, c’est à dire de bouée.


Dans un premier temps on étudie son évolution en considérant l’effet de la tension su-
perficielle. On ne prend pas en compte l’effet du poids et on suppose qu’au cours de son
évolution il conserve toujours une forme de tore. On rappelle que l’énergie associée à la
tension superficielle est E = νA où A est l’aire de la surface libre et ν le coefficient de
tension superficielle. On suppose que le rayon r du grand cercle qui passe par les centres
des sections est bien plus grand que le rayon des sections R du tore.

— Déterminer l’évolution du système. Interpréter physiquement.

— Quel est le comportement aux temps courts ? Et aux temps longs ?

— Comment peut on le stabiliser avec une rotation initiale ?

Aide 1 : Le volume et l’aire du tore


p sont√respectivement 2π 2 rR2 et 4π 2 Rr. p
√ √
Aide 2 : La fonction f (x) = −8/3 1 − x(1 + x/2) est telle que f 0 (x) = 1/ 1 − x.
Eléments de solution


On va se rendre compte que c’est semblable à un problème de force centrale en 1/ r.

1) On passe forcément par la conservation de l’énergie car c’est le seul aspect connu de
la tension superficielle et de plus le problème n’a qu’un degré de liberté. On a affaire à
un fluide incompressible et donc la conservation du fluide va nous donner une relation.
Comme le volume est ∝ R2 r on obtient
r
ri
R = Ri . (1)
r

Etant donné le rapport d’aspect du problème (r  R), les vitesses sont essentiellement
radiales, c’est à dire selon ~er en coordonnées cylindriques naturelles. On a donc (en posant
ρ = 1)
1
Ec = V ṙ2 = π 2 R2 rṙ2 . (2)
2
Et pour l’énergie de tension superficielle

Ep = νA = ν4π 2 Rr . (3)

Donc sans dissipation visqueuse (fluide parfait) la conservation donne


4ν 4ν
ṙ2 = − . (4)
Ri R
On a commencé le mouvement sans vitesse initiale pour fixer la constante. On remarque
qu’en dérivant et avec (1) on a
ν
r̈ = − √ (5)
Ri rri

C’est une force en 1/ r. On reformule ṙ
√ s r
2 ν r
ṙ = − √ 1− (6)
Ri ri

En utilisant l’aide on peut intégrer pour trouver (en posant x = r/ri )

3 x3/2 9ν 2
1− x− = t . (7)
4 4 16Ri ri2

La solution est implicite mais c’est déjà bien.

2) Aux temps courts on peut poser

x=1− (8)
et on trouve
νt2 νt2
' , ⇒ r = ri − . (9)
2Ri ri2 2Ri ri
C’est une parabole, donc le mouvement d’une chute sous une force constante. En com-
parant avec (5) on voit que ça revient à approximer la force par la force initiale évaluée
en r = ri . Aux temps longs on va s’effondrer en boule et on peut estimer le temps pour
avoir x = 0 qui est p
tcollapse = 4/3 Ri /νri (10)
3) Si on met une rotation initiale Ωi = vi /ri où vi est la vitesse orthoradiale initiale sur
le grand cercle. L’énergie cinétique par volume est augmentée de (encore une fois c’est
approché mais justifié par le rapport d’aspect)
1
δEc /V = v 2 . (11)
2
Le moment cinétique par unité de masse est environ L/V ' rv. La force est centrale
donc conservation de L et il vient
rv = ri vi . (12)
En répétant l’analyse on a une barrière de potentiel comme dans le problème à deux
corps habituel. i.e. √
2 4ν r (vi ri )2
ṙ + √ + = Cte . (13)
Ri ri r2
La force totale effective par unité de masse est

ν (vi ri )2
− √ + (14)
Ri ri r r3

et elle s’annule si  2/5


vi2 Ri
rstable = ri . (15)
ν
Le potentiel effectif peut être lu sur (13). On peut le tracer pour bien voir la stabilité de
l’équilibre. On peut enfin si on a encore du temps discuter ce qui se passe si le système
n’est pas initialement lancé à rstable . On va avoir une série de contractions et dilatations
entre un rmin et un rmax . Comment raffiner ? En incluant la viscosité. Les oscillations
vont s’amortir sans modifier le moment cinétique et donc on va tendre vers le rayon de
stabilité étant donnée une rotation initiale. Pour le voir il faut discuter symmétrie du
problème pour voir que les forces de viscosité ne peuvent pas être orthoradiales et donc
modifier le moment cinétique.
Trou noir 1

Le 10 avril 2019 la collaboration Event Horizon Telescope a publié la première photo d’un trou
noir, celui au centre de la galaxie M87. Ce trou noir a une masse d’environ 6.5 × 109 MS où
MS ' 2×1030 kg est la masse du soleil. Mais la masse d’un trou noir n’est pas nécessairement si
élevé. Un trou noir se caractérise par une attreaction gravitationnelle si forte dans l’espace qui
l’entoure que rien, ni même la lumière ne peut s’en échapper, une fois à l’intérieur de “l’horizon”
du trou noir. Ce que l’on voit sur la photo est de la matière chaude à l’extérieur de l’horizon
du trou noir.

a) Rappeler comment on calcule la vitesse de libération (p.ex. pour une fusée) sur terre.

b) Dans le cadre de la mécanique classique, calculer la vitesse de libération pour un objet


astronomique (étoile, planète,...), supposé à symétrie sphérique, de masse m et rayon R. Quelles
conditions trouve-t-on alors sur m et R pour que la vitesse de libération soit supérieue à la
vitesse de la lumière? Que devient cette condition si m = MS et quelle conclusion sur la masse
volumique (à tort supposée uniforme) d’un tel trou noir ? Comparer avec différentes masses
volumiques que vous connaissez.

c) Toujours dans le cadre de la mécanique classique, et en traitant à tort la lumière comme des
corpuscules ayant une masse et allant à la vitesse de la lumière, comment nous apparaı̂trait
alors une région du ciel dans laquelle se trouverait un trou noir, en supposant qu’il n’y ait pas
de gaz ou autre matière à proximité.

d) Maintenant on suppose qu’à un instant initial un tel trou noir de rayon RS se trouve entouré
à l’extérieur d’une couche sphérique d’hélium de masse volumique ρ0 et température initiale
RS
T0 = 3 K uniformes, située autour de r = 100RS et d’une épaisseur h = 10
. Décrire ce qui se
passe et commenter si cela pourra ressembler à la photo.

1
a) On obtient l’énergie potentielle gravitationnelle en intégrant la force de Newton : Ep (r) =
− GMr m . On aura donc une énergie totale zéro à r = ∞ si m 2
2
v − GM m
r
= 0 et donc pour la
2 2GM
vitesse de libération v = r
. Tout ceci suppose bien sûr une masse volumique à symétrie
sphérique.
2GM
b) Pour v = c on a r = c2
. Pour une masse M donné ce rayon s’appelle le rayon de
Schwarzschild:
2GM
RS =
. (1)
c2
C’est aussi la valeur correcte du rayon de l’horizon d’un trou noir en relativité générale. Pour
une masse solaire on trouve RS = 3 km, largement inférieur au rayon du soleil. On appréciera
l’incroyable coincidence entre le nom de Karl Schwarschild et la signification littérale : schwarz
= noir, schild = bouclier !
On aura donc un trou noir si le rayon de l’objet céleste r est inférieur à RS , soit pour la
masse volumique
M M M c6 3c6
ρ= 4π 3 > 4π 3 = 4π = . (2)
3
r 3
RS 3
8G3 M 3 32πG3 M 2
On voit donc que la densité critique est d’autant plus grande que la masse est petite. Pour
M = Msoleil et RS = 3 km on trouve par exemple ρ > 2 × 1016 cmg 3 A titre de comparaison, pour
g
la matière nucléaire on a, en prenant un noyau m ' 10−23 g et r ' 10−15 m, soit ρ ' 1016 cm3

c) Sans rentrer dans les calculs, il est clair qu’il existent deux rayons, le rayon de capture Rc
et le rayon de Schwarzschild RS avec Rc > RS . En effet, le rayon de capture correspond au
rayon maximal jusqu’auquel un ”corpuscule de lumière” peut s’approcher sans être capturé et
ensuite spiraler vers le le centre. Clairement, pour cela il suffit d’une attraction plus faible que
pour contraindre un ”corpuscule de lumière” à faire demi-tour sur une trajectoire radiale. Donc
Rc > RS . Les trajectoires pour lesquelles la distance de plus grande approche est supérieure à
Rc subissent juste une déviation vers le centre. Si on considère l’ensemble des rayon lumineux
qui proviennent des étoiles derrière le trou noir, on a un fort effet de lentille (gravitationnelle).
En particulier on ne verra pas de rayons qui proviennent du disque apparent du trou noir, et
par l’effet de lentille ce disque paraitra plus grand que Rc . On verra donc bien une tache noire
dans le ciel.
ρ
d) Le gaz (parfait) a une pression “thermique” de p = R T
µ B
où ρ est toujours la masse
−3
volumique, µ = 4 × 10 kg la masse molaire et RB la constante des gaz. Sous les forces
gravitationnelles, ce disque va subir une force radiale (vers le centre du trou noir) par unité de
GM ρ(4πr2 h) GM ρh
surface de r2 ×4πr2
= r2
. Ceci est l’équivalent d’une pression

GM ρh ρ GM hµ RS h c2 µ
pgrav = 2
= RB Tgrav , Tgrav = = . (3)
r µ RB r2 r2 2RB

2
c2 µ RS h
Mais 2RB
' 2 × 1013 K tandis que r2
= 10−5 , soit Tgrav ' 108 K. Evidemment ceci n’est
pas la pression thermique dans le gaz, mais l’équivalent d’une pression radiale qui va faire
contracter cette couche de gaz vers des r de plus en plus petits. Ce calcul montre que les forces
gravitationnelles sont largement supérieures aux forces thermiques. La contraction résulte en
une compression de la couche de gaz qui va donc augmenter sa (vraie) température thermique.
En passant de r = 100RS à r = 2RS (et en supposant h constant), le volume sera divisé
par 2500 et la tempŕature augmentera par un facteur du même ordre pour atteindre quelques
milliers de Kelvin résultant en une émission de lumière. Si on prend 3 K × 2500, on atteint
7500 K, très proche de la température de la surface du soleil et donc l’émission de lumière avec
un spectre très proche. D’où la photo.

3
Cylindre et dipole magnétique

Un cylindre de métal d’axe vertical, de longueur L, rayon R et d’épaisseur e, avec


L  R  e peut se déplacer librement selon son axe. Il a une masse m et subit également
l’action de son poids. La conductivité du métal du cylindre est σ. On place un dipôle
magnétique au coeur de ce cylindre. Le dipôle est fixe et son moment est orienté vertica-
lement (M ~ = M~ez ). On cherche à déterminer les équations du mouvement du cylindre
et à les résoudre. Le champ magnétique créé par un dipôle magnétique est (voir schéma)
en coordonnées sphériques

B~ = Br~er + Bθ~eθ
µ0 M 2 cos θ
Br =
4π r3
µ0 M sin θ
Bθ =
4π r3
Calculer le courant induit dans un élément infinitésimal du cylindre défini par la varia-
tion infinitésimale dθ (voir schéma).

En déduire la force exercée sur cet élément, puis celle sur tout le cylindre

Etablir et résoudre les équations du mouvement (On note V la vitesse verticale du cy-
lindre).

Aide : 0 sin6 θ cos2 θdθ = 5π/128.
Eléments de solution

1) On doit d’abord calculer le flux à travers un disque situé à une hauteur z et corres-
pondant donc à une distance radiale r et un angle θ tels que
p
r = R2 + z 2 (1)
R
sin θ = (2)
r
z
cos θ = (3)
r
R
tan θ = . (4)
z
(5)

On peut bien évidemment calculer le flux à travers le disque de hauteur z par rapport au
dipôle, mais cela fait intervenir à la fois du Br et du Bθ car tous deux ont une composante
~ = 0, on peut calculer le flux à travers une calotte de sphère
selon ~ez . En utilisant divB
s’appuyant sur ce disque. Dans ce cas, seul Br peut contribuer au flux. L’élement de
surface en coordonnées sphériques adaptées à cette calotte est d2 S = r2 dθ sin θdφ. Il
vient donc (en intégrant trivialement sur φ)
Z Z
µ0 M
Φ = 2π Br (r)r2 sin θdθ = cos θ sin θdθ (6)
r
µ0 M µ0 M
= sin2 θ = sin3 θ . (7)
2r 2R
La loi de Lenz va donner le champ électromoteur dans un élément de hauteur dz. Cet
élement à une conductance
σedz
dG = (8)
2πR
et on pourra écrire
dI = −dGΦ̇ . (9)
Ce courant va faire subir une force de Laplace et donc il faut le calculer et donc calculer
Φ̇.
3µ0 M
Φ̇ = sin2 θ cos θθ̇ . (10)
2R
Il faut relier la vitesse du cylindre V à θ̇. En utilisant tan θ = R/z il vient

dθ sin2 θ
=− (11)
dz R
et donc
sin2 θ
θ̇ = − V. (12)
R
Ainsi la variation du flux à cause du mouvement du cylindre (ou mouvement du dipole
à cylindre fixe)
3µ0 M
Φ̇ = −V sin4 θ cos θ . (13)
2R2
On aurait pu arriver à ce résultat en intégrant les forces de Lorentz (~v ∧ B).~ Cela offre
~ ~
deux voies de calcul. Il faut juste utiliser ~v = V ~ez puis j = σ~v ∧ B. Si on choisit cette
voie il suffit de voir que,
~ez ∧ ~er = sin θ~eφ ~ez ∧ ~eθ = cos θ~eφ (14)
et utiliser r3 = R3 sin3 θ pour n’exprimer les résultats qu’en fonction de (R, θ). Au final,
peu importe la méthode, on a (en utilisant également (11))
1 3µ0 M
dI = σe sin2 θ cos θV dθ . (15)
2πR 2R
2) On calcule les forces de Laplace sur cette tranche de cylindre. On met les signes en
sachant qu’on va pas inventer le mouvement perpétuel et que l’on doit avoir un freinage,
et on trouve
~ · ~ez = −2πRdI(Br sin θ + Bθ cos θ)
dFz = 2πRI(~eφ ∧ B) (16)
 
9 µ0 M 2
= − σeV cos2 θ sin6 θdθ . (17)
8π R2
Pour un cylindre très long on intègre θ de 0 à π et la force totale est (en utilisant l’aide)
 
45 µ0 M 2
Fz = − σeVz . (18)
1024 R2
3) C’est une force de frottement qu’on peut remettre sous la forme
m
Fz = − Vz (19)
τ
où m est la masse totale du cylindre. On obtient alors
1
z̈ + ż − g = 0 . (20)
τ
Passé le régime transitoire on trouve la vitesse constante de chute qui est gτ . En pratique
c’est en gardant le cylindre fixe mais en faisant tomber le dipole magnétique que l’on
peut réaliser cette expérience de freinage par induction, car il sera beaucoup plus léger,
et subira la même force.
Pour aller plus loin, on peut faire un bilan énergétique car l’énergie cinétique restant
constante, la baisse de l’énergie potentielle doit se retrouver en dissipation dans les cou-
rants. En effet la puissance dissipée est
|dP | = |edI| = |Φ̇dI| = |dG(Φ̇)2 | (21)
et en utilisant les expressions précédentes, on voit que c’est égal en norme à V dFz et une
fois intégré sur tout le cylindre c’est aussi égal à |mgV | en régime permanent qui est la
puissance de la gravité.
Gravitation maximale

On considère un point dans l’espace P entouré par une substance maléable de volume V et
de masse volumique constante ρ. On cherche quelle forme doit prendre la substance pour
obtenir un champ gravitationnel dont le module est maximal en P .

1) Montrer que P est nécessairement sur la surface de la substance.

2) Montrer que si la forme est une sphère, alors, ce n’est pas la forme qui donne le module
du champ de gravitation maximal en P .

3) Trouver une équation pour cette forme. Comprarer la autant que possible à une shpère de
même volume.
C’est un exercice de réflexion. Il y a beaucoup d’idées possibles à essayer qui montreront la
connaissance du cours, même si l’exercice en tant que tel n’est pas trouvé complètement.

1) Question facile: Si P n’est pas sur la surface, on considère une petite sphère autour de
P (centre de cette sphère) complètement incluse dans le volume V . Cette sphère donne un
champ de gravité nul en P . Donc, on a perdu au moins le volume ce cette sphère pour le
champ de gravité de la substance de volume total V en P . Ce n’est donc pas le bon endroit
pour placer P : P doit être sur la surface.

2) Si la forme de la substance est une sphère de cenre O, de rayon R et P est sur la surface,
le module du champ de gravité en P est immédiat: g(P ) = GM/R2 (avec M = 43 πR3 ρ) et
ce champ est radial. Maintenant, imaginons que l’on prenne un tout petit volume de masse
dm au point diamétralement opposé à P et que l’on déplace dm/2 de 900 et dm/2 de −900 .
Avant déplacement, dm contribuait au module du champ comme:
Gdm
gini (P ) =
4R2
Après déplacement, on le champ en P est toujours radial mais sa valeur est diminuée de gini
et augmentée de:
Gdm
gpost (P ) = √ .
2 2R2
Donc, le module du champ en P a augmenté de:
Gdm Gdm
√ − >0
2 2R 2 4R2
Le module champ de gravitation en P est donc plus grand après modification qu’avant, donc
la sphère ne peut pas être la forme qui donne le module du champ maximal en P .

3) Une fois le raisonnement du 2) acquis, le plus dur est fait. Ce que l’on comprend, c’est
que si jamais il y a un endroit sur la surface de la forme avec un petit élément dm que l’on
peut déplacer ailleurs sur la surface et qui donnerait un module du champ de gravitation plus
grand en P , alors c’est ce qu’il faut faire. Au final, on va aboutir à une configuration où tous
les petits éléments dm sur la surface vont donner le même module du champ en P (sinon, on
continue les déplacements). Voyons ce que cela donne techniquement: on repère un point M
de la surface (masse dm) par sa distance r par rapport à P et son angle θ par rapport à l’axe
du champ de gravitation maximal (en module) en P . Le raisonnement précédent donne:

Gdm
cos(θ) = const.
r2
L’équation de la surface de la forme est donc:

r2 = C cos(θ).

C > 0. Arriver jusque là, c’est déjà très bien. Pour voir un peu mieux à quoi cela correspond
(en particulier par rapport à une sphère), il faut écrire cette équation différemment. On prend
P comme origine des corrdonnées, l’axe x comme axe du champ de gravitation maximal (en
module). Il y a nécessairement une symétrie de rotation par rapport à (x) de la forme. Donc,
on a juste besoin de repérer un point M de la forme par un couple (x, y) et:
x
x2 + y 2 = C p .
x2 + y 2

Soit:
y 2 = C 2/3 x2/3 − x2 .

Rappelons que le point (0, 0) est P et y 2 > 0 tant que x < C. Donc, la forme commence

en (0, 0), ensuite y > 0 puis y redevient nul pour x = C. On peut d’ailleur voir aussi que

dy/dx = ∞ pour x = 0 et x = C. La forme est donc assez proche d’une sphère de diamètre

C et, au vu du raisonnement du 2), on se doute que c’est une forme proche d’une sphère
que l’on aurait un peu contractée suivant (x). Pour comprendre un peu mieux cette idée, le
plus simple est de calculer le volume de la forme que l’on sait égal à V :
Z √
C
4π 3/2
V = πy 2 dx = C .
0 15
Ensuite, il est facile de voir qu’une sphère de même volume V aurait un un diamètre un peu
plus grand, ce qui confirme le point ci-dessus.
MQ - Oscillateur anharmonique

On rappelle qu’en mécanique quantique, pour un oscillateur harmonique à une dimension, c’est
m 2 2
à dire dans un potentiel V (x) = 2
ω x, l’énergie de l’état fondamental et la fonction d’onde
de l’état stationnaire correspondant sont
~ω mω 2 
E0 = , ϕ0 (x) = A exp − x . (1)
2 2~

a) Discuter la signification de ϕ0 et déterminer la constante A.

b) On change maintenant le potentiel V en y ajoutant une “petite” modification

m2 ω 3 4
δV (x) = γ x . (2)
2~
Justifier cette façon d’écrire δV et préciser ce que veut dire “petit”. Comment change l’énergie
de l’état fondamental en première approximation?

1
a) La fonction d’onde ϕ0 est solution de l’équation de Schrödinger stationnaire:
~2 00 m ~ω
− ϕ0 (x) + ω 2 x2 ϕ0 (x) = ϕ0 (x) . (3)
2m 2 2
Son module au carré |ϕ0 (x)|2 représente une densité de probabilité de présence (soit dx |ϕ0 (x)|2
la probabilité de trouver la particule entre x et x + dx), d’où la condition de normalisation
Z ∞ Z ∞ r Z ∞ r
2 2 mω 2  2 ~ −ξ 2 2 ~π
1= dx |ϕ0 (x)| = |A| dx exp − x = |A| dξe = |A| , (4)
−∞ −∞ ~ mω −∞ mω
et donc, à un choix de phase irrelevant près.
 mω 1/4
A= . (5)

On notera le changement de variables
r

ξ= x (6)
~
effectué pour calculer l’intégrale. La variable ξ est sans unité. Le potentiel peut alors s’écrire
~ω 2
V (x) = 2
ξ .

b) Le terme supplémentaire δV est de même δV (x) = γ ~ω


2
ξ 4 . On voit donc que γ est un
paramètre numérique (sans unité) et la condition δV petit devant V sera alors exprimé par
γ  1.

Dans la suite il ne s’agit pas de résoudre Schrödinger avec V + δV mais d’amener le candidat
à raisonner sur l’interprétation de la probabilité de présence.

En x, l’énergie potentielle de la particule est augmenté de δV et en pondérant cette augmenta-


tion par la probabilité de présence on obtient
Z ∞
δE ' dx |ϕ0 (x)|2 δV (x) . (7)
−∞

La suite est le calcul de l’intégrale:


r Z ∞ Z ∞
~ −ξ 2 ~ω 4 ~ω 1 2
δE = |A| 2
dξe γ ξ = γ √ dξe−ξ ξ 4 . (8)
mω −∞ 2 2 π −∞
R∞ 2 ∂2
R∞ 2 ∂2
pπ p
Pour l’intégrale on a −∞ dξe−αξ ξ 4 = ∂α 2 −∞
dξe−αξ = ∂α 2 α
= 43 απ5 , puis en posant
R∞ √
α = 1 on a −∞ dξe−ξ ξ 4 = 43 π. Finalement δE = γ ~ω
2
2
× 43 et
~ω  3 
E0 = 1+ γ . (9)
2 4

2
Corde secouée

Une corde inextensible de masse linéique ρ et de longueur L est suspendue par l’une de ses
extrémités et libre à l’autre bout. Elle pend verticalement (axe (z)). On s’intéresse aux petits
mouvements suivant la perpendiculaire à cet axe lorsque l’on secoue (faiblement) la corde
suivant un axe horizontal (x).

Dans ces conditions, quelle est la plus petite fréquence pour une onde stationnaire qui peut
exister sur cette corde?

Indications:
On pourra étudier un développement en série entière des solutions stationnaires de l’équation
obtenue pour des petites oscillations suivant (x) de la corde, écrite avec des variables adimen-
sionnées.

Si besoin, on admettra que l’équation:

1 − x + x2 /2!2 − x3 /3!2 + x4 /4!2 − x5 /5!2 + ... = 0

a comme plus petite racine: xmin ' 3/2.


Le début de l’exercice repose sur le cours sur les cordes vibrantes, avec une variation qui est
que l’on considère la corde vibrante comme suspendue. La fin de l’exercice est plus subtile.

La tension le long de la corde (en z) est: T (z) = ρg(L − z), c’est le poids de la corde suituée
sous z. Alors, l’équation des petites oscillations suivant (x) s’écrit:

∂ 2x ∂ ∂x
2
= g [(L − z) ].
∂t ∂z ∂z
Soit:
∂ 2x ∂ 2x ∂x
2
= g(L − z) 2
−g .
∂t ∂z ∂z
On ne s’intéresse qu’aux ondes stationnaires, donc ici, de la forme: x(z, t) = cos(ωt)f (z), ce
qui donne:
z ∂ 2f ∂f
−ω 2 f = gL(1 − ) 2 −g .
L ∂z ∂z
L’énoncé indique de se placer en variable adimensionnée. On multiplie alors cette équation
par L, puis on introduit: ξ = z/L. L’équation devient:

∂ 2f ∂f
−(ω 2 L/g)f = (1 − ξ) 2
− .
∂ξ ∂ξ
Soit, avec Λ = ω 2 L/g:
∂ 2f ∂f
(1 − ξ) 2
− + Λf = 0.
∂ξ ∂ξ
C’est l’équation pour les petites oscillations horizontales (stationnaires) de la corde. Les
conditions aux limites sont que la corde est fixée en z = 0, soit f (ξ = 0) = 0 et que les
mouvements sont d’amplitudes finies en z = L, soit f (ξ = 1) < ∞. C’est tout ce que l’on
sait.
Idéalement, on cherche à résoudre cette équation, puis appliquer les conditions aux limites,
ce qui doit donner les solutions en ondes stationnaires avec des conditions sur les fréquences
possibles, et donc que Λ ne peut prendre que quelques valeurs dénombrables. Pour répondre
à la question de l’énoncé, on doit alors trouver la fréquence la plus petite. Pour aller dans
cette direction, il faut suivre l’indication de l’énoncé et chercher un développement en série
P
entière pour f (ξ). Naturellement, on a tendance à écrire: f (ξ) = n≥0 an ξ n (avec a0 = 0),
puis l’équation pour f (.) donne une relation entre an+2 , an+1 et an . Alors on est bloqué car
on ne pourra jamais aller dans la direction voulue avec cela.

P
Arrivé à ce point, il faudra fournir une indication. Il ne faut pas poser: f (ξ) = n≥0 an ξ n
P
mais f (ξ) = n≥0 an (1 − ξ)n , avec f (ξ = 1) fini et
X
f (0) = an = 0.
n≥0

Dans ce cas, la relation de récurrence issue de l’équation pour f (.) devient:

(n + 1)nan+1 + (n + 1)an+1 + Λan = 0.

Que l’on peut facilement résoudre: an = −(Λ/n2 )an−1 , soit: an = (−Λ)n /(n!2 )a0 . Ensuite,
l’équation f (0) = 0 devient:
X (−Λ)n
= 0.
n≥0
n!2

En suivant l’indication donnée dans l’énoncé, on trouve alors la plus petite valeur de Λ
2
possible solution de cette relation, soit: Λmin = 1.5 = ωmin L/g. On obtient alors la
plus petite fréquence possible pour une onde stationnaire dans les conditions du montage:
1
p
νmin = 2π 1.5g/L, ce qui conclut.
Axions en présence de champ magnétique intense.

~ =
On considère une onde EM, plane progressive et harmonique, polarisée linéairement (E
E(x, t)~ez ) qui se propage dans le vide suivant l’axe (x) dans le sens des x croissants. Elle
arrive dans un domaine (entre x = −L et x = 0) où règne un champ magnétique constant
~ 0 = B0~ez (plus grand que 1010 T, que l’on peut rencontrer proche de certaines
très intense B
étoiles à neutrons). On admet que dans ces conditions, il peut y avoir un couplage entre
l’onde EM et B ~ 0 , ce qui va engendrer une densité volumique de courant: ~j = − η B~0 ∂φ , où
µ0 c ∂t
φ = φ(x, t) = φ0 ei(ωt−kx) est le champ qui est non nul du fait du couplage et η un paramètre
qui traduit l’intensité du couplage. Il vérifie l’équation de propagation modifiée:

1 ∂ 2φ ∂ 2φ η ~ ~
2 2
− 2 =− E.B0 . (1)
c ∂t ∂x µ0 c
~ le champ électrique de l’onde EM. En particulier, en dehors du domaine où B0 est
Avec E
non nul, φ(x, t) vérifie une équation de propagation libre à la vitesse c. De plus φ(x, t) et
∂φ(x, t)/∂x sont continus en x = −L et x = 0.

1 ∂2E~ ~
∂2E
1) Entre x = −L et x = 0, montrer que: c2 ∂t2
− ∂x2
n’est pas nul. En utilisant également
l’équation (1), en déduire la relation de dispersion dans la région de champ B0 non nul.

2) On considère que le champ φ (dit champ d’Axions) une fois produit dans le domaine
−L ≤ x ≤ 0 s’échappe de cette région, dans le sens des x croissants pour x > 0 et dans le
sens des x décroissants pour x < −L. Montrer que |φ| est proportionnel à |E|B0 pour x > 0
et x < −L.
(L’intérêt est que si |φ| est assez grand, alors il pourrait donner des effets observables).

Indication: on rappelle que la solution de l’équation (1) est une solution particulière plus une
solution de l’équation homogène (sans second membre).
C’est un exercice qui propose une variation sur le cours des équations de Maxwell, très cer-
tainement pas faite en classe. La question 1) est proche du cours, une fois que l’on sait écrire
les équations de Maxwell. Pour la question 2), il faut avoir bien compris l’énoncé et aussi
ne pas se tromper dans l’utilisation de la notation complexe pour les ondes qui se propagent
dans un sens ou un autre. Ici, on utilise la convention e−ikx pour une propagation dans le
sens des x croissants. Dans son ensemble, la question 2) est la variation par rapport au cours
que propose l’exercice. Attention, je ne demande pas de résoudre le système à la fin de la
question.

1) Dans le domaine où B0 est non nul, en utilisant la densité volumique de courant donnée dans
l’énoncé et en utilisant le fait que le champ électrique de l’onde EM et le champ magnétique
B~ 0 sont tous les deux suivant (z), on trouve facilement que:

1 ∂ 2E ∂ 2E η
2 2
− 2
= B0 φ̈.
c ∂t ∂x c
On voit donc un premier point très important: si φ(x, t) est proportionnel à ei(ωt−kx) , donc φ̈
aussi, alors E est également de cette forme. On en déduit donc que:
ω2 η
(k 2 − 2
)E = B0 (−ω 2 )φ.
c c
De même, l’équation de l’énoncé pour φ(.) donne:
ω2 η
(k 2 − 2
)φ = − EB0 .
c µ0 c
En combinant ces deux équations, on obtient la relation de dispersion, c’est-à-dire la relation
en k et ω:
ω2 2 η2 2 2
(k 2 − ) = B ω .
c2 µ 0 c2 0

2) On considère ici que le champ φ(.) a été produit dans la région de champ B0 intense. Il
vérifie dans cette région les relations de la question 1. Il s’échappe alors dans le sens des x
décroissants pour x < −L, soit: φ(x < −L, t) = c1 ei(ωt+k0 x) . La pulsation ω est la même
que pour la question 1 car on doit avoir continuité de φ(.) en x = −L et ceci quelque soit
t. Par contre, le vecteur d’onde k0 est différent car pour x < −L, φ(.) vérifie une équation
de propagation (standard) à la vitesse c, et donc: k0 = ω/c. Egalement le champ φ(.)
s’échappe dans la région des x croissants pour x > 0, soit: φ(x > 0, t) = c2 ei(ωt−k0 x) . Dans les
expressions des champs φ(.) ci-dessus, c1 et c2 sont des constantes (complexes) à déterminer
avec les conditions aux limites en x = −L et x = 0. Mais avant cela, il faut trouver φ(.)
dans la région de champ B0 intense, en particulier en fonction de E et B0 . Il faut résoudre
l’équation (1). La question 1 fournit la solution particulière φp (.):
η 1
φp (x, t) = − ω2
EB0 ei(ωt−kx) .
µ0 c (k 2 − c2
)
Avec k solution de la relation de dispersion de la question 1. La solution sans second membre
φh (.) est simple mais il faut faire attention car on peut avoir une propagation dans le sens des
x croissants ou décroissants, soit:

φh (x, t) = c3 ei(ωt−k0 x) + c4 ei(ωt+k0 x) .

Avec c3 et c4 des constantes (complexes) et k0 = ω/c. La solution générale pour φ(.) dans le
domaine −L < x < 0 est donc:
η 1
φ(x, t) = − ω2
EB0 ei(ωt−kx) + c3 ei(ωt−k0 x) + c4 ei(ωt+k0 x)
µ0 c (k 2 − c2
)

Maitenant, il faut utiliser les relations de continuité en x = −L et x = 0 pour φ(x, t) et


∂φ(x, t)/∂x pour déterminer les 4 constantes. C’est bien car il y a 4 équations. Le problème
est complètement posé. L’énoncé ne demande pas de résoudre ce système mais seulement de
montrer que |c1 | et |c2 | sont proportionnels à |E|B0 . C’est évident une fois le système posé
car les égalités doivent être vraies pour tout |E| et B0 , cela veut dire que l’on doit pouvoir
simplifier par |E| et B0 chaque égalité, ce qui implique que |c1 | et |c2 | sont bien proportionnels
à |E|B0 , ce qui conclut.

Note: on remarque qu’il y a une petite différence par rapport aux exercices de MQ où l’on applique
des relations de continuité avec une variable de plus que le nombre d’équations, cette variable étant
liée au flux incident (libre). Ici, ce n’est pas le cas car φ(.) est produit dans la région de champ
intense et donc on n’envoie pas un flux (libre) de φ(.) vers quelque chose.
Bille sur escalier

On considère un grand escalier descendant, dont les marches numérotées de n = 1 à n = N


(depuis le haut) ont toutes la même hauteur b et des largeurs différentes ξn a avec ξ1 = 1. Une
bille en acier glisse sans frottement sur le plancher avec une vitesse v0 vers l’escalier et tombe
sur la première marche de telle façon qu’elle rebondit exactement au milieu de la marche. Puis
elle rebondit sur la deuxième, puis la troisième marche, toujours en leurs milieux, et ainsi de
suite. Déterminer les valeurs des ξn puis exprimer le temps qu’il faut à la bille pour atteindre
la dernière marche si N = 100. On évaluera ce temps avec une précision d’au moins 1% si
a = b = 10 cm.

1
Un exercice de chute libre, donc a priori facile ...
La bille va décrire un bout d’une parabole. Prenons des coordonnées x pour l’horizontale et z
pour la verticale, dirigé vers le bas, et x = z = 0 au bord du plancher, soit au point où la bille
commence sa chute libre. En ce point la tangente à la parabole est horizontale et comme la
parabole doit arriver à x = a/2 et z = b, on a forcément pour cette première parabole
 2
z 2x
n=1 : = . (1)
b a
Ensuite, le rebond se fait comme une réflexion sur un miroir. La bille décrit donc l’image
miroir par rapport à x = a/2 de sa précédente trajectoire, et remonte jusqu’à z = 0 avant de
recommencer une nouvelle descente. L’équation de cette parabole est
 2
z 2(x − a)
n=2 : = . (2)
b a
√ √
Elle atteint alors la deuxième marche en z = 2b pour 2 = 2(x2a−a) , soit x2 = a(1 + 2
2
). Mais

a est la largeur de la première marche et a 22 doit être la moitié de la largeur de la deuxième
√ √
marche. Celle-ci a donc une largeur 2a, soit ξ2 = 2. Puis la bille rebondit et décrit sa

troisième parabole qui culmine au dessus du bout de la deuxième marche en x = a(1 + 2).
Donc √ !2
z 2(x − a − 2a)
n=3 : = . (3)
b a
√ √
3
Elle touche la troisième marche en z = 3 en x3 = a(1 + 2 + 2
) et la troisième marche doit
√ √
avoir une largeur 3a, soit ξ3 = 3. Et ainsi de suite. Donc

ξn = n , (4)

et la largeur totale de l’escalier est


N
X Z N
√ √ 2a
Ltot = a n'a dn n = N 3/2 , (5)
n=1 0 3
√ 1
où l’erreur et au plus de l’ordre de N , soit une erreur relative de N
ce qui est 1% pour
N = 100.
Dans tout le mouvement la vitesse selon x reste inchangée et égale à v0 . Le temps pour
parcourir l’escalier est donc simplement
Ltot 2a
T = ' N 3/2 . (6)
v0 3v0
g 2 g
Il reste à déterminer v0 . Pour la première parabole on a x = v0 t et z = 2
t = 2v02
x2 . On
2
q
identifie 2vg 2 = a4b2 , soit v02 = ga
8b
ou a
v0
= 8b
g
. Donc
0
s
4 2b 3/2
T ' N ' 190 s . (7)
3 g

2
Ondes EM en champs intenses

On considère 3 ondes EM (issues de 3 sources différentes) toutes polarisées linéairement suivant


(x). Les champs électriques de ces 3 ondes sont de la forme:

~ i (z, t) = Ai (z)ei(ωi t−ki z)~ex + c.c.


E

avec Ai (z) pratiquement constant sur des échelles de quelques longueurs d’onde. Elles ne sont
pas de même fréquence. En principe, il ne devrait pas y avoir d’interaction entre ces 3 ondes.
On va admettre ici que si les champs électriques sont suffisamment intenses, alors, il faut
modifier les équations de Maxwell en rajoutant un vecteur densité de courant de la forme:
3 ~ 3
~j = γ ∂|E| ~ex ,
c2 ∂t
où γ est une constante (dimensionnée). On va chercher à montrer que ces 3 ondes peuvent
alors effectivement interagir, et ceci en observant qu’une 4è onde peut apparaitre dans la zone
où se rencontrent les 3 ondes.

En écrivant le champ électrique de la 4è onde (que l’on suppose également polarisée suivant
(x)) sous la forme:
~ 4 (z, t) = A4 (z)ei(ω4 t−k4 z)~ex + c.c.
E

avec k4 = ω4 /c et où A4 (z) est pratiquement constant sur des échelles de quelques λ4 , montrer
que dA4 /dz n’est pas nul.
On supposera que ω1 + ω2 = ω3 + ω4 et ~k1 + ~k2 = ~k3 + ~k4 .
C’est un exercice qui décrit un phénomène physique qui peut surprendre.

On introduit le vecteur densité de courant de l’énoncé dans Maxwell-Ampère, puis avec


Maxwell-Faraday, on obtient immédiatement:

γ3 2 3 1
∂z2 E = µ0 2
∂t E + 2 ∂t2 E. (1)
c c
Avec ici, E = E(z, t) qui est la somme de tous les champs électriques en jeu, soit:
i=4
X i=4
X
E(z, t) = Ei (z, t) = Ai (z)ei(ωi t−ki z) + c.c. (2)
i=1 i=1

Dans l’équation (1), il y a beaucoup de termes mais on ne va garder que les terme en eiω4 t .
On obtient alors:
1 2 γ3
∂z2 [A4 (z)ei(−k4 z) ] + ω A 4 (z) = µ 0 (−ω42 )6(A1 (z)A2 (z)Ā3 (z)). (3)
c2 4 c2
En utilisant l’équation de conservation donnée dans l’énoncé et où Ā3 est le complexe conjugué
de A3 ... On n’a pas précisé dans l’énoncé que les Ai étaient nécessairement réels, mais on peut
le supposer: cela n’a pas d’importance ici. Ensuite, l’énoncé indique que d2 A4 /dz 2 donne un
terme négligeable dans l’équation (3) et c’est fini:

dA4 (z)/dz = −iγ 3 k4 6A1 (z)A2 (z)Ā3 (z).


Corde suspendue

Une corde inextensible de masse linéique µ et longueur L = 10 m est suspendue par l’une de ses
extrémités. On impose à ce point d’attache un mouvement horizontal x(t) = x0 cos ωt avec x0 = 5 cm
avec ω = 10 s−1 . Déterminer les positions x(t, z) de la corde à la hauteur z.

1
On prend z dirigé vers le bas avec l’attache en z = 0. L’approximation des petits x consiste a
2
supposer que dx
dz
 1, en particulier l’angle α(z) que fait la corde avec la verticale à la cote z est
dx(z)
petit et sin α(z) ' tan α(z) = dz
. On pourra aussi confondre zmax avec la longueur L. La tension
de la corde à z est donnée par le poids de la corde “en dessous”, soit T (z) = µ(L − z)g. On regarde
un élément de la corde entre z et z + dz. En z la tension est −T (z)(sin α(z), cos α(z)) et à z + dz
elle est T (z + dz)(sin α(z + dz), cos α(z + dz)). On projette sur ~ex et on a le PFD
d2 x d 
dzµ 2 = T (z + dz) sin α(z + dz) − T (z) sin α(z) = dz T (z) sin α(z)
dt dz
dx(z) d2 x(z)
= −µg dz + µg(L − z) dz (1)
dz dz 2
soit
d2 x dx d2 x
= −g + g(L − z) . (2)
dt2 dz dz 2
Pour des petits z (z  L) on pourra remplacer (L − z) par L et chercher une solution de la forme
 
x(t, z) = x0 Re eiωt−ikz . Pour zq
= 0 on a bien le mouvement imposé x0 cos ωt. L’éq (2) donne alors
i ω2
2 2
ω = igk + gLk , soit k = 2L ± gL − 4L1 2 , soit
" s #
z ω2 1
x(t, z) ' x0 e 2L cos ωt − − 2 z tant que L − z ' L . (3)
gL 4L

Le choix du signe devant la racine indique que les ondes se propagent du point d’attache vers le bout
libre, et pas l’inverse. (Cependant, en réalité, il y a une sorte de réflexion au bout libre et on aura bien
une superposition des deux solutions avec les deux signes.) L’amplitude des oscillations augmente
exponentiellement avec z,qmais cette solution n’est approximativement valable que si (L − z) ' L,
ω2
disons jusqu’à z ' L4 . Ici gL − 4L1 2 ' 1 m−1 et les variations de l’argument du cos avec z sont donc
bien plus rapides que la variation de la tension ∼ (L − z). On obtient alors une solution approchée
en remplaçant simplement L par L − z:
" s #
z ω2 1
x(t, z) ' x0 e 2(L−z) cos ωt − − z tant que L − z ≥ 1 m−1 . (4)
g(L − z) 4(L − z)2

On observe toujours l’amplification exponentielle de l’amplitude lorsque z augmente. Au voisinage


de z ' L on pose ξ = L − z et (2) devient
d2 x dx d2 x
= g + gξ . (5)
dt2 dξ dξ 2
En fait, (5) est strictement équivalent à (2). L’ansatz x(t, ξ) = eiωt f (ξ) conduit à
ω2
ξf 00 + f 0 + f =0, (6)
g
ce qui est une équation différentielle connue (éq de Bessel). La solution exacte est donc donnée
en termes de la bonne combinaison linéaire des deux solutions de (6) (fonctions de Bessel) f (ξ) =
c1 f1 (ξ) + c2 f2 (ξ). Les constantes ci sont déterminées telles que f (L) = x0 et par une deuxieme
dx
condition au bout libre, en ξ = 0, à savoir dξ
= 0.

2
Raidissement d’une vague

On considère une vague qui se propage à la surface d’un fluide parfait et incompressible. On
admet que les variations de hauteur de fluide (c’est-à-dire le mouvement de la vague) peuvent
être décrites par l’équation:
∂δh ∂δh
+ [A + Bδh] = 0. (1)
∂t ∂x
Avec A e B des constantes positives. C’est une équation de propagation avec une vitesse de
propagation U (δh) = A + Bδh qui dépend de δh.

1) Montrer que: δh(.) = F (x − U (δh)t) est solution de l’équation (1) avec F (.) une fonction
suffisamment régulière.

2) On note G la fonction qui est l’inverse de F , soit: G(.) = F −1 (.) et x − U (δh)t = G(δh).
Montrer que G0 (δh) = dG/d(δh) vérifie la même équation que δh (donc l’équation (1)). En
déduire que la vague va se raidir à un moment donné.
La question 1) et le début de la question 2) sont simples. Il s’agit de reproduire des caculs qui
sont dans le cours sur les ondes progressives, avec la petite variation qu’il faut faire attention
car U (δh) se dérive aussi. La fin de la question 2) demande d’avoir compris ce qu’est une
dérivée en suivant le mouvement (dérivée particulaire), ce qui demande de réfléchir un peu.

1) le résultat est immédiat.

2) On cherche à montrer dans un premier temps que:

∂G0 ∂G0
+ [A + Bδh] = 0.
∂t ∂x
∂G0
On doit donc calculer par exemple: ∂t
. Mais G0 (.) est une fonction de δh, donc:

∂G0 ∂G0 ∂δh


= .
∂t ∂δh ∂t
∂G0
Même principe pour ∂x
, ce qui donne le résultat. Notons que c’est également vrai pour G(.)
qui est aussi une fonction de δh mais pour G c’est évident car: G(δh) = x−U (δh)t. Ce que l’on
a montré, c’est que la dérivée temporelle de G0 en suivant le mouvement avec la vitesse U (δh)
est nulle, soit: dG0 /dt = 0. Alors, si l’on suit le mouvement: G0 (δh)[t = T ] = G0 (δh)[t = 0].
Il faut voir ce que cela signifie. Jusqu’ici, on a suivi l’énoncé avec des questions sur les
ondes progressives. Ici, l’énoncé indique qu’il faut utiliser l’invariance de G0 en suivant le
mouvement (relation ci-dessus). La meilleure stratégie est donc de voir ce que cela donne.
On a : G(δh) = x − U (δh)t. Donc, en particulier:

∂G ∂δh ∂δh
=1−B t = G0 (δh) .
∂x ∂x ∂x
Donc:
1 − B ∂δh t 1
G0 (δh) = ∂x
∂δh
= ∂δh
− Bt.
∂x ∂x

L’invariance de G0 dans le temps en suivant le mouvement donne:


1 1
∂δh
= ∂δh
− BT.
|
∂x t=0
|
∂x t=T
∂δh
Alors, pour une certaine valeur de x à l’instant initial (pour lequel |
∂x t=0
< 0), on peut
trouver un temps T tel que:
1
∂δh
= −BT.
|
∂x t=0
Et donc, à cette même position dans la vague, en ayant suivi son mouvement, on obtient:
1 ∂δh
∂δh
|
= 0, soit |
∂x t=T
= ∞, la vague se raidit, ce qui conclut.
∂x t=T
Echange d’un solvant

On considère une boı̂te séparée en deux par une membrane fixe semi-perméable (figure).

Le membrane sépare les 2 milieux représentés sur la figure: (1) étant un milieu dilué constitué
d’un solvant et d’un soluté et (2) constitué du même solvant seul. Seul le solvant peut traverser
la membrane. Cette dernière est d’épaisseur e et elle est trouée par des fentes de hauteur b
et de largeur L (dans le plan perpendiculaire à la figure, L est donc aussi la largeur de la
boite). Le solvant s’échange entre les 2 milieux par ces fentes. On suppose que P1 > P2 (la
pression dans (1) est la plus grande) et que le solvant est un fluide Newtonien de viscosité η.
On suppose de plus que le volume molaire du solvant et du soluté sont identiques.

Sans utiliser la notion de potentiel chimique (vu en cours de chimie), déterminer l’expression
de ∆P = P1 − P2 lorsqu’un régime permanent s’établit (on montrera en particulier que ∆P
est proportionnelle à η).

Note: si besoin, les bords extérieurs de la boı̂te peuvent être un peu mobiles
L’exercice combine 2 exercices classiques, probablement faits en classe: un problème de dif-
fusion particulaire (dans un cas simple) et un problème d’écoulement de Poiseuille plan. Ce
qui est compliqué ici c’est qu’il faut rapprocher ces 2 histoires. Donc, si c’est fait, c’est un
excellent point même si les calculs ne sont pas conduits très loin. Sinon, il y a beaucoup de
choses à dire qui relèvent du cours.

Donc, le point à comprendre, c’est que le régime permanent ici, c’est l’équilibre entre les
2 courants possibles: celui de diffusion particulaire du solvant entre (1) et (2) et celui de
l’écoulement de Poiseuille plan par chaque fente (∆P 6= 0). L’égalité (en module) entre les
vecteurs densité de courant de ces 2 phénomènes est la seule manière d’avoir un équilibre (ou
un régime permanent). Cela veut dire évidemment que nécessairement, dans les conditions
du montage, la densité particulaire de solvant n’est pas la même dans (1) et (2) -sans que
ce soit besoin de le préciser dans l’énoncé-. Egalement, il serait préférable que les vecteurs
densité de courant soient de sens opposés.

On va toujours raisonner pour une fente. On oriente (x) de (1) vers (2). On commence
par exemple par le courant de diffusion du solvant: Jx = −D(nsolv,2 − nsolv,1 )/e où nsolv,1
est le nombre de particules de solvant dans (1) (Nsolv,1 ) divisé par V1 . En introduisant, la
concentration particulaire de soluté (petite devant celle du solvant), puis l’égalité des volumes
molaires (énoncé), on arrive facilement à:

nsolv,2 /V2 nsolute


Jx = −D .
e nsolv,1

C’est un courant négatif, qui tend à pousser la membrane vers la gauche. Et en termes
d’unité, Jx multiplié par une surface est la variation par seconde d’un nombre de particules
(de solvant).

Ensuite, on passe au courant de Poiseuille plan. L’énoncé donne ∆P = P1 − P2 > 0, ce qui


va provoquer un courant de Poiseuille vers la droite. Donc, tout va bien. L’établissement du
profil des vitesses entre 2 fentes rectangulaires est fait en cours. Dans le montage de l’exercice,
cela donne:
P2 − P 1 2
v(y) = (y − b2 /4)
2ηe
R b/2
Le débit volumique correspondant est: Dv = −b/2 v(y)dyL et vecteur densité de courant de
Dv
Poiseuille (densité de particules de solvant de (1) multipliée par la vitesse moyenne Lb
):

nsolv,1 Dv nsolv,1 ∆P 2
Jx0 = = b > 0.
V1 Lb V1 12ηe

Il ne reste plus qu’à écrire: |Jx | = |Jx0 | et c’est fini. On trouve:

12ηD nsolute
∆P = .
b2 nsolv,1

Si on prend des valeurs numériques raisonnables: D = 10−10 SI, η = 0.01 SI et b ∼ 1 µm avec


un rapport soluté/solvant de 1/10, on trouve: ∆P = 1 atm, ce qui est compatible avec une
pression osmotique.
Corde qui tombe

On considère une corde de longueur l et de masse linéique µ suspendue par une de ses
extrémités à un plafond, l’autre bout étant libre. A l’instant initial, on amène le bout li-
bre de la corde au niveau du plafond, à coté de l’attache. La corde est alors pliée en 2 et
l’ensemble est au repos. On lache alors le bout libre et a corde va se déplier. On suppose que
ce mouvement s’effectue sans dissipation de chaleur.

Combien de temps va mettre la corde à se déplier complètement? Commenter.

On donne: s
Z 1
1−x
dx ' 1.2
0 x(2 − x)
On prend l’axe (z) vertical avec z = 0 au niveau du plafond et sens des z > 0 vers le bas.
Lorsque le bout libre de la corde est à la cote z, le pli de la corde est en l/2 + z/2. La partie
de la corde qui tombe à la vitesse ż (le mouvement de chaque point de la partie qui tombe
est le même) a donc pour longueur: l/2 + z/2 − z = l/2 − z/2. L’énergie cinétique de la corde
est donc:
1
Ec = µ(l − z)ż 2 .
4
C’est bien l’Ec de toute la corde puisque l’autre coté de la corde est au repos.

Pour l’énergie potentielle, il faut faire attention. On prend Ep (t = 0) = 0. Alors, lorsque le


bout libre est en z, Ep < 0 (puisqu’un opérateur devrait fournir de l’énergie pour remettre la
corde comme elle était au départ). La partie de la corde qui était entre z = 0 et z s’est répartie
des 2 cotés entre les cotes: z = l/2 + z/2 et z = l/2 (qui a pour barycentre: (l + z/2)/2).
Donc:
Ep = −µgz/2[(l + z/2)/2 − z/4] − −µgz/2[(l + z/2)/2 − 3z/4].

Soit:
µgz
Ep = − (2l − z).
4
Donc:
1 µgz
E = µ(l − z)ż 2 − (2l − z) = 0,
4 4
car 0 est la valeur de l’énergie à l’instant initial et E est conservée. On suppose qu’il n’y a
pas de chaleur dissipée. Ce qui donne:

z(2l − z)
z̈ = g + g > g.
(l − z)2

Z Z Z s s
l
dz l 2l
T = dt = dz/|ż| = q ' 1.2 < .
0 gz(2l−z) g g
l−z

Commentaire: le fait que z̈ soit plus grand que g vient du fait qu’il y a une tension T au
niveau du pli de la corde:
l−z l−z
µ z̈ = T + µ g.
2 2
Ce qui donne: T = ... Attention: T = ∞ en z = l car il n’y a plus conservation de l’énergie à
ce niveau. La tension de la corde au niveau de l’attache du plafond est:
l+z
T 0 = T + µg .
2
Force d’adhérence

On considère un disque en verre circulaire de rayon a posé à plat sur une plaque de verre
avec entre les deux une mince couche d’eau, considérée comme un fluide Newtonien (donc
incompressible) de masse volumique ρ, viscosité η et d’épaisseur h. On tire le disque vers le
haut avec une force F , donc h = h(t) et dh/dt > 0. On sait que c’est difficile en pratique et on
va comprendre pourquoi ici. On se place en coordonnées cylindriques (r, θ, z) avec (z) orienté
vers le haut et on admet que la force F provoque un écoulement dans la mince couche d’eau
avec le champ des vitesses: ~u = ur (r, z, t)~er + uz (r, z, t)~ez . On suppose toujours h(t) << a.
On suppose également que le nombre de Reynolds (horizontal) associé à cet écoulement est
tel que: ρ|u|a/η << a2 /h2 (dit écoulement rampant).

1) Justifier que dans les conditions du montage l’équation du mouvement radiale pour un petit
∂P 2
élément de fluide se réduit à: ∂r
= η ∂∂zu2r .
Qu’en est-il pour l’équation du mouvement suivant z?

2) Montrer que la pression dans la couche de fluide s’écrit:

6η dh
P (r) − P (a) = 3dt (r2 − a2 ).
h
dh
En déduire la force F pour a = 10 cm, dt
= 0.02 mm/s et h(0) = 0.02 mm (η = 10−3 Pa.s).

1 ∂(rur ) ∂uz
Formulaire: Pour ~u = ur (r, z, t)~er + uz (r, z, t)~ez , div(~u) = r ∂r
+ ∂z
.
C’est un exercice dans lequel on comprend qu’il faut négliger les termes instationnaire et
convectif par rapport au terme de viscosité. Cela est fait dans le cours dans un cas plus
simple qui utilise uniquement la comparaison du nombre de Reynolds avec l’unité. Ici, on
part tout de suite avec une variation de ce qui est fait en classe, mais le principe est le même.
De plus, on voit dans l’énoncé que le terme de viscosité ne comprend qu’une dérivée seconde
suivant z, on comprend donc qu’il faut utiliser la condition h << a pour ne conserver que les
termes dominants. Il y a beaucoup d’exercices qui utilisent cette technique mais le candidat
n’est pas tenu de les avoir faits. L’énoncé donne toutes les clefs dans la question 1) pour
avancer. La question 2) utilise uniquement ce qui est fait en classe, en particulier pour les
conditions aux limites.

1) Le début de l’exercice utilise ce qui est vu dans différents exercices utilisant les couches
limites ou le nombre de Reynolds, ici que le mouvement radial de l’écoulement s’effectue
sur une distance caractéristique a alors que le mouvement suivant z s’effectue sur h avec
h << a. Ceci va permettre d’estimer les différents termes qui entrent dans les équations avec
des dérivées suivant r ou z. Déjà, l’équation de continuité, ici div(~u) = 0 donne: uz /h ≡
ur /a et donc, en ordre de grandeur uz << ur . Ensuite, on voit facilement que ρ∂ur /∂t et
~
ρ(~v .grad)u 2 2 2
r sont du même ordre de grandeur, soit: ρur /a ' ρu /a mais ρu /a << ηu/h , ce
2

~
qui correspond à la condition de l’énoncé: ρua/η << a2 /h2 . Les termes ρ∂ur /∂t et ρ(~v .grad)ur

sont donc négligeables devant le terme de viscosité dans l’équation du mouvement radiale pour
un petit élément de fluide. En effet, dans le terme de viscosité, le Laplacien est complètement
dominé par ∂ 2 /∂z 2 ≡ 1/h2 . On en déduit que l’équation du mouvement radiale se réduit à:
∂P 2
∂r
= η ∂∂zu2r . Pour l’équation du mouvement suivant z, le mêmme type de raisonnement donne
une équation très simple: ∂P/∂z = 0, ce qui conclut.

2) Comme P (.) est indépendant de z, on peut intégrer l’équation du mouvement radiale pour
obtenir ur (.). On obtient:
1 ∂P 2
ur (r, z, t) = (z − zh(t)).
2η ∂r
En z = 0 et z = h(t), on a bien ur qui est nulle. Ensuite, la seule équation qui reste à utiliser
pour avancer, c’est l’équation de continuité. De plus, on voit dans le résultat donné dans
l’énoncé qu’il y a un dh/dt, soit uz en z = h(t), or la seule manière pour faire apparaitre
un terme avec uz , c’est d’utiliser l’équation de continuité. Il n’y a donc pas d’autre option
possible:
1 ∂(rur ) ∂uz
+ = 0.
r ∂r ∂z
Avec ur (r, z, t) trouvé plus haut. On doit donc injecter l’expression pour ur dans l’équation
de continuité mais il faut également intégrer cette équation sur z entre z = 0 et z = h(t) pour
faire apparaitre un terme en uz (z = h(t)) = dh/dt (et uz (z = 0) = 0). On écrit:
Z h(t) !
1 2 1 ∂(r ∂P
∂r
) dh
dz (z − zh(t)) + = 0.
0 2η r ∂r dt
Soit:
1 h(t)3 1 ∂(r ∂P
∂r
) dh
(− )( )+ = 0.
2η 6 r ∂r dt
Ce qui donne:
1 ∂(r ∂P
∂r
) 12η dh
dt
)=
r ∂r h(t)3
On peut alors intégrer 2 fois suivant r cette expression, en sachant qu’il n’y a pas de divergence
∂P
pour ∂r
en r = 0. On obtient facilement le résultat proposé dans l’énoncé:
3η dh
dt
(r2 − a2 ).
P (r) − P (a) =
h3
Avec P (a) la pression atmosphérique en dehors de la petite couche d’eau. On voit que
P (r) − P (a) < 0, donc la couche d’eau tend à attirer le disque en verre, c’est pour cela qu’il
faut tirer avec une force: Z a
F =− (P (r) − P (a))2πrdr.
0
On trouve:
3πη dhdt 4
F = a.
h3
On observe en particulier le terme en 1/h3 , c’est pour cela qu’il est très difficile de décoller
le disque. L’AN donne F ' 2.3 103 N, soit l’équivalent d’une masse de 230 kg suspendue
à une poulie. Il est également préférable de vérifier que la condition: ρ|u|a/η << a2 /h2 est
dh
satisfaite. avec |u| la vitesse horizontale qui s’obtient comme: uz h/a ' dt
h/a et on voit que
la condition est bien vérifiée: 104 << 108 /4.
Ondes surface

On considère un bassin infini rempli d’eau jusqu’à une hauteur h. On considère l’eau
comme un fluide parfait incompressible et le fond du bassin est horizontal.

— Lorsque l’eau du bassin est au repos, quelle est la pression dans l’eau ? Montrer
qu’en général, on a ∆P = 0.

— Lorsqu’il y a de petites ondes de surface, c’est-à-dire des petites vagues, quelle est
la pression dans l’eau ? Quelle est la vitesse de l’eau à la surface et sous la surface ?
Eléments de solution

C’est le champ de vitesse et les trajectoires des particules que cet exercice cherche à
déterminer.

Dériver les équations des ondes de surface quand on passe par le potentiel est un classique.
On rappelle la dérivation habituelle et on donne plus loin celle suggérée par le sujet. La
coordonnée verticale z a pour origine le fond du bassin. On cherche des déformations
(z = h + s) de surface de type
s = s0 ei(ωt−kx) . (1)
On linéarise donc on néglige l’accélération convective. Et donc le champ de vitesses
est irrotationnel (on trouve que c’est une constante et donc constante nulle pour des
~
oscillations) et dérive d’un potentiel ~v = ∇φ. L’incompressibilité donne ∆φ = 0 et on
cherche des solutions de la forme

Φ(x, z, t) = ϕ(z)ei(ωt−kx) (2)

ce qui est juste une séparation de variable où on devine ce qu’on va trouver. On a alors

ϕ00 (z) − k 2 ϕ = 0 . (3)

La solution est
ϕ = a1 ekz + a2 e−kz . (4)
On doit avoir une condition de vitesse verticale nulle au fond (or vz = ∂z φ) donc

a1 = a2 . (5)

La condition de surface est ∂t s = vz (en fait la condition générale est (24) mais on
linéarise) donc en ne gardant que l’ordre le plus bas il vient

k[a1 ekh − a2 e−kh ] = iωs0 (6)

donc
iωs0
a1 = a2 = (7)
2k sinh(kh)
et on trouve
iωs0 cosh(kz) i(ωt−kx)
Φ= e . (8)
k sinh(kh)
Il faut encore la relation ω(k). A partir d’Euler on trouve que

P
∂t Φ + + gs (9)
ρ
est une constante globale (on trouve que son gradient est nul donc plus précisément c’est
une fonction globale du temps). En prenant la dérivée convective selon les particules à
la surface, et en négligeant les ordres deux ça donne

∂tt Φ + g∂t s = ∂tt Φ + g∂z Φ = 0 (10)

qu’on évalue à z = h car c’est déjà d’ordre un, et alors il vient la relation de dispersion

ω 2 = gk tanh(kh) . (11)

On remplace dans Φ et en prenant le gradient on a le champ de vitesses.


cosh(kz)
vx = ωs0 cos(ωt − kx) , (12)
sinh(kh)
sinh(kz)
vz = −ωs0 sin(ωt − kx) . (13)
sinh(kh)

1) La méthode guidée implicitement passe directement par la pression. En effet à la


première question on demande de rappeler l’équilibre hydrostatique

Pstat = Patm − ρg(z − h) . (14)

De l’incompressibilité div~v = 0 et en négligeant le terme convectif dans l’équation d’Euler


il vient juste ∆P = 0. On vérifie au passage que c’est le cas pour la pression statique
ci-dessus.

2) Dans le cas avec des vagues, on va utiliser ∆P = 0. La pression va être légèrement


modifiée autour de la pression statique donc on pose

P = Pstat + δP ⇒ ∆δP = 0 . (15)

La réponse à une onde de surface (1) va être de la même forme et donc

δP = δP (z)ei(ωt−kx) , (16)

et la condition de Laplacien donne

δP 00 (z) = k 2 δP (z) . (17)

La solution est
δP (z) = b1 ekz + b2 e−kz . (18)
Avec l’équation d’Euler on obtient alors la vitesse
k
vx = δP (19)
ωρ
ik  kz 
vz = b1 e − b2 e−kz ei(ωt−kx) (20)
ωρ
La condition de vitesse verticale nulle au fond du bassin donne directement

b1 = b2 . (21)

La condition de pression atmosphérique sur la surface donne


 
b1 ek(h+s) + e−k(h+s) + Patm − ρgs0 = Patm . (22)

Dans les exponentielles on peut prendre h au lieu de h + s car c’est déjà un terme d’ordre
un et on a alors
ρgs0
b1 = . (23)
2 cosh(kh)
On a déjà utilisé l’équation de conservation qui donne l’absence de divergence de la
vitesse, et l’équation d’Euler. On a dit que la vitesse au fond est bien compatible avec la
forme du fond. Il reste à dire que la vitesse en surface est compatible avec la forme de la
surface. En général on trouverait

∂t s + vx ∂x s = vz , (24)

mais en linéarisant c’est seulement

∂t s = iωs = vz . (25)

En injectant dans (20) est en ne gardant que les termes d’ordre un on a alors

−ρω 2 s0 = −2kb1 sinh(kh) ⇒ ω 2 = kg tanh(kh) . (26)

On a donc ω(k) et en remplaçant dans (19) (20) la forme complète de la vitesse qui est
bien évidemment la même que dans la méthode basée sur le potentiel de vitesse (12).

Les trajectoire des particules sont normalement calculées avec

ẋ = vx (x, z, t) ż = vz (x, z, t) . (27)

c’est à dire qu’on a un passage du champ Eulerien de vitesse à de la formulation Lagran-


gienne où on suit une particule. Mais la différence est d’ordre deux car les déplacement,
tout comme la vitesse, sont d’ordre un. On peut donc prendre une intégrale sur le temps
très brutalement à partir du champ de vitesses pour avoir les trajectoires
cosh(kz)
x = s0 sin(ωt − kx) (28)
sinh(kh)
sinh(kz)
z = s0 cos(ωt − kx) (29)
sinh(kh)
A la surface, et si le bassin est très profond devant la longueur d’onde, les rapports des
fonctions hyperboliques sont 1 et on a des cercles de rayon s0 . C’est intéressant de voir
dans quel sens ça tourne en fonction du sens de propagation de l’onde, ici vers les x
croissants. On peut tracer le mouvement d’un canard à la surface. Quand il est au plus
haut il va dans le sens de l’onde et inversement quand il est dans le creux de la vague. Si
on ne surfe pas la vague on fait des ronds dans l’eau !

En revanche, tout au fond du bassin, si kz  1 on a des mouvements elliptiques et


l’amplitude du mouvement vertical est kz celle du mouvement horizontal.

Si maintenant le bassin est peu profond par rapport à la longueur d’onde on a


s0
x = sin(ωt − kx) (30)
kh
s0 z
z = cos(ωt − kx) (31)
h
et c’est partout que le mouvement vertical est en kz de celui du mouvement vertical. On
a toujours des ellipses.

En question subsidiaire on peut s’intéresser à la pression qui est

cosh(kz)
δP = ρgs . (32)
cosh(kh)

A la surface c’est donc comme si on avait une pression hydrostatique modulée par la sur-
épaisseur ou sous-épaisseur d’eau. Mais en profondeur on a une redistribution horizontale
de ces fluctuations de pression hydrostatique (c’est une interprétation qui a ses limites
puisque cela n’a pas de sens de parler d’hydrostatique ici) et cela tend donc à les effacer
une fois que kz  kh, donc lorsqu’on descend à des profondeurs plus grandes que la
longueur d’onde.
Disque de gaz dans l’espace

On supposera qu’une agglomération de gaz intergalactique puisse être décrite en bonne ap-
proximation, pendant un intervalle de temps suffisamment long, comme un disque de rayon
R ' 1011 m et de hauteur 2h ' 107 m, de masse volumique ρ constante et homogène. A t = 0
une fusée se trouve sur l’axe du disque à z0 = 108 m “au dessus” du centre avec une vitesse
nulle lorsque ses réacteurs tombent en panne.
a) Décrire le mouvement qu’effectuera la fusée. En particulier, calculer la période.
b) Montrer comment l’observation de la période permet à un astronome terrestre (qui ne connait
que la masse totale M et le rayon R du disque, ainsi que z0 ) de déterminer h.
c) Discuter les hypothèses faite sur la forme du nuage intergalactique.

1
Il s’agit d’abord de voir l’analogie avec le champ élecrique créé par une plaque uniformément
chargée. Donc, au-dessus et en dessous du disque le champ gravitstionnel est uniforme et
“vertical”. Ceci est valable tant la hauteur z satisfait |z|  R. Si on écrit pour la force
gravitationnelle sur une masse m en ~r : F~ (~r) = −m∇V
~ grav (~r) le potentiel gravitationnel Vgrav
satisfait à ∆Vgrav = 4πGρ, alors que ∆Vel = −ρel /0 . Le théorème de Gauss est donc
Z Z
~
∇Vgrav · dS = 4πG d3 xρ = 4πGM (vol) . (1)
S vol

~ grav est vertical ∼ ~ez et qu’il doit s’annuler en


Par la symétrie du problème, on sait que ∇V
z = 0. En choisissant la surface de Gauss comme un cylindre d’axe parallèle à l’axe du disque,
avec une base à z = 0 et l’autre en z on trouve facilement
 
z2 h
Vgrav (z) = 4πGρ |z| ≤ h , Vgrav (z) = 4πGρ h |z| − |z| > h (2)
2 2
L’énergie potentielle d’une masse m en z est donc Epot = mV (z). La conservation de l’énergie
de la fusée donne
1 2
ż + Vgrav (z) = Vgrav (z0 ) . (3)
2
Tant que z > h cela donne une chute libre avec accélération constante, soit ż 2 q
= 8πGρh(z0 − z)
z0 −h
ou z̈ = −4πGρh, soit z(t) = z0 − 2πGρht2 et la fusée atteint z = h à t1 = 2πGρh avec une
p
vitesse ż = − 8πGρh(z0 − h). Ensuite la fusée se trouve dans un potentiel harmonique avec
ω 2 = 4πGρ, soit z(t) = A sin ω(t − t1 ) + B cos ω(t − t1 ). Les conditions initiales donnent B = h
p
et ωA = ż(t1 ), soit A = − 2h(z0 − h):
p
z(t) = h cos ω(t − t1 ) − 2h(z0 − h) sin ω(t − t1 ) . (4)

La fusée arrive en z = 0 à t2 avec


s r r s !
h 1 2(z0 − h) h
tan ω(t2 − t1 ) = ⇒ t2 = + arctan . (5)
2(z0 − h) 4πGρ h 2(z0 − h)
Ensuite, le mouvement est symétrique pour z < 0 pour atteindre z = −z0 , puis la fusée remonte
jusqu’a z0 . Au total, le mouvement est périodique avec une période T = 4t2 . Avec les données
√ 2
numériques, h  z0 et on peut faire un DL ( 1 − x = 1− x2 − x8 +O(x3 ) et arctan x = x+O(x3 ))
r   r  
z0 h2 8z0 R2 h2
T '4 1+ 2 = 1+ 2 (6)
2πGρh 8z0 GM 8z0
où on a utilisé M = πR2 hρ. Si on connait M , R et z0 , l’observation précise de la période
h2
permet donc de déterminer h (effet en 8z02
' 3 × 10−4 ).
Les hypothèses faite sur la forme en disque avec un ρ constant sont hautement non-réalistes.
Selon le temps qui reste, on peut alors discuter comment ce disque, soumis à son propre champ
gravitationnel pourrait être stabilisé par la pression thermodynamique du gaz à température
non-nulle et par la rotation du disque autour de son axe.

2
Transfert thermique barreau

On considère une barre solide infinie et de section S, de capacité thermique c, de masse


volumique ρ et de conductivité thermique λ.

— Etablir l’équation de la chaleur. Introduire le coefficient de diffusion D.

T0 ` 2
— Vérifier que T (x, t) = √
2 πDt
e−x /(4Dt) est solution où ` est une longueur.

— A quelle profil de température initial cela correspond-il ? Quel est l’injection d’éner-
gie thermique faut-il réaliser à t = 0 pour obtenir ce profil initial ?

— On prépare maintenant le barreau à t = 0 avec la condition initiale T = T0 pour


−L < x < 0 et T = 0 partout ailleurs. Quelle est l’évolution ultérieure de la
température ?
√ √
— Discuter cette solution quand L  4Dt puis quand L  4Dt.

R∞ 2 R∞ −x2 dx
Aide. On posera si nécessaire la fonction E(x) = √1 e−z dz et on rappelle
π x −∞ e =

π.
Eléments de solution

Au centre du raisonnement est le principe de superposition autorisé par la linéarité de


l’équation.

~ et on trouve
1) On fait un bilan avec le flux ~j = −λ∇T

ρc∂t T = λ∆T = λ∂xx T (1)

et on pose D = λ/(ρc).

2) La vérification prend une ligne pour la dérivée temporelle et deux lignes pour la double
dérivée spatiale.

3) On a une Gaussienne qui s’élargit. A t = 0 c’est nul partout sauf en x = 0 où c’est


infini. On a donc injecté beaucoup d’énergie thermique en une zone infiniment petite. On
peut calculer cette chaleur qui est
Z
Q0 = ρcST (x)dx = ρcS`T0 . (2)

C’est comme si on avait injecté l’énergie nécessaire pour chauffer à T0 une tranche de
longueur ` du barreau, sauf que cette énergie a en fait été injectée dans un espace infini-
ment petit centré en x = 0.

4) Si on a bien interprété la question précédente, alors on applique le principle de super-


position. On considère que sur un élément de longueur dx on a apporté

δQ(x) = ρcST0 (x)dx (3)

et que ça diffuse ensuite selon la solution trouvée. La différence c’est que ces différentes
injections sont comprises entre −L et 0. On a une convolution de la solution précédente
(celle associée à un Dirac initial) avec la condition initiale en superposant et on trouve
Z 0
T0 2 /(4Dt)
T (x, t) = √ e−(x−y) dy . (4)
2 πDt −L

On peut reformuler
√ avec la fonction E proposée en faisant un changement de variable
Z = (x − y)/ 4Dt et on trouve
Z (L+x)/√4Dt
1 2
T (x, t) = T0 √ √ e−Z dZ (5)
π x/ 4Dt
    
x L+x
= T0 E √ −E √ (6)
4Dt 4Dt

5) Si L  4Dt alors la seconde contribution est négligeable. Donc on trouve
 
x
T (x, t) ' T0 E √ . (7)
4Dt
En x = 0, comme E(0) = 1/2 (trouvable avec l’aide), on a

T (0, t) = T0 /2 . (8)

Physiquement si L est très grand on a la même chose que lorsque l’on met deux barreaux
semi-infinis en contact. La température au point de contact est alors la température
moyenne (ici la moyenne entre 0 et T0 ). On peut calculer le flux (intégré sur la surface)
2 √
et on trouve en utilisant E 0 (x) = −e−x / π

T0
Φx (x = 0) = Sλ √ . (9)
2 πDt
On a donc retrouvé la solution d’un barreau semi infini avec un bord maintenu à une
température T0 /2. La partie x < 0 assure de toujours subvenir au flux nécessaire pour
maintenir T0 /2 à x = 0.
Si en revanche L est petit devant la longueur de diffusion (ou vu plus physiquement
lorsqu’on attend longtemps pour que la longueur de diffusion soit plus grande que L)
alors on peut faire un DL de (5) qu’on peut voir comme f (x) − f (x + L) et on trouve

T0 L −x2 /(4Dt)
T (x, t) ' √ e . (10)
2 πDt
On retrouve la première solution où l’énergie avait été concentrée sur une longueur infi-
nitésimale. En fait il vaut mieux faire ce DL autour de x + L/2 car le terme d’ordre deux
s’annule donc le résultat est plus précis. On trouve alors
T0 L −(x+L/2)2 /(4Dt) √
T (x, t) ' √ e + O((L/ Dt))3 ) . (11)
2 πDt
Là c’est plus clair. La diffusion a effacé l’étalement initial de la température et tout se
passe comme si on avait concentré l’injection au milieu du segment [−L, 0].
Vitesse de l’effet dominos

On considère une ligne de dominos alignés verticalement sur une table. Il sont séparés
les uns des autres par des espaces de longueur s, ont une largeur D et une hauteur h.
On suppose que la largeur est négligeable devant la hauteur (D  h), que les dominos
ne peuvent que basculer mais pas glisser sur la table. On cherche à étudier la dynamique
de l’avancement du front de renversement des dominos.

— Lorsque les dominos sont renversés les uns sur les autres, relier l’angle de bascu-
lement (par rapport à la verticale) d’un domino en fonction de celui de ses voisins.

— En supposant s  h, adopter une description continue de l’angle de basculement


et calculer sa forme θ(x), où x est la coordonnée selon la direction de l’alignement
des dominos, en fonction de la longueur totale ` de dominos ayant basculé.

— Calculer l’énergie mécanique des dominos en fonction de la vitesse v d’avancement


du front de renversement des dominos.

— En supposant l’énergie mécanique constante en déduire la vitesse asymptotique.


Discuter les hypothèses et trouver une amélioration de la description de l’effet do-
mino.

Aide : le moment d’inertie d’un solide rectangulaire (de côtés a et b et de masse m)


associé à une rotation autour d’un de ses coins est I = 31 m(a2 + b2 ).
Eléments de solution

1) Il existe au moins trois méthodes et toutes font appel à de la trigonométrie de base. Le


plus immédiat consiste à utiliser le fait que dans un triangle les rapports du sinus d’un
angle avec la longueur du côté opposé sont constants. Une autre méthode plus analytique
consiste à exprimer la position du point de contact de deux dominos, une fois pour un
domino et une fois pour l’autre. Une troisième méthode consiste à exprimer de deux
\ est droit. On a AC = s cos θn+1
manières différentes AC (voir figure) où l’angle ACB
mais aussi AC = h sin(θn − θn+1 ) donc

h sin(θn − θn+1 ) = s cos θn+1 . (1)

2) Le passage au continu peut se faire à partir de la relation précédente ou en réitérant


la méthode mais pour des infinitésimaux. θn devient θ(x). On obtient

dθ cos θ
=− . (2)
dx h
On voit que l’angle varie essentiellement sur une longueur de l’ordre de h. Bien à l’arrière
du front, les dominos restent figés dans une position horizontale, complètement tombés
les uns sur les autres. On va vérifier par la résolution. La résolution procède facilement
par le changement de variable u = sin θ (suggéré aux candidats) et en se souvenant que
arctanh0 (x) = 1/(1 − x2 ). On obtient en effet
x
arctanh(u) = − + Cte (3)
h
et la constante est déterminée par le fait qu’au front d’onde localisé à une distance `, les
dominos sont verticaux et θ(`) = 0. Il vient
1
θ = arcsin[tanh((` − x)/h)] , cos θ = (4)
[cosh((` − x)/h)]2
On peut tracer cos θ et on voit bien que pour une distance  h à l’arrière du front, c’est
très petit et donc θ ' π/2, c’est-à-dire les dominos sont complètement allongés.
P R`
3) La somme sur les dominos est remplacée par une intégrale ( i → 0 s−1 dx). L’énergie
potentielle est simplement
Z
mgh dx
Ep = (cos θ − 1) (5)
2 s
 Z 
mgh h ` dx 1 `
= − (6)
2 s 0 h cosh((` − x)/h) s

Et quand `  x le premier terme est négligeable. Essentiellement l’energie potentielle est


réduite par le fait que tous les dominos jusqu’à ` sont tombés sur le table.
Pour l’énergie cinétique il faut θ̇(x). Il faut bien comprendre que θ est en fait une
fonction de (x, `) qui d’ailleurs est en fait une fonction de x − `, donc

∂θ ∂θ v
θ̇ = v = − v = cos θ , (7)
∂` ∂x h

car v = `˙ est la vitesse d’avancement du front. Le moment d’inertie de chaque domino


est I = mh2 /3 et donc
Z
I dx
Ec = θ̇2 (8)
2 s
Z
mhv 2 1 dx
= 2 (9)
6s cosh ((` − x)/h) h
mhv 2
= tanh(`/h) (10)
6s
4) Cette hypothèse est discutable car les chocs ne sont pas élastiques en vrai, et puis le
glissement des dominos les uns sur les autres apporte aussi sont lot de forces de Coulomb.
L’autre hypothèse extrême qui consiste à dire que les dominos restent collés sans rebondir
est aussi incorrecte quand on voit les films de dominos au ralenti. Bref, la conservation
de l’énergie mécanique est un cas extrême peu réaliste. En utilisant que tanh(x) → 1
quand x → ∞ et en gardant le terme dominant dans Ep on trouve
p
v → 3g` . (11)

Comme seulement les dominos sur une longueur typique h basculent et donc ont de
l’energie cinétique, mais que l’énergie potentielle apportée est linéaire en ` car provenant
de la bascule de tous les dominos déjà tombés, la vitesse augmente.

NB) En réalité, l’énergie n’est en fait pas conservée, ni d’ailleurs l’impulsion, ni


d’ailleurs le moment cinétique, voir arXiv:physics/0401018.
MQ - réflexion - transmission harmomique

Un jet d’électrons d’énergie E > 0 se propage selon l’axe x dans une région de potentiel nul
pour x < −a.  x = −a et x = a il rencontre une région de potentiel (d’énergie potentielle)
 2 Entre
x
V (x) = V0 a2 − 1 avec V0 > 0. Pour x > 0 on a de nouveau un potentiel nul. On supposera

~2 (E + V0 )2
= . (1)
2ma2 V0

m 2 2
a) En posant V0 = 2
ω a quelle interprétation pouvez-vous donner à cette condition (1) ?

b) Déterminer autant que possible la fonction d’onde dans les différentes régions, et discuter
comment obtenir les coefficients de réflexion et de transmission.

1
Un exercice pas facile qui demandera de guider le candidat à certaines étapes.
m 2 2 m 2 2
a) En posant V0 = 2
ω a on a V (x) = 2
ω (x − a2 ) ce qui montre déjà que ω est la fréquence
(pulsation) de l’oscillateur classique correspondant. Observons aussi que V0 est la profondeur
maximale du puit harmonique et que E + V0 est l’énergie des électrons mesurée depuis le
minimum de ce puit harmonique. La condition (1) devient alors
~2 ~2 m 2 2  ~ω 2 ~ω
(E + V0 )2 = V0 = ω a = ⇒ E + V0 = . (2)
2ma2 2ma2 2 2 2
~ω m 2 2
Mais 2
est l’énergie du niveau fondamental dans un potentiel harmonique V (x) = 2
ω x,
mesurée depuis son minimum, si −∞ < x < ∞. Ici, le choix de E correspond donc précisément
à l’énergie de cet état lié. Mais comme E > 0, ici il n’y a pas “d’état lié à l’intérieur du puits”.
Comme on verra après, ce choix de E implique que la solution de Schrödinger entre −a et a sera
bien donnée par la fonction d’onde du fondamental de l’oscillateur harmonique correspondant.

b) On a donc 3 régions, avec

région 1 : x < −a : V (x) = 0 ,


m
région 2 : −a < x < a : V (x) = ω 2 (x2 − a2 ) ,
2
région 3 : a<x : V (x) = 0 (3)

et l’équation de Schrödinger stationnaire est


~2 d2 φ
− + V (x)φ(x) = Eφ(x) . (4)
2m dx2
~2 k2
On posera comme d’habitude E = 2m
et on trouve pour les régions 1 et 3 :

φ1 (x) = Aeikx + Be−ikx , φ3 (x) = F eikx , (5)

sans terme e−ikx dans la région 3, vu qu’il n’y a pas de flux venant de +∞. Dans la région 2
on a
~2 00 m 2 2 m ~ω
− φ2 + ω x φ2 = (E + ω 2 a2 )φ2 = (E + V0 )φ2 = φ2 , (6)
2m 2 2 2
où on a utilisé la relation (2). On peut simplifier un peu pour obtenir

−φ002 + α2 x2 φ2 = αφ2 , α= . (7)
~
A ce stade on peut laisser le candidat chercher un peu pour trouver la solution, mais assez
2
rapidement il conviendra de suggérer d’essayer une solution de la forme eβx .
2
Si on essaie donc φ2 = Ceβx on trouve φ02 = 2βxφ2 et φ002 = (4β 2 x2 + 2β)φ2 ce qui donne
α2 α
−4β 2 x2 − 2β + α2 x2 = α ⇒ β2 = et β = − . (8)
4 2

2
2 /2 2 /(2~)
Les deux équations sont compatibles et donnent β = −α/2, soit φ2 = Ce−αx = Ce−mωx .
Le fait que les deux équations sont compatibles vient évidemment du choix particulier de E.
2
Autrement dit, avec un choix différent de E, la simple fonction eβx n’aurait pas été solution.
Lorsque l’on cherche les états liés dans un potentiel donné (à une dimension) et lorsque’on
a trouvé une solution normalisable on se contente de celle-ci. Mais ici, la question de la nor-
malisabilité ne se pose pas et on a besoin de la deuxième solution, linéairement indépendante.
On peut encore laisser le candidat chercher un peu, mais il faudra très probablement suggérer
2 /2
assez vite de chercher une deuxième solution de la forme φ(x) = f (x)e−αx .
−αx2 /2
 −αx2 /2
Insérons donc dans (7) l’ansatz φ2 (x) = f (x)e . On a (f e )00 = f (α2 x2 − α) −
 2 2
2αxf 0 + f 00 e−αx /2 et comme e−αx /2 est déja solution de (7) on obtient pour f une équation
plus simple :
f 00 = 2αxf 0 . (9)

Evidemment, f = C est solution et redonne juste notre solution précédente. Mais on a aussi la
2 Rx 2
nouvelle solution f 0 (x) = Deαx et donc f (x) = D −a dzeαz + C. Donc, la solution générale
pour φ2 est Z x
2 /2 2 /2 2
φ2 (x) = Ce−αx + DF (x) , F (x) = e−αx dz eαz . (10)
−a
2 /2
Notons que D/C est homogène à 1/a. Notons également que F (−a) = 0 et F 0 (−a) = Deαa
2 /2
et que F (a) peut être calculé en fonction de a et F 0 (a) = −αaF (a) + Deαa .
On peut maintenant imposer la continuité de φ et de φ0 en x = −a ce qui donne:
2 /2
Ae−ika + Beika = Ce−αa (11)
2 /2 2 /2
ik(Ae−ika − Beika ) = Cαae−αa + Deαa . (12)

Continuité en x = a donne
2 /2
F eika = Ce−αa + DF (a) (13)
−αa2 /2 αa2 /2
ikF eika = −Cαae + De − αaDF (a) . (14)

On a donc 4 équations pour 5 inconnues, et on peut exprimer B, C, D et F en fonction de A.


Les coefficients de réflexion et transmission sont alors donnés par

|B|2 |F |2
R= , T = . (15)
|A|2 |A|2

3
Petits mouvements dans une 1/2-sphère en rotation

On considère une 1/2-sphère de rayon R qui tourne autour de son axe de révolution à une
vitesse angulaire constante Ω. Une particule (*) de masse m se déplace sur la surface intérieure
de la 1/2-sphère, sous l’influence de la gravité. De plus, la particule est soumise à une force
de frottement: −k~vrel où ~vrel est la vitesse de m relativement à la 1/2-sphère.

On s’intéresse aux petits mouvements de la particule proche du bas de la 1/2-sphère. Pour


cela, on introduit le référentiel galiléen cartésien (x, y, z) avec son origine en bas de la 1/2-
sphère et l’axe (Oz) suivant l’axe de rotation. On se placera dans l’approimation |x| << R
et |y| << R.

A quelle(s) condition(s) ce mouvement est périodique? Discuter sa stabilité.

(*) Ici, particule = petite bille suffisamment petite pour négliger tout effet dû à sa taille.
On comprend bien ce qui se passe. Si k = 0 (pas de frottement) alors, le fait que 1/2-
sphère tourne ou pas n’a pas d’importance. Le mouvement de la particule est un mouvement
de petites oscillations autour de O avec la pulsation ω 2 = g/R. A un moment, il faudra
retrouver ce résultat évident.

~ ):
Maintenant, si k 6= 0, les équations du problème ne sont pas compliquées (avec ~r = OM

m~r¨ = m~g + N
~ − k~vrel .

Et la vitesse de m est: ~r˙ = ~vrel + Ω


~ ∧ ~r. Il n’y a pas besoin d’écrire les vecteurs en coor-
données sphériques (énoncé) car on s’intéresse uniquement à ce qui se passe proche de O. En
particulier: N ' mg, puis on projette l’équation du mouvement sur (x) et (y):

mẍ = Nx − kvrel,x .

mÿ = Ny − kvrel,y .

En utilisant la figure, on obtient facilement: Nx = −mgx/R et Ny = −mgy/R puis (vrel,x , vrel,y )


s’obtient facilement à partir de l’expression de la vitesse plus haut. Au final, on trouve:
x
mẍ = −mg − k ẋ − kΩy.
R
y
mÿ = −mg − k ẏ + kΩx.
R
A partir de là, c’est pratiquement fini. On cherche à quelles conditions un mouvement
périodique est possible, soit: x = x0 eiωt et y = y0 eiωt . Ce qui va donner un système de
2 équations pour (x0 , y0 ) qui admet une solution non nulle seulement si:

(−ω 2 + g/R − ikω/m)2 + (kΩ/m)2 = 0.


p
Soit: −ω 2 +g/R−ikω/m = ±ikΩ/m. Ce qui n’est possible que si les 2 conditions: ω = g/R
p
et Ω = ±ω si k 6= 0. Et si k = 0, seulement ω = g/R reste. Ce qui conclut.
Condensateur avec diélectrique tenu par un ressort

Dans un milieu homogène, appelé “diélectrique linéaire”, la loi de Gauss est modifiée en remplaçant
Q Q
0
par 0 r
où r est appelé la permittivité relative du milieu.

a) Déterminer la capacité d’un condensateur plan plongé dans un milieu diélectrique linéaire ho-
mogène de permittivité relative r > 1. Déterminer l’énergie contenue dans un tel condensateur
lorsqu’on applique une différence de potentiel V0 aux deux plaques.

b) On considère un condensateur formé de deux plaques rectangulaires de dimensions a et b, parallèles


et séparées d’une distance e. Un diélectrique rigide de permittivité relative r , d’épaisseur e et de
même dimensions a et b que les plaques peut coulisser entre les plaques sans frottement. Il est relié
d’un côté à un ressort de raideur k et de longueur à vide telle que la force de rappel soit nulle quand
le diélectrique est juste à l’entrée du condensateur. On charge le condensateur avec une charge Q
et on l’isole. Quelle(s) est/sont la/les position(s) d’équilibre du diélectrique? Peut-il y avoir des
ocillations et si oui, avec quelle fréquence? On pourra éventuellement supposer que r − 1 est petit,
puis regarder l’AN avec r − 1 = 10−1 , Q = 10−5 C, a = b = 20 cm, e = 5 mm et k = 10 N m−1 .

1
a) On fait l’argument habituel des deux plaques 1 et 2, parallèles, de normale ~ez , chargées de Q et
−Q. Soient leurs dimensions a × b et leur séparation e. La première produit un champ E ~ 1 = E1~ez
au dessus et E~ 1 = −E1~ez en dessous. Par Gauss modifié on a 2E1 ab = Q et idem pour le champ
0 r
produit par la deuxième plaque. En dehors des plaques les champs s’annulent et entre les plaques,
ils s’ajoutent, donc
~ =E
~1 + E
~2 = 2 × Q
E . (1)
2ab0 r
La différence de potentiel est V = E × e et donc la capacité
Q ab
C= = 0 r . (2)
V e
Pour déterminer l’énergie stockée, on peut calculer le travail qu’il faut fournir pour changer Q →
Q+δQ, en prenant une charge δQ de la plaque avec charge −Q et en la transportant jusqu’à la plaque
1
de charge Q. Le travail est EδQe et l’énergie change donc de E → E + δE avec δE = EdδQ = C
QδQ.
En intégrant de 0 à Q on trouve
Q2 C
E= = V2 , (3)
2C 2
comme dans le vide, la seule modification est le r contenu dans l’expression de la capacité. Pour
l’énergie volumique entre les plaques on trouve
 2
E 0 r V 0 r ~ 2
u= = = E . (4)
abe 2 e 2
En dehors du condensateur, le champ électrique s’annule. La présence du diélectrique en dehors
du condensateur n’a donc pas d’influence, et les résultats précédents restent valables même si le
diélectrique ne remplit que l’espace entre les plaques.

b) Pour résoudre cet exercice il faut raisonner sur l’énergie totale. Si le candidat essaie d’expliciter
les forces, il conviendra de lui suggérer, après un temps convenable, une approche énergétique.

Appelons x la coordonnée qui donne la position du bord du diélectrique qui est entré dans le conden-
sateur, telle que x = 0 correspond au cas du diélectrique juste à l’entrée du condensateur. La force de
rappel du ressort est donc F~ = −kx~ex . Soit a la longueur selon x. Le condensateur correspond donc
à deux condensateurs en parallèle, l’un avec dimensions x et b et permittivité  = r 0 et l’autre avec
i Si
dimensions a − x et b et permittivité 0 (ou air ' 0 ). La capacité de chacun est donné par Ci = e
où Si est l’aire de la partie des plaques sans ou avec le diélectrique: C1 = r e0 xb et C2 = 0 (a−x)b
e
.
Q 2
b0
Pour l’ensemble C = C1 + C2 = e
[r x + (a − x)]. Ensuite on peut raisonner sur l’énergie: Eel = 2C
et Emec = 12 kx2 . L’énergie potentielle du système est la somme des deux:

eQ2 1 1
Epot (x) = + kx2 . (5)
2b0 (r − 1)x + a 2

2
La position d’équilibre correspond au minimum de Epot :

0 eQ2
Epot (xe ) = 0 ⇒ η a3 = xe [(r − 1)xe + a]2 , η≡ (r − 1) . (6)
2b0 ka3
Q2
On vérifiera que le paramètre η est bien sans dimension car 0 a2
et ak sont bien homogènes à une
force. L’équation (2) est cubique. Elle a donc 3 solutions mais seules les solutions avec 0 ≤ xe ≤ a
ont un sens. Pour un équilibre stable la deuxième dérivée doit être positive. On a, en se servant de
(2)  
00 2xe
Epot (xe ) =k +1 . (7)
(r − 1)xe + a
Toute position d’équilibre avec xe ≥ 0 est donc stable.

Si r − 1  1 on peut approximer [(r − 1)xe + a)]2 ' a2 et la solution de (2) est xe ' η a, à
condition, bien sur, que η ≤ 1. (Evidemment, r − 1  1 n’implique pas forcément η  1.) Les deux
autres “solutions” sont à xe = O( ra−1 ) ce qui est en dehors de l’intervalle [0, a] autorisé. On a alors
00
Epot (xe ) ' k[1 + 2xae ] ' k[1 + 2η] > 0. Si m est la masse du diélectrique, on aura des oscillations avec
r r
00
Epot (xe ) k(1 + η)
ω= ' . (8)
m m
1
Pour l’AN, on a 0 = = 8, 8 × 10−12 C2 J−1 m−1 et on trouve η ' 0, 2. On estime la masse
c2 µ0 q
k
m = ρ abe à 0, 4 kg (avec un ρ = 2g cm−3 ) et donc m = 5 s−1 et ω = 5, 5 s−1 .

3
Jet sur un coin

On considère un jet d’eau bidimensionnel impactant un coin d’angle intérieur 2α. On note
L sa dimension transverse. Loin en amont, la vitesse du jet est U0 et son épaisseur h0 ,
symétrique par rapport au coin. On considère que l’écoulement est stationnaire, que le fluide
est incompressible et de masse volumique ρ. On négligera les effets de la gravité.

1) En supposant le fluide parfait, déterminer la force exercée par le jet sur le coin.

2) On ne suppose plus le fluide parfait. Ceci étant, on suppose que les effets de viscosité sont
non-négligeables uniquement dans une couche (limite) d’épaisseur δ(x) << h0 au niveau des
bords du coin (x repérant la progression le long d’un bord dans la direction du fluide). Donc,
au-delà de l’épaisseur δ(x), le fluide est considéré comme parfait.

En déduire la force du jet sur le coin jusqu’à une extension du coin x = a, le long de chaque
bord du coin. Cette expression fait intervenir en particulier δ(a).

Indication: le cours indique l’expression de l’épaisseur d’une couche limite sur une dimension
√ q
ηx
caractéristique L: δ = L/ Re, ce qui se traduit ici par: δ(x) = ρU .
Question 1: cours. Question 2: variation par rapport au cours. La notion de couche limite
est dans le cours, en particulier la formule rappelée dans l’énoncé. Simplement, on en fait ici
une utilisation qui n’est pas directement dans le cours.

1) Le fluide est incompressible, on a donc conservation du débit volumique. On prend une


section transverse (h0 L) de fluide en amont et les deux sections transverses à une position x
sur chaque bord: hL (le jet reste symétrique de part et d’autre du coin):

U0 h0 L = 2U hL.

De plus, en utilisant le théorème de Bernoulli le long de la surface libre: U = U0 , donc


h = h0 /2.

Ensuite, pour calculer la force extérieure exercée sur le fluide, on peut considérer un volume
donné de fluide qui se déplace en fonction du temps. Voir cours (systèmes ouverts). On
considère le volume entre une section h0 L en amont et deux sections hL à une position x
le long de chaque bord. La variation de quantité de mouvement pour tout le fluide dans ce
volume est la force extérieure exercée sur le fluide contenu dans ce volume. Donc ici, la force
exercée par le coin sur le fluide dans ce volume plus les forces de pression (qui s’exercent sur
le volume de contrôle).

Les seules contributions non nulles sont celles qui viennent des sections transverses. On obtient
en projetant sur la direction du jet initial (~eX ):
Z
d
[ ρ~v (.) dV ] = −ρU02 h0 L + ρU 2 hL(2 cos α) = ρU02 h0 L(cos α − 1)~eX .
dt V
C’est la force de l’extérieur sur le fluide. Les forces de pression sur le fluide s’expriment comme:
R
− pd2 S,~ où l’intégrale est prise sur la surface (fermée) qui délimite le volume de contrôle.
Mais la pression est partout égale à p0 , donc cette contribution est nulle. Finalement, la force
du coin sur le fluide est:
F~ = ρU02 h0 L(1 − cos α)~eX .
Elle est bien orientée suivant ~eX . De plus, elle est bien nulle pour α = 0 et maximale pour
α = π/2.
2) On va d’abord écrire la conservation du débit volumique dans les nouvelles conditions d’un
fluide non parfait. On suit le fluide le long d’un bord du coin, circulant le long de x. On
prend l’axe (y) perpendiculaire à (x). En y = 0 la vitesse du fluide est donc nulle (condition
à la limite pour un fluide visqueux), puis cette vitesse augmente de y = 0 à y = δ(x). Au-delà
de δ(x) le fluide est à nouveau parfait et sa vitesse est U0 (question 1)). Donc, la vitesse dans
le fluide visqueux peut s’écrire au 1er ordre:
y
vx (y) = U0 .
δ(x)

La conservation du débit volumique s’écrit donc (pour un bord):


Z δ(x) Z h(x)
h0
Q1 = L vx (y)dy + U0 dy = L U0 .
0 δ(x) 2

Attention: h(x) est égal à h0 /2 uniquement pour x = 0. Ensuite, h(.) doit être une fonction
de x pour que la conservation du débit volumique soit possible. On obtient:
1
h(x) = (h0 + δ(x)).
2

On reprend le même raisonnement qu’en 1) pour déterminer la force du coin (jusqu’à la


position x) sur le fluide:
Z h(x)
F = ρU02 h0 − 2 cos α vx2 (y)dy.
0
R h(x)
Il reste à calculer: 0
vx2 (y)dy, en utilisant la méthode ci-dessus. On obtient:
Z h(x)
1 1
vx2 (y)dy = U02 h0 (1 − δ(x)).
0 2 3

Donc la force du coin (jusqu’à la position a sur chaque bord) sur le fluide devient:

1 δ(a)
F = ρU02 h0 [1 − cos α + cos α ].
3 h0
Ce qui conclut.
Tube capillaire

De l’eau se trouve dans un très grand bac. Un tube cylindrique, d’axe vertical, en verre fin, de
rayon intérieur R, et ouvert aux deux bouts, se trouve avec l’un de ses bouts juste au dessus
de la surface de l’eau. A t = 0 en enfonce rapidement (mais sans faire de vagues) le bout du
tube cylindrique de h0 = 0, 1 mm dans l’eau. On observe que l’eau monte à l’intérieur du
tube cylindrique d’une certaine hauteur h(t) (mesurée à partir du bout du tube). La hauteur
atteinte à l’équilibre et la façon dont cette hauteur est atteinte dépendent du rayon intérieur R
du tube. La figure montre h(t) pour différentes valeurs de R allant de 0,2 mm à 0,5 mm.

Expliquer qualitativement ce qui se passe, puis décrire quantitativement le phéomène.

1
Equilibre statique à t = ∞ : Il y a une tension de surface A entre l’eau et le verre qui résulte
en un gain d’énergie proportionnel à la surface de contact, soit ∆E = −2πRhA. Si cette notion
n’est pas familière on donnera ce résultat. On a donc une force F~A = −∇(∆E)
~ = 2πRA~uz
dirigée ver le haut. Il y a le poids de l’eau qui est montée, proportionnel au volume, soit
πR2 (h − h0 ), et dirigée vers le bas. A l’équilibre cela donne (h → h∞ ), avec ρ la masse
volumique de l’eau,
2A
2πRA = πR2 (h∞ − h0 )ρg ⇒ h∞ − h0 = . (1)
ρg R
Regardons comment on pourrait obtenir cela à partir de Navier-Stokes. A l’équilibre on a ~v = 0
et on a seulement
~ − ρg~uz + . . . ,
0 = −∇p (2)

où + . . . représente d’autres forces. Ici on a la force due à la tension de surface, mais elle
n’est pas volumique, ni même proportionnelle à la suface de contact. En effet, c’est l’énergie
qui est proportionnelle à la surface de contact, et la force qui en est moins le gradient est
seulement proportionnelle à la circumférence du cylindre. Il faut alors intégrer (2) sur le volume
de la colonne d’eau et ajouter la force F~A = 2πRA~uz , ce qui donne (avec p(h) = patm et
p(0) = patm + ρgh0 )
 2A 
0 = πR2 − p(h) + p(h0 ) − ρgh + 2πRA = πR2 ρg(h0 − h) + , (3)
R
où h = h∞ , et on retrouve (1).
Comme h0  h∞ on pourra négliger h0 et on obtient en bonne approximation
2A
h∞ ' , (4)
ρg R
2A
soit h∞ R = ρg
est constant. On voit assez bien sur les figures que c’est bien le cas: on a
(en mm) (R, h∞ ) = (0.2, 7.2), , (0.3, 4.8), (0.4, 3.6) et (0.5, 2.8) , soit pour les produits Rh∞ :
14, 4 mm2 à chaque fois. Avec ρ = 103 kg m−3 et g = 9.8 m s−2 on trouve A ' 0.7 kg m−1 .

Approche de l’équilibre aux petits t : Au début, la force résultante est non-nulle et accélère
l’eau vers le haut. A cause de son inertie elle va monter plus haut que la position d’équilibre,
puis redescendre et osciller autour de la hauteur h∞ . Ces oscillations sont amorties par les
forces de viscosité. Ces dernières sont proportionnelles à la surface de contact avec le verre
(∼ R), alors que l’inertie est proportionnelle aux volume (∼ R2 ). Les oscillations seront donc
plus amorties pour des petites valeurs de R, ce qui est effectivement ce qu’on observe.
Essayons d’être plus quantitatif. Lorsque l’eau monte dans le tube, on pourrait essayer de
décrire le champ des vitesses dans le tube. A cause de la viscosité la vitesse sera plus grande

2
au centre et diminue vers le bord, ce qui impliquera que l’eau qui monte plus vite au milieu,
une fois arrivée en surface, s’étale sur toute la section du tube. Mais “loin” de la surface on a
~v = v(t, r)~uz avec v(t, r) independant de z à cause de l’incompressibilité qui donne ∇~ · ~v = 0,
R R
donc ∂v
∂z
= 0. Le flux à travers une section πR2 est donc 2π 0 dr r v(r, t) ce qui sera égal à
πR2 ḣ(t) :
Z R
2
ḣ(t) = 2 dr r v(r, t) . (5)
R 0

Evidemment, pour v indépendant de r on trouve simplement ḣ = v.


Si on ne tient pas compte des forces de tension de surface, l’équation de Navier-Stokes
projetée sur ~uz est
dv ∂p
= −ρg − ρ + η∆v . (6)
dt ∂z
Comme avant, les forces de tension de surface ne sont pas localisables et le plus simple est
d’intégrer l’équation (6) sur tout le volume de l’eau dans le tube et d’y ajouter la force due à
la tension de surface qui est 2πRA selon ~uz . On utilisant aussi (5) on obtient:
Z R
2 2 2

ρ πhR ḧ = −πhR ρg − πR p(h) − p(0) + η2πh dr r ∆v + 2πRA . (7)
0

En utilisant comme avant p(h) − p(h0 ) = −ρgh0 et en simplifiant par πR2 on a


Z
2A 2h R
ρhḧ = gρ(h0 − h) + +η 2 dr r ∆v . (8)
R R 0

Pour v = ḣ = 0 on retrouve évidemment (3) et (1). Le Laplacien en coordonnées polaires est



∆v(r) = 1r ∂r

r ∂v
∂r
et donc
Z R
∂v  R ∂v
dr r ∆v = r = R (R) . (9)
0 ∂r 0 ∂r
(C’est normal que l’intégrale du Laplacien sur le volume donne un terme de bord: c’est le
R
théorème de Gauss qui relie d3 x∆v au flux de ∂v ~e à travers le cylindre.) On a donc
∂r r

2A 2h ∂v
ρhḧ = gρ(h0 − h) + +η (R) . (10)
R R ∂r
∂v
Que vaut ∂r
(R) ? Avec un flux régulier dans un tube on aurait un profil des vitesses v(r) =
r2

v(0) 1 − R2 et donc ∂v
∂r
(R) = − 2v(0)
R
. On relie v(0) à ḣ en évaluant (5) pour ce profil de vitesses
v(0)
et on obtient ḣ = 2
. Donc ∂r (R) = − 4Rḣ . En insérant cette valeur dans (10) on a
∂v
finalement
l’équation différentielle pour h:

2A 8hḣ
ρhḧ = gρ(h0 − h) + −η 2 . (11)
R R
3
On voit déjà l’essentiel: plus R est petit, plus le terme de viscosité est important, en accord avec
notre discussion qualitative ci-dessus et en accord avec les figures. Cette équation différentielle
est non-linéaire, mais à l’approche de la solution d’équilibre on peut la linéariser en posant
h(t) = h∞ + (t), où h∞ est solution de (1). On obtient:

2 g ρR2
¨ + ˙ + ω 2  = 0 , ω2 = , τ= . (12)
τ h∞ 4η
On a bien un oscillateur amorti avec fréquence ω et temps d’amortissement τ . Sur la figure
g
pour R = 0.5 mm on estime 7 T ' 2 s, soit ω ' 20 s−1 (alors que ω 2 = h∞
donne ω ' 55 s−1 –
où est l’erreur ???) et τ ' 0.8 s pour R = 0.5 mm, ce qui donne une estimation de la viscosité
de l’eau à η ' 0.3 × 10−3 N sm−2 , alors que η à température ambiante vaut 10−3 N sm−2 , soit
encore une erreur par un facteur 3...

4
Optique géométrique

On considère une onde électromagnétique dont le champ électrique prend la forme


~ =E
E ~ 0 e(iωt−iϕ(x,y,z)) . (1)

On suppose que la phase ϕ varie beaucoup plus rapidement que l’amplitude, donc on
~  |∇
suppose que |∇ϕ| ~E~ 0 |/|E
~ 0 | (Par exemple |∂x ϕ|  |∂x E0 x |/|E0 x |). On suppose éga-
lement que les dérivées secondes de ϕ sont très petites devant ses dérivées premières et
~ ∇ϕ|
donc |∇ ~ ~ 2 (par exemple |∂xx ϕ|  |∂x ϕ|2 ). On veut montrer que sous ces
 |∇ϕ|
~
hypothèses, ∇ϕ/| ~
∇ϕ| décrit la direction de propagation d’un rayon lumineux.

~ est orthogonal à E
— Montrer que dans un milieu électriquement neutre ∇ϕ ~ 0 et que
~
le vecteur de Poynting est orienté selon ∇ϕ.

~ · ∇ϕ~ = ω2
— Montrer que dans le vide ∇ϕ c2
. Comment cette relation est modifiée dans
un milieu d’indice optique n ?

~ = ~k = k~u où ~u est un vecteur unitaire. Montrer que dans un milieu


— On pose ∇ϕ
~ u=∇
d’indice n variable (~u · ∇)~ ~ ln n − ~u(~u · ∇
~ ln n).

— Justifier qualitativement pourquoi on peut considérer que ~u décrit la direction d’un


rayon lumineux.

— Quel doit être le gradient vertical d’indice optique dans l’atmosphère pour qu’un
rayon lumineux puisse faire le tour de la Terre ?
Eléments de solution

Un exercice sur l’approximation eikonale.

~ = 0 donne lorsqu’on néglige les dérivées de l’amplitude


1) divE
~ ·E
∇ϕ ~ = 0. (2)

~ E
Ensuite rot ~ donne toujours en négligeant les dérivées de l’amplitude
~ = −∂t B

~ ∧E
∇ϕ ~ = ωB
~ (3)

si on écrit le champ magnétique comme on a écrit le champ électrique. On a donc


~ ∝E
Π ~ ∧B
~ ∝ ∇ϕ
~ . (4)
~
L’énergie s’écoule selon ∇ϕ.

~ = −∆E
2) On fait comme d’habitude rotrotE ~ si ρe = 0 et sans courants et on trouve
 
1 ~ = 0.
∆ − 2 ∂tt E (5)
c

En négligeant les dérivées secondes de ϕ devant le carré des dérivées premières, et en


négligeant la variation de l’amplitude, il vient
2
∇ϕ ~ =ω .
~ · ∇ϕ (6)
c2
Dans un milieu on aura plutôt
2 2
∇ϕ ~ =ω n .
~ · ∇ϕ (7)
c2

3) On part de (7) et on en prend le ∇. ~ D’abord sous la forme k 2 = ω 2 n2 /c2 on a


simplement
2
~ = ω n ∇n
k ∇k ~ ~ = ω ∇n
∇k ~ . (8)
c2 c
Il faut voir qu’on avait négligé les dérivées secondes de ϕ, mais seulement par rapport
aux dérivées premières ou par rapport à celles de l’amplitude, mais pas par rapport à
l’indice optique. Donc là on va garder les dérivées secondes. Pour faire le calcul on peut
soit être à l’aise avec les opérateurs, soit passer en coordonnées Cartésiennes. Il va falloir
faire commuter les dérivées secondes pour obtenir
  2
~ ∇ϕ
∇ ~ · ∇ϕ
~ = 2(~k · ∇) ~ = 2 ω n ∇n
~ ∇ϕ ~ . (9)
c2
Or
(~k · ∇) ~ = ~u(~k · ∇)k
~ ∇ϕ ~ + k(~k · ∇)~
~ u (10)
donc en utilisant (8) et après division de (9) par k 2 il vient
~ u=∇
(~u · ∇)~ ~ ln n − ~u(~u · ∇
~ ln n) (11)

4) On peut évoquer le sens d’écoulement de l’énergie selon ~u. On peut aussi faire le lien
avec les photons
E = ~ω = pc p = ~k ~ .
k = |∇ϕ| (12)
On peut également rajouter que le fait de négliger les gradients de l’amplitude revient
à décrire un rayon lumineux dont la taille transverse (normalement supposée infiniment
petite car c’est une ligne) est beaucoup plus grande que la longueur d’onde 2π/k. On a
donc une échelle mésoscopique qui est la taille du rayon lumineux. A l’intérieur du rayon
lumineux, son étendue mésoscopique étant suffisamment grande, tout se passe comme
dans une onde plane, mais seulement localement et la direction de cette onde plane évo-
lue.

5) On applique la loi de réfraction (11). Le gradient est vertical selon ~er en coordonnées
polaires et ~u = ~eθ . Donc le membre de gauche donne

∂θ~eθ ~er
=− (13)
RT RT

et le membre de droite donne (∂r ln n)~er . Il faut donc un gradient


n 1
∂r n = − '− (14)
RT RT
Trou noir 2

Le 10 avril 2019 la collaboration Event Horizon Telescope a publié la première photo d’un trou
noir, celui au centre de la galaxie M87. Ce trou noir a une masse d’environ 6.5 × 109 MS où
MS ' 2 × 1030 kg est la masse du soleil. Mais la masse d’un trou noir n’est pas nécessairement
si élevé.
Selon la relativité générale (classique, c’est à dire non quantique), un trou noir se caractérise
par une attraction gravitationnelle si forte dans l’espace qui l’entoure que rien, ni même la
lumière ne peut s’en échapper, une fois à l’intérieur de “l’horizon” du trou noir.
Cependant, en combinant relativité (gravitation) et mécanique quantique, Hawking a démon-
tré que, par un effet quantique relativiste, un trou noir doit néanmoins émettre un rayonnement
de corps noir caractérisé par la température de Hawking TH . Cela veut dire que le trou noir
se comporte comme un corps noir qui émet un rayonnement de corps noir à cette température
TH . Dans la suite on admettra simplement ce résultat.

a) Rappeler la loi de Stefan-Boltzmann. Que peut-on conclure pour la puissance rayonnée par
le soleil ? A quelle perte de masse par an cela correspond-il ? Commentaire ? On donne
σ = 5, 67 × 10−8 W m−2 K−4 .

b) Hawking trouve que pour un trou noir (sphérique) de masse M la température de Hawking
1 1
est TH ∼ M
. Déterminer la relation correcte, à un facteur numérique près (qui s’avère être 8π
).

c) La taille d’un trou noir est caractérisé par son rayon de Schwarzschild RH ∼ M . Comme
ci-dessus, trouver la bonne relation entre RH et M à un facteur numérique près (qui sera 2).

d) En considérant le trou noir de masse M comme un corps noir sphérique de rayon RH et de


température TH , calculer la puissance rayonnée. Quelle masse devrait avoir le trou noir pour
que cette puissance soit comparable à celle du soleil ? A quelle perte de masse correspond ce
rayonnement? Calculer la durée de vie d’un trou noir par cette “évaporation de Hawking”.

1
a) On a j = σT 4 où j est le flux d’énergie, i.e. la puissance par unité de surface. Le rayon du
soleil est environ 100 fois celui de la terre, soit RS ' 7 × 108 m. On trouve pour la puissance
totale rayonnée:
PS = 4πRs2 σTS4 , (1)
b
où on prendra TS = 6000 K. (Si besoin, on peut discuter la loi de Wien λmax = T
avec
b b
b ' 2, 9 × 10−3 K m qui donne pour la lumière jaune T = λmax
= 600 nm
' 5000 K.)
On trouve alors, avec T = 6000 K, PS ' 4, 5 × 10 W. Par la relation E = mc2 bien connue
26

dm
de tous, cette puissance correspond à une perte de masse par unité de temps de dt
= − cP2 '
5 × 109 kg s−1 ' 1, 5 × 1017 kg an−1 . Ceci est très faible devant MS et il faudrait donc de l’ordre
de 1013 années pour que toute la masse soit rayonnée. Evidemment la durée de vie d’une étoile
comme le soleil est plus de l’ordre 109 années, temps au bout duquel tout l’hydrogène sera
transformé en Hélium (et autres).

b) C’est un exercice sur les unités. D’abord, il faut voir que le phénomène fait intervenir,
d’après l’énoncé, la relativité (vitesse de la lumière c), la gravitation (constante de Newton G)
et la mécanique quantique (~). Donc, on pose

A
TH = , A = ~α cβ Gγ . (2)
8πkB M
En effet, ~, c ou G ne font intervenir le Kelvin, et donc la seule façon de relier les grandeurs
A
est de relier l’énergie kB TH à 8πM
.
Il faut ensuite tout exprimer en kg, m et s et on aura 3 équations pour 3 inconnues:

[~] = J s = N m s = kg m2 s−1 , [c] = m s−1 ,


[G] = N m2 kg −2 = m3 s−2 kg −1 , [kB TH ] = J = kg m2 s−2 (3)

On trouve alors le système de 3 équations (pour les kg) 1 = α−γ−1, (pour les m) 2 = 2α+β+3γ
et (pour les s) −2 = −α − β − 2γ, soit α = 1, β = 3, γ = −1 et donc

~c3
TH = . (4)
8πkB GM

c) Posons de même, RH = 2BM avec [B] = m kg −1 et on voit simplement que G/c2 a la bonne
unité, soit
2GM ~c
RH = ⇒ TH = . (5)
c2 4πkB RH
d) Comme pour le soleil, on écrit pour la puissance totale rayonnée par le trou noir
 4
2 σ ~c 1
P = 4πRH σTH4 = 2
(6)
(4π)3 kB RH

2
~c
On a kB = Rgaz /NAvogadro ' 8
6×1023
J K−1 ' 1, 3 × 10−23 J K−1 , ~ ' 10−34 J s et donc kB
'
2, 3 × 10−3 K m. On trouve alors
 2
−22 m
P ' 8 × 10 W . (7)
RH

Pour un trou noir de la masse du soleil on a RH ≡ RS ' 3 km et donc P ' 10−28 W, à comparer
avec les 4, 5 × 1026 W rayonnés par le soleil. Ce ne sera donc qu’au bout de 1028+26+13 =
1067 années qu’une fraction considérable de la masse de ce trou noir aura été rayonnée. En
revanche, un trou noir qui sera 1030 fois plus petit (soit d’un kg environ) rayonnera une fraction
considérable de sa masse en rayonnement de Hawking en seulement 10 millions d’années. Et
plus la masse diminue, plus TH augmente et plus il rayonnera fort car en termes de M on peut
réécrire (6) comme
 4  4
σ ~c c4 c2 σc2 ~c
P = ≡ γ où γ = (8)
(4π)3 kB 4G2 M 2 M2 4(4π)3 G2 kB

On peut même calculer la durée qu’il faut pour que le trou noir ait rayonné toute sa masse.
γ
Comme avant on a − dM
dt
= P
c2
= M2
ou −M 2 dM = γdt, soit intégré entre 0 et tf et entre
Mi = M et Mf = 0 on obtient la durée de vie tf comme

M 3 = 3γ tf . (9)

3
Transfert thermique cercle

On considère un solide en forme de tore, c’est à dire en forme de bouée et dont les
sections ont une aire S. On suppose que le rayon des sections de ce tore est négligeable
devant le rayon principal du tore et on décrira ce solide par un cercle de rayon R. La capa-
cité thermique est notée c, la conductivité thermique est notée λ et la masse volumique ρ.

— Etablir l’équation du transfert de la chaleur en fonction de θ en coordonnées po-


laires.
2 2
— Vérifier que T (θ, t) = 2√TπDt
0`
e−R θ /(4Dt) avec D = λ/ρc est solution de cette équa-
tion. Pourquoi ne peut-on pas l’utiliser ?

— On cherche une solution sous la forme T (θ, t) = g(t)f (θ). Déterminer les fonc-
tions f et g.

— Initialement une moitié de cercle est à la température T0 tandis que l’autre moitié
est à température nulle. Trouver une solution pour l’évolution de la température.
Quelle est son évolution asymptotique au temps longs ?

— On part d’une température nulle dans tout le tore, et on injecte une quantité
d’énergie Q en un point donné. Quelle est l’évolution de la température et son
profil asymptotique ?

La fonction fa (θ) définie sur [−π, π] et ayant pour valeur π/a sur −a/2 ≤ θ ≤ a/2 et 0
ailleurs est égale à la série

1 X sin(an/2)
fa (θ) = + cos(nθ) (1)
2 an/2
n=1
Eléments de solution

1) On peut faire une mise en équation par bilan. Il faut juste voir que si on prend une
coordonnée s sur le cercle ∂s = 1/R∂θ . Ou alors on connait le gradient et la divergence
en polaire en voyant que l’hypothèse de sections très petites implique un gradient de
température orthoradial. On trouve
λ
ρc∂t T = λ∂ss T = ∂θθ T . (2)
R2
On introduit un coefficient de diffusion D = λ/ρc.

2) On vérifie sans problème que c’est solution de l’équation en dérivant. Mais cela n’est
pas compatible avec la topologie du problème. Pourtant on a bien T (−π, t) = T (π, t) !
Mais c’est au niveau des flux que ça va pas. Il ne sont pas continus en −π et en π vu
qu’ils s’éloignent tout deux de la valeur θ = 0.

3) Il faut donc passer à autre chose. On passe par une méthode de séparation des variables.
Si on cherche T (θ, t) = f (θ)g(t), alors
g 0 (t) D f 00 (θ)
= 2 . (3)
g(t) R f (θ)
Ces combinaisons doivent être des constantes. On a donc
f 00 (θ)
= Cte . (4)
f (θ)
Les solutions sont des cos et sin (on élimine les cosh et sinh) mais les conditions pério-
diques imposent cos(nθ) et sin(nθ). On a alors
g 0 (t) n2 D
=− 2 (5)
g(t) R
donc
2 Dt/R2
g(t) ∝ e−n . (6)
On trouve donc les solutions de type
2 Dt/R2 2 Dt/R2
e−n cos(nθ) e−n sin(nθ) (7)
4) Par superposition une somme de telles solutions est une solution. La condition initiale
suggérée est la fonction de l’aide évaluée en a = π et multipliée par T0 et c’est donc

X sin(nπ/2)
T0
T (θ, 0) = + T0 cos(nθ) (8)
2 nπ/2
n=1

T0 2T0 X (−1)p
= + cos[(2p + 1)θ] (9)
2 π (2p + 1)
p=0
et la solution générale à tous temps est

T0 2T0 X (−1)p 2 2
T (θ, t) = + cos[(2p + 1)θ]e−(2p+1) Dt/R . (10)
2 π (2p + 1)
p=0

C’est le terme p = 1 qui a la décroissance la moins rapide. C’est logique car ce mode a les
gradients les plus petits (et les distances entre les points hauts et bas sont les plus grands).

5) On peut refaire le raisonnement précédent mais a < π c’est à dire avec une zone
chaude plus petite que la zone froide initiale. La fonction donnée dans l’aide à la même
intégrale même si on varie a. Quand on avait a = π on avait injecté Q = RπρcST0 . Si on
injecte cette même énergie sur une longueur Ra au lieu de Rπ on a la condition initiale

X sin(an/2)
T0
T (θ, 0) = + T0 cos(nθ) (11)
2 an/2
n=1

avec
Q
T0 = . (12)
RπρcS
Dans la limite d’une injection d’énergie en un seul point du cercle a → 1 (et on uti-
lise sinc(a) → 1) alors la solution générale est (en rajoutant les facteurs d’évolution
temporelle)

X ∞
T0 2 2
T (θ, 0) = + T0 cos(nθ)e−n Dt/R (13)
2
n=1

Pareil c’est le mode le plus bas n = 1 qui guide la forme asymptotique. Pour Dt > R,
donc quand la longueur de diffusion a dépassé le rayon du cercle, alors on entre dans ce
régime asymptotique.
Induction sphère tournante

~ = Ω~ez ,
On considère une sphère conductrice (conductivité σ), de rayon R, en rotation Ω
~
et plongée dans une champ magnétique uniforme B = B0~ex .

Quelle est la distribution de courants dans la sphère lorsque celle-ci est devenue sta-
tionnaire ?

Aide : La loi d’Ohm pour un conducteur en mouvement dans un champ électrique et un


~ + ~v ∧ B).
champ magnétique est ~j = σ(E ~
Eléments de solution

En régime stationnaire, toutes les dérivées temporelles des champs sont nulles et donc
~ = 0.
rotE (1)

Tout d’abord le champ magnétique extérieur est irrotationnel, et donc il dérive d’un
potentiel
B ~
~ = gradΦ Φ = B0 x . (2)
~ D’abord il faut ~v .
En sphérique x = r sin θ cos φ. On calcule ~v ∧ B.
~ ∧ ~r = r sin θΩ~eφ .
~v = Ω (3)

En utilisant ~eφ = cos φ~ey − sin φ~ex , il vient

~ = −r sin θ cos φΩB0~ez = −xΩB0~ez .


~v ∧ B (4)

On regarde maintenant la loi d’Ohm.


~ + ~v ∧ B)
~j = σ(E ~ (5)

En stationnaire la conservation de la charge implique

div~j = 0 . (6)

Aussi sur le bord de la sphère en r = R, on a jr = 0. Regardons d’abord cette deuxième


condition. Puisque

~ · ~er
(~v ∧ B) = −R sin θ cos θΩB0 cos φ , (7)
r=R

on en déduit que

−∂r V |r=R = Er (r = R) = R sin θ cos θΩB0 cos φ . (8)

Regardons la condition de conservation du courant. D’abord


~ = 0.
div(~v ∧ B) (9)

C’est évident en utilisant l’expression Cartésienne mais pas facile sinon ! Donc on déduit
que divE~ = 0 et donc avec E ~
~ = −gradV , on a ∆V = 0. Là on peut deviner un peu la
solution à partir de (8). On suppose que la condition de courant qui n’est pas selon ~er est
valable partout et donc (8) est valable pour tout r. On vérifiera après que ça fonctionne.
On trouve alors
1 1
V = − r2 ΩB0 sin θ cos θ cos φ = − ΩB0 zx . (10)
2 2
~ = ΩB0 (x~ez + z~ex ) .
E (11)
2
On peut vérifier immédiatement sur cette forme que c’est sans divergence et que ∆V = 0
en utilisant les expressions Cartésiennes. La densité de courant est

~j = σΩB0 (z~ex − x~ez ) . (12)


2
Les lignes de courant sont des cercles de normale ~ey . Ces courants induisent un champ
magnétique qui va induire un champ électrique et donc d’autres courants etc... La pre-
mière itération peut se démarrer. On a le champ additionnel

~ B
rot ~ 1 = µ0~j . (13)

Pour résoudre on peut reprendre le rot de cette relation en utilisant que la divergence
d’un champ magnétique est nulle et la fameuse relation rotrot. On trouve
~ 1 = µ0 rot~j = µ0 σΩB0~ey .
−∆B (14)

~ 1 = (ax2 + by 2 + cz 2 )~ey et en utilisant (13) on a facilement


On tente alors B

~ 1 = − 1 µ0 σΩB0 (x2 + z 2 )~ey .


B (15)
4
On suit alors la même procédure mais c’est beaucoup plus compliqué !

NB) Voir American Journal of Physics, 86, 443 (2018).


Oscillations dans un bassin

Un bassin fait 25 m de long et 15 m de large. Le fond du bassin est horizontal. La profondeur


de l’eau au repos est de 2 m. Déterminer les fréquences des petites oscillations de la surface de
l’eau.

1
Dans un premier temps il faut montrer qu’on peut travailler avec un potentiel qui satisfait
l’équation de Laplace avec des conditions aux bords du bassin (conditions de Neumann), puis
trouver les solutions. Ensuite il faudra exploiter Euler pour obtenir une condition de bord à la
surface libre, ce qui donnera la relation de dispersion. On guidera ou ne guidera pas le candidat
...
On appelera a, b et H la longueur, largeur et hauteur au repos et le fond du bassin sera
~ · ~v = 0 et en particulier ρ
à z = 0. L’eau sera considérée comme incompressible et donc ∇
est constant. De plus, pour de petites oscillations le champ des vitesses est de rotationel nul,
~ ∧ ~v = 0, comme on montre facilement: Pour de petites oscillations on a Euler linéarisé

∂~v ∂~v 
ρ ~ − ρ~g
= −∇p ⇒ ~ p + gz .
= −∇ (1)
∂t ∂t ρ
~ v)
~ ∧
On conclut que 0 = ∇ ∂~v
= ∂(∇∧~
et donc que ∇ ~ ∧ ~v est une constante indépendante du
∂t ∂t
temps qui doit être nulle pour des mouvements oscillatoires. ∇ ~ ∧ ~v = 0 implique que ~v dérive
d’un potentiel φ et ∇~ · ~v = 0 que ce potentiel satisfait l’équation de Laplace:

~
~v = ∇φ , ∆φ = 0 . (2)

Les conditions aux bords d’absence de vitesse normale deviennent


∂φ ∂φ ∂φ ∂φ ∂φ
= = = = =0. (3)
∂x x=0 ∂x x=a ∂y y=0 ∂y y=b ∂z z=0
L’équation de Laplace (2), ces conditions de bord (3), et le fait qu’on cherche des oscillations
suggèrent de chercher des solutions à variables cartésiennes séparées de la forme

φ(x, y, z, t) = f (x)g(y)h(z)e−iωt . (4)

On pourra aussi bien travailler directement avec la partie réelle cos ωt. Une fois les solutions
obtenues, on pourra les superposer. L’équation de Laplace ∆φ = 0 donne
f 00 g 00 h00 f 00 g 00 h00
+ + =0 ⇒ = −α2 , = −β 2 , = α2 + β 2 ≡ γ 2 . (5)
f g h f g h

Le choix des signes des constantes −α2 et −β 2 est le seul qui permet de satisfaire les deux
conditions de bord pour f et les deux pour g. On obtient
r
nπ mπ n2 m 2
f (x) = A cos x, g(y) = B cos y, h(z) = C cosh γnm z , γnm = π + 2 , (6)
a b a2 b
avec n et m des entiers positifs ou nuls. Evidemment, l’équation de Laplace ne nous apprend
rien sur les fréquences ω.

2
Jusqu’ici on n’a pas tenu compte de la forme des forces de gravité; seule l’incompressibilité
et l’approximation linéaires ont été utilisées. Il reste encore à exprimer la condition qu’à la
surface libre on a p = p0 . Evidemment, cette dernière condition sera cruciale et devra faire
intervenir la gravité g, absente des équations (2) et (3). Considérons un élément de fluide à la
~
surface dont l’équation est z = s(x, y, t). L’équation d’Euler linéarisée (1) donne avec ~v = ∇φ
 ∂φ p 
~
∇ + + gz = 0 . (7)
∂t ρ
∂φ
On conclut que + ρp + gz = ψ(t) est indépendant des coordonnées. On peut toujours changer
∂t
Rt
φ par une fonction de t uniquement sans changer ~v et en remplaçant φ → φ0 = φ − ψ on
arrive à
∂φ p
+ + gz = 0 , (8)
∂t ρ
(où on écrit de nouveau φ au lieu de φ0 ). Pour un élément de fluide à la surface cela devient
∂φ  p0
x, y, s(x, y, t), t + + gs(x, y, t) = 0 , (9)
∂t ρ
En prenant une dérivée par rapport au temps on a
∂ 2φ ∂ 2 φ ∂s ∂s
2
+ +0+g =0, (10)
∂t ∂t∂z ∂t ∂t
∂s ∂φ
toujours pris en z = s(x, y, t). Mais ∂t
= vz |z=s = |
∂z z=s
et l’équation précédente devient
∂ 2φ ∂ 2 φ ∂φ ∂φ
+ + g =0 , en z = s(x, y, t) . (11)
∂t2 ∂t∂z ∂z ∂z
∂ 2 φ ∂φ ∂vz
Le deuxième terme est ∂t∂z ∂z
= v
∂t z
et, dans l’approximation des petites oscillations, doit
être négligé devant les deux autres termes, en accord avec l’utilisation de l’équation d’Euler
sans le terme ~v · ∇~~ v . De plus, on peut évaluer (11) en z = H au lieu de z = s(x, y, t), car
pour les petites oscillations s − H est petit. (On peut avoir l’impression que remplacer s par
H revient à supprimer l’existence même des oscillations, mais ce n’est évidemment pas le cas,
∂φ
car la non-nullité du terme = vz en z = H est le signe d’une surface non plane. Plus
∂z
Rt
mathématiquement, on a φ(x, y, s) − φ(x, y, H) ' ∂φ
∂z
(x, y, H) (s − H) et s − H ' dtvz ' vωz ,
vz2
soit φ(x, y, s) − φ(x, y, H) ' ω
. Donc, à des termes d’ordre v 2 près on pourra remplacer
φ(x, y, s) par φ(x, y, H).) Finalement, la condition (11) implique pour la solution (4) et (6):
 
− ω 2 cosh(γnm z) + gγnm sinh(γnm z) z=H = 0 , (12)

soit la relation de dispersion

ω 2 ≡ ωnm
2
= g γnm tanh(γnm H) . (13)

3
Double oscillateur

On considère deux particules ponctuelles de même masse m se déplaçant uniquement selon


l’axe x dans un potentiel harmonique “extérieur”, à une seule dimension, c’est à dire tel que
leur énergie potentielle est
m 2 2
ω (x1 + x22 ) ,
Vext (x1 , x2 ) = (1)
2
où x1 et x2 sont les coordonnées des deux particles. A cette énergie potentielle “extérieure”
s’ajoute une énergie potentielle d’interaction entre les deux particules
m 2
Vint (x1 , x2 ) = α (x1 − x2 )2 . (2)
4
a) Décrire le mouvement classique de ces deux particules. (On supposera quelles peuvent se
traverser mutuellement sans aucune incidence sur leur mouvement.)

Maintenant on essaie de décrire le même système en mécanique quantique par une fonction
d’onde ψ(t, x1 , x2 ).

b) Quelle interprétation proposez-vous pour cette fonction d’onde dépendant des coordonnées
des deux particules ?

c) Quelle équation de Schrödinger pourrait-on écrire pour cette fonction d’onde ? Proposez une
manière pour déterminer les niveaux d’énergie de ce système de deux particules quantiques.

1
Il faudra, bien sûr, guider le candidat si besoin, après le temps de réflexion personnelle nécessaire.
a) On écrit x21 + x22 = 12 (x1 + x2 )2 + 12 (x1 − x2 )2 et donc
 2
2 x1 + x2 m
Vtot (x1 , x2 ) = Vext (x1 , x2 ) + Vint (x1 , x2 ) = mω + (ω 2 + α2 )(x1 − x2 )2
2 4
 
x1 + x2
≡ Vcm + Vrel (x1 − x2 ) . (3)
2
On voit donc apparaitre deux oscillateurs découplés, l’un avec masse 2m pour la coordonnée du

barycentre X = x1 +x
2
2
, et un autre avec masse m2 et fréquence Ω = ω 2 + α2 , pour la coordonnée
relative x = x1 − x2 . Il est donc clair que X et x décrivent séparément des oscillations avec
leurs fréquences respectives. Montrons-le : Le PFD donne
∂Vtot 1 0 0 ∂Vtot 1 0 0
mẍ1 = − = − Vcm − Vrel , mẍ2 = − = − Vcm + Vrel (4)
∂x1 2 ∂x2 2
et donc
0 m 0 m
2mẌ = −Vcm (X) = −2mω 2 X ẍ = −Vrel
, (x) = − (ω 2 + α2 )x , (5)
2 2
m
où on voit la masse totale 2m et la masse relative 2 .

b) L’analogie avec une particule à deux dimensions nous semble évidente. Mais il convien-
dra d’interagir et de guider le candidat si cette interprétation ne lui vient pas tout seul :
|ψ(t, x1 , x2 )|2 dx1 dx2 est la probabilité de trouver la première particule en x1 à dx1 près et
la deuxième particule en x2 à dx2 près. (Evidemment on considère les deux particules comme
discernables.)

c) L’équation de Schrödinger à proposer doit se baser sur le principe de correspondence E = ~ω


et l’expérience avec l’équation pour une seule particule. Il s’agira donc de simplement remplacer
~ ∂ 2 2
le potentiel à une particule par le potentiel écrit ci-dessus et de comprendre que le terme − 2m ∂x2
~2 k2
qui donne pour une particule libre 2m
= Ecin doit ici être remplacé par l’énergie cinétique des
deux particules. Donc
∂ h ~2 ∂ 2 ∂2  i
i~ ψ(t, x1 , x2 ) = − + + Vtot (x1 , x2 ) ψ(t, x1 , x2 ) . (6)
∂t 2m ∂x21 ∂x22
On se satisfera aussi par l’équation stationnaire obtenu par ψ(t, x1 , x2 ) = e−iEt ϕ(x1 , x2 ) :
h ~2 ∂ 2 ∂2  i
E ϕ(x1 , x2 ) = − + + V tot (x 1 , x2 ) ϕ(x1 , x2 ) . (7)
2m ∂x21 ∂x22
Le traitement du problème classique suggère de passer en variables X et x. Pour le potentiel
c’est déjà fait avant. Pour l’opérateur différentiel il suffit d’écrire proprement (probablement ils
ne savent pas faire et c’est hors programme - donc il faudra les aider ou donner cette formule)
∂ ∂X ∂ ∂x ∂ 1 ∂ ∂ ∂ ∂X ∂ ∂x ∂ 1 ∂ ∂
= + = + , = + = − . (8)
∂x1 ∂x1 ∂X ∂x1 ∂x 2 ∂X ∂x ∂x2 ∂x2 ∂X ∂x2 ∂x 2 ∂X ∂x

2
Puis on prend les carrés et on somme :

∂2 ∂2 1 ∂2 ∂2
+ = + 2 . (9)
∂x21 ∂x22 2 ∂X 2 ∂x2

Schrödinger devient alors avec ϕ(x1 , x2 ) ≡ φ(X, x)


h ~2 1 ∂ 2 ∂2  i
E φ(X, x) = − + 2 + V cm (X) + V rel (x) φ(X, x) , (10)
2m 2 ∂X 2 ∂x2
et le tour est joué : on raconte le coup habituel de la séparation des variables : φ(X, x) =
φcm (X)φrel (x) et
h ~2 ∂ 2 i h ~2 ∂ 2 i
− +Vcm (X) φcm (X) = cm φcm (X) , − +Vrel (x) φrel (x) = rel φrel (x) , (11)
4m ∂X 2 m ∂x2
avec E = cm + rel . Les énergies d’un seul oscillateur à une dimension sont au programme et

donc cm = ~ω(n1 + 21 ) et rel = ~ ω 2 + α2 (n2 + 21 ), soit

1 √ 1
E = ~ω(n1 + ) + ~ ω 2 + α2 (n2 + ) . (12)
2 2
Dans la limite où l’interaction entre les deux particules disparait (α = 0) on retrouve évidemment
E = ~ω(n1 + n2 + 1) comme il faut.

3
Conditions de l’effet dominos

On considère une ligne de dominos alignés verticalement sur une table. Il sont séparés
d’un espace de longueur s, ont une largeur D et une hauteur h. On suppose que les
dominos ne peuvent que basculer mais pas glisser sur la table.

— Lorsque les dominos sont renversés les uns sur les autres, relier l’angle de bascu-
lement (par rapport à la verticale) d’un domino en fonction de celui de ses voisins.

— Quel est l’angle final des dominos une fois qu’ils sont tous tombés ?

— A quelles conditions sur la forme des dominos ce renversement de tous les domi-
nos est possible ?

On considère deux dominos placés verticalement l’un à côte de l’autre et séparés par un
espace de longueur s. Le premier domino a une largeur D et une hauteur h mais le second
est proportionnellement r fois plus grand.

A quelle condition le premier domino peut-il entrainer la chute du second domino ?

Aide : le moment d’inertie d’un solide rectangulaire (de côtés a et b et de masse m)


associé à une rotation autour d’un de ses coins est I = 31 m(a2 + b2 ).
Eléments de solution

\ est droit. On a AB =
1) On exprime de deux manières différentes AB où l’angle ABC
(D + s) cos θn+1 − D mais aussi AB = h sin(θn − θn+1 ) donc

h sin(θn − θn+1 ) = (D + s) cos θn+1 − D . (1)

D’autres méthodes fonctionnent aussi bien et toutes nécessitent simplement de la trigo-


nométrie élémentaire.

2) Lorsque les dominos sont tous renversés θn+1 = θn et donc


D
cos θf = . (2)
D+s
3) Il faut voir si c’est énergétiquement possible et donc si l’énergie potentielle finale est
plus basse ou pas. On trouve facilement pour un domino en regardant où est son centre
de masse
mgh
Ep = [cos θ − 1 + (D/h) sin θ] (3)
2
On évalue Ep < 0 pour θf et on trouve la condition
Dp
s> 2sD + s2 (4)
h
qui se récrit
2D3
s> . (5)
h2 − D2
On peut discuter la possible annulation du dénominateur. Clairement il faut au moins
h > D. En fait les dominos doivent se toucher donc il faut au moins aussi avoir h > s+D.
Plus précisément, ils doivent même pouvoir reposer l’un sur l’autre donc la contrainte
est un peu plus faible et p
h > (s + D)2 − D2 (6)
qui donne p
s < −D + h2 + D2 . (7)
On note que (5) et (7) donnent l’encadrement (trop serré ça bascule pas assez et ça fait
que monter, trop espacé ça ne se touche plus en tombant)
2 p
< −1 + α2 + 1 (8)
α2 − 1
avec α = h/D. Avec un peu d’algebre, cela donne la condition

α2 (α2 − 3) > 0 (9)



donc il faut h > 3D.

4) On donne juste les grandes lignes à suivre. On peut supposer qu’on met le domino
1 juste sur sa diagonale, en position instable, prêt à tomber. On démarre donc avec un
angle initial tel que
h D
cos θ1,i = √ sin θ1,i = √ (10)
h2 + D2 h2
+ D2
Le choc a lieu quand sin θ1 = s/h. Par conservation de l’énergie totale, et en utilisant
la différence de Ep entre temps initial et impact on peut connaître l’énergie cinétique à
l’impact (en utilisant le rappel sur le moment d’inertie et 1/2Iω 2 )
1
Ec = m(h2 + D2 )ω12 . (11)
6
Ensuite ce qu’il faut voir c’est que les deux dominos sont liés (hypothèse donnée aux
candidats à ce stade) et donc les angles sont reliés par un raisonnement similaire à (1).
On trouve
h sin(θ1 − θ2 ) = (rD + s) cos θ2 − rD . (12)
Les variations de θ1 sont liés à celles de θ2 par différenciation

dθ1 h cos(θ1 − θ2 ) − (rD + s) sin θ2


= (13)
dθ2 h cos(θ1 − θ2 )

En particulier à l’impact θ2 = 0 et
dθ1
= 1. (14)
dθ2
Les rotations sont les mêmes ! Pour finir on peut argumenter que le moment cinétique
total est conservé. En effet, l’impact correspond à une force horizontale entre les deux
solides dans le cas idéalisé. Et donc le bras de levier de cette force de 2 sur 1 ou 1 sur
2 est le même, c’est la hauteur du point d’impact. La variation du moment cinétique de
1 est donc opposée à celle de 2. En utilisant cette conservation on peut en déduire la
vitesse de rotation commune ω juste après l’impact

I1 ω1 = (I1 + I2 )ω. (15)


Il faut donc calculer le rapport I2 /I1 sachant que le domino 2 et le meme que le 1 mais
avec une dilatation d’un facteur r. La masse va comme r3 (ou r2 si la dimension transverse
est gardée semblable) et donc I2 /I1 = r5 (ou r4 ). Une fois déduit ω, il suffit de regarder
l’énergie cinétique totale
1
Ec = (I1 + I2 )ω 2 (16)
2
et vérifier que ça suffira à aller jusqu’au point d’énergie potentielle maximale. Cette éner-
gie potentielle totale est la somme de (3) pour 1 et 2, mais avec l’angle θ1 qui est asservi
à l’angle θ2 via (12). Le point d’énergie potentielle maximale est un peu avant le moment
où 2 à son énergie potentielle maximale (sa diagonale verticale) car quand 2 monte, 1
descend.

NB) Voir arXiv:1301.0615 pour une résolution complète.


Ressort massif

Un ressort de constante de raideur K et de longueur au repos L a en plus une masse par unité
de longueur ρ. Ce ressort est suspendu verticalement, fixé à son extrémité supérieure et laissé
libre sans masse accrochée à son extrémité inférieure. Dans le champ de gravité uniforme g,
il s’étire donc du fait de son propre poids: c’est-à-dire qu’un point du ressort situé en z (avec
g = 0) se retrouve en z + s(z) (s(z) > 0).

1) Déterminer s(z) lorsque l’équilibre est atteint.

On coupe alors le champ de gravité (cela peut se faire par exemple en mettant le système
dans un ascenseur en chute libre).

2) Déterminer alors autant que possible s(z, t) sachant que la condition initiale est s(z, 0) =
s(z) de la question 1).
Pour la question 1), il est possible qu’un exercice de ce type ait été traité en classe. Il y a
plusieurs manières de faire, j’en propose une ici. La question 2) est plus subtile. Ici, j’attends
de voir comment le candidat prend des initiatives, propose des idées sans qu’il soit utile de
trouver la solution complète pour que je considère que c’est très bien.

L’axe (z) est l’axe vertical, parallèle au champ de gravité. Le point clef de l’exercice est de
bien s’imaginer la situation physique. Chaque petit morceau du ressort dz a un poids non nul
dmg et ce petit poids va donc tirer le ressort vers le bas. Et donc chaque morceau peut être
considérer comme une petite masse dm et un petit bout de ressort de longueur dz. C’est-à-
dire que le problème est exactement équivalent à considérer une suite de petite masses dm
attachées les unes aux autres par des petits ressorts de longueur dz de constante de raideur
kdz . Ce problème est alors très classique.

1) Dans un premier temps, on repère chaque position des petites masses par un indice discret
i et on note si l’allongement du ressort à cette position. On considère N masses le long
du ressort telles que dm = M/N où M = ρL est la masse du ressort. Aussi dz = L/N .
L’équation du mouvement de la masse dm = ρL/N située à la position i (position si (t) entre
2 autres petites masses) s’écrit:

dms¨i = dmg + kdz (si+1 + si−1 − 2si ).

Pour calculer kdz à partir de K la constante de raideur du ressort dans son ensemble, on peut
juste réaliser qu’il est beaucoup plus difficile de détendre un ressort tout petit qu’un grand.
Donc, kdz doit être de la forme: kdz = KL/dz. L’équation du mouvement devient:
KL
ρs¨i = ρg + (si+1 + si−1 − 2si ).
dz 2
Mais à la limite où N devient très grand, donc dz très petit, on peut écrire:
∂ 2 s(z, t)
si+1 + si−1 − 2si = s(z + dz, t) − s(z − dz, t) − 2s(z, t) = dz 2
∂z 2
Au final, l’équation du mouvement devient:
KL ∂ 2 s(z, t)
s̈(z, t) = g + .
ρ ∂z 2
On cherche s(z) lorsqu’il n’y a plus de variation temporelle (en régime permanent), donc:

∂ 2 s(z) gρ
2
=− .
∂z KL
Il nous faut donc 2 conditions aux limites pour répondre à la question. On a évidemment:
s(0) = 0. De plus, pour la dernière masse en bas du ressort, l’équilibre donne (dans le cas
discret): kdz (sN − sN −1 ) = dmg. Donc, en passant la limite continue:
KL 0
s (z)|z=L dz = ρdzg.
dz
Soit: s0 (z)|z=L = 0. On en déduit s(z) pour toutes les valeurs de z le long du ressort:

ρg z2
s(z) = (z − ).
K 2L

2) On coupe le champ de gravité, soit g = 0 dans l’équation du mouvement, qui devient:

KL ∂ 2 s(z, t)
s̈(z, t) = . (1)
ρ ∂z 2
Cette fois-ci, on n’est plus en régime permanent (à l’équilibre) car, ayant coupé g, le ressort
ne peut plus rester dans l’état de la question 1). On part donc de la condition initiale
ρg z2
s(z, t = 0) = K
(z − 2L ) (que l’on note s0 (z) dans la suite) puis s(z, t) varie suivant l’équation
(1). Cette équation est une équation de propagation (ici pour la fonction s(z, t)) avec la
p
vitesse de propagation: v = KL/ρ.

On peut écrire la solution générale de ce type déquation sous la forme de 2 ondes progressives
l’une allant dans le sens des z croissants et l’autre dans le sens des z décroissants, soit:

s(z, t) = f (z − vt) + g(z + vt).

La condition initiale donne: s(z, 0) = s0 (z) = f (z) + g(z). De plus, le champ de gravité étant
coupé à l’instant initial, on a: ∂s(z, 0)/∂t = 0, donc: f 0 (z) = g 0 (z) = s00 (z)/2. On en déduit:
f (z) = s0 (z)/2 + C et g(z) = s0 (z)/2 − C, où C est une constante. D’où:
1
s(z, t) = (s0 (z − vt) + s0 (z + vt)),
2
pour laquelle la constante s’élimine. Il reste à voir ce que donne les conditions aux limites en
ds
z = 0 et z = L qui sont comme en 1): s(0, t) = 0 et dz
(L, t) = 0. En z = L, il n’y a pas de
problème. Par contre, en z = 0, il semble qu’il y ait une difficulté. Si on prend simplement
1
l’expression fonctionnelle de s0 (.) plus haut et on calcule: (s (−vt)
2 0
+ s0 (vt)), alors cette
expression n’est pas nulle quelque soit t. Le point clef ici, c’est que physiquement s0 (−vt) a
un sens si on prend la fonction s0 (.) et on remplace l’argument par −vt mais physiquement,
cette opération n’a pas de sens. Ceci car s0 (z) de la question 1) est définie uniquement
pour z > 0. Donc, s0 (−vt) ne correspond à rien de physique. Finalement, le fait que la
fonction 21 (s0 (−vt) + s0 (vt)) ne soit pas identiquement nulle n’est pas un problème, car elle ne
correspond pas à s(0, t) pour le ressort. Une fois cette observation faite, on comprend que la
relation plus haut qui donnait: f (z) = s0 (z)/2 (si on prend la constante nulle) n’est valable
que pour z > 0. Pour z < 0 on est libre de choisir f (z) différemment et on peut parfaitement
choisir: f (z < 0) = −s0 (|z|)/2. Et, avec ce choix, la condition s(0, t) = 0 est bien vérifiée.

Si jamais un candidat souhaite utiliser comme approche une solution sous forme d’onde sta-
tionnaire, ce qui va donner les modes propres après vérification des conditions aux limites,
c’est bien aussi... Jusque là, c’est possible dans le cadre du programme. Mais il faut préciser
que ce n’est pas l’approche demandée car elle se heurte à la propriété hors programme qui
est nécessaire pour cette approche: que la solution générale de l’équation de propagation est
alors la somme discrète de tous les modes propres. De plus, trouver les coefficients de cette
somme qui sont compatibles avec la condition initiale est également hors programme.

C’est pourquoi, je pense qu’une fois que l’équation (1) est établie, et une fois qu’il est observé
que c’est une équation de propagation, il faut préciser au candidat que l’approche demandée
est celle des ondes progressives.
Un appontage rudimentaire

On considère un système de freinage de secours pour l’atterrissage d’un avion. Il s’agit


d’une chaîne parfaitement flexible de masse linéique λ = 10kg/m et disposé en longueur
dans le sens de l’atterrissage (voir schéma). L’avion peut accrocher l’extrémité de cette
chaîne et l’entrainer dans son mouvement, si bien qu’une partie de la chaîne est en mou-
vement et le reste est immobile. On décrira l’avion par une masse ponctuelle de masse
M = 104 kg et on note V0 = 100 m/s la vitesse de l’avion au moment où il accroche la
chaîne. On suppose que la chaîne glisse sans frottement sur le tarmac.

— Quelles sont les équations régissant le mouvement de l’avion atterrissant en sup-


posant qu’il n’y a aucune force de frottement sur l’avion ? Les résoudre. L’avion
peut-il s’arrêter uniquement avec ce dispositif ?

— On suppose de plus que l’avion atterrit sur le ventre et est soumis à un frottement
solide de coefficient f = 0.5. En combien de temps s’arrête-t-il ? Et sur quelle
distance ?

— Comment évolue l’énergie cinétique du système ?

— Réaliser un modèle microscopique de la chaîne pour rendre compte de ce qui se


passe.
Eléments de solution

L’énergie cinétique n’est pas conservée, si bien que la dernière question initie une discus-
sion à ce sujet (une fois cette propriété établie) permettant de faire le lien avec d’autres
parties du programmes. On peut considérer qu’il s’agit d’un exercice sur les systèmes
ouverts.

1) Lorsque l’avion a avancé de x, il a entrainé une masse λx/2 (voir schéma) et force
cette partie de la chaîne à aller à sa vitesse ẋ. Donc le PFD sans forces extérieures donne

(M + λx/2)ẍ + λẋ2 /2 = 0 . (1)

En effet, la masse en mouvement est (M + λx/2) et se déplace à une vitesse ẋ. On peut
le voir soit comme un système ouvert, soit comme une système fermé avec une partie à
vitesse nulle et une partie à vitesse ẋ. Eq. (1) est plus clair si écrit clairement comme la
conservation de la quantité de mouvement, et
d
[(M + λx/2)ẋ] = 0 . (2)
dt
Pour résoudre il suffit de trouver la constante

(M + λx/2)ẋ = M V0 . (3)

Il vient alors
(M + λx/2)dx = M V0 dt (4)
et donc (on met l’origine des temps à l’accrochage de la chaîne)
" r #
2M V0 tλ
M x + λx2 /4 = M V0 t , x= −1 + 1 + . (5)
λ M

On trouve que la distance ne fait que croitre avec le temps même si la vitesse décroit
puisque
V0 V0
ẋ = q = . (6)
1 + VM
0 tλ 1 + λx/(2M )

2) On peut refaire la même chose avec un frottement f M g.

d
[(M + λx/2)ẋ] = −f M g , (7)
dt
(M + λx/2)ẋ = M V0 − f M gt . (8)
On voit déjà qu’on s’arrête quand
V0
tstop = (9)
fg
donc en 20s avec les données du problème (un atterrissage sur le ventre violent...). On
ne trouve aucune dépendance en M ni en λ et cela appelle un commentaire. C’est le
même résultat que sans le cable, car il est immobile au début et à la fin. En revanche, la
distance d’arrêt va dépendre de λ donc on a quand même un mouvement différent. On
intègre en
M x + λx2 /4 = M V0 t − f M gt2 /2 , (10)
et on trouve alors " r #
2M λ
x= −1 + 1 + (V0 t − f gt2 /2) . (11)
λ M
Et on peut utiliser tstop pour avoir la distance d’arrêt.
 s 
2M  2
λV0 
xstop = −1 + 1 + . (12)
λ 2f M g

3) Le travail des forces de frottement est simplement W = −f M gx et donc le travail


total est Wtot = −f M gxstop . L’énergie cinétique est
1
Ec = (M + λx/2)ẋ2 . (13)
2
Son évolution est avec (7)

Ėc = ẋ(M ẍ + λxẍ/2 + λẋ2 /4) = ẋ(−f M g − λẋ2 /4) (14)

Le premier terme de cette dernière expression est la puissance des forces de frottements,
et s’il n’y avait que lui on aurait ∆Ec = −f M gxstop = Wtot . Mais on a le second terme
qui indique qu’on a perdu de l’énergie cinétique ailleurs. C’est d’ailleurs pour ça qu’on a
une distance d’arrêt plus courte. On a moins besoin du travail des forces de frottement
car il y a quelque part d’autres forces non conservatives qui ont agi. On remarque que
dans le cas initial sans frottement f = 0 on a Ėc = −λẋ3 /4. Ca ne suffit pas à arrêter
l’avion tout seul mais l’énergie cinétique est réduite.

4) On peut voir le cable ou la chaîne comme une série d’oscillateurs harmoniques couplés.
Quand on rajoute un nouvel élement du cable, c’est comme si on rajoutait un nouvel
oscillateur à la chaîne, mais sans lui donner la même condition initiale de vitesse que les
autres qui ont déjà été entrainés. Si on modélise le problème de cette manière on peut
passer au continu en notant ϕ(x) le champ de déplacement des oscillateurs. Si on ne met
aucune dissipation alors le passage au continu donnera

∂tt ϕ − c2 ∂xx ϕ = 0 , (15)

et notre rajout d’élément revient à mettre une condition initiale ϕ̇ < 0 juste au début de
la chaîne (dans le référentiel en mouvement avec la partie mobile de la chaîne). On va
avoir une onde sonore longitudinale qui va se propager. Mais une fois arrivée sur l’avion
elle va rebondir en arrière. C’est là qu’on voit qu’il faut rajouter quelque-chose car si la
perturbation initiale a fait un aller-retour, on n’a pas réussi à forcer l’élement rajouté à
suivre la vitesse ẋ.
Pour rendre compte du phénomène, une modification simple est de rajouter un frot-
tement. On peut imaginer que le cable ou la chaîne est comme une chaîne d’oscillateurs,
mais qui frottent sur qqch de manière visqueuse et lui transfèrent de la quantité de
mouvement. On aurait donc plutôt

∂tt ϕ + γ∂t ϕ − c2 ∂xx ϕ = 0 . (16)

et une onde amortie. L’amortissement conservant la quantité de mouvement totale, l’onde


s’arrête quand tous les éléments ont la même vitesse correspondant à cette conservation.
Au passage on a perdu de l’énergie cinétique dans ces frottements visqueux et cela a
chauffé le cable. Cela fonctionne si c  ẋ et donc si l’onde sonore a la temps de distri-
buer et uniformiser la quantité de mouvement. La vitesse du son dans les solides étant
plus grande que dans l’air, on est tranquilles dans le cas de l’avion qui ne va pas se poser
à Mach 3...

Une autre direction pour décrire le phénomène est d’invoquer directement une ana-
logie entre l’onde longitudinale dans le solide et les ondes acoustiques qui se dissipent
par viscosité. La discontinuité de vitesse sur l’élément qu’on rajoute induit un saut de
pression (et donc de tension), qui se propage mais se dissipe par viscosité. Avec cette
approche, l’équation pour la vitesse dans la chaîne vu comme une cable solide satisfait
(la même dérivation que d’habitude mais en gardant la viscosité)
η
∂tt v − c2 ∂xx v = ∂xx ∂t v . (17)
ρ
On peut étudier la relation de dispersion pour voir comment une condition initiale de
vitesse non-uniforme se dissipe.
MQ-Coulomb

On considère l’atome d’hydrogène en mécanique quantique. Écrire l’énergie potentielle de


l’électron dans le champ électrostatique du noyau que l’on considérera comme immobile. Dis-
cuter les symétries du problème et proposer une équation de Schrödinger stationnaire que l’on
pourra écrire pour la fonction d’onde ψ(~r). Discuter l’interprétation de ψ.

On peut montrer que lorsque l’on cherche des solutions à symétrie sphérique, on peut poser
u(r) 1
ψ(~r) = r


et alors u(r) satisfait une équation de Scrödinger “à une dimension” dans la seule
variable r, qui est l’analogue de l’équation unidimensionnelle habituelle (sauf que r ≥ 0). Ecrire
cette équation. Discuter l’interprétation de u(r) et en particulier la condition de normalisation
correspondante.

On cherche à déterminer l’énergie d’un état lié, à symétrie sphérique. Pour cela on essayera
une fonction u de la forme
u(r) = rα e−r/r0 (1)

dépendant d’un paramètre α. Déterminer dans quel cas ceci est une solution de l’équation de
Schrödinger et quelle est l’énergie correspondante ? Conclusion ?

1
e
Le potentiel électrostatique du noyau est φ(~r) = 4π0 r
et donc l’énergie potentielle de l’électron
dans ce champ est (−e)φ(~r), soit

e2
Epot (~r) = − . (2)
4π0 r
On a évidemment une symétrie sphérique. Mais cela n’implique pas que les solutions de
Schrödinger ont forcément aussi cette symétrie sphérique, seulement que l’on peut trouver
certaines solutions qui ont cette symétrie. L’équation de Schrödinger est (attention à 3 dimen-
sions l’équation n’est pas au programme, mais connu de presque tous les candidats ! Si besoin
la faire deviner en généralisant à partir d’une dimension.)
h ~2 e2 i
− ∆− ψ = Eψ , (3)
2m 4π0 r

où ∆ est le Laplacien. Discuter que |ψ(~r)|2 dxdydz = |ψ(r, θ, ϕ)|2 r2 dr sin θdθ dϕ est la proba.
R(r) 1 d2
Pour une solution ψ(~r) = √

ne dépendant que de r on a ∆R(r) = r dr2
(rR(r)) et
R(r) u(r)
Schrödinger pour ψ = √

= √
r 4π
devient effectivement pour u
h ~2 d2 e2 i
− − u(r) = Eu(r) , (4)
2m dr2 4π0 r
exactement comme à une dimension, sauf qu’ici r ≥ 0. La probabiité de trouver l’électron entre
r et r + dr est |ψ(r)|2 r2 dr4π = |u(r)|2 dr, ce qui explique le facteur √1 dans la normalisation

de u.
Essayons donc

u(r) = rα er /r0 ,
 
0 α−1 rα −r/r0
u (r) = αr − e ,
r0
   
00 α−2 rα−1 rα −r/r0 α(α − 1) 2α 1
u (r) = α(α − 1)r − 2α + 2 e = − + u(r) . (5)
r0 r0 r2 r0 r r02

On voit donc que u satisfait l’équation de Schrödinger (4) si et seulement si

~2 2α e2 ~2
α(α − 1) = 0 , = , − =E . (6)
2m r0 4π0 2m r02

La première égalité donne α = 0 ou α = 1. α = 0 est exclu par la deuxième égalité. Donc :


 2
4π0 ~2 m e2
α=1 , r0 = , E=− . (7)
me2 2 4π0 ~

On voit que r0 est le rayon de Bohr (53 pm) et E l’énergie de l’état fondamental de l’atome
d’hydrogène (−13, 6 eV).

2
Fabrication de barbapapa

Vous avez déjà vu une machine à fabriquer de la barbapapa ? Cela peut se modéliser
par à une passoire remplie de sucre liquide (car chauffé) et qui tourne à toute vitesse
autour d’un axe vertical. On veut modéliser le phénomène de formation des fibres de
sucre en faisant des hypothèses simplificatrices. La passoire peut être décrite par un pot
cylindrique de rayon r0 = 0.1 m (percé de nombreux trous sur son côté, mais pas au
fond) qui tourne autour de l’axe vertical à Ω = 100 rad/s. La viscosité du sucre chauffé
est µ = 1 Pa.s et sa densité la même que l’eau. Le sucre sort de chaque trou à une vitesse
radiale de u0 = 0.01 m/s. On cherche à décrire la fibre qui se forme par le fluide qui sort
d’un trou.

Estimer les différentes forces en jeu et montrer que la force d’inertie d’entrainement est
la plus importante.

Résoudre alors en ne considérant que cette force d’entrainement et vérifier le domaine de


validité des approximations.

Comment évolue le rayon de la fibre dans le régime stationnaire en fonction de la distance


à la passoire ?

En prenant en compte également la force de Coriolis (mais toujours en négligeant la


viscosité et le poids) trouver l’expression de la norme de la vitesse du sucre liquide dans
le référentiel tournant en fonction du temps de vol. On négligera aussi la vitesse radiale
initiale u0 . On utilisera que dans le référentiel Galiléen, la trajectoire du sucre est très
simple. Trouver également la forme de la fibre.
Eléments de solution

Ce sujet mélange mécanique (changement de référentiel éventuellement) pour une des-


cription balistique, mais également fluides visqueux pour comprendre les limitations de
cette approche.

1) On compare les accélérations associées à chaque force. Tout d’abord la force d’inertie
d’entrainement est d’ordre r0 Ω2 = 103 m/s2 soit cent fois plus que la gravité qu’on né-
gligera. Pour la force de Coriolis, on utilise qu’on a un peu de vitesse radiale en sortant
des trous donc Ωu0 = 1 m/s2 . Et pour la viscosité, il faut estimer µ∆v/ρ. Si on prend
∆v ' u0 /r02 ca donne 1 m/s2 . Mais on sait que le fluide accélère et donc on peut l’es-
timer avec la vitesse radiale Ωr0 plutôt que u0 et on trouve alors la même expression
que la force d’inertie et donc le même ordre de grandeur 103 m/s2 . On va voir que loin
de la passoire la force de Coriolis ne doit pas être négligée, mais pour ce qui est de la
viscosité on verra que c’est en sortie de passoire que l’on ne peut pas la négliger et c’est
la limitation de la description.

2) Avec ces hypothèses, dans le référentiel tournant on a θ̇ = 0 (coordonnées cylindriques)


et donc simplement
r̈ = Ωr . (1)
On considère que le temps d’éjection d’une particule dont on décrit la trajectoire est
te = 0. Les conditions initiales sont donc r(t = 0) = r0 et ṙ(t = 0) = u0 et on trouve
u0
r = r0 cosh(Ωt) + sinh(Ωt) (2)

v = u0 cosh(Ωt) + Ωr0 sinh(Ωt) . (3)

On va maintenant utiliser ces résultats pour calculer les termes négligés. On cherche des
ordres de grandeur donc on néglige u0 dans les expressions précédentes. Le rapport force
de Coriolis sur force centrifuge est alors

Ω2 r/(Ωv) = tanh(Ωt) . (4)

On voit qu’il faut Ωt  1 et donc r − r0  r0 . Il faut rester proche du bord de la passoire


(proche en comparaison au rayon de la passoire). Pour calculer les effects visqueux et il
faut trouver t(r) pour avoir v(r). On passe d’une description Lagrangienne à Eulerienne.
On va estimer ∆v = v 00 (r) puisque dans les fibres les vitesses sont essentiellement le long
de la fibre et notre fibre est radiale dans notre solution. On peut développer en Ωt les
solutions précédentes, mais cette fois-ci sans négliger u0 . On trouve
1
r = r0 + u0 t + Ω2 r0 t2 (5)
2
v = Ω2 r0 t + u0 (6)
Physiquement ça revient à considérer une accélération d’inertie constante et donc à rester
proche de la passoire. On aurait d’ailleurs pu faire ce choix dès le début en prenant une
accélération d’inertie constante Ω2 r0 pour trouver ce mouvement de chute libre dans la
direction radiale.
On cherche t(r) pour avoir le temps de vol de cette particule en fonction de la distance
dans la fibre p
−u0 + u20 + 2Ω2 r0 (r − r0 )
t= (7)
Ω2 r0
et donc q
v(r) = u0 1 + 2Ω2 r0 (r − r0 )/u20 (8)
En dérivant deux fois
C2
v 00 (r) = −u0 C ≡ Ω2 r0 /u20 . (9)
[1 + 2C(r − r0 )]3/2

Quand C(r−r0 )  1, donc proche de la sortie, voire à la sortie même, on trouve une accé-
lération de viscosité immense. Clairement ça bloque la solution balistique. En revanche si
C(r − r0 )  1 alors le terme de viscosité devient négligeable. Par exemple on trouve alors
µ/ρv 00 (r = 2r0 ) ' 102 m/s2 . Proche de la passoire on ne peut pas négliger la viscosité,
et l’élongation de la fibre n’est pas aussi efficace que l’approche balistique le laisse penser.

3) Une fois qu’on a v(r), alors en régime stationnaire la conservation de sucre donne que
a2 (r)v(r) est constant ou a(r) est le rayon de la fibre (il suffit de faire un bilan sur un
volume). On a donc
 −1/4
2Ω2 r0 (r − r0 )
a(r) = a0 1 + (10)
u20
Et on trouve qu’en r = 2r0 on a réduit d’un facteur (20)3/2 ' 100.

4) La résolution du problème peut se faire en prenant en compte Coriolis. Mais pour


cela on ne se place pas dans le référentiel tournant mais dans le référentiel de la fête
foraine. Après éjection à te une particule de sucre suit une trajectoire rectiligne, de
vitesse constante donc sa vitesse et position à t/te sont

~v (t) = ~vi (te ) (11)


~r(t) = ~vi (te )(t − te ) + ~ri (te ) . (12)

Les candidats peuvent être guidés en leur demandant d’exprimer les résultat en fonction
de ~ri , et ainsi les amener à écrire les deux équations précédentes.
On peut ne pas négliger la vitesse radiale initiale u0 mais rapidement c’est compliqué
même si ça reste toujours possible analytiquement. L’énoncé suggère donc de négliger u0 .
On a alors simplement la vitesse initiale reliée à la position d’éjection par
~ ∧ ~ri (te ) .
~vi (te ) = Ω (13)
En régime stationnaire la vitesse s’exprime en fonction de la vitesse dans le référentiel
tournant ~vR0 . Cette vitesse correspond à l’écoulement du fluide dans la direction de la
fibre. En effet en régime stationnaire, et dans le référentiel tournant, la fibre ne bouge
pas, mais le fluide se déplace à l’intérieur, comme dans un tuyau courbé mais fixe. La
vitesse dans le référentiel Galiléen est reliée à celle dans le référentiel tournant par
~ ∧ ~r + ~vR0 .
~v = Ω (14)

Il vient en utilisant (11) pour remplacer ~v et ~r (et le fait que ~ri est horizontal alors que
~ est vertical)

~ ∧ [Ω
~vR0 (t) = −Ω ~ ∧ ~ri (te )](t − te ) = Ω2~ri (te )τ , (15)
où on a posé le temps de vol
τ = t − te . (16)
On peut discuter ce résultat pas très intuitif. En effet, on peut être surpris que la norme
de la vitesse ne fasse qu’augmenter mais c’est essentiellement parce que c’est une vitesse
relative, et plus le temps de vol est grand, plus on est éloigné du centre de rotation et
donc plus la vitesse est grande par rapport au référentiel tournant, parce que justement
la vitesse associée à cette rotation est Ωr et r → Ωr0 τ .
La norme de (15) nous donne la vitesse d’écoulement dans la fibre. Si s est une
coordonnée curviligne (suggérer d’utiliser cette coordonnée ?) décrivant la fibre (s = 0
sur la passoire) alors pour une particule émise à te
ds 1
V ≡ |~vR0 | = = Ω2 r0 τ s = Ω2 r0 τ 2 . (17)
dt 2
On relie ainsi le temps de vol τ à la distance atteinte le long de la fibre s.
On peut ensuite inverser et trouver τ (s), c’est-à-dire le temps de vol en fonction de
la distance le long de la fibre, et on trouve
r
2s
τ (s) = , (18)
Ω2 r0
et donc ainsi avoir V (s). p
V (s) = 2Ω2 r0 s . (19)
En fait tout se passe comme si on avait une accélération constante Ω2 r0 égale à l’accé-
lération centrifuge initiale, mais dirigée toujours dans le sens de la fibre, qui elle prend
une forme courbée !

On peut ensuite trouver la forme de cette fibre. On travaille donc à un instant donné
et on veut la dépendance en s. Avant on travaillait en regardant un particule émise à te
et on regardait son évolution. Là on va prendre un t fixe, par exemple t = 0. On peut
alors relier le moment d’éjection te au temps de vol avec (16) et ainsi trouver le moment
d’éjection d’une particule à une coordonnée curviligne s :

te (s) = −τ (s) . (20)


On peut alors remplacer dans (12)
~ ∧ ~ri (−τ (s)) + ~ri (−τ (s))
~r(s) = τ (s) Ω (21)

On utilise la position d’éjection en fonction du moment d’éjection

~ri (te ) = r0 cos(Ωte )~ex + r0 sin(Ωte )~ey . (22)


p
Avec Ωte = −Ωτ (s) = − 2s/r0 , on trouve alors
h p p p i
~r(s) = r0 cos( 2s/r0 ) + 2s/r0 sin( 2s/r0 ) ~ex
h p p p i
+ r0 − sin( 2s/r0 ) + 2s/r0 cos( 2s/r0 ) ~ey (23)

ou en polaire (en reconnaissant r0 cos(b − a) et r0 sin(b − a) ci-dessus)


p
r(s) = r0 1 + 2s/r0 (24)
p p
θ(s) = arctan( 2s/r0 ) − 2s/r0 (25)

ce qui fait une jolie spirale. C’est surprenant, mais dans ce cadre d’approximation (pas
de viscosité, pas de frottement de l’air), la forme de la spirale ne dépend pas de Ω.

En fait on réalise qu’on peut le voir si on fait une construction géométrique depuis le
début. Il faut voir que les particules ont toujours la même vitesse avant et après l’éjection
(en norme). Soit elles tournent car elles ne sont pas encore éjectées, sont elles vont tout
droit une fois qu’elles ont été éjectées, mais c’est toujours avec la vitesse Ωr0 . Finalement
la forme de la courbe est la même que l’extrémité d’un fil autour d’une bobine que
l’on déroule mais en maintenant toujours la partie déroulée tendue. Cette construction
géométrique donne
~ ∧ ~ri (−τ ) + ~ri (−τ )
~r(τ ) = τ Ω (26)
où on utilise (22), c’est-à-dire

~r(τ ) = r0 [cos(Ωτ ) + Ωτ sin(Ωτ )] ~ex


+ r0 [− sin(Ωτ ) + Ωτ cos(Ωτ )] ~ey . (27)

Ensuite on fait
ds/dτ = |d~r/dτ | = r0 Ω2 τ . (28)
et on trouve s(τ ) puis τ (s) comme précédemment. On remplace alors dans (27) et on
retrouve (23). On a donc deux approches. Soit on calcule d’abord la vitesse dans le
référentiel tournant pour trouver ds/dτ et on calcule ensuite la forme de la fibre, soit on
calcule la forme de la fibre avec cet argument géométrique de la bobine qu’on déroule, et
ensuite une fois ~r(τ ) connu, on en déduit la vitesse et donc ds/dτ puis s(τ ) et enfin τ (s)
pour avoir la trajectoire en fonction de la longueur de fibre.
Effondrement gravitationnel

On considère un nuage de gaz sans pression de forme sphérique et de masse volumique


constante ρ0 .

En combien de temps s’effondre-t-il ?

Aide : On donne l’intégrale Z 1


dx π
p = (1)
0 1/x − 1 2

On considère maintenant un nuage de gaz avec pression non nulle, et remplissant tout
l’espace uniformément. On cherche les conditions pour que des perturbations de cette
configuration soient instables.

— On considère une distribution de masse volumique ρ(x, y, z). En utilisant une ana-
logie avec l’électrostatique, justifier que la force gravitationnelle exercée par cette
distribution volumique de masse sur une masse ponctuelle m que l’on rajoute peut
être mise sous la forme f~ = mG~ où divG~ = −4πGρ.

— Etablir les équations de propagation du son lorsque l’effet gravitationnel du gaz est
pris en compte.

— En considérant une onde sonore plane dans un milieu homogène et infini, trouver
à quelle condition celle-ci est instable.
Eléments de solution

Il s’agit de l’instabilité de Rayleigh-Jeans.



1) On peut d’abord faire une analyse dimensionnelle avec ρ0 et G et voir que 1/ ρ0 G est
un temps. Pour un nuage à symétrie sphérique on considère la dynamique d’un point situé
à r du centre (initialement r0 ). La force ressentie est due à la masse M (r0 ) = (4/3)πρ0 r03
et
GM (r0 )
r̈ = − (2)
r2
On multiple par ṙ pour intégrer et (avec la constante mise en prenant aucune vitesse
initiale)
1 2 GM (r0 ) GM (r0 ) GM (r0 )  r0 
ṙ = − = −1 (3)
2 r r0 r0 r
On récrit (il faut mettre un signe moins car quand le temps augmente le rayon diminue,
et à la fin on sait qu’on doit avoir un temps d’effondrement positif)
s
2GM (r0 ) dr
dt = − p . (4)
r0 r0 /r − 1

En utilisant l’aide on a donc pour l’effondrement total (r va de r = r0 à r = 0 et on


a juste un petit changement de variable x = r/r0 à faire et inverser les bornes)
s
2GM (r0 ) π
3 teff = , (5)
r0 2

et donc r r
1 3π 1
teff = . (6)
4 2 Gρ0
L’analyse dimensionnelle donne le bon ordre de grandeur et là on a le résultat complet.

2) la force est en 1/r2 et le potentiel en 1/r dans les deux cas. Pour l’électrostatique la
constante est 1/(4π0 ) et pour la gravité c’est G. Pour l’électrostatique c’est la distribu-
tion de charge électrique et pour la gravité la distribution de masses. Mais on a aussi un
signe différent car deux charges électriques de même signe se repoussent, alors que pour
deux masses c’est le contraire. Donc si en électrostatique on a

~ = 4π ρelec
divE (7)
4π0
alors avec cette analogie on a
~ = −4πGρmass .
divG (8)
3) On refait la dérivation habituelle de l’équation de propagation du son pour ρ = ρ0 + δρ
mais dans l’équation d’Euler on a une force volumique en plus de −gradP ~ ~ A
qui est ρG.
l’ordre le plus bas on ne peut avoir G~ car le gaz est homogène et il n’y a aucune direction
préférée. Donc on a ce terme qui apparait uniquement lorsqu’on considère la linéarisation
et divG~ = 4πGδρ. La linéarisation d’Euler donne

ρ0 ∂t~v = −c2 gradδρ + ρ0 G~ . (9)

Comme d’habitude on prend la divergence (on utilise (8)) et on remplace dans la dérivée
de l’équation de continuité pour trouver

∂tt δρ − c2 ∆δρ = −divρ0 G~ = 4πGρ0 δρ . (10)

4) On considère une onde plane selon x de mode k, donc en ∝ ei(ωt−kx) avec ω = kc et il


vient la relation de dispersion

ω 2 = c2 k 2 − 4πGρ0 . (11)

Si ρ0 est assez grand, ou plutôt si k est suffisamment petit, on va avoir ω 2 < 0 et donc
un mode croissant. La longueur d’onde de Jeans, au dessus laquelle on a instabilité est
telle que
 2
2 2π
c = 4πGρ0 . (12)
λ
On peut estimer pour l’atmosphère. On peut calculer la vitesse du son avec c2 = γ RT M
où M est la masse molaire et donc c2 ' (340m/s)2 . La masse volumique de l’air est
1.2kg/m3 , et donc λ ' 67000 km. L’atmosphère est stable à Rayleigh-Jeans. Pour l’ef-
fondrement des étoiles c’est beaucoup plus grand car ρ0 est beaucoup plus faible dans les
nuages moléculaires.

Comment interpréter la compétition√ entre les deux termes de (10) ? On a compétition


entre le temps d’effondrement en 1/ Gρ0 et le temps correspondant à la propagation des
ondes sur une longueur d’onde 1/(kc). Si le son est beaucoup plus rapide, il va effacer
les surdensités. S’il est plus lent que l’effondrement alors il n’arrive pas à contrebalancer
à temps l’effondrement.
Oscillations d’un piston dans un tube hemi-circulaire


On considère un tube de section circulaire s, coudé en demi cercle de rayon a (a  s), rempli
d’un gaz parfait. Le tube est fermé aux deux extrémités. Un piston de masse m et d’épaisseur
négligeable devant a peut coulisser sans frottement à l’intérieur du tube. La température du
gaz est maintenue constante à T . De chaque côté du piston on a la même quantité n de gaz.

a) En fonction des différentes données, déterminer la ou les positions d’équilibre du piston. En


particulier, on identifiera un régime de basse température et un régime de haute température
et on étudiera plus en détail la transition entre les deux régimes.

b) Déterminer la fréquence des petites oscillations autour de ces positions d’équilibre. Discuter
la validité de l’hypothèse d’une température constante. Par quoi peut elle être remplacée et
comment cela changerait les fréquences d’oscillation?

1
On répère la position du piston par l’angle φ compté depuis le haut du tube (φ = 0 est le
nRT nRT
point symétrique). Les pressions de chaque côté du piston sont alors pg = Vg
= as( π2 +φ)
et
nRT nRT
pd = Vd
= as( π2 −φ)
. Le PFD du piston s’écrit alors

nRT  1 1 
maφ̈ = mg sin φ + π − π . (1)
a 2
+φ 2
−φ

a) Les positions d’équilibre sont données par φ̈ = 0, soit


 1 1  nRT
sin φ = C π − π ≡ f (φ) , C= (2)
2
−φ 2
+φ mga

On peut essayer une résolution graphique. La courbe orange est sin φ, la verte est est f (φ) pour
C = 2 et la bleue est f (φ) pour C = 0.3. On voit que pour C grand (haute température) il n’y

-1.5 -1.0 -0.5 0.5 1.0 1.5

-1

-2

-3

a qu’une solution φ = φ0 = 0, mais que pour C petit (basse température) on a deux solutions
φ = ±φ1 avec f1 6= 0. Plus précisément, tant que C > C∗ il n’y aura que la solution φ0 = 0,
mais pour C < C∗ on aura en plus les deux solutions ±f1 avec f1 qui augmente au fur et à
mesure que C∗ − C diminue. La transition de une à trois solutions a lieu lorsque la dérivée en
φ = 0 de f (φ) est égale à la dérivée de sin φ, c’est à dire égal à 1. On a

8C π2 π 2 mga
f 0 (0) = ⇒ C∗ = ⇔ T∗ = . (3)
π2 8 8nR
Physiquement, dans la phase haute température T > T∗ , les pressions sont plus importantes
que le poids et le piston est maintenu en haut. Dans la phase basse température T < T∗ , le
poids l’emporte et le piston se stabilisera à droite ou à gauche quand le gaz sera suffisamment
comprimé pour que la pression égale de nouveau le poids. C’est donc φ0 qui devient instable
et les solutions ±φ1 stables.
π
Clairement, pour C → 0 on a φ1 → 2
. Pour C légèrement inférieure à C∗ on est dans
la région de la “transition de phase” où φ1 est petit et on peut développer f (φ) autour de
φ = 0. Comme sin φ, cette fonction est impaire et il faut pousser le DL jusqu’à l’ordre 3:

2
f (φ) = C( π82 φ + 32 3
π4
φ) + O(φ5 ) = C
C∗
φ + C
2C∗2
φ3 + O(φ5 ). L’équation (2) s’écrit alors (sin φ =
φ3
φ− 6
+ O(φ5 ))
h C C  2i

φ (C∗ − C) − + φ + O(φ5 ) = 0 . (4)
6 2C∗
On voit de nouveau que pour C > C∗ , seul φ = 0 est solution, alors que pour C < C∗ on a 3
solutions. Dans ce dernier cas, les solutions stables seront données par
s r r
6C∗ (C∗ − C) 3 1/2 3
±φ1 = ± 2
'± (C∗ − C) ' ± (T∗ − T )1/2 . (5)
C∗ + 3C) 2C∗ 2T∗
(Remarque : On retrouve l’exposant critique 1/2 caractéristique d’une transition de phase avec
brisure spontanée d’une symétrie Z2 .)

b) Le PDF (1) peut être réécrit comme


a
φ̈ = sin φ − f (φ) . (6)
g
Dans la phase de haute température il convient de développer autour de la solution stable
φ = φ0 = 0 et on trouve
a C C∗ − C T∗ − T g T − T∗
φ̈ = φ − φ= φ= φ ⇒ φ̈ + ω 2 φ = 0 , ω2 = . (7)
g C∗ C∗ T∗ a T∗
Pour T → T∗ on a ω 2 → 0 et lorsqu’on passe dans la phase basse température on aurait ω 2 < 0
ce qui est le signe que φ = 0 est devenu instable.
Dans la phase de basse température il faut développer autour de φ1 (ou −φ1 ). Soit φ =
φ1 + δφ et rappelons que sin φ1 − f (φ1 ) = 0. On trouve alors (cos δφ ' 1 et sin δφ ' δφ)
a¨    
δφ = cos φ1 − f 0 (φ1 ) δφ ⇒ δφ ¨ + ω 2 δφ = 0 , ω 2 = g cos φ1 − f 0 (φ1 ) . (8)
g a
Pour T proche de T∗ , mais toujours inférieur à T∗ , on sait que φ1 est petit (cf eq (5)) et on peut
approximer
g g C  g T∗ − T
ω 2 ' (1 − f 0 (0)) = 1 − = . (9)
a a C∗ a T∗
Proche de la température critique T∗ la fréquence des oscillations est donc faible et l’hypothèse
d’une température maintenue constante est réaliste. Par contre lorsque la fréquence est plus
élevée (amplitude fois fréquence n’est plus petit devant la vitesse du son dans le gaz) une
hypothése adiabatique parait plus raisonnable et il faudrait alors remplacer pV = nRT par
pV γ = p0 V0γ . Le PFD deviendrait alors
1 1 
maφ̈ = mg sin φ + p0 s π − π . (10)
( 2 + φ)γ ( 2 − φ)γ
ou
a po s  1 1 
φ̈ = sin φ − fe(φ) , fe(φ) = − . (11)
g mg ( π2 − φ)γ ( π2 + φ)γ
Le reste de l’analyse est similaire avec fe à la place de f .

3
Fluide parfait

On considère un fluide parfait de masse volumique ρ(t, x, y, z). On note le champ des vitesses
dans le fluide: ~v (t, x, y, z). On néglige l’effet de la gravité.

1) Montrer que si l’énergie interne de ce fluide est constante, donc son énergie volumique
1
totale (énergie cinétique plus énergie interne) est de la forme 2
ρv 2 à une constante près,
alors: div(~v ) = 0.

2) On ne suppose plus l’énergie interne du fluide constante, montrer alors que nécessairement
l’entropie (par unité de masse) s vérifie: ds/dt = 0, où la dérivée par rapport au temps est
faite en suivant le mouvement (tout mouvement dans le fluide est isentropique).

On rappelle que l’énergie interne d’un fluide vérifie la relation: dU = T dS − pdV .


L’exercice demande de manipuler l’équation d’Euler. Il faut donc être au point sur le cours
des fluides parfaits (équation d’Euler et équation de conservation de la masse) et avoir utilisé
ces relations dans des exercices pour s’y habituer. La question 1) est simple, proche du cours.
La question 2) est plus technique.

1) L’énoncé suggère une approche énergétique. On considère un volume fixe V de fluide et


R
on écrit le bilan énergétique pour ce volume. Soit l’énergie du fluide dans V : V 21 ρv 2 dV plus
une constante mais qui n’a pas d’importance car on va regarder sa variation dans le temps.
R
La quantité dtd ( V 12 ρv 2 dV ), peut s’écrire de différentes manières. Déjà, comme le volume est
fixe: Z Z
d 1 2 ∂ 1 2
( ρv dV ) = ( ρv )dV. (1)
dt V 2 V ∂t 2

Et cela sécrit aussi en fonction du flux entrant ou sortant (d’énergie volumique) plus la
puissance des forces de pression à la surface S du volume V :
Z Z Z
d 1 2 1 ~ − p~v dS.
~
( ρv dV ) = − ( ρv 2 )~v dS (2)
dt V 2 S 2 S

C’est l’égalité entre (1) et (2) en utilisant l’équation d’Euler qui va donner le résultat. Comme
l’égalité entre (1) et (2) est vraie pour tout volume V , on obtient:

∂ 1 2 1
( ρv ) + div(( ρv 2 )~v ) + div(p~v ) = 0.
∂t 2 2
Soit:
1 ∂~v 1 1 ~
(−div(ρ~v ))v 2 + ρ~v + (div(ρ~v ))v 2 + ρ~v .grad(v 2
) + div(p~v ) = 0.
2 ∂t 2 2
Ce qui donne, en utilisant l’équation d’Euler (projetée sur ~v ), div(~v ) = 0. Ce qui conclut.

2) On commence comme pour la question 1) sauf que dans (1) et (2), il faut remplacer 12 ρv 2
par 21 ρv 2 + ρ, où ρ est l’énergie interne volumique et donc  l’énergie interne par unité de
masse. En reprenant le raisonnement de 1), on obtient:

∂ρ
+ div(ρ~v ) + pdiv(~v ) = 0.
∂t
Il faut montrer que cette relation est équivalente à la conservation de l’entropie (massique)
en suivant le mouvement. On l’écrit sous la forme:
d(ρ)
+ ρdiv(~v ) + pdiv(~v ) = 0. (3)
dt
De plus: d = T ds + pdρ/ρ2 . On rappelle que  est l’énergie interne massique (ou par unité
de masse) du fluide. Donc, d(ρ)/dt = dρ/dt + ρT ds/dt + p/ρ dρ/dt. Alors (3) devient:

dρ/dt + ρT ds/dt + p/ρ dρ/dt + ρdiv(~v ) + pdiv(~v ) = 0.

~
On écrit alors les termes en dρ/dt comme ∂ρ/∂t + ~v .grad(ρ) et, après simplifications, on
obtient: ds/dt = 0, ce qui conclut.
Dissolution du sel dans l’eau

Un grand récipient sphérique de rayon R = 10 m est rempli d’eau. A t = 0 on introduit au


centre du récipient un cristal d’un certain sel, faiblement soluble dans l’eau. Pour simplifier,
on supposera ce cristal de forme sphérique, de rayon a0 = 1 cm. Déterminer comment la
concentration en sel dans l’eau varie en fonction du temps et de la position. On établira d’abord
les équations, puis on tentera de les résoudra en procédant à des approximations raisonnables.

1
L’énoncé est volontairement un peu vague. On ne sait pas de quel sel il s’agit à part qu’il est
faiblement soluble dans l’eau. En se basant sur une masse volumique de l’ordre de 2 g cm−3
on a environ 8 g de sel pour environ 4 × 106 litres d’eau, soit 2 × 10−6 g l−1 . La concentration
saturante de NaCl est environ 360 g par litre d’eau! Pour un sel très faiblement soluble tel
AgCl, la solubilité est de 2 × 10−3 g l−1 . Donc même pour un sel faiblement soluble, on pourra
supposer qu’à la fin tout le sel est dissous.
On a une symétrie sphérique. (C’est bien pour cela qu’on a pris un cristal de forme
sphérique!) On aura donc des concentrations n(t, r) et des flux ~j(t, ~r) = j(t, r)~ur . On aura
d’une part un processus de diffusion et d’autre part une source au centre qui injecte du sel dans
l’eau au fur et à mesure que le cristal se dissout. Le flux à la surface du cristal (nombre de
grammes dissous par unité de temps et de surface) dépendra de la différence entre la concen-
tration en sel dissous dans l’eau à la surface du cristal n(t, a) et de la concentration saturante
nsat :

~j(t, ~r)|r=a = γ nsat − n(t, a) ~ur (1)

où γ est une certaine constante dont l’unité est m s−1 . Mais ici a = a(t) ≤ a0 est le rayon du
cristal à l’instant t. En intégrant sur la surface du cristal on obtient la quantité (masse) de sel
injecté dans l’eau par unité de temps au centre du récipient:
Z

ṁ(t) = − dS~ · ~j(t, ~r)|r=a(t) = −4πa(t)2 γ nsat − n(t, a(t)) . (2)
r=a(t)

~ · ~j = 0 et la loi de Fick ~j = −D∇n,


Pour r > a(t) on a l’équation de conservation ṅ + ∇ ~ soit

∂  2 ∂n 
ṅ = D∆n ⇒ ṅ = D r . (3)
∂r ∂r
(Le laplacien en coordonées sphérique n’est pas exigible, mais peut être obtenu en utilisant le
~ · ~j = 0 avec une couche entre r et r + dr.) L’équation (1) donne
théorème de Gauss sur ṅ + ∇
~ donc la condition de bord
~j en r = a et par continuité cela donne D∇n,

 ∂n γ 
~
D∇n(t, a) = γ n(t, a) − nsat ~ur ⇒ (t, a) = n(t, a) − nsat , a ≡ a(t) , (4)
∂r D
ainsi que
n(t, r) = 0 , t ≤ 0 , r > a0 . (5)

Finalement il faut tenir compte comment a(t) diminue. Soit ρ la masse volumique du cristal
4π 3
de sel. Alors 3
a (t)ρ = m(t), soit en prenant la dérivée et en utilisant (2) :
γ 
ȧ(t) = − nsat − n(t, a(t)) . (6)
ρ

2
On a un système d’équations couplées compliqué, à savoir (3) avec conditions aux bords
et initiales (4) et (5), avec a(t) déterminé implicitement par (6). Afin de pouvoir avancer, on
supposera par exemple qu’on ne s’intéresse qu’à des temps suffisamment courts pour que le
cristal ne change pas considérablement de taille de façon à approximer

a(t) ' a0 . (7)

Ensuite, à la surface du cristal il y a une compétition entre l’augmentation de n par la dissolution


du cristal (γ) et la diffusion (D) qui va diminuer n. Puisqu’on est dans le cas d’un sel faiblement
soluble on pourra supposer que la diffusion évacue efficacement le sel vers les r croissant et qu’il
n’y a pas d’accumulation de sel dissous au voisinage du cristal. Plus précisément on supposera
D ∂n
que aγ
 1. Cela implique que, très vite pour t ≥ 0, on a n(t, a) constant et donc ∂r
(t, a)
constant. Mais ce régime stationnaire s’établit aussi pour des r > a jusqu’à un certain rmax (t)
qui constitue une sorte de front de la propagation de la concentration du sel. Donc pour
r < rmax l’éq. (3) donne (r2 n0 (r))0 = 0 soit

B
n(r) = , a0 < r < rmax (t) . (8)
r
γ

La condition de bord (4) donne alors − aB2 = D
B
a0
− nsat , soit
0

a0 nsat nsat a0
B= ⇒ n(r) = , a0 < r < rmax (t) . (9)
1 + aD0 γ 1 + aD0 γ r

On peut estimer rmax en calculant la masse totale du sel dissous à l’intérieur de rmax :
Z rmax
2 nsat a0  rmax
2
a20 
m(r < rmax ) = 4π dr r n(r) = 4π − , (10)
a0 1 + aD0 γ 2 2

ce qui doit être égal, d’après (2), à


 a0 Dnsat
4πa20 γ nsat − n(a0 )) × t = 4π ×t , (11)
1 + aD0 γ

soit
2
rmax (t) = a20 + 2Dt . (12)

3
Oscillateur sur corde vibrante

On considère une corde de longueur 2L, de masse linéique λ et tendue avec une tension
T entre deux parois solides. Au milieu de cette corde est attaché un ressort de raideur K
et au bout de ce ressort est fixée une masse m qui peut se mouvoir selon une direction
orthogonale à la direction de la corde au repos (voir schéma). Il n’y a pas de gravité.

— A quelles fréquences le système peut-il osciller en régime stationnaire ?


Eléments de solution

On note z(x) la position de la corde et on prend l’origine des x au milieu si bien que
z(−L) = z(L) = 0 par fixation aux parois. On note Z la coordonnée de la masse par
rapport à la position à vide (dans la configuration ou la corde aussi est au repos). Je
note z0 la position du point central de la corde, là où le ressort est attaché. L’équation
du mouvement de Z est
K
Z̈ + (Z − z0 ) = 0 . (1)
m
On note la fréquence propre du système ressort et masse ω02 = K/m. En régime oscil-
latoire stationnaire de fréquence ω on va avoir (bien sûr on peut détailler ou alors se
souvenir du cours de sup d’électrocinétique et aller vite)

ω02 ω2
Z= z0 (z0 − Z) = − z0 . (2)
ω02 − ω2 ω02 − ω2

En régime stationnaire
p on doit avoir des ondes stationnaires dans la corde. La vitesse des
ondes est c = T /µ et un mode d’onde est relié à la fréquence par k = ω/c. Une onde
stationnaire à gauche est (avec la condition z(−L) = 0)

zgauche (x, t) = A sin(k(L + x))eiωt (3)

et le symétrique à droite est

zdroite (x, t) = B sin(k(L − x))eiωt . (4)

On a A = B car ces fonctions doivent être égale en x = 0 (sauf si kL = nπ, cas envisagé
brièvement à part ci-après).
On fait une petite pause pour vérifier ce qui se passe quand on n’accroche pas un
ressort. Dans ce cas, la corde doit être continue et les dérivées doivent égales. Comme en
x = 0 elles valent

∂x zgauche (x = 0, t) = Ak cos(kL)eiωt (5)


iωt
∂x zdroite (x = 0, t) = −Bk cos(kL)e (6)

on en déduirait sans ressort ni masse que kL = π/2 + nπ car A = B (sauf le cas où on


avait pris kL = nπ plus haut et alors on a besoin de A = −B. On voit que kL = nπ est
un peu spécial car on a un noeud au milieu donc on peut bien accrocher un ressort si on
veut, la corde ne le fera pas bouger.

Mais on a un ressort. Donc au point d’attache, on va avoir un discontinuité de tangente


en général résultant en une force nette due à la tension et qui est équilibrée par la tension
du ressort (le point d’attache est sans masse). On suppose qu’on n’a pas kL = nπ et donc
pas de noeud au point d’attache (ou plutôt on a déjà identifié que ce sont des oscillations
possibles) et ainsi A = B comme vu plus haut. Dans la corde, la force de la droite sur la
gauche du fil étant ' T ∂x z, le bilan sur le point d’attache donne
T [∂x zdroite (x = 0, t) − ∂x zgauche (x = 0, t)] = K(z0 − Z) . (7)
soit  
ω2
2kT A cos(kL) = K A sin(kL) (8)
ω02 − ω 2
et on peut simplifier par A. Cette condition ne dépend pas de l’amplitude. On utilise
k = ω/c et on définit ωc = c/L. Cela se récrit
 ωcm
ω02 − ω 2 cos(ω/ωc ) = ω02 sin(ω/ωc ) (9)
2T
ou encore 
m ω
ω02 tan(ω/ωc ) = ω02 − ω 2 . (10)
2µL ωc
On voit en traçant la gauche et la droite qu’il va y avoir plein de solutions possibles. Si
K = 0 alors ω0 = 0 alors il faut ω = (π/2 + nπ)ωc ou kL = (π/2 + nπ). On retrouve que
les oscillations habituelles avec un ventre au milieu vont bien. Le cas qui est intéressant
à discuter c’est quand ω0  ωc mais pas nul. Alors, on peut développer la tangente en
regardant une solution proche de ω0 . On trouve
m ω2 ω2
= 1 − (11)
2µL ωc2 ω02
et donc
1 1 1 1
2
= 2 + . (12)
mω mω0 2µL ωc2
On se souvient que deux ressorts en série, ont les inverses de leur raideurs qui s’ajoutent.
1
Ici tout se passe comme si on avait ajouté l’inverse de la raideur du ressort 1/K = mω 2
0
1 1
avec l’inverse de la raideur de la corde 2µL ωc2
se comportant comme un ressort car ne
développant pas plein de noeuds. On note que 2µL est la masse de la corde. Bien sur
comme on est dans l’hypothèse ω0  ωc , on n’a qu’une faible modification et ω est
proche de ω0 . Le point d’attache du ressort n’étant pas stable on a une légère baisse
de la fréquence de régime stationnaire. On voit aussi sur (2) que l’amplitude du point
d’attache est beaucoup plus petite que l’amplitude de l’allongement du ressort dans ce
cas.

On peut aussi questionner le cas ω0  ω. On a alors approximativement


2µL ωc
tan(ω/ωc ) = (13)
m ω
On va avoir plein de solutions qui vont tendre vers ω → nπωc ou kL → nπ mais en
approchant par le supérieur. On a quasiment un noeud au point d’attache. La masse est
petite et elle fait une petite cassure de pente au point d’attache mais rien de méchant,
et la raideur du ressort va faire que l’allongement du ressort reste toujours très petit. On
doit retrouver les oscillations stationnaires d’une corde avec une perle massive accrochée
au milieu.
Stabilité d’un fluide 2D

On considère un fluide idéal incompressible (de masse volumique ρ constante) à 2 dimensions


qui occupe un grand espace Σ (dans le plan xy). On suppose que son champ des vitesses
~v (x, y, t) est régulier à l’instant initial t = t0 , c’est-à-dire que ni ~v (.), ni aucune de ses dérivées
spatiales ne divergent en aucun point de Σ. De plus, on suppose que la condition aux limites
suivante est toujours vraie: le champ des vitesses et ses dérivées spatiales tendent vers 0 aux
bornes de Σ.

La condition sur la réguarité de ~v (.) à t = t0 signifie en particulier:


Z
|~v (x, y, t = t0 )|2 dS < ∞ (1)
Σ
Z h i
~
|grad(v 2 ~
x )| + |grad(vy )|
2
(x, y, t = t0 )dS < ∞ (2)
Σ

Montrer qu’alors (1) est vraie pour tout t ≥ t0 .

Si besoin, on utilisera (avec ω ~ v ):


~ = rot~

~ v = 1 grad(v
(~v .grad)~ ~ 2
~ ∧ ~v .
)+ω
2

Note: on pourrait montrer également que (2) est vraie pour tout t ≥ t0 et conclure que la régularité
du champ des vitesses du fluide est assurée à long terme si elle est vraie à l’instant intitial: ni ~v (.),
ni aucune de ses dérivées spatiales ne divergent en aucun point de Σ quel que soit t.
Il y a une manière ’force brute’ de faire l’exercice, on va commencer par cette approche. C’est
la plus naturelle: il faut manipuler l’équation d’Euler et la relation de conservation de la
masse (dans le cas incompressible, div ~v = 0). La difficulté de l’exercice est qu’il ne faut
pas se perdre dans les calculs et penser à intégrer par partie des intégrales autant de fois que
nécessaire.

Montrer la validité de (1) quel que soit t est très simple. On montre que (Σ est fixe en fonction
du temps par définition):

Z Z
d 2 ∂
|~v (x, y, t)| dS = 2 ~v (x, y, t) ~v (x, y, t)dS = 0 (3)
dt Σ Σ ∂t


Mais ∂t
~v (x, y, t) s’obtient par l’équation d’Euler:
∂ ~ v − grad(p/ρ)
~
~v = −(~v .grad)~
∂t

Alors, l’intégrale dans (3) devient:


Z Z h i
d 2 ~ v ] + ~v (x, y, t).grad(p/ρ)
~
|~v (x, y, t)| dS = −2 ~v (x, y, t)[(~v .grad)~ dS
dt Σ Σ

Le deuxième terme de l’intégrale ci-dessus donne 0 après une intégration par partie, en utilisant
div ~v = 0, et le fait que les vitesses sont nulles aux bornes du domaine Σ. Pour le premier
~ v ] comme:
terme de l’intégrale, on peut par exemple écrire ~v (x, y, t)[(~v .grad)~

~ v ] = ~v (x, y, t)[ 1 grad(v


~v (x, y, t)[(~v .grad)~ ~ 2 1 ~
~ ∧ ~v ] = ~v (x, y, t)grad(v
)+ω 2 1
) = div(v 2~v )
2 2 2
Donc, ce terme s’intègre simplement et donne 0 du fait que les vitesses sont nulles aux bornes
du domaine. Au final, on a bien trouvé que la relation (1) est vraie pour tout t.

Montrer que la relation (2) est vraie pour tout t, donc:


Z h i
d ~ 2 ~ 2
|grad(vx )| + |grad(vy )| (x, y, t)dS = 0
dt Σ
demande plus de calculs. Mais sur le principe, c’est exactement comme plus haut: une suite
d’intégrations par partie avec l’utilisation de l’équation d’Euler et de la conservation de la
masse.

Pour montrer que (2) est vraie pour tout t, on peut aussi utiliser une astuce. (2) traduit
uniquement le fait que les dérivées spatiales du champ des vitesses sont finies à l’instant
initial en tout point de Σ. On peut traduire cette propriété différemment: (2) est équivalente
au fait que la vorticité est finie l’instant initial en tout point, avec en plus la condition sur la
conservation de la masse. Alors, les 4 termes: ∂x vx , ∂y vx , ∂x vy , ∂y vy doivent être finis en tout
point de Σ (à l’instant initial). Mais, div ~v = 0 est vraie pour tout t et de plus (en 2 D):

d
ω = 0.
dt
Cela veut dire que si la vorticité est finie à l’instant initial, alors, elle le reste pour tout t. Ce
qui conclut.
Télésiège

On considère un télésiège initialement au repos et qui se met en mouvement à t = 0. On


veut étudier la dynamique de l’un de ses sièges. On modélise le cable sur lequel est fixé
ce siège par une très longue tige solide mais de masse négligeable, et inclinée par rapport
à l’horizontale d’un angle α. Cette tige ne peut se déplacer que dans sa direction et on
note V la norme de sa vitesse. Le siège est modélisé par un pendule (fil de longueur ` et
masse m) fixé à cette tige.

— Déterminer les équations du mouvement en fonction de l’accélération du cable V̇ .

— On suppose que à t = 0− le siège est au repos et le cable ne bouge pas. Un moteur


actionne rapidement le cable afin de lui faire acquérir une vitesse montante V0
à t = 0+ et qui restera constante par la suite. Quelle est la vitesse du siège à
t = 0+ ? Quel est le travail fourni par ce moteur entre 0− et 0+ ?

— Quelle est l’évolution ultérieure de la position du siège ?

— Quelle est la puissance fournie par le moteur aux temps t > 0+ ? Quelle est sa
moyenne ?

— Si le siège est modélisé par une tige solide pouvant tourner librement autour de
son point de fixation, comment la modélisation est-elle modifiée ?
Eléments de solution

On peut traiter le problème dans le référentiel attaché au point de fixation O qui n’est
pas Galiléen entre t = 0− et t = 0+ à cause de son accélération brusque et on utilisera
alors les forces d’inertie pour une translation rectiligne non-uniforme.

1) Dans le référentiel attaché au point de fixation O (et défini avec ses axes en translation
par rapport au référentiel de la piste noire), l’accélération d’inertie est juste l’accélération
de ce point O notée ~aO (et sa norme V̇ ). On note ~vM/O la vitesse du point matériel M (le
bout du pendule) par rapport à O et on utilise une notation similaire pour l’accélération.
On a les lois de composition

~vM = ~vM/O + ~vO (1)


~aM = ~aM/O + ~aO , (2)

où à gauche ce sont les vitesses et accélérations dans le référentiel de la station. Le PFD


donne (on peut parler de forces d’inertie si on préfère)

T~
~aM/O = + ~g − ~aO . (3)
m
~ = r~er , on a simplement
En prenant des coordonnées polaires pour OM

`θ̈ = −g sin θ + V̇ cos(α + θ) (4)


T
−`θ̇2 = g cos θ + V̇ sin(α + θ) − (5)
m
où T est la tension dans le fil.

2) On va intégrer (4) entre 0− et 0+ . L’accélération ~aO est nulle sauf entre 0− et 0+ et


on a alors Z 0+
V0 = V̇ dt . (6)
0−
Pendant tout ce temps infinitésimal, rien n’a le temps de bouger et seules les vitesses
changent. On a donc θ = 0 si bien que l’intégration donne


` θ̇ + = V0 cos α . (7)
t=0

On se pose la question de la vitesse de la masse au bout du pendule car on va en


avoir besoin ensuite pour résoudre. En effet ça va nous servir de condition initiale. La
composition (1) donne
~vM = ~vO + `θ̇~eθ . (8)
Or la vitesse du cable est

~vO |t=0+ = V0 [− cos(α + θ)~eθ − sin(α + θ)~er ] (9)


donc à l’instant initial on a (en définissant ~ez associé à l’axe vertical)

~vM |t=0+ = V0 sin α~ez . (10)

On a encore utilisé que de 0− à 0+ , θ = 0. On va donc pouvoir connaître l’énergie


cinétique à 0+ et le travail fourni est donc
1
W = mV02 sin2 α . (11)
2
On remarque que si α = 0, on n’a pas mis en mouvement le pendule donc aucune énergie
cinétique et aucun travail. La rotation impulsée est juste due à la vitesse initiale du point
d’attache sans que la masse du pendule ne bouge. En revanche pour α 6= 0 on a impulsé
une vitesse verticale à la masse du pendule. On pourrait retrouver le travail fourni en
utilisant que
P = T~ · ~vM = T |~vO | sin(α + θ) (12)
et en intégrant de 0− à 0+ une fois exprimée la tension T avec (5). On voit alors que
le terme venant de la gravité est un infinitésimal, tout comme celui venant de `θ̇2 et on
trouve alors (toujours en prenant θ = 0)
Z 0+ Z 0+
1
W = Pdt = m~aO · ~vO sin2 αdt = mV02 sin2 α (13)
0− 0− 2

Mais c’est quand même moins facile.

3) Maintenant que l’on connait θ̇ à t = 0+ on a juste le mouvement libre d’un pendule


dans un référentiel Galiléen car le cable a une vitesse constante. Pour des petits angles
on a juste r 
V0 cos(α) g
θ= √ sin t (14)
g` `
qui satisfait bien la condition initiale (7) et θ(t = 0+ ) = 0.

4) Pour la puissance fournie il faut utiliser (12). Ca s’obtient en regardant le travail


fourni par la tension sur la masse (cf. (12)), ou alors de manière équivalente au travail
fourni par les forces de contact sur le point d’attache, car seule la composante selon le
cable de ces forces de contact joue. Il faut donc déterminer T . On peut intégrer (4) (en
ayant multiplié par θ̇) ou regarder la conservation de l’énergie cinétique du pendule (c’est
pareil) pour trouver

`θ̇2 = 2g(cos θ − 1) + ` θ̇2 + (15)
t=0
et donc remplacer dans (5)

mV02 cos2 α
T = gm(3 cos θ − 2) + . (16)
`
En particulier initialement on a T > mg. On est écrasé au fond du siège au redémarrage.
Sur un vrai télésiège le cable n’est pas un solide, et donc ça fait un petit mouvement
d’oscillation vers le bas du point d’attache car on augmente la composante normale au fil
T cos(α + θ) = T cos α et donc la forme du fil (sa discontinuité de tangente et ensuite ça
oscille de bas en haut et c’est assez désagréable) doit s’ajuster à cette nouvelle contrainte
normale. La puissance fournie par le moteur est alors donnée par (12) donc

P = T V0 sin(α + θ) . (17)

Pour calculer sa moyenne c’est un peu traitre. En effet il faut développer mais à l’ordre
2 dans les petits angles. On trouve pour l’expression ci-dessus

sin(α + θ) ' sin α(1 − θ2 /2) + θ cos α (18)

et pour la tension  
3 V 2 cos2 α
T ' mg 1 − θ2 + 0 . (19)
2 g`
Ensuite on utilise que hθi = 0 et hθ2 i = θmax
2 /2 pour trouver

hPi = mgV0 sin α . (20)

La puissance moyenne fournie est évidemment celle nécessaire pour augmenter l’énergie
potentielle moyenne.

5) On peut refaire avec une tige avec un moment d’inertie I = m`2 /3 le même exercice.
Le théorème du moment cinétique donne
 
` `
I θ̈ = m − g sin θ + V̇ cos(α + θ) . (21)
2 2

On fait pareil que pour le pendule, on intègre de 0− à 0+ en prenant θ = 0 pour avoir la


vitesse de rotation à t = 0+ . On trouve
3

` θ̇ + = V0 cos α . (22)
t=0 2
On peut calculer la vitesse du point bas du pendule en utilisant (8) et on trouve alors
1
~vM = V0 cos α~eθ . (23)
2
Dans le cas simple α = 0, on voit que cette impulsion donnée vers la gauche fait déplacer
le point bas vers la droite. Quand le télésiège redémarre, on se sent donc reculer un peu,
en plus du tassement au fond du siège. On peut faire la même chose pour le barycentre
et on voit qu’il avance.
On peut aussi recalculer le travail fourni par le moteur de 0− à 0+ . On peut simple-
ment faire une intégrale le long de la tige en sommant les énergies cinétiques. En effet la
composante horizontale de la vitesse pour un point de distance r au point d’attache, est
 
3r
V0 cos α −1 (24)
2`

et la composante verticale est la même que dans le cas du pendule (10). On intègre les
énergies cinétiques des éléments dm = m/`dr de 0 à ` et on trouve
1 1
W = mV02 cos2 α + mV02 sin2 α . (25)
8 2
En plus de ce qu’il faut apporter pour imprimer un mouvement vertical, il faut apporter
ce qu’il faut pour la mise en rotation, ce que l’on n’avait pas dans le cas du pendule avec
masse ponctuelle. On aurait aussi pu faire un calcul direct en obtenant l’expression de
l’énergie mécanique (cinétique et potentielle), et on calculant sa dérivée temporelle.
Corde sur une table

Une corde de longueur L (de masse linéique σ) est enroulée sur elle-même pour former un tas
sur une table, excepté pour une toute petite partie de la corde qui pend en passant par un
trou fait dans la table. La corde est alors laissée libre et commence à glisser à travers le trou.
Le mouvement de la corde se fait sans frottement. Dans ce montage, on admet que seule la
partie de la corde qui arrive au niveau du trou et en dessous est en mouvement. La partie de
la corde qui est posée en tas sur la table ne bouge pas (sauf quand elle arrive au niveau du
trou).

Déterminer l’énergie cinétique de la corde lorsqu’elle quitte complètement la table (toute la


corde a passé le trou). En déduire qu’il y a une quantité de chaleur Q > 0 dissipée lorsque la
corde quitte complètement la table dans ce montage.
Le première partie de la question consiste à faire un bilan de quantité de mouvement pour un
système mécanique de masse variable. C’est donc proche du cours qui traite le cas de la fusée
ou d’une avalanche. Il y a une petite variation par rapport au cours car, ici, la résolution de
l’équation du mouvement demande de réfléchir un peu. La deuxième partie utilise ce résultat
pour une observation un peu subtile (comme souvent en mécanique).

On note x = x(t) la partie de la corde qui pend en-dessous du trou. Sa masse est donc σx et
l’équation du mouvement pour ce morceau de corde s’écrit:
d(σxẋ)
= σxg.
dt
Soit:
xẍ + ẋ2 = xg.
Pour trouver l’énergie cinétique de la corde lorsqu’elle a complètement passé le trou, il faut
obtenir le temps T tel que: x(T ) = L puis la vitesse ẋ pour t = T . Il faut donc résoudre
l’équation ci-dessus. C’est une variation par rapport au cours car, ici, il faut garder un sens
physique pour trouver simplement la solution. On sait que le cas physique de ce problème est
un mouvement de chute libre dans le champ de gravité g. Donc la solution ne peut être que de
la forme: x(t) = αgt2 avec α un nombre sans dimension. On obtient alors: 2α2 +4α2 = α, soit
α = 1/6. Donc: x(t) = gt2 /6 et L = gT 2 /6. La vitesse de la corde lorsqu’elle a complètement
quitté la table est alors: r
2gL
v= .
3
Et son énergie cinétique:
σL2 g
Ec = .
3
On peut maintenant faire un bilan dénergie. Lorsque la corde est complètement en tas sur la
table (sa vitesse est nulle), son énergie mécanique peut s’écrire comme:

E1 = (σL)gh,

où h est la hauteur du centre de gravité du tas de corde par rapport à la table que l’on
prend comme origine. Lorsque la corde lorsqu’elle a complètement passé le trou, son énergie
potentielle est: (σL)g(−L/2) et son énergie mécanique:

σL2 g
E2 = − σgL2 /2.
3
En suivant l’énoncé, on écrit que l’énergie se conserve et on cherche la quantité de chaleur Q
telle que:
σL2 g σL2 g
(σL)gh = − σgL2 /2 + Q = − + Q.
3 6
Soit: Q = σLg(L/6 + h) > 0. Ce qui conclut.
L’énoncé ne demande pas d’explication. Ce n’est pas simple. Si on avait résolu le cas de la
corde initialement tendue sur la table, il n’y aurait pas de dissipation de chaleur. Et dans ce
cas, tous les points de la corde ont la même vitesse à un instant donné. Le fait que Q > 0
dans le montage de l’énoncé est donc clairement une conséquence du fait que la partie de la
corde sur la table est immobile alors que la partie qui pend ne l’est pas. Ce qui se passe, c’est
qu’il y a un mouvement très soudain des éléments de la corde qui passent d’une vitesse nulle
à non nulle au niveau du trou et c’est cet effet qui provoque une dissipation d’énergie sous
forme de chaleur.
Anneau et boules

Deux petites boules percées, en acier, de masses m et rayon a, peuvent coulisser sans frottement
sur un anneau (sans le quitter). L’anneau, de masse M et rayon R  a, est posé verticalement
sur le sol. Initialement les deux boules sont au sommet de l’anneau. Elles sont lachées à t = 0
et commencent à glisser lelong de l’anneau, l’une à gauche, l’autre à droite.

a) Quelle est la valeur minimale m∗ de m à partir de laquelle l’anneau s’élèvera du sol ?

b) Si m < m∗ , au bout de quel temps T les deux boules arrivent en bas de l’anneau? On fera
l’AN pour R = 1 m et a = 1 cm.

1
~ = N~ur de l’anneau.
Les forces sur chaque boule sont son poids −mg~uz et la réaction normale N
Le PFD donne alors (avec r = R constant)

mRθ̈ = mg sin θ , −mRθ̇2 = −mg cos θ + N. (1)

La première équation est intégrée pour donner la conservation de l’énergie, soit


2g
θ̇2 = (1 − cos θ) . (2)
R
En reportant dans la deuxème équation (1) on trouve

N = mg(3 cos θ − 2) ⇒ ~ = mg(3 cos θ − 2)~ur .


N (3)

Ceci est valable pour chacune des deux boules avec leur ~ur respectifs. Chaque boule exerce
donc une force −N ~ sur l’anneau. Les composantes horizontales s’annulent et les composantes
verticales s’ajoutent, soit une force −2mg(3 cos θ − 2) cos θ ~uz . On voit que cette force peut être
dirigée selon +~uz uniquement si cos θ ≥ 0. En ajoutant le poids de l’anneau on trouve la force
totale qu’exerce l’anneau sur le sol:
h i
F~a−>s = −g M + m 2(3 cos θ − 2) cos θ ~uz . (4)

a) L’anneau commencera à se soulèver du sol lorsque cette force s’annule pour la première fois,
soit pour θ = θ∗ avec r
1 1 3M
cos θ∗ = + 1− , (5)
3 3 2m
3M
ce qui est réel à condition que 2m
< 1 ou m ≥ m∗ avec

3
m∗ = M . (6)
2
Si cette condition est satisfaite, l’anneau se soulève dès que θ = θ∗ . Mais ensuite, pour θ > θ∗ ,
l’analyse ci-dessus cesse d’être valable car l’anneau ne constitue plus un référentiel inertiel.

b) Dans le cas m < m∗ l’anneau ne se soulève jamais et l’analyse des forces et accélérations
ci dessus reste valables tout le temps. D’après (2) on a θ̇2 = 4g
R
sin2 2θ , soit θ̇ = 2ω sin 2θ avec
pg R
ω = R ou sindθθ = 2ωdt et donc (si besoin, on donnera la formule sindθθ = 2 ln tan 4θ )
2 2

Z Z Z
1 θ
dθ 1 θ
dθ 1 θ
dθ 1 θ θ0 
(t − t0 ) = θ
= θ θ
= θ θ
= ln tan − ln tan .
2ω θ0 sin 2 4ω θ0 sin 4 cos 4 4ω θ0 tan 4 cos2 4
ω 4 4
(7)
Evidemment, si on pose θ0 = 0 on trouve T = tf − t0 = ∞ : une boule posée exactement en
a
θ0 = 0 y restera indéfiniment. Mais ici les boules ont un rayon a et donc θ0 ' tan θ0 = R
. De

2
2
même, θf = π − Ra . On a ln tan 4R
a a
' ln 4R et ln tan( π4 − a
4R
) ' ln 1 = 0 à des termes O( Ra 2 )
près. On trouve alors s
1 a R 4R
T ' − ln = ln . (8)
ω 4R g a
q
R
Pour l’AN on trouve ' 0.3 s et ln 4R
g a
= ln 400 ' 6, soit T ' 1.8 s. A comparer avec une
q
chute libre sur une hauteur 2R : 2R = g2 T 2 qui donne T = 2 Rg ' 0.6 s.

3
Pluie dans un nuage

On décrit un nuage comme un ensemble de gouttelettes distribuées uniformément dans un


large volume. A l’intérieur de ce nuage, on suppose qu’il pleut, c’est-à-dire que des gouttes de
pluie tombent à travers le nuage. On suppose que cet ensemble de gouttes de pluie représente
une distribution uniforme de gouttes et reste toujours sphérique.

On considère que quand une goutte de pluie rencontre une gouttelette du nuage sur son
chemin, la masse de la goutte de pluie augmente de la valeur de la masse de la gouttelette
rencontrée.

Déterminer la vitesse de l’ensemble sphérique des gouttes de pluie.


Exercice de mécanique sur un système à masse variable, donc qui ressemble au cours sur la
fusée.

Il y a une relation évidente entre la masse des gouttes de pluie M (t) et le rayon de cet
ensemble de gouttes R(t) (toujours de forme sphérique): M (t) = 43 πR(t)3 ρ, où ρ est la masse
volumique des gouttes. Cette masse M (t) va augmenter du fait que cet ensemble tombe dans
le nuage en agrégeant les gouttelettes du nuage lui-même. On note ρ0 la masse volumique du
nuage: comme le nuage est certainement moins dense que l’ensemble des gouttes de pluie,
c’est assez logique de prendre une masse volumique ρ0 6= ρ qui sera donc plus petite que ρ. Il
faut trouver quel est le volume balayé par les gouttes de pluie pendant dt en supposant qu’à
cet instant l’ensemble des gouttes de pluie ont par exemple une vitesse v (suivant z). Chaque
petite couronne de la demi-sphère inférieure de l’ensemble des gouttes balaye le volume vdt
fois l’aire de la petite couronne projetée sur la surface horizontale (perpendiculaire à (z)). Au
final, cela représente le volume: πR(t)2 vdt, où R(t) est le rayon de l’ensemble des gouttes. Et
donc, l’augmentation de la masse des gouttes est:

dM = πR(t)2 vdt ρ0 .

Il reste à écrire l’équation du mouvement pour la boule de gouttes de pluie et le problème


sera complètement posé.
Ṁ v + M v̇ = M g.

Soit (après division par πR(t)2 ):


4 4
ρ0 v 2 + R(t)ρv̇ = R(t)ρg.
3 3
Avec les relations pour M (t) et dM (t), on obtient aussi:
ρ
v = 4Ṙ 0 .
ρ
Ce qui permet d’obtenir une équation pour R(t) seul, puis on pourra ensuite déduire v (via
l’équation ci-dessus):
ρ0
12Ṙ2 + 4RR̈ = gR.
ρ
On n’a pas une équation toute simple qui lie R̈ avec g mais on a quand même un mouvement
0
de chute libre, donc on doit chercher une solution de la forme: R(t) = a ρρ gt2 et déterminer a
pour que cette forme soit solution de l’équation pour R(t). On trouve a = 1/56, donc:
g
v = t.
7
Ce qui conclut.
Il va y avoir une augmentation de température de l’ensemble des gouttes de pluie. Ce n’est
pas évident à montrer. Mais il se trouve que cette augmentation de température restera faible
sur les distances raisonnables et donc le nuage de pluie ne va pas se mettre à s’évaporer auquel
cas le problème posé n’aurait pas de sens.
Courant dans cylindre

a) Un cylindre en cuivre de rayon R = 1 cm et longueur L = 1 m est parcouru par un courant


I = 100 A. En régime stationnaire, déterminer autant que possible les champs électrique et
magnétique. Si le cylindre est thermoisolé, déterminer l’élévation de sa température après une
heure.

b) Faire un bilan local d’énergie.

La masse molaire du cuivre est de 63 g mol−1 , sa masse volumique de ρ = 9 g cm−3 , sa con-


ductivité électrique de σ = 6 × 107 A V−1 m−1 et sa capacité thermique massique de cm =
380 J kg−1 K−1 .

1
a) En première approximation le courant est uniforme, ~j = I
~e .
πR2 z
Pour avoir ce courant, on a
forcément un champ électrique externe

~ 0 = 1 ~j = I ~ez ,
E (1)
σ πR2 σ
IL
ce qui correspond à une tension entre les deux bords du cylindre de V0 = E0 L = πR2 σ
'
−3
5 × 10 V.
Le courant ~j crée un champ magnétique B ~ = µ0 I r 2 ~eθ . Ce champ exerce une force
2πR
volumique sur les porteurs de charge du courant f~ = ~j ∧ B
~ = µ0 I 2 2r 4 ~er . Cette force conduit
2π R
à une séparation des charges positives et négatives et donc une création d’un champ électrique
~ H + ~j ∧ B
(de Hall) selon ~er qui équilibrera cette force: ρE ~ = 0. Avec ρ = ne on a

2
~ H = − µ0 I
E
r
~er (2)
ne 2π 2 R4
et donc une différence de potentiel (de Hall) entre le centre r = 0 et le bord r = R de
RR
VH = 0 E ~ · d~l = − 2µ0 I 2 2 .
4π n e R
On peut leur demander d’estimer VH en prenant un électron libre par atome de cuivre. Avec
9 mol
63 cm3
cela donne environ n = 1029 m−3 et donc ne ' 1010 C m−3 . On a µ0 = 4π10−7 kg m C−2
10−9 kg m2 s−2
et donc VH ' π C
' 3 10−10 J
C
= 3 10−10 V, soit VVH0 ' 10−7 .
= 3 10−10 Ws
C
~ H . La puissance volumique
Pour calculer l’effet Joule, on peut complètement négliger E
dissipée est donc
2 2
~ 0 · ~j = j = I
u̇ = E . (3)
σ π 2 σR4
La capacité thermique volumique est cvol = ρcm et u̇ = ρcm dT
dt
, soit

u̇ I2
∆T = ∆t = 2 4 ∆t ' 5 × 10−4 K s−1 ∆t . (4)
ρcm π σR ρcm
Cela fait une augmentation de la température de 2 K environ en une heure.

b) L’énergie dissipée est amenée par le champ électromagnétique. Le flux est donné par le
vecteur de Poynting
2 2
~ = 1E
P ~ = − I r ~er = − I
~0 ∧ B ~r (5)
µ0 2π 2 R4 σ 2π 2 R4 σ
Le flux d’énergie est radial et pas selon ~ez comme on pourrait croire naivement ! On a (dans le
~ · ~r = 2 et donc
plan) ∇
2
~ ·P
∇ ~ =− I ⇒ u̇ = −∇ ~ ·P ~ , (6)
π 2 R4 σ
ce qui est le bilan souhaité car −∇ ~ ·P ~ correspond bien au flux entrant dans un volume
élémentaire.

2
Fluide entre 2 plaques, gradient de pression oscillant

On considère un fluide Newtonien (ρ, η) emplissant l’espace entre 2 plaques immobiles situées
en y = ±a. Le fluide est illimité dans les directions (x) et (z). On considère les écoulements
dans le plan (xOy) invariants par translation suivant (z) et dont les lignes de courant sont
suivant (x). On suppose qu’un opérateur applique un gradient de pression suivant (x) de la
forme: −G0 cos(ωt).

Déterminer le champ des vitesses dans le fluide dans la limite des basses fréquences (on
p
montrera que cela correspond à: a << δ avec δ = 2η/(ρω)), puis dans la limite des hautes
fréquences. Commenter.
Exercice qui rentre dans la grande catégorie: fluide entre 2 plaques, ce qui en fait un exercice
proche du cours. Il faut faire les calculs correctement. La compréhension physique de δ doit
sortir du calcul dans le cas a >> δ. Si jamais un candidat comprend ce que signifie δ par un
raisonnement de type nombre de Reynolds puis couche limite (avant les calculs), et utilise cela
pour court-circuiter tous les calculs, c’est bien aussi. Mais l’approche suggérée par l’énoncé
est plus simple et plus linéaire: on fait les calculs de type cours et on essaie à comprendre.

On a immédiatement que: ~v = u(y, t)~ex . u(.) n’a pas de dépendance en z (invariance) et


pas de dépendance en x (équation de conservation de la masse). L’ équation du mouvement
suivant (y) pour un petit volume de fluide donne: P = P (x, t), puis l’équation suivant (x):

∂u ∂ 2u ∂P
ρ −η 2 =− .
∂t ∂y ∂x
Et donc:
∂u ∂ 2u ∂P
ρ −η 2 =− = −G(t),
∂t ∂y ∂x
où la fonction G(t) est celle donnée dans l’énoncé: Re(−G0 e−iωt ) (en notations complexes).
On cherche alors le champ des vitesses sous la forme: Re(u0 (y)e−iωt ) avec u0 (.) qui vérifie:

d2 u0 (y)
η + iωρu0 (y) = −G0 .
dy 2
u0 (.) est donc de la forme: −G0 /(iωρ) plus la solution (u00 (.)) de l’équation sans second
membre:
d2 u00 (y)
η + iωρu00 (y) = 0.
dy 2

La suite est simple. On cherche u00 (.) sous la forme ery puis on impose les conditions aux
limites en y = ±a à la solution complète u0 (.). On obtient facilement:

r = ±(1 − i)/δ ≡ ±k,


p
avec: δ = 2η/(ρω). Donc, la solution complète est:

u0 (y) = −G0 /(iωρ) + Aeky + Be−ky .


Les constantes A et B sont déterminées par les conditions aux limites: u0 (y = ±a) = 0.

Dans la limite |ka| << 1 (a << δ), on trouve:

−G0 a2
u(y, t) = (1 − y 2 /a2 ) cos(ωt).
ωρ δ 2
C’est un écoulement de type Poiseuille, en phase avec les variations temporelles du gradient
de pression.

Dans la limite |ka| >> 1, lorsque a est grand devant δ, on trouve que u0 (y) est pratiquement
constant entre les 2 plaques sauf lorsque y est proche de ±a à δ près, auquel cas la vitesse
tend rapidement vers 0 sur l’épaisseur δ. Cette épaisseur δ est donc la couche limite en dehors
de laquelle les effets de viscosités sont négligeables. Le fluide se comporte alors comme un
fluide parfait soumis à une densité volumique de force −G(t). Les effets de la viscosité se font
alors sentir dans la couche limite, et c’est donc dans cette couche que la vitesse du fluide va
tendre rapidement vers 0 au niveau des plaques immobiles. Evidemment, ceci n’a de sens que
lorsque a >> δ.
Lorsque l’on est dans le premier cas étudié, soit a << δ, on est par conséquent dans une
configuration où les effets de viscosité se manifestent dans tout le fluide. C’est pourquoi, on
a retrouvé un écoulement de type Poiseuille.
Particule chargée dans une onde électromagnétique

On veut caractériser le mouvement d’une particule chargée dans une onde électromagné-
tique de basse fréquence.

— Soit un champ électrique uniforme E ~ = E~ex et un champ magnétique uniforme


~ = B~ey . On néglige le poids. Quel est le mouvement d’une particule de charge q
B
et de masse m initialement au repos à t = 0 ?

— Pourquoi ne peut-on utiliser le résultat précédent pour étudier le mouvement d’une


particule chargée dans une onde plane progressive monochromatique de fréquence
Ω si qE/(mc)  Ω

— On suppose que la fréquence de l’onde vérifie qE/(mc) < Ω. Déterminer le mou-


vement de particule chargée. On pourra effectuer un développement en qE/(mcΩ)
et on choisira des conditions initiales pour cette particule qui impliquent un mou-
vement confiné.
Eléments de solution

1) Les équations du mouvement donnent

q q ż
ẍ = E− B (1)
m m
ÿ = 0 (2)
q ẋ
z̈ = B. (3)
m
On pose ω = qB/m la fréquence synchrotron. La résolution se fait facilement en posant
Z = ẋ + iż et on trouve
q
Ż = iωZ + E . (4)
m
On choisit de partir sans vitesse initiale et à l’origine du système de coordonnées et on
trouve (en y c’est évident c’est constant)

E
ẋ = sin(ωt) (5)
B
E
ż = [1 − cos(ωt)] , (6)
B

E
x = [1 − cos(ωt)] (7)

E sin(ωt)
z = [t − ]. (8)
B ω
Là il est intéressant de tracer et d’interpréter. Malgré le champ selon ~ey on a une dé-
rive selon z. La particule accélère selon x, et le champ magnétique la fait tourner jusqu’à
rebiquer et rebrousser chemin en x si bien que le champ électrique la ralentit et elle finit
par s’arrêter et ça recommence. On a une cycloide. Ou plutôt un mouvement de rotation
sur un cercle de rayon R = E/(Bω) avec la coordonnée z du centre du cercle va comme
(E/B)t On a donc comme pour un roulement sans glissement vzcentre = Rω.

2) Si ce champ croisé est en fait une onde plane de très basse fréquence alors E = Bc et la
vitesse de dérive est celle de la lumière. Pire que ça, comme dans une cycloide la vitesse
sur le bord du cercle est double quand le point est en haut, on a même une vitesse à 2c
quand ωt = π. C’est donc impossible et de plus l’expression de la force de Lorentz utilisée
était non relativiste. On pourrait choisir une condition initiale telle qu’on n’aurait pas
cette dérive, mais on aurait toujours le mouvement sur un cercle à la vitesse constante c.

3) si qE/(mc) < Ω donc ω < Ω alors ce mouvement de cycloide n’a jamais le temps
de se mettre en place. On peut donc étudier le mouvement. En effet, si ω  Ω alors
la vitesse typique atteinte est qE/(mΩ) = cω/Ω (l’accélération pendant une période de
l’onde électromagnétique). Donc on peut faire le rapport de l’accélération de Lorentz sur
l’accélération du champ électrique et c’est de l’ordre

(cω/Ω)(B/E) = ω/Ω . (9)

Donc considérer ω/Ω petit revient à considérer que l’effet principal vient du champ
électrique et l’effet secondaire du champ magnétique.
On considère une onde plane

Ex = E cos[Ω(t − z/c)] (10)


By = E/c cos[Ω(t − z/c)] (11)

et on va d’abord la regarder en z = 0. Clairement on va faire une erreur car la particule


va se déplacer selon z un peu mais ça va venir aux ordres supérieurs. On va faire une
résolution itérative et on pose
qE
α = ω/Ω = , (12)
mcΩ
comme le petit paramètre.
(0) (1)
On cherche une solution de la forme vx = vx + vx + . . . et pareil pour vz et les
positions.
A l’ordre le plus bas de (1) on trouve (on choisit judicieusement les conditions initiales
pour éviter toute dérive) juste l’effet du champ électrique

vx(0) = cα sin(Ωt) (13)


vz(0) = 0 (14)

(toujours en mettant des conditions initiales sans problème de dérive)



x(0) = − cos(Ωt) (15)

z (0) = 0 (16)

et donc pas de déplacement selon z. Ensuite quand on regarde l’ordre suivant il vient

v̇x(1) = 0 (17)
qBy (0)
v̇z(1) = v . (18)
m x
On a donc une correction uniquement sur vz . On trouve (avec les conditions initiales qui
vont bien)
qE qE
v̇z(1) = sin(Ωt) cos(Ωt) (19)
mc mΩ
et donc
c
vz(1) = − α2 cos(2Ωt) , (20)
4
c
z (1) = − α2 sin(2Ωt) . (21)
8Ω
Pour continuer il faudrait commencer à pendre en compte que l’onde ne peut plus être
évaluée en z = 0 et c’est un effet d’ordre α3 . En tout cas le mouvement est jusqu’à l’ordre
α2

x = − cos(Ωt) (22)

c
z = − α2 sin(2Ωt) (23)
8Ω
On est plus du tout sur le mouvement d’un cercle mais on a une forme de 8. On quitte
˙ = nq/mE.
le régime linéaire et on n’a plus pour le plasma une relation de type~j ~ Si on
met des nombres on voit qu’en pratique il faut des champs très intenses pour commencer
à avoir α non négligeable ou alors des basses fréquences. Mais on trouve quand même

q/(me c) ' 500Cs/kg/m , (24)

donc à basse fréquence (par exemple Ω = 1 kHz) il faudrait E = 1V /m pour avoir un α


non négligeable.
Puissance EM dissipée dans une barre

On considère un barreau métallique fixe de longueur h = 1 m, de rayon R = 2 cm et


de conductivité γ = 5 107 S/m, soumis à un champ magnétique spatialement uniforme,
variable dans le temps à fréquence f = 50 Hz et de direction parallèle à l’axe du cylindre (z):
~ = B0 cos(2πf t)~ez . On considérera que R << h et on négligera les effets de bord.
B

Déterminer l’énergie EM moyenne dans le temps dissipée dans le barreau. On montrera en


particulier qu’il est nécessaire d’exprimer celle-ci en conservant le(s) terme(s) en f 4 . Estimer
alors l’erreur en % que l’on fait sur le calcul de l’énergie dissipée?
C’est un exercice classique d’EM dans lequel on arrive très rapidement à une réponse qui ne
correspond pas à ce qui est demandé dans l’énoncé. Cette partie est la partie proche du cours.
Ensuite, il faut réfléchir.

~ r, t)
Pour trouver la puissance dissipée dans le barreau, il faut déjà trouver ~j(~r, t) donc E(~
dans le barreau, non nul du fait de l’existence d’un champ magnétique variable dans le temps
(équation de Maxwell-Faraday). Du fait de l’invariance du problème par rotation autour de
l’axe (z) et de l’invariance par translation suivant z (on néglige les effets de bord -énoncé-), on
~ r, t) = E(r,
a: E(~ ~ t). Comme: rot ~ E~ = −∂ B/∂t
~ est parallèle à z on voit que Ez = 0 et tout plan
contenant l’axe (z) contenant B,~ E ~ est perpendiculaire à ce plan, donc: E(~
~ r, t) = E(r, t)~eθ .
On en déduit donc E~ en intégrant l’équation de Maxwell-Faraday sur une surface du cylindre
(de rayon r) perpendicualire à (z). Soit:

~ = − r Ḃ~eθ .
E
2
~ puis l’énergie EM dissipée en moyenne dans le barreau:
On a alors: ~j = γ E
Z R
r2 γπhR4 B02
E= 2πhrdrγ < Ḃ 2 >= (2πf )2 .
0 4 16

Ce qui ne conclut pas car cette énergie est en f 2 et ne contient pas de terme(s) en f 4 comme
demandé dans l’énoncé.

Ce que l’on a oublié, c’est que le champ E~ plus haut va lui même contribuer à modifier le
~ B
champ magnétique initial via l’équation de Maxwel-Ampère: rot ~ = µ0~j = µ0 γ E.
~ On peut
~
vérifier facilement que le terme en c12 ∂ E/∂t est négligeable. On va noter cette contribution
~ 0
au champ magnétique B . On doit donc trouver B ~ 0 tel que:
h r i
~ ~ 0
rotB = µ0 γ − Ḃ~eθ .
2
~ 0 est aussi suivant (z), mais atten-
Avec: Ḃ = −B0 (2πf ) sin(2πf t). On voit facilement que B
~ 0 est nul en dehors du barreau et donc nul aussi en R. On trouve (il faut utiliser un
tion B
petit rectangle situé dans un plan contenant (z) ayant une longueur suivant (z) à la cote r et
l’autre à la cote R):
~ 0 (~r, t) = µ0 γ (r2 − R2 )Ḃ~ez .
B
4
~ 0 va donc produire un champ E
Et ce champ B ~ 0 tel que:

~ E
rot ~ 0 = −∂ B
~ 0 /∂t.

Une fois calculé ce champ E ~ 0 , on pourra calculer à nouveau l’énergie EM dissipée en utilisant
cette fois-ci le champ électrique: E ~ +E ~ 0 et voir ce que cela change. On revient plus loin sur
la nécessité d’itérations supplémentaires. Pour trouver E ~ 0 , on procède comme plu haut mais
en faisant attention à l’intégration sur la surface:
Z
0 µ0 γ 1 r 02 µ0 γ
E (r, t) = − B̈ (r − R2 )r0 dr0 = − B̈[r3 /4 − R2 r/2]
4 r 0 4
On peut alors recalculer l’énergie dissipée dans le barreau:
Z R
E= 2πhrdrγ < (E(r, t) + E 0 (r, t))2 >
0

On obtient:  
γπhR4 B02 2 µ0 γ(2πf )R2 2
E =[ (2πf ) ] 1 + α( )
16 4
Avec α = 11/24. On peut se tromper dans le calcul de ce coefficient, ce n’est pas très
grave. Ceci étant, la forme de l’expression est intéressante: on trouve la valeur plus haut (en
f 2 ) et un terme correctif, qui donne un terme en f 4 . Le terme correctif dans les crochets:
2
( µ0 γ(2πf
4
)R 2
) ' 0.36 avec les valeurs numériques de l’énoncé: il n’est donc pas négligeable et
c’est pour cela qu’il fallait pousser le calcul (au moins) jusque là. On peut faire apparaitre
l’épaisseur de peau pour ce montage mais ce n’est pas demandé. Si on voulait aller plus loin
dans les itérations, en suivant la méthode ci-dessus, on obtiendrait une nouvelle correction
2
en ( µ0 γ(2πf
4
)R 4
) ' 0.12 qui n’est pas encore négligeable par rapport à l’unité. Cela veut dire
que l’on a obtenu une expression pour l’énergie dissipée dans le barreau qui contient bien un
terme en f 4 , comme demandé, mais l’erreur sur ce calcul est encore d’environ 10 %.
Propagation d’une onde électromagnétique dans un gaz

Etudier la propagation d’une onde électromagnétique de fréquence ω dans de l’hélium à tempé-


rature ambiante et pression atmosphérique. On considérera que tous les électrons sont élastique-
ment liés avec une même force de rappel que l’on estimera. Déterminer la fréquence caracté-
ristique de cette liaison élastique. Déterminer la relation de dispersion de l’onde ainsi que
l’indice de réfraction de l’hélium dans ce modèle. Comparer les résultats avec ceux pour la
propagation d’une onde électromagnétique dans un plasma.

1
Ecrivons la force de rappel Fel = mω02 r où m est la masse de l’électron et r sa distance au
m 2 2
noyau. Cela donne une énergie potentielle ω r .
2 0 0
L’énergie cinétique est du même ordre de
grandeur (on pourra éventuellement le démontrer pour des mouvements classiques circulaires
dans ce potentiel élastique; pour cela il suffit d’égaler la force centrifuge et la force de rappel).
Donc mω02 r02 ' Eliaison ' 25eV. (On a Z = 2 ce qui double l’énergie de liaison du premier
électron par rapport à l’hydrogène.) Pour r0 on prendra le rayon de Bohr r0 ' 5 × 10−11 m, et
la masse de l’électron est donnée par mc2 ' 5 × 105 eV. Donc
 ω r 2  ω r 2 1 c
0 0 0 0
mc2 ' Eliaison ⇒ ' × 10−4 ⇒ ω0 ' 0, 7 × 10−2 ' 4 × 1016 s−1 (1)
c c 2 r0
Il faudra comparer ce ω0 à ω. Pour une onde em, une fréquence de ce ω0 correspond à une
λ c r0
longueur d’onde de 2π
= ω0
' 0,7×10−2
, soit λ ' 103 r0 ' 5 × 10−8 m = 50 nm. On est bien dans
l’ultraviolet.
L’équation pour la position classique ~xe d’un électron est

m ~x¨e = −mω02~xe + q E(~


~ xe ) + q ~x˙ e ∧ B(~
~ xe ) . (2)

Pour une onde em on a |B| ' |E|


c
et comme |~x˙ e | ' ω0 r0 ' 10−2 c on a |~x˙ e ∧ B|
~ ' 10−2 |E|
~
et on pourra négliger le terms magnétique dans la force. Une autre difficulté vient de ce que
~ x) = E
l’onde plane E(~ ~ 0 e−i(ωt−~k·~x) est à prendre à la position de l’électron ~x = ~xe , ce qui
fait de (2) une équation différentielle non-linéaire compliqué. Pour éviter cela, on se placera
~
dans un régime de λ  r0 de telle sorte que eik·~x peut être supposé constant sur l’échelle
d’un atome. On remplacera alors ~x par ~xn , la position du noyau. Cette condition implique
c
ω ' ck  r0
' 1, 4 × 102 ω0 . On ne regardera donc que des ondes em avec ω  100 ω0 , disons
18 −1
ω ≤ 10 s . L’équation (2) est alors simplement une équation pour un oscillateur forcé:
~
~ 0 eik·~xn
−iωt qE
~xe (t) = ~x(0)
e e , ~x(0)
e = (3)
m(ω02 − ω 2 )

Ce déplacement des électrons (il y en a 2 par atome et n atomes par unité de volume) génère
un courant volumique associé

~ ik·~xn ~
ωp2 2nq 2
~j = 2nq~x˙ e = −2iωnq q E0 e e−iωt ⇒ ~j(t, ~x) = −0 iω ~ ~x) ,
E(t, ωp2 = .
m(ω02 − ω 2 ) (ω02 − ω 2 ) 0 m
(4)
où on a défini la fréquence plasma par ωp . Ce courant est à insérer dans l’équation de Maxwell

~     ~
∂E ωp2 1 ωp2 ∂E
~ ~
rotB = µ0 j + µ0 0 = −iωµ0 0 1 + 2 ~
E = 2 1+ 2 . (5)
∂t ω0 − ω 2 c ω0 − ω 2 ∂t

2
On obtient l’équation d’onde comme d’habitude en utilisant rot(rotE) ~ = grad(divE)−∆
~ ~ avec
E
~ = ρ/0 = 0 car à l’échelle des variations spatiales de E
divE ~ on a un gaz globalement neutre.

~   ~
ωp2 ∂ 2E
∆E ~ = rot ∂ B = ∂ rotB
~ = −rot(rotE) ~ = 1 1+ . (6)
∂t ∂t c2 ω02 − ω 2 ∂t2

On obtient la relation de dispersion et l’indice de réfraction n


  s
2
ω p ωp2
c2~k 2 = ω 2 1 + 2 , n = 1 + . (7)
ω0 − ω 2 ω02 − ω 2

L’indice n (et donc aussi ~k) devient imaginaire si ω02 < ω 2 < ω02 + ωp2 . Estimons ωp2 . n est
de l’ordre de 6 × 1023 par 22 l, soit environ 3 × 1025 m−3 . On a la constante de structure fine
q2
1
137
' 0 ~c
et mc2 ' 5 × 105 eV, ~ ' 10−34 Js et 1 eV ' 1, 6 × 10−19 J. A la fin on trouve
14 −1
ωp ' 10 s , soit ωp  ω0 et l’interval des ω qui donne des n et k imaginaires est donc très
petit.
Pour discuter le cas d’unq plasma avec des électrons libres il suffit de poser ω0 = 0 et on a
ω2
alors c2~k 2 = ω 2 − ωp2 et n = 1 − ωp2 . Dans un plasma il y a propagation (k réel) uniquement
si ω > ωp .

3
Passoire tournante

On considère une passoire qui tourne autour de son axe vertical de symétrie. On modé-
lise cette passoire par un cylindre de rayon R, et dans le disque de base du cylindre sont
disposés un très grand nombre N de trous. Chaque trou est de forme cylindrique, de
rayon a et de hauteur h  a, tous deux très petits devant R. On veut savoir en combien
de temps la passoire (qui est initialement remplie d’eau) va se vider.

En régime stationnaire, et en négligeant les pertes par les trous de la passoire, quel est le
profil du fluide et sa distribution de pression ?

Quel est le débit dans chaque trou en fonction de la pression au fond de la passoire ?

Etablir les équations décrivant la passoire qui se vide. Discuter leur domaine de validité.
Eléments de solution

Il est crucial de voir que l’on a un écoulement de Poiseuille dans les trous de la passoire
(donc il ne s’agit pas d’une passoire de cuisine...).

1) L’équilibre hydrostatique donne

~ 1
gradP = −ρg~ez + Ω2~er . (1)
2
Et donc avec P0 la pression de l’atmosphère le profil de pression est
1
P = P0 − ρg(z − H) + ρΩ2 r2 . (2)
2
où H est la hauteur du point le plus bas de la surface libre. La surface libre est une
équipotentielle et donc
1
z = H + Ω2 r 2 . (3)
2g
Au passage, on calcule que le volume d’eau est
Z Z
R 2 Ω2
V = 2πrdrdz = 2πz(r)rdr = πR2 (zmin + ) (4)
4g

~
2) L’écoulement de Poiseuille est classique et correspond à un régime de Stokes gradP =
µ∆~v
∂P 1
= µ ∂r (r∂r vr ) . (5)
∂z r
Après intégration et conditions aux limites

r2 − a2 ∂P
vr = (6)
4µ ∂z
et le débit par section est alors
πa4 ∂P
D=− . (7)
8µ ∂z
Le signe moins parce que le fluide coule vers la direction de perte de charge (pression
décroissante). Dans le cas présent la pression à l’entrée des trous est
1
Ptrous = P0 + ρgH + ρΩ2 r2 , (8)
2
donc le gradient est  
∂P 1 1 2 2
=− ρgzmin + ρΩ r . (9)
∂z h 2
3) Le nombre de trous par unité de surface est N/(πR2 ). On fait une description continue
des trous et on intègre pour avoir la débit total. On obtient
Z
N
Dtot = 2πrdrD . (10)
πR2
On trouve en prenant que D dépend de r puisque selon la position du trou, le gradient
de pression ne sera pas le même,
 
2 H
Dtot = πR +U (11)
τ

avec le temps caractéristique τ et la vitesse U

1 N a4 ρg a4 N Ω2 ρ
= U= . (12)
τ 8µhR2 32µh
Or le débit correspond à une baisse du volume donc
∂V
V̇ = Ḣ = πR2 Ḣ = −Dtot . (13)
∂H
Donc on a
Ḣ = −H/τ − U (14)
dont la solution est

H = Hi exp(−t/τ ) + U τ [exp(−t/τ ) − 1] . (15)

On a H = 0 quand  
Hi + U τ
t = τ ln . (16)

En effet grâce au terme de forçage de rotation on n’a pas une simple décroissance expo-
nentielle car cela pousse toujours même quand H → 0. On peut aussi discuter la validité
de la modélisation. En effet il faut que le temps d’établissement du régime permanent
pour le profil soit bien plus court que le temps τ pour que la description hydrostatique
reste valable. Il faut que les trous de passoire soient assez longs pour que Poiseuille (va-
lable pour un cylindre infini) soit une bonne description.

4) Question subsidiaire : que se passe-t-il quand H = 0 ? En effet il reste du fluide sur


les bords à cause du profil parabolique...
Mise en rotation d’un cylindre rempli d’eau

Un tube cylindrique vertical de rayon a = 2 cm et de longeur L = 1 m, fermé à l’une de ses extrémités,


est rempli presque entièrement d’eau. Le cylindre tourne autour de son axe avec une vitesse angulaire
variable ω(t).

Discuter les équations qui déterminent la vitesse ~v (t, ~r) d’un élément d’eau et la pression dans l’eau,
ainsi que le profil de la surface de l’eau en fonction du temps. Discuter les différents régimes en
fonction de la valeur de ω̇/ω. (On pourra considérer par exemple ω = ω0 eαt .)

1
On ne peut négliger la viscosité ici, sinon l’eau resterait immobile. On rappelle que pour l’eau
la viscosité dynamique est η ' 10−3 kg m−1 s−1 , la masse volumique ρ = 103 kg m−3 , la dimension
spatiale caractéristique ici est a = 2 × 10−2 m. Donc
η
' 2, 5 × 10−3 s−1 . (1)
ρa2
On verra que l’inverse fournit une échelle de temps caractéristique de 400 s qui est le temps qu’il faut
pour que la viscosité soit suffissisante pour entrainer toutes les couches d’eau. Pour ω̇/ω < 10−3 s−1
on aura sensiblement la même distribution de vitesses que pour une situation stationnaire avec un
ω constant. Par contre pour ω̇/ω > 10−2 s−1 , seules les couches près du bord suivent la rotation
imprégnée par le bord, alors que les couches plus près du centre ont un mouvement de rotation
correspondant aux ω(t0 ) aux temps antérieurs.
Par les symétries on s’attend à ce que ~v (t, ~x) soit dirigé selon ~eθ et indépendent de θ et z, au moins
si on néglige les effets de bord (a  L). Certes, la surface libre de l’eau n’est pas plane et avec un
ω variable cette suface se modifie au cours du temps et on a donc forcément une petite composante
de ~v selon ~ez . Nous allons faire l’hypothèse qu’elle est négligeable. De même, la pression ne dépend
que de r et z :
~v (t, ~x) = v(t, r)~eθ , p = p(t, r, z) . (2)
Dans la suite on n’écrira pas la dépendence en t explicitement. On traitera l’eau comme incompress-
~ · ~v = 0, ce qui est satisfait par notre ansatz (2). Pour écrire l’équation
ible, ρ est constant, et donc ∇
de Navier Stokes  
∂~v ~ ~ + η∆~v − ρg~ez
ρ + ~v · ∇~v = −∇p (3)
∂t
on a besoin de ∆~v en coordonnées! cylindriques (à fournir au candidat qui le demande). Alternative-
−y
ment, on peut écrire ~v = v(r)
r
et calculer le Laplacien en cartésien. On trouve
x
   v 0 (r) v(r)

00
∆ v(r)~eθ = v (r) + − 2 ~eθ . (4)
r r
~ = v(r)∇θ =
Ensuite, ~v · ∇ v(r) ∂
et ∂
~e = −~er , donc
r ∂θ ∂θ θ
2
~ v = − v(r) ~er
~v · ∇~ (5)
r
Posons d’abord p(r, z) = pe(r, z) − ρgz afin de tenir compte de l’effet du poids. On obtient alors la
même équation que (3) mais avec pe et sans le dernier terme. On substitue (4) et (5) et projette sur
~er et ~eθ pour obtenir
 
∂~v 00 v 0 (r) v(r) v(r)2 ∂ pe(r)
ρ = η v (r) + − 2 , ρ = . (6)
∂t r r r ∂r

2
La deuxième éq donne pe une fois v(r) connu:
Z r
v(r0 )2
p(r) = pe(0) + ρ dr0 , (7)
0 r0
Rr 0 2
ce qui montre, bien sur, que la pression augmente avec r. On a alors p(r, z) = pe(0)−ρgz+ρ 0
dr0 v(rr0 ) .
La surface libre est déterminée par p(r, z) = patm , son profil est donné par
Z
1 r 0 v(r0 )2
h(r) = dr . (8)
g 0 r0
Finalement il reste à résoudre la première équation (6) avec les conditions aux bords

v(t, r = 0) = 0 , v(t, r = a) = a ω(t) . (9)


ηω
Ces conditions de bords permettent d’éstimer le membre de droite de la première éq. (6) à η aω
a2
= a
et
ω̇ η −3 −1
le membre de gauche à ρaω̇. Le membre de gauche sera alors négligeable si ω
 ρa2
' 2, 5 × 10 s .
00 v 0 (r) v(r)
Dans ce cas, une bonne approximation est la solution quasi-statique de v (r) + r − r2
= 0 qui
correspond à l’eau tournant avec le cylindre, soit v(r) = rω(t). Dans ce cas (8) redonne le résultat
ω2 r2
bien connu ∆h(r) = 2g
.
ω̇
Si ω
est plus grand il faut résoudre l’équation aux dérivées partielles. On essaie la séparation des
variables en posant v(t, r) = g(t)f (r) ce qui conduit à
η ηC
ġ = C g ⇒ g(t) = g(0)e ρ t , (10)
ρ
f0 f
f 00 + − 2 = C f (11)
r r
g(0)f (a) ηC
La condition de bord aω(t) = v(t, a) = g(t)f (a) donne ω(t) = a
e ρ t. C’est donc convenable
pour ω(t) = ω0 eαt avec
ρα
C= , g(0)f (a) = aω0 (12)
η
Ensuite, l’équation (11) pour f admet forcément une solution (réelle) avec f (0) = 0 et f (a) > 0. Les
condition de bords sont alors satisfaites et on a une solution. En fait, la solution pour f est connue

puisque (11) est l’équation de Bessel d’ordre 1: f (r) = Im J1 (i Cr). Finalement, la solution s’écrit
f (r)
v(t, r) = aω(t) . (13)
f (a)
√ q √
ρα2 α
f (a) fait intervenir l’argument Ca = η
. Si ce nombre est très petit, on a f (r) = Im J1 (i Cr) '
r et on retrouve la solution quasi-stationnaire. Par contre si ce nombre est par exemple ègal à 2, on
f (r) 3
a f (a)
' 0.62 ar + 0.31 ar 3 pour 0 ≤ r ≤ a et on voit bien que la vitesse angulaire près de r = 0 est
maintenant nettement inférieure à ω(t).

3
Tubes (lignes) de vortex

On considère un fluide idéal incompressible de champ des vitesses ~v (x, y, z, t). On note
ω ~ v ) sa vorticité. Ce fluide a une distribution spatiale particulière: sa vor-
~ (x, y, z, t) = rot(~
ticité est localisée seulement dans deux tubes C1 et C2 qui se croisent comme sur la figure
ci-dessous et la vorticité est nulle partout ailleurs. De plus, chaque tube C1 et C2 est un tube
(dit vortex) tel que ω
~ dans le tube est parallèle à l’axe du tube en chaque point. Egalement,
on considère que ces tubes sont assez fins, assez proches de lignes.
R
1) Pour chaque tube, on définit: KS = ω ~ où l’intégrale est prise sur S qui est une section
~ ·dS
S
du tube donc parpendiculaire à ω ~ est parallèle à ω
~ (et dS ~ ). Montrer que KS est constante le
long d’un tube. On note K1 et K2 les valeurs pour les tubes C1 et C2 respectivement.
R
2) Déterminer l’intégrale ~ · ~v prise sur l’ensemble du fluide en fonction de K1 et K2 .
dV ω
A priori, K1 et K2 peuvent dépendre du temps. On admet que ce n’est pas le cas. Quelle est
alors la généralisation du résultat de la question 2). Au vu de l’exercice, y-a-t’il des hypothèses
à cette généralisation?

C2 passe à l’intétieur de la boucle formée par C1


1) Résultat immédiat avec la propriété: div(~ω ) = 0.
R
2) On décompose l’intégrale: dV ω~ · ~v sur chaque tube 1 et 2:
Z Z Z
~ · ~v = dV1 ω
dV ω ~ · ~v + dV2 ω
~ · ~v .

Mais dV1 = dS1~n1 d~l1 et:


Z Z Z I
dV1 ω ~ · ~v = dS1~n1 d~l1 ω ~ ][d~l1~v ] = K1 d~l1~v .
~ · ~v = [dS1~n1 ω

Car l’intégrale sur le volume complet veut dire que l’on fait tout le tour du tube. Mais:
I Z
d~l1~v = ~
rot~v dS
S1

où l’intégrale est prise sur la surface enlacée par la ligne C1 (on rappelle que l’on a pris les
tubes assez fins, ce qui est utile ici). On poursuit:
Z Z
~
rot~v dS = ω ~ = K2 .
~ dS
S1 S1

Car le seul domaine de S1 pour lequel la vorticité est non nulle est celui du tube C2 . Finale-
ment: Z
~ · ~v = K1 K2 .
dV1 ω

C’est exactement la même chose pour le tube 2 et donc:


Z
dV ω~ · ~v = 2K1 K2 .

Ce qui répond à la question.


R
Comme on admet que K1 et K2 ne dépendent pas du temps, cela veut dire que: ~ · ~v est
dV ω
aussi indépendant du temps. Ce sera vrai dès lors que la vorticité forme des boucles dans le
volume du fluide. C’est assez restrictif, mais c’est ce que donne l’exercice.
Pendule double

On considère un pendule double constitué d’une deuxieme masse m2 accrochée par un fil de
longueur l2 à une première masse m1 elle même attachée par un fil de longueur l1 au plafond.
Etablir les équations du mouvement et les analyser. Pour simplifier on supposera que tous les
mouvements sont dans le même plan x − O − z.

1
D’abord, remarquons que rien n’oblige ni le premier pendule ni le deuxième pendule d’avoir un
mouvement dans un plan, ni d’avoir les deux mouvements plans dans le même plan. Mais c’est
déjà assez compliqué, même en supposant les deux mouvements dans le même plan...

On prendra l’axe des z vers le bas et l’axe x horizontal. On appelera θ1 l’angle entre le
premier fil et la verticale tel que θ1 augmente avec x, et θ2 l’angle entre le deuxième fil et la
verticale, compté dans le même sense, voir la figure.

On appelera T1 et T2 les modules des tensions des fils 1 et 2. Pour la première masse on écrit
le PFD en coordonnée polaires et on projette sur ~eθ (ce qui élimine T1 ):

m1 l1 θ̈1 = −m1 g sin θ1 − T2 sin(θ1 − θ2 ) . (1)

Pour la deuxième masse on écrit le PFD en coordonnées cartésiennes:

m2 ẍ2 = −T2 sin θ2 , m2 z̈2 = m2 g − T2 cos θ2 . (2)

Il faut ensuite exprimer x2 et z2 en termes des θi et éliminer T2 . On aura alors deux équations
différentielles pour les deux fonctions θ1 et θ2 . On a

x2 = l1 sin θ1 + l2 sin θ2 , z2 = l1 cos θ1 + l2 cos θ2 , (3)

ce qui entraine

ẍ2 = l1 cos θ1 θ̈1 − l1 sin θ1 θ̇12 + l2 cos θ2 θ̈2 − l2 sin θ2 θ̇22


z̈2 = −l1 sin θ1 θ̈1 − l1 cos θ1 θ̇12 − l2 sin θ2 θ̈2 − l2 cos θ2 θ̇22 . (4)

Afin d’éliminer T2 dans (2) on prend la combinaison

cos θ2 ẍ2 − sin θ2 z̈1 = −g sin θ2 (5)

2
En utilsant (4), ceci devient

l1 cos(θ2 − θ1 ) θ̈1 + sin(θ2 − θ1 ) θ̇12 + l2 θ̈2 = −g sin θ2 . (6)

Il reste à éliminer T2 de l’équation (1). En développant sin(θ1 − θ2 ) = sin θ1 cos θ2 − cos θ1 sin θ2
et en utilisant les deux équations (2) on obtient

m1 l1 θ̈1 + (m1 + m2 )g sin θ1 = m2 (sin θ1 z̈2 − cos θ1 ẍ2 ) , (7)

et en utilisant de nouveau (4):


 m1 + m2
l2 cos(θ1 − θ2 ) θ̈2 + sin(θ1 − θ2 ) θ̇22 + rl1 θ̈1 = −rg sin θ1 où r = . (8)
m2
Les deux équations (4) et (6) sont donc

l1 cos(θ2 − θ1 ) θ̈1 + sin(θ2 − θ1 ) θ̇12 + l2 θ̈2 = −g sin θ2 ,

l2 cos(θ1 − θ2 ) θ̈2 + sin(θ1 − θ2 ) θ̇22 + rl1 θ̈1 = −rg sin θ1 (9)

On remarquera la symétrie entre les deux équations sous l’échange 1 ↔ 2, aux facteurs r près.
Pour tenter de résoudre ces équations, on peut regarder différents cas limite: petits angles
et/ou une masse beaucoup plus grande que l’autre. Commençons par les petits angles. On a
alors cos → 1 et sin α → α et on néglige tous les termes d’ordre θ2 . Les éqs (6) et (8) deviennent

l1 θ̈1 + l2 θ̈2 = −gθ2 , l2 θ̈2 + rl1 θ̈1 = −rgθ1 . (10)

Un cas limite simple est m2  m1 . Dans ce cas, la masse 2 “tire” beaucoup plus fort que la
masse 1 et les deux fils sont quasi-aligés. On devrait donc trouver θ1 ' θ2 . En effet, dans ce cas
r ' 1 et la différence des deux eqs (9) donne θ1 ' θ2 ≡ θ, et donc (l1 + l2 )θ̈ = −gθ : on retrouve
un pendule de longueur l1 + l2 . Mais on devrait retrouver le pendule de longueur l = l1 + l2
lorsque r = 1 même si les angles ne sont pas petits. Clairement, pour r = 1, θ1 = θ2 = θ est
une solution de (6) et (8) si (l1 + l2 )θ̈ = −g sin θ, mais il n’est pas évident que ce soit la seule
possibilité dans ce cas.

3
Systéme mécanique tournant

On considère le montage suivant: un losange ABCD peut tourner autour d’un axe vertical
(z) tel que le point A est fixé sur l’axe et le point C peut se déplacer le long de l’axe. Deux
masses m sont situées en B et D et un ressort de raideur k et de longueur au repos nulle relie
les points A et C. Voir figure.

Les tiges et le ressort sont de masse négligeable. De plus, on suppose que le système tourne
très vite autour de (z).

Déterminer les positions (trajectoire) d’équilibre du montage et les petits mouvements autour
de ces positions (s’ils sont possibles).

Indication: on pourra utiliser comme système [le losange (tiges et masses)].


On prend comme système le losange (tiges et masses) (énoncé). Les forces extérieures qui
s’exercent sur lui sont: les poids de chaque masse m et la force de rappel du ressort. On
peut remarquer que le moment suivant (z) de chacune de ces forces est nul. Donc, le mo-
ment cinétique suivant (z) pour le système (losange) est conservé, de même que son énergie
mécanique. On obtient facilement:

Lz = 2ma2 sin2 θ φ̇ = cte = L.

Avec φ̇ la vitesse angulaire de rotation autour de (z). Puis:


1
E = ma2 [θ̇2 + sin2 θ φ̇2 ] + 2mga(1 − cos θ) + k[2a cos θ]2 = cte.
2
Le premier terme est le terme d’énergie cinétique et les 2 derniers représentent l’énergie poten-
tielle. Il faut faire attention ici car ce qui compte pour l’équilibre, c’est l’énergie potentielle
effective (comme dans le cas des systèmes gravitationnels), c’est pourquoi on va procéder
comme suit.

On doit déjà trouver une équation pour la variable θ seule. On remplace φ̇ dans la conservation
de l’énergie par sa valeur tirée de la conservation de Lz et on dérive par rapport au temps.
On obtient une équation différentielle (équation du mouvement) pour θ de la forme:

θ̈ + F (θ) = 0.

On cherche alors les positions (trajectoire) d’équilibre. On obtient ces positions avec la con-
dition: θ̈ = 0 à l’équilibre, soit F (θ = θequil ) = 0. Par ailleurs, on suppose (énoncé) que le
système tourne très vite, donc la position d’équilibre doit être proche de θ = π/2. On cherche
donc une valeur de θ telle que θ = π/2 − γ avec γ petit, telle que: F (θ) = 0.

Après calculs simples, on trouve:


2mga
γequil = .
L2 /(2ma2 )
+ 4ka2
Comme le système tourne vite, L2 /(ma2 ) >> ka2 et on a bien γequil << 1 et:
2mga
θequil = π/2 − .
L2 /(2ma2 ) + 4ka2
Cela veut dire que la trajectoire: θ = θequil donne l’équilibre du système tournant (vite).

On cherche alors les petites variations δθ (telles que θ = θequil −γequil −δθ) autour de l’équilibre,
on obtient:
L2 /(2ma2 ) + 4ka2
δ θ̈ + δθ = 0.
2ma2
L2 /(2ma2 )+4ka2
Donc, des petits mouvements sont possibles avec: ω 2 = 2ma2
> 0. Cet équilibre est
donc stable.
Gravité anneau

On considère un anneau de masse M , que l’on modélise par un cercle de rayon R de


centre O dans le plan Oxy, et de masse linéique λ = M/(2πR). On cherche à caractériser
le mouvement d’une particule ponctuelle dans le champ de gravité généré par cet anneau.

— Quelles sont les équations du mouvement d’une particule placée initialement proche
du centre de l’anneau ?

— Le moment cinétique de cette particule est-il conservé ?

— Quel est le mouvement d’une particule placée initialement sur l’axe de symmétrie
Oz à l’infini ?

— On rajoute une masse MO en O. Jusqu’à quelle distance peut-on place une parti-
cule pour qu’elle reste proche de O ?

— On considère une configuration de trois anneaux de rayon R de même masse, de


même centre et respectivement dans les plans Oxy, Oyz et Oxz. Quel est le mou-
vement d’une particule proche du centre ?
Eléments de solution

1) On utilise que le potentiel gravitationnel est en 1/r ce qui se trouve une fois le problème
reformulé comme un problème d’électrostatique, et on somme les potentiels créés par
chaque élément. Z 2π
λRdφ
V (~r) = −G (1)
0 ~
|~r − ξ(φ)|
où ~r est la position du point où on calcule le champ et ξ~ est la position d’un élément de
masse du cercle. Ici il faut choisir entre des coordonnées cylindriques ou des coordonnées
sphériques. On peut faire l’un ou l’autre. Les cylindriques sont plus simples mais il faut
bien distinguer l’angle φ pour la position dans le cercle de l’angle θ associé au point où
on regarde le potentiel. On a simplement
~
ξ(φ) = R~er (φ) = R(cos φ~ex + sin φ~ey ) (2)
~r = r~er + z~ez = r(cos θ~ex + sin θ~ey ) + z~ez . (3)

Etant données les symétries, on regarde un point dans le plan Oxz et donc θ = 0. Il vient
p p
~ = (r − R cos φ)2 + (R sin φ)2 + z 2 = r2 + R2 + z 2 − 2rR cos φ .
|~r − ξ| (4)

On fait un développement en r  R et z  R. En utilisant


r
p r2 + z 2 r
r2 + R2 + z 2 − 2rR cos φ − R 1 + − 2 cos φ (5)
R2 R
on développe
 2
R r r2 + z 2 3 2r cos φ
p ' 1 + cos φ − + . (6)
r2 + R2 + z 2 − 2rR cos φ R 2R2 8 R

On peut donc calculer (1) en utilisant que la moyenne de cos est nulle et celle de cos2
donne 1/2. Et on trouve
 
GM r2 z2
V '− 1+ − . (7)
R 4R2 2R2

Si on le veut en Cartésien on utilise juste r2 = x2 + y 2 . En voyant que ce potentiel


diminue quand r augmente à z fixe, on réalise bien l’instabilité qui pousse à s’écarter
vers le cercle. On trouve l’accélération gravitationnelle
GM  r 
~a = ~
e r − z~
e z . (8)
R3 2
On peut discuter qualitativement que c’est comme un ressort sur l’axe z mais en revanche
répulsif dans la direction radiale. Si on veut les équations du mouvement il faut juste
utiliser la formule générale de l’accélération en cylindrique. Mais on ne va pas résoudre
(sauf pour la coordonnée z ou c’est facile).

~ en coordonnées cylindriques est r~er + z~ez donc la force n’est pas cen-
2) Le vecteur OM
trale.

3) Sur l’axe, r = 0, et on peut calculer le potentiel et la force exactement. Le potentiel


est
GM
V = −√ (9)
z 2 + R2
et la conservation de l’énergie mécanique donne

ż 2 GM
=√ . (10)
2 z 2 + R2
On peut pas résoudre simplement, mais pour z  R on retrouve l’équation harmonique
comme vu précédemment.

4) Avec cette masse au centre, on rajoute une accélération en

GM
~a2 = − (r~er + z~ez ) . (11)
(r2 + z 2 )3/2

On peut regarder en utilisant le développement proche du centre de l’accélération due


à l’anneau jusqu’à quel rayon la force radiale va vers le centre plutôt que repousse vers
l’extérieur. C’est quand
 
2 2M2 2/3 2
r < R − z2 (12)
M
5) Le potentiel proche du centre est (7) et donc en Cartésien

∝ x2 + y 2 − 2z 2 . (13)

Si on rajoute les deux autres cercles dans les plans Oxz et Oyz, on va sommer

(x2 + y 2 − 2z 2 ) + (y 2 + z 2 − 2x2 ) + (z 2 + x2 − 2y 2 ) = 0 . (14)

Le potentiel est-il constant ? Non, mais il faudrait développer à des ordres plus élevés
le potentiel car les termes les plus bas seront tous en x4 , y 4 , z 4 , x2 y 2 , x2 z 2 , . . . . Et si on
continue à rajouter des cercles on peut éliminer également ces termes et ne conserver
que des x6 , y 6 , . . . . En poussant jusqu’au bout, et en rajoutant plein de cercles, on aura
formé une sphère creuse, et là le potentiel sera constant avec aucune force à l’intérieur
grâce au théorème de Gauss.
Induction cylindre tournant

On considère un cylindre infini conducteur (conductivité σ), de rayon R, en rotation


~ = Ω~ez autour de son axe de symétrie et plongé dans un champ magnétique B
Ω ~ = B0~ex .

Quelle est la distribution de courants dans le cylindre lorsque celle-ci est devenue station-
naire ?

Quelles sont les forces et moments subis par le cylindre ?

Aide : La loi d’Ohm pour un conducteur en mouvement dans un champ électrique et un


~ + ~v ∧ B).
champ magnétique est ~j = σ(E ~
Eléments de solution

On utilise les coordonnées cylindriques naturellement associées avec le cylindre.

1) En régime stationnaire, toutes les dérivées temporelles des champs sont nulles et donc
par exemple
~ = 0 , div~j = 0 .
rotE (1)
~ D’abord il faut ~v .
On calcule ~v ∧ B.
~ ∧ r~er = rΩ~eφ .
~v = Ω (2)

En utilisant ~eφ = cos θ~ey − sin θ~ex = x~ey − y~ex , il vient

~ = −r cos θΩB0~ez = −xΩB0~ez .


~v ∧ B (3)

On regarde maintenant la loi d’Ohm


~ + ~v ∧ B)
~j = σ(E ~ (4)

En stationnaire la conservation de la charge implique

div~j = 0 . (5)

Or on peut vérifier facilement (quand on passe aux coordonnées Cartésiennes) que


~ = 0,
div~v ∧ B (6)

donc on pourrait se passer d’un champ électrique. Comme sans champ électrique cela fait
un courant selon ~ez , les lignes de courant ne sortent pas du conducteur donc c’est pas la
peine d’en rajouter et on conclut qu’on n’a pas de champ électrique. On a simplement

~j = −σΩB0 x~ez . (7)

2) Deux directions s’ouvrent. Calculer les forces volumique et le moment des forces par
unité de longueur du cylindre. Ou alors voir si les courants induisent un autre champ
magnétique et éventuellement un autre champ électrique etc... et raffiner la solution.
Pour les forces c’est juste
f~ = ~j ∧ B
~ = −σΩB 2 x~ey ,
0 (8)
et la moyenne sur une section est nulle. Pour le moment par unité de longueur c’est
Z
T~ = − σΩB02 x2 rdddθ~ez (9)

πR2
= −σΩB02 ~ez (10)
4
On a bien un freinage.

Maintenant on essaie de raffiner en calculant le champ à l’intérieur créé par les cou-
rants. On cherche donc
~ B
rot ~ 1 = µ0~j . (11)
On ne peut avoir de dépendance en z par symétrie. On peut prendre le rot de cette
équation (en utilisant qu’un champ magnétique est sans divergence) pour avoir une idée.
On trouve
−∆B ~ 1 = µ0 rot~j = µ0 σΩB0~ey . (12)
On a donc l’intuition que B ~ 1 est selon ~ey , et la fonction devant doit être de Laplacien
constant (et ne pas dépendre de z d’après les symétries du problème), donc on essaie x2
et y 2 avec possiblement un soupçon de xy. On trouve que pour satisfaire (11) c’est x2
qu’il faut et
2
~ 1 = − x σµ0 ΩB0~ey
B (13)
2
On a
2
~ 1 = x y σµ0 Ω2 B0~ez .
~v ∧ B (14)
2
Et donc on trouve ~j1 = σ~v ∧ B~ 1.
On peut aussi continuer et chercher le correctif au correctif
~ 2 = µ0~j1 .
rotB (15)

On voit qu’on a rot~j1 =∝ x2~ex + ∝ xy~ey et donc ∆B ~ 2 est du même genre. On a alors By
qui est une combinaison de x y et xy , mais pour B2 seul By ∝ x3 y convient. Le même
3 3

genre de raisonnement donne Bx ∝ x4 . En le faisant correctement on trouve


 
~ σ 2 µ20 Ω2 B0 1 4 3
B2 = − x ~ex + x y~ey (16)
6 4

On peut vérifier que les forces sont toujours de moyenne nulle dans les sections.
Figure 1 – Lévitateur Alcon fonctionnant avec deux cylindres tournants (rotations de
sens opposés) et un aimant en lévitation au dessus.
Onde EM entre 2 plans

On considère une onde électromagnétique (EM) plane harmonique qui se propage dans le vide
suivant l’axe (x) avec un vecteur d’onde ~k = (kx , ky ) et polarisée linéeairement suivant l’axe
(z). Cette onde entre dans l’espace vide confiné entre 2 plans conducteurs d’équations: y = 0
et y = a. On admet que dans les conditions du montage, le champ électrique est nul au
niveau des plans conducteurs et que, pour le champ magnétique, By = 0 au niveau des plans
conducteurs.

Dans ces conditions, trouver quelle est la forme complète du champ électrique associé à l’onde
entre les plans conducteurs. Montrer sur cet exemple que la vitesse de propagation de l’énergie
EM est égale à la vitesse de groupe.
C’est un exercice proche du cours. La variation (ou la difficulté) est qu’il faut être précis dans
les calculs.

Pour le début de la question, il faut savoir que le champ électrique d’une onde EM dans
le vide, de vecteur d’onde ~k = (kx , ky ) et polarisée linéairement suivant l’axe (z), s’écrit en
notations complexes:

~
E(x, y, z, t) = E0 exp[i(kx x + ky y − ωt)]~ez .

Avec k = ω/c. Entre les 2 plans y = 0 et y = a, il y a une onde de cette forme de vecteur
d’onde (kx , ky ) et une onde réfléchie de vecteur d’onde (kx , −ky ). Le champ électrique entre
les 2 plans devient donc:

~
E(x, y, z, t) = (A exp[i(kx x + ky y − ωt)] + B exp[i(kx x − ky y − ωt)])~ez .

Avec A et B des constantes (complexes) à déterminer à partir des conditions aux limites, ici
le champ électrique est nul au niveau des plans conducteurs. On trouve alors:

~ πy
E(x, y, z, t) = 2iA sin( ) exp[i(kx x − ωt)]~ez .
a
Le champ magnétique se trouve alors avec la relation de Maxwell-Faraday.
∂Ez π πy
= 2iA cos( ) exp[i(kx x − ωt)] = iωBx .
∂y a a
∂Ez πy
− = 2Akx sin( ) exp[i(kx x − ωt)] = iωBy .
∂x a
Soit:
π πy kx πy
(Bx , By ) = (2A cos( ) exp[i(kx x − ωt)], 2A sin( ) exp[i(kx x − ωt)]).
ωa a iω a
Ce qui conclut pour le début de la question. On observe bien que By = 0 au niveau des
plans conducteurs, comme indiqué dans l’énoncé. On ne demande bien évidemment pas de
commenter le fait que Bx ne soit pas nul (hors programme).
Pour répondre à la fin de la question, on doit déjà trouver la vitesse de groupe. Avec la forme
du champ EM entre les plans, on obtient la relation de dispersion:
π2
+ kx2 = ω 2 /c2 .
a2
Soit, une vitesse de groupe:
r
(πc/a)2
vg = dω/dkx = c 1 −
ω2
Ensuite, pour étudier la vitesse de propagation de l’énergie EM, on doit calculer la valeur
moyenne dans le temps (sur une période T = 2π/ω) du vecteur de Poynting, soit (attention
à la notation complexe):
1
< P~ >T =< Re[E ~ ∧B~ ∗ /µ0 ] >T .
2
Et la moyenne temporelle de la densité volumique d’énergie EM (attention à la notation
complexe):
1 ~ E~ ∗ /2 + B.
~ B~ ∗ /(2µ0 )] >T .
< u >T =< Re[0 E.
2
On obtient:
|A|2 πy
< P~ >T = 2 kx sin2 ( )~ex .
µ0 ω a
 2 2

2 kx c 2 πy π 2 c2 /a2 2 πy
< u >T = 0 |A| (1 + 2 ) sin ( ) + cos ( ) .
ω a ω2 a
En suivant ce qui est appris dans le cours, on sait alors que l’énergie EM qui traverse une
surface S perpendiculaire au plans y = 0 et y = a (et de longueur L arbitraire suivant (z)) et
pendant T vaut: Z
T < P~ >T dS~ex .
S
Si l’énergie EM se propage à la vitesse v, cette quantité peut aussi s’écrire comme:
Z
vT < u >T dS.
S
Pour trouver v, on doit ensuite écrire l’égalité entre ces 2 expressions. Soit:
Z Z
vT < u >T dS = T < P~ >T dS~ex .
S S
Les calculs sont simples, on trouve que v = vg , ce qui conclut.
Un calcul du même type est fait en cours dans le cas d’une onde plane progressive harmonique
dans le vide, auquel cas: v = vg = c.

Vous aimerez peut-être aussi