Next Article in Journal
Path Planning Strategy for a Manipulator Based on a Heuristically Constructed Network
Next Article in Special Issue
Kinematic Characteristics Analysis of Orbiting Scroll and Structural Optimization of Oldham’s Coupling in Scroll Compression
Previous Article in Journal
A Calculation Method for Tooth Wear Depth Based on the Finite Element Method That Considers the Dynamic Mesh Force
Previous Article in Special Issue
A Review of Recent Research and Application Progress in Screw Machines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of the Influence of Thermophysical Properties of a Droplet Liquid on the Work Processes of a Two-Stage Piston Hybrid Power Machine

1
Hydromechanics and Machines Department, Omsk State Technical University, 644050 Omsk, Russia
2
Foreign Languages Department, Omsk State Technical University, 644050 Omsk, Russia
*
Author to whom correspondence should be addressed.
Machines 2022, 10(2), 70; https://doi.org/10.3390/machines10020070
Submission received: 5 December 2021 / Revised: 16 January 2022 / Accepted: 17 January 2022 / Published: 19 January 2022
(This article belongs to the Special Issue Advances in Positive Displacement Compressors)

Abstract

:
The article considers the influence of the main thermophysical properties of the power liquid on the working processes and energy characteristics of a new highly efficient two-stage piston hybrid power machine designed to compress gas to medium and high pressures. Mathematical model of the working processes of the machine under study with a pro-filed working chamber of stage 2 has been developed. To verify the developed mathematical model of working processes, a prototype and a stand for its study were created. The visualization of the movement of the liquid piston, measurement of the instant pressure in the working chamber of the machine, energy and consumption characteristics verified the developed mathematical model. As a result of a numerical experiment, we determined that there is an optimal value of the dynamic viscosity of the power liquid, which depends on the angular velocity of the crankshaft and is in the range from 0.005 to 0.02 Pa∙s, with large values referring to large values of the angular velocity. With a decrease in the density of the power liquid, the work of friction forces decreases and the total indicator efficiency increases. The findings of the research can be used when choosing coolants in reciprocating compressors with a liquid piston.

1. Introduction

The ever-growing demand for new technology prototypes stimulates the development of research and development works, mostly in energy-intensive branches of technology. These branches include pump and compressor construction. More than 20% of the power used is spent on the drive of pumping units [1,2]. As a result, the task of increasing the efficiency of recently produced pumping and compressor units is important and urgent.
An increase in the total efficiency of a piston compressor is connected with an increase in indicator efficiency, mechanical efficiency, and transmission efficiency. Indicator efficiency can be increased in the following ways: by developing the cooling of the compressed gas [2], a decrease in the leakage of the compressed gas through the leaks of the cylinder-piston group and valves [3,4,5], a decrease in losses of work in suction and discharge by reducing the gas and valve closure oscillations [3,4,5]. An increase in the displacement ratio of a reciprocating compressor is due to a decrease in the leakage of the compressed gas, a decrease in dead space, and an increase in the density of the intake gas due to an increase in pressure (resonant or mechanical boost), as well as a decrease in the temperature of the intake gas [6]. It should be noted that the most in the indicator efficiency can be achieved by approaching compression to isothermal by intensive cooling of the compressed gas [7,8]. Currently, there are three main types of cooling used in compressors: air, water and coolant injection.
From a thermodynamic point of view, the coolant injection is the most efficient due to the developed heat exchange surface and absence of thermal resistance between the compressed gas and the coolant [9,10]. However, significant energy costs connected with liquid spraying and its separation, as well as operational difficulties, significantly slow down the widespread introduction of this type of cooling.
Air cooling is simple for practical implementation, however, due to its low efficiency, mostly it is used in general-purpose compressors with a discharge pressure of up to 1.0 MPa of low capacity.
The majority of medium and high pressure piston compressors of small, medium, and high capacity are water-cooled [11].
As a result, the idea of combining a piston compressor and a piston pump into a single unit called the “piston hybrid power machine” (PHPM) appeared. The schematic diagram of the cross-head PHPM is shown in Figure 1.
The PHPM upper part is the compressor section, and the PHPM lower part is the pumping section. One piston is used to compress and move gas and liquid; there is liquid in the gap between the piston and the cylinder, which eliminates the leakage of the compressed gas. In the first case, the pumping section can be independent and be considered a pump for pumping liquid. In this case, the compressor is cooled along the way without energy costs. In the second case, the pumping section can be considered a pump for pumping liquid through the cooling jacket. Then, the cost of driving the pump section must be added to the cost of energy to drive the compressor.
When combining the pump and compressor into a single unit, we achieve intensive cooling of the compressed gas, elimination of leaks in the cylinder–piston group, significant reduction or elimination of dead space, and reduction of friction forces in the cylinder–piston group due to the regulation of liquid friction.
All the aspects result in improved performance and overall compressor efficiency.
Combining with a compressor is useful for a pump. When combined, the suction cavitation margin increased due to the flow of liquid from the gap in the cylinder–piston group; the work supplied in the cycle decreases due to the utilization of the compressed gas heat [12].
Multistage reciprocating compressors are used to compress gas to medium and high pressures. The number of stages varies from 3 to 4, depending on the type of cooling. When using air cooling, the degree of pressure increase in the stage is within 5–6, and when using water cooling, it can be increased to 10–12. When using coolant injection, its value can be raised up to 20, even for gases with a high adiabatic index, such as helium.
One of the promising schemes for multistage gas compression to high pressures is its preliminary compression in a reciprocating compressor and final compression in the cylinder using a liquid piston [13,14].
The use of a liquid piston regulates its intensive cooling, eliminates leaks in the cylinder-piston group, and the size of the dead space. These factors lead to an increase in the indicator efficiency and the compressor delivery ratio.
Thus, this gas compression scheme needs a compressor for preliminary gas compression and a pump for its final compression. This condition is satisfied by the two-stage PHPM shown in Figure 2. The PHPM stage 1 is a single-stage crosshead power machine with piston 1 located in cylinder 2. Piston 1 divides the working volume of cylinder 2 into compressor chamber 5 and pump chamber 13. The pumping chamber is connected by pipeline 14 with stage 2 of the machine, which has liquid volume 16 (liquid piston) and gas volume 17.
The compressor section of stage 1 has suction chamber 6 and discharge chamber 7 connected through suction valve 3 and discharge valve 4 to working chamber 5. The discharge chamber 7 of stage 1 is connected by pipeline 8 to suction chamber 9 of stage 2. Suction chamber 9 is connected by suction valve 11 of stage 2 of the compressor section to working chamber 17, and the latter is connected to the pressure chamber of stage 2 by pressure valve 12. Compressed gas is supplied from discharge chamber 10. There are heat exchangers 15 and 18 to remove heat from the preliminary compressed gas and the power liquid.
The two-stage PHPM works as follows: During the upward stroke of piston 1, gas is compressed and pumped from working chamber 5 through discharge valve 4 into pumping chamber 7. At the same time, pumping chamber 13 increases and the liquid from chamber 16 of stage 2 enters chamber 13 through connecting pipeline 14. The preliminary compressed gas enters suction chamber 9 through connecting pipeline 8, and through suction valve 11 enters gas chamber 17 of stage 2, where the gas pressure decreases.
During the downward stroke of piston 1, gas is sucked from suction chamber 6 through suction valve 3 into working chamber 5. The volume of pumping chamber 13 decreases, thus the liquid enters chamber 16 through pipeline 14, its volume increases, and volume of gas chamber 17 decreases.
A decrease in chamber 17 leads to an increase in its pressure. When the discharge pressure is reached, discharge valve 12 opens and gas enters discharge chamber 10.
When liquid in the pumping section of stage 1 exceeds the volume of the working chamber of stage 2, the excess liquid enters the gas cap installed on the connecting pipeline. In the considered machine, there is no gas cap and the volume of stage 1 pumping section is equal to the working volume in stage 2 compressor section.
Thus, the two-stage PHPM provides intensive cooling in stages 1 and 2, elimination of leaks in their cylinder-piston groups, elimination of dead space in the stages of the machine, which provides a high indicator efficiency and compressor delivery coefficient at high weight and dimensional characteristics.
Completed works on the study of reciprocating compressors with a liquid piston [15,16,17] were aimed at studying the processes of heat transfer and its effect on the compression efficiency. In [18], we analyzed the effect of inserts from a porous medium on the operation of a liquid reciprocating compressor. The influence of the speed of the liquid piston on the working processes of the compressor was studied in [19]. The authors of [20] considered the influence of the main size and operational parameters on the dynamics of liquid movement and work processes in a reciprocating compressor with a liquid piston.
The references review shows that mostly the influence of size, operational, and thermal parameters on the operation of a liquid piston was studied. However, studies on the influence of the main thermophysical properties of a liquid on the dynamics of the movement of a liquid piston could not be found.
A significant influence on the dynamics of the motion of a droplet liquid, which ensures gas compression in stage 2 of the compressor, is exerted by its main thermophysical properties: density, viscosity, coefficient of volumetric compression, and coefficient of thermal expansion. Taking into account that the effective pressure drop on the liquid is not more than 10 MPa, and the modulus of elasticity of the droplet liquids is about 2000 MPa, the effect of the compressibility of the liquid on the working processes in the machine can be neglected. The performed theoretical studies established that, in the range of the investigated pressures (up to 15 MPa), the compressibility of a liquid in a two-stage PHPM under the action of gas forces, inertial forces, and resistance forces does not exceed 1%, and the decrease in the instant velocity of the liquid piston is within 5%. Thus, we assume that the liquid is incompressible.
It should be noted that the temperature coefficient βT for dripping liquids is higher than the volumetric compression ratio βp, but its absolute value remains small. So, for water, it is 0.63 × 10−4 K−1. It should be noted that low absolute βp and βT even if they are changed several times will not have a significant effect on the working processes of the machine.
The density ρw and the dynamic viscosity of the liquid μw are different. Values ρw and μw change significantly both when changing thermodynamic parameters, primarily temperature, and when using various liquids. So, for example, the dynamic viscosity of water at T = 293K is μw = 1004 × 10−3 Pa·s, but the dynamic viscosity of industrial oil M-20A at T = 313K is μw = 26.7 × 10−3 Pa·s, i.e., 26 times greater. We observe similar trends with the density of the liquid.
This research is devoted to the analysis of the influence of the density and viscosity of the working fluid on the working processes and integral characteristics (indicative power and pressure loss in connecting communications) of a two-stage PHPM. The issue of designing and manufacturing a prototype of a piston hybrid power machine with a liquid piston, designed to compress gas to medium and high pressures, is considered; the physical picture of the working processes in the machine under study is studied. With the developed mathematical model of the working processes of the prototype, the influence of the main thermophysical properties of the working fluid: dynamic viscosity and density on the working processes and energy characteristics of the new two-stage PHPM will be considered.

