Pressure Control of High Pressure Tubing ()

Shilin Hu^{1}, Yu Quan^{2*}

^{1}College of Science, Yanbian University, Yanji, China.

^{2}College of Economics and Management, Yanbian University, Yanji, China.

**DOI: **10.4236/oalib.1106351
PDF HTML XML
102
Downloads
291
Views
Citations

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.

Share and Cite:

Hu, S. and Quan, Y. (2020) Pressure Control of High Pressure Tubing. *Open Access Library Journal*, **7**, 1-8. doi: 10.4236/oalib.1106351.

1. Question Restatement

Fuel entry and ejection of high-pressure fuel pipes are the basis of the work of many fuel engines. Figure 1 shows the working principle of a high-pressure fuel system. Fuel enters the high-pressure fuel pipe from A through the high-pressure fuel pump and is then ejected from nozzle B. The intermittent working process of fuel entry and ejection will cause the pressure in the high-pressure fuel pipe to change, which will cause the deviation of the amount of fuel injected, thus affecting the working efficiency of the engine.

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 Figure 2.

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?

2. Question Analysis

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.

3. Assumptions and Explanations about Model Building

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.

Figure 1. Schematic diagram of high-pressure oil pipe.

Figure 2. Schematic diagram of fuel injection rate.

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.

4. Modeling and Solving of Questions

4.1. Determination of Fuel Pressure and Density

According to the change of fuel and the change of density is proportional to the proportional coefficient of $\frac{E}{\rho}$ , it can be seen that the pressure and density of fuel meet the differential equation:

$\frac{\text{d}P}{\text{d}\rho}=\frac{E}{\rho}$ (5-1)

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

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

$E=0.0001{P}^{3}-0.0011{P}^{2}+5.47449P+153.9$ (5-2)

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

$\frac{\text{d}P}{\text{d}\rho}=\frac{0.0001{P}^{3}-0.0011{P}^{2}+5.47449P+153.9}{\rho}$

