1 – Exercices : 12 - Conduction - Convection - Rayonnement Sciences Physiques MP 2017-2018

Exercices : 12 - Conduction - Convection - Rayonnement

A. Régime stationnaire
1. Température d’interface et régime stationnaire
On met en contact, suivant leur surface commune, d’aire S, deux conducteurs thermiques limités par des plans
parallèles. En régime stationnaire, l’ensemble des deux conducteurs, de même épaisseur e, se comporte comme
un système dont l’état ne dépend que de la seule coordonnée spatiale z le long de l’axe perpendiculaire à leur
plan. En outre, les températures des faces des deux conducteurs qui ne sont pas en contact sont maintenues
aux valeurs T1 = 293 K et T2 = 373 K respectivement. On désigne par λ1 et λ2 les conductivités thermiques des
deux corps.
1. Quelle est l’expression de la résistance thermique de chaque conducteur en fonction de e, S et de sa
conductivité thermique ? En déduire la résistance thermique Rth de l’ensemble des deux conducteurs
placés en série.
2. En s’appuyant sur l’analogie avec la loi d’Ohm, montrer que la température Ti à l’interface est telle que :
Ti − T1 = α(T2 − T1 ) où α est une quantité que l’on exprimera en fonction des résistances thermiques
Rth1 et Rth2 des deux conducteurs. En déduire Ti en fonction de T1 , T2 , λ1 et λ2 .
3. Application : Calculer Ti pour un conducteur organique comme le corps humain, (λ1 = 0, 5 W ·m−1 ·K−1 )
en contact avec du bois (λ2 = 0, 2 W · m−1 · K−1 ) puis en contact avec du cuivre (λ2 = 390 W · m−1 · K−1 ).
1 e 1 e R2 λ1 T1 +λ2 T2
Réponses : R1 = λ1 S , R2 = λ2 S , Rtot = ( λ11 + 1 e
λ2 ) S , T2 − Ti = R1 +R2 (T2 − T1 ), Ti = λ1 +λ2 , Ti = 43 ◦ C,
Ti ≃ 100 ◦ C.

2. Résistance thermique cylindrique, sphérique

On considère un manchon cylindrique de conductivité λ, de hauteur H, de rayon intérieur r1 et de rayon extérieur
r2 . La paroi intérieure est portée à la température T1 et la paroi extérieure à T2 > T1 . Le régime permanent
indépendant du temps est établi.
1. Représenter les lignes de densité de courant de transfert thermique. Réfléchir aux invariances et aux
symétries.
2. En déduire la forme de dépendance en fonction de r du vecteur densité de courant de transfert thermique
et de la température T (r).
3. Calculer la résistance thermique équivalente de ce manchon.
4. Reprendre la même étude pour une coquille sphérique.

Réponses : T (r) = (T2 − T1 ) lnlnrr/r 1

2 /r1
+ T1 , ~jcond = − λr lnT2r−T 1
2 /r1
~er , R = ln r2 /r1
λ2πH , T (r) = r1 r2 1
r2 −r1 [ r1 − 1r ](T2 − T1 ) + T1 ,
~jcond = − λ2 1 2 (T2 − T1 )~er , R = 2 1 .
r r r −r
r r2 −r1 4πλr1 r2

3. Ailettes de refroidissement
Pour éviter l’échauffement d’un appareil dû à l’effet Joule, on munit son boı̂tier d’ailettes de refroidissement
métalliques. Chaque ailette est parallélépipédique, de dimensions a = 2, 0 mm (épaisseur), b = 10 cm (largeur)
et c = 20 cm (longueur). On pourra admettre que a est négligeable devant b. En fonctionnement, le boı̂tier
de l’appareil M sera maintenu à la température TM = 60 ◦ C. L’air extérieur, qui circule, est de température
constante et uniforme TA = 20 ◦ C, sauf au voisinage immédiat de l’ailette, entourée d’une couche limite d’air
thermiquement peu conductrice dont la température reste localement voisine de celle de la surface de l’ailette.
Dans l’ailette, on admettra que le transfert thermique, de type conductif, est monodimensionnel dans la direction
de l’axe Ox. Il obéit à la loi de Fourier, la conductivité thermique étant λ = 16 W · m−1 · K−1 . On note T (x)
la température de l’ailette à l’abscisse x. Il existe aussi un transfert thermique de l’ailette vers l’air ambiant, à
travers la couche limite. Le flux thermique au niveau d’une surface dS de l’élément de l’ailette de longueur dx
est de la forme :
dP = h(T (x) − TA )dS
où h = 150SI est un coefficient uniforme et constant.
1. Expliquer la loi de Fourier et donner l’unité du coefficient h dans le système international.
2. Écrire le bilan des transferts d’énergie pour
r la tranche d’ailette comprise entre les abscisses x et x + dx,
λa
en régime permanent. On posera : L = et on donnera la valeur numérique de L ainsi que son unité.
2h
En déduire l’équation différentielle dont T (x) est la solution.
3. Résoudre cette équation différentielle pour déterminer l’expression de T (x). On vérifiera que L ≪ c et
on pourra considérer c comme infini pour simplifier.

