Abstract

The ongoing development of modern electronic systems leads to smaller, more powerful devices that are expected to operate in complex environments. Due to this, advanced thermal management technologies are required to meet the growing demand, especially in space where two-phase thermal systems are limited by the absence of gravity. Electrohydrodynamic (EHD) and dielectrophoretic (DEP) forces can be used to sustain stable liquid film flow boiling in the absence of gravity, which is otherwise impractical due to the lack of a required buoyancy force to initiate bubble departure. EHD is a phenomenon that is represented by the interaction between electric fields and fluid flow. The DEP force is characterized by its ability to act on liquid/vapor interfaces due to a high gradient of electrical permittivity. This study investigates the heat transfer characteristics of EHD conduction pumping driven liquid film flow boiling coupled with DEP vapor extraction during a microgravity parabolic flight and on the ground. The results of this study show that EHD and DEP raise the critical heat flux, lower heater surface temperature, and successfully sustain boiling in both microgravity and on the ground with low power consumption. Additionally, the heat transfer data captured in terrestrial, microgravity, and 1.8 g conditions compare well, indicating that combining these mechanisms can provide thermal enhancement independent of gravity. This study provides fundamental understanding of electrically driven liquid film flow boiling in the presence of phase change, paving the way toward developing next-generation heat transport devices for space and terrestrial applications.

Introduction

Electrohydrodynamics (EHDs) characterizes the interaction between electric fields and fluid flow. This phenomenon can induce fluid flows, create mixing, and enhance heat transfer. As a pumping mechanism, EHD has shown promise due to its low power consumption, rapid response, geometric simplicity, low weight, and heat transport capabilities. It also involves no mechanical parts, avoiding issues with vibration and fatigue, which is especially advantageous for space applications.

The EHD body force that describes this phenomenon is described in the following [1]:
(1)

where E represents the electric field, ρe represents the net charge density, ρ is the fluid density, and ε is the fluid electrical permittivity. The first term in this expression represents the Coulomb force density, which exists when free charges are present in the fluid under an electric field. The second term in the expression represents the dielectrophoretic (DEP) force density, which exists, for example, when there is a liquid/vapor interface or a large temperature gradient corresponding to a gradient in the electrical permittivity of the medium. The third term in the expression is the electrostriction term, which is negligible when the working fluid is an incompressible liquid film. Both the Coulomb force density and DEP force density are the driving mechanisms in this study.

Electrohydrodynamic pumping can be separated into three main forms: ion drag pumping, induction pumping, and conduction pumping [2]. This study utilizes EHD conduction pumping, which uses ionic dissociation to generate free charges in the fluid [3]. In a dielectric fluid, the rate of dissociation and recombination of these free charges exists at an equilibrium. If acted upon by a strong enough electric field, the rate of dissociation exceeds the rate of recombination. Via this mechanism, nonequilibrium charge layers, comprised of these free charges, form near the electrodes, which are referred to as heterocharge layers. As charges are attracted to these layers, bulk fluid motion occurs. If the electrodes are designed with an asymmetrical geometry, a net Coulomb force is induced in the desired direction, and therefore, a fluid net flow can be produced. Figure 1 [4] shows these heterocharge layers along with the induced net fluid flow direction. The smaller arrows in Fig. 1 near the heterocharge layers represent the attraction of charges. Further analysis into the physical fundamentals of EHD conduction phenomenon can be found from the efforts of Vázquez et al. [5], Yazdani et al. [6], and Talmor et al. [7].

Fig. 1
Heterocharge layers formed over electrodes via EHD conduction [4]
Fig. 1
Heterocharge layers formed over electrodes via EHD conduction [4]
Close modal

EHD conduction pumping has been explored in macro-and meso scales as well as microscale [8]. Additionally, pumps in both rigid and flexible configurations have been explored, with the inclusion of particle image velocimetry and particle tracking [912]. The results of these studies have shown that EHD conduction pumping, in both rigid and flexible configurations, can generate strong net flows upwards of 20 cm/s of dielectric fluids in the presence of electric fields.

The use of DEP forces as a method to enhance boiling heat transfer has been studied previously through the efforts of Ogata and Yabe [13], Yagoobi et al. [14,15], and Darabi et al. [16]. Additional studies on enhancing nucleate boiling via a nonuniform electrical field were conducted by Kano and Takahashi [17], including a study using HFE 7100 [18]. Their work consisted of a microscale electrode design where heights from the surface of the heater ranged from 100 to 300 μm. It was found that the maximum critical heat flux (CHF) value was approximately three times higher than the observed CHF under the influence of no electric field. Recently, Garivalis et al. [19] explored the enhancement of critical heat flux in microgravity using microstructured surfaces and applied electric fields. The CHF in this study was found to increase by over 50% when utilizing the combination of electric fields with enhanced surfaces. Patel et al. [20] explored the impact of gravity on EHD conduction pumping and its effect on heat transfer enhancement. When activated, EHD pumping lowered the heater surface temperature at higher heat fluxes in microgravity. However, no DEP mechanism was involved in this study.

The coupling of the EHD conduction pumping mechanism and the DEP extraction force on boiling heat transfer enhancement in the absence of gravity is explored for the first time in this work. It is important to note that liquid pool boiling has been successfully explored in the absence of gravity by Stephan et al. [21] and Grassi et al. [22], the latter of which used an electric field to reduce bubble diameter and enhance heat transfer in pool boiling conditions. Despite this, it is still considered impractical, especially at higher heat fluxes. This has been confirmed by studies from Henry and Kim [23] and Dhir et al. [24], which showed that under microgravity conditions, pool boiling is dominated by a large primary vapor bubble formation, causing dryout. Additionally, previous studies [2124] focused on pool boiling, and this study explores liquid film flow boiling. The working fluid in this study is Novec 7100 (C4F9OCH3). This fluid, also known as HFE 7100, was selected for its favorable dielectric properties, and its thermophysical and electric properties are tabulated in Table 1 [25,26].

