Optimization Design and Analysis of Bionic Friction Reducing Nozzle in Oil Shale High-Pressure Jet Mining

: The borehole hydraulic mining method has unique advantages for underground oil shale exploitation. Breaking rock with a high-pressure water jet is a crucial step to ensure the smooth implementation of borehole hydraulic mining in oil shale. The hydraulic performance of the nozzle determines the efﬁciency and quality of high-pressure water jet technology. To obtain a superior hydraulic performance nozzle, based on the bionic non-smooth theory, a circular groove was selected as the bionic unit to design a bionic straight cone nozzle. The structural parameters of the circular groove include the groove depth, width, and slot pitch. The optimization objective was to minimize the pressure drop, where the ﬂuid has the least resistance. A genetic algorithm was used to optimize the structural parameters of the circular grooves in the inlet and outlet sections of the bionic straight cone nozzle. The optimal structural parameters of the nozzle were as follows: the inlet diameter was 15 mm, the inlet length was 20 mm, the outlet diameter was 4 mm, the length-to-diameter ratio was 3, and the contraction angle was 30 ◦ . In addition, in the inlet section, the groove width, slot pitch, and groove depth were 3.9 mm, 5.2 mm, and 5.5 mm, respectively, and the number of circular grooves was 2. Moreover, in the outlet section, the groove width, slot pitch, and groove depth were 2.25 mm, 3 mm, and 5.5 mm, respectively, and the number of circular grooves was 2. The CFD numerical simulation results showed that under the same numerical simulation conditions, compared with the conventional straight cone nozzle, the bionic straight cone nozzle velocity increase rate could reach 13.45%. The research results can provide scientiﬁc and valuable references for borehole hydraulic mining of high-pressure water jets in oil shale drilling.


