Simplified Coarse-Grained Dynamic Model for Real Gases

A simplified model is proposed for an easy understanding of the coarsegrained technique and for achieving a first approximation to the behavior of gases. A mole of a gas substance, within a cubic container, is represented by six particles symmetrically moving. The impacts of particles on container walls, the inter-particle collisions, as well as the volume of particles and the inter-particle attractive forces, obeying a Lennard-Jones curve, are taken into account. Thanks to the symmetry, the problem is reduced to the nonlinear dynamic analysis of a SDOF oscillator, which is numerically solved by a stepby-step time integration algorithm. Five applications of proposed model, on Carbon Dioxide, are presented: 1) Ideal gas in STP conditions. 2) Real gas in STP conditions. 3) Condensation for small molar volume. 4) Critical point. 5) Iso-kinetic energy curves and iso-therms in the critical point region. Results of the proposed model are compared with test data and results of the Van der Waals model for real gases.


Introduction
Recently, in Computational Chemistry, the coarse-grained molecular dynamics technique is often used, by which millions of molecules are represented by a few hundred particles [1]- [13].For example, if a mole of a gas substance is simulated by a thousand particles, by use of Avogadro number, it is noticed that every par-handling of chemical problems becomes possible and usually a satisfactory approximation to observed behavior is achieved.
If we consider an amount of a gas substance, represented by a few particles, first in a large container (Figure 1(a)) and then in a small container (Figure 1(b)), the following two observations can be made [14]: 1.In the large container of Figure 1(a), the volume of particles is not significant, compared with the volume of the container.On the contrary, in the small container of Figure 1(b), the volume of particles is significant.
2. In the large container of Figure 1(a), the mean distance between a couple of particles is large, so, as well known from Physical Chemistry [14] and described by Lennard-Jones curve [15], the inter-particle attractive forces result small up to negligible.On the contrary, in the small container, the mean distance, between a couple of particles, is small, so the inter-particle attractive forces exhibit significant values.
For the above two reasons, for a quite large molar volume, as in Figure 1(a), the volumes of particles and the inter-particle attractive forces can be ignored.
So, the amount of gas substance under consideration obeys the ideal gas laws.
On the contrary, for a small molar volume, as in Figure 1(b), the particles volumes and the inter-particle attractive forces must be taken into account.That is, we have a real gas, which significantly deviates from the behavior of ideal gases.[16], by taking into account the molecular volumes and the inter-molecular attractive forces, developed a semi-rational, semi-empirical model, which is simple and exhibits a satisfactory approximation to the observed behavior of real gases.

J. D. van der Waals
The kinetic behavior of gases, ideal and real ones, is often described by the Maxwell-Boltzmann stochastic model [17].The stochastic models are accurate but complicated.On the other hand, they obey some required symmetries.And it is recognized [18] [19] that, alternatively to a stochastic model, a symmetric deterministic model can be used, which is much simpler, but usually exhibits satisfactory approximation to test data.Also, a very coarse-grained model can be used, that is consisting of very few, very large particles.This is similar to the concept of fundamental vibration mode of structural dynamics.Where there exist a lot of high vibration modes (with small periods and usually small amplitudes, too), which are of negligible interest, but complicate the computation and require a very small time steplength.Whereas, the fundamental vibration mode is the simplest mode and, at the same time, the most representative of the dynamic behavior of the structure.
In the present work, such a symmetric deterministic model is proposed for real gases, which is very coarse-grained, that is, it consists of very few -very large particles, and is compared to corresponding test data [14] [20], as well as to results of the Van der Waals model [14] [16] for real gases.
In the recent literature on the coarse-grained technique [1]- [13], advanced models are proposed, for the accurate description of the behavior of real materials, which can be used in the Design.Whereas, the proposed here simplified model, aims to an easy understanding, of the coarse-grained technique, by researchers of other than Chemistry fields and to a first approximation to the behavior of gases.

