1. Introduction
Cavitation has received more and more attention during the last decade due to the potential for its implementation in many applications, such as biomedical ones [
1], energy harvesting [
2], surface cleaning, and water treatment [
3,
4,
5]. Cavitation is associated with a process involving the nucleation, growth, and implosion of bubbles filled with vapor or gas. Such cavities are generated in a fluid when the static pressure drops below the vapor saturation pressure of the fluid. These bubbles implode violently when encountering a higher pressure region, resulting in high temperature spots. The static pressure drop can be obtained using acoustic waves generated by ultrasonic transducers (acoustic cavitation) or by suddenly reducing the hydrodynamic pressure using geometrical changes, typically by increasing the velocity of the fluid or by boundary layer separation (hydrodynamic cavitation) [
6,
7]. Hydrodynamic cavitation is known to be more energy efficient because of rather simple configuration and construction of related systems, low energy requirements, and easier scaling up possibility for industrial applications [
8,
9,
10].
Hydrodynamic cavitation can be observed when the fluid passes through a physical constriction, such as an orifice plate [
11] or a venturi [
12]. Orifices are one of the most commonly used geometries as hydrodynamic cavitation generators. Despite the extensive use of orifices, the functions and interactions among different designs, and operating parameters of the resulting cavitation behavior, have not yet been sufficiently understood. The geometry of devices, thermophysical properties of the working fluids, surface configurations, operating temperatures and pressures, and contaminant gases in the fluid play important roles in the cavitation phenomenon.
The geometry of the flow restrictive element is one of the most important parameters in the hydrodynamic cavitation phenomenon. There are numerous studies on geometrical parameters of cavitation generators (flow restrictive elements), such as venturis and orifice plates, in the literature. For instance, He et al. studied cavitating flows in a diesel injector nozzle and showed that hydrodynamic cavitation was beneficial in the separation process of sprays. It was concluded that the radius of curvature (ROC) of the entrance hole was the most crucial parameter for a nozzle in affecting internal cavitating flow characteristics [
13]. Rudolf et al. investigated hydrodynamic cavitation in single and multi-hole orifices. Their experimental results illustrated that using multi-hole orifices was more effective than single-hole orifices due to their low energy dissipation. Moreover, based on their study, multi-hole orifices lead to a significant drop of loss coefficient in cavitating and non-cavitating regimes [
14]. Dong et al. experimentally studied the effects of number, arrangement, and ratio of the holes in an orifice plate. They proved that cavitation number was minimal for a diagonal structure, while it was maximal for a radial structure. Besides, the turbulent kinetic energy had the maximum value near the edges, and minimum values were seen at the axial lines [
15]. Ranade et al. numerically studied the influence of critical geometrical parameters such as orifice thickness, hole inlet sharpness, and wall angle on cavitation flow patterns. According to this study, orifice geometry had a significant influence on cavitation inception. In contrast with the rounded edged orifice, the sharp edged orifice generated more intense cavitating flows. In addition, the required pressure for cavitation inception was increased by 30%–40% with angled orifice walls in comparison to a straight throat section, and as a result, angled orifice walls were undesirable for intensified cavitation [
16].
Cavitating flows in micro scale differ from those in macro scale. For instance, experimental investigation of cavitating flows in micro-orifices with a rectangular cross-section showed that cavitating flow patterns in micro-orifices are different from macroscale orifices. It was demonstrated that incipient and choking cavitation numbers, and cavitating flow patterns of micro-orifices are heavily influenced by the size scale [
17]. In this regard, Jin et al. numerically studied the effect of the ratio between the length and diameter of micro-orifices. They reported that the vapor cavity increased with an increase in length to diameter
(l/d) ratio, and there existed an optimum value for
l/d of a micro-orifice based on both cavitation intensity and flow rate [
18]. When sufficiently many cavitation bubbles collapse, they release a high amount of energy downstream of the flow channel. This released energy could be utilized to provide the required power for some applications, such as water treatment, surface cleaning, biomedical applications, and energy harvesting. Energy harvesting with hydrodynamic cavitation flows could, thus, be realized in a cost-effective and environmentally-friendly fashion. Surface topology and geometry of the device are among the most significant parameters, which lead to an increase in the intensity of cavitation flows while causing earlier inception. Moreover, the lower inlet pressure for the inception and intensified cavitating flows will result in a higher amount of energy from the collapse of bubbles and improved performance. In this regard, Ghorbani et al. focused on the effects of surface roughness of microfluidic devices on hydrodynamic cavitation. Their results showed that the roughness on the channel surface led to intensified cavitating flows in comparison with a smooth channel over the same range of flow rate. They also suggested that longer microchannels were suitable for energy harvesting applications because of the penetration of the emerging twin cavities [
19]. It should be noted that the entire downstream region of microchannel could be filled with vapor and large vapor bubbles by causing twin cavities, which could be exploited in energy harvesting.
Hydrodynamic cavitation is directly affected by thermophysical properties of the working fluid. Therefore, it is crucial to consider this key parameter in microfluidic devices. There are a few studies addressing the application of other working fluids and its effect on flow patterns. For example, Mossaz et al. [
20] reported that the mixture of isopropanol and water could change the cavitation inception from turbulent regime to laminar regime. Low percentage of isopropanol in water is likely to increase the dynamic viscosity so that cavitation inception corresponding to laminar flow could be more easily achieved. Hydrodynamic cavitation of water, ethanol, and refrigerant liquids in micro venturis was investigated by Mishra et al. [
21]. They compared the results of R-123 to those of water and observed that lower surface tension of this fluid enhanced bubble formation. In a subsequent study, different types of liquids were utilized to obtain flow patterns at low upstream pressures for similar devices [
22]. A rapid transition from inception to choked flow condition, beyond which the flow rate no longer increased with upstream pressure, was observed. Ghorbani et al. [
23] used poly(vinyl alcohol) (PVA) microbubbles in water inside microchannels with different geometries and compared the results with those of water. They could show that the devices working with PVA had lower cavitation number in comparison with the case of pure water under the same working conditions. The experimental results from this study also demonstrated that the impact pressure upon bubble collapse was more for the case of PVA. Li et al. [
24] investigated the influence of thermodynamic effects in the evaporation and condensation coefficient of the cavitation model using an accurate model based on the Rayleigh–Plesset equation. They developed a computational fluid dynamics (CFD) model for simulation of two airfoil cavitation. Their results indicated that the fluid temperature had an effect on the cavitation intensity around the hydrofoil, and the cavitation became more intense with the increase in the temperature. In order to study the cavitation phenomenon and energy generation from hydrodynamic cavitation, Gevari et al. [
25] fabricated micro orifices on chips and tested them with different working fluids, including water and a perfluoropentane (PFC5) droplet-water suspension. Based on their results, cavitation inception occurred earlier for the droplet-water suspension. In the case of droplet-water suspension, supercavitation condition appeared sooner, and as a result, the input energy of the system was minimized, which suggested higher energy efficiency.
Energy harvesting with hydrodynamic cavitation has attracted considerable attention during recent years. However, only few attempts have been made regarding the design and optimization of cavitation on chip devices for the purpose of energy harvesting. In this study, the effects of two key parameters in the design and optimization, namely, working fluid and roughness elements inside the microfluidic device, were investigated. In microscale systems, the surface tension force becomes more dominant over body forces such as gravity because of increasing surface to volume ratio. Consequently, surface roughness significantly affects the flow behavior in microfluidic devices. Furthermore, surface roughness is a vital parameter for multiphase flows, where the inception of the secondary phase is from the surface (i.e., boiling, cavitation, condensation). The surface roughness and cavity size are critical parameters in bubble nucleation. In phase change processes, bubble nucleation and bubble departure frequency increase with surface roughness up to a critical size, beyond which surface roughness has no considerable effect on the nucleation process. In this study, the microfluidic devices with a single orifice were fabricated by conventional microfabrication methods. The microfluidic devices were tested at different inlet pressures, and the cavitating flow patterns were visualized. The performances of the microfluidic devices were evaluated, and wall roughness elements were added to the originally fabricated plain surface microfluidic device. The effect of these elements on cavitation flow patterns was studied, and the optimum design area was obtained. Changing the working fluid to intensify cavitating flows in the extension area was the last step in optimization. Finally, an optimum design was suggested for microfluidic devices for achieving facile and intensified cavitating flows.
3. Results and Discussion
In the fabricated microfluidic devices, the heat dissipation upon cavitating bubble collapse was meant to be harvested and converted into electricity by a thermoelectric module installed at the end wall of the extension, where bubbles were targeted to. As a result, intense cavitating flows at lower inlet pressures were to increase the efficiency of the energy harvesting system. In other words, inception and intensified cavitating flows in the microfluidic device were desired to be at lower pressures. As a result, lower inception and cavitation development pressures of cavitating flows constitute the optimization goals of this study. The optimization process was performed in three phases, which are explained in the next sections.
3.1. Phase One: Initial Design and Analysis
The inlet and extension width of the initial design, as mentioned before, was 900 µm, and the width of the microchannel was 152 µm. The widths and lengths of the microfluidic devices were kept constant during the optimization process. Cavitation number, Equation (4), as an indicator to characterize the cavitation phenomena, is considered here as an optimization parameter along with the inlet pressure, at which cavitating flow incepts:
where
P is the inlet pressure,
Pvap is the saturation vapor pressure of the working fluid,
is the fluid density, and
V is the characteristic velocity of the fluid in the microfluidic device, which is calculated at the beginning of the nozzle based on the volumetric flow rate of the system (flow rate/cross sectional area).
Flow pattern images obtained from the high speed camera during the experiments are also compared at different inlet pressures to indicate the cavitation intensity. Achieving more intensity of the generated cavitating flows is of interest. In addition, reaching developed cavitating flows at lower inlet pressures raises the efficiency in possible applications. In conclusion, more intense cavitating flows at a lower inlet pressure are in favor of this optimization.
The initially designed microfluidic device (device 0) was a plain surface micro orifice with the geometrical features shown in
Table 1. The device was installed on the chip holder sandwich (package), and the inlet pressure was gradually increased. The flow patterns were recorded simultaneously during the experiment. For the plain wall configuration, no inception of cavitation was detected over the entire inlet pressure range. The microfluidic devices are durable to inlet pressures of about 7584 kPa. This microfluidic device needed to be optimized.
3.2. Phase Two: Wall Roughness
Procedure: The inlet pressure is gradually increased and is kept constant every 350 kPa for a few minutes to make sure that the transient state of the flow pattern passes, and steady state is achieved. The inlet pressure is increased until inception of cavitating flow could be captured by the visualization system of the experimental setup. The cavitation numbers at three different inlet pressures are calculated and compared after the experiments. It should be noted that the inlet pressure is considered as the pressure term in cavitation number (Equation (4)), and the saturation vapor pressure of water at 20 °C, 2333 Pa, is the other pressure term in this equation. The velocity of the fluid is also calculated using the volumetric flow rate at the beginning of the nozzle, where the velocity is maximum.
The density was 998.2 kg/m
3 in this analysis.
Figure 5 shows the cavitation number and the inlet pressure of each microfluidic device and the corresponding flow patterns. The cavitation number in microfluidic devices first decreased and then increased in all the tests. According to Equation (4), the decrease in the value of the cavitation number is as a result of the velocity increase with the inlet pressure. This trend continues until the microfluidic channel flow is choked. Beyond this condition, the increase in the inlet pressure does not change the velocity of the working fluid. Thus, the cavitation number increases, which could be considered as the indicator of the choked flow. The microfluidic devices were tested until this trend was observed for all the cases.
As shown in
Figure 5, the extension inception (the onset of cavitating flow in the extension zone of the microfluidic devices) in device A happens at 2.06 MPa, while the same flow pattern is seen at a higher inlet pressure of 3.10 MPa in device C. The difference between these two devices is the total length of the roughness elements. According to this observation, shorter total length of the roughness elements in device A exhibits a better performance compared to device C. In addition, the nozzle inception (the inception of cavitating flow at the beginning of the nozzle) happens at 4.13 MPa, while it never happens in device C. This trend also confirms the better performance of device A compared to device C. On the other hand, the performance of device B resembles the performance of device C in the extension. However, the nozzle inception was recorded at 5.17 MPa in device B, which suggests a relatively better performance compared to device C with no nozzle inception. For the first three microfluidic devices (devices A, B, and C in
Table 1), considering the inlet pressure for inception of cavitating flows in the extension section of the microfluidic devices as the optimization goal, it can be concluded that the optimum total length of the wall roughness elements is one third of the total length of the microchannel. Furthermore, the shorter the total length of the lateral roughness elements, the better the performance acquired from the device will be. However, comparing the performances of devices D, E, and F (
Table 1) leads to a different conclusion. The difference between this group of the microfluidic devices and the previous three devices is the height of the roughness elements, which is ten times larger. Device F with the longest total length of the roughness elements exhibits the best performance among the tested microfluidic devices. The extension and nozzle inception happens at 1.86 MPa, while the inceptions in the other devices happen at 2.06 MPa and 2.48 MPa, respectively.
The results of the experimental observations could be interpreted in two scenarios. In the first scenario, the microfluidic devices are divided in two groups. The members of each group have the same roughness element height but different total lengths. In the second scenario, on the other hand, the microfluidic devices are divided in three groups, in which the lengths of the roughness elements are the same but their heights are different. Based on the first scenario, the shorter total lengths of roughness elements for the tall elements, device A, and longer total length for short elements, device F, are in favor of this optimization process (
Figure 6). On the other hand, from the second scenario perspective (roughness height), comparing devices A and D, B and E, and C and F, it can be concluded that the shorter elements lead to a better performance, and this conclusion is consistent for all three groups in this scenario. The difference between these pairs of devices lies in the height of the roughness elements. For instance, the roughness height in device C is more than in device F, and the inception happens at 3.10 MPa, while the inception inlet pressure for device F is 1.86 MPa. The other two pairs of the devices follow the same trend in their results. The experimental results for the nozzle inception are also consistent with the extension results.
After both scenarios, devices C and E could be named as the worst cases for all six microfluidic devices, which will be the interest of the next phase of optimization. The inception of cavitating flows in the extension happens at 3.10 MPa in device C, and the intensity of cavitating flows in the extension is not very high at the highest inlet pressure of 6.20 MPa (
Figure 5). In device E, on the other hand, the inception happens at 2.48 MPa, and cavitating flows are not intense in the extension at 6.20 MPa (
Figure 5).
The analysis results from the first scenario mark device A as the optimum design, which has a shorter total length (
LR) and taller roughness elements (
HR) on the wall, whereas the second scenario suggests device F as the optimum design, which has the longest total length and the shortest roughness elements. The contradiction between the analysis results shows an interaction between the design parameters. A full factorial design on the microfluidic device considering both length and height of the roughness elements shows this interaction in
Figure 7. The steep angle between the parameter lines in the interaction profile bolsters this claim.
3.3. Phase Three: Working Fluid Change to Ethanol
We aimed use the microfluidic device we were working on in an energy harvesting system. As a result, a larger number of bubbles in the extension area was desired so that the thermal energy could be harvested from the collapsing bubbles. The fabricated and tested microfluidic devices in the previous phases did not exhibit good performances in this regard. The extension areas in neither of the cases showed an intense presence of cavitating bubbles, even at high inlet pressures. As mentioned in
Section 2.2, the tensile strength of the working fluid plays a critical role in bubble generation. Equation (1) shows the tensile strength of the working fluid, and Equation (2) expresses the required energy to initiate the phase change. Decreasing the surface tension of the working fluid leads to a lower tensile strength and lower energy requirement. The surface tension of water is 72 mN/m at 25 °C, while the surface tension for ethanol is 21.78 mN/m at the same temperature. The lower surface tension of ethanol makes it a good candidate for the optimization of the microfluidic devices; 100% ethanol was used to run new sets of tests on the fabricated microfluidic devices from the previous phase of the optimization. In the second phase of the optimization, the worst performing devices in terms of the cavitation inception were devices C and E.
In the case of ethanol as the working fluid, the intensity of cavitating flows in the extension is significantly higher in comparison with the cases with water. The images from the high speed camera show that the number of the bubbles generated in the extension section is significantly higher than that of the case of water at the same inlet pressure. This observation was due to the lower surface tension of ethanol in comparison with water. Apart from the inception pressure difference between ethanol and water, the physics of extension inception are also different. The inception in experiments with water happens gradually with a weak bubble cloud in the extension, while the bubble cloud in the case of ethanol appears suddenly and intensely. In addition, the lower surface tension of ethanol increases its tendency to form bubbly flow rather than a bubble cloud, which could be seen in the experiments with water (
Figure 5).
Figure 8 shows cavitating flows in the extension of devices C and E working with ethanol.
The inception pressure in device C working with water is 3103 kPa, while this value drops to 1793 kPa for ethanol. The same trend is visible for device E with inception pressure difference of 414 kPa. The reason for this observation could be explained by flow rate expressed as Equation (5) [
28]:
where the width of inlet and nozzle are denoted as
W and
w, respectively;
H is the depth of the etched channels in the silicon substrate; and
Cd is the discharge coefficient of the device.
The flowrate, at which cavitation incepts, is formulated with the geometrical parameters and the thermophysical properties of the working fluid. The vapor pressure of ethanol at 25 °C is 7.83 kPa, and the density is 789 kg/m3, while the vapor pressure and the density of water are 3.17 kPa and 998.2 kg/m3, respectively. The higher vapor pressure of ethanol in comparison with water decreases the critical flowrate, at which cavitation incepts, which explains the earlier inception of cavitating flows.
Since the proposed energy harvesting system uses the released heat energy from the collapsing bubbles to generate electricity, ethanol could be nominated as a better working fluid. In this phase of the optimization, the working fluid was changed so that the fluid properties would be in favor of that application. Comparing the flow patterns of ethanol and water, more bubbles in the extension are visible for ethanol. On the other hand, the generated bubbly flow and larger bubbles in the extension release more energy upon collapse. Equation (6) can be employed to find the potential energy of the bubbles generated in the microfluidic device:
where
is the static pressure of the surrounding. As can be seen, the potential energy of bubbles (
) is in direct relation with the third power of the radius of the bubble. More than half of this energy is converted to heat after collapse of the bubbles, which is the target of the proposed system [
29].
Table 2 summarizes and compares the design and optimization phases in this study. The initially designed microfluidic device did not show a good performance in terms of cavitation generation. In Phase II, wall roughness elements with different heights (
HR) and lengths (
LR) were employed in the nozzle section. The change in the tensile strength of the working fluid facilitates cavitation inception in the extension. The worst cases (performances) in this phase were chosen for the last phase. Ethanol as the working fluid was nominated for Phase III so that the thermophysical differences with water could work in the favor of this optimization goal. The results indicate that the early inception of cavitating flows and its intensification can be achieved by adding lateral roughness elements to the microfluidic devices. Changing the working fluid from water to ethanol also serves for intensification of cavitating flows. Thus, more cavitating bubbles can be generated in the extension, and more energy harvesting from the collapsing bubbles could be attained.
4. Conclusions
In this study, single orifice microfluidic devices were designed and fabricated on silicon and were anodically bonded to glass. The fabrication of the devices in silicon ensures their resistance to high pressure flowing flows without having any mechanical failure or geometrical deformation. An experimental setup was designed and assembled in order to visualize cavitating flow patterns and to measure pressures at different points. The initially designed microfluidic device was equipped with six different lateral wall roughness element configurations, and the performances of the configurations were studied experimentally. Cavitation number was used as the major parameter to provide consistent comparisons during the tests. Two scenarios were considered to analyze the results. In the first scenario, the microfluidic devices were divided in two groups with the same roughness element heights but with different total lengths of roughness elements. The results showed that shorter total length leads to an earlier and more intense cavitating flows in the extension for taller roughness elements. In addition, longer total length of the lateral roughness elements increased the performance for shorter roughness elements. In the second scenario, on the other hand, the devices were divided into three groups, while the total lengths of the roughness elements were the same and their heights were different. The experimental results confirmed that shorter roughness elements exhibited a better performance regardless of the total length of the roughness elements. Since the conclusions from two scenarios did not match, a parametric effect study was also performed to examine the interactions among the design parameters.
In the last phase of optimization, the worst cases from the previous phase were selected for the tests with ethanol as the working fluid instead of water. The results showed more intense cavitating flows at low inlet pressures for these devices due to the thermophysical differences between the working fluids.
The employment of wall roughness elements inside the microchannel could facilitate cavitation flows. Although the initially designed microfluidic device did not generate cavitating bubbles at any inlet pressure, all the secondarily designed and fabricated devices led to cavitating flows in their extension. Ethanol, due to its thermophysical properties, could considerably enhance the performances of the microfluidic devices. Thus, based on the results from the present study, ethanol could be considered a strong candidate for the working fluid to enhance the intensity of cavitating flows. Further research on the effect of ethanol on other microfluidic channels with different configurations will shed more light on this finding.
It is evident that the geometry of roughness elements could significantly alter the performance of the microfluidic devices. As a result, we suggest testing different geometries as a future research direction.