Introduction
Oil shale is an unconventional oil and gas resource, as coal, oil, and natural gas provide non-renewable fossil energy [1,2]. Such resources worldwide are abundant, equivalent to more than 400 billion tons of synthetic shale oil, three times the world's crude reserves [3,4]. The borehole hydraulic mining method has unique advantages for underground oil shale exploitation [5]. As shown in Figure 1, when borehole hydraulic mining technology is used to mine underground oil shale, first, drilling rigs must be used to drill holes that are 100~1000 mm in diameter and pass through oil shale layers. Second, the high-pressure water gun is inserted into the borehole, and the underground oil shale bed is broken into several small pieces by the high-pressure water jet from the nozzle connected to the end of the ejector. Then, it is stripped from the parent rock to form a solid-liquid two-phase flowable slurry at the bottom of the borehole. Finally, pneumatic or hydraulic lifting devices are used to extract slurry to the surface [6]. Borehole hydraulic mining technology has the advantages of low cost, short cycle, reliable use, easy operation, low labor intensity, and little environmental pollution [7]. High-pressure water jet broken rock is crucial to ensure the smooth implementation of oil shale borehole hydraulic mining [8,9]. The entire underground oil shale deposit has to be broken into small pieces of rock and stripped off. A series of follow-up work can be carried out [10,11]. The hydraulic performance of the nozzle determines the efficiency and quality of high-pressure water jet technology [12,13]. Thus, it is necessary to study the hydraulic performance of the nozzle. Due to the mechanical properties of oil shale, it is most easily broken when the water jet is parallel to the bedding direction of the oil shale [14,15]. The primary flow resistance of the straight cone nozzle appears in the outlet section. When drilling hydraulic mining oil shale, a high-pressure water jet often uses a straight cone nozzle because the straight cone nozzle can improve the jet quality [16,17]. The structural parameters of the straight-cone nozzle are the key factors affecting the spray performance of the nozzle [18,19]. The transition from the convergence section to the outlet section must be optimized. It is necessary to reduce the outlet section length as much as possible while satisfying the beam-gathering ability. All the optimization work involves increasing rock breaking effect, saving hydraulic energy, and improving energy conversion efficiency [20,21]. Wen et al. [22] chose a circular annular groove as a bionic unit and arranged the circular annular groove evenly on the inner cavity surface of a straight cone nozzle. The original smooth flow channel inside the nozzle became a bionic non-smooth surface structure, effectively improving the flow channelʹs hydraulic characteristics inside the nozzle. However, how the size of the circular annular groove affects the hydraulic characteristics of the internal flow channel of the nozzle has not been studied.
In summary, based on the existing pressure drop formula of the variable cross-section of a circular tube, this study deduced the theoretical formula for the pressure drop calculation of a bionic circular tube (BCT). Using a genetic algorithm, combined with the pressure drop calculation formula of a BCT, the minimum pressure drop was taken as the optimization goal. The circular annular groove structure parameters of the inlet and outlet sections of the bionic straight cone nozzle (BSCN) were optimized. The optimal structural parameters of the BSCN were obtained by using genetic algorithm optimization analysis and CFD numerical simulation results. The drag reduction characteristics of BSCN were further analyzed. The research results can provide scientific and valuable references for borehole hydraulic mining of high-pressure water jets in oil shale drilling. At the same High-pressure water jet broken rock is crucial to ensure the smooth implementation of oil shale borehole hydraulic mining [8,9]. The entire underground oil shale deposit has to be broken into small pieces of rock and stripped off. A series of follow-up work can be carried out [10,11]. The hydraulic performance of the nozzle determines the efficiency and quality of high-pressure water jet technology [12,13]. Thus, it is necessary to study the hydraulic performance of the nozzle. Due to the mechanical properties of oil shale, it is most easily broken when the water jet is parallel to the bedding direction of the oil shale [14,15]. The primary flow resistance of the straight cone nozzle appears in the outlet section. When drilling hydraulic mining oil shale, a high-pressure water jet often uses a straight cone nozzle because the straight cone nozzle can improve the jet quality [16,17]. The structural parameters of the straight-cone nozzle are the key factors affecting the spray performance of the nozzle [18,19]. The transition from the convergence section to the outlet section must be optimized. It is necessary to reduce the outlet section length as much as possible while satisfying the beam-gathering ability. All the optimization work involves increasing rock breaking effect, saving hydraulic energy, and improving energy conversion efficiency [20,21]. Wen et al. [22] chose a circular annular groove as a bionic unit and arranged the circular annular groove evenly on the inner cavity surface of a straight cone nozzle. The original smooth flow channel inside the nozzle became a bionic non-smooth surface structure, effectively improving the flow channel's hydraulic characteristics inside the nozzle. However, how the size of the circular annular groove affects the hydraulic characteristics of the internal flow channel of the nozzle has not been studied.
In summary, based on the existing pressure drop formula of the variable cross-section of a circular tube, this study deduced the theoretical formula for the pressure drop calculation of a bionic circular tube (BCT). Using a genetic algorithm, combined with the pressure drop calculation formula of a BCT, the minimum pressure drop was taken as the optimization goal. The circular annular groove structure parameters of the inlet and outlet sections of the bionic straight cone nozzle (BSCN) were optimized. The optimal structural parameters of the BSCN were obtained by using genetic algorithm optimization analysis and CFD numerical simulation results. The drag reduction characteristics of BSCN were further analyzed. The research results can provide scientific and valuable references for borehole hydraulic mining of high-pressure water jets in oil shale drilling. At the same time, it provides a scientific and valuable reference for improving the hydraulic performance of conventional nozzles. The bionic non-smooth surface drag reduction technology is based on the principle of bionics [23]. Changing the morphological characteristics of the object surface will change the kinematic and dynamic characteristics of the fluid flowing along the wall [24]. Research on bionic non-smooth surface-drag-reduction technology mainly focused on different rib surfaces [25,26]. Currently, two main drag-reduction mechanisms of bionic non-smooth surfaces are recognized by many scholars, namely, the 'second vortex group' theory and the 'bulgy height' theory [27,28]. The 'second vortex group' theory holds that a small secondary vortex is generated due to the interaction between the backward-rotating flow vortex and the tip of the groove. The existence of the groove restricts the movement of the flow direction vortex, which makes the wall suddenly weak, and ultimately reduces the wall friction [29,30]. The 'bulgy height' theory argues that the existence of grooves increases the thickness of the viscous bottom in the boundary layer, resulting in a decrease in friction, which was confirmed by many researchers [31].