Table 1

Thermophysical and electrical properties of Novec 7100 [23,24]

Liquid Range at 1 atm−135 °C to 61 °C
Liquid Density (kg/m3)2.2690T[K]+2157.737
Vapor Pressure (Pa)exp(22.4153641.9T[K])
Kinematic Viscosity (mm2/s)85.44*exp(1859.59T[K]0.009152T[K]+332097T[K]2)
Relative Permittivity0.026T[°C]+7.3969
Electric Conductivity (S/m)2.2907*109+1.20114*108*exp(2(T[°C]18.570620.295)2)
Liquid Range at 1 atm−135 °C to 61 °C
Liquid Density (kg/m3)2.2690T[K]+2157.737
Vapor Pressure (Pa)exp(22.4153641.9T[K])
Kinematic Viscosity (mm2/s)85.44*exp(1859.59T[K]0.009152T[K]+332097T[K]2)
Relative Permittivity0.026T[°C]+7.3969
Electric Conductivity (S/m)2.2907*109+1.20114*108*exp(2(T[°C]18.570620.295)2)

This work is an extension of a previous study carried out by Patel and Yagoobi [27], which explored the boiling heat transfer enhancement of EHD conduction and DEP mechanism together in a terrestrial environment only using HCFC-123 as the working fluid. The main difference between the previous study and the current study is that the sole and combined DEP and EHD mechanisms and their effects on liquid film flow boiling are investigated in different gravity conditions. Additional differences between these studies include the designs of the DEP electrode, EHD conduction pump, and heater, as well as the working fluid. Similar to the previous experiment, a radial EHD conduction pump drives dielectric fluid inward toward the center as a way to electrowet a heater. This heater generates vapor bubbles at the center of the experiment. A DEP electrode is installed 1.6 mm above the heater surface to provide a local nonuniform electric field, which enhances vapor bubble extraction during liquid film flow boiling via the DEP force. As the vapor is extracted from the heater surface by the DEP electrode, it travels to the periphery of the experimental chamber and condenses back onto the EHD conduction pump, where it travels as liquid back to the heater.

The boiling phenomenon of this study is explored in terrestrial and microgravity flight environments, and the results are compared. Combining DEP and EHD forces, this study explores enhanced liquid film flow boiling as a way of removing heat in space. DEP and EHD conduction mechanisms are enabling technologies for the next-generation cooling of high-powered electronics and thermal management in space.

Experimental Design

The experiment consists of several main components including a radial EHD conduction pump, a DEP electrode, and a resistance heater. The design of the radial EHD conduction pump used in this study is shown in Fig. 2. The characteristic length, denoted as d, represents the spacing between the high voltage and ground electrode of a selected pair. Yazdani and Yagoobi [6] determined the electrode sizing ratio and spacing between the neighboring pairs to be the ideal dimensions for EHD conduction via numerical optimization. The main difference between the previous EHD conduction pump design [20,27] and the one used in this study is the material of the substrate and electrodes. Here, the substrate material of the pump is silicon, and the electrodes are platinum. The DEP electrode design, shown in Fig. 3, is made from stainless steel and exhibits grate-like spacing, similar to the DEP electrode from [27], to produce a nonuniform diverging electric field. The design imposes a DEP force on the bubbles that form on the heater surface and allows the bubbles to pass through the spaces within the electrode. In this configuration, the heater surface acts as the ground electrode to the DEP electrode. Constructed via platinum vapor deposition, the 2.25 cm2 heater is located in the center of the chamber on an Ultem 1000 pedestal, which thermally isolates it. This heater is different from the design used in Patel [27], which used a circular copper heater encased in a Delrin housing. Finally, a ThermoCube recirculating chiller utilizes chilled water to provide a constant temperature of 15 °C below the EHD conduction pump disk. The EHD disk acts as the condenser in the system, allowing for the vapor generated by boiling to condense on the disk to subsequently be pumped back to the heater. A diagram showing top-down and cross-sectional views of this system is shown in Fig. 4.

Fig. 2
Schematic showing EHD conduction pump dimensions and electrode spacing
Fig. 2
Schematic showing EHD conduction pump dimensions and electrode spacing
Close modal
Fig. 3
Schematic showing DEP electrode dimensions (units are millimeters)
Fig. 3
Schematic showing DEP electrode dimensions (units are millimeters)
Close modal
Fig. 4
Diagram of EHD conduction pump, DEP electrode, and heater assembly
Fig. 4
Diagram of EHD conduction pump, DEP electrode, and heater assembly
Close modal

The experimental chamber utilizes an aluminum circular housing with NPT and Yor-Lok fittings for filling, vapor pressure measurement, electrical connections, and supplying chilled water. Additionally, two Stanford Research Systems PS350 high voltage power supplies are used to power the conduction pump and the DEP electrode, and an Agilent 6634B DC power supply is used to control power the heater. The complete flight rack with the experimental chamber and all associated equipment is shown in Fig. 5.

Fig. 5
Photo of fully assembled experiment aboard flight rack (top view)
Fig. 5
Photo of fully assembled experiment aboard flight rack (top view)
Close modal

