Next Article in Journal
An Improved Passivity-based Control of Electrostatic MEMS Device
Next Article in Special Issue
Patterning Biological Gels for 3D Cell Culture inside Microfluidic Devices by Local Surface Modification through Laminar Flow Patterning
Previous Article in Journal
High-Fidelity Harmonic Generation in Optical Micro-Resonators Using BFGS Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Measuring Surface and Interfacial Tension In Situ in Microdripping Mode for Electrohydrodynamic Applications

by
Karim I. Budhwani
1,2,*,
Gerald M. Pekmezi
2 and
Mohamed M. Selim
2
1
CerFlux, Inc., Birmingham, AL 35205, USA
2
School of Medicine and School of Engineering, University of Alabama at Birmingham (UAB), Birmingham, AL 35294, USA
*
Author to whom correspondence should be addressed.
Micromachines 2020, 11(7), 687; https://doi.org/10.3390/mi11070687
Submission received: 18 June 2020 / Revised: 6 July 2020 / Accepted: 7 July 2020 / Published: 16 July 2020
(This article belongs to the Special Issue Microfluidics: Tissue Chips and Microphysiological Systems)

Abstract

:
Walking on water is made possible, at least for tiny insects, by molecular interaction at the interfaces of dissimilar materials. Impact of these interactions—surface tension (SFT) and, more broadly, interfacial tension (IFT)—is particularly evident at micro and nano sizescales. Thus, implications of walking on water can be significant for SFT or IFT (S/IFT)-driven nanofabrication technologies, such as electrohydrodynamic atomization (EHDA), in developing next generation biomimetic microphysiological systems (MPS) and drug delivery systems (DDS). However, current methods for estimating S/IFT, based on sessile drops or new surface formation on a ring or plate, are unsuitable for integration with EHDA assemblies used in electrospinning and electrospraying. Here, we show an in situ method for estimating S/IFT specifically devised for EHDA applications using signal processing algorithms that correlate the frequency and periodicity of liquid dispensed in EHDA microdripping mode with numerical solutions from computational fluid dynamics (CFD). Estimated S/IFT was generally in agreement with published ranges for water–air, 70% ethanol–air, chloroform–air, and chloroform–water. SFT for solutions with surfactants decreased with increasing concentrations of surfactant, but at relatively higher than published values. This was anticipated, considering that established methods measure SFT at boundaries with asymmetrically high concentrations of surfactants which lower SFT.

Graphical Abstract

1. Introduction

