Durée: 5 HEURES

Aucun document 'est autorisé.

Pour les épreuves d'admissibilité, l'usage de calculatrices électroniques de poche à alimentation autonome, non imprimantes et sans document d'accompagnement, est autorisé, une seule à la fois étant admise sur la table ou le poste de travail, et aucun n'échange n'est autorisé entre les candidats.

L'épreuve se décompose en deux problèmes, à traiter de préférence dans l'ordre proposé. Ils concernent tous deux la modélisation mathématique des phénomènes de combustion et d'explosion.

Le premier problème se focalise plus particulièrement sur les réactions explosives. Leurs principales caractéristiques ont été étudiées au milieu du XIXeme siècle par Van’t Hoff et Bunsen. Elles présentent la particularité d'avoir un taux de réaction qui dépend de manière non-linéaire et forte de la température (par opposition aux réactions lentes dont la dépendance est faible). Cette propriété de forte dépendance est due à une grande énergie d'activation de la réaction de combustion; elle se couple en général à une forte puissance thermique de la réaction. Dans ces conditions, un phénomène d'autocatalyse thermique peut se produire, provoquer l'emballement de la réaction et une très forte augmentation de la température; ce phénomène est appelé explosion thermique. Divers travaux expérimentaux ont alors conduit N.N. Semenov à développer la théorie de l'explosion thermique dans les années 1920 pour le cas d'un réacteur homogène. Frank-Kamenetskii a ensuite repris cette théorie dans le cas non-homogène et proposé un traitement exhaustif dans un certain nombre de cas.

Le but du problème I est d'introduire ces théories et de saisir les mécanismes de base de couplage entre le terme source thermique issu de la réaction chimique, la conduction thermique et les pertes thermiques au bord des domaines considérés. Le but est d'obtenir les conditions critiques d'explosion, i.e. les valeurs limites de certains paramètres qui déterminent si le phénomène d'explosion a lieu ou non. Finalement, le dernier point que nous aborderons sera le couplage entre l'explosion thermique et la convection naturelle en présence du champ de gravitation. Ces études sont importantes et trouvent de nombreuses applications dans la détermination des conditions de sécurité pour le stockage de matières dangereuses.

Dans le problème , nous nous contentons de considérer des phénomènes explosifs localisés mais nous ne traitons pas leur propagation car il s'agit d'un problème essentiellement plus difficile qui nécessite de prendre en compte les effets de densité volumique de masse et de pression puisque des ondes de chocs ou de détonations peuvent se former. On s'intéresse par contre, dans le problème II, à un modèle simple de déflagration à pression constante qui décrit la propagation d'une flamme à faible vitesse par rapport à la vitesse du son. L'étude fondatrice dans ce domaine date de 1928 avec les travaux de Zeldovich sur un modèle que nous allons étudier.

Notation : Dans toute l'épreuve, on utilisera les notations suivantes pour les dérivées, désigne la dérivée temporelle à l'instant de la fonction et désignent respectivement la dérivée partielle par rapport au temps, la dérivée partielle par rapport à la variable spatiale et la dérivée seconde par rapport à la variable , au point du champ .

Note : Le premier problème requiert l'utilisation du théorème de Cauchy-Lipschitz. Il est rappelé dans l'énoncé.

Nous allons dans un premier temps considérer les équations qui régissent l'écoulement d'un mélange gazeux réactif compressible dans un domaine bi-dimensionnel . Le candidat ne devra pas être impressionné par la complexité du système car les quatre premières parties se concentrent sur des versions très simplifiées pour lesquelles le mouvement hydrodynamique du mélange est supposé inexistant (la vitesse du fluide est zéro, uniformément en temps et en espace).