Modeling and Optimization
The rib drag reduction method was used in many fields and was shown to positively reduce fluid resistance [32,33]. In this study, according to the typical L-shaped rib, we proposed a new type of bionic circular tube (BCT), where the schematic diagram is shown in Figure 2. time, it provides a scientific and valuable reference for improving the hydraulic performance of conventional nozzles.

BCT Model and Pressure Drop Calculation
The bionic non-smooth surface drag reduction technology is based on the principle of bionics [23]. Changing the morphological characteristics of the object surface will change the kinematic and dynamic characteristics of the fluid flowing along the wall [24]. Research on bionic non-smooth surface-drag-reduction technology mainly focused on different rib surfaces [25,26]. Currently, two main drag-reduction mechanisms of bionic nonsmooth surfaces are recognized by many scholars, namely, the 'second vortex group' theory and the 'bulgy height' theory [27,28]. The 'second vortex group' theory holds that a small secondary vortex is generated due to the interaction between the backward-rotating flow vortex and the tip of the groove. The existence of the groove restricts the movement of the flow direction vortex, which makes the wall suddenly weak, and ultimately reduces the wall friction [29,30]. The 'bulgy height' theory argues that the existence of grooves increases the thickness of the viscous bottom in the boundary layer, resulting in a decrease in friction, which was confirmed by many researchers [31].
The rib drag reduction method was used in many fields and was shown to positively reduce fluid resistance [32,33]. In this study, according to the typical L-shaped rib, we proposed a new type of bionic circular tube (BCT), where the schematic diagram is shown in Figure 2. As shown in Figure 2, the typical non-smooth surface bionic prototype mainly includes shield scales of shark skin, the mastoid structure of the lotus leaf, the spherical nodular structure of a humpback fin, and radial stripes of scallops. A BCT model was proposed based on the bionic non-smooth surface drag reduction theory. A circular annular groove was uniformly added to the internal flow channel surface of the conventional circular pipe according to certain rules as the bionic unit so that the internal flow channel of the conventional circular pipe was changed from the original smooth flow As shown in Figure 2, the typical non-smooth surface bionic prototype mainly includes shield scales of shark skin, the mastoid structure of the lotus leaf, the spherical nodular structure of a humpback fin, and radial stripes of scallops. A BCT model was proposed based on the bionic non-smooth surface drag reduction theory. A circular annular groove was uniformly added to the internal flow channel surface of the conventional circular pipe according to certain rules as the bionic unit so that the internal flow channel of the conventional circular pipe was changed from the original smooth flow channel to the bionic non-smooth surface structure. This could improve the kinematic and dynamic characteristics of the fluid flowing on the pipe wall to a certain extent.

Modeling for a Single-Type Variable Cross-Sectional Structure
A single type of circular tube with a variable cross-sectional structure has a reduced or increased diameter, as shown in Figure 3. In a Newtonian fluid, the structure has a mature formula for calculating the pressure drop [34]. Appl channel to the bionic non-smooth surface structure. This could improve the kinematic and dynamic characteristics of the fluid flowing on the pipe wall to a certain extent.