Proposed Model
A mole of a gas substance is considered, within a cubic container of side L   Obviously, all the six particles, as they move symmetrically, they pass simultaneously from the above three characteristic states.
The mutual repulsive forces F between a particle and a container wall (Figure 4(a)) are described by the repulsive part of a Lennard-Jones curve [14] [15] (Figure 4(b)).If the perpendicular distance  of the center of particle from container wall is quite large: , then 0, F = that is no force F is developed between the particle and the wall.
On the contrary, if  is quite small: 1.1225 2 , D <  a mutual repulsive force F is developed between particle and wall, given by the equation (Figure 4(b)): where 2 r D =  and the determination of force coefficient F 0 of Lennard-Jones curve is described in the following Section 4.1.
The mutual forces F, attractive or repulsive, between any couple of particles, with a distance  of their centers (Figure 5 For any distance  of the centers of two particles, their mutual force F is given by the following equation, corresponding to the Lennard-Jones curves of where r D =  and 0.7071 Finally, the left particle L acts on particle R, by an, always attractive, force F  , which is given by Figure 5 and Figure 6 and Equation (2), with 2 y =  .
So, the horizontal force on particle R, due to inter-particle action, is and the total horizontal force on particle R, due to inter-particle action and impact on wall is, and the acceleration of particle R, under consideration, is, at any instant,

Step-by-Step Algorithm
It has been described, in the previous Section 2, how the proposed model is reduced, thanks to the symmetric movement of its six particles, to the study of the movement of the single particle R (Figure 6).So, the problem is reduced to the nonlinear dynamic analysis of a SDOF (single degree of freedom) oscillator, which can be solved numerically by a step-by-step time integration algorithm.
For this purpose, the algorithm of trapezoidal rule (or Newmark algorithm of constant average acceleration) is chosen, combined with a predictor-corrector technique, with two corrections per step [21], which has been proved simple and effective.

Flow-Chart
The flow-chart of the proposed algorithm is shown in Figure 7 and is briefly described below.
First, the constant input data are read: particle mass m and diameter D, force coefficient o F of L-J (Lennard-Jones) curve, side L of cubic container, time step-length t ∆ of the algorithm.
The initial conditions are read: position y, temperature T and speed v of the particle R under consideration.The initial speed v results from initial temperature T, by a thermodynamic postulate, which will be described in following section 4.1.
The subroutine L-J (Lennard-Jones) is called, which, from the initial position y of the particle, determines the initial forces i F and w F , acting on it, and its initial acceleration Figure 7. Flow-chart of the proposed step-by-step time integration algorithm.
Within each step of the algorithm, first the steps counter n is increased by 1 and time t by t ∆ .
Then, the prediction is performed, which determines the predicted values , p p y v of state variables and the subroutine L-J is called, which, from the given p y determines the predicted acceleration p γ .
The first correction, by trapezoidal rule, determines the first corrections , y v of the state variables.The subroutine L-J, from 1  y , finds the first correction of acceleration 1 γ .
The second and final correction finds the final values of , y v , for present step, and the subroutine L-J, from y, determines the final forces i F , w F and accele- ration γ , for present step.
The output, of present step of algorithm, is printed: steps counter n, time t, position y and speed v of the particle, forces i F due to inter-particle action and w F due to impact on wall, instantaneous total kinetic energy, for all six particles , where M = 6 m.
At the end of step of algorithm, three summations are made: The present force i F , due to inter-particle action, is summed to i F ∑ .The present w F , due to impact on wall is summed to Then, if the first cycle of oscillation has not yet been completed, we continue with the next step of the algorithm.
When the first cycle of oscillation is completed, by returning to the initial state, if we continued the algorithm, everything would be repeated the same, with only a small algorithmic damping.So, the algorithm is interrupted and the global output data are printed, which are: , and mean (rms) speed ( ) which, for small molar volumes, results significantly lower than the initial speed.

Computer Program
Based on the step-by-step time integration algorithm, described in the previous The program is written in the version Force 2.0 of Fortran, whose compiler is free available, even in Internet cafés.

Applications
The proposed simplified coarse-grained dynamic model for real gases is applied on Carbon Dioxide (CO 2 ), which exhibits a particular behavior in Critical point region, as it condensates for rather high temperatures, slightly lower than 31 C 304 K T = =  .
From the next Section 4.1, it is apparent that, in order to calibrate the proposed model on other gases, the following data are required: molar mass, incompressibility limit of molar volume, as well as temperature, pressure and molar volume at the Critical Point.

Determination of Parameters
The numerical values of parameters of proposed model are determined below, which will be used in the following applications: 1) The mass of a particle is 6 0.044 kgr 6 0.007333 kgr m M = = = , where M = 0.044 kgr is the molar mass of Carbon Dioxide.
2) The diameter D of a particle is determined on the basis of criterion of incompressibility of closely-packed equal spherical particles, as shown in Figure 8.According to experimental evidence [14] [16] [20], the in-compressibility limit of molar volume, for Carbon Dioxide, is about V = 50 cm 3 , which corresponds to a cubic container with side L = 3.684 cm.In Figure 8, the spherical particles of proposed model are shown, closely -packed in such a small container.The inter-particle distances are 1.1225D and the particle-wall distances are 1.1225D/2, as, for smaller distances, mutual repulsive forces begin to develop (see Figure 4 and Figure 5).So, on the basis of configuration of Figure 8. Closely-packed particles of proposed model, in a small container, in the limit of in-compressibility of Carbon Dioxide, according to experimental evidence [14] [16] [20].
3) The force coefficient o F of Lennard-Jones curve (Equations ( 1) and ( 2 A maximum stiffness max K appears in two cases: inter-particle collision (Figure 5 and Figure 6 and Equation ( 3)) and particle-wall impact (Figure 4 and Equation ( 1)).By linearization of branch AK in the Lennard-Jones curves in However, the cost, from using a further shorter time step-length t ∆ of the algorithm, is negligible, as the computing time, for the first oscillation cycle of the model, is only a few seconds.So, a 4 10 msec t − ∆ = is chosen, much shorter than that required by the above accuracy criterion of the algorithm, so that to achieve more accuracy.
6) The initial position of the particle R under consideration is 0 4 y L = , where L side of cubic container, as mentioned in Section 2 and according to Figure 2 and Figure 6.8) The initial speed 0 v , of the particle under consideration, is obtained from the initial temperature T 0 , by the thermodynamic postulate: ( ) where R = 8.3144 Joules −1 •K −1 is the value of gas constant for ideal gases.However, for small molar volumes, through the oscillation of the particle, the speed v is significantly reduced, which implies mean values of R much smaller than the initial one, as will be shown in the applications.

