Méthode de la zone fondue

La thermodynamique, plus particulièrement les diagrammes de phases entre le composé à purifier et ses impuretés, aide à comprendre le fonctionnement du procédé et à prévoir son efficacité. Plusieurs comportements d'impuretés peuvent être distingués.

Une impureté qui abaisse le point de fusion implique, sur le diagramme de phases entre le composé à purifier et l'impureté, des solidus et liquidus décroissants avec augmentation croissante de l'impureté dans le système. Ainsi, à concentration donnée d'impureté dans la phase liquide, la phase solide qui recristallise présente une concentration plus faible en impureté que la zone fondue.

Ces impuretés sont donc déplacées dans le même sens que la zone fondue.

Une impureté qui augmente le point de fusion implique, sur le diagramme de phases entre le composé à purifier et l'impureté, des solidus et liquidus croissants avec augmentation croissante de l'impureté dans le système. Ainsi, à concentration donnée d'impureté dans la phase liquide, la phase solide qui recristallise présente une concentration plus importante en impureté que la zone fondue.

Ces impuretés sont donc déplacées dans le sens inverse du sens de déplacement de la zone fondue.

Il peut arriver que certaines impuretés n'abaissent, ni n'augmentent le point de fusion du composé à purifier. Ce cas de figure est relativement restreint : il peut par exemple concerner un système contenant deux énantiomères en solution solide totale et idéale, ou un solide à purifier dont la composition est celle d'un composé à fusion congruente.

Dans ce cas, à chaque déplacement élémentaire de zone fondue, la fraction recristallisée présente la même concentration en impureté que la zone. La fusion de zone devient ainsi une technique non-pertinente pour purifier le solide considéré.

Bien qu'il soit rare de trouver de tels cas de figure, il est possible de trouver des situations qui s'en rapprochent. Si une impureté abaisse ou augmente le point de fusion du composé à purifier, mais que les courbes de liquidus et de solidus sont très proches l'une de l'autre, la fusion de zone n'est pas la technique la plus adaptée à la purification du composé.

