Open access peer-reviewed chapter

Measurement of Transient Fluid Temperature in the Heat Exchangers

Written By

Magdalena Jaremkiewicz

Submitted: 13 April 2016 Reviewed: 09 September 2016 Published: 27 April 2017

DOI: 10.5772/65686

From the Edited Volume

Heat Exchangers - Design, Experiment and Simulation

Edited by S M Sohel Murshed and Manuel Matos Lopes

Chapter metrics overview

2,238 Chapter Downloads

View Full Metrics

Abstract

In the chapter, a method for measuring the transient temperature of the flowing fluid based on time temperature changes of the thermometer is described. In the presented method, the thermometer is considered as an inertial system of first and second order. To reduce the influence of random errors in the temperature measurement, the local polynomial approximation based on nine points is used. As a result, the first and second derivatives of a temperature, which indicate how the temperature of the thermometer varies over time, are determined very accurately. Next, the time constant is defined as a function of fluid velocity for sheathed thermocouples with different diameters. The applicability of the presented method is demonstrated on real data in the experiment. The air temperature is estimated from measurements carried out by the three thermocouples having different outer diameters when the air velocity varied in time. A comparison of the computed temperatures of air gives confidence to the accuracy of the presented method. The method presented in this chapter for measuring the transient temperature of the fluid can be used for the online monitoring of fluid temperature change with time.

Keywords

  • temperature measurement
  • transient conditions
  • first-order model
  • second-order model
  • time constant

1. Introduction

Under steady-state conditions, when the fluid temperature is constant, temperature measurement can be performed with the high degree of accuracy owing to the absence of damping and time lag. However, when fluid temperature varies rapidly, for example, during start-up, there is a significant difference between the actual and the measured temperature of the fluid. These differences occur because it takes time for heat to transfer through the heavy housing to the temperature sensor.

Most of the scientific publications concerning the measurement of temperature mainly discuss the problem of temperature measurement at steady state [19]. Only the step response of thermometers is studied to estimate the dynamic error of the temperature measurement. Few studies refer to the measurement of the transient fluid temperature, despite the high practical significance of the problem [1013].

An example is the measurement of transient temperature steam or flue gases in power plants, which is very difficult. Measured temperature differs significantly from the real temperature of the fluid, which is caused by massive thermowells of the thermometers and their low heat transfer coefficients. Some thermometers may have a time constant of 3 min or more, which makes the implementation of a single temperature measurement requiring about 15 min [13].On the other side, some designs of thermometers need more than one time constant to describe the unsteady response of temperature sensor inserted into the thermowell. Measuring the fluid temperature in a controlled process may require the knowledge of two or three time constants, which describe the transient response of the thermometer [14].

The problem of dynamic errors in temperature measurements becomes particularly important in superheated steam temperature control systems, which use injection coolers (spray attemperators). Due to the large inertia of the thermometer, the measurements of transient fluid temperature can be inaccurate, which causes the automatic control system of the superheated steam temperature not to work properly.

A similar problem occurs in measuring the exhaust gas temperature as the time constants of the thermometer and the time delay are large.

Advertisement

2. Mathematical models of the thermometers

Generally, thermometers are considered as elements with lumped thermal capacity. It is assumed that the temperature of the thermometer is only a function of time and temperature differences inside it are neglected.

Based on these assumptions, the mathematical model of the thermometer is the differential equation describing the inertial system of the first order [14]:

τdTtdt+Tt=TftE1

where τ—time constant of the thermometer in the first-order model in s, T(t)—measured temperature in °C and Tf(t)—fluid temperature in °C.

The time constant is expressed by formula:

τ=mtcthtAtE2

where mt—thermocouple mass in kg, ct—average-specific heat of the thermocouple in J/(kg⋅K), ht—heat transfer coefficient on the outer surface of the thermocouple in W/(m2⋅K) and At—outer surface area of the thermocouple in m2.

The ordinary differential Eq. (1) was solved for the initial condition:

T0=T0=0E3

where T0 signifies initial thermometer temperature in °C.

The initial problem, (Eqs. (1) and (3)), was solved using the Laplace transformation. The operator transmittance G(s) then assumes the following form:

Gs=TsTfs=1τs+1E4

where T(s)—Laplace transform of the thermometer temperature, Tf(s)—Laplace transform of the fluid temperature and s—complex variable.

For the step increase of the fluid temperature from T0 = 0°C to the constant value Ts, the Laplace transform of the fluid temperature assumes the form Tf(s) = Ts/s and the transmittance formula can be simplified to

TsTs=1sτs+1E5

After writing Eq. (5) in the form:

TsTs=1s1s+1τE6

it is easy to find the inverse Laplace transformation and determine the thermometer temperature as a function of time:

ut=TtT0TsT0=1exptτE7

where u(t) denotes unit-step response of the thermometer.

For structurally complex thermometers that measure the temperature of fluid under high pressure, the accuracy of the first-order model (1) is inadequate. The cross section of the temperature sensor placed in the massive housing is shown in Figure 1. This example will be analysed to show that the thermometers can be modelled as second-order inertial systems [14].

Figure 1.

Part of the cross section of the temperature sensor together with the housing [14].

Air gap can occur between the outer thermowell and the temperature sensor (Figure 1). The discussion assumes that the thermal capacity is neglected because of its small value.

Introducing the overall heat transfer coefficient uin referenced to the inner surface of the housing

1uin=1hin+1+Din/dDind4kw+Dind1hTE8

and accounting for the radiation heat transfer from the housing to the inner sensor, the heat balance equation for the sensor located within the housing assumes the form:

AρcdTdt=PinuinTwT+CTw4T4E9

where the symbol C denotes:

C=πdσ1εT+dDin1εw1E10

In Eqs. (8) and (9), the symbols hin and hT represent heat transfer coefficient on the inner surface of the housing and the outer surface of the thermocouple in W/(m2K), respectively, Din inner diameter of the housing in m, d outer diameter of the thermocouple in m, ko housing thermal conductivity in W/(m⋅K), A surface area of the thermocouple cross section in m2, ρ average density of the thermocouple in kg/m3, c average-specific heat of the thermocouple in J/(kg⋅K), Pin perimeter of the internal surface of the housing in m, Tw housing temperature in °C, σ = 5.67 × 10−8 W/(m2⋅K4) Stefan-Boltzmann constant and εw and εT emissivity of the inner housing and the thermocouple surface, respectively.

The convection and conduction of heat transfer between the fluid and the thermometer housing are characterised by the overall heat transfer coefficient uout referenced to the outer housing surface:

1uout=1hout+1+Dout/Din2δwλwE11

where hout is the heat transfer coefficient on the outer surface of the housing in W/(m2⋅K), Dout outer diameter of the housing in m, δw housing thickness in m and kw housing thermal conductivity in W/(m⋅K).

The formulas (8) and (11) for the overall heat transfer coefficients were derived using the basic principles of heat transfer [2, 4, 8]. The heat transfer equation for the housing (thermowell) can be written in the following form:

AwρwcwdTwdt=PoutuoutTfTwPinuinTwTCTw4T4E12

where Aw is the surface area of the housing cross section in m2, ρw average density of the housing in kg/m3, cw average-specific heat of the housing, J/(kg⋅K) and Pout perimeter of the external surface of the housing in m.

Further analysis of the heat exchange between the housing of the thermometer and the temperature sensor omits the heat transfer by radiation [14]. This is possible when the gap between the thermowell and the temperature sensor is filled with a non-transparent substance or if one of the two emissivities εw and εT is near to zero.

Transforming Eq. (9), we get:

Tw=AρcPinuindTdt+TE13

Substituting Eq. (13) into Eq. (12) yields:

AwρwcwAρcPinuinPoutuoutd2Tdt2+AwρwcwPoutuout1+PoutuoutAρcPinuinAwρwcw+AρcAwρwcwdTdt+T=TfE14

Introducing the following coefficients:

a2=AwρwcwAρcPinuinPoutuout,a1=AwρwcwPoutuout1+PoutuoutAρcPinuinAwρwcw+AρcAwρwcwE15

the ordinary differential equation of the second order (14) can be written in the form:

a2d2Tdt2+a1dTdt+T=TfE16