Modeling for a Single-Type Variable Cross-Sectional Structure
A single type of circular tube with a variable cross-sectional structure has a reduced or increased diameter, as shown in Figure 3. In a Newtonian fluid, the structure has a mature formula for calculating the pressure drop [34].  [35][36][37] (as shown in Figure 3a), the pressure drop calculation formula is: For the increased pipe diameter structure [35][36][37] (as shown in Figure 3b), the pressure drop calculation formula is: In these formulas: s p  and k p  -reduced or increased pipe diameter structure pres- When the fluid density is constant, the velocity is inversely proportional to the crosssectional area. The fluid flow through the contraction channel will accelerate and the flow through the expansion channel will decelerate. Combined with the continuous equation of fluid mechanics, the pressure drops of the pipe diameter reduction or increased structure were further simplified and expressed as: For the increased pipe diameter structure [35][36][37] (as shown in Figure 3b), the pressure drop calculation formula is: In these formulas: ∆p s and ∆p k -reduced or increased pipe diameter structure pressure drop (MPa), respectively; D and d-large diameter and small diameter (m), respectively; v D and v d -mean velocity of the large diameter and small diameter cross-sections (m/s), respectively; ρ-liquid density (kg/m 3 ).
When the fluid density is constant, the velocity is inversely proportional to the crosssectional area. The fluid flow through the contraction channel will accelerate and the flow through the expansion channel will decelerate. Combined with the continuous equation of fluid mechanics, the pressure drops of the pipe diameter reduction or increased structure were further simplified and expressed as: In these formulas: ∆p k1 and ∆p s1 -simplified increased and reduced pipe diameter structure pressure drops (MPa), respectively; d-trail diameter (m); h-groove depth (m); v D and v d -mean velocities of large diameter and small diameter cross-sections (m/s), respectively; A D and A d -the cross-sectional areas of the large and small paths (m 2 ), respectively; ρ-liquid density (kg/m 3 ).

Modeling for a Comprehensive-Type Variable Cross-Sectional Structure
The comprehensive type of circular tube with a variable cross-sectional structure is a reduction-increase or increase-reduction structure, as shown in Figure 4. For a Newtonian fluid, the structure has a related mature formula for calculating the pressure drop [38,39].
In these formulas:

Modeling for a Comprehensive-Type Variable Cross-Sectional Structure
The comprehensive type of circular tube with a variable cross-sectional structure is a reduction-increase or increase-reduction structure, as shown in Figure 4. For a Newtonian fluid, the structure has a related mature formula for calculating the pressure drop [38,39]. In Figure 4, and are the length of the reduction section and increase section (m), respectively. For the pipe diameter reduction-increase structure [35,40,41] (as shown in Figure 4a) In these formulas:  In Figure 4, L d and L D are the length of the reduction section and increase section (m), respectively. For the pipe diameter reduction-increase structure [35,40,41] (as shown in Figure 4a), the pressure drop calculation formula is In these formulas: ∆p sk is the reduction-increase pipe diameter structure pressure drop; K d is the influence coefficient of the pressure drop of the increased structure; and For the increase-reduction pipe diameter structure [35,40,41] (as shown in Figure 4b), the pressure drop calculation formula is In these formulas: ∆p ks is the pressure drop of the increase-reduction pipe diameter structure; K D is the influence coefficient of the pressure drop of the increased structure; and n D = L D /D, where 0.15 ≤ n D ≤ 0.37.
Combined with the geometric characteristics of the BCT, the pressure drop of the reduction-increase and increase-reduction pipe diameter structures was further simplified and expressed as In these formulas: ∆p sk1 -reduction-increase BCT pipe diameter structure pressure drop (MPa); ∆p ks1 -increase-reduction BCT pipe diameter structure pressure drop (MPa); s-slot pitch (m); d-trail diameter (m); h-groove depth (m); w-groove width (m); ρ-liquid density (kg/m 3 ).

