Next Article in Journal
Influence of SST in Low Latitudes on the Arctic Warming and Sea Ice
Previous Article in Journal
Walleye Pollock Gadus chalcogrammus, a Species with Continuous Range from the Norwegian Sea to Korea, Japan, and California: New Records from the Siberian Arctic
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pipeline Pressure Loss in Deep-Sea Hydraulic Systems Considering Pressure-Dependent Viscosity Change of Hydraulic Oil

1
College of Mechanical and Electrical Engineering, Central South University, Changsha 410083, China
2
State Key Laboratory of High Performance Complex Manufacturing, Central South University, Changsha 410083, China
3
State Key Laboratory of Exploitation and Utilization of Deep-Sea Mineral Resources, Changsha Research Institute of Mining & Metallurgy, Changsha 410012, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2021, 9(10), 1142; https://doi.org/10.3390/jmse9101142
Submission received: 27 September 2021 / Revised: 8 October 2021 / Accepted: 11 October 2021 / Published: 18 October 2021
(This article belongs to the Section Ocean Engineering)

Abstract

:
The high ambient pressure in deep-sea conditions greatly increases the viscosity of hydraulic oil and then the pipeline pressure loss in deep-sea hydraulic systems. Large pipeline pressure loss can lead to a further change of viscosity on the basis of the viscosity increase caused by the ambient pressure when the hydraulic oil flows through the pipeline. Therefore, the classic Poiseuille’s law can no longer accurately calculate the pipeline pressure loss in deep-sea conditions since it treats the viscosity as a constant. In this paper, based on laminar flow theory and the viscosity-pressure characteristics of hydraulic oil, a novel equation for pipeline pressure loss is proposed, in which the viscosity change when flowing through the pipeline is taken into account. A CFD (Computational Fluid Dynamics) model of a pipeline in the deep-sea hydraulic system has been established, and CFD simulations have been conducted to verify the correctness of the proposed equation for pipeline pressure loss. Theoretical analysis shows that the proposed novel equation for pipeline pressure loss is equivalent to the classic Poiseuille’s law when the pipeline pressure loss or the viscosity change is low. The research results in this paper can provide theoretical support for work efficiency optimization, load capacity improvement, and precise control of deep-sea operation equipment or deep-sea hydraulic systems.

1. Introduction