2. Theoretical Research

2.1. Determination of Basic Geometric Parameters

For a given capacity and number of revolutions, using well-known techniques [11], the diameter of stage 1 piston, stage 1 piston stroke, the main parameters of the self-acting valves, and the diameters of the interstage pipelines are determined. In our case, we coordinate the operation of stage 1 of the pump section and stage 2 of the compressor section. In this case, we determine the volume of the working chamber of stage 2 of the compressor section and the diameter of the stem in the pump section. In addition, to improve the operation of stage 2 of the PHPM, profiling the working chamber has been completed.

2.1.1. Determination of Working Chamber of Stage 2

We determine the amount of intake gas into stage 2 of the compressor, neglecting leaks and overflows of gas in stage 1 of the compressor, from the equality of the masses of gas entering stages 1 and 2. Assuming that the compressible gas follows the laws of an ideal gas, we obtain:
V s c 2 = V s c 1 T s c 2 T s c 1 P s c 1 P s c 2
Considering that V s c 1 = V h 1 / λ 1 , but V s c 2 = V h 2 / λ 2 Equation (1) is transformed to:
V h 2 = V h 1 λ 2 λ 1 T s c 2 T s c 1 P s c 1 P s c 2
Values T s c 2 and P s c 2 can be determined as T s c 2 = T s c 1 + Δ T H 0 , P s c 2 = P D 1 Δ P c .
In view of the above, Equation (2) is transformed to:
V h 2 = V h 1 λ 2 λ 1 ( 1 + Δ T H O T s c 1 ) ( 1 ε 1 Δ P c P s c 1 )  
Values λ 1 ,   λ 2 and Δ T H 0 , Δ P c are determined on well-known recommendations [11].
The working volume of the pumping section of stage 1 V c p 1 equals to the working volume of stage 2. Thus, the rod diameter is determined as:
d r = d 1 2 4 V h 2 π S h 1

2.1.2. Profiling of Stage 2 Working Chamber

When designing and operating compressor units with a liquid piston, a number of measures are used to improve the kinematics and integrity of the liquid piston by reducing foaming, crushing, and carrying part of the liquid through the discharge valve out.
For this purpose, the compressor chamber is profiled (the working chamber of stage 2 of the PHPM compressor section) and devices for stabilizing the movement of liquid of various designs are used, for example, made in the form of movable and stationary permeable solid inserts and elements.
The main reason for the liquid piston destruction is the effect on the liquid of inertial forces caused by the acceleration of the pump piston supplying liquid to the working chamber. In order to reduce inertial forces, it is necessary to bring the liquid velocity in the working chamber closer to a certain constant value, for example, the average piston speed in the pump.
To achieve this goal, it is necessary to profile the working chamber, i.e., change the flow area (F) along the length of the working chamber (S).
Currently, various profiling techniques have been developed based on different optimization criteria or their combination [21]. In practice, it is advisable to use working chambers obtained using simple bodies of revolution: a cylinder and a cone (see Figure 3).
In [21], after analyzing the working processes and operating features of stage 2 of the high-pressure compressor, it was proposed to replace the cylindrical chamber (position “a” Figure 3) with a combined one (position “b” Figure 3) having two conical parts and one cylindrical. Conducted theoretical studies [21] established that the working chamber, consisting of two truncated cones, provides a high heat exchange surface, good stabilization of the liquid piston speed with high ergonomic parameters.
Considering the above, we use the working chamber of stage 2 consisting of two truncated cones (see Figure 3, position “c”), its working volume can be determined as:
V h 2 = 2 S 1 1 3 π [ R 1 2 + R 1 R 2 + R 2 2 ]
The ratio of the height of the cone to the diameter of its base, in accordance with [15] is to be chosen in the range from 10 to 20.

2.2. Review of Mathematical Models and System of Accepted Assumptions

The two-stage piston hybrid power machine has stages 1 and 2 of compression and inter-stage connections for gas and liquid. Stage 1 of the investigated machine is a single-stage piston hybrid power machine. In [22], the mathematical model of working processes of a cross-head PHPM with a smooth seal is considered for the first time. The works [13,23] are devoted to the development of a mathematical model of the working processes of a cross-head PHPM with a gas chamber in the pump section. In these works, the simultaneous compression of liquid and gas in the working processes of a piston pump was considered for the first time.
The developed mathematical models are based on the following basic equations:
  • for the chambers of the compressor section: the first law of thermodynamics of a body of variable mass, the equation of conservation of mass, the equation of volume change, the equation of state, the equation of the dynamics of motion of the closure body of the self-acting valve;
  • for the working chamber of the pump section: the energy conservation equation written in the form of the Bernoulli equation; Hooke’s equation, volume conservation equation, mass conservation equation. In [24,25], a general approach was developed for calculating work processes in the working chamber of a compressor and a positive displacement pump.