First Application. STP Conditions. Ideal Gas
A mole of Carbon Dioxide is considered, within a cubic container of side L = 28.195cm, that is volume In this first application, point particles are assumed, that is with zero volume, and the inter-particle attractive forces are ignored.So, we have an ideal gas.This case is simple, so it will be solved by hand.
Within the first cycle of oscillation, the particle, starting from the position 0 4 y L = (Figure 6), goes to impact on wall at right, where the speed is reversed.Then, an inter-particle collision occurs at left, where the speed is again reversed and the particle returns to the initial position.So, the particle runs twice the distance L/2 (Figure 6), with the constant speed v = 393.5 m/sec and the period of oscillation is  At inter-particle collision and particle-wall impact, the impulse-momentum conservation equation can be written:  as was expected for an ideal gas in STP conditions.
As the speed is constant, the total kinetic energy is, at any instant, as is due for an ideal gas.That is, the proposed model describes accurately the behavior of an ideal gas.

Second Application. STP Conditions. Real Gas
The same input data, of the previous first application, in STP conditions, are again considered, that is a cubic container with side L = 28.195cm, thus volume , are taken into account.So, we have a real gas, and the proposed step-by-step time integration algorithm is used, in order to follow the oscillation of the particle.
In Figure 11(a), is shown the large container of side L = 28.195cm, with the particles of diameter D = 1.35 cm, in the initial conditions.In the Figures 11(b)-(f), are presented, within the first oscillation cycle, the variations, with respect to time t, of five quantities: b) position y of the particle.c) speed v. d) inter-particle force i F .e) particle-wall force w F .f) total kinetic energy  that is, it slightly deviates from the ideal gas value.
The mean 2 nd power of speed is So, the mean (rms) speed is that is, slightly larger than initial speed.
And the mean kinetic energy is

Third Application. Condensation for Small Molar Volume
A small cubic container with side L = 4.55 cm, thus molar volume  The application run by the proposed step-by-step algorithm, with

K E Mv =
In the present application, because of the low initial temperature, thus low initial speed and kinetic energy, too, the inter-particle attractive forces i F over- come the kinetic energy of the particle, thus preventing it from reaching to impact on the wall.So zero particle-wall forces 0 w F = and zero pressure P result, which mean that a liquid phase exists.Because of the zero particle-wall forces, 0 w F = , the sum of inter-particle forces results zero, , for equilibrium, as noted in the Figure 12(d).
The mean 2 nd power of speed results

Fourth Application. Critical Point
The same small cubic container of side L = 4.55 cm of previous application is considered, which implies a volume  The application run by the proposed step-by-step algorithm, with The first oscillation cycle was completed in 562 steps, so the period is 0.0562 msec τ = In Figures 13(b)-(f), the variations, with respect to time t, of five quantities, are presented: b. position y of the particle.c. speed v. d. inter-particle force i F .e.

K E Mv =
The mean inter-particle force results 8.500 kN 562 15.12 kN.The pressure on the wall is The mean 2 nd power of speed is The present value of R is much smaller than the value R = 8.3144 of ideal gases.
This will be discussed in the fifth application of next section 4.6.The present application is adapted to the Critical point by its initial temperature 31 C 304 K T = =  , which is the critical temperature of Carbon Dioxide, according to test data [14].In the fifth application of next Section 4.6.The Critical point of CO 2 will be determined in two different ways: By the group of isokinetic energy curves of Figure 14 and by the group of iso-therms of Figure 16.Both cases are close to the Critical point of present application.

Fifth Application. Iso-Kinetic Energy Curves in the Critical Point Region
For side of cubic container ranging from L = 4.0 cm up to 9.0 cm, with a step 0.2 cm L ∆ = , that is, volume V = L 3 ranging from 64 cm 3 up to 729 cm  With the help of this graph, the iso-kinetic energy curves of proposed model, of Figure 14, have been transformed to the iso-therms shown in Figure 16, for temperatures ranging from T = 250 K up to 400 K with a step 10 K T ∆ = .The above iso-therms of proposed model are compared with corresponding ones of the test data of Eastman-Rollefson [14] [20], in Figure 17, as well as with those of Van der Waals model [14] [16], in Figure 18.
It is observed, in the Figure 17 and Figure 18, that the proposed model better represents the wave-shaped isotherms of the Van der Waals model, in the Vapor region, than the horizontal linear isotherms of the test data by Eastman-Rollefson.Also, the proposed model approximates better the larger incompressibility limit of the molar volume given by the Van der Waals model (Figure 18), than the smaller one of the test data by Eastman-Rollefson (Figure 17).

Conclusions
A simplified coarse-grained dynamic model, for real gases, is proposed.Five applications of this model, on Carbon Dioxide, are presented: 1) In STP conditions, by ignoring particle volume and inter-particle attractive forces, the proposed model accurately represents the behavior of an ideal gas.2) Again in STP conditions, but taking into account the particles volume and inter-particle attractive forces, the proposed model slightly deviates from the behavior of an ideal gas, as was expected.3) For a small molar volume and a low initial temperature, the inter-particle attractive forces overcome the initial kinetic energy of particles and prevent them from reaching at impact with container wall.So, a Liquid phase exists, with zero pressure.5) Iso-kinetic energy curves have been determined, by the proposed model, in the Critical point region of Carbon Dioxide.By comparing these iso-kinetic energy curves to corresponding iso-therms of test data by Eastman-Rollefson [14] [20] and to those of Van der Waals model [14] [16], a variation of values of gas constant R, in Critical point region, is revealed, ranging from 2.85 up to 5.0 Joules mole −1 •K −1 .With the help of this variation of values of R, the iso-kinetic energy curves of proposed model are transformed to iso-therms, which are compared to corresponding ones of test data by Eastman-Rollefson, as well as to iso-therms of Van der Waals model.And a better agreement is achieved between the proposed model and the Van der Waals model, as shown in Figure 18, as regards a larger in-compressible molar volume and particularly the wave-shaped iso-therms in the Vapor region.
The above five numerical experiments show that the proposed simplified model can approximate the observed behavior of real gases.
In the present work, in order to achieve simplicity, the accuracy is reduced.
However, if a refined version of the proposed model, with more particles, is developed, the accuracy can be improved.

Figure 1 .
Figure 1.An amount of a gas substance, represented by a few particles.(a) In a large container; (b) In a small container.

(Figure 2 )
, represented by six equal spherical particles, each one with mass m = M/6, where M molar mass.A reference axes system Ixyz is considered, with origin I at the center of cube and the axes x, y, z parallel to the principal directions of cube.The centers of the six particles are located on the axes x, y, z, initially at the middles of distances of I from the centers of six faces of cube, that is they have the six particles move symmetrically.So, by considering the plane Iyz (Figure 3(a)), what happens in this plane, the same happens in the planes Ixy, Ixz, too.

Figure 2 .
Figure 2. Initial positions , , o o o x y z of the six particles of proposed model in the cubic

Figure 3 .
Figure 3.The three successive characteristic states of the particles of proposed model, within the plane Ixy (a) of the container; (b) Initial state; (c) Impacts of particles on container walls; (d) Inter-particle collisions in the central region of the container.The arrows represent instantaneous velocities of the particles.
(a)), are described by the Lennard-Jones curve of Figure 5(b).In Figure 5(c) is shown enlarged the attractive part of this Lennard-Jones curve, because of the significance of attractive forces.

Figure 4 .
Figure 4. (a) Perpendicular distance  of a particle from container wall and mutual repulsive force F between particle and wall; (b) Repulsive part of a Lennard-Jones curve describing the function ( ) F  .

Figure 5 .
Figure 5. (a) Distance  and mutual inter-particle force F, attractive (+) or repulsive (−), between a couple of particles; (b) Lennard-Jones curve describing the function ( ) F  ; (c) The Lennard-Jones curve with enlarged its attractive part.
Thanks to the symmetrical movement of all the six particles, it is enough to study the movement of only one particle, let choose that on right part of axis Iy of Figure2and name it R (right), as shown in Figure6.If the distance of this particle from container wall at right is quite small 1wall repulsive force is activated, according to Figure4and equation (1), and let call this force F w .At left of Figure 6, the particle R interacts with the other five particles of the model.The relative position of point R under consideration with respect to four of these particles, F, B, O, U (front, back, over, under) is symmetric.So, the horizontal resultant F 4 of the four equal forces (attractive or repulsive), by which the points F, B, O, U act on point R is, according to Figure 5 and Figure 6 and Equation (2):

Figure 6 .
Figure 6.The particle R under consideration is moving on axis Iy.At right, it reaches up to impact with container wall.At left, it interacts with all the other five particles: (F, B, O, U) and L.

2 )
Mean particle-wall impact force w as is due for global equilibrium.3)Pressure on wall 2 w P F L = , in Pascals = N/m 2 , which, divided by 101,325 N/m 2 , turns to atm units.4) Mean 2nd power of speed 2 section 3.1.and the flow-chart of Figure 7, a simple and short computer program has been developed, with only about 45 Fortran instructions for the main program of step-by-step algorithm and only about 25 Fortran instructions for the L-J (Lennard-Jones) subroutine, that is totally only about 70 Fortran instructions.

4 )
)) is determined by calibration of proposed model on the Critical point of Carbon Dioxide.For a cubic container with the critical volume[14] [16][20], various values of o F are tried, until to achieve a value of pressure, with satisfactory approximation to the experimental critical pressure of CO 2 , point is described in Section 4.5.For the side L of cubic container, in STP (standard temperature-pressure) conditions, the value 28.195 cm L = is chosen, which corresponds to a volume V = 22.414 liters.And, in the Critical point region of CO 2 , values of L ranging from 4.0cm up to 9.0 cm are used, which correspond to container volumes V = L 3 ranging from 64 cm 3 up to 729 cm 3 .5)The time step-length t ∆ of the proposed step-by-step integration algorithm can be determined on the basis of the accuracy criterion of the algorithm[21],

Figure 4 (Figure 9 .
Figure 4(b) and Figure 5(b), the stiffness of the above two cases can be determined on the basis of Figure 9(a) and Figure 9(b), respectively.It is observed that both give the same value of stiffness, represented by Figure 9(c), which is 6 2 max 2 675 kN 2 0.1225 816.4 10 kgr sec 0.1225 0.0135 cm o K F D × = = = × ×

7 )
Initial temperatures, ranging from 0 50 C 223 K T = − =  in liquid phase region up to 0 100 C 373 K T = =  in gas phase region, are used.

In Figure 10
Figures 10(b)-(f), for the first oscillation cycle, the variations, with respect to time t, of five quantities, are presented: (b).Position y of the particle.(c).Speed v. (d).Inter-particle impulse i F dt . (e).Particle-wall impulse w F dt . (f).Total ki- netic energy, for all six particles,

Figure 10 . 2 .
Figure 10.First application.STP conditions.Ideal gas.(a) Large container of side L = 28.195cm with point-particles in the initial state.In the following diagrams, variations of five quantities, with respect to time t; (b) Position y of the particle; (c) Speed v; (d) Inter-particle impulse i F dt ; (e) Particle-wall impulse w F dt ; (f) Total kinetic energy 2 . . 1 2 K E Mv = .
∞ (tend to infinity).How- ever, as everyone, of the two above impulses, occurs once in an oscillation cycle, we can obtain the finite mean values of forces i F , w F , by simply dividing the above impulses by the period τ:

F
are opposite, as is due for equilibrium, and are noted in Figure 10(d) and Figure 10(e), respectively.The pressure on the wall is speed of the particle 0 393.5 m sec v = .However, now, the volume of particles, with diameter D = 1.35 cm, and the inter-particle attractive forces, described by a Lennard-Jones curve (Figure5, Equation (2)), with a force coefficient 0 675 kN F =

2 .
6312 steps of the algorithm have been performed, within the first oscillation cycle, with a time-steplength as is due for global equilibrium.The mean forces i F , w F are noted in the Figure 11(d) and Figure 11(e), respectively.The pressure on the wall is

Figure 11 . 2 .
Figure 11.Second application.STP conditions.Real gas.(a) Large container of side L = 28.195cm with particles of diameter D = 1.35 cm, in the initial state.In the following diagrams, variations of five quantities with respect to time t; (b) Position y of the particle; (c) Speed v; (d) Inter-particle force i F ; (e) Particle-wall force w F ; (f) Total kinetic en- ergy it slightly deviates from the corresponding value of ideal gas.The mean kinetic energy is noted on the diagram K.E.-t of Figure 11(f), for comparison.
, the same as in the Critical point of Carbon Dioxide, according to test data [14].And a low initial temperature 0 50 C 223 K T = − =  , which implies a low initial speed of the particle ( ) the small container, with the particles of diameter D=1.35cm, in the initial state, is shown in Figure 12(a).

Figure 12 . 2 .
Figure 12.Third application.Condensation for small molar volume.(a) Small container of side L = 4.55 cm with particles of diameter D = 1.35 cm, in the initial state.In the following diagrams, variations of four quantities, with respect to time t; (b) Position y of the particle; (c) Speed v; (d) Inter-particle force i F ; (e) Total kinetic energy 2 . . 1 2 K E Mv = .

e) total kinetic energy 2 .
oscillation cycle was completed in 475 steps, thus the period is 0.0475 msec.τ = In Figures 12(b)-(e), the variations, with respect to time t, of four quantities, are presented: b) position y of the particle.c) speed v. d) inter-particle force i F .

∑
thus the mean (rms) speed is 0 221.8 m sec 355.6 , v v = < = significantly smaller than the initial speed.Finally, the mean kinetic energy results on the diagram K.E.-t of Figure 12(e), for comparison.
In Figure13(a), a sketch of the above small container, of side L = 4.55 cm, is shown, with the particles of diameter D = 1.35 cm, in the initial state, with the speeds directed outwards.

Figure 13 . 2 .
Figure 13.Fourth application.Critical point.(a) The small container of side L = 4.55 cm with particles of diameter D = 1.35 cm, in the initial state.In the following diagrams, variations of five quantities with respect to time t; (b) Position y of the particle; (c) Speed v; (d) Inter-particle force i F ; (e) Particle-wall force w F ; (f) Total kinetic energy 2 . . 1 2 K E Mv = .

F
are noted on the diagrams F i − t, F w − t of Figure13(d) and Figure13(e), respectively.
on the diagram K.E.-t of Figure 13(f), for comparison.The above mean kinetic energy . .K E corresponds to a value of gas constant R = 2.856 Joules mole −1 •K −1 , as obtained by equating 1302 Joules 3 2 304 K. R = × 3 .And for initial temperature ranging from 0 50 C 223 K T = − =  up to 100 C 273 K =  , with a step 10 C 10 K T ∆ = =  , thus, for initial speed ranging from 0 355.6 m sec v = up to 459.8 m/sec, for every couple of L, T 0 , the corresponding pressure P and mean kinetic energy . .K E have been determined.Every point

Figure 15 .
Figure 15.Variation of gas constant R values in the Critical Point region of Carbon Dioxide, obtained by comparison of iso-. .K E curves of proposed model of Figure 14 with corresponding iso-therms of Eastman-Rollefson tests [14] [20] and those of Van der Waals model [14] [16].

Figure 16 .
Figure 16.Iso-therms of proposed model in Critical Point region of Carbon Dioxide, obtained from transformation of iso-. .K E curves of proposed model of Figure 14, by use of variation of gas constant R values of Figure 15.C.P.: Critical Point.

Figure 17 .
Figure 17.Comparison of iso-therms of proposed model to corresponding ones of Eastman-Rollefson tests [14] [20], in the Critical Point region of Carbon Dioxide.

Figure 18 .
Figure 18.Comparison of iso-therms of proposed model with corresponding ones of Van der Waals model [14] [16] in the Critical Point region of Carbon Dioxide.