Numerical simulation of pressure pulsation effects of a snubber in a CNG station for increasing measurement accuracy

Natural companies employ turbine flow meters to measure natural gas which delivered to Com- pressed Natural Gas stations. The stations utilize compressors to increase pressure. The compressor produces a flow pulsation, which affects the accuracy of the measurement. The main aim of this arti- cle is to decrease the compressor effects on measurement accuracy by utilizing a snubber between the turbine flow meter and the reciprocating compressors. For this aim, numerical modeling has been built to simulate natural gas flow through a snubber. The effects of various snubber parameters on pressure pulsation have been investigated. The parameters included snubber volume to the minimum volume ratio, the ratio of height to diameter, outlet pipe length, and the existence and non-existence a buffer. The Ansys Fluent has been used for numerical modeling with transient analysis. Results show that in H/D value of 3, the maximum reduction in the percentage of pressure pulsation drop is about 47% and increasing the outlet pipe length to the 10 times of initial length causes a decrease of about 83% in pressure pulsations. Besides, for the ratio of snubber volume to the minimum volume from 1 to 16.7, the amplitude of pressure pulsations decreases from 4.1% to 0.25%.


Introduction
The utilization of NG (Natural Gas) as a fuel for vehicles has been underway for decades. One of the forms of NG used all over the world is CNG (Compressed Natural Gas). CNG technology has been developed and applied for decades and is becoming an established applicable technology. Methane is the main element (97%) of Natural Gas. This gas has no color or odor and is lighter compared with air. With ever growing concerns about environmental pollution, the high cost of petroleum fuels, the low cost of natural gas compared to conventional fuels and the low cost of the process for preparing natural gas for industrial application have led to a strong desire for different countries to use natural gas as a fuel for automobiles. If Natural Gas is properly used as a clean fuel in cars, the amount of harmful pollutant emissions will be greatly reduced. Thus development of CNG stations will facilitate the usage of Natural Gas and will lead to lower environmental pollution. The most important component of a CNG station is a reciprocating compressor. In a CNG station, the pressure of the gas must be increased from approximately 0.5 MPa to about 20-25 MPa. A multi-stage compressor is employed to achieve this pressure increment. The largest share of the current and initial costs of the stations depends on the reciprocating compressor (Farzaneh-Gord, Niazmand, Deymi-Dashtebayaz, & Rahbari, 2015).
In reciprocating compressors, pulsations can occur due to the opening and closing of the valves in the inlet and outlet flow lines. If vibrations do not reduce to a standard level, they may cause various issues. In general, these vibrations should be eliminated for reasons of providing a laminar flow into and out of the compressor, preventing excessive loading on compressors and reducing overall compressor vibrations (Yu, Na, & Kim, 2002). Mujic, Kovacevic, Stosic, and Smith (2008) stated that the pulsations of gas in discharge and suction chambers are the main reasons for noise creation in screw compressors. The effects of compressor working conditions and the geometric features of the mentioned phenomena were analyzed in the study. An area function was distinguished as the crucial factor affecting the pulsations of the gas. In addition, it was concluded that the amplitude could be decreased by port shape optimization. Wu, Xing, Peng, and Shu (2004) investigated unsteady equations with periodic flow changes at the discharge port of a twin-screw compressor. In their study, one-dimensional unsteady equations, by considering the impacts of gas and pipe heat transfer in addition to friction for gas flow, was derived to explain the pulsation of pressure in discharge mode. The two-step Lax-Wendroff scheme was used in order to obtain the results of the problem. They concluded that difference in pressure between upstream of the discharge port and discharge line some way downstream has the key role in the pulsation of discharge pressure. The lowest fluctuation happened in the cases where there was the same value for discharge pressure and the pressure belonged to the built-in volume ratio. In the case of constant discharge pressure, the amplitude of pressure pulsation increased with the rotational speed. Furthermore, accurate natural gas measurement is an important issue, especially in residential natural gas consumption usage levels (Farzaneh-Gord, Parvizi, Arabkoohsar, Machado, & Koury, 2015). In CNG stations, Turbine gas flow meters are used extensively to measure the volume of NG entering the site. These types of flow meters are very sensitive to oscillating flows. When the reciprocating compressor is operated, due to the operation of the piston and the valves, a pulse is produced at the suction line and downstream of the flow meter. It results in false pulses in flow meters and, ultimately, errors in measurement.
Experimental and theoretical research has been performed to determine the pulsation effects on turbine flow meters' measurement accuracy. Atkinson (1992) solved the turbine rotor motion equation and employed the magnetic pickup registering the passing of a rotor blade to compute the real volume flow. This approach should only be employed if the amplitude of the pulsations is detectable by the turbine signal. By the rapid reduction of the amplitude of pulsation in the turbine meter signal due to an increase in frequency, it was difficult to forecast the real flow in the cases of pulsations with high frequency. Lee, Cheesewright, and Clark (2004) developed another software tool called the 'Watchdog System'. In addition, this system utilized the equation of the motion of the rotor, but an accelerometer is also used to measure the flow noise of a valve or bend. Assuming incompressible and quasi-steady flow and negligible friction force, a relationship can be obtained between the pulsations of velocity and the error of measurement (Lee et al., 2004). Experiments have been conducted in the past by Lee et al. (2004), Jungowski and Weiss (1996) and McKee (1992). Vetter and Schweinfurther (1987), Zhao, Xiao, Wang, Luo, and Cao (2018) and Hou, Zhang, Li, and Zhang (2017) investigated pressure pulsation in a reciprocating pump. Vetter and Schweinfurther (1987) presented a calculation approach to model precisely different pump installations. Comparison of calculated and experimental data shows a good agreement and provides a validation of the computational model. Notzon (1996) discussed different methods of pulsation dampening in high pressure reciprocating pumps and described ways to control and optimize the effectiveness of damping measures by numerical pulsation analysis. Reciprocating pumps are often used as pressure generators in high-pressure process technology. Modern numerical simulation methods are able to determine pressure pulsations in high pressure plants with high accuracy. Therefore, a prediction of pressure pulsations in the planning stage becomes possible and measures for improvements can be taken in time. Commonly used simulation programs are based on a numerical solution of the one-dimensional classical fluid mechanics equation by the method of characteristics. Zhou, Liu, Wang, Liu, and Zhao (2017) numerically investigated an unsteady multiphase flow in the whole passage of a Francis hydraulic turbine. In their research, pulsation of the pressure was forecast and compared with the data experimentally obtained. Afterwards, the relationship between the pulsation of pressure and air admission was analyzed. According to the numerical data, air admission from a spindle hole reduced the difference in the pressure in the draft tube horizontal part, which in turn reduced the amplitude of low-frequency pulsation of pressure in the draft tube; the nitration of rotor-stator between the runner and air inlet resulted in an increment in the pressure pulsation of the blade-frequency in front of the runner. Jeong et al. (2008) explained the numerical investigation of gas flow moving via a snubber in a hydrogen gas compressor unit. In their study, 30 snubber models were generated by changing the height of the snubber angle of buffer. Afterwards, the models were simulated with the real working operation of a hydrogen gas compressing system. A Star CD, which had transient analysis and used a k-ε / high Reynolds number transient model, was used as a Computational Fluid Dynamic (CFD) code package. The research was carried out by comparing pressure loss and the pulsation of pressure, since the factors noted are the objective functions of optimizing the snubber. Numerical outcomes showed that pressure loss increases with the increment in the volume of the snubber. To the contrary, however, the pressure pulsation is decreased. By determining the angle of buffer as the adjusted variable, the lowest pressure loss was obtained at 30°; however, the trend of the pulsation of pressure was escalating. Akbar, Shim, and Yi (2006) and Akbar, Shim, Yi, Chung, and Jeong (2006) experimentally investigated the decrease in pressure pulsation in the compressing system using a snubber with and without buffer. The experiment measured pressure at the inlet and outlet of the snubber. It used an air compressor as a model of a hydrogen one. Snubbers with and without a buffer were utilized in order to reach comprehensive comparison between both of the systems noted. An analysis by employing a Fast Fourier Transform (FFT) approach was carried out to verify the frequency of operating pressure. They observed that pulsation was increasing proportionally in the cases where there was an increment in the driver motor frequency.
During investigating different aspects of technology, CFD has many engineering and medical applications (Alizadeh et al., 2018;Ghalandari, Mirzadeh Koohshahi, Mohamadian, Shamshirband, & Chau, 2019;Ramezanizadeh, Alhuyi Nazari, Ahmadi, & Açıkkalp, 2018). CFD-based simulation of utilizing natural gas in a dual-fueled diesel engine was done on a dual-fueled constant-speed engine (Akbarian et al., 2018). The outcomes of the study revealed that under different loads and pilot to gaseous fuel ratios, the nitrous oxides and particle materials emissions in the dual mode were lower than the corresponding values of the diesel engine. Faizollahzadeh Ardabili et al. (2018) reviewed a computational intelligence method for modeling hydrogen production. Jiang (2002, 2004) developed a 3-D numerical model for pollutant transport, which was on the basis of orthogonal curvilinear and sigma coordinate systems in horizontal and vertical directions, respectively. A simple open boundary condition was considered in the model noted. Furthermore, Ramezanizadeh, Alhuyi Nazari, Ahmadi, and Chau (2019) numerically investigated a thermosyphon-based heat exchanger and compared its performance with other kinds of heat exchangers.
The aim of the pulsation analysis is to prevent difficulties such as: an unfavorable level of fluctuations and vibrations and cyclic pipe stresses leading to fatigue failure; increased consumption of power and loss of capacity; hammering of compressor valves, relief valves and check valves due to oscillations; and inaccuracies of devices used for flow measurement. Most of the studies in this area have been carried out to investigate the effect of snubber applications in the output of the compressor; a few studies have been conducted to investigate the effect of snubber application in the suction line of the compressor. The aim of this article is to propose a snubber at the suction line of the compressor of a CNG station to decrease pressure pulsation, and finally eliminate or decrease compressor effects on the measuring of the NG flow rate. To achieve this aim, 3D numerical investigation is conducted in ANYSYS FLUENT for a snubber between a turbine flow meter and a reciprocating compressor. This investigation studies the effect of buffer application on pressure specifications. Significant pressure specifications such as pressure drop and pressure pulsation are also studied in order to optimize the snubber. It should be noted that the best snubber has minimum pressure loss and pressure pulsation.