Ce type d'écoulement est régi par un système d'équations aux dérivées partielles pour lequel le vecteur des variables est , où désigne la transposition et où est la densité volumique de masse partielle de la espèce, variant de 1 à , la densité volumique de masse totale, et les composantes de la vitesse du fluide dans les directions (horizontale) et (verticale) et enfin , la température du fluide. Ces quantités dépendent donc du temps et des coordonnées d'espace : . Il n'est pas nécessaire de détailler ce système d'équations, mais il faut souligner qu'il doit être accompagné par des relations thermodynamiques et une loi d'état. Cette loi d'état relie la pression (qui apparaît dans les équations pour la vitesse à travers son gradient) à la densité volumique de masse totale, à la température et à la composition du mélange gazeux. On peut alors écrire une équation sur le champ de température qui prend la forme :

où est la capacité thermique massique à pression constante, représente la puissance des forces de viscosité, , le vecteur densité volumique de courant thermique et désigne sa divergence, , la puissance thermique (terme source) due aux réactions chimiques.

La pression vérifie la loi d'état des mélanges de gaz parfaits :

où est la masse molaire de la espèce, la masse molaire moyenne du mélange, la constante universelle des gaz, la fraction massique de l'espèce et on a ,
par définition de . Cette loi est la simple généralisation de la loi d'état des gaz parfaits à un mélange gazeux.

Le système de base que nous venons d'évoquer est souvent utilisé dans l'approximation de Oberbeck-Boussinesq qui considère que la pression et la densité volumique de masse sont quasiment constantes et (le signe signifie "très petit devant" et indique une séparation d'échelles. Ce genre d'hypothèse implique par exemple que , i.e. que cette expression peut être approchée formellement par 1 , ou encore que les variations de pression ou de densité volumique de masse par rapport à un état de référence ou sont négligeables).
I. 1 On suppose que le mélange gazeux est initialement à la température , à la densité et à la pression , uniformément en espace. On néglige les différences de masse molaire, c'est-à-dire . Montrer que la loi d'état (I.2) se réduit alors à la loi d'état d'un gaz parfait simple .
I. 2 Montrer, en utilisant cette loi d'état, que l'approximation de Oberbeck-Boussinesq est valable pour de faibles variations de température, c'est-à-dire .

Avant de revenir sur les équations dans la dernière partie du problème, nous allons commencer par comprendre les mécanismes de base sur des modèles simplifiés.

Le but de cette partie est de faire apparaître, dans la cadre d'une modélisation minimale, l'échelle de temps du problème. On se place dans le cadre d'une chimie simple, c'est-à-dire, d'une chimie constituée d'une unique réaction globale irréversible "Fuel Produits". Le degré d'avancement de cette réaction est donné par la variable fraction massique de fuel : ; elle vaut 1 à l'instant initial quand le fuel n'a pas encore commencé à réagir et 0 quand tout le fuel présent dans le réacteur a été consommé en produits et en source thermique. La température du réacteur à l'instant initial est notée et évolue du fait de la réaction chimique : on note la température à l'instant . L'évolution de la fraction massique et de la température est décrite par un système de deux équations différentielles ordinaires. Elles font intervenir la vitesse de la réaction qui dépend de manière non-linéaire de la température. Ce modèle, donné pour décrire l'évolution du réacteur, est valable pour et ; il s'écrit :

est le facteur de fréquence (homogène à l'inverse d'un temps), est l'énergie d'activation de la réaction, la température du mélange lorsque tout le fuel est consommé; cette température est associée à une enthalpie de réaction, à la capacité thermique du milieu et à une quantité de fuel présente dans le réacteur. On suppose dans ce problème que la réaction est exothermique et donc que . Par ailleurs et sont des constantes
strictement positives. On ne prend en compte, dans ce modèle, que le chauffage du milieu par la réaction et l'on suppose que le système est isolé thermiquement de l'extérieur; on parle alors d'explosion thermique adiabatique.

Note : On donne ici un théorème utile pour la suite sous la forme d'une adaptation du théorème d'existence et unicité de solution de Cauchy-Lipschitz. On considère un système d'équations différentielles ordinaires autonome:

avec comme donnée initiale du problème de Cauchy et .
THÉORÈME : On considère un ensemble compact de sous la forme d'un produit d'intervalles fermés , avec . On suppose que applique dans , que est , Lipschitz sur et que appartient à . On fait par ailleurs l'hypothèse que dans le cas où la donnée initiale appartient à , la frontière de , le vecteur est soit nul (point singulier) soit non dirigé vers l'extérieur du domaine . Alors il existe une unique solution maximale , du système (II.3) de régularité sur un intervalle de temps , avec éventuellement , telle que . Dans le cas où est fini, on a , quand .

COROLLAIRE: Sous les hypothèses du théorème précédent, si deux solutions et du système (II.3) sont égales en un instant , même si (par exemple si ou si est un temps d'existence maximal fini), alors ces deux solutions sont égales sur l'intersection des intervalles de temps d'existence des deux solutions maximales.
L'idée sous-jacente est que sous des hypothèses de régularité de , deux solutions distinctes du système (II.3) ne peuvent pas "se croiser" dans le domaine .
II. 1 Le but de cette question est de montrer l'existence d'une solution pour tout temps et d'obtenir des renseignements sur la structure qualitative de la solution du système (II.1-II.2).
II.1.a Montrer que, du système (II.1-II.2), on peut déduire une équation sur une quantité , combinaison de et , qui se conserve au cours du temps. Donner la valeur de .
II.1.b En déduire une équation autonome sur la variable

et expliciter la fonction .
II.1.c En appliquant le théorème donné ci-dessus, montrer qu'il existe une unique solution maximale sur l'intervalle de temps . Montrer que la fonction est monotone strictement décroissante, puis que . En déduire que est strictement positive (Indication : supposer que est fini, montrer que la seule possibilité est que quand , considérer la solution constante égale à zéro et conclure en utilisant le corollaire ci-dessus).
II.1.d Montrer que est monotone croissante bornée supérieurement par .
II.1.e Déduire de II.1.c que la limite de en est 0 .
II.1.f Représenter graphiquement le comportement qualitatif de la solution.
II.1.g Proposer une équation différentielle autonome décrivant l'évolution de la température , ne faisant pas intervenir , sous la forme et donner une expression explicite pour la fonction .

La structure de la non-linéarité dans cette dernière équation sur la température ne permet pas une résolution analytique. Nous allons donc utiliser deux caractéristiques des réactions de combustion, introduites sous la forme de deux hypothèses, et expliciter le comportement qualitatif du système.
[H1] La réaction satisfait l'hypothèse des grandes énergies d'activation si est très grand devant 1 , ce qui implique :

La réaction satisfait l'hypothèse d'une grande enthalpie de réaction si et donc :

où est appelée la température de Frank-Kamenetskii et est notée .
II.1.g On se place dans le cas où . Représenter sur deux graphiques les fonctions pour et , et pour . (Afin de pouvoir comparer les deux non-linéarités, on divisera la fonction par sa valeur maximum sur l'intervalle ). En comparant les valeurs de pour et pour , dans les deux situations et en observant les graphiques, expliquer ce que signifie "une forte non-linéarité dépendant de l'énergie d'activation".
II. 2 Nous allons maintenant étudier un modèle approché pour le système (II.1-II.2).
II.2.a Adimensionner l'équation vérifiée par la température et écrire pour cela une équation sur la température réduite . Cela revient à utiliser comme température de référence la température initiale et comme échelle, la température de Frank-Kamenetskii.
II.2.b Montrer que, sous les hypothèses et , et tant que peut être approchée par une fonction , solution de l'équation différentielle :

avec la condition initiale . Le paramètre est appelé temps d'induction et l'on donnera son expression explicite en fonction de et .
II.2.c Dans l'équation précédente, on a remplacé sous l'hypothèse [ ], l'argument de l'exponentielle par . Reprendre cette approximation dans les variables dimensionnées et expliquer pourquoi on l'appelle "linéarisation de Frank-Kamenetskii".
II.2.d On pose ; déduire de (II.5) l'équation différentielle sur . Résoudre l'équation correspondante et donner une formule explicite pour .
II.2.e Décrire le comportement "explosif" de cette solution. Tracer l'évolution temporelle de la solution , pour , et ceci selon deux échelles différentes pour l'axe des ordonnées : la première est telle que le maximum de l'axe des ordonnées est 5 , et la seconde telle que le maximum est (on prendra égal à 50 pour l'exemple). Décrire le comportement aux deux échelles et commenter.
II.2.f On appelle le temps pour lequel atteint . Exprimer la différence relative entre et , ainsi que le pourcentage de que représente .
II. 3 Il s'agit de montrer que l'équation simplifiée (II.5) est un "bon modèle" et d'en déduire le comportement qualitatif général de la solution de (II.1-II.2).
II.3.a Montrer que la condition de validité n'est par contraignante en utilisant la fin de la question II. 2 et en déduire le comportement du système de départ sur l'intervalle .
II.3.b Calculer la température correspondant au taux de réaction maximum, en valeur absolue, pour le système (II.1-II.2), sachant que le taux de réaction est le second membre de l'équation sur et vaut donc . Comment se situe-t-elle par rapport à ?
II.3.c Quel scénario peut-on déduire sur l'évolution du système (II.1-II.2) dans la mesure où le taux de réaction maximum est atteint dans un proche voisinage de ?
II.3.d Caractériser l'évolution du système dans le plan des variables dimensionnées et (appelé aussi plan de phase). Tracer par ailleurs les fonctions et . Décomposer l'évolution du système en trois étapes que l'on précisera.
II.3.e Que se passe-t-il pour la consommation de fuel (caractérisée par la fraction massique de fuel ) dans l'intervalle ? Interpréter l'hypothèse [H2].
II.3.f Que devient l'évolution du système quand l'énergie d'activation devient de plus en plus grande?
II.3.g On choisit une application numérique avec . On donnera la valeur de la température de Frank-Kamenetskii et on refera les diagrammes de la question II.3.d en marquant les points particuliers suivants tels que

Commenter. (Remarque : les données numériques ne permettent pas de calculer explicitement; les temps seront donc exprimés en pourcentage de )

Dans la partie précédente, l'absence de perte thermique au bord du réacteur engendre systématiquement l'explosion thermique car la température ne peut qu'augmenter et atteindre (température de l'ordre de 3000 K pour les réactions de combustion). Dans ces conditions, le modèle cesse d'être valide quand on arrive à de fortes températures dans des délais courts, car d'autres phénomènes rentrent alors en jeu (création d'ondes de pression, etc...). Nous levons ici l'hypothèse d'adiabaticité, mais conservons l'homogénéité spatiale. Dans un modèle plus réaliste que dans la partie précédente, les pertes thermiques à la frontière peuvent empêcher l'explosion si elles compensent le terme source issu de la réaction chimique. Il s'agit de caractériser les conditions d'explosion, c'est-à-dire dans quels cas les pertes thermiques à la frontière empêchent l'explosion et dans quel cas elles ne suffisent pas à limiter la croissance de la température du système.

Un modèle simple de perte, dans le cas d'un réacteur fermé, est de fixer la température à la paroi à la température de réference . Le réacteur peut alors être modélisé par le nouveau système, pour :

où est un temps de transfert à la paroi et où l'on prend les mêmes données initiales que dans la partie II. On suppose que les deux hypothèses [H1] et [H2] sont vérifiées.
III. 1 Nous nous intéressons brièvement au comportement qualitatif général des solutions du système pour et .
III.1.a Appliquer le théorème donné dans l'énoncé et montrer l'existence d'une solution maximale sur l'interval .
III.1.b Montrer que sur l'intervalle (indication : écrire une équation .
III.1.c Montrer que la quantité introduite dans la partie II n'est plus conservée au cours du temps par (III.1-III.2) en présence de pertes thermiques. Montrer qu'elle est majorée par la valeur obtenue pour dans le cas adiabatique de la partie II.
III.1.d Montrer que .
III.1.e Montrer, en utilisant le même type d'argument que dans la partie II, que et que (Indication : utiliser et montrer que la seule possibilité, si est fini, est que et quand . Construire alors une solution telle que et qui passe par le point ( ) en ).
III.1.f Déduire des questions qui précèdent que tend vers zéro à l'infini et par conséquent que tend vers à l'infini.
III.1.g Donner le comportement qualitatif du système. Tracer la forme de et .
III.1.h Adimensionner les équations (III.1-III.2) en utilisant le même adimensionnement que dans la partie II, et en prenant comme temps adimensionné (donner en fonction de et ).
III. 2 Comme la quantité n'est plus conservée, on ne peut plus se ramener à une seule équation et il subsiste alors le problème de la présence de la fraction massique dans l'équation sur qui décrit l'influence de la consommation de fuel sur la vitesse de la réaction chimique.
III.2.a Montrer que, sous les hypothèses et , l'équation sur peut être approchée par :

où et donner l'expression de .
III.2.b Sur quel intervalle de temps peut-on considérer que l'équation précédente est représentative de la physique du système (III.1-III.2) ? Quelle hypothèse implique son utilisation sur la consommation de fuel dans cet intervalle?
III.2.c Représenter sur un diagramme les graphes des deux fonctions et en fonction de pour diverses valeurs de . Etudier les points stationnaires de l'équation (III.3) définis par , en fonction du rapport Proposer trois scénarios. On donnera en particulier la valeur critique de ainsi que la température du point stationnaire correspondant. Interpréter physiquement les trois cas.
III.2.d Relier , en variables dimensionnées, à la température de Frank-Kamenetskii.
III.2.e Classifier et caractériser les diverses dynamiques possibles pour l'équation (III.3) en fonction de rapport . On étudiera en particulier les limites et .
III.2.f Que dire de la stabilité des points stationnaires quand ils existent (un point stationnaire est dit stable pour l'équation différentielle (III.3) si l'évolution d'une solution, dont la donnée initiale est prise dans un voisinage du point, est la convergence vers cet état stationnaire et instable dans le cas contraire) ?
III. 3 Parmi les divers scénarios proposés, certains correspondent au comportement du système (III.1-III.2) et d'autres non. Nous allons étudier ce point dans cette sous-partie.
III.3.a Montrer que cette correspondance est effective dans le cas et faire le lien avec la partie II.
III.3.b Pour , montrer que l'existence d'un état stationnaire vers lequel converge la solution de (III.3) implique une borne sur la température réduite pour le système (III.1-III.2) et garantit la non-explosion.
III.3.c Dans le cas où il y a explosion pour l'équation (III.3), c'est-à-dire pour , on estime , le temps d'induction non-adiabatique, par:

Justifier la borne d'intégration en utilisant les hypothèses.
III.3.c On cherche la limite de ce temps quand tend vers par valeur inférieure. Pour cela on pose et et l'on suppose et donc
(cette dernière estimation, que l'on admettra, vient du fait que lorsque , la température s'établit en un temps de l'ordre de 1 autour de et reste dans ce voisinage pendant une grande partie de la dynamique). Donner une estimation de en fonction de , pour petit et montrer que tend vers l'infini quand tend vers 0 (indication : on écrira l'intégrale en fonction de et , puis on remplacera par puisque supposé petit et enfin on utilisera l'égalité ). En déduire la dynamique de l'explosion quand l'on s'approche des conditions critiques par valeur inférieure, et la représenter sur un diagramme .
III.3.d Que conclure sur le modèle (III.3) au voisinage des conditions critiques en terme de l'hypothèse : "on néglige la consommation de fuel et l'on suppose sur
III. 4 En fait, la restriction sur l'applicabilité du modèle quand on s'approche des conditions critiques n'a d'effet que dans un très petit voisinage inférieur de .
III.4.a Présenter un diagramme des diverses évolutions temporelles possibles de la température dans le cas général du système (III.1-III.2).
III.4.b Donner l'évolution dans le plan de phase ( ) et faire le lien avec la partie II, afin de synthétiser les scénarios possibles.
III.4.c Conclure cette étude.

Nous allons maintenant considérer le cas inhomogène pour lequel les champs de température et fraction massique dépendent non seulement du temps, mais aussi de la variable d'espace. On se place dans le cadre particulier d'une bande infinie dans la direction horizontale et de hauteur dans la direction . On suppose le problème homogène dans la direction ce qui revient donc à travailler sur un problème monodimensionnel présenté sur la Figure 1.

On considère alors le système de deux équations aux dérivées partielles suivant :

où désigne la dérivée partielle seconde en la variable d'espace, la dérivée partielle en temps, avec les conditions initiales et , pour tout et les conditions aux limites et . La condition aux limites sur impose la non pénétration du fuel dans la paroi du tube et le fait d'imposer la température au bord induit des pertes thermiques. On peut alors montrer qu'il existe une solution régulière de ce système pour , telle que , et telle que la fonction est concave pour tout .

IV.1.a Ecrire une équation d'évolution sur la moyenne de sur l'intervalle et sur , celle de .
IV.1.b En utilisant les conditions aux limites, le signe de , et donc le signe du gradient de aux limites en et , évaluer les champs de température et de fraction massique à la limite .
IV.1.c Faire le lien avec le modèle homogène de la partie III.
IV. 2 On introduit un temps de diffusion , ainsi que , son rapport au temps d'induction introduit dans la partie II. En utilisant comme échelle de longueur, en notant , et défini comme précédemment, donner la forme adimensionnée de l'équation aux dérivées partielles sur lorsque l'on néglige la consommation de fuel et que l'on prend comme nous l'avons fait auparavant. La mettre sous la forme

où on explicitera et justifiera le calcul de .
IV. 3 Donner une expression analytique implicite de l'état stationnaire associé à l'équation (IV.3) (indication : remarquer que l'équation différentielle sur en admet une intégrale première en multipliant par la dérivée, utiliser la symétrie et la concavité pour caractériser en au centre du réacteur. La seconde intégration se fait par changement de variable adéquat que l'on justifiera) ainsi qu'une équation du type

reliant à la température maximum au centre.
IV. 4 Tracer et donner la forme de cette fonction ; montrer qu'elle admet un maximum et un paramètre critique associé. En déduire, comme dans la partie III, trois scénarios, sur l'existence ou non de profils de température réduite stationnaires, suivant si ou .
IV. 5 Dans le cas , combien existe-t-il de profils stationnaires? En se référant à la partie III, que peut-on prédire sur leur stabilité?

Nous n'envisageons pas dans cette partie une étude complète des phénomènes de couplage entre l'explosion thermique dans une bande, comme dans le cas précédent, et la convection naturelle qui va se déclencher lorsque la température au centre dépasse un seuil de stabilité. Ce genre d'étude demande un ensemble de simulations numériques sur ordinateur et une investigation paramétrique. Dans ce contexte, le modèle à utiliser est un couplage entre une équation sur la température en variables adimensionnées :

où est la variable adimensionnée dans la direction , et les équations de Navier-Stokes pour un fluide incompressible sur ( ) avec un terme source, dans l'équation de la vitesse verticale adimensionnée , qui inclue les variations de densité volumique de masse dues au faibles variations de température.
V. 1 Expliquer, à la lumière des parties précédentes, pourquoi un modèle qui néglige la consommation du fuel, suppose de faibles variations de température et considère l'équation (V.1) permet de rendre compte d'une grande partie des phénomènes d'explosion pour un modèle plus complet quand les hypothèses et sont vérifiées.

Dans ce second problème, nous allons aborder la théorie de la propagation de flammes planes laminaires. Ces objets monodimensionnel peuvent se voir, dans certains cas simples, comme des structures auto-similaires (c'est-à-dire dont le profil spatial est invariant au cours du temps) se propageant à vitesse .

Nous allons aborder ces structures que l'on appelle "ondes progressives" pour des systèmes d'équations aux dérivées partielles, en dimension un d'espace, qui modélisent le couplage instationnaire entre des phénomènes de conduction thermique et des phénomènes de réaction chimique. Nous allons considérer une chimie simple du type "Fuel Produits"; l'évolution physique peut alors être décrite par celles de la fraction massique de fuel (qui donne le degré d'avancement de la réaction) et de la température :

pour et . On rappelle que désigne la dérivée partielle seconde en la variable d'espace la dérivée partielle en temps; désigne le champ de fraction massique de fuel et le champ de température. Par ailleurs, est le facteur de fréquence homogène à l'inverse d'un temps; est une fonction monotone de qui vaut 0 en . La température du mélange où tout est brûlé de manière adiabatique est notée , comme dans le Problème I. On suppose dans ce problème que la réaction est exothermique et donc que . Par ailleurs est une constante strictement positive. On remarquera que le coefficient de diffusion de masse et conduction thermique sont pris égaux pour simplifier l'analyse menée dans ce problème.

Nous allons donc chercher des profils et et une vitesse ( fait partie des inconnues du problème), tels que et et soient solution de (I.1-I.2). Ils nous faut aussi imposer des conditions aux limites et c'est sur les comportements à l'infini que nous les fixons en supposant que , et . Les gaz frais sont donc du "côté " et les gaz brûlés du "côté ". La fraction massique vaut 1 quand la réaction n'a pas eu lieu et 0 quand tout le fuel a été brûlé. On suppose que et . La flamme se déplace à la vitesse inconnue strictement positive, des gaz brûlés vers les gaz frais (Figure 2)

II. 1 De la même manière que dans le problème I, on peut se ramener à une seule équation sur la température; c'est l'objet de cette question.
II.1.a Dans un premier temps, on se donne une vitesse de flamme positive. Ecrire le système de deux équations différentielles ordinaires satisfaites par les deux fonctions et .
II.1.b On note . Montrer que cette quantité est constante, c'est-à-dire qu'elle ne dépend ni du temps ni de l'espace. Donner sa valeur. (Indication : résoudre l'équation différentielle satisfaite par en utilisant les conditions à l'infini pour déterminer les constantes).
II.1.c Montrer alors que le problème se réduit à une équation sur du type

et donner une expression explicite pour la non-linéarité .
II. 2 On suppose que l'équation sur ci-dessus admet une solution régulière . Montrer que cette solution est monotone en (indication : supposer l'existence d'un minimum et montrer qu'il y a une contradiction en étudiant la convexité de la solution).
II. 3 Nous allons maintenant faire quelques hypothèses sur la non-linéarité afin de montrer l'existence et surtout l'unicité du profil de température et de la vitesse . On suppose donc qu'il existe tel que, pour pour et , , avec et fini. Une telle fonction est représentée sur la Figure 3.

On va réduire l'étude d'une équation du second ordre à celle d'un système de deux équations à deux inconnues du premier ordre et traiter l'existence dans le plan de phase . Par soucis de simplicité, on suppose que .
II.3.a On note , la dérivée de . Ecrire le système de deux équations différentielles du premier ordre sur le couple de nouvelles variables . Le plan est appelé le plan de phase.
II.3.b Montrer que l'on peut reposer ce système dans le plan ( ) avec et sous la forme

avec et les conditions aux limites , dont la solution (courbe intégrale) permet de joindre les deux états stationnaires et comme on le montre sur la Figure 4.
II.3.c On cherche dans un premier temps , la pente de la solution près du point . On peut effectuer un développement de Taylor et supposer et , avec . En utilisant le fait que vérifie l'équation (II.1), déterminer en fonction de et de .
II.3.d On se place sur l'intervalle . On peut alors considérer qu'à une vitesse donnée, on peut trouver un profil qui vérifie (II.1). Montrer que la solution sur cet intervalle est une fonction décroissante de , donner l'expression de cette solution pour ainsi que la valeur de en fonction de .
II.3. e Résoudre l'équation dans l'intervalle . Quel est le sens de variation de cette solution sur en fonction de ? En déduire l'existence d'une unique vitesse telle qu'il existe une courbe intégrale continue qui joigne et dans le plan de phase . On note cette solution.