FR 2853404 A1 20041008 FR 0304284 A 20030407 i La présente invention est relative à un procédé de détermination de paramètres de fonctionnement d'une installation de refroidissement thermique d'articles. L'invention s'applique notamment aux installations de surgélation d'articles alimentaires. Les installations connues de surgélation comportent, par exemple, une enceinte ou tunnel, de surgélation traversé de part en part par un convoyeur à bandes sur lequel sont déposés des articles à congeler, le convoyeur circulant en continu ou séquentiellement au travers du tunnel de surgélation. Une enceinte cryogénique utilise un fluide inerte à basse température qui échange de la chaleur directement par contact avec les articles à surgeler. De manière classique, une enceinte cryogénique utilise soit la neige carbonique (-800C), l'air liquide ou l'azote liquide (-1960C) comme vecteur de froid. La neige carbonique permet le transport de produits frais ou surgelés sans 15 crainte de rupture de la chaîne de froid. L'azote et l'air liquide permettent soit la surgélation individualisée de produits alimentaires, soit le durcissement de produits fragiles, déformables ou collants (du type crème glacée...). Les paramètres de fonctionnement des installations forment des recettes créées de manière expérimentale. Une recette mémorise les paramètres de 20 réglage d'une installation pour une production donnée. Il n'existe pas aujourd'hui de capteur capable de mesurer en continu et sans contact, la température interne d'un article, de sorte que la détermination des recettes passe par des procédés de détermination de paramètres de fonctionnement qui sont ensuite ajustés à l'aide de tests réalisés en situation, les25 quels tests sont le plus souvent destructifs. Les procédés de détermination de recettes comprennent généralement une étape de détermination d'une consigne portant sur la température de sortie des articles suivie d'une étape de détermination de paramètres de fonctionnement initiaux puis la réitération d'un cycle de test comprenant une étape de 30 prédiction de la température de sortie des articles, une étape de comparaison de la température prédite avec la consigne et, en cas de différence, une étape de modification des paramètres de fonctionnement. Le cycle de test est renouvelé jusqu'à ce que les paramètres de fonctionnement permettent d'obtenir une température prédite sensiblement proche de la température de consigne. Si l'on examine le système constitué de l'enceinte et de la charge d'ar5 ticles, plusieurs paramètres peuvent influer sur la température de l'article en sortie: le débit de production qui, pour un taux de chargement donné, implique une variation du temps de séjour dans l'enceinte, le débit du fluide qui agit sur le profil des températures, la température d'entrée des articles, le profil convectif de l'enceinte, et le taux de chargement. Le système est donc un système multivariables et les procédé de détermination de paramètres existants, et notamment l'étape de prédiction de la température des articles en sortie, ne prennent pas en compte ces éléments. Dans les étapes de prédiction de l'état de la technique, afin de traiter des systèmes monovariables, on a été amené à considérer le profil convectif et le 15 taux de chargement comme constants, et à fixer le débit de production, la température des articles en entrée de l'enceinte de surgélation ainsi que d'autres paramètres de fonctionnement des installations. De ce fait, les recettes déterminées par les procédés de détermination existants sont relativement imprécises et nécessitent la mise en oeuvre d'une 20 production suivie de tests destructifs sur les articles. La présente invention vise à remédier à ce problème en définissant un procédé de détermination de paramètres de fonctionnement, précis et aisément mis en couvre. A cet effet, l'invention a pour objet un procédé de détermination des 25 paramètres de fonctionnement d'une installation de refroidissement thermique d'articles, comprenant une enceinte traversée par lesdits articles d'une entrée à une sortie et utilisant un fluide de refroidissement, le procédé comprenant: - une étape de détermination d'une consigne de température d'articles en sortie de ladite enceinte; - une étape de détermination de paramètres de fonctionnement initiaux pour ladite installation; et - un cycle de test des paramètres de fonctionnement comprenant: - une étape de prédiction de la température d'articles en sortie de ladite enceinte; - une étape de comparaison de la température de consigne avec la température prédite; et - si ladite étape de comparaison révèle un écart supérieur à un seuil prédéterminé, une étape de modifications de paramètres de fonctionnement de ladite installation et une réitération du cycle de test, caractérisé en ce que ladite étape de prédiction est réalisée à partir de paramètres de fonctionnement de ladite enceinte, de caractéristiques thermodynamiques et physiques de ladite enceinte et de caractéristiques thermodynamiques et physiques desdits articles. Suivant d'autres caractéristiques: - ladite étape de prédiction comporte une étape de prédiction du comportement de ladite enceinte fondée sur la résolution de bilans thermiques sur des tranches élémentaires du volume de ladite enceinte, réalisée au moins à partir de caractéristiques thermodynamiques dudit fluide de refroidissement et de 15 caractéristiques thermodynamiques et physiques de ladite enceinte; - ladite étape de prédiction du comportement de ladite enceinte est réalisée en outre, à partir de paramètres de fonctionnement de ladite installation; - lesdits paramètres de fonctionnement de ladite installation représentent au moins l'un des éléments choisis dans le groupe consistant en: - la vitesse d'un convoyeur de transport desdits articles au travers de ladite enceinte; le taux de chargement; et - la ventilation de l'atmosphère de ladite enceinte; - ladite étape de prédiction comporte une étape de prédiction du com25 portement desdits articles fondée sur la résolution de l'équation de conservation de la chaleur discrétisée et appliquée à un réseau de points spatiaux et temporels constituant un maillage desdits articles, réalisée au moins à partir de caractéristiques thermodynamiques et physiques desdits articles; - ladite étape de prédiction du comportement desdits articles est réali30 sée en outre, à partir de paramètres de fonctionnement de ladite installation; - lesdits paramètres de fonctionnement de ladite installation comprennent la température desdits articles en entrée de ladite enceinte; - ladite étape de prédiction du comportement desdits articles est optimisée par des calculs de modification dudit maillage desdits articles selon des suites mathématiques; - ladite étape de prédiction du comportement desdits articles est opti5 misée par suppression des calculs de prédiction pour des points spatiaux et temporels dudit maillage desdits articles pour lesquels les variations d'enthalpie sont inférieures à un seuil prédéterminé; - ladite étape de prédiction de la température desdits articles en sortie de ladite enceinte est fondée sur ladite étape de prédiction du comportement de 10 ladite enceinte ainsi que sur ladite étape de prédiction du comportement desdits articles; - ladite étape de modification des paramètres de fonctionnement comprend une étape de modification manuelle d'au moins une partie des paramètres de fonctionnement; - ladite étape de modification des paramètres de fonctionnement comprend la modification automatique d'au moins une partie desdits paramètres de fonctionnement; et - ladite étape de modification des paramètres de fonctionnement comprend la modification d'au moins un des paramètres choisis par le groupe consis20 tant en: - le débit dudit fluide de refroidissement; - le temps de séjour desdits articles dans ladite enceinte; - le débit de gaz extrait de ladite enceinte; - la mise en vitesse des gaz; - la recirculation des gaz; et - la balance entre les entrées d'air et les sorties de gaz. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés, sur lesquels: - la Fig.1 représente un schéma synoptique illustrant une installation de refroidissement; - la Fig.2 est un organigramme général du procédé de l'invention; - la Fig.3 illustre la modélisation numérique des articles à traiter; - la Fig.4 illustre la modélisation numérique de l'enceinte de refroidissement; et - la Fig.5 représente l'organigramme détaillé du cycle de test du procédé de l'invention. Sur la figure 1, on a représenté une installation classique de traitement d'articles alimentaires, pour laquelle des paramètres de fonctionnement sont déterminés par un procédé selon l'invention. Cette installation comprend une enceinte ou tunnel cryogénique 2, de type classique, permettant la congélation d'articles alimentaires P par leur mise 10 en présence avec un fluide cryogénique 4 véhiculé par une ligne d'alimentation 5, à partir d'une source quelconque. Par exemple, l'enceinte 2 a une forme de parallélépipède rectangle. Comme cela est dit précédemment, le fluide cryogénique 4 utilisé peut être par exemple, de la neige carbonique ou de l'azote liquide et est injecté en un 15 ou plusieurs endroits de l'enceinte 2. Cette enceinte 2 est associée à un convoyeur 6 de type classique, permettant l'introduction des articles P dans l'enceinte 2 et leur extraction et fonctionnant soit de manière séquentielle, soit de manière continue. L'installation présente plusieurs paramètres de fonctionnement, à sa20 voir le profil des températures dans l'enceinte, le temps de séjour des articles P dans l'enceinte 2 ou la vitesse de déroulement du convoyeur 6, et la température d'entrée des articles P. L'installation comporte enfin des moyens de commande 12 de la quantité de fluide cryogénique 4 injectée dans l'enceinte 2. Ces moyens 12 comportent des moyens de pilotage 14 du débit de fluide cryogénique 4. Par exemple, les moyens de pilotage 14 sont constitués par des systèmes d'électrovannes ou des vannes proportionnelles de type classique, disposés sur la ligne d'alimentation 5 en fluide cryogénique 4. Avantageusement, l'installation comporte également un système de 30 ventilation de gaz contrôlant les flux de gaz et la ventilation de l'atmosphère de l'enceinte 2. Par exemple, ce système est composé de ventilateurs spécifiques permettant la mise en vitesse des gaz, de ventilateurs contrôlant la recirculation des gaz et d'une combinaison entre des ventilateurs et des portes mobiles contrôlant la balance entre les entrées d'air et les sorties de gaz. En référence à la figure 2, on va maintenant décrire un organigramme général du procédé de détermination des paramètres de fonctionnement selon l'invention. Ce procédé débute par une étape 16 de saisie d'une consigne portant sur la température de sortie des articles après refroidissement thermique. L'étape 16 est suivie d'une étape 18 de détermination des paramètres de fonctionnement initiaux. Les paramètres déterminés au cours de cette étape 10 18 sont des paramètres connus tels que les caractéristiques mécaniques de l'enceinte 2 ou encore les caractéristiques physiques et thermodynamiques des articles P et des paramètres variables, tels que les paramètres de fonctionnement de l'installation qui sont fixés arbitrairement. Le procédé comporte ensuite une étape 20 de prédiction de la tempé15 rature des articles P en sortie de l'enceinte 2. Cette étape 20 comporte une étape 22 de prédiction du comportement de l'enceinte 2 et une étape 24 de prédiction du comportement des articles P. L'étape 22 de prédiction du comportement de l'enceinte 2, permet de prédire par calcul ainsi que cela est décrit plus loin en référence à la figure 4, le 20 profil théorique des températures du fluide cryogénique à l'intérieur de l'enceinte 2. Les résultats délivrés par l'étape 22 dépendent des caractéristiques thermodynamiques du fluide cryogénique 4, des caractéristiques convectives de l'enceinte 2 ainsi que les caractéristiques des moyens d'injection du fluide cryo25 génique 4, des caractéristiques du système de ventilation et des caractéristiques physiques de l'enceinte 2. L'étape 22 prend également en compte des paramètres de fonctionnement de l'installation comme la vitesse du convoyeur 6. L'étape 24 de prédiction du comportent des articles P permet de dé30 terminer par calcul, ainsi que cela est décrit plus loin en référence à la figure 3, des variations d'enthalpie des articles P en fonction de leur environnement extérieur et de leur température initiale. Ainsi, les résultats délivrés par l'étape 24 de prédiction du comportement des articles P, dépendent de leurs caractéristiques physiques et thermodynamiques. Les étapes 22 de prédiction du comportement de l'enceinte 24 et de 5 prédiction du comportement des articles P, sont couplées l'une à l'autre ainsi que cela est décrit plus en détail en référence à la figure 5, afin de délivrer une température théorique des articles P en sortie de l'enceinte 2. Ainsi, la prédiction réalisée au cours de l'étape 20 de prédiction de la température des articles P en sortie de l'enceinte 2 prend en compte les caracté10 ristiques thermodynamiques et physiques de l'enceinte 2 et des articles P, ainsi que les paramètres de fonctionnement de l'installation. De ce fait, la détermination de la température des articles P en sortie de l'enceinte 2 est dynamique, paramétrable et d'une grande précision. L'étape 20 de prédiction de la température des articles P en sortie de 15 l'enceinte 2 est suivie d'une étape 26 de comparaison de la température prédite lors de l'étape 20, avec la consigne de température déterminée lors de l'étape 16. Par exemple, cette étape 26 de comparaison permet la prise en compte d'un intervalle de tolérance de l'ordre de quelques degrés autour de la consigne de température de sortie déterminée lors de l'étape 16. Quand l'écart entre la température prédite et la consigne est supérieur à un seuil prédéterminé, l'étape 26 de comparaison est suivie d'une étape 28 de modification de paramètres de fonctionnement de l'installation. Par exemple, les paramètres de fonctionnement modifiés au cours de cette étape 28, comprennent la vitesse du convoyeur 6 et les paramètres interve25 nant dans la détermination du profil théorique des températures du fluide 4 dans l'enceinte 2, soit par exemple le débit du fluide cryogénique 4, le contrôle de la ventilation à l'intérieur de l'enceinte 2 et le taux de chargement du convoyeur 6. Les modifications des paramètres de fonctionnement en cours de l'étape 28 peuvent être réalisées directement par un opérateur ou être faites de 30 manière automatique par un calculateur prenant en compte des limites maximale et minimale pour chacun des paramètres. De plus, les modifications peuvent affecter un ou plusieurs paramètre à chaque fois. Il est également possible de définir un ordre de modification des paramètres de fonctionnement pour tenter d'atteindre la température de consigne en ne modifiant qu'un paramètre à la fois. Si la consigne ne peut être atteinte en modifiant un premier paramètre entre des valeurs limites. Ce paramètre est fixé 5 sur une valeur limite ou une valeur moyenne et aux itérations suivantes, le paramètre suivant dans la liste est modifié. A l'issue de l'étape 28 de modification des paramètres de fonctionnement, le procédé reprend à l'étape 20 de prédiction de la température des articles P en sortie de l'enceinte 2, le procédé formant ainsi un cycle de test des paramè10 tres de fonctionnement de l'installation, comprenant l'étape 20 de prédiction de la température des articles en sortie, l'étape 26 de comparaison de la température prédite avec la température de consigne et l'étape 28 de modification des paramètres de fonctionnement. Ce cycle de test, désigné par la référence numérique générale 29, est 15 répété jusqu'à ce que la comparaison effectuée lors de l'étape 26 entre la température prédite et la température de consigne révèle à un écart inférieur à une valeur de seuil prédéterminée. Le cycle 29 est alors interrompu et l'étape 26 est suivie d'une étape 30 d'enregistrement des derniers paramètres de fonctionnement testés, qui forment 20 ainsi une recette. La prédiction réalisée au cours de l'étape 20 étant d'une grande précision et prenant en compte des paramètres de fonctionnement de l'enceinte 2, des caractéristiques thermodynamiques et physiques de l'enceinte 2 et des caractéristiques thermodynamiques et physiques de l'article P, la recette détermi25 née par le procédé de l'invention est précise et proche du fonctionnement réel. De plus, une telle recette est facilement adaptée à des modifications des conditions de fonctionnement. Par exemple, si un paramètre de fonctionnement tel que la température d'entrée des articles dans l'enceinte 2 est modifié, la recette peut être corrig30 ée par l'exécution du procédé de l'invention en utilisant les paramètres actuels de recette comme paramètres de fonctionnement initiaux lors de l'étape 18, l'exécution du procédé de l'invention permettant de déterminer rapidement et simplement les corrections à apporter aux paramètres de fonctionnement pour obtenir la bonne température de sortie des articles P. En référence aux figures 3, 4 et 5, on va maintenant décrire plus en détail l'étape 20 de prédiction de la température des articles P en sortie de l'enceinte 2. Sur la figure 3, on a représenté un exemple de maillage d'un article 5 alimentaire P tel que mis en oeuvre au cours de l'étape 24 de prédiction du comportement des articles P. La prise en compte, dans le procédé de détermination de paramètres de fonctionnement, des caractéristiques thermodynamiques et physiques des articles P lors de l'étape 24 de prédiction de leur comportement, est basée sur 10 une modélisation des articles P à laquelle est appliquée l'équation de conservation de la chaleur discrétisée. En effet, l'équation de la conservation de la chaleur ne peut pas être résolue en tout point de l'espace et en tout instant, par une fonction intégrale simple. La méthode employée consiste à discrétiser cette équation de sorte que l'on ne la résout plus que sur des points spatiaux et temporels appelés noeuds et désignés par la référence numérique générale 32. Après définition d'un maillage de l'article P, on applique l'équation de conservation de la chaleur à chacun des noeuds 32. On obtient ainsi un système d'équation qu'il faut résoudre pour connaître l'état thermique de l'article P dans le temps et dans l'espace. a(ia)+ a )+ a (iaT)=Pa(C.T) aFx ax ay ay ( az âz at X, Y et Z sont des axes définissant un repère spatial orthonormé autour de l'article P. T est la température de l'article P exprimée en kelvin (K), et C 25 sa chaleur spécifique exprimée en watt par kilogramme et par kelvin (W/(Kg*K) ). Les articles alimentaires P que l'on surgèle sont généralement constitués de corps différents. Cela signifie que le changement de phase s'accompagne d'une variation de température et que l'équation de la conservation de la chaleur peut tou30 jours s'appliquer. En revanche, si l'on est amené à traiter un corps pur, l'équation n'est plus continue. Dans ce cas, on simplifie le problème en modifiant la table d'enthalpie du corps pur pour que le changement de phase engendre une faible variation de température. La discrétisation est faite grâce à la méthode mathématique des différences finies en régime variable. De manière connue, celle-ci peut être effectuée de deux manières. La première, la discrétisation implicite, a l'avantage d'être stable quelle que soit la configuration spatiale et temporelle. A un instant donné, elle permet de déterminer la température d'un noeud 32 en fonction de la température des noeuds voisins au même instant. Elle implique cependant des conditions aux 10 frontières constantes et une résolution matricielle du système d'équation formé par chacun des noeuds 32. La seconde, la discrétisation explicite, permet de déterminer directement le température d'un noeud 32 à un instant T+AT d'après les conditions à l'instant T. Le résultat est immédiat, en revanche, il faut choisir un pas de temps 15 adapté de manière à éviter l'instabilité du modèle. La première méthode est préconisée dans le cas o l'on cherche à obtenir principalement la température de surface d'un article, ce qui correspond à l'opération communément appelée opération de " croûtage ". La seconde est préconisée lorsque l'on veut faire de la surgélation et connaître la température à 20 coeur d'un article. Le maillage de l'article P est un problème crucial. De cela dépend directement la simplicité du traitement ultérieur et la précision des résultats. Un nombre de noeuds important amène une grande précision dans le résultat mais impose un temps de calcul élevé. Il s'agit de trouver un compromis 25 entre la précision et le temps de calcul. Par exemple, pour le cas d'un produit alimentaire de dimensions extérieures 100 x 60 x 10 mm avec un maillage régulier tous les millimètres, il faut plus de 17000 noeuds et autant d'équations pour définir le comportement de l'article P. Dans la version décrite de l'invention, l'étape 24 comporte des calculs qui permettent d'optimiser le maillage qui est effectué pour la prédiction du comportement des produits P. Par exemple, dans le cas du croûtage, on surveille plus particulièrement la solidification d'une fine épaisseur de la peau du produit par changement de phase. Il faut donc un maillage dense en périphérie et plus large à coeur. Pour ne pas saisir manuellement les coordonnées de chacun des 5 noeuds et pour conserver des relations simples entre les noeuds et faciliter le traitement, une solution consiste à distribuer des noeuds dans chaque direction de l'espace à l'aide par exemple d'une progression géométrique, comme cela est représenté sur la figure 3. Par exemple, sur l'axe X, on distribue les noeuds de la manière sui10 vante: on considère Ax la valeur du premier terme qui correspond à l'abscisse du premier noeud, et r la raison, différente de 1, de la suite géométrique mis en oeuvre. La valeur du nème terme est: Ax*rnl1, cela correspond à la position sur l'axe X du nème noeud. La somme des n premiers termes est: sa1-rn S = Ax-+1F 1-r La figure 3 représente le positionnement des noeuds selon ce maillage sur un article P parallélépipédique o l'on a imposé une condition de parité sur le nombre de noeuds de manière à simplifier la résolution. On obtient ainsi une valeur correspondant à la dimension du produit P selon l'axe X: LP = Ax+ AX(l + r Ax = X(+ 1 (1- r)) 1+r l(R avec R 1+r Avec I qui correspond à l'abscisse du noeud central sur cette longueur: /L1--i Ln L'imprécision sur l'axe X s'exprime alors de la manière suivante: Lni- L"- fIR1 -/X2 - Dans le cas des opérations de surgélation o l'on cherche à déterminer la température à coeur du produit, un terme correctif peut être inséré dans les formules. Dans le cas de l'axe X, on insère le terme correctif suivant: 1 0 Une autre méthode d'optimisation possible consiste à diminuer le temps de traitement en omettant certains calculs. En effet, sur chaque noeud, on additionne le flux thermique aux six faces de son volume élémentaire. Cependant, il existe des zones o les effets thermiques sont assimilables à des problèmes monodimensionnels. 1 5 Pour exploiter cette particularité, on décompose le traitement en sommant non plus globalement les flux thermiques sur chaque face, mais selon chaque direction. Sur chaque direction, on résout, pour un pas de temps AT, les équations aux noeuds en allant de la frontière vers le coeur, jusqu'à ce que la variation d'enthalpie soit considérée comme étant négligeable car inférieure à un 20 seuil prédéterminé. En effectuant cette opération dans chaque direction, on définit un volume du produit P englobant tous les noeuds pour lesquels les variations d'enthalpie seront négligeables, et donc pour lesquels aucun calcul ne sera fait. On peut ainsi économiser du temps de calcul, surtout dans les pre25 miers instants de l'échange. Dans le cas o l'article P est de forme complexe, on peut le décomposer en un ensemble de formes élémentaires auxquelles on applique le maillage défini ci-dessus ou tout autre maillage adapté à la forme de l'article P. Sur la figure 4, on a représenté schématiquement l'enceinte de traite30 ment 2 des articles alimentaires P telle que modélisée pour la mise en oeuvre de l'étape 22 de prédiction de son comportement. La prise en compte, dans le procédé de détermination des conditions de fonctionnement, des caractéristiques thermodynamiques et physiques de l'enceinte 2 lors de l'étape 22 de prédiction du comportement de l'enceinte 2, est basée sur sa modélisation sous forme de tranches élémentaires. Comme cela a été décrit précédemment en référence à la figure 1, l'enceinte 2 de refroidissement est associée à un convoyeur 6. Elle est alimentée 5 en fluide cryogénique 4 via à une ligne d'alimentation 5. L'enceinte 2 est assimilée à un parallélépipède rectangle. Pour déterminer le profil théorique des températures du fluide 4, la méthode mise en oeuvre pour la prédiction du comportement de l'enceinte 2 consiste à effectuer une succession de bilans thermiques locaux. Dans ce but, on considère une modélisation du système thermodynamique de l'enceinte 2, en régime établi, sous la forme de tranches élémentaires 34, à 34., perpendiculaires à sa longueur. La somme de ces tranches élémentaires 34, à 34., représente le volume interne de l'enceinte 2. Pour chaque tranche élémentaire 34, à 34, on effectue le bilan des 15 transferts thermiques afin de déterminer l'enthalpie du fluide 4 et donc sa température. Ce bilan doit tenir compte: - des déperditions thermiques avec l'extérieur de l'enceinte 2; - du liquide cryogénique 4 injecté dans les zones de pulvérisation; et 20 - des échanges entre les produits P et le fluide 4. Dans le cas de la tranche 34j de l'enceinte 2 de dimensions L*l*h, le bilan thermique se traduit par l'équation suivante: H +2K T(1 +h JT Am, T feO,)hX -A Spray (i,*(H fe(i)-H fLiq)+f tp(IH pe(i)- H ps(() lf m fe(A)+ m fspr (i) Dans cette équation: Hfs(i) correspond à l'enthalpie du fluide cryogénique 4 en sortie de la tranche élémentaire 34j exprimé en joules par kilogramme (J/Kg); Hîfe() correspond à l'enthalpie du fluide cryogénique 4 en entrée de la tranche élémentaire 34j exprimée en joules par kilogramme (J/Kg); H fi correspond à l'enthalpie liquide du fluide cryogénique 4 injecté 30 exprimée en joules par kilogramme (J/Kg); Hpe(i) correspond à l'enthalpie de l'article P en entrée de la tranche 34j exprimée en joules par kilogramme (J/Kg); H ps(0 correspond à l'enthalpiede l'article P en sortie de la tranche 34j exprimée en joules par kilogramme (J/Kg); KT correspond au coefficient d'échange thermique de l'enceinte 2 avec l'extérieur exprimé en watt par mètre carré et par kelvin (W/(m2K)); ff'l"Spray(i) correspond au débit massique de fluide cryogénique 4 vaporisé dans la tranche 34j exprimé en kilogramme par seconde (Kg/s); mfef(i) correspond au débit massique de fluide cryogénique 4 entrant dans la tranche 34j exprimé en kilogramme par seconde (Kg/s); mh correspond au débit massique de produits à traiter exprimé en ki10 logramme par seconde (Kg/s); TAmb correspond à la température ambiante exprimée en kelvin; et Tfe(o correspond à la température du fluide cryogénique 4 en entrée de la tranche 34j exprimée en kelvin. Sur la figure 4 on a également représenté les flux thermiques: 15 1imshpray(i)Hziq est représenté par la lettre A; /nfe(i) Hfe(0) est représenté par la lettre B; fi(i)Hf,(i) est représenté par la lettre C; îhpHpe(,, est représenté par la lettre D; et hpHps(,, est représenté par la lettre E avec mfs(i)=-fe(i) +Afspray(i) Par expérience, on sait que dans certaines conditions de fonctionnement (débit de production trop faible ou température du fluide cryogénique 4 trop basse), le liquide cryogénique 4 injecté n'est que partiellement vaporisé et une fraction du liquide s'écoule vers l'entrée de l'enceinte 2. Si l'on souhaite prendre en compte ce phénomène, il est préférable de résoudre les bilans locaux en commençant par la tranche élémentaire située en sortie de l'enceinte 2. On fait donc les calculs dans le sens inverse du trajet des articles P, soit selon l'axe X tel que représenté sur la figure 4. Dans ce sens en effet, on peut reporter la fraction de liquide non vaporisé dans la suivante et ainsi de suite jusqu'à atteindre les zones de ventilation o les débits injectés sont nuls et o les surplus liquides sont vaporisés. Pour déterminer la fraction de liquide cryogénique 4 non vaporisé dans 5 une tranche élémentaire, on désigne une enthalpie du fluide limite, endessous de laquelle apparaîtra un titre liquide. Cela revient à fixer une température de fluide gazeux minimale dans l'enceinte 2. Le titre liquide non vaporisé en sortie de la tranche élémentaire 34 10 correspond à XL(i) et s'exprime sous la forme suivante: mlfLiq(i) XL(i)= en fs(i) Si on simplifie les calculs en considérant que l'enthalpie de cette fraction liquide est sensiblement égale à l'enthalpie du fluide cryogénique 4 injecté, on obtient l'expression du titre du liquide suivante: (ô fim HAs)) 1 5 XL= W fLim -H JLi Dans cette équation, Hfiim correspond à l'enthalpie limite de formation d'un titre liquide dans une tranche élémentaire de l'enceinte 2. Sur la figure 5, on a représenté le détail du cycle de test 29 et notamment de l'étape 20 de prédiction de la température des articles P en sortie de 20 l'enceinte 2. Afin de pouvoir réaliser une prédiction de la température des articles P en sortie de l'enceinte 2, le procédé de refroidissement fait intervenir l'étape 22 de prédiction du comportement de l'enceinte 2 et l'étape 24 de prédiction du comportement des articles P. On commence par réaliser l'étape 22 de prédiction du comportement de l'enceinte 2, lors d'une étape 40. Cette étape 40 délivre les pertes thermiques 42 par tranche élémentaire, qui sont réintroduites dans l'étape 22. Après avoir réitéré cette opération un certain nombre de fois, on obtient les pertes thermiques totales 44 ainsi que le profil 46 des températures du fluide 4 dans l'enceinte 2, lors de l'étape 40. Pour calculer le bilan thermique de chaque tranche, l'étape 22 requiert 5 les variations enthalpiques des articles P. De fait, lors de la première itération, le profil des températures du fluide 4 dans l'enceinte 2 ne pouvant pas être calculé, il est fixé arbitrairement. On met ensuite en couvre l'étape 24 de prédiction du comportement des articles P, lors d'une étape 50. Cette étape 50 délivre l'enthalpie 52 de 10 l'article P en sortie de l'enceinte 2, soit sa température. Eventuellement, l'étape 24 de prédiction du comportement des articles P délivre également les variations enthalpiques 54 d'un article P pour chaque tranche élémentaire de l'enceinte 2. Dans ce cas, cette information est retournée à l'étape 22 de prédiction de comportement de l'enceinte 2, qui l'insère 15 dans le bilan thermique de chaque tranche élémentaire. L'enthalpie 52 de l'article P en sortie de l'enceinte 2 ainsi que le profil 46 des températures du fluide 4 dans l'enceinte 2 et les pertes thermiques totales 44, sont mis en relation afin de déterminer le débit total du fluide, à l'étape 60. Eventuellement, on obtient également le débit 62 injecté dans chaque 20 tranche élémentaire. Dans ce cas, cette information est retournée à l'étape 22 de prédiction du comportement de l'enceinte 2, qui l'insère dans le bilan thermique de chaque tranche élémentaire. On vérifie ensuite si le profil des températures du fluide 4 dans l'enceinte 2 est stable, à l'étape 70. Par exemple, le profil des températures du fluide est considéré comme stable s'il répond deux fois de suite au critère suivant: Hps(lIk)H ps(lkl) Au premier passage, le profil est considéré comme instable. Tant que le profil est considéré comme instable, on retourne à l'étape et on recommence la succession d'opérations permettant de définir un profil. Une fois qu'un profil stable a été obtenu, on vérifie lors de l'étape 26, si la consigne déterminée lors de l'étape 26 et portant sur la température des articles P en sortie de l'enceinte 2 a été atteinte. Si la consigne a été atteinte, les paramètres de fonctionnement abou5 tissant au dernier profil des températures du fluide 4 à l'intérieur de l'enceinte 2 sont enregistrés lors de l'étape 30 et forment une recette. Si la consigne n'a pas été atteinte, à l'étape 28 les paramètres de fonctionnement de l'installation sont modifiés. Par exemple, cette modification comporte une correction 102 sur le débit du fluide 4 avant de réitérer l'algorithme. 10 Eventuellement, la modification comporte une correction 104 directement sur les paramètres de fonctionnement conditionnant le profil des températures du fluide 4, qui est utilisée à l'étape 22 de prédiction du comportement de l'enceinte 2. Dans cet exemple, la prédiction de la température des articles en sortie de l'enceinte est utilisée pour déterminer le débit du fluide 4 injecté dans une 15 enceinte cryogénique 2. On peut de la même manière influer sur le temps de séjour des articles P dans l'enceinte 2 en modifiant la vitesse du convoyeur 6 ou les temps d'arrêt dans le cas d'un convoyeur séquentiel. On peut également agir sur la vitesse d'extraction des gaz ou sur le taux de chargement, ou sur une combinaison de 20 ces paramètres. De même, il est possible d'influer sur la balance entre les entrées d'air et les sorties de gaz, la vitesse d'extraction des gaz, la mise en vitesse des gaz, ou encore la recirculation des gaz en influant sur les éléments de contrôle de ces paramètres. Par ailleurs, si on met en oeuvre le procédé de l'invention pour une installation disposant de capteurs sans contact de la température des articles en sortie, par exemple des capteurs basés sur le rayonnement thermique ou l'image infrarouge, ou encore sur une mesure par thermométrie micro-onde (TMO), tel que le capteur décrit dans le brevet FR-A-2 771 552, les résultats délivrés lors de 30 l'étape de prédiction de la température des articles en sortie de l'enceinte peuvent alors être recoupés avec les mesures délivrées par ces capteurs. Dans ce cas, on utilise l'une ou l'autre des informations pour vérifier l'autre. Dans un autre mode de réalisation, l'information délivrée par le capteur est utilisée pour corriger la prédiction. Bien qu'un mode particulier de réalisation ait été décrit, il n'est pas considéré comme limitatif de la portée de la présente invention. En variante, le procédé de refroidissement de l'invention peut aussi être appliqué dans une installation de froid mécanique présentant un dispositif d'échange de chaleur indirect. L'invention a été décrite dans le cas du refroidissement d'articles alimentaires, cependant on peut l'appliquer également à d'autres types d'articles 10 notamment métalliques. De plus, le terme refroidissement couvre également les systèmes visant au maintien et au contrôle d'une température inférieure à la température initiale d'un article. Le procédé de l'invention peut être mis en oeuvre à l'aide par exemple 15 d'un programme exécuté sur un ordinateur ou encore toute autre solution logicielle et/ou matérielle adaptée. REVENDICATIONS 1. Procédé de détermination des paramètres de fonctionnement d'une installation de refroidissement thermique d'articles (P), comprenant une enceinte (2) traversée par lesdits articles (P) d'une entrée à une sortie et utilisant un fluide de refroidissement (4), le procédé comprenant: une étape (16) de détermination d'une consigne de température d'articles (P) en sortie de ladite enceinte (2); - une étape (18) de détermination de paramètres de fonctionnement initiaux pour ladite installation; et - un cycle (29) de test des paramètres de fonctionnement comprenant: - une étape (20) de prédiction de la température d'articles (P) en sortie de ladite enceinte (2); - une étape (26) de comparaison de la température de consigne 15 avec la température prédite; et - si ladite étape (26) de comparaison révèle un écart supérieur à un seuil prédéterminé, une étape (28) de modifications de paramètres de fonctionnement de ladite installation (2) et une réitération du cycle de test, caractérisé en ce que ladite étape (20) de prédiction est réalisée à par20 tir de paramètres de fonctionnement de ladite enceinte (2), de caractéristiques thermodynamiques et physiques de ladite enceinte (2) et de caractéristiques thermodynamiques et physiques desdits articles (P). 2. Procédé selon la revendication 1, caractérisé en ce que ladite étape (20) de prédiction comporte une étape (22) de prédiction du comportement de 25 ladite enceinte (2) fondée sur la résolution de bilans thermiques sur des tranches élémentaires du volume de ladite enceinte (2), réalisée au moins à partir de caractéristiques thermodynamiques dudit fluide de refroidissement (4) et de caractéristiques thermodynamiques et physiques de ladite enceinte (2). 3. Procédé selon la revendication 2, caractérisé en ce que ladite étape 30 (22) de prédiction du comportement de ladite enceinte (2) est réalisée en outre, à partir de paramètres de fonctionnement de ladite installation. 4. Procédé selon la revendication 3, caractérisé en ce que lesdits paramètres de fonctionnement de ladite installation représentent au moins l'un des éléments choisis dans le groupe consistant en: - la vitesse d'un convoyeur (6) de transport desdits articles (P) au travers de ladite enceinte (2); - le taux de chargement; et - la ventilation de l'atmosphère de ladite enceinte (2). 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que ladite étape (20) de prédiction comporte une étape (24) de prédiction du comportement desdits articles (P) fondée sur la résolution de l'équation de conservation de la chaleur discrétisée et appliquée à un réseau de points spatiaux et temporels constituant un maillage desdits articles (P), réalisée au moins à 10 partir de caractéristiques thermodynamiques et physiques desdits articles (P). 6. Procédé selon la revendication 5, caractérisé en ce que ladite étape (24) de prédiction du comportement desdits articles (P) est réalisée en outre, à partir de paramètres de fonctionnement de ladite installation. 7. Procédé selon la revendication 6, caractérisé en ce que lesdits pa15 ramètres de fonctionnement de ladite installation comprennent la température desdits articles (P) en entrée de ladite enceinte (2). 8. Procédé selon l'une quelconque des revendications 5 à 7, caractérisé en ce que ladite étape (24) de prédiction du comportement desdits articles (P) est optimisée par des calculs de modification dudit maillage desdits articles (P) 20 selon des suites mathématiques. 9. Procédé selon l'une quelconque des revendications 5 à 8, caractérisé en ce que ladite étape (24) de prédiction du comportement desdits articles (P) est optimisée par suppression des calculs de prédiction pour des points spatiaux et temporels dudit maillage desdits articles (P) pour lesquels les variations 25 d'enthalpie sont inférieures à un seuil prédéterminé. 10. Procédé selon les revendications 2 et 5 prises ensembles, caractérisée en ce que ladite étape (20) de prédiction de la température desdits articles (P) en sortie de ladite enceinte (2) est fondée sur ladite étape (22) de prédiction du comportement de ladite enceinte (2) ainsi que sur ladite étape (24) de prédic30 tion du comportement desdits articles (P). 11. Procédé selon l'une quelconque des revendications 1 à 10, caractérisé en ce que ladite étape (28) de modification des paramètres de fonctionnement comprend une étape de modification manuelle d'au moins une partie des paramètres de fonctionnement. 12. Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce que ladite étape (28) de modification des paramètres de fonctionnement comprend la modification automatique d'au moins une partie desdits paramètres de fonctionnement. 13. Procédé selon l'une quelconque des revendications 1 à 12, caractérisé en ce que ladite étape (28) de modification des paramètres de fonctionnement comprend la modification d'au moins un des paramètres choisis par le groupe consistant en: - le débit dudit fluide de refroidissement (4); - le temps de séjour desdits articles (P) dans ladite enceinte (2); - le débit de gaz extrait de ladite enceinte (2); - la mise en vitesse des gaz; - la recirculation des gaz; et - la balance entre les entrées d'air et les sorties de gaz.