With one in two men and one in three women poised to be diagnosed with cancer [1], disease modeling with biomimetic microphysiological systems (MPS) and targeted drug delivery systems (DDS) is of keen interest in micro and nanoscale processing. EHDA fabrication technologies, such as electrospinning (ESp) and electrospraying (ESy), offer promising pathways in both MPS and DDS applications. ESp is well established and continues to build momentum in tissue engineering but is eclipsed by soft lithography in MPS for investigating disease progression and treatment evaluation. However, the porosity of scaffolds, which is integral to in vivo barrier and interface functions, is either entirely absent in MPS, like lab-chip systems, or introduces considerable cost, complexity, and an unrealistic uniformity in pore geometry. Nanofibrous porous scaffolds produced using ESp would be more suitable for recapitulating the in vivo tissue microenvironment in MPS, particularly for modeling organ–capillary transport, the air–liquid interface, and tumor progression [2]. Mechanisms to improve the throughput of nanofibrous scaffolds continue to gain momentum [3]. Similarly, ESy can be used to manufacture micro and nanocontainers [4,5,6] that reduce systemic toxicity and increase efficacy by delivering therapeutics only to the target site. While emulsification is the established method for manufacturing gas-filled microbubbles, used clinically as ultrasound contrast agents [7] (UCA), ESy offers considerable advantages for fabricating drug-loaded core-shell devices; specifically, in terms of higher encapsulation efficiency, near monodisperse distribution, narrow size distribution [8,9], and versatile morphology options including core-shell nanofibers [10,11]. Combining the two ideas, EHDA can also be used to produce bioactive nanofibrous scaffolds for tissue engineering applications that encapsulate growth factors or genetic material and are capable of programmed delivery of these to engineered microtissue over time [12,13]. Aside from tissue engineering and biomedical applications, EHDA is integral to a wide range of technologies, including filtration, renewable and green energy, and electrospray thrusters for space propulsion applications. New multiplexing and planar arrays [14], and mechanisms for optimizing these [15], continue to emerge to meet the increasing throughput demands of such applications. Other limiting factors, such as clean room fabrication, are being overcome with additive manufacturing capabilities [16].
However, the very principle governing Taylor cone formation and jet breakup in EHDA—interplay among S/IFT, gravity, mechanical forces, and electric fields—adds a dimension of instability and unpredictability [9,17] to the process. Variations in S/IFT from alterations in material composition or concentration, equipment, or environment, can dramatically alter output, efficiency, or both. Consider variations in temperature (energy of particles at the interface), humidity (saturation and concentration gradients), gravity, and other forces (centrifugal, centripetal, electromagnetic, mechanical) driving particles at material interfaces closer or further apart. Accounting for these in situ could assist in extemporaneous tuning of parameters to compensate for such variances. For instance, applied electric fields can be reduced in response to lower than expected S/IFT from changes in concentration of added pharmaceutical surfactants. Additionally, identification of S/IFT values asynchronously would be useful in the analysis and post-processing phase of EHDA.
It is informative to briefly review the merits of the EHDA microdripping mode as a starting point for developing computational modeling. EHDA processing regimes as a function of applied voltage (V) and flow rate (Q) are loosely illustrated in Figure 1. The Microdripping mode is generally observed for combinations of very low flow rates (Q = Qlow) and very low, even zero, applied voltage (V = Vmin) [18,19,20,21,22,23]. Due to the relatively low influence of mechanical and electrical forces in this mode, drops are dispensed almost solely due to the interplay between S/IFT and gravity. With electromechanical stresses muted, at Vmin = 0, this mode is particularly suited for in situ measurement of S/IFT. Moreover, as one would intuitively expect, particle production frequency in this mode is relatively very low—from near zero Hz to low kHz. Another interesting consequence of the virtually nonexistent influence of Vmin is that there is no cone formation at the tip of the nozzle. Instead, as the solution flows through the assembly to the nozzle tip, the drop at the tip continues to grow primarily as a function of SFT, capillary diameter, and Qlow until the weight (gravity) of the drop overcomes SFT. At this point, one part of the drop from the nozzle tip is dispensed, while the remaining solution relaxes to regain its original smaller shape at the tip of the nozzle. This mode is also characterized with near monodisperse particle sizes.
Many refinements and alternatives have been developed for Worthington’s classic Pendant Drop [24], from reducing measurement times to accounting for dynamic elements [25,26,27]. However, established S/IFT measurement methods, including du Noüy ring and Wilhelmy plate, are not suited for EHDA integration. Typically, EHDA assemblies include pumps that dispense solutions at controlled rates through precision diameter nozzles, from which the solution is spun or sprayed. This setup precludes rings or plates, which are central to the du Noüy and Wilhelmy methods. Similarly, kinetics and forces involved in spinning-drop and bubble-pressure tensiometry [28] are incompatible with EHDA, where the solution is neither in a state of perturbation nor in equilibrium. On the other hand, pendant or dispensed drops in Worthington’s pendant drop, stalagmometric, and drop volume methods, seem better suited for EHDA integration. These methods derive SFT from the shape of pendant drops, average weight of collected drops, and residual drop-volume on the nozzle tip after the drop is dispensed, respectively. Bridging elements from these with statistical analysis, digital signal processing (DSP) algorithms and computational fluid dynamics (CFD) simulations [29,30], we devised a drop kinetics (or “drop-kicks”) tensiometry platform that seamlessly integrates with EHDA configuration, as shown in Figure 2A,B. Our approach, which is notably different from the du Noüy ring and Wilhelmy plate methods, offers considerable advantages including form-factor, interface characteristics, flow kinetics that are germane to EHDA and the capacity for IFT determination, which is particularly relevant for coaxial EHDA applications. Our approach also differs from pendant drop in terms of flow kinetics (dispensed vs. sessile drops) and computation (periodicity and frequency vs. radii of surface curvature). The use of drop frequency measurements enables us to compute IFT under dynamic conditions, where inertial effects are not negligible and the equilibrium assumption of the Young–Laplace equation is not valid. This precludes the use of the Pendant Drop method, which is typically considered appropriate for drop frequencies up to 1 Hz [31].

2. Materials and Methods

EHDA setup modifications include the addition of a computer-connected charge-coupled device (CCD) camera, a light source, and a dark mat to absorb scattered light. For IFT, the lighter fluid is placed in a conical collector, with the nozzle-tip positioned just under the meniscus, while the heavier fluid is dispensed from the nozzle. Video images are captured (Figure 2C) at 200+ frames per second (fps) as fluid is dispensed in microdripping [20] mode—low flowrate (Q) and 0 kV applied voltage (V)—to reduce variables and multiphysics complexity. Nozzle diameter and flow velocity were chosen to maximize droplet frequency, while maintaining a sufficiently low Weber number; the Weber number is the ratio of inertial force to SFT. This is to ensure that SFT has a significantly higher impact, compared to fluid viscosity, on flow kinetics. Custom MATLAB® (2017a, MathWorks, Inc., Natick, MA, USA) modules were developed to derive time domain signals from the flux in the columnar sum of pixel intensities of light reflected from the dispensed drops in the captured frames. The algorithm then analyzes this signal to determine periodicity and frequency (Figure 2D).
If the signal is periodic, S/IFT is determined by cross-referencing it against periodicities and frequencies computed using the same algorithms previously run on volume-fraction output (instead of reflected light intensity) from a coupled level-set volume of fluid (VOF) model simulation.

2.1. Computational Model

The axisymmetric computational domain consists of a 2D cross section of stainless steel 14 g needle (1.600/2.108 mm inner/outer diameter). In order to improve performance, only the needle tip was included in the domain. Further to this, it was assumed that the needle tip has a smooth inner and outer surface with perfect circular inner and outer diameters, as shown in Figure 3A, with a non-adaptive mesh. Applying levels of mesh adaptation functions (Figure 3B) had no significant impact, indicating mesh-independence.
The parameters used in the simulations are summarized in Table 1. Maximum iterations per time step were set to 10 for all simulations. Output from simulations was analyzed using custom Python modules for post-processing speed, consistency, and accuracy. Please refer to the Supplementary Material section for Ansys® Fluent® (18.0, Ansys, Inc., Canonsburg, PA, USA) project runtime parameters.