La vitesse de déplacement de la zone fondue présente une importance notable pour l'efficacité du procédé. L'efficacité maximale théorique de la technique est obtenue pour une vitesse de déplacement qui tend vers zéro. À l'inverse, plus la vitesse est rapide, moins l'efficacité est bonne (jusqu'à devenir nulle).

Coupe transversale schématique d'un lingot

On considère un lingot cylindrique droit de solide impur présentant une teneur massique ${\displaystyle {C}_{0}}$ en impureté. Sa longueur vaut ${\displaystyle l}$ et sa section droite, supposée constante, ${\displaystyle {S}_{0}}$. Le lingot est raffiné avec une zone cylindrique droite de longueur ${\displaystyle e}$ maintenue constante, à vitesse de translation constante (voir schéma ci-dessus). On désigne par :

• ${\displaystyle x}$ : la position du front de solidification dans le lingot (${\displaystyle {0}\leqslant {x}\leqslant {l}}$).
• ${\displaystyle {C}_{{\text{S}},i}\left(x\right)}$ : la teneur massique fractionnaire en impureté dans la tranche de solide redéposée à la position ${\displaystyle x}$ après le passage de zone fondue de rang ${\displaystyle i}$ (${\displaystyle 0\leqslant {C}_{{\text{S}},i}\left(x\right)\leqslant 1{\text{$ ; }}i\geqslant 1} ).
• ${\displaystyle {C}_{{\text{L}},i}\left(x\right)}$ : la teneur massique fractionnaire en impureté dans la zone fondue lorsque le front de solidification se trouve à la position ${\displaystyle x}$ au passage de rang ${\displaystyle i}$ (${\displaystyle 0\leqslant {C}_{{\text{L}},i}\left(x\right)\leqslant 1}$).
• ${\displaystyle {\rho }_{\text{S}}}$ et ${\displaystyle {\rho }_{\text{L}}}$ : les masses volumiques des phases solide et liquide (fondue).

On cherche à exprimer ${\displaystyle {C}_{{\text{S}},i}\left(x\right)}$, la teneur du lingot en impureté quels que soient la position dans le lingot et le rang du passage de zone fondue. Pour cela, on introduit une grandeur utile : le coefficient de ségrégation de l’impureté, qu’on peut supposer dépendre, dans un cas général, de la composition de la zone fondue (donc, de la position de la zone et du rang du passage), et qu’on note par conséquent ${\displaystyle {k}_{i}\left(x\right)}$. Ce coefficient traduit le rapport entre la composition de la tranche de solide déposée au front de solidification de la zone à un instant donné, et celle de la zone : ${\displaystyle {k}_{i}\left(x\right)={\frac {{C}_{{\text{S}},i}\left(x\right)}{{C}_{{\text{L}},i}\left(x\right)}}}$.

Pour résoudre ce problème, il est nécessaire de réaliser deux bilans de matière distincts : en effet, tant que le front de fusion n’a pas atteint l’extrémité finale du lingot (${\displaystyle 0\leqslant x\leqslant l-e}$), la zone subit une entrée de matière permanente puisqu’elle incorpore une tranche de solide au front de fusion à chacun de ses déplacements élémentaires. Ceci n’est en revanche plus vrai lorsque le front de fusion a atteint la longueur du lingot (${\displaystyle l-e<x<l}$).

À chaque déplacement élémentaire ${\displaystyle {\mathit {dx}}}$ de la zone, la masse d’impureté incorporée à la zone au front de fusion est ${\displaystyle {S}_{0}\cdot {\rho }_{\text{S}}\cdot {C}_{{\text{S}},i-1}\left(x+e\right)\cdot {\mathit {dx}}}$. La masse d’impureté rejetée dans la tranche de solide redéposée au front de solidification est ${\displaystyle {S}_{0}\cdot {\rho }_{\text{S}}\cdot {C}_{{\text{S}},i}\left(x\right)\cdot {\mathit {dx}}}$. La masse d’impureté accumulée dans la zone s’écrit quant à elle ${\displaystyle {S}_{0}\cdot {\rho }_{\text{L}}\cdot {\mathit {dC}}_{{\text{L}},i}\left(x\right)\cdot e}$. En supposant que les densités des phases liquide et solide sont identiques, le principe de conservation de la matière implique l’égalité suivante :

${\displaystyle {\mathit {dC}}_{{\text{L}},i}\left(x\right)\cdot e=\left({C}_{{\text{S}},i-1}\left(x+e\right)-{C}_{{\text{S}},i}\left(x\right)\right)\cdot {\mathit {dx}}}$

On peut alors utiliser la définition du coefficient de ségrégation, ${\displaystyle {k}_{i}\left(x\right)={\frac {{C}_{{\text{S}},i}\left(x\right)}{{C}_{{\text{L}},i}\left(x\right)}}}$, pour réexprimer ${\displaystyle {\mathit {dC}}_{{\text{L}},i}\left(x\right)}$ en fonction de ${\displaystyle {C}_{{\text{S}},i}\left(x\right)}$. On obtient alors ${\displaystyle {\mathit {dC}}_{{\text{L}},i}\left(x\right)={\frac {{\mathit {dC}}_{{\text{S}},i}\left(x\right)}{{k}_{i}\left(x\right)}}}$, ce qui aboutit à l’équation différentielle suivante :

${\displaystyle {\frac {{\mathit {dC}}_{{\text{S}},i}\left(x\right)}{\mathit {dx}}}={\frac {{k}_{i}\left(x\right)}{e}}\cdot \left({C}_{{\text{S}},i-1}\left(x+e\right)-{C}_{{\text{S}},i}\left(x\right)\right)}$

Pour chaque passage de zone de rang ${\displaystyle i}$, la condition initiale est donnée par :

${\displaystyle {C}_{{\text{S}},i}\left(x=0\right)={\frac {{k}_{i}\left(x=0\right)}{e}}\cdot {\underset {0}{\overset {e}{\int }}}{C}_{{\text{S}},i-1}\left(x\right)\cdot {\mathit {dx}}}$

Si le coefficient de ségrégation et la longueur de zone sont supposés constants tout le long du procédé (on note ainsi ${\displaystyle {k}_{i}\left(x\right)=k}$), il existe alors des expressions analytiques de la condition initiale et du profil de distribution de l’impureté dans la région ${\displaystyle 0\leqslant x\leqslant l-e}$ pour le premier passage de zone fondue :

${\displaystyle {C}_{{\text{S}},1}\left(x=0\right)={\frac {k}{e}}\cdot {C}_{0}}$

${\displaystyle {C}_{{\text{S}},1}\left(0\leqslant x\leqslant l-e\right)={C}_{0}\cdot \left(1-\left(1-k\right)\cdot \exp \left(-{\frac {k\cdot x}{e}}\right)\right)}$

Lorsque le front de fusion atteint le haut du lingot, il n’y a plus d’entrée de matière dans la zone et la longueur de cette dernière diminue au fur et à mesure qu’elle avance (${\displaystyle e=l-x}$). À chaque déplacement ${\displaystyle {\mathit {dx}}}$ de la zone, une tranche de liquide contenant ${\displaystyle {S}_{0}\cdot {\rho }_{\text{L}}\cdot {C}_{{\text{L}},i}\left(x\right)\cdot {\mathit {dx}}}$ d’impureté en masse se scinde en deux parties : l’une s’incorpore à la tranche de solide déposée au front de solidification et contient ${\displaystyle {S}_{0}\cdot {\rho }_{\text{S}}\cdot {C}_{{\text{S}},i}\left(x\right)\cdot {\mathit {dx}}}$ en masse d’impureté, et l’autre contribue à l’accumulation d’impureté dans la zone à hauteur d’une masse de ${\displaystyle {S}_{0}\cdot {\rho }_{\text{L}}\cdot {\mathit {dC}}_{{\text{L}},i}\left(x\right)\cdot \left(l-x\right)}$. En supposant que ${\displaystyle {\rho }_{\text{L}}={\rho }_{\text{S}}}$, la conservation de la matière donne donc l’égalité qui suit :

${\displaystyle \left({C}_{{\text{L}},i}\left(x\right)-{C}_{{\text{S}},i}\left(x\right)\right)\cdot {\mathit {dx}}=\left(l-x\right)\cdot {\mathit {dC}}_{{\text{L}},i}\left(x\right)}$

En utilisant la définition du coefficient de ségrégation, cette équation devient :

${\displaystyle {\frac {{\mathit {dC}}_{{\text{S}},i}\left(x\right)}{\mathit {dx}}}={\frac {{C}_{{\text{S}},i}\left(x\right)\cdot \left(1-{k}_{i}\left(x\right)\right)}{l-x}}}$

La condition initiale de cette équation est la valeur de la teneur du lingot en impureté en ${\displaystyle x=l-e}$, qu’on calcule à l’aide de l’équation différentielle établie dans la partie précédente. Si le coefficient de ségrégation est supposé constant, l’équation présente une solution analytique valable quel que soit le rang du passage de zone :

${\displaystyle {C}_{{\text{S}},i}\left(x\right)={C}_{{\text{S}},i}\left(l-e\right)\cdot {\left({\frac {l-x}{e}}\right)}^{k-1}}$

L’une des méthodes consistant à résoudre numériquement les équations différentielles permettant de déterminer le profil de distribution d’une impureté après raffinage par zone fondue est de discrétiser le lingot en intervalles ${\displaystyle \mathrm {\Delta } x}$ les plus petits possibles.

Pour chaque passage de rang ${\displaystyle i}$, la teneur en impureté du lingot en ${\displaystyle x=0}$ se calcule grâce à l’équation suivante :

${\displaystyle {C}_{{\text{S}},i}\left(x=0\right)={\frac {{k}_{i}\left(x\right)\cdot \mathrm {\Delta } x}{e}}\cdot \sum _{0}^{e-\mathrm {\Delta } x}{C}_{{\text{S}},i-1}\left(x\right)}$

Pour le premier passage de zone fondue, ${\displaystyle {C}_{{\text{S}},i-1}\left(x\right)={C}_{0}}$.

Pour chaque passage, la composition de la tranche déposée en ${\displaystyle x+\mathrm {\Delta } x}$ vaut :

${\displaystyle {C}_{{\text{S}},i}\left(x+\mathrm {\Delta } x\right)={C}_{{\text{S}},i}\left(x\right)+{\frac {{k}_{i}\left(x\right)\cdot \mathrm {\Delta } x}{e}}\cdot \left({C}_{{\text{S}},i-1}\left(x+e\right)-{C}_{{\text{S}},i}\left(x\right)\right){\text{ pour }}0<x\leqslant l-e-\mathrm {\Delta } x}$

${\displaystyle {C}_{{\text{S}},i}\left(x+\mathrm {\Delta } x\right)={C}_{{\text{S}},i}\left(x\right)+{\frac {\mathrm {\Delta } x}{l-x}}\cdot {C}_{{\text{S}},i}\left(x\right)\cdot \left(1-{k}_{i}\left(x\right)\right){\text{ pour }}l-e\leqslant x\leqslant l-2\times \mathrm {\Delta } x}$

Ces calculs peuvent être mis en œuvre grâce à un simple tableur, ou un logiciel de calcul type Scilab/Matlab/GNU Octave.