An Agilent 34970A data acquisition system is used to collect data during the experiment. Applied voltage is monitored along with the corresponding current for the EHD pump, DEP electrode, and heater. Temperatures in the experiment are monitored in multiple locations. The heater has a platinum resistance temperature detector (RTD) embedded on its surface, as well as a thermistor for redundancy. The surface temperatures of the EHD conduction pump disk are also monitored with platinum resistance thermometers (PRTs) at five locations labeled in Fig. 2. These PRTs are covered in epoxy to prevent electrical interference due to proximity to the EHD electrodes. Another RTD probe is used to measure the vapor temperature. The ambient temperature in the aircraft and in the lab, depending on the gravitational scenario, as well as the chamber outer wall and experiment housing base plate temperatures, are recorded. The base plate houses the majority of the sensors and the experiment itself. In order to capture the gravitational aspects of the experiment, an accelerometer is installed on the baseplate. LabVIEW is used to control and monitor the experiment and has an automatic experiment shut down procedure, should the heater temperature exceed 70 °C. Visualization data at 1080p/60fps were obtained with a camcorder mounted directly above the experiment overlooking the glass window installed on the top of the test chamber.

Prior to charging, the chamber is vacuumed and then filled with a liquid film of HFE 7100 up to a height of 2 mm. The chamber is then degassed using a vacuum pump to the saturation pressure that corresponds to the operating temperature of the experiment. This was done to remove any noncondensable gasses that may enter when charging with fluid. The system loop is as follows; the working fluid is boiled at the heater surface. The DEP electrode extracts the vapor from the heater surface. The vapor condenses on the chilled silicon EHD conduction pump disk. The EHD conduction pump then drives the condensed liquid film back to the heater. In a typical thin-film boiling system, buoyancy/gravity drives the wetting of the heater and removal of vapor bubbles. The addition of EHD and DEP forces allows the loop to be sustained in the absence of gravity, where such a buoyancy force/gravity is negligible or not present.

The electric potential applied to the EHD pump was ramped up from 0 to 1500 V in 500 V increments. The incremental increase in voltage is required due to the sensitivity of forming the heterocharge layers on the conduction pump electrodes. Based on the spacing between electrodes, the voltage was limited to 1500 V to avoid dielectric breakdown of the fluid. The DEP electrode was either kept at 0 V or energized at 2000 V. Due to the nature of the parabolic microgravity flight, voltages of the EHD conduction pump and the DEP electrode were changed only during longer segments where the plane was leveled out. This is due to electrical steady-state taking several minutes to be reached. The microgravity flight consists of a series of parabolic maneuvers in which those onboard experience accelerations ranging from hypergravity (1.6 g to 1.8 g) on the lower vertex, to microgravity (≈0 g) at the higher vertex. The microgravity increments typically last between 15-20 s, and the hypergravity increments typically lasted 20-30 s. The flight campaign for this study consisted of a total of 120 flight parabolas conducted over 4 consecutive days. During the microgravity flights, the heater power was operated from 0 W/cm2 to around 7 W/cm2. The heater power was changed between parabolas to develop boiling curves from 0 W/cm2 to dryout values for various cases of EHD and DEP voltages.

In addition to the microgravity flight, the experiment was further tested terrestrially on the ground using an updated heater design, which has an area of 1.0 cm2, as well as a smaller DEP electrode to accommodate the heater. The purpose of this testing was to obtain full boiling curves and determine CHF with the enhancement of the EHD conduction pump and DEP mechanism in a 1 g environment. CHF for this ground testing occurred when the applied heat flux caused the heater surface temperature to rise very rapidly over the course of a few seconds. This rapid rise in temperature is caused by the boiling over the heater being strong enough to overcome the replenishment of cooler liquid HFE 7100. The CHF was increased from 9.76 W/cm2 under no electric fields to 15.63 W/cm2 with the application of electric fields.

Results and Discussions

Terrestrial Results.

Figure 6 shows the heat flux as a function of temperature difference between the heater surface and vapor for the cases of baseline (no EHD/DEP), EHD only, DEP only, and combined EHD and DEP cases for data taken on the ground. The ground data show that EHD conduction pumping enhances heat transfer by providing an electrically-driven forced convection over the heater surface. This is attributed to an increased flowrate of liquid over the heater and local liquid circulation near the heater due to the EHD conduction pumping. This is confirmed by video taken during testing. This also results in an increased critical heat flux. However, the enhancements are not significant because the pressure generated by the EHD conduction pump is only slightly higher than that already provided by hydrostatic pressure at 15 °C. This is expected since the hydrostatic pressure provided by the 2 mm fluid film is approximately 30 Pa, while the EHD pump is able to provide pressure generation of approximately 50 Pa at 1.5 kV. Therefore, naturally, in microgravity, the pump has more impact, as there is no hydrostatic pressure in the fluid. It is important to recognize that the role of the EHD conduction pumping is essential in the absence of gravity to electrowet the heater surface to maintain a fluid film for boiling.

Fig. 6
Heat flux versus temperature difference between heater surface and vapor showing electric field enhancement in the presence of gravity (g = 9.81 m/s2). Rohsenow correlation is included for comparison.
Fig. 6
Heat flux versus temperature difference between heater surface and vapor showing electric field enhancement in the presence of gravity (g = 9.81 m/s2). Rohsenow correlation is included for comparison.
Close modal

