Validation of numerical simulations and experiments on impulse characteristics induced by self-excited oscillation

The high-frequency pulse flow, equivalent to the natural frequency of rocks, is generated by a self-excited oscillating cavity to achieve resonance rock-breaking. The flow field and oscillating mechanism of the self-excited oscillating cavity were simulated using the large eddy simulation method of Computational Fluid Dynamics (CFD). A field-scale testing apparatus was developed to investigate the impulse characteristics and verify the simulation results. The results show that the fluid at the outlet at the tool is deflected due to the pulse oscillation of the fluid. The size and shape of low-pressure vortices constantly change, leading to periodic changes in fluid impedance within the oscillating cavity. The impulse frequency reaches its highest point when the length–diameter ratio is 0.67. As the length–diameter ratio increases, the tool pressure loss also increases. Regarding the cavity thickness, the impulse frequency of the oscillating cavity initially decreases, then increases, and finally decreases again. Moreover, both the impulse frequency and pressure loss increase with an increase in displacement. The numerical simulation findings align with the experimental results, thus confirming the validity of the theoretical model. This research provides theoretical guidance for the practical application of resonance rock-breaking technology.


List of symbols d
Distance to the nearest wall (m) D Diameter of the oscillating cavity (m) L Length of the oscillating cavity (m) i Coordinate direction in the three-dimensional rectangular coordinate system (1, 2, 3) j Direction of the velocity component in the three-dimensional rectangular coordinate system (1, 2, 3) L s Mixing length of the grid (m) P Static pressure (Pa) S ij Rotation of subgrid turbulence tensor (-) t Time taken to flow through the cell (s) T Thickness of the oscillating cavity (m) u i Velocity component in the i direction (m/s) u j Velocity component in the j direction (m/s) V Volume of the cell (m 3 ) ρ Fluid density (kg/m 3 ) μ Dynamic viscosity of the fluid (Pa s) µ t Subgrid turbulent viscous force (N) κ Constant (-) In the oil and gas industry, drilling operations in deep and hard formations are becoming increasingly common.However, enhancing rock-breaking efficiency in these formations poses a significant challenge [1][2][3][4] .The target rocks for deep drilling include volcanic rock, granite, conglomerate, etc., which are characterized by high abrasiveness, hardness, and uniaxial compressive strength exceeding 200 MPa.Resulting in fast bit wear and low rock-breaking efficiency are observed [5][6][7] .Therefore, there is an urgent need to explore new and efficient rock-breaking technologies.One such method is resonance rock-breaking technology, which has been proposed to improve rock-breaking efficiency in deep and hard formations 8,9 .By applying the rock's periodic vibration, resonance can be achieved when the applied vibration frequency matches the natural frequency of the rock 10 .This phenomenon significantly reduces the rock's failure strength and improves rock-breaking efficiency 11 .The self-excited oscillating induced by hydraulic power is capable of generating high-frequency pulsating jets.The self-excited oscillator is widely utilized for surface cleaning and blockage removal in producing oil wells due to its simple structure, compact size, and absence of additional excitation sources [12][13][14] .The self-excited oscillating pulse jet's high pulse pressure and resonance effect can effectively reduce the coal rock breaking threshold pressure and significantly improve the coal rock breaking efficiency 15 .When the pulse frequency of the waterjet is closer to the natural frequency of the rock mass, the resonance effect generated can intensify the degree of rock failure.When the frequency of the pulsed water jet is within a specific range, the rock-breaking impact of the pulsed water jet is significantly better 16 .
Numerous studies have been conducted to examine the properties of waterjet flow and resonance effect, and their effectiveness in breaking rocks.The impact capacity of the waterjet flow and the relaxation of coal can be quantified by evaluating the degree of disturbance and coal relaxation, respectively 17 .Furthermore, a theoretical model has been developed and experimentally validated to determine the self-propelled force of a multi-orifice nozzle 18 .The use of a supercritical CO 2 jet has shown a gradual increase in rock-breaking efficiency with higher jet pressure 19 .Additionally, Dehkhoda and Hood 20 , by controlling the number of impacts, have shown that the pulsation frequency directly influences the failure zone on the impact surface and that internal cracks tend to propagate in the direction of the principal stress.When the pressure was below 50 MPa and the traverse speed exceeded 600 mm/min, the water jet resulted in minimal damage to the coating 21 .Raghu Prasad and Vidya Sagar 22 similarly developed a correlation between fracture energy and acoustic emission energy.Experimental findings have also indicated that increasing the jet pressure results in a decrease in residual stress at the stainless steel surface while increasing the cutting speed elevates the residual stress 23 .Furthermore, both impact depth and rock-breaking volume have been observed to increase with higher water inlet velocity 24 .The incorporation of circular grooves in the inlet and outlet sections has proven effective in reducing fluid resistance 25 .Moreover, the erosion depth has been found to increase proportionally with pump pressure; however, the erosion depth of a convergent-straight nozzle is significantly lower than that of a convergent-divergent nozzle 26 .The center impact pressure of the circular nozzle was found to be highest when the inlet pressure and the target distance were small 27 .The pressure at the nozzle part is low due to the contraction of the thin nozzle, while the velocity is high 28 .In the case of hydro-cutting technology, experimental results showed that the operating time and the distance of the water jet have significant effects on the effectiveness of coal face mining 29 .The frequency changes in an oscillating cavity are primarily caused by the variation of vortices.The development of vortex structures is consistent with the frequency changes 30 .By applying shock vibration and pulse jet, the stress condition of the rock bottom is improved, the cleaning of rock debris at the bottom hole is strengthened, and the rock breaking efficiency is increased 31 .When the impact frequency is close to the resonance frequency of the rock, the displacement response amplitude is the most significant 32 .It is necessary to comprehensively consider the influence of pulse frequency and pressure amplitude on the resonance effect to obtain a better rock-breaking effect 33 .It is concluded that high-frequency and low-amplitude harmonic, dynamic loads can improve the dynamic characteristics of the resonance impact drilling system, thereby improving the drilling efficiency 34 .
Numerous investigations have been carried out to examine the structure of the jet flow field and identify the key factors influencing rock-breaking efficiency, with a primary focus on direct waterjet impingement on the rock 35 .Moreover, several studies have been conducted to achieve impact rock-breaking using mechanical percussion tools, which often possess complex structures and generate lower impacting frequencies 36 .In this paper, combining the advantages of resonance rock breaking and self-excited oscillating, it is proposed to use a self-excited oscillation cavity to generate high-frequency impulse waterjet which can be equivalent to the natural frequency of hard rock (more than 200 Hz) 2,3,37 .Subsequently, the resonance between the impulse jet and the rock will be generated to achieve the resonance rock-breaking of hard rock.This approach significantly differs from rock-breaking by waterjet or mechanical percussion methods.This study proposes a pulse jet that uses self-excited oscillation to break rocks through resonance impact, which can significantly reduce the rockbreaking threshold strength and improve efficiency compared to high-pressure water jet and mechanical impact methods.The tool has minimal activities and sealing parts, resulting in a longer tool life.Previous studies have verified the effectiveness of jet impact resonance 15,16,33 .In this study, the flow field and impulse mechanism of the self-excited oscillating cavity were simulated using the large eddy simulation method of Computational Fluid Dynamics (CFD).An experimental bench for a full-scale self-excited oscillating tool was established to analyze the effects of length-diameter ratio, cavity thickness, and displacement on the cavity's impulse characteristics through experiments and numerical simulations.

