1. Introduction
Thermal hydraulics is the phenomenon whereby the temperature and flow fields interact with each other. Examples of the strong coupling between temperature and flow fields include natural convection flow and heat transfer. The specific applications are air conditioning systems in buildings [
1], thermal storage tanks [
2,
3], and passive safety systems in nuclear power plants [
4]. Therefore, the measuring of temperature and velocity fields is crucial to understanding the thermal hydraulics in these applications.
Numerous measurement methods have been developed. Thermocouples for temperature and hot wire anemometers [
5] for flow velocity are widely used because they are simple and accurate. On the other hand, immersing these sensors in the measurement domain disturbs the temperature and flow fields, resulting in only a few selected points being measured [
6]. Two optical methods are combined to measure temperature and flow fields: Hanwook Park [
7] and Stig Grafsrønningen et al. [
8] used PIV (Particle Image Velocimetry) and Laser Induced Fluorescence with phosphor tracer particles, and In Mei Sou et al. [
9] used PIV and Thermography. Although their optical methods explored temperature and velocity fields, measurement calibration was complex due to significant changes in the refractive index and for optical visibility through windows, three-dimensional shapes, and large temperature gradients [
8]. In addition, when tracer particles were mixed into the measurement, the tracking of the particles became excessive. When the circulation force was weak, the tracer might have settled out due to the long time measuring and natural convection. Another measurement method is distribution measurement using optical fibers. Lomperski et al. [
10] measured the temperature distribution, although the velocity distribution measurement has not been thoroughly investigated. The simultaneous measurement of temperature and velocity has not been investigated.
We have developed a new temperature and velocity measurement with a high temporal and spatial resolution to overcome these drawbacks using an optical fiber thermometer. The new method can measure three-dimensional temperature and flow velocity distributions even in closed systems without significantly disturbing them.
Section 2 describes the temperature and velocity measurement methods using optical fibers.
Section 3 validates the measurement performance using the airflow experiment results.
Section 4 describes the condition setting of this measurement method. Finally,
Section 5 discusses an applied example of the new method developed in this study.
2. Measurement Method
This section describes the newly developed measurement method. The temperature was measured by an optical fiber measurement method. The flow velocity was measured by the temperature change, similar to with a hot wire anemometer. Simultaneous temperature and velocity measurements could be obtained by combining the above two measurements.
Figure 1 shows the schematic of the temperature and velocity measurement devices. For the flow velocity measurement, the following modifications were made to the standard optical fiber thermometer system: (1) coating the optical fiber with gold, and (2) adding devices for pulsed Joule heating of the gold-coated fiber (see the red-dashed boxes in
Figure 1).
2.1. Temperature Measurement Method
In this study, the Optical Distributed Sensor Interrogator B (ODiSI-B) from Luna Innovation Inc. was used as the optical fiber measurement device. The measurement device consisted of a tunable laser device, a detector, and an optical fiber sensor, as shown in
Figure 1. The frequency shift amount
of the Rayleigh backscattered light was detected by ODiSI-B based on the optical frequency domain reflectometry method [
11]. The measurement performance of ODiSI-B can measure the amount of frequency shift at a measurement frequency of 100 Hz and a measurement interval of 2.56 mm. The sensor part was calibrated for position and set to 2 m. The fiber diameter was approximately 0.16 mm (+0.002/−0.006).
The frequency shift amount
is a function of the temperature difference
and the strain
of the optical fiber from a reference state. The relationship is given by the following equation [
6,
10],
where
and
are proportional constants, which are determined in the calibration process. If the fiber strain is kept unchanged, the temperature can be estimated from Equation (1) and
. Since the values of
are obtained at nearly continuous spatial points along the fiber, the temperature can also be measured at the same spatial points.
Lomperski et al. [
10] and Sang et al. [
12] demonstrated the ability to measure temperature distributions by using optical fibers under high radiation conditions in spent fuel pools of nuclear reactor facilities. Kim et al. [
13] demonstrated the applicability of the measurement method at the operating range of 30–70 °C. Weathered et al. [
14] and Arora et al. [
15] applied the measurement method at high temperatures (up to 600 °C) in sodium and molten salt, respectively.
2.2. Velocity Measurement Method
Velocity across an optical fiber is estimated based on the principle of hot wire anemometry [
16]. The gold-coated optical fiber was electrically heated for velocity measurement by a DC power supply with a function generator signal. The temperature of the optical fiber was measured as described in the previous section.
Figure 2 shows the temperature change of the optical fiber when pulsed Joule heating was applied to the fiber. The temperature rose from the ambient fluid temperature
during the heating period and returned to
during the non-heating period.
If the heat conduction in the optical fiber is neglected, the temperature decay process during the non-heating period is expressed by the following energy conservation equation:
where
is the heat transfer coefficient the between ambient fluid and the optical fiber,
are the density, specific heat capacity, temperature, and radius of the optical fiber, respectively. The general solution of Equation (2) is presented by:
where
is a constant that depends on the applied heat, and
Since depends on the flow velocity, also depends on the flow velocity from Equation (4). Thus, the flow velocity can be estimated from the value, which is obtained by fitting the measured temperature with Equation (3).
2.3. Temperature and Velocity Simultaneous Measurement Method
As shown in
Section 2.2, the flow velocity was estimated by measuring the temperature change during the temperature decay process after pulsed Joule heating. After the temperature decay process, the temperature approached the ambient fluid temperature
, which is the temperature measurement of our method. Thus, one temperature value and one velocity value were measured per pulse of Joule heating.
3. Calibration and Measurement Performance
This section describes the calibration for the temperature and velocity measurement methods presented in
Section 2 and the measurement performance.
3.1. Temperature Measurement
The values of the constants in Equation (1) depend on the type of fiber and the temperature range. Calibration experiments were conducted by measuring the water temperature in which the gold-coated optical fiber was immersed to determine the value of
by comparing it with the reference temperature obtained from the Type-K thermocouple (measurement error of ±2.5 °C). The temperature varied in the range of 20–90 °C. As a result,
was obtained, and thus the temperature was estimated by the measured
and the following equation,
where
is the base temperature.
Lomperski et al. evaluated the measurement error of optical fiber thermometers [
6,
10,
15]. The measurement error of the optical fiber thermometer was reported to be ±3.38 °C [
15]. In our experiments, the measurement error was also found to be
. Here, the measurement error of the reference thermocouple was ±2.5 °C and that of the gold-coated fiber was ±2.8 °C.
3.2. Velocity Measurement
Figure 3 shows the experimental setup for the velocity calibration. The gold-coated optical fiber sensor was installed in a 0.1 m diameter pipe. A honeycomb structure was installed at the inlet, and a small fan (SANYO DENKI CO., LTD., Tokyo, Japan) was installed at the outlet to stabilize the airflow. The hot wire anemometer (KANOMAX Inc., Shimizu Suita City, Japan) was used to measure the wind velocity as a reference, and the measurement error was greater than ±3% of the displayed value or 0.02 m/s.
A small fan was used to generate a steady flow for the velocity calibration, and the flow velocity was measured with the hot wire anemometer. At the same time, the time constant
of the temperature decay process after Joule heating was measured as described in
Section 2.2.
Figure 4 shows the relationship between the reference velocity
measured by the hot wire anemometer and the time constant
. For the experimental range shown in
Table 1, the flow velocity
is accurately approximated by the following expression for powers of
,
The above multiplier relationship can be understood from the fact that Equation (3) gives
and the
number is given by the power of the Re number proportional to the flow velocity. For example, the Zukauskas model [
17] for flow around a cylinder shows
, giving
, which is consistent with Equation (6).
Figure 5 shows the results of the comparison between the flow velocity
measured with the optical fiber using Equation (6), and the reference velocity
measured with the hot wire anemometer. The green line is the standard deviation, and the flow velocity can be measured by the optical fiber with approximately 10% deviation from the hot wire anemometer. From the above, the temperature decay process after heating depends on the wind velocity, and the prospect of velocity–current measurement using τ was obtained.
3.3. Simultaneous Temperature and Velocity Measurements
The simultaneous measurement method described in
Section 2.3 was applied to the airflow experiment shown in
Figure 3. The reference temperature was provided by the same K-type thermocouple used in the temperature calibration, and the reference velocity was provided by the hot wire anemometer used in the velocity calibration. Regarding experimental conditions, the electric power was 21 kW/m
2 with heating for 0.2 s, with the output turned on every 2 s. The experiment started when the fan started, and the wind velocity stabilized at 1 m/s. As shown in
Figure 6, the temperature and velocity measured once per pulse matched the reference values from the thermocouple and the hot wire anemometer. The response time was 2 s, sufficient for the measurement prospect.
The measurement errors of the simultaneous measurements are summarized in
Table 2. Although the velocity measurement error was the same as the value obtained from the velocity calibration, the temperature measurement error was smaller than the value obtained from the temperature calibration. This is because the temperature in the simultaneous measurements was evaluated by fitting the temperature change during the temperature decay process.
4. Discussion
This section discusses the setting of the measurement parameters and the effect of wind velocity on the fiber.
4.1. Setting of Measurement Parameters
It was necessary to set the power and the heating duration time of the Joule heating applied to the fiber to perform the velocity measurement described in
Section 2.2. From the point of view of minimizing the effect on the fluid to be measured, it was desirable to use low power and a short time for the Joule heating. On the other hand, from the point of view of maintaining the accuracy of the fit when determining
, it was desirable to increase the power or heating duration time of the Joule heating to increase the amount of heat applied to the fiber and to increase the temperature response.
Figure 7 shows the fiber’s temperature response as the power varied between 7, 11, 21, 32, and 42 kW/m
2 with a heating duration of 0.2 s for the Joule heating. The smaller the power applied to the fiber, the smaller the temperature response to pulse heating and the larger the temperature noise effect. As shown in
Table 3, the coefficient of determination R
2 for obtaining τ was degraded under low power conditions (7 kW/m
2), so it was set to 11 kW/m
2 or higher in the experiments in this paper. Since the fiber temperature returns to the ambient temperature faster at lower power, a lower power should be selected when the interval between pulse heating, i.e., the interval between velocity measurements, is small.
In this paper, to confirm the temperature response to pulse heating with high accuracy, the heating duration was set to 0.2 s, which was sufficiently longer than the 0.01 s interval between temperature measurements of the optical fiber. However, from the point of view of shortening the interval between velocity measurements, the heating duration should be shorter.
4.2. Effect of Wind Velocity on Fiber Measurement
Equation (1) indicates that the frequency shift of the optical fiber is affected by stress. However, the temperature calibration equation in Equation (5) was calculated under the assumption that the optical fiber is not under stress. Therefore, we checked the effect of the stress received from the airflow using the experimental apparatus shown in
Figure 3. The amount of frequency shift when the wind velocity across the optical fiber was varied from 0.0 (stagnant) to 15.0 m/s without heating is shown in
Figure 8. The frequency shift was less than 1 GHz (equivalent to 0.67 K) regardless of the flow velocity. This indicates that the stress caused by the wind velocity did not affect the temperature measurement.
5. Application Example of Optic Fiber Velocity Measurement System
This section introduces the characteristics of using optical fiber: distribution measurements and time variation of velocity distribution.
5.1. Measurement of Wind Velocity Distribution
The calibration equation obtained from a single measurement point on the optical fiber was applied to the entire fiber to measure the velocity distribution. Three measurements were taken each time at wind velocities of 1 m/s, 2 m/s, and 3 m/s in the experimental setup shown in
Figure 3. A range of 70 mm was measured because the area near the electrode could not be accurately measured.
Figure 9 shows the wind velocity and hot wire anemometer with the optical fiber. The optical fiber could measure velocities at 2.56 mm intervals, with a tendency for higher velocities in the center of the cylinder and slower velocities near the walls. As a trend in wind velocity, there was a slight variation in the optical fiber measurements as the velocity increased. However, compared to the three hot wire anemometers, the measurement was within ±5%, which is less than the ±10% velocity measurement error of the optical fiber.
From the above, the optical fiber velocity meter could measure the velocity distribution within the measurement error of a hot wire anemometer in the measurement range. This is an advantage not found in conventional measuring instruments. On the other hand, the measurement error of this instrument tended to be larger than that of the hot wire velocimeters. This can be attributed to several factors. The main factors considered were the uncertainty in the temperature conversion equation of thermocouples, the uncertainty in calculating the decay time constant, the uncertainty in the velocity calibration equation, and the uncertainty in the gold layer thickness.
5.2. Measurement of Wind Velocity Transients
Figure 10 shows the change in velocity distribution when the fan changed the wind direction from the left side to the right side of the optical fiber in 10 s using the experimental airflow apparatus and the conditions shown in
Figure 3 in
Section 3.2. When the fan was fixed at the center position of the fiber using a hot wire anemometer, the wind velocity was 2.5 m/s. On the other hand, the velocity distribution using the optical fiber was smaller than 2.5 m/s overall because the fan was moved in the wind direction. However, the velocity distribution was observed to move from left to right.
From the above, the optical fiber anemometer could capture velocity changes on the order of seconds. Therefore, we qualitatively obtained the prospect that the anemometer can be applied as a measurement device for transient changes.
6. Conclusions
Using an optical fiber thermometer, a new high-resolution, multi-point temperature and velocity measurement method was developed. The scope of application of optical fiber velocimeters and their applicability to distribution, temperature, and velocity measurements were investigated. For the velocity measurement of optical fibers, we applied the principle of hot wire anemometry to temporarily increase the temperature of an optical fiber by pulsed Joule heating of a gold-coated optical fiber. Then, we investigated the velocity measurement method from the temperature time series data of the optical fiber. As a result, the velocity measurement method focusing on the temperature decay process was more effective in the experimental range of 0.2–6.8 m/s, with an error margin of approximately 10% compared to hot wire velocimetry. These results showed that the velocity measurement method focusing on the temperature decay process was adequate for the present experimental range.
The optical fiber velocity measurement method focusing on the temperature decay process confirmed that simultaneous temperature and velocity measurement could be performed every 2 s. Furthermore, the velocity distribution could be measured at 2.56 mm intervals. As a transient change in the velocity distribution, it could be qualitatively captured for wind velocities changing every second, giving us the prospect of applying this measurement system to transient changes.
Future research will focus on its applicability to different fluids, velocity and temperature range limits, and transient temperature and velocity changes.
Author Contributions
M.S.: Conceptualization, Methodology, Data curation, Validation, Writing—Original draft preparation. M.F.: Methodology, Writing—Review and editing. All authors have read and agreed to the published version of the manuscript.
Funding
This research received no external funding.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
Not applicable.
Acknowledgments
The authors are grateful for many constructive suggestions made by Naofumi Tsukamoto and Yasuhiro Masuhara of NRA, Japan. Thanks are extended to Kouji Suzuki, Kenji Kusunoki of Techno Service Co., Ltd., and Maki Kurashima of Waseda University for helping in the experiments for this study.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Van Hooff, T.; Blocken, B. Full-Scale Measurements of Indoor Environmental Conditions and Natural Ventilation in a Large Semi-enclosed Stadium: Possibilities and Limitations for CFD Validation. J. Wind Eng. Ind. Aerodyn. 2012, 104–106, 330–341. [Google Scholar] [CrossRef]
- Altuntop, N.; Arslan, M.; Ozceyhan, V.; Kanoglu, M. Effect of Obstacles on Thermal Stratification in Hot Water Storage Tanks. Appl. Therm. Eng. 2005, 25, 2285–2298. [Google Scholar] [CrossRef]
- Shin, M.; Kim, H.; Jang, D.; Lee, S.; Lee, Y.; Yoon, H. Numerical and Experimental Study on the Design of a Stratified Thermal Storage System. Appl. Therm. Eng. 2004, 24, 17–27. [Google Scholar] [CrossRef]
- International Atomic Energy Agency. Natural Circulation Phenomena and Modelling for Advanced Water Cooled Reactors, IAEA-TECDOC-1677; IAEA: Vienna, Austria, 2012. [Google Scholar]
- Bruun, H.H. Hot Wire Anemometry; Oxford University Press: Oxford, UK, 1995. [Google Scholar]
- Holler, D.; Vaghetto, R.; Hassan, Y. High-Resolution Wall Temperature Measurements with Distributed Fiber Optic Sensors. Int. J. Therm. Sci. 2019, 145, 106042. [Google Scholar] [CrossRef]
- Park, H.; Park, J.; Jung, S.Y. Measurements of Velocity and Temperature Fields in Natural Convective Flows. Int. J. Heat Mass Transf. 2019, 139, 293–302. [Google Scholar] [CrossRef]
- Grafsrønningen, S.; Jensen, A. Simultaneous PIV/LIF Measurements of a Transitional Buoyant Plume Above a Horizontal Cylinder. Int. J. Heat Mass Transf. 2012, 55, 4195–4206. [Google Scholar] [CrossRef]
- In Sou, I.M.; Allen, J.S., III; Layman, C.N.; Ray, C. A Synchronized Particle Image Velocimetry and Infrared Thermography Technique Applied to an Acoustic Streaming Flow. Exp. Fluids 2011, 51, 1201–1208. [Google Scholar] [CrossRef]
- Lomperski, S.; Gerardi, C.; Pointer, W.D. Fiberoptic Distributed Temperature Sensor Mapping of a Jet-Mixing Flow Field. Exp. Fluids 2015, 56, 55. [Google Scholar] [CrossRef]
- Soller, B.J.; Gifford, D.; Wolfe, M.; Froggatt, M. High Resolution Optical Frequency Domain Reflectometry for Characterization of Components and Assemblies. Opt. Express 2005, 13, 666–674. [Google Scholar] [CrossRef]
- Sang, A.K.; Froggatt, M.E.; Gifford, D.K.; Kreger, S.T.; Dickerson, B.D. One Centimeter Spatial Resolution Temperature Measurements in a Nuclear Reactor Using Rayleigh Scatter in Optical Fiber. IEEE Sens. J. 2008, 8, 1375–1380. [Google Scholar] [CrossRef]
- Kim, R.; Park, C.H.; Yoo, W.J.; Moon, J.H. Development and Characterization of a Fiber-Optic Monitoring System for the Key Environment Variables of the Spent Nuclear Fuel Pool at a Nuclear Power Plant. Ann. Nucl. Energy 2017, 99, 183–192. [Google Scholar] [CrossRef]
- Weathered, M.; Anderson, M. Development of optical fiber instrumentation for use in sodium cooled fast reactors. Proc. Nureth 2015, 16, 1050–1059. [Google Scholar]
- Arora, O.; Lancaster, B.; Yang, S.R.; Vaghetto, R.; Hassan, Y. Advanced Flow and Temperature Measurements in a Forced Convection Molten Salt Test Loop. Ann. Nucl. Energy 2021, 159, 108269. [Google Scholar] [CrossRef]
- Lundström, H. Investigation of Heat Transfer from Thin Wires in Air and a New Method for Temperature Correction of Hot-Wire Anemometers. Exp. Therm. Fluid Sci. 2021, 128, 110403. [Google Scholar] [CrossRef]
- Žukauskas, A. Heat Transfer from Tubes in Crossflow. Adv. Heat Transf. 1972, 8, 93–160. [Google Scholar] [CrossRef]
| Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).