JR Seigne Clemenceau Nantes

Sciences Physiques MP 2017-2018 Exercices : 12 - Conduction - Convection - Rayonnement – 2

4. Donner l’expression de la puissance thermique dP sortant de la surface latérale dS de la tranche d’ailette

comprise entre les abscisses x et x + dx. En déduire l’expression de la puissance thermique totale P
évacuée par l’ailette, faire l’application numérique.
5. Exprimer et calculer la puissance thermique transmise du boı̂tier de l’appareil M à l’ailette en x = 0.
Conclure.
6. Combien faudrait-il fixer d’ailettes sur le boı̂tier pour évacuer un flux thermique total de 0, 9 kW ? La taille
de chaque ailette peut-elle être réduite sans changer notablement l’ensemble des résultats précédents ? Si
oui, expliquer comment et pourquoi.
2
Réponses : L = √ 1, 03 cm ; ddxT2 − LT2 = − TLA2 , T = TA + (TM − TA ) exp − Lx ; P = λ ab
L TM ≃ 100 W ; la puissance
évacuée est en b a, on peut réduire c tout en conservant c ≫ L.
4. Isolation et pertes thermiques
Un tube cylindrique de rayon a, à la température T1 = 304 K, est séparé de l’extérieur, à la température à
la température T2 = 275 K, par une gaine cylindrique d’épaisseur e, constituée d’un matériau de conductivité
thermique λ = 0, 9 W · m−1 · K−1 .
1. Représenter les lignes de densité de courant de transfert thermique. Réfléchir aux symétries et aux
invariances.
2. Quelle est l’expression de la résistance thermique de la gaine sur une longueur l ?
3. Trouver la puissance thermique qui traverse la gaine sur une longueur l. Application numérique : a =
20 cm, e = 4 cm et l = 1 m.
On augmente l’isolation du tube au moyen d’une couche supplémentaire cylindrique, d’épaisseur e′ ,
constituée d’un matériau isolant de conductivité thermique λ′ = 0, 03 W · m−1 · K−1 .
4. Quelle doit être la valeur de e′ pour que les pertes thermiques soient divisées par 10 ?

5. Homéothermie
Une sphère de rayon a est maintenue en permanence à la température T1 , dans un milieu fluide qui, à grande
distance de la sphère, est à la température T0 < T1 . La conductivité thermique du fluide est notée λ.
On néglige toute discontinuité de température (transfert pariétal parfait) à la surface de la sphère.
Le problème est étudié en régime permanent.
1. Expliciter la puissance P thermique produite par la sphère.
2. On donne T1 = 310 K, T0 = 280 K et a = 25 cm (modélisation d’un animal homéotherme avec un
rapport surface/volume comparable à celui d’un être humain). Calculer P si le fluide est de l’air (λ =
2, 6 × 10−2 W · m−1 · K−1 ).
3. L’homéothermie est-elle plus aisée pour un petit animal ou pour un gros ?
2
Réponses : div ~j = r12 ∂(r∂r j) = 0, j = rα2 , T = βr + γ, T (r) = (T1 − T0 ) ar + T0 , P = j(r = a)4πa2 = λ(T1 − T0 )4πa,
P = 2, 45 W ; plus a grand, plus P est grande, plus facile pour un petit animal.
6. Transferts entropiques
−−→ T
Dans un milieu conducteur thermique de conductivité λ, on pose ~jS = −λgrad ln où T0 est une constante
T0
positive quelconque.
On note aussi s l’entropie massique et ρ la masse volumique du milieu étudié.
Enfin, l’irréversibilité des transferts thermiques impose la création d’une entropie par unité de temps et de
volume de matériau ṡc .
1. Relier ~jS , ρ, s et ṡc .
2. On considère une transformation isobare et on note cp la capacité thermique massique correspondante. ρ
et cp sont supposés constants. On étudie une situation unidimensionnelle en régime permanent. Exprimer
ṡc , examiner son signe et conclure.
λa2
Réponses : div ~jS + ρ ∂s
∂t = ṡc ;
T
T0 = ax + b en régime stationnaire, ṡc = (ax+b)2 > 0, irréversible.
7. Effet Joule et effet Thomson
Sir William Thomson (Lord Kelvin) découvre en  l’effet maintenant appelé Joule-Kelvin ou effet
Thomson. Il résulte du passage d’un courant électrique dans un conducteur où règne un gradient de température.
La puissance Thomson absorbée - comme dans un puits où de l’énergie se perdrait - lors du passage d’un courant
d’intensité I dans un tronçon de fil dont les extrémités sont portées à des températures différant de dT s’écrit
dPT = hIdT . Le coefficient h s’exprime en V · K−1 . Il est appelé coefficient de Thomson du conducteur. Cet
effet dépend du sens du passage du courant. On convient de compter h > 0 si le passage d’un courant dans le
sens des températures croissantes s’accompagne d’une absorption d’énergie. C’est le cas du cuivre.

