Process of determining at least one thermal property of a fluid under investigation with a thermal sensor

Abstract

A Process of determining at least one thermal property of a fluid under investigation with a thermal sensor. The thermal sensor has at least a first sensor element that is heated to provide heat to the fluid under investigation. The first or a second sensor that can sense the temperature of the fluid under investigation, wherein the process is characterized by the following steps: a) Providing a calibrated reduced order model which is calibrated with one or more thermal properties of at least a second and a third fluid; b) (Applying an amount of heat to the fluid under investigation by the first sensor element and) measuring the temperature T.sub.sens at the first and/or second sensor element said fluid under investigation; and c) Extracting one or more thermal property of the fluid under investigation by applying the said temperature T.sub.sens to said calibrated reduced order model.

Claims

1. A process of determining at least one thermal property of a fluid under investigation with a thermal sensor by using a calibrated reduced order model, the thermal sensor having at least a first sensor element that is heated to provide heat to the fluid under investigation, wherein the first or a second sensor element can sense the temperature of the fluid under investigation, wherein the process comprises the following steps: a) providing a calibrated reduced order model which is calibrated with one or more thermal properties of at least a second and a third fluid; b) applying an amount of heat to the fluid under investigation using the first sensor element, and measuring the temperature at the first and/or second sensor element of said fluid under investigation; and c) extracting one or more thermal properties of the fluid under investigation by applying the said temperature to said calibrated reduced order model.

2. The process according to claim 1, wherein: the step of providing a calibrated reduced order model includes a generation of a reduced order model from a 3D model of the sensor by applying one model order reduction method.

3. The process according to claim 1, wherein: the step of providing a calibrated reduced order model comprises steps that are applied before the start of the process of extracting at least one thermal property.

4. The process according to claim 1, wherein: the reduced order model is a Proper Orthogonal Decomposition Model.

5. The process according to claim 1, wherein: the thermal sensor operates with a harmonic excitation, and the reduced order model is fitted according to the following equation:
minT.sub.sens,modelT.sub.sens,measured.sup.2, wherein a measured amplitude and a measured phase of T.sub.sens,measured measured by the first or second sensor element, are the entry values for the aforementioned equation, while the output values of the aforementioned equation are thermal properties, wherein T.sub.sens,model is provided by the calibrated reduced order model.

6. The process according to claim 1, wherein: the thermal sensor operates with a 3-omega method or transient method instead of a harmonic excitation.

7. The process according to claim 1, wherein: the thermal sensor determines a thermal conductivity and a volumetric heat capacity of the fluid under investigation.

8. The process according to claim 1, wherein: the step of providing a calibrated reduced order model according to step a) comprises the step of: a1) providing a set of data representing a 3D-model of the thermal sensor.

9. The process according to claim 8, wherein: the step of providing a calibrated reduced order model according to step a) comprises the step of: a2) comprising the adaptation of the set of data of the 3D-model to a reduced order model.

10. A thermal sensor for determining at least one thermal property of a fluid under investigation, having at least: a first sensor element that is heated to provide heat to the fluid under investigation, wherein said first or a second sensor can sense the temperature of the fluid under investigation; and an evaluation unit that is adapted to determine at least one thermal property according to claim 1.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) Hereinafter the invention is explained more in detail with some drawings

(2) FIG. 1 is a direct parametric ROM-approach (not part of the invention);

(3) FIG. 2 is an inverse paramentric ROM-approach;

(4) FIG. 3 is a graphical representation of the thermal conductivity and volumetric heat capacity calculation;

(5) FIG. 4 is a ROM-classes according to the state of the art;

(6) FIG. 5 is a parametric ROM extraction process for receiving a calibrated ROM; and

(7) FIG. 6 is a schematical presentation of a thermal sensor.

DETAILED DISCUSSION IN CONJUNCTION WITH THE DRAWINGS

(8) The current invention relates to a process for the determination of at least one thermal property by thermal sensor. The thermal sensor is for example able to determine the thermal conductivity (k), the volumetric heat capacity (c.sub.p) and/or the thermal diffusivity (a).

(9) It is known in the art that sensor has a model for the calculated of the end physical quantity from the direct measured quantity. An example of an analytical model which is already known in the art for a thermal sensor which consist of a hot wire sensor element can be found in E. Yusibani and P. Woodfield, A procedure for application of the three-omega method to measurement of gas thermal conductivity, Journal of Thermal Science and Technology, vol. 4, no. 1, pp. 146-158, 2009.

(10) In the case of said thermal property sensor the direct measured quantity is the temperature of the temperature sensors, which is derived from the measured resistance or the measured voltage by sensor elements of the thermal property sensor. However, if there is no 1-1 (monotonic) relation between the direct measured quantity and the end measured quantity, respectively the desired properties, normal models often fail.

(11) In the case of a thermal property sensor this non-monotonic behavior could be observed for the end measured quantities, which are the desired thermal properties, such as k and c.sub.p. The direct measured quantities, that means the temperatures, are coupled with both of the thermal properties (k, c.sub.p) and the real physical behavior cannot be described accurately with an analytical modelling approach or with an equivalent electrical circuit thermal model approach.

