^{1}

^{2}

This article focuses on how to control the pressure in the high-pressure tubing. First, an overall analysis of the system is performed to determine the pressure balance conditions of the high-pressure fuel pipe; that is, the mass of fuel flowing out of the high-pressure fuel pipe is equal to the mass of fuel flowing in. This article uses Excel to fit and organize the data in all data processing. When solving differential equations of pressure and density, the fourth-order Runge-Kutta formula is used, and the numerical solution is obtained by MATLAB2017a.

Fuel entry and ejection of high-pressure fuel pipes are the basis of the work of many fuel engines.

The inner cavity length of a certain type of high-pressure oil pipe is 500 mm; the inner diameter is 10 mm; the diameter of the small hole at the oil supply inlet A is 1.4 mm; the length of the oil supply time is controlled by the check valve switch, and the check valve needs to be opened every time. Turn off for 10 ms. The fuel injector works 10 times per second, and the fuel injection time is 2.4 ms during each operation. The rate of fuel injection from the fuel injection nozzle B when the fuel injector is working is shown in

This article raises the following two questions:

1) The pressure provided by the high-pressure oil pump at the inlet A is constant at 160 MPa, and the initial pressure in the high-pressure oil pipe is 100 MPa. If the pressure in the high-pressure oil pipe is to be stabilized at about 100 MPa as much as possible, how should the opening time of the check valve be adjusted?

2) If the pressure in the high-pressure oil pipe is to be increased from 100 MPa to 150 MPa and stabilized at 150 MPa after an adjustment process of about 2 s, how should the duration of the check valve be adjusted?

The process of fuel pumping into the high-pressure fuel pipe is constantly changing. It is undoubtedly difficult and complicated to try to achieve dynamic balance and analyze it every moment, which is not conducive to the establishment of mathematical models. In order to simplify the model, this dynamic process can be divided into a superposition of many small segments. In order to maintain the pressure of the high-pressure oil pipe stable, the quality of the oil pumped by the high-pressure oil pump should be equal to the quality of the oil sprayed by the fuel injection nozzle. By analyzing and calculating the amount of oil sprayed by the fuel injection nozzle within one working cycle (100 ms), the amount of oil that the high-pressure oil pump needs to pump in the same period of time can be obtained, and the time when the check valve opens is derived.

1) The pressure of the high-pressure oil pipe is only restricted by the conditions specified in this article, and we will not discuss other influencing factors.

2) The model stipulates that each coefficient is for a general general situation, and we will not discuss other situations.

3) In the analysis of mathematics and physical equations, we set the measurement error to be within the scope of the problem.

4) Ignore fuel leakage, cavitation and temperature changes.

5) At the same moment, the pressure in each point of the high-pressure oil pipe is equal.

6) When calculating the pressure balance, ignore the slight pressure changes in the high-pressure tubing.

7) Regardless of the propagation speed of pressure waves.

Appendix shows the explanation of symbols used in this article.

According to the change of fuel and the change of density is proportional to the proportional coefficient of E ρ , it can be seen that the pressure and density of fuel meet the differential equation:

d P d ρ = E ρ (5-1)

Processing the given data, the fitted curve and correlation are as follows:

Fitting the equation of elastic modulus E with respect to pressure P:

E = 0.0001 P 3 − 0.0011 P 2 + 5.47449 P + 153.9 (5-2)

Take (5-2) into (5-1) to get the differential equation that pressure P and density ρ satisfy:

d P d ρ = 0.0001 P 3 − 0.0011 P 2 + 5.47449 P + 153.9 ρ