Pressure Drop Calculation
According to the superposition of a pressure drop, the pressure drop calculation formula of a BCT is extended as follows: In this formula: ∆p-BCT pressure drop (MPa); ∆p sk1 -reduction-increase BCT pipe diameter structure pressure drop (MPa); ∆p ks1 -increase-reduction BCT pipe diameter structure pressure drop (MPa). As shown in Figure 5, the inlet section L1 and outlet section L3 of the conventional straight cone nozzle were circular tubes. The BSCN model was proposed based on the bionic non-smooth surface drag reduction theory, which can effectively improve the kinematic and dynamic characteristics of fluid flow along the pipe wall to a certain extent.

Optimization Process
The genetic algorithm is a bionic algorithm for solving complex global optimization based on Darwin's natural evolution theory and genetic variation theory [42]. It draws on natural selection and natural genetic mechanisms in the biological world. It conducts a random search in the solution space based on probability theory and finally finds the op- As shown in Figure 5, the inlet section L1 and outlet section L3 of the conventional straight cone nozzle were circular tubes. The BSCN model was proposed based on the bionic non-smooth surface drag reduction theory, which can effectively improve the kinematic and dynamic characteristics of fluid flow along the pipe wall to a certain extent.

Optimization Process
The genetic algorithm is a bionic algorithm for solving complex global optimization based on Darwin's natural evolution theory and genetic variation theory [42]. It draws on natural selection and natural genetic mechanisms in the biological world. It conducts a random search in the solution space based on probability theory and finally finds the optimal solution to the problem [43,44]. A genetic algorithm is a mature optimization algorithm with unique advantages for solving multi-variable optimization problems. It can save the time and cost of experimentation and provide more accurate results. Figure 6 shows the optimization process of the structural parameters of the BSCN. The BSCN optimization variables' parameter combinations were converted to the genetic algorithm to identify the search space for processing, and chromosome coding was realized for the parameter data. The constraint function created each generated initial population within the range of values. In the optimization process of genetic algorithms, fitness is an important indicator used to evaluate individual performance. This study selected the objective function as the fitness function to evaluate the fitness. As shown in Figure 5, the inlet section L1 and outlet section L3 of the conventional straight cone nozzle were circular tubes. The BSCN model was proposed based on the bionic non-smooth surface drag reduction theory, which can effectively improve the kinematic and dynamic characteristics of fluid flow along the pipe wall to a certain extent.

Optimization Process
The genetic algorithm is a bionic algorithm for solving complex global optimization based on Darwin's natural evolution theory and genetic variation theory [42]. It draws on natural selection and natural genetic mechanisms in the biological world. It conducts a random search in the solution space based on probability theory and finally finds the optimal solution to the problem [43,44]. A genetic algorithm is a mature optimization algorithm with unique advantages for solving multi-variable optimization problems. It can save the time and cost of experimentation and provide more accurate results. Figure 6 shows the optimization process of the structural parameters of the BSCN. The BSCN optimization variablesʹ parameter combinations were converted to the genetic algorithm to identify the search space for processing, and chromosome coding was realized for the parameter data. The constraint function created each generated initial population within the range of values. In the optimization process of genetic algorithms, fitness is an important indicator used to evaluate individual performance. This study selected the objective function as the fitness function to evaluate the fitness.

Objective Function and Constraint Function
To further study the structural parameters of the groove of the BSCN, the pressure drop formula of the BCT was taken as the objective function for the inlet section L1 and the outlet section L3. As shown in Figure 7, a pressure drop calculation formula of the BCT was obtained by selecting a circular groove element:

Objective Function and Constraint Function
To further study the structural parameters of the groove of the BSCN, the pressure drop formula of the BCT was taken as the objective function for the inlet section L1 and the outlet section L3. As shown in Figure 7, a pressure drop calculation formula of the BCT was obtained by selecting a circular groove element:   The example Formula (18) was further simplified to the objective function: (19) In these formulas: d-trail diameter (m); w-groove width (m); s-slot pitch (m); h-groove depth (m); ρ-liquid density (kg/m 3 ); v d -flow rate at the inlet (m/s).