CNG station equipment and the application of a snubber
CNG stations receive fuel via a local utility line at a pressure lower than that used for vehicle fueling. The gas is compressed to a higher pressure for vehicle fueling. According to the storage capacity, the size of the compressors and dispensing rate, CNG stations are of three types: fast-fill; time-fill; and combination-fill. The CNG fueling equipment includes compressors, a dryer, storage, dispensers and flow meter. As mentioned, one of the problems of such a fueling station is the error occurring in measuring the gas flow due to the oscillating flows. To eliminate or decrease these pulsations, the solution is the application of a snubber in these stations. Figure 1 schematically shows a snubber. To improve the performance of the snubber, a metal plate called a buffer is placed inside the snubber. It is also shown how to connect this buffer to the snubber.

Minimum volume of the snubber
In this project, the snubber volume of a CNG station is investigated for 250 psi (17.1 bar) inlet gas pressure. The minimum volume of a snubber in a suction line of a compressor is calculated according to API standard 618. This standard covers the lowest needs for reciprocating  compressors and their drivers for utilization in the gas, petroleum and chemical industries. One of the important subjects of the standard is pulsation analysis. This standard proposes a relation for a minimum volume of a snubber in the suction line of a compressor as below: In Equation (1), V s , γ , T s , M and PD are the minimum volume of a snubber (m 3 ), compressibility factor, gas temperature at suction line of compressor (K), molar mass and displacement volume after one revolution of the crankshaft. The maximum allowed pulsation of a snubber is calculated as below: In Equation (3), P L is mean pressure of line. Thus according to the line pressure of 17.1 bar, the maximum allowed pulsation equals 1.6% of line pressure. The minimum volume is presented in Table 1.