Various deep-sea operation equipment, such as human occupied vehicles (HOVs) [1], remotely operated vehicles (ROVs) [2], submarine mining vehicles [3,4], autonomous underwater vehicles (AUVs) [5], and submarine cable trenchers [6], have been successfully designed and developed. In addition, with them, humans have carried out much scientific exploration and resource development in the deep-sea environment. These high-power deep-sea operation equipment are usually powered by a motor-driven hydraulic system [7,8], and the power is provided by the mother ship on the sea surface through an umbilical cable, as shown in Figure 1. The loss in each power transmission link is also indicated in Figure 1.
The deep-sea hydraulic system is equipped with a pressure compensator [9,10,11], which introduces the seawater ambient pressure that increases with depth into the hydraulic system. Then, the internal and external pressures of the hydraulic system tend to be balanced. Therefore, the components in the deep-sea hydraulic system do not need to be designed with a special pressure-resistant structure, and the mature technology in the land hydraulic system can be transplanted and applied. As a result, the development and design costs of the deep-sea hydraulic system can be greatly reduced.
However, the viscosity of hydraulic oil increases exponentially with pressure. Hydraulic oil is the energy transfer medium of the hydraulic system. Therefore, the increase in its viscosity would cause the loss of each transmission link in the deep-sea hydraulic system to increase significantly compared with those in land conditions. Furthermore, the working efficiency, load driving capability, and precise control of deep-sea operation equipment are seriously affected. Some scholars have researched the influence on deep-sea hydraulic systems caused by the viscosity increase of hydraulic oil.
Cao et al. [12] established a deep-sea hydraulic power unit model based on the linear varying parameters modeling method, which considered the compressibility and viscosity of hydraulic oil affected by ambient pressure in the deep-sea environment. Based on the established model, simulation analysis was carried out to predict the dynamic performance of the hydraulic power unit, which was verified by experimental tests. Tian et al. [13] developed a mathematical model of the pressure control variable displacement pump, which is the power source of the deep-sea hydraulic manipulator. In the model, the increased viscosity of the hydraulic oil caused by the ambient pressure of the seawater was considered. The influence of the viscosity on the dynamic characteristics and performance of the pump was simulated and corresponding experiments were conducted. Tian et al. [14] tested the change in hydraulic oil viscosity with pressure and established a nonlinear model of an underwater manipulator that took into account the increase in hydraulic oil viscosity caused by ambient pressure. Then, the influence of ambient pressure on the dynamic performance of the manipulator was studied, and the online test under the ambient pressure of 115 MPa was carried out.
In the performance analysis of the land hydraulic system, the dynamic process and pressure loss of the connecting pipeline are ignored to simplify the model [15,16]. Under the effect of the high ambient pressure in the deep sea, the viscosity of hydraulic oil increases significantly, and then the pipeline pressure loss cannot be ignored directly. Deep-sea operation equipment usually has multiple actuators in concert, and many actuators are distributed far away from the hydraulic pump. As a result, there are numerous long-distance pipelines in the deep-sea hydraulic system. Therefore, the pipeline pressure loss has a significant impact on the operation of deep-sea operation equipment. The flow in the pipeline is generally a typical Poiseuille flow, and many scholars have carried out related studies on the Poiseuille flow or pipeline pressure loss.
Savić et al. [17] developed a numerical model for determining the pressure losses within flat pipelines. In this numerical model, the influence of pressure and temperature on the density and viscosity of the oil was taken into account. Ganji et al. [18] utilized the Arrhenius model to describe the temperature-viscosity characteristics of the fluid and researched the steady-state Hagen-Poiseuille flow in a circular pipe. Housiadas and Georgiou [19] assumed that density and the viscosity of the fluid have a linear relationship with pressure, and then proposed a novel solution for Poiseuille flows. Based on a Runge–Kutta–Fehlberg integration scheme with shooting technique, Makinde et al. [20] investigated a magnetohydrodynamic Couette-Poiseuille flow between two parallel plates in a rotating permeable channel, in which the temperature-viscosity characteristics of the fluids was taken into account. Lee et al. [21] studied how surface roughness affects the flow characteristics, and the result showed that the roll mode for the Couette-Poiseuille flow over a rough wall can be inhibited by the surface roughness. Hong-Xiang [22] calculated the pipeline pressure loss in the deep-sea hydraulic system, taking into account the increase in viscosity caused by high pressure and low temperature in deep-sea conditions. Then, a variable gain control algorithm of the underwater manipulator based on the pressure loss was developed, which was verified through simulation and experimental testing. With numerical methods, Luo et al. [23] studied the propagation of Taylor vortices in Taylor-Couette-Poiseuille flow considering the effects of an abruptly contracting and expanding annular gap. Li and Wu [24] researched the deformation and leakage mechanisms at hydraulic clearance fit in the deep-sea extreme environment, in which the viscosity-pressure characteristics of hydraulic oil was considered. With the Chebyshev collocation and the Galerkin methods, B.M. and Shivakumara [25] researched how the uniform vertical throughflow affects the stability of Poiseuille flow in a Newtonian fluid-saturated Brinkman porous medium.
Generally, the pressure loss in the pipeline can be determined by Poiseuille’s law.
Δ p = 128 η L Q π D 4 ,
where Δ p is the pressure loss, η is the viscosity, L and D are the length and the inner diameter of the pipeline, respectively, and Q is the flow rate.
The viscosity-pressure characteristics of hydraulic oil can be fitted by the Barus formula [26], which is expressed as Equation (2).
η = η 0 e α p o ,
where η 0 is the initial viscosity of the oil at standard atmospheric pressure, α is the viscosity-pressure index, and p o is the outlet pressure of the pipeline. Due to the function of the pressure compensator, for an inflow pipeline that connects the hydraulic pump and the hydraulic actuator, p o is the sum of the ambient pressure and the load pressure. As for the return pipeline that connects the hydraulic actuator and the oil tank, p o equals ambient pressure.
With Equations (1) and (2), the pressure loss considering the viscosity increase caused by the extremely high ambient pressure in the deep sea is obtained as follows.
Δ p = 128 η 0 L Q e α p o π D 4 .
The fluid viscosity in Equation (1) is a constant value. In Equation (3), the viscosity takes into account the change caused by ambient pressure, but it is still a constant value when the hydraulic oil flows through the pipeline. When in land conditions or shallow sea conditions, the ambient pressure is tiny, and so are the outlet pressure of the pipeline p o and the pressure loss calculated by Equation (3). Therefore, the viscosity of the hydraulic oil changes minimally, and the calculation using Equation (3) is still reasonable.
However, when working in deep-sea conditions, the viscosity of hydraulic oil is greatly increased compared with land working conditions, so the pipeline pressure loss also increases significantly, which can exceed 10 MPa. Such a huge internal loss in the system seriously reduces the working efficiency of deep-sea operation equipment or deep-sea hydraulic systems. Generally, the maximum output pressure of the hydraulic pump usually does not exceed 30 MPa. Therefore, the pipeline pressure loss of up to 10 MPa reduces the load capacity of the deep-sea hydraulic system by an incredible 33%. In addition, proportional valves or servo valves are usually used in hydraulic systems to achieve high-precision control. For a proportional valve or servo valve with a given control signal, the flow rate that flows through it to the hydraulic actuator is related to the pressure difference between the inlet and outlet of the valve [14,16]. When the deep-sea hydraulic system works at different depths, the viscosity of the hydraulic oil changes and so does the pipeline pressure loss, which results in the variations of the pressure difference between the inlet and outlet of the valve and, furthermore, the flow rate to the hydraulic actuator. As a consequence of this series of changes, the hydraulic actuators have different output responses at different underwater depths, even with the same control signal, which is very unfavorable for the precise control. The above analysis shows that pipeline pressure loss has a significant influence on the work efficiency, load capacity, and precise control of deep-sea operation equipment or deep-sea hydraulic systems.
Besides, a pressure loss of 10 MPa can cause a significant change in the viscosity of the hydraulic oil when flowing through the pipeline. Taking α as 2.2 × 10 8 Pa 1 , the change percentage in hydraulic oil viscosity caused by a pressure loss of 10 MPa can be calculated using Equation (2), which is
η 0 e α ( p o + Δ p ) η 0 e α p o η 0 e α p o × 100 % = e α Δ p 1 × 100 % = e 2.2 × 10 8 P a 1 × 10 MPa 1 × 100 % = 24.61 % .
Equation (4) indicates that the viscosity of the hydraulic oil flowing through the pipeline changes by as much as 24.61%, which means that the classical Poiseuille’s law, in which the viscosity is treated as a constant value, can no longer accurately calculate the pipeline pressure loss in the deep-sea hydraulic system.
However, few studies involved pipeline pressure loss in deep-sea hydraulic systems. In addition, the corresponding research only considered the increase in viscosity caused by the ambient pressure of the deep-sea environment, instead of the viscosity change during the flowing process, as in the research of Tian et al. [14] and Hong-Xiang [22]. Therefore, it is necessary to conduct an in-depth study on the pipeline pressure loss in the deep-sea hydraulic system, in which the viscosity change of the hydraulic oil during the flowing process in the pipeline should be taken into account.
In this paper, based on laminar flow theory and the viscosity-pressure characteristics of hydraulic oil, a novel equation for pipeline pressure loss is derived, in which fluid viscosity is variable. Then, a CFD model of a pipeline in the deep-sea hydraulic system is established, and simulations are carried out to verify the correctness of the proposed novel equation for pipeline pressure loss. Finally, the proposed novel equation for pipeline pressure loss is analyzed and discussed, and it is also compared with the classic Poiseuille’s law. The research results in this paper can provide theoretical support for work efficiency optimization, load capacity improvement, and precise control of deep-sea operation equipment or deep-sea hydraulic systems.

2. Theory Analysis