The DEP extraction mechanism significantly enhances the heat transfer. Higher heat fluxes are reached with lower heater surface superheat. For example, at the heat flux of 9.76 W/cm2, the heater surface temperature is lowered by 18.65 °C compared to the baseline case. The case of EHD and DEP mechanisms combined further enhance heat transfer as shown in Fig. 6. This is clear as the CHF reaches its maximum value of 15.63 W/cm2 compared to only 9.76 W/cm2 for the baseline case. Considering the combined power consumption (maximum of 400 mW) of the DEP electrode and EHD conduction pump, this result demonstrates a significant improvement in heat removal at the cost of low power. For qualitative comparison, the Rohsenow correlation for boiling is estimated below and presented in Fig. 6.

The heat transfer characteristics of nucleate boiling can be estimated using the correlation developed by Rohsenow [28]
(2)

Here, Csf and r are constants that are specific to the system and the working fluid properties. Equation (2) has been rearranged from the original equation given by Ref. [28] to solve for ΔT as a function of heat flux, q″. Fluid properties from Ref. [25] along with the constants (Csf = 0.0036, r = 0.33) [29] were used to find these values of ΔT. These temperature values were calculated for the case of a terrestrial environment (i.e., 1 g = 9.81 m/s2). The 1 g prediction based on this correlation is compared to the experimental ground data and is plotted in Fig. 6.

The trend qualitatively agrees with the data in the baseline case. Patel et al. [27] conducted a similar Rohsenow analysis and also found that the correlation qualitatively aligned with their baseline (no electric field present) data. Comparing the data obtained in [27] with this study reveals similar trends, but a direct comparison is not practical, as the fluids used were different in each study (HFE 7100 versus HCFC-123). Again, this comparison is qualitative, as Rohsenow is used for pool boiling, and this study evaluates liquid film flow boiling.

Figure 7 shows the resultant heat transfer coefficient as a function of temperature difference between the heater surface and vapor for the same cases on the ground given in Fig. 6. Activating the DEP electrode significantly raises the heat transfer coefficient, by as much as three times when compared to the case of no EHD or DEP influence. The enhanced heat transfer coefficient with the DEP mechanism activated remains constant at lower ΔT values below the ONB. At these values, the DEP force stirs the liquid film (in the absence of vapor bubbles). This is due to the electric permittivity gradient that is expected to be present within the liquid film height because of the presence of a temperature gradient. Further enhancement is observed with the combined operation of EHD pumping and DEP extraction force.

Fig. 7
Heat transfer coefficient versus temperature difference between heater surface and vapor in the presence of gravity (g = 9.81 m/s2)
Fig. 7
Heat transfer coefficient versus temperature difference between heater surface and vapor in the presence of gravity (g = 9.81 m/s2)
Close modal

Hysteresis Results.

Figure 8 illustrates the impact of presence of EHD conduction pumping and DEP vapor extraction mechanism on hysteresis, on the ground. As seen in the figure, in the absence of EHD pumping and DEP mechanism, there is a clear hysteresis present when comparing increasing versus decreasing heat flux curves. Furthermore, the onset of nucleate boiling (ONB) occurs at a ΔT of 16.6 °C, corresponding to 1.75 W/cm2.

Fig. 8
Comparison of heat flux versus temperature difference for baseline, EHD = 1500V and DEP = 2000V showing hysteresis
Fig. 8
Comparison of heat flux versus temperature difference for baseline, EHD = 1500V and DEP = 2000V showing hysteresis
Close modal

In the presence of EHD conduction pumping or DEP mechanism, as depicted in Fig. 8, the hysteresis almost gets eliminated. This is primarily due to the fact that EHD conduction pumping effectively wets the heater surface as it provides forced convection, and the DEP mechanism extracts bubbles as they get formed. The DEP mechanism also raises the ONB from 1.75 W/cm2 to 2.95 W/cm2, which was visually observed during testing. This delay in the ONB with the DEP mechanism activated can be attributed to additional heat transfer taking place before boiling. As mentioned before, the DEP force stirs the liquid film due to the presence of an electric permittivity gradient in the liquid film height.

Gravity Condition Results.

The compiled flight results are shown in Figs. 911(d). Due to the nature of parabolic flights and the short time allocated between changes in gravity conditions, inertial forces result in fluid avalanching between transitions. Therefore, the results correspond to unsteady state situations. One the other hand, the ground data were obtained under steady-state. For reference, the data show the time required to reach steady-state on the ground was approximately 10-15 s. The data from the parabolic flights were taken over different, but consecutive days, which could also contribute to challenges in comparing data points. Despite the above issues with testing in parabolic flights, the results show that the EHD conduction pumping, and, to much extent, the DEP mechanism enhanced the heat transfer. To obtain full boiling curves and determine the true impact of these mechanisms in the absence of gravity at steady-state, a steady microgravity environment, such as the International Space Station, is needed.

Fig. 9
Heat flux versus temperature difference between heater surface and vapor in the absence of gravity (g = 0 m/s2)
Fig. 9
Heat flux versus temperature difference between heater surface and vapor in the absence of gravity (g = 0 m/s2)
Close modal

