Next Article in Journal
Comparison of Membrane-Based Treatment Methods for the Removal of Micro-Pollutants from Reclaimed Water
Previous Article in Journal
Influence of Seawater Salinity on the Survival, Growth, Development and Neonate Production of Scapholeberis mucronata (O. F. Müller) (Crustacea: Cladocera)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Investigation on Influence of Number of Bubbles on Laser-Induced Microjet

1
Graduate School of Engineering, Tokyo University of Agriculture and Technology, Tokyo 184-8588, Japan
2
Institute of Engineering, Tokyo University of Agriculture and Technology, Tokyo 184-8588, Japan
*
Author to whom correspondence should be addressed.
Water 2022, 14(22), 3707; https://doi.org/10.3390/w14223707
Submission received: 3 October 2022 / Revised: 7 November 2022 / Accepted: 15 November 2022 / Published: 16 November 2022
(This article belongs to the Section Hydraulics and Hydrodynamics)

Abstract

:
In this study, the impact of the number of bubbles on the velocity of laser-induced microjet is numerically investigated, focusing on the pressure wave propagation generated by multiple laser-induced bubbles. First, we show that the microjet velocity increases with the increasing impulse of the pressure wave propagating to the meniscus direction. This result indicates that it is possible to study the structure of the pressure field generated from bubbles to investigate the effect on microjet generation. In addition, it is found that the microjet is weakened with the increase in the number of bubbles. Next, we show that the propagation of the pressure waves has two types. The first type is propagating from a bubble to a meniscus. The second type is propagating round trip between nearby bubbles or by the bubble itself. Finally, we explain the reason for the decrease in the microjet velocity with the increasing number of bubbles by an expansion history of the bubbles, which depends on their interaction with the pressure waves. These results could help to design not only laser-induced microjet generation but also devices that use laser-induced bubbles generated in a microchannel.

1. Introduction

Tagawa et al. proposed a device for generating a microjet by focusing pulsed laser on a liquid-filled capillary [1], and Figure 1 shows schematics of the laser-induced microjet. This microjet has a diameter of up to 30 μm and travels at a speed of up to 850 m/s. The jet formation from the meniscus has been studied by several research groups. Antkowiak et al. studied the jet generation from the meniscus of a test tube filled with liquid by impacting the tube bottom on a hard surface [2]. Jet formation by bubble collapse on a free surface has been studied [3,4,5]. However, in these methods, compared with the method proposed by Tagawa et al. [1], the jet scale is larger, and the jet speed is about one order of magnitude smaller. High-speed microjets are very useful for practical engineering applications, such as needle-free drug injection [6,7,8,9], inkjet printing [10], and nano-droplet generation [11]. Therefore, considerable interest has been focused on microjet generation in the field of microfluidic devices.
The generation process of laser-induced microjets is as follows: a high-pressure laser-induced bubble is generated by focusing a pulsed laser in a liquid-filled capillary (Figure 1a). The pressure wave generated by the laser-induced bubble propagates to the meniscus (Figure 1b). When the pressure wave interacts with the meniscus, the focusing effect of the flow occurring at the curved meniscus causes high-speed microjet generation (Figure 1c). Note that the pressure wave propagation is not a simple process. A complex pressure field is formed inside the microtube through multiple reflections on the tube wall and other bubbles, and the microjet velocity also has complex behavior responding to the bubble conditions.
It has frequently been reported that multiple bubbles are often generated by focusing the laser in a liquid-filled capillary. Tagawa et al. experimentally investigated the jet velocity in a case in which multiple bubbles are generated. Additionally, they especially investigated the relationship between the approximated total bubble volume and jet velocity by changing the laser focusing position [1]. In their study, the bubble volume was geometrically approximated from the optical path of the laser, meaning that the volume increment is correlated with the number of bubbles. The experimental results [1] showed that the jet velocity increased with the approximated bubble volume, which indicates that the bubble volume and number of bubbles significantly affect the jet velocity. However, the mechanisms through which these parameters affect the jet velocity are not clear. To efficiently generate a microjet with less laser energy, it is necessary to investigate the effects of the laser-induced bubble volume and the number of bubbles on the jet velocity in detail and clarify the mechanism. In this study, we mainly focused on the effect of the number of bubbles.
It is important to clarify the characteristics of the pressure waves propagating from laser-induced bubbles to understand the microjet formation. Peters et al. numerically simulated laser-induced jets using the boundary integral method and showed that the summation of the pressure impulses decides the microjet velocity [12]. In addition, using a hydrophone and a high-speed camera, Hayasaka et al. experimentally showed that the jet velocity is roughly proportional to the pressure impulse [13].
There have been many studies on the pressure waves and pressure fields generated by laser-induced bubbles [14,15,16,17,18,19,20,21,22,23]. Kodama et al. investigated the pressure waves produced by an argon fluoride excimer laser, ruby laser, and shock tubes [15]. Moreover, Sankin et al. experimentally investigated the interactions between the pressure waves generated by a single laser-induced bubble and lithotripter shock waves [16]. Klaseboor et al. investigated the interaction of a single laser-induced bubble with lithotripter shock waves using the boundary element method [17]. Quinto-Su et al. investigated a cavitation bubble cloud between two laser-induced bubbles [19]. They also investigated the time evolution of the bubbles and the pressure waves in the cavitation cloud through experiments and simulations. Hsiao et al. simulated the dynamics of two bubbles and the pressure field around bubbles in a channel using a boundary element method [20]. Beig simulated two bubbles near a wall and showed the pressure field around two bubbles [22]. Kyriazis et al. numerically investigated the effects of the initial pressure of a single bubble and meniscus shape on microjet velocity [23]. However, few investigations have focused on the structure of the pressure waves propagating in the meniscus direction from multiple bubbles inside a microchannel.
In this study, we numerically simulate the microjet generation and the pressure waves propagating from multiple bubbles. First, we confirm a relationship between impulse and microjet velocity by simulating the case with and without the meniscus as in the previous studies [12,13]. In addition, we evaluate the effect of the number of bubbles on the microjet velocity and the impulse. The simulation is conducted for various numbers of bubbles with fixed total initial energy and bubble volume. Next, we investigate the detailed structure of pressure waves propagating in the meniscus direction and the pressure-time variation is explained at a specific point based on the pressure wave propagation path. Finally, we discuss the reason for the changing microjet velocity with the number of bubbles using the expansion history of the bubbles. To reduce the computational cost, we simulate the pressure wave propagation in a two-dimensional microchannel in this study. In a microtube, the focusing effect of the pressure waves occurs because the waves generated at the center of the axis are reflected on the surrounding tube wall and return to the tube axis. In the two-dimensional microfluidic channel, the focusing effect becomes weaker, and the quantitative discrepancy appears in the pressure impulse strength and jet velocity. However, the basic process of pressure field generation such as wave reflection on the wall and bubbles can be taken into account even in the two-dimensional microchannel, and the objective of this study is to make a quantitative discussion and understand the basic phenomena of pressure field generation.

2. Computational Methods

2.1. Governing Equations

The boundary integral method used in the previous study accurately predicts the microjet velocity [12]. However, we consider it unsuitable for our research purpose because pressure wave propagation needs to be investigated. In this study, the underwater pressure wave propagation is numerically simulated using the 5-equation model [24,25], which is one of the compressible multi-phase fluid models. The 5-equation model is widely used for compressive multiphase flow simulations such as bubble-shock interaction, bubble collapses, and underwater explosions [22,26,27,28,29,30]. This model can be written as follows:
α i ρ i t +   ·   α i ρ i u = 0 ,
ρ u t +   ·   ρ u u + p I = 0 ,
ρ u t +   ·   ρ u u + p I = 0 ,
α g t +   ·   u α g = α g + K   ·   u   ,
where α is the volume fraction, ρ is the density, u = u , v , w denote the velocity components (Cartesian coordinates), e is the total specific energy, p is the pressure, and the subscript i     g , l indicates gas or liquid. Equation (1) shows that the mass conservation laws are defined both for the gas and liquid phases. Here, K is
K = α g α l ρ l c l 2 ρ g c g 2 α l ρ g c g 2 + α g ρ l c l 2
where c is the speed of sound. A term K · u in Equation (4) follows from an asymptotic expansion of the total-disequilibrium model of Baer and Nunziato [31]. The 5-equation model can be corrected to capture bubble collapse accurately by adding term K · u proposed by Kapila et al. [24,32,33]. On the other hand, it is also known that the term K · u induces numerical instability in strong compression or expansion in mixture regions due to its non-conservative nature [33,34]. The following relationships are held for the physical quantities.
1 = α g + α l ,
ρ = α g ρ g + α l ρ l ,
ρ e = ρ ϵ + 1 2 ρ u 2 + v 2 + w 2 ,
ρ ϵ = α g ρ g ϵ g + α l ρ l ϵ l ,
p = p g = p l ,
u = u g = u l ,
where ϵ is the internal specific energy. The closure of the system is obtained through the specification of an equation of state (EoS). To account for fluid properties, various EoS have been proposed considering molecular effects (See Métayer and Saurel [35]). To close the system, we introduced the Stiffened gas EoS (SG EoS):
p i = γ i 1 ρ i ϵ i γ i p r e f , i ,
where γ is the specific heat ratio, and p r e f is the reference pressure. By assuming that the pressures of each phase are equal at the interface, Equation (12) can be replaced by [22,25,36]
p = ρ ϵ η ξ ,
ξ = i g , l α i γ i 1 ,
η = i g , l α i γ i p r e f , i γ i 1 .
All the physical quantities were nondimensionalized as follows:
x = x ˜ R ,
ρ = ρ ˜ ρ r e f ,
u = u ˜ c r e f ,
p = p ˜ ρ r e f c r e f 2 ,
ρ e = ρ ˜ e ˜ ρ r e f c r e f 2 ,
where x is the position vector. The tilde (~) indicates a dimensional quantity. The reference values are set as R = 35   μ m , ρ r e f = 1.204   kg / m 3 , and c r e f = 341.0   m / s .
In this study, we disregard the viscosity, cavitation, and surface tension. The previous study by Peter et al. also ignored cavitation and viscosity, but fortunately, it was reported that the microjet velocity was successfully reproduced accurately compared to the experiment [12]. In addition, Hayasaka et al. showed that the surface tension did not hinder the microjet formation [13]. On the other hand, the surface tension and viscosity may have some quantitative influence on the expansion motion of bubbles, and cavitation bubbles may have some impact on pressure wave propagation. However, in this study, the pressure propagation is analyzed and discussed in the simple condition and settings to obtain basic knowledge.

2.2. Simulation Settings and Conditions

In this study, we numerically simulate the pressure wave propagation from multiple high-pressure bubbles (the number of bubbles N     1 , 2 , 3 ) inside a two-dimensional microchannel. To evaluate microjet velocity and impulse, we set up two simulation conditions.
Figure 2 shows the simulation settings of a microchannel with the meniscus (w-meniscus setting) for evaluating the microjet velocity. To capture the bubble interface and meniscus, the grid spacing is set to 0.04 in the range from   x 6.2 ,   30 . In this case, we do not apply the term K · u   due to strong expansion in the meniscus. In some studies simulating an underwater explosion with strong expansion waves, it was reported that the term K · u should not be applied for stable computation [27,30].
Figure 3 shows the simulation settings for a microchannel without the meniscus (w/o-meniscus setting) for evaluating impulse. The computational grid of the w/o-meniscus setting is uniform in the range from x 12 ,   12 ,   and the grid spacing is 0.04. In the regions of x < 12 and x < 12 , the computational grid is proportionally stretched in the directions of the left and right boundaries. In this study, the pressure wave impulse is evaluated at the inspection surface. The inspection surface for investigating the pressure impulse is the cross-section of the tube at 12 from the bubble center in the x-direction. The grid dependence is identified in Appendix B.
To measure the time history of the pressure, the pressure probe is placed at (12, 0). In this case, we evaluate the impulse as the following equation.
I = 0 T S p d s d t = Δ t 0 ,   T j cell   indexes p j s j t ,
where t is the time interval, j is the index of the cell on the inspection surface, s is the area of the cell face.
The bubble radius is set to 1 when the bubble number is 1. The parameters of the SG EoS is set as in Equation (22) following Shukla et al. [26] in the nondimensional value. The initial conditions are set as in Equation (23).
γ g ,   γ l ,   p r e f , g ,   p r e f , l = 1.4 ,   6.12 ,   0 ,   2450 ,
ρ g , ρ l , u , v , w , α g , p = 83.3 ,   0 ,   0 ,   0 ,   0 ,   1 ,   571.4 Laser-induced bubble 0 ,   833 ,   0 ,   0 ,   0 ,   0 ,   0.7143 Liquid 1 ,   0 ,   0 ,   0 ,   0 ,   1 ,   0.7143    Air .
These conditions are the same for all the computational cases. The bubble pressure is based on the laser energy of 36 [ μ J], which is in the same order of magnitude as that of the pulsed-laser devices frequently used in the experiments [6,7,8,9,13,37]. To keep the total internal energy of the bubbles constant among the three cases, the sum of the bubble volumes is set to be constant. From this restriction, the radius of the bubble r is determined using the following equation:
r = R 2 N .
The bubble radius complies with the typical spot size of a pulsed laser [37]. The bubbles are placed at equidistant intervals, as shown in Table 1.
Usually, a bubble interface is smeared for computational stability [22,28,29,32,33]. The smeared volume fraction is given a hyperbolic tangent function [32].
α g = 1 2 1 tanh l x r h ,
where x is the position of the cell center, and l x is the distance of the bubble center. The variable h is the characteristic length of the corresponding computational cell.
Only the lower side of the microchannel is computed by the implemented symmetric boundary condition. The slip condition is implemented to the microchannel wall, and the outflow boundary condition using zero-order extrapolation is implemented at the left and right sides. The outflow boundary is set far enough away to stretch the grid so that the reflections are not a problem.
The time interval is fixed at t = 1.0 × 10 4   [-] from the Courant–Friedrichs–Lewy condition, and the Courant number was set as 0.2.

2.3. Numerical Schemes

The 5-equation model is discretized using the finite volume method. The numerical fluxes are computed using the Harten–Lax–van Leer contact approximate Riemann solver extended for the 5-equation model [38], and we use the second-order total variation diminishing Runge–Kutta for the time stepping. The following reconstruction methods are used depending on the settings for reasonable computation.
  • ρ - THINC-MUSCL3
To stabilize the simulation of the w-meniscus setting, we use the third-order monotonic upstream-centered scheme for conservation laws [39] with the minmod limiter [40] (MUSCL3) and   ρ -THINC [29]. ρ - THINC is one of the variations of the tangent of the hyperbola interface capturing method (THINC) [41] to sharpen the volume fraction and density on the interfaces. This coupling reconstruction method ( ρ -THINC-MUSCL3) captures the interface and stabilizes the simulation.
  • MP-WENO5-JS
To accurately capture the pressure waves propagating in the w/o-meniscus setting, we use the five-order monotonicity preserving weighted essentially non-oscillatory schemes with the indicator developed by Jiang and Shu (MP-WENO5-JS) [42,43,44] for a reconstruction.
The validation of the computational solver is described in Appendix A.

3. Result and Discussion

3.1. Microjet Evolution and Time History of Impulse

In this section, we show the relationship between impulse and microjet velocity. In addition, we also discuss the time evolution of microjets and the time history of the impulse. In this case, we used the ρ -THINC-MUSCL3 reconstruction to stabilize the simulations.
Figure 4 shows the time evolution of the microjet by the single bubble case. As reported by Tagawa et al. [1], it is confirmed that the microjet is generated from the center of the meniscus. Figure 5 shows the time history of the microjet velocity and microjet formations at the final state. Figure 5a shows that the microjet velocity at the final phase decreases as the number of bubbles increases. In addition, the final position of the microjet head also differs with the difference in the time history of the microjet velocity (Figure 5b–d).
Next, the impulse acting on the inspection surface is discussed. Figure 6 shows the time history of the impulse calculated by Equation (21). Figure 6 shows that the case with the three bubbles has the highest impulse until the time of 3.5. After that, the case with the two bubbles has the highest impulse in the time range from 3.5 to 4.5. After t = 4.5 , the case with the single bubble comes to have the highest impulse. As a summary of the above-mentioned results, the impulse increases with more multiple bubbles in the initial stage. However, after that, the impulse with the single bubble comes to exceed that with more bubbles.
As a summary of the results in this section, the pressure impulse at the final phase decreases as the number of bubbles increases. This tendency is the same as in the microjet velocity. Thus, we conclude that the microjet velocity has a strong correlation with the impulse of pressure propagating from the bubble to the meniscus. These results indicate that it is possible to discuss the effect of bubble conditions on the microjet velocity by investigating the structure of pressure field in the w/o-meniscus setting.

3.2. Pressure Field Evolution around a Single Bubble

In this section, we investigate the structure of the pressure field generated from the single bubble in the w/o-meniscus setting. We used the MP-WENO5-JS reconstruction to accurately capture the pressure waves.
Figure 7 shows the pressure field evolution generated by a single laser-induced bubble. At first, a spherical pressure wave (compression wave) is generated by a high-pressure laser-induced bubble (Figure 7a) and propagates in the liquid (Figure 7b). The compression wave then reaches the microchannel wall and is reflected as a compression wave (Figure 7c). Next, the pressure wave from the wall is reflected at the bubble interface and becomes an expansion wave (Figure 7d). After that, the expansion wave is reflected at the wall (Figure 7e) and is also reflected as a compression wave at the bubble interface (Figure 7f). As a result of this series of wave reflections, a complex pressure field is formed inside the microchannel. In addition, the pressure waves propagating repeatedly between the bubble and wall can be observed.
Figure 8 shows the time history of the pressure at the probe position. As observed in the figure, the time history of the pressure has periodic oscillations with the period T . In addition, the pressure peaks are confirmed at a unique timing t n in the initial phase.
It can be considered that the time-pressure history is generated by the following process: Figure 9 shows the path of the pressure wave from the bubble to the pressure probe. There are two types of dominant pressure wave paths reaching the pressure probe. The first path is by the reflected pressure wave not only on the wall but also on the bubble (Figure 9a). The pressure wave generated by the high-pressure bubble returns to the bubble through the reflection at the wall and then propagates to the pressure probe through the reflection on the bubble. The pressure wave then repeatedly goes back and forth between the bubble and wall. Therefore, the pressure wave propagation to the pressure probe is also periodically repeated. The second path is by the pressure wave propagation from the bubble to the probe through multiple reflections at the wall (Figure 9b). The pressure wave generated by the high-pressure bubble reaches the probe directly or through more than one reflection on the wall. It is considered that the periodical oscillation in the time-pressure history shown in Figure 8 is due to the first pressure path and that the pressure peaks seen in the initial phase come from the second pressure path, from which the pressure wave arrives at the probe through a variety of paths. This results in the delay of the pressure peak depending on the propagation paths.
To verify the above-mentioned pressure oscillation mechanism, first, the periodic time of the pressure oscillation described in Figure 9a is calculated using a simple theory. The period of the pressure peaks T   is expected to be calculated from the time taken by the pressure wave propagating over the distance L p , which is the round-trip distance between the bubble and wall, as follows:
T = L p c l = D c l
where c l is the sound speed of the liquid, and D = 14 is the width of the channel. In this study, c l is 4.8 by a dimensionless quantity. Next, the delay time of the pressure peak coming from the second pressure wave path is estimated by considering the length of the pressure path L n , as depicted by the dashed line in Figure 9b. The delay time t n for the pressure wave to arrive at the probe could be calculated using the following equation:
t 1 = L 1 c l = L c l t 2 = L 2 c l = L 2 + D 2 c l t 3 = L 3 c l = L 2 + 2 D 2 c l t n = L 2 + ( n 1 D 2 c l
where L   is the length between the bubble and the inspection surface. The results of the calculations of Equations (26) and (27) are summarized in Table 2. The results in Equations (26) and (27) are consistent with the simulation results. However, the equations evaluate a longer time than that of the simulation results. The sound speed in the liquid is changed in the range from 4.8 to 5.0 by the pressure wave propagating from the laser-induced bubble. The pressure wave propagated faster than the speed of sound considered in the equations. In addition, the bubble radius is not considered in Equations (26) and (27). Therefore, the equations evaluated a longer time.
Finally, we summarize the structure of pressure field generated from the single laser-induced bubble in a microchannel. We can see two types of pressure wave propagation. The first type is propagation from the bubble to the meniscus, and the second type is reciprocating propagation around the bubble. These two types of pressure waves generate the structure of the pressure field.

3.3. Pressure Field Evolution around Multi-Bubbles

In this section, the results in the case of multiple bubbles are discussed. The MP-WENO5-JS reconstruction and w/o-meniscus setting are used as in Section 3.2.
Figure 10 shows a comparison of the time evolution of the pressure field among the single bubble, two bubbles and three bubbles. Figure 10a shows that a pressure wave overlap occurred more quickly with increase in the number of bubbles. When the pressure waves overlapped, a strong pressure area is generated. As shown in Figure 10b, a pressure area higher than p = 226 is generated by the overlap of two pressure waves in the case of the two bubbles, and a stronger pressure area over p = 455 is generated by the overlap of three pressure waves in the case of the three bubbles. However, in the case of the single bubble, the pressure waves did not overlap, so the maximum pressure in Figure 10c is about p = 161. In this manner, with the increase in the number of bubbles, multiple pressure waves overlap more quickly, and as a result of the quick pressure-wave overlap, stronger pressure waves arrive at the inspection surface more quickly.
Figure 11 shows the time history of the pressure at the probe. As shown in Figure 11, the pressure oscillation with period T can be seen even in the case of multiple bubbles. Note that period T differs depending on the number of bubbles.
Equations (26) and (27) can be extended to multiple bubbles. Figure 12 shows the pressure wave path to the probe in the case of multiple bubbles. First, the equation for the period T (Equation (26)) is extended. In the case of the multiple bubbles, the round-trip distance of the pressure wave in Equation (26) is turned into the distance between the bubbles or twice the distance between the bubble and the wall. The propagation distance L p of the pressure wave can be calculated using the following equation:
L p = D N
and the period T can be calculated using L p through the following equation:
T = L p c l = D N c l .
Next, the equation for the delay time (Equation (27)) is extended. As in the case of the single bubble, the dashed line in Figure 12b is the path on which the pressure wave propagates from the bubble to the probe. The delay time t n can be calculated using the following equations:
t 1 = L 2 + s 2 c l t 2 = L 2 + D s 2 c l t 3 = L 2 + D + s 2 c l t n = L 2 + s 2 c l   If   n = 1 L 2 + n 1 D s 2 c l   If   n   is   an   even   number , L 2 + n 2 D + s 2 c l   If   n   is   an   odd   number ,
where s is the distance from the channel center to the bubble. The results of the calculations of Equation (29) and (30) for multiple bubbles are shown in Table 3 and Table 4. For the delay time, Equation (30) is calculated for each bubble, and the three least times are shown. There is a little discrepancy between the analytical and simulation results due to not considering the changing sound speed and the bubble radius. However, the error is less than 10%, and the analytical results are consistent with the simulation results. These results indicate that we can see two types of pressure wave propagation (propagating from the bubble to the meniscus and reciprocating around a bubble) even with multiple bubbles.
Table 5 shows a comparison of the maximum pressure values and the arrival times of the maximum pressure peak at the probe. The maximum pressure increases with the increase in the number of bubbles. Moreover, the maximum pressure arrival time becomes faster with the increase in the number of bubbles. This result is caused by the larger number of overlapping pressure waves when the number of bubbles increases, as confirmed in Figure 10. In Figure 6, we can see that the impulse increases with increasing bubbles in the initial stage. The increased impulse with a larger number of bubbles at the initial stage is concluded to be due to the pressure wave overlap effect.

3.4. Discussion

All figures and tables should be cited in the main text as Figure 1, Table 1, etc.
The force applied to the surrounding fluid by the expansion of the bubbles is expected to turn into an impulse. In this section, based on the time history of the total bubble volume, the discussion focuses on the reason why the impulse increased with the decrease in the number of bubbles (Figure 13).
Figure 13 shows that the case with the three bubbles had the largest total volume at t 0.0 ,   1.2 and that the case with the two bubbles then had the highest total volume at t 1.2 ,   2.4 . However, the volume in the case with the single bubble became the largest eventually. The time history of the bubble volume shows the similar trend to the impulse-time history shown in Figure 6. Figure 13 b shows a zoom-up view of the timing with the change in the expansion gradient of the bubble. The timing of the change in the bubble expansion gradient coincides with the timing at which the bubble first interacted with the compression pressure wave, which inhibited the expansion of the bubble. Table 6 summarizes the times when the bubble first contacted the compression pressure wave. In the setup of this study, the interval between the bubbles become shorter with the increase in the number of bubbles. Therefore, as the number of bubbles increases, the interference between the pressure waves and bubbles took place more quickly, which resulted in a smaller bubble expansion, i.e., the impulse decreased.
Tagawa et al. conducted experiments varying the objective lens magnification used to focus the laser with constant laser pulse energy and confirmed that the number of bubbles generated initially reduces with higher lens magnification, which leads to higher pressure impulse [45]. The results of this experiment can be interpreted as the generation of a larger impulse with smaller number of bubbles, and it is consistent with the results obtained in this study. This study clearly shows that the larger pressure impulse with the smaller number of bubbles can be explained to be due to the inhibition of the bubble expansion in the larger number of bubbles, as shown in Figure 13.

4. Conclusions

In this study, the propagation of the pressure waves generated by multiple laser-induced bubbles inside a microchannel filled with liquid was numerically simulated using the 5-equation model. First, we investigated a proportional relationship between impulse and microjet velocity. The result showed that the final pressure impulse and microjet velocity decreased with increasing the number of bubbles. It was confirmed that the pressure impulse has a strong correlation with the jet velocity; therefore, it was concluded that the effect of bubble conditions on the microjet velocity could be addressed by investigating the structure of the pressure field even without a meniscus. Next, we investigated the time evolution of the pressure field and considered the path of the pressure wave in the case with a single bubble. The results showed that the propagation of the pressure waves has two types: from a bubble to a meniscus and reciprocating around a bubble. These two types of propagation were observed also in the case of multiple bubbles. Furthermore, by focusing on the pressure peak observed at the probe position, the maximum pressure increased with the increase in the number of bubbles, and the time of reaching the peak pressure decreased. Finally, the time history of bubble volume was investigated. The reason for the decrease in the impulse with the increase in the number of bubbles can be explained by the expansion history of the bubbles. The increase in the number of bubbles inhibited the early expansion of the bubbles and consequently led to the impulse reduction. It was suggested that the interaction between the bubbles and the pressure wave caused by the placement of multiple bubbles negatively affected the impulse at the inspection surface.

Author Contributions

T.I.: software, investigation, formal analysis, visualization, data curation, writing (original draft). H.N.: software, validation, resources, funding acquisition, project administration, writing (review and editing). Y.T.: conceptualization, supervision, writing (review and editing). All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by AMED under Grant Number JP22he0422016j0001, Japan Society for the Promotion of Science (Grant No. 20H00223, 20H00222, and 20K20972), Japan Science and Technology Agency PRESTO (Grant No. JPMJPR21O5).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors have no conflict to disclose.

Appendix A. Solver Varidation

We conducted the solver validation using a one-dimensional gas–liquid Riemann problem adopted by Cocchi et al. [46]. The initial conditions are given by
ρ l , ρ g , u , p = 0 ,   1.241 ,   0 ,   2.753 x 1 ,   0 0.991 ,   0 ,   0 ,   3.059 × 10 4 x 0 , 1
α l = 1 1 + exp x 0.5 h 1
where h is cell spacing and x is the cell center position. The properties of SG EoS are set to γ l , p r e f , l , γ g ,   p r e f , g = 5.5 ,   1.505 ,   1.4 ,   0 . It is evolved with the time step Δ t = 2.5 × 10 4 , and the constant grid spacing h = 0.0025 following Garrick et al. [47]. In Figure A1, the results show that the numerical solutions correctly identify with the exact one obtained analytically.
Figure A1. Result of Gas–liquid Riemann problem. The exact and numerical solutions are compared at t = 0.2 .
Figure A1. Result of Gas–liquid Riemann problem. The exact and numerical solutions are compared at t = 0.2 .
Water 14 03707 g0a1

Appendix B. Grid Independence Study

The grid independence was investigated using ρ -THINC-MUSCL3, which is the lower spatial accuracy reconstruction method in this study. We used the w/o-meniscus setting and compared the impulses in each grid width h = 0.028   (fine grid), h = 0.04   (intermediate grid), and h = 0.056 (coarse grid). Figure A2 shows that the coarse grid differs significantly from other grids. On the other hand, the difference between the fine grid and the middle grid is 9.1% and, therefore, we selected the intermediate grid.
Figure A2. Result of the grid independence study using ρ -THINC-MUSCL3 on w/o-meniscus setting.
Figure A2. Result of the grid independence study using ρ -THINC-MUSCL3 on w/o-meniscus setting.
Water 14 03707 g0a2

References

  1. Tagawa, Y.; Oudalov, M.; Visser, C.W.; Peters, I.R.; Meer, D.; Sun, C.; Prosperetti, A.; Lohse, D. Highly focused supersonic microjets. Phys. Rev. X 2012, 2, 31002. [Google Scholar] [CrossRef] [Green Version]
  2. Antkowiak, A.; Bremond, N.; Le Dizès, S.; Villermaux, E. Short-term dynamics of a density interface following an impact. J. Fluid Mech. 2007, 577, 241–250. [Google Scholar] [CrossRef] [Green Version]
  3. Zeff, B.W.; Kleber, B.; Fineberg, J.; Lathrop, D.P. Singularity dynamics in curvature collapse and jet eruption on a fluid surface. Nature 2000, 403, 401–404. [Google Scholar] [CrossRef]
  4. Duchemin, L.; Popinet, S.; Josserand, C.; Zaleski, S. Jet formation in bubbles bursting at a free surface. Phys. Fluids 2002, 14, 3000. [Google Scholar] [CrossRef]
  5. Bergmann, R.; van der Meer, D.; Gekle, S.; Lohse, D. Controlled impact of a disk on a water surface: Cavity dynamics. J. Fluid Mech. 2009, 633, 381–409. [Google Scholar] [CrossRef] [Green Version]
  6. Tagawa, Y.; Oudalov, N.; El Ghalbzouri, A.; Sun, C.; Lohse, D. Needle-free injection into skin and soft matter with highly focused microjets. Lab Chip 2013, 13, 1357–1363. [Google Scholar] [CrossRef] [Green Version]
  7. Kiyama, A.; Endo, N.; Kawamoto, S.; Katsuta, C.; Oida, K.; Tanaka, A.; Tagawa, Y. Visualization of penetration of a high-speed focused microjet into gel and animal skin. J. Vis. 2019, 22, 449–457. [Google Scholar] [CrossRef] [Green Version]
  8. Krizek, J.; Delrot, P.; Moser, C. Repetitive regime of highly focused liquid microjets for needle-free injection. Sci. Rep. 2020, 10, 5067. [Google Scholar] [CrossRef] [Green Version]
  9. Miyazaki, Y.; Usawa, M.; Kawai, S.; Yee, J.; Muto, M.; Tagawa, Y. Dynamic mechanical interaction between injection liquid and human tissue simulant induced by needle-free injection of a highly focused microjet. Sci. Rep. 2021, 11, 14544. [Google Scholar] [CrossRef]
  10. Delrot, P.; Modestino, M.A.; Gallaire, F.; Psaltis, D.; Moser, C. Inkjet printing of viscous monodisperse microdroplets by laser-induced flow focusing. Phys. Rev. Appl. 2016, 6, 24003. [Google Scholar] [CrossRef]
  11. Xiong, S.; Chin, L.K.; Ando, K.; Tandiono, T.; Liu, A.Q.; Ohl, C.D. Droplet generation via a single bubble transformation in a nanofluidic channel. Lab Chip 2015, 15, 1451–1457. [Google Scholar] [CrossRef] [PubMed]
  12. Peter, I.R.; Tagawa, Y.; Oudalov, N.; Sun, C.; Prosperetti, A.; Lohse, D.; van der Meer, D. Highly focused supersonic microjets: Numerical simulations. J. Fluid Mech. 2013, 719, 587–605. [Google Scholar] [CrossRef] [Green Version]
  13. Hayasaka, K.; Kiyama, A.; Tagawa, Y. Effects of pressure impulse and peak pressure of a shock wave on microjet velocity in a microchannel. Microfluid Nanofluid 2017, 21, 116. [Google Scholar] [CrossRef] [Green Version]
  14. Vogel, A.; Noack, J.; Nahen, K.; Theisen, D.; Busch, S.; Parliz, U.; Hammer, D.X.; Noojin, G.D.; Rockwell, B.A.; Birngruber, R. Energy balance of optical breakdown in water at nanosecond to femtosecond time scales. Appl. Phys. B 1999, 68, 271–280. [Google Scholar] [CrossRef]
  15. Kodama, T.; Hamblin, M.R.; Doukas, A.G. Cytoplasmic molecular delivery with shock waves: Importance of impulse. Biophys. J. 2000, 79, 1821–1832. [Google Scholar] [CrossRef] [Green Version]
  16. Sankin, G.N.; Simmons, W.N.; Zhu, S.L.; Zhong, P. Shock wave interaction with laser-generated single bubbles. Phys. Rev. Lett. 2005, 95, 34501. [Google Scholar] [CrossRef]
  17. Klaseboor, E.; Fong, S.W.; Turangan, C.K.; Khoo, B.C.; Szeri, A.J.; Calvisi, M.L.; Sankin, G.N.; Zhong, P. Interaction of lithotripter shockwaves with single inertial cavitation bubbles. J. Fluid Mech. 2007, 593, 33–56. [Google Scholar] [CrossRef]
  18. Ando, K.; Liu, A.-Q.; Ohl, C.-D. Homogeneous nucleation in water in microfluidic channels. Phys. Rev. Lett. 2012, 109, 44501. [Google Scholar] [CrossRef]
  19. Quinto-Su, P.A.; Ando, K. Nucleating bubble clouds with a pair of laser-induced shocks and bubbles. J. Fluid Mech. 2013, 733, R3. [Google Scholar] [CrossRef]
  20. Hsiao, C.-T.; Choi, J.-K.; Singh, S.; Chahine, G.L.; Hay, T.A.; Ilinskii, Y.A.; Zabolotskaya, E.A.; Hamilton, M.F.; Sankin, G.; Yuan, F.; et al. Modelling single- and tandem-bubble dynamics between two parallel plates for biomedical applications. J. Fluid Mech. 2013, 716, 137–170. [Google Scholar] [CrossRef]
  21. Supponen, O.; Obreschkow, D.; Kobel, P.; Tinguely, M.; Dorsaz, N.; Farhat, M. Shock waves from non-spherical cavitation bubbles. Phys. Rev. Fluids 2017, 2, 93601. [Google Scholar] [CrossRef] [Green Version]
  22. Beig, S.A. A Computational Study of the Inertial Collapse of Gas Bubbles Near a Rigid Surface. Ph.D. Thesis, University of Michigan, Ann Arbor, MI, USA, 2018. [Google Scholar]
  23. Kyriazis, N.; Koukouvinis, P.; Gavaises, M. Numerical investigations on bubble-induced jetting and shock wave focusing: Application on a needle-free injection. Proc. R. Soc. A Math. Phys. Eng. Sci. 2019, 475, 20180548. [Google Scholar] [CrossRef] [PubMed]
  24. Kapila, A.K.; Menikoff, R.; Bdzil, J.B.; Son, S.F.; Stewart, D.S. Two-phase modeling of deflagration-to-detonation transition in granular materials: Reduced equations. Phys. Fluids 2001, 13, 3002–3024. [Google Scholar] [CrossRef]
  25. Allaire, G.; Clerc, S.; Kokh, S. A five-equation model for the simulation of interfaces between compressible fluids. J. Comput. Phys. 2002, 181, 577–616. [Google Scholar] [CrossRef] [Green Version]
  26. Shukla, R.K.; Pantano, C.; Freund, J.B. An interface capturing method for the simulation of multi-phase compressible flows. J. Comput. Phys. 2010, 229, 7411–7439. [Google Scholar] [CrossRef]
  27. Daramizadeh, A.; Ansari, M.R. Numerical simulation of underwater explosion near air–water free surface using a five-equation reduced model. Ocean Eng. 2015, 110, 25–35. [Google Scholar] [CrossRef]
  28. Rasthofer, U.; Wermelinger, F.; Hadijdoukas, P.; Koumoutsakos, P. Large Scale Simulation of Cloud Cavitation Collapse. Procedia Comput. Sci. 2017, 108, 1763–1772. [Google Scholar] [CrossRef]
  29. Garrick, D.P.; Hagen, W.A.; Regele, J.D. An interface capturing scheme for modeling atomization in compressible flows. J. Comput. Phys. 2017, 344, 260–280. [Google Scholar] [CrossRef]
  30. Yu, J.; Liu, J.; He, B.; Li, H.; Xie, T.; Pei, D. Numerical research of water jet characteristics in underwater explosion based on compressible multicomponent flows. Ocean Eng. 2021, 242, 110135. [Google Scholar] [CrossRef]
  31. Baer, M.R.; Nunziato, J.W. A two-phase mixture theory for the deflagration-to-detonation transition (ddt) in reactive granular materials. Int. J. Multiph. Flow 1986, 12, 861–889. [Google Scholar] [CrossRef]
  32. Tiwari, A.; Freund, J.B.; Pantano, C. A diffuse interface model with immiscibility preservation. J. Comput. Phys. 2013, 252, 290–309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Schmidmayer, K.; Bryngelson, S.H.; Colonius, T. An assessment of multicomponent flow models and interface capturing schemes for spherical bubble dynamics. J. Comput. Phys. 2020, 402, 109080. [Google Scholar] [CrossRef] [Green Version]
  34. Saurel, R.; Petitpas, F.; Berry, R.A. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixtures. J. Comput. Phys. 2009, 228, 1678–1712. [Google Scholar] [CrossRef]
  35. Métayer, O.L.; Saurel, R. The Noble-Abel Stiffened-Gas equation of state. Phys. Fluids 2016, 28, 46102. [Google Scholar] [CrossRef] [Green Version]
  36. Saurel, R.; Métayer, O.L.; Massoni, J.; Gavrilyuk, S. Shock jump relations for multiphase mixtures with stiff mechanical relaxation. Shock Waves 2007, 16, 209–232. [Google Scholar] [CrossRef]
  37. Gutiérrez-Hernández, U.J.; Colle, F.D.; Ohl, C.-D.; Quinto-Su, P.A. Transient time-delay focusing of shock waves in thin liquids. J. Fluid. Mech. 2021, 910, A27. [Google Scholar] [CrossRef]
  38. Tian, B.; Toro, E.F.; Castro, C.E. A path-conservative method for a five-equation model of two-phase flow with an HLLC-type Riemann solver. Comput. Fluids 2011, 46, 122–132. [Google Scholar] [CrossRef]
  39. van Leer, B. Towards the ultimate conservative difference scheme. V. A second-order sequel to Godunov’s method. J. Comput. Phys. 1979, 32, 101–136. [Google Scholar] [CrossRef]
  40. Roe, P.L. Characteristic-based schemes for the Euler equations. Annu. Rev. Fluid Mech. 1986, 18, 337–365. [Google Scholar] [CrossRef]
  41. Skyue, K.M.; Xiao, F. An Eulerian interface sharpening algorithm for compressible two-phase flow: The algebraic THINC approach. J. Comput. Phys. 2014, 268, 326–354. [Google Scholar]
  42. Jiang, G.-S.; Shu, C.-W. Efficient Implementation of Weighted ENO Schemes. J. Comput. Phys. 1996, 126, 202–228. [Google Scholar] [CrossRef]
  43. Balsara, D.S.; Shu, C.-W. Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy. J. Comput. Phys. 2000, 160, 405–452. [Google Scholar] [CrossRef] [Green Version]
  44. Coralic, V.; Colonius, T. Finite-volume WENO scheme for viscous compressible multicomponent flows. J. Comput. Phys. 2014, 274, 95–121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Tagawa, Y.; Yamamoto, S.; Hayasaka, K.; Kameda, M. On pressure impulse of a laser-induced underwater shock wave. J. Fluid Mech. 2016, 808, 5–18. [Google Scholar] [CrossRef] [Green Version]
  46. Cocchi, J.P.; Saurel, R.; Loraud, J.C. Treatment of interface problems with Godunov-type schemes. Shock Waves 1996, 5, 347–357. [Google Scholar] [CrossRef]
  47. Garrick, D.P.; Owkes, M.; Regele, J.D. A finite-volume HLLC-based scheme for compressible interfacial flows with surface tnsion. J. Comput. Phys. 2007, 339, 46–67. [Google Scholar] [CrossRef]
Figure 1. Illustration of the generation of microjets by focusing a pulsed laser in liquid-filled capillary. (a) Laser-induced bubbles is created by a focused laser pulse. (b) High-pressure laser-induced bubbles induce pressure waves toward the meniscus. (c) Microjet is generated due to kinematic focusing on the meniscus.
Figure 1. Illustration of the generation of microjets by focusing a pulsed laser in liquid-filled capillary. (a) Laser-induced bubbles is created by a focused laser pulse. (b) High-pressure laser-induced bubbles induce pressure waves toward the meniscus. (c) Microjet is generated due to kinematic focusing on the meniscus.
Water 14 03707 g001
Figure 2. Simulation settings for the microchannel with meniscus (w-meniscus setting). This figure shows the setting for the case of two bubbles.
Figure 2. Simulation settings for the microchannel with meniscus (w-meniscus setting). This figure shows the setting for the case of two bubbles.
Water 14 03707 g002
Figure 3. Simulation settings for the microchannel without meniscus and multiple laser-induced bubbles (w/o-meniscus setting). This figure shows the setting for the case of two bubbles.
Figure 3. Simulation settings for the microchannel without meniscus and multiple laser-induced bubbles (w/o-meniscus setting). This figure shows the setting for the case of two bubbles.
Water 14 03707 g003
Figure 4. Time evolution of microjet by the single bubble. Interfaces are shown by the density gradient.
Figure 4. Time evolution of microjet by the single bubble. Interfaces are shown by the density gradient.
Water 14 03707 g004
Figure 5. Time history of microjet velocity and microjet formations at the final state. The orange dashed line in (bd) shows the position of the microjet head in the single bubble case. Interfaces are shown by the density gradient.
Figure 5. Time history of microjet velocity and microjet formations at the final state. The orange dashed line in (bd) shows the position of the microjet head in the single bubble case. Interfaces are shown by the density gradient.
Water 14 03707 g005
Figure 6. Time history of the impulse. This result is obtained from the simulation result of the w/o-meniscus setting.
Figure 6. Time history of the impulse. This result is obtained from the simulation result of the w/o-meniscus setting.
Water 14 03707 g006
Figure 7. Pressure field evolution generated by the single laser-induced bubble at (a) t = 0.0, (b) t = 0.9, (c) t = 1.8, (d) t = 3.3, (e) t = 4.5, and (f) t = 5.4; numerical schlieren at the upper side and pressure distribution at the lower side of each figure.
Figure 7. Pressure field evolution generated by the single laser-induced bubble at (a) t = 0.0, (b) t = 0.9, (c) t = 1.8, (d) t = 3.3, (e) t = 4.5, and (f) t = 5.4; numerical schlieren at the upper side and pressure distribution at the lower side of each figure.
Water 14 03707 g007aWater 14 03707 g007b
Figure 8. Time history of the pressure at the pressure probe in the case of the single bubble. The blue dashed line shows the arrival time of the pressure wave.
Figure 8. Time history of the pressure at the pressure probe in the case of the single bubble. The blue dashed line shows the arrival time of the pressure wave.
Water 14 03707 g008
Figure 9. Schematics of the path of the pressure wave from the bubble to the pressure probe point. (a) Path of the pressure wave through the wave reflection between the bubble and wall; (b) Path of the pressure wave to the pressure probe through the wave reflection at the wall.
Figure 9. Schematics of the path of the pressure wave from the bubble to the pressure probe point. (a) Path of the pressure wave through the wave reflection between the bubble and wall; (b) Path of the pressure wave to the pressure probe through the wave reflection at the wall.
Water 14 03707 g009
Figure 10. Time evolution of the numerical schlieren at the upper side and the pressure distribution at the lower side at (a) t = 0.25, (b) t = 1.0, and (c) t = 2.25. The red dash line shows the inspection surface.
Figure 10. Time evolution of the numerical schlieren at the upper side and the pressure distribution at the lower side at (a) t = 0.25, (b) t = 1.0, and (c) t = 2.25. The red dash line shows the inspection surface.
Water 14 03707 g010aWater 14 03707 g010b
Figure 11. Pressure-time history measured at the pressure probe.
Figure 11. Pressure-time history measured at the pressure probe.
Water 14 03707 g011
Figure 12. Path of the pressure wave from the bubble in the case of multiple bubbles. The path of the pressure wave is the same as the length shown by the dashed line. (a) Path of the pressure wave between the bubbles; (b) Path of the pressure wave to the inspection surface.
Figure 12. Path of the pressure wave from the bubble in the case of multiple bubbles. The path of the pressure wave is the same as the length shown by the dashed line. (a) Path of the pressure wave between the bubbles; (b) Path of the pressure wave to the inspection surface.
Water 14 03707 g012
Figure 13. Time history of the total bubble volume.
Figure 13. Time history of the total bubble volume.
Water 14 03707 g013
Table 1. Bubble positions of each case.
Table 1. Bubble positions of each case.
Number of Bubbles123
Positions x , y [-] 0 , 0 0 , 3.5 ,   0 , 3.5 0 , 4.67 ,   0 , 0 ,   0 , 4.67
Table 2. Comparison between the simulation and analytical results.
Table 2. Comparison between the simulation and analytical results.
Times [-] T t 1 t 2 t 3
Simulation Results (Figure 4)2.92.63.76.2
Results of Equations (26) and (27)2.92.53.96.4
Table 3. Comparison between the simulation and analytical results for the case of the two bubbles.
Table 3. Comparison between the simulation and analytical results for the case of the two bubbles.
Times [-] T t 1 t 2 t 3
Simulation Results (Figure 4)1.42.63.24.3
Results of Equation (26) and Equation (27)1.52.63.34.5
Table 4. Comparison between the simulation and analytical results for the case of the three bubbles.
Table 4. Comparison between the simulation and analytical results for the case of the three bubbles.
Times [-] T t 1 t 2 t 3
Simulation Results (Figure 4)0.92.42.63.0
Results of Equation (26) and Equation (27)1.02.52.63.1
Table 5. Comparison of the maximum pressure values and maximum pressure arrival times.
Table 5. Comparison of the maximum pressure values and maximum pressure arrival times.
Number of Bubbles1 2 3
Maximum Pressure [-]2.92.63.7
Time at the Maximum Pressure [-]2.92.53.9
Table 6. Times when the pressure waves first contacted the bubble surfaces.
Table 6. Times when the pressure waves first contacted the bubble surfaces.
Number of Bubbles123
Time [-]2.51.00.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ishikawa, T.; Nishida, H.; Tagawa, Y. Numerical Investigation on Influence of Number of Bubbles on Laser-Induced Microjet. Water 2022, 14, 3707. https://doi.org/10.3390/w14223707

AMA Style

Ishikawa T, Nishida H, Tagawa Y. Numerical Investigation on Influence of Number of Bubbles on Laser-Induced Microjet. Water. 2022; 14(22):3707. https://doi.org/10.3390/w14223707

Chicago/Turabian Style

Ishikawa, Tatsumasa, Hiroyuki Nishida, and Yoshiyuki Tagawa. 2022. "Numerical Investigation on Influence of Number of Bubbles on Laser-Induced Microjet" Water 14, no. 22: 3707. https://doi.org/10.3390/w14223707

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