Laminar flow has less energy loss than turbulent flow. Therefore, in the design of the hydraulic system, according to the volume flow rate and the set maximum average flow velocity, the pipe diameter can then be determined so that the hydraulic oil flows in a laminar flow state in the pipeline [27]. Therefore, in this paper, the laminar fluid domain in the circular pipeline is used for analysis, and the analysis model is shown as Figure 2. A micro-cylinder with a length of d l and a radius of r located on the axis of the pipeline is selected as the research object. The outlet pressure of the micro-cylinder is p, the pressure loss when flowing through the micro-cylinder is d p , and the flow velocity on the surface of the micro-cylinder is v.
The pipeline model in Figure 2 is placed horizontally. The diameter of the pipeline usually does not exceed 100 mm. Based on the hydraulic oil density of 850 kg / m 3 , the pressure difference between the upper and lower hydraulic oil in the pipeline due to gravity is
850 kg / m 3 × 9 . 8 N / kg × 100 mm = 8 . 33 × 1 0 4 MPa ,
which is very low. Therefore, for the horizontal model shown in Figure 2, gravity can be ignored.
When the pipeline is placed vertically, the pressure difference due to gravity in a pipeline with a length of 5 m is
850 kg / m 3 × 9 . 8 N / kg × 5 m = 4 . 165 × 1 0 2 MPa ,
which is far less than the pressure loss due to the viscosity of the hydraulic oil. Therefore, the results derived from the horizontal pipeline model in this paper are also applicable to vertical pipelines. Furthermore, it is obviously applicable to inclined pipes, as well.
The influence of temperature on viscosity is not considered in this analysis, and it is assumed that the hydraulic system is operating at a design condition of 40 °C.
When in the state of laminar flow, the fluids flow in the axial direction without interfering with each other. Therefore, in the radial direction, the forces of the micro-cylinder are balanced. In the axial direction, the micro-cylinder is subjected to the forces caused by pressure at both ends and the viscous friction force caused by the rest of the fluid on its surface, which are expressed as follows.
F i = π r 2 ( p + d p ) , F o = π r 2 p , F f = A · η d v d r = 2 π r η · d l · d v d r .
where F i and F o are the forces caused by pressure at the inlet and outlet, respectively, F f is the viscous friction force, and A is the side area of the micro-cylinder.
When the hydraulic oil in the pipeline is in the state of laminar flow, the axial forces of the micro-cylinder are also balanced. Then,
F i = F o + F f ,
or
π r 2 ( p + d p ) = π r 2 p + 2 π r η · d l · d v d r ,
or
d v = 1 2 η d p d l r d r .
Substituting the viscosity-pressure characteristics expressed by Equation (2) into Equation (8) and replacing p o with p, we have
d v = 1 2 η 0 e α p d p d l r d r .
The inner radius of the pipeline is R, and the flow velocity here is 0. With these boundary conditions, we integrate and process Equation (9) as follows.
0 v d v = R r 1 2 η 0 e α p d p d l r d r = 1 2 η 0 e α p d p d l R r r d r ,
or
v = R 2 r 2 4 η 0 e α p d p d l ,
or
4 v η 0 R 2 r 2 d l = e α p d p .
The pipeline length is L and the outlet pressure of the pipeline is p o . The inlet pressure is the sum of the outlet pressure p o and the pressure loss Δ p . With these boundary conditions, we integrate and process Equation (12) as follows.
L 0 4 v η 0 R 2 r 2 d l = 4 v η 0 R 2 r 2 L 0 d l = p o + Δ p p o e α p d p ,
or
4 v η 0 L R 2 r 2 = 1 α e α p o e α p o + Δ p .
The flow velocity v can be obtained as follows by transforming Equation (14).
v = R 2 r 2 4 α η 0 L e α p o e α p o + Δ p .
According to Equation (15), the flow velocity v in the pipeline shows a quadratic parabolic distribution with the radius r, and the maximum flow velocity v m a x , which is expressed as follows, occurs when r is 0, or on the axis.
v m a x = R 2 4 α η 0 L e α p o e α p o + Δ p .
Based on Equation (15), the total flow rate of hydraulic oil in the pipeline can be derived. Figure 3 shows the pipeline cross-section. A micro-ring with radius r and width d r in the cross-section is taken as the research object.
The flow rate through this micro-ring is
d Q = v d A = 2 π v r d r ,
where d A is the area of the micro-ring.
Substituting Equation (15) into Equation (17), we have
d Q = π 2 α η 0 L e α p o e α p o + Δ p R 2 r 2 r d r .
Then, the flow rate Q can be deduced as follows by integrating Equation (18).
Q = 0 R π 2 α η 0 L e α p o e α p o + Δ p R 2 r 2 r d r = π 2 α η 0 L e α p o e α p o + Δ p 0 R R 2 r 2 r d r = π 2 α η 0 L e α p o e α p o + Δ p 0 R R 2 r d r 0 R r 3 d r = π 2 α η 0 L e α p o e α p o + Δ p 1 2 R 4 1 4 R 4 = π R 4 8 α η 0 L e α p o e α p o + Δ p = π D 4 128 α η 0 L e α p o e α p o + Δ p .
Dividing the flow rate Q by the flow area π R 2 , the average flow rate v a v e can be obtained as follows.
v a v e = Q π R 2 = 1 π R 2 π R 4 8 α η 0 L e α p o e α p o + Δ p = R 2 8 α η 0 L e α p o e α p o + Δ p .
According to Equations (16) and (20), the average flow velocity v a v e is half of the maximum flow velocity v m a x , namely
v a v e = v m a x 2 .
In the textbooks or handbooks of the hydraulic system [27,28,29,30], when analyzing the Poiseuille flow with a constant viscosity in a circular pipeline, the law expressed by Equation (21) is also described. Therefore, it can be known that the change in viscosity when flowing through the pipeline does not affect the velocity distribution.
With Equations (16), (20) and (21), Equation (15) describing the velocity distribution can be rewritten as follows.
v = 2 Q π R 4 R 2 r 2 .
Finally, the pressure loss considering variable viscosity can be obtained as follows by transforming Equation (19).
Δ p = 1 α ln ( 1 128 α η 0 L Q e α p o π D 4 ) .
It can be known from Equation (23) that the pressure loss considering variable viscosity is related to the viscosity-pressure index α and initial viscosity η 0 of hydraulic oil, the length L and the inner diameter D of the pipeline, the flow rate Q, and the outlet pressure p o .
In order to verify the correctness of the proposed novel equation for pipeline pressure loss, a CFD [31,32] model of an inflow pipeline in the deep-sea hydraulic system is established and simulations are carried out in the following subsections.

