Calculation of the Incompressible Viscous Fluid Flow in Piston Seals of Piston Hybrid Power Machines

The article considers the calculation of the flow of a viscous incompressible fluid in piston seals of piston hybrid power machines. The most widely used and effective seals are considered: a smooth gap seal and a step-type gap seal, and—based on the references analyses—the initial boundary conditions for their calculation are determined. The laminar and turbulent flows in gap seals are calculated based on the well-known analytical relationships, experimental studies and flow models, including the k-ε, Menter’s Shear StressTransport (SST) and Reynolds Stress (RSM) turbulence models. The effectiveness of using each model to determine average velocities, flow rates, and velocity plots in the cross section of a gap seal, as well as the adequacy of the description of known physical laws, is estimated. The results proved that the RSM turbulence model is better for the gap seals of different types under different modes of motion.


Introduction
Considerable energy is spent on the drive of piston compressors and piston pumps, including radial piston and axial piston, considering their widespread use in the automotive, machine-tool, aviation and other industries. As a result, the task of increasing their efficiency and effectiveness is not new, but still relevant and timely.
One of the main ways to increase dramatically the efficiency and economy of gas compression is to combine a piston compressor with a piston pump and to obtain a piston hybrid power machine (PHPM) in which gas and liquid are compressed and moved by changing the volume of the working chambers (see Figure 1). Due to the intensive cooling of the suction and compressed gas and the elimination of leaks and leakage losses, the indicator efficiency increases and work is reduced; by eliminating dead space in the compressor section by filling it with fluid, the feed coefficient and productivity of the compressor section are increased; due to the fluid in the gap seal, the work of friction forces in the cylinder-piston group decreases; due to leakage of liquid into the working cavity of the pump section from the gap in the piston seal, the cavitation reserve increases [1]. friction forces in the cylinder-piston group decreases; due to leakage of liquid into the working cavity of the pump section from the gap in the piston seal, the cavitation reserve increases [1].
Thus, the piston seal through which the fluid flows from the pump section to the compressor section and back is a key factor providing all positive effects in the PHPM.
Over the years, gap seals have been heavily used in hydraulic, oil, and fuel systems, as well as in pumps, hydraulic motors, distribution and control devices, and power cylinders. Thus, the piston seal through which the fluid flows from the pump section to the compressor section and back is a key factor providing all positive effects in the PHPM.
Over the years, gap seals have been heavily used in hydraulic, oil, and fuel systems, as well as in pumps, hydraulic motors, distribution and control devices, and power cylinders.
Its widespread use is due to constructive simplicity and efficiency. A gap seal is a gap between two mating parts in the simplest case. There is fluid flow observed in the gap seal under the action of a pressure differential between the sealing cavity and the low-pressure cavity. There will be two flow modes observed: laminar and turbulent, depending on the pressure drop, physical properties of the fluid and the geometrical dimensions of the gap seal.
In the general case, the gap between the two mating parts can be straight, expanding, tapering and of another shape, with and without labyrinth grooves.
Analytical solutions have been obtained for the distribution of the diagrams of velocities, average velocity and flow rate for the laminar flow in the simplest slots. So, formulas for determining the flow rate were obtained for flat and annular slots, called Hagen-Poiseuille formulas [2]. In the general case, the Navier-Stokes equation and the continuity equation are used to describe the laminar flow in the gaps [3,4]. Due to the fact that solving complete Navier-Stokes equations causes great mathematical difficulties, the equations are simplified by introducing the following assumptions: fluid flow is steady, isothermal between solid impermeable walls; -the mass terms of the equation of motion are negligible (they should be taken into account only when the fluid flows on a free surface or when the fluid density is inhomogeneous); -inertia forces are also negligible [4,5] (inertia forces can be neglected if the reduced Reynolds where v is determining speed; ν is kinematic viscosity; l is seal length; h is slot height); -∂ 2 v x ∂x 2 value is less then ∂ 2 v x ∂y 2 value and it can be neglected; -transverse pressure gradient is negligible; -all members containing z coordinate (coordinate in the transverse direction) and derivatives with respect to z are discarded.
As a result of the assumptions made, the system of equations is significantly simplified and can be used to obtain analytical solutions.
Analytical calculation is impossible in a turbulent fluid flow, which was not an overcome problem for their calculation with a Reynolds number more than 1000.
Currently, due to the development of computer technology and numerical methods of analysis, two-parameter semi-empirical turbulence models using Reynolds averaging and turbulent viscosity are widely used. These models include k − ε, k − ω and combinations thereof, for example SST [6][7][8][9][10][11][12][13]. There are also many other turbulence models based on an additional tensor in the Reynolds equations, which includes so-called Reynolds stresses. A detailed analysis of turbulence models will be given in the article below when considering turbulent flow.
It should be noted that all developed models require experimental verification. The article will also present the results of experimental studies.
Thus, the main objectives of this article can be formulated as follows: 1.
To analyze gap seals used in hybrid power machines and select the most structurally simple from them and effective in operation.

