Accurate measurement of unsteady state fluid temperature

In this paper, two accurate methods for determining the transient fluid temperature were presented. Measurements were conducted for boiling water since its temperature is known. At the beginning the thermometers are at the ambient temperature and next they are immediately immersed into saturated water. The measurements were carried out with two thermometers of different construction but with the same housing outer diameter equal to 15 mm. One of them is a K-type industrial thermometer widely available commercially. The temperature indicated by the thermometer was corrected considering the thermometers as the first or second order inertia devices. The new design of a thermometer was proposed and also used to measure the temperature of boiling water. Its characteristic feature is a cylinder-shaped housing with the sheath thermocouple located in its center. The temperature of the fluid was determined based on measurements taken in the axis of the solid cylindrical element (housing) using the inverse space marching method. Measurements of the transient temperature of the air flowing through the wind tunnel using the same thermometers were also carried out. The proposed measurement technique provides more accurate results compared with measurements using industrial thermometers in conjunction with simple temperature correction using the inertial thermometer model of the first or second order. By comparing the results, it was demonstrated that the new thermometer allows obtaining the fluid temperature much faster and with higher accuracy in comparison to the industrial thermometer. Accurate measurements of the fast changing fluid temperature are possible due to the low inertia thermometer and fast space marching method applied for solving the inverse heat conduction problem.

comparison to the industrial thermometer. Accurate measurements of the fast changing fluid temperature are possible due to the low inertia thermometer and fast space marching method applied for solving the inverse heat conduction problem.

Abbreviations c(T) Specific heat of the thermometer material (J/(kgK)) c pa
Specific heat of the air at constant pressure (J/(kgK)) CS Control volume surface d Outer diameter of the thermometer, d = 2r out (m) h Heat transfer coefficient (W/(m 2 K)) k a Air thermal conductivity (W/(mK)) k(T) Thermal conductivity of the thermometer material (W/(mK)) n Normal unit surface vector directed towards the outside Nu a Nusselt number, Nu a = hd/k a Pr a Prandtl number, Pr a = c pa μ a /k ȧ q Heat flux vectoṙ q v Energy generation rate per unit volume (uniform within the body) (W/m 3

) Q i,j
Heat flow rate between the nodes i and j (W) Re a Reynolds number, Re a = wd/ν a r Radius (m) r out Outer radius of the thermometer (m) S Surface area (m 2 ) t Time (s) T(t) Temperature (°C) T 0 Initial temperature (°C) T f (t) Fluid temperature (°C) T i (t) Temperature at the node located inside the cell under consideration (°C) u (t) Response to the unit step Abstract In this paper, two accurate methods for determining the transient fluid temperature were presented. Measurements were conducted for boiling water since its temperature is known. At the beginning the thermometers are at the ambient temperature and next they are immediately immersed into saturated water. The measurements were carried out with two thermometers of different construction but with the same housing outer diameter equal to 15 mm. One of them is a K-type industrial thermometer widely available commercially. The temperature indicated by the thermometer was corrected considering the thermometers as the first or second order inertia devices. The new design of a thermometer was proposed and also used to measure the temperature of boiling water. Its characteristic feature is a cylinder-shaped housing with the sheath thermocouple located in its center. The temperature of the fluid was determined based on measurements taken in the axis of the solid cylindrical element (housing) using the inverse space marching method. Measurements of the transient temperature of the air flowing through the wind tunnel using the same thermometers were also carried out. The proposed measurement technique provides more accurate results compared with measurements using industrial thermometers in conjunction with simple temperature correction using the inertial thermometer model of the first or second order. By comparing the results, it was demonstrated that the new thermometer allows obtaining the fluid temperature much faster and with higher accuracy in