JR Seigne Clemenceau Nantes

3 – Exercices : 12 - Conduction - Convection - Rayonnement Sciences Physiques MP 2017-2018

1. On néglige dans un premier temps l’effet Thomson pour se consacrer à l’étude de l’effet Joule. Une
barre conductrice en cuivre calorifugée de longueur L, de section S, de conductivité électrique γ et de
conductivité thermique λ, est parcourue par un courant d’intensité I uniformément réparti. Les tempé-
ratures imposées aux extrémités sont T1 en x = 0 et T2 en x = L. La masse volumique du cuivre est µ, sa
capacité thermique massique est c. La température T (x, t) est identique en tout point d’abscisse x. Quelle
est la puissance Joule fournie à la barre entre les sections x et x + dx ? Réaliser un bilan énergétique
pour une section comprise entre x et x + dx. Trouver l’équation différentielle vérifiée par T (x, t).
2. On considère le régime permanent T (x). Montrer que l’équation différentielle vérifiée par T (x) est :
d2 T
= −KI 2 . Calculer K pour la barre de cuivre avec λ = 400 W · m−1 · K−1 , γ = 6 × 107 S · m−1 et
dx2
S = 2 mm2 . Déterminer T (x). À quelle condition la fonction T (x) passe-t-elle par un maximum entre
x = 0 et x = L ? On suppose que T2 − T1 = 100 K. Déterminer la valeur minimale I1 que doit posséder
l’intensité I pour qu’un maximum de température existe entre les extrémités du fil. La longueur est
L = 1 m.
3. On prend en compte maintenant l’effet Thomson et on se place en régime permanent dans la situation
où le courant I et les températures T1 et T2 satisfont les conditions I < I1 et T2 > T1 . Le courant
I circule dans le sens des températures croissantes. Trouver l’équation différentielle à laquelle obéit
la distribution de température T (x). Pour quelle valeur I2 de l’intensité I obtient-on un gradient de
dT T2 − T1
température uniforme = = Cte ? Pour le cuivre h = 2, 2 × 10−6 V · K−1 . Calculer I2 en
dx L
utilisant les données précédentes.

8. Amélioration de la conductivité thermique d’un liquide

On rappelle l’expression, pour une grandeur scalaire F (r, θ, ϕ) exprimée en coordonnées sphériques, du laplacien :

∂2F
   
1 ∂ 2 ∂F 1 ∂ ∂F 1
∆F = 2 r + 2 sin θ + 2 2
r ∂r ∂r r sin θ ∂θ ∂θ r sin θ ∂ϕ2
1. Un milieu conducteur thermique de conductivité λ est étudié en régime permanent. On étudie un cylindre
de centre O, d’axe Oz, de hauteur L, de surface de base circulaire de rayon L/2, entièrement empli
de ce matériau (cf. fig. 1 à gauche) ; les surfaces de base du cylindre sont portées aux températures
T (z = −L/2) = T1 et T (z = +L/2) = T2 .

L
2 z z
b b b

Figure 1 – Amélioration de la conductance thermique

Déterminer dans ce milieu la température T (z). Calculer en particulier le flux thermique Φ à travers le
cylindre étudié. On négligera tout effet de bord.
2. Dans le milieu précédent, on dispose dans le plan z = 0 une sphère, de rayon a ≪ L, formée d’un matériau
de conductivité thermique λ′ ≫ λ, de sorte que sa température soit considérée comme uniforme égale à
T1 + T2
(cf. fig. 1 à droite). On cherche alors la température dans le milieu environnant, en coordonnées
2
T1 + T2 T2 − T1
sphériques d’axe Oz, sous la forme T (r, θ) = + cos θf (r).
2 L
(a) Établir l’équation différentielle vérifiée par f (r).
(b) Résoudre cette équation ; on pourra chercher des solutions de la forme f (r) = Aru . Montrer que la
solution correspond à une combinaison linéaire des deux formes trouvées. On prendra garde à vérifier
L
les conditions aux limites pour r = a ∀θ et pour r = en θ = 0 par exemple sans oublier que a ≪ L.
2
~
(c) Exprimer le flux thermique conductif jc en tout point extérieur à la sphère.
(d) Calculer le flux thermique Φ′ à travers le disque de rayon L/2, situé dans le plan z = 0. Exprimer Φ′
en fonction de Φ, a et L.

JR Seigne Clemenceau Nantes

Sciences Physiques MP 2017-2018 Exercices : 12 - Conduction - Convection - Rayonnement – 4

3. Un liquide homogène de conductivité thermique λ contient, à raison de n par unité de volume, des
sphères de très forte conductivité thermique et de rayon a. En admettant que ces sphères soient bien
dispersées dans le milieu liquide, montrer que la conductivité thermique du milieu semble améliorée par
l’introduction de ces sphères ; calculer la conductivité thermique équivalente λe en fonction de λ, n et a.
 
