Canada

Zero Emission Vehicles (ZEVs) are necessary to help reduce the emissions in the transportation sector which is responsible for 40% of overall greenhouse gas emissions. There are two types of ZEVs Battery Electric Vehicles (BEVs) and Fuel Cell Electric Vehicles (FCEVs) Commercial Success of BEVs has been challenging thus far also due to limited range and very long charging duration. FCEVs using H2 infrastructure with SAE J2601 and J2799 standards can be consistently fuelled in a safe manner fast and resulting in a range similar to conventional vehicles. Specifically fuelling with SAE J2601 with the SAE J2799 enables FCEVs to fill with hydrogen in 3-5 minutes and to achieve a high State of Charge (SOC) resulting in 300+ mile range without exceeding the safety storage limits. Standardized H2 therefore gives an advantage to the customer over electric charging. SAE created this H2 fuelling protocol based on modelling laboratory and field tests. These SAE standards enable the first generation of commercial FCEVs and H2 stations to achieve a customer acceptable fueling similar to today's experience. This report details the advantages of hydrogen and the validation of H2 fuelling for the SAE standards.

Buoyancy effects and nozzle geometry can have a significant impact on turbulent jet dispersion. This work was motivated by applications involving hydrogen. Using helium as an experimental proxy buoyant horizontal jets issuing from a round orifice on the side wall of a circular tube were analyzed experimentally using particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF) techniques simultaneously to provide instantaneous and time-averaged flow fields of velocity and concentration. Effects of buoyancy and asymmetry on the resulting flow structure were studied over a range of Reynolds numbers and gas densities. Significant differences were found between the centreline trajectory spreading rate and velocity decay of conventional horizontal round axisymmetric jets issuing through flat plates and the pipeline leak-representative jets considered in the present study. The realistic pipeline jets were always asymmetric and found to deflect about the jet axis in the near field. In the far field it was found that the realistic pipeline leak geometry causes buoyancy effects to dominate much sooner than expected compared to horizontal round jets issuing through flat plates.

In order to facilitate the application of hydrogen energy and ensure its safety the compressed hydrogen storage tank on board needs to be full of hydrogen gas within 3 minutes. Therefore to meet this requirement the effects of refueling parameters on the filling time need to be investigated urgently. For the purpose of solving this issue a novel analytical solution of filling time is obtained from a lumped parameter model in this paper. According to the equation of state for real gas and dimensionless numbers Nu and Re the function relationships between the filling time and the refueling parameters are presented. These parameters include initial temperature initial pressure inflow temperature final temperature and final pressure. These equations are used to fit the reference data the results of fitting show good agreement. Then the values of fitting parameters are further utilized so as to verify the validity of these formulas. We believe this study can contribute to control the hydrogen filling time and ensure the safety during fast filling process.

With the simple structure and quick refuelling process the compressed hydrogen storage system is currently widely used. However thermal effects during charging-discharging cycle may induce temperature change in storage tank which has significant impact on the performance of hydrogen storage and the safety of hydrogen storage tank. To address this issue we once propose a single zone lumped parameter model to obtain the analytical solution of hydrogen temperature and use the analytical solution to estimate the hydrogen temperature but the effect of the tank wall is ignored. For better description of the heat transfer characteristics of the tank wall a dual zone (hydrogen gas and tank wall) lumped parameter model will be considered for widely representation of the reference (experimental or simulated) data. Now we extend the single zone model to the dual zone model which uses two different temperatures for gas zone and wall zone. The dual zone model contains two coupled differential equations. To solve them and obtain the solution we use the method of decoupling the coupled differential equations and coupling the solutions of the decoupled differential equations. The steps of the method include: (1) Decoupling of coupled differential equations; (2) Solving decoupled differential equations; (3) Coupling of solutions of differential equations; (4) Solving coupled algebraic equations. Herein three cases are taken into consideration: constant inflow/outflow temperature variable inflow/outflow temperature and constant inflow temperature and variable outflow temperature. The corresponding approximate analytical solutions of hydrogen temperature and wall temperature can be obtained. The hydrogen pressure can be calculated from the hydrogen temperature and the hydrogen mass using the equation of state for ideal gas. Besides the two coupled differential equations can also be solved numerically and the simulated solution can also be obtained. This study will help to set up a formula based approach of refuelling protocol for gaseous hydrogen vehicles.

In this paper a vertical turbulent round jet of helium was studied numerically using the PHOENICS software package. The flow was assumed to be steady incompressible and turbulent. The jet discharge Froude number was 14000 and the turbulent Schmidt number was 0.7. The incompressible Reynolds average Navier-Stokes equations and helium transport equation expressed in 2-D axisymmetric domain were applied to model the underlying helium release. The k-e RNG turbulence model was used for the calculations of the corresponding turbulent viscosity diffusivity velocity and concentration fields in the domain. The simulation results are compared with the experimental measurements from the earlier published studies on helium jets in non-buoyant jet region (NBJ) intermediate region (I) and buoyant plume region (BP). The numerical results show that the radial profiles of mean velocity and mean concentration are consistent with the empirical data scaled by the effective diameter and density-ratio dependence. The mean velocity and concentration fields along the axis of the jet agree with the decay laws correlated from the previous experiments. The discrepancy between the numerical and experimental data is within 10% proving that the current CFD model for gas release and dispersion is robust accurate and reliable and that the CFD technique can be used as an alternative to the experiments with similar helium jets. The authors believe that the current CFD model is well validated through this study and can be further extended to predict similar hydrogen releases and dispersion if the model is properly applied with hydrogen properties.

