FR 3048509 A1 20170908 FR 1651774 A 20160302 1 La présente invention concerne un procédé et un dispositif de mesure de propriétés thermo-physiques d'un milieu. Une méthode connue pour mesurer les propriétés thermo-physiques d'un milieu est la méthode dite du fil chaud modulé, ou « méthode 3w », qui consiste à exciter à une fréquence w un fil métallique de faible diamètre, de l'ordre de quelques dizaines de la- mètres, dont la longueur peut varier de quelques millimètres à quelques dizaines de millimètres et dont la résistivité électrique en am -1 et le coefficient de température en K-1 sont connus. Une telle méthode est notamment décrite dans l'article de C. Dames and G. Chen, « 1, 2, and 3 omega methods for measurements of thermal properties », Review of Scientific Instruments, n°76, 124902, 2005. L'excitation thermique par effet Joule apparait à la fréquence double 2e) et génère une variation de température du fil, et par conséquent de la résistance du fil. La tension, produit du courant de fréquence w et de la résistance à 2e) se compose majoritairement d'un terme dit « parasite » à lio et d'une petite composante utile à 3e). Pour éliminer la composante majoritaire, le fil est placé dans un pont dont on équilibre les deux branches et un voltmètre à détection synchrone effectue la mesure de la composante à 3e). Cette mesure permet alors de déterminer, par comparaison avec une référence, la conductivité thermique du milieu dans lequel se trouve le fil, quelle que soit la nature du milieu, par exemple gaz, liquide ou solide. Lorsque l'on souhaite effectuer des mesures en continu sur une large gamme de températures, la résistance du fil ne cesse d'évoluer avec la température du milieu ambiant et l'équilibrage du pont doit être réalisé en continu. Cette méthode, dite du « lock-in », a une capacité limitée à séparer correctement la composante 3w en présence d'autres harmoniques. En particulier, l'amplitude du fondamental à 1e) ne doit pas dépasser de plus d'un facteur 10 l'amplitude de l'harmonique à 3e), et pour une bonne précision de la mesure de la phase, les deux amplitudes devraient être proches. Avec un pont avec équilibrage à température ambiante, le déséquilibre avec la montée en température peut résulter en un rapport d'amplitudes 1 e)/3w égal à 1000. Un filtrage, notamment par un filtre étroit (« notch filter » en anglais), peut être utilisé pour diminuer la composante parasite, mais cela affecte la phase de la composante utile. Des constantes de temps supérieures à au moins trente fois la période de modulation peuvent être utilisées, mais le temps d'acquisition du signal est alors fortement augmenté. 3048509 2 Le principe de construire un pont actif avec un seul fil et un transistor FET à la place du potentiomètre est connu pour assurer automatiquement l'équilibre du pont quand la résistance varie en fonction de la température, grâce à une réaction. Mais un tel réglage automatique du pont génère des fluctuations qui gênent encore plus la mesure des 5 composantes à 3w de très faible amplitude. Par ailleurs, les méthodes de mesure des propriétés thermiques en régime stationnaire ne sont pas adaptées aux fluides du fait de la présence de courant de convection naturelle lors de l'établissement du gradient thermique. En régime non stationnaire, la méthode du fil chaud, décrite notamment dans le livre de Lomas, C.G. 10 « Fundamentals of hot wire anemometry », Cambridge Univ. Press, 2011, et les méthodes en régime périodique, décrites notamment dans l'article de P. Bhattacharya et al. «Characterization of the temperature oscillation technique to measure the thermal conductivity of fluids », Int. Journal of Heat and Mass Transfer, n°49, pages 2950-2956, 2006, sont plus adaptées et couramment utilisées. 15 Une autre méthode utilisant les oscillations thermiques et la diffusion d'ondes thermiques a été proposée par A. J. Angstrom en 1863, et a ensuite été utilisée pour la mesure de la diffusivité thermique de fluides, comme expliqué dans l'article d'A. Santucci, L. Verdini, and P.G. Verdini. «Data-acquisition system for measurement of thermal diffusivity and propagation properties of thermal waves by a non-steady-state method », 20 Review of Scientific Instruments, n°57, pages 1627-1632, 1986. La conductivité thermique, la diffusivité thermique ou les deux paramètres à la fois sont mesurés en géométrie cylindrique soumise à un flux thermique longitudinal périodique. Des dispositifs à double fil chaud, un long, l'autre court, parcourus par le même courant, sont également connus, dans le but de compenser par soustraction les effets 25 de bords dus à la longueur finie du fil. Il subsiste par conséquent un besoin pour améliorer encore les méthodes de mesure de propriétés thermo-physiques d'un milieu, afin de remédier aux inconvénients énoncés ci-dessus. L'invention a pour but de répondre à ce besoin et elle y parvient, selon l'un de 30 ses aspects, grâce à un procédé de mesure de propriétés thermo-physiques d'un milieu dans lequel sont placées une sonde de mesure et une sonde de compensation, la sonde de mesure étant parcourue par un courant et la sonde de compensation par un courant de même 3048509 3 fréquence et d'amplitude proportionnelle à celle du courant parcourant la sonde de mesure, procédé dans lequel un signal contenant une composante harmonique utile pour la mesure des propriétés thermo-physiques du milieu est généré à partir d'une amplification différentielle des tensions aux bornes de la sonde de compensation et de la sonde de 5 mesure, avec un gain choisi pour diminuer l'amplitude de la première composante harmonique au sein du signal résultant. La composante harmonique utile est de préférence la troisième harmonique, de fréquence 3w. L'invention porte encore, selon un autre de ses aspects, sur un dispositif de 10 mesure de propriétés thermo-physiques d'un milieu, notamment pour la mise en oeuvre du procédé de mesure défini ci-dessus, comportant une sonde de mesure et une sonde de compensation placées dans le milieu à mesurer, la sonde de mesure étant parcourue par un courant et la sonde de compensation par un courant de même fréquence et d'amplitude proportionnelle à celle du courant parcourant la sonde de mesure, le dispositif étant 15 configuré de telle sorte qu'un signal contenant une composante harmonique utile pour la mesure des propriétés thermo-physiques du milieu soit généré à partir d'une amplification différentielle des tensions aux bornes de la sonde de compensation et de la sonde de mesure, avec un gain choisi pour diminuer l'amplitude de la première composante harmonique au sein du signal résultant. 20 L'avantage de l'invention est d'une part d'exploiter la grande précision et sensibilité du traitement du signal en régime modulé périodique et d'autre part d'offrir une méthode et une solution technique implémentée sur une sonde à double fil, pour la compensation de la dérive en température de l'équilibre du pont de mesure. En l'absence de cet équilibrage, le traitement de type lock-in du signal utile est affecté par des erreurs de 25 mesure systématiques, comme expliqué précédemment. L'invention permet d'obtenir un «pont actif» amélioré et équilibré en continu lors des variations de température car les deux sondes, correspondant chacune à une branche du pont, sont parcourues par des courants parfaitement connus et subissent les mêmes variations thermiques. 30 L'invention permet la séparation des harmoniques grâce à l'asymétrie des courants et l'amplification définie sur l'une des deux sondes. 3048509 4 L'invention permet d'obtenir des courbes stables des propriétés thermiques qui varient en fonction de la température. Le dispositif de mesure selon l'invention est beaucoup plus simple à mettre en oeuvre et à utiliser que les dispositifs connus, ce qui se traduit par un coût de construction 5 réduit. Une composante continue peut être ajoutée. Dans ce cas, la composante harmonique utile est la deuxième harmonique, de fréquence 2(o. Génération de la composante harmonique utile Le dispositif selon l'invention comporte de préférence un amplificateur, de 10 type « lock-in », comportant un générateur interne générant une tension de référence Vie, les amplitudes des courants parcourant les sondes étant de préférence proportionnelles à ladite tension de référence Via. Le rapport N entre le courant parcourant la sonde de compensation et celui I. parcourant la sonde de mesure est avantageusement inférieur à 1. Ce rapport peut être 15 choisi en fonction de la précision souhaitée pour les mesures. Les sondes peuvent être parcourues par des courants provenant de générateurs de courant, configurés de telle sorte que les amplitudes des courants soient différentes et leurs fréquences identiques. Dans une variante, les sondes sont montées en série avec des résistances 20 d'ajustement de la tension, choisies de telle sorte que, une même tension étant appliquée à leurs bornes, les amplitudes des courants soient différentes et leurs fréquences identiques. Les coefficients de variation avec la température des résistances internes des sondes sont avantageusement identiques. Les résistances internes des sondes peuvent être identiques. 25 Les calculs suivants sont établis dans l'hypothèse où aucune composante continue n'est ajoutée. L'amplitude de la tension aux bornes de la sonde de mesure, mesurée en un point A, peut s'exprimer comme suit : V A=[Vh,±V3,,(T)]A=11,,RA(T)+ K(T)113,,RA(T), 30 où le courant 13. généré par la non-linéarité de la résistance de la sonde dépend du courant d'excitation Iia : I (T) = K (T)I , avec K(T) le coefficient de non-linéarité, et II. = Im le courant parcourant la sonde de mesure. 3048509 5 L'amplitude de la tension aux bornes de la sonde de compensation, mesurée en un point B', peut s'exprimer comme suit: VB, =[171,,±V 3,,(T)]B= NI u,RB(T) + K (T )(NI 1,, )3 RB (T ) , avec N11,) = le courant parcourant la sonde de compensation. 5 La tension VB' aux bornes de la sonde de compensation est avantageusement amplifiée d'un gain G, par passage dans un amplificateur opérationnel, l'amplitude de la tension amplifiée, mesurée en un point B, s'exprimant de préférence comme suit : VB =GV B, =GNI u,RB(T)+ GK(T)(N11)3 RB(T), où G est le gain de l'amplificateur opérationnel, avec G > 1. 10 Le gain G d'amplification de la tension aux bornes de la sonde de compensation est avantageusement égal à l'inverse du rapport N entre le courant parcourant la sonde de compensation et celui parcourant la sonde de mesure. Cela permet l'annulation de la première composante harmonique à 1w. L'amplification différentielle des tensions aux bornes de la sonde de 15 compensation et de la sonde de mesure, relevée entre le point A et le point B, après amplification par le gain G, peut s'exprimer comme suit : VA-VB = (1- GN)11R(T)+ (1 - GN3 )K (T)IR(T) , en considérant que les résistances des deux sondes sont égales, soit RA(/)=RB(/)=R(/). Le gain G d'amplification de la tension aux bornes de la sonde de 20 compensation étant égal à l'inverse du rapport N entre le courant parcourant la sonde de compensation et celui parcourant la sonde de mesure, l'amplification différentielle s'exprime avantageusement comme suit : VA-VB = (1- N2 )K (T)113'R(T) = (1- N2 )V3' (T) . La première composante harmonique à 1w est ainsi complètement supprimée 25 ou largement diminuée, tandis que la composante harmonique à 3w a peu perdu en amplitude, par exemple seulement 1 % avec G = 1/N = 10. La fréquence w peut être comprise entre 0.5 Hz et 10 Hz, mieux entre 1 Hz et 5 Hz, notamment pour un fil de 40 jam de diamètre. 30 3048509 6 Lien entre composante harmonique utile et température La température du milieu peut être calculée en fonction de la tension aux bornes d'une des deux sondes. Une sonde de température séparée n'est ainsi pas nécessaire. 5 L'évolution de la température d'une sonde peut s'exprimer en fonction de la puissance dissipée : T(t) = 1 f P (t)dt = 1 ° R ° t sin(2wt) 2 , avec m= pSL la masse de la mc mc sonde, c la capacité thermique massique en J.Kg-1.1(-1, et Ro la résistance électrique de la sonde à température ambiante en absence du courant. La puissance dissipée est fonction du courant transmettant l'excitation électrique à la sonde : P(t)=12 (OR° avec /(t)=Tocos(ot). 10 Par effet Joule, la température T(t) peut avoir une composante harmonique à 2w avec une amplitude égale à : T2.-2/2effRo /(4wmc). La résistance de la sonde peut osciller à la fréquence 2w : R(t) = Ro[l+r T 2.cos(2 ot- 9)] . Par la loi d'Ohm, en fonction de la température de la sonde, la tension aux 15 bornes de la résistance de la sonde peut s'exprimer comme suit : T7(t) = I(t)R(t) = IoRo{cos(wt)+(1/2) rT2w[cos(ot-9)+cos(3ot-9)] 1, où r est le coefficient de variation de la résistance avec la température de la sonde. Le terme à 3 w apparaît à cause du mixage du courant alternatif à la fréquence w avec l'oscillation de la résistance à 2w. Cette équation montre que la composante à 3w est 20 proportionnelle à T2co, l'oscillation de la température à la pulsation 2w. Cela sert de base aux mesures de propriétés thermo-physiques par la méthode du fil chaud à 3w. Le premier terme de l'équation ci-dessus est bien plus grand que le dernier terme, et peut être annulé en pratique en plaçant la sonde dans un pont de Wheatstone. Le ratio entre la troisième composante harmonique et le fondamental de la 25 tension peut s'exprimer ainsi : V.30,47t. (1/4)rTpp rI2 ef f Ro /(4wmc), avec Tpp = 2T2co l'oscillation crête-à-crête de la température de la sonde dans l'air. Ce ratio est proportionnel à la puissance dissipée. Dans le cas d'une excitation de la sonde uniquement par courant alternatif, la 30 température de la sonde oscille à 2w, alors que la tension aux bornes de la sonde a 3048509 7 uniquement comme composantes celles à 10 et à 3o, comme décrit dans l'article de Lefèvre S. et Volz S., « 3o-scanning thermal microscope », Review of Scientific Instruments, 2005, 76 (3), 033701. Dans le cas où la sonde est excitée à la fois par un courant alternatif et par un courant continu, la tension résultante aux bornes de la sonde a 5 trois composantes harmoniques, 1w, 2w et 3w, comme mentionné dans l'article de Dames C. et Chen G., « 1, 2 and 3w methods for measurement of thermal properties », Review of Scientific Instruments, no. 76, 124902, 2005. Diffusivité thermique et conductivité thermique Le lien entre la température T2 d'une sonde et la diffusivité thermique et la 10 conductivité thermique du milieu peut s'exprimer comme suit : 7.c Pli ( + ln as r Pli ( 1.2594r 7r Pli ( ln a . T = 0(f)= = ln + / - 27eks y 2 27eks 1-ts 4 ) 47ck 4.9828 r2 .2f 2 . 0 ) avec 0(f) la réponse transitoire de la température, a la diffusivité thermique, k la conductivité thermique, ro le rayon du fil de la sonde, P/l la puissance dissipée par unité de longueur, C = exp(y) = 1,7810, et bis la longueur de diffusion thermique dans la direction 15 radiale ,us = la/(7r.2fill/ 2 . En utilisant la partie réelle de la réponse transitoire de la température, Pli( 1 , on peut obtenir la conductivité thermique Re[0 (f )] = ln a - ln(2f)+ ln 47rk 4.9828 ro ) en utilisant la pente de la partie réelle de 0(f ) en fonction de ln(!). En utilisant la phase 9 de la réponse transitoire de la température, 20 - cot(9) = ln(2f)- ln a + ln(4.9828 r02) , on peut obtenir la diffusivité thermique à partir de cette phase à une fréquence fixée, par exemple f = 1 Hz. La définition de la diffusivité thermique est a = k/(pc), avec p la masse volumique et c la capacité thermique massique en J.Kg-1.K-1. Leur produit correspond à la capacité thermique volumique en J.m-3.1(-1 : pc=k/a. 25 Grâce à l'invention, il est possible d'obtenir des valeurs indépendantes, simultanées et absolues de la diffusivité thermique et de la conductivité thermique d'un milieu, ainsi que de la capacité thermique volumique. 3048509 8 Sondes Les sondes peuvent être jetables, à usage unique ou non. Les sondes peuvent comporter un fil, notamment de longueur comprise entre 0.1 mm et 10 cm, par exemple égale à 4 cm, et notamment de diamètre compris entre 5 iam 5 et 200 iam, notamment égal à 40 iam. Le fil est avantageusement métallique, étant notamment réalisé en platine, tantale, nickel ou d'alliages à faible coefficient de dilatation thermique, par exemple un alliage de fer, de nickel et de cobalt. Le fil est de préférence lisse, ayant une surface parfaite. 10 Il peut être recouvert par une fine couche isolante électrique d'oxyde, de polymère ou de résine époxy, adaptée pour des mesures dans des liquides conducteurs électriques. Le fil est de préférence réalisé en un matériau ayant un coefficient de variation de la résistance avec la température élevé, notamment égal environ à 0.005 K-1. 15 Description détaillée L'invention pourra être mieux comprise à la lecture de la description détaillée qui va suivre, d'exemples de mise en oeuvre non limitatifs de celle-ci, et à l'examen du dessin annexé, sur lequel : - la figure 1 représente un dispositif de mesure de propriétés thermo- 20 physiques d'un milieu selon l'invention, et - les figures 2 à 4 représentent des chronogrammes de signaux obtenus en appliquant le procédé selon l'invention. Un exemple de dispositif 1 destiné à la mise en oeuvre du procédé de mesure de propriétés thermo-physiques d'un milieu selon l'invention, est représenté à la figure 1. Une 25 sonde de mesure 2 et une sonde de compensation 3 sont immergées dans un milieu à mesurer, à la même température T. La sonde de mesure 2 est parcourue par un courant I. de pulsation w et la sonde de compensation 3 est parcourue par un courant de même fréquence w et d'amplitude proportionnelle à celle du courant L. De préférence et comme dans l'exemple 30 décrit, le rapport N entre le courant parcourant la sonde de compensation 3 et celui In, parcourant la sonde de mesure 2 est inférieur à 1 : N = 1. Comme expliqué précédemment et comme visible à la figure 1, l'amplification 30 différentielle des tensions aux bornes de la sonde de compensation 3 et de la sonde de mesure 2, relevée entre le point A où est mesurée la tension aux bornes de la sonde de mesure 2 et le point B, après amplification par le gain G, s'exprime comme suit : 3048509 10 VA-VB=(1-N2)K(T)113'R(T)= (1- N2)V3'(T) . La première composante harmonique à 1w est ainsi complètement supprimée. Comme décrit précédemment, la température du milieu peut être calculée en fonction de la tension aux bornes d'une des deux sondes. 5 Les figures 2 à 4 représentent l'évolution de signaux générés selon l'invention : l'amplitude de la fréquence fondamentale à 1e) est représentée à la figure 2, avec w = 1 Hz, et l'amplitude et la phase de la troisième composante harmonique à ko sont représentées aux figures 3 et 4. Les chronogrammes représentent les courbes obtenues selon l'invention, et celles obtenues, selon la méthode connue, avec un pont passif contenant trois résistances 10 et une seule sonde comportant un fil résistif. De préférence et comme dans l'exemple décrit, dans le dispositif 1 selon l'invention, les sondes de mesure 2 et de compensation 3 sont identiques et comportent chacune un fil, par exemple de longueur égale à 4 cm et de diamètre égal à 40 p.m. Le fil est avantageusement métallique, étant réalisé en nickel dans l'exemple considéré, ayant 15 une résistance interne de 2,2 ohm. Dans l'exemple considéré, la ou les deux sondes ont été immergées dans de l'eau à température ambiante. Le pont a été équilibré en s'assurant d'une amplitude égale à 0 V entre les points A et B. Puis la ou les deux sondes ont été immergées dans de l'eau à 60°C et simultanément l'enregistrement des courbes a démarré. Avec le temps, la 20 température de l'eau diminue jusqu'à la température ambiante, et les courbes tendent vers zéro. La figure 2 montre le déséquilibre initial de tension, à t = 0 s, qui est de 18% pour le dispositif selon l'art antérieur et 1,9% pour le dispositif selon l'invention. L'efficacité de ce dernier correspond au facteur de réduction de la dérive thermique : 18/1,9 25 = 9,5. Cela montre que la dérive de la première composante harmonique est suffisamment faible pour ne pas perturber la mesure de la troisième composante harmonique utile. Dans l'exemple considéré, le dispositif selon l'invention a réduit de 10 fois la dérive en température. Cette réduction est suffisante pour permettre de réaliser des mesures exploitables de la composante harmonique utile, sans l'intervention de l'opérateur pour 30 réajuster le dispositif à plus haute température. Dans une variante de dispositif selon l'invention, non représentée, utilisant des résistances ajustables, le facteur de réduction de la dérive thermique a été d'environ 100. 3048509 11 Les figures. 3 et 4 permettent de déterminer les courbes de variation de la conductivité thermique et de la diffusivité thermique de l'eau avec la température. La courbe obtenue selon l'art antérieur montre l'erreur causée par le déséquilibre du pont. Cette oscillation est un phénomène connu sous le nom de "aliasing". Il dépend des réglages 5 du dispositif et peut être réduit dans certains cas, comme expliqué précédemment, mais en augmentant le temps de réponse d'une manière souvent inacceptable, avec une constante de temps de 10 s ou plus. Dans le cas de mesures relatives, la précision des résultats est de l'ordre de 0,1% sur la conductivité thermique k et de 0,3% sur la diffusivité thermique a. Les 10 instruments à fil chaud connus du commerce avec excitation impulsionnelle ne répondent pas à ces exigences. L'invention permet d'établir un suivi temporel de la température, d'effectuer une mesure de température à l'aide d'une des deux sondes, sans nécessiter l'utilisation d'un dispositif tiers, ou de mesurer des débits. 15 L'invention n'est pas limitée aux exemples qui viennent d'être décrits. L'invention peut notamment être utilisée également dans d'autres applications que pour des liquides, et sous d'autres géométries que celle d'un fil, par exemple pour des tracés de dépôts métalliques sur des matériaux solides. L'invention peut être appliquée à l'anémométrie. 20 L'invention peut être également utilisée dans des applications de microscope à sonde thermique locale, utilisant la correspondance entre la température locale et la conductivité thermique d'une interface. Dans cette technique, l'utilisation de courant alternatif pour l'excitation résulte en une meilleur sensibilité, stabilité et résolution qu'en utilisant un courant continu. L'invention a été appliquée à un alliage nickel-titane à 25 mémoire de forme, microstructure par fraisage et implantation par un faisceau d'ions de gallium. De très faibles variations de signal, 1.75% en amplitude et 0.60° en phase ont pu être observées, à cause de la variation de la conductivité thermique pendant la transition de phase martensite - austénite de l'alliage, comme décrit dans l'article de M. Chirtoc et al. « Temperature-dependent quantitative 3omega scanning thermal microscopy: Local 30 thermal conductivity changes in NiTi microstructures induced by martensite-austenite phase transition », Review of Scientific Instruments no. 79, 093703, 2008. REVENDICATIONS1. Procédé de mesure de propriétés thermo-physiques d'un milieu dans lequel sont placées une sonde de mesure (2) et une sonde de compensation (3), la sonde de mesure (2) étant parcourue par un courant (I.) et la sonde de compensation (3) par un courant (L) de même fréquence et d'amplitude proportionnelle à celle du courant (I.) parcourant la sonde de mesure (2), procédé dans lequel un signal contenant une composante harmonique utile pour la mesure des propriétés thermo-physiques du milieu est généré à partir d'une amplification différentielle des tensions aux bornes de la sonde de compensation (3) et de la sonde de mesure (2), avec un gain (G) choisi pour diminuer l'amplitude de la première composante harmonique au sein du signal résultant. 2. Procédé selon la revendication 1, dans lequel le rapport (N) entre le courant (1,) parcourant la sonde de compensation (3) et celui (I.) parcourant la sonde de mesure (2) est inférieur à 1. 3. Procédé selon la revendication 1 ou 2, dans lequel les sondes de mesure (2) et de compensation (3) sont parcourues par des courants (I., Ic) provenant de générateurs de courant (4, 5), configurés de telle sorte que les amplitudes des courants soient différentes et leurs fréquences (w) identiques. 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel les coefficients de variation avec la température des résistances internes des sondes de mesure (2) et de compensation (3) sont identiques. 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel les résistances internes des sondes de mesure (2) et de compensation (3) sont identiques. 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'amplitude (VA) de la tension aux bornes de la sonde de mesure (2), mesurée en un point A, s'exprime comme suit : V A=[Vh,±V3,,(T)]A=11,,RA(T)+ K(T)I13,,RA(T), où le courant (13.) généré par la non-linéarité de la résistance de la sonde (2) dépend du courant d'excitation (Ii.) : I3p (T) = K(T)I13,,, avec (K(I)) le coefficient de non-linéarité, et = Im le courant parcourant la sonde de mesure (2). 3048509 13 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'amplitude (VB') de la tension aux bornes de la sonde de compensation (3), mesurée en un point B', s'exprime comme suit: VB, =[Vu,±V 3,,(T)]B= NI u,RB (T ) + K (T )(NI 1,, )3 RB (T ) , 5 avec NI 1,) = le courant parcourant la sonde de compensation (3). 8. Procédé selon l'une quelconque des revendications précédentes, dans lequel la tension (VB') aux bornes de la sonde de compensation (3) est amplifiée d'un gain (G), par passage dans un amplificateur opérationnel (7), l'amplitude (Vs) de la tension amplifiée, mesurée en un point B, s'exprimant de préférence comme suit : 10 VB =GV B, =GNI u,RB (T ) + G K (T )(NI 1,, )3 RB (T ) où le gain (G) d'amplification est supérieur à 1. 9. Procédé selon la revendication précédente, dans lequel le gain (G) d'amplification de la tension aux bornes de la sonde de compensation (3) est égal à l'inverse du rapport (N) entre le courant (1,) parcourant la sonde de compensation (3) et 15 celui (I.) parcourant la sonde de mesure (2). 10. Procédé selon la revendication précédente, dans lequel l'amplification différentielle des tensions aux bornes de la sonde de compensation (3) et de la sonde de mesure (2), relevée entre le point A et le point B, après amplification par le gain (G), s'exprime comme suit : V A-V B = (1- N2)K(T)I13,,R(T)= (1- N2)V (T) . 20 11. Procédé selon l'une quelconque des revendications précédentes, dans lequel la température (T) du milieu est calculée en fonction de la tension aux bornes d'une des deux sondes (2, 3), la tension aux bornes de la résistance d'une sonde s'exprimant comme suit en fonction de la température : V(t) I(t)R(t) = oRo{cos(ot)+(112) rT2w[cos(ot-9)+cos(3 ot-9)ll , 25 où r est le coefficient de variation de la résistance avec la température de la sonde. 12. Procédé selon l'une quelconque des revendications précédentes, dans lequel les sondes (2, 3) comportent un fil, notamment de longueur comprise entre 0.1 mm et 10 cm, et notamment de diamètre compris entre 5 lm et 200 ùm. 13. Procédé selon la revendication précédente, dans lequel le fil est 30 métallique, étant notamment réalisé en platine, tantale, nickel ou d'alliages à faible coefficient de dilatation thermique, par exemple un alliage de fer, nickel et cobalt. 3048509 14 14. Procédé selon les deux revendications immédiatement précédentes, dans lequel le fil est recouvert par une fine couche isolante électrique d'oxyde, de polymère ou de résine époxy, adaptée pour des mesures dans des liquides conducteurs électriques. 5 15. Dispositif (1) de mesure de propriétés thermo-physiques d'un milieu, notamment pour la mise en oeuvre du procédé de mesure selon l'une quelconque des revendications précédentes, comportant une sonde de mesure (2) et une sonde de compensation (3) placées dans le milieu à mesurer, la sonde de mesure (2) étant parcourue par un courant (I.) et la sonde de compensation (3) par un courant (1,) de même fréquence 10 et d'amplitude proportionnelle à celle du courant (I.) parcourant la sonde de mesure (2), le dispositif étant configuré de telle sorte qu'un signal contenant une composante harmonique utile pour la mesure des propriétés thermo-physiques du milieu soit généré à partir d'une amplification différentielle des tensions aux bornes de la sonde de compensation (3) et de la sonde de mesure (2), avec un gain (G) choisi pour diminuer 15 l'amplitude de la première composante harmonique au sein du signal résultant. 16. Dispositif selon la revendication précédente, comportant un amplificateur (6) comportant un générateur interne générant une tension de référence (Vie), les amplitudes des courants (I., 1,) parcourant les sondes de mesure (2) et de compensation (3) étant de préférence proportionnelles à ladite tension de référence (Vie).