Réponses : T (z) = T2 −T
L z +
1 T1 +T2
2 , Φ = λ(T1 −T4
2 )πL
, ∆T = a1 ∂T d
∂t = 0 d’où dr r
2 df (r)
dr = 2f (r), u(u + 1) = 2
donc u = 1 ou u = −2, u = 1 correspond à la première situation puisque qu’alors z = r cos θ et u = −2,
3
T (r, θ) = T1 +T
2
2
+ T2 −T B 3
L (Ar + r 2 ) cos θ, f (r = a) = 0 impose B = −Aa et l’autre
1
h condition A(1 − 8a L3 ) = i1
3 3 3 3
et donc A ≃ 1 + 8a3 , T (r, θ) = T1 +T2 + A T2 −T1 (r − a2 ) cos θ, ~jc = Aλ T1 −T2 (1 + 2a3 ) cos θ~er − (1 − a3 )~eθ ,
L 2 L r L r r
T1 −T2 8a3 T1 −T2 2
a3 3
′ 2
Φ = Φ1 + Φ2 , Φ1 = λ L 3πa (1 + L3 ) pour la sphère de rayon a et Φ2 = λ L π[ L4 + L − 3a ](1 + 8a
2
L3 )
L2 a3 8a3 4a3 8a3 3
12a3
pour le reste, Φ′ = λ T1 −T ′
L π[ 4 + L ](1 + L3 ), Φ = Φ[1 + L3 ](1 + L3 )
2
Φ′ ≃ Φ(1 + 12a
L3 ), λe ≃ λ(1 + Ntot L3 )
3
avec Ntot = n πL 3
4 d’où λe ≃ λ(1 + 3πna ).

B. Régime dépendant du temps

9. Mise en équilibre thermique, analogie
On considère la conduction thermique entre deux sphères de rayons respectifs R1 et R2 avec R1 < R2 . Entre
ces sphères l’espace est occupé par un matériau homogène et isotrope de conductivité thermique λ supposée
constante. Les sphères sont portées respectivement aux températures T1 et T2 < T1 . Le régime est supposé
stationnaire.
1. Calculer en fonction de R1 , R2 et λ la résistance thermique Rth entre les deux sphères.

Les deux sphères ont une même capacité calorifique C et ont une grande conductivité thermique de sorte
qu’à chaque instant on peut considérer les températures T1 et T2 comme uniformes. On désignera par T01
et T02 les températures initiales des sphères. On définira une constante de temps τ qui fixe l’évolution
des températures. On supposera que l’ensemble est isolé thermiquement avec le milieu extérieur.
2. Déterminer les équations qui déterminent les évolutions temporelles des températures T1 et T2 .
3. Quelles analogies peut-on faire ?

Réponses : div ~j = 0 conduit à r2 dT α R1 R2

dr = α, T = − r + β, condition aux limites α = R2 −R1 (T2 − T1 ), jcond =
4πλ RR21−R
R2
1
R2 −R1
(T1 − T2 ) et Rth = λ1 4πR 1 R2
, Rth C dT dT2
dt + T1 = T2 et Rth C dt + T2 = T1 , τ = Rth C, T1 =
1 T01 +T02
2 +
T01 +T02 2t T01 +T02 T01 +T02 2t
2 exp − τ et T 2 = 2 − 2 exp − τ , électricité circuit RC avec τ = RC.

10. Estimation de l’âge de la Terre par Lord Kelvin

On néglige la sphéricité et les sources radioactives de la planète, mais on ne se place pas en régime permanent.
On admet que la température dépend de t et de la profondeur z comptée positivement. Elle vérifie l’équation
de diffusion :
∂T ∂2T
ρcp =λ 2
∂t ∂z
où ρ est la masse volumique, cp la capacité thermique massique à pression constante et λ la conductivité
thermique.
1. Démontrer l’équation différentielle vérifiée par q (puissance surfacique) :

∂q ∂2q
=D 2
∂t ∂z
dans la quelle on notera D la diffusivité thermique D = λ/ρcp .

Au milieu du xixe siècle, Lord Kelvin a imaginé que la Terre avait été formée à une température élevée
T1 uniforme à la date t = 0. Il a proposé d’autre part qu’à cette même date, sa surface avait été soumise
instantanément à une température TS . Depuis ce temps-là, la planète se refroidirait. Lord Kelvin a
modélisé le refroidissement pour en déduire l’âge de la Terre. La densité de flux thermique est donc une
fonction de la profondeur et du temps q(z, t).
2. Dans l’hypothèse de Lord Kelvin, quelle doit être la valeur de la densité de flux thermique en z = 0
lorsque t tend vers zéro et lorsqu’il tend vers l’infini ? Quelle doit être la valeur de la densité de flux
thermique à une profondeur z non nulle lorsque t tend vers zéro et lorsqu’il tend vers l’infini ?

JR Seigne Clemenceau Nantes

5 – Exercices : 12 - Conduction - Convection - Rayonnement Sciences Physiques MP 2017-2018