The initial conditions are:

T0=T0=0,dTtdtt=0=vT=0E17

Equations (16) and (17) are solved using the Laplace transformation. The operator transmittance G(s) is as follows:

Gs=TsTfs=1τ1s+1τ2s+1E18

The time constants τ1 and τ2 in Eq. (18) are:

τ1,2=2a2a1±a124a2E19

The differential Eq. (16) can be written in the following form:

τ1τ2d2Tdt2+τ1+τ2dTdt+T=TfE20

Equation (20) is solved for a step change in fluid temperature from T0 = 0°C to the constant value Tf. Laplace transform of the constant fluid temperature Tf is Tf(s) = Ts/s, and the operator transmittance assumes the following form:

TsTs=1sτ1s+1τ2s+1E21

Eq. (21) can be written in another form:

TsTs=1s+τ1τ2τ11s+1τ1τ2τ2τ11s+1τ2E22

Making inverse Laplace transformation of Eq. (22), the thermometer temperature as a function of time is obtained:

ut=TtT0TfT0=1+τ1τ2τ1exptτ1τ2τ2τ1exptτ2E23

If we assume in Eq. (23) with τ2 = 0, then we obtain Eq. (7) with τ = τ1.

A time delay in the time response of the first order inertial system does not appear. Measuring the temperature of the fluid at high pressure requires the use of thermometers with the massive housings. In such cases, there is a time delay between the temperature indicated by the sensor and the changing temperature of the fluid. The inertial system of second order is suitable to describe the response behaviour with a time delay [13].

Advertisement

3. Identification of time constants

The time constant τ in Eq. (7) or time constants τ1 and τ2 in Eq. (23) can be determined on the basis of experimental data. For this purpose, the method of least squares is proposed to use [13, 15]. Finding the minimum of the function S

S=i=1Numtiuti2=minE24

allows to determine values of the time constants. In Eq. (24), u(ti) denotes the approximating function given by Eq. (7), and N is the number of conducted measurements (ti, um(ti)). The sum of the squares of the differences of the measured values um(ti) and the fitted values u(ti) is minimised. When the time constant τ or time constants τ1 and τ2 are determined, their values can be substituted into Eq. (24) to calculate Smin.

The uncertainties of the calculated time constant τ or time constants τ1 and τ2 are calculated using the mean square error [13, 1618]:

SN=SminNmE25

where m denotes the amount of time constants (i.e. m = 1 for Eq. (7) and m = 2 for Eq. (7) and m = 2 for Eq. (23)). Based on the determined mean square error SN, which is an approximation of the standard deviation, the uncertainties in the determined time constants can be calculated using the formulas given in the TableCurve 2D software [18].

Advertisement

4. Determination of transient fluid temperature

The fluid temperature can be determined on the basis of measured histories of the thermometer temperature T(t) and known time constant τ, using Eq. (1), or time constants τ1 and τ2, using Eq. (20). Temperature of the thermometer T(t) and its first- and second-order time derivatives can be smoothed using the formulas [3, 13]:

Tt=169363ft4Δt+42ft3Δt++117ft2Δt+162ftΔt+177ft+162ft+Δt++117ft+2Δt+42ft+3Δt63ft+4ΔtE26
Tt=dTtdt=11188Δt86ft4Δt142ft3Δt193ft2Δt126ftΔt+126ft+Δt++193ft+2Δt+142ft+3Δt86ft+4ΔtE27
Tt=d2Ttdt2=1462Δt228ft4Δt+7ft3Δt8ft2Δt17ftΔt20ft17ft+Δt8ft+2Δt+7ft+3Δt+28ft+4ΔtE28

Symbol f(t) in Eqs. (26)–(28) denotes the temperature measured by the thermometer, and Δt is a time step. Application of nine-point digital filter allows eliminating the influence of random errors of the measured temperature T(t) on the determined temperature of the fluid Tf(t). If the measured temperature is not disturbed by significant errors, the first and second derivatives of the temperature can be estimated using the central difference method [15]:

Tt=ft+ΔtftΔt2ΔtE29
Tt=ft+Δt2ft+ftΔtΔt2E30

