Numerical Simulation on Fast Filling of Hydrogen for Composite Storage Cylinders

Author(s):  
Yanlei Liu ◽  
Jinyang Zheng ◽  
Ping Xu ◽  
Yongzhi Zhao ◽  
Lei Li ◽  
...  

Compressed gaseous hydrogen storage is widely used in hydrogen filling stations and fuel cell automobiles. As an effective filling is required to complete in a short time about 3 to 15 minutes, a high filling velocity is necessary which may lead to temperature rise and heterogeneity inside cylinders. In this research, a numerical method is proposed based on the finite volume method to investigate the thermal behavior such as temperature rise and distribution in the process of fast filling for the 15 L cylinder that designed by us. By numerical simulation, temperature, velocity and pressure distribution in the filling process are shown intuitively. Also, temperature and velocity gradients in the axial direction are obtained. In the simulation, parameters such as flow rates of fast filling are considered. Thus, curves of the temperature rise in the cases of different flow rates of fast filling are given. Furthermore, the limit of mass filling rate for the 15 L cylinder is obtained. The numerical results are compared with the existing experimental results and show high consistency.

2008 ◽  
Vol 33-37 ◽  
pp. 1303-1306
Author(s):  
Hajime Adachi ◽  
Kazuhiro Suga ◽  
Masanori Hayase ◽  
Shigeru Aoki

A numerical simulation of corrosion in a tube is performed with the solution velocity effect taken into account. A two dimensional tube, the cross-section of which is widening or narrowing with increase in distance, is considered. The velocity distribution in the tube is calculated with the Finite Volume Method (Open FOAM), and the derivatives of velocity with respect to the distance from the tube wall is determined at any location of the tube. The corrosion rate of the tube wall is estimated under the assumption that the corrosion rate depends on the velocity gradient, i.e. , it is estimated by solving the Laplace equation under the boundary conditions given with the polarization curves measured under various velocity gradients. The Boundary Element Method (3D-CAFE) is used to solve the Laplace equation. It is shown that the distribution of corrosion rate, including the maximum corrosion rate and its location, is different between the widening and narrowing tubes, even if the average velocities in the two tubes are equal.


2019 ◽  
Vol 14 (2) ◽  
pp. 82-88
Author(s):  
M.V. Alekseev ◽  
I.S. Vozhakov ◽  
S.I. Lezhnin

A numerical simulation of the process of the outflow of gas under pressure into a closed container partially filled with liquid was carried out. For comparative theoretical analysis, an asymptotic model was used with assumptions about the adiabaticity of the gas outflow process and the ideality of the liquid during the oscillatory one-dimensional motion of the liquid column. In this case, the motion of the liquid column and the evolution of pressure in the gas are determined by the equation of dynamics and the balance of enthalpy. Numerical simulation was performed in the OpenFOAM package using the fluid volume method (VOF method) and the standard k-e turbulence model. The evolution of the fields of volumetric gas content, velocity, and pressure during the flow of gas from the high-pressure chamber into a closed channel filled with liquid in the presence of a ”gas blanket“ at the upper end of the channel is obtained. It was shown that the dynamics of pulsations in the gas cavity that occurs when the gas flows into the closed region substantially depends on the physical properties of the liquid in the volume, especially the density. Numerical modeling showed that the injection of gas into water occurs in the form of a jet outflow of gas, and for the outflow into liquid lead, a gas slug is formed at the bottom of the channel. Satisfactory agreement was obtained between the numerical calculation and the calculation according to the asymptotic model for pressure pulsations in a gas projectile in liquid lead. For water, the results of calculations using the asymptotic model give a significant difference from the results of numerical calculations. In all cases, the velocity of the medium obtained by numerical simulation and when using the asymptotic model differ by an order of magnitude or more.


2021 ◽  
Vol 11 (11) ◽  
pp. 4990
Author(s):  
Boris Benderskiy ◽  
Peter Frankovský ◽  
Alena Chernova