2.2. Signal Processing and Spectral Analysis

The fast frame video captured using LabView® (2016, National Instruments, Austin, TX, USA) was processed using a custom MATLAB® program to confirm the periodicity of dispensed droplets and calculate droplet frequencies from the input signal comprising of a series of image frames. Briefly, the algorithm comprised of transforming 2D image frames to a 1D array of columnar intensity sums as a means to rapidly pinpoint the droplet’s position and size in each frame. The input signal, which is essentially an array of 2D image frames, is thus converted into an array of 1D arrays, which forms a convenient numerical proxy for the original fast frame video captured by the LabView® module. A visual representation of this array of 1D arrays is shown in the first panel in Figure 2D. Next, spectral density and period are computed by the algorithm using Welch’s method. Spectral density of a signal is simply the power of the signal at different frequencies. The algorithm starts by converting the input signal from the time to frequency domain to obtain the periodogram/spectrum. Compared to the standard periodogram spectrum and Bartlett’s method, Welch’s method reduces noise by compromising (reducing) frequency resolution. This is advantageous when the signal contains noise from imperfect or finite data and where noise reduction is desired. A visual representation of the Welch spectral density computed is shown in the second panel in Figure 2D. The algorithm also performs a confirmatory Fast Fourier-Transform (FFT) analysis. It is important to note that FFT is “conjugate symmetric” with a 2-sided spectrum function of both positive and negative frequencies (−Fs/2 to +Fs/2). The algorithm computes a 1-sided spectrum (with twice the amplitude) based on a 2-sided spectrum and signal length equal to the original signal. Absolute values of the complex-valued FFT are computed to extract the 1-sided magnitude and filter out phase data. A visual representation of the FFT is shown in the third panel in Figure 2D. Please refer to Supplementary Material section for LabView® and MATLAB® source files.

2.3. Experimental

Minor changes were made to the experimental setup described in earlier studies. Briefly, a high precision syringe pump was connected to a stainless steel 14 g needle (1.600/2.108 mm inner/outer diameter; purchased from ramé-hart instrument co.) using 1/16” chemical resistant autoclavable Tygon® tubing (Saint-Gobain Corporation, La Défense, Courbevoie, France). For flow rates below 500 µL/min, a BeeHive® MD-1020 (Bioanalytical Systems, Inc., Lafayette, IN, USA) syringe pump was used, while a KDS-410 (kdScientific®, Holliston, MA, USA) was used for flow rates ranging from 500 µL/min to 20 mL/min. A high-speed Teledyne Dalsa Genie® CR-GM00-H6401 camera (300 fps, Teledyne DALSA, Waterloo, Ontario, Canada) fitted with extension tubes and a 10x macro lens and connected to a computer system running custom LabView (LabView® 2016, Vision Acquisition Software® 2017, and Vision Development Module® 2016) and MATLAB (MathWorks® MATLAB® 2017a) programs, was positioned 40 mm away from the needle tip, normal to which, a focused LED light source was also positioned 40 mm away from the needle. As noted earlier, to ensure a low Weber number, parameters such as nozzle diameter, flow rate, density, and viscosity must be appropriately configured. Droplet frequencies near 1–5 Hz minimize computational resources and also control experimental variance. In addition, care must be exercised to manage conditions such as atmospheric saturation, nozzle-plane leveling, and surrounding vibration. Tubing and apparatus must be examined to ensure that no contaminants are inadvertently introduced in the system. Furthermore, pumps must be calibrated and fine-controlled to ensure smooth and accurate dispensing of solution, which in turn, ensures a near-constant flow velocity. Ethanol (70% v/v) was purchased from Ricca Chemical Company (Arlington, TX, USA). Saline (0.9% NaCl with trace HCl or NaOH for pH adjustment) was purchased from Hospira (Lake Forest, IL, USA). Molecular biology grade chloroform (with 0.75% Ethanol preservative) and non-ionic surfactant Tween® 20 (Polyoxyethylene-20-sorbitan Monolaurate) were purchased from Fisher Scientific (Fair Lawn, NJ, USA). All experiments were at least performed in triplicate. At least 3 measurements were made for all experimental observations.

3. Results

Frequencies for discrete S/IFT from CFD solutions were interpolated to obtain reference values at finer S/IFT intervals (included in the Supplementary Material). The resulting Chloroform–Air SFT and Chloroform–Water IFT curves were approximate linear functions of “drop-kicks” frequency. However, 70 v/v % Ethanol–Air and Water–Air SFT (Figure 4A) were curvilinear. Water–air SFT sported a parabolic feature between 70 and 75 dyne/cm with a virtual flatline between 72 and 72.8 dyne/cm. The flatline was resolved by replacing it with two points at 72.5 ± 0.3 dyne/cm with simple frequency averages for either side. It was also observed that water drops were dispensed slightly faster in simulation. A straight-line adjustment (labeled “Adj” in Figure 4A) of −4% for water, and +1.5% for other cases, was applied.
Figure 4B shows S/IFT obtained by cross-referencing in situ recorded “drop-kicks” frequencies with those from CFD solutions along with corresponding published values [32,33] (dashed lines).
Finally, for solutions with surface-active additives or surfactants, which are widely used in pharmaceutical applications, we expected our SFT estimates to be higher. This is because it is energetically favorable for surfactants to orientate toward the surface, which in turn, lowers the SFT at the interface. The plate/ring setup, thus, makes it so that SFT is measured at boundaries with asymmetrically high concentrations of surfactants. As shown in Figure 5, we reported a reduction in SFT with increasing concentrations of Tween® 20 aqueous solutions similar to published profiles but, as anticipated, at relatively higher values of SFT [34,35,36].