Figure 9 presents the results of the microgravity sections of the parabolic flights. The most important trend to notice in this plot is the effect of the DEP electrode on the liquid film flow boiling heat transfer. The DEP electrode's capability to extract vapor raises the heat flux to its highest value of 7 W/cm2 before dryout. This value is considered the dryout point since the heater surface temperature rises very rapidly over a few seconds. It should be noted that dryout is not necessarily the same as CHF, as the experiment is not at steady-state, and the nature of the microgravity flight is not controllable. This increase in heat flux is 1.5 W/cm2 or a 27.3% increase compared to the dryout point of the baseline case. Furthermore, this trend exists for lower ranges of heat flux (0–4 W/cm2), indicating that high heat flux values are not necessary to capture the heat transfer enhancement that these EHD mechanisms can provide in microgravity. Additionally, combining the EHD conduction pump with the DEP electrode provides further enhancement. The combined effect of these mechanisms reached higher heat fluxes for a given temperature difference compared to baseline and EHD pumping-only cases. This is shown and labeled in Fig. 9 as a difference of 15 °C exists for the same heat flux value between the combined EHD/DEP case and the baseline case. This result demonstrates that using EHD and DEP together under unsteady state, a heated surface experiencing liquid film flow boiling in the absence of gravity can reach higher heat flux values and still maintain a lower surface temperature when compared to traditional baseline cases where EHD and DEP are not present. Without gravity, the role of EHD conduction pumping is critical in providing the liquid to the heater surface. Regardless of the transient environment of the flights, the DEP mechanism shows significant enhancement of heat transfer raising the CHF to 7 W/cm2.

Similar to the microgravity heat flux results, a set of results in 1.6 g to 1.8 g ranges was also plotted in the same fashion. Figure 10 highlights these results and shows that the DEP electrode continues to dominate the heat transfer enhancement, especially toward dryout values. The combined effect of the EHD conduction pump and DEP electrode also provides enhancement, especially at lower heat flux values, and results in a 9 °C drop in temperature for the same heat flux value when compared to the baseline test case. This result indicates that these mechanisms could provide heat transfer enhancement, even under harsh conditions associated with parabolic flights.

Fig. 10
Heat flux versus temperature difference between heater surface and vapor in the presence of hypergravity (g = 1.8 × 9.81 m/s2)
Fig. 10
Heat flux versus temperature difference between heater surface and vapor in the presence of hypergravity (g = 1.8 × 9.81 m/s2)
Close modal

Figures 11(a)11(d) directly compare the ground, microgravity, and 1.8 g data for different test cases of baseline, EHD, DEP, and combined EHD with DEP. These figures demonstrate that these mechanisms provide very similar enhancement in microgravity as compared to the ground, indicating that they could be effective not only for terrestrial applications, but also in space. It should be noted that the ground data reach significantly higher than the microgravity or 1.8 g data as a result of the updated hardware mentioned earlier. The heat flux values in the microgravity data are observed to plateau as this is associated with the dryout point. To elaborate, nearing dryout, the heat flux remains constant, while the superheat (ΔT) increases rapidly. Thus, a plateau is observed in these plots. Also indicated by these plots is the significance of the DEP electrode as a method of thermal management, reaching higher heat fluxes than those of EHD pumping only or baseline. Additionally, for a given value of heat flux, the heater surface temperature is decreased the most under the influence of the DEP electrode in all gravitational conditions (unsteady state) as compared to the baseline and EHD pumping only cases, up to a difference of 12 °C.

Fig. 11
(a) Heat flux versus temperature difference between liquid and vapor for ground (g = 9.81 m/s2), microgravity (g = 0 m/s2), and 1.8 g (g = 1.8 × 9.81 m/s2) conditions at baseline condition (0 V EHD, 0 V DEP)l; (b) Heat flux versus temperature difference between liquid and vapor for ground (g = 9.81 m/s2), microgravity (g = 0 m/s2), and 1.8 g (g = 1.8 × 9.81 m/s2) conditions at EHD = 1500V, DEP = 0V; (c) Heat flux versus temperature difference between liquid and vapor for ground (g = 9.81 m/s2), microgravity (g = 0 m/s2), and 1.8 g (g = 1.8 × 9.81 m/s2) conditions at EHD = 0V, DEP = 2000V; and (d) Heat flux versus temperature difference between liquid and vapor ground (g = 9.81 m/s2), microgravity (g = 0 m/s2), and 1.8 g (g = 1.8 × 9.81 m/s2) conditions at EHD = 1500V, DEP = 2000V
Fig. 11
(a) Heat flux versus temperature difference between liquid and vapor for ground (g = 9.81 m/s2), microgravity (g = 0 m/s2), and 1.8 g (g = 1.8 × 9.81 m/s2) conditions at baseline condition (0 V EHD, 0 V DEP)l; (b) Heat flux versus temperature difference between liquid and vapor for ground (g = 9.81 m/s2), microgravity (g = 0 m/s2), and 1.8 g (g = 1.8 × 9.81 m/s2) conditions at EHD = 1500V, DEP = 0V; (c) Heat flux versus temperature difference between liquid and vapor for ground (g = 9.81 m/s2), microgravity (g = 0 m/s2), and 1.8 g (g = 1.8 × 9.81 m/s2) conditions at EHD = 0V, DEP = 2000V; and (d) Heat flux versus temperature difference between liquid and vapor ground (g = 9.81 m/s2), microgravity (g = 0 m/s2), and 1.8 g (g = 1.8 × 9.81 m/s2) conditions at EHD = 1500V, DEP = 2000V
Close modal

Dielectrophoretic Force Analysis.