Experimental
In the current work, a field-scale experimental apparatus was developed and used to carry out the impulse characteristics experiments induced by self-excited oscillation.

Experimental apparatus
The experimental setup for the self-excited oscillating tool is shown in Fig. 2. The setup consists primarily of three systems: the fluid supply system, the experimental bench system, and the measuring system.The fluid supply system comprises a water tank, a fluid supply line, a high-pressure pump, and a high-pressure line.The water tank, with a volume of 15 m 3 , is the liquid storage container.The fluid supply line, with an inner diameter of 101.6 mm, is used for transporting the low-pressure fluid to the high-pressure pump.The high-pressure pump used is a plunger-type pump with a power of 320 kW, capable of producing a maximum pressure of 35 MPa.The highpressure line, with a diameter of 50.8 mm, is used to transport the high-pressure water flow into the experimental bench system.The experimental bench system includes an experimental bench and a self-excited oscillating tool.The experimental bench, weighing 3.5 tons, is used to secure the self-excited oscillating tool and install the test sensor.The self-excited oscillating tool, with an outer diameter of 178 mm and a length of 0.73 m, is utilized to generate a high-frequency oscillating jet.The self-excited oscillating tool is fixed to the experimental bench by three screw hoops.The measuring system consists of a vibration sensor, a pressure sensor, a flow sensor, and a computer.The vibration sensor is capable of testing the pulse frequency of the self-excited oscillating tool, with a test frequency range of 1-5000 Hz and a sensitivity of 52.35 mV/g.The vibration sensor is a kind of piezoelectric acceleration sensor.The ceramic crystal induces the amount of charge, and the vibration frequency is obtained by subsequent integral processing.The uncertainty of conventional vibration sensors is mainly affected by the temperature and electromagnetic environment.Still, the piezoelectric acceleration sensor used in this experiment is not affected by external electromagnetic and temperature, and the measurement stability and accuracy are high.The pressure sensor, with a measuring range of 0-60 MPa and an accuracy of ± 0.1%, is used to measure the pressure of the self-excited oscillating tool.The uncertainty of the pressure sensor is mainly affected by the temperature and the conveying medium.The temperature of this experiment is the average temperature, and the conveying medium is clean water without impurities, so the measurement stability and accuracy can meet the experimental requirements.The flow sensors are used to measure the fluid flow rate of the self-excited oscillating tool, with a range of 0-150 m 3 /h and an accuracy of 0.5 level.The uncertainty of the flow sensor is mainly affected by the temperature, conveying medium, and vibration.The temperature of this experiment is normal temperature, and the conveying medium is clean water without impurities.The vibration of the flow sensor is small as it is fixed on the experimental bench with the screws.So, the measurement stability and accuracy of the flow sensor can meet the experimental requirements.The computer is used for receiving and processing the sensor signals.The flow direction of the fluid, described by the red arrows, is depicted in Fig. 2a.The tap water stored in the water tank is inhaled by the high-pressure pump through the water supply pipeline.Then, water is pressured by the high-pressure pump and transported to the flow rate and pressure sensor through the high-pressure pipe.Finally, the high-pressure water flow arrives and ejects from the self-excited oscillating tool.