Calculation of drop in pressure pulsations
The main objective of this article is to decrease the pressure pulsations at the inlet of the snubber. Thus three equations are defined as below (Jeong et al., 2008): In these equations, ΔP, P in , P out , P max and P min , are the pressure difference between the inlet and outlet of the snubber, pressure at the inlet of the snubber, pressure at the outlet of the snubber, and the maximum and minimum pressure of the cycle, respectively. In addition, A, A red , A in and A out are the percentage of pressure pulsation amplitude related to the line pressure, percentage of pressure pulsation drop, pressure pulsation amplitude of inlet and pressure pulsation amplitude of outlet. Decrease in pulsations means a decrease in A (Akbar, Shim, Yi, Lee, et al., 2006;Jeong et al., 2008).

Governing equation
One of the most significant choices in CFD is the verdict of an applicable turbulence model. The main aim of any turbulence model is providing an approach for determining the impact of turbulence fluctuations on the mean flow field. In the present study, the k-ε model has been implemented. The k-ε model is the most conventional and validated turbulence model, applicable for different flow fields, which indicates its appropriateness. It is generally applicable for free-shear layer flows with relatively average pressure gradients in addition to the confined flows where the Reynolds shear stresses are most important. A differential form of continuity equation is as follows (Jeong et al., 2008): Generally, conservative equations of turbulent fluid flow are written as: A conservative equation of Energy is written as (Jeong et al., 2008): λ is Thermal Conductivity in Equation (9). These equations, called RANS, are nearly identical to the original Navier-Stokes equations written for the mean flow variables. These equations are explicit and no assumptions have been applied to achieve them. Stress tensor is defined by Equation (10) and components of −ρu i u j are known as Reynolds stress (Jeong et al., 2008).