It is known that ρ = 0.85 mg / mm 3 when P = 100 MPa . To solve the above equations, the fourth-order Runge-Kutta formula (RK) method is used. First discretize the equation, generally using the difference quotient to approximate the derivative [

{ d y d x = f ( x , y ) a ≤ x ≤ b y ( a ) = y 0 (5-3)

Let h n = x n + 1 − x n be called the step size from x n to x n + 1 . Unless otherwise specified, we always take the step size as a constant h. If the forward difference quotient y ( x n + 1 ) − y ( n ) h is substituted for y ′ ( x n ) into the differential equation, then:

y ( x n + 1 ) − y ( n ) h ≈ f ( x n , y ( x n ) ) ( n = 0 , 1 , ⋯ , N − 1 )

Simplify:

y ( x n + 1 ) ≈ y ( x n ) + h f ( x n , y ( x n ) )

If the approximate value y n of y ( x n ) is substituted into the right end of the above formula, the result is taken as the approximate value of y ( x n + 1 ) , which is recorded as y n + 1 :

{ y n + 1 = y n + h f ( x n , y n ) ( n = 0 , 1 , ⋯ , N − 1 ) y 0 = y ( a )

According to the median theorem of differential equations:

y ( x n + 1 ) − y ( x n ) h = y ′ ( x n + θ h ) , 0 < θ < 1

Note that the equation y ` = f ( x , y ) can be obtained:

y ( x n + 1 ) = y ( x n ) + h f ( x n + θ h , y ( x n + θ h ) )

Let K ¯ = f ( x n + θ h , y ( x n + θ h ) ) be called the average slope in the interval [ x n , x n + 1 ] . Taking a few more points in the interval [ x n , x n + 1 ] and using their slope-weighted average as K ¯ , you can construct a calculation formula with better accuracy, the fourth-order Runge-Kutta formula (RK) takes 4 points to construct the following formula [

{ y n + 1 = y n + h ( λ 1 k 1 + λ 2 k 2 + λ 3 k 3 + λ 4 k 4 ) k 1 = f ( x n , y n ) k 2 = f ( x n + α 1 h , y n + β 1 h k 1 ) k 3 = f ( x n + α 2 h , y n + β 2 h k 1 + β 3 h k 2 ) k 4 = f ( x n + α 3 h , y n + β 4 h k 1 + β 5 h k 2 + β 6 h k 3 )

Among them, there are 13 undetermined coefficients λ i , α j , β m , which are derived and calculated. Eleven equations with local errors y ( x n + 1 ) − y ( x n ) = O ( h 5 ) are obtained. Taking simple λ i , α j , β m that satisfy these equations [

{ y n + 1 = h 6 ( k 1 + 2 k 2 + 2 k 3 + k 4 ) k 1 = f ( x n , y n ) k 2 = f ( x n + h 2 , y n + h k 1 2 ) k 3 = f ( x n + h 2 , y n + h k 2 2 ) k 4 = f ( x n + h , y n + h k 3 )

Use MATLAB to calculate the numerical solution of the equation. After writing the data into an Excel file, fit the data to obtain a function image with a correlation coefficient of 1:

When P = 160 MPa , the fuel density is 0.8711 mg/mm^{3}; when P = 150 MPa , ρ = 0.8679 mg / mm 3 . In order to maintain the stability of the pressure in the high-pressure fuel pipe, the fuel mass flowing into the high-pressure fuel pipe at port A should satisfy m i n = m o u t . The calculation of the mass of fuel flowing into port A and the mass of fuel ejected from the fuel injector during a working cycle of the fuel injector is as follows:

{ m i n = Q i n × ρ 160 × t t o t a l Q = C A 2 Δ P ρ m o u t = 44 × ρ 100

Simultaneous calculation of the above equations:

{ m o u t = m i n = 37.4 mg t t o t a l = 2.8 ms

Within 100 ms, the number of inflows is 100 t + 10 , t is the time for each opening, and t t o t a l is the total time for the A port one-way valve to open during one working cycle of the nozzle port, so:

100 t + 10 × t = t t o t a l

Therefore, the duration of each valve opening is t = 0.288 ms .

In order to make P in the high-pressure tubing from 100 MPa to 150 MPa in 2 s, the pressure task that needs to be increased can be equally divided into each working cycle of the nozzle port in 2 s, so the opening time of the check valve is changed every 100 ms to The pressure of the high-pressure tubing rises by 2.5 MPa within 100 ms. Let P n be the high pressure tubing pressure.

0 - 100 ms:

{ m i n = Q × ρ × t t o t a l ( m i n − m o u t ) / V = ρ 102.5 − ρ 100 m o u t = 44 × ρ 100 m i n = ρ 102.5 − ρ 100 × V + m o u t

where V is the volume of the high-pressure tubing, V = 500 × 5 2 × π = 12500 π mm 3 , P 1 = 100 MPa , simultaneous calculation of the above equation:

t t o t a l _ 1 = ρ 102.5 − ρ 100 × V + 44 × ρ 100 0.85 × 0.49 π 2 × ( 160 − P 1 ) × 0.8711

Therefore, the check valve opening time is:

t 1 = 10 t t o t a l _ 1 100 − t t o t a l _ 1

100 - 200 ms: Similarly, we can get:

t t o t a l _ 2 = ρ 105 − ρ 102.5 × V + 44 × ρ 102.5 0.85 × 0.49 π 2 × ( 160 − P 2 ) × 0.8711

P 2 = 102.5 MPa

t 2 = 10 t t o t a l _ 2 100 − t t o t a l _ 2

200 - 300 ms:

…

1900 - 2000 ms:

t t o t a l _ 20 = ρ 150 − ρ 147.5 × V + 44 × ρ 147.5 0.85 × 0.49 π 2 × ( 160 − P 20 ) × 0.8711

t 20 = 10 t t o t a l _ 20 100 − t t o t a l _ 20

With this iteration, MATLAB can be used to calculate the adjustment time of the check valve every 100 ms within the adjustment time of 2 s. In order to stabilize the adjusted pressure at 150 MPa, the mass flowing out of the nozzle is

Total adjustment time 2 s | |||
---|---|---|---|

Work cycle (100 ms) | Check valve opening time (ms) | Work cycle (100 ms) | Check valve opening time (ms) |

1 | 0.599997543 | 11 | 0.768145842 |

2 | 0.611289684 | 12 | 0.796051254 |

3 | 0.623415249 | 13 | 0.827681804 |

4 | 0.636477259 | 14 | 0.863911514 |

5 | 0.650597192 | 15 | 0.905926936 |

6 | 0.665919426 | 16 | 0.955386087 |

7 | 0.68261705 | 17 | 1.014688202 |

8 | 0.700899583 | 18 | 1.087458477 |

9 | 0.721023391 | 19 | 1.179483494 |

10 | 0.743305952 | 20 | 1.300690386 |

equal to the mass flowing in port A within a working cycle, and the opening time of the one-way valve at this time is t, then within 100 ms:

{ m o u t = m i n m o u t = 44 × ρ 150 m i n = C A 2 Δ P ρ 160 × ρ 160 × t t o t a l

The calculation is:

t t o t a l = 6.992 ms , t = 0.752 ms .

The adjustment scheme of the check valve thus obtained is:

In order to stabilize the pressure of the high-pressure tubing at 100 MPa, the opening time of the check valve every 100 ms should be t = 0.288 ms.

In order to increase the pressure of the high-pressure tubing to 150 MPa, when the adjustment time is 2 s, the opening time of the check valve changed every 100 ms is shown in

In order to stabilize the pressure of the high-pressure tubing at 150 MPa, the opening time of the check valve every 100 ms should be t = 0.752 ms.

This model achieves reasonable assumptions and is close to actual problems. Compared with other complex and computationally intensive models, it has the advantages of simple and easy to understand and simple programming.

The authors declare no conflicts of interest regarding the publication of this paper.

Hu, S.L. and Quan, Y. (2020) Pressure Control of High Pressure Tubing. Open Access Library Journal, 7: e6351. https://doi.org/10.4236/oalib.1106351