3. CFD Model and Settings

Due to the symmetry of the circular pipeline, a quarter of the fluid domain in the pipeline is established in the CFX software to save computing resources, as shown in Figure 4a. The inner diameter D and length L of the pipeline are set as 4 mm and 2 m, respectively, indicating that it is a significant slender pipeline. In order to obtain better mesh quality, the fluid domain is divided into three regions containing a cuboid with a square bottom for mesh generation, as shown in Figure 4b. For the three regions on the end face, each edge or curve is divided into 20 equal parts. In addition, the model is divided into 4000 equal parts along the axial direction to reduce the aspect ratio of the elements and so as to improve the quality. Then, the number of generated hexahedral elements and nodes are 4,800,000 and 5,045,261, respectively.
The fluid is set as a 22# hydraulic oil, whose kinematic viscosity is 22cSt at a standard atmospheric pressure and a working temperature of 40 °C. Its density is set as 850 kg/m 3 . Then, its initial dynamic viscosity η 0 is 22 cSt × 850 kg / m 3 = 0.0187 Pa · s .
A user-defined expression is created to import the viscosity–pressure characteristics expressed by Equation (2). For hydraulic oils, the viscosity–pressure index α is within 1.5 × 10 8 Pa 1 to 3.5 × 10 8 Pa 1 [27,29,30]. In this CFD model, set α as 2.2 × 10 8 Pa 1 , which is the viscosity-pressure index of the HFD fire-resistant hydraulic oil [30].
The flow rate Q is set as 5 L/min, and its equivalent mass flow rate is
5 L / min × 850 kg / m 3 = 0.07083333 kg · s .
Then, the inlet of the fluid domain is set as a mass flow rate inlet with one-quarter of the calculated mass flow rate since it is a quarter model. The ambient pressure at the bottom of the Mariana Trench is about 110 MPa, and, if the load pressure is 15 MPa, the outlet pressure p o of the fluid domain is set as 125 MPa. The cut plane and the cylindrical surface are set as the symmetry boundary and the wall boundary, respectively.
According to the aforementioned parameter settings and the hydraulic oil viscosity of 125 MPa at the outlet of the pipeline, the Reynolds number is 77.08, indicating laminar flow. Then, the analysis type is set as steady-state and the turbulence type is set as none, which means laminar flow. The minimum and maximum iterations are 50 and 100, respectively. The convergence criteria is an RMS (Root Mean Square) type residual with a target value of 1 × 10 4 .
To improve the precision and reliability of the simulation results, a mesh independence study was conducted based on the aforementioned simulation settings. Table 1 lists the number of total mesh elements in six cases. How the pipeline pressure loss, the average flow velocity at the outlet, and the calculation time change with the number of elements are shown in Figure 5. According to Table 1 and Figure 5, the pressure loss and the average flow velocity can be considered stable when the number of mesh elements exceeds 4,800,000. The calculation time increases with the number of mesh elements and is about 1200 s or about 20 min in the case with 4,800,000 elements. Therefore, the CFD simulation of a model with 4,800,000 elements has satisfied calculation accuracy and fast calculation speed at the same time.
Based on the established CFD model, the pipeline flow simulation of the deep-sea hydraulic system is carried out. At the same time, a simulation that does not consider the change in hydraulic oil viscosity, which is based on Equation (3), is also carried out for comparison.

4. CFD Results

Dynamic viscosity, pressure distribution with variable viscosity, and pressure distribution with constant viscosity are shown in Figure 6, Figure 7 and Figure 8, respectively. Figure 6 shows that the established CFD successfully takes into account the viscosity change during the flow process. Calculated with the viscosities of the inlet and outlet in Figure 6, the viscosity of the hydraulic oil changes by 20.79% when flowing through the pipeline. Comparing Figure 7 and Figure 8, it can be known that the pressure loss considering the viscosity change during the flow process is greater than in the case where the viscosity is constant.
The simulated pressure loss can be obtained by reading the average value of the inlet pressure and subtracting the set outlet pressure of 125 MPa from it, as shown in Table 2. The theoretical values calculated according to Equations (3) and (23) are also listed.
From the data in Table 2, it can be known that the theoretical results are in good agreement with the simulation results. According to the setting in Section 3, the load pressure is 15 MPa. Therefore, it can be calculated that the load capacity has dropped by 39.64%, which indicates that the operating performance of the deep-sea hydraulic system is greatly affected. In addition, the pressure loss that takes the viscosity change into account is 9.6% larger than the pressure loss in which the viscosity is constant. This means, in the deep-sea environment, the viscosity change during the flow process in the pipeline has a significant impact on the pipeline pressure loss.
The flow velocity distribution is shown in Figure 9, in which the theoretical distribution curve according to Equation (22) is also plotted. Figure 9 shows that, regardless of whether the viscosity change during the flow process is considered, the flow velocity distribution is in good agreement with each other, which is consistent with the theoretical analysis in Section 2. The velocity distribution at the inlet or at a position less than 0.01 m away from the inlet is quite different from the theoretical distribution. This is caused by the uniform input of fluid at the flow inlet with an average flow velocity. The difference only exists in positions less than 0.01 m away from the inlet, whose range is extremely low compared to the pipeline length of 2 m. As for a position more than 0.01 m away from the inlet, the velocity distribution is in good agreement with the theoretical distribution.
In order to more comprehensively verify the proposed novel equation for pipeline pressure loss expressed by Equation (23), more CFD simulation calculations with variable input parameters are carried out, whose inputs and pressure losses are all listed in Table 3 and Table 4. According to the data in Table 3 and Table 4, all simulated pressure losses are in good agreement with the theoretical pressure losses calculated by Equation (23), and the maximum error is −1.80%. All the CFD results justify the correctness of the proposed novel equation for pipeline pressure loss in this paper.

5. Discussions