2.
To determine the basic geometric and operational parameters for the selected types of gap seals.

3.
To determine the flow patterns in them based on existing recommendations.

4.
To calculate the laminar and turbulent flow modes. To determine the most suitable turbulence models for the turbulent flow mode and conduct a comparative analysis of the effectiveness of their application based on a comparison with the results of experimental studies and known physical concepts.
To provide recommendations on the applicability of existing turbulence models for calculating fluid flows in gap seals of hybrid power machines.

Analysis of the Main Types of Piston Seals Used in PHPM
The classic piston seal connecting the compressor and pump sections is a smooth gap of constant clearance along the length formed by the piston and cylinder [14,15], while the piston has a concentric or eccentric arrangement. When using a smooth gap seal, the flow rate of the fluid in the forward and reverse direction is the same. In accordance with the previously performed researches [16,17], the PHPM can operate in three modes, with the main mode being the priority flow of fluid from the pump fluid to the compressor. In this case, the piston seals must have different hydraulic resistance in the forward and return directions. When the fluid moves from the pump section to the compressor, the hydraulic resistance must be less than when the fluid moves from the compressor to the pump.
This effect can be achieved in various constructive ways. [18][19][20][21][22]. We consider the main ones. There was considered a gap seal (see Figure 2), which has protrusions on the piston body and is profiled to provide difference in the flow rates of the fluid in the forward and return directions [19,23]. The multidimensional nonlinear optimization of the basic geometric parameters establishes their optimal values and achieve a ratio of fluid flow G p /G c = 1.3 ÷ 1.4, where G p is mass flow of fluid per cycle from the pump section to the compressor; G c is mass flow of fluid per cycle from the compressor section to the pump. It should be noted that this ratio is small and it was necessary to increase it. This was achieved in subsequent studies.

Analysis of the Main Types of Piston Seals Used in PHPM
The classic piston seal connecting the compressor and pump sections is a smooth gap of constant clearance along the length formed by the piston and cylinder [14,15], while the piston has a concentric or eccentric arrangement. When using a smooth gap seal, the flow rate of the fluid in the forward and reverse direction is the same. In accordance with the previously performed researches [16,17], the PHPM can operate in three modes, with the main mode being the priority flow of fluid from the pump fluid to the compressor. In this case, the piston seals must have different hydraulic resistance in the forward and return directions. When the fluid moves from the pump section to the compressor, the hydraulic resistance must be less than when the fluid moves from the compressor to the pump. This effect can be achieved in various constructive ways. [18][19][20][21][22]. We consider the main ones.
There was considered a gap seal (see Figure 2), which has protrusions on the piston body and is profiled to provide difference in the flow rates of the fluid in the forward and return directions [19,23]. The multidimensional nonlinear optimization of the basic geometric parameters establishes their optimal values and achieve a ratio of fluid flow Gp/Gc = 1.3 ÷ 1.4, where Gp is mass flow of fluid per cycle from the pump section to the compressor; Gc is mass flow of fluid per cycle from the compressor section to the pump. It should be noted that this ratio is small and it was necessary to increase it. This was achieved in subsequent studies. A piston seal of a step-type form has been studied (see Figure 3). It consists of a gap with small δ1 and large δ2 clearances and it works as follows: when the piston moves from the top dead center (TDC) to the bottom dead center (BDC), the fluid is compressed and injected in the pump section and the gas expands and sucks in the compressor. Piston seal length with low clearance 1 decreases and with a large clearance 2 increases, and in general, the hydraulic resistance of the gap along the piston is constantly decreasing. During the piston stroke from BDC to TDC, gas compression and injection are in the compressor section, and expansion and absorption are in the pump section. The fluid from the compressor section through the gap seal enters the pump section. It should be noted that in this case, the hydraulic resistance of the gap increases due to the fact that the length of the gap with a small clearance δ1 increases, and the length of the gap with a large clearance δ2 decreases. Thus, for a  A piston seal of a step-type form has been studied (see Figure 3). It consists of a gap with small δ 1 and large δ 2 clearances and it works as follows: when the piston moves from the top dead center (TDC) to the bottom dead center (BDC), the fluid is compressed and injected in the pump section and the gas expands and sucks in the compressor. Piston seal length with low clearance 1 decreases and with a large clearance 2 increases, and in general, the hydraulic resistance of the gap along the piston is constantly decreasing. During the piston stroke from BDC to TDC, gas compression and injection are in the compressor section, and expansion and absorption are in the pump section. The fluid from the compressor section through the gap seal enters the pump section. It should be noted that in this case, the hydraulic resistance of the gap increases due to the fact that the length of the gap with a small clearance δ 1 increases, and the length of the gap with a large clearance δ 2 decreases. Thus, for a given type of gap seal, the G p /G c value can be 4 ÷ 5. Even better results can be expected from a gap seal, which does not change the length along the piston, but the clearance [22]. However, studies of the PHPM with such gap seals were not carried out. given type of gap seal, the Gp/Gc value can be 4 ÷ 5. Even better results can be expected from a gap seal, which does not change the length along the piston, but the clearance [22]. However, studies of the PHPM with such gap seals were not carried out. Thus, to study the fluid flow in a piston seal, it is advisable to consider two types of gap seals: a smooth gap with a constant clearance and a gap with a stepwise change in the clearance.