3. Vérifier que la solution proposée par Lord Kelvin :

A z2
q(z, t) = − √ exp −
Dt 4Dt
où t est le temps écoulé depuis la formation de la Terre est bien la bonne. Dessiner schématiquement la
valeur absolue de la densité de flux thermique en fonction de la profondeur pour deux époques différentes.
4. Les paramètres du problème sont T1 − TS , λ, ρ et cp . On suppose que A s’exprime par :

1
A = √ (T1 − TS )α λβ ργ cδp
π
Déterminer par analyse dimensionnelle, les valeurs des exposants de cette loi.
∂T
5. Exprimer la valeur du gradient thermique en surface de la Terre . Lord Kelvin a admis que T1 − TS
∂z
−6 2 −1
était de l’ordre de 1000 à 2000 K et que D est proche de 10 m · s . Sachant que l’augmentation de
température mesurée dans les mines indiquait un gradient proche de 30 K · km−1 , quel âge de la Terre
Lord Kelvin a-t-il déduit de son modèle ?
6. Que pensez-vous de l’estimation précédente ? Quel est le ou les ingrédients que Lord Kelvin n’aurait pas
dû négliger ?

Réponses : q(z = 0, t → 0) → ∞, q(z = 0, t → ∞) → 0, q(z 6= 0, t → 0) → 0, q(z 6= 0, t → ∞) → 0 ; A en

W · m−1 , α = 1, β = 1, γ = 0, δ = 0 ; pour 1500 K : 80 millions d’années ; trop faible, radioactivité à l’intérieur
de la Terre.
11. Régime transitoire et série de Fourier
Un solide (C) la forme d’un cylindre droit à base circulaire de hauteur L, de rayon R est constitué d’un matériau
homogène et isotrope de masse volumique µ, de capacité thermique massique c et de conductivité thermique
λ supposées constantes. T désigne la température du cylindre. On appelle x la direction parallèle à l’axe du
cylindre et on suppose que T ne dépend que de x et de t. On place les deux faces extrêmes de (C) (en x = 0 et
x = L) en contact avec deux sources de chaleur (S ′ ) et (S ′′ ) de température respectives T ′ et T ′′ et on empêche
λ
tout transfert thermique par la face latérale. On posera a = . Pour t < 0 on a T ′ = T ′′ = T1 . A t = 0 on
µc
change les sources (S ′ ) et (S ′′ ) et pour t > 0 on a T ′ = T ′′ = T0 .
1. Établir l’équation de diffusion thermique dans le cylindre.
2. Donner en le justifiant la fonction T (x, t) juste avant t = 0.
On s’intéresse désormais à la fonction T (x, t) pour t > 0 et on pose θ = T − T0 .
3. Quelles sont les conditions aux limites pour θ en x = 0 et x = L ? Quelles sont les conditions initiales à
t = 0 (en fonction de x) ?
4. On cherche θ(x, t) sous la forme θ = f (t)g(x).
d2 g
— Montrer que g(x) est solution de l’équation différentielle 2 = αg où a est une constante indéterminée
dx
à ce stade des calculs.
— À l’aide des conditions aux limites, montrer que −α est positive (on posera −α = k 2 ) et ne peut
prendre que certaines valeurs dépendant d’un entier n.
— À quelle équation différentielle obéit f (t) ?
— Montrer que la solution la plus générale que l’on peut obtenir par cette méthode est :
∞
X t
θ(x, t) = Bn exp − sin kn x
n=1
τn
en donnant les expressions de τn et kn en fonction de n, L et a.
5. À l’aide des conditions initiales montrer que le calcul des coefficients Bn se ramène au calcul des coef-
ficients de Fourier d’une fonction g ∗ (x) dont on précisera la parité et la période. Calculer Bn et donner
l’expression θ(x, t).
6. Calculer le rapport rn entre l’amplitude d’un terme quelconque du développement de θ et l’amplitude
du premier terme et montrer qu’à partir d’un instant t1 dont on donnera un ordre de grandeur on peut
garder uniquement le premier terme. Donner l’allure de θ(x, t) pour t ≪ t1 , pour t > t1 puis pour t → ∞.
À partir de quel instant t2 a-t-on :

|T (x, t) − T0 |
< 10−2
T0

JR Seigne Clemenceau Nantes

Sciences Physiques MP 2017-2018 Exercices : 12 - Conduction - Convection - Rayonnement – 6

Données numériques : L = 1 m ; R = 2 cm ; µ = 9000 kg · m−3 ; c = 400 J · kg−1 · K−1 ; λ = 400 W · m−1 · K−1 ;
T1 = 370 K et T0 = 300 K.