It should be noted that the developed models simulate the working processes in the compressor section, pump section, piston seal, as well as the flow of gas and liquid in the connecting pipelines.
When calculating the gas flow in the connecting pipelines of the compressors, we use a system of equations including the continuity equation, the equation of motion, and the equation for the conservation of internal and kinetic energy. The equation of state for an ideal gas is used as the equation of state. One of the effective methods for the numerical solution of such a system of partial differential equations as for single-phase flow or two-phase flow is the “large particles” method [7].
Currently, when modeling work processes in machine chambers, both distributed-parameter and lumped-parameter models are used. Taking into account the small linear dimensions of the working chambers of the displacement machines [11], as well as significant mathematical difficulties and implementation time, the researchers use models with lumped parameters. As a result, the use of mathematical models with lumped parameters in the above works is verified and useful. It seems reasonable to use a mathematical model with lumped parameters when developing a mathematical model of a two-stage PHPM.
When considering the flow of gas and liquid in the connecting pipes of the machine, a one-dimensional flow is considered either in an unsteady or in a quasi-stationary setting [26]. If the pipeline is simple, without local resistance, and a mathematical model with distributed parameters is used to describe the working processes, we use a non-stationary one-dimensional flow model. If there are complex pipelines with closure and control valves, as well as when using a mathematical model with lumped parameters to describe thermogasdynamic processes in machine chambers, we use a quasi-stationary model for calculating the flow of gas and liquid in connecting pipes.
System of accepted basic assumptions:
  • taking into account the above, we accept a mathematical model with lumped parameters to simulate the working processes of a two-stage PHPM;
  • taking into account that the mathematical models of the working processes of single-stage cross-head and crosshead-free PHPMs are well developed, we consider the working processes in stage 2 of the PHPM and the liquid movement from the pump section of stage 1 through the connecting pipeline to stage 2;
  • the thermodynamic gas parameters in stage 2 suction chamber are constant and equal to their rated values;
  • the motion of the liquid at each moment of time is steady and stationary, i.e., the motion of the liquid is quasi-stationary;
  • the working liquid has constant temperature and follows Newton’s laws;
  • the thermophysical properties of the liquid are constant;
  • coefficients of friction along the length and local coefficients of resistance obtained for stationary conditions are also valid for unsteady flow;
  • the bottom of the liquid piston is a plane parallel to the ground;
  • there is no droplet liquid in the compressed gas in stage 2 of the PHPM;
  • compressed gas follows the laws that are true for an ideal gas;
  • the simulated processes in a compressible gas are equilibrium and reversible;
  • the change in kinetic and potential energy for a compressible gas can be neglected.
  • the movement of the closure body of a self-acting valve is considered as the movement of a material point;
  • the walls of the pipeline and working chambers are considered absolutely rigid.

2.3. Mathematical Model of the Working Processes of a Two-Stage PHPM

The working chamber of stage 2 consists of two truncated cones (see Figure 4), having the same geometric dimensions and located symmetrically about the axis passing through the center of the working chamber and perpendicular to the vector of the liquid velocity. Thus, we have a working chamber consisting of a diffuser and confuser parts, when the liquid moves in both forward and reverse directions.
The system of differential equations describing the change in the thermodynamic parameters of gas in stage 2 of the PHPM and the motion of liquid includes: the first law of thermodynamics of a variable mass body for an open thermodynamic system; the equation of conservation of mass, the dynamics of motion of the closure body of the self-acting valve, the equation of state for a compressible gas, as well as the equation of conservation of energy and the equation of continuity for the liquid phase, can be written as:
d U g = d Q g p g d V g + i = 1 N 1 i n i d M n i i = 1 N 2 i 0 d M 0 i                          
d M g = i = 1 N 1 d M n i i = 1 N 2 d M 0 i      
m r d v d 2 h c d τ 2 = F i          
d p g = ( k 1 ) ( d U g V g U g V g 2 d V g )
d T g = 1 k ( p g d V g + V g d p g ) T g d M g    
d z w + d p w ρ w g + d ( v w 2 2 g ) + d ( h l + h ζ + h i n ) = 0
d Q w = d ( V w f w ) = 0
When determining the elementary amount of heat removed from the compressed gas, the Newton–Richmann hypothesis was used:
d Q g = α w i = 3 N 3 = 3 F n i ( T ¯ w T g ) d τ
The total heat transfer surface during gas compression can be defined as:
i = 1 N 3 = 3 F n i = F p + F b + F v p
Piston surface area Fp and valve plate Fvp can be defined as F p = π ( R 1 + z 21 t g α ) 2 , F v p = π R 1 2 .
The side surface of the working chamber (see Figure 4) can be determined as
F b = π ( R 2 + R 1 ) · S 1 c o s α + π ( R 2 + R 1 + z 21 · t g α ) ( S 1 z 21 ) c o s α
Heat transfer coefficient α w , in general for a reciprocating compressor is determined on the experimental studies. In most practical cases, an empirical ratio is used to determine the heat transfer coefficient in the cylinder of a single-acting compressor N u = f ( A ,   B , R e ) [11]. The piston speed is used as the determining speed in defining the Reynolds number. In our case, we use the average mass velocity of the liquid in the working chamber of stage 2 as the determining speed.
The average integral value of the surface temperature of the working chamber can be determined as:
T ¯ w = T ¯ v p F v p + F b T ¯ b + F p · T ¯ w i = 1 N 3 F n i
When determining the mass flows of the attached (dMni) and detached (dM0i) through the gas distribution bodies, the San Venant–Wanzel equation was used for a one-dimensional isentropic flow in a quasi-stationary setting. Changing working chamber dVg was limited by a change in the volume in the pumping section of stage 1 in the absence of dropping liquid compressibility.
When integrating Bernoulli Equation (11), two sections were used: the first section (I-I) was aligned with the lower boundary of the mechanical piston of stage 1, and the second section (II-II) was aligned with the surface of the liquid piston.
The first coordinate can be determined as:
z 1 = S l d + S h 2 [ ( 1 c o s φ ) + λ m 4 ( 1 c o s 2 φ ) ]
The second coordinate can be determined as:
  • the height of the surface of the liquid piston is less than the height of the truncated cone (S1), then z 2 = z 20 + z 21 . Value z21 is determined from the solution of the following nonlinear volume conservation equation
    1 3 π z 21 [ R 1 2 + R 1 · ( R 1 + z 21 t g α ) + ( R 1 + z 21 t g α ) 2 ] = { S h S h 2 [ ( 1 c o s φ ) + λ m 4 ( 1 c o s 2 φ ) ] } · F 1
  • the surface of the liquid piston is greater than the height of the truncated cone (S1), then z 2 = z 20 + S 1 + z 22 . Value z22 is determined from the equation of conservation of the liquid volume
    1 3 π S 1 ( R 1 2 + R 1 · R 2 + R 2 2 ) + 1 / 3 π z 22 · [ R 2 2 + R 2 ( R 2 z 22 t g α ) + ( R 2 z 22 t g α ) ] = { S h S h / 2 [ ( 1 c o s φ ) + λ m / 4 ( 1 c o s 2 φ ) ] } · F 1  
The liquid velocity in the first section V1 is determined by the kinematics of stage 1 piston drive mechanism, and the liquid velocity in the second section V2 is determined based on the continuity equation. Total head losses consist of head losses along the length Δ h l , head loss on local resistance Δ h ζ and inertial head losses Δ h i n ( Δ h Σ = Δ h l + Δ h ζ + Δ h i n ).
We determine the values included in the total head losses Δ h Σ .
Head losses along the length ( Δ h l ) are determined from the head losses along the length in the pumping chamber of stage 1, head losses along the length in pipeline 3 and head losses along the length in the profiled working chamber of stage 2 Δ h l = Δ h l 1 + Δ h l pp + Δ h l 2 .
Head losses along the length Δ h l 1 and Δ h l p p can be determined based on the Darcy–Weisbach law as follows Δ h l 1 = λ c 1 ( z 1 / d 1 ) ( v 1 2 / 2 g ) , Δ h l pp = λ c p p ( l pp / d pp ) ( v pp 2 / 2 g ) .
Length friction coefficients λ c 1 , λ c p p , λ c 2 in general, are a function of the Reynolds number and relative roughness and is determined on the known recommendations.
The compression chamber of stage 2 consists of a diffuser and a confuser, and the head losses can be determined in them as follows:
  • for the diffusor
    Δ h l 2 = λ c 2 8 sin α ( 1 1 n g 2 ) v 21 2 2 g
    where n g = ( R 1 + z 21 t g α ) 2 / R 1 2 is diffuser expansion ratio; v 21 = v 1 F 1 / ( π R 1 2 ) is liquid velocity in the initial section of the diffuser.
  • for the confusor
    Δ h l 2 = λ c 2 8 sin α ( 1 1 n K 2 ) v 22 2 2 g
    where n k = R 2 2 / ( R 2 z 22 t g α ) 2 is the degree of narrowing of the confuser; v 22 = v 1 F 1 / ( π [ R 2 z 22 t g α ] ) is liquid velocity in confuser.
The head losses on local resistances ( Δ h ζ ) depend on the direction of the liquid flow, as when moving forward the local resistance is a sudden narrowing, and with the reverse liquid movement—there is a sudden expansion.
As a consequence, when determining local resistances, it is necessary to take into account the direction of the piston movement.
All local resistances can be divided into four types:
  • local resistance when connecting stage 1 pumping section to connecting pipeline 3—sudden narrowing ( ζ 1 );
  • local resistance when flow turns ( ζ 2 );
  • local resistance when connecting pipeline 3 to the working chamber of stage 2—sudden expansion ( ζ 3 );
  • constant expansion in the diffuser part of the working chamber ζ g and constant narrowing in the confuser part of the working chamber ζ k .