To estimate the electric field distribution of the DEP electrode, an analysis in COMSOL Multiphysics was carried out. This analysis followed the same procedure as Ref [27], however, with updated geometry and different fluid properties to reflect the experiment presented in this paper. Solving this in COMSOL indicated that the DEP electrode would induce a strong electric field distribution over the heater. This electric field distribution was determined in the absence of vapor bubbles as an approximation. The magnitude of the electric field reached a maximum value of 5.0 × 106 V/m, mostly near the edges of the heater and DEP electrode. In addition to calculating the electric field distribution, DEP force calculations were performed to directly compare the force induced by the electrodes to buoyancy forces. This is important as it shows how the DEP electrode can significantly increase the force a vapor bubble experiences as opposed to without it. The DEP force [29] acting on a vapor bubble with a radius a is as follows [30]:
(3)
Here, ε1 and ε2 represent the permittivities of two different mediums. In the case of this study, the liquid medium is HFE 7100 which has a dielectric constant of 7.4 [26], and a permittivity of ε1 = 6.55196 × 10−11 F/m. For the vapor bubble, the permittivity of free space is utilized, such that ε2 = 8.854 × 10−12 F/m. Since the DEP force is directly proportional to the gradient of the electric field magnitude squared, a nonuniform dielectric field will result in a DEP force that acts upon the vapor phase. More specifically this results in a diverging electric field where the vapor phase (i.e., the bubble) will move toward the lower electric field. A strong diverging electric field is made possible by designing the DEP electrode to be larger than the heater and by applying a high enough electric potential. To determine the expected bubble radius for an imposed heat flux and temperature difference, the following correlation [31] was used
(4)
Fluid properties for HFE 7100 are supplied in Ref. [8], and the acceleration due to gravity is g =9.81 m/s2. Values for heat flux and temperature difference were selected as q″ = 2 W/cm2 and ΔT =15 °C. These values represented nucleate boiling for the baseline ground tests shown in Fig. 6. All of these inputs result in a bubble diameter of roughly 0.2 mm, giving a radius of approximately 0.1 mm. Therefore, this radius was chosen for the DEP force and buoyancy force calculations. Taking advantage of the DEP phenomenon in Eq. (3) allows for vapor bubbles to be extracted at forces much higher than those induced by buoyancy forces alone. To demonstrate, the buoyancy force on a bubble of radius a is calculated using the following equation [32]
(5)

This results in a buoyancy force of Fbuoyancy = 6.02 × 10−8 N. The maximum DEP force on a vapor bubble was calculated to reach maximum force of 55.0 × 10−8 N, resulting in a ratio of DEP force to buoyancy force of 9.14. This indicates that the DEP mechanism can effectively extract bubbles away from the heater surface. In addition to providing vapor extraction in the absence of gravity, this also indicates that on the ground, the DEP electrode will still be a superior vapor extraction method compared to the natural buoyancy forces induced by gravity as illustrated in Fig. 6.

Error Analyses.

Experimental uncertainty in all measurements and the derived values are compiled and represented in Table 2. The horizontal and vertical error bars depicted in the data graphs correspond to two standard deviations. Overall, the derived error contributions come from the ΔT values (±0.7 °C), the heat flux values (±3.1%), and the heat transfer coefficient values (±14%).

Table 2

Maximum systematic uncertainty of measured and derived values

MeasurementMaximum Uncertainty
Temperature±0.5 °C
Absolute pressure±40 Pa
Voltage (EHD)±2.5 VDC
Current (EHD)±2.5 μA
Voltage (DEP)±2.5 VDC
Current (DEP)±2.5 μA
Voltage (Heater)±30 mVDC
Current (Heater)±2 mA
MeasurementMaximum Uncertainty
Temperature±0.5 °C
Absolute pressure±40 Pa
Voltage (EHD)±2.5 VDC
Current (EHD)±2.5 μA
Voltage (DEP)±2.5 VDC
Current (DEP)±2.5 μA
Voltage (Heater)±30 mVDC
Current (Heater)±2 mA
Derived valueMaximum uncertainty
Heat transfer coefficient, h±14%
Heat flux, q±3.1%
ΔT±0.7 °C
Derived valueMaximum uncertainty
Heat transfer coefficient, h±14%
Heat flux, q±3.1%
ΔT±0.7 °C

Conclusions

This study investigated the effect of EHD conduction pumping and the DEP mechanism under various gravitational conditions. The testing conditions aboard the parabolic flight corresponded to unsteady state. Regardless of this, both mechanisms illustrated that they enhance heat transfer, even during harsh conditions associated with parabolic flights, especially when combined. Additionally, they also allow for sustainable liquid film boiling in the absence of gravity. It was also observed that the presence of electric fields from the EHD and DEP mechanisms eliminated boiling hysteresis. Throughout all operating conditions, the maximum current drawn by the EHD and DEP was less than 0.3 mA with the average being substantially lower. This resulted in less than at most 600 mW of power drawn, with the average being much lower, highlighting this technology's advantage of low power consumption while enabling two-phase heat transfer in microgravity. This study paves the way for future implementation of DEP/EHD-driven two-phase heat transport devices into space and aeronautical applications.

Acknowledgment

The work was supported by the Biological and Physical Sciences Division in the Science Mission Directorate at NASA HQ under NASA Grant NNX16AT09G. The work of the first author was supported by the NASA Space Technology Research Fellowship NSTRF Grant 80NSSC19K1170.

Funding Data

  • National Aeronautics and Space Administration (NASA) (Funder ID: 10.13039/100000104).

Nomenclature

a =

bubble radius (m)

cp =

specific heat (kJ/kg K)

Csf =

surface-fluid constant for Rohsenow correlation

d =

characteristic length (m)

E =

electric field (V/m)

Ee =

electric field vector (V/m)

Fbuoyancy =

buoyancy force (N)

FDEP =

DEP force (N)

fEHD =

electrohydrodynamic body force (N)

g =

gravity (9.81 m/s2)

hfg =

heat of vaporization (kJ/kg)

kl =

liquid thermal conductivity (W/m K)

Prl =

liquid Prandtl number

q″ =

heat flux (W/m2)

r =

system-specific constant for Rohsenow correlation

T =

temperature (K or °C)

V =

electric potential (V)

ΔT =