Computational fluid dynamics is used to simulate the release of high pressure Hydrogen from a reservoir with an exit of increasing diameter. Abel-Noble real gas equation of state is used to accurately simulate this high pressure release. Parallel processing based on Message Passing Interface for domain decomposition is employed to decrease the solution time. The release exit boundary is increased in time to simulate a scenario when the exit area increases during the release. All nodes and elements are moved accordingly at each time step to maintain the quality of the mesh. Different speeds of increasing diameter are investigated to see the impact on this unsteady flow.

The effect of vertical surface on the extent of high pressure unignited jets of both hydrogen and methane is studied using computer fluid dynamics simulations performed with FLACS Hydrogen. Results for constant flow rate through a 6.35 mm round leak orifice from 100 barg 250 barg 400 barg 550 barg and 700 barg compressed gas systems are presented for vertical jets. To quantify the effect of the surface on the jet the jet exit is positioned at various distances from the surface ranging from 0.029 m to 12 m. Free jets simulations are performed for comparison purposes.

A buoyant round vertical hydrogen jet is investigated using Large Eddy Simulations at low Mach number (M = 0.3). The influence of the transient concentration fields on the extent of the gas envelope with concentrations within the flammability limits is analyzed and their structure are characterized. The transient flammable region has a complex structure that extends up to 30% beyond the time-averaged flammable volume with high concentration pockets that persist sufficiently long for potential ignition. Safety envelopes devised on the basis of simplified time-averaged simulations would need to include a correction factor that accounts for transient incursions of high flammability concentrations.

During the early phase of the transient process following a hydrogen leak into the atmosphere a contact surface appears separating air heated by the leading shock from hydrogen cooled by expansion. Locally the interface is approximately planar. Diffusion leads to a temperature decrease on the air side and an increase in the hydrogen-filled region and mass diffusion of hydrogen into air and of air into hydrogen potentially resulting in ignition. This process was analyzed by Li ˜nan and Crespo [1] for unity Lewis number and Li ˜nan and Williams [2] for Lewis number less than unity. We included in the analysis the effect of a slow expansion [3 4] leading to a slow drop in temperature which occurs in transient jets. Chemistry being very temperature-sensitive the reaction rate peaks close to the hot side of the interface where only a small fuel concentration present close to the warm air-rich side which depends crucially upon the fuel Lewis number. For Lewis number unity the fuel concentration due to diffusion is comparable to the rate of consumption by chemistry. If the Lewis number is less than unity diffusion brings in more fuel than temperature-controlled chemistry consumes. For a Lewis number greater than unity diffusion is not strong enough to bring in as much fuel as chemistry would burn; combustion is controlled by fuel diffusion. If the temperature drop due to expansion associated with the multidimensional jet does not lower significantly the reaction rate up to that point analysis shows that ignition in the jet takes place. For fuel Lewis number greater than unity chemistry does not lead to a defined explosion so that eventually expansion will affect the process; ignition does not take place [3 4]. In the current paper these results are extended to consider multistep chemical kinetics but for otherwise similar assumptions. High activation energy is no longer applicable. Instead results are obtained in the short time limit still as a perturbation superimposed to the self-similar solution to the chemically frozen diffusion solution. In that approximation the initiation step which consumes fuel and oxidant is taken to be slow compared with steps that consume one of the reactants and an intermediate species. The formulation leads to a two point boundary value problem for set of coupled rate equations plus an energy equation for perturbations. These equations are linear with variable co-effcients. The coupled problem is solved numerically using a split algorithm in which chemical reaction is solved for frozen diffusion while diffusion is solved for frozen chemistry. At each time step the still coupled linear problem is solved exactly by projecting onto the eigenmodes of the stiff matrix so that the solution is unaffected by stiffness. Since in the short time limit temperature is only affected at the perturbation level the matrix depends only on the similarity variable x t but it is otherwise time-independent. As a result determination of the eigenvalues and eigenvectors is only done once (using Maple) for the entire range of discretized values of the similarity variable. The diffusion problem consists of a set of independent equations for each species. Each of these is solved using orthogonal decomposition onto Hermite polynomials for the homogeneous part plus a particular solution proportional to time for the non-homogeneous (source) terms. That approach can be implemented for different kinetic schemes.

The heat transfer analysis in the filling process of compressed on-board hydrogen storage tank has been the focus of hydrogen storage research. The initial conditions mass flow rate and heat transfer coefficient have certain influence on the hydrogen filling performance. In this paper the effects of mass flow rate and heat transfer coefficient on hydrogen filling performance are mainly studied. A thermodynamic model of the compressed hydrogen storage tank was established by Matlab/Simulink. This 0D model is utilized to predict the hydrogen temperature hydrogen pressure tank wall temperature and SOC (State of Charge) during filling process. Comparing the simulated results with the experimental data the practicability of the model can be verified. The simulated results have certain meaning for improving the hydrogenation parameters in real filling process. And the model has a great significance to the study of hydrogen filling and purification.