C. Diffusion de matière
12. Neutronique
On étudie la diffusion de neutrons dans la matière fissile d’un réacteur à neutrons rapides. On appelle, en tout
point ~r du matériau, ρ (r, t) le nombre de neutrons par unité de volume et ~v (r, t) la vitesse de ces neutrons. Du
fait des collisions multiples des neutrons, on admet que v = k~v k est  uniforme  dans le milieu. On rappelle qu’en
1 ∂ 2 ∂f 1 ∂ 2 (rf )
coordonnées sphériques avec invariance en θ et ϕ, on a ∆f = 2 r = .
r ∂r ∂r r ∂r2
−−→
1. La loi de diffusion des neutrons dans le matériau fissile s’écrit ρ (r, t) ~v (r, t) = −Dv grad ρ (r, t). Préciser
le signe, l’unité de mesure et le nom du coefficient D.
2. On suppose d’abord qu’aucun neutron n’est absorbé ou émis dans le milieu. En déduire l’équation aux
dérivées partielles satisfaite par ρ (r, t).
3. Du fait des absorptions des neutrons et des réactions de fission qui en produisent,le nombre
 de neutrons
dρ
cédés au milieu par unité de volume et par unité de temps se met sous la forme = Σvρ.
dt m
Préciser l’unité de Σ et le signe de Σ selon que la fission l’emporte sur l’absorption ou vice-versa.
Écrire alors l’équation aux dérivées partielles satisfaite par ρ ; on supposera ici Σ > 0.
4. Le milieu fissile est sphérique de rayon R ; on écrira donc, compte tenu de la symétrie, que ρ ne dépend
que de la distance r au centre O de la sphère et du temps t.
(a) Déterminer ρ(r) en régime permanent. On justifiera les conditions aux limites utilisées.
Montrer que ce régime permanent n’existe que si R < Rc et exprimer la valeur critique Rc .
(b) Si R < Rc , ou R > Rc , montrer qu’il existe des solutions de la forme ρ(r, t) = exp (γt) ϕ(r) et
interpréter le résultat.

Réponses : D en m longueur caractéristique, D > 0 car la diffusion s’effectue des zones les plus concentrées vers
1 ∂ρ
les moins concentrées, ∆ρ = Dv 2 −1
∂t , a = Dv diffusivité en m ·s , Σ > 0 la fission l’emporte, Dv∆ρ+ Σvρ = ∂ρ ∂t ,
d2 X
Σ en m−1 , dr 2 + = 0 avec X = rρ, ρr=0 = ρ0 bornée et maximale dρ
Σ
DX dr = 0, tous les neutrons créés
q r=0 q
R Σ d2 Y Σv−γ
sortent vρr=R 4πR2 = Σv 0 ρ(r)4πr2 dr, ρ(r) = ρ0 sinc D r, ρ(r) > 0 d’où 0 < R < π D
R
Σ , dr 2 + Dv Y = 0
sh αr
avec Y = rϕ, si γ < 0 les neutrons disparaissent, si γ < σv alors ϕ(r) en sinc sinon en αr si γ > 0 alors le
nombre de neutrons diverge : explosion.

13. Colonie de bactéries

On s’intéresse à la croissance d’une colonie de bactéries par division et diffusion. La densité (nombre de bactéries
par unité de longueur) n(x, t) ne dépend que d’une coordonnée x et du temps t.
1. On suppose que la loi d’évolution des bactéries pour la diffusion est identique à la loi de diffusion
thermique unidimensionnelle. Rappeler l’expression de cette loi, on notera D le coefficient de diffusion
encore appelé diffusivité. Sachant que chaque bactérie change aléatoirement de sens tous les T = 1 s,
donner une estimation de D lorsque la vitesse des bactéries est v = 10 µm· s−1 . À l’aide d’un raisonnement
dimensionnel, trouver le temps au bout duquel une petite colonie remplit un récipient de longueur L =
10 cm, commenter brièvement le résultat obtenu.
2. Une bactérie se divise en deux en moyenne tous les τ = 1200 s. Faire un bilan de nombre de bactéries sur
une tranche comprise entre x et x+dx et un intervalle de temps compris entre t et t+dt. Trouver l’équation
de diffusion modifiée à laquelle satisfait la densité de bactérie. Quelles sont les solutions stationnaires de
l’équation obtenue ? Ont-elles un sens physique ? Déterminer l’évolution d’un profil de densité indépendant
de x et valant n(x, t = 0) = n0 initialement.
3. Tous les τ = 1200 s en moyenne, une bactérie meurt avec une probabilité proportionnelle au nombre
moyen de ses voisins dans un intervalle de longueur a. Modifier le bilan de particules établi à la question
précédente en tenant compte de la mort des bactéries et montrer que :

∂n ∂2n
= D 2 + d1 n − d2 n2
∂t ∂x
où d1 et d2 sont des coefficients dont on précisera les expressions, les dimensions et les valeurs numériques
(on prend a = 10 µm). Quelles sont les solutions stationnaires et uniformes de l’équation précédente ? À
quoi correspondent-elles physiquement ? On les notes n1 et n2 avec n1 < n2 .

JR Seigne Clemenceau Nantes

7 – Exercices : 12 - Conduction - Convection - Rayonnement Sciences Physiques MP 2017-2018

