Fluid surface tension evaluation using capillary wave measurement with optical coherence tomography

The surface tension of biological fluids is an important parameter because the mechanical properties of fluids are closely linked with hematological diseases and other pathophysiologies. Capillary waves are associated with fluid mechanical properties. Here, we propose a method that utilizes the acoustic radiation force (ARF) to generate propagating waves and optical coherence tomography (OCT) to measure the wave motion. This ARF-OCT method is capable of evaluating the surface tension of fluids, water and porcine whole blood in this study, based on the dispersion relation of capillary waves. Two-dimensional Fourier transforms were used to decompose frequency components of wave motion images to obtain a k-space representation and estimate the wave phase velocity. The phase velocity of capillary waves was obtained from the experimental results and compared to theoretical calculations. The surface tensions of water and porcine whole blood were determined from the experimental results. We first report that capillary waves measured with OCT can be a new promising modality for measuring the surface tension of fluids. The proposed method could be used to differentiate actual pathologic fluids or blood from those taken from healthy subjects and as a biomarker in future biomedical applications.

Surface tension is derived from cohesive forces of molecules pulled equally in every direction by neighboring liquid molecules so that the net force is zero (equilibrium). 1 Surface tension properties, defined as energy per unit area (dyn/cm) or force per unit length (N/m), 2 of biological liquids are closely linked with various physiological processes. [3][4][5] Pathological processes such as lung carcinogenesis 6 influence the surface tension of alveolar sacs, 7 and blood clot formation is associated with the surface tension of blood. 8 Therefore, the determination of the surface tension of biological fluids could be beneficial for understanding the pathological steps and physiological conditions. [3][4][5] The DuNoüy ring method and the Wilhelmy plate method are two classical techniques for measuring the surface tension at a liquid-air interface. 1,2,9 However, extra correction factors (Harkins and Jordan) and large volumes of studied solutions are limitations of these methods. 2 The droplet method is also a well-known technique to measure the surface tension, but the limitation to this method is related to the concentration of the surfactants due to the exhibition of relatively weak surface activity. 10 Kalantarian et al. proposed an image analysis method named axisymmetric drop shape analysis (ADSA) to estimate the surface tension. 11 However, a high image quality is required to be able to define the boundary of drops. Laser-based optical methods such as optical traps 12 and CO 2 laser 13 were proposed to estimate the surface tension of protein condensates in the μN/m scale and soft glass fibers at a high temperature (at least 1000 ○ C) for industrial applications, respectively. In addition, a He-Ne laser was used to measure the light reflection from the liquid surface to determine the contact angle and surface tension. 14 For microscopy-based methods, Bertocchi et al. determined the surface curvature of drops discriminated by backreflected light to determine the surface tension. 15 For both laserbased and microscopy-based methods, a delicate calibration process will always be necessary because the surface curvature is a critical parameter. 15 Optical coherence tomography (OCT) has been used to evaluate the mechanical properties of solid biological materials due to high resolution, the direct measurement of the wave propagation, and noncontact operation. [16][17][18][19] However, it has not been widely studied to investigate the rheological properties in liquids. Capillarywave techniques are being developed to obtain more insight into surface properties. [20][21][22][23] Our previous report demonstrates that capillary waves can be generated by an acoustic radiation force (ARF) and monitored by OCT. 24 In this study, we first report that OCT could be a new modality to evaluate the mechanical property, surface tension, of liquids based on the characterization of the dispersion of capillary waves in frequency dependent rheology. The perturbation of capillary waves can be completely monitored due to a relatively large field-of-view of OCT compared with laser-based methods. 20,23 In addition, we characterize biological fluid, porcine whole blood, in this paper. Acoustic radiation forces were exerted on the fluid-air interface to generate a capillary wave, which can be imaged similar to watching ripples on a water surface created by dropping a stone in a pond, assuming that the depth of the water is smaller than half of a wavelength. 25 The theoretical calculations of the surface tension were compared with our experimental results.
In Navier-Stokes theory, the surface tension is associated with the phase velocities of capillary waves Cp, fluid depth d, and fluid density ρ. 4,20,26 A capillary wave is a surface wave propagating along the interface of fluids between a liquid and air. In this study, we consider a fluid-air free surface with small amplitude waves and a fluid that is assumed to be of uniform depth, incompressible, and irrotational. The wave propagation in the x-direction with a phase velocity Cp can be described by the following equation: 25 where g is the gravity, k is the wavenumber (2π/λ), d is the depth of the fluid, ω is the angular frequency, σ is the surface tension, and Δρ is the mass density difference between two media. ρ and ρ ′ are the mass densities of the liquid and air, respectively. The Eötvös (Eo) number (or Bond number) is a factor to characterize the regime between the gravity and capillary waves, where λ capillary = √ σ/Δρg is the capillary length and Δρ is the difference in densities between the liquid and air. For a liquid-air interface, if λ is much smaller than 2πλ capillary (Eo ≪ 1), the wave is classified as a capillary wave. Figure 1 illustrates the theoretical calculation of phase velocities Cp of water and porcine whole blood for capillary waves with Eq. (1). We assume surface tension values of 73 mN/m for water 27 and 55 mN/m for whole blood. 4 The mass densities of water and whole blood are 1000 kg/m 3 and 1060 kg/m 3 , respectively. 28,29 The mass density of air is 1.2 kg/m 3 , and the gravity is 9.81 m/s 2 . The fluid depth d is assumed to be 3 mm. By substituting these parameters into Eq. (1), the theoretical phase velocity of capillary waves in whole blood (red curve) and water (blue curve) can be obtained, as shown in Fig. 1. Therefore, the critical wavelength λc in water and in whole blood was evaluated to be 17.25 mm (blue curve in Fig. 1) and 14.45 mm (red curve in Fig. 1), respectively. The corresponding minimum phase velocity C p(min) in water and in whole blood was 0.224 m/s and 0.209 m/s, respectively. Once the wavelength is smaller than the critical wavelength, the waves will enter the capillary regime and the gravity is considered to be negligible. 26,30 The dispersion relation in (1) can be written as Due to the liquid-air interface, the first term on the right side of Eq. (3) is approximately equal to 1 and Δρ is approximately 1000 kg/m 3 . The deep water and shallow water regimes are defined as the regimes where d is larger than 0.5λ and smaller than 0.05λ, respectively. 25 In our case, the depth of fluids was considered to be 3-5 mm and λ was approximately 6 mm (as will be presented later); therefore, the capillary waves are considered to be in the deep water (d > 0.5λ) regime in the most cases, and tanh(kd) = tanh(2πd/λ) is close to 1. Capillary waves in the shallow water regime are not considered in this study because the depth d of the liquid would need to reach to a few hundred micrometers. We rewrite Eq. (3), and the surface tension based on experimental parameters σexp can be evaluated by the following equation: In practice, kexp can be obtained by the k-space calculated by twodimensional Fourier transform (2D-FT) 31 of capillary wave motions in the temporal space monitored by OCT. The capillary wave phase velocity Cp_exp and the surface tension σexp in water and in whole blood can be estimated according to Eqs. (1) and (4), respectively. The details are described in the following paragraphs.
The whole blood preparation is followed by a standard protocol. 32 The porcine whole blood with the K3-EDTA anticoagulant was purchased from LAMPIRE Biological Laboratories, Inc. (Pipersville, PA, USA). The 1 l bottle with porcine whole blood was inverted 5 times gently and then 4 ml of blood was transferred to the Petri dish for the experiment. The bottom of the Petri dish was replaced with a Mylar film. The thickness of the Mylar film is only 100 μm, and the film is acoustically transparent than the plastic of the Petri dish. 33 Therefore, the acoustic power will be only weakly attenuated. The tap water from our laboratory was used in this study. The depth of fluids was around 3-5 mm. The temperature of water and whole blood was 21 ○ C (room temperature). The capillary waves on water and porcine whole blood were recorded by using an acoustic radiation force optical coherence tomography (ARF-OCT) system. It was composed of a spectral domain OCT (SD-OCT) scanner with 1300 nm central wavelength (TEL320C1, Thorlabs, Inc., Newton, NJ, USA), a 7.5 MHz focused transducer (ISO703HR, Valpey-Fisher, Hopkinton, MA, USA), three function generators (33250A, Agilent, Santa Clara, CA, USA), and a radiofrequency (RF) amplifier (240L, Electronics and Innovation, Ltd., Rochester, NY, USA), as illustrated in Fig. 2. Three function generators were employed to manage the whole system. Acoustic radiation force was initially generated by function generator 2 and was amplified to 50 dB by a radiofrequency (RF) power amplifier to drive the transducer. Function generator 1 was used to trigger excitation and provide synchronization for the whole system. Function generator 3 was set to control the OCT scan rate at 10 kHz with square pulse trains. More details about the excitation methods can be found in a previous article. 24 The OCT system is capable of producing 13 μm lateral resolution and 3.6 mm of penetration depth within 10 × 10 mm 2 field of view (FOV) to detect the motion of particles on the surface. The 3.6 mm is provided by Thorlabs, Inc., which is the maximum imaging depth in air that the system displays. The penetration depth is governed by the level of optical scattering and absorption for a given sample. In the whole blood, this depth was on the order of 0.3 mm. To be able to capture the propagation of capillary waves on water and whole blood, the M-B scan mode with 10 kHz scan rate was used to track dynamic processes in the space-time domain. The customized acquisition was set to obtain the data at 100 lateral positions with 0.1 mm interval within a 10 mm FOV and 500 axial scans at 10 kHz scan rate at each position. The acoustic radiation force generated with a 2 ms excitation signal repeated every 50 ms was transmitted to excite the fluids to generate capillary waves. Our previous research reports that the double sideband suppressed carrier amplitude modulation (DSB-SC AM) technique was a useful tool for characterizing the viscoelastic properties of a tissue mimicking a phantom material. 34 Therefore, it was adopted in the study to create the capillary waves on water and whole blood. In this study, ten cases were collected for water and 20 cases were collected for blood.
The capillary wave motions were determined by the autocorrelation method. 35 Before autocorrelation, a median filter was employed to remove noise from an original OCT wave motion image. The data in each column of the wave motion image were upsampled by five times by spline interpolation so that the maximum of the motion signal and the time of that peak can be found at each spatial location. After the autocorrelation calculation, spatiotemporal wave motion images were reconstructed, as shown in Fig. 3(a) for water and Fig. 3(b) for the porcine whole blood, with white dashed pentagons. A 2D-FT was used to decompose multiple frequencies from wave motion images of water and whole blood to obtain a k-space illustrated in Figs. 3(c) and 3(d), respectively. 31 To determine the surface tension of water and whole blood, wavenumbers kexp and capillary wave phase velocities Cp_exp need to be determined. The maximum energy peak in the k-space, as shown in Fig. 3(c) for water and Fig. 3(d) for whole blood, was selected to determine the corresponding kexp and frequency coordinates to determine Cp_exp. Figures 3(c) and 3(d) illustrate an example of a wavenumber and its corresponding dispersion frequency in capillary wave on water and whole blood. The capillary wave phase velocities can be determined by a single dispersion frequency divided by its associated wavenumber. Figure 3(e) demonstrates the experimental and theoretical results of the phase velocity in the capillary wave on water and whole blood.
All the parameters that were assumed to construct the curves in Fig. 1 are assumed here. For the theoretical calculation, we assumed the values of surface tension of 73 mN/m for water 27 and 55 mN/m for whole blood. 4 The mass densities of water and whole blood are 1000 kg/m 3 and 1060 kg/m 3 , respectively. 28,29 The mass density of air is 1.2 kg/m 3 and gravity is 9.81 m/s 2 . The fluid depth d is assumed to ARTICLE scitation.org/journal/adv be 3 mm. According to Eq. (1), the theoretical phase velocities with various wavelengths from 0 mm to 20 mm were calculated and illustrated in Fig. 1. Based on Fig. 1, the theoretical capillary phase velocity is 0.2891 m/s for water at the measured wavelength of 6.1 mm (or wavenumber of 165.07 m −1 ) and 0.2262 m/s for whole blood at the measured wavelength of 8.6 mm (or wavenumber of 116.23 m −1 ). These wavelengths were identified from the energy center of the kspace, as they represent the primary wave propagating on the surface of the fluid. For the experimental calculation, the experimental phase velocity is calculated by the frequency divided by the wavenumber at the energy center of a k-space. According to Fig. 3(c), the experimental phase velocity is 0.2899 m/s for water and 0.2201 m/s for whole blood according to Fig. 3(d). The accuracy of the phase velocity is 99.72% for water and 97.3% for whole blood. This accuracy evaluation is valid with respect to the reported literature values for the surface tension and density and for the measured wavelengths.
We performed ten replicate measurements in water and 20 replicate measurements in whole blood. The reason that more measurements were made in the blood is possible inhomogeneity in the blood compared to the water (for testing the system) and to produce more reliable statistical results for the biological case. The phase velocities in experimental and theoretical calculations were 0.28 ± 0.0026 m/s and 0.28 ± 0.0011 m/s for water, respectively, and 0.22 ± 0.002 m/s and 0.22 ± 0.0008 m/s for porcine whole blood, respectively. The standard deviations in the theoretical results in Fig. 3(e) are due to slight shifting of energy center in each measured k-space realization. In this study, we observed very close agreement and high accuracy of phase velocities between experimental and theoretical results in both water and whole blood.
According to Eq. (4), the surface tension is directly proportional to C 2 p_exp and inversely proportional to k 3 exp . Figure 4 clearly exhibits that the surface tension in water is larger than that in whole blood, which matches our results that the phase velocities in water are more than those in whole blood, as shown in Fig. 3(d). Figure 4 shows that the mean values of surface tension are 71.26 ± 1.22 mN/m for water and 51.14 ± 1.34 mN/m for whole blood. A well-known work published by Vargaftik et al. reported that the surface tension of water at 20 ○ C is 72.75 mN/m. 36 The surface tension of distilled water at 22 ○ C is 72.45 mN/m. 4 Takamura et al. described that the salinity at 20 ○ C is approximately 73 dyn/cm. 27 On the other hand, the surface tension of human whole blood at 22 ○ C was reported to be approximately 55 mN/m. 4 Our results have very good agreement with previous studies.
In summary, we demonstrate that the capillary waves-based method ARF-OCT might be a promising tool to evaluate the surface tension of biological fluids. The acoustic radiation force is a useful noncontact method to generate capillary waves on the surface of fluids, which can be recorded by OCT. The capillary wave phase velocities from the theory and experiments revealed significant agreement. Human biological fluids such as mucus, cerebrospinal fluids, endocrine glands, and alveolar lining fluids contain considerable surfactants and proteins. 37 Various molecular weight surfactants control the surface tension in tissues of the human body. The surface tension behavior is an important physicochemical feature affected by various disorders. 37,38 Changes in the surface tension behavior of biological fluids are associated with certain disease and its severity, such as respiratory distress syndrome, glomerulonephritis, chronic bronchitis, and neurosyphilis. 37 The ARF-OCT would be a promising tool to evaluate the surface tension in various biological fluids and hematological diseases and a useful biomarker in the future biological applications.