Basic Structure Parameters
Combined with the structural parameters of the conventional straight cone nozzle, the structural parameters of the BSCN were determined. As shown in Table 1, the BSCM mainly increased the three structural parameters of groove depth, groove width, and slot pitch compared with the conventional straight cone nozzle.  (19) of the BCT, the inlet velocity of the inlet section L1 of the bionic straight cone nozzle was 1 m/s, the inlet velocity of the outlet section L3 of the BSCN was 14 m/s, and the liquid density was 1 kg/m 3 . The inlet section L1 and outlet section L3 of the BSCN were regarded as two BCTs, and the minimum pressure drop of the BCT was taken as the optimization objective. In the range of 0.1 mm-14.5 mm in groove depth, the genetic algorithm was used to analyze and solve the pressure drop example in Formula (19), and the groove width and slot pitch of the structural parameters of the BSCN were obtained.

Optimization Results
The inlet section L1 and outlet section L3 of the BSCN were regarded as two BCTs. The inlet velocity of the inlet section L1 of the BSCN was 1 m/s, the inlet velocity of the outlet section L3 of the BSCN was 14 m/s, and the liquid density was 1 kg/m 3 . In the groove depth range 0.1 mm-14.5 mm, the genetic algorithm was used to analyze and solve Formula (19) for the pressure drop example and the groove width, and the slot pitch of the structural parameters of a circular groove unit of the BSCN was obtained.
As shown in Figure 8, in the groove depth range of 0.1 mm~14.5 mm, the minimum pressure drop of the BCT was taken as the optimization objective, and the genetic algorithm was used to solve the objective function. The values of the groove width and slot pitch could be obtained. The values of the groove width and slot pitch spacing obtained showed a linear change trend. mula (19) for the pressure drop example and the groove width, and the slot pitch of the structural parameters of a circular groove unit of the BSCN was obtained.
As shown in Figure 8, in the groove depth range of 0.1 mm ~ 14.5 mm, the minimum pressure drop of the BCT was taken as the optimization objective, and the genetic algorithm was used to solve the objective function. The values of the groove width and slot pitch could be obtained. The values of the groove width and slot pitch spacing obtained showed a linear change trend.

Determining the Number of BSCN Grooves
Based on the optimization results of the genetic algorithm, the number of grooves in the inlet and outlet sections of the BSCN was determined by combining the terms in the following Formula (21), as shown in Figure 9. The models of the BSCN and conventional straight cone nozzle were established.

Determining the Number of BSCN Grooves
Based on the optimization results of the genetic algorithm, the number of grooves in the inlet and outlet sections of the BSCN was determined by combining the terms in the following Formula (21), as shown in Figure 9. The models of the BSCN and conventional straight cone nozzle were established.
In this formula, n1 and n3 are the maximum positive integers, respectively; w1 and s1 are the groove width and slot pitch of the inlet section of the BSCN, respectively; and L1 is the length of the inlet section of the BSCN. w3 and s3 are the groove width and slot pitch of the outlet section of the BSCN, respectively, and L3 is the length of the outlet section of the BSCN.

Basic Settings for Numerical Simulation
COMSOL Multiphysics software was used for the simulation. The inlet velocity was 1 m/s and the outlet was under a static pressure condition. The fluid was incompressible, the wall material of the model was set to ordinary steel, the turbulence model (k − ξ) was selected, the mesh size was selected to be more refined, and the calibration was fluid dynamics. The maximum unit size was 0.91 mm, the minimum unit size was 0.026 mm, and the boundary layer number was set to 8. The solver selected was the steady-state solver. The velocity in each model was analyzed using CFD numerical simulation.

