These Anne Lenoble
These Anne Lenoble
THESE
présentée pour obtenir le grade de
Discipline : Mécanique-Energétique
par
Anne LENOBLE
THESE
présentée pour obtenir le grade de
Discipline : Mécanique-Energétique
par
Anne LENOBLE
Remerciements
Je tiens tout d’abord à remercier les membres du jury pour avoir accepter de juger ce travail, et tout
spécialement les rapporteurs.
Je remercie également le groupe Saint-Gobain Recherche, à travers Paul-Henri Guering, et l’ADEME,
pour avoir financer ce travail de thèse.
J’exprime toute ma gratitude à Fabrice Onofri, pour avoir diriger ce travail de thèse.
J’adresse également tout mes remerciements aux personnes qui m’ont aidée au cours de ce travail, par
leur travail ou leurs conseils, et notamment :
– Hervé Bultynck, pour mon encadrement à Saint-Gobain Recherche durant les deux dernières années
de thèse.
– Nicolas Marsault, Arnaud d’Aubigny et Pierre Deleplace, de Saint-Gobain Vetrotex-Chambéry, pour
leurs conseils avisés et pour la mise en route de la filière.
– Ouamar Rahli, et mes parents, pour la relecture de ce manuscrit.
– Fabrice Rigollet, pour les tests de caractérisation des fibres par Microscopie à Force Atomique.
– Scott A. Schaub (Mission Research Corporation, USA) et Amir Naqwui (Univ. Minesota, USA),
pour les validations dans le cas d’une fibre homogène.
– Stefan Radev (Académie des Sciences de Bulgarie) ainsi que Richard Saurel et Eric Daniel pour
leurs conseils sur l’étude de l’écoulement du verre.
– Roland Faure, Jean-Claude Morellini, Jean Lequéré et Paul Cervetti, ingénieurs et techniciens du
laboratoire, pour leur contribution à la réalisation des parties les plus techniques de ce travail.
Introduction générale 11
3
4 Table des matières
Bibliographie 160
6 Table des matières
Nomenclature
7
8 Table des matières
Lin
σD résolution moyenne de FIBS à 1 − σD , PDRS linéaire
n
σD résolution moyenne de FIBS à 1 − σD , PDRS polynomiale (n > 1)
FFT
σD résolution moyenne de FIBS à 1 − σD , PDRS filtrée par FFT
Abréviations :
CSDW Fonction interspectre à fenêtre glissante
F IBS Interféromètre : FIBER SIZER
HIREDI Diffractomètre : HIGH RESOLUTION DIFFRACTOMETER
P DRS Relation phase-Diamètre (Phase-Diameter RelationShip)
M EB Microscopie Electronique à Balayage
T LM Théorie de Lorenz-Mie
Contexte de l’étude
Les «fibres» sont des objets qui tendent vers l’unidimensionnalité et qui relient la matière à différentes
échelles, notamment le monde microscopique (par leur section) au monde macroscopique (par leur lon-
gueur).
Cette propriété d’échelle confère à la matière, ainsi mise en forme, des propriétés spécifiques. La
Nature utilise largement les «fibres», les matériaux «fibreux», et ceci, pour essentiellement trois types
de fonctions : pour l’isolation thermique (fourrures, plumages...), pour le transport de la matière et
de l’information (radicelles, systèmes vasculaires, nerfs...), et pour leurs propriétés mécaniques (fibres
musculaires, fibres ligneuses,...)
Par mimétisme, très certainement, l’industrie utilise également depuis longtemps les fibres et ceci, à
quelques exceptions près, pour les mêmes fonctions que la Nature : pour l’isolation thermique (tissus,
fourrure synthétique, laine de roche,...), pour le transport de la matière (fluides moléculaires, électrons...)
et de l’immatériel (photons), et pour améliorer les propriétés mécaniques des matériaux (fibres de verre,
tresses métalliques, nano tubes de carbone...).
L’utilisation toujours croissante des fibres dans l’industrie concourt au développement de nombreux
travaux de recherche sur leur procédé de fabrication et sur leur utilisation. Ainsi durant ces dernières
années, au sein du laboratoire, l’IUSTI-UMR CNRS n◦ 6595, différentes études mettant en jeux des
milieux fibreux ont été conduites : compaction d’enchevêtrements de fibres, sédimentation de suspensions
de fibres, hydrodynamique et transfert thermique dans les milieux fibreux ou dans les capillaires, les
micro-canaux...
Ce travail de thèse porte, quant à lui, sur les fibres de verre utilisées pour les applications de renfort,
et plus spécifiquement, sur le développement des outils nécessaires à l’étude du domaine de stabilité d’un
procédé industriel de fabrication de ces fibres. Ce travail a été effectué en collaboration avec le Groupe
Saint-Gobain (Saint-Gobain Recherche & Saint-Gobain Vetrotex) et a été soutenu financièrement par
l’ADEME (Agence de l’Environnement et de la Maı̂trise de l’Energie).
12 Introduction générale
Fig. 2 – Analyse par Microscopie électronique des éléments présents dans le verre-E : Silicium, Oxygène,
Aluminium, Calcium, Magnésium, Bore, Sodium...
La problématique de l’étude
Bien que le procédé industriel de fabrication de ces fibres de renforcement soit opérationnel, le groupe
Saint-Gobain souhaite l’optimiser en vue notamment de : i) augmenter la tirée des filières (i.e. une vitesse
de bobinage maximale) pour accroitre la productivité du procédé ; ii) obtenir un diamètre de fibre le plus
constant possible (propriétés spécifiques du produit constantes) ; iii) réduire l’occurrence des fibres creuses
qui posent problèmes pour les applications électroniques ; iv) réduire l’occurrence de la ”casse” des fibres
en cours de fibrage. En effet, lors d’une casse, pendant les quelques minutes nécessaires aux opérateurs
pour redémarrer le procédé, le verre continue à s’écouler du four. Du fait de l’ensimage, ce verre ne
peut être refondu dans le bain de la composition initiale. Il est donc placé en décharge. La casse induit
donc des pertes à plusieurs niveaux : matière première, énergie, transport et stockage des déchets. En
France, la quantité de déchets verriers ainsi produits s’élève à plusieurs milliers de tonnes par an pour
une production annuelle de plus de 590 000 tonnes de fibres de verre3 .
L’optimisation du procédé de fibrage nécessite d’étudier son domaine de stabilité. Cependant, les
mécanismes physiques qui peuvent conduire à sa déstabilisation sont a priori nombreux et variés :
- instabilités inertielles et capillaires qui se développent le long du jet visqueux,
- interaction rayonnement/hydrodynamique, renforcée par la nature semi transparente du verre et les
hautes températures nécessaires au fibrage du verre ; perturbations radiatives et convectives externes
(casse d’autres fibres, turbulences induites par le bobinage...),
- inhomogénéités dans la composition de verre (température/viscosité, bulles de gaz...),
1 sous
forme d’oxydes : Silice, Potasse, Alumine...
2 La
figure 7.1 illustre le principe de fonctionnement d’une filière de laboratoire, mono-téton, et l’instrumentation mise
en place pour son étude.
3 en 1998
Introduction générale 13
Organisation du mémoire
Conformément à la problématique de cette étude, ce manuscrit se compose de deux parties. La
première partie présente les modèles de diffusion de la lumière et les travaux expérimentaux qui ont
été nécessaires au développement d’un granulomètre laser capable de suivre en temps réel le diamètre des
fibres en cours de production. La seconde partie présente nos travaux de modélisation de l’hydrodyna-
mique du procédé de fibrage ainsi que les résultats expérimentaux préliminaires obtenus sur la stabilité
d’une filière de laboratoire.
Au terme de ce manuscrit, nous dresserons un bilan des principaux résultats obtenus et prolongerons
la réflexion en ouvrant des perspectives possibles à ce travail.
Première partie
15
Chapitre 1
Dans ce chapitre nous présentons le formalisme et les algorithmes numériques qui ont été spécifiquement
développés au cours de ce travail, pour simuler les propriétés de diffusion de la lumière d’un cylindre mul-
ticouche.
Les résultats théoriques et numériques de ce chapitre seront utilisés pour simuler la réponse de l’in-
terféromètre (Chapitre 2) et du diffractomètre (Chapitre 4) à la mesure des fibres de renfort.
Z
Onde Y
Plane D/2
eϕ
er
ez
i
Eϕ
r
ϕ
k ext m1
Ez
i Z X
m2
...
mj
...
mL
mext
ω rL ... r j ... r2 r1
λ 0 /m ext
1.1 Théorie
On considère un cylindre infini de section circulaire, composé de L-couches concentriques avec
1/2
j = 1, 2, . . . , L − 1, L. Chaque couche est caractérisée par son indice de réfraction complexe mj = (εj ) ,
par sa constante de perméabilité magnétique µj et par son rayon rj , voir la figure 1.1.
Ce cylindre, placé dans un milieu d’indice réel mext , est éclairé par une onde plane harmonique de
vecteur d’onde k = kex = (2πmext /λ) ex . La géométrie du modèle et les conventions sont les mêmes que
celles introduites par Bohren et Huffman pour un cylindre homogène [10].
Le problème se ramène à la résolution de l’équation d’onde (∇2 Ψ + k 2 Ψ = 0) pour les conditions
limites du problème : décroissance des champs diffusés en 1/r (champ lointain) ; les champs internes,
incidents et diffusés, sont de dimension finie au centre de la particule ; la continuité des composantes
tangentielles des champs (électriques et magnétiques) à la traversée de la surface externe de la particule
et de chaque couche interne j.
17
18 Chapitre 1. Diffusion électromagnétique par un cylindre multicouche
Les solutions de l’ équation 1.3 sont les composantes des harmoniques cylindriques [10] :
(p) (p)
Mn = −kZn (ρ) einϕ eϕ
(p) (p) (1.4)
Nn = kZn (ρ) einϕ ez
où les Jn (ρ) sont les fonctions de Bessel[1] d’ordre entier (n ∈ Z) et du premier ordre ; les Yn (ρ) sont les
fonctions de Bessel du second ordre (ou encore fonction de Neumann) ; les Hn (ρ) sont les fonctions de
Hankel correspondantes.
∞
(1)
EiT M = En Nn
n=−∞
∞ (1.6)
−ikext (1)
HiT M = ωµext En Mn
n=−∞
avec pour le vecteur d’onde de l’onde incidente : kext = 2πmext /λ0 où λ0 est la longueur d’onde de
n
l’onde incidente «dans l’air». En = E0 (−i) /kext où E0 est l’amplitude du champ incident [10]. Les
composantes des harmoniques se réduisent alors à :
(1)
Mn (mext xL ) = −kext Jn (mext xL ) eϕ
(1) (1.7)
Nn (mext xL ) = kext Jn (mext xL ) ez
Ce qui donne pour les composantes des champs électriques et magnétiques à la surface de la particule :
1 les dérivées correspondantes sont notées : Z (p)
2Y
n (ρ → 0) → −∞
1.1. Théorie 19
∞
Ezi (mext xL ) = kext En Jn (mext xL )
n=−∞
Hzi (mext xL ) =0
(1.8)
Eϕi (mext xL ) =0
i(kext )2
∞
Hϕi (mext xL ) = ωµext En Jn (mext xL )
n=−∞
∞
(3)
EsT M = − En bnI Nn
n=−∞
∞ (1.9)
ikext (3)
HsT M = ωµext En bnI Mn
n=−∞
En introduisant le paramètre de taille du cylindre : xL = 2πrL /λ0 , les composantes des harmoniques
s’écrivent sur la surface externe :
(3) (1)
Mn (mext xL ) = −kext Hn (mext xL ) eϕ
(3) (1) (1.10)
Nn (mext xL ) = kext Hn (mext xL ) ez
Ce qui donne pour les composantes tangentielles des champs diffusés :
∞
(1)
Ezs = −kext En bnI Hn (mext xL )
n=−∞
Hzs = 0
(1.11)
Eϕs = 0
−i(kext )2
∞
(1)
Hϕs = ωµext En bnI Hn (mext xL )
n=−∞
Ces dernières expressions, pour le champ électrique et le champ magnétique, sont strictement iden-
tiques à celles obtenues pour un cylindre homogène [10]. Seule l’expression du coefficient de diffusion
externe bnI diffère (voir plus loin).
(j) (j)
Calcul des champs internes : ET M , HT M
Nous utilisons ici la méthode de résolution par récurrence introduite par Wu et al. [109] pour le cas
de la sphère stratifiée éclairée par une onde plane et qui a été étendue par Onofri et al. [71] pour une
onde incidente de forme arbitraire [37].
Pour une onde incidente TM, les expressions du champ électrique et du champ magnétique internes
à la couche j sont de la forme :
∞
(j) (j) (1) (j) (2)
ET M = En cn Nn − en Nn
n=−∞
∞ (1.12)
(j) (j) (1) (j) (2)
HT M ikext
= ωµext
En c n M n − e n M n
n=−∞
(j) (j)
Les coefficients de diffusion internes cn , en sont a priori non nuls. Les composantes des harmoniques
s’écrivent avec kj = 2πmj /λ0 et ρ = mj xj :
(p) (p)
Mn (ρ) = −kj Zn (ρ) eϕ
(p) (p) (1.13)
Nn (ρ) = kj Zn (ρ) ez
En fait, dans le cas du cylindre homogène, pour la même raison que pour le champ incident, la fonction
de Bessel du second ordre, Yn , ne peut être solution de l’Eq. (1.3). Dans le cas du cylindre multicouche,
3I (r) ∝ Hn (r) H̄n (r) → α/r, r→∞
20 Chapitre 1. Diffusion électromagnétique par un cylindre multicouche
la condition de dimension finie n’est critique que pour la première couche (qui est en fait un cylindre
homogène), j = 1. Pour les autres couches, j = 2 . . . L, cette restriction n’a pas lieu d’être. C’est pourquoi
on doit exprimer les champs internes, pour chaque couche j, comme une combinaison linéaire des deux
(1)
fonctions de Bessel, Jn , Yn , avec ρ = mj xj et pour le cas particulier de la première couche : en = 0. On
obtient ainsi :
∞
(j) (j) (j)
En,z = kj En cn Jn (mj xj ) − en Yn (mj xj )
n=−∞
(j)
Hn,z = 0
(j) (1.14)
En,ϕ = 0
∞
(j) i(kj )2 (j) (j)
Hn,ϕ = ωµj En cn Jn (mj xj ) − en Yn (mj xj )
n=−∞
Avec pour le cas particulier de la couche externe du cylindre, j=L, ρ = mL xL
∞
(L) (L) (L)
En,z = kL En cn Jn (mL xL ) − en Yn (mL xL )
n=−∞
(L)
Hn,z = 0
(L) (1.15)
En,ϕ = 0
∞
(L) i(kL )2 (L) (L)
Hn,ϕ = ωµL En cn Jn (mL xL ) − en Yn (mL xL )
n=−∞
(j) (j−1)
en en
mj c(j)
n Jn (mj xj−1 ) − Y (mj xj−1 ) =
(j) n
mj−1 c(j−1)
n Jn (mj−1 xj−1 ) − (j−1)
Yn (mj−1 xj−1 )
cn cn
(1.17)
(j) (j−1)
(mj )2 (j) en (mj−1 )2 (j−1) en
cn Jn (mj xj−1 ) − (j) Yn (mj xj−1 ) = cn Jn (mj−1 xj−1 ) − (j−1) Yn (mj−1 xj−1 )
µj cn µj−1 cn
(1.18)
Plutôt que de calculer les 2n (L − 1) coefficients internes on peut se limiter au calcul des n (L − 1)
(1)
quotients suivants ( avec Bn = 0) :
(j)
en
(j)
Bn(j) = (1.19)
cn
Ce qui permet, avec le calcul du rapport Eq. (1.18)/Eq.(1.17), d’éliminer une série d’inconnues :
(j) (j−1)
Jn (mj xj−1 ) − Bn Yn (mj xj−1 ) µj mj−1 Jn (mj−1 xj−1 ) − Bn Yn (mj−1 xj−1 )
(j)
= (j−1)
(1.20)
Jn (mj xj−1 ) − Bn Yn (mj xj−1 ) µj−1 mj Jn (mj−1 xj−1 ) − Bn Yn (mj−1 xj−1 )
1.1. Théorie 21
(j−1)
(j−1) Jn (mj−1 xj−1 ) − Bn Yn (mj−1 xj−1 )
Hb,n = (j−1)
(1.21)
Jn (mj−1 xj−1 ) − Bn Yn (mj−1 xj−1 )
µj mj−1 (j−1)
Jn (mj xj−1 ) − Bn(j) Yn (mj xj−1 ) = Hb,n Jn (mj xj−1 ) − Bn(j) Yn (mj xj−1 ) (1.22)
µj−1 mj
(j)
De cette dernière expression on peut tirer la forme récurrente du quotient Bn :
(j−1)
mj−1 µj Jn (mj xj−1 ) Hb,n − mj µj−1 Jn (mj xj−1 )
Bn(j) = (j−1)
(1.23)
mj−1 µj Yn (mj xj−1 ) Hb,n − mj µj−1 Yn (mj xj−1 )
En remplaçant dans l’équation (1.24) les champs par leur expression sur la surface externe, et en
simplifiant, on obtient le système suivant :
i (kext )
2
i (k )2
(L)
En Jn (mext xL ) − bnI Hn(1) (mext xL ) = c(L)
L
n En Jn (mL xL ) − Bn Yn (mL xL ) (1.25)
ωµext ωµL
kext En Jn (mext xL ) − bnI Hn(1) (mext xL ) = kL c(L) (L)
n En Jn (mL xL ) − Bn Jn (mL xL ) (1.26)
Le calcul du rapport des Eqs. (1.25)/ (1.26) permet d’éliminer les inconnues suivantes : le coefficient
(L)
de diffusion interne cn et la constante En , d’où :
(1) (L)
kext Jn (mext xL ) − bnI Hn (mext xL ) kL Jn (mL xL ) − Bn Yn (mL xL ) kL L
= = H (1.27)
µext Jn (mext xL ) − bnI Hn(1) (mext xL ) µL Jn (mL xL ) − Bn(L) Jn (mL xL ) µL b,n
L
Dans le second membre de l’expression précédente on reconnaı̂t l’expression du rapport Hb,n introduite
L
par l’Eq. (1.21). En remplaçant Hb,n par son expression et en simplifiant on obtient :
kext µL Jn (mext xL ) − kext µL bnI Hn(1) (mext xL ) = kL µext Hb,n
L
Jn (mext xL ) − kL µext Hb,n
L
bnI Hn(1) (mext xL )
(1.28)
De cette dernière équation on peut tirer l’expression du coefficient bnI :
L
mL µext Hb,n (mext xL ) Jn (mext xL ) − mext µL Jn (mext xL )
bnI = (1) (1)
(1.29)
ext xL ) Hn (mext xL ) − mext µL Hn
L (m
mL µext Hb,n (mext xL )
22 Chapitre 1. Diffusion électromagnétique par un cylindre multicouche
∞
(1)
EiT E = −i En Mn
n=−∞
∞ (1.30)
kext (1)
HiT E = ωµext En Nn
n=−∞
soit sur la surface externe de la particule et en remplaçant par les composantes des harmoniques cylin-
driques :
Ezi (mext xL ) = 0
(kext )2
∞
Hzi (mext xL ) = ωµext En Jn (mext xL )
n=−∞
∞
(1.31)
Eϕi (mext xL ) = ikext En Jn (mext xL )
n=−∞
Hϕi (mext xL ) = 0
- Champs diffusés :
∞
(3)
EsT E = i En anII Mn
n=−∞
∞ (1.32)
(3)
HsT E = − ωµ
kext
ext
En anII Nn
n=−∞
soit à la surface de la particule et en remplaçant par les composantes des harmoniques cylindriques :
Ezs = 0
2
∞
(1)
Hzs = − (kext )
ωµext En anII Hn (mext xL )
n=−∞
∞ (1.33)
(1)
Eϕs = −ikext En anII Hn (mext xL )
n=−∞
Hϕs = 0
- Champs internes :
∞
(j) (j) (1) (j) (2)
ET M = En dn Mn − fn Mn
n=−∞
∞ (1.34)
(j) (j) (1) (j) (2)
HT M ikext
= ωµext
E n d n N n − f n N n
n=−∞
(1)
la condition de «dimension finie» des champs impose pour la première couche : fn = 0.
Soit pour les champs internes à la couche j et en remplaçant par les composantes des harmoniques
cylindriques :
(j)
En,z = 0
∞
(j) kj (j) (j)
Hn,z = ωµj En dn Jn (mj xj ) − fn Yn (mj xj )
n=−∞
∞ (1.35)
(j) (j) (j)
En,ϕ = kj En dn Jn (mj xj ) − fn Yn (mj xj )
n=−∞
(j)
Hn,ϕ = 0
1.2. Formulation pour le calcul numérique 23
Au final la condition de continuité, équation (1.24), permet d’obtenir la forme du coefficient de diffusion
externe pour une onde Transversale Electrique (TE), anII :
L
mext µL Ha,n (mext xL ) Jn (mext xL ) − mL µext Jn (mext xL )
anII = (1) (1)
(1.40)
ext xL ) Hn (mext xL ) − mL µext Hn
L (m
mext µL Ha,n (mext xL )
Jn (ρ) Jn (ρ)
Dn(0) (ρ) = (1)
, Dn(x) (ρ) = (1.42)
Hn (ρ) Yn (ρ)
Il convient ensuite de faire apparaı̂tre ces dérivées et ces rapports dans les équations précédentes. Pour
les champs TM, par exemple, on peut reformuler l’équation (1.23) comme suit :
(j−1) J (m x )
Jn (mj xj−1 ) mj−1 µj Hb,n − mj µj−1 Jn (mj xj−1 )
n j j−1
Bn(j) = (1.43)
Yn (mj xj−1 ) mj−1 µj H (j−1) − mj µj−1 Yn (mj xj−1 )
b,n Yn (mj xj−1 )
24 Chapitre 1. Diffusion électromagnétique par un cylindre multicouche
(j)
En procédant de même pour An on obtient :
(j−1) (1)
mj µj−1 Ha,n − mj−1 µj Dn (mj xj−1 )
A(j) (x)
n = Dn (mj xj−1 ) (j−1) (2)
(1.45)
mj µj−1 Ha,n − mj−1 µj Dn (mj xj−1 )
(j)
Pour le rapport Hb,n , on cherche à faire apparaı̂tre les dérivées définies par l’équation(1.41) et les
rapports explicités par l’équation(1.42) en partant de l’expression :
(j)
(j) Jn (mj xj ) − Bn Yn (mj xj )
Hb,n = (j)
(1.46)
Jn (mj xj ) − Bn Yn (mj xj )
(j)
(j) Jn (mj xj ) Bn Yn (mj xj )
Hb,n = (j)
− (j)
(1.47)
Jn (mj xj ) − Bn Yn (mj xj ) Jn (mj xj ) − Bn Yn (mj xj )
Jn (mj xj ) (j) Yn (mj xj )
(j) Jn (mj xj ) Jn (mj xj ) Bn Yn (mj xj )
Hb,n = Jn (mj xj ) (j)
− Jn (mj xj ) (j)
(1.48)
Yn (mj xj ) − Bn Yn (mj xj ) − Bn
Yn (mj xj )
avec au final :
(L) (1)
mext µL Hb,n − mL µext Dn (mext xL )
anII = Dn(0) (mext xL ) (L) (3)
(1.52)
mext µL Hb,n − mL µext Dn (mext xL )
(L) (1)
mL µext Hb,n − mext µL Dn (mext xL )
bnI = Dn(0) (mext xL ) (L) (3)
(1.53)
mL µext Hb,n − mext µL Dn (mext xL )
Cette formulation est formellement identique à celle obtenue pour la sphère multicouche [109, 71] en
faisant les correspondances équations (1.88-1.89) avec µj = 1, j = 1 . . . L + 1 et mext = 1.
Une première version du code fortran Diff MCyl a mis en évidence une instabilité du calcul et no-
(x)
tamment au niveau du calcul du rapport Dn et ceci, même pour des particules de faible diamètre et de
quelques couches.
Dans le paragraphe suivant, nous reformulons les équations précédentes en transposant la modification
introduite par Wu et al. [110] pour le cas de la sphère multicouche.
1.2. Formulation pour le calcul numérique 25
Yn (mj xj )
βn(j) = Bn(j) (1.55)
Jn (mj xj )
et l’on obtient :
Jn (mj xj ) (j) Jn (mj xj ) Yn (mj xj )
(j) Jn (mj xj ) − βn Yn (mj xj ) Jn (mj xj )
Hb,n = (j)
(1.56)
1 − βn
avec au final :
(1) (j) (2)
(j) Dn (mj xj ) − βn Dn (mj xj )
Hb,n = (j)
(1.57)
1 − βn
(x)
dans cette dernière équation, le rapport Dn a été éliminé implicitement. On opère de même pour la
(j)
fonction βn :
(j−1) J (m x )
Yn (mj xj ) Jn (mj xj−1 ) mj−1 µj Hb,n − mj µj−1 Jn (mj xj−1 )
n j j−1
βn(j) = (1.58)
Jn (mj xj ) Yn (mj xj−1 ) mj−1 µj H (j−1) − mj µj−1 Yn (mj xj−1 )
b,n Yn (mj xj−1 )
(x) (j)
On introduit alors le rapport de deux fonctions Dn à travers la nouvelle fonction Rn :
(x)
Yn (mj xj ) Jn (mj xj−1 ) Dn (mj xj−1 )
Rn(j) = = (x)
(1.59)
Jn (mj xj ) Yn (mj xj−1 ) Dn (mj xj )
Au final on obtient :
(j−1) (1)
mj−1 µj Hb,n − mj µj−1 Dn (mj xj−1 )
βn(j) = Rn(j) (j−1) (2)
(1.60)
mj−1 µj Hb,n − mj µj−1 Dn (mj xj−1 )
En procédant de même pour une onde TE, on obtient :
(1) (j) (2)
(j) Dn (mj xj ) − αn Dn (mj xj )
Ha,n = (j)
(1.61)
1 − αn
(j−1) (1)
mj µj−1 Ha,n − mj−1 µj Dn (mj xj−1 )
α(j) (j)
n = Rn (j−1) (2)
(1.62)
mj µj−1 Ha,n − mj−1 µj Dn (mj xj−1 )
Ce profil est pris en compte dans le calcul en dimensionnant ce profil : xj = 2π (r1 /rL ) / (λ0 /mext ).
On pose par ailleurs mL+1 = mext , mext ∈ R.
De manière classique maintenant [10] (voir section 1.5.2), comme pour une particule homogène (sphère
ou cylindre), les séries infinies peuvent être tronquées a partir de :
1/3
nstop = xL + 4xL + 2 (1.64)
Cette dérivée est calculée par récurrence descendante a partir de la relation suivante :
n 1
Dn(1) (ρ) = − (1) (1.65)
ρ Dn+1 (ρ) + n+1
ρ
(1)
avec n = nmax , nmax − 1 . . . 0. Pour la valeur initiale on pose : Dn=max (ρ) = 0 + 0i, avec :
ρ = mj xj−1 , j = 2...L + 1
(1.67)
ρ = mj xj , j = 1...L
(2) (3)
Dérivées : Dn , Dn
(2) (3)
Les dérivées logarithmiques Dn et Dn sont calculées par récurrence ascendante à partir de la relation
(n = 0, 1, · · · nstop − 1) :
(p) 1 n+1
Dn+1 (ρ) = (p)
− (1.68)
n/ρ − Dn (ρ) ρ
(2)
Dans le cas de Dn , ces calculs sont effectués pour chaque couche (ρ = mj xj ) et avec comme valeurs
initiales :
(2) Y (ρ) −Y1 (ρ)
D0 (ρ) = 0 = (1.69)
Y0 (ρ) Y0 (ρ)
(3)
Dans le cas de Dn , ces calculs sont effectués uniquement pour l’interface externe (ρ = mext xL ) et avec
comme valeur initiale :
(1)
(3) H (mext xL ) J1 (mext xL ) + iY1 (mext xL )
Dn=0 (mext xL ) = 0 =− (1.70)
H0 (mext xL ) J0 (mext xL ) + iY0 (mext xL )
En pratique, pour initialiser les récurrences, les fonctions de Bessel J0 , J1 , Y0 et Y1 sont calculées à
l’aide des routines proposées par la référence [16]. Des détails sur les relations de récurrence entre les
fonctions de Bessel et leurs dérivées sont disponibles dans les références [1, 10]. Les prédictions de ces
routines ont été comparées avec celles des routines de Numerical Recipies [90], les données tabulées de
la référence [1] et des calculs effectués avec MAPLE et MATHCAD. Il est important de noter qu’un
manque significatif de précision a été constaté pour les routines de Numerical Recepies et les calculs sous
MATHCAD.
1.2. Formulation pour le calcul numérique 27
(0)
Rapport : Dn
Ce rapport de fonctions de Bessel peut être calculé à partir de la relation de récurrence ascendante :
(1)
(0) n/ρ − Dn (ρ)
Dn+1 (ρ) = Dn(0) (ρ) (3)
(1.71)
n/ρ − Dn (ρ)
Il est calculé uniquement pour l’interface externe (ρ = mext xL ) à partir de la valeur initiale :
(x) (j)
Rapports : Dn , Rn
(j)
En remplaçant les dérivées à l’ordre n + 1 en fonction de leur valeur à l’ordre n, Rn peut se calculer
par récurrence ascendante avec n = 1, . . . nstop − 1 et j = 2, 3, . . . , L :
(1) (2)
n
− Dn (mj xj−1 ) mnj xj − Dn (mj xj )
(j) mj xj−1
Rn+1 = Rn(j) (1.73)
(2) (1)
n
mj xj−1 − Dn (mj xj−1 ) mnj xj − Dn (mj xj )
(x)
(j) D0 (mj xj−1 ) J0 (mj xj−1 ) Y0 (mj xj )
R0 = (x)
= (1.74)
Dn=0 (mj xj ) Y0 (mj xj−1 ) J0 (mj xj )
1.2.5 Récurrence
La reformulation analytique n◦ 2 s’est révélée beaucoup plus stable que la première formulation. Dans
ce qui suit, nous présentons l’algorithme développé, pour le calcul des coefficients externes de diffusion,
à partir de la reformulation n◦ 2.
Algorithme :
Toutes les dérivées logarithmiques et autres rapports des fonctions de Bessel cylindriques sont cal-
culés au préalable. Les coefficients externes de diffusion peuvent alors être déterminés par récurrence
ascendante :
α(1) (1)
n = βn = 0.0 + i0.0 (1.75)
(1) (1)
Ha,n = Hb,n = Dn(1) (m1 x1 ) (1.76)
- Etape n◦ 2/Couche n◦ 2 :
(1) (1)
(2) (2) m2 µ1 Ha,n −m1 µ2 Dn (m2 x1 )
αn = Rn (1) (2)
m2 µ1 Ha,n −m1 µ2 Dn (m2 x1 )
(1) (1) (1.77)
(2) (2) m1 µ2 Hb,n −m2 µ1 Dn (m2 x1 )
βn = Rn (1) (2)
m1 µ2 Hb,n −m2 µ1 Dn (m2 x1 )
(1)
(2) Dn (m2 x2 )−α(2) (2)
n Dn (m2 x2 )
Ha,n = (2)
1−αn
(2) (1)
Dn (2) (2)
(m2 x2 )−βn Dn (m2 x2 )
(1.78)
Hb,n = (2)
1−βn
- Etape n◦ j/Couche n◦ j :
j−1 (1)
(j) (j) mj µj−1 Ha,n −mj−1 µj Dn (mj xj−1 )
αn = Rn j−1 (2)
mj µj−1 Ha,n −mj−1 µj Dn (mj xj−1 )
(j)
j−1 (1)
(j) mj−1 µj Hb,n −mj µj−1 Dn (mj xj−1 )
(1.79)
βn = Rn j−1 (2)
mj−1 µj Hb,n −mj µj−1 Dn (mj xj−1 )
28 Chapitre 1. Diffusion électromagnétique par un cylindre multicouche
(1)
j Dn (mj xj )−α(j) (2)
n Dn (mj xj )
Ha,n = (j)
1−αn
j (1)
Dn (j) (2)
(mj xj )−βn Dn (mj xj )
(1.80)
Hb,n = (j)
1−βn
- Dernière étape/Etape n◦ L/Couche n◦ L :
L−1 (1)
(L) (j) mL µL−1 Ha,n −mL−1 µL Dn (mL xL−1 )
αn = Rn L−1 (2)
mL µL−1 Ha,n −mL−1 µL Dn (mL xL−1 )
(L)
L−1 (1)
(j) mL−1 µL Hb,n −mL µL−1 Dn (mL xL−1 )
(1.81)
βn = Rn L−1 (2)
mL−1 µL Hb,n −mL µL−1 Dn (mL xL−1 )
(1)
(L) Dn (mL xL )−α(L) (2)
n Dn (mL xL )
Ha,n = (L)
1−αn
(L) (1)
Dn (L) (2)
(mL xL )−βn Dn (mL xL )
(1.82)
Hb,n = (L)
1−βn
(L) (L)
Au final, les Eqs. (1.52) et (1.53) et l’évaluation des fonctions Ha,n et Hb,n , permettent de déterminer
les deux coefficients externes de diffusion : anII et bnI .
2 2 ∞
2 2 ∞
Qsca
I = x |b0I | + 2 |b0I | , QextI = x Re b 0I + 2 b 0I
n=1 n=1 (1.83)
2 2 ∞
2 2
∞
Qsca
II = x |a0II | + 2 |a0II | , QI = x Re a0II + 2
ext
a0II
n=1 n=1
1.4 Validations
1.4.1 Cas limite du cylindre homogène
(L) (L)
Dans le cas d’une couche unique nous avons : Hb,n = Ha,n = Dn1 (mL xL ). Pour un cylindre de
paramètre de taille xL = x placé dans un milieu défini par mext = m et µL = µext , l’équation (1.54)
devient :
J (mx) J (x)
Jn (x) m Jnn (mx) − Jnn (x)
bnI = (1) (1) (1.86)
Hn (x) mL Jn (mx) − Hn(1) (mx)
Jn (mx) Hn (mx)
1.4. Validations 29
Cette dernière expression est strictement identique à celle donnée par la référence [10] dans le cas du
cylindre homogène. On trouve le même accord pour anII .
Sphère Cylindre
ψn → Jn
(1.88)
χn → Yn
(1)
ξn → Hn
Sphère Cylindre
Dn1 = ψn
ψ Jn
→ Jn
n
Dn2 = χn
χ Yn
→ Yn
n
(1) (1.89)
Dn3 = ξn
ξ Hn
→ (1)
n Hn
(0)
Dn = ψ χn
n
→ Jn
(1)
Hn
(x)
Dn = ψ n
χn → Jn
Yn
En utilisant ces correspondances, il est facile de montrer que les expressions des coefficients de diffu-
sion externe TM et TE du cylindre multicouche, équations (1.52-1.53), sont effectivement formellement
identiques à celles obtenues pour la sphère [109, 71].
La figure 1.2 présente une simulation du profil d’indice du coeur d’une fibre multimode Corning à
gradient d’indice et de diamètre de coeur : 85 µm ; indice du milieu extérieur mext = 14 . La figure 1.3
présente l’évolution des deux sections efficaces lorsque l’on raffine la discrétisation du profil (augmentation
du nombre de couches) : L = 2 → 5000. L’épaisseur d’une couche dr est également donnée en fraction
de la longueur d’onde moyenne λ/dr = 2 → 50. On constate que les sections convergent rapidement vers
une même valeur. Pour un incrément en dr = λ/10 la différence est inférieure à 0.1%, pour dr = λ/40,
la différence est de moins de 0.006%. La figure 1.4 montre les diagrammes de diffusion calculés pour une
fibre optique multimode Corning SX+ complète (coeur et cladding), placée dans l’air.
La stabilité des algorithmes développés est telle que, en discrétisant finement le profil d’indice, on
peut simuler les propriétés de diffusion d’un cylindre à profil d’indice continu.
4 Le cladding n’est donc pas simuler ici, voir plus loin pour sa prise en compte.
30 Chapitre 1. Diffusion électromagnétique par un cylindre multicouche
1,51
<m>=1.48663
1,50 mmax=1.5 <λ>=λ0/<m>=0,425661
Max{dm/dr}=0.00188
1,49
1,47
1,46 mmin=1.46
1,45
1,0
-50 -40 -30 -20 -10 0 10 20 30 40 50
Rayon [µm]
Fig. 1.2 – Profil d’indice d’une fibre optique multimode «Corning» à gradient d’indice avec 82µm de diamètre
de coeur.
Fig. 1.3 – Evolution des coefficients de diffusion en fonction du nombre de couches (raffinement de la
discrétisation du profil.)
Fig. 1.4 – Profil d’indice d’une fibre multimode à gradient (Corning type Sx+) : coeur et cladding, et évolutions
des diagrammes de diffusion correspondants. Le profil d’indice du coeur a été discrétisé en 3000 couches.
32 Chapitre 1. Diffusion électromagnétique par un cylindre multicouche
• Pour une fibre absorbante (i.e. métallique par exemple), figure 1.6 à droite, les diagrammes ont
une structure beaucoup plus simple. C’est la transparence des fibres de verre-E, dans le visible, qui est à
l’origine de la complexité des diagrammes obtenus.
• La figure 1.7 présente les même évolutions que précédemment mais pour la zone de diffraction
(diffusion avant : 0 < θ < 15◦ ) et la polarisation perpendiculaire. L’ intensité est normalisée par la valeur
obtenue à 0◦ (1er lobe de diffraction) ; avec à gauche, une fibre de verre-E homogène et à droite, une fibre
métallique.
Là encore, il est clair que la transparence des fibres de verre induit une grande complexité des diagrammes
de diffusion.
Fig. 1.5 – Evolution angulaire de l’intensité diffusée par une fibre de verre-E, homogène, éclairée sous incidence
normale. En abscisse, l’angle de diffusion évolue de 150 à 180◦ de gauche à droite. En ordonnée, le diamètre
de la fibre évolue de 0.02 à 42 micromètres de bas en haut. Les courbes de niveau représentent le logarithme de
l’intensité diffusée, normalisée par la valeur maximale, pour la polarisation perpendiculaire.
34 Chapitre 1. Diffusion électromagnétique par un cylindre multicouche
Fig. 1.6 – Idem figure précédente excepté : à gauche, pour la polarisation parallèle et, à droite, pour une fibre
métallique.
1.5. Sur la complexité des processus de diffusion 35
Fig. 1.7 – Idem figures précédentes excepté que l’angle de diffusion évolue de 0 à 15◦ (≈ zone de diffraction), la
polarisation est perpendiculaire et l’intensité est normalisée par la valeur obtenue à 0◦ (pic de diffraction) ; avec
à gauche, une fibre de verre-E homogène, et à droite, une fibre métallique.
36 Chapitre 1. Diffusion électromagnétique par un cylindre multicouche
Fig. 1.8 – Schéma de principe sur ”les ondes de surface” : le bord externe d’une fibre de diamètre 10.7235µm est
tracé en gras, l’ordre d’expansion des champs est schématisé par la stratification en grisé.
Fig. 1.9 – Comparaison de l’intensité diffusée calculée pour la série complète et lorsque l’on pose égale à
0 le coefficient de diffusion externe b61 .
imaginaire
0.0
0.6
0.4 -0.2
0.001µm
0.2 zoom
0.0
-0.2
Re(b61)
-0.4 0.03µm
Im(b61)
-0.6
10.1 10.2 10.3 10.4 10.5 10.6 10.7 10.8
Diamètre [µm]
Fig. 1.10 – Evolution de la partie réelle et imaginaire du coefficient b61 .
1.6. Conclusion 37
• Comme la figure 1.9, la figure 1.11 présente l’évolution de l’intensité du champ diffusé en fonction
du diamètre de la fibre. Mais, dans ce dernier cas, les séries sont tronquées de façon encore plus brutale :
a) en haut, la série est complète ; b) au milieu, tous les termes correspondant aux ondes de surfaces sont
annulés bn = 0.0 + 0.0i, n >= x ; c) en bas, tous les termes qui correspondent à des rayons qui impactent
sur la fibre sont annulés, bn = 0.0 + 0.0i, n < x. Pour des raisons de lisibilité de la figure, l’intensité du
cas c) a été multipliée par un facteur 10. Les évolutions obtenues sont cohérentes :
- Le cas b) montre une évolution, périodique à basse fréquence, caractéristique des ’basses fréquences’.
- Le cas c) semble rassembler l’essentiel des contributions des ondes de surface.
En fait, même visuellement, l’addition des évolutions b) et c) ne correspond que qualitativement à
l’évolution a). On peut avancer comme explication à cela que les termes bn tels que n = int(x) ± 1
rendent compte très vraisembablement à la fois des ondes externes et des ondes internes près de la sur-
face. La troncature que nous avons utilisé plus haut, dissocie trop ”violament” des processus qui sont
nécessairement couplés près de la surface (i.e. continuité des composantes tangentielles des champs).
Une analyse par la méthode de décomposition de Debye, plus rigoureuse, peut également aider à la
compré- hension de ces mécanismes [68].
• Toujours pour les mêmes paramètres, la figure 1.12 présente l’évolution de l’intensité du champ
diffusé en fonction du diamètre pour : i) une fibre homogène et différents indices de réfraction : m =
1.555, 1.55499, 1.5549; ii) une fibre avec un gradient d’indice, dont le profil linéaire m = 1.55, . . . , 1.555
est discrétisé en 1000 couches. Il semble qu’une variation légère de l’indice ne contribue qu’à faire glisser,
ou décaler, la structure de résonance, sans en changer fondamentalement la forme ou l’amplitude. Les
MDRs sont quant à elles beaucoup plus sensibles aux gradients internes, ce qui se comprend aisément.
• La figure 1.13 montre l’évolution de l’intensité du champ diffusé par une fibre à coeur, en fonction de
son diamètre et de la taille relative du coeur : 0, 10, 30, 60 et 90 % du rayon externe de la fibre. L’indice
de réfraction du coeur est pris égal à 1. On constate que, bien que les évolutions à basse fréquence soient
assez différentes, la position des résonances semble ne pas trop varier pour un coeur atteignant jusqu’à
60% du rayon externe de la fibre. Ici encore, ce résultat s’accorde avec la nature des résonances de surface.
1.6 Conclusion
Dans ce chapitre nous avons résolu le problème de la diffusion d’une onde électromagnétique plane,
harmonique, de polarisation linéaire, par un cylindre de section circulaire, d’extension infinie, composé de
L couches concentriques d’indice de réfraction complexe. Pour ce faire, nous avons utilisé une méthode
de séparation des variables (i.e. théorie de Lorenz-Mie) et une méthode de récurrence pour le calcul des
champs internes. L’onde incidente est d’incidence normale. Les algorithmes numériques et les codes de
calcul développés sont très stables. Les propriétés de diffusion de cylindres de grand paramètre de taille
(α > 1000), composés de plusieurs milliers de couches, peuvent ainsi être calculées. Ces travaux seront
largement utilisés pour simuler la réponse de l’interféromètre (Chapitre 2) et du diffractomètre (Chapitre
4) lors de la mesure des fibres de renfort.
De façon plus générale, ces travaux permettent également de simuler les propriétés de diffusion de
”particules” cylindriques rencontrées dans les laboratoires et l’industrie :
– Homogènes : fibres de verre, polymères, carbone, métalliques ; jets liquides axi-symétriques...
– A saut d’indice : fibres creuses (UV, IR, ”air-line”), fibres optiques à saut d’indice (coeur-cladding),
jets liquides enrobés, capillaires...,
– Stratifiées en indice (gradient discret ou continu) : fibres optiques , jets liquides et capillaires soumis
à un choc thermique, en cours d’évaporation, avec mélange diffusif, etc.
Ils s’appliquent également à d’autres domaines de longueurs d’ondes (ondes radar, infra-rouge loin-
tain...) et pourraient être généralisés dans le cas d’un éclairement par un faisceau gaussien [73, 36].
38 Chapitre 1. Diffusion électromagnétique par un cylindre multicouche
Fig. 1.11 – Calcul de l’intensité diffusée par une fibre lorsque l’expansion des champs est complète,
tronquée vers le bas ou vers le haut.
Fig. 1.12 – Effet de la partie réelle de l’indice de réfraction sur la structure de résonance (ondes de
surface) de l’intensité diffusée par une fibre homogène et une fibre à gradient d’indice.
Fig. 1.13 – Mise en évidence de l’effet d’un coeur d’air (bulle de gaz) sur la structure de résonance (ondes
de surface) de l’intensité diffusée par une fibre.
Chapitre 2
Modélisation de la réponse de
l’interféromètre ”FIBS”
Nous venons de résouldre le problème de la diffusion d’une onde plane par un cylindre multicouches.
Dans ce chapitre, le modèle ainsi développé va nous permettre de simuler la réponse de l’interféromètre
afin de l’optimiser, d’étudier sa sensibilité à divers paramètre et effets.
2.1 Introduction
Contraintes pour la mesure du diamètre des fibres de renfort
Nous avons établi au début de cette étude, d’un commun accord avec Saint-Gobain Recherche, qu’un
système de mesure et de contrôle en ligne, de la granulométrie des fibres de renfort, doit satisfaire aux
contraintes suivantes : i) La résolution sur le diamètre des fibres doit être telle que σD < 0.5µm pour
des fibres dans la gamme de taille D = 5 − 30µm ; ii) La résolution longitudinale/spatiale du système
doit pouvoir atteindre 1mm pour une fibre défilant à 50m/s, c’est-à-dire un taux d’acquisition de 50 000
mesures/s ; iii) Le système de mesure doit être compact afin de pouvoir être utilisé sur une filière de
laboratoire ou industrielle (place disponible très réduite). iv) La mesure doit être effectuée à grande
distance de la fibre pour ne pas perturber les conditions de fibrage.
Nos recherches bibliographiques [23, 51] sur les techniques de mesure du diamètre des fibres (optique,
polymère, renfort...), ont montré qu’il n’existait pas sur le marché1 de système de mesure capable de
satisfaire à toutes les contraintes précèdentes.
39
40 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
un traitement des données réduit au calcul d’une FFT... De ce fait, et compte tenu de notre expérience
dans le développement de cette technique, nous avons opté pour le développement d’un système de mesure
des fibres de renfort qui utilise le principe de l’IPD. La condition préalable à la réalisation expérimentale
d’un tel système est de modéliser théoriquement, et d’optimiser numériquement, sa réponse lors de la
mesure des fibres de renfort.
êz E//1 E1 E//s Es
êz
δ1 z δa P
êy1 E 1 E s
k1.êx=-k k.ê r=k êϕ
1 1
P
Da
E1
k1 Ω ψa>0 2
êx1
x α O
E2 α ψb<0
ϕ êz 1
k2 r êϕ Ω
êy1
êr Db
y
VZ
On considère ici que les dimensions du volume de mesure, formé par le croisement des deux faisceaux,
sont très supérieures au diamètre de la fibre. On peut alors utiliser les résultats établis dans le chapitre
précédent, pour des ondes planes. De manière classique, on suppose également que les faisceaux laser sont
décalés en fréquence (de la quantité νs ) par des cellules de Bragg ou par un réseau de diffraction tournant
[26].
On peut alors calculer le champ total diffusé par la fibre sur la surface d’un détecteur Da , d’angle
de polariseur δa , placé dans le plan (Oxy) à une grande distance r de la fibre, et dont l’axe optique fait
un angle ”d’élévation” ψa avec l’axe −x, voir la figure 2.1. Dans ces conditions, le détecteur collecte la
lumière diffusée aux angles de diffusion : faisceau n◦ 1, θ1 = ψa + α et pour le faisceau n◦ 2, θ1 = ψa − α :
Es (ψa ) = Es1 (ψa + α) + Es2 (ψa − α) ez + [E⊥s1 (ψa + α) + E⊥s2 (ψa − α)] eϕ (2.1)
Classiquement, ν1 /ν2 ≈ 1 avec pour la fréquence moyenne ν = (ν1 + ν2 ) /2. Le vecteur de Poynting total
2.2. Modélisation de la réponse de l’interféromètre 41
s’écrit alors :
k 2 2
Ss = E//s1 + E//s2 + |E⊥s1 |2 + |E⊥s2 |2 + 2Re E//s1 E//s2
∗ ∗
+ E⊥s1 E⊥s2 (2.2)
2µ0 ν
On peut le reformuler sous la forme suivante :
k
Ss = (G + Re {Hexp [−iωD (t)]}) (2.3)
2µ0 ω
avec 2 2
G = E//s1 + E//s2 + |E⊥s1 | + |E⊥s2 |
2 2
(2.4)
∗ ∗
H = 2 E//s1 E//s2 + E⊥s1 E⊥s2 (2.5)
νD = (k1 − k2 ) · r + νs = −4VY (t) sin (α) /λ + νs (2.6)
Dans les équations (2.4) et (2.5) les expressions des champs peuvent être explicitées à partir des
résultats du chapitre précèdent. Pour le faisceau 1 et pour le détecteur Da , on a par exemple :
3πi 2
E//s1 = exp cos2 δa T1 (ψa + α) E1 cosδ1 + sinδa cosδa T2 (ψa + α) E⊥1 sinδ1
4 πkr
(2.7)
3πi 2
E⊥s1 = exp sinδa cosδa T1 (ψa + α) E1 cosδ1 + sin2 δa T2 (ψa + α) E⊥1 sinδ1
4 πkr
Dans l’équation (2.6) le produit scalaire (k1 − k2 ) .r a été réécrit sous une forme qui dépend explicite-
ment du temps. Quand la fibre se déplace dans le plan (Oxy), le produit scalaire varie proportionnellement
à la composante de vitesse Vy (t) de la fibre. Lors du fibrage, toute variation de la «fréquence hétérodyne
Doppler» νD pourra être associée à une fluctuation de la position de la fibre suivant l’axe y. De plus,
du fait du décalage en fréquence des faisceaux, νs = (ν1 − ν2 ) > 0, la lumière diffusée par la fibre est
modulée dans le temps même lorsqu’elle est fixe dans le volume de mesure.
Remarque : Il peut être utile d’analyser le degré de polarisation linéaire de la lumière diffusée par la fibre
(cas d’une fibre biréfringente : § 2.5.5). Pour ce faire, une partie de la lumière qui est diffusée par la fibre
doit traverser un polariseur linéaire avant d’être focalisée sur l’élément photo détecteur. Sur la figure 2.1,
pour le détecteur Da , l’axe optique du polariseur P fait un angle δa avec ez . Au final, le champ diffusé et
collecté s’écrit après la traversée du polariseur linéaire :
Es cos2 δa sinδa cosδa Es
= (2.8)
E⊥s P sinδa cosδa sin2 δa E⊥s
Pour un détecteur réel, l’équation (2.3) doit être intégrée sur l’angle solide Ω de collection du détecteur.
Sous incidence normale, le calcul se réduit à une intégration sur le domaine angulaire [ψa − Ω/2, ψa + Ω/2].
Si la réponse du détecteur est linéaire (i.e. signal électrique proportionnel au signal optique), le signal
électrique produit par ce détecteur est de la forme 2 :
I (t) = |Ss (t)|dψ (2.9)
Ω
Fig. 2.2 – Exemple de simulations des caractéristiques des signaux Doppler produits par des fibres de
verre-E homogènes : déphasage, visibilités et piedestal (intensité diffusée) des deux signaux.
Nous introduisons à présent les quantités caractéristiques des signaux Doppler : le Piedestal P, la
Visibilité V et la phase φ[62] :
P = k GΩ / (2µ0 ω) (2.12)
V = 2 Hi 2Ω + Hr 2Ω / GΩ (2.13)
φ = tan−1 (Hi Ω / Hr Ω ) (2.14)
Finalement, l’équation (2.10) peut être mise sous la forme caractéristique d’un signal Doppler laser :
La phase φ est une phase absolue, elle ne peut pas être mesurée directement. C’est pourquoi, dans le
principe de la technique de l’IPD, deux détecteurs Da et Db ou plus (voir également [79, 82]) sont utilisés
pour mesurer la différence de phase ∆φab = φa − φb entre deux signaux Doppler Ia (t) et Ib (t). Ces
derniers sont obtenus à l’aide de deux détecteurs placés à des angles de diffusion différents ψa et ψb :
Ia (t) = Pa [1 + Va cos (2πνD t + φa )]
(2.16)
Ib (t) = Pb [1 + Vb cos (2πνD t + φb )]
La dépendance de la différence de phase avec le diamètre de la particule diffusante, φab = φa − φb , sera
appelée PDRS3 dans ce qui suit. D’après le principe de l’IPD, le diamètre de la particule (la fibre dans
3 Phase Diameter RelationShip
2.2. Modélisation de la réponse de l’interféromètre 43
le cas présent) est déduit de la mesure de la différence de phase ∆φab et de la connaissance «théorique/
numérique» de la PDRS. La figure 2.2 montre, à titre d’exemple, l’évolution du déphasage, de la visibilité
et du piedestal des signaux Ia et Ib , pour différents diamètres de la fibre. La fibre est ici homogène et
son diamètre varie entre D = 0.01µm et D = 70µm . Le déphasage φab évolue de manière plus ou moins
linéaire, par morceaux. Cette linéarité chute de concert avec la visibilité des signaux. Le premier saut du
déphasage (modulo 2π) apparaı̂t pour D ≈ 23µm. La visibilité évolue de manière beaucoup plus complexe
( i.e. ≈ fonction de Bessel). Elle passe par un minimum pour D ≈ 44µm. Le piedestal évolue de manière
plus ou moins linéaire. Ces évolutions ont été calculées pour les paramètres optiques suivants : α = 1.11◦ ,
Ω/2 = 5.15◦ , m = 1.555, ψa = 162.5◦ , ψb = −166.73◦, λ = 0.6328µm et une polarisation perpendiculaire
(système F IBS).
La validation des modèles précédents a été réalisée en 3 étapes. Premièrement, nous avons comparé les
valeurs numériques obtenues pour un cylindre homogène avec celles obtenues pour une sphère [62, 71]. Les
résultats étaient cohérents, compte tenu des différences prévisibles. Deuxièmement, nous avons comparé
nos résultats avec ceux de la littérature (fibres homogènes et d’indice inférieur) : ceux de Schaub et
al. [94], basés sur la Théorie de Lorenz-Mie, et ceux de Mignon et al., basés sur l’optique Géométrique
[58]. Nos résultats concordaient parfaitement avec ceux de Schaub et al. L’accord avec les prédictions
de l’optique géométrique n’était que très qualitatif, voir à ce propos le § 2.2.2. Au final, différents tests
expérimentaux ont démontré la validité de nos modèles de diffusion (§ 3.4).
p=2 θ2
τ ,
τ
r
h
,
τ , p=1
, τ τ
y1 r h τ θ1
,
τ τ
,
τ
L
p=0 θ0
Fig. 2.3 – Optique géométrique : parcours des rayons lumineux à travers une particule réfractante.
diffusion θp du rayon p, par rapport à la direction incidente, s’écrit [105] :
L’énergie véhiculée par les différents rayons peut être calculée au moyen des coefficients de Fresnel. Ces
coefficients en amplitude s’écrivent, pour les deux polarisations, perpendiculaire et parallèle :
sin τ − m sin τ m sin τ − sin τ
r⊥ = r = (2.19)
sin τ + m sin τ m sin τ + sin τ
44 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
Fig. 2.4 – Optique géométrique : suivi des rayons qui intéragissent avec une fibre dans le plan.
I0 D2 2⊥ γ I0 D2 2 γ
I⊥ (θ) = I (θ) = (2.20)
r2 r2
où r est la distance des détecteurs au point de mesure ; I0 est l’intensité de référence de l’onde incidente ; γ
est un facteur qui rend compte de la dispersion de l’énergie du rayon incident et qui est liée à la courbure
locale de la surface de la particule, avec γ = sinτ /| dθ dτ |. La figure 2.4 montre le résultat du suivi d’une
centaine de rayons incidents, parallèles, qui impactent sous incidence normale sur une fibre de verre-E de
diamètre 100µm. Dans cette simulation, pour des raisons de clarté de la figure, le parcours des rayons n’a
été suivi que jusqu’à la deuxième réflexion interne (i.e. p = 3). Différents phénomènes sont remarquables :
point de focalisation vers l’avant, rayons ayant subi une seule réflexion spéculaire (p = 0), présence de
caustiques à l’intérieur de la particule, zones de localisation du premier4 (p = 2) et second arc-en-ciel
(p = 3)... La figure 2.5 montre l’évolution de l’intensité de ces rayons en fonction de la distance qu’ils ont
parcourue à l’ intérieur de la particule. On constate que leur intensité décroı̂t rapidement mais qu’il existe
aussi, des périodicités, des plages de distances ou ”bandes” interdites, qui sont à l’origine de certaines
des périodicités observées dans les diagrammes de diffusion et notamment des MDRs (voir les sections
1.5.1).
4 Les optiques de collection de l’interféromètre FIBS sont situées dans cette région, voir le chapitre suivant
5 On néglige ici la dépendance temporelle et le déphasage dû à la distance particule-détecteur
2.2. Modélisation de la réponse de l’interféromètre 45
Intensité diffusée
0.6
0.5
réfraction, p=2
0.4 1er Arc-en-Ciel
0.3 réfraction, p=3
2eme Arc-en-Ciel...
0.2
0.15
0.1
0.06 0.04...
0.0
0 1 2 3 4 5 6 7 8
Bandes interdites... Parcours des rayons dans la fibre [en équivalent diamètre]
Fig. 2.5 – Optique géométrique : évolution de l’intensité des différents rayons diffusés en fonction de la
distance géométrique qu’ils ont parcourue à l’intérieur de la fibre.
]g 400
eD
[es 300
ah
P 200
100
0
0 10 20 30 40 50
Diamètre[µm]
Fig. 2.6 – Comparaison de la relation phase diamètre prédite par la théorie électromagnétique et par
l’optique géométrique (uniquement les rayons p = 2), pour ψa = −ψb = 165◦ , α = 0.855◦ et Ω = 3.13◦ .
où ϕ1 et ϕ2 représentent le déphasage des rayons réfractés (par rapport à un rayon non dévié par la
particule) du fait de leur parcours optique à l’intérieur de la particule (i.e. mL1 et mL2 ). La figure 2.6
compare la relation phase diamètre prédite par la théorie électromagnétique et l’optique géométrique
(uniquement les rayons p = 2), pour ψa = ψb = 165◦ , α = 0.855◦ et Ω = 3.13◦ . On constate que l’optique
géométrique ne prédit que très approximativement la pente de la relation phase-diamètre.
Pour tenter de prédire les oscillations haute fréquence de la PDRS, avec l’optique géométrique, il
conviendrait de sommer les contributions de tous les rayons qui sont susceptibles d’être diffusés dans la
direction ψa (i.e. pour p = 0, 1, 2...∞) et de prendre en compte les effets dus aux points focaux. Cependant,
il est important de noter que cette approche sera limitée par les limites mêmes de l’optique géométrique
[105]. Des effets purements ondulatoires comme la diffraction, le phénomène comme d’arc-en-ciel ou les
ondes de surface... ne seront jamais, ou que très grossièrement, décrits par l’optique géométrique.
La figure 2.7 présente l’évolution de la fréquence des oscillations de la PDRS, calculée avec la théorie
électromagnétique, pour différentes valeurs du rapport λ/m et pour la géométrie IPD correspondant à la
figure 2.6. Tous les points sont alignés suivant une droite de pente ≈ 1/2. Ces oscillations proviennent en
fait de l’interférence entre les rayons p ≥ 2. L’amplitude des oscillations est maximale lorsque ces rayons
sont déphasés entre eux de la quantité : ≈ n2 λm0 (la fibre se comporte comme une cavité résonnante).
46 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
Fig. 2.7 – Evolution de la fréquence des oscillations (hautes fréquences) de la PDRS, en fonction du
paramètre : longueur d’onde/indice réfraction.
ET (y, z, t) = E+1 ej[k+1 .r−(ω+ωs )t] + E−1 ej[k−1 .r−(ω−ωs )t] + E0 ej[k0 .r−ωt] (2.22)
2.3. Principe du fonctionnement en mode N-faisceaux 47
avec pour les vecteurs d’ondes respectifs exprimés dans le repère (Oyz) :
0 −sinα sinα
k0 = k k1 = k k−1 = k (2.23)
1 cosα cosα
Fig. 2.9 – Simulation de la distribution d’intensité dans un volume de mesure, type Doppler laser, formé
par l’interférence de trois faisceaux cohérents.
La figure 2.9 montre la distribution d’intensité dans le plan (Oyz), obtenue avec l’équation (2.25)
pour α = 1.11◦ et λ = 0.6328µm (i.e. les paramètres de l’interféromètre FIBS). Nous avons posé ici,
t = 0 et E0 = E1 . On constate que les franges sont toujours rectilignes mais qu’elles présentent une
modulation d’amplitude suivant l’axe optique (elles semblent apparaı̂tre, disparaı̂tre...). L’analyse de la
figure précédente et de l’expression (2.25) montre que le volume de mesure présente désormais trois
(1) (2)
périodicités spatiales : deux suivant l’axe y, avec des interfranges de δy = π/A et δy = 2π/A ; et une
périodicité suivant l’axe optique z avec une modulation spatiale de pas δz = 2π/B. Pour les paramètres
(1) (2)
de la figure 2.9 nous obtenons : δz ≈ 843µm, δy ≈ 8.15µm et δy ≈ 16.3µm.
(1)
Le réseau de franges de pas δy est identique à celui produit par un système à deux faisceaux (−1, +1)
(2)
d’angle total 2α. Le réseau de franges de pas δy est dû à la superposition des deux réseaux de franges
produits par les faisceaux (−1, 0) et (0, +1). Ce dernier réseau correspond approximativement à celui qui
serait produit par des faisceaux (−1, +1) qui se croisent avec un angle total α. La modulation suivant
l’axe optique ne peut pas être reproduite avec un système classique à deux faisceaux.
La figure 2.10 a) illustre l’évolution du profil d’intensité correspondant à la figure 2.9 en fonction de
y et pour z = 0, 0.2δz , 0.4δz , · · · , 1.0δz . En fait, ces évolutions simulent la forme d’un signal Doppler qui
serait produit par une fibre qui se déplace suivant y 6 . Les signaux Doppler sont maintenant bi-fréquence.
(1) (2)
La figure 2.10 b) présente l’évolution du profil d’intensité suivant l’axe z, pour y = δy et y = δy . On
observe que les modulations suivant z, des deux réseaux de franges, sont de même période mais déphasées
spatialement de 2δz .
6 ou qui est immobile si les franges défilent avec t > 0
48 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
La figure 2.10 c) montre l’évolution du rapport de modulation des deux réseaux de franges :
(1)
I δy , z, t
R (z, t) = R (∀y, z, t) = (2.27)
(2)
I δy , z, t
En guise de conclusion sur cette technique d’interférométrie à trois faisceaux, nous dirons qu’elle
présente deux avantages essentiels : elle permet d’étendre la dynamique de mesure d’un facteur 2 et
de mesurer la position de la fibre le long de l’axe optique. Sa résolution sur z est contrôlée par l’angle
entre les faisceaux, ce qui rend très aisé la production d’intervalles de modulation de seulement quelques
dizaines de micromètres (i.e. 2δz ≈ 80µm pour α = 5◦ ). Cette technique originale devrait trouver d’autres
applications (caractérisation des couches limites, des écoulements en micro canaux...) en permettant la
mesure ”instantanée” de profils de vitesse.
Fig. 2.10 – Profils d’intensité dans un volume de mesure formé par trois faisceaux cohérents et évolution
du rapport de modulation suivant z, t = 0.
Fig. 2.11 – Evolution des relations phase-diamètre en mode 3-faisceaux cohérents : (a, b), calcul rigoureux ;
(c) calcul simplifié et (d) différence entre les deux méthodes de calcul.
50 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
De manière générale, et dans le cas particulier de la figure 2.12, les meilleurs résultats sont obtenus
avec la PDRS basée sur un filtrage FFT et par polynôme d’ordre élevé (N.B. n = 9 donne les même
résultats que la méthode de filtrage).
Cette procédure d’optimisation, des considérations sur la dynamique du système et sur les contraintes
expérimentales 9 , nous ont permis de fixer : pour le demi-angle entre les faisceaux, la valeur α = 1.11◦ ;
et pour l’ouverture des détecteurs, la valeur Ω/2 = 5.1◦ . Ces valeurs correspondent à celles utilisées par
le système FIBS et sont imposées dans ce qui suit.
• Configuration A, polarisation parallèle, les deux détecteurs sont localisés en diffusion avant, de
part et d’autre de l’axe optique. L’optimum est obtenu pour ψa = 75.8◦ et ψb = −72.8◦. Dans ce cas, le
Lin
processus de diffusion principal est la réfraction du premier ordre (p=1) et l’on a σab−D = 0.091µm avec
◦
b = 22.09 /µm.
• Configuration B, un détecteur est localisé dans la région du premier arc-en-ciel ψa = 163.0◦ (p=2,
une reflexion interne) et le second est localisé vers l’avant ψb = −33.0◦ (p=1). Cette configuration est
Lin
particulièrement encombrante et présente un résolution faible avec σab−D = 0.37µm et b = 19.33◦ /µm.
• Configuration C, les deux détecteurs sont légèrement vers l’arrière, de façon plus ou moins
symétrique par rapport à l’axe optique. Ces régions angulaires correspondent à la ”zone noire d’Alexan-
der” (p=0, réflexion spéculaire dominante), ψa = 125.6◦ et ψb = −115.3◦. Cette configuration est celle
Lin
qui offre la moins bonne résolution : σab−D = 0.75µm avec b = −15.27◦/µm.
• Configuration D, les deux détecteurs sont localisés vers l’arrière, dans la région de l’angle d’arc-en-
ciel, de façon plus ou moins symétrique. Cette configuration est de loin la plus compacte avec ψa = 162.5◦
et ψb = −166.7◦, et elle offre une résolution acceptable : σab−D Lin
= 0.37µm, σab−D Lin
= 0.54µm avec
◦
b = 16.75 /µm.
Fig. 2.12 – Evolution de la PDRS du système FIBS, configuration (D) optimisée, et des PDRS filtrées :
linaire, polynôme d’ordre n = 9 et par FFT.
Fig. 2.13 – Schéma de principe des quatre configurations optiques ”intéressantes”. Les configurations A
et D ont été testées expérimentalement.
52 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
180
165
(D)
150 (B) +0.37
+0.37
135
(C)
120 -0.75
105
Perp.
σ ab-D
ψ1 90
0.75
0.40
0.25
[deg] 75 0.10
-0.91
60 (A) [µ m]
+0.091
45
Paral.
30
15
0
0 15 30 45 60 75 90 105 120 135 150 165 180
−ψ2 [deg]
Fig. 2.14 – Simulation des isovaleurs de l’erreur estimée en fonction de l’angle d’élévation de chacun des
détecteurs (ψa = ψ1 , ψb = ψ2 ), pour les polarisations perpendiculaire et parallèle.
42 0 .4
0
0.30 0.35
39
Diametre superieur [ µ m]
36
0.26
33
0.22
30
0.18
σ FFT
0.30
27
ab-D
24 0.60
0 .3
0.50
0.35
21 0.40
0.35
18
0 .4 0
0.26 0.30
15 0.26
0.22
0.
12 0.18
50
[ µ m]
0.
9
55
0.60
6
6 9 12 15 18 21 24 27 30 33 36 39 42
Diametre inferieur [ µ m]
Fig. 2.15 – Simulation des isovaleurs de l’erreur estimée en fonction de la gamme des diamètres
considérés. Cette carte a été calculée pour le système FIBS.
2.5. Sensibilité à divers paramètres et effets physiques 53
800
260
700
600 255
400 N
Φ (D)
17.0 17.1 Trans
Zoom Φ (D)
300
200
N Trans
100*(Φ (D)- Φ (D))
100
0
σ=0.28°
-100
10 15 20 25 30 35 40
D [µm]
Fig. 2.16 – Système FIBS : influence sur la PDRS de la fonction de transmission anglulaire des optiques
de collection.
Du fait des résonances, la résolution attendue sur les diamètres dépend également de la gamme des
tailles considérées10. La Figure 2.15 présente une carte des iso valeurs de l’erreur estimée σab−D
NF F T
pour
une configuration optique de type D11 , en fonction des valeurs du diamètre minimum et du diamètre
maximum. La largeur maximale de la distribution est de D = 5−45µm. A titre d’exemple, on peut déduire
FFT
de cette carte que l’erreur moyenne sur les diamètres est égale à σab−D = 0.30µm, 0.26µm, 0.30µm et
0.35µm pour des fibres dont le diamètre évolue linéairement sur les intervalles [Dmin , Dmax ] ∼ : [12, 14],
[12, 17], [12, 30], [12, 42] [µm, µm].
cohérents !).
54 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
Fig. 2.17 – Schéma de principe très simplifié sur l’influence des angles d’inclinaison de la fibre, γ et τ ,
sur la forme adoptée par le cône de diffusion.
hauteur angulaire réelle des détecteurs de l’interféromètre FIBS est de seulement ±0.7◦ (i.e. les fibres
inclinées de plus de 2γ = 0.7◦ ne peuvent donc pas être détectées).
Pour des fibres métalliques, tout ce passe comme si le ”diamètre” de la fibre diminuait avec l’angle
d’inclinaison γ. Une régression sur les variations de la phase montre que cette décroissance est en cosinus
de l’angle γ. Pour les fibres de verre, les tendances sont les mêmes. Cependant, du fait des MDRs, les
évolutions sont plus complexes. Pour des angles d’inclinaison inférieurs à γ ≈ 1◦ , les variations de la
phase sont extrêmement faibles. C’est pourquoi, pour l’interferomètre FIBS, on peut négliger l’influence
de cet angle d’inclinaison.
La figure 2.17 schématise également les effets de l’angle d’inclinaison τ (la fibre est dans le plan Oyz).
Dans ce cas, le cône de lumière difusé par la fibre prend la forme d’une nappe plane, comme pour la
figure 2.1, mais cette nappe est inclinée par rapport au plan des faisceaux Oxy. Expérimentalement, nous
avons constaté que cet angle d’inclinaison pouvait être facilement détecté en comparant la visibilité des
signaux mesurés à celle prédite pour une fibre non inclinée (§ 3.3). C’est pourquoi, nous pouvons négliger
cet effet d’inclinaison dans ce qui suit.
Ces angles d’inclinaison n’ont pas de répercution significative sur la résolution du système F IBS. En
revanche, lorsqu’ils sont excessifs, du fait de fortes vibrations de la fibre par exemple, ils entraı̂nent une
forte chute du taux de validation du système.
Le verre utilisé pour produire les fibres de renfort est un produit de qualité ”industrielle”. La com-
position du verre, de même que la température de la fibre au point de mesure, peuvent être sujets à des
modifications ou à des variations importantes. Ces dernières peuvent influer sur l’indice de réfraction du
verre et donc, sur la mesure du diamètre.
Deux options sont alors possibles : soit on néglige l’influence de ces effets sur la réponse du système
FIBS, soit on recalcule la PDRS pour le nouvel indice supposé de la fibre. La figure 2.19 simule l’évolution
Lin
de l’erreur σab−D pour ces deux cas de figure. Au premier ordre, les différences sont minimes. La réponse
moyenne de l’interféromètre est très peu modifiée pour des indices compris entre m = 1.52 et m = 1.56.
L’erreur augmente significativement en dehors de ce domaine. Ce dernier point s’explique par le fait que
les détecteurs de l’interféromètre ont été placés dans la région de l’angle d’arc-en-ciel produit par des
fibres d’indice m = 1.555 et dont le diamètre est inférieur à D < 35µm.
Des variations de la température de la fibre, au point de mesure, n’auront donc que très peu d’influence
sur la réponse moyenne du système FIBS. En effet, une variation de température de ∆T = 100◦ C
n’entraı̂ne qu’une variation de ∆m ≈ ±0.0001 de l’indice13 .
Fig. 2.18 – Simulation de l’influence de l’angle d’inclinaison γ sur les propriétés du signal Doppler Ia
produit par des fibres de diamètre D = 14µm, 17µm, 20µm : métalliques et constituées de verre-E, pour
Ω/2 = 5.15◦ , ψa = 165◦ , α = 1.11◦ .
Fig. 2.19 – Sensibilité du système FIBS à une variation de l’indice de réfraction de la fibre mesurée.
Corpus et al. [18] ont trouvé un bon accord entre les prédictions de leur modèle et des données
expérimentales obtenues pour des fibres composées d’un verre NBS-170 ”soda-lime-silica”. A partir de
ce travail, on peut dériver une formule analytique pour la dépendance de l’indice de réfraction avec le
diamètre de la fibre (taux de refroidissement) :
βT (θT + 2) ln (D) + αT (θT + 2) − 1
m (D)T = (2.28)
βT (θT − 1) ln (D) + αT (θT − 1) − 1
Pour estimer les valeurs des coefficients αT , βT et θT pour le verre-E, nous avons utilisé les données
expérimentales de Clementin-de-Leusse [19]. Une procédure de moindres carrés a été utilisée pour ”fitter”
les données expérimentales avec l’équation théorique précédente. Ce qui nous a permis d’obtenir les valeurs
suivantes : αT = 119.68047, βT = −82.49331 et θT = 3.10733. La figure 2.20 a) montre l’évolution de
l’indice de réfraction prédite par l’Eq. 2.28. Il varie ainsi de m = 1.5407 → 1.55573 pour une fibre dont
le diamètre évolue de D = 5 → 45µm.
1.556 1.5550
a) b)
1.554 1.55573 1.55497
Indice de réfraction, partie réelle [-]
1.5549
1.550 1.5547
1.544 1.5544
1.5543
1.542 1.54071 1.55421
1.5542
5 10 15 20 25 30 35 40 45 5 10 15 20 25 30 35 40 45
Diametre, D [µm] Diametre, D [µm]
Fig. 2.20 – Evolution de l’indice de réfraction du matériau de la fibre : a) en fonction de la tension et b)
en fonction du taux de refroidissement.
L’influence du taux de refroidissement sur la réponse du système FIBS peut être évaluée à l’aide de
la figure 2.12. Cette figure illustre la PDRS obtenue lorsque l’équation 2.28 est utilisée ΦTab (D), ou non
ΦN T N
ab (D) (indice de réfraction constant, m=1.5550). La différence entre les deux réponses,Φab (D)−Φab (D),
est également présentée. Pour la gamme des tailles, D = 5 − 43µm, l’écart-type de la différence de phase
◦
ΦTab (D) − ΦNab (D) est de 3.99 . Ce qui signifie que, si l’on néglige la dépendance de l’indice de réfraction
Lin
avec le taux de refroidissement, on ajoute une erreur moyenne de ∆σab−D = 0.24µm à l’estimation du
diamètre moyen de la fibre.
La dépendance de l’indice de réfraction avec le diamètre semble avoir une forte influence sur la
structure des résonances : un décalage en fréquence (qui semble être un décalage de phase) est observé
pour les diamètres D < 12µm et D > 30µm, voir la figure 2.12. Pour de faibles diamètres, cette évolution
s’explique par la forte évolution de l’indice correspondant. Pour les grands diamètres, on peut avancer que
l’augmentation de l’épaisseur optique des particules compense la décroissance de l’évolution de l’indice
avec le diamètre. Pour les diamètres intermédiaires, les différences sont plus faibles. Une explication
plausible est que les structures de résonances (des deux détecteurs) sont quasiment ”en phase” pour ces
diamètres.
2.5. Sensibilité à divers paramètres et effets physiques 57
En conclusion, pour simuler la réponse de l’interféromètre FIBS, nous prendrons désormais en compte
la dépendance de l’indice avec le diamètre, dû au taux de refroidissement.
où ε⊥ et ε sont deux constants diélectriques (égales lorsque la fibre est coupée) et DX , DY et DZ les
trois composantes du vecteur déplacement électrique.
En négligeant l’absorption du verre dans le domaine visible, l’indice extraordinaire de la fibre m et
son indice ordinaire, m⊥ , peuvent être simplement reliés aux constantes diélectriques par :
1/2
m⊥ = (µε⊥ /µ0 ε0 )1/2 , m = µε /µ0 ε0 (2.30)
Pour de faibles contraintes on peut considérer que la biréfringence de la fibre m − m⊥ est proportion-
nelle à la tension appliquée, ce qui permet d’écrire :
m⊥ σ (T0 ) F
m − m⊥ = −1.29 104 (2.31)
D2
où F est la tension d’étirement de la fibre, en grammes ; σ (T0 ) le coefficient de stress optique, en MPa −1 ;
D représente le diamètre de la fibre, en micromètres. σ (T0 ) et F dépendent des conditions de fibrage et
principalement de la température de fibrage qui influent fortement sur la viscosité du verre. Dans le cas
du fibrage du verre E, une température de fibrage de T0 = 1200◦C, le coefficient de stress optique vaut
σ (T ) ≈ 6.1×10−7M P a−1 [41]. En utilisant l’Eq. (2.31) on trouve que la différence entre les deux indices,
0
m − m⊥ , est de 2.4 × 10−4 et 4.3 × 10−5 pour des fibres de D = 5µm et D = 10µm. Théoriquement,
la prise en compte de cette biréfringence est réalisée en remarquant que (§ 2.2.1) le champ électrique
diffusé Es dépend uniquement du coefficient externe de diffusion bnI et, E⊥s dépend uniquement du
coefficient externe de diffusion anII . C’est pourquoi, la prédiction des propriétés de diffusion d’une fibre
biréfringente uniaxe peut être obtenue au moyen des équations introduites dans le chapitre 1 et dans la
section 2.2.1, en posant bnI λ, D, m = m et anII (λ, D, m = m⊥ ).
La figure 2.20 b) présente l’évolution de l’indice extraordinaire m en fonction du diamètre (i.e. m⊥ =
cte = 1.555 − 0.0i). Elle indique que l’effet de la tension est plus important pour les fibres de faible
diamètre(D < 15µm). Le changement de l’indice de réfraction est de l’ordre de ∆m = +0.0008 pour des
fibres dans la gamme D = 5 − 45µm.
L’influence de la tension de fibrage sur la PDRS du système FIBS peut être estimée à partir de la figure
2.22. Cette figure présente l’évolution de la différence entre : la PDRS calculée en utilisant l’équation 2.35
(fibre sous tension, i.e. durant le procédé de fibrage) et la PDRS calculée sans tension (fibre coupée, i.e.
mesure en laboratoire) : ΦTab (D) − ΦN ab (D). Pour cette gamme de diamètres, D = 5 − 7µm, on observe
des pics importants. La différence de phase varie entre +24◦ et -45◦ . Ces pics apparaissent pour des
couples très particuliers (indice, diamètre). L’écart type moyen total est de 3.29◦ . Ce qui signifie que,
en négligeant les effets de la tension sur les propriétés optiques de la fibre, on surestime l’écart-type
58 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
Ε//1 z
ε//
x
ε O
Ε ε
1
y
Faisceau
Laser
F
Fig. 2.21 – Géométrie du modèle de diffusion par un cylindre biréfringent.
30
]
g
e
25 Φab (D)-Φab N (D) )
F
-
)
D
(
10
F ba
5
Φ
°
9
e
s
a
h
0 .2
3
=
p
σ
e
d
-5
e
-10
c
n
e
r
e
f
f
i
-46
D
5,0 5,2 5,4 5,6 5,8 6,0 6,2 6,4 6,6 6,8 7,0
Diametre, D [µm]
Fig. 2.22 – Comparaison de l’ évolution de la PDRS du système FIBS lorsque l’on prend ou non en
compte la dépendance de l’indice de réfraction avec la tension (biréfringence uniaxe).
Lin
moyen sur le diamètre de près de ∆σab−D = 0.20µm. Cette surestimation n’est plus, en moyenne, pour
D = 5 − 43µm, que de ∆σab−D = 0.06µm (i.e. écart-type moyen de la différence de phase : 0.94◦ ).
Lin
Expérimentalement, on constate que les bulles observées au niveau du bulbe sont évacuées avec le verre
et donc ”fibrées”... Elle conduisent de ce fait à la formation de fibres creuses. Ces fibres creuses posent
différents problèmes et notamment pour les applications de renforcement des cartes électroniques. En
effet, les fibres creuses peuvent induire la formation de ponts électriques (court-circuits). Les fibres
creuses ont également un indice de réfraction ”effectif” très différent de celui des fibres homogènes. Cette
différence d’indice peut être problématique pour les applications qui nécessitent des fibres ”invisibles”
(i.e. index-matching). On peut également penser que ces bulles perturbent l’écoulement au sein du bulbe,
qu’elles fragilisent le fil de verre et donc au final, qu’elles induisent des ”casses”. C’est pourquoi, si l’on
cherche à minimiser le taux de formation/présence de bulles, il faut pouvoir les détecter, les compter, et
éventuellement les caractériser (section, longueur d’étirement...) .
Simulations
Si l’on considère une bulle initialement sphérique, de diamètre Ds , entrainée par l’écoulement de verre
fondu à travers le téton, on peut envisager deux cas de figure : i) le débit diphasique (gaz+verre) est
conservé et la section totale de la fibre ne change pas, ii) seul le débit de verre est conservé. La figure
2.26 illustre les profils auxquels on peut s’attendre une fois la fibre figée14 . Nous avons alors :
– i) le diamètre externe De reste constant et le diamètre interne Di évolue ”d’une certaine façon”.
Dans ce qui suit nous imposerons arbitrairement une évolution sinusoı̈dale du profil (voir plus loin).
– ii) De et Di évoluent simultanément. A nouveau deux cas se présentent : on impose arbitrairement
le profil d’évolution de Di ou celui de De (ces deux profils sont reliés par la conservation du débit)
avec pour une fibre de diamètre initial, D0 :
2 2
De (t) = D0 + Di (t) Di (t) = De2 (t) − D02 (2.32)
La figure 2.25 montre l’évolution simulée du profil d’une fibre creuse lorsque :
– Profil n◦ 1 : Le débit diphasique est constant (i.e. diamètre externe est constant, le diamètre interne
varie sinusoı̈dalement.) ;
– Profil n◦ 2 : Le débit de verre est constant avec un profil interne sinusoı̈dal ;
– Profil n◦ 3 : Le débit de verre est constant avec un profil externe sinusoı̈dal.
Fig. 2.23 – Photographie du bulbe lors du passage d’une bulle de gaz (en bas, sur le côté droit), qui va
entraı̂ner la formation d’une fibre creuse.
Pour chacun des profils, 9 cas (i.e. bulles) sont considérés. On augmente sucessivement de 11% le diamètre
maximal. Par exemple, pour le Profil n◦ 1, Di évolue de 0.11D0 à 0.99D0 . Le diamètre initial de la fibre
homogène est pris égal à D0 = 8.24µm. En abscisse on trouve le paramètre d’itération, qui est assimilable
à un temps. La dilatation ou la contraction de l’échelle de ce paramètre permet donc de rendre compte
14 Nous ne disposons d’aucune information sur la nature réelle du profil ou la loi de conservation suivie par le débit
diphasique. Cependant, des observations qualitatives, sur l’écoulement de jets visqueux isothermes, semblent indiquer qu’au
passage d’une bulle le diamètre externe du jet augmente et que la bulle adopte un profil type sinoı̈dal (une seule période).
60 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
de vitesses de fibrage différentes et/ou de bulles initiales de plus ou moins grande dimension. En fait, les
observations expérimentales montrent que les bulles ont un rapport d’aspect beaucoup plus important
(i.e. plus allongées/étirées) [77].
L’intégration du profil de chaque bulle permet de déduire le diamètre de la bulle sphérique correspondante.
A titre indicatif, pour la 5eme bulle du Profil n◦ 2, on trouve Ds ≈ 106, 226, 285µm en faisant l’hypothèse
qu’elle est étirée sur une longueur de 0.1, 1 ou 2m.
Fig. 2.24 – Schéma de principe sur l’évolution du profil de la fibre creuse lorsque l’on utilise une hypothèse
de : (a) conservation du débit diphasique (diamètre externe fixe) ; (b) conservation du débit de verre
(diamètre externe et diamètre interne variables.)
-8
8 Profil n°2: Débit verre constant, profil interne sinusoidale
4
0
D0=8.24 µm
-4
-8
8 Profil n°3: Débit verre constant, profil externe sinusoidale
4
0
-4
D0=8.24 µm
-8
0 2000 4000 6000 8000 10000 12000 14000 16000 18000
Evolution "temporelle" [U.A]
Fig. 2.25 – Simulation de différents profils de fibres creuses (i.e. bulles).
450
360
270
180
90
0
630 Profil n°3
540
450
360
270
180
90
0
0 2000 4000 6000 8000 10000 12000 14000 16000 18000
Evolution "temporelle" [U.A]
(a)
2.0 Profil n°1
1.5
1.0
0.5
Intensité des signaux, Normalisée
0.0
2.0 Profil n°2
1.5
1.0
0.5
0.0
2.0 Profil n°3
1.5
1.0
0.5
0.0
0 2000 4000 6000 8000 10000 12000 14000 16000 18000
Evolution "temporelle" [U.A]
(b)
Fig. 2.26 – Simulation de la réponse de l’interferometre FIBS pour les différents profils de fibres creuses
(i.e. bulles). (a) Déphasage et (b) Intensité des signaux Doppler.
62 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
profils, les fluctuations du déphasage ne deviennent très importantes que pour des fibres creuses dont le
diamètre interne est supérieur à 0.4-0.5 fois le diamètre externe. Le déphasage pour le profil n◦ 3 est lui
très sensible à la présence d’une bulle (même ”petite”).
La détection de fluctuations de ce type, dans une série temporelle, doit donc permettre de détecter
les fibres creuses.
L y Valve
Chopper Air
supply
2a
2-fix points
-F x
winding z
Vf : drawing velocity
Une autre technique bien connue, pour la mesure de la tension, est celle de la corde vibrante. Cette
technique, valable pour une fibre de longueur L et aux extrémités fixes, consiste à mettre en vibration
transverse la fibre puis à mesurer sa fréquence de vibration :
15 la composition du verre...
2.6. Mesure de la tension de fibrage par la méthode de la corde vibrante 63
sinα dy (t)
νD (t) = νs − 4π (2.37)
λ dt
V, φ sont indépendants de la position de la fibre dans le volume de mesure alors que, pour des fibres
telles que D ωy , le piédestal dépend fortement de cette position et peut être modélisé simplement :
P (y) = P (0) I (y). Si la position transverse de la fibre est fixe dans le volume de mesure, nous avons
νD = νs et P (y) = cte, et l’on peut déduire le diamètre de la fibre de manière classique (§ 2.2.1).
Le principe de la mesure de la tension repose sur l’utilisation du profil gaussien d’intensité dans le
volume de mesure et notamment suivant y. Suivant nos hypothèses, lorsque la fibre vibre, sa position dans
le plan xy est donnée au temps t par y (t) = A1 cos (2πν1 t) + y0 (avec x = 0). A1 > 0 et y0 représentent
son déplacement maximal et sa position moyenne suivant l’axe y. Ce qui donne au final pour l’équation
du signal Doppler produit par cette fibre :
La dépendance temporelle de ces signaux est plutôt complexe. En fait, nous n’avons pas pu obtenir une
expression analytique pour la transformée de Fourier : S (ν) = T F [S (t)]. Dans l’équation (2.38) on peut
cependant remarquer que le premier terme ne dépend pas du temps ; le second terme correspond à des
évolutions à haute fréquence (du fait de la présence du décalage en fréquence νs ) dont l’influence devrait
être négligeable sur le spectre basse fréquence des signaux(ν1 νs ) ; le troisième terme correspond à des
évolutions à basse fréquence, autour de ν1 . Pour ce troisième terme, il est cependant difficile d’évaluer,
intuitivement, la contribution du terme exponentiel. C’est pourquoi, la figure 2.29 présente les résultats
du calcul numérique de S (ν) pour trois signaux Doppler simulés à l’aide de l’équation (2.38) et pour les
paramètres suivants : a) α = 0◦ et y0 = ωy /2 ; b) α = 0◦ et y0 = 0 ; c) α = 1.10◦ et y0 = ωy /2, avec
pour les trois cas : A1 = ωy /2, P (0) = V = 1, φ = 0 et ν1 = 200Hz. Il est à noter que, pour calculer
ces transformées de Fourier, nous avons utilisé les mêmes paramètres de traitement que pour analyser les
signaux expérimentaux, à savoir : νs = 90kHz, une fréquence d’échantillonnage de νdigit = 2048Hz et
64 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
N = 1024 échantillons. On remarquera également que, pour la mesure de la tension, les signaux Doppler
sont largement sous échantillonnés : ν1 νdigit νD . Les cas pour lesquels α = 0◦ correspondent
à un système optique à un seul faisceau (i.e. il n’y a pas de ”franges” dans le volume de mesure). Le
pic correspondant, dans le spectre de Fourier, est très facilement identifiable. Cependant, si la position
moyenne de la fibre se situe au centre du volume de mesure y0 = 0, la détermination de la fréquence
ν1 est biaisée (la fréquence mesurée est double de celle attendue), voir la figure 2.29 b. En fait, pour les
mesures de tension, on doit s’assurer que les oscillations de la fibre ne se produisent que dans une zone
où l’évolution du piédestal est monotone, figure 2.28. La position moyenne optimale correspond au points
d’inflexion de P (y), à savoir : y0 = ±ωy /2 avec A1 ωy /2 . Si cette condition est remplie, l’interféromètre
FIBS (i.e. α > 0) permet de déterminer sans biais cette fréquence de vibration. La contribution de la
fréquence de décalage νs se réduit alors à un bruit de fond qui ne pénalise pas la mesure, voir la figure
2.29 c).
Fig. 2.28 – Schéma de principe pour la mesure de tension : la fibre vibre suivant le profil gaussien
d’intensité de la sonde optique.
Expérimentalement, différentes façons d’exciter la fibre ont été testées (perturbations mécaniques,
acoustiques...). La méthode la plus simple s’est révélée être une technique de perturbation par un jet
d’air pulsé qui se propage suivant l’axe y. Dans notre montage, la vitesse du jet d’air, au niveau de
la fibre, est de l’ordre de 2m/s pour une durée de ≈ 50ms. Pour les plus petites fibres, D < 15µm,
l’efficacité de cette méthode décroı̂t rapidement. Il y a deux raisons à cela : la surface d’interaction
air/fibre décroı̂t et la vitesse de bobinage croı̂t (les perturbations sont plus rapidement «convectées» en
dehors de la zone d’étude L). Une façon de résoudre ce problème consiste à perturber la fibre avec un jet
d’air modulé en fréquence. Ceci peut être réalisé en hachant le jet d’air, initialement continu, à l’aide d’un
choppeur (optique dans notre cas). Le moteur du chopper est alors alimenté par une rampe de tension afin
d’augmenter progressivement sa vitesse de rotation et donc, la fréquence de modulation du jet d’air. La
fréquence d’oscillation de la fibre est alors associée à la fréquence d’excitation pour laquelle l’amplitude de
modulation des signaux Doppler est maximale (i.e. le pic le plus énergétique dans le spectre de Fourier).
La mesure simultanée du diamètre et de la tension de la fibre peut donc être réalisée à partir de la
séquence suivante :
- (1) Le diamètre de la fibre est mesuré, classiquement, à partir de la mesure du déphasage entre deux
signaux Doppler. La fréquence de décalage choisie ici est de νD = νs = 90kHz.
- (2) La vitesse de bobinage est mesurée par l’interféromètre au moyen de l’analyse fréquentielle du
signal de sortie du bobinoir (équipé d’une roue codeuse).
- (3a) D 15µm, un jet d’air pulsé est généré.
- (3b) D < 15µm, le jet d’air est mis en route et la rampe en fréquence du chopper démarre. Cette
rampe balaie un domaine de fréquence autour de la fréquence d’oscillation attendue. Cette rampe dure
approximativement ≈ 60s.
- (4) L’interféromètre numérise le signal de sortie d’un des deux détecteurs sur une période de 0.25 à
2s. Ce temps d’acquisition, plutôt long, est nécessaire si l’on veut déterminer avec une bonne précision
la fréquence de vibration de la fibre (basses fréquences).
- (4a) La numérisation démarre avec un retard de l’ordre 200ms, afin de s’assurer que le mode dominant
est bien le mode fondamental. En effet, les autres harmoniques s’atténuent rapidement.
- (4b) La numérisation est effectuée de manière périodique, par blocs, tout au long de la rampe du
chopper. Les spectres des signaux Doppler sont alors cumulés pour obtenir une fonction de transfert.
- (5) Cette procédure peut être répétée avec une cadence maximale de l’ordre de : une mesure simul-
tanée par seconde pour D 15µm et une mesure simultanée par minute dans le cas ou D < 15µm
16 Remarque : avec la procédure d’inversion par régression linéaire, ces grandeurs sont liées linéairement aux moments
statistiques des phases mesurées : Φ̄0 et σΦ0 . Pour une PDRS théorique linéarisée ΦT = aDT + b, on obtient en effet
D̄ = (Φ0 − b)/a) et σD = σΦ0 /a.
66 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
σDT ) et la loi de probabilité des fluctuations du diamètre de la fibre est gaussienne (voir la figure 7.5) :
2
1 1 D̄ − D
P (D) = √ exp − (2.39)
σD 2π 2 σ̄D
On définit un estimateur de type ”moindres carrés”17, où κ est un facteur de pondération qui permet
de privilègier la moyenne ou bien l’écart-type :
2 2
Γ D̄T , σDT = Φ̄ D̄T , σDT − Φ̄0 + κ σΦ D̄T , σDT − σΦ0 (2.40)
Fig. 2.30 – Séries temporelles centrées sur D̄ = 15µm : a) évolutions temporelles ; b) distributions des
diamètres, pour 3 cas : 1Hz, 0.5Hz et 0.1Hz.
Fig. 2.31 – Différence entre le diamètre moyen théorique et celui obtenu par inversion avec la méthode
linéaire et la méthode des distributions : a) Sans erreur de phase ; b) avec une erreur de phase de +1◦
et/ou une erreur sur l’écart-type de la phase de +0.5◦.
Pour tester cette méthode d’inversion, nous avons utilisé une série expérimentale filtrée passe-bas à
1Hz, 0.5Hz et 0.1Hz, figure 2.30. Avec le filtre à 1Hz les fluctuations du diamètre sont importantes et
17 L’ajout de moments d’ordre supérieur n’apporte pas d’avantages significatifs
2.7. Perspectives sur l’amélioration de l’inversion des phases 67
Fig. 2.32 – Différence entre l’écart-type du diamètre théorique et celui obtenu par inversion avec la
méthode linéaire et la méthode des distributions : a) sans erreur de phase et b) avec une erreur de phase
de +1◦ et/ou une erreur sur l’écart-type de la phase de +0.5◦ .
les distributions sont plutôt de forme gaussienne. Avec le filtre à 0.1Hz, l’évolution du diamètre apparaı̂t
”quasi-périodique” et la distribution des diamètres s’écarte significativement de la forme gaussienne (ce
qui affaibli notre hypothèse...) 18 . Ces séries temporelles sont ensuite ”décalées” de manière à obtenir un
diamètre moyen centré sur D̄ = 10µm, 12.5µm, 15µm, 17.5µm, 20µm et 22.5µm. On simule ainsi 3 × 6
séries temporelles pour lesquelles on peut calculer les moments statistiques réels. La figure 2.30 présente :
a) les séries centrées sur D̄ ≈ 15µm et b) les distributions correspondantes.
La figure 2.31 compare les résultats obtenus, sous la forme |D̄reel − D̄T |, lorsque la phase simulée
est inversée : a) directement, et b) pour simuler une erreur systématiqueon ajoute un décalage de phase
de +1◦ et/ou un décalage sur l’écart-type de 0.5◦ . La figure 2.32 présente les résultats obtenus pour les
écart-types, |σDreel − σDT |.
La méthode des distributions donne de meilleurs résultats que la méthode d’inversion linéaire, excepté
pour quelques rares cas pour lesquels la distribution réelle des diamètres diffère trop d’une distribution
gaussienne.
où les DT (i)19 représentent les diamètres théoriques associés à la phase φ(t) par la PDRS de l’appareil.
On peut maintenant calculer, à chaque instant, la probabilité d’occurence du diamètre DT (i) dans la
série expérimentale.
18 N.B.La PDF d’une variable harmonique n’est absolument pas gaussienne !
19 La PDRS est calculée pour i = 1, 2...N diamètres. L’incrément sur le diamètre conditionne nécessairement la résolution
de cette méthode d’inversion.
68 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
Diamètre [ m]
Cône initial
Temps [s]
Fig. 2.33 – Schéma de principe : formes du ”cône” limitant la variation maximale autorisée pour le
diamètre autour du temps t.
On pose, comme autre hypothèse, que le diamètre de la fibre évolue de manière continue (i.e. la résolution
temporelle du système de mesure doit être suffisante). On introduit alors la fonction de Densité de
Probabilité Jointe (DPJ) suivante :
!
j=+n/2
DP J {t, DT (i)i=1,2...N } = P {t + j∆t, DT (i)i=1,2...N } (2.43)
j=−n/2
où ∆t est l’intervalle de temps entre deux points de mesure consécutifs et n, le nombre de points de
mesure, ou voisins, pris en considération pour le calcul de la DPJ au temps t.
Au final, au temps t, le diamètre le plus probable pour la fibre est pris égal au diamètre théorique DT (i)
tel que la fonction de DPJ est maximale :
d’améliorer significativement la résolution sur les diamètres mesurés. Ceci, en inversant tout ou en partie
les contributions des MDRs. Elles présentent cependant deux inconvénients importants : i) les calculs sont
très lourds (surtout la DPJ), ce qui ne permet d’utiliser ces techniques que pour le post-traitement des
données brutes ; ii) ces méthodes d’inversion nécessitent une connaissance précise de la PDRS théorique
de l’appareil. Elles requièrent donc une connaissance très précise des paramètres optiques de l’appareil et
la prise en compte des effets de la tension de fibrage, du taux de refroidissement... L’analyse d’échantillons
de fibres par MEB (§ 3.4.4) devait nous aider à déterminer plus précisément ces paramètres. Cela n’a pas
été le cas, bien au contraire ! Nous avons donc dû stopper le développement de ces techniques d’inversion
afin de développer un appareil de référence : le diffractomètre haute résolution...
Fig. 2.34 – Comparaison des séries temporelles obtenues avec la méthode d’inversion linéaire, la DPJ, et
la série réelle (l’évolution obtenue avec la DPJ est pratiquement confondue avec l’évolution réelle).
Fig. 2.35 – Ecart-type entre la valeur théorique du diamètre et celle déduite par inversion : avec la méthode
linéaire et la méthode des DPJ en fonction de l’écart-type σφ introduit sur les phases.
revient à déterminer les paramètres optiques pour lesquels un processus de diffusion est dominant et pour
lesquels les oscillations sont en opposition ”de phase” (pour chacun des détecteurs).
Une autre solution possible, pour lisser ces oscillations, est d’éclairer la fibre avec une source de lumière
à spectre plus large. En effet, nous avons vu que la périodicité des oscillations dépend fortement de λ,
on peut donc s’attendre à ce qu’un spectre plus large de longueurs d’onde permette de moyenner ces
oscillations.
Afin de tester cette solution, nous avons calculé la réponse du système de type D pour 1000 longueurs
d’onde dans la gamme : λ ∈ [0.59, 0.69]µm. Pour simuler l’effet de la largeur spectrale de la source nous
avons alors simplement sommé de manière incohérente les différentes PDRS22 . La figure 2.36 illustre les
différents cas considérés : les sources émettent un spectre de forme gaussienne, centré sur la longueur
d’onde d’un laser HeNe λ = 623.8nm, et dont la largeur spectrale à mi-hauteur ∆λF W HM , varie de 0 à
28nm.
La figure 2.37 montre l’évolution des PDRSs correspondantes. On constate que plus la largeur spectrale
de la source augmente plus la PDRS est lissée. Dans le cas limite, ∆λF W HM = 28nm, les oscillations ont
complètement disparu...
La figure 2.38 présente, pour chaque source, l’évolution de l’écart-type entre la PDRS ainsi obtenue
et celle obtenue par filtrage FFT des oscillations : σφF F T .
L’efficacité de cette solution décroı̂t rapidement avec la largeur spectrale de la source. Cependant, par
rapport à une source parfaitement monochromatique, ∆λF W HM = 0, l’écart-type σφF F T décroı̂t d’un
facteur 10 lorsque la largeur spectrale de la source atteint 11.5nm (i.e. 1.8% de la longueur d’onde). La
FFT
résolution sur le diamètre de la fibre évolue dans le même sens : source monochromatique σD ≈ 0.32µm,
source de largeur spectrale ∆λF W HM = 11.5 : σD FFT
≈ 0.032µm. Dans ce dernier cas, cette solution
permet d’augmenter d’un ordre de grandeur la résolution sur le diamètre.
En fait, plus que la largeur spectrale absolue de la source, c’est sa largeur spectrale relative à la
périodicité des oscillations de la PDRS qui est importante. Ainsi, à indice constant (hypothèse forte),
on peut s’attendre à des résultats similaires pour une source centrée sur 380nm et de largeur spectrale
∆λF W HM = 6.9nm, ou bien une source centrée sur 780nm et de largeur ∆λF W HM = 14.2nm...
Fig. 2.37 – Evolution de la PDRS en fonction de la largeur spectrale de la source. Pour les besoins de la
figure, les PDRS ont été décalées artificiellement.
Fig. 2.38 – Evolution de l’écart-type de l’amplitude des oscillations des PDRS par rapport à la PDRS
correspondante filtrée par FFT, en fonction de la largeur spectrale de la source.
d’une largeur spectrale supérieure à quelques nanomètres implique de se tourner vers des sources du type
diodes électroluminescentes, ou lampes spectrales. Il est bien évident que l’emploi de ces sources pose de
gros problèmes : divergence/collimation, puissance, polarisation...
- La longueur de cohérence (i.e. cohérence temporelle) de la source doit être suffisante pour obtenir un
volume de mesure dont les dimensions ne soient pas trop petites (i.e. du fait des mouvements transverses
de la fibre). De plus, cette longueur de cohérence doit être, a priori, supérieure au diamètre de la fibre. La
longueur de cohérence d’une source polychromatique est estimée, classiquement, à partir de la relation
suivante :
λ2
L0 ≈ (2.45)
∆λF W HM
Une simple application numérique montre que, pour ∆λF W HM = 11.5nm et λ = 632.8nm, on obtient
une longueur de cohérence extrémement faible, L0 ≈ 35µm ; ce qui peut rendre très problématique
l’utilisation d’une source à large spectre.
72 Chapitre 2. Modélisation de la réponse de l’interféromètre ”FIBS”
Chapitre 3
Développement de l’interféromètre
FIBS
La figure 3.1 présente une photographie1 du premier interféromètre développé [23]. Cette version
fonctionnait avec des optiques de collection placées en diffusion avant (i.e. configuration de type A §2.4.2).
Ce système, bien que doté d’une très bonne résolution, a été jugé trop encombrant pour mener des études
sur site industriel. C’est pourquoi, nous avons opté pour le développement d’un nouvel interféromètre
(figure 3.2), plus compact, basé sur une configuration de type D. Nous détaillons dans ce qui suit les
caractéristiques techniques de ce système.
Fig. 3.1 – Interférométre de type A lors de son implantation sur le site de Saint-Gobain Recherche, à
Aubervilliers.
3.1 Caractéristiques
3.1.1 Générales
L’interféromètre Phase Doppler FIBS (i.e. FIBRE SIZER) se compose essentiellement de trois éléments :
une tête optique, des composants électroniques d’alimentation et de contrôle, et un ordinateur équipé
d’une carte de numérisation.
La tête optique renferme le système d’émission des faisceaux laser (§ 3.1.2), deux optiques de collection
(§ 3.1.3) et une partie de l’électronique du système (§ 3.1.4), voir la figure 3.2. Le schéma de principe
correspondant est présenté par la figure 3.3.
Elle prend la forme d’un boı̂tier rectangulaire de dimensions L∗P ∗H = 330mm∗330mm∗100mm et d’un
poids d’environ 9kg. Sa face avant (A) est fermée par un hublot (H), en verre borosilicate, transparent
à la longueur d’onde laser et absorbant dans l’infrarouge. Ce hublot est lui-même fixé sur une plaque en
aluminum ajourée qui joue le rôle d’un masque optique. Cet ensemble (hublot+masque) permet de laisser
1 lors de son implantation sur le site d’Aubervilliers de Saint-Gobain Recherche
73
74 Chapitre 3. Développement de l’interféromètre FIBS
passer les faisceaux émis par le système ainsi qu’une partie de la lumière rétrodiffusée par la fibre, tout
en limitant l’entrée du rayonnement de sources parasites comme la lumière naturelle ou le rayonnement
thermique du four. La face arrière du boı̂tier (AR) permet le passage des câbles électriques. Ce boı̂tier est
lui-même fixé sur deux platines de translation et deux platines de contrôle de l’inclinaison dans le plan
horizontal. Ces déplacements permettent de positionner avec précision le volume de mesure sur la fibre à
analyser. Pour des mesures sur site industriel, l’ensemble du système peut être placé sur un trépied photo
muni de réglages de position supplémentaires, voir la figure 3.8.
2λf1
ω1 = (3.1)
πωL
- Le réseau de diffraction par transmission (R) est un disque de verre sur lequel trois réseaux concen-
triques ont été gravés. Ces trois ”pistes” p = 1, 2, 3 ont des pas ( ∆1 = 4µm, ∆2 = 8µm, ∆3 = 16µm)
et un nombre de traits (N1 = 23483, N2 = 12527 et N3 = 6627) différents. Les faisceaux diffractés par le
réseau, autres que les ordres ±1 et 0, sont bloqués (S). L’angle de diffraction des ordres ±1 s’écrit :
λ
θp = ±sin−1 (3.2)
∆p
avec θ1 = 9.10◦ , θ2 = 4.54◦ et θ3 = 2.27◦ .
Le réseau est fixé sur l’arbre d’un moteur électrique (M ). La vitesse de rotation de ce dernier, Rmot ,
permet de contrôler la fréquence de rotation du réseau et, en conséquence, celle du décalage entre les
faisceaux νs (elle se ramène à la fréquence Doppler des signaux dans le cas d’une fibre sans mouvement
suivant l’axe y) :
νD = νs = Np ∗ Rmot (3.3)
3.1. Caractéristiques 75
(A) (H)
L4
L5 L3
(L)
(D)
(B)
(T) (+1)
(0)
(P)
(-1)
L2 (APDs)
(F1) L1
Z (F2)
X
Y (S) (R) (M) (AR)
(BNCs) (//)
Fig. 3.3 – Schéma de principe de la tête de l’interféromètre FIBS.
- Le moteur est lui-même fixé à une platine de translation. Ce dispositif permet de choisir la piste du
réseau utilisée pour la production des faisceaux diffractés.
- Les faisceaux ±1 et 0 sont collimatés par une lentille plan-convexe (L2 ) de focale f2 = 100mm, avec
un diamètre au col, dans le plan image de cette dernière, de : ω2 ≈ 1040µm.
- Trois miroirs, m2 , m3 , m4 sont utilisés pour réaliser un parcours optique de longueur f2 + f3 . Ce
montage ”focale à focale”, bien qu’encombrant, permet d’assurer, avec une bonne précision, que le col
des faisceaux est bien situé dans le volume de mesure final.
- En fin de parcours, juste avant la lentille de transmission L3 , les faisceaux passent à travers une
lame demi-onde (D). Cette lame permet de contrôler l’angle de polarisation des faisceaux (et notamment
les effets de polarisation dus aux miroirs).
- La lentille de transmission L3 , une achromate de focale f3 = 400mm, forme le volume de mesure par
focalisation des trois faisceaux ±1 et 0. Le volume de mesure est formé à 340mm de la face avant de la
tête optique et son diamètre à 1/e2 est de : 2ω0 ≈ 620µm. Ce dernier dépend comme suit des paramètres
optiques précèdents :
2λf3 f1
ω0 = (3.4)
πω0 f2
Après la lentille de transmission, les faisceaux diffractés 0, ±1 ont un demi-angle d’intersection :
θp f 2
αp /2 = (3.5)
f3
avec α1 = 1.147◦, α2 = 0.568◦ et α3 = 0.283◦. Pour les fibres de renfort, nous utilisons exclusivement la
piste p = 1 soit α1 = 1.147◦ .
distribution centrale correspond à l’ordre 0. Sa faible intensité est due au réseau de diffraction qui a été
optimisé par le constructeur pour concentrer 80% de l’énergie initiale dans les faisceaux ±1.
Cette différence d’intensité pourrait être minimisée, voire annulée, par l’emploi d’un réseau de diffrac-
tion optimisé à cet effet. Des réseaux de ce type existent sur le marché, mais pas sous la forme de réseaux
tournants.
Fig. 3.4 – Profil d’intensité des faisceaux laser incidents suivant Y avec X = 0, Z ≈ −30mm, mode 3
faisceaux.
C(a) (b)
Fig. 3.5 – Profil d’intensité du volume de mesure suivant Y, X=Z=0. (a) en mode 2 faisceaux et (b) en
mode 3 faisceaux.
(F2 ). Cet ensemble est lui même fixé sur une platine de micro-translation (P ) à 2 axes (≈ X, Y ). Cette
platine est le seul élément qui permette de compenser les problèmes éventuels de décentrage des optiques
du système. Ce nouvel ensemble est fixé sur un rail optique qui autorise des déplacements le long de
l’axe optique ≈ Z et donc, l’ajustement de la position de la fente avec celle de la tache de focalisation
des deux achromates. La fente optique a une largeur de 200µm et une hauteur de 3mm. Le filtre spatial
ainsi réalisé permet de limiter les dimensions du champ observé par les optiques de collection et donc,
les dimensions effectives du volume de mesure. Ce filtrage permet de renforcer l’hypothèse de ”détecteur
lointain”, chapitre 2.
La fibre optique, F1 ou F2 , est une fibre optique polymère de 1mm de coeur, 2m de longueur et
d’ouverture numérique d’environ N.A. ≈ 0.5. Ces fibres transmettent vers les diodes Photo-
Avalanche (AP D) la lumière collectée par les optiques de collection. Il est à noter que l’utilisation de
fibres optiques simplifie grandement la caractérisation des optiques et les réglages.
0.8
Intensité transmisse
Da Db
0.6
< -5.09° 5.11° > < -5.15° 5.14° >
0.4
a=162.50° b=166.73°
0.2
0.0
-168 -166 -164 -162 -160 -158 162 164 166 168 170 172
[deg] [deg]
Fig. 3.6 – Intensité transmise par les optiques de collection : achromates et fibres optiques.
De ces mesures on peut ainsi déduire :
- La fonction de transmission angulaire T r(ψ), introduite dans le chapitre 2, que l’on utilise pour les
simulations numériques.
- L’angle d’ouverture équivalent de chaque détecteur, 2Ωa = 10.20◦, 2Ωb = 10.29◦
- Et par intégration, l’angle moyen d’élévation de chaque détecteur ψa = +162.50◦, ψb = −166.73◦
La figure 3.7 montre, pour l’optique Db , les profils d’intensité correspondants : (a) suivant X et (b)
suivant Z. On peut ainsi déterminer la dimension latérale du champ observé : Lx ≈ 540µm et la dimension
longitudinale Lz ≈ 2600µm. Ces dimensions sont ”compatibles” 4 avec les dimensions de l’image projetée
du couple (filtre spatial, fibre optique) : sur X ≈ 420µm et sur Z ≈ 2040µm.
Le fait que 2ω0 < Lz permet de conclure que, dans le système FIBS, c’est le diamètre des faisceaux qui
impose la résolution longitudinale (i.e. axe de la fibre) de la mesure. La résolution suivant l’axe optique
est déterminée par les capacités du mode 3-faisceaux cohérents.
C(a) (b)
Fig. 3.7 – Profil d’intensité correspondant au ”champ d’observation” de l’optique Db suivant : (a) X et
(b) Z (axe de la fibre).
3.1.4 Electronique
Composants
- Le laser est alimenté en haute tension par son alimentation d’origine, placée à l’extérieur de la tête
optique.
- Le réseau de diffraction est mis en mouvement par un mini moteur électrique à courant continu,
sans broches, équipé d’un réducteur planétaire sans jeu et d’une roue codeuse. Sa tension d’alimentation,
comprise entre 0.5-12V, permet de faire varier sa vitesse de rotation et celle du réseau dans la plage
≈ 0.5 − 5tr/s. Le contrôle de la tension d’alimentation du moteur est assuré par un boı̂tier de commande
et d’alimentation (BCA) externe à la tête optique. Ce boitier est relié à la tête obtique par un connecteur
parallèle (//), voir les figures 3.3- 3.8
- Les diodes photoavalanches (APDs) sont équipées d’un étage de préamplification, de bande passante :
[0.1 − 10M hz]. Elles sont placées à l’intérieur de la tête optique. Leur alimentation électrique est controlée
à partir du boı̂tier (BCA). Le gain des deux APDs est réglable manuellement. Chaque APD restitue
un signal de tension proportionnel au signal optique. Ce signal électrique sort de la tête optique par
l’intermédiaire de câbles BNC blindés, (BN Cs) .
- Le boı̂tier de commande et d’alimentation (BCA) a été réalisé au laboratoire. Il assure l’alimentation
des APDs et le contrôle de la vitesse de rotation du moteur du réseau. Ce boı̂tier est commandé via le
port parallèle du PC.
Carte de numérisation et PC
- Les signaux de tension issus des APDs sont envoyés vers un boı̂tier de connexion et de synchroni-
sation (BCS). Ce boı̂tier permet le conditionnement les signaux Doppler avant leur numérisation par
une carte PCI implantée dans un PC. Ce boı̂tier permet également de recueillir des signaux logiques de
synchronisation (i.e. TTL, externes/internes).
4 Dimensions
physiques de la fente (200µm × 3mm) qui sont a multiplier par le grandissement de l’optique de collection
(×2) et un facteur de projection dû à l’angle d’élévation (1/cosψ2 ≈ 1.04)
3.2. Principe du traitement des signaux Doppler 79
Fig. 3.8 – Les principaux éléments de l’interféromètre Fibs : a) Alimentation et contrôle du réseau, b)
PC et carte de numérisation, c) Alimentation du laser, d) Tête optique montée sur trépied, e) Boı̂tier de
connexion et synchronisation et f ) Tête optique implantée sur la filère mono-téton.
- La carte de numérisation, 12 bits, a un taux d’échantillonnage maximale de 2.5 M/s pour deux voies,
1.25 M/s pour quatre voies...
- L’ordinateur utilisé est un PC, blindé, équipé de 128M o de mémoire et d’un processeur type P III
cadencé à 350M Hz.
Calcul du déphasage :
L’extraction de l’évolution temporelle du déphasage est réalisée à l’aide du calcul de la fonction inter-
spectre des deux signaux, Gab :
Gab (ν) = Cab (ν) + jQab (ν) (3.6)
80 Chapitre 3. Développement de l’interféromètre FIBS
où Cab (ν) correspond au spectre de coı̈ncidence des deux signaux et où Qab (ν) correspond à leur spectre
en quadrature. En notation exponentielle nous avons :
où Gab est le module de Gab (ν), et θab (ν) son angle de phase, avec :
1
2
Gab (ν) = (Cab (ν) + Q2ab (ν)) 2 (3.8)
Qab (ν)
θab (ν) = tg −1 (3.9)
Cab (ν)
Dans le spectre en fréquence du module de Gab la fréquence Doppler commune aux deux signaux est
identifiée par un pic caractéristique de forme ”gaussienne”. A partir de là et avec l’équation (3.9), on
peut déterminer le déphasage entre les signaux. Dans le cas du mode trois faisceaux on observe deux pics
caractéristiques et l’on détermine les deux déphasages correspondants. En pratique, les signaux Doppler
sont numérisés, Ia (n∆t) et Ib (n∆t) :
2∆t ∗
Gab (νk ) = [i (νk ).ib (νk )] (3.10)
N a
k
νk = (3.11)
N ∆t
où,
– n = 0, 1, 2, . . . , N − 1 et k = 0, 1, 2, . . . , N − 1
– ia (νk ) et ib (νk ) sont les transformées discrètes des deux séries, Ia (n∆t) et Ib (n∆t).
– N correspond au nombre d’échantillons compris dans la fenêtre d’analyse.
– ∆t est égal à l’inverse de la fréquence d’échantillonnage des signaux.
Fig. 3.9 – Principe de l’extraction de l’évolution temporelle du déphasage avec la fonction interspectre à
fenêtre glissante.
Analyse temporelle :
- L’évolution temporelle du déphasage est déduite du calcul de la fonction interspectre : sur une fenêtre
d’analyse de N échantillons, que l’on fait glisser généralement de N/2 échantillons à chaque pas d’analyse
(résolution temporelle ≈ N ∆t/2), voir la figure 3.9.
- La résolution sur la mesure du déphasage dépend principalement de la largeur temporelle des fenêtres
d’analyse (soit de N, à fréquence d’échantillonnage constante) [66]. La résolution temporelle/spatiale de
la mesure décroı̂t avec la largeur des fenêtres d’analyse.
- Une Transformée de Fourier Rapide est utilisée pour le calcul des transformées discrètes ia (νk ) et
ib (νk ). N doit être une puissance de 2.
- La résolution en fréquence est améliorée par l’utilisation d’un fittage des pics caractéristiques par
une fonction gaussienne.
- En mode 3 faisceaux, on extrait les déphasages correspondant aux deux fréquences Doppler.
Acquisition
Le système dispose de trois modes d’acquisition des signaux Doppler (signaux de tension issus des
APDs) :
- Le mode Temps Réel réalise : i) l’acquisition simultanée de signaux sur une durée réduite à une seule
fenêtre d’analyse, ii) le traitement de ces deux signaux et iii) l’affichage de la mesure correspondante. Ce
mode permet de suivre en temps réel l’évolution du diamètre de la fibre sans conduire à la saturation des
disques. Le nombre de mesures réalisées est de l’ordre de 80 à 200 mesures/s, suivant le type de données
et les statistiques affichées. La résolution temporelle/spatiale de ce mode d’acquisition est cependant
limitée (100 à 250mm pour une vitesse de défilement de Vf = 20m/s). Il est à noter que ce mode est très
utile pour la phase de réglage de l’interféromètre.
- Le mode Continu DMA permet de remplir la mémoire du PC jusqu’à ”saturation”. Les données sont
ensuite post-traitées. Ce mode permet d’atteindre une très haute résolution temporelle/spatiale. Dans ce
mode, la carte de numérisation permet d’atteindre un taux d’échantillonnage maximal de 2.5M ech/s par
voie. En pratique, ce mode d’acquisition permet d’atteindre une résolution de 1.3mm sur une fibre défilant
à 50m/s, c’est à dire ≈ 50000 mesures/s. Le nombre d’échantillons que l’on peut acquérir est cependant
limité par la mémoire du PC à 2 × 21M ech. A titre indicatif, le temps nécessaire au post-traitement de
toutes ces données est inférieur à 16s.
- Le mode Continu sur Disque gère l’acquisition en continue des données sur le disque dur du PC.
Les données sont ensuite post-traitées. Dans ce mode, le taux de transfert des données vers le disque est
limité à 2 × 0.25M ech/s. En pratique, ce mode d’acquisition permet d’atteindre une résolution de 5mm
sur une fibre défilant à 20m/s et cela, sur plusieurs dizaines de kilomètres.
Chacun de ces modes peut-être déclenché manuellement (à partir du PC) ou à partir d’un signal de
synchronisation externe (TTL) relié au boı̂tier (BCS). En début et en fin d’acquisition le système peut
également produire un signal de synchronisation sortant. Ces fonctionnalités permettent de synchroniser
le système avec un événement particulier ou d’autres appareils de mesure...
Le système peut également numériser de manière synchrone deux signaux analogiques supplémentaires :
signal de vitesse du bobinoire, mesure de la température du bulbe, etc.
Les paramètres entrant en compte dans la qualité du traitement des signaux Doppler sont très nom-
breux. Afin de simplifier le réglage du système nous avons développé deux algorithmes d’optimisation
de ces paramètres en fonction des paramètres du procédé de fibrage. Ainsi, en fonction de la résolution
spatiale souhaitée, ou du nombre de mesures par seconde et de la vitesse de bobinage, le système ajuste
automatiquement : la fréquence de rotation du réseau tournant, la fréquence d’échantillonnage et la
banque passante des filtres, la largeur et le recouvrement des fenêtres d’analyse, le gain de la carte de
numérisation, etc.
Affichage et sauvegardes
Dans sa version actuelle le logiciel dispose de 5 écrans d’affichage principaux :
- Acquisition des données : sur cet écran on visualise les signaux Doppler, l’évolution temporelle du
diamètre de la fibre, l’histogramme correspondant, ainsi que différentes statistiques et aides aux réglages,
et notamment, le positionnement de la fibre par le mode 3 faisceaux, voir la figure 3.10.
82 Chapitre 3. Développement de l’interféromètre FIBS
- Séries temporelles : sur cet écran on visualise l’évolution temporelle du diamètre de la fibre, de la
fréquence Doppler des signaux (mouvements transverses de la fibre ou instabilités du réseau), du rapport
signal sur bruit, de la visibilité ou de l’inclinaison de la fibre. Ces trois dernières quantités sont des aides
au contrôle du bon fonctionnement du système.
- Paramètres du système : cet écran rassemble tous les paramètres de réglages non optiques du
système. Il regroupe également toutes les fonctionnalités de sauvegarde ou d’importation des données
et paramètres, voir la figure 3.10.
- Calcul de Diffusion : cet écran permet de calculer la réponse théorique de l’appareil (PDRS), en
fonction des différents paramètres optiques du système et des propriétés de la fibre.
- Mesure de la Tension : cet écran gère l’acquisition et le traitement des signaux Doppler pour la
mesure en ligne de la tension de fibrage.
3.4 Validations
Valider un nouveau système de mesure est une tache particulièrement ardue lorsque l’on ne dispose
pas de mesures de référence.
Dans ce qui suit, nous détaillons le principe et comparons les résultats obtenus avec différentes tech-
niques ”classiques” de mesure du diamètre des fibres. Ces techniques se sont révélées peu fiables et peu
précises.
La valeur théorique du débit de verre peut-être estimée en faisant un bilan hydrodynamique [41] :
πR04 ρ(∆p + ρgL0 ) h
Q= = Q 0
3πR0 µ 0
8µ0 L0 1 + 8L0
Fig. 3.10 – Copie de deux écrans du logiciel FIBS : en haut, l’écran pour l’acquisition et en bas, l’écran
des paramètres.
84 Chapitre 3. Développement de l’interféromètre FIBS
Fig. 3.11 – Copie de l’écran du logiciel FIBS pour le calcul de la réponse théorique de l’interféromètre à
l’aide du modèle électromagnétique.
Au total, on estime que pour une expérience de durée t[mn], l’évolution relative du diamètre de la
fibre est de :
∆D
≈ t[s]1.10−5
/s + 0.8% ≈ t[mn]0.06% + 0.8% (3.16)
D
Exemple : entre le début et la fin d’une expérience, d’une durée t = 10mn, le diamètre moyen d’une
fibre de D ≈ 20µm diminue de ≈ 0.3µm si le creuset n’est pas réapprovisionné en verre (voir le § 7.1).
Fig. 3.12 – Evolution de la tirée de la filière avec la température selon différentes méthodes de calcul.
La masse M de la bobine est connue à une incertitude ∆M = 5.10−5 kg, sachant que M = Qt =
1.10−5t kg/s. La densité ρ est de 2600kg/m3 avec une incertitude, d’après Gupta, de ∆ρ = 30kg/m3 (
d’après nos mesures elle est plutôt de l’ordre de ∆ρ = 100kg/m3). L’influence de la variation du rayon
du bobinoire Rc a déjà été prise en compte. L’incertitude sur la fréquence de rotation, Fc est estimée à
1%. Celle sur le temps, t, est de l’ordre de 1s. Ce qui donne au final :
∆D 300% 5%
≈ 2.5% + ≈ 2.5% + (3.19)
D t[s] t[min]
La figure 3.13 présente l’évolution de cette erreur relative en fonction du temps. Il existe un durée
d’expérience optimale pour minimiser l’incertitude sur le diamètre : t = 6 − 8mn pour une incertitude de
≈ 4.5%.
8
Dérive du diamètre au cours du temps
Erreur estimée pour la mesure de pesée
6 Somme des deux erreurs
]%
[d 4
/d
∆
0
0 2 4 6 8 10 12 14
Temps [mn]
Fig. 3.13 – Evolution temporelle de l’incertitude sur le diamètre déduit par pesée.
86 Chapitre 3. Développement de l’interféromètre FIBS
∆D b∆T
= (3.20)
D 2.(T − Tf ig )2
Une autre façon d’estimer cette incertitude est de déterminer les diamètres minimums et maximums de
la fibre pour différentes conditions de fibrage, voir la figure 3.14. L’incertitude relative sur le diamètre
diminue avec la température de fibrage. Les deux méthodes donnent des résultats légèrement différents,
par la suite, on utilisera les résultats obtenus par la seconde méthode (maximums et minimums).
La figure 3.16 compare les résultats obtenus avec l’interféromètre, la pesée, l’estimation hydrodyna-
mique et la mesure filamentaire (voir plus loin). Les mesures sont plus ou moins en accord. L’estimation
hydrodynamique semble sous-estimer le diamètre de la fibre formé. La mesure filamentaire semble le
sur-estimer.
11,5
] ∆D/D = 2.8+22exp(-T/200)
%
[
11,0
é 2
∆D/D = 3538*10/(2(T-766) )
im
10,5
ts
ee 10,0
rt
è 9,5
m
ai 9,0
de
lr 8,5
us
ed
8,0
uti 7,5
tr
ec 7,0
nI 6,5
Fig. 3.15 – Comparaison pour différentes températures des mesures obtenues par interféromètrie laser et
par estimation hydrodynamique.
Fig. 3.16 – Comparaison des mesures obtenues par interféromètrie laser, pesée, estimation hydrodyna-
mique et microscopie optique (mesure filamentaire).
échantillon de fibre on obtient des variations importantes du diamètre (0.1µm à 1µm). Nous les attribuons
à la mesure par MEB et non pas à des variations du diamètre de la fibre !
En fait, nous avons été très déçu des résultats obtenus par MEB : les résultats se font attendre
plusieurs semaines à plusieurs mois... et semblent peu fiables7 .
Cette technique permet cependant de visualiser l’état de surface des fibres8 . La figure 3.20 présentes
quatre photographies de fibres atypiques. a) la surface de cette fibre présente des variations de contraste
importantes. Nous n’avons pas d’explication très claire à ce dernier phénomène (métallisation non ho-
mogène, un début de corrosion due à l’humidité de l’air ?). b) la surface de cette fibre est très irrégulière
7 Outre les problèmes de délai, de représentativité des échantillons... il faut savoir que, le facteur de grandissement des
et semble contaminée par quelques particules (i.e. cette fibre était stockée dans du papier kraft depuis
plus de deux ans)9 . c) et d) le diamètre de cette fibre semble varier brusquement. Sur une centaine
d’échantillons analysés, c’est la seule fois où nous avons observé ce phénomène. Il ne nous est pas possible
de dire s’il s’agit d’une variation réelle du diamètre de la fibre ou un effet optique du à la remontée de
l’adhésif carbone 10 .
La figure 3.18 compare les résultats obtenus avec cette méthode, l’interféromètre et le diffractomètre11
pour différents échantillons de fibres coupées. L’accord global est bon. On constate cependant que pour
certains diamètres la réponse de l’interféromètre semble davantage ”fluctuer” (i.e. la procédure actuelle
d’inversion des phases de l’interféromètre néglige les contributions des MDRs).
La figure 3.19 compare les résultats obtenus par MEB, sur des échantillons de fibres coupées, avec
la réponse moyenne (et l’écart-type) de l’interféromètre ontenue sur plusieurs centaines de mètres de
fibre, l’estimation hydrodynamique. L’accord FIBS-MEB est à mieux que 5% sur toute la gamme. Les
contributions des MDRs sont ici ”lissées” (du fait de la largeur de la distribution) la moyenne obtenue avec
l’interféromètre fluctue donc moins. Le diamètre prédit avec l’estimation hydrodynamique est relativement
correct, excepté quelques écarts. Cet ecart est surtout dû aux bulles.
Fig. 3.18 – Comparaison des mesures obtenues par interféromètrie, diffractométrie et microscopie
électronique pour des échantillons de fibres coupées.
Fig. 3.19 – Comparaison des mesures obtenues avec l’interféromètre et le modèle hydrodynamique, sur
des fibres en cours de production (t = 4 min), et celles obtenues par microscopie électronique sur des
échantillons de fibres coupées.
3.4.7 Conclusion
Les prédictions des techniques classiques : pesée, estimation hydrodynamique, microscopie optique
et électronique ont pu être comparées entre-elles et avec celles de l’interféromètre et du diffractomètre.
Les techniques classiques se sont révélées lourdes d’emploi et peu fiables comparées à celles que nous
avons développé. Cependant, nous avons validé la réponse de l’interférométre FIBS et du diffractomètre
HIREDI, que ce soit pour une fibre coupée ou en cours de production.
Globalement, nous dirons que l’incertitude de mesure sur le diamètre est de l’ordre de : ≈ 8% avec
l’estimation hydrodynamique et ≈ 5% avec la pesée. La résolution de l’interféromètre est conforme à ce
que nous avions prédit numériquement.
Fig. 3.20 – Microscopie électronique : photographies de fibres atypiques... a) la surface de cette fibre
présente des variations de contraste importantes ; b) la surface de cette fibre, stockée depuis plus de deux
ans, est très irrégulière ; c) et d) Fibre dont le diamètre semble varier brusquement.
Fig. 3.21 – Microscopie optique : image de la coupe polie d’une nappe de fibres.
Chapitre 4
4.1 Introduction
Nous avons vu dans le chapitre précédent que les méthodes classiques de mesure du diamètre des fibres
se sont révélées trop peu précises pour valider la réponse de l’interféromètre FIBS. Nous avons donc dû
développer un autre système de mesure de référence doté d’une bonne résolution spatiale et temporelle.
Ce système, nécessairement optique, est un diffractomètre laser à haute résolution (HIgh REsolution
DIffractometer, HIREDI) [76].
On trouve dans la littérature de nombreux travaux sur la granulométrie par diffraction de particules
métalliques, absorbantes, ou d’un ensemble de particules transparentes ou non (suspension, pulvérisa-
tion, . . .). Ces travaux contournent en fait un problème difficile : la prise en compte des Résonances
Morphologico-Dépendantes (MDRs). En effet, les MDRs n’existent pas, ou sont très fortement atténuées,
pour de telles particules (voir le cas des fibres métalliques, figure 1.7) ou bien elles sont lissées, moyennées
(distribution granulométrique). Nombre de ces travaux utilisent la théorie de Fraunhofer pour inverser
les diagrammes de diffraction. Souvent, parce qu’elle est simple et rapide d’emploi, et dans certains cas,
parce qu’elle permet de prédire relativement simplement la diffusion par des particules de forme complexe
(particules irrégulières, cristaux,. . .). Au final la résolution obtenue avec ces systèmes, en laboratoire, sur
des fibres de verre, dans la gamme D = 5 − 30µm, n’excède pas 2 à 3 µm [54, 49, 84].
Dans ce qui suit, notre choix a donc été d’utiliser la théorie de Lorenz-Mie (avec prise en compte
des effets de biréfringence, taux de refroidissement. . .) pour calculer les diagrammes théoriques des fibres
91
92 Chapitre 4. Diffractomètre Haute Résolution HIREDI
à mesurer et, d’inverser les MDRs en définissant un estimateur de corrélation entre les diagrammes
théoriques et expérimentaux. Dans un souci de concision, nous mettrons principalement l’accent sur les
originalités de notre travail, sans revenir sur le principe bien connu de la diffractométrie.
Y
(13)
Ff (10)
(3) (2) (1)
Z
X
O
- (10) Le point de focalisation du faisceau direct. Il est localisé à quelques centaines de micromètres
de la zone sensible du CCD. Ceci permet d’éviter la saturation du capteur. En phase de réglage, après
diminution de l’efficacité du couplage (2), la tâche de focalisation du faisceau direct peut être observée
en déplaçant la caméra suivant Y à l’aide de (13).
- (11) Une caméra linéaire numérique. La fréquence de ligne (i.e. d’acquisition) de cette caméra 8 bits
peut être ajustée dans la gamme : 70-43000 Hz. Elle dispose de 1024 pixels carrés (∆y = ∆z = 10 µm).
Le signal numérique issu de cette caméra est transféré dans la mémoire d’un PC via par une carte vidéo
numérique au format LDVS.
- (12) La ligne CCD. Elle se situe dans le plan image de la lentille de Fourier (9).
- (13) Deux micro-déplacements. Ils permettent d’ajuster la position de la caméra dans le plan (YOZ).
La figure 4.3 présente un enregistrement des figures de diffraction successives, enregistrées à une
cadence de 1.024 kHz, lorsque la fibre est soumise à une rampe de vitesse croissante. On remarquera
la complexité des évolutions, de même que la présence d’un bruit périodique (réflexions parasites des
optiques) 2 . Ces dernières remarques justifient, en partie, pourquoi l’intensité des lobes de diffraction,
même relative, n’est pas utilisée pour l’inversion.
Le passage de bulles (ou fibres creuses) se caractérise par une brusque variation de l’intensité diffusée,
figure 4.3 a. Ces passages sont aisément détectables, à l’oeil, ou par un logiciel de traitement d’images,
en recherchant les discontinuités.
Toute évolution du diamètre de la fibre est également aisément repérable. La figure 4.3 b met l’accent
sur une évolution ”rapide” et la figure 4.3 c), sur une évolution plus ”lente” 3 . La figure 4.9 présente le
résultat de l’analyse de ces diagrammes.
En principe, pour un diamètre donné, le nombre des extremums trouvés pour un diagramme expérimental
+
devrait être égal à celui obtenu pour le diagramme théorique correspondant, i.e. kE (n) = kT+ (i) et
− −
kE (n) = kT (i). En pratique, en raison du bruit optique et électronique, et du fait du peu de contraste de
+
certains extremums, ces nombres sont souvent différents même s’ils restent proches : kE (n) ≈ kT+ (i) , kE
−
(n) ≈
kT− (i).
S’ils sont très différents, c’est que le diagramme associé au diamètre théorique DT (i) est très différent
du diagramme expérimental mesuré pour le diamètre DE (n), i.e. DT (i) = DE (n). On introduit alors un
coefficient de discrimination δk (i) tel que :
" + −
0, pour kE (n) − kT+ (i) > σk (i) ou
k (n) − k − (i) > σk (i)
E T
δk (i, n) = (4.3)
1, pour kE +
(n) − kT+ (i) σk (i) et k − (n) − k − (i) σk (i)
E T
où σk (i) est un paramètre de contrôle qui quantifie l’écart maximum toléré, entre le nombre d’extré-
mums d’un diagramme théorique et d’un diagramme expérimental, pour que ceux-ci correspondent à une
fibre ayant les mêmes caractéristiques.
Pour une large plage de variation du diamètre des fibres, le nombre d’extrémums correspondant varie
beaucoup (i.e. 0 kT+ 9, pour 5 D 35µm). Ce nombre est plus important pour les ”grosses”
fibres, ceci peu générer un biais de la procédure d’inversion, au détriment des fibres les plus petites. Afin
d’éviter ce biais, on introduit la fonction de pondération ou de «normalisation» suivante :
δk (i, n)
∆K (i, n) = (4.4)
kT+ (i) kT− (i) + 1
Cette dernière est maximale pour les petites fibres et décroı̂t avec leur taille. Elle est nulle si le nombre
d’extrémums trouvés pour le diagramme expérimental n, diffère trop de celui du diagramme théorique, i.
L’écart angulaire absolu entre la position théorique et la position expérimentale de chacun des
extrémums s’écrit : + +
∆θ+ i, n, jT+ , jE
+
= θT i, jT − θE+
n, jE+
(4.5)
∆θ− i, n, jT− , jE
−
= θT− i, jT− − θE
− −
n, jE
4 assez lourds
4.3. Méthode d’inversion des diagrammes de diffraction 95
q=1-11±
Zoom s
Fig. 4.4 – Evolution de la position angulaire des extrémums, θT+ (i), θT− (i), en fonction du diamètre de la
fibre (i.e. analyse de la figure 1.7 gauche).
Pour chaque couple d’extrema (Expérience/Théorie) on peut alors associer un estimateur de ”corrélation” :
σθ (i)
∆P + i, n, jT+ , jE
+
= ∆θ + (i,n,jT+ +
,jE )+1
−
− −
σθ (i) (4.6)
∆P i, n, jT , jE = ∆θ − (i,n,jT− −
,jE )+1
où σθ (i) 1 est un second paramètre de contrôle. Il représente l’écart angulaire maximum toléré,
entre la position d’un extremum théorique et d’un extremum expérimental, pour que ceux-ci puissent être
considérés comme identiques. ∆P + i, n, jT+ , jE +
et ∆P − i, n, jT− , jE
−
décroissent linéairement lorsque
l’écart angulaire croı̂t et, il sont inférieurs à 1 si l’écart angulaire dépasse la valeur du paramètre de
contrôle correspondant.
Le diamètre de la fibre, DE (n) , n = 1, 2 · · · nE , est alors posé égal au diamètre théorique donnant un
estimateur de corrélation maximum :
Afin de réduire les temps de calcul, si l’on connaı̂t a priori la plage de variation du diamètre de la
fibre, Dmin D Dmax , il est souhaitable de limiter le nombre de P (i, n) à l’intervalle :
P (i = imin , imin + 1 · · · imax , n) où DT min = D (imin ) et DT max = D (imax ).
La figure 4.5 présente une évolution typique de P (i, n), obtenue pour une fibre de diamètre DE ≈
DT = 23.38µm. L’estimateur P (i, n) est infiniment petit en dehors de l’intervalle considéré. De plus,
on constate que l’amplitude du maximum est très supérieure à celle des maximums locaux. Ces derniers
correspondent à des diamètres de fibres, DT = 22.39 et DT = 24.32µm, qui donnent des diagrammes de
diffraction très «semblables» à celui d’une fibre de DT = 23.38µm (périodicités).
1
Diam ètreleplusprobable
23.38µm
facteur:80
0.1
P(i,n)
Spectressim ilaires...
22.39µm
24.32µm
0.01
1E-3
22.25 22.50 22.75 23.00
23.25 23.50 23.75 24.00 24.25 24.50
Diametre[µm]
Fig. 4.5 – Exemple d’évolution de l’estimateur P (i, n).
cherche des extremums requière également de définir la largeur (et le type) de la fenêtre de lissage des
diagrammes5 , la distance minimale entre deux extremums...
Au terme d’une analyse paramétrique nous avons déterminé, une fois pour toute, les valeurs optimales
pour ces paramètres. Pour les fibres de verre-E et le système HIREDI, dans sa version actuelle, elles sont
données par la figure 4.7. Ces valeurs sont données en fonction du nombre de maximums du diagramme
correspondant (et non pas du diamètre de la fibre, que l’on cherche à déterminer...), avec : Filtering :
la largeur totale de la fenêtre de lissage en centièmes de degré, Width : la distance angulaire minimum
entre deux extremums consécutifs, Neighb. correspond à σk et Sigma correspond à σθ .
En mode auto adaptatif, le logiciel HIREDI adapte automatiquement ses filtres et les autres pa-
ramètres d’analyse, en consultant cette table, pour chaque diagramme mesuré.
Fig. 4.6 – Ecran principal du logiciel HIREDI. On voit une partie des diagrammes en cours de traitement,
un diagramme brut et sa version filtrée, diverses données sur le nombre d’extrémums détectés, le diamètre
mesuré...
100 Chapitre 4. Diffractomètre Haute Résolution HIREDI
l’ordre de σD ≈ 0.02µm, contre σD ≈ 0.3−0.4µm pour l’interféromètre (suivant la gamme des diamètres).
Cette différence vient de ce que le diffractomètre inverse les contributions des MDRs, ce qui n’est pas le
cas de l’interféromètre.
- Les diagrammes obtenus sont très... ”visuels”. Les perturbations, tout comme les bulles, sont aisément
détectables, même à l’oeil nu, alors que les signaux Doppler de l’interféromètre ne sont pas interprétables
visuellement.
Avantages : FIBS/HIREDI
- L’interféromètre est beaucoup plus compact et dispose de tous les réglages de position et d’inclinaison.
Il ne nécessite qu’un seul accès optique. La mesure est effectuée à plus grande distance de la fibre (340 mm
au lieu de moins de 50mm pour le diffractomètre). Il est donc mieux adapté à des mesures sur site
industriel.
- L’interféromètre permet l’enregistrement de séries expérimentales beaucoup plus longues et avec une
meilleure résolution temporelle. Le traitement des données brutes est beaucoup plus rapide (traitement
de signaux 1D au lieu de signaux 2D). Il est donc mieux adapté aux études statistiques, sur de grands
échantillons.
Perspectives de développement du système HIREDI
Le système HIREDI a été conçu pour réaliser des expériences de validation, sur une filière de la-
boratoire. Sa configuration opto-mécanique actuelle n’est pas adaptée à la réalisation de campagnes de
mesures sur une filière industrielle. Cependant, nous disposons à présent de tous les outils nécessaires
pour envisager une miniaturisation de ce système.
La procédure d’inversion, originale, fonctionne parfaitement. Si bien que d’autres diagnostics sont
envisageables avec ce système :
- Mesure de l’indice de réfraction des fibres. Des tests préliminaires ont montré que la résolution du
système est telle que l’on pouvait déduire l’indice m en plus du diamètre D (en entrant dans le calcul des
P (i, n) une inconnue supplémentaire : l’indice)...
- La caractérisation des fibres creuses ou bulles est envisageable si elles sont axi-symétriques (en entrant
dans le calcul des P (i, n) une inconnue supplémentaire : le diamètre interne)... Les tests préliminaires
n’ont pas été concluants.
- La mesure d’autres gammes de tailles (en changeant la lentille de Fourier) ou d’autres types de fibres
(polymère, carbone, métallique,...) ne pose a priori pas de difficultés.
15
Diffractometre
3.1
Interferometre
14
13
Interferometre (i.e. MDRs)
3.0
12
+
11
Bulle
2.9
10
9
2.8
Bulles
Temps [s]
8
27.5
27.0
26.5
26.0
25.5
7
6
Zoom n˚2
Zoom n˚1
5
Filtrée
Brute
4
5.2
Bulle
3
Diffractometre seul:
2
5.0
1
0
9
30
27
24
21
18
15
12
Diametre [µm]
Fig. 4.9 – Comparaison des résultats obtenus avec l’interféromètre FIBS et le diffractomètre HIREDI
lorsque le procédé est soumis à une rampe de la vitesse de fibrage.
Deuxième partie
Etude de l’hydrodynamique du
procédé
103
Chapitre 5
Depuis les années 60, de nombreuses études ont été publiées sur la modélisation du fibrage, du verre
[65, 32, 22, 15, 85, 112, 43, 42, 91] ou des polymères [39, 46, 47, 48, 12]. On peut trouver une synthèse
d’une partie de ces travaux dans les ouvrages suivants : [11, 114].
Vers la fin des années 80, Gupta [11] a proposé un modèle hydrodynamique pour le fibrage des fibres
de renfort, qui fait encore référence. Ce modèle, non-isotherme et instationnaire, prend en compte de
nombreux effets : la forte thermo-dépendance de la viscosité du verre, les transferts radiatifs, la trainée
aérodynamique de la fibre, les effets de la tension superficielle... Cependant, comme c’est le cas dans la
plupart des travaux sur le fibrage, ce modèle suppose l’axi-symétrie du jet de verre et que ce dernier peut
être considéré comme étant 1D (voir le chapitre suivant).
Le taux de refroidissement :
Lors du fibrage du verre, à haute température, les transferts thermiques par convection et par rayon-
nement sont très importants2 . Le taux de refroidissement des fibres de renfort peut en effet atteindre
−100 000◦ /s. Cependant, les études expérimentales [53] montrent qu’augmenter le taux de refroidisse-
ment des fibres (i.e. échanges thermiques entre l’air et la fibre) stabilise le procédé.
D’après Glicksman [34, 33], les transferts radiatifs sont dominants dans la partie supérieure du jet
z < R0 , et ils deviennent rapidement négligeables devant les transferts convectifs, en aval du jet z > 3R0 .
Cet auteur a également introduit un nombre de Nusselt qui rend compte de la stratification en température
(gradient) de l’air ambiant près de la surface du jet.
Pour les transferts radiatifs, les modèles utilisent pratiquement toujours un coefficient d’échange qui
prend la forme non linéaire de la loi Stephan-Boltzman (voir par exemple Gupta [11]).
Cependant, plus récemment, Song. et al.[98] ont proposé d’utiliser un modèle radiatif à deux dimensions
où la trajectoire des photons ”thermiques” est prise en compte par une méthode des ordonnées discrètes.
Ce modèle introduit une dépendance de la distance parcourue par les photons avec l’absorptivité spectrale
du verre. Ceci permet de simuler les effets ”optiques” induits par la courbure du profil du jet de verre
sur la propagation du rayonnement thermique. Le tableau 5.1 donne quelques ordres de grandeurs des
distances qui sont ainsi parcourues.
Tab. 5.1 – Evolution de la distance parcourue Λ (à 1/e2 ) par les photons en fonction de la partie imagi-
naire de l’indice du verre E (l’absorptivité : K = 2πk/λ) pour différentes longueurs d’onde. La température
correspondante du corps noir (CN) est donnée à titre indicatif.
Effets capillaires :
Donnelly et al [25] ont étudié l’influence de la tension superficielle sur la stabilité du procédé. Cette
étude a été développée par Geyling et al[32], avec la prise en compte des effets visqueux. Entov et Yarin [28]
2 Dansle cas du fibrage des polymères, les transferts radiatifs sont négligeables, mais les réactions de polymérisation son
exothermiques...
5.2. Les différentes méthodes de résolution 107
ont également étudié l’influence des forces de tension superficielle en déterminant un taux de croissance
pour les petites perturbations. Pour les fibres de renfort, les études concluent généralement que les effets
capillaires sont négligeables par rapport aux effets inertiels3 .
Bechtel et al. ont étudié, théoriquement et numériquement, les effets de torsion sur la stabilité du procédé
[7]. Effets que nous avons nous-mêmes observés expérimentalement [74].
3 pourla filière mono-téton, T0 = 1140 − 1250◦ C, nous avons (Re/W e)z=0 < 0.28 avec W e(z = 0) = 2R0 V02 ρ/σ
4 Fibresoptiques.
5 Equations de Bilan (Eq.), Conditions Limites (C.L.)
108 Chapitre 5. Bibliographie sur la modélisation du fibrage
La méthode analytique
Eq. + C.L. −→ Simplication −→ Résolution analytique
Cette méthode nécessite l’identification des différentes contributions mécaniques et énergétiques pour
pouvoir simplifier les équations de continuité. Elle permet donc d’avoir une bonne compréhension des
phénomènes physiques.
Elle a été utilisée dans le cas du fibrage isotherme, puis non-isotherme, ainsi que pour étudier la stabilité
du procédé soumis à de petites perturbations du rayon [97]. radiale de
Tab. 5.2 – Ordre de grandeur du paramètre de stabilité de Shah et Pearson pour différentes températures
de fibrage.
Fig. 5.2 – Positions de nos points de mesure dans le diagramme de stabilité théorique de Shah et Pearson.
Modélisation physique de
l’écoulement
Dans ce chapitre, nous détaillons les hypothèses et le formalisme du modèle ”physique” que nous
avons développé pour simuler l’hydrodynamique de l’écoulement du jet de verre en sortie du téton (avant
la zone de figeage). Les prédictions de ce modèle sont ensuite comparées à celles obtenues avec d’autres
modèles hydrodynamiques et des solutions analytiques, ainsi que des données expérimentales.
Fig. 6.1 – Schéma de principe de l’écoulement du verre et des repères utilisés pour sa modélisation.
La figure 7.26 présente des images du jet de verre produit par une filière mono-téton, pour différentes
températures de fibrage. On oberseve que le jet de verre est axi-symétrique et que, sous l’action de la
tension imposée par la vitesse du bobinoir, il se contracte fortement pour former, en aval, une fibre de
verre ∼ 1D dont le diamètre varie de D = 5 à 40µm.
112 Chapitre 6. Modélisation physique de l’écoulement
Le transfert de chaleur q est plus complexe à définir. On peut le décomposer en la somme de deux
termes, un terme radiatif q r , et un terme conductif q c , soit q = q c + q r .
– en r = R(t, z), à l’interface, certains travaux négligent la friction de l’air et la pression atmosphérique.
Seule la force de capillarité [57] est prise en compte :
1 1
(T.n).s = 0 ; (T.n).n = σ +
Rzz Rnθ
– en t = 0, on s’intéresse aux solutions stationnaires du problème :
R(0, z) = Rt0 (z) ; u(0, z, r) = ut0 (z, r) ; w(0, z, r) = wt0 (z, r) ; T (0, z, r) = Tt0 (z, r)
– La viscosité dynamique ne dépend que de la température qui elle-même ne dépend que de la coor-
donnée axiale z, µ(T ) = µ(z).
– Il est à noter que nous avons négligé le rayonnement volumique interne qui est, d’après Glicksman
[34, 33], très faible par rapport au rayonnement surfacique et à la conduction longitudinale, pour
des températures inférieures à 1100◦ C. Le terme de rayonnement devient q r = σ(T ∗4 − T 4 ).
Les équations utilisées dans notre modèle sont définies ci-dessous Eqs. (6.9-6.11). Des détails sur
l’obtention de ces équations sont apportés dans l’annexe A.
∂S
+ Sw = 0 (6.9)
∂t
∂wS ∂w2 S ∂ ∂w √ √
ρ + = 3Sµ + σ πS − 2 πSτa + ρgS (6.10)
∂t ∂z ∂z ∂z
∂T S ∂T wS √
ρcp + = λair πN u(Ta − T ) + 2 πSσ̃(T ∗4 − T 4 )
∂t ∂z
2
∂w ∂ S∂T ∂(Sqzr )
+3µS +λ − (6.11)
∂z ∂z ∂z ∂z
où L est une longueur introduite par Ph. Demont [21]. Elle est mesurée entre le fond de filière et le
bas du téton ; zc est une longueur caractéristique ajustable zc ∼ 3mm.
114 Chapitre 6. Modélisation physique de l’écoulement
– Nous utilisons la loi d’évolution de la conductivité thermique de l’air introduite par D. Lemmonier
[19] :
λair (z) = 2, 92.10−4[Tf ilm ]0,8
– Pour la viscosité cinématique de l’air [50] :
1,72
Tf ilm
νair (z) = 1.58.10−5
300
– Pour le nombre de Nusselt, nous prenons celui introduit par Glicksman [33].
1 3 4νair z
N u(z) = 4, 3 − 3 avec ε = max 3, ln
ε ε w0 R02
Caractéristiques prises pour le verre-E :
– Pour le couplage mécanique/thermique, nous prenons une viscosité qui dépend de la température
selon la loi expérimentale de Tamman-Vogel-Fulcher :
bµ
log(µ) = aµ + (6.12)
T − Tf ig
où µ est exprimé en Poises, et T en Kelvin, aµ , bµ , Tf ig varient avec le verre utilisé avec pour un
verre-E ”classique” : aµ = −1.65, bµ = 3140 et Tf ig = 790 + 273◦ K,
– La chaleur spécifique est prise comme constante tant que la température du verre est supérieure à
la température de transition vitreuse T g (T g = 731 + 273◦K soit µ = 1013 P oise) :
cp = 1238J/kg/K
ρ = 2450 − 2600kg/m3
λ = 1, 8W/m/K
6.2.3 Normalisation
Les équations 6.9 à 6.11 sont adimensionnées de manière à réduire les ordres de grandeur entre les
différents termes des équations et à mettre en évidence les paramètres de similitude.
∗ Choix des paramètres caractéristiques :
– la longueur, R0 ,
– la température, T0 ,
– la viscosité, µ0 ,
– le temps, R0 /w0 ,
– la masse volumique, ρ.
∗ Nombres caractéristiques sans dimensions :
– le nombre de Reynolds, e = ρw0 R0 /(3µ0 ),
6.2. Le modèle développé 115
avec Θf ig = Tf ig /T0 .
∗ Conditions extérieures sans dimension :
– la température ambiante près de la fibre, Θa ,
– la température de l’environnement radiatif, Θ∗ ,
∗ Grandeurs adimensionnées des différents termes de l’équation :
– le coefficient d’échange thermique (convection libre), H = (λf N u)/(ρcp w0 R0 ),
– le coefficient d’échange radiatif, E = (2σ̃T03 (x))/(ρcp w0 ),
– la tension de surface, Ω = σ/(3µ0 w0 ),
– le terme de convection forcée, C1 = e,
– le cisaillement de l’air, C2 = 2R0 τa /(3µ0 w0 ),
– la gravité, C3 = ρgR02 /(3µ0 w0 ),
– la diffusion visqueuse, D1 = 3µ0 w0 /(ρcp T0 R0 ),
– la conduction, D2 = λ/(ρcp w0 R0 ),
– la transfert radiatif suivant z, D3 = σT03 /(ρcp w0 ).
6.2.4 Discrétisation
Les équations (6.13-6.15) sont discrétisées par la méthode des différences finies. La discrétisation de
l’équation de conservation de la masse se fait au pas d’espace i − 1/2 et au pas de temps n + 1/2, celle
116 Chapitre 6. Modélisation physique de l’écoulement
de la quantité de mouvement au pas d’espace i et au pas de temps n + 1/2, et celle de l’énergie au pas
d’espace i et au pas de temps n + 1/2. Le pas d’espace est égal à ∆x avec ∆x = L/(N bR0 ), Nb étant le
nombre de points du maillage. L’erreur spatiale est d’ordre 2. Le pas de temps sera défini ultérieurement.
Les schémas sont présentés en détail dans l’annexe B. On considère deux schémas distincts : le schéma
implicite et le schéma semi-implicite.
- le schéma implicite :
On se place au temps n+1. La stabilité du schéma est bonne mais l’erreur temporelle du schéma est
d’ordre 1.
- Le schéma semi-implicite :
De manière à augmenter la précision temporelle, nous avons choisi de développer un schéma de Cranck-
Nicolson linéarisé. On pose X n+1/2 = (1 − β)X n + βX n+1 avec β compris entre 0 et 1, dans le schéma
semi-implicite on pose, β = 0.5. Pour le schéma implicite on pose β = 1, pour un schéma explicite β = 0,
ce dernier cas pose des problèmes de stabilité. L’erreur temporelle du schéma est d’ordre 2.
Fig. 6.2 – Comparaison des résultats analytiques et numériques obtenus pour un jet isotherme lorsque :
cas 1, la viscosité est dominante ; cas 2, la viscosité et les forces d’inertie sont dominantes.
6.3. Validations et exemples de résultats 117
Fig. 6.3 – Comparaison des profils de température obtenus avec Modfib-substitution, Modeft et les points
expérimentaux de l’expérience n◦ 1 de Clementin-de-Leusse.
- Expérience n◦ 2/Modeft/GY/Modfib
La figure 6.4 compare le profil de température du jet obtenu avec ModFib-substitution et ModFib-
Newton Raphson, les codes Modeft et Gy (Saint-Gobain Recherche), des données expérimentales de C.
Clémentin-de Leusse [19]. Conditions aux limites prisent pour ce calcul : T0 = 1184◦C, w0 = 1, 6.10−3m/s,
wL = 22, 4m/s.
Les prédictions de notre modèle sont ”correctes”. Il semble cependant que la version ModFib-Newton
Raphson donne de meilleurs résultats. En ce qui concerne le profil du bulbe, figure 6.5, il semble que
notre modèle sur-estime la contraction du jet (pour ce cas, le code Modeft ne congerge pas)2 .
1 Obtenues par spectroscopie Brillouin, précision attendue ≈ ±5◦ C
2 Rigoureusement, les modèles 1D comme Modeft, Gy et Modfib, ne sont valides que pour des distances z > 3R0 ...
L’intérêt de ces comparaisons est donc limité... Cependant, elles nous permettent de comparer nos résultats avec ceux
d’études antérieures.
118 Chapitre 6. Modélisation physique de l’écoulement
Fig. 6.4 – Comparaison des profils de température obtenus avec Modeft, Gy et Modfib et les points
expérimentaux de l’expérience n◦ 2 de Clementin-de-Leusse.
Fig. 6.5 – Evolution du rayon de la fibre dans un cas ”défavorable” : cas 11 [19].
• Temps de calculs
Dans le tableau ci-suivant, nous comparons les temps de calcul de ModFib3 , pour un cas stationnaire et
les deux méthodes de résolution.
La méthode de Newton-Raphson nécessite des temps de calcul beaucoup plus longs. Cependant, si
pour un cas stationnaire les temps de calculs nécessaires avec la méthode de Newton-Raphson ne sont
par prohibitifs, ils le deviennent pour un cas instationnaire (le calcul du cas présenté dans la figure 6.8
nécessite déjà 8 heures avec la méthode de substitution). C’est pourquoi, dans ce qui suit, les études en
instationnaire ont été effectuées avec la version de ModFib qui utilise la méthode de substitution.
Perturbation convective
On considère ici que le jet de verre est soumis à une perturbation convective qui prend la forme d’un
écoulement d’air axisymétrique dirigé perpendiculairement à l’axe de la fibre. Ce flux d’air va entraı̂ner
une modification des échanges convectifs entre la fibre et l’air ambiant. Pour simuler numériquement cette
3 Les calculs ont été effectués sur un Pentium II à 266Mhz et avec 64Mo de RAM.
4 Alors que notre modèle le suppose constant.
120 Chapitre 6. Modélisation physique de l’écoulement
perturbation, nous prenons pour le nombre de Nusselt une loi d’évolution du type N u(t) = C(t)Re(t))en(t)
où C(t) et n(t) dépendent du nombre de Reynolds de l’écoulement d’air [100]. L’évolution temporelle du
2
profil de vitesse de l’écoulement d’air est supposé gaussien avec vair = vmax e−(t−t0) /σ .
La figure 6.7, présente l’évolution du diamètre de la fibre formée dans ces conditions, avec vmax = 0.05m/s
√
et ≈ 4 σ = 0.4s. Le flux d’air entraı̂ne un fort accroissement du diamètre de la fibre (i.e. le bulbe rétrécit)
puis, après relaxation, le diamètre terminal de la fibre revient à sa valeur initiale. Cette évolution est
assez conforme avec nos observations expérimentales [77].
Fig. 6.7 – Evolution du diamètre de la fibre produite lorsque le jet est soumis à une perturbation convective.
Fig. 6.8 – Perturbation radiative : évolution du diamètre de la fibre produite lorsqu’on augmente bruta-
lement, à t=0, la température de rayonnement du milieu ambiant.
Perturbation radiative
On simule ici l’effet d’une perturbation radiative externe sur le jet de verre. La perturbation prend
la forme d’une augmentation de la température de rayonnement du milieu ambiant5 . Cette perturbation
est maintenue [52].
Les figures 6.8 à 6.11 présentent les évolutions obtenues pour : le diamètre de la fibre produite, le profil
radiale du jet de verre et son profil longitudinal de température, la tension de fibrage dûe à la viscosité.
On peut tenter d’interpréter les résultats obtenus en se rappelant que, dans notre modèle, le débit
de la filière est posé comme constant. L’augmentation de la température du milieu ambiant induit un
allongement du bulbe. L’accumulation de verre, au niveau du bulbe, implique alors une diminution du
diamètre de la fibre produite. Une fois que le nouveau régime est établi, le diamètre de la fibre reprend
sa valeur initiale.
5 dû fait de la casse d’une autre fibre par exemple...
6.3. Validations et exemples de résultats 121
Fig. 6.9 – Perturbation radiative : évolution du profil axial du jet : avant et après l’application de la
perturbation.
Fig. 6.10 – Perturbation radiative : évolution du profil axial de température du jet de verre.
Fig. 7.1 – Dispositif expérimental mis en place pour l’étude de la stabilité du procédé.
123
124 Chapitre 7. Etude expérimentale de la stabilité du procédé
– Le creuset céramique est alimenté en verre à l’aide de billes dont la masse est de l’ordre de 15g. Bien
que la filière soit dotée d’un système d’alimentation automatique en billes, nous n’avons pas utilisé
ce système. En effet, il n’est pas synchronisable et il génère un bruit électrique/électromagnétique
qui perturbe les autres systèmes de mesure.
– La température du bain de verre est contrôlée par un automate programmable. La régulation est
basée sur une mesure de la température du bain par thermo-couple.
– La hauteur h du bain de verre est contrôlée en immergeant une tige d’acier inoxidable dans le
creuset. La hauteur de ”mouillage” de la tige nous permet d’en déduire la hauteur de verre dans le
creuset. De la mesure de cette hauteur et de la température du bulbe on peut déduire la tirée (i.e.
débit de la filière).
– La température du bulbe est contrôlée à l’aide d’un pyromètre à dispersion de filament (lecture
”manuelle” uniquement). La mesure de luminance est effectuée pour une longueur d’onde de λ ≈
0.6µm et pour un angle d’inclinaison dans le plan vertical d’approximativement 45◦ . Ce pyromètre
a été calibré à l’aide d’un corps noir et sa mesure a été comparée avec celle du système PyroVerre
de Saint-Gobain Recherche (mesure à λ ≈ 5µm). La figure 7.26 indique la zone du bulbe où la
luminance est mesurée.
– Les mouvements d’oscillation de la fibre sont atténués au moyen de deux crayons en graphite (i.e.
embarrage).
– La vitesse de fibrage Vf est déduite de la mesure de la vitesse de rotation du bobinoir. Cette dernière
est elle-même déduite de l’analyse de la fréquence du signal produit par une roue codeuse, fixée à
l’arbre du moteur du bobinoir, à l’aide d’un fréquence-mètre ou du système F IBS (utilisation de
la troisième voie de numérisation et analyse par CSDW).
– Le bobinoir est fixé sur un profilé motorisé (moteur pas à pas). Ceci permet d’assurer une meilleure
homogénéité de l’enroulement de la fibre sur le bobinoir. Ce système de chariotage est contrôlé
par un PC. La vitesse de chariotage est typiquement de 10mm/s pour une course de 20mm (i.e.
fréquence : 2Hz).
– L’interféromètre FIBS est positionné d’un côté de la fibre (voir le chapitre 3). Il est contrôlé par
un PC. Ce dernier numérise les signaux Doppler mais aussi le signal de vitesse du bobinoir, il est
également synchronisé avec le diffractomètre et le système d’imagerie du bulbe.
– Le diffractomètre HIREDI est positionné de part et d’autre de la fibre, voir le chapitre 4.
– Un système d’imagerie vidéo permet d’enregistrer la forme du bulbe au cours du temps, voir la
figure 7.26.
– Un boı̂tier de synchronisation permet de générer un signal de déclenchement vers les trois systèmes
de mesure et donc, la réalisation de mesures synchrones : Diamètre interféromètre, Diamètre dif-
fractomètre, Images du Bulbe et Vitesse d’étirage.
– L’alimentation en verre, la mesure de la hauteur du bain de verre et de la température du bulbe
sont réalisées manuellement...
60
50
Vitesse [m/s]
40
30
20
10
Fig. 7.3 – Photographie du bulbe à faible température de fibrage : il n’est pas axisymétrique, on observe
également le passage d’une bulle.
La figure 7.4 montre une évolution temporelle typique obtenue pour les paramètres de fibrage suivants :
(T0 , Vf ) = (1187, 28.0) [◦ C, m/s] et h = 139mm4. Le diamètre moyen de cette fibre est centré sur
D = 13.87µm.
En bas de la figure 7.4 on trouve également la densité spectrale de puissance (DSP) correspondant à
la série temporelle. Cette DSP a été calculée par une méthode classique de moyenne de spectres. Son
évolution indique que les ”fluctuations” ou ”oscillations” du diamètre de la fibre sont quasi-périodiques,
de fréquence moyenne : νosc. ≈ 0.75Hz.
La distribution obtenu pour les diamètres est présentée par la figure 7.5. Elle est approximativement de
forme gaussienne 5 . L’amplitude des fluctuations pour une valeur de l’écart-type (1 − σD ), ou 2 valeurs
(2 − σD ) atteint ≈ 0.29µm et ≈ 0.58µm.
La figure 7.6 montre la distribution correspondante du taux de fluctuation du diamètre, en fonction du
temps : dD(t)/dt, et en fonction de la longueur de fibre : dD(L)/dL6 . Pour un interval de confiance de
2 − σD ce taux atteint ≈ 3.8µm/s et ≈ 0.14µm/m, soit une variation relative de ≈ 27%/s et ≈ 1%/m.
A travers cet exemple on constate que, même pour des conditions stationnaires, le diamètre de fibrage
3 Sauf précision, ces mesures ont été obtenues avec le diffractomètre HIREDI, en effet le taux de fluctuation du diamètre
des fibres est de l’ordre de grandeur de l’écart-type de la réponse de l’interféromètre FIBS.
4 Pour tous les résultats présentés ici nous avons essayé de maintenir constante cette dernière valeur.
5 Pour des temps d’acquisition plus courts, du fait de la quasi-périodicité des oscillations, elle diffère significativement de
la forme gaussienne.
6 Nous avons simplement utilisé le fait que dL = V dt
f
126 Chapitre 7. Etude expérimentale de la stabilité du procédé
Fig. 7.4 – Régime stationnaire : évolution temporelle du diamètre et Densité Spectrale de Puissance
(DSP) ; Diamètre moyen et écart-type : D = 13.87µm et σD = 0.29µm, conditions de fibrage : (T0 , Vf ) =
(1187, 28.0) [◦ C, m/s] et h = 139mm.
fluctue de manière assez importante. Dans ce qui suit, nous nous intéressons successivement à l’évolution
de la valeur moyenne puis à l’écart-type du diamètre des fibres produites, pour différentes conditions de
fibrage. Nous assimilerons l’amplitude moyenne des fluctuations du diamètre à l’écart-type de la série
temporelle.
Fig. 7.5 – Distribution statistique du diamètre de la fibre produite pour (T0 , Vf ) = (1187, 28.0) [◦ C, m/s]
et h = 139mm.
Fig. 7.6 – Distribution statistique du taux de fluctuation spatiale et temporelle du diamètre de la fibre
produite pour (T0 , Vf ) = (1187, 28.0) [◦ C, m/s] et h = 139mm.
Fig. 7.8 – Evolution du diamètre moyen en fonction de la vitesse de fibrage et pour trois températures de
fibrage.
7.2. Régimes stationnaires 129
Fig. 7.10 – Ecart-type absolu du diamètre en fonction de la température de fibrage et du diamètre moyen.
130 Chapitre 7. Etude expérimentale de la stabilité du procédé
Fig. 7.12 – Ecart-type relatif du diamètre en fonction de la température de fibrage et du diamètre moyen.
7.2. Régimes stationnaires 131
Fig. 7.13 – Ecart-type absolu du diamètre en fonction de la vitesse de fibrage et du diamètre moyen.
Fig. 7.14 – Ecart-type relatif du diamètre en fonction de la vitesse de fibrage et du diamètre moyen.
132 Chapitre 7. Etude expérimentale de la stabilité du procédé
Fig. 7.15 – Variation de la fréquence temporelle des fluctuations du diamètre en fonction de la température
de fibrage (pour différentes vitesses de fibrage).
Oscillations du bulbe
La figure 7.16 compare l’évolution du diamètre de la fibre et les variations relatives de la hauteur du
centroı̈de de l’image du bulbe9 , pour un régime stationnaire (T0 , Vf ) = (1170◦C, 9.0 m/s). Le diamètre de
la fibre oscille10 , tout comme la hauteur du centroı̈de du bulbe. Ces fluctuations ne semblent pas corrélées...
Les Densités Spectrales de Puissance (DSP) montrent que les variations de hauteur du centröide se font
principalement à des fréquences autour de 0.50 − 0.68Hz alors que, pour le diamètre, la fréquence des
oscillations est plutôt de 0.60 − 0.78Hz. On peut s’étonner de trouver des gammes de fréquences si proches
et pourtant distinctes... On remarquera cependant qu’elles se recouvrent en partie et que leur largeur
est la même : 0.18 Hz. Si, comme nous l’intuitons, les fluctuations du diamètre de la fibre sont dues aux
oscillations du bulbe, il semble que les fréquences générées par le bulbe soient ”filtrées” : puisque l’on
passe de à 0.50 − 0.68 → 0.60 − 0.78Hz.
Dans l’état actuel, il nous est difficile d’établir l’origine de ces oscillations car plusieurs mécanismes
couplés entrent probablement en jeux. Sans être exhaustifs on peut remarquer que :
- En amont du téton, les forts gradients thermiques induisent probablement des recirculations et donc,
des fluctuations de la viscosité du verre, des variations du diamètre de la fibre formée. Avec les trans-
ferts conductifs et convectifs, on s’attend à ce que les perturbations correspondantes soient de très
basse fréquence. Cependant, dans cette zone de l’écoulement, les transferts radiatifs sont généralement
considérés comme prépondérants. Les perturbations ”thermiques” peuvent donc se propager ”instan-
tanément”, en amont comme en aval, et sur des distances importantes11 .
- Les gradients thermiques sont encore plus marqués dans la zone d’écoulement aval : la zone du
bulbe. Les dimensions du bulbe sont relativement faibles, l’inertie thermique de ce volume fluide est donc
d’autant réduite. Les effets évoqués précédemment devraient donc être plus prononcés dans cette zone.
A ces effets pourraient s’ajouter les effets capillaires liés à la présence de l’interface air-verre (bien que,
généralement, on les considère comme négligeables...).
9 Dans l’exemple de la figure 7.16, le centroı̈de est indiqué par une petite croix.
10 On remarquera également les nombreuses signatures du passage de bulles : brusques variations du ”diamètre” mesuré.
Nous avons dénombré pas moins de 20 bulles en 32 s !
11 Le libre parcours moyen des photons ”thermiques” de longueur d’onde λ = 2.2µm (i.e. ≈ 1045◦ C) est de ≈ 330mm...
7.2. Régimes stationnaires 133
- Les vibrations mécaniques de la fibre formée pourraient également influer sur l’écoulement visqueux
amont, de manière directe ; ou à travers le rayonnement thermique (modification de la forme du guide
d’onde)...
ont été effectuées sur moins de 100 mesures. Les barres d’erreurs correspondent aux écart-types obtenus :
pour le diamètre (instabilité naturelle du procédé et résolution de l’interféromètre) et pour la fréquence
(fluctuations du diamètre et de la mesure de ν1 ).
En utilisant l’équation 2.35 et les points expérimentaux de la figure 7.17, on peut déterminer l’évolution
correspondante de la tension d’étirage, voir la figure 7.18. Les barres d’erreur sur la tension rendent compte
des fluctuations du diamètre de la fibre et de l’écart-type sur les autres grandeurs mesurées : ν1 et Vf
(résolution ≈ 1%). Dans les figures 7.17 et 7.18, la fréquence de vibration et la tension semblent suivre une
loi du type exponentielle décroissante. Ce résultat est en accord avec les résultats que l’on peut obtenir
à partir de considérations sur hydrodynamique du procédé [34].
Les résultats expérimentaux donnent des tensions relativement importantes. Ceci s’explique par la re-
lative faiblesse de la température de fibrage et par le fait que l’embarrage introduit, très certainement,
une tension additionnelle (par frottement). En utilisant alternativement, deux puis trois crayons pour
l’embarrage, on doit pouvoir mesurer cette tension additionnelle.
Fig. 7.19 – FIBS : (a) Evolution temporelle et (b) statistique du diamètre mesuré lorsque la vitesse de
fibrage est diminuée par paliers successifs.
On peut forcer d’avantage encore le procédé en imposant des rampes sur la vitesse de fibrage, continues,
et avec une des accélérations plus élevées. Expérimentalement, ceci peut se faire en freinant énergiquement
le bobinoir.
Les figures 7.20, 7.21 et 7.22, présentent les résultats obtenus pour trois décélérations γ = −5.4m/s2,
−7.4m/s2, −20.0m/s2 et pour, respectivement, des températures de fibrage de 1170◦C, 1200◦C et
1220◦C. Les évolutions tracées sont : le diamètre mesuré avec l’interféromètre et le diamètre prédit
numériquement12 , la vitesse de fibrage et l’accélération associée, le débit instantané déduit de la mesure
du diamètre et de la vitesse instantanée de fibrage.
On remarque que les effets induits par les rampes de la vitesse de fibrage sont assez faibles. Pour la
rampe n◦ 3, le débit instantané de la filière oscille. En fait, le diamètre de la fibre oscille autour de la
valeur prédite par le code ModFib. Or, ce dernier est basé sur une hypothèse de débit constant pour la
filière. Le procédé semble donc toujours être dans un régime quasi-stationnaire.
La figure 7.23 présente, pour une rampe de vitesse très comparable à la rampe n◦ 3, l’évolution du
diamètre de la fibre produite, du débit de verre associé, de la vitesse de fibrage et de deux dimensions
du bulbe. La dimension latérale du bulbe LD est mesurée à 1.17mm de l’orifice de sortie du téton. La
dimension longitudinale du bulbe LZ correspond à la distance entre l’orifice de sortie du téton et la zone
sombre, caractéristique, en bas du bulbe. Les images correspondantes du bulbe sont présentées par la
figure 7.24.
Fig. 7.20 – Evolution du diamètre et du débit instantané lorsque la vitesse de fibrage est diminuée rapi-
dement, la décélération moyenne est ici de γ ≈ −5.4m.s−2 .
Fig. 7.21 – Idem figure précédente pour une décélération moyenne de γ ≈ −7.4m.s−2 .
Fig. 7.22 – Idem figure précédente pour une décélération moyenne de γ ≈ −20.0m.s−2 .
7.3. Régimes transitoires 137
Fig. 7.23 – Comparaison de l’évolution de la forme du bulbe (dimensions latérale et axiale) avec le
diamètre de fibrage et le débit instantané. Rampe de vitesse d’accélération maximale : γ = −20.7m/s et
T0 = 1220◦C.
138 Chapitre 7. Etude expérimentale de la stabilité du procédé
Fig. 7.24 – Evolution de la forme du bulbe pour la rampe n◦ 3 de vitesse. Les prises de vues sont espacées
de ≈ 1/10s.
7.3. Régimes transitoires 139
Les prises de vues de la figure 7.24 sont espacées de ≈ 1/10s et l’instant de prise de vue est déduit
du numéro de l’image par l’expression : Timage = −12.36 + 0.108 ∗ F rame. Ces données ont été obtenues
de manière synchrone, pour T0 = 1220◦C, et pour une rampe de vitesse de décélération maximale :
γ = −20.7m/s2.
Au regard de la figure 7.23, les corrélations entre la forme du bulbe, ses dimensions, et les autres grandeurs
mesurées, ne sont pas évidentes...13 . LZ semble légèrement osciller, comme nous l’avions déjà observé pour
les régimes stationnaires.
On remarquera qu’à partir du temps t ≈ 3.7s, le diamètre de la fibre n’est plus mesuré. En fait, comme
le montre la figure 7.24 pour Frame ≥ 154, il est manifeste que la fibre a cassé : le verre tend à former
une goutte de verre qui va finir par tomber.
Fig. 7.25 – Evolution du diamètre de la fibre produite lorsque l’on impose une rampe de température, avec
pour i) ∆T = +12◦ C/mn, Vf = 25m/s ; ii) ∆T = −37.5◦C/min et Vf = 20m/s.
augmente de (∆T = +12◦ C/min). L’évolution correspondante de la forme du bulbe est présentée par la
figure 7.26, pour T0 = 1145 −→ 1245◦C.
Les valeurs obtenues pour le diamètre de la fibre sont similaires à celles obtenues pour des conditions
stationnaires, voir la figure 7.7. Nous sommes donc, ici encore, dans un régime quasi-stationnaire.
7.4 Conclusion
Dans ce chapitre nous avons réalisé une étude préliminaire, sur la stabilité de la filière mono-téton en
fonction des deux paramètres principaux de contrôle du procédé : la température et la vitesse de fibrage
avec (T0 , Vf ) ∈ (1145 − 1245, 5 − 65) [◦ C, m/s]. Sur ce domaine, le diamètre de la fibre varie entre 10µm
et 35µm.
Pour des conditions de fibrage stationnaires, le diamètre de la fibre fluctue de manière assez im-
portante. En assimilant l’amplitude des fluctuations à l’écart-type des séries temporelles, on trouve que
l’amplitude des fluctuations est comprise entre ±1.8 et ±3.7% (i.e. ±3.6 et ±7.4% pour 2 − σD ). Il
existe deux domaines de température pour lesquels l’amplitude relative des fluctuations est minimale :
i) T < 1175◦ C et ii) T > 1230◦C. Le minimum absolu σD /D ≈ 1.8% est observé pour le domaine
(T0 , Vf ) ∈ (1145 − 1175, 18 − 45) [◦ C, m/s] et des fibres dont le diamètre est inférieur à ≈ 15µm. Ces
fluctuations, ou ”oscillations”, sont quasi-périodiques. Leur fréquence augmente de 0.5Hz à 0.9Hz pour
des températures allant de 1145◦C à 1245◦C. Pour une fibre de D ≈ 14µm produite pour (1187, 28.0)
[◦ C, m/s], le taux des fluctuations temporelle et spatiale du diamètre atteint ±1.9µm/s et ±0.07µm/m à
1 − σD . Nous intuitons que ces perturbations trouvent leur origine au niveau du bulbe. Ce dernier oscille
également, à des fréquences proches de celles mesurées pour le diamètre de la fibre. L’origine réelle de
l’instabilité du bulbe nous est cependant inconnue.
Nous avons également étudié la réponse du procédé à des conditions transitoires de fibrage : rampes
de la vitesse ou de la température de fibrage. Nous avons constaté que, à chaque instant, le diamètre de
la fibre ainsi produite est très similaire à celui que l’on obtient pour le régime stationnaire correspondant.
En fait, tout se passe comme si le procédé était toujours dans un régime quasi-stationnaire (même pour
des décélérations de la vitesse de fibrage de près de −20m/s2 ).
La réponse du procédé à des perturbations singulières (perturbation convective, bulles dans la compo-
sition...) n’a pas été abordée dans ce manuscrit, faute de temps. Cependant, différents tests expérimentaux
ont montré que l’interféromètre FIBS, tout comme le diffractomètre HIREDI, pouvaient être utilisés pour
étudier la réponse du procédé à ces diverses perturbations [77, 76].
Conclusion générale
Au terme de ce travail, différents outils théoriques, numériques et expérimentaux ont été développés
pour permettre l’étude des conditions de stabilité, mais également le contrôle en ligne, d’un procédé de
fabrication de fibres de verre utilisées pour des applications de renfort. Ce travail a été réalisé au sein
de l’UMR CNRS 6595-IUSTI, à Marseille, et il a été financé par le groupe Saint-Gobain Recherche et
l’ADEME. Il concilie de ce fait des aspects de recherche fondamentale, une problématique industrielle et
une problématique environnementale. Ce travail revêt également un caractère pluridisciplinaire, de part
la physique du procédé et les outils mis en place pour son étude : hydrodynamique et thermique des
jets visqueux, électromagnétisme, instrumentation scientifique. . . De manière synthétique, les résultats
obtenus peuvent être classés en trois catégories :
Instrumentation scientifique :
Conclusions :
Un Interféromètre phase Doppler (FIBS) a été développé spécifiquement pour ce travail [74, 75, 80, 77].
Ce système, par ses fonctionnalités, son encombrement et sa résolution, répond au cahier des charges défini
au début de cette étude avec le groupe Saint-Gobain Recherche. Il a été optimisé pour la caractérisation
des fibres de verre-E (indice élevé, m ≈ 1.555). Sa résolution est de l’ordre de σD ≈ 0.35µm sur la
gamme de tailles D = 5 − 42µm avec un taux d’acquisition qui peut atteindre ≈ 50kHz (i.e. une mesure
tous les millimètres pour Vf = 50m/s). Ce système a été testé sur une filière mono-téton implantée
au laboratoire et sur deux filières semi-industrielles implantées sur le site de Saint-Gobain Recherche à
Aubervilliers. Entre autres originalités, cet interféromètre utilise trois faisceaux cohérents pour former le
volume de mesure [78, 79, 82], ce qui permet d’étendre sa dynamique sur les tailles mesurables mais aussi
14 Pour de grand paramètres de taille et plusieurs milliers de couches
141
142 Chapitre 7. Etude expérimentale de la stabilité du procédé
de déterminer la position de la fibre le long de l’axe optique. Ce système permet également de déterminer
la tension de fibrage, en ligne, au moyen de la technique de la «corde vibrante» [81].
Un diffractomètre haute résolution (HIREDI) a été développé afin de valider la réponse de l’in-
terféromètre [76]. Ce système et son algorithme original de traitement permettent d’inverser les contribu-
tions des Résonances Morphologico-Dépendantes (MDRs). Sa résolution est de l’ordre de σD ≈ 0.02µm
pour une gamme de tailles de D = 8 − 40µm avec un taux maximal d’acquisition de 43kHz. Cependant,
du fait de son encombrement et de la distance de mesure15 , ce diffractomètre est difficilement utilisable
sur site industriel, dans sa version actuelle.
Perspectives :
Ces deux systèmes optiques peuvent être utilisés pour caractériser d’autres types de fibres : polymère,
carbone, métalliques, etc. d’autres gammes de tailles.
Dans leur version actuelle, ils ne permettent pas de mesurer le diamètre de fibres d’orientation incon-
nue. Cependant, ils peuvent détecter les fibres inclinées et donc permettre de déterminer la granulométrie
de fibres pour des angles d’inclinaison donnés.
La résolution de l’interféromètre laser pourrait être augmentée en déconvoluant, ou en atténuant, les
contributions des MDRs. Plusieurs solutions ont été proposées : déconvolution des distributions, utilisation
de probabilités jointes, utilisation d’une source de plus grande largeur spectrale. Il conviendrait d’explorer
plus en détail ces différentes solutions.
Le diffractomètre laser pourrait être miniaturisé afin de permettre des mesures sur site, il présente en
effet une résolution élevée.
effets instationnaires (la fibre se comporte comme un guide d’onde pour les photons thermiques émis par
le four). Une solution possible, et à moindre coût, consisterait à développer un modèle radiatif 2D, du
type Monte Carlo, qu’il conviendrait de coupler de manière itérative à notre modèle d’écoulement.
Suite aux résultats expérimentaux obtenus, il nous est apparu que la filière mono-téton n’est pas suf-
fisamment contrôlée et instrumentée pour que l’on puisse développer beaucoup plus en avant les études
de stabilité. Pour les études avenir, il conviendrait de modifier cette filière afin notamment : i) de mieux
contrôler son débit (par ajout en continu d’une poudre de verre préchauffée) ; ii) d’améliorer l’isola-
tion du bulbe vis-à-vis des perturbations externes (vibrations mécaniques induites par le transformateur
électrique, convection, turbulence induite par le bobinage...) ; iii) de mesurer avec plus de précision la
hauteur de bain de verre et la température du bulbe.
144 Chapitre 7. Etude expérimentale de la stabilité du procédé
Conclusion générale
Au terme de ce travail, différents outils théoriques, numériques et expérimentaux ont été développés
pour permettre l’étude des conditions de stabilité, mais également le contrôle en ligne, d’un procédé de
fabrication de fibres de verre utilisées pour des applications de renfort. Ce travail a été réalisé au sein
de l’UMR CNRS 6595-IUSTI, à Marseille, et il a été financé par le groupe Saint-Gobain Recherche et
l’ADEME. Il concilie de ce fait des aspects de recherche fondamentale, une problématique industrielle et
une problématique environnementale. Ce travail revêt également un caractère pluridisciplinaire, de part
la physique du procédé et les outils mis en place pour son étude : hydrodynamique et thermique des
jets visqueux, électromagnétisme, instrumentation scientifique. . . De manière synthétique, les résultats
obtenus peuvent être classés en trois catégories :
Instrumentation scientifique :
Conclusions :
Un Interféromètre phase Doppler (FIBS) a été développé spécifiquement pour ce travail [74, 75, 80, 77].
Ce système, par ses fonctionnalités, son encombrement et sa résolution, répond au cahier des charges défini
au début de cette étude avec le groupe Saint-Gobain Recherche. Il a été optimisé pour la caractérisation
des fibres de verre-E (indice élevé, m ≈ 1.555). Sa résolution est de l’ordre de σD ≈ 0.35µm sur la
gamme de tailles D = 5 − 42µm avec un taux d’acquisition qui peut atteindre ≈ 50kHz (i.e. une mesure
tous les millimètres pour Vf = 50m/s). Ce système a été testé sur une filière mono-téton implantée
au laboratoire et sur deux filières semi-industrielles implantées sur le site de Saint-Gobain Recherche à
Aubervilliers. Entre autres originalités, cet interféromètre utilise trois faisceaux cohérents pour former le
volume de mesure [78, 79, 82], ce qui permet d’étendre sa dynamique sur les tailles mesurables mais aussi
14 Pour de grand paramètres de taille et plusieurs milliers de couches
141
142 Chapitre 7. Etude expérimentale de la stabilité du procédé
de déterminer la position de la fibre le long de l’axe optique. Ce système permet également de déterminer
la tension de fibrage, en ligne, au moyen de la technique de la «corde vibrante» [81].
Un diffractomètre haute résolution (HIREDI) a été développé afin de valider la réponse de l’in-
terféromètre [76]. Ce système et son algorithme original de traitement permettent d’inverser les contribu-
tions des Résonances Morphologico-Dépendantes (MDRs). Sa résolution est de l’ordre de σD ≈ 0.02µm
pour une gamme de tailles de D = 8 − 40µm avec un taux maximal d’acquisition de 43kHz. Cependant,
du fait de son encombrement et de la distance de mesure15 , ce diffractomètre est difficilement utilisable
sur site industriel, dans sa version actuelle.
Perspectives :
Ces deux systèmes optiques peuvent être utilisés pour caractériser d’autres types de fibres : polymère,
carbone, métalliques, etc. d’autres gammes de tailles.
Dans leur version actuelle, ils ne permettent pas de mesurer le diamètre de fibres d’orientation incon-
nue. Cependant, ils peuvent détecter les fibres inclinées et donc permettre de déterminer la granulométrie
de fibres pour des angles d’inclinaison donnés.
La résolution de l’interféromètre laser pourrait être augmentée en déconvoluant, ou en atténuant, les
contributions des MDRs. Plusieurs solutions ont été proposées : déconvolution des distributions, utilisation
de probabilités jointes, utilisation d’une source de plus grande largeur spectrale. Il conviendrait d’explorer
plus en détail ces différentes solutions.
Le diffractomètre laser pourrait être miniaturisé afin de permettre des mesures sur site, il présente en
effet une résolution élevée.
effets instationnaires (la fibre se comporte comme un guide d’onde pour les photons thermiques émis par
le four). Une solution possible, et à moindre coût, consisterait à développer un modèle radiatif 2D, du
type Monte Carlo, qu’il conviendrait de coupler de manière itérative à notre modèle d’écoulement.
Suite aux résultats expérimentaux obtenus, il nous est apparu que la filière mono-téton n’est pas suf-
fisamment contrôlée et instrumentée pour que l’on puisse développer beaucoup plus en avant les études
de stabilité. Pour les études avenir, il conviendrait de modifier cette filière afin notamment : i) de mieux
contrôler son débit (par ajout en continu d’une poudre de verre préchauffée) ; ii) d’améliorer l’isola-
tion du bulbe vis-à-vis des perturbations externes (vibrations mécaniques induites par le transformateur
électrique, convection, turbulence induite par le bobinage...) ; iii) de mesurer avec plus de précision la
hauteur de bain de verre et la température du bulbe.
144 Chapitre 7. Etude expérimentale de la stabilité du procédé
Annexe A
∂ρ 1 ∂(ρru) ∂(ρw)
+ + = 0 (A.1)
∂t r ∂r ∂z
∂(ρu) 1 ∂(ρru2 ) ∂(ρwu) 1 ∂ (r(τrr − p)) ∂τzr ρu2 τθθ − p
+ + = + + + (A.2)
∂t r ∂r ∂z r ∂r ∂z r r
∂(ρw) 1 ∂(ρruw) ∂(ρw2 ) ∂(τzz − p) 1 ∂(rτrz )
+ + = + + ρg (A.3)
∂t r ∂r ∂z ∂z
r ∂r
∂(ρT cp ) 1 ∂(ρruT cp ) ∂(ρwT cp ) 1 ∂(rqr ) ∂qz
+ + = Φ− + (A.4)
∂t r ∂r ∂z r ∂r ∂z
Les différentes composantes du tenseur des contraintes et le flux de chaleur (6.6,6.8) deviennent :
∂u 2 1 ∂(ru) ∂w
τrr = +2µ − µ +
∂r 3 r ∂r ∂z
u 2 1 ∂(ru) ∂w
τθθ = 2µ − µ +
r 3 r ∂r ∂z
∂w 2 1 ∂(ru) ∂w
τzz = 2µ − µ +
∂z 3 r ∂r ∂z
∂w ∂u
τrz = τrz = µ +
∂r ∂z
2 2
∂u
2 u 2 ∂w ∂w ∂u
Φ = 2µ + + +µ +
∂r r ∂z ∂r ∂z
(A.5)
On peut écrire les équations précédentes (A.1-A.4) sous une forme du type ’transport-diffusion’. Cette
forme sera utilisée ultérieurement pour faciliter l’intégration du système :
Avec φ remplaçant respectivement les quantités ρ, u, w et T dans les équations de A.1 à A.4, Az la
diffusion suivant z, Ar la diffusion suivant r, et ρU les forces conservatives
146 Annexe A. Formulation Intégrale du système
De ces projections, on peut alors déduire τn la contrainte sur la surface libre et ses projections dans
le repère (n, s) tangent à la surface de la fibre :
D’après les expressions précédentes A.8 et A.9, exprimant la contrainte de pression τnn (normale à la
surface libre), on peut déduire l’expression de la pression dans la fibre :
1 1
p = τnn + pa + σ +
Rzr Rnθ
2 −1/2 2
∂R ∂R ∂R 1 1
= 1+ τrr − 2 τrz + τzz + pa + σ + (A.11)
∂z ∂z ∂z Rzr Rnθ
q = qc + qr
qnc = h(T − Ta )
qnr = σ̃(T 4 − T ∗4 ) (A.12)
On calcule le coefficient d’échange convectif h à l’aide du nombre de Nusselt N u par la loi suivante :
hL 2Rh
Nu = = (A.13)
λair λair
147
Intégration
R R
∂φ ∂ ∂R
dr = φdr − φ
0 ∂t ∂t 0 ∂t
Ce théorème permet d’intégrer sur la section de la fibre les deux membres de l’équation ’transport
diffusion’ (A.6), référencés par les chiffres (I) et (II) :
R(z,t)
(I)dS = 2π (I)rdr
S(z,t) 0
S S S
∂ρφ ∂ρwφ 1 ∂ρruφ
= dS + dS + dS
0 ∂t 0 ∂z 0 r ∂r
s s
∂ ∂ ∂R ∂R
= ρφdS + ρwφdS − 2πRρφ + w + 2πRρφu
∂t 0 ∂z 0 ∂t ∂z
s s
∂ ∂
= ρφdS + ρwφdS (A.14)
∂t 0 ∂z 0
R(z,t)
(II)dS = 2π (II)rdr
S(z,t) 0
S S S
∂Az 1 ∂rAr
= dS + dS + ρU dS
0 ∂z 0 r ∂r 0
s S
∂R ∂
= 2πR Ar − Az + Az dS + ρU dS
∂z ∂z 0 0
s S
∂
= n.AδS + Az dS + ρU dS (A.15)
∂z 0 0
2
avec : δS = 2πR 1 + ∂R
∂z et A = (Ar , Az )
En remplaçant les expressions de (I) et de (II) dans les équations précédentes A.14 et A.15, par les
différents termes des équations de bilan (A.1-A.4) on obtient la formulation intégrale du modèle :
S S
∂ ∂
ρdS + ρwdS = 0 (A.16)
∂t 0 ∂z 0
S S S S
∂ ∂ ∂ τθθ
ρudS + ρwudS = τrz dS + τnr δS + dS (A.17)
∂t 0 ∂z 0 ∂z 0 0 r
S S S S
∂ ∂ ∂
ρwdS + ρw2 dS = τzz dS + τnz δS + ρgdS (A.18)
∂t 0 ∂z 0 ∂z 0 0
S S S S
∂ ∂ ∂
ρcp T dS + ρcp wT dS = − qz dS − q.nδS + ΦdS (A.19)
∂t 0 ∂z 0 ∂z 0 0
148 Annexe A. Formulation Intégrale du système
Simplifications
Si on considère maintenant le fluide incompressible (ρ = constante), le système unidimensionnel
( ∂φ
∂r 1) et une fibre de faible rayon de courbure ( ∂R
∂z 1), les expressions A.5 et A.11 deviennent :
1 ∂ru ∂w r ∂w
∇.v = + =0 ⇒ u=−
r ∂r ∂z 2 ∂z
σ ∂w
p = pa + − µ
R ∂z
σ
τrr − p = −pa −
R
σ
τθθ − p = −pa −
R
σ ∂w
τzz − p = −pa − + 3µ
R ∂z
2
µr ∂w
τrz = − 1
2 ∂z
∂R σ ∂σ
τnz δS = 2πR pa + − τa +
∂z R ∂z
2 2
2 2
∂w µr ∂ w
φ = 3µ −
∂z 4 ∂z 2
qn δs = 2πR(h(T − Ta ) + σ̃(T 4 − T ∗4 ))
En remplaçant ces expressions dans les équations A.16 à A.19, on obtient le système d’équations du
modèle qu’il nous faut résoudre :
∂S
+ Sw = 0 (A.20)
∂t
∂wS ∂w2 S ∂ σ ∂w
ρ + = S −pa − + 3µ
∂t ∂z ∂z R ∂z
∂R σ ∂σ
+2πR pa + − τa + + ρgS
∂z R ∂z
∂ ∂w ∂S ∂R ∂ σS
= 3Sµ − pa + 2πR pa + ρgS −
∂z ∂z ∂z ∂z ∂z R
∂R σ ∂σ
+2πR + − τa
∂z R ∂z
∂ ∂w √ √
= 3Sµ + σ πS − 2 πSτa + ρgS (A.21)
∂z ∂z
2 2
∂cp T S ∂cp T wS ∂qz S ∂w µπR4 ∂ 2 w
ρ + = − + 3µS − − 2πRqn
∂t ∂z ∂z ∂z 16 ∂z 2
2
√ ∂w
= λair πN u(Ta − T ) + 2 πSσ̃(T ∗4 − T 4 ) + 3µS
∂z
∂ S∂T ∂(Sqzr )
+λ −
∂z ∂z ∂z
(A.22)
Il est à noter qu’à l’intérieur de la fibre les échanges convectifs sont définis par la relation suivante :
∂(qzc S) ∂ ∂T
= −λ S (A.23)
∂z ∂z ∂z
Annexe B
Les équations de bilan sont discrétisées avec la méthode des différences finies. Nous utilisons un schéma
temporel purement implicite du premier ordre et une discrétisation spatiale du second ordre. Pour les
équations du mouvement et de l’énergie, nous avons choisi de prendre deux maillages distincts qui se
chevauchent (Fig.B.1) d’une demi-maille.
Fig. B.1 – Schématisation du maillage utilisé pour la discrétisation des équations bilan adimensionnées.
De manière à augmenter la précision temporelle, nous avons choisi de développer un schéma de Cranck-
Nicolson linéarisé, X n+1/2 = (1 − β)X n + βX n+1 avec β compris entre 0 et 1. L’erreur de calcul à chaque
itération donne la consistance du schéma.
On part des équations de bilan adimensionnées suivantes :
∂A ∂(W A)
+ = 0 (B.1)
∂t ∂x
√ ∂W ∂W ∂ ∂W √
−C3A + C2 A + C1 A + AW = MA +Ω A (B.2)
∂t ∂x ∂x ∂x
∂Θ ∂Θ √
A + AW = H(Θa − Θ) + E A(Θ∗4 − Θ4 ) (B.3)
∂t ∂x
Avec, comme conditions aux limites en x = 0, A0 = W0 = Θ0 = 1, et en x = L, WN b+1 = WL
La conservation de la masse
La discrétisation de l’équation de conservation de la masse (eq. B.1) se fait au pas d’espace i − 1/2 et
au pas de temps n + 1/2, ce qui donne :
150 Annexe B. Discrétisation des équations
pour 1 < i ≤ N b :
(An+1
i − Ani + An+1
i−1 − Ai−1 )
n
0 =
∆t
(1 − β)(Win Ani − Wi−1
n
Ani−1 ) + β(Win+1 An+1
i − Wi−1
n+1 n+1
Ai−1 )
+ 2 (B.4)
∆x
pour i = 1 :
(An+1
1 − An1 ) (1 − β)(W1n An1 − 1) + β(W1n+1 An+1
1 − 1)
0= +2 (B.5)
∆t ∆x
L’erreur est de
2 3 2 2 2 3
1 ∆t ∂ A 1 ∆x ∂ ∂A 1 ∆x ∂ WA
E1 = − 3
− 2
−
6 2 ∂t 2 2 ∂x ∂t 6 2 ∂x3
2 2
∆t ∂ ∂W A
− + o(∆x4 , ∆t4 , ∆x2 ∆t2 ) (B.6)
2 ∂t2 ∂x
La quantité de mouvement
La discrétisation de l’équation de conservation de la quantité de mouvement (eq. B.2) se fait au pas
d’espace i et au pas de temps n + 1/2, ce qui donne :
pour 1 < i < N b :
n
2(1 − β)(C2 An
i − C3A i ) + 2β(C2 An+1
i − C3An+1
i )
# &
n+1 n n n+1 n+1
n n+1 Wi − Win n n Wi+1 − Wi−1 n+1 n+1 Wi+1 − Wi−1
+C1 2((1 − β)Ai + βAi ) + (1 − β)Ai Wi + βAi Wi
∆t ∆x ∆x
n
(1 − β) n Wi+1 − Win Win − Wi−1
n n
Wi+1 − 2Win + Wi−1
n
= Mi+1 An
i+1
n
− Mi−1 An
i−1 + Min An
i
∆x ∆x ∆x ∆x
β n+1 n+1
+ Ω An i+1 − Ω Ani−1 + Ω Ai+1 − Ω Ai−1
∆x
n+1
n+1 n+1 Wi+1 − Win+1 n+1
n+1 n+1 Wi
n+1
− Wi−1 n+1
Wi+1 − 2Win+1 + Wi−1
n+1
+ Mi+1 Ai+1 − Mi−1 Ai−1 + Min+1 An+1
i
∆x ∆x ∆x
(B.7)
pour i = 1 :
n
2(1 − β)(C2 An 1 − C3A1 ) + 2β(C2 An+1
1 − C3An+1
1 )
n+1 n n n+1
n+1 W1 − W1 n W2 − 1 n+1 W2 −1
+C1 2((1 − β)An 1 + βA 1 ) + (1 − β)A n
1 W 1 + βA n+1
1 W 1
∆t ∆x ∆x
(1 − β) n
n n W2 − W1
n n
W1 − 1 n n
n n W2 − 2W1 + 1
= M2 A2 − + M1 A1 + Ω An 2 −Ω
∆x ∆x ∆x ∆x
β W2n+1 − W1n+1 W1n+1 − 1 n+1 n+1 W2
n+1
− 2W1n+1 + 1
+ M2n+1 An+1
2 − + Mi A 1 + Ω An+1
2 −Ω
∆x ∆x ∆x ∆x
(B.8)
pour i = N b :
n n+1
2(1 − β)(C2 An Nb − C3ANb ) + 2β(C2 An+1
Nb − C3ANb ) + C1
# n+1
&
n+1 n n
n n+1 WNb − WNb n n WL − WNb−1 n+1 n+1 WL − WNb−1
2((1 − β)ANb + βANb ) + (1 − β)ANb WNb + βANb WNb
∆t ∆x ∆x
n n n
(1 − β) n 7WL − 8WNb + WNb−1 WL − WNb−1
= MNb AnNb
n
+ (−4MNb−1 An n n
Nb−1 + 3MNb−2 ANb−2 )
2∆x 2∆x ∆x
β
+ 3Ω An Nb − 4Ω An Nb−1 + Ω AnNb−2 + 3Ω An+1
Nb − 4Ω An+1
Nb−1 + Ω An+1
Nb−2
∆x
n+1 n+1 n+1
n 7WL − 8WNb + WNb−1 WL − WNb−1
+ MNb An+1
Nb
n+1
+ (−4MNb−1 An+1
Nb−1 + 3M n+1
A n+1
Nb−2 Nb−2 ) (B.9)
2∆x 2∆x
151
L’erreur est de
2 3 2 3 2 2
1 ∆t ∂ W 1 ∆x ∂ WA ∆t ∂ ∂W
E1 = −C1A − C1AW − C1 AW
6 2 ∂t3 6 2 ∂x3 2 ∂t2 ∂x
2 2
∆t ∂ ∂W A
− + o(∆x4 , ∆t4 , ∆x2 ∆t2 ) (B.10)
2 ∂t2 ∂x
La conservation de l’énergie
La discrétisation de l’équation de l’énergie (eq. B.3) se fait au pas d’espace i et au pas de temps
n + 1/2, ce qui donne,
pour 1 < i < N b :
(B.11)
pour i = 1 :
(B.12)
pour i = N b :
N b − ΘN b
Θn+1 n
((1 − β)AnN b + βAn+1
Nb )
∆t
3Θ − 4ΘnN b1 + ΘnN b−2 n+1 3ΘN b − 4ΘN b1 + ΘN b−2
n n+1 n+1 n+1
+ (1 − β)(AnN b WNn b ) N b + β(An+1
N b WN b )
2∆x 2∆x
4∗ ∗4 4n+1
= H Θa − ((1 − β)ΘN b + βΘN b )) + E (1 − β) AN b (Θ − ΘN b ) + β AN b (Θ − ΘN b )
n n+1 n n4 n+1
(B.13)
L’erreur est de
2 2
2 ∂ Θ 2 ∂ Θ
EΘ = Θ3 (∆t) 2
− 2Θ3 (∆x) 2
+ o(∆x4 , ∆t4 , ∆x2 ∆t2 )
∂t # ∂x
2 2 2 3 2 2 &
∆x ∂ A ∂Θ 1 ∆t ∂ Θ ∆x ∂ ∂Θ
ErrΘ = − +A − −
2 ∂x2 ∂t 6 2 ∂t3 2 ∂x2 ∂t
# 2 2 &
2 2 √
1 ∆t ∂ Θ ∆x ∂ Θ
+H 2
+ 2
− E AEΘ + o(∆x4 , ∆t4 , ∆x2 ∆t2 )
6 2 ∂t 2 ∂x
152 Annexe B. Discrétisation des équations
Annexe C
Jacobienne du système
A partir des équations précédemment définies, on calcul les différents termes de la jacobienne, qui
sont pour 1 < i < N b :
∂F 1
= −2βWi−1
n
∂An+1
i−1
∂F 1
= 2βWin
∂An+1
i
∂F 1
n+1 = −2βAni−1
∂Wi−1
∂F 1
= 2βAni
∂Win+1
# &
∂F 2 ∆xΩ
= −β Mi−1
n
(Win − Wi−1
n
)+ n
∂An+1
i−1 2 Ai−1
# # &&
∂F 2 C2
= β Min (Wi+1
n
− 2Win + Wi−1
n
) − C1∆xWin (Wi+1
n
− Wi−1
n
) + ∆x2 2C3 − n
∂An+1
i Ai
# &
∂F 2 ∆xΩ
= β n
Mi+1 n
(Wi+1 − Win ) + n
∂An+1
i+1 2 Ai
∂F 2 n n
n+1 = β Mi−1 Ai−1 + Min Ani + C1∆xAni Win
∂Wi−1
∂F 2 n n
= −β Mi+1 Ai+1 + 2Min Ani + Mi−1
n
Ani−1 + C1∆xAni (Wi+1
n
− Wi−1
n
)
∂Win+1
∂F 2 n n
n+1 = β Mi+1 Ai+1 + Min Ani − C1∆xAni Win
∂Wi+1
∂F 2 ∂M
= −βAni−1 (Win − Wi−1
n
)
∂Θn+1
i−1
∂Θ
∂F 2 ∂M
= βAni (Wi+1
n
− 2Win + Wi−1
n+1
)
∂Θn+1
i
∂Θ
∂F 2 ∂M
= βAni+1 (Wi+1
n
− Win )
∂Θn+1
i+1
∂Θ
# &
4n 4∗
∂F 3 E(Θ − Θ )
= β Win (Θni+1 − Θni−1 ) + i
∆x
∂An+1
i 2 Ani
∂F 3
= βAni (Θni+1 − Θni−1 )
∂Win+1
∂F 3
= −βAni Win
∂Θn+1
i−1
∂F 3
= β(H + 4E Ani Θ3n
j )∆x
∂Θn+1
i
∂F 3
= βAni Win
∂Θn+1
i+1
157
en i = 1 on a :
∂F 1
= 2βW1n
∂An+1
1
∂F 1
= 2βAn1
∂W1n+1
# # &&
∂F 2 2 C2
= β M1n (W2n − 2W1n + 1) − C1∆xW1n (W2n − 1) + ∆x 2C3 − n
∂An+1
1 A1
# &
∂F 2 ∆xΩ
= β M2n (W2n − W1n ) + n
∂An+1
2 2 A1
∂F 2
= −β (M2n An2 + 2M1n An1 + 1 + C1∆xAn1 (W2n − 1))
∂W1n+1
∂F 2
= β (M2n An2 + M1n An1 − C1∆xAn1 W1n )
∂W2n+1
∂F 2 ∂M
= βAn1 (W2n − 2W1n + 1)
∂Θn+1
1 ∂Θ
∂F 2 ∂M
= βAn2 (W2n − W1n )
∂Θn+1
2 ∂Θ
# &
∂F 3 E(Θ4n 4∗
1 −Θ )
= β W1n (Θn2 − 1) + ∆x
∂An+1
1 2 An1
∂F 3
= βAn1 (Θn2 − 1)
∂W1n+1
∂F 3
= β(H + 4E An1 Θ3n
1 )∆x
∂Θn+1
1
∂F 3
= βAn1 W1n
∂Θn+1
2
158 Annexe C. Jacobienne du système
en i = N b on a :
∂F 1
= −2βWNn b−1
∂An+1
N b−1
∂F 1
= 2βWNn b
∂An+1
Nb
∂F 1
= −2βAnN b−1
∂WNn+1
b−1
∂F 1
= 2βAnN b
∂WNn+1
b
# &
∂F 2 ∆xΩ 1 n
= β n + MN b−2 (WL − WNn b−1 )
∂An+1
N b−2 AN b−2 2
# &
∂F 2 4∆xΩ
−β 2MN b−1 (WL − WN b−1 ) + n
n n
=
∂An+1
N b−1 AN b−1
# &
∂F 2 7W L − 8W n
+ W n
3∆xΩ
= β MN n
b
N b N b−1
+ n − C1∆xWNn b (WL − WNn b−1 )
∂An+1
Nb
2 AN b
# &
C2
+β∆x2 2C3 − n
AN b
b−1 AN b−1 − MN b−2 AN b−2 + MN b AN b
n n n n n n
∂F 2 4MN
= β + C1∆xAnN b WNn b
∂WNn+1
b−1
2
∂F 2
= −β 4MN
n
b AN b + C1∆xAN b (WL − WN b−1 )
n n n
∂WNn+1
b
∂F 2 1 n ∂M
= βA (WL − WNn b−1 )
∂Θn+1
N b−2
2 N b−2 ∂Θ
∂F 2 ∂M
= −2βAnN b−1 (WL − WNn b−1 )
∂Θn+1
N b−1
∂Θ
18
valeur analytique
16
valeur numérique
14
12
Vitesse[m/s]
10
0
0,00 0,01 0,02 0,03 0,04 0,05
z[m]
Fig. D.1 – Profil de vitesse avec dissipation visqueuse dominante : M odf ib/solution analytique.
12
Vitesse[m/s]
10
0
0,00 0,01 0,02 0,03 0,04 0,05
z[m]
Fig. D.2 – Idem figure précédente avec forces d’inertie incluses.
160 Annexe D. Validations pour un jet isotherme
Bibliographie
[1] Milton Abramowitz and Irene A. Stegun. Handbook of mathematical Functions. Dover inc., NY, 1965.
[2] Y. Aizu, J. Domnick, F. Durst, G. Gréhan, F. Onofri, H.H. Qiu, M. Sommerfeld, T.-H. Xu, and M. Ziema.
A new generation of Phase Doppler instruments for particle velocity, size and concentration measurements.
Part. Part. Syst. Charat., 2 :43–54, 1994.
[3] H. E. Allbrecht, N. Damaschke, M. Borys, and C. Tropea. Laser Doppler and Phase Doppler Measurement
Techniques. Berlin Springer, 2003.
[4] F.R. Archibald and A. G. Emslie. The vibration of a string having a uniform motion along its length. Appl.
Mech., 25 :347–348, 1954.
[5] W.D. Bachalo and M.J. Houser. Phase/Doppler spray analyzer for simultaneous measurements of drop size
and velocity distributions. Optical Engineering, 23 :583–590, 1984.
[6] S. Bateson. Critical study of the optical and mechanical properties of glass fibers. Journal of Applied
Physics, 29(1) :13–21, 1958.
[7] S.E. Bechtel, K.D. Bolinger, J.Z. Cao, and M.G. Forest. Torsional effects in high-order viscoelastic thin-
filament models. SIAM J. Appl. Math., 55(1) :58–99, 1995.
[8] R. Beerkens. Chemical equilibrium reaction as driving forces for growth of gas bubbles during refining.
Glastechnische Berichte, 63 :222–242, 1990.
[9] A. Bergonzoni and A.J. DiCresce. The phenomenon of draw resonance in polymeric melts parti qualitative
view. Polymer Engineering and Science, pages 45–49, 1966.
[10] C.F. Bohren and D. R. Huffman. Absorption and Scattering of Light by Small Particles. John Wiley &
Sons, NY, 1998.
[11] A.R. Bunsell, editor. Fiber Reinforcements for composite materials. ELSEVIER, 1988.
[12] W.F. Busse. Mechanical structure in polymer melts . (I) mesurements of melt strength and elasticity.
Journal of polymer science, 5 :1249–1259, 1967.
[13] J.C. Chang and M.M. Denn. An experimental study of isothermal spinning of a newtonian and viscoelastic
liquid. Journal of Non-Newtonian Fluid Mechanics, 5 :369–385, 1979.
[14] J.C. Chang, M.M. Denn, and F.T. Geyling. Effects of inertia, surface tension, and gravity on the stability
of isothermal drawing of newtonian fluids. Ind. Eng. Chem. Fundam., 20 :147–149, 1981.
[15] J.C. Chang, M.M Morton, and S. Kase. Dynamic simulation of low-speed melt spinning. Ind. Eng. Chem.
Fundam., 21 :13–17, 1982.
[16] S-C. Chang, J. Jin J-M. Jin, and S. Zhang. Computations of Special Functions. Published by John Wiley
& Sons, 1996.
[17] P. Chylek, J.T. Kiehl, M.K.W. KO, and A. Ashkin. Surface waves in light scattering by spherical and non
spherical particles, pages 153–164. Plenum Press, NY, in lights scattering by irregularly shaped particles
(ed. d.w. schuerman) edition, 1980.
[18] J.M. Corpus and P.K. Gupta. Diameter dependence of the refrective index of melt-drawn glass fibers. J.
Am. Ceram. Soc., 76(5) :1390–1392, 1993.
[19] C. Clementin de Leusse. Etude de fibres de verre E. Caractérisation expérimentale par spectrométrie de
Brillouin et test de modèle de fribrage. PhD thesis, Université de Monpellier II, Jul. 2000.
[20] Y. Demay and F. Agassant. Experimental study of draw resonance in fiber spinning. Journal of Non-
Newtonian Fluid Mechanics, 18 :187–198, 1985.
[21] P. Demont. Etude de base fibrage : modélisation du fibrage, et application à la conception des filières.
rapport d’avancement, DDI R12 1985.
[22] M.M. Denn. Continous drawing of liquids to form fibers. Ann. Rev. Fluid Mech., 12 :365–387, 1980.
[23] A. Déroche. Conception d’un interféromètre phase-Doppler. Rapport de DEA-IUSTI-UMR CNRS n◦ 6595,
Univ. de Provence, Marseille, France, 1999.
[24] G.J. Donnelly and C.B. Weinberger. Stability of isothermal fiber spinning of a newtonian fluid. Ind. and
Eng. Chem. Fundamental, 14(4) :334–337, 1975.
161
162 Bibliographie
[25] R.J. Donnerly and W. Glaberson. Experiments on the capillary instability of a liquid jet. Proc. Roy. Soc.
A, 290 :547–556, 1966.
[26] F. Durst, A. Melling, and J.H. Whitelaw. Principles and practice of laser-Doppler anemometry. Academic
Press, London, 1981.
[27] F. Durst and M. Zaré. Laser Doppler measurements in two-phase flows. In Proceedings of LDA-Symposium,
Copenhagen, pages 403–429, 1975.
[28] V.M. Entov and A.L. Yarin. The dynamics of thin liquid jets in air. J. Fluid Mech., 140 :91–111, 1984.
[29] R.J. Fisher and M.M. Denn. Finite-amplitude stability and draw resonance in isothermal melt spinning.
Chemical Engeneering Science, 30 :1129–1134, 1975.
[30] R.J. Fisher and M.M. Denn. A theory of isothermal melt spinning and draw resonance. AIChE Journal,
22(2) :236–246, 1976.
[31] D. Gelder. The stability of fiber drawing process. Ind. Eng. Chem. Fudam., 10 :534–535, 1971.
[32] F.T. Geyling. Basic fluid-dynamic considerations in the drawing of optical fibers. The Bell Systeme Technical
Journal, 55(8) :1011–1056, 1976.
[33] L.R. Glicksman. The colling of glass fibres. Glass Technology, 9 :131–138, 1968.
[34] L.R. Glicksman. The dynamics of heated free jet of variable viscosity at low Reynolds numbers. Journal
of Basic of Basic Engineering, 90 :343–354, 1968.
[35] P. Gospodinov and A.L. Yarin. Draw resonance of optical microcapillaries in non-isothermal drawing. Int.
J. Multiphase Flow, 23(3) :967–976, 1997.
[36] G. Gouesbet and G. Gréhan. Interaction between shaped beams and an infinite cylinder, including a
discussion of gaussian beams. Part. Part. Syst. Charact., 11 :299–308, 1994.
[37] G. Gouesbet, B. Maheu, and G. Gréhan. Light scattering from a sphere arbitrarily located in a Gaussian
beam, using a Bromwich formulation. J.O.S.A. A., 5(9) :1427–1443, 1988.
[38] G. Gréhan, G. Gouesbet, A. Naqwi, and F. Durst. Particle trajectory effects in Phase Doppler systems :
Computations and experiments. Part. Part. Syst. Charact., 10 :332–338, 1993.
[39] R.M. Griffith. Velocity, temperature and concentration distributions during fiber spinning. I&EC Funda-
mental, 3(3) :245–250, 1964.
[40] A.K. Gupta, K.R Srinivasan, and P Krishna Kumar. Glass fiber reinforced polypropylene/epdm blends. ii
mechanical properties and morphology. Journal of Applied Polymer Science, 43 :451–462, 1991.
[41] G. Gupta. Fiber Reinforcements for composite materials, chapter 2, pages 19–71. A.R. Bunsell (edited by
ELSEVIER), 1988.
[42] G. Gupta and W.W. Schultz. Non-isothermal flows of newtonian slender glass fibers. Int. J. Non-Linear
Mechanics, 33(1) :151–163, 1998.
[43] G.K. Gupta, W.W. Schultz, E.M. Arruda, and X. Lu. Nonisothermal model of glass fiber drawing stability.
Rheologica Acta, 35(6) :584–596, 1996.
[44] I. Gurwich, N. Shiloah, and M. Kleiman. The recursive algorithm for electromagnatic scattering by tilted
infinite circular multilayered cylinder. J. of Quantitative Spectroscopy and Rad. Transfer, 63 :217–229, 1999.
[45] J. Huang and P.K. Gupta. Enthalpy relaxation in thin glass fibers. Journal of Non-Cristalline Solids,
151 :175–181, 1992.
[46] H.M. Hulburt and S. Katz. Some problems in particle technology. Chemical Engeneering Science, 19 :555–
574, 1964.
[47] S. Kase and T. Matsuo. Studies on melt spinning. (I). fundamental equations of the dynamics of melt
spinning. Journal of Polymer Science : Part A, 3 :2541–2554, 1965.
[48] S. Kase and T. Matsuo. Studies on melt spinning. (II). steady-state and transient solution of fundamental
equations compared with experimental results. Journal of Applied Polymer Science, 11 :251–287, 1967.
[49] D. Lebrun, A. Belaid, C. Ozkul, K-F. Ren, and G. Gréhan. Opt. Eng., 35 :946–950, 1996.
[50] D. Lemmonier. Modélisation dynamique et thermique du fibrage du verre. rapport ENSMA, 1996.
[51] A. Lenoble. Etude expérimentale et modélisation physique des modes d’instabilité d’un procédé de fibrage
du verre, 2000. stage de D.E.A. de l’Université Paris XI-Orsay - IUSTI- UMR CNRS 6595-Univ. Provence,
Marseille, France.
[52] A. Lenoble, F. Onofri, and S. Radev. Modélisation physique du procédé de fibrage des fibres de renforcement.
In Congrès Français de Thermique, Vittel (France), 2002.
[53] X.L. Luo and R.I. Tanner. A computeur study of film blowing. Polymer Engeneering and Science,
25(10) :620–629, 1985.
[54] L.J. Lynch and N. Thomas. Optical diffraction profiles of single fibers. Textile reseach journal, 41 :568–572,
1971.
Bibliographie 163
[55] U. Manasse, T. Wriedt, and K. Bauckhage. Reconstruction of real size distributions hidden in Phase Doppler
Anemometry results obtained from droplets of inhomogeneous liquids. Part. Part. Syst. Charat., 11 :84–90,
1994.
[56] M.A. Matovitch and J.R.A. Pearson. Spinning a molten threadline - stability. I&EC Fundamentals,
8(4) :605–609, 1969.
[57] S. Middleman. Modeling Axisymetric Flows -Dynamics of Films, Jets, and Drops, chapter 1, pages 2–5.
Academic Press, 1995.
[58] H. Mignon. Anémométrie Phase Doppler et Particules Non Sphériques : Cas des Cylindres et des Ellipsoı̈des.
PhD thesis, La Faculté des Sciences de l’Université de Rouen, janvier 1997.
[59] H. Mignon, G. Gréhan, G. Gouesbet, T. H. Xu, and C. Tropea. Mesurement of cylindrical particles with
phase Doppler anemometry. Applied Optics, 35(25) :5180–5190, 1996.
[60] A. Napolitano and E.G. Hawkins. Viscosity of a standard soda-lime-silica glass. A. Physics and Chemistry,
68A(5) :439–448, 1964.
[61] A. Naqwi. Sizing of irregular particles using a phase doppler system. Part. Part. Syst. Charact., 8 :343–349,
1996.
[62] A. Naqwi and F. Durst. Analysis of the laser light scattering interoferometric devices for the in-line diagnostic
of moving particles. Appl. Opt., 32 :4003–4018, 1993.
[63] L. Nẽmec. The behaviour of bubbles in glass melts. part 1 bubble size controlled by diffusion. Glass
Technology, 21(3) :134–138, 1980.
[64] L. Nẽmec and S. Laurentova. The main technological characteristics of glass melting zone with respect to
the course of the melting process : (III) the main characteristics holding till complete sand dissolution and
refining. Silicaty, 33 :105–122, 1989.
[65] R.E. Nickell, R.I. Tanner, and B. Caswell. The solution of viscous incompressible jet and free-surface flows
using finite element methods. J. Fluid Mech., 65(part 1) :189–206, 1974.
[66] F. Onofri. Prise en compte des dimensions finies des faisceaux d’éclairage en granulométrie optique :
Anémométrie Phase Doppler, Applications aux Diagnostics des Milieux Diphasiques. PhD thesis, Université
de Rouen, URA CNRS 230 (FRANCE), 1995.
[67] F. Onofri. Etude numérique et expérimentale de la sensibilité de l’Anéomométrie phase Doppler à l’état de
surface des particules détectées. In 7ème Congrès Francophone de Vélocimétrie Laser, 19-22 Sept. Marseille
(France), pages 335–342, 2000.
[68] F. Onofri. Calcul des séries de Debye pour la sphère et le cylindre homogènes. Rapport interne, UMR
CNRS 6595-IUSTI, Univ. Provence, Marseille, France., Nov. 2003.
[69] F. Onofri, L. Bergougnoux, J.-L. Firpo, and J. Misguich-Ripault. Size, velocity, and concentration in
suspension measurements of spherical droplets and cylindrical jets. Applied Optics, 38(n◦ 21) :4681–4690,
1999.
[70] F. Onofri, D. Blondel, G. Gréhan, and G. Gouesbet. On the optical diagnosis and sizing of spherical coated
and multilayered particles with phase-Dopple anemometry. Particle &Particle Systems Characterization,
13 :104–111, 1996.
[71] F. Onofri, G. Gréhan G., and G. Gouesbet. Electromagnetic scattering from a multilayered sphere located
in an arbitrary beam. Appl. Opt., 34(30) :7113–7124, 1995.
[72] F. Onofri, T. Girasole, G. Gréhan, G. Gouesbet, G. Brenn, J. Domnick, T.-H. Xu, and C. Tropea. Phase-
Doppler anemometry with the dual burst technique for measurement of refractive index and absorption
coefficient simultaneously with size and velocity. Particle and Particle Systems Characterization, vol. 13 :pp.
112–124, 1996.
[73] F. Onofri, G. Gréhan, and G. Gouesbet. Electromagnetic scattering from a multilayered sphere located in
an arbitrary beam. Applied Optics, 34(30) :7113–7124, 1995.
[74] F. Onofri and A. Lenoble. Sizing of single fibers under torsional stress. In 6th International Congress on
Optical Particle Characterisation, Brighton, U.K., 2001. paper 7.7.
[75] F. Onofri, A. Lenoble, H. Bultinck, P.-H. Guering, and N. Marsault. Interferometric sizing og single-
axis birefingent glass fiber. In 10 th International Symposium on Applications of laser Technique to Fluid
Mechanics - Lisbon Portugal, 2002.
[76] F. Onofri, A. Lenoble, H. Bultynck H., and P.H. Guering. High-resolution diffractometry for the on-line
sizing of small transparent fibres. Opt. Com., 234 :184, 2004.
[77] F. Onofri, A. Lenoble, and N. Marsault H. Bultynck, P.H. Guering. Interferometric sizing of single-axis
birefringent glass fiber. Part. Part. Syst. Charact., 20 :171–182, 2003.
[78] F. Onofri, A. Lenoble, and S. Radev. Superimposed non interfering probes to extend the phase Doppler
anemometry size range. In 10 th International Symposium on Applications of laser Technique to Fluid
Mechanics - Lisbon Portugal, 2000. paper 20.1, 8 pages.
164 Bibliographie
[79] F. Onofri, A. Lenoble, and S. Radev. Superimposed non interfering probes to extend the phase Doppler
anemometry capabilities. Applied Optics, 41(18) :3590–1600, 2002.
[80] F. Onofri, A. Lenoble, S. Radev, and H. Bultinck. Contrôle temps réel du diamètre de fibre de verre en
cours de production. In 3 ème Colloque Francophone de Méthodes et Techniques Optiques pour l’Industrie,
pages 63–68, Saint Aubin de Médoc, France, 2002. paper.
[81] F. Onofri, A. Lenoble, S. Radev, and P-H. Guering. Optical measurement of the drawing tension of small
glass fibres. Submitted for publication to Meas. Sci. Tech., 2003.
[82] F. Onofri, A. Lenoble, and Radev S. 3-coherent beams phase Doppler and laser Doppler techniques. sub-
mitted for publication to Opt. Com., 2004.
[83] F. Onofri, H. Mignon, G. Gouesbet, and G. Gréhan. On the extension of phase-Doppler anemometry to the
sizing of spherical multilayered particles and cylindrical particles. In 4 th Optical Particle Sizing Conference
Partec, pages 275–284, Nüremberg, Germany, 1995.
[84] C. Ozkul, D. Lebrun, D. Allano, A. Abdelghani-Idrissi, and A. Leduc. Opt. Eng., 1991 :1855–1861, 1991.
[85] T.C. Papanastasiou, C.W. Macosko, L.E. Scriven, and Z. Chen. Fiber spinning of vicoelestic liquid. AIChe
Journal, 33(5) :834–842, 1987.
[86] E.A. Patterson, W. Ji, and Z.F. Wang. On image analysis for birefringence measurements in photoelasticity.
Optics and Lasers in Engineering, 28 :17–36, 1997.
[87] J.R.A. Pearson. Instability in non-newtoniean flow. Ann. Rev. Fl. Mech., 8 :163–181, 1976.
[88] J.R.A. Pearson, Y.T. Sahh, and R.D. Mhaskar. On the stability of fiber spinning of freezing fluids. Ind.
Eng. Chem. Fundam., 15 :31–37, 1976.
[89] J.R.A. Pearson and Y.T. Shah. On the stability of isothermal and non-isothermal fiber spinning of power-law
fluids. Ind. Eng. Chem. Fundam., 13(2) :134, 1974.
[90] W.H. Press, B.P. Flannery, S.A. Teukolsky, and W.T. Vetterling. Numerical Recipes : The art of scientific
computing. Cambridge University Press, 1986.
[91] R. Rajamanickam, S.M. Hansen, and S. Jayaraman. Studies on fiber-process-structure-property relationships
in air-jet spinning part ii : Model development. Journal of Textile Institute, 89(2) :243–265, 1998.
[92] J. Rheims, H. Dahl, T. Wriedt, and K. Bauckage. Sizing coated spheres with Phase Doppler Anemometry.
In 4th Optical particle sizing conference-Partec95, Nurnberg, Germany, pages 409–418, 1995.
[93] T. Roi, O. Seidel, G. Nölle, and D. Höhne. Modelling of bubble behaviour in melts. In 2nd int seminar on
mathematical simulation in the glass melting, pages 30–56, 1993.
[94] S.A. Schaub, A. A. Naqwi, and F.L. Harding. Design of a phase/Doppler light-scattering system for me-
surement of small-diameter glass fibers during fiberglass manufacturing. Applied Optics, 37(3) :573–585,
1998.
[95] Y.T. Shah and J.R.A. Pearson. On the stability of nonisothermal fiber spinning. Ind. Eng. Chem. Fundam.,
11 :145–149, 1972.
[96] Y.T. Shah and J.R.A. Pearson. On the stability of nonisothermal fiber spinning - general case. Ind. Eng.
Chem. Fundam., 11 :150–153, 1972.
[97] Y.T. Shah and J.R.A. Pearson. Stability of fiber spinning of power-law fluids. Polymer Engeneering and
Science, 12(3) :219–223, 1972.
[98] M. Song, K.S. Ball, and T.L. Bergman. A model for radiative cooling of a semitransparent molten glass jet.
Journal of Heat Transfert, 120 :931–938, 1998.
[99] R.L. Stover, C.W. Tobias, and M.M. Denn. Bubble coalescence dynamics. AIChE Journal, 43(10) :2385–
2392, 1997.
[100] H. Tadrist. Les outils du thermicien. UMR CNRS 6595-Univ. Provence, Marseille, France., 1997.
[101] C. Toutiri. Modélisation dynamique et thermique d’une fibre de verre en cours d’étirage. PhD thesis,
LET-UMR CNRS 6608-ENSMA-Université de Poitiers, Poitiers, France., janvier 1999.
[102] C. Tropea, T-H. Xu, F. Onofri, G. Grehan, and P. Haugen. Dual-mode phase Doppler anemometer. Part.
and Part. Syst. Charact., 13 :165–170, 1996.
[103] A. Ungan. Three dimension numerical simulation of fining process in glass melting furnaces : Mathematical
formulation. Glastechnische Berichte, 63 :19–28, 1990.
[104] A. Ungan and B. Balkanli. Numerical simulation of refining bubble behavior in glass melting furnace.
ASME, 3 :15–21, 1992.
[105] H.C van de Hulst. Light Scattering by Small Particles. Dover Publications, Inc. New York, 1957.
[106] R. Viskanta. Review of three-dimensional mathematical modeling of glass melting. Journal of Non Crys-
talline Solids, 177 :347–362, 1994.
[107] J.L. White and J.F. Roman. Instabilities and failure in elongational flow and melt spinning of fibers. Journal
of Applied Polymer Science, 22 :3057–3074, 1978.
Bibliographie 165
[108] Z. Wu and L. Guo. Electromagnetic scattering from a multilayered cylinder arbitrarily located in a gaussian
beam, a new recursive algorithms. Progress in Electromagnatics Researc, 18 :317–333, 1998.
[109] Z-S. Wu and Y.P. Wang. Electromagnetic scattering from multilayered sphere : recursive algorithms. Radio
Science, 26(6) :1393–1401, 1991.
[110] Z.S. Wu, L.X Guo, K.F. Ren, and G. Gouesbet et G. Gréhan. Improved algorithms for electromagnetic
scattering of plane waves and shaped beams by multilayered spheres. Appl. Opt., 36(21) :5188–5198, 1997.
[111] T.H Xu and C. Tropea. Improving the performance of two-component Phase Doppler Anemometry. Meas.
Sci. Techn., 5 :969–975, 1994.
[112] A. Yarin, Vl. Rusinov, P. Gospodinov, and S. Radev. Quasi one-dimensional model of drawing of glass
microcapillaries and approximate solutions. Journal of Theoretical and Applied Mechanics, (3) :55–62,
1989.
[113] A.L. Yarin, P. Gospodinov, O. Gottieb, and M.D. Graham. Newtonnian glass fiber drawing : Chaotic
variation of the cross section radius. Physics of Fluids, 11(11) :3201–3208, 1999.
[114] A. Ziabicki. Fundamentals of fibre formation. New-York : Wiley, 1976.
Titre : Caractérisation optique et étude de la stabilité d’un procédé de fibrage du verre.
Résumé : Dans ce travail de thèse différents outils théoriques, numériques et expérimentaux, ont été
développés afin d’étudier les conditions de stabilité d’un procédé de fibrage du verre-E.
Un interféromètre phase Doppler et un diffractomètre haute résolution ont été mis au point pour
mesurer, en temps réel, l’évolution du diamètre d’une fibre (D = 5 − 42µm) en sortie de filière (vitesses :
Vf = 5 − 65m/s). L’interféromètre permet également de mesurer la tension de fibrage. Les résolutions
obtenues avec ces systèmes (respectivement de 0.35µm et 0.02µm) l’ont été grâce au développement d’un
modèle rigoureux de diffusion de la lumière par un cylindre multicouche et en prenant en compte des
effets optiques propres à ce procédé : biréfringence uni-axe induite par la tension de fibrage, dépendance
de l’indice de réfraction avec les conditions de refroidissement, fibres creuses...
Le procédé de fibrage a été modélisé par l’étirage d’un jet visqueux à haute température, axisymétrique
et 1D. Ce modèle hydrodynamique, physique, permet notamment de prédire l’évolution du profil axiale
de température du jet, son profil de contraction... pour des régimes stationnaires et non stationnaires.
Les travaux expérimentaux réalisés sur une filière mono-téton ont montré que la stabilité du procédé
est maximale pour des températures de fibrage de T0 = 1145 − 1175◦ C et des fibres dont le diamètre est
inférieur à D ≈ 15µm. Cependant, même dans ces conditions, le diamètre de la fibre produite fluctue
périodiquement avec une amplitude moyenne de l’ordre de σD /D̄ ≈ 1.8%. La fréquence de ces oscillations
croı̂t avec la température de fibrage : νosc = 0.5 → 0.9Hz pour T0 = 1145 → 1250◦C. Le taux moyen des
fluctuations du diamètre est de l’ordre de dσD /dt ≈ 1.9µm/s et dσD /dL ≈ 0.07µm/m (pour Vf = 20m/s).
Mots clefs : Fibres de renforcement, hydrodynamique, instabilités, déchets verriers, tension, granu-
lométrie, interférométrie, diffractométrie, théorie de Lorenz-Mie, cylindre multicouche, biréfringence uni-
axe.
Title : Optical characterization and study of the stability of a glass fibre drawing process
Abstract : In this PhD work several theoretical, numerical and experimental tools are developed to
study the stability of a E-glass fibre drawing process.
A phase Doppler interferometer and a high resolution diffractometer have been developed to determine
the diameter of a fibre (D = 5 − 42µm) under drawing conditions (VelocityVf = 5 − 65m/s). These
systems allow high resolutions measurements (0.35µm and 0.02µm respectively) thanks to the derivation
of a rigorous model for the light scattering by a multilayered cylinder. This last model takes also into
account various particular properties of these fibres : single-axis birefringence induced by the drawing
stress, the refractive index dependence with the cooling rate, hollow fibres...
The process is modelled by the drawing of a one dimensional and axi symmetrical viscous jet at
high temperature. This model allows predicting various quantities such as the fibre’s temperature and
contraction profiles... under stationary or non stationary conditions.
Preliminary experimental results, performed on a single fibre drawing bench, have shown that the
stability of the drawing process is maximum for temperatures in the range T0 = 1145 − 1175◦C and for
fibres with a diameter bellow D ≈ 15µm. Nevertheless, even under such conditions, the fibre’s diameter
fluctuates periodically with a mean amplitude of σD /D̄ ≈ 1.8%. The frequency of theses oscillations in-
creases νosc = 0.5 → 0.9Hz with the temperature T0 = 1145 → 1250◦C. The mean rate of the diameter’s
fluctuations is about dσD /dt ≈ 1.9µm/s and dσD /dL ≈ 0.07µm/m (for Vf = 20m/s).
Key words : Reinforcement glass fibres, hydrodynamic, instabilities, glass wastes, sizing, tension, inter-
ferometry, diffractometry, Lorenz-Mie theory, multilayered cylinder, single-axis birefringence.