Theoretical Research
The analysis of theoretical and experimental studies of PHPM is to determine the most characteristic values of geometric and operational parameters: 1. The suction pressure in the compressor and pump section was always assumed to be the same and equal to 0.1 MPa.
2. The discharge pressure in the compressor and pump sections was in the range from 0.3 to 1.0 MPa. To conduct theoretical research and to compare the flow rates of fluid in the forward and return directions, we assume that Pdp = Pdc = 1МPа (where Pdp is discharge pressure in the pump section; Pdc discharge pressure in the compressor section).
3. The piston diameter of the studied PHPM was from 0.04 to 0.06 m. We assume that the piston diameter is 0.04 m in these studies. 4. The length of the piston actually determines the length of the piston seal. It varies in different PHPM designs. We assume that the length of the piston is 0.04 m in these studies.
5. The magnitude of the radial clearance is one of the most important parameters. It was shown in [24] that the rational values for a smooth gap seal for a crosshead free PHPM are (20 ÷ 30) μm. It was proved in [25] that rational clearances δ1 and δ2 are (30 ÷ 55) μm and (80 ÷ 110) μm for a steptype piston seal in a crosshead PHPM respectively. In this paper, we assume that δ1 is 25 μm, and δ2 is 50 μm, i.e., those values are close to rational.

Laminar Fluid Flow
For the steady flow of incompressible viscous fluid, we obtain the equations for the distribution of fluid in the annular gap, the average velocity and flow rate for the fixed and moving walls with a concentric and eccentric arrangement of the piston from the equilibrium of the acting forces [24]. Thus, to study the fluid flow in a piston seal, it is advisable to consider two types of gap seals: a smooth gap with a constant clearance and a gap with a stepwise change in the clearance.

Theoretical Research
The analysis of theoretical and experimental studies of PHPM is to determine the most characteristic values of geometric and operational parameters: The suction pressure in the compressor and pump section was always assumed to be the same and equal to 0.1 MPa.

2.
The discharge pressure in the compressor and pump sections was in the range from 0.3 to 1.0 MPa. To conduct theoretical research and to compare the flow rates of fluid in the forward and return directions, we assume that P dp = P dc = 1 MPa (where P dp is discharge pressure in the pump section; P dc discharge pressure in the compressor section).

3.
The piston diameter of the studied PHPM was from 0.04 to 0.06 m. We assume that the piston diameter is 0.04 m in these studies.

4.
The length of the piston actually determines the length of the piston seal. It varies in different PHPM designs. We assume that the length of the piston is 0.04 m in these studies.

5.
The magnitude of the radial clearance is one of the most important parameters. It was shown in [24] that the rational values for a smooth gap seal for a crosshead free PHPM are (20 ÷ 30) µm.

Laminar Fluid Flow
For the steady flow of incompressible viscous fluid, we obtain the equations for the distribution of fluid in the annular gap, the average velocity and flow rate for the fixed and moving walls with a concentric and eccentric arrangement of the piston from the equilibrium of the acting forces [24].
The fluid flow rate without taking into account the frictional motion and for the concentric arrangement of the piston in the cylinder is determined as: where d is the cylinder diameter; µ is dynamic fluid viscosity; ∆p is piston differential pressure; δ is radial clearance. The average velocity is determined as: Reynolds number is defined as: It should be noted that the critical number Re cr the laminar flow mode is carried out in a smooth gap is within 1000 ÷ 1200.
Using Equations (2) and (3), the value of the radial clearance the laminar mode of fluid flow will be observed is defined as: For accepted values n and ∆p, and also considering that the studied fluid is water, and at T = 293 K = 998.2 kg/m 3 and µ = 10 −2 Pa s, as well as considering that Re cr = 1000, we obtain that δ cr ≤ 80.5 µm, laminar mode will be realized with radial clearance values up to 80.5 µm. To calculate the step-type gap seal in [25], the equality of the combined costs is used Q 1 and Q 2 through its narrow and wide parts, while it is believed that the local resistance due to the step-type change in the radial clearance is negligible.
Using equality of expenses Q 1 and Q 2 , pressure P3 is written in the section of step-type clearance as where 1 and 2 are seal lengths with clearance δ1 and δ2.
The fluid flow rate in this case is determined as: Reynolds number is determined as: Values 1 and 2 are interconnected by the following equality 1 + 2 = n , where n is the piston length.