Numerical simulation
The study aims to investigate the impulse characteristics and mechanism generated by the self-excited oscillating cavity and to achieve a higher impulse frequency exceeding 200 Hz.In this paper, the flow field and impulse characteristics of the self-excited oscillating cavity were simulated using a large eddy simulation method of computational fluid dynamics.The analysis focused on evaluating the effects of length-diameter ratio, cavity thickness, and displacement on the impulse characteristics of the self-excited oscillating cavity.

Governing equations
The large eddy simulation (LES) is a significant numerical simulation method in fluid dynamics that has been developed in recent decades.In LES, the flow field is filtered to remove vortices smaller than a certain scale, and only the large-scale vortices are calculated.The solution for this vortex is obtained by solving the additional equation.Typically, the grid-scale is used as the filter scale.The advantage of the large eddy simulation method is that it is more computationally efficient and requires fewer system resources compared to directly solving the Reynolds-averaged Navier-Stokes (RANS) equation or the Direct Numerical Simulation (DNS) equation.Additionally, LES provides greater accuracy compared to turbulence modeling methods 38 .
For this study, the continuity equation and the N-S equation are adopted as the basic governing equations for the large eddy simulation: In the three-dimensional rectangular coordinate system, the coordinate direction is denoted by i (i = 1, 2, 3).The velocity component direction in the three-dimensional rectangular coordinate system is denoted by j (j = 1, 2, 3) u i and u j are the velocity in the direction of x i and x j respectively; t is the time taken to flow through the cell; ρ is fluid density; μ is the dynamic viscosity of the fluid; p is the static pressure.
During the calculation, the N-S equation undergoes a filtering process.The filter function is used to decompose the flow variables into large scale quantities f (x, t) and small scale quantities f ′ (x, t): (1) ∂ρ ∂t where G(x, x ′ ) is the given kernel function, which is called filter function.Equation ( 3) is substituted into N-S equation and continuity equation: The large eddy simulation equation satisfying large scale quantity f ′ (x, t) can be obtained: where τ ij is turbulence subgrid stress.If τ ij is the subgrid tension, the model must be closed: where µ t is subgrid turbulent viscous force; S ij is the rotation of subgrid turbulence tensor; L s is the mixing length of the grid; κ is a constant; d is the distance to the nearest wall; V is the volume of the cell [39][40][41] .