It is known that $\rho =0.85\text{\hspace{0.17em}}\text{mg}/{\text{mm}}^{\text{3}}$ when $P=100\text{\hspace{0.17em}}\text{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 [1]. Let the differential equation be:

$\{\begin{array}{l}\frac{\text{d}y}{\text{d}x}=f\left(x,y\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}a\le x\le b\\ y\left(a\right)={y}_{0}\end{array}$ (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 $\frac{y\left({x}_{n+1}\right)-y\left(n\right)}{h}$ is substituted for ${y}^{\prime}\left({x}_{n}\right)$ into the differential equation, then:

$\frac{y\left({x}_{n+1}\right)-y\left(n\right)}{h}\approx f\left({x}_{n},y\left({x}_{n}\right)\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(n=0,1,\cdots ,N-1\right)$

Figure 3. The relationship between E and P.

Simplify:

$y\left({x}_{n+1}\right)\approx y\left({x}_{n}\right)+hf\left({x}_{n},y\left({x}_{n}\right)\right)$

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

$\{\begin{array}{l}{y}_{n+1}={y}_{n}+hf\left({x}_{n},{y}_{n}\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(n=0,1,\cdots ,N-1\right)\\ {y}_{0}=y(\; a\; )\end{array}$

According to the median theorem of differential equations:

$\frac{y\left({x}_{n+1}\right)-y\left({x}_{n}\right)}{h}={y}^{\prime}\left({x}_{n}+\theta h\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}0<\theta <1$

Note that the equation ${y}^{`}=f\left(x,y\right)$ can be obtained:

$y\left({x}_{n+1}\right)=y\left({x}_{n}\right)+hf\left({x}_{n}+\theta h,y\left({x}_{n}+\theta h\right)\right)$

Let $\stackrel{\xaf}{K}=f\left({x}_{n}+\theta h,y\left({x}_{n}+\theta h\right)\right)$ be called the average slope in the interval $\left[{x}_{n},{x}_{n+1}\right]$ . Taking a few more points in the interval $\left[{x}_{n},{x}_{n+1}\right]$ and using their slope-weighted average as $\stackrel{\xaf}{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 [2]:

$\{\begin{array}{l}{y}_{n+1}={y}_{n}+h\left({\lambda}_{1}{k}_{1}+{\lambda}_{2}{k}_{2}+{\lambda}_{3}{k}_{3}+{\lambda}_{4}{k}_{4}\right)\\ {k}_{1}=f\left({x}_{n},{y}_{n}\right)\\ {k}_{2}=f\left({x}_{n}+{\alpha}_{1}h,{y}_{n}+{\beta}_{1}h{k}_{1}\right)\\ {k}_{3}=f\left({x}_{n}+{\alpha}_{2}h,{y}_{n}+{\beta}_{2}h{k}_{1}+{\beta}_{3}h{k}_{2}\right)\\ {k}_{4}=f\left({x}_{n}+{\alpha}_{3}h,{y}_{n}+{\beta}_{4}h{k}_{1}+{\beta}_{5}h{k}_{2}+{\beta}_{6}h{k}_{3}\right)\end{array}$

Among them, there are 13 undetermined coefficients ${\lambda}_{i},{\alpha}_{j},{\beta}_{m}$ , which are derived and calculated. Eleven equations with local errors $y\left({x}_{n+1}\right)-y\left({x}_{n}\right)=O\left({h}^{5}\right)$ are obtained. Taking simple ${\lambda}_{i},{\alpha}_{j},{\beta}_{m}$ that satisfy these equations [3], we can get:

$\{\begin{array}{l}{y}_{n+1}=\frac{h}{6}\left({k}_{1}+2{k}_{2}+2{k}_{3}+{k}_{4}\right)\\ {k}_{1}=f\left({x}_{n},{y}_{n}\right)\\ {k}_{2}=f\left({x}_{n}+\frac{h}{2},{y}_{n}+\frac{h{k}_{1}}{2}\right)\\ {k}_{3}=f\left({x}_{n}+\frac{h}{2},{y}_{n}+\frac{h{k}_{2}}{2}\right)\\ {k}_{4}=f\left({x}_{n}+h,{y}_{n}+h{k}_{3}\right)\end{array}$

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: Figure 4.

4.2. Calculation of Question (1)

When
$P=160\text{\hspace{0.17em}}\text{MPa}$ , the fuel density is 0.8711 mg/mm^{3}; when
$P=150\text{\hspace{0.17em}}\text{MPa}$ ,
$\rho =0.8679\text{\hspace{0.17em}}\text{mg}/{\text{mm}}^{\text{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}_{in}={m}_{out}$ . 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:

$\{\begin{array}{l}{m}_{in}={Q}_{in}\times {\rho}_{160}\times {t}_{total}\\ Q=CA\sqrt{\frac{2\Delta P}{\rho}}\\ {m}_{out}=44\times {\rho}_{100}\end{array}$

Simultaneous calculation of the above equations:

$\{\begin{array}{l}{m}_{out}={m}_{in}=37.4\text{\hspace{0.17em}}\text{mg}\\ {t}_{total}=2.8\text{\hspace{0.17em}}\text{ms}\end{array}$

Within 100 ms, the number of inflows is $\frac{100}{t+10}$ , t is the time for each opening, and ${t}_{total}$ is the total time for the A port one-way valve to open during one working cycle of the nozzle port, so:

Figure 4. Relationship between ρ and P.

$\frac{100}{t+10}\times t={t}_{total}$

Therefore, the duration of each valve opening is $t=0.288\text{\hspace{0.17em}}\text{ms}$ .

4.3. Calculation of Question (2)

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:

$\{\begin{array}{l}{m}_{in}=Q\times \rho \times {t}_{total}\\ \left({m}_{in}-{m}_{out}\right)/V={\rho}_{102.5}-{\rho}_{100}\\ {m}_{out}=44\times {\rho}_{100}\\ {m}_{in}={\rho}_{102.5}-{\rho}_{100}\times V+{m}_{out}\end{array}$

where V is the volume of the high-pressure tubing, $V=500\times {5}^{2}\times \pi =12500\pi \text{\hspace{0.17em}}{\text{mm}}^{3}$ , ${P}_{1}=100\text{\hspace{0.17em}}\text{MPa}$ , simultaneous calculation of the above equation:

${t}_{total\text{\_}1}=\frac{{\rho}_{102.5}-{\rho}_{100}\times V+44\times {\rho}_{100}}{0.85\times 0.49\pi \sqrt{2\times \left(160-{P}_{1}\right)\times 0.8711}}$

Therefore, the check valve opening time is:

${t}_{1}=\frac{10{t}_{total\text{\_}1}}{100-{t}_{total\text{\_}1}}$

100 - 200 ms: Similarly, we can get:

${t}_{total\text{\_}2}=\frac{{\rho}_{105}-{\rho}_{102.5}\times V+44\times {\rho}_{102.5}}{0.85\times 0.49\pi \sqrt{2\times \left(160-{P}_{2}\right)\times 0.8711}}$

${P}_{2}=102.5\text{\hspace{0.17em}}\text{MPa}$

${t}_{2}=\frac{10{t}_{total\text{\_}2}}{100-{t}_{total\text{\_}2}}$

200 - 300 ms:

…

1900 - 2000 ms:

${t}_{total\text{\_}20}=\frac{{\rho}_{150}-{\rho}_{147.5}\times V+44\times {\rho}_{147.5}}{0.85\times 0.49\pi \sqrt{2\times \left(160-{P}_{20}\right)\times 0.8711}}$

${t}_{20}=\frac{10{t}_{total\text{\_}20}}{100-{t}_{total\text{\_}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

Table 1. Adjustment table of check valve opening time every 100 ms.

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:

$\{\begin{array}{l}{m}_{out}={m}_{in}\\ {m}_{out}=44\times {\rho}_{150}\\ {m}_{in}=CA\sqrt{\frac{2\Delta P}{{\rho}_{160}}}\times {\rho}_{160}\times {t}_{total}\end{array}$

The calculation is:

${t}_{total}=6.992\text{\hspace{0.17em}}\text{ms}$ , $t=0.752\text{\hspace{0.17em}}\text{ms}$ .

4.4. Solution to Questions

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 Table 1.

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.

5. Evaluation of the Model

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.

Appendix

Explanation of Symbols Used in This Article

NOTES

*Corresponding Author: Yu Quan.

Conflicts of Interest

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

[1] |
Breit, D., Feireisl, E. and Hofmanová, M. (2020) Solution Semiflow to the Isentropic Euler System. Archive for Rational Mechanics and Analysis, 235, 167-194.
https://doi.org/10.1007/s00205-019-01420-6 |

[2] |
Buium, A. and Previato, E. (2017) Arithmetic Euler Top. Joumal of Number Theory, 173, 37-63 https://doi.org/10.1016/j.jnt.2016.09.024 |

[3] |
Milici, C., Machado, J.T. and Dr?g?nescu, G. (2018) On the Fractional Comuspirals. Communications in Nonlinear Science and Numerical Simulation, 67, 100-107.
https://doi.org/10.1016/j.cnsns.2018.07.004 |

Journals Menu

Copyright © 2021 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.