Geometry, modeling and grid system
Problem geometry includes an inlet pipe of 2 in diameter and 1 m length, and an outlet pipe of 2 in diameter and different lengths. In addition, a cylinder exists between these two pipes. We study this cylinder in two cases, with and without a buffer. In order to find the suitable snubber, many geometries and many parameters are studied including the effects of snubber volume to the minimum volume ratio, the ratio of height to diameter (H/D), the outlet pipe length and usage or not to use a buffer in the snubber. Figure 2 shows the geometry of the equipment and Figure 3 presents the geometric parameters and the location of symbols in the geometry. In Figure 3, plane 1 is observed at the inlet of snubber cylinder. To measure pressure pulsation, this plate is located at a distance of 10 cm from the cylinder wall. Then it is required to generate a suitable mesh. High mesh quality is necessary to solve the problem. Course  mesh produces large numerical errors in the solution while very fine mesh can be unacceptably expensive. Thus finding the right mesh that meets the accuracy and be run in an acceptable time is very important. To check the mesh quality, the parameter of Skewness is used. The value of this parameter for each cell is between 0 and 1. A smaller value of Skewness presents a higher mesh quality. For 3D geometries, it should not exceed 0.75. Figure 4 depicts the mesh of geometry for two cases, with and without a buffer. The mesh includes the structured and unstructured meshes. Far from the wall, unstructured mesh with pyramid shaped elements and near the wall of the pipe, cylinder and buffer, hexahedral structured mesh is considered. To decrease mesh number and iteration, symmetry of xy plane is used in the modeling of the problem.