In sudden narrowing, ζ 1 can be determined as:
ζ 1 = 0.5 ( 1 n ζ )  
where n ζ = ω 2 / ω 1 is the ratio of the areas of flow sections after narrowing and before narrowing. We assume that ζ 1 = 0.5 , as ω 1 .
When turning the pipe, the value ζ 2 is determined as [21]:
ζ 2 = 0.05 + 0.2 d pp R pp  
where R pp is bending radius of the pipeline (see Figure 4).
Local resistances during sudden expansion are determined by the Borda formula. Local resistances in the diffuser and confuser parts are determined according to [26].
Inertial head losses Δ h i n consist of inertial head losses in stage 1 of the pump section in the connecting pipeline and in stage 2 Δhinhin1hinpphin2.
In stage 1 of the pump section and in the connecting pipeline, the determination of the inertial head losses does not depend on the direction of liquid movement. Their value can be defined as Δ h in 1 = a n g z 1 ; Δ h i n pp = a pp g l cpp .
We consider the inertial head losses in stage 2 of the PHPM.
When the liquid piston moves up:
  • At z2i < S1 the inertial head losses are determined as:
    Δ h i n 2 = 1 g R 1 R i a i t g α d R  
    where R i = R 1 + z 21 t g α ; a i = F 1 f i a n F 1 v 1 f i 2 ω d f i d ϕ ; f i = π R i 2 .
  • At z2i > S1
    Δ h i n 2 = 1 g R 1 R 2 a i t g α d R + 1 g R 2 R i a i t g α d R  
    where, R i = R 2 z 22 t g α .
When the piston moves down, Δhin2 is defined similarly.

3. Experimental Studies

The main purpose of the experiment is to check the performance of the two-stage PHPM, to gain new knowledge about its work, and to verify the developed mathematical model.

3.1. Subject of Research, Stand, and Measuring Equipment

A schematic diagram of a two-stage PHPM was previously shown in Figure 4. To exclude the influence of gas pulsations from the first stage, a 25 L receiver was additionally installed between the compressor stages.
The developed prototype of the two-stage PHPM had basic dimensions, Table 1.
The basic geometrical dimensions of the second stage are shown in Figure 5a.
The profiling of the working chamber of the second stage was carried out by installing the insert in the cylinder. The insert consisted of two cones connected at their vertices. Some accesses are made on the upper cone to ensure the passage of gas to the suction and discharge ports, and the base of the lower cone is made conical to improve the flow of liquid around it. The 3D view of the redesigned second stage is shown in Figure 5b. A photograph of the finished layout of the modernized second stage is shown in Figure 5c. Lobe valves were used as gas distribution bodies in the compressor section. To visualize the dynamics of the movement of the liquid piston, the second stage of the PHPM was made of Plexiglas, and the value of the discharge pressure was limited.
Table 2 shows the main dimensions of self-acting valves in stage 2 compressor section.
Working bodies used during experimental research are air, mineral oil (MGE-46V).
The prototype was tested on the developed stand; its general view is shown in Figure 6. A general view of the prototype installed on the stand is shown in Figure 7a. An image of the prototype on the stand is shown in Figure 7b.
To change smoothly the rotational speed of the drive shaft of the machine, a volumetric hydraulic drive is used, which includes a variable displacement axial piston pump model 313.3.56.804 with a discharge pressure 6.3 MPa and an axial piston hydraulic motor model 310.3.56.01.03.V.U.
Measurement of static pressures is performed by manometers, instantaneous pressures are measured with strain gauges.
Suction and discharge temperatures are measured by TW-N PT100 sensors, the temperature of the liquid in the tanks is controlled by DTS-035-50M.B3 sensors.
The air flow rate at the suction of the first stage of the compressor is measured by the “Vector-04” model flow meters; we used the PSE530-M5-l sensors to measure the suction and discharge pressures.
The hydraulic and pneumatic lines of the stand are equipped with devices for liquid and gas purification, safety devices, and control equipment to operate the tested machine in different modes.

3.2. Verification of the Mathematical Model of Work Processes

One of the main tasks of the experimental research of the developed machine is to verify the mathematical model. The mathematical model verification was carried out at different values of the discharge pressure in stage 2, the suction pressure in stage 2 (discharge pressure in stage 1) and the crankshaft revolutions. An estimate of the experimental error and a comparison of the results of theoretical and experimental studies were published in [27].
It should be noted that the developed mathematical model describes the friction and inertia forces in moving liquid, as the data of indicator diagrams obtained experimentally and theoretically match; the maximum and minimum values of instantaneous pressures match as well.
Thus, the discrepancy in determining the instantaneous pressure in the machine cavities is within 10–15%.

4. Results of a Numerical Experiment and Their Discussion

To carry out a numerical experiment, we choose a two-stage piston hybrid power machine with a profiled working chamber shown in Figure 8. The main geometrical and operational parameters are shown in Table 3.
As independent variables, we take the above-stated coefficient of dynamic viscosity and density of the liquid.
The dynamic viscosity coefficient varied from 0.0005 Pa to 0.04 Pa, and the density of the liquid varied from 600 kg/m3 to 1400 kg/m3.
In carrying out the numerical experiment, we used the classical design with fractional replicas.
All response functions can be divided into two large groups. The first group is made up of functions that determine the flow of liquid in the machine: Δ p ¯ 1 = Δ p 1 / p Σ ; Δ p ¯ 2 = Δ p 2 / p Σ ; Δ p ¯ 3 = p 3 / p Σ , where Δ p Σ = Δ p 1 + Δ p 2 + Δ p 3 ;   Δ p 1 = 0 2 π | z 2 z 1 | ρ w g d φ ; Δ p 2 = 0 2 π ρ w | ( α 2 V 2 2 ) / 2 ( α 1 V 1 2 ) / 2 | d φ ; Δ p 3 = 0 2 π ρ w g | Δ h l + Δ h ζ + Δ h i n | d φ .
Δ p ¯ 1 determines the ratio of the geometric head to the sum of the geometric head, velocity head, and total head losses per cycle. Δ p ¯ 2 determines the ratio of the velocity head to the sum of the geometric head, velocity head, and total head losses per cycle. Δ p ¯ 3 determines the ratio of total head losses to the sum of geometric head, velocity head, and total head losses per cycle. Δ h ¯ l = ( 0 2 π | Δ h l | d φ ) / Δ h Σ –the ratio of head losses to hydraulic resistance along the length to the total head losses per cycle; Δ h ¯ ζ = ( 0 2 π | Δ h ζ | d φ ) / h Σ –the ratio of head losses to local resistances to total head losses per cycle; Δ h ¯ i n Σ = ( 0 2 π | Δ h i n Σ | d φ ) / Δ h Σ –ratio of inertial losses, taken in modulus, to total head losses per cycle, where Δ h Σ = 0 2 π ( | Δ h l | + | Δ h ζ | + | Δ h i n Σ | ) d φ .
The second group is made of functions determining the efficiency of stage 2 of the PHPM.
Indicator work of the cycle in stage 2 is L i n d = V 2 d p .
Technical work to overcome the frictional forces of the liquid in the machine L f r = ρ w g ( | Δ h l | + | Δ h ζ | ) d V .
The ratio of the work of friction forces to indicator work L ¯ f r = L f r / L i n d .
Indicator isothermal efficiency of stage 2 of PHPM is:
η ind . is = G c p R T s c ln ( p d / p s c ) L i n d
Full indicator efficiency of stage 2 of PHPM:
η ind . f = G c p R T s c ln ( p d / p s c ) L i n d + L f r

4.1. Analysis of the Influence of the Coefficient of Dynamic Viscosity