Turbulent Fluid Flow
In general case to describe the unsteady flow of viscous compressible fluid, the following system of equations can be used, including the continuity equation, the equation of motion in the form of the Navier-Stokes equation, the energy conservation equation: The Navier-Stokes Equation (10) is Newton's second law written for an infinitely small volume of fluid taking into account pressure forces and viscous friction forces. The left side of the equation includes local and convective acceleration, and the right side includes the forces of pressure and viscous friction, respectively grad(p) and Div(τ). Thermodynamic pressure p is determined as p = 1 3 σ xx + σ yy + σ zz , and the stress tensor τ is determined on the basis of the generalized Newton's law of friction: Equation (11) is an energy conservation equation for an infinitely small volume, which is written taking into account the braking heat expressed in terms of the braking enthalpy H, which is associated with static enthalpy as follows: To determine steady flow and incompressible fluid, taking into account constant temperature and, accordingly, constant thermophysical properties, the system of Equations (9)-(11) is transformed to the following form: -Navier-Stokes equation: where V is velocity vector.
The system of Equations (13) and (14) in the general case has no analytical solution and can be solved only numerically. Currently, three main methods are used to solve such systems of equations: the finite difference method, the finite volume method, and the finite element method. The finite volume method is most widely used due to the high computational stability. The idea of this method is to replace infinitely small volumes of fluid with volumes of finite sizes-cells. Equations (13) and (14) are written for each created cell, and the totality of all cells is called the computational grid. The theory of turbulent motions proves that turbulent fluid motions consist of many vortex formations having different geometric dimensions [26,27]. As a result, the cell sizes used for the numerical solution of the equations should be much smaller than the smallest turbulent formation. This condition determines the number of cells required to solve the problem. As a rule, these are tens and hundreds of millions of cells. The cell size of the computational grid is chosen much larger than the size of small-scale turbulent formations, as existing computers have insufficient performance to solve similar problems of fluid mechanics and it is necessary to reduce the number of cells [28]. In this way the number of cells in the grid can be significantly reduced. In this case, a certain number of turbulent structures will fall inside each of the grid cells, and, therefore, their characteristics cannot be determined on the basis of Equations (13) and (14), i.e., additional equations are required to calculate turbulence characteristics. To date, several empirical methods for calculating turbulence characteristics have been developed. Reynolds offered the most widely used method which consists in averaging variables.
The proposed averaging of Reynolds does not solve the problem, because it adds six additional unknowns to the system of equations, which are called "Reynolds stress". Using the Boussinesq hypothesis simplifies significantly the Navier-Stokes Equation (14) and reduces the number of unknowns characterizing turbulence from six to one, turbulent viscosity. Turbulent viscosity is not a constant value, like molecular viscosity, but depends both on time and on spatial coordinates, and its value is calculated for each cell of the design string at each step of integration over time. The calculation of the turbulent viscosity is carried out using a semi-empirical model of turbulence. At present, a large number of such models have been created; however, two-parameter semi-empirical turbulence models are most widely used. This is due to the fact that these models are the most stable, although they are less accurate than a number of other models.
The analysis of existing turbulence models revealed three models; their use is expedient for calculating the flow in a gap seal of a piston hybrid power machine.

Turbulence Model k − ε
Turbulence model k − ε is one of the most widely used models of turbulence and it consists of two equations: kinetic energy of turbulence k and dissipation of kinetic energy ε. The main idea underlying this model is to describe the processes of nucleation, transfer, and dissipation of turbulence. To implement this idea, the following differential transport equations k and ε were used: In these equations k is kinetic energy of turbulence; ε is turbulence energy dissipation; P k is generation (generation) of turbulence energy, which is calculated by the following dependence for an incompressible fluid: where P kb and P εb are terms that take into account the effects of gravity which are usually neglected. Joint solution of Equations (15) and (16) is to calculate the distribution k and ε over the computational domain at a given point in time, then the value of turbulent viscosity for each computational cell is calculated. A number of empirical constants are used in the turbulence model k − ε, its value is determined by calibrating the model according to a series of field experiments. References prove that the model k − ε calibrated for developed turbulent fluid flows and can mainly be used to describe the free flow of fluid at a sufficient distance from the solid walls.
The main disadvantages of the model include: 1.
Instability of work with flows having a large pressure gradient and a powerful degree of turbulence. 2.
This model gives an overestimated value of kinetic energy at high speeds of movement of the working fluid.
3.4. Turbulence Model k − ω k − ω model was developed to eliminate k − ε model flaws. It was developed to calculate fluid motion taking into account the presence of solid walls. Turbulent viscosity in the model k − ε is determined as follows: It should be noted that the k − ω model equation has a similar physical meaning and structure with equations k − ε models and differ in closure constants and terms describing the generation and dissipation of variables k, ω and ε.
Despite all the positive properties, k − ω model of turbulence, it is unacceptable in a number of cases, for example, when describing the developed turbulent motion away from solid walls and having direct contact with them. In this case, high accuracy is required for calculating the turbulence characteristics both near and far from solid walls. An example of such cases is the problem of separation of the boundary layer from a solid wall. In these cases when there is separation of the boundary layer, k − ε turbulence models are poorly described as well as k − ω models.