Introduction
The temperature measurement is an important and widely used measurement in industrial processes [1,2]. Various techniques are used to measure the transient temperature of the fluid and the solid surface. One of the most commonly used measurement techniques is based on the use of liquid crystals. Heiland et al. [3] studied experimentally flow and temperature field measurements in a small vertical gap using liquid crystals. Thermochromic liquid crystals particles are dispersed in the fluid flow. The crystal particles change their reflected colour with temperature. Stephan and Fuchs measured local heat temperature fluctuations in wall and fluid in nucleate boiling systems using unencapsulated thermochromic liquid crystals with high spatial and temporal resolution [4]. The instantaneous shape of a growing bubble and the temperature of the liquid crystals beneath the bubble were observed with a high-speed camera. The method of liquid crystals is very useful in the laboratory. Under steady-state conditions when the fluid temperature is constant, there is no damping and time lag and temperature measurement can be made with high accuracy. When the fluid temperature changes over time, the response time of thermometer can be very long, depending on the type of thermometer and the outer diameter of its housing [1,2,[5][6][7][8].
Even if thermocouples or resistance thermometers are thin, dynamic temperature measurement errors can be large [9,10]. Modeling the thermocouple as an inertial object of the first or second order dynamic temperature measurement errors can be greatly reduced [10]. The techniques presented in the paper [10] take into account the dependence of the thermometer time constant on the fluid velocity.
Thick thermowell must be used in an industrial environment when the fluid temperature is measured at high pressure and high flow rate [5][6][7][8]. Thermometer housing must be robust because of the risk of damage to the thermometer due to fatigue. Vibrations of the thermometer pocket are caused by Karman vortices, which appear behind the thermometer [6]. During time variations of the fluid temperature, the thermal inertia of the housing has a significant influence on the measured values. Temperature sensor readings are delayed and attenuated compared to changes in temperature of the fluid. Thermocouples or resistance sensors are located in a housing containing ceramic insulating elements. The increase in the time constant of the thermometer is also due to a large air gap between the sensor and the housing. Accurate measurement of transient fluid temperature is required for proper operation of temperature control systems. Also during the start-up of large steam boilers dynamic errors of temperature measurement should be very small to carry out the optimum heating or cooling of thick-walled pressure elements [11,12].
Accurate measurement of transient fluid temperature is also important in power plants, refineries, and chemical plants. Proper measurement of fluid temperature can prevent overheating of materials from which pipes and pressure elements are made, as well as improve the quality of products in the petrochemical and chemical industries.
This paper presents a new measurement technique that enables significant improvement in the dynamics of the industrial thermometer. By using a thermometer of the new design, as well as a new data reduction method, the thermometer delay, and time constant can be reduced several times.

Techniques for determining transient fluid temperature
In this section, two methods for indirect measurements of transient fluid temperature are presented. First, the thermometer is approximated by the first or second order inertial models and then the inverse heat conduction problem is solved to determine fluid temperature.

Considering thermometer to be a first and second order inertia device
Usually, the thermometer is modelled as an element with lumped thermal capacity. In this way, it is assumed that the temperature of the thermometer is only the function of time, and temperature differences occurring within the thermometer are neglected. The temperature changes of the thermometer in time T(t) have been described by an ordinary first order differential equation (first order thermometer model) [1,2] (1) For thermometers with a complex structure used for measuring the temperature of the fluid under high pressure, the accuracy of the first order model (1) is inadequate. To take into account that the temperature sensor readings are delayed when a transient fluid temperature is measured with heavy thermometers, the second order model will be used. In this case, thermometer can be modelled using the ordinary differential equation of the second order which can be written in the form [8,9]: Equations (1) and (2) were solved for the step increase of the fluid temperature from T 0 = 0 °C to the constant value T f . The initial conditions are: The initial problem (2)-(4) was solved using the Laplace transformation and the temperature as a function of time is expressed by the following formula In the same way, the initial problem expressed by Eqs. (1) and (3) was solved. For the first order model the thermometer response for a unit step fluid temperature change is determined by the simple expression: The time constant τ in Eq. (6) and time constants τ 1 and τ 2 in Eq. (5) will be estimated from experimental data.
The fluid temperature can be determined based on the measured history of the thermometer temperature T(t) and known time constant τ using Eq. (1) or time constants τ 1 and τ 2 when using Eq. (2).

Inverse marching method
An inverse space marching method is used to determine the temperature distribution inside the thermometer and the fluid temperature. The transient heat conduction equation for a cylindrical wall is [8] Based on the temperature T 1 (t) measured at the cylinder axis r = 0 the temperature distribution inside the cylinder including the outer surface will be determined. The unknown fluid temperature T f (t) appears in the convective boundary condition at the outer surface of the cylinder The solution of the inverse problem (7,8) is obtained by using the Finite Volume Method (FVM) [13].
The Finite Volume Method (FVM), also called control volume method, is a universal and effective method for solving heat conduction problems. Assuming that the temperature field is two-dimensional and thermal properties c, ρ, λ and power density of internal heat sources q v are temperature dependent, then heat conduction equation is as follows Once Eq. (10) is integrated over the control volume (CV), the following equation for a single cell (CV) is obtained: Applying the Green theorem to the first term on the right-hand-side of Eq. (11) gives where CS is the control volume surface, while n a normal unit surface vector directed to the outside of the control volume. If ΔV i designates the volume of a control cell, the individual terms in Eq. (12) can be approximated in the following way: where Q i,j is the heat flow that is transferred from the node j located inside the neighbouring cell and the node i.
The node i is located inside the analyzed region. Given that heat transfer takes place with the n control volumes adjacent to the analyzed control domain, and taking into account Eqs. (14) and (15), Eq. (13) can be written as: The thermometer housing is modelled as a solid cylindrical without an inner hole, divided into four control volumes (Fig. 1). The temperature sensor is placed at the axis of the cylinder. Temperatures in the nodes 2, 3 and 4 are calculated based only on the temperature measured at node 1 solving the inverse heat conduction problem.
Applying the FVM, the energy balance equations can be written for all nodes which are located at the centres of the control volumes. The energy balance equations are written for the nodes: 1, 2, and 3 of the control volumes. Taking into account that the temperature T 1 (t) at the node 1 is measured, the temperature at the node 2 is determined form the energy balance equation for the node 1. In a similar way, the temperature at the node 3 is obtained from the energy balance equation for the node 2, and subsequently the temperature at the node 4 from the energy balance equation for the node 3. The thickness of the control volume adjacent to the outer surface of the cylinder is equal to Δr/2. Such arrangement of nodes ensures equal distance Δr between all nodes. Furthermore, the temperature of the thermometer surface can be determined with higher accuracy, if the node 4 is located on the thermometer surface. The location of the node 4 on the cylinder boundary facilitates the formulation of the boundary condition (9) in a discrete form. The temperatures at the nodes 2, 3, and 4 are as follows: where the symbol Δr = r out /3. The energy balance equation for the node 4 is The surface temperature T f is determined sequentially. Calculating the time derivatives dT 1 /dt and using Eq. (17), the temperature at the node 2 is obtained at first. Next, the temperatures at the nodes T 3 and T 4 are calculated in sequence using Eqs. (18) and (19).
Values of thermal conductivity k(T i ) for temperature T i at the i node (where i = 2, 3, 4) are unknown. It is necessary to execute several iterations to obtain temperature T i . It was assumed that the starting value of the thermal conductivity for the first iteration is The iteration continues until where the small value ε is equal to 0.000 01 K.
The value of ε must be small enough so that the temperature dependent properties of the thermometer are taken into consideration. Iterative calculations were carried out at different values of tolerance ε. When the values of ε equal (17) Having determined temperatures T 2 , T 3 , and T 4 , the fluid temperature is calculated using Eq. (21).

Computational validation of the inverse method
The method developed in the paper was used to determine transient fluid temperature. At first, direct heat conduction problem was solved for two cases: stepwise and linear temperature increase. The temperature T 1 (t) at the centre of the solid cylinder r = 0 for sudden temperature increase was calculated using an analytical method [13]. In the case of a linear rise of fluid temperature, the temperature at the centre of the cylinder was determined by Duhamel integral using the solution of the direct heat conduction problem for the unit step change in temperature of the fluid.
The calculations were carried out for a thermometer with an outer diameter of 7 mm made of 1.4541 steel (austenitic steel), with the following physical properties: ρ = 7.9·10 3 kg/m 3 , k = 18 W/(m K), c = 500 J/(kg K). The heat transfer coefficient at the outer surface is h = 2000 W/(m 2 K). The cylinder temperature was calculated (direct problem) for the stepwise increase in fluid temperature from 20 to 100 °C and linear rise from 0 to 170 °C with temperature change rate equal: v T = 0.33333 K/s. Finally, the inverse problem was solved using T 1 (t) as the input data. The time step size for the inverse analysis was determined from the following condition [13]:

Experimental study
Experimental studies included measurement of water and transient air temperature. Two thermometers were used, both with a housing diameter of 15 mm, but of different construction. Measurements were carried out using the developed thermometer and an industrial thermometer. Time constants of the industrial thermometer were determined while the water or air temperature was measured. In the first experiment, the test thermometer was suddenly immersed in boiling water to determine its time constants. Next, the time constant of the thermometer was determined in air flowing at different velocities.
The first gauge is an industrial thermometer widely available commercially. The housing of thermometer having an inner diameter of 9 mm is made of 1H18N9T austenitic steel. The chemical composition of the steel is as follows: 0.1 % C, 18 % Cr, 9 % Ni, 2 % Mn, 0.8 % Si, 0.045 % P, and 0.03 % S. The temperature sensor is a K type thermocouple (NiCr-NiAl). Two wires with a diameter of 0.5 mm each form a thermocouple. Hot junction of  the thermocouple is not grounded. Six ceramic cylinders of 6 mm diameter and 50 mm length are inserted into the thermometer pocket. There are four holes with a diameter of 1 mm in the cross-section of the each ceramic cylinder, and only through two holes thermocouple wires are conducted. There is a gap of 1.5 mm between ceramic cylinders and the inner surface of the thermometer thermowell.
The second thermometer of the new design has a cylindrical casing of external diameter 15 mm, inside which a hole with a diameter of 1.1 mm is drilled. The outer diameter of the sheath thermocouple inserted into the inner hole is 1 mm. The thermometer presented in this paper is characterized by a very small air gap between the thermocouple sheath and the surface of the hole drilled axially throughout the solid cylinder. The thermocouple of K type with the grounded junction is used to measure temperature. The thermowell is made of 1H18N9T stainless steel.
The temperature sensors used in the experiment have been calibrated in the temperature range from 25 to 120 °C in a calibration oven to obtain their characteristics.

Water temperature measurement
In the experiments, the saturated water temperature was measured using the developed and industrial thermometer. Thermometers at room temperature were suddenly immersed in boiling water. In addition, the water temperature was measured by a bare sheath thermocouple having a diameter of 1.5 mm. The K type sheath thermocouple was used to measure the water temperature continuously during the experiment. The time constant of the thermocouple is low due to its small diameter (Fig. 4).
The temperature indicated by the thermometer sensors are shown in Fig. 5.  The measurement of transient temperature by the industrial thermometer has a very big dynamic error. To reduce this error, mathematical model of the first and second order was used to calculate water temperature. For this purpose, the time constant τ was determined for the model of the first order and for the second order model two time constants τ 1 and τ 2 were determined. Based on water temperature measurement the time constant τ was found to be: τ = 67.156 ± 0.484 s. The time constants for the thermometer considered as the second order model are: τ 1 = 5.360 ± 0.334 s and τ 2 = 61.399 ± 0.445 s.
The 95 % uncertainties in the calculated time constants were estimated using the variance propagation rule formulated by Gauss [14]. The mean square error s N of the measured temperature was assumed equal for all measuring time points and calculated using the formula: where m is the number of time constants and N denotes the number of time points in a measurement series. The symbol S min designates the minimum sum of squares of differences between measured and calculated temperatures.
The fluid temperature T f (t) was determined using Eqs. (1) and (2). To eliminate the random errors inherent to actual temperature measurements and to determine the first and second order derivatives of measured temperature with respect to time, a 9 point digital filter was used [15]. The results of the calculations are shown in Fig. 5a.
The proposed thermometer had much better dynamic properties when water temperature was measured in comparison with a conventional thermometer. The temperature measurement made by the thermometer with the housing proposed is burdened with much smaller dynamic error compared to the measurement made by the industrial thermometer (Fig. 5b). Even for quicker response of the thermometer and more accurate measurements, an inverse marching method was applied. Equations (17)- (19) and (21) were used to obtain the water temperature. Calculations were conducted for constant physical properties of the thermometer housing: specific heat c = 500 J/ (kg K), density ρ = 7 900 kg/m 3 and thermal conductivity k = 18 W/(m K). Heat transfer coefficient on the outer surface of the thermometer was h = 2 381 W/(m 2 K). The results of the calculations are shown in Fig. 5b.
The fluid temperatures obtained by using a conventional thermometer in conjunction with the first and second order thermometer model and by using developed thermometer in combination with the inverse method are compared with the reference temperature measured by a thermocouple with the outer diameter of 1.5 mm (Figs. 5 and 6).
The analysis of the results depicted in Figs. 5 and 6 show that the technique proposed for measuring the fluid temperature by the new thermometer is much more accurate. The conventional thermometer has a large inertia and even temperature correction using the model of the first or second order does not give satisfactory results.

Air temperature measurement
At first, the time constants of the industrial thermometer τ for various air velocities w were determined in an open benchtop wind tunnel (Fig. 7). The WT4401-S benchtop wind tunnel is designed to give uniform flow rate over a 100 mm × 100 mm cross section (WT4401-S & WT4401-D Benchtop Wind Tunnels). Based on the thermometer response for a step change in air temperature, the time constants for different air velocities were determined using the first order thermometer model. The time constant of a conventional thermometer that was estimated experimentally is shown in Fig. 8 as a function of air velocity. The experimental data were approximated by the least squares method. The experimental data very well approximates the following function ( Fig. 8): where τ is expressed in s, and w in m/s. (1), 2-temperature difference between the reference temperature and fluid temperature calculated by Eq. (2), 3-temperature difference between the reference temperature and fluid temperature calculated using the inverse marching method Next, the transient temperature in the wind tunnel ( Fig. 9) was measured. Two thermometers described in Subsection 4 were used, both with a housing diameter of 15 mm.
First, the measurements were made using the industrial thermometer and then using a newly designed thermometer. The air in the tunnel was heated by a heat exchanger, and its velocity was also varied. The thermometers measured the temperature behind the heat exchanger. The air velocity was measured by the vane anemometer. Both the velocity and temperature data were collected using the Ahlborn ALMEMO 5990-0 data acquisition system. Also, the calibrated K type sheath thermocouple with an outer diameter of 0.5 mm was used to measure continuously the air temperature behind a heat exchanger (Fig. 9).
By the sudden drop in air velocity of 2.5-1.1 m/s, the air temperature after the heat exchanger increased from 48 to 55 °C (Fig. 10). Indications of a thermocouple located inside the industrial thermometer are significantly delayed and damped with respect to changes in air temperature. Similar behavior is observed for the thermocouple inside of the proposed thermometer (Fig. 11). As a result of a decrease in air velocity of 2.5-1.0 m/s the air temperature after heat exchanger rose from 47 to 54 °C. The air-site heat transfer coefficient on the outer surface of the housing of the thermometer is small. With such a small heat transfer coefficient heating of the commercial thermometer and the developed thermometer is very slow. From the comparison of the results of measurements shown in the Figs. 10 and 11 can be seen that there are no major differences in the time of heating of the thermocouple inside the thermometer.  [16]: 1-test chamber with opening for the thermometer; 2-differential pressure measurement; 3-data acquisition system, 4-opening for the thermometer insertion, 5-restrictive plate for adjusting velocity measurement range  Based on measurements, the transient temperature of the air was determined assuming that the industrial thermometer is an inertial object of the first order. Calculating the time constant from Eq. (25) and using Eq. (1) the air temperature was estimated. Because the air velocity varied in time, the time constant of the thermometer is also a function of time. The computations were performed for the sampling time Δt = 6 s. The results of the calculations are shown in Fig. 10.
Then, the transient air temperature was measured using a new thermometer. To determine the temperature of the air more accurate data reducing method was used. The physical properties of the thermowell material are the same as in the case of measuring the water temperature.
Different correlations may be used to calculate the heat transfer coefficient on the outer surface of the thermometer [18,19]. In this paper, the heat transfer coefficient on the outer surface of the thermometer housing in the Eq. (21) was calculated from the relationship: (26) Nu a = 1.3Re 0.5 a Pr 0.31 a , 10 2 ≤ Re a ≤ 10 4 where: Nu a = hd k a , Re a = wd ν a and Pr a = c pa µ a k a . The inverse marching method given by Eq. (17)- (19) and (21) was applied for determining the air temperature. The following physical properties of the air were used in the evaluation of the heat transfer coefficient: Prandtl number Pr a = 0.698, thermal conductivity k a = 0.0283 W/(m K), and kinematic viscosity ν a = 17.95·10 −6 m 2 /s. The calculations were made with the time step of 6 s.
An inspection of results depicted in Fig. 11 shows that the proposed measuring technique gives the air temperature that is very close to the temperature measured using a thin thermocouple.

Conclusions
Both methods presented in the paper allow for the accurate measurements of transient fluid temperature but are suitable for different applications. The first method is suitable for measuring gas temperature, e.g. the temperature of the flue gas or air when the temperature is measured using thin bare thermocouples. The second technique based on the solution of the inverse heat conduction problem is more appropriate for measuring the fluid temperature at high pressure and temperature, e.g. water or steam in boilers, pipelines or turbines. By comparing the measurement results, it was demonstrated that the proposed measuring technique allows obtaining the real fluid temperature much faster compared to industrial thermometers. Accurate measurement of fast fluid temperature changes is possible due to the designed low inertia thermometer and fast space marching method for solving inverse heat conduction problem. The fluid temperature is calculated on the basis of temperature indicated by the thin thermocouple located at the axis of the cylindrical thermometer.
The results presented in the paper are very important from a practical point of view. Using the inverse method for transient fluid temperature measurement proposed in the paper, the start-up of steam boilers and steam or gas turbines can be carried out in the right way. Using the method developed for measuring the temperature of the fluid can greatly improve the quality of temperature control of superheated steam, flue gases, etc.
The measuring technique is simple, yet accurate. Assuming that the heat transfer coefficient can be calculated from available in literature heat transfer correlations, the fluid temperature can be calculated in on-line mode. The advantage of the inverse method is the independence of the results from the initial temperature distribution in the thermometer. This is very important because if the previous measurement data are lost, temperature measurement can be carried further after picking the first four points. The experimental test, which has been carried out, illustrates the good accuracy of the method and its computational efficiency.