The CFD calculation results in the last section fully prove the correctness of the proposed novel equation for pipeline pressure loss. In this section, this novel equation will be further analyzed and discussed.
According to Equation (23), the pressure loss considering variable viscosity is related to hydraulic oil properties, pipeline dimensions, and working conditions. Because the expression contains logarithmic and exponential operations, it becomes difficult to directly analyze the relationship between pressure loss and various parameters. Therefore, based on Equation (23) and the set values in Section 3, curves are drawn to show how the pressure loss varies with each parameter. In addition, the increment percentage of pressure loss compared to classical Poiseuille’s law, which can be calculated with Equations (3) and (23), is plotted. In addition, so is the change percentage of viscosity calculated by Equation (4).

5.1. Influence of Hydraulic Oil Properties

It can be known from Figure 10 and Figure 11 that the pressure loss increases with the increase in the viscosity-pressure index or the initial viscosity. The greater the viscosity index is, the faster the pressure loss increases, while the pressure loss changes approximately linearly with the initial viscosity. The increment percentage of pressure loss compared to classical Poiseuille’s law and the change percentage of viscosity show a similar pattern of change.
According to the pressure loss changes shown in Figure 10 and Figure 11, for deep-sea hydraulic systems, it is recommended to use hydraulic oil with a low viscosity-pressure index and a low initial viscosity, which can reduce pipeline pressure loss, thereby improving work efficiency and load capacity. In addition, comparing Figure 10 and Figure 11, it can be known that the viscosity-pressure index has a greater impact on pressure loss, so hydraulic oils with a low viscosity-pressure index should be given priority.
According to the percentage curves in Figure 10 and Figure 11, it can be known that, when the viscosity of the hydraulic oil flowing through the pipeline changes minimally, the pressure loss calculated by Equation (23) is less different from the classic Poiseuille’s law. When the viscosity changes greatly, the difference between the two calculation methods becomes larger. According to Figure 10, when the viscosity-pressure index is 3 × 10 8 Pa 1 , the viscosity changes by 172.30%, and the difference in pressure loss calculated by the two methods reaches 58.31%.

5.2. Influence of Pipeline Dimensions

According to Figure 12, the pressure loss increases approximately linearly with the increase in the pipeline length, and so do the two percentage curves. As for the influence of pipeline inner diameter that is shown in Figure 13, the pressure loss decreases sharply as the pipeline inner diameter increases when the pipeline inner diameter is less than 12 mm. When the pipeline inner diameter is greater than 12 mm, the pressure loss changes minimally. The two percentage curves in Figure 13 also show a similar pattern of change.
Then, for deep-sea hydraulic systems, it is recommended to use pipelines with short lengths and large inner diameters to reduce pipeline pressure loss, thereby improving work efficiency and load capacity. In addition, according to Figure 13, an excessively large pipeline inner diameter contributes minimally to reducing pipeline pressure loss but increases the weight and space of the deep-sea hydraulic system.
The percentage curves in Figure 12 and Figure 13 also indicate that, when the viscosity of the hydraulic oil changes minimally, the pressure loss calculated by Equation (23) is less different from the classic Poiseuille’s law.

5.3. Influence of Working Conditions

The maximum outlet pressure of 140 MPa in Figure 14 is obtained from the ambient pressure of 110 MPa at the bottom of the Mariana Trench plus the load pressure of 30 MPa. It can be known from Figure 14 and Figure 15 that the pressure loss increases with the increase in the outlet pressure or the flow rate. The greater the outlet pressure is, the faster the pressure loss increases, while the pressure loss changes approximately linearly with the flow rate. The percentage curves show a similar pattern of change.
According to Figure 14 and Figure 15, for deep-sea hydraulic systems, attention should be paid to the working conditions of large working flow rate or large pipeline outlet pressure, in which large pipeline pressure loss would occur and, thereby, reduce the work efficiency and load capacity of the deep-sea hydraulic system. Besides, by comparing Figure 14 and Figure 15, it can be known that the pipeline outlet pressure has a greater impact on pressure loss. For the inflow pipeline which connects the hydraulic pump and the hydraulic actuator, the pipeline outlet pressure is the sum of the load pressure and the ambient pressure in the deep sea. As for the return pipeline that connects the hydraulic actuator and the oil tank, pipeline outlet pressure equals ambient pressure. Then, it can be known that the extremely high ambient pressure in the deep-sea environment has a significant impact on the working efficiency and load capacity of the deep-sea hydraulic system.
The percentage curves in Figure 14 and Figure 15 also indicate that, when the viscosity of the hydraulic oil changes minimally, the pressure loss calculated by Equation (23) is less different from the classic Poiseuille’s law.
The temperature also has a great influence on the viscosity of hydraulic oil. However, because the hydraulic system has a large number of heat sources, and the low-temperature seawater with a large specific heat capacity is sufficient, the thermal field distribution of deep-sea hydraulic pipelines becomes very complicated. Therefore, the influence of temperature is not considered in the theoretical analysis and CFD simulation in this paper, and it is assumed that the hydraulic system is operating at a design condition of 40 °C. The following only briefly discusses the influence of temperature.
The viscosity of hydraulic oil increases greatly in low-temperature conditions. Due to sufficient seawater with a large specific heat capacity and low temperature, the heat dissipation conditions of deep-sea hydraulic systems are much better than those in land conditions. When the deep-sea hydraulic system has not been started or is in a standby state for a long time, the temperature of the hydraulic oil will tend to be the same as the seawater temperature, which is only about 2–4 °C. As a result, the viscosity of hydraulic oil greatly increases. This can significantly increase the pipeline pressure loss when the deep-sea hydraulic system is cold started or restored from a long-term standby state compared with the normal working state. Therefore, the work efficiency and the load capacity of deep-sea hydraulic systems are greatly reduced.

5.4. Relationship between the Novel Equation and the Classic Poiseuille’s Law