Boundary conditions
The boundary conditions include constant velocity at the inlet, pressure pulsation at the outlet, symmetry at xy plane for the inlet and outlet pipes and cylinder wall. Based on the type of CNG station with an inlet pressure of 17.1 bar and according to Table 2, the volumetric flow rate under normal condition is 950 m 3 /h by the measure of a turbine flow meter. According to the inlet pipe diameter of 2 in, the equality of mass flow rate at normal condition and real condition, inlet velocity is considered to be 8.42m/s. It should be noted that under normal condition, pressure is 1 atm and temperature is 15°C; under real condition, the pressure is 17.1 bar and temperature is  25°C. The values of velocity and density of the methane under real and normal conditions are presented in Table 2.
The pressure outlet boundary condition is selected for an outlet boundary which is connected to the compressor inlet or the source of pressure pulsations. To enter the amount of pressure pulsations as sinusoidal functions of time, UDFs are attached to FLUENT ( Figure 5). It should be noted that the base pressure is 250 psi (17.1 bar). The problem is solved for three pressure pulsations: 5, 10 and 25 psi (equivalents to 2%, 4% and 10% of line pressure). Also according to the rotational speed of the compressor crankshaft which is equal to 1200rpm, each compressor cycle occurs in 0.05 s. Figure 5 shows boundary condition at the outlet of snubber for amplitude of pressure pulsation of 4% of the line pressure. The boundary condition of inner cylinder wall and pipes is defined as a rigid wall with no-slip conditions. It means that the adjacent gas has no velocity. This boundary condition is employed for continuity, momentum, turbulence kinetic energy and turbulence dissipation rate.

Solver settings
Fluent settings, Model Solver and the solution method used in the simulation are presented in Table 3. A second order implicit unsteady solver was used in the simulation.

Verification of the computational method
For mesh dependency, seven mesh grids of 1383, 5676, 10,535, 50,580, 106,032, 222,234 and 417,362 are employed. In this study for a case with no buffer, the percentage of pressure pulsation amplitude of geometry is studied for minimum volume, output length of 1 m and height to diameter ratio of 2. As Figure 6 depicts, the variation of the percentage of pressure pulsation amplitude is   Figure 7 shows comparison between the results of this paper and the experimental study of Akbar, Shim, Yi, Chung, et al. (2006) and CFD work of Jeong et al. (2008). In this figure, the same trend between CFD and the experimental results is observed. The mean pressures of each case were compared for three motor frequencies. Good agreement can be seen with a maximum of only 0.43% of the difference.

Results and discussion
In the current study, the effect of different geometrical parameters on inlet pressure pulsation have been   investigated. The effects of these parameters are discussed in this section.

Height to diameter ratio
API standard 618 describes a relationship for the minimum volume of a snubber, while it does not express a relationship for the ratio of height to diameter of the cylinder. One of the effective parameters in determining the optimum design of the snubber is the ratio of height to diameter of the cylinder (H/D). For H/D values of 1, 2, 3 and 4, study is performed at the minimum snubber volume, pressure pulsation of 10 psi (4% of line pressure), and exit length of 1 m with no buffer. Height, diameter and the ratio of height to diameter of the cylinder are presented in Table 4. For H/D values of 1, 2, 3 and 4, study is performed at the minimum constant snubber volume, pressure pulsation of 10 psi (4% of line pressure), and outlet pipe length of 1 m with no buffer. In API Standard 618, it is stated that for a snubber with a single cylinder, the ratio of height to diameter should not exceed more than 4. Thus the maximum value of H/D is set to 4. Figures 8 and 9 state that the maximum reduction in percentage of pressure pulsation drop is about 47%, which takes place with a H/D value of 3.
For plane 1, inlet and outlet of the snubber, Figures 10 and 11 present pressure and velocity variation resulting from pressure pulsations of the outlet within 0.3 s of the compressor start for H/D = 3. Time phase difference between the inlet and outlet of the snubber is produced due to the fact that time is running out when pulsations reach the inlet and outlet of the snubber. As is clear, the pressure pulsations decrease with the presence of the snubber. In Figure 11, it is seen that the velocity is high in some moments due to the high pressure difference between the inlet and outlet of the snubber. In some moments the velocity approximates to 0 and this is while the pressure difference between inlet and outlet is insignificant. According to Figure 10, when pressure at the inlet of the snubber is at the maximum value of 1770 kPa, pressure at the outlet is at the minimum value of 116 kPa. Based on these results, from the compressor start-up (from the first cycle to the second cycle), the peak pressure pulsation decreases; after the second cycle, we witness the constant pulsations and equality of the minimum and maximum values of pressure in each cycle. In comparison with outlet velocity amplitude, smaller velocity amplitude is produced at the inlet.
For six working cycles of a compressor, Figures 12 and  13 indicate the variation of pressure and velocity on plane 1, located 10 cm from the snubber inlet and 0.3 s from the start of pulsation at the outlet. Excepting the H/D value of 1, for other ratios, there is little difference between the pulsations for both pressure and velocity. Figures 14 and 15 present contours of pressure and velocity for different values of H/D. These contours are for 0.3 s. Figure 14 shows comparison between H/D values of 1, 2, 3 and 4. Pressure difference between the   inlet and outlet of the snubber decreases while for higher values of H/D, no difference is observed. According to Figure 15, the same trend occurs for velocity difference. It means velocity difference decreases for H/D value of 1 and 2 while for H/D values of 2, 3 and 4, no remarkable velocity difference is observed in the velocity field.