Turbulence Models SST
Menter developed the Shear StressTransport (SST) turbulence model that includes both k − ε and k − ω turbulence models to increase the accuracy of the description of turbulent flows in which there is a probability of boundary layer separation. This became possible due to the high similarity of these models' equations. With a certain combination of empirical constants SST turbulence model can be reduced to k − ε model or to k − ω model.
SST turbulence model includes two differential equations describing the transfer of turbulence energy k and the intensity of its dissipation ω. The model uses two functions F 1 and F 2 -these are the transition functions between k − ε and k − ω turbulence models. Far from a solid wall, these functions turn into 0, while the SST turbulence model comes down to a standard model k − ε. Function F 1 and F 2 tend to 1 inside the boundary layer, which corresponds k − ω model. A smooth transition from one model to another takes place in the space between the boundary layer and the developed turbulence model. All empirical constants of SST turbulence models are calculated by recalculating the corresponding constants k − ε and k − ω models using the function F1 followed by addiction: where a is SST constant of turbulence model; a 1 and a 2 are constant k − ε and k − ω models respectively. For the successful application of the SST turbulence model, it is necessary to choose the right cell size directly near the solid wall. The SST model does not use the transport equation for turbulent shear stresses, but calculates the shear stresses themselves using the Boussinesq hypothesis based on the application of the concept of turbulent viscosity. The value of turbulent viscosity is determined similarly to the model k − ω, but with the introduction of a limiter. This representation of turbulent viscosity is to solve problems with the appearance of separated flow flows with changes in the sign of the pressure gradient near the wall.
There are many other turbulence models, in particular, a group of RSM models based on the additional tensor in the Reynolds equations which includes the so-called Reynolds stresses or turbulent stresses. These are RSM, Explicit Algebraic Reynolds Stress Model (EARSM) and other models. The main advantage of such models is the inclusion of turbulence anisotropy which takes into account the turbulence characteristics in different directions of space. The disadvantages include a large number of calibration empirical coefficients which significantly reduce the versatility of this model and the inability to apply it without changing the coefficients to other turbulent flows. The analysis of references, based on material costs, decision time and quality of the results allowed focusing on three main models: k − ε, SST and RSM.
It should be noted that a complex of experimental studies was carried out in [2] to determine the flow rate of liquid through smooth slots in a turbulent flow mode, its results were processed and fluid flow formulas were obtained for concentric and eccentric arrangement of the piston.
For the concentric piston arrangement in the cylinder we have: where ν = µ/ρ is kinematic viscosity. For a step-type seal using the equality of volumetric flow rates in the narrow and wide parts of the piston seal we obtain the following Equation for pressure p 3 : and flow: Using Equations (20) and (22), we can evaluate the reliability of the turbulence models considered above.

Data Analysis
A numerical experiment will be conducted on the applicability of the turbulence models discussed above, as well as the Laminar model, for calculating smooth and step-type gap seals in laminar and turbulent modes in a piston hybrid power machine. We calculate the flows in gap seals using the ANSYS CFX application software package which includes turbulence model data.

Implementation Features
Assuming that the problem under consideration is axisymmetric, a gap seal segment with an angular width was considered 0.01 0 to reduce the calculation time and reduce the number of cells. When determining the total flow rate in a gap seal, the flow rate determined in this segment is multiplied by 36,000.
The segments of gap seals with constructed meshes are shown in Figures 4 and 5. The construction of the finite elements was carried out automatically in a smooth gap seal, and it was carried out automatically with manual correction in a step-type gap seal at the points of rotation on 90 0 .
Using Equations (20) and (22), we can evaluate the reliability of the turbulence models considered above.

Data Analysis
A numerical experiment will be conducted on the applicability of the turbulence models discussed above, as well as the Laminar model, for calculating smooth and step-type gap seals in laminar and turbulent modes in a piston hybrid power machine. We calculate the flows in gap seals using the ANSYS CFX application software package which includes turbulence model data.

Implementation Features
Assuming that the problem under consideration is axisymmetric, a gap seal segment with an angular width was considered 0.01 0 to reduce the calculation time and reduce the number of cells. When determining the total flow rate in a gap seal, the flow rate determined in this segment is multiplied by 36,000.
The segments of gap seals with constructed meshes are shown in Figures 4 and 5. The construction of the finite elements was carried out automatically in a smooth gap seal, and it was carried out automatically with manual correction in a step-type gap seal at the points of rotation on 90 0 .  We conduct a sequential review of gap seals.

Laminar Flow
Diagrams of velocity distribution along the cross section of the gap are presented in Figure 6. We conduct a sequential review of gap seals.

Laminar Flow
Diagrams of velocity distribution along the cross section of the gap are presented in Figure 6. We conduct a sequential review of gap seals.

Laminar Flow
Diagrams of velocity distribution along the cross section of the gap are presented in Figure 6. Analyzing the presented data, we have the following conclusions: X, mm V, m/s Analyzing the presented data, we have the following conclusions: 1.
The velocity distribution diagrams of the RSM and SST turbulence models are closest to the theoretical velocity distribution determined in accordance with [4,5]. The discrepancy in determining the maximum speed between them is approximately equal to 1 m/s, i.e., 14.2%.