According to the discussions in the previous subsections, when the pressure loss or the viscosity change is minimal, the pressure loss calculated by Equation (23) is almost the same as the result of the classic Poiseuille’s law. In this subsection, the relationship between them will be discussed from the perspective of theoretical analysis.
The transformation of Equation (23), which is the proposed novel equation for pipeline pressure loss, is as follows
α Δ p = ln ( 1 128 α η 0 L Q e α p o π D 4 ) .
Taking the natural logarithm of Equation (24), we have
e α Δ p = 1 128 α η 0 L Q e α p o π D 4 .
According to Equations (2) and (4), the term e α Δ p indicates the degree of change in hydraulic oil viscosity. When the pressure loss Δ p is tiny or the viscosity changes minimally, we have
e α Δ p 1 α Δ p .
Substituting Equation (26) into Equation (25), we have
1 α Δ p 1 128 α η 0 L Q e α p o π D 4 ,
or
Δ p 128 η 0 L Q e α p o π D 4 .
The expression of Equation (28) is exactly the same as that of Equation (3), which is obtained based on classical Poiseuille’s law.
In addition, the proposed novel equation can be expanded by the Taylor series. The Taylor series of the logarithmic function is
ln ( 1 + x ) = x x 2 2 + x 3 3 .
Then,
Δ p = 1 α ln ( 1 128 α η 0 L Q e α p o π D 4 ) = 1 α ln 1 + 128 α η 0 L Q e α p o π D 4 = 1 α 128 α η 0 L Q e α p o π D 4 128 α η 0 L Q e α p o π D 4 2 2 + 128 α η 0 L Q e α p o π D 4 3 3 = 128 η 0 L Q e α p o π D 4 + 1 α 128 α η 0 L Q e α p o π D 4 2 2 1 α 128 α η 0 L Q e α p o π D 4 3 3 + 1 α × .
The first term in Equation (30) is a one-time term, which is also exactly the same as Equation (3).
The derivation of the above two methods shows that the proposed novel equation for the pipeline pressure loss is equivalent to the classic Poiseuille’s law when the pipeline pressure loss or the viscosity change is minimal. Then, we can conclude that the classic Poiseuille’s law is a simplified expression of the novel proposed equation for pipeline pressure loss in the case of small viscosity changes. In other words, the novel pressure loss equation is an extension of the classic Poiseuille’s law.

6. Conclusions

(1.) Based on laminar flow theory and the viscosity-pressure characteristics of hydraulic oil, a novel equation for pipeline pressure loss is proposed, in which the viscosity change when flowing through the pipeline is taken into account. The pipeline pressure loss is related to the viscosity-pressure index and the initial viscosity of the hydraulic oil, the length and inner diameter of the pipeline, the flow rate, and the outlet pressure.
(2.) Theoretical analysis shows that the proposed novel equation for pipeline pressure loss is equivalent to the classic Poiseuille’s law when the pipeline pressure loss or the viscosity change is minimal, which means the novel pressure loss equation is an extension of the classic Poiseuille’s law.
(3.) The larger the viscosity-pressure index or the initial viscosity of the hydraulic oil, the length of the pipeline, the flow rate, or the outlet pressure is, the greater the deviation between the pressure loss calculated by the proposed novel equation and the classic Poiseuille’s law. In addition, the smaller the inner diameter of the pipeline is, the greater the deviation is. The difference in pressure loss calculated by the two methods reaches 58.31% when the viscosity-pressure index of the hydraulic oil is 3 × 10 8 Pa 1 .
(4.) A CFD model of a pipeline in the deep-sea hydraulic system is established, and CFD simulations are conducted. The maximum error between the pressure loss calculated by the novel equation and the one calculated by the CFD simulation is only −1.80%. The fluid velocity distribution obtained by CFD simulation is consistent with the theoretical analysis. All the CFD results justify the correctness of the proposed novel equation for pipeline pressure loss.

Author Contributions

The paper emerged from the long-term cooperation of the authors. Conceptualization, L.L.; Formal analysis, J.-B.W.; Funding acquisition, L.L.; Methodology, J.-B.W.; Project administration, L.L.; Software, J.-B.W.; Validation, J.-B.W.; Writing—original draft, J.-B.W.; Writing—review and editing, L.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Hunan Provincial Science and Technology Department (Grant No. 2019SK2271, 2020GK1020).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available in a publicly accessible repository.

Acknowledgments

The authors would like to thank the editors and anonymous reviewers for their careful work and thoughtful suggestions that have helped improve this paper substantially.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AUVsAutonomous underwater vehicles
CFDComputational Fluid Dynamics
HOVsHuman occupied vehicles
RMSRoot Mean Square
ROVsRemotely operated Vehicles