This paper considers the issues of numerical modeling of nonstationary spatial gas dynamics in the pre-nozzle volume of the combustion chamber of a power plant with a cylindrical slot channel at the power plant of the mass supply surface. The numerical simulation for spatial objects is based on the solution conjugate problem of heat exchange by the control volume method in the open integrated platform for numerical simulation of continuum mechanics problems (openFoam). The calculation results for gas-dynamic and thermal processes in the power plant with a four-nozzle cover are presented. The analysis of gas-dynamic parameters and thermal flows near the nozzle cover, depending on the canal geometry, is given. The topological features of the flow structure and thermophysical parameters near the nozzle cap were studied. For the first time, the transformation of topological features of the flow structure in the pre-nozzle volume at changes in the mass channel’s geometry is revealed, described, and analyzed. The dependence of the Nusselt number in the central point of stagnation on the time of the power plants operation is revealed.


2010 ◽  
Vol 35 (15) ◽  
pp. 8092-8100 ◽  
Author(s):  
Lei Zhao ◽  
Yanlei Liu ◽  
Jian Yang ◽  
Yongzhi Zhao ◽  
Jinyang Zheng ◽  
...  

2014 ◽  
Vol 5 (1) ◽  
pp. 30-44 ◽  
Author(s):  
Qing-rui Meng

Purpose – The purpose of this paper is to reveal the temperature rise characteristics of the disc and pads under different load types. Design/methodology/approach – Evolutions of the disc and pads temperature under a stable, gradual changing and sine-wave contact pressures widely used at present are analyzed numerically by using ANSYS software. Findings – The results show that during the loading process, the temperature increases most rapidly under a stable contact pressure, most slowly under a gradual changing contact pressure; the disc temperature rise curves expose saw-shaped character, the closer it is to the friction surface, the more serious the fluctuations will be, the pads temperature rise curves are rather smooth; temperature gradient in the axial direction is higher than that in the other two directions under all of the three types of contact pressure and shows a sine-wave variation under a sine-wave contact pressure. Originality/value – It indicates that a gradual changing contact pressure should be adopted preferentially in practical application. The simulation results of this work provide theoretical basis for load simulation.


Author(s):  
Keyong Cheng ◽  
Xiulan Huai ◽  
Jun Cai ◽  
Zhixiong Guo

In the present study, numerical simulation is carried out for impingement/effusion cooling on the leading edge of a turbine blade similar to an experimental model tested previously. The k-ε turbulence model is used, and simulation parameters are set in accordance with the experimental conditions, including temperature ratio, blowing ratio, and Reynolds number of the main stream. The accuracy and reliability of the simulation is verified by the experimental data, and the influence of various factors on fluid flow and heat transfer is analyzed in detail. The results indicate that the blowing ratio is one critical factor which affects the cooling effectiveness. The greater the blowing ratio is, the higher the cooling effectiveness is. In addition, a staggered-holes arrangement is numerically studied and compared with a line-holes arrangement. The results show that the staggered-holes arrangement has a lower temperature on the outer surface of the leading edge and has improved the cooling effectiveness.


2015 ◽  
Vol 719-720 ◽  
pp. 46-49 ◽  
Author(s):  
Ginka Ranga Janardhana ◽  
Mani Senthil Kumar ◽  
B. Dhanasekar

The plasma cutting technology has been emerged as a developing technology which finds tremendous potential in fabrication and metal cutting industries. Thus for the cutting operation, the electrode inside the plasma torch plays a vital role for the plasma arc generation. The temperature of the arc is very high and at the electrode is around 3500°C. The cutting torch requires proper cooling system in order to prevent the electrode from quick wear due to the existence of high thermal gradient. The presented work aimed to study the impact of three coolants propylene glycol, ethylene glycol and de-ionized water flow over the electrode life. The experimental setups were arranged to study the heat transfer capabilities of the three coolants for different flow values and aimed to achieve the optimal flow rates for the efficient heat removal. The electrode life test trials were conducted to measure the electrode life for the flow values of three coolants in the temperature rise test. The optimal flow rates arrived from temperature rise test and the electrode life measured from life test are compared for the three coolant cases considered.


Author(s):  
Merouane Habib ◽  
Senouci Mohammed

In this paper, we investigate the no-reacting swirling flow by using the numerical simulation based to the unsteady Reynolds-averaged Navier-Stokes approach. The numerical simulation was realized by using a computational fluid dynamics CFD code. The governing equations are solved by using the finite volume method with two classical models of turbulence K-epsilon and Shear Stress K-ω. The objective of this paper is therefore to evaluate the performance of the two models in predicting the recirculation zones in a swirled turbulent flow. The current models are validated by comparing the numerical results of the axial, radial and tangential velocities to the experimental data from literature.


Sign in / Sign up

Export Citation Format

Share Document