Analysis of Simulation Results
As shown in Figure 10, based on CFD numerical simulation results, the velocity along the L3 axis of each BSCN and the conventional straight cone nozzle was extracted, and then the average value was taken. The average velocity along the central L3 axis of each BSCN was subtracted from the average velocity on the L3 central axis of the conventional straight cone nozzle. The obtained results were divided by the average velocity along the central axis of L3 in the outlet of the conventional straight cone nozzle, and finally, the velocity increase rate of the BSCN was obtained. Then, the curves of the velocity rise rate of each BSCN were drawn, as shown in Figure 11.  The average velocity along the central L3 axis of each BSCN was subtracted from the average velocity on the L3 central axis of the conventional straight cone nozzle. The obtained results were divided by the average velocity along the central axis of L3 in the outlet of the conventional straight cone nozzle, and finally, the velocity increase rate of the BSCN was obtained. Then, the curves of the velocity rise rate of each BSCN were drawn, as shown in Figure 11.
In this formula: η 0 -the rate of velocity increase; V bn -average velocity along the axis of the BSCN (m/s); V cn -average velocity along the central axis of the conventional straight cone nozzle (m/s).
It can be observed from Figure 11 that the velocity increase rate of the BSCN was more significant than zero after the structural parameters were optimized. It was further noted that the resistance against the fluid flowing through the inner wall was reduced by setting circular grooves in the inlet and outlet sections of the conventional straight cone nozzle, and the drag reduction effect of the outlet section of the BSCN was better than that of the inlet section.
In the optimized BSCN model, in the L1 inlet section, the groove width was 3.9 mm, the slot pitch was 5.2 mm, the groove depth was 5.5 mm, and the number of grooves was 2.
In the L3 outlet section, the groove width was 2.25 mm, the slot pitch was 3 mm, the groove depth was 5.5 mm, and the number of grooves was 2. The maximum velocity increase rate was 13.45%, and the drag reduction effect was the best. BSCN was obtained. Then, the curves of the velocity rise rate of each BSCN were drawn, as shown in Figure 11.    Figure 12 is the velocity vector diagram of the BSCN with the best drag reduction effect and a velocity increase rate of 13.45%. In Figure 12, each circular annular groove contained a swirl opposite the direction of the internal flow. Appl It can be observed from Figure 11 that the velocity increase rate of the BSCN was more significant than zero after the structural parameters were optimized. It was further noted that the resistance against the fluid flowing through the inner wall was reduced by setting circular grooves in the inlet and outlet sections of the conventional straight cone nozzle, and the drag reduction effect of the outlet section of the BSCN was better than that of the inlet section.

Analysis of the Effect
In the optimized BSCN model, in the L1 inlet section, the groove width was 3.9 mm, the slot pitch was 5.2 mm, the groove depth was 5.5 mm, and the number of grooves was 2. In the L3 outlet section, the groove width was 2.25 mm, the slot pitch was 3 mm, the groove depth was 5.5 mm, and the number of grooves was 2. The maximum velocity increase rate was 13.45%, and the drag reduction effect was the best. Figure 12 is the velocity vector diagram of the BSCN with the best drag reduction effect and a velocity increase rate of 13.45%. In Figure 12, each circular annular groove contained a swirl opposite the direction of the internal flow. As shown in Figure 13, the BSCN and the conventional straight cone nozzle were divided perpendicular to the axis. The first dividing line was at the left end of the first circular groove, the second dividing line was located at the middle line of the first circular groove slot pitch in the inlet section, the third dividing line was located at the right end of L1, the fourth dividing line was located at the left end of L3, the fifth dividing line was located at the middle line of the first circular groove slot pitch in the outlet section, and the sixth dividing line was located at the right end of L3. As shown in Figure 13, the BSCN and the conventional straight cone nozzle were divided perpendicular to the axis. The first dividing line was at the left end of the first circular groove, the second dividing line was located at the middle line of the first circular groove slot pitch in the inlet section, the third dividing line was located at the right end of L1, the fourth dividing line was located at the left end of L3, the fifth dividing line was located at the middle line of the first circular groove slot pitch in the outlet section, and the sixth dividing line was located at the right end of L3. Appl It can be observed from Figure 11 that the velocity increase rate of the BSCN was more significant than zero after the structural parameters were optimized. It was further noted that the resistance against the fluid flowing through the inner wall was reduced by setting circular grooves in the inlet and outlet sections of the conventional straight cone nozzle, and the drag reduction effect of the outlet section of the BSCN was better than that of the inlet section.