2.
Speed distribution by model k − ε is significantly better than laminar model. The discrepancy in determining the maximum speed model k − ε is 2 m/s, i.e., 28.4%. The discrepancy in determining the maximum speed of the laminar model is 3.5 m/s, i.e., 50%. Table 1 presents the average fluid velocities in the gap and the Reynolds numbers determined by different flow models. The presented results lead to the following conclusions: 1. The Reynolds numbers determined by different models are in the range from 120 to 210, which allows supposing that there is a laminar flow mode in the gap seal. It should be noted that this conclusion could be made from the previously considered results (velocity diagrams).

2.
The maximum average velocities of 4.234 and 4.174 m/s are observed using the SST and RSM turbulence models. These average velocities exceed significantly (almost by 0.5 m/s) the average velocities determined by calculation for the laminar fluid flow in the gap [3], and according to the turbulence model k − ε.

3.
The lowest average velocity of 2.404 m/s is shown by the laminar fluid flow model which is almost 1.3 m/s less than the theoretically calculated. 4. k − ε model shows the closest to theoretically calculated values by the definition of average velocity. The discrepancy in determining the average velocity, compared with the theoretically determined one, is 0.033 m/s or 0.88%, i.e., less than 1%. The discrepancy in the definition of average velocities is equivalent to the discrepancy in the definition of volumetric costs presented in Table 2.
Thus, summarizing the analysis of the laminar flow in the smooth gap of the PHPM, we can conclude that it is preferable to use SST and RSM to determine the maximum values of velocities, and the k − ε model to determine average velocities and costs.

Turbulent Flow
The turbulent flow, in contrast to the laminar flow, is characterized by the presence of vortices and intensive mixing of the fluid in the transverse direction. To obtain a turbulent flow in the gap seal while maintaining the differential pressure and length of the gap seal, it is necessary to increase the gap in the gap seal. In accordance with previous calculations, the gap should be more than 80.5 microns. We take 100 µm for accuracy, because if the gap is more than 100 µm, then the operation of the gap seal becomes ineffective and, accordingly, the study is not of scientific interest.
The velocity distribution when using different flow models is close to the laminar one, then the laminar model was considered together with turbulent models. Figure 7 shows cross-sectional velocity plots for the models under consideration. 2. The lowest maximum speed is observed when using the laminar model. It is about 5 m/s, which is more than two times less than the maximum speed obtained by the RSM model. From the presented data, we can draw the following conclusions: 1.
The maximum values of velocities are almost doubled and amount to (12 ÷ 13) m/s. The highest maximum velocity is observed in the RSM model. The differences between the maximum values of the fluid flow rate, determined by the RSM and SST models, are about 1 m/s, while they practically coincided in laminar flow.

2.
The lowest maximum speed is observed when using the laminar model. It is about 5 m/s, which is more than two times less than the maximum speed obtained by the RSM model. Table 3 presents the average velocities and Reynolds numbers for different models of fluid flow determined by the results of experimental studies [2]. Analyzing the presented results, we have: 1.
The average velocities calculated from the considered flow models are significantly less than the average velocity determined experimentally. The average velocity determined by the SST model is the closest value. However, the discrepancy is 1.83 m/s, which is 15.67%.

2.
The lowest average velocity of 4.66 m/s is shown by the laminar model, which is less than two times less than the velocity determined experimentally.

3.
It should be noted that the Re number determined by the average velocity for all turbulence models is less than 1000. Only the Re number determined by the average velocity determined on the basis of experimental data is more critical than the number Re cr (Re cr = 1000) and is 1160. This suggests that we are generally in the transition area between the laminar and turbulent flow modes.
The volumetric flow rate of the fluid (see Table 4), determined by the SST turbulence model, as well as the average velocity, is closest to the flow rate determined from experimental studies. Summing up the study of turbulent flow in a smooth gap, we can conclude that the results obtained using the SST and RSM turbulence models are closest to the experimental results and, therefore, their application in this case is preferable.

Laminar Flow
There is considered the flow of a viscous fluid in a step-type gap seal with a change in the gap from 25 to 50 µm in the forward and reverse directions. The step is organized from the left border of the gap seal at a distance cm = yn /2.0 = 20 MM, i.e. half its length. Reducing the gap by half the length of the gap seal leads to a decrease in the overall average gap in the gap seal, which leads to a decrease in the velocity and flow rate of the fluid in the gap seal. During the analysis, we will consider the diagrams of instantaneous velocities and the average fluid velocity in two sections: with a gap of 25 µm and with a gap of 50 µm. Figure 8 shows the velocity plots in the cross section with a gap of 50 µm obtained using various flow models.
The maximum velocity in a gap seal of about 1.5 m/s is achieved using the SST and laminar turbulence models. The maximum velocities in the remaining flow models are approximately equal and less by about 0.2 m/s. Figure 8 shows the velocity plots in the cross section with a gap of 50 μm obtained using various flow models.
The maximum velocity in a gap seal of about 1.5 m/s is achieved using the SST and laminar turbulence models. The maximum velocities in the remaining flow models are approximately equal and less by about 0.2 m/s.  Table 5 presents the mean velocities and Reynolds numbers for fluid flow in a gap seal. The results presented in the Table prove that the minimum discrepancy in average velocity and, accordingly, in flow rate between theoretically calculated by formula [2]       Table 5 presents the mean velocities and Reynolds numbers for fluid flow in a gap seal. The results presented in the Table 5 prove that the minimum discrepancy in average velocity and, accordingly, in flow rate between theoretically calculated by formula [2] υ = ∆p 2µl n δ 2 4 − y 2 and model, provides the RSM turbulence model, then k − ε model and the SST model. The greatest discrepancy is observed when using the laminar model. Table 5. Values of average velocities and Reynolds numbers for fluid flow in a step-type gap seals in the forward direction determined by different flow models.