4. On s’intéresse au bord de la colonie et on cherche des solutions de l’équation différentielle de la question

3 sous la forme n(x, t) = f (x − ct), telle que :

lim f (x) = n2 lim f (x) = n1

x→−∞ x→+∞

Que signifie ce choix ? Montrer que l’équation différentielle qui détermine f se met sous la forme :

e1 f ′′ = e2 f ′ + e3 f + e4 f 2

où les ei sont des coefficients dont on donnera l’expression. Reconnaı̂tre dans cette équation le mouvement
d’une particule de masse D soumise à une force de frottement visqueux et à un potentiel à préciser. Faire
un Zbilan des travaux des forces entre les positions n2 et n1 et en déduire une expression de c en fonction
de f ′2 (u)du. Dans quel sens la colonie se déplace-t-elle ? Un récipient se remplirait-il plus vite ou moins
vite que dans la question 1 ?
2 2
Réponses : ∂∂xn2 = D 1 ∂n 2
∂t , D en m · s
−1
, D ≃ v 2 T = 10 × 10−10 m2 · s−1 , ∆t = LD = 108 s donc trois ans,
∂2n
c’est trop long, la diffusion n’est pas le seul mécanisme, ∂n n
∂t = D ∂x2 + τ , solutions sinusoı̈dales, solutions
stationnaires sans sens physique car la division entraı̂ne forcément une augmentation du nombre de particules,
solutions indépendantes de x n(t) = n0 exp τt diverge cohérent avec l’observation précédente, Nmoy ≃ na,
∂2n
probabilité en na ∂n n a 2 1
τ , ∂t = D ∂x2 + τ − τ n , d1 = τ = 8, 3 × 10 s , d2 = τa = 8, 3 × 10−9 m · s−1 , solutions
−4 −1
1
stationnaires et uniformes n1 = 0 aucune bactérie, n2 = a la mortalité est exactement compensée par la
df d2 f
division, source de bactérie en −∞ qui meurent toutes avant d’arriver en +∞, ξ = x − ct, f ′ = dξ et f ′′ = dξ 2 ,
2
′2 ′2
D ddξf2 = − 1c df 1 a 2 1 2 a 3 1
dξ − τ f + τ f , Epot = 2τ f − 3τ f , ∆Ec + ∆Epot = Wfrott avec ∆Ec = 2 D(f n2 − f n1 ) = 0 d’où
2
c = 6τ a2 f ′ (ξ)dξ, plus vite car avec la solution f (x − ct) la propagation s’effectue dans le sens x croissant
R

avec des durées proportionnelles aux distances.

D. Soleil - Terre
14. Dilution du rayonnement solaire
La Terre et le Soleil sont considérés comme des corps noirs aux températures respectives de 300 K et 5800 K.
Leurs rayons sont de 6 400 km et 700 000 km et ils sont distants de 150 millions de kilomètres.
1. Calculer les puissances rayonnées par la Terre et le Soleil.
2. Calculer l’angle solide Ω sous lequel la Terre est vue depuis le Soleil. Le rayonnement du Soleil étant
supposé isotrope, en déduire la puissance reçue par la Terre due au rayonnement solaire. Conclure.

πR2T
Réponses : 7, 5 × 1016 W et 4 × 1026 W ; Ω = d2 , 2 × 1017 W, situation de quasi équilibre (imprécision des
calculs).

15. Température d’équilibre de la Terre et effet de serre

Le tableau 1 donne quelques valeurs numériques et constantes astronomiques.

Grandeur physique Symbole Valeur numérique

Constante de Stefan σ 5, 67 × 10−8 W · m−2 · K−4
Température de surface du Soleil TS 6000 K
Distance de la Terre au Soleil a 150 × 106 km
Rayon solaire RS 700 000 km

Table 1 – Tableau de valeurs numériques

1. Dans un premier modèle, on considère que la Terre réfléchit la portion A (A porte le nom d’albédo)
du rayonnement solaire. Le Soleil et la Terre sont assimilés à des corps noirs de températures TS et T1
uniformes.
Déterminer la température T1 à l’équilibre. On donne A = 0, 35. Faire l’application numérique.
2. On tient compte de plus des propriétés de la troposphère, couche d’épaisseur e faible devant le rayon
terrestre RT , de température T . La troposphère rayonne comme un corps noir, mais elle absorbe une
fraction α du rayonnement solaire. La température de la Terre est maintenant T2 .
RS
(a) Exprimer l’équation décrivant l’équilibre thermique de la Terre, reliant A, α, T , T2 , TS et γ = .
a

JR Seigne Clemenceau Nantes

Sciences Physiques MP 2017-2018 Exercices : 12 - Conduction - Convection - Rayonnement – 8

(b) Exprimer l’équation décrivant l’équilibre thermique de la troposphère, reliant A, α, T , T2 , TS et γ.

(c) Exprimer la nouvelle température de la Terre en fonction de celle déterminée précédemment et de α.
(d) Faire l’application numérique si α = A. Commenter.