Analysis of the Effect
In the optimized BSCN model, in the L1 inlet section, the groove width was 3.9 mm, the slot pitch was 5.2 mm, the groove depth was 5.5 mm, and the number of grooves was 2. In the L3 outlet section, the groove width was 2.25 mm, the slot pitch was 3 mm, the groove depth was 5.5 mm, and the number of grooves was 2. The maximum velocity increase rate was 13.45%, and the drag reduction effect was the best. Figure 12 is the velocity vector diagram of the BSCN with the best drag reduction effect and a velocity increase rate of 13.45%. In Figure 12, each circular annular groove contained a swirl opposite the direction of the internal flow. As shown in Figure 13, the BSCN and the conventional straight cone nozzle were divided perpendicular to the axis. The first dividing line was at the left end of the first circular groove, the second dividing line was located at the middle line of the first circular groove slot pitch in the inlet section, the third dividing line was located at the right end of L1, the fourth dividing line was located at the left end of L3, the fifth dividing line was located at the middle line of the first circular groove slot pitch in the outlet section, and the sixth dividing line was located at the right end of L3. As shown in Figure 14, based on the results of the CFD numerical simulation, the value of the velocity at each dividing line was extracted, and the curve of the change in the velocity at each dividing line with the distance of the dividing line was plotted. Line01/line1 to line06/line6 represent the positions shown in Figure 13. As shown in Figure 14, based on the results of the CFD numerical simulation, the value of the velocity at each dividing line was extracted, and the curve of the change in the velocity at each dividing line with the distance of the dividing line was plotted. Line01/line1 to line06/line6 represent the positions shown in Figure 13. It can be concluded from Figure 14a that in the inlet section, the velocity of the BSCN at different positions perpendicular to the axis was greater than that of the conventional straight cone nozzle at a distance of 2 mm from the wall. It can be seen from Figure 14b that in the outlet section, at the position of 0.75 mm away from the wall, the velocity of the BSCN at different positions perpendicular to the axis was greater than that of the conventional straight cone nozzle. In the outlet section, compared with the conventional straight cone nozzle, the BSCN improved the fluid velocity and had a more vital constraint on the fluid.

Conclusions
(1) Under the same numerical simulation conditions, compared with the conventional straight cone nozzle, the velocity increase rate of the BSCN could reach 13.45%. This can provide valuable scientific references for the hydraulic mining of high-pressure water jets in oil shale drilling. (2) CFD numerical simulation results showed that the circular groove in the inlet and outlet sections could effectively reduce the fluid resistance. The drag reduction effect of the circular groove arranged in the outlet section of BSCN was better than that arranged in the inlet section. (3) Based on the results of genetic algorithm optimization analysis and CFD numerical simulation, the optimal structural parameters of the BSCN were as follows: Fixed values were inlet diameter D = 15 mm, inlet length L1 = 20 mm, outlet diameter d = 4 mm, L3/d = 3, and contraction angle θ = 30 • . Variable values for the inlet section of L1, namely, the groove width, slot pitch, and groove depth, were 3.9 mm, 5.2 mm, and 5.5 mm, respectively, and the number of circular grooves was 2. In the outlet section of L3, the groove width, slot pitch, and groove depth were 2.25 mm, 3 mm, and 5.5 mm, respectively, and the number of circular grooves was 2.
Future work is planned to manufacture the real object of the BSCN and conduct a series of experiments. At the same time, in addition to the bionic unit used in this study, we can also try to find different types of bionic units to explore nozzles with better hydraulic performance.