(12) The current invention proposes the use of a parametric reduced order model for the determination of the desired property values in a thermal property sensor, especially in a transmitter of such a thermal property sensor. This model can provide a far more accurate calculation of the end quantity, respectively the desired properties, such as k and c.sub.p, as the model can take into consideration and compensate effects that cannot be modeled by an analytical model.

(13) Common modern mathematical models of real-life processes used in numerical simulations are often complex and have a large size dimension and need therefore a computing capacity. Compared to these models, a reduced order model has lower the computational complexity of such problems, for example, in simulations of large-scale dynamical systems and control systems. By a reduction of the model's associated state space dimension or degrees of freedom, an approximation to the original model is computed. This reduced order model (ROM) can then be evaluated with slightly lower accuracy but in significantly less time.

(14) An approach for model order reduction is projection-based reduction. The following methods were comprised in this class:

(15) Proper Orthogonal Decomposition (POD-Model); Balanced Truncation; Approximate Balancing; Reduced Basis Method; Matrix Interpolation; Transfer Function Interpolation; Piecewise Tangential Interpolation; Loewner Framework and (Empirical) Cross Gramian. Some more reduced order models are classified in FIG. 4.

(16) In general all these methods can be applied for the modelling of a thermal property sensor and for the determination of the desired properties.

(17) The reduced order model has very low dimensions. In the current case a reduced order model with a matrices with dimensions lower than 3030 can been used. Models of this size can be handled by embedded system micro-processors which are present in the transmitter of said thermal property sensor.

(18) Usually a parametric reduced order model (ROM) is extracted from a respective finite element method (FEM) model and is used for the calculation of the system response, this is the direct or forward problem.

(19) In FIG. 1 the direct parametric reduced order modelling approach is shown. As it is presented in FIG. 1, for a thermal sensor with harmonic excitation (temperature oscillating technique), the two input parameters are the thermal conductivity (k) and the volumetric heat capacity (c.sub.p), solving the reduced order model for specific values of the parameters, the temperature amplitude and phase of the temperature sensor can be calculated. The volumetric heat capacity is the product between the density and the specific heat capacity (c.sub.p).

(20) For the better understating, the parametric reduced order model is a numerical model in form of matrices and/or the matrices have coefficients which consists of the parameters of the system. In other words, is a system of equations of Ordinary Differential Equations (ODEs) or Algebraic Equations (AEs) in matrix form. The system is the combination of the sensor, including sensor geometrics, and the fluid.

(21) Although the parametric reduced order modelling approach has been used for parametric analysis of sensors and their optimization, no use of the parametric reduced order model in the inverse problem has been found in literature for applications of fluid thermal property sensors. Therefore the use of a reduced order model by an evaluation unit of a thermal property sensor for the determination of said properties is so far not known in the art. Especially the reduced order modeling can be used in a sensor which is fully described in WO2015/074833 A1. Reference is made to the construction of the sensors disclosed therein, which can be also applied for a thermal property sensor in the context with the current invention.

(22) The proposed approach of the current invention can preferably but not exclusively be applied for a sensor under harmonic excitation (temperature oscillating technique) and can be summarized in the diagram of FIG. 2.

(23) In case of a thermal sensor with harmonic excitation (temperature oscillating technique), the amplitude and the phase on the temperature sensor are measured and then the model is fitted according to equation (1), cited below, on the measured value of temperature T=|T|e.sup.i with parameters of the thermal fluid properties, especially k and c.sub.p. We have to mention here that the parametric reduced order model has to be calibrated first in known fluid, where k, c.sub.p are known.
minT.sub.sens,modelT.sub.sens,measured.sup.2(1)

(24) FIG. 2 shows schematically an inverse parametric reduced order modelling (ROM) approach for a sensor with harmonic excitation (temperature oscillating technique). A least square fitting algorithm or another approach can be used complementary.

(25) Another option for the extraction of thermal properties (k, c.sub.p) from the reduced order model avoiding a fitting algorithm would be the pre-calculation of an amplitude and a phase for all of the values of thermal conductivity (k) and volumetric heat capacity (c.sub.p). Then during the measurement process the combination of measured temperature amplitude and phase |T|, provides the combination of the thermal properties k and c.sub.p. A graphical presentation of this approach is given in FIG. 3, where the cross point of amplitude |T| and phase iso-curves is used for the extraction of the thermal conductivity and the volumetric heat capacity of the fluid.

(26) Applicability of the method is possible in every thermal property sensor independent of the measurement principle that it uses. For example temperature oscillating technique, transient method, 3, steady state (only thermal conductivity measurement).

(27) Different parametric reduced order modelling techniques can be used for the modelling of a fluid thermal sensor. For example, Reduced Basis method, Krylov-Sub-Space methods, Truncated Balanced Realization and the Proper Orthogonal Decomposition (POD) and other. However the POD-model is the preferred modeling technique.