Advertisement

5. Examples of application of the method for a step change in temperature of the fluid

The first example demonstrates the use of the described method of transient fluid temperature measurement in the case of thin sheathed thermocouple application. During the experiment, the thermocouple with outer diameter 1.5 mm at the ambient temperature was suddenly immersed into hot water at saturation temperature. The results are presented in Figure 2. The time step Δt of the measured temperature was 1.162 s. The measured history of the temperature was approximated by Eq. (7) in order to determine the time constant τ of the thermocouple. In the calculations, TableCurve 2D code [18] was used. The calculations of the time constant and the uncertainty at the confidence level of 95% gave the following results: τ = 1.54 ± 0.09 s.

Figure 2.

Fluid and thermometer temperature changes determined from the first-order equation (1) for the sheathed thermocouple with outer diameter 1.5 mm.

Next, the transient fluid temperature Tf was determined using Eq. (1) together with Eqs. (26) and (27). The first-order time derivative dT/dt in Eq. (1) was also calculated using the central difference quotient of Eq. (29). The results obtained show that the use of a first-order model for thin thermometers is sufficient (Figure 2). The results also indicate that the central difference approximation of the time derivative in Eq. (1) leads to less accurate results, since it is more sensitive to random errors in the experimental data.

Another example shows the application of the method for measuring transient fluid temperature with an industrial thermometer with massive thermowell and its complex construction (Figure 3). As in the previous example, the thermometer at the ambient temperature was suddenly immersed into hot water at a temperature of about 100°C. To compare the two methods of determining the unknown fluid temperature (using the first-and second-order model) for this thermometer, Eqs. (7) and (23) were used as the functions approximating the transient response of the thermometer. The following values with the 95% confidence uncertainty were obtained: τ = 14.07 ± 0.39 s, τ1 = 3.0 ± 0.165 s and τ2 = 10.90 ± 0.2 s.

Next, the fluid temperature changes were determined from Eq. (1) for the first-order model and from Eq. (20) for the second-order model (Figure 4).

The analysis of the results presented in Figure 4 indicates that the second-order model delivers more accurate results than the first-order model.

Figure 3.

Industrial thermometer and its dimensions: D = 18 mm, l = 65 mm and L = 140 mm.

Figure 4.

Fluid and industrial thermometer temperature changes determined from the first-order equation (1) and the second-order equation (20).

Advertisement

6. Experimental determination of time constant as a function of fluid velocity

The method presented in this chapter can also be used to determine the temperature of the flowing fluid. In this case, the time constant of the thermometer τ should be determined as a function of fluid velocity w. After specifying function τ(w), it should be substituted into Eq. (1).

Research has shown that for thin sheathed thermocouples with a diameter ranging from 0.5 to 6 mm, the equation of time constant as a function of fluid (air) velocity is as follows [15]:

τ=1a+bwE31

where τ is the time constant of the thermometer in s, a constant in 1/s, b constant in (m·s)−1/2 and w fluid velocity in m/s.

As an example, the constants a and b were determined for K-type sheathed thermocouples with grounded hot junction with the outer diameter of 0.5, 1.0, 1.5 and 3.0 mm.

The thermocouple time constant τ for various air velocities w was determined in an open benchtop wind tunnel (Figure 5). The WT4401-S benchtop wind tunnel is designed to give uniform flow rate over a 100 mm × 100 mm cross section [19].

Figure 5.

Wind tunnel used for determining thermocouple time constant—overall view [19]: (1) test chamber with an opening for the thermometer, (2) differential pressure measurement, (3) data acquisition system and (4) opening for the thermometer insertion.

The calculated time constants τ of the sheathed thermocouples with the outer diameter of 0.5, 1.0, 1.5 and 3.0 mm for different velocities of the air are shown in Figure 6. The experimental data collected for all thermocouples, as presented in Figure 6, were approximated by the least squares method. The best estimates for the constants a and b, with uncertainty at the 95% confidence level, for the tested thermocouples are [15]:

  • dt = 0.5 mm

    a = 0.004337 ± 0.000622 1/s and b = 0.022239 ± 0.001103 (m·s)−1/2

  • dt = 1.0 mm

    a = 0.020974 ± 0.006372 1/s and b = 0.103870 ± 0.011240 (m·s)−1/2

  • dt = 1.5 mm

    a = 0.040425 ± 0.003301 1/s and b = 0.056850 ± 0.004479 (m·s)−1/2

  • dt = 3.0 mm

    a = 0.128220 ± 0.035716 1/s and b = 0.220641 ± 0.051122 (m·s)−1/2

The time constant of the thermocouple τ strongly depends on the heat transfer coefficient ht on the outer thermometer surface, which results from Eq. (2). The heat transfer coefficient is a function of Nusselt number, and this, in turn, is a function of fluid (air) velocity [20].

Figure 6.

Time constants τ of sheathed thermocouple with outer diameters of 0.5, 1.0, 1.5 and 3.0 mm as a function of air velocity w with 95% confidence intervals.

When the velocity of the fluid and its temperature varies over time, time constant as a function of the velocity formulated by Eq. (31) can be used in Eq. (1) to determine the transient fluid temperature based on measurements made using a sheathed thermocouple.

To show how the described method can improve the temperature readings, the experimental measurements were presented [21]. The temperature of the flowing air in an open wind tunnel (Figure 7) was measured by K-type sheathed thermocouples with outer diameters of 0.5, 1.0 and 1.5 mm. During measurements, the temperature and velocity of air flowing through the tunnel were changed. The thermocouples were placed in the tunnel behind the heat exchanger and very close to each other (i.e. they measured the same temperature, but did not influence each other). The air velocity was measured by the vane anemometer FV A915 S220. Both the temperature and velocity data were collected using the Ahlborn ALMEMO 5990-0 data acquisition system.

Figure 7.

Diagram of an open wind tunnel [22]: (A) heat exchanger, (B) fan, (C) chamber, (D) air channel, (E) water outlet pipe and (F) hot water feeding pipe.

The comparison of the computed temperatures with the measured temperatures, when the time constants of the thermocouples are known, shows that the above method provides decent results (Figures 8 and 9).

The time histories of temperature obtained from calculations are very similar, especially for thermocouples with the outer diameter of 0.5 and 1.0 mm. In the small degree, the temperature calculated using the measurements with the thermocouple with the outer diameter of 1.5 mm deviates from them. This difference is due to the large time constant of the sheathed thermocouple with the outer diameter of 1.5 mm.

Figure 8.

Temperature of the air measured by the thermocouples with outer diameters of 0.5, 1.0 and 1.5 mm and temperature calculated by Eq. (1) when the velocity of the air was constant.

Figure 9.

Temperature of the air measured by the thermocouples with outer diameters of 0.5, 1.0 and 1.5 mm and temperature calculated by Eq. (1) when the velocity of the air was changed.

Advertisement

7. Conclusions

Both methods of transient fluid temperature measurement presented in the chapter can be used online.

The first method, in which a mathematical model of the thermometer is first-order differential equation, is the most suitable for thermometers with very small time constants. In such cases, the delay of thermometer indication is small compared to changes in fluid temperature. In turn, the industrial thermometers, which are designed to measure the temperature of the fluid of high pressure, are characterised by a considerable delay of indications in reference to the actual changes of fluid temperature. For such thermometers, the second-order thermometer model, allowing for modelling the signal delay, is more adequate [14].

The method described in this chapter is the most suitable for measuring the transient temperature of gases, such as air or exhaust gases. This is due to the fact that the time constant depends on the heat transfer coefficient on the outer surface of the thermometer and in turn on the Reynolds and Prandtl numbers. For gases such as air or exhaust gas, Prandtl number varies slightly in a wide range of temperatures. However, during measurement of the transient steam temperature, the value of Prandtl number varies significantly depending on the pressure and temperature [23]. In this case, the inverse marching method described in [24] is more appropriate for correction of the dynamic error.

