Thèse APPROCHE DÉTERMINISTE DU SÉCHAGE DES AVIVÉS ...

Chapitre III
Figure III.10: Représentation de la méthode du gradient conjugé dans un espace de dimension 2.
(a) Descente en plusieurs itérations vers le minimum de f ; (b) Zoom sur une itération : recherche suivant
une direction perpendiculaire au contour f ( X k ) du minimum local X k + 1 , à partir duquel une nouvelle
direction de recherche orthogonale est construite.[Numerical Recipes In Fortran 77 (1986-1992)].
2.2 Ajout du module mécanique au code TransPore
A partir du champ d’humidité et de température dans l’épaisseur d’une planche, un module, écrit en
Fortran, assemble la matrice A et le vecteur B et résout le système. Ce programme est appelé par
Transpore (programme principal) à chaque pas de temps (Fig.III.11).
Le temps de calcul du code a été allongé par l’ajout de ce programme. Néanmoins, l’ensemble
conserve son aspect interactif grâce à des durées de calculs faibles, restant de l’ordre de quelques
secondes. Pour l’illustrer, le code TransPore simule le séchage complet d’une planche de 50 mm
d’épaisseur, discrétisée en 100 contrôles de volume en 2.5 secondes (Pentium IV, 1,3 GHz). L’ajout
du sous-programme de mécanique rallonge ce temps de calcul d’environ 2 secondes. La méthode du
gradient biconjugué converge dans ce cas après seulement 10 itérations.
Lorsque les conditions de transferts, de part et d’autre de la planche sont symétriques, seul le
séchage d’une demi-épaisseur de planche est simulé, afin de réduire le nombre de volumes de
contrôle, et ainsi réduire le temps de calcul du code. Les conditions aux limites sont modifiées :
- inexistence des transferts externes sur l’une des faces,
- rayon de courbure imposé et tendant vers l’infini.
Page-110-