(28) FIG. 3 shows a graphical representation of thermal conductivity and volumetric heat capacity calculation from the measured amplitude and phase. Amplitude solid lines and phase dashed lines. The cross point of the iso-curves corresponds to the thermal properties k,c.sub.p.

(29) For the further explanation of the use of a parametric reduced order model according to the invention a preferred embodiment is described by the following example:

(30) In the following analysis the thermal property sensor operates using the temperature oscillating technique (harmonic excitation). The thermal property sensor consists of a heater, as a first sensor element, which operates using a harmonic excitation
{dot over (Q)}=|{dot over (Q)}|e.sup.it [W/m.sup.3],
wherein {dot over (Q)} is defined as Energy generation in the volume of the heater
|{dot over (Q)}| is defined as Amplitude of the Energy generation
e.sup.it is defined as exponential form of an harmonic excitation and a temperature sensor, as a second sensor element. The heater and the temperature sensor can also be included in only one first sensor element, which is can operate between two operation modes, a heating mode and a temperature sensing mode.

(31) The temperature on the temperature sensor
T.sub.sens=|T.sub.sens|e.sup.it+,
wherein
T.sub.sens is the temperature by the temperature sensor element of the thermal property sensor in a fluid;
|T.sub.sens| is the amplitude of the aforementioned measured temperature T.sub.sens; and
e.sup.it+ is the exponential form of the harmonic response of the temperature of the sensor including the phase ().

(32) In FIG. 5 a parametric reduced order model extraction process according to an advantageous embodiment is shown. In FIG. 5 the parametric reduced order model extraction process is described. The process begins with a 3D CAD model of the sensor. The parametric reduced order model (ROM) is extracted following one of the reduced order model methods. The extracted reduced order model is parametrised and a simple form can be seen in equation 2 and 3. As can be seen in these equations the fluid thermal properties, which are scalar quantities, are parameter of fluid matrices. It could be a gas or a liquid. The K.sub.r,fluid and Mr.sub.r, fluid are the fluid matrices with dimensions N.sub.rN.sub.r extracted from the ROM. Similarly, the K.sub.solid and M.sub.solid are the matrices which correspond to sensor (solid part) of the domain which are extracted from the ROM. In a more general frame, we can say that these matrices contain the information of the fluid and solid domains respectively. The vector {dot over (Q)}.sub.r(N.sub.r1) corresponds to the energy generation in the domain, while the vector N.sub.r,sens (N.sub.r1) corresponds to the domain of the temperature sensor and is used for the calculation of the averaged temperature. Finally, T.sub.r is the state variable vector (temperature vector). The temperature in the temperature sensor can be calculated from equation 3.
k.sub.fluidK.sub.r,fluidT.sub.r+K.sub.r,solidT.sub.r+i(c.sub.p).sub.fluidM.sub.r,fluidT.sub.r+iM.sub.r,solidT.sub.r={dot over (Q)}.sub.r(2)
T.sub.sens=N.sub.r,sens.sup.HT.sub.r(3) N.sub.r,sens.sup.H is the conjugated transposed vector of N.sub.r,sens.

(33) The second step of the ROM extraction is, according to FIG. 5, the calibration of the parametric reduced order model. The reduce order model can also be calibrated by adding calibration factors in equation 2 and 3. This can eliminate any discrepancies between the model and the real sensor response. During the calibration process the calibration parameters extracted by fitting the model to measured data of a fluid with known properties.

(34) Finally, this calibrated parametric reduced order model can be used inversely for the calculation of the thermal properties (k,c.sub.p) from the direct measured quantities, the temperature of the temperature sensor (T.sub.sens) or sensors in the current case.

(35) According to the current invention it is proposed the use of a parametric reduced order model as a replacement of the analytical models or the equivalent electrical circuit thermal models, which has been used so far in a fluid thermal property sensors with scope the calculation of fluid thermal properties, especially the thermal conductivity and/or the thermal diffusivity as well as the volumetric heat capacity. The reduced order models have low dimensions, which reduces significantly the computation effort. For the extraction of the parametric reduced order model any reduction method, which is applicable in this field can be used, but the parametric POD-model is preferred among the other models.

(36) As mentioned before proposed invention can be used for thermal property sensor independent of the working principle. For example, temperature oscillation technique, 3 method, steady state method, transient method and other. The said invention could also be used in other thermal sensors.

(37) A simplified thermal sensor 1 is shown in FIG. 6. This sensor has a first and a second sensor element 2 and 3, which are arranged in a free space of the sensor. The first sensor element 2 provides heat to a fluid of investigation with unknown properties which surrounds the sensor element, displayed in the plan view of FIG. 6. The second sensor element 3 can measure the temperature of the sensor, including an amplitude and a phase. The measured signals are transferred to an evaluation unit 4, having a microprocessor and a data storage unit. The evaluation unit (4) is adapted to use the reduced order model to determine the thermal properties of the fluid under investigation.