Brevet FR2965922A1 — publié le 20120413

L'invention concerne l'identification de la conductivité k et/ou de la capacité thermique C d'une paroi plane, comme par exemple mais non exclusivement, une paroi de bâtiment. Dans la construction, la connaissance de l'un et/ou l'autre de ces deux paramètres est essentielle dans la mesure où elle permet de calculer de manière précise les déperditions thermiques susceptibles de se produire à travers la paroi étudiée. Une telle étude peut alors être exploitée pour définir une stratégie de rénovation du bâtiment visant à réduire la consommation énergétique nécessaire pour le chauffer. Rappelons l'importance des enjeux énergétiques liés au chauffage des bâtiments : En France, le parc de bâtiments est responsable de 40 % des consommations d'énergie. Dans ce document, la capacité thermique C (en J.K-» désigne la quantité de chaleur que peut emmagasiner la paroi considérée ; elle est également égale au produit entre la masse volumique et la chaleur massique. Une paroi dont la capacité thermique est élevée est souvent qualifiée de paroi d'inertie importante ou « paroi lourde ». Par ailleurs, la conductivité k de la paroi (en W.m-'.K» est égale à la quantité de chaleur transférée par unité de surface et par unité de temps sous un gradient de température de 1°C par mètre. En outre, le flux thermique surfacique sur chacune des faces de la paroi est contrôlé par la conductivité de la paroi, le flux thermique surfacique étant la quantité de chaleur par unité de temps et par unité de surface (en W.m-Z) qu'échange une face de la paroi avec l'atmosphère ou les objets environnants. Il est égal à la variation du champ de température selon l'axe perpendiculaire à la paroi multipliée par la conductivité k. METHODES ACTUELLES Les méthodes connues pour l'identification de la conductivité k et/ou de la capacité thermique C d'une paroi plane reposent sur l'observation de la réponse en température de la paroi suite à une sollicitation thermique. Cette sollicitation peut être imposée ou non. Cela permet donc de distinguer les méthodes dites actives, dans lesquelles une quantité de chaleur connue est communiquée à la paroi, des méthodes passives dans lesquelles la sollicitation provient de l'environnement extérieur et n'est pas contrôlée. Les méthodes actives nécessitent une parfaite maîtrise des conditions environnementales et sont de ce fait difficilement applicables in situ. Par ailleurs, une partie des méthodes d'identification des propriétés thermiques d'une paroi sont des méthodes dites 'en régime permanent'. Ces méthodes se caractérisent par la mesure simultanée d'un flux de chaleur traversant l'échantillon et d'une différence de température ; elles permettent de déterminer la conductivité k de la paroi mais sont insuffisantes pour déterminer sa capacité thermique C. Parmi ces méthodes figurent la méthode de la « plaque chaude gardée » ou la méthode du « fil chaud », couramment employées dans l'industrie. Y figure également l'essai normalisé « Hot Box », qui utilise une méthode active en régime permanent pour la détermination de la résistance d'une paroi, et est largement utilisé aujourd'hui comme méthode de certification des matériaux de construction. Ce dernier essai n'est cependant envisageable que dans des conditions expérimentales de laboratoire. A contrario, la méthode Flash et ses variantes sont des méthodes de mesure en régime transitoire. Ces méthodes consistent à imposer un échange de chaleur variable (échelon, Dirac, sinusoïde, ...) sur l'une des faces et à observer la réponse sur l'autre face. Ces méthodes permettent d'obtenir les caractéristiques dynamiques de la paroi et notamment les coefficients k et C. Cependant, ces méthodes sont également difficilement utilisables pour des mesures sur site. En outre, certaines méthodes permettant la détermination des coefficients C et k nécessitent de mesurer la température en différents points à l'intérieur même de la paroi. De telles méthodes, étant destructives, sont en général inenvisageables pour des mesures sur site. Aussi, il n'existe actuellement aucun dispositif ni aucune méthode utilisable in situ, de mise en oeuvre simple qui permette de déterminer la conductivité k et/ou la capacité thermique C d'une paroi. Un premier objectif de l'invention est donc de proposer un dispositif d'identification de la conductivité k et/ou de la capacité thermique C d'une paroi plane s'étendant entre une première et une deuxième faces, le dispositif comportant : un premier capteur de température apte à être placé sur la première face pour acquérir des mesures de la température de celle-ci ; un deuxième capteur de température apte à être placé sur la deuxième face pour acquérir des mesures de la température de celle-ci ; un système d'acquisition de données, apte à acquérir périodiquement des mesures fournies par lesdits capteurs ; et des moyens de traitement des données aptes à déterminer le coefficient de conductivité k et/ou la capacité thermique C de la paroi ; dispositif qui soit utilisable in situ, simple à mettre en oeuvre, et dont l'utilisation ne nécessite pas de porter atteinte à l'intégrité de la paroi. Cet objectif est atteint grâce au fait que le dispositif comporte en outre un troisième capteur, qui est un capteur de température apte à être disposé au voisinage de la première face pour fournir une mesure de la température, ou un capteur de flux thermique surfacique, apte à être placé sur la première face de la paroi pour mesurer le flux thermique surfacique de la première face ; et que les moyens de traitement déterminent le coefficient de conductivité k et/ou la capacité thermique C, en calculant l'évolution des températures dans la paroi pendant un intervalle de temps, à partir notamment des mesures fournies par ledit troisième capteur pendant ledit intervalle, et en retenant pour le coefficient de conductivité k et/ou la capacité thermique C, la ou les valeurs pour lesquelles, durant ledit intervalle, les différences entre les valeurs de température calculées pour les première et deuxième faces et les températures mesurées respectivement par les premier et deuxième capteurs sont minimales. L'avantage du dispositif selon l'invention est qu'il n'est nécessaire de réaliser que trois mesures différentes : les températures sur les deux faces de la paroi, ainsi que la température ambiante du côté de la première face (température au voisinage de la première face) ou le flux thermique surfacique échangé par la première face. Dans le cas de parois de bâtiments, en considérant que la première face correspond à la face intérieure, il suffit donc par exemple que le troisième capteur mesure la température intérieure, ou le flux thermique surfacique sur la face intérieure de la paroi, qui est en général d'accès relativement facile. La mesure du flux thermique surfacique sur la face extérieure de la paroi n'est en revanche pas nécessaire, parce que celui-ci est estimé grâce aux moyens de traitement des données. Comme on peut le constater, les mesures prévues selon l'invention restent simples à réaliser. En outre, ces mesures ne sont pas destructives et ne nécessitent pas la mesure de température à l'intérieur même de la paroi. En ce qui concerne le matériel nécessaire pour la mise en oeuvre de l'invention, avantageusement, celui-ci est relativement simple et courant : les premiers et deuxièmes capteurs peuvent être des capteurs de température usuels. Le premier et/ou le deuxième capteur peuvent ainsi fonctionner par mesure du rayonnement infrarouge, ou encore être des thermocouples, des thermistances, etc. On comprend par ailleurs que, pour que le dispositif produise un résultat pleinement satisfaisant, il est préférable que l'intervalle de temps soit suffisamment long pour que pendant celui-ci, la température varie sur l'une et/ou l'autre des faces de la paroi. Les moyens de traitement sont en général un ordinateur qui exécute un programme qui exploite les mesures fournies par les capteurs grâce au système d'acquisition de données. Le coefficient de conductivité k et/ou la capacité thermique C sont déterminés de manière à faire coïncider les températures calculées et mesurées sur les deux faces de la paroi. Dans un mode de réalisation préféré de l'invention, pour déterminer les valeurs minimales des différences entre les températures calculées et mesurées sur les deux faces de la paroi, les moyens de traitement exécutent un programme informatique qui minimise un critère d'erreur. Dans ce programme, en général, la détermination du coefficient de conductivité k et/ou de la capacité thermique C est faite via un algorithme itératif qui permet de minimiser le critère d'erreur. Dans un mode de réalisation, les moyens de traitement sont aptes à exécuter : - une phase d'initialisation durant laquelle des valeurs approchées des températures dans la paroi pendant ledit intervalle de temps sont calculées ; et - une phase de calcul durant laquelle le coefficient de conductivité k et/ou la capacité thermique C sont calculés en minimisant lesdites différences, seules étant utilisées dans cette phase de calcul les mesures acquises par lesdits premier, deuxième et troisième capteurs. La phase d'initialisation est celle au cours de laquelle, pour la première fois, de premières valeurs approchées des températures dans la paroi sont calculées, les paramètres k et/ou C n'étant pas modifiés pendant cette phase. Ce calcul prend en compte les différentes mesures fournies par les capteurs, et en particulier la mesure du flux thermique surfacique de la première face de la paroi, ou la mesure de température au voisinage de la première face de la paroi. Ensuite, au cours d'une deuxième phase dite de calcul, les moyens de traitement exécutent une partie du programme qui vise à identifier les valeurs minimales des différences entre les températures mesurées et calculées, sur les deux faces. C'est pendant cette phase que les valeurs précises des grandeurs k et C sont calculées. Avantageusement, pour cette deuxième phase, seules sont utilisées les mesures provenant des trois capteurs, à l'exclusion d'autres mesures fournies par des capteurs. Grâce à cela, la détermination des coefficients k et/ou C peut être réalisée en utilisant seulement des grandeurs relativement simples à mesurer ou obtenir (températures sur les deux faces, et mesure du flux thermique surfacique ou mesure de température ambiante sur la première face). Dans un mode de réalisation, l'architecture du programme exécuté pour déterminer le coefficient de conductivité k et/ou la capacité thermique C, est la suivante. Les moyens de traitement comportent : - un module de simulation, apte à calculer les valeurs de température de la paroi pendant ledit intervalle, en fonction notamment de la conductivité k et/ou de la capacité thermique C de la paroi, - un module d'adaptation, apte à déterminer une ou des valeurs modifiée(s) d'au moins une variable d'entrée du module de simulation, dont notamment le coefficient de conductivité k et/ou la capacité thermique C, de manière à réduire ou minimiser les différences précitées. Le module de simulation peut également calculer les valeurs de température 25 de la paroi en prenant en compte comme paramètre, entre autres, le flux thermique surfacique (noté OL) sur la deuxième face de celle-ci. Le module d'adaptation est une partie de programme qui vise à modifier, en général en plusieurs étapes de manière itérative, les valeurs des paramètres utilisés par le module de simulation pour calculer les températures dans la paroi. Ces 30 modifications sont réalisées de manière à réduire ou à minimiser la valeur du critère d'erreur. Les variables que le module d'adaptation est apte à modifier comprennent la conductivité k et/ou la capacité thermique C de la paroi, et éventuellement un ou plusieurs des paramètres suivants : flux thermique surfacique de la deuxième face, champ initial de température dans la paroi, et un ou des coefficient(s) d'échange thermique liant le flux thermique surfacique de la première et la deuxième face respectivement à la température en surface de la paroi. De nombreux algorithmes peuvent être utilisés pour la minimisation du critère d'erreur. Il est possible par exemple d'utiliser l'algorithme de Levenberg-Marquardt. Dans un mode de réalisation, à chaque itération, le module de simulation détermine les valeurs de la température dans la paroi ; puis le module d'adaptation détermine la ou les valeur(s) modifiée(s) de ladite ou desdites variable(s) modifiée(s), qui minimisent ou du moins réduisent la valeur du critère d'erreur. Différents modes de réalisation de l'invention envisageables se distinguent par la manière dont le flux thermique surfacique, sur la première et/ou la deuxième face de la paroi, est mesuré et/ou exploité. Dans un mode de réalisation, le flux thermique surfacique sur la deuxième face de la paroi est mesuré. Dans ce cas, en plus des trois premiers capteurs, le dispositif comprend un quatrième capteur de mesure apte à être disposé au voisinage de ladite deuxième face pour fournir des valeurs approchées du flux surfacique échangé par ladite deuxième face ; et durant la phase d'initialisation, les valeurs approchées des températures dans la paroi sont calculées à partir notamment de valeurs approchées du flux thermique surfacique mesurées par ledit quatrième capteur pendant ledit intervalle de temps. Grâce à cette estimation du flux thermique surfacique échangé par la deuxième face, le calcul des valeurs de température dans la paroi (notamment celui effectué dans une phase d'initialisation) est plus rapide et plus précis. Cela étant, que sa valeur soit ou non mesurée (ou estimée à l'aide de mesures provenant d'un capteur), la valeur du flux thermique surfacique de la deuxième face peut être déterminée par les moyens de traitement dans certains modes de réalisation. Pour cela, les moyens de traitement peuvent avoir recours à l'équation suivante : O(L, t) (S3) + a F~( L ( 0, T] dans laquelle aL est un coefficient d'échange thermique de la deuxième face et OL est le flux thermique surfacique de la deuxième face. Dans le programme informatique exécuté par les moyens de traitement, les mesures fournies par le troisième capteur servent à prendre en compte les échanges thermiques réalisés par la première face de la paroi. Dans un mode de réalisation, les moyens de traitement calculent les valeurs de température dans la paroi en prenant en compte l'équation suivante relative aux échanges thermiques sur la première face : v~0(0, t) (S2) -éz, I- cxoH(0, t,) = cxoT(t,) [0, T] dans laquelle : a0 est un coefficient d'échange thermique ; 00, qui est égal au produit a0 x T(t), est le flux thermique surfacique, de la première 15 face ; et T(t) est la température mesurée par le troisième capteur, exprimée en fonction du temps. Une contrainte à prendre en compte lors de la réalisation des mesures tient au fait que celles-ci doivent être faites pendant un intervalle de temps durant lequel 20 se produit une variation de température. Bien souvent, il n'est pas nécessaire d'utiliser une source de chaleur contrôlée ; le dispositif ou la méthode selon l'invention peut être mis(e) en oeuvre en exploitant les variations de température affectant naturellement la paroi au cours du temps. Cependant, pour gagner du temps il peut être préférable d'utiliser une source de chaleur. Dans ce cas, le 25 dispositif comporte en outre au moins une source de chaleur, apte à être disposée en regard d'une face de ladite paroi et à solliciter celle-ci par convection ou par rayonnement, modifiant ainsi le flux thermique surfacique de ladite face. Il s'agit ici d'une source de chaleur dont on peut régler la position, et dans une certaine mesure la puissance rayonnée ; par exemple, un convecteur électrique. La source de chaleur 30 est utilisée en général pour provoquer une soudaine élévation de température sur un côté de la paroi. La réaction de la paroi suite à ce changement permet alors de déterminer relativement rapidement les valeurs du coefficient de conductivité k et/ou de la capacité thermique C.10 Un deuxième objectif de l'invention est de proposer une méthode d'identification de la conductivité k et/ou de la capacité thermique C d'une paroi plane s'étendant entre une première et une deuxième face, méthode qui soit utilisable in situ, qui soit simple à mettre en oeuvre, et ne nécessite pas de porter atteinte à l'intégrité de la paroi. Cet objectif est atteint grâce au fait que la méthode comporte les étapes suivantes : a) on dispose : sur une première face de la paroi, un premier capteur de température, afin d'acquérir des mesures de la température de la première face ; - sur une deuxième face de la paroi, un deuxième capteur (A2) de température, afin d'acquérir des mesures de la température de la deuxième face ; - un troisième capteur, qui est soit un capteur de température disposé au voisinage de la première face pour fournir une mesure de la température, soit un capteur de flux thermique surfacique placé sur la première face de la paroi pour mesurer le flux thermique surfacique (1)0 de cette première face ; b) pendant un intervalle de temps, on acquiert périodiquement les températures des première et deuxième faces de la paroi, ainsi que les mesures fournies par le troisième capteur; puis, c) on détermine le coefficient de conductivité k et/ou la capacité thermique C, en calculant les températures dans la paroi pendant ledit intervalle de temps, à partir notamment des mesures fournies par ledit troisième capteur pendant ledit intervalle ; et en retenant pour le coefficient de conductivité k et/ou la capacité thermique C, la valeur pour laquelle ou les valeurs pour lesquelles, durant ledit intervalle, les différences entre les valeurs de température calculées pour les première et deuxième faces et les températures mesurées respectivement par les premier et deuxième capteurs sont minimales. Dans un mode de réalisation de cette méthode, l'étape c) comporte : c1) une phase d'initialisation pendant laquelle on détermine des valeurs approchées des températures dans la paroi pendant ledit intervalle ; puis c2) une phase de calcul durant laquelle on détermine le coefficient de conductivité k et/ou la capacité thermique C en minimisant lesdites différences, seules étant utilisées dans ce calcul les mesures acquises par lesdits premier, deuxième et troisième capteurs. Dans un mode de réalisation, pendant l'intervalle de temps des mesures, la première face de la paroi est protégée des rayonnements thermiques par un écran opaque aux radiations. Par ailleurs, lorsque les variations de température naturellement causées par l'environnement de la paroi (soleil, etc...) sont suffisamment rapides, dans un mode de réalisation, pendant l'intervalle de temps des mesures, la paroi n'est pas soumise à l'action d'une source de chaleur contrôlée. La mise en oeuvre de la méthode est dans ce cas particulièrement simple. Dans un autre mode de réalisation au contraire, durant au moins une partie de l'intervalle de temps des mesures et/ou pendant une période précédant cet intervalle, on dispose une source de chaleur contrôlée, en regard d'une face de ladite paroi pour modifier le flux thermique surfacique de cette face. Dans un mode de réalisation, pendant au moins une partie de l'intervalle de temps, le flux thermique surfacique sur la deuxième face de la paroi est mesuré au moyen d'un capteur de flux thermique. L'invention sera bien comprise et ses avantages apparaîtront mieux à la lecture de la description détaillée qui suit, d'un mode de réalisation représenté à titre d'exemple non limitatif. La description se réfère aux dessins annexés, sur lesquels : - la figure 1A est une coupe d'une paroi faisant apparaître la position des capteurs d'un premier dispositif selon l'invention ; - les figures 1B et 1C sont des coupes d'une paroi, représentées respectivement perpendiculairement à l'axe de la paroi (fig.1B) et en perspective (fig.1C), et faisant apparaître la position des capteurs d'un deuxième dispositif selon l'invention ; - les figures 2 et 3 sont des représentations schématiques montrant deux dispositifs selon l'invention, dans deux modes de réalisation différents ; - la figure 4 est une coupe respectivement d'une paroi homogène et d'une paroi multicouche ; - la figure 5 est un diagramme montrant une comparaison entre deux courbes d'évolution du flux de température en fonction du temps, sur la deuxième face de la paroi : la première représentant l'évolution réelle du flux, et la deuxième représentant l'évolution du flux telle qu'estimée grâce à l'invention ; - les figures 6 et 7 montrent la convergence des paramètres de conductivité thermique (k) et de capacité thermique (C), en fonction du nombre d'itérations de l'algorithme de minimisation ; et - la figure 8 est un diagramme représentant de manière schématique le fonctionnement de l'algorithme mis en oeuvre dans les dispositifs d'identification présentés par les figures 1 à 3. De manière générale, l'identification des paramètres k et C selon l'invention comporte trois phases : a) La mise en place du dispositif ; une phase b) d'acquisition des mesures, et une phase c) de traitement des mesures permettant l'identification du ou des coefficients cherchés. Ces différentes phases vont successivement être présentées. Mise en place du dispositif de mesure L'acquisition des mesures est réalisée au moyen de dispositifs d'identification de la conductivité k et/ou de la capacité thermique C d'une paroi. Ces dispositifs (référence 10) sont présentés dans différents modes de réalisation par les figures 1 (1A à 1C) à 4. La figure 1 présente un dispositif 10 disposé sur une paroi 12, dans un premier mode de réalisation. La paroi 12 est une paroi plane comme par exemple une paroi d'un bâtiment. Cette paroi s'étend entre une première face FO et une deuxième face FL, qui s'étendent respectivement dans les plans (x=0) et (x=L). Le dispositif 10 comporte trois capteurs : - un premier capteur AO de température placé sur la première face FO de la paroi pour acquérir des mesures de la température (MM» de celle-ci ; - un deuxième capteur AL de température placé sur la deuxième face FL de la paroi pour acquérir des mesures de la température (8d2(t)) de celle-ci ; - un troisième capteur Aamb, disposé au voisinage de la première face de la paroi. Dans l'exemple présenté, le capteur Aamb est un capteur de température ; il aurait été possible d'utiliser par exemple un capteur de flux thermique. Chacun des trois capteurs indiqués précédemment peut être un capteur simple, fournissant directement sa propre mesure, ou peut être un ensemble de capteurs individuels 30, qui sont exploités ensemble de manière à fournir une mesure agrégée unique, obtenue à partir des mesures fournies par les capteurs individuels. En l'occurrence, le capteur AO comme le capteur AL comporte trois capteurs individuels 30, comme cela est représenté sur les figures 2 et 3. Le dispositif 10 comporte en outre un système 16 d'acquisition de données (fig.2,3), apte à acquérir périodiquement des mesures fournies par lesdits capteurs, et des moyens 14 de traitement des données ainsi acquises. Ces moyens de traitement comportent en général un ordinateur, par exemple de type PC et de préférence un PC portable, qui exécute un programme de traitement des données de mesure. Le système d'acquisition des données peut être une carte d'acquisition de données ou un système externe 16 comme représenté en figures 3 et 4. Le PC 14 (en tant que moyen de traitement) permet d'exécuter le programme de traitement des données de mesure. Grâce à ce logiciel, le PC 14 permet de calculer l'évolution des températures (6(x,t)) dans la paroi, en fonction de la position (x) suivant l'axe perpendiculaire à la paroi et du temps (t), comme cela sera détaillé par la suite. Dans un deuxième mode de réalisation (figures 1B et 1C), un dispositif 10 comporte en outre un capteur de flux thermique Ae, disposé au voisinage de la deuxième face de la paroi. Celui-ci est placé sur la paroi extérieure FL au voisinage du capteur de température AL, de manière à fournir une évaluation du flux thermique surfacique de la face externe FL. Dans ce deuxième mode de réalisation, le dispositif 10 comporte de plus un écran 20 de protection contre les rayonnements incidents parasites. Celui-ci a la forme d'une plaque plane carrée, de dimensions environ d x d m2, d étant la longueur du côté de l'écran 20. Pour que l'écran 20 apporte une protection satisfaisante, il convient qu'il soit opaque aux rayonnements, que sa face dirigée du côté de la paroi 12 soit non-émissive, qu'il soit positionné près de la paroi, et que sa dimension d soit nettement supérieure à l'épaisseur de la paroi. On peut avoir par exemple, d > 5L, et e Grâce à l'écran 20, opaque aux radiations, pendant la durée des mesures la première face de la paroi est protégée des rayonnements thermiques. Cela évite des apports de chaleur sous forme de rayonnement reçu par la première face, et permet de s'affranchir de l'utilisation d'un capteur de flux sur cette première face. Le flux thermique surfacique sur la première face peut ainsi être estimé simplement à partir de la température de surface et de la température ambiante, sur le premier côté de la paroi. Les figures 2 et 3 présentent deux modes de réalisation du dispositif 10, présentant quelques différences qui vont être détaillées plus loin. Malgré ces différences, les mêmes références numériques sont attribuées aux éléments identiques ou similaires dans les deux modes de réalisation. Une première différence réside dans le mode de transmission des mesures entre les capteurs A0, AL et le système d'acquisition de données 16. Dans le mode de réalisation de la figure 2, chacun des capteurs individuels 30 et 40 composant le capteur mis en oeuvre (A0, AL, Aamb) communique directement par voie hertzienne avec le système d'acquisition des données 16, relié pour sa part par câble au PC 16. En revanche, dans le mode de réalisation de la figure 3, chacun des capteurs individuels composant les capteurs mis en oeuvre (A0, AL, Aamb) est relié par fil avec un boîtier de communication 18. Il y a donc deux boîtiers 18, chacun étant placé sur une face de la paroi et étant relié aux capteurs individuels disposés sur ou au voisinage de cette face de la paroi. Les boîtiers 18 communiquent par voie hertzienne avec le système d'acquisition des données 16, relié par câble au PC 16. Une deuxième différence entre les modes de mise en oeuvre de l'invention correspondant aux figures 2 et 3 réside dans le fait que la méthode présentée par la figure 2 est une méthode passive, sans mise en place par l'opérateur d'une source de chaleur, alors que la méthode présentée par la figure 3 est active et suppose la mise en place par l'opérateur d'une source de chaleur 32. Ainsi, dans le mode de mise en oeuvre de la figure 2, pendant la durée des mesures, la paroi n'est pas soumise à l'action d'une source de chaleur contrôlée. La chaleur reçue par la paroi est le résultat des fluctuations de température de l'air environnant les deux faces de la paroi et du rayonnement infrarouge reçu par la deuxième face de celle-ci. Inversement, dans le mode de mise en oeuvre de la figure 3, durant les mesures (ou juste avant, ou au début, des mesures), la source de chaleur 32 est placée en regard de la deuxième face de la paroi de manière à modifier le flux thermique surfacique reçu par cette face. Les mesures peuvent ainsi être effectuées plus rapidement que dans le cas précédent. Un délai d'une heure peut être suffisant. La source de chaleur peut être placée en regard d'une face ou l'autre de la paroi. Acquisition des mesures Une fois le dispositif en place, les mesures sont effectuées. Pendant un intervalle de temps on acquiert une série de mesures de températures à l'aide des capteurs A0, AL et Aamb, ainsi que des mesures du flux thermique sur la deuxième face, à l'aide du capteur Ae (ces mesures étant optionnelles). L'intervalle de temps choisi est suffisamment long pour que des variations significatives de sollicitations thermiques de la paroi aient lieu pendant celui-ci. Typiquement, un intervalle de 24 heures peut être choisi. Une fois les mesures effectuées, l'étape c) de traitement des données est réalisée : Traitement informatique des données La figure 8 représente de manière simplifiée les traitements informatiques réalisés pour déterminer les coefficients k et/ou C. Ces traitements comportent deux modules principaux : un module de simulation 60 et un module d'adaptation 70. Ces modules reposent sur une modélisation des phénomènes physiques mis en oeuvre, qui va d'abord être présentée : Le module 60 de simulation numérique est fondé sur un modèle numérique de l'évolution des températures dans la paroi. Il permet de déterminer le champ de températures 0(x,t) dans la paroi, à partir de la mesure du troisième capteur Aamb et d'un vecteur U ((DL(t), 0i(x), k, C, a0, aL) de paramètres, fourni en entrée au module 60. Dans le mode de réalisation présenté, ce vecteur U comprend les paramètres suivants :,échange thermique de la deuxième face OL(t), champ initial de température Oi(x) dans la paroi, conductivité k et capacité thermique C de la paroi, et les coefficients d'échange thermique a0, aL mis en ceuvre dans les équations d'échange thermique sur les deux faces, citées précédemment. (On notera que dans un autre mode de réalisation, si au lieu de mesurer la température ambiante sur la première face, on mesurait le flux thermique surfacique reçu par la première face, alors, le vecteur U pourrait inclure en outre le flux thermique sur la première face (Mn. .2965922 14 Le modèle numérique mis en oeuvre repose sur une représentation simplifiée de la paroi : On considère que la paroi est plane et présente des surfaces parallèles FO et FL disposées dans les plans (x=0) et (x=L). On suppose que son épaisseur est faible 5 par rapport aux deux autres dimensions, de telle sorte que les échanges thermiques à l'intérieur de la paroi ont lieu dans la direction X de l'épaisseur, perpendiculaire aux faces FO et FL. La première surface est en contact avec l'air à la température T(t) et échange de la chaleur par convection avec celui-ci. Selon le modèle numérique (et les 10 capteurs) utilisé(s), on peut supposer que la première surface échange ou non de la chaleur sous forme de rayonnement. La chaleur échangée se traduit par un flux thermique surfacique (D0. La deuxième surface est en contact avec l'air et échange de la chaleur par convection avec celui-ci. Elle échange également de la chaleur sous forme de 15 rayonnement avec l'environnement (soleil, bâtiments environnants, ciel, sol,...) ce qui se traduit par un flux thermique surfacique c13L. Les échanges de chaleur ont lieu de manière dynamique entre chacune des deux surfaces et les environnements extérieur et intérieur. Dans le mode de réalisation présenté, ces échanges sont pris en compte dans 20 les calculs de la manière suivante : Tout d'abord, en ce qui concerne la deuxième face, le flux thermique surfacique OL(t) échangé est utilisé comme paramètre pour la détermination du champ de température 0(x,t). Les contraintes quant aux échanges thermiques au niveau de la deuxième face sont prises en compte via l'équation citée précédemment 25 (S3) C L, + LO(L = t,{t} [Q,T] Cette équation signifie que le flux thermique qui est échangé par la paroi en x=L 30 (membre de gauche) est égal au flux thermique échangé avec l'extérieur (cDL(t), membre de droite). Le cas échéant, après avoir déterminé l'évolution des températures dans la paroi, cette équation permet de calculer la valeur du flux thermique surfacique de la deuxième face. D'autre part, en ce qui concerne la première face, le troisième capteur (Aamb) fournit des mesures de la température au voisinage de la première face. A l'aide des informations fournies par ce capteur, ainsi que par le capteur A0, les moyens de traitement prennent en compte l'échange thermique sur la première face et calculent les valeurs de température (6(x,t)) en prenant en compte l'équation suivante, citée précédemment : 00(0t) (S2) k; o z, + aoO(O, t) = a - T (t) [0,T] Par ailleurs, dans la détermination du modèle du champ de température (8(x,t)), un champ initial de température (8i(x) dans la paroi à l'instant initial est utilisé comme paramètre. En résumé, l'évolution de la température dans la paroi est modélisée par un système d'équations (1) aux dérivées partielles de type parabolique : (S1) 00(0, t) (S2) - (:.:toO«), t) = c~~oT(t) [0, T] 0x (S3) 09(L, t) c;H(L, t) = (t) [0, T] (S4) °G 0) _ i (x) 25 30 Dans ce modèle (1) : La première équation (S1) dite équation de la chaleur régit l'évolution de la température (8(x,t)) à l'intérieur de la paroi. Cette évolution dépend essentiellement de la conductivité k et de la capacité thermique C de la paroi. Les équations S2 et S3 (1) correspondent aux conditions d'échange thermique sur les faces FO et FL. L'équation S4 impose que le champ de température dans la paroi à l'instant initial t=0 des mesures, soit égal à un champ de température donné 0i(x). D'autres formulations peuvent être adoptées notamment pour les équations S2 et S3. L'équation S3 présente un modèle d'échange de chaleur générique sur la paroi FL. En effet, le flux thermique surfacique échangé OL(t) peut être absolument quelconque. L'équation S2 présente un modèle d'échange de température simplifié sur la paroi F0, selon l'hypothèse que cette face ne reçoit pas de chaleur sous forme de rayonnement. A la place de l'équation S2, il serait possible pour la paroi FO d'utiliser une équation similaire à l'équation S3, en remplaçant simplement les références à la paroi FL, par des références à la paroi F0. De manière connue en soi, le système d'équations (1) est résolu grâce à une discrétisation des inconnues et des paramètres d'entrée et d'un algorithme d'intégration numérique, en supposant connu l'ensemble des paramètres inconnus U. La méthode des éléments finis ou la méthode des différences finies est utilisée pour l'approximation dans l'espace ce qui résulte à un système matriciel dépendant du temps. Le domaine temporel est discrétisé en pas de temps de longueur égale, et le champ de température est obtenu, pour chaque pas de temps, grâce à l'algorithme implicite d'Euler, en utilisant les valeurs du champ de température au pas de temps précédent. Par ailleurs, dans un mode de réalisation non présenté, le modèle pourrait prendre en compte en outre deux équations supplémentaires, qui traduisent la contrainte d'égalité entre les température 0(0,t) et 0(L,t) telles que fournies par le modèle pour les deux surfaces de la paroi, et les températures Od1(t), 0d2(t) effectivement mesurées sur l'une et l'autre de ces surfaces. Dans le mode de réalisation présenté inversement, cette contrainte d'égalité entre les températures calculées et mesurées sur les deux faces de la paroi n'est pas prise en compte parmi les équations du système (S) comme une contrainte du modèle, à l'instar des équations S2 à S4. En effet, cette contrainte d'égalité est prise en compte par un critère d'erreur J, que la méthode vise à minimiser. Lorsque le minimum est atteint, la contrainte d'égalité est alors satisfaite ou quasiment satisfaite. Après avoir présenté le modèle numérique sur lequel reposent les calculs effectués, l'algorithme mis en oeuvre à l'aide des modules 60 et 70 va maintenant être précisé : Le module de simulation numérique permet de calculer des valeurs du champ de température, à partir des variables du vecteur U et de la mesure fournie par le troisème capteur Aamb. Le programme de calcul est itératif. Lors de la première itération, des valeurs arbitraires par défaut sont utilisées comme vecteur U, et permettent un premier calcul des températures 0(x,t) dans la paroi. Ce calcul du champ de température 0 fournit en particulier les valeurs de température sur les deux faces, 0(0,t) et 0(L,t). Ces valeurs, ainsi que les mesures de températures Odi(t), 0d2(t) sont fournies à un comparateur 65, qui détermine la valeur d'un critère d'erreur J. Ce critère d'erreur est donné par l'équation : J= (e(xi, t) - o (t))2 dt dans laquelle on utilise la convention suivante : xi = 0 et x2 = L. Le critère J mesure l'écart entre les valeurs du modèle de température 0(O,t) et 0(L,t) respectivement sur les première et deuxième faces et les températures 20 Odi(t), 0d2(t) mesurées respectivement par les premier et deuxième capteurs. La minimisation du critère J est faite en faisant varier les valeurs des différents paramètres du vecteur U, au moyen du module d'adaptation 70. Ce module permet de déterminer un nouveau jeu de valeurs des paramètres du vecteur U, de 25 manière à réduire la valeur du critère d'erreur J, dans l'objectif de minimiser celui-ci. En sortie, le module d'adaptation 70 fournit ainsi de nouvelles valeurs pour les variables suivantes (ou seulement certaines d'entre elles, en fonction du mode de mise en oeuvre) : la conductivité (k) et/ou la capacité thermique (C) de la paroi, le flux thermique surfacique de la deuxième face ((DL(t)), le champ initial de 30 température Oi(x) dans la paroi, et un ou des coefficient(s) (aO,aL) d'échange thermique liant le flux thermique surfacique de la première et la deuxième faces respectivement à la température en surface de la paroi. Pour minimiser le critère J, les modules 60 et 70 de simulation numérique et d'adaptation sont mis en oeuvre en utilisant l'algorithme de Levenberg-Marquardt. Le modèle présenté précédemment comporte à la fois une inconnue principale, le champ de température 8(x,t) dans la paroi, et un ensemble de paramètres inconnus, regroupés en un vecteur de paramètres U égal à { cl)L(t), 6i(x), k, C, a0, aL}. La méthode proposée permet d'identifier l'ensemble de ces inconnues, en déterminant des valeurs de A(x,t) et de U qui font coïncider sensiblement les températures mesurées avec celles provenant de la simulation, c'est-à-dire les températures de surface calculées par le module de simulation 60. Avantageusement, cette identification est réalisée simplement à partir des mesures de température edi(t), 8d2(t) sur les deux faces de la paroi, et de la mesure de la température de l'air dans le voisinage de la première face ou de la quantité de chaleur totale reçue par cette face. La détermination des paramètres inconnus U est faite en utilisant un algorithme de minimisation itératif. Celui-ci consiste à chaque pas de temps à fixer une valeur pour l'ensemble des paramètres inconnus U, calculer la valeur du champ de température au niveau de chaque face résultant de cet ensemble de paramètres, et modifier en conséquence les valeurs de l'ensemble de paramètres U jusqu'à ce que la fonction J soit minimale. Le problème posé étant fortement non-linéaire, la méthode détaillée ici recourt à la linéarisation et notamment à l'algorithme de Levenberg-Marquardt. Cet algorithme apporte une vitesse de convergence plus élevée que d'autres algorithmes, tout en restant de complexité raisonnable. Plus précisément, le programme informatique est fondé principalement sur l'algorithme de Levenberg-Marquardt, et a recours à la méthode du gradient conjugué en tant que sous-processus. Suivant cette méthode, à chaque itération n, ayant déterminé le champ de température 8(x,t) dans la paroi en fonction du temps pour le vecteur Un (Un désignant la valeur attribuée au vecteur U lors de l'itération n), on détermine le vecteur Un+1 en minimisant un critère d'erreur linéarisé J'(U). Ce critère d'erreur linéarisé J' défini par linéarisation du critère J peut s'exprimer par l'équation (2) suivante (2) T`(ir..) = 2 J (H(xj, t) + 88(r.i, t) e(t))2 Olt -I u..r,.11v _, Jc) où 80 = 80(8U) est la solution du système (1) linéarisé autour d'un point Un = { (DL(t), 0i(x), kn,Cn, ao,n, ai.,n, }- (on note SU = { (1)L(t), 0i(x), Sk, SC, SaO, SaL}). Le dernier terme sur la droite dans l'équation (2) est un terme de régularisation de type Tikhonov autour du point de régularisation Un, et c est un paramètre de valeur proche de zéro. Il - Il représente une mesure des valeurs de l'ensemble U, est peut être donné par la racine carrée de la somme des carrés des valeurs des paramètres, intégrés dans le domaine de définition si les paramètres sont des fonctions d'une variable. D'autre part, 80 = 80(8U) est la solution du système (3) obtenu par linéarisation du système (1) et qui s'écrit : 060 0 060 060 . . 00 rl Ot O 1' -(C5l ()f + (~~a~ «5k - h,lt l ox 060(0) 00(0» / -kn + ao ,ibO(0) (nk - k,) - «Sao - aom.)(8(0) (3) Ox; Oz ( car>8(L) - 00(L) n +CkL,r,ôo(L) - _«5k kii) i)x - (f~Q'f, ~kL,°n.)O(L) (Îti 60(x, 0) = 0 L'algorithme itératif consiste donc à procéder à un premier niveau d'itérations (suivies par une variable n). Lors de chacune de ces itérations, on détermine le champ de température 0(x,t) dans la paroi en fonction du temps pour le vecteur Un, puis on détermine le vecteur Un+1 qui minimise le critère d'erreur linéarisé Y(Un). Pour cette dernière opération de minimisation, qui fournit Un+l, l'algorithme fait appel à un sous-processus de minimisation, lui-même itératif en fonction d'une variable I, et qui utilise la méthode du gradient conjugué. Ainsi dans ce mode de réalisation, la minimisation du critère d'erreur linéarisé est réalisée par une méthode de gradient conjugué. Ce sous-processus nécessite de calculer à chaque itération (sur 1) le gradient de la fonctionnelle J'(U). Ce gradient est calculé en calculant l'état adjoint p, solution du problème (4) suivant : YOp t) Off} c---Kr Kr + Mer,. Op«)) (4) -k + aop«..)) - 0 o x k° Op(L) + (1t11"}(L) = Q O x 15 où l'on note Sxj la fonction Dirac qui vaut 1 en xj et 0 ailleurs. Dans l'algorithme de gradient conjugué, la fonction de descente est une fonction de 0 et de p, notée 20 F(8,p). Récapitulatif Au final, l'algorithme d'identification simultanée du champ de température dans la paroi, et du vecteur de paramètres U, qui vaut { (DL(t), 0i(x) , k,C, a0, aL}, est 25 conduit de la manière suivante (Uni désigne la valeur attribuée au vecteur U lors de l'itération n): - Initialisation : Uol - Boucle sur n (algorithme de Levenberg-Marquardt) : A chaque itération 30 (n>=1) : o Initialisation :1=0, Uno=Un_11+l o Boucle sur 1 (algorithme du gradient conjugué) : A chaque itération (I>1) :10 ^ Calculer 0, solution de (1) : on détermine le champ de température correspondant le mieux au vecteur Un courant ; ^ Calculer 80 solution de (3) ^ Calculer p solution de (4) ^ Calculer la direction de descente d = F(0,p) ^ Calculer le pas optimal p tel que J'( Uni+ p d)= min r J'(Uni+ r d) ^ Mettre à jour les variables : Uni+l=Uni +p d : on détermine le vecteur Uni+1 optimisé en fonction du champ de température identifié précédemment ; ^ Si II Uni+1- Uni 11 REVENDICATIONS1. Dispositif (10) d'identification de la conductivité k et/ou de la capacité thermique C d'une paroi plane s'étendant entre une première (FO) et une deuxième (FL) face (x=O,x=L), le dispositif comportant : un premier capteur (AO) de température apte à être placé sur la première face pour acquérir des mesures de la température (0d1(t)) de celle-ci ; un deuxième capteur (AL) de température apte à être placé sur la deuxième face pour acquérir des mesures de la température (0d2(t)) de celle-ci ; un système d'acquisition de données, apte à acquérir périodiquement des mesures fournies par lesdits capteurs ; et des moyens de traitement des données aptes à déterminer le coefficient de conductivité k et/ou la capacité thermique C de la paroi, le dispositif se caractérisant en ce qu'il comporte en outre un troisième capteur (Aan,b), qui est un capteur de température apte à être disposé au voisinage de la première face pour fournir une mesure de la température ou un capteur de flux thermique surfacique apte à être placé sur la première face de la paroi pour mesurer le flux thermique surfacique (00) de la première face ; et que les moyens de traitement déterminent le coefficient de conductivité k et/ou la capacité thermique C, en calculant l'évolution des températures dans la paroi pendant un intervalle de temps, à partir notamment des mesures fournies par ledit troisième capteur pendant ledit intervalle, et en retenant pour le coefficient de conductivité k et/ou la capacité thermique C, la ou les valeurs pour lesquelles, durant ledit intervalle, les différences entre les valeurs de température (0(0,t), 0(L,t)) calculées pour les première et deuxième faces et les températures (0d1(t), Od2(t)) mesurées respectivement par les premier et deuxième capteurs sont minimales. 2. Dispositif selon la revendication 1, dans lequel les moyens de traitement sont 30 aptes à exécuter : - une phase dinitialisation durant laquelle des valeurs approchées des températures dans la paroi pendant ledit intervalle de temps sont calculées ; - une phase de calcul durant laquelle le coefficient de conductivité k et/ou la capacité thermique C sont calculés en minimisant lesdites différences, seules 35 étant utilisées dans cette phase les mesures acquises par lesdits premier, deuxième et troisième capteurs. 3. Dispositif selon la revendication 1 ou 2, dans lequel pour déterminer le coefficient de conductivité k et/ou la capacité thermique C, les moyens de 40 traitement comportent - un module de simulation, apte à calculer les valeurs de température de la paroi (0(x,t)) pendant ledit intervalle, en fonction notamment de la conductivité (k) et/ou de la capacité thermique (C) de la paroi ; et- un module d'adaptation, apte à déterminer une ou des valeurs modifiée(s) d'au moins une variable d'entrée du module de simulation, dont notamment le coefficient de conductivité k et/ou la capacité thermique C, de manière à réduire ou minimiser lesdites différences. 4. Dispositif selon la revendication 3, dans lequel les variables que le module d'adaptation est apte à modifier comprennent la conductivité (k) et/ou la capacité thermique (C) de la paroi, et éventuellement un ou plusieurs des paramètres suivants : flux thermique surfacique de la deuxième face ($:DL(t)), champ initial de température 8i(x) dans la paroi, et un ou des coefficient(s) (aO,aL) d'échange thermique liant le flux thermique surfacique de la première et/ou la deuxième face à la température à l'intérieur ou en surface de la paroi. 5. Dispositif selon l'une quelconque des revendications 2 à 4, dans lequel les 15 moyens de traitement sont aptes à exécuter : - une phase d initialisation durant laquelle des valeurs approchées des températures dans la paroi pendant ledit intervalle de temps sont calculées ; - une phase de calcul durant laquelle le coefficient de conductivité k et/ou la capacité thermique C sont calculés en minimisant lesdites différences, seules 20 étant utilisées dans cette phase les mesures acquises par lesdits premier, deuxième et troisième capteurs ; le dispositif comprend un quatrième capteur (Ae) apte à être disposé au voisinage de ladite deuxième face pour fournir des valeurs approchées du flux surfacique (00 échangé par ladite deuxième face ; et 25 durant la phase d initialisation, les valeurs approchées des températures dans la paroi sont calculées à partir notamment de valeurs approchées du flux thermique surfacique (00 mesurées par ledit quatrième capteur pendant ledit intervalle de temps. 30 6. Dispositif selon l'une quelconque des revendications 1 à 5, dans lequel pour déterminer les valeurs minimales desdites différences, les moyens de traitement exécutent un programme informatique qui minimise un critère d'erreur 0). 7. Dispositif selon l'une quelconque des revendications 1 à 6, comportant en 35 outre au moins une source de chaleur (32), apte à être disposée en regard d'une face de ladite paroi et à solliciter celle-ci par convection ou par rayonnement, modifiant ainsi le flux thermique surfacique ((DL) de ladite face. 8. Méthode d'identification de la conductivité k et/ou de la capacité thermique C 40 d'une paroi plane s'étendant entre une première (FO) et une deuxième (FL) face (x=O,x=L), la méthode comportant les étapes suivantes : a) on dispose :sur une première face de la paroi, un premier capteur (Al) de température, afin d'acquérir des mesures de la température (8d1(t)) de ladite première face ; sur une deuxième face de la paroi, un deuxième capteur (A2) de température, afin d'acquérir des mesures de la température (8d2(t)) de ladite deuxième face ; un troisième capteur (Aan,b), qui est soit un capteur de température disposé au voisinage de la première face pour fournir une mesure de la température, soit un capteur de flux thermique surfacique placé sur la première face de la paroi pour mesurer le flux thermique surfacique ((D0) de cette première face ; b) pendant un intervalle de temps, on acquiert périodiquement les températures des première et deuxième faces de la paroi, ainsi que les mesures fournies par le troisième capteur; puis, c) on détermine le coefficient de conductivité k et/ou la capacité thermique C, en calculant les températures dans la paroi pendant ledit intervalle de temps, à partir notamment des mesures fournies par ledit troisième capteur pendant ledit intervalle ; et en retenant pour le coefficient de conductivité k et/ou la capacité thermique C, la valeur pour laquelle ou les valeurs pour lesquelles, durant ledit intervalle, les différences entre les valeurs de température (8(0,t), 8(L,t)) calculées pour les première et deuxième faces et les températures (MW), 8d2(t)) mesurées respectivement par les premier et deuxième capteurs sont minimales. 9. Méthode selon la revendication 8, dans laquelle l'étape c) comporte : c1) une phase d initialisation pendant laquelle on détermine des valeurs approchées des températures dans la paroi pendant ledit intervalle ; puis c2) une phase de calcul durant laquelle on détermine le coefficient de conductivité k et/ou la capacité thermique C en minimisant lesdites différences, seules étant utilisées dans ce calcul les mesures acquises par lesdits premier, deuxième et troisième capteurs. 10. Méthode selon la revendication 8 ou 9, dans laquelle pendant ledit intervalle de temps, la première face de la paroi est protégée des rayonnements thermiques par un écran opaque aux radiations. 11. Méthode selon l'une quelconque des revendications 8 à 10, dans laquelle pendant ledit intervalle de temps, la paroi n'est pas soumise à l'action d'une source de chaleur contrôlée. 12. Méthode selon l'une quelconque des revendications 8 à 10, dans laquelle durant au moins une partie dudit intervalle de temps et/ou pendant une période précédant ledit intervalle, on dispose une source de chaleur contrôlée, en regard40d'une face de ladite paroi pour modifier le flux thermique surfacique de ladite face. 13. Méthode selon l'une quelconque des revendications 8 à 12, dans laquelle pendant au moins une partie dudit intervalle de temps, le flux thermique surfacique sur la deuxième face de la paroi est mesuré au moyen d'un capteur de flux thermique.