4. Discussion

The level-set method [37] is an interface-tracking method for computing two-phase flows with topologically complex interfaces, wherein the interface is captured and tracked by the level-set function. Because this function is smooth and continuous, its spatial gradients can be accurately calculated, which in turn, produces accurate estimates of interface curvature and SFT caused by the curvature. However, this method is found to be deficient in volume conservation [38]. On the other hand, VOF applied to a fixed Eulerian mesh [39] is naturally volume-conserving, as it computes and tracks volume-fraction of a particular phase in each cell rather than the interface itself. In this model, a single set of momentum equations is shared by the fluids, and the volume-fraction of each fluid in each computational cell is tracked throughout the domain by solving a continuity equation for the volume-fraction of one or more phases. The weakness of the VOF method lies in calculating spatial derivatives, since the VOF function (volume fraction of a particular phase) is discontinuous across the interface. To overcome the deficiencies of both, a coupled level-set VOF [40] approach was adopted. Setting interface Γ as zero level-set, the level-set function φ is defined as a signed distance to interface φ {x,t} and can be expressed as:
φ { x , t } = { + | d | if   x the primary phase 0 if   x Γ | d | if   x the secondary phase
where d is distance from the interface. Similarly, evolution of the level-set function is:
t ( φ ) + · ( φ v ) = 0
where v is the flow velocity.
Since SFT acts to minimize surface area of the interface, it gives rise to effects such as pressure discontinuity at the interface and capillary effects at adhesive walls. The continuum surface force (CSF) model [41,42] results in a source-term in the momentum equation from the addition of SFT to the VOF calculation.
t ( ρ v ) + · ( ρ v v ) = p + · τ + ρ g F s f
where F s f is the force from SFT effects, ρ is volume-averaged-density, p is the pressure, and τ is the deviatoric stress tensor. For two phases, the force from SFT effects becomes:
F s f = σ i j ρ k i α i 1 2 ( ρ i + ρ j )
where σ i j is the coefficient of SFT between primary phase i and phase j, ki is local mean interface curvature for i, and α i is volume-fraction of i. Equation (4) shows that SFT source-term for a cell is proportional to the average density in the cell. For two-phase systems, volume-fraction-averaged density becomes:
ρ = a j ρ j + ( 1 a j ) ρ i
Although, it is established that temperature plays an important role in S/IFT, temperature was explicitly measured and controlled on the experimental side. Given the temperature ranges expected in a controlled laboratory setting, on the computational side, temperature was an implicit parameter with negligible effect on viscosity and density of fluid. Numerical solutions were obtained for these governing equations over the domain shown in Figure 3A. Applying levels of mesh adaptation functions (Figure 3B) had no significant impact (Figure 3C), indicating mesh-independence. For preliminary evaluation of this approach and its applicability to a broad range of fluids, simulations and tests were run for: (i) Chloroform (CHCl3), an organic solvent with high density but low dynamic viscosity; (ii) 70 v/v % ethanol (etOH), with low density but high viscosity; (iii) water.
The “drop-kicks” frequency and periodicity of signals, derived from both CFD models and CCD videos from experiments, were determined by first computing an average drop rate. This was used as bracketing input for DSP algorithms, including Welch’s spectral density (WSD), to determine the fundamental frequency, harmonics, and power at those frequencies (Figure 2D). Periodicity is an important indicator of the requisite uniformity of volume dispensed in each drop. Solutions with aperiodic (no uniformity) or quasiperiodic signals, where uniform volume is dispensed over a series of drops, are excluded from the scope of this work. Runtime parameters and scripts encoding algorithms for analyzing signals from both CFD and CCD are included under the Supplementary Material.

5. Conclusions and Outlook

Tissue engineering applications such as MPS—for modeling disease progression and evaluating efficacies of new treatments—can be enhanced by replacing soft lithography-fabricated membranes with biomimetic nanofibrous scaffolds produced using ESp. With improved configuration control of ESp, MPS parameters could be refined to enable control and capabilities which are simply not possible in animal models—capabilities such as isolating, tweaking, and observing the impact of specific variables in a humanized model of target physiological and pathophysiological conditions. This could substantially reduce the number of animals used in preclinical studies. Similarly, EHDA can make a marked improvement in developing targeted and on-demand DDS. However, despite its demonstrated superiority over competing methods, adoption of EHDA in these and other nanomedicine applications remains dwarfed due to its weakness in reproducibility and scale.
In the words of Melcher and Taylor regarding EHDA research, from five decades ago, “the center of attention in almost any discussion is the lack of reproducibility in experiments and the inadequacies of theoretical models… Yet the foundations of fluid mechanics are formed from work that relates carefully designed experiments to analytical models”. The note regarding lack of reproducibility and of theoretical models rings true even today. Addressing these challenges is not trivial, given the number of material, process, and environmental variables, the number and complex geometries of interfaces—which vary dramatically across facilities and applications further exacerbating matters—in EHDA assemblies, surface-active components, the complexity of interactions underlying S/IFT, and the ensuing vagaries. Ironically, the very principle—interplay among S/IFT, gravity, mechanical forces, and electric fields—bestowing upon EHDA its characteristic advantages also leaves it vulnerable to the Achilles heel of instability and unpredictability. While gravity can be assumed constant and applied voltage and mechanical force can be equipment controlled, variations in S/IFT can dramatically diminish stability, alter resulting properties, or both. Our approach for estimating S/IFT relates “carefully designed experiments” to numerical and DSP models—while making minimal changes to the assembly—to innately account for conditions specific to the EHDA site, materials, and equipment.
In comparing the du Noüy ring with their preferred Wilhelmy plate, Krüss GmbH noted that “rings are for fingers — plates are for surface tension”, which we would update to “rings are for fingers, plates are for dinner — drop kinetics are for surface tension”. Dry humor notwithstanding, as evidenced in Moore’s Law [43], nanotechnology has led to the exponential growth of in silico capabilities; better hardware supporting increasingly demanding software. It is thus, appropriate to complete this loop by applying advances in computational capacity, including computational modeling and DSP, toward furthering nanotechnology; specifically, nanofabrication technologies such as EHDA in developing new theranostic devices and biomimetic microphysiological systems that enhance both the quantity and quality of life of people around the world. We hope that ours is but one step toward building that feedback loop.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-666X/11/7/687/s1.

Author Contributions

Conceptualization, K.I.B., G.M.P. and M.M.S.; Data curation, K.I.B., G.M.P. and M.M.S.; Formal analysis, K.I.B., G.M.P. and M.M.S.; Funding acquisition, K.I.B.; Investigation, K.I.B., G.M.P. and M.M.S.; Methodology, K.I.B., G.M.P. and M.M.S.; Resources, K.I.B.; Software, K.I.B., G.M.P. and M.M.S.; Writing—original draft, K.I.B., G.M.P. and M.M.S.; Writing—review and editing, K.I.B. and G.M.P. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by grants from elixir international and the National Science Foundation (NSF EPS 1158862).

Acknowledgments

We thank J.M. Rogers, H. Zhang, D.J. Buchsbaum, M. Bolding, M.N. Saleh, V. Thomas, D.L. Littlefield, P.G. Oliver, H. Budhwani, C.A. Monroe, J.C. Sellers, M. Li, P. Sethu, and G. Walcott for their insight and expertise and the O’Neal Comprehensive Cancer Center (OCCC) for use of equipment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sweeney, S.M. AACR Cancer Progress Report 2015. Am. Assoc. Cancer Res. 2014, 126, 8–22. [Google Scholar]
  2. Budhwani, K.I.; Oliver, P.G.; Buchsbaum, D.J.; Thomas, V. Novel Biomimetic Microphysiological Systems for Tissue Regeneration and Disease Modeling. In Advances in Experimental Medicine and Biology; Springer Science and Business Media LLC: Berlin/Heidelberg, Germany, 2018; Volume 1077, pp. 87–113. [Google Scholar]
  3. Srivastava, Y.; Marquez, M.; Thorsen, T. Multijet Electrospinning of Conducting Nanofibers from Microfluidic Manifolds. J. Appl. Polym. Sci. 2007, 106, 3171–3178. [Google Scholar] [CrossRef]
  4. Niu, C.; Wang, Z.; Lu, G.; Krupka, T.M.; Sun, Y.; You, Y.; Song, W.; Ran, H.; Li, P.; Zheng, Y. Doxorubicin Loaded Superparamagnetic PLGA-Iron Oxide Multifunctional Microbubbles for Dual-Mode US/MR Imaging and Therapy of Metastasis in Lymph Nodes. Biomaterials 2013, 34, 2307–2317. [Google Scholar] [CrossRef]
  5. Mura, S.; Nicolas, J.; Couvreur, P. Stimuli-Responsive Nanocarriers for Drug Delivery. Nat. Mater. 2013, 12, 991–1003. [Google Scholar] [CrossRef] [PubMed]
  6. Budhwani, K.I.; Dettmann, M.A.; Saleh, M.N.; Thomas, V. Nano and Microbubble Systems for On-Demand Cancer Drug Delivery. Curr. Nanosci. 2018, 14, 33–41. [Google Scholar] [CrossRef]
  7. Feinstein, S.B. The Evolution of Contrast Ultrasound. J. Am. Coll. Cardiol. 2016, 67, 2516–2518. [Google Scholar] [CrossRef] [PubMed]
  8. Xie, J.; Jiang, J.; Davoodi, P.; Srinivasan, M.P.; Wang, C.H. Electrohydrodynamic Atomization: A Two-Decade Effort to Produce and Process Micro-/Nanoparticulate Materials. Chem. Eng. Sci. 2015, 125, 32–57. [Google Scholar] [CrossRef] [Green Version]
  9. Bock, N.; Woodruff, M.A.; Hutmacher, D.W.; Dargaville, T.R. Electrospraying, a Reproducible Method for Production of Polymeric Microspheres for Biomedical Applications. Polymers 2011, 3, 131–149. [Google Scholar] [CrossRef] [Green Version]
  10. Sill, T.J.; Von Recum, H.A. Electrospinning: Applications in Drug Delivery and Tissue Engineering. Biomaterials 2008, 29, 1989–2006. [Google Scholar] [CrossRef]
  11. Budhwani, K.I.; Thomas, V.; Sethu, P. Lab-on-a-Brane: Nanofibrous Polymer Membranes to Recreate Organ–Capillary Interfaces. J. Micro. Micro. 2016, 26, 035013. [Google Scholar] [CrossRef]
  12. Ji, W.; Sun, Y.; Yang, F.; Van Den Beucken, J.J.J.P.; Fan, M.; Chen, Z.; Jansen, J.A. Bioactive Electrospun Scaffolds Delivering Growth Factors and Genes for Tissue Engineering Applications. Pharm. Res. 2011, 28, 1259–1272. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Liu, C.; Wang, C.; Zhao, Q.; Li, X.; Xu, F.; Yao, X.; Wang, M. Incorporation and Release of Dual Growth Factors for Nerve Tissue Engineering Using Nanofibrous Bicomponent Scaffolds. Biomed. Mater. 2018, 13, 044107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Deng, W.; Klemic, J.F.; Li, X.; Reed, M.A.; Gomez, A. Increase of Electrospray Throughput Using Multiplexed Microfabricated Sources for the Scalable Generation of Monodisperse Droplets. J. Aerosol Sci. 2006, 37, 696–714. [Google Scholar] [CrossRef]
  15. Krpoun, R.; Smith, K.L.; Stark, J.P.W.; Shea, H.R. Tailoring the Hydraulic Impedance of Out-of-Plane Micromachined Electrospray Sources with Integrated Electrodes. Appl. Phys. Lett. 2009, 94, 5–7. [Google Scholar] [CrossRef]
  16. García-López, E.; Olvera-Trejo, D.; Velásquez-García, L.F. 3D Printed Multiplexed Electrospinning Sources for Large-Scale Production of Aligned Nanofiber Mats with Small Diameter Spread. Nanotechnology 2017, 28, 425302. [Google Scholar] [CrossRef]
  17. Morad, M.R.; Rajabi, A.; Razavi, M.; Sereshkeh, S.R.P. A Very Stable High Throughput Taylor Cone-Jet in Electrohydrodynamics. Sci. Rep. 2016, 6, 38509. [Google Scholar] [CrossRef] [PubMed]
  18. Pareta, R.; Edirisinghe, M.J. A Novel Method for the Preparation of Starch Films and Coatings. Carbohydr. Polym. 2006, 63, 425–431. [Google Scholar] [CrossRef]
  19. Luo, C.J.; Loh, S.; Stride, E.; Edirisinghe, M. Electrospraying and Electrospinning of Chocolate Suspensions. Food Bioprocess Technol. 2012, 5, 2285–2300. [Google Scholar] [CrossRef]
  20. Castillo-Orozco, E.; Kar, A.; Kumar, R. Electrospray Mode Transition of Microdroplets with Semiconductor Nanoparticle Suspension /639/166/988 /639/925/927/351 /123 Article. Sci. Rep. 2017, 7, 1–10. [Google Scholar] [CrossRef]
  21. Jaworek, A.; Krupa, A. Main Modes of Electrohydrodynamic Spraying of Liquids. In Proceedings of the Third International Conference on Multiphase Flow, ICMF’98, Lyon, France, 8–12 June 1998; pp. 8–12. [Google Scholar]
  22. Jaworek, A.; Sobczyk, A.T. Electrospraying Route to Nanotechnology: An Overview. J. Electrostat. 2008, 66, 197–219. [Google Scholar] [CrossRef]
  23. Nemes, P.; Marginean, I.; Vertes, A. Spraying Mode Effect on Droplet Formation and Ion Chemistry in Electrosprays. Anal. Chem. 2007, 79, 3105–3116. [Google Scholar] [CrossRef]
  24. Worthington, A.M. On Pendent Drops. Proc. R. Soc. Lond. 1881, 32, 362–377. [Google Scholar]
  25. Gunde, R.; Kumar, A.; Lehnert-Batar, S.; Mader, R.; Windhab, E.J. Measurement of the Surface and Interfacial Tension from Maximum Volume of a Pendant Drop. J. Colloid. Interf. Sci. 2001, 244, 113–122. [Google Scholar] [CrossRef]
  26. Zuo, Y.Y.; Ding, M.; Bateni, A.; Hoorfar, M.; Neumann, A.W. Improvement of Interfacial Tension Measurement Using a Captive Bubble in Conjunction with Axisymmetric Drop Shape Analysis (ADSA). Colloids Surf. A Physicochem. Eng. Asp. 2004, 250, 233–246. [Google Scholar] [CrossRef]
  27. Hyde, A.; Phan, C.; Ingram, G. Determining Liquid-Liquid Interfacial Tension from a Submerged Meniscus. Colloids Surf. A Physicochem. Eng. Asp. 2014, 459, 267–273. [Google Scholar] [CrossRef]
  28. Kovalchuk, N.M.; Roumpea, E.; Nowak, E.; Chinaud, M.; Angeli, P.; Simmons, M.J.H. Effect of Surfactant on Emulsification in Microchannels. Chem. Eng. Sci. 2018, 176, 139–152. [Google Scholar] [CrossRef]
  29. Batchelor, G.K. An Introduction to Fluid Dynamics; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  30. Boyd, S.; Vandenberghe, L. Convex Optimization; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  31. Berry, J.D.; Neeson, M.J.; Dagastine, R.R.; Chan, D.Y.C.; Tabor, R.F. Measurement of Surface and Interfacial Tension Using Pendant Drop Tensiometry. J. Colloid Interface Sci. 2015, 454, 226–237. [Google Scholar] [CrossRef]
  32. Kadlec, P.H.S. Properties of Ethanol and Ethanol-Water Solutions-Tables and Equations. Sugar Ind. 2010, 135, 607–613. [Google Scholar]
  33. Bahramian, A.; Danesh, A. Prediction of Liquid-Liquid Interfacial Tension in Multi-Component Systems. Fluid Phase Equilib. 2004, 221, 197–205. [Google Scholar] [CrossRef]
  34. Patist, A.; Bhagwat, S.S.; Penfield, K.W.; Aikens, P.; Shah, D.O. On the Measurement of Critical Micelle Concentrations of Pure and Technical-Grade Nonionic Surfactants. J. Surfactants Deterg. 2000, 3, 53–58. [Google Scholar] [CrossRef]
  35. Rehman, N.; Ullah, H.; Alam, S.; Jan, A.K.; Khan, S.W.; Tariq, M. Surface and Thermodynamic Study of Micellization of Non Ionic Surfactant / Diblock Copolymer System as Revealed by Surface Tension and Conductivity. J. Mater. Environ. Sci. 2017, 8, 1161–1167. [Google Scholar]
  36. Bąk, A.; Podgórska, W. Interfacial and Surface Tensions of Toluene/Water and Air/Water Systems with Nonionic Surfactants Tween 20 and Tween 80. Colloids Surfaces A Physicochem. Eng. Asp. 2016, 504, 414–425. [Google Scholar] [CrossRef]
  37. Osher, S.; Sethian, J.A. Fronts Propagating with Curvature-Dependent Speed: Algorithms Based on Hamilton-Jacobi Formulations. J. Comput. Phys. 1988, 79, 12–49. [Google Scholar] [CrossRef] [Green Version]
  38. Olsson, E.; Kreiss, G.; Zahedi, S. A Conservative Level Set Method for Two Phase Flow II. J. Comput. Phys. 2007, 225, 785–807. [Google Scholar] [CrossRef]
  39. Hirt, C.W.; Nichols, B.D. Volume of Fluid Methods for the Dynamics of Free Boundaries. J. Comput. Phys. 1981, 39, 201–225. [Google Scholar] [CrossRef]
  40. Cao, Z.; Sun, D.; Wei, J.; Yu, B. A Coupled Volume-of-Fluid and Level Set Method Based on Multi-Dimensional Advection for Unstructured Triangular Meshes. Chem. Eng. Sci. 2018, 176, 560–579. [Google Scholar] [CrossRef]
  41. Brackbill, J.U.; Kothe, D.B.; Zemach, C. A Continuum Method for Modeling Surface Tension. J. Comput. Phys. 1992, 100, 335–354. [Google Scholar] [CrossRef]
  42. Shaikh, J.; Sharma, A.; Bhardwaj, R. On Comparison of the Sharp-Interface and Diffuse-Interface Level Set Methods for 2D Capillary or/and Gravity Induced Flows. Chem. Eng. Sci. 2018, 176, 77–95. [Google Scholar] [CrossRef]
  43. Waldrop, M.M. The Chips Are down for Moore’s Law. Nature 2016, 530, 144–147. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Approximate regimes of electrohydrodynamic atomization (EHDA) modes shown as a function of applied voltage (V) and flow rate (Q).
Figure 1. Approximate regimes of electrohydrodynamic atomization (EHDA) modes shown as a function of applied voltage (V) and flow rate (Q).
Micromachines 11 00687 g001
Figure 2. Schematic illustration of the modified EHDA setup for in situ (A) surface tension (SFT) and (B) interfacial tension (IFT) measurements. In both cases, a light-emitting diode (LED) light source is pointed slightly below the tip of the EHDA nozzle, while a charge-coupled device (CCD) camera captures images at a minimum 200 frames per second (fps). (C) An image of a dispensed drop captured by the camera. (D) Three panels showing the signal derived from captured image sequence (video) in the time domain, where dS/dt is the derivative of signal intensity with respect to time, which highlights peak intensities at drop locations, Welch’s power spectral density (WSD) and fast Fourier-transform (FFT) in the frequency domain. Both WSD and FFT indicate that the signal is not only periodic, given that almost all the power is in the fundamental frequency and its multiples (harmonics), but that the fundamental frequency is 2.13 Hz, which is close to the average drop rate.
Figure 2. Schematic illustration of the modified EHDA setup for in situ (A) surface tension (SFT) and (B) interfacial tension (IFT) measurements. In both cases, a light-emitting diode (LED) light source is pointed slightly below the tip of the EHDA nozzle, while a charge-coupled device (CCD) camera captures images at a minimum 200 frames per second (fps). (C) An image of a dispensed drop captured by the camera. (D) Three panels showing the signal derived from captured image sequence (video) in the time domain, where dS/dt is the derivative of signal intensity with respect to time, which highlights peak intensities at drop locations, Welch’s power spectral density (WSD) and fast Fourier-transform (FFT) in the frequency domain. Both WSD and FFT indicate that the signal is not only periodic, given that almost all the power is in the fundamental frequency and its multiples (harmonics), but that the fundamental frequency is 2.13 Hz, which is close to the average drop rate.
Micromachines 11 00687 g002
Figure 3. Computational fluid dynamics (CFD) domain, adaptive mesh, and drop kinetics simulation. (A) Domain showing tip of the nozzle and non-adaptive mesh. (B) Mesh adaptation function. (C) Volume fractions from a detaching drop in two-phase simulations of dispensed drops.
Figure 3. Computational fluid dynamics (CFD) domain, adaptive mesh, and drop kinetics simulation. (A) Domain showing tip of the nozzle and non-adaptive mesh. (B) Mesh adaptation function. (C) Volume fractions from a detaching drop in two-phase simulations of dispensed drops.
Micromachines 11 00687 g003
Figure 4. Estimating SFT and IFT. (A) Frequencies from CFD numerical solutions for discrete S/IFT values were interpolated, labeled “CFD,” and adjusted, “Adj”, −4% for Water–Air and +1.5% for other interfaces. Bold dashed lines indicate published [32,33] range of S/IFT for respective interfaces. (B) S/IFT derived by cross-referencing in situ recorded “drop-kicks” frequencies against corresponding “Adj” frequencies from CFD numerical solutions. Aside from a few outliers, resulting S/IFT were in agreement with published ranges, particularly for Chloroform–Water IFT. Flowrate (Q) for CHCl3–Air SFT was 3 mL/min, all others were at 5 mL/min.
Figure 4. Estimating SFT and IFT. (A) Frequencies from CFD numerical solutions for discrete S/IFT values were interpolated, labeled “CFD,” and adjusted, “Adj”, −4% for Water–Air and +1.5% for other interfaces. Bold dashed lines indicate published [32,33] range of S/IFT for respective interfaces. (B) S/IFT derived by cross-referencing in situ recorded “drop-kicks” frequencies against corresponding “Adj” frequencies from CFD numerical solutions. Aside from a few outliers, resulting S/IFT were in agreement with published ranges, particularly for Chloroform–Water IFT. Flowrate (Q) for CHCl3–Air SFT was 3 mL/min, all others were at 5 mL/min.
Micromachines 11 00687 g004
Figure 5. Effect of surfactants and additives on SFT. (A) Frequencies from CFD numerical solutions for a wider range of discrete Water–Air SFT values interpolated, labeled “CFD,” and adjusted, “Adj,” −4%. (B) “Drop-kicks” frequencies cross-referenced against Water–Air “Adj” frequencies to derive SFT for Water–Air and various micromolar concentrations of aqueous Tween® 20 solutions. Bold dashed lines indicate published ranges [34,35,36] of SFT. SFT decreased with increasing concentrations of Tween® 20 but at relatively higher than published values. This was anticipated considering that established methods measure SFT at boundaries with asymmetrically high concentrations of surfactants which lower SFT.
Figure 5. Effect of surfactants and additives on SFT. (A) Frequencies from CFD numerical solutions for a wider range of discrete Water–Air SFT values interpolated, labeled “CFD,” and adjusted, “Adj,” −4%. (B) “Drop-kicks” frequencies cross-referenced against Water–Air “Adj” frequencies to derive SFT for Water–Air and various micromolar concentrations of aqueous Tween® 20 solutions. Bold dashed lines indicate published ranges [34,35,36] of SFT. SFT decreased with increasing concentrations of Tween® 20 but at relatively higher than published values. This was anticipated considering that established methods measure SFT at boundaries with asymmetrically high concentrations of surfactants which lower SFT.
Micromachines 11 00687 g005
Table 1. Summary of simulation parameters.
Table 1. Summary of simulation parameters.
SimulationMaterialDensity (kg/m3)Viscosity (kg/m·s)Velocity (m/s)Time Step (s)
Water-AirWater998.20.0010030.0835 × 10−6
Air1.2251.7894 × 10−5
Ethanol-Air70% Ethanol8800.00250.04151 × 10−6
Air1.2251.7894 × 10−5
Chloroform-AirChloroform14900.0005630.04155 × 10−6
Air1.2251.7894 × 10−5
Chloroform-WaterChloroform14900.0005630.04151 × 10−6
Water998.20.001003
Air1.2251.7894 × 10−5

Share and Cite

MDPI and ACS Style

Budhwani, K.I.; Pekmezi, G.M.; Selim, M.M. Measuring Surface and Interfacial Tension In Situ in Microdripping Mode for Electrohydrodynamic Applications. Micromachines 2020, 11, 687. https://doi.org/10.3390/mi11070687

AMA Style

Budhwani KI, Pekmezi GM, Selim MM. Measuring Surface and Interfacial Tension In Situ in Microdripping Mode for Electrohydrodynamic Applications. Micromachines. 2020; 11(7):687. https://doi.org/10.3390/mi11070687

Chicago/Turabian Style

Budhwani, Karim I., Gerald M. Pekmezi, and Mohamed M. Selim. 2020. "Measuring Surface and Interfacial Tension In Situ in Microdripping Mode for Electrohydrodynamic Applications" Micromachines 11, no. 7: 687. https://doi.org/10.3390/mi11070687

APA Style

Budhwani, K. I., Pekmezi, G. M., & Selim, M. M. (2020). Measuring Surface and Interfacial Tension In Situ in Microdripping Mode for Electrohydrodynamic Applications. Micromachines, 11(7), 687. https://doi.org/10.3390/mi11070687

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop