CC BY
7
Determination of Flow Characteristics in Technological Processes with Controlled Pressure
T. Reader, V.A. Tenenev, A.A. Chernova
Kalashnikov Izhevsk State Technical University, Studencheskaya str., 7, Izhevsk 426069, Russia
Received 22.04.2020
Accepted for publication 02.09.2020
Abstract
Assessment of the environmental and economic performance of a safety valve requires information about the flow of the substance through the valve when it is actuated. The goal of this paper was to determine the flow rate of the discharged substance and the mass flow rate of the substance entering the safety valve when it is actuated.
Proposed a mathematical model to describe the processes occurring in the valve. The model includes a system of differential equations describing the physical laws of conservation in the internal volume of the valve and differential equations, which link the value of gas flow through the valve with the pressure and the amount of movement of the shut-off disk. Used a modified method by S.K. Godunov to solve gas-dynamic equations.
Established that the determination of the flow and power characteristics of the valve requires the preliminary construction of a mathematical model of the safety valve operation. Based on this, proposed a method for determining the flow rate of the discharged substance and the mass rate of the substance entering the safety valve when it is actuated.
Obtained the flow characteristics of the valves under review and the dynamics of movement of the shut-off disc of the valve, as well as the dependence of the pressure change on the opening time of the valve. Comparison of the calculated values with available experimental data gives good agreement of results (no more than 5.6 % for a gas flow rate, under 10 % for the movement of the valve and change the arrival of gas in time using the standard deviation function of the flow characteristics of 0.6 %), confirms the correctness of the defined mathematical model, used numerical schemes and algorithms, as well as the proposed method and recoverability of the arrival of gas in a pressure-time curve.
Keywords: safety valve, flow measurement, pressure monitoring, numerical simulation. DOI: 10.21122/2220-9506-2020-11-3-204-211
Адрес для переписки:
А.А. Чернова
Ижевский государственный технический университет, ул. Студенческая, 7, г. Ижевск 426063, Россия e-mail: alicaaa@gmail.com
Для цитирования:
T. Reader, V.A. Tenenev, A.A. Chernova.
Determination of Flow Characteristics in Technological Processes
with Controlled Pressure.
Приборы и методы измерений.
2020. - Т. 11, № 3. - С. 204-211.
DOI: 10.21122/2220-9506-2020-11-3-204-211
Address for correspondence:
A.A. Chernova
Kalashnikov Izhevsk State Technical University, Studencheskaya str., 7, Izhevsk 426069, Russia e-mail: alicaaa@gmail.com For citation:
T. Reader, V.A. Tenenev, A.A. Chernova.
Determination of Flow Characteristics in Technological Processes
with Controlled Pressure.
Devices and Methods of Measurements.
2020, vol. 11, no. 3, pp. 204-211.
DOI: 10.21122/2220-9506-2020-11-3-204-211
Определение расходных характеристик в технологических процессах с контролируемым давлением
Т. Редер, В.А. Тененев, А.А. Чернова
Ижевский государственный технический университет имени М.Т. Калашникова, ул. Студенческая, 7, г. Ижевск 426063, Россия
Поступила 22.04.2020 Принята к печати 02.09.2020
Оценка экологичности и экономичности работы предохранительного клапана требует информации о расходе вещества через клапан при его срабатывании. Целью данной работы являлось определение величины расхода сбрасываемого вещества и массовой скорости поступления вещества в предохранительный клапан при его срабатывании.
Для описания процессов, протекающих в клапане, предложена математическая модель, включающая систему дифференциальных уравнений, описывающих физические законы сохранения во внутреннем объёме клапана и дифференциальные уравнения, связывающие величину расхода газа через клапан с давлением и величиной перемещения запорного диска. Для решения газодинамических уравнений применялся модифицированный метод С.К. Годунова.
Установлено, что определение расходной и силовой характеристик клапана требует предварительного построения математической модели функционирования предохранительного клапана. На основании чего предложена методика определения величины расхода сбрасываемого вещества и массовой скорости поступления вещества в предохранительный клапан при его срабатывании.
Получены расходные характеристики рассматриваемых клапанов и динамика перемещения запорного диска клапана, а также зависимость изменения давления от времени открытия клапана. Сопоставление расчётных значений с имеющимися экспериментальными данными даёт хорошее совпадение (не более 5,6 % для расхода газа, менее 10 % для перемещения клапана и изменения прихода газа во времени при среднеквадратичных отклонениях функции расходных характеристик 0,6 %) результатов, подтверждает корректность сформулированной математической модели, используемых численных схем и алгоритмов, предложенной методики и возможность восстановления прихода газа по кривой давление-время.
Ключевые слова: предохранительный клапан, измерение расходных характеристик, контроль давления, численное моделирование.
БОТ: 10.21122/2220-9506-2020-11-3-204-211
Адрес для переписки:
А.А. Чернова
Ижевский государственный технический университет, ул. Студенческая, 7, г. Ижевск 426063, Россия e-mail: alicaaa@gmail.com
Address for correspondence:
A.A. Chernova
Kalashnikov Izhevsk State Technical University, Studencheskaya str., 7, Izhevsk 426069, Russia e-mail: alicaaa@gmail.com
Для цитирования:
T. Reader, V.A. Tenenev, A.A. Chernova.
Determination of Flow Characteristics in Technological Processes
with Controlled Pressure.
Приборы и методы измерений.
2020. - Т. 11, № 3. - С. 204-211.
DOI: 10.21122/2220-9506-2020-11-3-204-211
For citation:
T. Reader, V.A. Tenenev, A.A. Chernova.
Determination of Flow Characteristics in Technological Processes
with Controlled Pressure.
Devices and Methods of Measurements.
2020, vol. 11, no. 3, pp. 204-211.
DOI: 10.21122/2220-9506-2020-11-3-204-211
Introduction
Chemical and technological processes, oil and gas pipeline systems, and high-pressure devices require the use of pressure level monitoring. The safety valve refers to the class of pipeline fittings designed to protect against mechanical destruction of equipment and pipelines, which is caused by overpressure, using automatically relief of excess working medium from systems and pressure vessels in excess of the set pressure. The valve should also ensure that the discharge of the medium is stopped when the operating pressure is restored [1-3]. Dangerous overpressure may occur in the system due to foreign factors (improper operation of equipment, heat transfer from exterior sources, incorrectly assembled thermal and mechanical circuit), as well as due to internal physical processes caused by an event beyond the normal operation.
In this case, the safety valve refers to quickoperation devices [4] that require the use of specialized quick-operation measuring instruments [5], the range of operation of which is quite limited [6]. However, we need to find the values of the discharge flow rate to assess the economic and environmental impacts. The determination of the mass rate of the substance input due to the process deviation from the normal technological mode is required to eliminate the causes of abnormal operation.
If there is a record of measurement of the current pressure in the tank, the value of gas arrival, GP, and gas flow through the valve, Gv, can be determined based on the mathematical model of the safety valve operation. Therefore, the goal of this paper was to determine the flow rate of the discharged substance and the mass flow rate of the substance entering the safety valve when it is actuated.
Mathematical model of gas arrival calculation
To obtain the dependence between the flow characteristics and the pressure in the tank, let us consider the mass balance equation:
wdP=Gp- G •
dp _ kRT(0) ^ p \ k
dt
W
Gp - G )
p(0)
(1)
where k, R is the adiabatic exponent and the gas constant of the working medium; 7X0), />(0) are the initial values of temperature and pressure in the tank.
The gas flow through the valve is a function dependant on the pressure and the value of the movement, X, of the shut-off disc: Gv = Gv(p, X). We can obtain this dependence either experimentally or by calculation. We propose to find this dependence from the numerical solution of the equations of the mathematical model of valve operation [7, 8] with the verification using the results of experimental measurements for a specific type ofvalve. Equation (1) gives us the expression for the gas arrival:
dp W
Gp = Gv (X, p)) +
k-i
'p(0) ^ ^
dt kRT(0)
(2)
Expression (2) includes the value ofthe movement of the valve disc, X(t). We find this dependence from the solution of the disk movement equation [7].
d ^ dX
m^ = Ff - Fs , ^ =
(3)
The movement of the disk in the axial direction is determined by the action of the force from the gas, Ff, and the elastic force of the spring, Fs = Ks(X + X0). Here: Ks is the spring stiffness coefficient; X0 is the initial spring compression (preload): n is the speed of the disk movement. Initial conditions: X(0) = 0, n(0) = 0. The height of the disk lift is limited by the value Xk.
For the numerical solution of the disk movement equations, we introduce a difference grid
{ < 'i < ••• < tn < ••• < tq = tk, K = tn+l - tn} and a two-step
difference scheme [8]:
To calculate the movement of the disk, we should define the dependence of the force acting on the disk from the gas on the movement and pressure, Fr(X, p)• To calculate dependencies, Gv = Gv(p, X), Fc(X, p), we use the following mathematical model
Xn
(( ( X", p(tn ))-KsX 0 ) + KdX" +
Xn - X"-1 X"
ml--1--
h h
к h + m
s " h
where: p is the gas density; W is the operation capacity of the tank; t is the time.
In the adiabatic approximation, we obtain an equation that links the pressure change, p, to the flow characteristics, GP, Gv:
Mathematical model of safety valve operation
Processes occurring in the valve are described by a mathematical model in the form of a system of
k-i
differential equations describing the physical laws of conservation in the internal volume of the valve. For a gas safety valve, we consider the processes of internal gas dynamics within the framework of the viscous compressible gas model. The change in the thermodynamic parameters of the gas in the tank is subject to the equation of state: p = pRT.
The continuity equation is written down general form:
in
dP + pVU = 0,
dt K
(4)
where U is the gas flow rate vector. The momentum conservation equations and the energy conservation equation have the following form:
p dU = -Vp + DivP;
dt
d i— dt
CT +
U_
2
2 Л
= V(U ) + Vq,
(5)
(6)
in a cylindrical coordinate system, we construct a difference grid (Figure 1b) in the 9 = const plane using the complex method of boundary elements [13]. For equations written down in a Cartesian coordinate system, the difference grid is unstructured. The gradients of variables included in the tensor components P are calculated in the middle of each face through the values of variables in the surrounding control volumes, as described in [14]. For time integration, we use a two-step Runge-Kutta scheme with second-order accuracy.
where P is the tensor of viscous stresses; Cv is the specific heat capacity of the gas at a constant volume; q is the heat flow vector.
The system of equations (4)-(6) is supplemented by equations for the transfer of the kinetic energy of turbulence and the rate of turbulence dissipation [9]. The axial component of the force acting on the disk from the gas is determined by the integral on the disk surface, S,
d-
Ff = j PdS.
(7)
The calculation area (Figure 1a) is divided into two parts.
The first part is axially-symmetric and it uses a cylindrical coordinate system (x, r, 9). The second part uses a rectangular coordinate system (x, y, z). We use the control volume method for the numerical solution of a system of gas-dynamic equations. We determine gas parameters at the boundaries of control volumes using the S.K. Godunov method [10]. To increase the order of approximation of the Godunov difference method, we use the MUSCL (Monotone Upwind Schemes for Conservation Laws) scheme. In accordance with this scheme, we define the values of gas-dynamic parameters for solving the problem of discontinuity decay using extrapolation with a limiter [11, 12]. For equations written down
Figure 1 - Valve diagram (a), difference grid (b), and flow structure (c)
The gas-dynamic equations (4)-(7) are solved numerically together with the disk movement equations (3). The difference grid is adapted to the disk movement.
We consider two designs of the LESER safety valve, 2J3 and 441. The gas-dynamic force and flow through the valve are written down as: Ff(X,p) = py(X), Gv(p, X) = py(X, where y(X), Y(X) functions are characteristics of the valve type. The calculated and experimental dependencies are shown in Figure 2.
b
a
d
c
300
g200 150 100
experiment 441
-calculate 441
experiment 2J3 —calculate 2.T3
6
X, mm
10
12
Figure 2 - Valve characteristics
The calculated dependencies, y(X), are approximated by 5th-degree polynomials,
5
Y (X) = with coefficients, pk, in Table 1.
Table 1
Polynomial coefficients for
Valve ß*
441 182.5704 2.2713 2.4635 -0.3243 0.0155 -0.0002
2J3 123.8281 9.6135 0.1595 -0.3951 0.0570 -0.0024
Polynomial coefficients for y(X)
Valve ß*
441 0 0.0345 -0.00071
2J3 0 0.0292 -0.001
Valve X0, mm Ks, N/m L, m Dt, m
441 0.017 54700 1.2 0.05
2J3 0.025 26000 1.2 0.05
serves as a gas source for the vessel tank with pressure control. The intensity of the gas arrival is regulated by the operator during the test (regulator, Rx). One of the goals of the tests was to check the stability of the valves with a long pipe between the tank and the valve. A pipe with a length of L and a diameter of Dt connects the tank to the valve. Pressure and temperature are measured at the inlet and at the end of the pipe. The movement of the rod, X", is measured on the valve. The pressure is measured using a P900 series load cell with an error of 0.2 %. The force is measured using a U2B sensor with an accuracy class of 0.2. The movement of the disk is controlled by a WA standard displacement inductive transducer with an error of 1 %. Sensor readings are recorded in 0.000417 seconds.
The flow characteristics of valves obtained by calculation are approximated by 2nd-degree polynomials with coefficients from Table 2.
Table 2
The results of tests processing
The scheme for testing safety valves in accordance with the ASME PTC 25-2014 standard is shown in Figure 3.
The characteristics of the tested LESER 441 and 2J3 valves are shown in Table 3.
Table 3
Characteristics of valves
The tests were carried out in the air. The vessel0 tank has a large volume under high pressure and
Figure 3 - Safety valve testing scheme: vessel0 - a tank with a large volume under high pressure, which is the source of gas for the vessel tank; vessel - a pressure-controlled tank with operating parameters P, T, W; R1 - gas flow regulator; Gp - gas arrival; P1, T1 - gas pressure and temperature at the outlet of the tank with pressure control; L - length of a pipe with a diameter, Dt, connecting the tank to the valve; P2, T2 - gas pressure and temperature at the inlet to the valve; Gv - gas flow through the safety valve
The use of formula (2) for calculating the gas arrival value based on the results of experimental measurements of pressure, p(t) , and disk movement, X(t ), requires calculating the derivative —.
dt
Differentiation of functions based on the results of experimental measurements is an incorrect operation. Figure 4 shows the measured dependence, p(t) (red line).
The prese nce of acoustic oscillations and measurement errors leads to the occurrence of high-frequency oscillations. Therefore, the results of measuring all parameters were smoothed using smoothing cubic splines. The value of the derivative was determined by the value of the relevant spline coefficient. The smoothed dependency, p(t), is shown in Figure 4 with a black line.
Figure 4 - Pressure versus time (valve 441)
The time change of the value of air arrival in the vessel tank, which is calculated using the formula (2), equations (3), and pre-calculated functions, y(X), Y(X), is shown in Figure 5.
Figure 5 - Dependence of gas intake on time
The scheme of the valve tests under review provided for measuring the movement of the valve disc, as well as the pressure and temperature in two sections of the pipe (Figure 1). Based on this information, we can determine the value of gas flow through the valve during testing and compare it with the calculated dependence, Gv(p, X) = py(X). If p2, T2 are the pressure and temperature at the end of the pipe before the valve, we determine the mass flow of gas through the pipe using the following expression:
D2
P2
Gexp =9 2 2 (P - P2 )RT^
where 92 is the flow coefficient related to a sharp narrowing at the inlet to the pipe from the tank (Figure 6) and friction in the pipe. We determine the flow coefficient when solving the gas-dynamic problem in an axially-symmetric setting in the calculated area (Figure 6).
The value of the flow coefficient, 92 = 0.814. The measured movement of the valve disc is shown in Figure 7.
Figure 6 - Fragment of the vessel - pipe design area
Figure 7 - Moving the valve disc
The calculated and experimental flow values corresponding to these values, X, are shown in Figure 8.
Figure 8 - Calculated and experimental dependencies for gas flow
The experimental and calculated, Y(X) = 0.0345X - 0.00071X2 , flow characteristics of the 441 valve are shown in Figure 9.
The average square deviation is 0.006. We may conclude from the presented comparison of the results of calculations and experiments that the design characteristics of the valve provide the determination of gas arrival in the controlled tank only through the pressure-time dependence.
Let us review the test results of the 2J3 valve. The time change in pressure is shown in Figure 10.
0,4 i
experiment —calculate
4 6 X, mm
Figure 9 - Calculated and experimental characteristics
10 12
flow
Figure 10 - Time dependence of pressure (2J3 valve)
During this time, the three-time actuation of the valve occurs, as shown in Figure 11.
Figure 11 - Measured and calculated movement of the 2J3 valve disc
The time interval under review corresponds to the dependence of the gas arrival on time, which is obtained by the calculation and confirmed in experiments (Figure 12).
o —experiment —calculate
The comparison of the change in the value of the disk movement obtained in the experiment and by the calculated method (solution of equations (3) using the power characteristic of the valve 2J3, y(X) = 123.8281 + 9.6135X + 0.159X2 - 0.3951X3 + + 0.0570X4 - 0.0024X5), which is shown in Figure 11, demonstrates that the power characteristic of the valve was obtained correctly.
Figure 12 - Change of gas arrival in time (valve 2J3)
Conclusion
We have shown that in order to find the arrival of a gas mass in a controlled technological tank, which leads to increase in pressure, it is required to have the following information: the pressure-time dependence recorded in the technological process; calculated flow and power characteristics of the valve. We confirmed during our work that the determination of the flow and power characteristics of the valve requires the preliminary construction of a mathematical model of the safety valve operation. Based on this, we proposed a method for determining the flow rate of the discharged substance and the mass rate of the substance entering the safety valve when it is actuated.
Following the numerical simulation, we obtained the flow characteristics of the valves under review and the dynamics of movement of the shut-off disc of the valve, as well as the dependence of the pressure change on the opening time of the valve. Comparison of the calculated values with available experimental data gives good agreement of results (no more than 5.6 % for a gas flow rate, under 10 % for the movement of the valve and change the arrival of gas in time using the standard deviation function of the flow characteristics of 0.6 %), confirms the correctness of the defined mathematical model, used numerical schemes and algorithms, as well as the proposed method in general and recoverability of the arrival of gas in a pressure-time curve in particular.
Acknowledgments
We have performed the paper within the framework of the Research Activity of the Udmurt Federal Research Center of the Ural Branch of the Russian Academy of Science "Artificial Intelligence in the Development, Training, and Support of Expert Systems for the Representation and Use of Knowledge in Natural-Scientific, Technical, and Socio-Humanitarian Fields" AAAA19-119092690104-4
References
1. Faizov A.R., Churakova S.K., Pruchay V.S. Simulation mode of emission from the distillation column. Bashkir chemistry journal, 2013, no. 2, pp. 52-54.
2. Beune A., Kuerten J.G.M., Schmidt J. Numerical calculation and experimental validation of safety valve flows at pressures up to 600 bar. AIChE Journal, 2011, vol. 57, pp. 3285-3298.
3. Bernad S.I., Susan-Resiga R. Numerical model for cavitational flow in hydraulic poppet valves. Modelling and Simulation in Engineering, 2012, vol. 2012, pp. 1-10. DOI: https://doi.org/10.1155/2012/742162
4. Tusyuk S.K., Belyanskaya E.S. Diagnostic internet - testing first - year students in the tula state university. Izvestiya TulGU [News of the Tula state university. Technical sciences], 2012, no. 12-2, pp. 262-267 (in Russian).
5. Boronenko M.P. Overview of application of high-velocity television measuring systems in a physical experiment. Yugra state university, 2014, no. 2(33), pp. 43-55.
6. Fyodorov S.A. Embedded software of high speed IP camera for optical diagnostic of the fleeting processes. Proceedings of TUSUR, no. 1, 2 part [25], pp. 225-229.
7. Song X., Cui L., Cao M., Cao W., Park Yo., Dempster W.M. A CFD analysis of the dynamics of a direct-operated safety relief valve mounted on a pressure vessel. Energy Conversion and Management, 2014, vol. 81, pp. 407-419.
DOI: 10.1016/j.enconman.2014.02.021
8. Reader T., Tenenev V.A., Paklina N.V. Numerical 3D Simulation of Safety Valve Gas Dynamics. Vestnik Izhevsk state technical University named after M.T. Kalashnikov [Bulletin of Kalashnikov ISTU], 2018, vol. 21, no. 4, pp. 174-181 (in Russian).
DOI: 10.22213/2410-9304-2018-2-28-40
9. Volkov K.N., Emelyanov V.N. Large-eddy simulation in the calculation of turbulent flows. Moscow. 2008, 368 p.
10. Godunov S.K., Kulikov I.M. Computation of discontinuous solutions of fluid dynamics equations with entropy nondecrease guarantee. Computational Mathematics and Mathematical Physics, 2013, vol. 53, no. 8, pp. 1179-1182. DOI: 10.1134/S0965542514060086
11. Smirnova N.S. Comparison of flux splitting schemes for numerical solution of the compressible Euler equations. Proceedings of MIPT, 2018, vol. 10, no. 1, pp. 122-141.
12. Wesseling Pieter Dr. Principles of computational fluid dynamics. Springer series in computational mathematics, ISSN 0179-3632; 29. Mathematics Subject Classification, 1991: 76M,65M, XII, 644 p.
DOI: 10.1007/978-3-642-05146-3
13. Tenenev V.A., Gorokhov M.M., Rusyak I.G. Numerical study of particle Gorenje in a two-phase flow. Mathematical modeling, 1997, vol. 9, no. 5, pp. 87-96.
14. Wesseling P., Segal A., Kassel C.G.M. Computing Flows on General Three-Dimensional Nonsmooth Staggered Grids. Journal of Computational Physics, 1999, vol. 149, pp. 333-362.
DOI: 10.1007/978-94-017-1564-5 2