It should be noted that when analyzing the effect of liquid viscosity on the working processes of a two-stage PHPM, the total speed of liquid movement in the machine, including the pump section of the right-hand stage, the connecting pipeline and stage 2, as well as the liquid velocity only in the connecting pipeline, is of great importance, since it is the main source of hydraulic losses.
In the first case, the change in the total speed will be achieved by changing the angular speed of rotation, and in the second case—by changing the diameter of the connecting pipeline.
With a decrease in the dynamic viscosity coefficient μw at a constant density of the medium, the Reynolds number increases, the flow turbulization increases, leading to an increase in the hydraulic friction coefficient and an increase in the relative hydraulic losses along the length Δ h ¯ l and the value of the relative pressure losses per cycle Δ p ¯ 3 (see Figure 9 and Figure 10).
It should be noted that a sharp increase in Δ h ¯ l and Δ p ¯ 3 is observed from the value of the coefficient of dynamic viscosity equal to 0.0008 Pa∙s. Up to this value, we observe a linear drop in these values in the range of variation of μw from 0.02 Pa∙s to 0.0008 Pa∙s, which is fully consistent with the theory of laminar liquid motion. Against the sharp increase of Δ h ¯ l and Δ p ¯ 3 , there is a sharp decrease in Δ p ¯ 1 and Δ h ¯ i n Σ with a smooth decrease in Δ p ¯ 2 .
In the range of μw variation from 0.02 Pa∙s to 0.005 Pa∙s, a decrease in the work of friction forces L f r is observed and with a further decrease in μw the value L f r increases (with a small local maximum) to 6.17 J, which is 14.89% (see Figure 11) from the indicator work of compression: the maximum value of the total indicator efficiency is 73.2% at μw = 0.005 Pa∙s and its value drops to 64.2% at μw = 0.0005 Pa∙s, i.e., by almost 10%. Values L ind and η ind . is remain constant when changing μw.
Thus, the minimum L ind and η ind . f does not coincide with the minimum Δ h ¯ l and Δ p ¯ 3 and is shifted towards large values of μw.

4.2. Influence of the Angular Velocity of the Crankshaft

With an increase in the angular velocity of the crankshaft, the speed of liquid movement in the stages of the machine and in the connecting pipeline increases. This leads to an increase in Δ p ¯ 3 and Δ h ¯ l while maintaining the previously described physical picture. The onset of flow turbulization and a sharp increase in the hydraulic friction coefficient along the length is shifted towards large values of μw.
Thus, if at nrev = 300 rpm, minimum Δ h ¯ l is achieved at μw = 0.005 Pa∙s, then at nrev = 500 rpm minimum Δ h ¯ l is achieved at μw = 0.010 Pa∙s, and at nrev = 750 rpm minimum Δ h ¯ l is achieved at μw = 0.02 Pa∙s (see Figure 12).
With an increase in nrev and an increase in Δ h ¯ l there is an increase in L f r and the total indicator efficiency decreases. Maximum η ind . f shifts towards larger values of μw while simultaneously decreasing it in absolute terms. The maximum η ind . f and minimum L f r points coincide with the minimum points Δ h ¯ l . Absolute decrease in the maximum value of the indicator efficiency is 0.041 or 5.93%.
It should be noted that at the number of revolutions nrev = 750 rpm and μw = 0.0005 Pa∙s, the value of the work of the friction forces L f r slightly exceeds the work for gas compression in the second stage of the machine L ind . is and η ind . f drops to 0.362, i.e., decreases by almost half.

4.3. Influence of the Diameter of the Connecting Pipe

It should be noted that the connecting pipeline is of decisive importance in the balance of total hydraulic losses along the length. This led to its consideration.
A change in the diameter of the connecting pipeline leads to a change in the speed in it and, therefore, the physical picture of working processes and changes in energy characteristics will be similar to that considered. As a result, we consider the working processes of the machine under study and its energy characteristics at constant μw and nrev and a variable value of the pipeline diameter (dpp).
With an increase in the diameter of the pipeline (dpp), velocity of the liquid in it decreases, which leads to a decrease in the hydraulic resistance and an increase in Δ p ¯ 3 with a simultaneous increase in Δ p ¯ 1 and Δ p ¯ 2 .
It should be noted that when dpp changes from 0.014 m to 0.022 m, there is a linear drop in Δ p ¯ 3 from 0.9924 to 0.981 (see Figure 13).
With a decrease in dpp, there is a drop in Δ h ¯ l and Δ h ¯ ζ with an increase in the relative inertial losses Δ h ¯ i n Σ . The curve h ¯ l = f ( d p p ) has a parabolic character with a sharp increase with decreasing dpp, but the curve h ¯ ζ = f ( d p p ) is close to linear (see Figure 14).
With an increase in dpp and a decrease in hydraulic losses along the length and on local resistances there is a decrease in the work of friction forces L f r and an increase in the total indicator efficiency. With an increase in dpp, the value L f r / L c o m p r decreases from 8.8% to 1%, when increasing dpp from 0.014 m to 0.022 m, but the value of the total indicator efficiency increases from 0.674 to 0.726, i.e., more than 0.05, which is almost 10%.

4.4. Analysis of the Effect of Liquid Density

To analyze the influence of the density of the coolant on the working processes and characteristics of the machine, we take: dpp = 0.018 m; nrev = 500 rpm; μw = 0.01 Pa∙s.
With a change in ρw values Δ p ¯ 1 and Δ p ¯ 2 have maximum at ρw = 1200 kg/m3, although these maxima are negligible and can be neglected, which allows us to conclude that Δ p ¯ 1 , Δ p ¯ 2 , and Δ p ¯ 3 remain constant when ρw change in the range from 600 kg/m3 to 1400 kg/m3.
In the range of density variation, we observe a minimum of h ¯ l at ρw = 1200 kg/m3 and maximum h ¯ l . It should be noted that Δ h ¯ ζ changes insignificantly, while h ¯ l decreases almost twofold. The influence of the density of the working fluid is very complicated. With a change in density, the Re number changes at a constant μw, which leads to a change in the value of the coefficient of friction along the length, and the values of the hydrostatic and velocity heads change. The maximum relative inertial losses are achieved at ρw =1200 kg/m3 and make up the biggest share of the total inertial losses—92%.
With an increase in the density of the liquid the work of friction forces L f r increases almost linearly from ρw =600 kg/m3 to 1100 kg/m3, with a further increase in ρw the value L f r changes parabolically (see Figure 15). With a decrease in density due to the fact that L f r decreases, there is an increase in the total indicator efficiency. Its value increases from 0.7036 at ρw = 1400 kg/m3 to 0.722 at ρw = 600 kg/m3, i.e., by almost 2% (see Figure 15).

5. Conclusions

In the manuscript, a mathematical model of the working processes of a two-stage piston hybrid power machine designed to compress gas to medium and high pressures has been developed and validated. A numerical experiment was carried out to study the influence of the viscosity and density of the working fluid on the working processes in the machine under study.
As a result of the numerical experiment, we determined:
(1)
There is a value of the optimal value μw, which provides the maximum value of the total indicator efficiency and bypassing the work of the forces of liquid friction. The optimal μw is in the range from 0.005 Pa∙s to 0.02 Pa∙s, with a change in the number of revolutions of the crankshaft from 300 rpm to 750 rpm.
(2)
With an increase in the angular velocity of the crankshaft to obtain maximum efficiency, it is necessary to increase the dynamic viscosity coefficient. The use of water as a coolant is not effective, because as the temperature rises, it decreases and at the operating temperature T = 323 K; its value is μw = 0.0005494 Pa∙s, which corresponds to the minimum indicator efficiency and the maximum work of friction forces. An increase in the diameter of the connecting pipeline leads to a decrease in the work of friction forces and an increase in the indicator efficiency.
(3)
With a decrease in the density of the medium, the total indicator isothermal efficiency increases and the work of friction forces L f r decreases. The minimum relative hydraulic losses along the length and the maximum Δ p ¯ 1 and Δ p ¯ 2 correspond to ρw = 1200 kg/m3.
(4)
The analysis showed that industrial oil I-12A, which has ρw = 880 kg/m3 and μw/ρw = (13÷17)∙10−6 m2/s can be used as power liquid in the machines.

Author Contributions

Conceptualization, V.S. and E.N.; methodology, V.S.; software, E.P.; validation, V.S. and E.P.; formal analysis, V.S.; investigation, V.S., E.P., E.N., and I.B.; writing—original draft preparation, V.S. and I.B.; writing—review and editing, I.B. and E.P.; visualization, E.P. and E.N.; supervision, V.S.; project administration, V.S. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by the Russian Science Foundation Grant No. 22-29-00399, https://rscf.ru/project/22-29-00399/ (Accessed date 16 January 2022).

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