References

  1. Chu, Z.; Chen, Y.; Zhu, D.; Zhang, M. Observer-based fault detection for magnetic coupling underwater thrusters with applications in jiaolong HOV. Ocean Eng. 2020, 210, 107570. [Google Scholar] [CrossRef]
  2. Capocci, R.; Dooly, G.; Omerdić, E.; Coleman, J.; Newe, T.; Toal, D. Inspection-Class Remotely Operated Vehicles—A Review. J. Mar. Sci. Eng. 2017, 5, 13. [Google Scholar] [CrossRef]
  3. Dai, Y.; Xue, C.; Su, Q. An Integrated Dynamic Model and Optimized Fuzzy Controller for Path Tracking of Deep-Sea Mining Vehicle. J. Mar. Sci. Eng. 2021, 9, 249. [Google Scholar] [CrossRef]
  4. Dai, Y.; Ma, F.; Zhu, X.; Liu, H.; Huang, Z.; Xie, Y. Mechanical Tests and Numerical Simulations for Mining Seafloor Massive Sulfides. J. Mar. Sci. Eng. 2019, 7, 252. [Google Scholar] [CrossRef] [Green Version]
  5. Li, D.; Du, L. AUV Trajectory Tracking Models and Control Strategies: A Review. J. Mar. Sci. Eng. 2021, 9, 1020. [Google Scholar] [CrossRef]
  6. Vu, M.T.; Choi, H.S.; Kim, J.Y.; Tran, N.H. A study on an underwater tracked vehicle with a ladder trencher. Ocean Eng. 2016, 127, 90–102. [Google Scholar] [CrossRef]
  7. Huo, X.; Wang, X.; Ge, T. Impulse control method for hydraulic propulsion system used in 3500m work-class ROV. Appl. Ocean Res. 2016, 60, 75–83. [Google Scholar] [CrossRef]
  8. Liu, Y.; Wu, D.; Li, D.; Deng, Y. Applications and Research Progress of Hydraulic Technology in Deep Sea. J. Mech. Eng. 2018, 54, 14–23. [Google Scholar] [CrossRef]
  9. Wang, F.; Chen, Y. Dynamic characteristics of pressure compensator in underwater hydraulic system. IEEE ASME Trans. Mechatron. 2014, 19, 777–787. [Google Scholar] [CrossRef]
  10. Li, S.; Du, X.; Zhang, L.; Chen, K.; Wang, S. Study on dynamic characteristics of underwater pressure compensator considering nonlinearity. Mech. Sci. 2020, 11, 183–192. [Google Scholar] [CrossRef]
  11. Wu, J.B.; Li, L.; Wei, W. Research on dynamic characteristics of pressure compensator for deep-sea hydraulic system. Proc. Inst. Mechan. Eng. Part M J. Eng. Maritime Environ. 2021. Available online: https://journals.sagepub.com/doi/abs/10.1177/14750902211033265 (accessed on 1 October 2021). [CrossRef]
  12. Cao, X.p.; Ye, M.; Deng, B.; Zhang, C.; Yu, Z. LVP modeling and dynamic characteristics prediction of a hydraulic power unit in deep-sea. China Ocean Eng. 2013, 27, 17–32. [Google Scholar] [CrossRef]
  13. Tian, Q.; Zhang, Q.; Huo, L.; Zhang, Y.; Du, L.; Sun, J. Modeling and experimental validation of the hydraulic power source for deep-sea manipulator. In Proceedings of the OCEANS 2017, Anchorage, AK, USA, 18–21 September 2017; pp. 1–6. [Google Scholar]
  14. Tian, Q.; Zhang, Q.; Chen, Y.; Huo, L.; Li, S.; Wang, C.; Bai, Y.; Du, L. Influence of Ambient Pressure on Performance of a Deep-sea Hydraulic Manipulator. In Proceedings of the OCEANS 2019, Marseille, France, 17–20 June 2019. [Google Scholar] [CrossRef]
  15. Xin, Y.; Yao, Y.; Liu, M.; Zhang, W. Establishment of valve control mathematical model of hydraulic differential cylinder. In Proceedings of the 2nd International Conference on Advanced Engineering Materials and Technology, AEMT 2012, Zhuhai, China, 6–8 July 2012; Trans Tech Publications: Stafa-Zurich, Switzerland, 2012; pp. 1124–1131. [Google Scholar] [CrossRef]
  16. Litong, C. Hydraulic Control Systems; Tsinghua University Press: Beijing, China, 2014. (In Chinese) [Google Scholar]
  17. Savić, V.; Knežević, D.; Lovrec, D.; Jocanović, M.; Karanović, V. Determination of pressure losses in hydraulic pipeline systems by considering temperature and pressure. Stroj. Vestnik J. Mech. Eng. 2009, 55, 237–243. [Google Scholar]
  18. Ganji, D.; Ashory Nezhad, H.; Hasanpour, A. Effect of variable viscosity and viscous dissipation on the Hagen-Poiseuille flow and entropy generation. Numer. Methods Part. Differ. Equ. 2011, 27, 529–540. [Google Scholar] [CrossRef]
  19. Housiadas, K.D.; Georgiou, G.C. New analytical solutions for weakly compressible Newtonian Poiseuille flows with pressure-dependent viscosity. Int. J. Eng. Sci. 2016, 107, 13–27. [Google Scholar] [CrossRef]
  20. Makinde, O.; Iskander, T.; Mabood, F.; Khan, W.; Tshehla, M. MHD Couette-Poiseuille flow of variable viscosity nanofluids in a rotating permeable channel with Hall effects. J. Mol. Liq. 2016, 221, 778–787. [Google Scholar] [CrossRef]
  21. Lee, Y.M.; Kim, J.H.; Lee, J.H. Direct numerical simulation of a turbulent Couette-Poiseuille flow with a rod-roughened wall. Phys. Fluids 2018, 30, 105101. [Google Scholar] [CrossRef]
  22. Qiu, H. Research on the Compensation of Viscous Pressure Characteristics and Variable Gain Control of Deepwater Hydraulic Manipulator. Master’s Thesis, Zhejiang University, Hangzhou China, 2018. (In Chinese). [Google Scholar]
  23. Luo, G.; Yao, Z.; Shen, H. Mechanism of pressure oscillation in Taylor-Couette-Poiseuille flow with abruptly contracting and expanding annular gap. Phys. Fluids 2019, 31, 075105. [Google Scholar] [CrossRef]
  24. Li, L.; Wu, J.B. Deformation and leakage mechanisms at hydraulic clearance fit in deep-sea extreme environment. Phys. Fluids 2020, 32, 067115. [Google Scholar] [CrossRef]
  25. Shankar, B.M.; Shivakumara, I.S. Stability of porous-Poiseuille flow with uniform vertical throughflow: High accurate solution. Phys. Fluids 2020, 32, 044101. [Google Scholar] [CrossRef] [Green Version]
  26. Barus, C. Isothermals, isopiestics, and isometrics relative to viscosity. Am. J. Sci. 1893, 45, 87–96. [Google Scholar] [CrossRef]
  27. Jiwei, W. Hydraulic Transmission, 2nd ed.; China Machine Press: Beijing, China, 2007. (In Chinese) [Google Scholar]
  28. Jinchun, S. Practical Handbook of Hydraulic Technology; China Electric Power Press: Beijing, China, 2010. (In Chinese) [Google Scholar]
  29. En, M.; Sumin, L. Hydraulic and Fluid Power Transmission; Tsinghua University Press: Beijing, China, 2015. (In Chinese) [Google Scholar]
  30. Vacca, A.; Franzoni, G. Hydraulic Fluid Power: Fundamentals, Applications, and Circuit Design; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2021. [Google Scholar]
  31. Martins, N.M.C.; Soares, A.K.; Ramos, H.M.; Covas, D.I.C. CFD modeling of transient flow in pressurized pipes. Comput. Fluids 2016, 126, 129–140. [Google Scholar] [CrossRef]
  32. Martins, N.M.C.; Brunone, B.; Meniconi, S.; Ramos, H.M.; Covas, D. Efficient CFD model for transient laminar flow modeling: Pressure wave propagation and velocity profile changes. ASME J. Fluids Eng. 2018, 140, 011102. [Google Scholar] [CrossRef]
