1. Introduction
The current work employs a photothermal imaging technique recently developed by the authors [
1]. This procedure is particularly useful for detecting convective heat flow in liquids arising from localized heating by a laser beam. In particular, thermal starting plumes in pure water and salt water were imaged and characterized in spite of the well-known fact that water is considered the worst of the common liquids for thermal lensing methodologies [
2,
3]. Thermal starting plumes are the initial part of the formation of a thermal plume originating from a heat source [
4].
Of particular interest is the behavior of thermal starting plumes near the temperature of maximum density. Under these conditions, the buoyant force nearly matches the weight force. While there is potential for thermal pluming due to local heating, a plumehead does not rise rapidly or at all. It may, in fact, descend. These cases will be referred to here as ‘frustrated’ thermal starting plumes.
The temperature of maximum density in both pure and salt water has a critical effect on the coefficient of thermal expansion. In pure water, this temperature occurs near 4 °C. At this specific temperature, water molecules are arranged in a manner that maximizes their packing efficiency, resulting in the highest density achievable under normal atmospheric pressure. As the temperature decreases below 4 °C, the water molecules form a more ordered structure because of hydrogen bonding, resulting in a decrease in density.
One motivation for this work is to contribute to the understanding of temperature-driven convection near the temperature of maximum density in fresh water lakes during fall and spring turnover [
5]. In temperate lakes, the phenomena of fall and spring turnover play a pivotal role in the maintenance of aquatic ecosystems. As the air temperature decreases in fall, the surface water cools, increases in density, and begins to descend, initiating a mixing process, which continues until the water column reaches a uniform temperature of approximately 4 °C. In the spring, the ice that formed on the surface over winter begins to melt, mirroring the thorough mixing that occurred in the fall. These seasonal turnovers are critical for redistributing oxygen from the surface to the lake’s depths, essential for the survival of various aquatic organisms, and for cycling nutrients from the bottom to the surface.
An additional motivation for this work is the importance of thermal plumes, both natural and human-made, in oceans and other bodies of water [
6,
7,
8,
9,
10,
11]. Understanding thermal plumes is crucial for comprehending various oceanic processes, including ocean circulation, heat transport, and their impact on climate systems [
12]. Specific to pluming, mid-ocean ridges are a source of heat in the form of line plumes [
12]. They also spawn volcanoes and hydrothermal vents. These offer places of distinctive chemical and biological activity [
13]. To that end, the simulated seawater Instant Ocean
® was included in the study in addition to sodium chloride solutions.
In the case of salt water, the presence of dissolved ions, primarily sodium (Na
+) and chloride (Cl
−) ions, alters the temperature of maximum density. This change occurs because the presence of salt disrupts the hydrogen-bonding network among water molecules, preventing them from forming the same dense structure observed in pure water. The temperature of maximum density is a linearly decreasing function of salinity,
S, whose slope is greater in magnitude (more negative) than that for the melting point, which is also a decreasing function of salinity [
14]. The temperature of maximum density reaches the melting point at a salinity of about
[
14] (the salinity units reported throughout this paper are grams of salt per kilogram of solution). For a recent thorough discussion of the importance of the effects and fluid mechanical treatment related to the frustrated thermal starting plumes discussed in the current paper, see the 2022 Ph.D. Thesis of Grace [
15].
Following this Introduction, the various physical principles that underlie both the technique used in this work and the behavior of pure and salt water under the conditions of the reported experiments are discussed along with a semi-empirical model. This is followed by a description of the materials and methods. The results are then presented and discussed. The paper culminates with a conclusion and prospectus on future work.
3. Materials and Methods
In this section, the experimental setup is detailed, the sample preparation is described, and the data analysis procedure is outlined.
3.1. Experimental Setup
The experimental setup, as illustrated in
Figure 3, includes an excitation laser, the QSL103A Q-Switched Picosecond Microchip Laser System by Thorlabs, Newton, NJ, USA emitting light at a 1030 nm wavelength. This laser generates pulses approximately 500 ps long at a repetition rate of approximately 9 kHz, and its beam predominantly follows a Gaussian spatial distribution. Note that, given the timescales involved in the photothermal processes (on the order of seconds), the laser’s operation resembles that of a continuous wave source.
The power reaching the sample, which was largely unabsorbed, was measured using the PM100D Digital Optical Power and Energy Meter from Thorlabs, equipped with an S140C powerhead. This meter has an uncertainty of about 7% at 1030 nm, as specified in the product datasheet. Additionally, it was observed that approximately 7% of the laser’s power was reflected off the front surface of the cuvette. Therefore, the power values reported in the subsequent
Section 4 were adjusted to account for 93% of the initial measurements directly from the laser. A computer-controlled solenoid-actuated shutter enabled the recording of the background (shutter closed), active (shutter open), and recovery (shutter closed) phases during the experiment.
The probe beam was generated by a low-power diode laser from HiLetgo, operating at a 650 nm wavelength. This laser was placed within a passive aluminum heat sink and connected to a potentiometer, allowing for the adjustable irradiation of the sample. While the emitted powers had minimal impact on the sample, this adjustability was crucial for fine-tuning the camera’s exposure settings.
The imaging system utilized an Arducam Camera equipped with a 5-megapixel Omnivision OV5647 sensor. Although capable of operating at a resolution of 1080 p, the camera was configured to function at pixels to achieve a higher frame rate, which was advantageous for capturing dynamic processes in the experiment with greater temporal resolution. The recorded video was encoded in raw .h264 format. Because this format lacks timestamp information in individual frames, time calibration was obtained by recording ten, 10.0 s videos and averaging the total number of frames collected. This process established a time calibration of seconds between successive frames. Monitoring the software’s timing of the shutter opening revealed an estimated timing jitter of fewer than 2 frames. Space calibration of the video data was achieved by introducing a 1.25 mm obstruction at the sample’s position, creating a shadow visible to the camera. Due to the inherent diffraction effects, a subjective assessment was made to identify which pixels corresponded to the obstruction’s edges. The spatial calibration value was 3.54 µm per pixel.
A Raspberry Pi 4 microcomputer was responsible for coordinating the operation of the shutter mechanism and acquiring data from the camera. The Raspberry Pi 4 ran custom-designed software tailored specifically for this experiment (
https://github.com/ulnessd (accessed on 14 December 2023)).
3.2. Sample Preparation
The sample was placed in a 1 cm × 1 cm fluorimeter cell, which was set into a brass jacket built in-house. Ethylene glycol flowing through the jacket was heated/cooled using an FTS model RS33AL10 heater/chiller. Sample temperatures compared to reservoir temperatures were calibrated from 0 °C to 70 °C. The heating/chilling was regulated to ±0.5 °C. For temperatures below 10 °C, argon flowed over the cuvette to prevent condensation.
The water used in this work was obtained from an in-house reverse-osmosis system. The sodium chloride was obtained from Fisher Scientific and used without further purification to prepare salinity values of
,
,
, and
. Instant Ocean
® is a fairly complete seawater substitute used to house marine animals. The composition of Instant Ocean
® is listed in Table 2 of Reference [
53]. There is considerable agreement with the ion concentrations compared to actual seawater reported in Table I of Wiesenburg and Little [
14]. In particular, the top six ions by composition are chloride, sodium, sulfate, magnesium, calcium, and potassium, and they are represented in Instant Ocean
® within 1–2% of that for seawater. These six ions constitute 99.35% of the ions in real ocean water. Missing from Instant Ocean
® are the many trace elements found in the oceans (see Table II of Reference [
14]. Instant Ocean
® was prepared at
.
3.3. Data Collection
For each sample and temperature, 20 runs were recorded sequentially with either a 30 s (≥20 °C) or 90 s (≤15 °C) rest interval between runs, i.e., a given run began its background phase of 30 s or 90 s after the recovery phase of the previous run. During data analysis, careful observation for systematic trends across the 20-run set was made to ensure that the sample was not starting each run at a warmer temperature than the previous.
3.4. Data Analysis
The data analysis in this experiment was conducted using a set of custom-written programs
https://github.com/ulnessd (accessed on 14 December 2023). The data processing involved three main steps. First, the video data were processed. Second, dynamic information related to the starting plumes was extracted and combined with relevant physical data from the existing literature. Finally, the dynamic and physical information were plotted and interpreted.
The video-processing phase followed a structured procedure. Initially, the 20 individual .h264 files were averaged and converted into a single .avi file, which is referred to as the raw sample video and serves as the initial data for the given sample liquid. The top row of
Figure 4 shows an example of this for water at
°C. There, single frames are shown at times
s,
s, and
s. During the first 2 s, the laser was blocked by the shutter; therefore, the images show the raw visual field at
s,
s, and
s after the light first reached the sample. In the image, the laser position is about half-way up the frame and appears as a dark horizontal region. It is dark relative to the spatial background because the heated volume had a lower index of refraction, causing light to diverge from this region. Although somewhat difficult to see, there are two regions that are slightly brighter than the spatial background directly above and below the heated region. Finally, the curved structure at the bottom of the frame is the edge of the probe beam and is of no physical relevance.
Subsequently, a new .avi video was generated by averaging the first 120 frames (equivalent to 1.23 s) of the temporal background from the raw sample video and subtracting this averaged background from all frames in the raw sample video file. The absolute value of the subtracted difference then becomes an individual frame in this second video file referred to as the background-subtracted video. The bottom row of
Figure 4 shows the precise frame in the subtracted video as its corresponding raw video. This process enhances the visualization of both the transient and stationary components of the material’s response. Specifically, the transient feature corresponds to the thermal starting plume, while the stationary component represents the region where the laser interacts with the sample, referred to as the plume source.
For static visualization of the video data in the form of a single figure, the averaged background-subtracted video was compressed into a heatmap. The heatmap image was created by integrating each row in a given frame to produce one individual column in the heatmap. The presentation and analysis of the heatmap representations make up the bulk of the
Section 4. All the video files can be found at
https://www.darinulness.com/research/data/frustrated (accessed on 14 December 2023).
3.5. Artificial Intelligence
It is appropriate to acknowledge the extensive use of OpenAI’s ChatGPT 4.0 in assisting with the development of the programs to acquire and analyze the data [
54]. Use in writing was minimal and limited to paraphrasing previous writings from the authors and producing references in the correct format. Sundry use as background and/or support of the authors included collecting and organizing the papers both used and not used based on the AI’s understanding of the abstracts.
4. Results and Discussion
Results for pure water, an artificial seawater, and several salinities of salt water at a variety of temperatures are presented and discussed.
4.1. Pure Water
As mentioned at the end of the previous section, the results are compressed into heatmaps. By way of example, the heatmap for pure water at 70 °C is shown in
Figure 5. This example offers an illustrative representation of thermal plume dynamics within which several distinct features are apparent. Notably, there is a visually pronounced rising plumehead and source, which manifest as prominent regions of elevated temperature. In contrast, the conduit is less discernible, primarily due to a less abrupt temperature gradient when compared to the leading edge of the plumehead and the laser-heated region. The excitation laser is positioned along the horizontal green line. This line serves as the reference for kinetics traces in subsequent analyses. Beneath this, one observes a significant feature resulting from the thermal lens effect, causing light to bend outward from the cylindrical region. Furthermore, the vertical green line indicates the point of plumehead formation, or the time to form the plumehead,
. Upon shutter closure, the heated fluid within the laser region rapidly ascends, marking the initiation of a thermal ending plume. The ascension angle,
, is denoted by an angled green line, while curly brackets delineate integration regions encompassing asymmetry arising from convection effects. This heatmap provides valuable insights into the complex dynamics of thermal plumes under the specified conditions.
4.1.1. Heatmap Representation of Thermal Starting Plumes
Figure 6 shows heatmaps for pure water at varying temperatures ranging from 11 °C to 70 °C. In this examination of the heatmaps, distinct temperature-dependent features become apparent. As illustrated in
Figure 6, the thermal plume dynamics exhibit notable variations with changing temperature. A key observation is the diminishing time required for plumehead formation and the accelerated ascent of the plumehead as the temperature increases. While the signal-to-noise ratio exhibits a decline at lower temperatures, it remains clearly perceptible above the background level.
In contrast, the investigation of the heatmaps for temperatures ranging from 0 °C to 10 °C reveals a more rapidly changing set of characteristics (
Figure 7). At these lower temperatures, the signal-to-noise ratio exhibits significant degradation, leading to reduced clarity in the observed phenomena.
Figure 8 shows a plot of the total integrated signal for each of the temperatures studied.
Despite the challenging noise conditions at low temperatures, certain features remain discernible. In particular, a rising plumehead is evident, persisting down to 7 °C. A central observation to the main theme of this paper is the rapid decrease in the ascension angle of the thermal ending plume in contrast to the nearly constant ascension angle observed in the higher temperature range depicted in
Figure 6 (the ascension angle is plotted in
Figure 9). Furthermore, indirect evidence of the transition from ascending frustrated starting thermal plumes to descending frustrated thermal starting plumes is discernible through an assessment of the relative intensities of light accumulation above and below the laser position, as observed in
Figure 6.
Serving as a counterexample, the heatmap for silicone oil is shown in the bottom right panel of
Figure 7. One sees here a similar case to the fully frustrated thermal starting plume of water at about 2 °C, albeit a much stronger signal. Convection via thermal pluming also does not appear for silicone oil, but for a different physical reason. The silicone oil is indeed heated in the laser region, and that region is indeed less dense than the ambient fluid. However, because of the very high viscosity of silicone oil, nearly 10,000-times that of water, the viscous drag force effectively stops ascension. For the frustrated thermal starting plumes of water, there is the potential for thermal starting plumes in terms of viscosity, but because the buoyant force matches the weight force, there is no net movement regardless of viscosity.
The transition observed at low temperatures, characterized by a shift from a bright central line (representing the laser position) to a dark central line, warrants the following discussion. A comparison between the heatmaps for 6 °C and 0 °C highlights this transition. This phenomenon arises due to the nature of the heatmaps, which depict the absolute difference from the background. For low temperatures, the thermal lens effect causes the deviation of light from the laser path, resulting in reduced light intensity within the laser region. Given that the background signal is heavily averaged (over 120 frames), its contribution to noise in the heatmap is relatively diminished compared to the individual frames from which it is subtracted. Consequently, noise fluctuations in the laser region exhibit a lower magnitude in comparison to the rest of the visual field, manifesting as a darker central line. Conversely, noise becomes more pronounced in the regions directly above and below the laser position, resulting in brighter lines in the heatmap.
4.1.2. Semi-Empirical Model
The total signal strength as a function of temperature and Equation (
1) can be viewed in the context of a very simplified chemical equilibrium model (dots representing hydrogen bonding):
Using the local mode/combination band model discussed above, the signal strength is determined by the density of ‘free’ (not involved in a hydrogen bond) O–H,
[
55], where
Y is the percent population of free O–H. Unfortunately, while in principle, it is possible to construct a chemical equilibrium model to obtain thermodynamic information, in practice, these thermodynamic values as fit parameters carry too much variance, rendering them useless. Instead, one must settle for a semi-empirical model that fits the lensing-corrected signal strength (
), discussed below, to the Hill equation [
56]. These data reveal the ultimate effect of extended water structures on heat transport. Studies of the extended network of water near and below freezing have been performed by Magazú et al. in the late 1980s through the mid-1990s [
57,
58,
59] and are relevant to the free OH model of Fang and Swofford [
50].
Reconsidering Equation (
1), and hiding the constant laser power, but exposing the source of the absorption, one expresses
where
is the molar absorptivity of free O–H and the constant laser power has been absorbed in the proportionality. Rearranging gives
The measured data for
are collected in
Figure 8. As discussed in
Section 3.4, Data Analysis, each datum represents a 20-run average. Individual error bars are not generated because the lens-corrected total integrated signal was determined from the averaged data set. However, the relative standard deviation is estimated at 2.25% based on an approximate 10% standard deviation of the individual runs, and thus scaled by a factor of
. These data are fit to the Hill equation with Hill coefficient
:
The Hill coefficient is an empirical measure of cooperativity/anti-cooperativity with
representing cooperativity and
representing anti-cooperativity (
recovers a neutral Langmuir behavior). The visual quality of the fit is reasonable for values of between
and
, but
appeared to be the best. The blue and green curves in
Figure 8 show the fits of
for pure water and Instant Ocean
®, respectively. It is also important to note that one must not read too much into the actual values of the parameters as the Hill equation is purely empirical [
56]. This fact, combined with the chemical model above, is the rationale for calling the model semi-empirical.
Although the Hill equation is used to fit cooperative ligand binding systems, the analogy here is that the ligand concentration is replaced with temperature. Temperature is serving as a surrogate for bound and unbound waters per Equation (
2).
The ‘cooperativity’ stems from the various states of hydrogen bonding that water can be in as indicated in the simplified chemical reaction of Equation (
2). At low temperatures, the equilibrium shifts to the right, quenching the signal, because the O–H stretch is both diminished and shifted [
55]. Conversely, at high temperatures, the equilibrium shifts to the left, liberating O–H bonds and, thus, enhancing the signal. This behavior is readily apparent upon inspection of
Figure 6 and
Figure 7.
4.1.3. Asymmetry in Steady-State and Thermal Ending Plumes
There is a visually apparent asymmetry as one views integration region 1 from
Figure 5. For ascending thermal starting plumes, there is a brighter streak below the central streak of the laser region. The converse is true for descending plumes in which the streak above the laser line is brighter. The reason for this is that the conduit in either case acts to reduce the spatial temperature gradient, thereby diminishing the thermal lensing effect.
There is a way to better discern the asymmetry, which is by plotting the integrated signal strength within integration region 1 from
Figure 5. The graphs for these values for temperatures near the temperature of maximum density are collected in
Figure 10. The graph of
Figure 10 shows a comparison between temperatures of 0 °C and 5 °C. At 0 °C, the peak on the positive side of zero (the laser position) exceeds that on the negative side. Visual scrutiny of the asymmetry observed in the peaks on the left and right sides suggests the occurrence of the fully frustrated thermal starting plume to transpire between 1 °C and 2 °C. The authors estimate the fully frustrated plume occurs at
°C.
Advantageously, this allows for an estimate of the laser-induced temperature increase. By utilizing the empirical formula of reference [
14] and solving
for the temperature in the laser region,
, and subtracting the temperature of the bulk,
, give
°C.
The left panel of
Figure 11 illustrates the thermal ending plumes across a temperature range from
°C to
°C. The observations reveal a transition from an ascending frustrated thermal starting plume at
°C to a fully frustrated thermal starting plume occurring near
°C, followed by a descending frustrated thermal starting plume at
°C. To quantify these observations, quadratic functions were fit to the troughs within a narrow range near their minima. The minima extracted from these fits were then plotted against temperature, as shown in the right-hand graph. The trend is remarkably linear, giving insight into the temperature-dependent behavior of thermal starting plumes near the temperature of maximum density.
The kinetics traces at the position of the laser (green horizontal line in
Figure 5) are shown for water in the left graph of
Figure 12. As temperature is reduced, the peak in the traces grows smaller and wider and occurs later. This is consistent with the reduced signal due to increasing hydrogen bonding, as well as an increase in viscosity as temperature decreases. The time of the peak is the time to form the plumehead [
21].
The time to form the plumehead is plotted for both water (blue) and Instant Ocean
® (green) in the right graph of
Figure 12. Notably, water follows what appears to be a
behavior, whereas Instant Ocean
® appears more linear with temperature. This is a result of the softening of the effects of hydrogen bonding with temperature for the high salt condition. This follows from the semi-empirical model discussed above and
Figure 8.
4.2. Salt Water
As noted in the Introduction, water is a poor liquid for thermal-lens-based methods due to its low magnitude of
and its high thermal conductivity [
2,
3]. As discussed in
Section 2.2.1, hydrogen bonding inordinately quenches absorption into overtone bands of O–H stretching and presumably shifts them out the laser excitation energy as well. The addition of salt can lead to a decrease in hydrogen bonding. This is the case for sodium chloride with a salinity-dependent enhancement [
3], evident in
Figure 8, where the Hill equation fit shows a significantly softer sigmoidal curve with lower amplitude than for water.
The salinity study was primarily focused on temperatures near that of maximum density. For all salinities, data were collected at temperatures of , 10, 8, 6, 4, 2, and 0 °C. Additionally, °C was used for all samples, except pure water and . The Instant Ocean® data were recorded in descending steps of 5 °C from 70 °C to 10 °C in addition to the temperatures mentioned above. Finally, temperatures of 25 °C and 30 °C were used for the and samples.
In the analysis of the heatmaps for a temperature of 0 °C across the range of salinities investigated, one sees the expected behavior of the thermal starting plumes. Specifically, one sees both a descending thermal ending plume for pure water and
and a fully frustrated thermal ending plume for
. Ascending thermal ending plumes are seen for
,
, and Instant Ocean
®. Moreover, asymmetry in integration region 1 changes from stronger above the laser region to stronger below it with increasing salinity. These observations are presented in
Figure 13.
For salinities below 23, the temperature of maximum density surpasses the melting point. This characteristic gives rise to a unique behavior in the thermal plumes [
14]. Specifically, when salinity is set at 0 and 5, descending frustrated thermal starting plumes become evident.
Figure 11 shows that, for
, fully frustrated thermal starting plumes occur at
°C, and for
, they occur at
°C. Conversely, the scenario changes for a salinity level of 14.5, where the system exhibits a fully frustrated thermal starting plume. These values suggest laser-induced heating of
for
and
for
.
In contrast, for salinities of 25.1 and 35.0 and Instant Ocean® (), ascending frustrated thermal starting plumes dominate the observed behavior. Under these conditions, all thermal starting plumes are less dense than the ambient fluid.
5. Conclusions
The findings presented in this study provide insight into the nature of thermal plume dynamics, particularly for pure and saltwater systems near their respective temperatures of maximum density. Analysis employing the photothermal imaging technique has revealed critical insights into the behavior of frustrated thermal starting plumes. These plumes, characterized by their almost equivalent buoyant and weight forces near the temperature of maximum density, exhibit unique patterns of movement, ranging from stagnation to descending motion, a stark contrast to the typical ascending thermal plumes seen in all other common liquids.
In pure water, the occurrence of maximum density at approximately 4 °C plays a unique role in thermal starting plume dynamics compared to other molecular liquids. This study has demonstrated that the transition from ascending to descending convection, as seen at about 2.3 °C, exemplifies the more complex interplay of temperature, density, and buoyancy in pure water compared to other molecular liquids. The saltwater scenarios add another layer of complexity. The presence of Na+ and Cl− ions disrupts the hydrogen bonding among water molecules, leading to a shift in the temperature of maximum density. For low salinities, there remains a significant-enough gap between the melting point and the temperature of maximum density. Ascending, descending, and fully frustrated thermal starting plumes were observed. For higher values of salinity, only ascending thermal starting plumes can occur. This understanding has potential implications for the study of natural and human-made thermal plumes in various water bodies.
In the future, investigators could (i) consider more extreme temperature and pressure situations, (ii) study the effect of dissolved organic matter, in particular relatively stronger absorbers such as alcohols and amines, (iii) study plumes in nonhomogeneous fluids such as layered salt solutions, (iv) develop protocols such as producing a thermal followed by a thermal starting plume [
60], and (v) employ computational fluid dynamics simulations.