NUMERICAL SIMULATION OF METHANE UNDER PHASE TRANSFER CONDITIONS WITH ITS REAL PROPERTIES

The technique and results of calculation of parameters of the flow of methane in single and two-phase state is discussed. The flow is described by the one-dimensional nonstationary equations of a continuous medium. Is taken into account the processes of friction, heat transfer and phase transitions. The results of the application of various schemes of numerical integration of the equations are discussed. In this paper the parameters of the flow of metal having generally two phases investigates. For a description of the thermodynamic properties of the working wide range of parameters an approximation of the experimental data suggested in the article [1] in a fluid equation of state is used. The essence of this approach lies in the fact that the experimental data for methane were approximated by the authors for the molar Helmholtz free energy, as well as approximating were chosen easily integrable and is easily differentiable functions. The parameters are determined by the current system of equations [2]: ρ ρ τ condensate u G x , (1) 2 (ρ ) ρ friction condensate u p w F G w x , (2) 2 2 ρ(ε 0.5 ) ρ ( 0.5 ) τ out condensate condensate u u h u Q G h x , (3) Where τ, , ρ, , , ε, , , , condensate friction out x u p h G F Q , respectively, time, coordinate, density, velocity, pressure, specific internal energy, specific enthalpy, mass flow rate of the phase transition from the condensation on the walls, or evaporation from the surface of the film of condensation on the walls of the channel, the force of friction and heat flux. In this case , , condensate friction out G F Q given on unit volume. a Corresponding author: bvborisov@tpu.ru DOI: 10.1051/ C © Owned by the authors, published by EDP Sciences, 201 / 0 0 (201 ) 201 epjconf EPJ Web of Conferences , 0 01 1 61 6 6 10 10 1 1 1

In this paper the parameters of the flow of metal having generally two phases investigates.For a description of the thermodynamic properties of the working wide range of parameters an approximation of the experimental data suggested in the article [1] in a fluid equation of state is used.The essence of this approach lies in the fact that the experimental data for methane were approximated by the authors for the molar Helmholtz free energy, as well as approximating were chosen easily integrable and is easily differentiable functions.
The parameters are determined by the current system of equations [2]: ρ(ε 0.5 ) ρ ( 0.5 ) τ Where τ, , ρ, , , ε, , , , The initial conditions are given a uniform initial distribution of density, velocity and internal energy as the rest gas, and for uniformly moving stream for whole canal.Border conditions: System ( 1) -( 4) together with the state equation is solved numerically.Using of equation of state is implemented as a procedure, an algorithm which allows for the input parameters of the density and specific internal energy, found from the system of equations ( 1) -( 3), determine the values of the appropriate temperature and then all the other thermodynamic parameters [1], used for numerical integration of ( 1) - (3).
Simulation is performed with taking into account of the hydraulic and thermal interaction with the walls of the tube.For this the numerical implementation of the control (finite) volume method has been selected [2].To determine the greatness of the original stream was used TVD scheme with approximation of derivatives against the flow and using the principle of the minimum of the derivative.The use of this approach was followed by the appearance of significant instability that is a consequence of a very strong dependence of the pressure and specific internal energy according to the equation of state.The introduction of artificial viscosity in the flow of the total momentum and the total enthalpy in the form proposed in 1944 and von Neumann J. Richtmyer R.D. [2] only partly solved the problem.Scheme proposed by Lax P.D. in 1954 [2] for numerical integration helped fully repay nonphysical oscillation of solutions.The very limited number of Courant was used.This scheme has a substantial viscosity, which allows it to avoid significant fluctuations nonphysical solutions.But some embodiments, the calculation results included an unjustified increase in the pressure in the steady stream through the channel of constant cross section with increasing density and fall of the speed, in some cases, when a two-phase flow simulation included on hydraulic losses and heat exchange with the wall of the channel for providing condensation conditions.
Example of a detailed calculation for the steady flow is shown in Figure 1.This the pressure behavior is explained due to on the one hand the strong diffusion properties of Lax's schemes, and on the other hand, the special features determining the pressure according to of the obtained values of density and internal energy [1].Authors have the experience the use for determine of gas-dynamic parameters in difficult areas and in conditions of high gradients of pressure due to monotone SK Godunov scheme of decay of a discontinuity [3,4].This method of determining the values of flows calculated on the channel border and at the borders of cells is well physically justified.For an ideal gas, this method is deeply studied and widely used.
As for implementation for a real gas of such an algorithm is difficult.Based on the proposal of SK Godunov [3] trusts that the flow does not contain strong discontinuities parameters.In this case the shock waves and the waves underpressure converted to sound waves.Dynamic relations on the discontinuities derived from the equations of continuity, momentum and energy (1-3) are determined by [3]: Density Ro , pressure P , velocity U and specific energy E at the interface between the cell is determined according to the following analysis of the decay of an arbitrary discontinuity options: . There is the supersonic flow "to right".Parameters at the cell border assumed equal the parameters of the left cell: . There is the supersonic flow "to left".Parameters at the cell border assumed equal the parameters of the cell right: .The density and internal energy is determined from (5) -(7) according to the "left" elementary wave: Option (D): and 0 CD U .The density and internal energy is determined from (5) -( 7) according to the "right" elementary wave: The solutions of elementary waves are used to implement the boundary conditions in accordance with the sign of the velocity at the boundary.It is believed that the channel is finished at borders of the large volumes, which are having a complete set of relevant thermodynamic parameters.For definiteness we consider the ratio of parameters on the left border.It is believed that the flow is carried out in the subsonic mode, and in accordance with the pressure at the corresponding pressure in the volume ( 0 P p ).Speed and density determined by the decision on the appropriate wave:  Thanks to the SK Godunov circuit application and the appropriate implementation of the boundary conditions is received significant impact to the calculation of the parameters (Fig. 3) obtained in the same option of that the respective calculation of Lax PD (Fig. 1).The presented method of calculation can be recommended for use in the process design of the special devices.

a/
, time, coordinate, density, velocity, pressure, specific internal energy, specific enthalpy, mass flow rate of the phase transition from the condensation on the walls, or evaporation from the surface of the film of condensation on the walls of the channel, the force of friction and heat flux.In this case , , Corresponding author: bvborisov@tpu.ruDOI: 10.1051/ C Owned by the authors, published by EDP Sciences, 201

Figure 1 .
Figure 1.Parameters in the two-phase flow field of methane produced by the scheme Lax PD.

Figure 2 .
Figure 2. Schemes options decay of an arbitrary discontinuity.
the two other options the speed and pressure of the contact discontinuity is determined according to by SK Godunov proposal[3 the specific enthalpy corresponds to a value in the adjacent volume and in case outflow ( 0 U) the specific internal energy's value is determined according to the following (5) -(7) the relation:

Figure 3 .
Figure 3. Parameters in the two-phase flow field of methane produced by the scheme Godunov SK.