The length of outlet pipe
The other parameter studied here is the outlet pipe length, L out , the distance of the connecting pipe between the snubber and the compressor. The lengths of 0.2, 0.5, 1.0, 2.0 and 5.0 m are considered. This case is done at minimum snubber volume, pressure pulsation of 10 psi (4% of line pressure) and H/D value of 2.5 with no buffer. The length of snubber outlet pipe, in fact, is the distance between the snubber and the compressor or the source of pulsations. The results were obtained for a snubber with a volume equal to the minimum volume calculated by standard API 618. The length of inlet pipe is 1 m, H/D is 2.5 and amplitude of pressure pulsation is 4% of line pressure. The aim of this section is to prove that by increasing the length of the snubber outlet pipe, the pressure pulsation effects on the flow meter can be decreased. The length of the snubber outlet pipe is designed according to the space of CNG station. While Figures 16 and 17      show that the greater the distance from the compressor, the lower the effect of pressure pulsations of compressor on the pressure pulsations of the snubber inlet, the amplitude of pressure pulsations of 3.2% of line pressure for outlet pipe length of 0.5 m reaches to 0.7% for outlet pipe of 5 m. It means that the drop in pressure pulsations will be about 83% for an outlet pipe of 5 m. The reason is the pressure drop, which damps the impact of these  pulsations by increasing the length of the pipe. This result also states that with a small snubber volume equal to the minimum volume with increasing distance from the compressor, pulsations can be significantly decreased. It should be noted that the best place to install a snubber at the station is immediately after the measurement unit and before the dryer of the compressor.
For outlet pipe length of 2 m, Figures 18 and 19 show the effects of pressure pulsation and velocity pulsation of the outlet on the pressure pulsation and velocity pulsation of inlet for plane 1 and in 0.4 s from the compressor start-up. At first cycles, the decreasing trend of the pressure pulsation at the inlet shows the effect of the application of the snubber. From the fourth cycle, the amplitude of pressure and velocity pulsations at the inlet and outlet reaches constant values, the maximum pressure of 1745.5 kPa and the minimum pressure of 1704.7 kPa. According to Figure 19, a maximum velocity of 9.25 m/s and a minimum velocity of 7.32 m/s are achieved at the inlet of the snubber. Thus amplitude of velocity pulsations is 0.9 m/s.
To compare the pressure pulsation and velocity pulsation over time at the inlet of the snubber for different values of outlet pipe length, Figures 20 and 21 are presented. These figures clearly show a decreasing trend of pulsation amplitude with increasing the outlet pipe length in 0.3 s from the compressor start-up; as output pipe length increases, the pulsation amplitude also decreases. The  results are not in the same time phase because with an increase in outlet length, it takes more time to affect pulsations on the input.