Réponses : (1 − A)σTS4 ( RaS )2 πRT2 = 4πRT2 σT14 , T1 = TS [ 1−A RS 2 1/4

4 ( a ) ] = 260 K, inférieure à la valeur constatée :
285 K ; σT2 4πRT = σT 4πRT + (1 − α)(1 − A)σTS ( a ) πRT , T24 = T 4 + (1 − α)T14 , α(1 − A)σTS4 ( RaS )2 πRT2 +
4 2 4 2 4 RS 2 2

σT24 4πRT2 = 2σT4 4πRT2 , αT14 + T24 = 2T 4, T2 = (2 − α)1/4 T1 = 295 K, plus agréable . . .

E. Objets terrestres en régime indépendant du temps

16. Plaque chauffée
On dispose d’une plaque de cuivre, de surface S = 12 cm2 , reliée à un circuit électrique. Il passe un courant
d’intensité I à travers la section s la plaque, on note que s < S. On réalise successivement 3 expériences.
Dans le première expérience, on éclaire un côté de la plaque pendant une durée ∆t1 = 140 s en l’absence de
courant (I = 0). On observe une élévation de température ∆θ1 = 11 ◦ C. Dans la deuxième expérience, la plaque
est plongée dans l’obscurité et on impose une intensité I = 2, 23 A pendant une durée ∆t2 = ∆t1 = 140 s.
On observe la même élévation de température que dans la première expérience ∆θ2 = ∆θ1 = 11 ◦ C. Dans la
troisième expérience, on se place dans les même conditions que dans la première expérience pendant une durée
de 70 s. Après ces 70 s, on place la plaque de cuivre à l’obscurité et on mesure, à l’aide d’un ohmmètre, la
résistance de la plaque et on trouve R = 0, 18 Ω.
1. Déterminer la moyenne de la densité de courant de rayonnement jray (ou encore nommé vecteur de
Poynting) reçu par la plaque.
2. Estimer la perte de masse du Soleil pour une année terrestre. On donne la distance Terre-Soleil d =
1, 5 × 1011 m, la masse MS = 2 × 1030 kg du Soleil et enfin le rayon du Soleil Rs = 7 × 108 m.

Réponses : dans la première expérience la puissance reçue est P1 = jray S = mc dθ dt si m est la capacité thermique
j S j S
massique du cuivre et m la masse de la plaque, on a θ = ray mc t + θ 0 , d’où ∆θ = raymc ∆t. Pendant la seconde
2
expérience, la puissance reçue est Pélec = RI 2 et on a ∆θ = RI mc ∆t. La mesure de la résistance électrique
dépend de la température mais la mesure est effectuée à la température qui correspond à la moitié de l’élévation
de la température pour minimiser l’influence de θ. La puissance électrique est Pélec = RI 2 = 0, 9 W. Cela
signifie que jray S = 0, 9 W. On en déduit que jray = 750 W · m−2 . La dilution du rayonnement solaire fait que
jray 4πd2 = PS représente la puissance émise par le Soleil. Par la formule d’Einstein, Ps = − dm 2 ∆m 2
dt c = − ∆t c .
16
On trouve ∆m = 7, 4 × 10 kg perdus par an, ce qui est relativement faible comparativement à la masse du
Soleil.

17. Effet de serre

On étudie l’effet de serre produit par l’interposition d’une vitre au-dessus d’une plaque qui reçoit le rayonnement
solaire. La plaque est noircie et assimilée à un corps noir. Le verre est supposé totalement transparent au
rayonnement solaire (sauf à la question 3) où l’on tient compte du rayonnement solaire réfléchi. La vitre est en
revanche totalement absorbante pour le rayonnement infrarouge émis par la plaque qui absorbe le rayonnement
solaire. On désigne par ϕS le flux solaire surfacique supposé arriver normalement à la vitre et à la plaque.
1. On suppose l’équilibre radiatif de la plaque et de la vitre. Écrire les équations exprimant ces équilibres
et en déduire la température T de la plaque. On donne : ϕS = 600 W · m−2 . Calculer T et la température
T1 de la vitre.
2. Reprendre la question précédente dans le cas de deux vitres puis de n vitres.
3. Chaque vitre réfléchit une fraction r de l’énergie solaire incidente. On néglige toujours l’absorption du
rayonnement solaire par les vitres. Montrer qu’il existe une valeur nm optimale du nombre de vitres.
Application numérique : r = 0, 08.

Réponses : ϕS + ϕV = ϕP et ϕP = 2ϕV d’où σT 4 = 2ϕS et σT14 = ϕS , ϕS + ϕV 1 = ϕP ϕP + ϕV 2 = 2ϕV 1 et

2ϕV 2 = ϕV 1 d’où σT 4 = 3ϕS , ensuite σT 4 = (n + 1)ϕS , cela diverge, σT 4 = (1 − r)n (n + 1)ϕS optimale pour
1
nm = − ln(1−r) − 1, nm ≃ 11.

JR Seigne Clemenceau Nantes