Corrigé de Laurent BEAU
Professeur de Sciences Physiques en Math Spé MP*Lycée Mohamed V. CASABLANCA
N’hésitez pas à me signaler des erreurs ou à me suggérer des commentaires ou des réponses plus "élégantes". Merci.Collisions nucléaires et fragmentation
Première partieAnalyse cinématique d'une collision
- Cinématique du problème à deux corps.
Nous noterons M1 et M2 les masses respectives et B1 et B2 les positions respectives des noyaux cible (indice 1) et projectile (indice 2)
- {RG=M1r1+M2r2M1+M2=A1r1+A2r2A1+A2r=B1B2=r2−r1 ⇒{r1=RG−A2A1+A2rr2=RG+A1A1+A2r
- D’après sa définition, la quantité de mouvement relative des deux noyaux s’écrit : p=A1A2A1+A2m(v2−v1)soit, en introduisant la masse réduite µ du système : μ=A1A2A1+A2m, on obtient :p=μ(v2−v1). On remarque que p ne dépend pas du référentiel dans lequel elle est calculée. En particulier : p∗=p.
-
La quantité de mouvement totale du système s’écrit : P=p1+p2
Le moment cinétique total du système par rapport à O est défini par : L=r1∧p1+r2∧p2
Dans le référentiel barycentrique : →ℓ=L∗=r1∧p∗1+r2∧p∗2
Dans le référentiel barycentrique : P∗=→0 donc p=p∗=p∗2=−p∗1. On en déduit : →ℓ=r∧p
En utilisant la composition des vitesses : p1=p∗1+M1vG et p2=p∗2+M1vG. Il vient alors :
L=→ℓ+(M1r1+M2r2)∧vG
ou encore L=→ℓ+RG∧P
résultat qui constitue le théorème de Kœnig relatif au moment cinétique.
- Energie cinétique
L’énergie cinétique totale Ec s’écrit dans le référentiel du laboratoire : Ec=p122M1+p222M2 qui s’écrit en composant les vitesses : Ec=E∗c+Ec(G) où Ec(G)=P22m(A1+A2) car P∗=→0.
Ceci constitue le théorème de Kœnig relatif à l’énergie cinétique. - Energie potentielle d'interaction coulombienne.
- Considérons une sphère uniformément chargée en volume. Le champ électrostatique E(M) créé par cette distribution de charge en M est radial car tout plan passant par M et par le centre de la sphère est un plan de symétrie.
De plus, la norme de ce champ ne dépend que de r (en coordonnées sphériques) car la distribution est invariante par rotation autour du centre de la sphère. On en déduit : E(M)=E(r)er.
Le flux de ce champ à travers une sphère de rayon r vaut donc : ∬E(M)dS=E(r)4πr2
La charge intérieure à cette sphère vaut : {Qint=Qsir>RQint=Q(rR)3sir<R
Le théorème de Gauss s’écrivant : ∬E(M)dS=Qintε0, on en tire :
{E=Q4πε0r2ersir>RE=Q4πε0rR3ersir<R
On en déduit alors le potentiel électrostatique : dV=−E.dl=−E(r)dr
En prenant V = 0 à l’infini , il vient :
{V=Q4πε0rsir>RV=−Q4πε0r22R3+V(0)sir<R
La distribution de charge étant uniforme, le potentiel au centre de la sphère se calcule simplement :V(O)=∭sph ˊe recharg ˊe eQ43πR3dτ4πε0r avec dτ=4πr2dr
⇒V(O)=38Qπε0R
Il vient alors :
{V=Q4πε0rsir>RV=3Q8πε0R(1−13(rR)2)sir<R -
L’énergie d’interaction électrostatique Eel s’écrit alors Eel=Q1V2(r>R)=Q2V1(r>R), c’est-à-dire :
Eel=Z1Z2e24πε0r
- Considérons une sphère uniformément chargée en volume. Le champ électrostatique E(M) créé par cette distribution de charge en M est radial car tout plan passant par M et par le centre de la sphère est un plan de symétrie.
Collision et évolution du système composite
- Energie du système avant le contact.
- L’énergie totale du système s’écrit :
E=E∗c+Ec(G)+Eel+Ub - Le système des deux noyaux étant isolé, sa quantité de mouvement est constante au cours du mouvement. Le centre de masse a donc un mouvement rectiligne et uniforme.
-
L’énergie cinétique du centre de masse s’écrit :
Ec(G)=P22m(A1+A2)=P24mA=cte car A1 = A2 = A
Cette énergie cinétique étant constante, calculons-la à l’instant initial :
Ec(G)=(p1(0))24mA=Elab2
La moitié de l’énergie cinétique initiale du noyau projectile correspond à l’énergie cinétique du centre de masse du système. -
Le système est isolé donc sa quantité de mouvement et son moment cinétique sont conservées.
Dans le référentiel barycentrique, le moment cinétique du système est : →ℓ=r∧p. En définissant le mobile équivalent (ou mobile fictif) M par un mobile de masse µ et de position GM=r, le moment cinétique du mobile fictif est également →ℓ=r∧p.
Si ce moment cinétique est nul à l’instant initial (paramètre d’impact nul) alors le mouvement du mobile équivalent dans le référentiel barycentrique est rectiligne.
Si ce moment cinétique est non nul alors le mouvement est plan, dans le plan perpendiculaire à →ℓ. - Calculons l’énergie cinétique barycentrique : E∗c=Ec(mobilefictif)
E∗c=12μ(drdt)2 avec drdt=˙rer+r˙θeθ ⇒E∗c=12μ(˙r2+r2˙θ2)
Le moment cinétique s’écrit : →ℓ=r∧p=rer∧μ(˙rer+r˙θeθ)=μr2˙θez ⇒˙θ2=(ℓμr2)2
On en déduit : E∗c=12μ(˙r2+(ℓμr)2) et donc l’énergie totale barycentrique :
E∗=12μ˙r2+ℓ22μr2+Eel+Ub
- L’énergie totale du système s’écrit :
- Barrière coulombienne.
- Si le nucléon test est un neutron, l’énergie potentielle totale de ce nucléon test est : Vtot(r)=Vnuc(r) car le neutron n’est pas chargé.
Si le nucléon test est un proton (particule de charge e), l’énergie potentielle totale de ce nucléon s’écrit : Vtot(r)=Vnuc(r)+Eel(r) où :
Si r > 0 : {Eel(r)=Ze24πε0rVnuc(r)=0 Si r < 0 : {Eel(r)=eV(r)=3Ze28πε0R(1−13(rR)2)Vnuc(r)=−V0
- Le noyau projectile étant constitué de Z protons et de N neutrons, il existe (par rapport au zéro d’énergie potentielle) entre les deux noyaux identiques(dont les centres sont distants de 2R au contact) une barrière d’énergie potentielle coulombienne Ucoul=Z.Ze24πε0(2R) : Ucoul=(Ze)24πε0(2R)≈73MeV
-
L’énergie totale barycentrique du système est :
E∗=cte=E∗c(t=0)+Ub=(12μ˙r2+ℓ22μr2)r=R+Ucoul+Ub
Or, à t = 0, l’énergie cinétique barycentrique du système vaut Elab/2 (cf. question 1c) de la 2ème partie)
D’où : E∗=Elab2+Ub
Au moment du contact : {˙r2>0r=2R⇒Elab2−Ucoul−ℓ22μ(2R)2>0
Le moment cinétique du mobile fictif vaut, à t = 0 :
ℓ=μbv(t=0)
et l’énergie cinétique barycentrique :
E∗c(t=0)=Elab2=12μv2(t=0)
On en déduit : ℓ22μ(2R)2=14(bR)2Elab2
La condition pour que les noyaux fusionnent s’écrit donc :
Elab2(1−(b2R)2)>Ucoul - Dans le référentiel barycentrique, les trajectoires des 2 noyaux sont hyperboliques lorsque la condition précédente n’est pas remplie (répulsion des deux noyaux).
-
L’énergie minimale permettant la fusion correspond à un paramètre d’impact nul (collision frontale) c’est-à-dire Eminlab=2Ucoul=146MeV
Prenons Elab>Eminlab ; la fusion est possible si ℓ22μ(2R)2<Elab−Eminlab2
La valeur maximale du moment cinétique est donc : ℓmax=2R√Am2(Elab−Eminlab) car μ=Am2
On prend ElabA=10MeV soit Elab=400MeV>Eminlab. La fusion est donc possible à condition que bR<2√1−EminlabElab≈1,6.
- Si le nucléon test est un neutron, l’énergie potentielle totale de ce nucléon test est : Vtot(r)=Vnuc(r) car le neutron n’est pas chargé.
- Compression.
- L’énergie de rotation du système composite est : Erot=l22μr2avec r = 2R soit : Erot=l28μR2
-
L’énergie totale barycentrique s’exprime sous la forme :
E∗=Erot+Ec,int+2AU(ρ)
L’énergie de rotation est considérée comme constante. La densité maximale est obtenue quand l’énergie cinétique interne est nulle (compression maximale). Ceci est également valable pour la densité minimale.
On a alors : U(ρmax)=E∗−Erot2A=U0+K18(ρmax−ρ0ρ0)2
On en tire :
ρmax=ρ0(1+√18K(E∗−Erot2A−U0))
Ensuite, le système subit une expansion jusqu’à une densité minimale donnée par : ρmax−ρ0=ρ0−ρmin c’est-à-dire :
ρmin=ρ0(1−√18K(E∗−Erot2A−U0))
Le système oscille ainsi autour de la position d’équilibre caractérisée par ρ = ρ0. -
Pour une collision frontale (paramètre d’impact nul), le composite ne tourne pas après la collision et son énergie cinétique de rotation est nulle.
L’énergie barycentrique du système est : E∗=Elab2+2AU0
On en déduit :
ρmax=ρ0(1+√18KElab4A) et ρmin=ρ0(1−√18KElab4A)
A.N : ρmax = 1,42.ρ0 = 0,22 fm–3 et ρmin = 0,58.ρ0 = 8,9.10–2 fm–3
Fragmentation du système composite formé
- Etude de la pression.
- dUT=TdS−pdV d’où à température nulle : p=−(∂(2AU(ρ))∂V)T=0
- Le volume du système composé de 2A nucléons s’exprime ainsi : V=2A/ρ. Donc dV=−2Aρ2dρ
On en tire : p=ρ2∂U∂ρ - p=ρ2K9(ρ−ρ0ρ20) avec l’expression (2) pour l’énergie interne.
- L’état d’équilibre défini par la densité ρ0 correspond à la pression nulle.
Etudions la stabilité de cet équilibre : ∂p∂ρ=K9ρρ0(3ρρ0−2)⇒∂p∂ρ|ρ=ρ0=K9>0. Une légère compression du composite (dρ > 0) entraîne une augmentation de pression donc une expansion du composite (retour vers l’état d’équilibre). L’équilibre est donc stable vis-à-vis des oscillations de densité.
La pression devient négative pour ρ > 0 ; ceci correspond aux forces de cohésions nucléaires responsables de la cohésion du noyau.
- Une équation d'état réaliste à température nulle.
- U=Ucin+Upot=Cρ2/3+t0ρ+t3ρ2.
Posons ρ=ρ0(1+ε) avec ε << 1. En développant U au 2ème ordre en ε, il vient :
U=(Cρ2/30+t0ρ0+t3ρ20)+ε(23Cρ2/30+t0ρ0+2t3ρ20)+ε2(t3ρ20−C9ρ2/30)
En identifiant à U=U0+K18ε2, on obtient :
{Cρ2/30+t0ρ0+t3ρ20=U023Cρ2/30+t0ρ0+2t3ρ20=0 d’où l’on tire :
t3=C3ρ−4/30−U0ρ−20=989MeV.fm6
t0=−43Cρ−1/30+2U0ρ−10=−396MeV.fm3 - p(ρ)=ρ2∂U(ρ)∂ρ=23Cρ5/3+t0ρ2+2t3ρ3
- ∂p∂ρ=0⇔59Cρ2/3+t0ρ+3t3ρ2=0 soit numériquement : 41,7ρ2/3−396ρ+2967ρ2=0 (avec le MeV comme unité d’énergie et le femtomètre comme unité de longueur).
On trouve alors deux valeurs de ρ non nulles qui vérifient cette équation : {ρ′s=1,2.10−3fm−3ρs=1,0.10−1fm−3
- U=Ucin+Upot=Cρ2/3+t0ρ+t3ρ2.
- L'équation d'état à température non nulle.
- L’équation d’état du gaz parfait s’écrit : pV=nRT=NNART soit pGP=ρkBT
La pression totale est donc : p=ρkBT+ρ2∂Upot∂ρ=ρkBT+t0ρ2+2t3ρ3 - (∂p∂ρ)T=0⇔kBT+2t0ρ+6t3ρ2=0. Cette équation du second degré admet 2 racines réelles si son discriminant réduit est strictement positif c’est-à-dire : Δ′=t20−6t3kBT>0 donc pour :
T<Tc avec Tc=t206t3kB
On a alors : ρ′s=−t06t3(1−√1−TTc) et ρ′s=−t06t3(1+√1−TTc) - kBTc=t206t3=26,4MeV et ρs(T=Tc)=−t06t3=6,7.10−2fm−3
- On peut écrire : (∂p∂ρ)T=6t3(ρ−ρs)(ρ−ρ′s) pour T < Tc.
Donc (∂p∂ρ)T<0 pour ρ′s<ρ<ρs et T < Tc.
Pour T > Tc, (∂p∂ρ)T>0 : le système composite se comporte « normalement » c’est-à-dire qu’une légère compression du composite (dρ > 0) entraîne une augmentation de pression donc une expansion du composite et vice versa, ceci quelle que soit la densité particulaire ; le système ne se fragmentera donc pas.
- L’équation d’état du gaz parfait s’écrit : pV=nRT=NNART soit pGP=ρkBT
- Equation d'état et fragmentation.
- L’équation (5) est une équation d’onde (équation de d’Alembert).
c2s=1χλ=−Vλ(∂p∂V)T=0 car (∂p∂V)T=0=1(∂V∂p)T=0
On en déduit : c2s=ρλ(∂p∂ρ)T=0 car V=2A/ρ - D’après la question 1d) de la troisième partie : c2s=K9λρ2ρ0(3ρρ0−2)
λ représente la masse volumique à l’équilibre soit : λ = mρ0 d’où : c2s=K9(ρρ0)2(3ρρ0−2) et au voisinage de ρ0, c2s≈K9m c’est-à-dire :
cs≈c√K9mc2≈c6=5.107m.s−1 - Si on admet que quel que soit ρ et quelle que soit T : c2s=ρλ(∂p∂ρ)T, on obtient en utilisant le résultat de la question 3a) de la troisième partie que c2s est du même signe que (∂p∂ρ)T. Donc c2s<0 quand ρ devient inférieure à ρs(T) au cours de l’expansion du système.
La solution de l’équation Δp−1c2s∂2p∂t2=0 diverge quand c2s<0, la pression devenant alors très importante : il y a ainsi possibilité de brisure du système composite (fragmentation).
- L’équation (5) est une équation d’onde (équation de d’Alembert).
Aucun commentaire:
Enregistrer un commentaire