References

  1. 1. Nicholas JV, White DR. Traceable Temperatures. An Introduction to Temperature Measurement and Calibration. 2nd ed. New York: John Wiley & Sons; 2001.
  2. 2. Michalski L, Eckersdorf K, Kucharski J, McGhee J. Temperature Measurement. 2nd ed. Chichester: John Wiley & Sons; 2001.
  3. 3. Wiśniewski S. Temperature Measurement in Engines and Theramal Facilities. Warszawa: WNT; 1983.
  4. 4. Taler J. Theory and Practice of Identification of Heat Transfer Processes. Wrocław: Zakład Narodowy imienia Ossolińskich; 1995.
  5. 5. Kabza Z, Kostryko K, Zator S, Łobzowski A, Szkolnikowski W. Room Climate Control. Warszawa: Agenda Wydawnicza, Pomiary Automatyka Kontrola; 2005.
  6. 6. Littler DJ et al. Instrumentation, Controls & Testing. Modern Power Station Practice. Oxford: Pergamon Press; 1991.
  7. 7. Childs PRN. Practical Temperature Measurement. Oxford: Buterworth-Heinemann; 2001.
  8. 8. Gerashchenko OA, Gordov AN, Lakh VI, Stadnyk BI, Yaryshev NA. Temperaturnye Izmereniya. Kiev: Naukova Dumka; 1984.
  9. 9. Han JC, Dutta S, Ekkad SV.Gas Turbine Heat Transfer and Cooling Technology. 2nd ed. Boca Raton: CRC Press; 2012.
  10. 10. Szekely V, Ress S, Poppe A, Torok S, Magyari D, Benedek Z, Torki K, Courtois B, Rencz M. New approaches in the transient thermal measurements. Microelectronics Journal. 2000;31(9):727–733.
  11. 11. Crocker DS, Parang M. Unsteady temperature measurement in an enclosed thermoconvectively heated air. International Communications in Heat and Mass Transfer. 2001;28(8):1015–1024.
  12. 12. Chau PC. Process Control: A First Course with MATLAB. Cambridge: Cambridge University Press; 2002.
  13. 13. Taler J, Jaremkiewicz M, Taler D, Sobota T. Fluid temperature measurement under transient conditions. Archives of Thermodynamics. 2009;30(3):75–88.
  14. 14. Jaremkiewicz M, Taler D, Sobota T. Measuring transient temperature of the medium in power engineering machines and installations. Applied Thermal Engineering. 2009;29(19):3374–3379.
  15. 15. Jaremkiewicz M. Reduction of dynamic error in measurements of transient fluid temperature. Archives of Thermodynamics. 2011;32(4):55–66.
  16. 16. ASME. Policy on reporting uncertainties in experimental measurements and results. Journal of Heat Transfer. 2000;122:411–413.
  17. 17. Moffat RJ. Describing the uncertainties in experimental results. Experimental Thermal and Fluid Science. 1988;1:3–17.
  18. 18. TableCurve 2D v. 5.1. Automated Curve Fitting & Equation Discovery. San Jose, California: Systat Software Inc.; 2002.
  19. 19. WT4401-S & WT4401-D Benchtop Wind Tunnels. Stamford: Omega; DOI: https://www.omega.com/manuals/manualpdf/M1776.pdf–Accessed July 15, 2016
  20. 20. Sanitjai S, Goldstein RJ. Forced convection heat transfer from a circular cylinder in crossflow to air and liquids. International Journal of Heat and Mass Transfer. 2004;47:4795–4805.
  21. 21. Jaremkiewicz M. Measurement of transient fluid temperature. In: Richard B. Hetnarski, editor. Encyclopedia of Thermal Stresses. Dordrecht: Springer; 2014. p. 2932–2938.
  22. 22. Taler D. The Dynamics of Tubular Heat Exchangers. Cracow: Uczelniane Wydawnictwa Naukowo-Dydaktyczne AGH; 2009.
  23. 23. Meyer CA, McClintock RB, Silvestri GJ, Spencer RC. ASME Steam Tables. New York: The American Society of Mechanical Engineers; 1993.
  24. 24. Jaremkiewicz M, Taler D, Sobota T. Measurement of transient fluid temperature. International Journal of Thermal Sciences. 2015;87:241–250.

Written By

Magdalena Jaremkiewicz

Submitted: 13 April 2016 Reviewed: 09 September 2016 Published: 27 April 2017