Turbulence Models
The Diagrams of the velocity distribution in the same section are presented in Figure 9, but when the fluid flows in the opposite direction for different turbulence models. When comparing the diagrams of the velocities of the forward and reverse flows, we can conclude that these diagrams qualitatively and quantitatively coincide.

12
  μ μ Diagrams of the velocity distribution in the same section are presented in Figure 9, but when the fluid flows in the opposite direction for different turbulence models. When comparing the diagrams of the velocities of the forward and reverse flows, we can conclude that these diagrams qualitatively and quantitatively coincide.  In general, it is known from the theory of hydraulic resistances that hydraulic losses during sudden expansion exceed hydraulic losses during sudden narrowing [4,5].
A significant parasitic vortex is formed with sudden expansion, which prevents the flow of the main stream. This phenomenon is well observed when using all flow models (for example, see Figures 10 and 11). With a sudden narrowing, a parasitic vortex is also formed, but its size is much smaller and the localization is different. It is localized at a wall with a large gap near a sudden narrowing (see Figures 12 and 13). Such a distribution of velocity vectors is characteristic of all flow models. In general, it is known from the theory of hydraulic resistances that hydraulic losses during sudden expansion exceed hydraulic losses during sudden narrowing [4,5].
A significant parasitic vortex is formed with sudden expansion, which prevents the flow of the main stream. This phenomenon is well observed when using all flow models (for example, see Figures 10 and 11). With a sudden narrowing, a parasitic vortex is also formed, but its size is much smaller and the localization is different. It is localized at a wall with a large gap near a sudden narrowing (see Figures 12 and 13). Such a distribution of velocity vectors is characteristic of all flow models.         Table 6 presents the results of average velocities and Reynolds numbers for different models of fluid flow and the known theoretical dependence. The results presented in Table 6 on qualitative and quantitative characteristics practically correspond to the results presented in Table 5. Table 6. Values of average velocities and Reynolds numbers during fluid flow in a step-type gap seals in the opposite direction determined by different flow models.  Table 6 presents the results of average velocities and Reynolds numbers for different models of fluid flow and the known theoretical dependence. The results presented in Table 6 on qualitative and quantitative characteristics practically correspond to the results presented in Table 5.
Comparing the results presented in Tables 5 and 6, we have the following conclusions: 1.
Results on average velocity and fluid flow in a gap seal obtained by model k − ε are closest to the results obtained from known theoretical dependencies [2,4,5].

2.
The RSM model is the only turbulence model that has lower velocity and lower flow rate during sudden expansion (forward flow) compared to sudden compression (reverse flow).

3.
Taking into account that the RSM model describes correctly qualitatively and quantitatively the fluid flow in a step-type gap seals, the RSM model is to be used when calculating a step seal in laminar mode.

Turbulent Flow
The flow of a viscous fluid in a step-type gap seal is considered when the gap changes from 50 to 100 µm. In this case, a fluid flow is realized with a sudden expansion. Analyzing the patterns of the distribution of modules and velocity vectors, we can conclude that when using the laminar, SST, and RSM models, two areas with a stagnant-flow areas can be seen (see Figures 14 and 15) and the curvature of the flow core.         For further analysis, we choose a section with a gap size of 100 μm. Figure 18 shows plots of the velocity distribution in the cross section of the gap seal. For further analysis, we choose a section with a gap size of 100 µm. Figure 18 shows plots of the velocity distribution in the cross section of the gap seal.  Table 7), we have the following conclusions:

X, mm
1. The smallest discrepancy between the results of experimental studies and flow models is observed when using the RSM model and it is 0.15 m/s, which is 4.3%.
2. Using the SST and laminar models, it is possible to determine the average fluid flow rate with approximately the same error of 0.39 m/s, i.e., 11.33%. The highest maximum speed (a little more than 6 m/s) is observed when using the laminar model, and the lowest maximum speed (about 3.5 m/s) is observed when using the k − ε turbulence model. This can also be seen in Figure 16. The RSM and SST flow models show approximately the same average values of maximum velocity and it is about 5 m/s.
Analyzing the values of average velocities in a gap seal obtained using different models of viscous fluid flow (see Table 7), we have the following conclusions:

1.
The smallest discrepancy between the results of experimental studies and flow models is observed when using the RSM model and it is 0.15 m/s, which is 4.3%.

2.
Using the SST and laminar models, it is possible to determine the average fluid flow rate with approximately the same error of 0.39 m/s, i.e., 11.33%. 3.
k − ε flow model has the worst results. The discrepancy in determining the average velocity is 0.64 m/s, which is 18.6%. In the reverse fluid flow, the same physical picture is observed in all models: there is separation of the flow and two circulation zones. The first significant circulation area is observed before the step-type change in the gap, and the second small circulation zone and flow separation are observed in the narrow part of the gap, immediately after the cross section (see Figures 19 and 20 In the reverse fluid flow, the same physical picture is observed in all models: there is separation of the flow and two circulation zones. The first significant circulation area is observed before the steptype change in the gap, and the second small circulation zone and flow separation are observed in the narrow part of the gap, immediately after the cross section (see Figures 19 and 20).  The velocity plots obtained in the section with a gap of 100 μm for different flow models are presented in Figure 21. The velocity plots obtained in the section with a gap of 100 µm for different flow models are presented in Figure 21. The velocity plots obtained in the section with a gap of 100 μm for different flow models are presented in Figure 21. For the reverse flow diagrams, the velocity distributions differ primarily, qualitatively and slightly quantitatively from the diagrams in the forward flow: For the reverse flow diagrams, the velocity distributions differ primarily, qualitatively and slightly quantitatively from the diagrams in the forward flow:

1.
The laminar flow model has the highest maximum speed, as in the case of a direct flow, but this value is less than in a direct flow.

2.
The SST flow model increased its maximum velocity compared to the direct flow, and the velocity distribution diagram almost coincided with the laminar velocity distribution diagram.

3.
Velocity distribution plots for k − ε and RSM turbulence models have a qualitative and quantitative coincidence in the forward and reverse flows.
Analyzing the results of calculations at average velocities in the reverse flow for different models (see Table 8), we have the following conclusions: The RSM turbulence model gives the closest result to the experimental values of flow rate and average velocity (see Tables 4 and 8). This is observed both in the forward and reverse directions.
It should be noted that in the opposite direction, this difference increases and it is almost 6%. 2.
The average velocities and flows obtained using the laminar and SST flow models are approximately the same in the opposite direction of the flow, as well as in the forward direction and they are equal to 3.77 m/s, the difference between the experimental value is 0.33 m/s or 9.6%. 3.
The largest discrepancy between the experimental results and the flow model in average velocity and flow in the forward and reverse directions is observed when using the k − ε model. Comparing the values of average velocities and flows in the forward and reverse directions, we can conclude that the RSM and k − ε models increase average velocities and costs during sudden compression, while the laminar and SST models reduce it. The increase in average velocities and costs correspond to the general theory of local resistances.
Summing up the above on the turbulent flow in a step-type compaction, we can conclude that the RSM turbulence model describes it most accurately qualitatively and quantitatively.

1.
As a result of the analysis of existing gap seals, two seals were selected: a simple annular gap with a constant gap and an annular gap with a stepwise change in the gap. These types of gap seals are structurally simple, as well as highly efficient from the point of view of ensuring a high mass flow ratio in the forward (from the pump section to the compressor) and reverse directions (from the compressor section to the pump).

2.
The analysis of the existing pressure drops acting on gap seals has established that the pressure drop is within 1.0 MPa, and the slot in the gap reaches 100 µm. With these parameters, both laminar and turbulent fluid flow modes can be implemented in hybrid power machines.

3.
The calculation of the flow of laminar and turbulent fluid flow in the gap seal in the forward and reverse directions. An analysis of existing turbulence models is carried out for a turbulent flow; the main models are selected that are widely used in the calculation of fluid flows, including seals and the features of their implementation in the ANSIS CFX package considered. The article provides calculations for the main models: laminar, k − ε, SST and RSM. A comparison of the velocity distribution diagrams in the gap for laminar and turbulent flows in the forward and reverse directions, average fluid velocities, and flow rates in a gap seal is carried out. The results of experimental studies are also compared. The results of the analysis on using various flow models for calculating gap seals of hybrid power machines of volumetric action are presented in Table 9. Each column presents two models showing the best result in the velocity distribution in the gap, average speed, flow rate and adequacy to known physical concepts.

4.
To calculate the piston seals of hybrid power machines proceeding from the results proposed in Table 9, RSM turbulence models should be used, as the most accurately determining the speed, flow rate and physical flow pattern.

5.
Further development of the presented studies is possible by considering the unsteady fluid flow in the gap seal and taking into account the inertia forces on the flow characteristics. This is very important because the pressure in the pump and compressor cavities vary significantly in the angle of rotation of the shaft. Table 9. The results of the analysis on using various flow models for calculating gap seals of piston hybrid power machines.

Smooth Gap Seal
Step-Type Gap Seal