The ratio of current volume to the minimum volume (CV)
The other parameter playing an important role is the ratio of snubber volume to the minimum volume (CV = V/V s ). A few cases are studied with CV ratios of 1, 2, 5, 10 and 16.7 and pressure pulsations of 5, 10 and 25 psi (2%, 4% and 10% of line pressure). These cases are done at the outlet pipe length of 1 m, and an H/D value of 2.5 with no buffer. Table 4 presents the ratio of height to diameter of the cylinder. Table 5 presents the geometric dimensions of a snubber cylinder for various CV ratios. Figure 22 shows the percentage of pressure pulsation amplitude. Increasing the snubber volume effectively controls pressure pulsations, so that for amplitude of pressure, pulsations decrease from 4.1% to 0.25% at the inlet for CV ratios of 1-16.7. According to Equation (3), the maximum allowed value for the pressure pulsations is 1.6% of line pressure.
Percentage of pressure pulsation drop resulting from the variation of snubber volume is shown in Figure 23. It is clear that for CV ratios greater than 5, 90% of pulsations     For different values of CV, Figure 25 shows that velocity as well as pressure has constant pulsations from the second cycle. In CV = 1, the maximum and minimum velocities are 10.3 and 6.2 m/s with an amplitude about 2 m/s, while in CV = 16.7, the maximum and minimum velocities are 8.43 and 8.14 m/s with an amplitude about 0.15 m/s. Figure 26 shows average velocity for pulsation amplitude of 2%, 4% and 10% of the line pressure. Here line pressure more than 10% is not considered, because  for these cases, the produced velocity at the outlet is very high and, in accordance with the Natural Gas standards, gas velocity in the pipelines should not exceed 20 m/s. For line pressure of 10%, the average velocity is about 68 m/s and this high velocity damage severely the valves and fittings. One way to slow down the velocity at the outlet is an increase in the distance between the snubber and the compressor or an increase in the outlet pipe diameter.

Buffer existence
Installation of a buffer or a simple blade in the snubber is a proposal to decrease snubber volume. The effect of buffer use in the snubber on pressure pulsations is studied. This parameter is examined for an H/D value of 2.5, amplitude of pressure pulsations 4% of the pressure line, outlet pipe length of 1 m, 4 buffer angles of 20°, 30°, 40°and 50°and 3 CV ratios of 2.5 and 10.  The geometric dimensions of the buffer are presented in Table 6.
As a proposal to reduce the snubber volume, the presence of a buffer or a simple blade inside the snubber is investigated. This parameter is examined for four angles (20°, 30°, 40°and 50°), the CV values of 2, 5 and 10, H/D value of 2.5, pulsation amplitude of 4% line pressure and an outlet pipe length of 1 m. Then a comparison with a no buffer case was performed. The application of a buffer in the snubber is provided by Jeong et al. (2008). To decrease pressure pulsations in the outlet of the compressor, they   used blades and achieved good results in damping pulsations. But as shown in Figures 27 and 28, when a snubber is installed at the compressor inlet, the utilization of a buffer with different angles or the lack of a buffer have little effect in reducing amplitude, especially for CV = 10.
For CV = 5, Figure 29 shows the variation of pressure pulsations at the inlet with and without a buffer. As is clear, the variation of buffer angle has no effect in reducing the pulsation amplitude, but in the model without the buffer, the time phase is earlier compared to the model with a buffer. This is due to the rotation of the gas flow behind the buffer for a while. Figures 30-33 show the stream lines and the velocity vector in a cycle, for four moments. At 0.2 s, pulsations start at the outlet and increase gradually until 0.2125 and reach the maximum pressure value. In the outlet pipe, flow direction is against the direction of inlet flow.  Then pressure reaches the minimum value at the outlet at 0.2375 s and again the flow continues in the direction of inlet flow.

