2019 PC Sujets Pour Oral de Physique
2019 PC Sujets Pour Oral de Physique
2019 PC Sujets Pour Oral de Physique
Filière PC
Quelle est l’évolution de la position de la barre suite à un petit déplacement vers le bas ?
Et celle du courant ?
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)
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
ω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.
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).
ρ0 ∂t vr = −∂r δP (2)
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.
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
2(kr0 )2
α= . (11)
1 + 2(kr0 )2
√
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)
3 x3/2 9ν 2
1− x− = t . (7)
4 4 16Ri ri2
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
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.
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
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).
Rπ
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 .
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~
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
mω
ξ= x (6)
~
effectué pour calculer l’intégrale. La variable ξ est sans unité. Le potentiel peut alors s’écrire
~ω 2
V (x) = 2
ξ .
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.
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.
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
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:
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
)
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
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)
2
Ondes EM en champs intenses
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.
γ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:
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
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.
∂G0 ∂G0
+ [A + Bδh] = 0.
∂t ∂x
∂G0
On doit donc calculer par exemple: ∂t
. Mais G0 (.) est une fonction de δh, donc:
∂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
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 à:
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).
nsolv,1 Dv nsolv,1 ∆P 2
Jx0 = = b > 0.
V1 Lb V1 12ηe
12ηD nsolute
∆P = .
b2 nsolv,1
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.
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.
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?
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
ce qui est juste une séparation de variable où on devine ce qu’on va trouver. On a alors
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
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
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)
δP = δP (z)ei(ωt−kx) , (16)
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)
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)
∂t s = iωs = vz . (25)
En injectant dans (20) est en ne gardant que les termes d’ordre un on a alors
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).
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
2
Transfert thermique barreau
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 ?
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
~ et on trouve
1) On fait un bilan avec le flux ~j = −λ∇T
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.
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.
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.
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
∂θ ∂θ v
θ̇ = v = − v = cos θ , (7)
∂` ∂x h
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.
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.
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
mω
−φ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)
|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.
(*) 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 .
mẍ = Nx − kvrel,x .
mÿ = Ny − kvrel,y .
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.
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.
U0 h0 L = 2U hL.
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)
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
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.
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
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 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 ?
— 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
~ E
Ensuite rot ~ donne toujours en négligeant les dérivées de l’amplitude
~ = −∂t B
~ ∧E
∇ϕ ~ = ωB
~ (3)
~ = −∆E
2) On fait comme d’habitude rotrotE ~ si ρe = 0 et sans courants et on trouve
1 ~ = 0.
∆ − 2 ∂tt E (5)
c
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
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).
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:
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 ρ.
— 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 ?
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)
div~j = 0 . (6)
on en déduit que
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
~ 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
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)
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)
ω 2 ≡ ωnm
2
= g γnm tanh(γnm H) . (13)
3
Double oscillateur
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.)
2
Puis on prend les carrés et on somme :
∂2 ∂2 1 ∂2 ∂2
+ = + 2 . (9)
∂x21 ∂x22 2 ∂X 2 ∂x2
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.
\ 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
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
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
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).
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:
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
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:
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 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 ?
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
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 )
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
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
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
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 peut montrer que lorsque l’on cherche des solutions à symétrie sphérique, on peut poser
u(r) 1
ψ(~r) = r
√
4π
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) = √
4π
ne dépendant que de r on a ∆R(r) = r dr2
(rR(r)) et
R(r) u(r)
Schrödinger pour ψ = √
4π
= √
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
4π
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
~2 2α e2 ~2
α(α − 1) = 0 , = , − =E . (6)
2m r0 4π0 2m r02
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.
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.
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
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.
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 :
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
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 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
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
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
ω 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.
√
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.
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
−φ
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
-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 .)
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).
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:
~
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
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)
∂ 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
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
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é.
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)
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
~ 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
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)~
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
— 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
T~
~aM/O = + ~g − ~aO . (3)
m
~ = r~er , on a simplement
En prenant des coordonnées polaires pour OM
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
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
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
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).
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.
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)
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 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.
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 .
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.
∂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:
−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.
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)
Bω
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
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
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
cα
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
~ 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
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
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
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
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)
τ
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
On a H = 0 quand
Hi + U τ
t = τ ln . (16)
Uτ
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.
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
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?
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 ω
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 ):
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
ce qui entraine
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
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
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).
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.
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
— Quelles sont les équations du mouvement d’une particule placée initialement proche
du centre de l’anneau ?
— 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 ?
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 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
~ en coordonnées cylindriques est r~er + z~ez donc la force n’est pas cen-
2) Le vecteur OM
trale.
ż 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.
GM
~a2 = − (r~er + z~ez ) . (11)
(r2 + z 2 )3/2
∝ x2 + y 2 − 2z 2 . (13)
Si on rajoute les deux autres cercles dans les plans Oxz et Oyz, on va sommer
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
Quelle est la distribution de courants dans le cylindre lorsque celle-ci est devenue station-
naire ?
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)
div~j = 0 . (5)
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
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
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.