temperature difference (K or °C)

ε =

electric permittivity (F/m)

μ =

dynamic viscosity (kg/m s)

ρl =

liquid density (kg/m3)

ρv =

vapor density (kg/m3)

σ =

surface tension (N/m)

References

1.
Melcher
,
J.
,
1981
,
Continuum Electromechanics
,
MIT Press
,
Cambridge, MA
.
2.
Yagoobi
,
J.
,
2005
, “
Electrohydrodynamic Pumping of Dielectric Liquids
,”
J. Electrostatics
,
63
, pp.
861
869
.10.1016/j.elstat.2005.03.047
3.
Atten
,
P.
, and
Seyed-Yagoobi
,
J.
,
2003
, “
Electrohydrodynamically Induced Dielectric Liquid Flow Through Pure Conduction in Point/Plane Geometry
,”
IEEE Trans. Dielectrics Electr. Insul.
,
10
(
1
), pp.
27
36
.10.1109/TDEI.2003.1176555
4.
Patel
,
V.
, and
Seyed-Yagoobi
,
J.
,
2014
, “
Recent Experimental Advances in Electrohydrodynamic Conduction Pumping Research
,”
IEEE Industry Application Society Annual Meeting
, Vancouver, BC, Canada, Oct. 5–9, pp.
1
10
.10.1109/IAS.2014.6978363
5.
Vázquez
,
P. A.
,
Talmor
,
M.
,
Seyed-Yagoobi
,
J.
,
Traoré
,
P.
, and
Yazdani
,
M.
,
2019
, “
In-Depth Description of Electrohydrodynamic Conduction Pumping of Dielectric Liquids: Physical Model and Regime Analysis
,”
Phys. Fluids
,
31
(
11
), p.
113601
.10.1063/1.5121164
6.
Yazdani
,
M.
, and
Seyed-Yagoobi
,
J.
,
2009
, “
Electrically Induced Dielectric Liquid Film Flow Based on Electric Conduction Phenomenon
,”
IEEE Trans. Dielectrics Electr. Insul.
,
16
(
3
), pp.
768
777
.10.1109/TDEI.2009.5128517
7.
Talmor
,
M.
, and
Seyed-Yagoobi
,
J.
,
2021
, “
Numerical Study of Micro-Scale EHD Conduction Pumping: The Effect of Pump Orientation and Flow Inertia on Heterocharge Layer Morphology and Flow Distribution Control
,”
J. Electrostatics
,
111
, p.
103548
.10.1016/j.elstat.2020.103548
8.
Pearson
,
M.
, and
Seyed-Yagoobi
,
J.
,
2013
, “
Electrohydrodynamic Conduction Driven Single- and Two-Phase Flow in Microchannels With Heat Transfer
,”
ASME J. Heat Transfer-Trans. ASME
,
135
(
10
), p.
101701
.10.1115/1.4007617
9.
Talmor
,
M.
,
Louste
,
C.
, and
Seyed-Yagoobi
,
J.
,
2018
, “
PIV Flow Field Measurements of Electrohydrodynamic Conduction Pumping
,”
Proc. of the 2018 Electrostatics Joint Conference Conf
, Boston, MA, June, pp.
1
15
.
10.
Louste
,
C.
,
Traoré
,
P.
,
Vazquez
,
P. A.
, and
Seyed Yagoobi
,
J.
,
2019
, “
PIV Flow Measurements of Conduction Pumping Flow Created by Nine Pairs of Asymmetric Surface Electrodes
,” 2019 IEEE 20th International Conference on Dielectric Liquids (
ICDL
), Rome, Italy, June 23–27, pp.
1
4
.10.1109/ICDL.2019.8796777
11.
Vázquez
,
P. A.
,
Seyed-Yagoobi
,
J.
,
Traoré
,
P.
, and
Louste
,
C.
,
2019
, “
EHD Pumping in Flexible Conic Nozzle
,”
IEEE 20th International Conference on Dielectric Liquids
(
ICDL
), Rome, Italy, June 23–27, pp.
1
4
.10.1109/ICDL.2019.8796552
12.
O'Connor
,
N.
,
Castaneda
,
A.
, and
Yagoobi
,
J.
,
2020
, “
Experimental Study of Flexible Electrohydrodynamic Conduction Pumping for Electronics Cooling
,”
J. Electron Packag.
,
142
(
4
), p.
041105
.10.1115/1.4047459
13.
Ogata
,
J.
, and
Yabe
,
A.
,
1993
, “
Augmentation of Boiling Heat Transfer by Utilizing the EHD Effect—EHD Behaviour of Boiling Bubbles and Heat Transfer Characteristics
,”
Int. J. Heat Mass Transfer
,
36
(
3
), pp.
783
791
.10.1016/0017-9310(93)80054-X
14.
Seyed-Yagoobi
,
J.
,
Geppert
,
C. A.
, and
Geppert
,
L. M.
,
1996
, “
Electrohydrodynamically Enhanced Heat Transfer in Pool Boiling
,”
ASME J. Heat Transfer-Trans. ASME
,
118
(
1
), pp.
233
237
.10.1115/1.2824048
15.
Seyed-Yagoob
,
J.
,
Hardesty
,
J. T.
,
Raghupathi
,
P.
, and
Bryan
,
J. E.
,
1997
, “
Experimental Study of Electrohydrodynamically Augmented Pool Boiling Heat Transfer on Smooth and Enhanced Tubes
,”
J. Electrost.
,
40–41
, pp.
597
602
.10.1016/S0304-3886(97)00109-5
16.
Darabi
,
J.
,
Ohadi
,
M. M.
, and
Dessiatoun
,
S. V.
,
2000
, “
Augmentation of Thin Falling-Film Evaporation on Horizontal Tubes Using an Applied Electric Field
,”
ASME J. Heat Transfer-Trans. ASME
,
122
(
2
), pp.
391
398
.10.1115/1.521478
17.
Kano
,
I.
, and
Takahashi
,
Y.
,
2013
, “
Effect of Electric Field Generated by Microsized Electrode on Pool Boiling
,”
IEEE Trans. Ind. Appl
,.,
49
(
6
), pp.
2382
2387
.10.1109/TIA.2013.2263213
18.
Kano
,
I.
,
Higuchi
,
Y.
, and
Chika
,
T.
,
2013
, “
Development of Boiling Type Cooling System Using Electrohydrodynamics Effect
,”
ASME J. Heat Transfer-Trans. ASME
,
135
(
9
), p.
091301
.10.1115/1.4024390
19.
Garivalis
,
A. I.
,
Manfredini
,
G.
,
Saccone
,
G.
,
Di Marco
,
P.
,
Kossolapov
,
A.
, and
Bucci
,
M.
,
2021
, “
Critical Heat Flux Enhancement in Microgravity Conditions Coupling Microstructured Surfaces and Electrostatic Field
,”
NPJ Microgravity
,
7
(
1
), p.
37
.10.1038/s41526-021-00167-3
20.
Patel
,
V. K.
,
Yagoobi
,
J.
,
Robinson
,
F.
, and
Didion
,
J. R.
,
2016
, “
Effect of Gravity on Electrohydrodynamic Conduction Driven Liquid Film Flow Boiling
,”
AIAA J. Thermophys. Heat Transfer
,
30
(
2
), pp.
429
437
.10.2514/1.T4696
21.
Schweizer
,
N.
, and
Stephan
,
P.
,
2009
, “
Experimental Study of Bubble Behavior and Local Heat Flux in Pool Boiling Under Variable Gravitational Conditions
,”
Multiphase Sci. Technol.
,
21
(
4
), pp.
329
350
.10.1615/MultScienTechn.v21.i4.40
22.
Di Marco
,
P.
, and
Grassi
,
W.
,
2009
, “
Effect of Force Fields on Pool Boiling Flow Patterns in Normal and Reduced Gravity
,”
Heat Mass Transfer
,
45
(
7
), pp.
959
966
.10.1007/s00231-007-0328-6
23.
Henry
,
C. D.
, and
Kim
,
J.
,
2004
, “
Heater Size, Subcooling, and Gravity Effects on Pool Boiling Heat Transfer
,”
Int. J. Heat Fluid Flow
,
25
(
2
), pp.
262
273
.10.1016/j.ijheatfluidflow.2003.11.019
24.
Dhir
,
V. K.
,
Warrier
,
G. R.
,
Aktinol
,
E.
,
Chao
,
D.
,
Eggers
,
J.
,
Sheredy
,
W.
, and
Booth
,
W.
,
2012
, “
Nucleate Pool Boiling Experiments (NPBX) on the International Space Station
,”
Microgravity Sci. Technol.
,
24
(
5
), pp.
307
325
.10.1007/s12217-012-9315-8
25.
3M Novec 7100 Engineeering Fluid
,” 3M, accessed Sept. 15, 2022, https://multimedia.3m.com/mws/media/1613377O/3mtm-novectm-7100-engineered-fluid-technical-data-sheet.pdf
26.
Nassar
,
M.
,
Vázquez
,
P. A.
,
Chauris
,
N.
,
Daaboul
,
M.
,
Michel
,
A.
, and
Louste
,
C.
,
2020
, “
Experimental Models of the Variation of HFE-7100 and HFE-7000 Electric Properties With Temperature
,”
IEEE Trans. Ind. Appl.
,
56
(
4
), pp.
4193
4199
.10.1109/TIA.2020.2990367
27.
Patel
,
V. K.
, and
Seyed-Yagoobi
,
J.
,
2017
, “
Combined Dielectrophoretic and Electrohydrodynamic Conduction Pumping for Enhancement of Liquid Film Flow Boiling
,”
ASME J. Heat Transfer-Trans. ASME
,
139
(
6
), p.
061502
.10.1115/1.4035709
28.
Rohsenow
,
W. M.
,
1952
, “
Method of Correlating Heat-Transfer Data for Surface Boiling of Liquids
,”
ASME Trans.
,
74
(
6
), pp.
969
975
.10.1115/1.4015984
29.
Manetti
,
L. L.
,
Moita
,
A. S. O. H.
,
de Souza
,
R. R.
, and
Cardoso
,
E. M.
,
2020
, “
Effect of Copper Foam Thickness on Pool Boiling Heat Transfer of HFE-7100
,”
Int. J. Heat Mass Transfer
,
152
, p.
119547
.10.1016/j.ijheatmasstransfer.2020.119547
30.
Pohl
,
H. A.
,
1978
, ”
Dielectrophoresis: The Behavior of Neutral Matter in Nonuniform Electric Fields
,
Cambridge University Press
,
Cambridge, UK
.
31.
Zuber
,
N.
,
1959
, “
Hydrodynamic Aspects of Boiling Heat Transfer
,”
U.S. AEC
,
Washington, DC
, Report No. AECU
-4439.
32.
Jones
,
T. B.
, and
Bliss
,
G. W.
,
1977
, “
Bubble Dielectrophoresis
,”
J. Appl. Phys.
,
48
(
4
), pp.
1412
1417
.10.1063/1.323806