Conclusions
Turbine flow meters are extensively used in CNG fueling stations to measure the volumetric flow rate of the inlet NG. These types of flow meters are very sensitive to oscillating flows. When a reciprocating compressor starts operation, due to piston and valve performance, oscillating flows are created in the suction line of the compressor (downstream of turbine flow meter). These pulsations make errors in measurement and are the source of difference between the gas supplier and gas consumer. In this study, the numerical simulation of a snubber was investigated. The snubber is located after the turbine flow meter and before the compressor. The final target is to eliminate  or reduce the effects of the compressor on the turbine flow meter in a CNG station and consequently to reduce measurement errors. To achieve this goal and study the effects of different geometrical parameters, governing equations were supplemented for compressible fluid. The boundary conditions are a non-slip condition on the walls, constant velocity at inlet and pressure pulsation at outlet. Continuity, momentum and turbulent equations were solved by ANSYS FLUENT software. Numerical simulation of the snubber was carried out for a CNG station with a gas inlet pressure of 17.1 bar and inlet and outlet pipe diameter of 2 in. In this paper, several parameters, including ratio of volume to minimum volume, distance between snubber and compressor, ratio of height to diameter of snubber, and utilization of a buffer in the cylinder are investigated to achieve the best design for snubber dimensions. Comparing the results with previous studies shows that the difference is less than 2%, which indicates the validity of simulation.
Results showed that for the ratio of volume to a minimum volume of 5, the pressure pulsation amplitude would decrease less than the maximum amount limit of 1.6% of pressure line. Furthermore, increasing the distance between snubber and compressor has significant effect on decreasing the pressure pulsation while the application of a buffer inside the snubber has no effect on a decrease in pressure pulsation. The maximum decrease in the percentage of pressure pulsation drop takes place in a height to diameter ratio of 3. The maximum reduction in the percentage of pressure pulsation drop is about 47%, which takes place with a H/D value of 3. Pressure difference between the inlet and outlet of the snubber decreases, while for higher values of H/D, no difference is observed. The same trend occurs for velocity difference, which decreases for H/D values of 1 and 2, while for H/D values of 2, 3 and 4, no remarkable velocity difference is observed in the velocity field. Based on the results, from the first cycle to the second cycle, the peak pressure pulsation decreases; after the second cycle, we witness the constant pulsations. In comparison with an amplitude of velocity pulsation of the outlet, a smaller amplitude of velocity pulsation is produced at the inlet. The greater the distance from the compressor, the lower the effect of pressure pulsations of the compressor on the pressure pulsations of the snubber inlet. The amplitude of pressure pulsations are 3.2% and 0.7% of the line pressure for outlet pipe length of 0.5 and 5 m. This means that the drop in pressure pulsations will be about 83% for an outlet pipe of 5 m in comparison with an outlet pipe length of 0.5 m. In addition, increasing the ratio of snubber volume to the minimum volume effectively controls pressure pulsations, so that for amplitude of pressure, pulsations decrease from 4.1% to 0.25% at the inlet for CV ratios of 1-16.7.

A
Percentage of pressure pulsation amplitude related to the line pressure (%) A in Pressure pulsation amplitude of inlet (%) A inlet Area of inlet entrance (m2) A red Percentage of pressure pulsation drop (%) A out Pressure pulsation amplitude of outlet (%) M Molar mass (kg/kmol) P1 Maximum allowed pulsation of snubber (%) PD Displacement volume (m3/rev) P in Pressure at the inlet of snubber (KPa) P out Pressure at the outlet of snubber (KPa) P max Maximum pressure of the cycle (KPa) P min Minimum pressure of the cycle (KPa) ΔP Pressure difference between inlet and outlet of snubber (KPa) P L Mean pressure of line (KPa) Q N Inlet mass flow at normal condition (Kg/s) T s Gas temperature at suction line of compressor (K) V R Fluid velocity at the inlet (m/s) V s Minimum volume of snubber () (m3) ρ N Density of methane under normal condition (kg/m3) ρ R Density of methane under real condition (kg/m3) γ Compressibility factor (-)

Disclosure statement
No potential conflict of interest was reported by the authors.