Figure 1. Power transmission block diagram of deep-sea operation equipment.
Figure 1. Power transmission block diagram of deep-sea operation equipment.
Jmse 09 01142 g001
Figure 2. Analysis model of laminar flow in a circular pipeline.
Figure 2. Analysis model of laminar flow in a circular pipeline.
Jmse 09 01142 g002
Figure 3. Pipeline cross-section and micro-ring.
Figure 3. Pipeline cross-section and micro-ring.
Jmse 09 01142 g003
Figure 4. CFD model of the fluid domain in the pipeline.
Figure 4. CFD model of the fluid domain in the pipeline.
Jmse 09 01142 g004
Figure 5. Pressure loss, average flow velocity, and calculation time versus total mesh elements.
Figure 5. Pressure loss, average flow velocity, and calculation time versus total mesh elements.
Jmse 09 01142 g005
Figure 6. Dynamic viscosity.
Figure 6. Dynamic viscosity.
Jmse 09 01142 g006
Figure 7. Pressure distribution with variable viscosity.
Figure 7. Pressure distribution with variable viscosity.
Jmse 09 01142 g007
Figure 8. Pressure distribution with constant viscosity.
Figure 8. Pressure distribution with constant viscosity.
Jmse 09 01142 g008
Figure 9. Velocity distributions at different positions.
Figure 9. Velocity distributions at different positions.
Jmse 09 01142 g009
Figure 10. Pressure loss versus viscosity-pressure index.
Figure 10. Pressure loss versus viscosity-pressure index.
Jmse 09 01142 g010
Figure 11. Pressure loss versus initial viscosity.
Figure 11. Pressure loss versus initial viscosity.
Jmse 09 01142 g011
Figure 12. Pressure loss versus pipeline length.
Figure 12. Pressure loss versus pipeline length.
Jmse 09 01142 g012
Figure 13. Pressure loss versus pipeline inner diameter.
Figure 13. Pressure loss versus pipeline inner diameter.
Jmse 09 01142 g013
Figure 14. Pressure loss versus outlet pressure of pipeline.
Figure 14. Pressure loss versus outlet pressure of pipeline.
Jmse 09 01142 g014
Figure 15. Pressure loss versus flow rate.
Figure 15. Pressure loss versus flow rate.
Jmse 09 01142 g015
Table 1. Sets of mesh elements.
Table 1. Sets of mesh elements.
Case123456
Number of mesh elements300,0001,200,0002,025,0004,800,0005,808,0007,500,000
Table 2. Pressure loss.
Table 2. Pressure loss.
Simulated ValueTheoretical ValueError
Pressure loss with variable viscosity8.53858 MPa8.50808 MPa−0.36%
Pressure loss with constant viscosity7.78795 MPa7.75926 Mpa−0.37%
Increment percentage9.64%9.65%-
Table 3. Variable inputs and pressure losses: part A.
Table 3. Variable inputs and pressure losses: part A.
ParameterSymbolUnitCase 1Case 2Case 3Case 4Case 5Case 6Case 7
Viscosity-pressure index α Pa 1 2.2 × 10 8 1.5 × 10 8 2.5 × 10 8 2.2 × 10 8 2.2 × 10 8 2.2 × 10 8 2.2 × 10 8
Initial kinematic viscosity-cSt22222232462222
Pipeline lengthLm222221.52.5
Pipeline inner diameterDmm4444444
Outlet pressure of the pipeline p o  MPa125125125125125125125
Flow rateQL/min5555555
Reynolds number--77.08184.9052.9852.9932.8677.0877.08
Simulated pressure loss Δ p  MPa8.538583.3420313.2961013.0044020.107006.2610510.93880
Theoretical pressure loss Δ p  MPa8.508083.3156413.2648812.9733220.067976.2271710.90953
Error--−0.36%−0.79%−0.23%−0.24%−0.19%−0.54%−0.27%
Table 4. Variable inputs and pressure losses: part B.
Table 4. Variable inputs and pressure losses: part B.
ParameterSymbolUnitCase 1Case 2Case 3Case 4Case 5Case 6
Viscosity-pressure index α Pa 1 2.2 × 10 8 1.5 × 10 8 2.5 × 10 8 2.2 × 10 8 2.2 × 10 8 2.2 × 10 8
Initial kinematic viscosity-cSt222222222222
Pipeline lengthLm222222
Pipeline inner diameterDmm684444
Outlet pressure of the pipeline p o  MPa1251258545125125
Flow rateQL/min555546
Reynolds number--51.3938.54185.83448.0261.6692.49
Simulated pressure loss Δ p  MPa1.562940.488553.364621.379816.6943010.46120
Theoretical pressure loss Δ p  MPa1.559130.487563.338031.354946.6742910.41898
Error--−0.24%−0.20%−0.79%−1.80%−0.30%−0.40%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, J.-B.; Li, L. Pipeline Pressure Loss in Deep-Sea Hydraulic Systems Considering Pressure-Dependent Viscosity Change of Hydraulic Oil. J. Mar. Sci. Eng. 2021, 9, 1142. https://doi.org/10.3390/jmse9101142

AMA Style

Wu J-B, Li L. Pipeline Pressure Loss in Deep-Sea Hydraulic Systems Considering Pressure-Dependent Viscosity Change of Hydraulic Oil. Journal of Marine Science and Engineering. 2021; 9(10):1142. https://doi.org/10.3390/jmse9101142

Chicago/Turabian Style

Wu, Jia-Bin, and Li Li. 2021. "Pipeline Pressure Loss in Deep-Sea Hydraulic Systems Considering Pressure-Dependent Viscosity Change of Hydraulic Oil" Journal of Marine Science and Engineering 9, no. 10: 1142. https://doi.org/10.3390/jmse9101142

APA Style

Wu, J. -B., & Li, L. (2021). Pipeline Pressure Loss in Deep-Sea Hydraulic Systems Considering Pressure-Dependent Viscosity Change of Hydraulic Oil. Journal of Marine Science and Engineering, 9(10), 1142. https://doi.org/10.3390/jmse9101142

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