The nomenclature of the paper is shown below:
α angle between the generatrix of the cone and the height;
α 1 coriolis coefficient for Section 1;
α 2 coriolis coefficient for Section 2;
α w heat transfer coefficient of compressible gas;
β p coefficient of volume compressibility;
α w coefficient of temperature expansion;
β p head loss due to inertial resistance in stage 1 of the pump section;
β T head loss due to inertial resistance in the working chamber of stage 2;
Δ h i n 1 head loss due to inertial resistance in the pipeline;
Δ h i n 2 head loss along the length in the cylinder of stage 1 pumping section;
Δ h i n p p head loss along the length in the working chamber of stage 2;
Δ h l 1 head loss along the length in connecting pipeline 3;
Δ h l 2 pressure loss in the interstage connecting of the compressor section;
Δ h l p p undercooling value of gas sucked into stage 2 of the compressor section;
Δpcthe ratio of the discharge pressure to the suction pressure of stage 1;
ΔTH0flow rates in stages 1 of the compressor section;
ε1flow rates in stages 2 of the compressor section;
λ1ratio of piston stroke to twice the length of the connecting rod;
λ2coefficient of friction along the length in the cylinder of stage 1;
λmcoefficient of friction along the length in the working chamber of stage 2;
λ c 1 coefficient of friction along the length in the connecting pipe;
λ c 2 dynamic viscosity;
λ c p p liquid density;
μ w sum of forces acting on the valve closure;
ρ w crankshaft angle;
Σ F i piston acceleration;
φacceleration in the connecting pipeline;
a n stage 1 cylinder diameter;
a p p elementary time interval;
d 1 connecting pipe diameter;
d τ elementary i:th separable mass of gas;
d p p elementary i:th attached mass of gas;
dM0irod diameter in PHPM stage 1;
dMnielementary heat input;
drstage 1 mechanical piston area;
dQgworking chamber lateral surface area;
F 1 heat exchange surface;
F b piston surface area;
F n i valve plate area;
F p cross-sectional area;
F v p gravity acceleration;
f w mass of gas supplied to the discharge line;
ghead loss on local resistance;
G c p degree of valve closure lift;
h ζ inertial head losses;
hchead loss along the length;
h i n specific enthalpy of the separated gas mass;
h l specific enthalpy of the added i:th mass of gas;
i0adiabatic exponent;
iniconnecting pipe length;
kreduced mass of the self-acting valve closure;
l c p p number of sources of separated and attached gas masses;
mrdvpressure in stage 1 of discharge chamber;
N1, N2gas pressure;
p D 1 gas pressure at the suction of stage 1 of the compressor section;
p g gas pressure at the suction of stage 2 of the compressor section;
psc1liquid pressure;
psc2liquid flow;
p w small radius of the truncated cone;
Q w large truncated cone radius;
R 1 truncated cone height;
R 2 piston stroke;
S 1 piston stroke in PHPM stage 1;
S h linear dead space;
Sh1average integral temperature of the lateral surface of the working chamber;
S l d temperature of the compressible gas;
T ¯ b suction gas temperature in stage 1 of the compressor section;
T g suction gas temperature in stage 2 of the compressor section;
Tsc1mean integral temperature of the valve plate surface;
Tsc2mean integral temperature of the surface of the compression chamber;
T ¯ v p total internal energy of the compressible gas;
T ¯ w mechanical piston speed in stage 1 pumping section;
U g gas volume;
V 1 the working volume of the pumping section of stage 1;
V g working volume of stage 1 of the compressor section;
V c p 1 working volume of stage 2 of the compressor section;
V h 1 liquid velocity in the connecting line;
V h 2 suction gas volume in stage 1 of the compressor section;
V p p suction gas volume in stage 2 of the compressor section;
V s c 1 liquid velocity;
V s c 2 coordinate of the gravity center of the first cross-section;
V w coordinate of the gravity center of the second cross-section;
z 1 distance from the axis of connecting pipeline 3 to the lower surface of the working chamber of stage 2;
z 2 coordinate of the surface of the liquid piston;
z 20 distance from the surface of O-O cross-section to the surface of the liquid piston;
z 21 the coordinate of the gravity centers of the cross-sections;

References

  1. Abdel-Hadi, A.; Salem, A.R.; Abbas, A.I.; Qandil, M.; Amano, R.S. Study of energy saving analysis for different industries. J. Energy Resour. Technol. Trans. ASME 2021, 143, 052101. [Google Scholar] [CrossRef]
  2. Silva, E.; Dutra, T. Piston trajectory optimization of a reciprocating compressor. Int. J. Refrig. 2021, 121, 159–167. [Google Scholar] [CrossRef]
  3. Yan, B.; Wieberdink, J.; Shirazi, F.; Li, P.Y.; Simon, T.W.; Van de Ven, J.D. Experimental study of heat transfer enhancement in a liquid piston compressor/expander using porous media inserts. Appl. Energy 2015, 154, 40–50. [Google Scholar] [CrossRef]
  4. Dindorf, R. Estimating Potential Energy Savings in Compressed Air Systems. Procedia Eng. 2012, 39, 1877–7058. [Google Scholar] [CrossRef] [Green Version]
  5. Ueno, K.; Hunter, K.S. Compressor Efficiency Definitions; University of Colorado: Boulder, CO, USA, 2003. [Google Scholar]
  6. Shcherba, V.E.; Paramonov, A.M.; Blinov, V.N.; Surikov, V.I.; Nosov, E.Y.; Tegzhanov, A.S. Comparative analysis of process of cooling of compressible gas in crosshead and crossheadless hybrid positive-displacement piston power machines. Chem. Pet. Eng. 2020, 55, 733–742. [Google Scholar] [CrossRef]
  7. Seshaiah, N.; Sahoo, R.K.; Sarangi, S.K. Theoretical and experimental studies on oil injected twin-screw air compressor when compressing different light and heavy gases. Appl. Therm. Eng. 2010, 30, 327–339. [Google Scholar] [CrossRef]
  8. Vittorini, D.; Cipollone, R. Energy saving potential in existing industrial compressors. Energy 2016, 102, 502–515. [Google Scholar] [CrossRef]
  9. Qv, D.; Dong, B.; Cao, L.; Ni, L.; Wang, J.; Shang, R.; Yao, Y. An experimental and theoretical study on an injection-assisted air-conditioner using R32 in the refrigeration cycle. Appl. Energy 2017, 185, 791–804. [Google Scholar] [CrossRef]
  10. Wang, X.; Hwang, Y.; Radermacher, R. Two-stage heat pump system with vapor-injected scroll compressor using R410A as a refrigerant. Int. J. Refrig. 2009, 32, 1442–1451. [Google Scholar] [CrossRef]
  11. Plastinin, P.I. Reciprocating Compressors; Tom 1: Moscow, Russia, 2006; 456p. [Google Scholar]
  12. Shcherba, V.E.; Bolshtyanskii, A.P.; Kaigorodov, S.Y.; Kuzeeva, D.A. Benefits of Integrating Displacement Pumps and Compressors. Russ. Eng. Res. 2016, 36, 174–178. [Google Scholar] [CrossRef]
  13. Zhang, C.; Li, P.Y.; van de Ven James, D. Design analysis of a liquid-piston compression chamber with application to compressed air energy storage. Appl. Therm. Eng. 2016, 101, 704–709. [Google Scholar] [CrossRef] [Green Version]
  14. Patil, V.C.; Acharya, P.; Ro, P.I. Experimental investigation of heat transfer in liquid piston compressor. Appl. Therm. Eng. 2019, 146, 169–179. [Google Scholar] [CrossRef]
  15. ToroK, A.; Petrescu, S.; Popescu, G. Quasi-Izothermal Compressors and Expanders with Liquid Piston. Available online: https://hal.archives-ouvertes.fr/hal-02461111/document (accessed on 16 January 2022).
  16. Piya, C.; Sircar, I.; Van de Ven James, D.; Olinger, D.J. Numerical Modeling of Liquid Piston Gas Compression. In Proceedings of the ASME International Mechanical Engineering Congress and Exposition, Lake Buena Vista, FL, USA, 13–19 November 2009; pp. 507–517. [Google Scholar]
  17. Kermani, N.A.; Rokni, M. Heat transfer analysis of liquid piston compressor for hydrogen applications. Int. J. Hydrog. Energy 2015, 40, 11522–11529. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, C.; Yan, B.; Wieberdink, J. Thermal analysis of a compressor for application to compressed air energy storage. Appl. Therm. Eng. 2014, 73, 1402–1411. [Google Scholar] [CrossRef]
  19. Saadat, M.; Li, P.Y.; Sirnon, T.W. Optimal trajectories for a liqu.id piston compressor/expander in a compressed air energy storage system with consideration of heat transfer and friction. In Proceedings of the American Control Conference (ACC), Montréal, QC, Canada, 27–29 June 2012; IEEE: Piscataway, NJ, USA, 2012; pp. 1800–1805. [Google Scholar]
  20. Shcherba, V.E.; Zanin, A.V.; Khrapskii, S.F. Fluid motion dynamics in a two-stage piston hybrid power machine. Chem. Pet. Eng. 2020, 56, 534–542. [Google Scholar] [CrossRef]
  21. Shcherba, V.E.; Shalai, V.V.; Pustovoi, N.; Zanin, A.V. On Profiling of the Working Cavity of the Second Stage in a Hybrid Piston Volumetric Power Machine. Chem. Pet. Eng. 2020, 56, 125–136. [Google Scholar] [CrossRef]
  22. Vimmr, J. Mathematical modelling of compressible inviscid fluid flow through a sealing gap in the screw compressor. Math. Comput. Simul. 2003, 61, 187–197. [Google Scholar] [CrossRef]
  23. James, D.; Van de Ven Perry, Y.L. Liquid piston gas compression. Appl. Energy 2009, 86, 2183–2191. [Google Scholar] [CrossRef]
  24. Xu, W.; Du, Z.; Wang, X.; Cai, M.; Jia, G.; Shi, Y. Isothermal piston gas compression for compressed air energy storage. Int. J. Heat Mass Transf. 2020, 155, 119779. [Google Scholar] [CrossRef]
  25. Bradshaw, C.R.; Groll, E.A.; Garimella, S.V. A comprehensive model of a miniature-scale linear compressor for electronics cooling. Int. J. Refrig. 2011, 34, 63–73. [Google Scholar] [CrossRef] [Green Version]
  26. Kundu, P.K.; Cohen, I.M.; Dowling, D.R. Fluid Mechanics; Academic Press: London, UK, 2015. [Google Scholar]
  27. Zanin, A.; Pavlyuchenko, E.A.; Shcherba, V.E. Numerical and experimental study on fluid compressibility in a two-stage reciprocating pump-compressor. Appl. Therm. Eng. 2021, 194, 117106. [Google Scholar] [CrossRef]
