Simulation Study on Dynamics of Hydraulic Turbines Used in Drilling Engineering

Hydroturbines have a very wide range of applications, which are commonly found in wind turbines, water turbines, aero engines, etc. ,is paper provided a detailed turbine design and a design method of turbine blade shape. Using the CFD (computational fluid dynamics) method, based on the realizable k-ε turbulence model and Euler multiphase flow model, the effects of different external loads, blade numbers, blade installation angles, and flow rates on the force condition of turbine and the influence of different solid contents, particle sizes, and densities on turbine performance were studied.,e simulation results show that, under the action of fluid, when the starting torque of turbine is larger than the external load, the turbine starts to move, the angular velocity increases until it remains constant, the absolute value of impact force decreases, and the impact torque decreases until it is equal to the external load; while the starting torque of turbine is smaller than the external load, the turbine stays still. ,e increase of the particle size, content, and density of the solid phase will lead to an increase in the torque and pressure drop of the turbine and ultimately leads to the increase of turbine input, output power, and efficiency.


Introduction
Hydroturbines have a very wide range of applications, which are commonly found in wind turbines, water turbines, aero engines, etc. [1,2]. Blade components also have important applications in drilling engineering and are often used as turbine power generators and turbine drilling tools as downhole power devices [3]. Underground turbine generator is one of the common underground power supplies. e downhole turbine generator acts as an axial flow impeller, and the turbine is rotated by the impact of the drilling fluid circulating in the downhole. At the same time, the magnetic coupling of the magnetic coupler drives the permanent magnet generator to generate electricity, thereby achieving the function of continuous power supply for downhole instruments. Turbine drilling tools [1,2,4] also belong to the typical blade-type axial flow machinery, which works by means of the interaction between the blades and the fluid medium and realizes the conversion of hydraulic energy and mechanical energy based on the principle of momentum moment and finally the speed and torque of the turbine shaft output.
With the increasing difficulty of drilling and production in oil and gas fields, downhole drilling tools have been developed rapidly, their power gradually increases, the working time becomes longer and longer, and harsh underground wells, horizontal wells, multibranch wells, LWD (logging while drilling) measuring instruments, and downhole intelligent tools such as mud loggers and rotary guides are indispensable [5,6]. ese instruments and tools require electrical energy to work. With the increase of measurement parameters and the complexity of downhole control mechanisms, the requirements for power supply are becoming higher and higher. erefore, it is very important to reliably provide continuous and stable electrical energy for underground instruments and tools.
In the purely mechanical structure, fluid pressure [7] is often used as a source of power, compared to the motor drive system; hydraulic system has higher power density [8]. e fluid power technology can achieve high power output with low system mass compared to other common powertrain technologies [9]. Hydraulic turbine [10][11][12] is the most common energy conversion device that converts hydraulic pressure into rotary motion, which is commonly used [2,[13][14][15][16][17][18][19][20][21][22][23][24] in wind energy, hydraulic turbines, hydraulic motors, etc. Ge et al. [15] studied the Reynolds number of airfoil aerodynamic performances and its influence on the optimal design of wind turbine rotor and discussed the maximum power coefficient; Hansen [16] developed and tested an airfoil optimization method for wind turbine applications, which can control the performance loss caused by leading edge contamination. In view of the failure of traditional fluid machinery in supersonic flow, Paniagua et al. [17] designed and analyzed advanced high-sonic axial flow turbines. Pereiras et al. [18] calculated the efficiency of the twin-turbine structure based on the results of the numerical model and paid attention to the influence of the turbine's torque and flow rate under the reverse mode conditions. Kwon et al. [19] used the Rayleigh-Ritz hypothesis model method to analyze a turbine generator, determined the resonant frequency and intensity of the rotating multipackage blade system excited by multiple nozzle forces, and studied the influence of different parameters such as the number of nozzles and damping coefficient on the system response. Müller et al. [20] analyzed the wind energy converter and discussed the modern transformation of this traction type energy converter for building integration. Lanzafame and Messina [21] improved the blade element momentum theory and established a mathematical model for the design of hydrodynamic wind turbines at low wind speed, which is used to optimize the performance of the rotor. e axial flow Wells turbine [22] is widely operated in most energy stations. e principle of operation of the Wells turbine is based on aerodynamics, and several variants of the Wells turbine have been proposed and developed [23,24]. e performance and shape of Wells turbines have been optimized by numerical [25] and experimental [26] methods for each of the sea states [27]. A novel twin-rotor radialinflow air turbine was recently proposed [28], which has two sets of rotor blades equipped on a common shaft that are axially offset from each other. Another new topology that uses twinunidirectional turbines has been proposed and studied [29].
is paper provided a detailed turbine design and design method of turbine blade shape. At present, the research on downhole turbines is static. In this paper, by using FLUENT software and sliding mesh method, UDF (user-defined function) file is compiled to simulate the dynamic performance of downhole turbine. e effects of external load, number of blades, flow rate, and blade installation angle on the hydrodynamic performance of downhole turbine are studied, and the start-up characteristics and movement process of downhole turbine are analyzed.

Modeling and Methods
In the drilling industry, downhole hydraulic turbine is often used as the power mechanism of downhole drilling tools. For example, as a part of the turbine generator, it provides power for the underground electronic mechanism. It can also be used as the rotor of turbodrill to drive other mechanical structures directly by using the rotating characteristics of the turbine. e downhole turbine can be divided into cascade turbine and spiral turbine according to its form. In this paper, the traditional cascade turbine is used to design the hydraulic turbine.

Turbine Cascade Flow Parameter
Design. Turbine motion is a dynamic process. e turbine is passively forced to rotate under the action of fluid. e speed of the turbine is related to the flow rate, turbine structure parameters, and the external load imposed on the turbine.
Turbine cascade flow parameter design refers to the calculation of the inlet and outlet liquid flow angle of the blade. Turbine stator and rotor blades have different flow angles at the inlet and outlet, and the law of liquid flow rate change in the turbine is also different, which will cause changes in performance parameters such as turbine output torque, speed, and power.
According to the unary flow theory [2], the absolute velocity, relative velocity, and circumferential velocity of the blade inlet and outlet under the nonimpact conditions constitute the blade inlet and outlet velocity triangle (as shown in Figure 1). e triangle of blade inlet and outlet speed is closely related to the characteristics of the turbine. According to the design parameters of the driving turbine and the geometrical dimensions of the turbine pair, the inlet and outlet speed triangles of the blades are determined, and the relationship between the angles of the triangles determines the parameters such as the flow angle.
Euler's formula for univariate flow theory is where T i , N i , and H i represent the torque, power, and conversion head of the single-stage turbine, respectively, Q, g, ρ, and R represent the inlet flow rate, gravity acceleration, working fluid density, and equivalent average radius of overflow, u is equivalent of the circumferential speed at the average flow radius, and C 1u and C 2u are the projections of the absolute inlet and outlet speeds C 1 and C 2 on the circumferential speed, respectively. In the speed triangle at the inlet and outlet of the rotor blade, there are the following equations: 2 Shock and Vibration where F is the cross-sectional area of the flow path perpendicular to the axial velocity, F � 2πBRφ; φ is the clogging factor affected by the thickness of the blade, which is generally 0.9. B is the radial width of the flow channel, and C Z is the axial flow velocity. Substituting (5) and (6) into (1), (2), and (3), respectively, we can obtain According to the import and export speed triangle, the fixed turbine, rotor inlet and outlet angles, and dimensionless coefficients have the following relationship: where C z is the axial velocity coefficient, C z is the axial flow velocity, and u is the peripheral speed, C z � C z /u; m a is the impact degree coefficient; C mu is the component of the average value of the absolute velocity of the entrance and exit in the direction of the circumferential velocity, m a � C mu /u; C u � C 1u − C 2u /u is circulation coefficient; as the circulation coefficient increases, the turbine output torque and turbine blade curvature also increase.

Blade Profile
Model. e turbine blade flow channel is constricted; that is, the flow channel from the blade inlet to the throat should be continuously and smoothly contracted, and the diffusion should be minimized from the throat to the trailing edge. e blade profile has a great influence on the flow field, mainly by changing the velocity and pressure distribution near the blade surface to change the flow field across the blade flow channel. It is generally expected that the blade surface speed and pressure distribution should change smoothly. e shape of the blade must ensure that the profile has a continuous curvature; that is, the profile of the pressure surface and the suction surface must have a continuous third derivative. e traditional blade pressure surface and suction surface curves are arc, parabola, double button, logarithmic spiral and hyperbola, and quadratic curve. ere are five main combinations of these curves when forming the turbine blade type: (1) completely composed of circular arcs, (2) composed of circular arcs and hyperbolic spirals, (3) composed of double curved lines, (4) composed of parabola, and (5) composed of parabola and arc. ese combined profiles often have discontinuous curvature points, which reduces turbine performance. In order to solve this problem, the blade pressure surface and suction surface profile should be in a single profile rather than in a combined profile and have a continuous Figure 1: Turbine cascade parameters. α 1 and α 2 : the outlet angle and inlet angle of turbine stator working fluid flow. β 1 and β 2 : inlet and outlet angles of turbine rotor working fluid flow. α 1k and α 2k : inlet structure angle and outlet structure angle of turbine stator. β 1k and β 2k : inlet structure angle and outlet structure angle of turbine rotor. c 1 : the absolute flow velocity at the inlet and outlet of the turbine rotor. c 2 : the absolute flow velocity at the turbine stator inlet and rotor outlet. c 1u : circumferential partial velocity of the absolute flow velocity at the inlet of the turbine rotor. c 2u : the circumferential component of the absolute flow velocity at the outlet of the turbine rotor. c 1z : the axial velocity component of the absolute flow velocity at the inlet of the turbine rotor. third-order derivative; high-order parameter splines, highorder polynomials, and high-order Bezier are the main selection objects.
In general, the shape of the blade is designed with a fifthdegree polynomial.
Establish the coordinates as shown in Figure 2. Let the pressure and suction surface profiles of the blade be y p � f(x) and y s � g(x), respectively; they have the following form: e first derivative is e second derivative is Let the coordinates, first and second derivatives of the first point on the pressure surface, be (x p1 , y p1 ), y p1 ′ and y p1 ″ ; let the coordinates, first and second derivatives of the last point on the pressure surface, be (x pn , y pn ), y pn ′ and y pn ″ ; let the coordinates, first and second derivatives of the first point on the suction side, be (x s1 , y s1 ), y s1 ′ and y s1 ″ ; let the coordinates, first and second derivatives of the last point on the suction side, be (x sn , y sn ), y sn ′ and y sn ″ . en, we can obtain the following two linear equations: Solve the linear equations (12) and (13) to obtain the blade pressure surface and suction surface profile.
Using the parameters in Table 1, combined with (12) and (13), we obtain that the pressure surface equation is e suction surface equation is

Fluid Motion Equations.
A flowing fluid follows the law of conservation of mass, momentum, and energy. However, it is very difficult to solve the velocity field and pressure field in complex flows accurately. e computational fluid dynamics (CFD) method can be used to solve the approximate solution to satisfy the engineering application. e accuracy of CFD methods has been verified in many fields [30][31][32]. is paper employed numerical simulation tools to demonstrate the results of the theoretical analysis and to study the variation of geometric parameters. e working medium of the turbine is liquid. Generally, the liquid is regarded as an incompressible fluid. e basic  e absolute motion continuity equation and motion in the cylindrical coordinate system are used for the stator, and the rotor adopts the relative motion continuity equation and motion equation under the cylindrical standard system.
Using the traditional numerical method to solve the fluid flow in the turbine stator and rotor across the blade flow channel is very complicated; generally, there is no analytical solution; as a result, CFD method is used.

Simulation Turbulence Model.
A flow simulation was carried out using the software of ANSYS FLUENT 19.0, which uses the finite-element numerical method for solving the Reynolds-averaged Navier-Stokes equations by means of the pressure-based solver. e whole 3D geometry model of the turbine was established using ANSYS Workbench, and grids were generated using the preprocessing software ICEM-CFD. e k-ε realizable model is selected. e transport equations of turbulent kinetic energy and dissipation rate are as follows: where k is the turbulent kinetic energy; ε is the turbulent dissipation rate; and σ k and σ ε are the Prandtl number corresponding to the turbulent flow energy k and the turbulent dissipation rate ε, σ k � 1.0 and σ ε � 1.2; x i and x j are the position coordinate component; u i is the velocity coordinate component; v is the molecular kinematic viscosity coefficient; G k is the average velocity gradient which causes the turbulent kinetic energy k generation term; G b is the buoyancy induced turbulent kinetic energy generation term; Y M is the effect of compressible turbulence fluctuation on the total dissipation rate; and μ t is the turbulent viscosity coefficient, C 1ε � 1.44 and C 2 � 1.9.
In order to better simulate the flow near the wall, a boundary layer needs to be generated when meshing, and the height of the first layer of the boundary layer is determined by the following formula: where where ρ is the fluid density; μ is the dynamic viscosity coefficient; U is the free flow velocity; L is the characteristic length; and y + is selected according to the standard wall function. e turbine simulation model and grid are shown in Figure 3. In the 3D model, the turbine wall is placed in the rotating domain, and the upper and lower parts are static domains.
e method of combining structured grid and unstructured grid is adopted for grid generation. e total number of grids is about 400000, skewness is less than 0.89, aspect ratio is less than 5 : 1, and growth rate is 1.13. e number and quality of grid meet the requirements of simulation accuracy.

Results and Discussion
e stress contour and streamline diagram of the turbine when the turbine speed is stable are shown in Figure 4. It is known from the stress distribution of the blade that the fluid pressure above the turbine is large, the turbine bears a large pressure on the upper part of the blade, and the vortex is formed under the fluid domain of the turbine. Conventional turbine mechanisms, such as wind turbines and water turbines, often use the momentum-blade theory to analyze the forces on the turbine. In the analysis, the momentum change caused by fluid impact is usually ignored; that is, the impact force is ignored. However, relevant researchers pointed out that when the turbine size becomes smaller, the force on the turbine due to fluid impact cannot be ignored. erefore, the turbine force is composed of three parts: impact force, pressure drag, and friction resistance. e impact force is caused by the change of fluid momentum, while the pressure drag is caused by the pressure difference between the front and back of the blade, and the friction resistance is caused by the viscosity of the fluid. Whether the turbine rotates or not depends on whether the torque of the fluid acting on the turbine is greater than the external load. When the torque generated by the fluid is less than the external load, the turbine is stationary; when the torque generated by the fluid is greater than the external load, the turbine starts to move. e following will analyze the influence of four factors, namely, external load, blade installation angle, flow rate, and blade number, on the movement and force of turbine. Shock and Vibration 5

Static Analysis.
e moment and impact force on the turbine when it is stationary are called starting torque and initial impact force. e force on the turbine at rest determines whether the turbine can overcome the load and then moves. e influence of different factors on starting torque and initial impact force of the turbine is shown in Figure 5(a). e curve in the graph is fitted nonlinearly. As shown in Figure 5(a), under the condition of keeping the turbine structure design and fluid parameters unchanged, the change of external load on the turbine will not affect the starting torque and initial impact force when the turbine is stationary. e fitting curve is a straight line with a value of 7 N·m.
As shown in Figure 5(b), the fitting results show that the change of the installation angle has four-time relationship with the starting torque and initial impact force. With the increase of the installation angle, the impact torque and initial impact force of the turbine have an increasing trend. When the installation angle of turbine blade increases, the fluid velocity and the angle of attack of the blade will change, which will lead to the change of momentum and flow field around the blade, resulting in the change of starting torque and initial impact force of turbine. e increase of flow rate means the increase of inlet flow rate. e starting torque and initial impact force obtained under different flow conditions are fitted. It is found that the relationship between the inlet flow rate and the starting torque and impact force of turbine is quadratic, as shown in Figure 5(c). With the increase of flow rate, namely, the increase of inlet flow rate, the torque of turbine increases and amplitude rises; meanwhile, the absolute value of impact force increases with the increase of velocity and the amplitude also rises.
As shown in Figure 5(d), the starting torque and initial impact force obtained under different blade numbers are    Shock and Vibration also fitted, and it is found that the number of blades is approximately quadratic with the starting torque but linearly with the impact force.

Dynamic
Analysis. Now, the overall analysis of the movement process of the turbine is carried out. e influence of external load, blade installation angle, flow rate, and blade number on the external torque, impact force, and angular velocity of the turbine movement process is shown in Figures 6-8 When the load of the turbine is large, the torque of the fluid acting on the turbine cannot overcome the external load. e turbine is stationary, the angular velocity is always zero, and the impact force and torque are maintained. When the torque generated by the fluid is greater than the load, the turbine starts to move, and the change of the angular velocity of the turbine is directly related to the impact torque of the turbine.
When the turbine starts to rotate from static state, its angular velocity increases continuously, which causes the angle of attack of the fluid acting on the turbine blade to change. On the one hand, it causes the flow field around the blade to change and the projection of the resistance and lift force on the blade to change in the vertical direction; on the other hand, it causes the momentum change when the fluid impacts the turbine blade, and finally impact force and hydraulic moment of the turbine change with the attack angle. e speed is constantly changing. When the torque generated by the fluid is equivalent to the external load, the total external moment of the turbine is zero and the turbine speed no longer changes. As a result, the flow field tends to be stable, and the angular velocity, impact torque, and impact force of the turbine are stable.
(1) Turbine Angular Velocity Variation. e influence of external load, blade installation angle, flow rate, and blade number on turbine angular velocity is shown in Figure 6. Installation angle (°) Starting torque Initial impact force Under different conditions, the change rule of turbine angular velocity with time is the same, turbine angular velocity increases with time, and the final speed is constant. When the external load is 8N·m, the external load is greater than the starting torque of the turbine, so the turbine stays still. e larger the installation angle of turbine, the larger the stable speed of turbine, and the longer the stability time; the larger the inlet flow, the larger the stable speed of turbine, and the shorter the stability time. However, the increase of the inlet flow leads to the increase of the erosion intensity of the drilling fluid to the turbine, and ultimately the service life of turbine is reduced. e more the number of blades, the higher the stability speed of turbine, and the shorter the stability time. e increase in the number of blades will affect the flow distribution. e more the blades, the more stable the fluid flow field. However, the high number of blades leads to an increase in rotation speed, and the pressure difference between the upper and lower turbines increases, which is not conducive to the flow of drilling fluid.
(2) Turbine Impact Force Variation. e influence of flow rate, external load, blade installation angle, and blade number on turbine impact force is shown in Figure 7. e impact force of turbine has the same change trend with time, which increases first and then stabilizes. e change of turbine impact force is closely related to the change of turbine rotation speed. With the increase of turbine rotation speed, the absolute value of impact force becomes smaller.  Shock and Vibration With the increase of flow rate, external load, installation angle, and number of blades, the absolute value of impact force decreases when it is stable. When the external load of the turbine is 8 N·m, the impact force of the turbine remains unchanged, because the starting torque of the turbine is smaller than the external load of the turbine; that is, the turbine stays still. In CFD analysis and calculation, the basic fluid parameters are density 1000 kg/m 3 , viscosity 10 mPa·s, and flow 30 L/s. e predicted mechanical characteristic curve is shown in Figure 9.  Figure 10 shows the streamline diagram of liquid phase and solid phase of the turbine. In the simulation, the solid phase and the liquid phase have the same inlet velocity, but when the turbine works, the solid phase and the liquid phase have different velocity distribution, and the flow velocity distribution range of the solid phase is smaller than that of the liquid phase. e pressure drop of the turbine depends on the change of fluid momentum caused by the impact of solid and liquid phases on turbine blades. Under different solid phase conditions, the flow field of solid-liquid two-phase caused by the turbine rotation is different, and the pressure drop and torque of the turbine will change. e effect of different solid content on the performance of the turbine is shown in Figure 11. With the increase of solid content, the proportion of solid phase in the solidliquid two-phase flow increases, and more solid particles in the flow field interact with the blades of the turbine. With the increase of solid content, the pressure drop and torque increase continuously, and then the input and output power increase. It can be seen from the figure that the efficiency curve of the turbine shows an upward trend, but its linear slope decreases gradually, which shows that the working efficiency of the turbine increases with the increase of solid content and increase amplitude of efficiency decreases. e effect of solid density on turbine performance is  shown in Figure 12. Although the inlet velocity is constant, the kinetic energy of particles increases with the increase of solid density. e impact of solid on turbine is more severe. It is shown that the torque, pressure drop, and input and output power of the turbine increase linearly with the increase of solid density. e efficiency increases with the increase of solid density, but the increase amplitude also decreases. e effect of solid particle size on turbine performance is shown in Figure 13. In the case of constant solid density   and solid content, the increase of solid particle size substantially leads to the decrease of the number of particles entering the fluid domain. Although the inlet velocity is constant, the kinetic energy of single particle increases with the increase of particle size. e change of solid particle size will cause the change of particle number and kinetic energy, which will lead to different results of two-phase flow impacting turbine blades. According to the figure, with the increase of solid particle size, the pressure drop and torque of turbine first increase; with the further increase of particle size, the pressure drop and torque of turbine tend to be stable. e change trend of turbine input and output power is the same as that of pressure drop and torque. e efficiency increases with the increase of solid particle size, and the increase amplitude decreases.

Conclusion
Turbine is a common mechanical structure in drilling industry, which is usually used in downhole turbogenerator, turbodrill, and so on. With the continuous development and innovation of computational fluid dynamics (CFD) theory and numerical methods, CFD method has been become an important means for the study of downhole turbine performance. At present, the research on downhole turbines is static. In this paper, by using FLUENT software and sliding mesh method, UDF file is compiled to simulate the dynamic performance of downhole turbine. e effects of external load, number of blades, flow rate, and blade installation angle on the hydrodynamic performance of downhole turbine are studied, and the start-up characteristics and movement process of downhole turbine are analyzed. e results show that when the start-up moment is less than the load, the turbine is stationary; when the start-up moment is greater than the load, the turbine starts to rotate. In the process of the turbine's rotation, the force moment of the turbine decreases, while the speed and impact force increase continuously. When the rotational velocity is stable, the torque produced by the fluid is balanced with the load, and the impact force is stable. Solid phase of drilling fluid will affect the performance of downhole turbogenerator. Based on Eulerian multiphase flow model and FLUENT software, three-dimensional numerical simulation of turbogenerator is carried out. e effects of solid content, particle size, and solid density of drilling fluid on the performance of downhole turbogenerator are studied. e simulation results show that, with the increase of solid phase content, the pressure drop, torque, input and output power, and efficiency of turbogenerator increase; with the increase of solid particle size, the pressure drop, torque, and input and output power of turbogenerator increase first and then stabilize, while the efficiency increases gradually and the increase decreases; with the increase of solid density, the pressure drop, torque, input and output power, and efficiency of turbogenerator increase linearly, and the efficiency increase decreases gradually.
Data Availability e data in the tables used to support the findings of this study are included within the article. e data in the figures used to support the findings of this study are available from the corresponding author (wjzym286@163.com) upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.