Geometric modeling and meshing
To achieve the desired frequency modulation of the oscillating impulsion, the key structural parameters of the self-excited oscillating cavity underwent optimization.The physical model of impulse characteristics of a selfexcited oscillating cavity is shown in Fig. 3.The continuous flow fluid enters the upper nozzle through the inlet, the jet from the upper nozzle enters the self-excited oscillating cavity, and the jet fluid continues to reach the lower nozzle.Due to the rebound and reflux effect of the wall at the lower nozzle, the continuous jet generated by the upper nozzle becomes a pulsed jet.Finally, the pulse flow is formed from the outlet after it is ejected from the lower nozzle.When the impact frequency of the pulsed jet is equivalent to the natural frequency of the rock, the rock begins to resonate, which can significantly reduce the rock-breaking threshold strength and improve the rock-breaking efficiency.Considering the 8-1/2″ well size, the space limitations in the downhole space, and the fluid up-flow, the outer diameter (D) of the oscillating cavity was set at 120 mm.The length-diameter ratio (L/D) varied within the range of 0.41 to 0.83.The cavity thickness (T) was adjusted between 20 and 70 mm, while the displacement ranged from 15 to 40 L/s.The cone angle of the downstream impact wall was set at 120°.Parameters for the impulse characteristics analysis of the self-excited oscillating cavity are depicted in Table 1.
The SolidWorks 2018 software was used to construct the three-dimensional model.
( The physical mode has meshed using the ANSYS ICEM software.Since large eddy simulation requires highquality mesh, a hexahedral structured mesh was adopted to improve the overall mesh quality.The meshing model is shown in Fig. 4, and the total grid amount is 348,267.

Model assumption and boundary condition
A model assumption was made, focusing solely on the high-frequency impulse characteristics and mechanism generated by the self-excited oscillating cavity, while disregarding the influence of temperature.The fluid flow simulation was carried out using ANSYS Fluent 18.0 software.The SIMPLEC technique was used to calculate the pressure field and velocity field, improving the convergence fitting.The control residual for each equation was maintained below 10 -5 or the number of iterations was less than 5000.The pressure term was discretized using the standardized discrete method, while the convection and diffusion terms were discretized using the secondorder upwind method.A turbulence intensity of 3% and a wall roughness of 10 μm were set 42 .The roughness constant was assigned the default value of 0.5, "velocity inlet" was designated as the inlet boundary, "reflect" was set for the wall boundary, and "pressure outlet" was specified for the outlet boundary 43 .

Results and discussion
By simulating and analyzing the impulse modulation characteristics of the self-excited oscillating cavity with specific parameters (length-diameter ratio of 0.67, cavity thickness of 60 mm, and flow rate of 30 L/s), we can uncover the distribution patterns of the velocity field, pressure field, and turbulent kinetic energy field.This analysis enables us to clarify the impulse modulation mechanism of the self-excited oscillating cavity.

Velocity field
The velocity profile of the self-excited oscillating cavity is shown in Fig. 5.The highest velocity is observed at the upper nozzle exit, reaching a value of 33.1 m/s as depicted in Fig. 5a.The interaction between the free shear jet at the upper nozzle outlet and the pressure disturbance wave generated by the collision feedback of the lower nozzle wall leads to the formation of shear circulation and fluid oscillation, as shown in Fig. 5b.At the outlet of the self-excited oscillating cavity, the fluid undergoes deflection due to the pulsating oscillation of the fluid.

Pressure field
The pressure profile of the self-excited oscillating cavity is shown in Fig. 6.The maximum negative pressure value within the oscillating cavity is observed at the outlet of the upper nozzle, and there exists a region of negative pressure inside the oscillating cavity.As the fluid enters the oscillating cavity, it undergoes acceleration, resulting in a high-speed jet directed into the cavity.Due to the diverging flow of the jet within the cavity, when it reaches the lower section of the oscillating cavity, a portion of the fluid collides with the lower impact wall, subsequently following upward along the oscillating cavity wall.It then collides with the upper impact wall, generating a vortex ring around the central jet beam within the oscillating cavity.This vortex ring, moving at high speed, forms a central area of low pressure in the cavity's center.The presence of these vortex rings exerts a damping effect on the fluid flowing through the middle region between the two vortex rings.

Turbulent kinetic energy
The turbulent kinetic energy profile of the self-excited oscillating cavity is depicted in Fig. 7.The region with the highest turbulent kinetic energy is observed on both sides of the lower nozzle and at the outlet of the tool.This indicates that the fluid in these areas exhibits intense turbulence.Of particular note is the presence of a strong impulse jet at the exit.The flow velocity is higher and there is significant interaction between the fluids on both sides of the lower nozzle, leading to larger turbulent kinetic energy.As a result of the self-excited oscillation chamber, the fluid at the exit of the tool produces an obvious pulse effect.This makes the fluid pulsed fluid, hence the turbulent kinetic energy is large.

Oscillation mechanism
The time evolution of the pressure in the self-excited oscillating cavity is shown in Fig. 8.A pair of symmetrical large-scale vortices is observed around the central jet area within the self-excited oscillating cavity, forming a low-pressure core area.By comparing the cloud images of the pressure field at different time intervals (0.001 to 0.003 s), it becomes apparent that the negative pressure vortex in the self-excited oscillating cavity initially forms near the impact wall of the cavity.It quickly develops into a symmetrical large-scale vortex.The vortex flows upstream along with the impact shear flow, moving periodically in a circumferential manner along the central axis and the cavity wall.The low-pressure vortices undergo a continuous cycle of occurrence, growth, collapse, and regeneration.Consequently, the size and shape of the vortices constantly change.These periodic changes in the fluid impedance within the oscillating cavity have an impact on the central main jet area.This, in turn,

Impulse frequency
The impulse frequency generated by the self-excited oscillating tool is shown in Fig. 9.In Fig. 9a, it is observed that before the high-pressure pump is activated, the velocity amplitude and frequency of the tool are zero.However, as the pump displacement increases, the velocity amplitude and frequency of the tool continuously rise.When the pump displacement reaches its maximum, the tool's velocity amplitude and frequency also reach   their maximum.As the pump displacement starts to decrease, the tool's velocity amplitude and frequency gradually decrease as well.When the pump is turned off, the velocity amplitude and frequency of the tool gradually decrease until reaching zero.By performing a Fourier transform on the original vibration waveform shown in Fig. 9a, the impact frequency graph is obtained which is shown in Fig. 9b. Figure 9b shows two main frequencies: one at 225 Hz and the other around 500 Hz.Since the amplitude of 500 Hz is small, it can be assumed that the frequency generated by the self-excited oscillating tool can reach 225 Hz.This is based on the assumption that the length-diameter ratio is 0.67, the cavity thickness is 60 mm, and the flow rate is 25 L/s.This frequency exceeds the frequency requirement for resonance rock breaking, which is set at (200 Hz).

Validation
Validation results for CFD simulation are shown in Table 2.The frequency generated by the self-excited oscillating tool experiments can reach 225 Hz (assuming a length-diameter ratio is 0.67, a cavity thickness of 60 mm, and a flow rate of 25 L/s).The frequency calculated by the simulation is 240 Hz (assuming a length-diameter ratio is 0.67, a cavity thickness of 60 mm, and a flow rate of 25 L/s).Previous studies have also shown that the self-excited oscillating cavity can generate high frequency pulse oscillating fluid 16,31 .The frequency obtained from the ground test is slightly lower than the result from the simulation, with a deviation of 6.67%, which falls within an acceptable range.The test results have confirmed the accuracy of the numerical model.This deviation  www.nature.com/scientificreports/could be attributed to the fluid pressure loss between the pressure gauge and the self-excited oscillating tool, resulting in a decreased frequency.

Effect of the oscillating cavity parameters
The effects of length-diameter ratio, cavity thickness, and displacement on impulse frequency and pressure loss of the self-excited oscillating cavity were analyzed, leading to the identification of optimal parameters.Since the experimental study on the influence of length-diameter ratio and cavity thickness on the pulse effect of the tool requires constant change of the tool's structural parameters and costly, time-consuming experiments, numerical simulation research was conducted solely for these parameters.Once the optimized parameters for length-diameter ratio and cavity thickness were obtained through numerical simulation, the effect of displacement on the tool's pulse effect was studied through simulation and experiments, with the experimental results serving as a verification for the numerical simulation outcomes.

Length-diameter ratio
Figure 10 illustrates the effects of the length-diameter ratio of the self-excited oscillating cavity on impulse frequency and pressure loss, considering a cavity thickness of 60 mm and a flow rate of 30 L/s.The impulse frequency initially increases and then decreases with the increase in the length-diameter ratio of the oscillating cavity.The highest frequency, reaching 280 Hz, is achieved when the length-diameter ratio is 0.67.The cavity length directly influences the development of shear layer instability and perturbation feedback, which in turn affects the natural frequency of the oscillating cavity.If the cavity length is too short, it becomes challenging to establish a stable shear layer.Consequently, if the cavity length is too long, the feedback disturbance frequency components become excessive, resulting in ineffective excitation and increased energy dissipation in the jet.Additionally, as the length-diameter ratio increases, the pressure loss of the oscillating cavity also increases.This can be attributed to the larger movement distance of the fluid within the oscillating cavity associated with a higher length-diameter ratio, leading to increased pressure loss.

Cavity thickness
Figure 11 presents the effects of cavity thickness on impulse frequency and pressure loss of the self-excited oscillating cavity, considering a length-diameter ratio of 0.67 and a flow rate of 30 L/s.The impulse frequency of the oscillating cavity initially decreases, then remains unchanged, and subsequently decreases again as the cavity thickness increases.The cavity thickness exerts a significant influence on the impulse frequency of the selfexcited oscillating cavity.When the cavity thickness increases (with a cavity diameter is less than 45 mm) the jet velocity within the oscillating cavity steadily decreases, resulting in reduced impact feedback from the jet and a subsequent decrease in the impulse frequency.The cavity thickness exerts a significant influence on the impulse frequency of the self-excited oscillating cavity.When the cavity thickness increases (with a cavity diameter of less than 45 mm), the jet velocity within the oscillating cavity steadily decreases, resulting in reduced impact feedback from the jet and a subsequent decrease in the impulse frequency.Once the cavity thickness exceeds 45 mm, the jet flow within the oscillating cavity stabilizes, leading to a constant impulse frequency.However, when the cavity thickness exceeds 60 mm, the vortex ring's velocity with the oscillating cavity becomes too low, resulting in a decrease in the oscillation frequency and consequently diminishing the effect of the self-excited oscillation impulse jet.The pressure loss of the oscillating cavity decreases as the cavity thickness increases.With an increase in cavity thickness, the flow velocity of the jet within the oscillating cavity gradually decreases, resulting in a continuous reduction in pressure loss.Considering both the output frequency and pressure loss, the optimal cavity thickness is determined to be 60 mm.www.nature.com/scientificreports/Displacement Figure 12 illustrates the effects of displacement on impulse frequency and pressure loss of the self-excited oscillating cavity, considering a cavity thickness of 60 mm and a length-diameter ratio of 0.67.As the displacement increases, the impulse frequency also increases.This is due to the rising velocities of the inlet, interior, and outlet of the self-excited oscillating cavity, along with the shortened period of the negative pressure vortex within the cavity.The feedback time of the lower nozzle wall is shorter, while the feedback frequency of the jet at the upper nozzle is higher, resulting in an elevated pulse frequency.Consequently, as the displacement increases, the pressure loss also increases.This is attributed to the rapidly increasing fluid velocity within the self-excited oscillating chamber, continuous intensification of fluid motion turbulence, and subsequent augmentation of fluid energy consumption, leading to a continuous rise in tool pressure loss.Therefore, if field conditions permit, maximizing the displacement can enhance the impulse frequency.

Conclusions
The flow field and impulse mechanism of the self-excited oscillating cavity were simulated using the large eddy simulation method.The influences of length-diameter ratio, cavity thickness, and displacement on the impulse characteristics of the self-excited oscillating cavity were analyzed, and the theoretical model was validated through experiments.
(1) The upper nozzle outlet exhibits the maximum velocity and negative pressure within the cavity, while the maximum turbulence kinetic energy is observed on both sides of the lower nozzle and at the cavity outlet.At the outlet of the self-excited oscillating cavity, the fluid undergoes deflection due to the pulse oscillation.impedance within the oscillating cavity, which, in turn, affects the central main jet area and changes the pressure and velocity at the outlet of the lower nozzle.(3) Experimental validation of the theoretical model was conducted on the self-excited oscillating tool.The numerical simulation results were found to be in agreement with the experimental findings.The impulse frequency of the oscillating cavity initially increases and then decreases with the increase in the lengthdiameter ratio.The highest frequency is observed when the length-diameter ratio is 0.67.As the length diameter increases, the pressure loss of the oscillating cavity also increases.When it comes to the cavity thickness, the impulse frequency initially decreases, remains unchanged, and then decreases again as the thickness increases.Conversely, the pressure loss decreases with an increase in cavity thickness.A recommended cavity thickness of 60 mm is suggested.Additionally, increasing the displacement leads to an increase in both impulse frequency and pressure loss.As the pump displacement increases, the velocity amplitude and frequency of the oscillating cavity tool continue to rise.(4) The method used in this paper, known as large eddy simulation, is typically employed for numerical simulations of turbulent flow with high Reynolds numbers.The theoretical basis of this method is that there are inertial subscale vortices in high Reynolds number turbulent flow.Therefore, it is not suitable for numerical simulation of low Reynolds number flows.In practical engineering, the complexity of on-site construction conditions necessitates modulation of the structure and construction parameters of the tool to achieve resonance rock breaking and to obtain a better rock-breaking effect.In future studies, it is essential to continue optimizing the tool structure and improving the energy utilization efficiency.

Figure 1 .
Figure 1.Structure and picture of self-excited oscillating tool.

Figure 2 .
Figure 2. (a) Flow chart and (b) image of the experimental apparatus.

3 .
Figure 3. Physical model of impulse characteristics of the self-excited oscillating cavity.

Figure 7 .
Figure 7. Turbulent kinetic energy profile of self-excited oscillating cavity.

Figure 9 .
Figure 9. Impulse frequency generated by the self-excited oscillating tool (a) original frequency diagram and (b) frequency diagram after Fourier change.

Figure 10 .
Figure 10.Effects of length-diameter ratio of the self-excited oscillating cavity on impulse frequency and pressure loss.

Figure 11 .
Figure 11.Effects of cavity thickness of the self-excited oscillating cavity on impulse frequency and pressure loss.

Figure 12 .
Figure 12.Effects of displacement of the self-excited oscillating cavity on impulse frequency and pressure loss.

Table 1 .
Parameters for the impulse characteristics analysis of self-excited oscillating cavity.

Table 2 .
Validation results for CFD simulation.