Figure 1. Construction diagram of PHPM (1—cylinder; 2—piston; 3, 4—compressor and pump cavity; 5, 6—discharge and suction liquid nozzles; 7—plunger).
Figure 1. Construction diagram of PHPM (1—cylinder; 2—piston; 3, 4—compressor and pump cavity; 5, 6—discharge and suction liquid nozzles; 7—plunger).
Machines 10 00070 g001
Figure 2. Schematic diagram of a two-stage PHPM without a cylinder cap.
Figure 2. Schematic diagram of a two-stage PHPM without a cylinder cap.
Machines 10 00070 g002
Figure 3. Different profiles of the working chambers of stage 2 of the PHPM ((a)—cylindrical; (b)—combined, made up of cones and cylinders; (c)—combined, made up of two truncated cones).
Figure 3. Different profiles of the working chambers of stage 2 of the PHPM ((a)—cylindrical; (b)—combined, made up of cones and cylinders; (c)—combined, made up of two truncated cones).
Machines 10 00070 g003
Figure 4. Schematic diagram of a two-stage PHPM with a profiled working chamber of stage 2.
Figure 4. Schematic diagram of a two-stage PHPM with a profiled working chamber of stage 2.
Machines 10 00070 g004
Figure 5. The profiling of the working chamber of the second stage: (a) Section of the working chamber of PHPM stage 2 with indication of the main size dimensions; (b) 3D image of the layout of the prototype of PHPM stage 2 (1. Suction line. 2. The suction valve of stage 2 of the compressor (shown conventionally). 3. Valve body. 4. Top flange. 5. Discharge window. 6. Plexiglas cylinder. 7. Liquid piston. 8. Bottom flange with tube. 9. Profiled insert); (c) Image of the PHPM stage 2 prototype.
Figure 5. The profiling of the working chamber of the second stage: (a) Section of the working chamber of PHPM stage 2 with indication of the main size dimensions; (b) 3D image of the layout of the prototype of PHPM stage 2 (1. Suction line. 2. The suction valve of stage 2 of the compressor (shown conventionally). 3. Valve body. 4. Top flange. 5. Discharge window. 6. Plexiglas cylinder. 7. Liquid piston. 8. Bottom flange with tube. 9. Profiled insert); (c) Image of the PHPM stage 2 prototype.
Machines 10 00070 g005
Figure 6. Hydropneumatic diagram of the experimental stand for two-stage PHPM studying. Legend: 1—fluid control valve; 2, 46—tank; 3—electromotor; 4, 54—pump; 5, 34, 45, 47—oil filter; 6, 20, 29—safety valve; 7, 15, 26, 35, 41, 55—pressure sensor; 8—air filter; 9, 24, 40, 51—ratemeter; 10, 14, 16, 31, 37, 44, 52, 56—temperature sensor; 11, 23, 28—manometer; 12—rotation sensor; 13—compressor; 17—oil-and-moisture separator; 18, 21, 25—globe valve; 19—receiver; 22—pressure switch gauge; 27—fluid-power motor; 30—heat exchanger; 32—throttle; 33, 39, 42, 43, 50—ball valve; 36—hydro-pneumatic accumulator; 38—measuring container; 48—fluid level gauge; 49—filler neck strainer; SV1-SV3—safety valve; GV1-GV3—globe valve; PSG1—pressure switch gauge; RC1—receiver; MN1-MN3—manometer; RM1-RM4—ratemeter; OS1—oil-and-moisture separator; TS1-TS8—temperature sensor; PS1-PS7—pressure sensor; CM1—compressor; AF1—air filter; F1-F4—oil filter; FN1—filler neck strainer; T1, T2—tank; P1, P2—pump; RS1—rotation sensor; FM1—fluid-power motor; FV1—fluid control valve; TT1—throttle; HE1—heat exchanger; BV1-BV5—ball valve; AC1—hydro-pneumatic accumulator; FL1—fluid level gauge; MC1—measuring container.
Figure 6. Hydropneumatic diagram of the experimental stand for two-stage PHPM studying. Legend: 1—fluid control valve; 2, 46—tank; 3—electromotor; 4, 54—pump; 5, 34, 45, 47—oil filter; 6, 20, 29—safety valve; 7, 15, 26, 35, 41, 55—pressure sensor; 8—air filter; 9, 24, 40, 51—ratemeter; 10, 14, 16, 31, 37, 44, 52, 56—temperature sensor; 11, 23, 28—manometer; 12—rotation sensor; 13—compressor; 17—oil-and-moisture separator; 18, 21, 25—globe valve; 19—receiver; 22—pressure switch gauge; 27—fluid-power motor; 30—heat exchanger; 32—throttle; 33, 39, 42, 43, 50—ball valve; 36—hydro-pneumatic accumulator; 38—measuring container; 48—fluid level gauge; 49—filler neck strainer; SV1-SV3—safety valve; GV1-GV3—globe valve; PSG1—pressure switch gauge; RC1—receiver; MN1-MN3—manometer; RM1-RM4—ratemeter; OS1—oil-and-moisture separator; TS1-TS8—temperature sensor; PS1-PS7—pressure sensor; CM1—compressor; AF1—air filter; F1-F4—oil filter; FN1—filler neck strainer; T1, T2—tank; P1, P2—pump; RS1—rotation sensor; FM1—fluid-power motor; FV1—fluid control valve; TT1—throttle; HE1—heat exchanger; BV1-BV5—ball valve; AC1—hydro-pneumatic accumulator; FL1—fluid level gauge; MC1—measuring container.
Machines 10 00070 g006
Figure 7. A general view of the prototype installed on the stand: (a) General view of two-stage PHPM with a gas cap (1—Gas cap; 2—Suction line; 3—Stage 1 of a two-stage PHPM; 4—Discharge line; 5—Manometer; 6—Inter -stage pipeline; 7—Stage 2 of two-staged PHPM); (b) Image of the prototype on the stand.
Figure 7. A general view of the prototype installed on the stand: (a) General view of two-stage PHPM with a gas cap (1—Gas cap; 2—Suction line; 3—Stage 1 of a two-stage PHPM; 4—Discharge line; 5—Manometer; 6—Inter -stage pipeline; 7—Stage 2 of two-staged PHPM); (b) Image of the prototype on the stand.
Machines 10 00070 g007
Figure 8. Two-stage piston hybrid power machine with a profiled working chamber: (a) objective design diagram of PHPM stage 2; (b) given design diagram of PHPM stage 2.
Figure 8. Two-stage piston hybrid power machine with a profiled working chamber: (a) objective design diagram of PHPM stage 2; (b) given design diagram of PHPM stage 2.
Machines 10 00070 g008
Figure 9. Dependence of the relative pressure losses along the length, to overcome local resistances and inertia forces on the dynamic viscosity coefficient (1- Δ h ¯ l ; 2- Δ h ¯ ζ ; 3- Δ h ¯ i n Σ ).
Figure 9. Dependence of the relative pressure losses along the length, to overcome local resistances and inertia forces on the dynamic viscosity coefficient (1- Δ h ¯ l ; 2- Δ h ¯ ζ ; 3- Δ h ¯ i n Σ ).
Machines 10 00070 g009
Figure 10. Dependence of relative pressure drops Δ p ¯ 1 ; Δ p ¯ 2 ; Δ p ¯ 3 , on the dynamic viscosity coefficient (1- Δ p ¯ 1 ; 2- Δ p ¯ 2 ; 3- Δ p ¯ 3 ).
Figure 10. Dependence of relative pressure drops Δ p ¯ 1 ; Δ p ¯ 2 ; Δ p ¯ 3 , on the dynamic viscosity coefficient (1- Δ p ¯ 1 ; 2- Δ p ¯ 2 ; 3- Δ p ¯ 3 ).
Machines 10 00070 g010
Figure 11. Dependence of indicator isothermal efficiency, total indicator efficiency and work of friction forces on the dynamic viscosity coefficient (1- η ind . is ; 2- η ind . f ; 3- L f r ).
Figure 11. Dependence of indicator isothermal efficiency, total indicator efficiency and work of friction forces on the dynamic viscosity coefficient (1- η ind . is ; 2- η ind . f ; 3- L f r ).
Machines 10 00070 g011
Figure 12. Dependence of the relative losses along the length on the coefficient of dynamic viscosity at different speeds of the crankshaft.
Figure 12. Dependence of the relative losses along the length on the coefficient of dynamic viscosity at different speeds of the crankshaft.
Machines 10 00070 g012
Figure 13. Dependence of relative pressure loss Δ p ¯ 1 ; Δ p ¯ 2 ; Δ p ¯ 3 , on the diameter of the supply pipeline (1- Δ p ¯ 1 ; 2- Δ p ¯ 2 ; 3- Δ p ¯ 3 ).
Figure 13. Dependence of relative pressure loss Δ p ¯ 1 ; Δ p ¯ 2 ; Δ p ¯ 3 , on the diameter of the supply pipeline (1- Δ p ¯ 1 ; 2- Δ p ¯ 2 ; 3- Δ p ¯ 3 ).
Machines 10 00070 g013
Figure 14. Dependence of the relative head losses along the length to overcome local resistances and inertial forces on the diameter of the supply pipeline (1- Δ h ¯ l ; 2- Δ h ¯ ζ ; 3- Δ h ¯ i n Σ ).
Figure 14. Dependence of the relative head losses along the length to overcome local resistances and inertial forces on the diameter of the supply pipeline (1- Δ h ¯ l ; 2- Δ h ¯ ζ ; 3- Δ h ¯ i n Σ ).
Machines 10 00070 g014
Figure 15. Dependence of the indicator isothermal efficiency, total indicator efficiency, and work of friction forces on the liquid density (1- η ind . is ; 2- η ind . f . ; 3- L f r ).
Figure 15. Dependence of the indicator isothermal efficiency, total indicator efficiency, and work of friction forces on the liquid density (1- η ind . is ; 2- η ind . f . ; 3- L f r ).
Machines 10 00070 g015
Table 1. Main dimensions PHPM.
Table 1. Main dimensions PHPM.
Indicator DescriptionValues
The diameter of the first stage pistondp = 0.050 m
Diameter of the rod of the first stagedr = 0.036 m
Full stroke of the 1st stage pistonSh = 0.050 m
Initial coordinate of the liquid in the second stageZ20 = 0.100 m
Diameter of the supply pipeline dsp = 0.016 m
The length of the supply pipelinelp = 0.5 m
Linear dead spaceSd = 0.0025 m
Turning radius of the supply pipelineRp = 0.5 m
Table 2. Size dimensions of self-acting valves of stage 2 compressor section.
Table 2. Size dimensions of self-acting valves of stage 2 compressor section.
Indicator DescriptionValues
Maximum lift height of the shut-off element of the suction valve in the compressor sectionhscvmax = 0.0015 m
The width of the passage in the seat of the suction valve of the compressor sectiondscvc = 0.023 m
Spring stiffness of the suction valve in the compressor sectioncsscv = 3200 N/m
Maximum lift height of the shut-off element of the discharge valve in the compressor sectionhdvmax = 0.0015 m
The width of the passage in the seat of the discharge valve of the compressor sectionddvc = 0.023 m
Spring stiffness of the discharge valve in the compressor sectioncsdv = 3200 n/m
Table 3. Main geometrical and operational parameters of two-stage PHPM.
Table 3. Main geometrical and operational parameters of two-stage PHPM.
Indicator DescriptionValues
The diameter of stage 1 pistond1 = 0.050 m
Diameter of the stem of stage 1dst = 0.036 m
Full stroke of stage 1 pistonSh1 = 0.050 m
The number of revolutions of the crankshaftnrev = 500 rpm
Suction pressurePsc2 = 5 × 105 Pa
Discharge pressurePd2 = 25 × 105 Pa
The average temperature of the walls of the surface of the working chamber of stage 2 T ¯ w = 330 K
Intake air temperatureTsc2 = 310 K
Initial coordinate of the liquid in stage 2Z20 = 0.100 m
Diameter of the supply pipelinedpp = 0.016 m
The length of the supply pipelinelpp = 0.5 m
Linear dead spaceSds = 0.0025 m
Turning radius of the supply pipelineRpp = 0.5 m
The radius of the conical insert in the design scheme of stage 2 of the PHPM (see Figure 8)R21 = 0.003 m, R11 = 0.015 m
The radii of the conical working chamber in the given design scheme of stage 2 of the PHPM:R1 = 0.0132 m, R2 = 0.01977 m
The angle of inclination of the generatrix of the truncated coneα = 0.235 rad
Truncated cone heightS1 = 0.0272 m
The maximum lifting height of the closure body of the discharge valvehdmax = 0.001 m
The maximum lifting height of the closure body of the suction valvehscmax = 0.001 m
The width of the passage in the saddle of the suction valvedscv = 0.09 m
The width of the passage in the saddle of the discharge valveddv = 0.09 m
Spring rate of the discharge valveCsprd = 16,280 N/m
Spring rate of the suction valveCsprsc = 5980 N/m
Maximum clearance between the saddle and the closure body of the suction valveδsc = 1 × 10−8 m
Minimum clearance between the saddle and the closure body of the discharge valveδd = 1 × 10−8 m
Roughness of the inner surface of the connecting pipelineK = 0.000005 m
Compressed gasair
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shcherba, V.; Pavlyuchenko, E.; Nosov, E.; Bulgakova, I. Analysis of the Influence of Thermophysical Properties of a Droplet Liquid on the Work Processes of a Two-Stage Piston Hybrid Power Machine. Machines 2022, 10, 70. https://doi.org/10.3390/machines10020070

AMA Style

Shcherba V, Pavlyuchenko E, Nosov E, Bulgakova I. Analysis of the Influence of Thermophysical Properties of a Droplet Liquid on the Work Processes of a Two-Stage Piston Hybrid Power Machine. Machines. 2022; 10(2):70. https://doi.org/10.3390/machines10020070

Chicago/Turabian Style

Shcherba, Victor, Evgeniy Pavlyuchenko, Evgeniy Nosov, and Irina Bulgakova. 2022. "Analysis of the Influence of Thermophysical Properties of a Droplet Liquid on the Work Processes of a Two-Stage Piston Hybrid Power Machine" Machines 10, no. 2: 70. https://doi.org/10.3390/machines10020070

APA Style

Shcherba, V., Pavlyuchenko, E., Nosov, E., & Bulgakova, I. (2022). Analysis of the Influence of Thermophysical Properties of a Droplet Liquid on the Work Processes of a Two-Stage Piston Hybrid Power Machine. Machines, 10(2), 70. https://doi.org/10.3390/machines10020070

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop