G. Spazzafumo Department of Industrial Engineering, University of Cassino, Italy

Similar documents
A new cryogenic high-pressure H2 test area: First results

ME ENGINEERING THERMODYNAMICS UNIT III QUESTION BANK SVCET

Compact liquefied gas expander technological advances

Teknologi Pemrosesan Gas (TKK 564) Instructor: Dr. Istadi ( )

Analysis of a Supercritical Hydrogen Liquefaction Cycle


OUTCOME 2 TUTORIAL 2 STEADY FLOW PLANT

Transportation in a Greenhouse Gas Constrained World

Heat Effects in Gas Systems

Estimation of Boil-off-Gas BOG from Refrigerated Vessels in Liquefied Natural Gas Plant

FLOW & HEAT TRANSFER IN A PACKED BED - TRANSIENT

ENVIRONMENT-FRIENDLY HYDROGEN GAS AS FUEL IN FUEL CELL AND ITS CHALLENGES

Review Questions for the FE Examination

a. The power required to drive the compressor; b. The inlet and output pipe cross-sectional area. [Ans: kw, m 2 ] [3.34, R. K.

Hydrogen oxygen steam generator integrating with renewable energy resource for electricity generation

The Future of the Hydrogen Economy: Bright or Bleak?

Modified Reverse-Brayton Cycles for Efficient Liquefaction of Natural Gas

Keywords: GTL, Fisher-Tropsch, natural gas, synthesis gas, steam reforming NCPO, auto-thermal reforming. UNESCO EOLSS

OVERALL EFFICIENCY CONSIDERATION OF PNEUMATIC SYSTEMS INCLUDING COMPRESSOR, DRYER, PIPE AND ACTUATOR

A Hydrogen Economy. Dr. Mazen Abualtayef. Environmental Engineering Department. Islamic University of Gaza, Palestine

Exergy in Processes. Flows and Destruction of Exergy

CONTROL VOLUME ANALYSIS USING ENERGY. By Ertanto Vetra

The DAΦNE Cryogenic System

Performance of a counterflow heat exchanger with heat loss through the wall at the cold end

HT A COMPUTATIONAL MODEL OF A PHASE CHANGE MATERIAL HEAT EXCHANGER IN A VAPOR COMPRESSION SYSTEM WITH A LARGE PULSED HEAT LOAD

STUDY OF CRYOGENIC CYCLES WITH ASPEN - HYSYS SIMULATIONS

LNG UNIT (ENGINEERING DESIGN GUIDELINE)

Waste Heat Recovery for Fuel Cell Electric Vehicle with Thermochemical Energy Storage

Article Thermodynamic Analysis of Three Compressed Air Energy Storage Systems: Conventional, Adiabatic, and Hydrogen-Fueled

Process HEAT PROGRESS REPORT. Lauren Ayers Sarah Laderman Aditi Verma Anonymous student

Thermodynamic analysis on post combustion CO 2 capture of natural gas fired power plant

Michigan State University DEPARTMENT OF CHEMICAL ENGINEERING AND MATERIALS SCIENCE. ChE 321: Thermodynamics Spring 2017

Low-Grade Waste Heat Recovery for Power Production using an Absorption-Rankine Cycle

Magnetically Coupled Submerged Cryogenic Pumps and Expanders for Ammonia Applications

ENVIRONMENT AND DEVELOPMENT-Vol.II- Capacity For Developing Innovative Production Technologies - Ying Chen

llilliiiillllllllilllllllllllllllllllllllil CERN AT, (CR, 94'36

Hydrogen Workshop for Fleet Operators

Fluid Mechanics, Heat Transfer, Thermodynamics. Design Project. Production of Ammonia

Smart Integration of Thermal Management Systems for Electronics Cooling

Session 4.1: Solid Storage Technology. Dr. N. Eigen (GKSS) 25 th 29 th September 2006 Ingolstadt. Session 1.2: Introductory Lectures. K.

A Study on Re-liquefaction Process of Boil-off Gas of LCO 2 Transfer Ship

ENERGY CARRIERS AND CONVERSION SYSTEMS Vol. II - Molten Carbonate Fuel Cells - Kouichi Takizawa

Chapter 9: Applications of the Laws of Thermodynamics

A new re-liquefaction system of MRS-F (Methane Refrigeration System Full re-liquefaction) for LNG Carriers

To Hydrogen or not to Hydrogen. Potential as a Ship Fuel. Dr. John Emmanuel Kokarakis. Emmanuel John Kokarakis University of Crete

Qualitative Phase Behavior and Vapor Liquid Equilibrium Core

NUCLEAR ENERGY MATERIALS AND REACTORS - Vol. II - Advanced Gas Cooled Reactors - Tim McKeen

Gas Generator ULTRA HIGH PURE - HIGH FLOW. Gas Generator

DESIGN ANALYSIS OF A REFRIGERATED WAREHOUSE USING LNG COLD ENERGY

SHRI RAMSWAROOP MEMORIAL COLLEGE OF ENGG. & MANAGEMENT B.Tech. [SEM IV (ME-41, 42,43 & 44)] QUIZ TEST-1 (Session: )

A MODEL FOR RESIDUAL STRESS AND PART WARPAGE PREDICTION IN MATERIAL EXTRUSION WITH APPLICATION TO POLYPROPYLENE. Atlanta, GA 30332

CONDITION ASSESSMENT AND LIFE EXTENSION. R.A. Chaplin Department of Chemical Engineering, University of New Brunswick, Canada

Revue des Energies Renouvelables Spécial ICT3-MENA Bou Ismail (2015) Numerical study of a single effect ejector-absorption cooling system

Thermodynamic Data. CO (g, 0 C, 1 atm) CO (g,100 C, 1 atm):

GTI Small-Scale Liquefier Technology. March 2013

Combined Mass and Energy Transients

Modeling of HTPEM Fuel Cell Start-Up Process by Using Comsol Multiphysics

FIREX Foam Cryogenic Target Development - Attempt of Residual Voids Reduction with Solid Hydrogen Refractive Index Measurement -

S.E. (Chemical) (First Semester) EXAMINATION, 2012 PROCESS CALCULATIONS (2008 PATTERN) Time : Three Hours Maximum Marks : 100

Heat Pump Efficiencies simulated in Aspen HYSYS and Aspen Plus

PRISM Membrane Systems for petrochemical applications... tell me more

Available online at ScienceDirect. Energy Procedia 75 (2015 )

ENGINEERING INFORMATION Hot water and steam service

POLLUTION FROM MOTOR VEHICLES

Kinetics of AlH 3 decomposition and subsequent Al oxidation / Hydrogen Jetfires

Improvement of distillation column efficiency by integration with organic Rankine power generation cycle. Introduction

Induction Heating. Jean Callebaut, Laborelec. 1 Introduction Physical principles Induction Installations... 5

Chapter 7 Mass Transfer

COLD NEUTRON SOURCE AT CMRR

Innovative Design of a Compact Reformer for PEMFC

Exergy Analysis of a Power Plant in Abu Dhabi (UAE)

ADVANCED PROCESS CONTROL QATAR GAS ONE YEAR EXPERIENCE

ELECTROTECHNOLOGIES FOR ENERGY END-USES: APPLICATION OF MVR TO FOOD INDUSTRIES

ENERGY RECOVERY IMPROVEMENT USING ORGANIC RANKINE CYCLE AT COVANTA S HAVERHILL FACILITY

Thermal Management of Densely-packed EV Battery Set

Reforming Natural Gas for CO 2 pre-combustion capture in Combined Cycle power plant

Fluid Mechanics, Heat Transfer, Fluid Mechanics Design Project. Production of Ethanol

Source Characterization of Ammonia Accidental Releases for Various Storage and Process Conditions

The Reduced Pumping Power Requirements from Increasing the Injection Well Fluid Density

Liquid-Flooded Ericsson Power Cycle

Linde Kryotechnik AG. Effects of oil contamination on cryogenic plants. Ionic Liquids.

Investigation on a three-stage hydrogen thermal compressor based on metal hydrides

OPTIMIZING A STEAM-METHANE REFORMER FOR HYDROGEN PRODUCTION USING A LUMPED PARAMETER MODEL

Secondary Systems: Steam System

INSULATION AGING & GREEN FOAM

Chem 355 Jasperse DISTILLATION

THE INFLUENCE OF THE AMOUNT OF FIBRE REINFORCEMENT ON THE AIR PERMEABILITY OF HIGH PERFORMANCE CONCRETE

ENERGY CARRIERS AND CONVERSION SYSTEMS Vol. I - Hydrogen Separation and Handling - Itsuki Uehara

ME 343 Exam 2 November 24, 2014

IAPWS Certified Research Need ICRN. Thermodynamic Properties of Humid Gases and CO 2 -Rich Mixtures

NOVEL SCHEME FOR SMALL SCALE LNG PRODUCTION in POLAND. W.H. Isalski

Modelling of CO 2 capture using Aspen Plus for EDF power plant, Krakow, Poland

Advanced combustion systems for tube manufacturing

Refrigeration Kylteknik

2291-6A. Joint ICTP-IAEA Course on Science and Technology of Supercritical Water Cooled Reactors. 27 June - 1 July, 2011

Connecting Coolers to Superconducting Magnets with a Thermal-Siphon Cooling Loop

MOLECULAR GATE TECHNOLOGY FOR (SMALLER SCALE) LNG PRETREATMENT

Chemistry of Petrochemical Processes

Fundamental Investigation Of Whole-Life Power Plant Performance For Enhanced Geothermal Systems

Cooling Ceiling Panel

Transcription:

STORAGE OF HYDROGEN G. Spazzafumo Department of Industrial Engineering, University of Cassino, Italy Keywords: Hydrogen, Storage, Gaseous Hydrogen, Liquid Hydrogen, Hydrogen Slush, Cylinder, Cryogenic Tank, Glass Microsphere, Metal Hydride, Liquid Hydride, Cryogenic Adsorption, Carbon nanotube, Graphite Nanofiber Contents 1. Introduction 2. Gaseous hydrogen 3. Liquid hydrogen 4. Slush hydrogen 5. Metal hydrides 6. Liquid hydrides 7. Cryogenic adsorption 8. Carbon nanostructures 9. Hydrogen producing systems 10. Comparison among storage systems 11. Conclusion Glossary Bibliography Biographical Sketch Summary The hydrogen storage systems currently used or under development have been classified and described with reference to process, energy density, safety and other relevant aspects. Taking into account the hydrogen application fields as an energy carrier and taking into account various parameters which are significant in relation to the applications, a comparison among the hydrogen storage system has been carried out. Different systems are suitable for different applications and therefore could be developed as commercial systems in the future. The most recent storage system, graphite nanofibers, if experimental data will be confirmed, seems the unique one to be suitable for almost all applications. 1. Introduction Hydrogen is an energy carrier, probably the main energy carrier of the future. This is due to its extreme versatility: on the one hand it can replace the traditional fuels; on the other hand it is easily convertible into electrical energy through fuel cells with high efficiencies even by means of small plants. Until a few years ago the need for an energy carrier like this was due to the efforts towards substitution of the traditional sources of energy with renewable sources (including the nuclear fusion) and was therefore futuristic. Now the development of the fuel cells for transport and a greater awareness for pol-

lution problems in the urban environment are opening up these prospects for the use of hydrogen. An energy carrier must, of course, be stored and distributed. These operations may be performed according to different processes. Each method of storage has particular features which make it come out more or less suitable for the use in a particular field and which may influence the choice in an even more decisive way than cost considerations. For example, essentially important aspects are safety, especially when people are exposed to hazard, or the weight and the overall dimensions, for movable applications. The main forms of hydrogen storage may be pure or combined, and in the latter case there can be a chemical bond or a physical bond. Table 1 shows a classification of the different options. 2. Gaseous Hydrogen Pure Gaseous in cylinders Gaseous in microspheres Liquid Slush Combined Metal hydrides Liquid hydrides Carbon nanotubes Graphite nanofibers Table 1: Hydrogen storage systems Gaseous hydrogen is the simplest form for storage; in fact molecular hydrogen is found in nature in that state. Usually gaseous hydrogen is stored in pressurised tanks or in cylinders with pressure of 20 MPa. There are no particular problems with the gas compression and with the subsequent reuse. The idea of recovering part of the energy used for compression can be exploited with the steady applications, while it must be rejected for movable applications because the resulting advantage is insignificant, when compared with the problems arising out of the system dimensions and weight. The need to reach such high storage pressures is justified because of the gas low density. So it is necessary that the containers are made with materials and dimensions such as to endure higher pressures. The use of composite materials make possible to construct tanks that are light and suitable for very high pressures, up to 69 MPa. They can be used for hydrogen transport because they allow us to multiply the transport capacity by a factor 1.5 and this can balance the high cost of such tanks. However, as the pressure increases, there is an advantage as regards the dimensions which decrease progressively, while the weight decreases only at the beginning, when the variation is caused by the diameter decrease, and then it increases, when the thickness increase prevails (Figure 1). The energy density of a gaseous hydrogen tank is in the range 1-2.5 MJ/l when storage pressure varies between 20 and 70 MPa. Variation in weight also depends strongly on the materials: energy density is in the range 2-10 MJ/kg.

Figure 1: Energy density for gaseous hydrogen in cylinders By cooling compressed hydrogen down to 80 K, which is the liquefaction temperature for nitrogen, it is possible to increase density by a factor of 3.5. This requires materials suitable for very low temperature and also requires technological development especially for composite materials. If hydrogen is produced through electrolytic cells at a pressure of 4 MPa and a temperature of 300 K, to compress it until the pressure of 20 MPa is reached, with compression efficiency of 87%, it is necessary to supply a work of compression of 2.99 MJ/kg, that is to say equivalent to 2% of the higher heating value of hydrogen. Regarding safety the best solution would be to make a tank with spherical shape, because there would be a uniform stress on the material. On the other hand, this shape is undoubtedly the least convenient from the point of view of the disposition inside the system, especially if it is on board. Again regarding safety it is also important to consider the possible pressurisation due to the heating and therefore the arrangement of appropriate safety valves. Otherwise gaseous hydrogen could be stored in hollow glass microspheres, with a diameter between 50 and 100 nm. This method is based on the dependence of permeability to hydrogen on temperature of glass. Hydrogen is put inside the microspheres at a pressure of 40-50 MPa and at temperature of 300-400 C. After the cooling at ambient temperature, the hydrogen diffusion through the glass walls decreases so that, even if the outside pressure decreases to the atmospheric value, inside the microspheres hydrogen remains at high pressure. Actually, at ambient temperature as well there is a slight loss of hydrogen, because the glass at that temperature is not quite impermeable. The phenomenon, often unwelcome, is influenced by the characteristics of the glass used to make the microspheres, as well as, of course, by the pressure and temperature values. As the molar percentage of forms of glass without reticulated structure increases, the

losses of hydrogen decrease, but it has negative consequences when it is necessary to store hydrogen because of the decreased permeability. During the filling procedures, to limit the number of microspheres which break, to the minimum values, it is important to handle them taking particular precautions. First of all it is necessary to make the sample of microspheres reach the chosen filling temperature, usually 350 C. Then we must raise gradually the pressure, so that the difference of pressure between the microspheres outside and inside is not too high, usually 20 MPa. Amid one increase and the other of the pressure, it is necessary to allow hydrogen to pass through the glass walls, keeping a constant pressure for a sufficiently long time. Finally, the temperature will be brought back to the ambient value. To obtain the release of stored hydrogen, it is necessary to make the microspheres reach a temperature which must increase at the increasing of the desired outflow speed. The maximum working pressure of the microspheres is limited by the glass mechanical resistance and by the dimensions of the microspheres. The quantity of hydrogen which can permeate the microspheres walls, during the filling and emptying procedures, depends fundamentally upon the balls dimensions, upon the composition of the glass, upon the gradient of pressure through the walls and upon the temperature. In the ideal case of microspheres which are exactly alike and packed in a compact way, and ignoring the thickness as regards the diameter, the available volume is about the 74% of the total space, so hydrogen is stored with a size reduction of 150-200 times (storage energy density 1.8-2.4 MJ/l) From a gravimetric point of view the mass of hydrogen that can be stored per unit mass of the sample is about 0.04 (storage energy density 5.6 MJ/kg). If Q is the volume of hydrogen per unit mass of the sample passing through the microspheres walls with a difference of pressure ΔP, we have: S Q=k Δ P s where k is the permeability of the microspheres with respect to hydrogen, S is the medium surface per unit mass and s is the microspheres medium thickness. If T is the absolute temperature and M the molar percentage of forms of glass without reticulated structure, then we have: 3600+165 M - -19-4 3 T ( ) k = 2.21 10 3.4 + 8 10 M T e and if d is the medium diameter and ρ m is the microspheres density and ρ v is the glass density, then:

6 d ρ S= and m s= 1-3 1- d ρ m 2 ρv. 3. Liquid Hydrogen In order to increase further the volumetric density of hydrogen storage it is possible to resort to liquefaction and storage in cryogenic containers. Liquid hydrogen has, in fact, a density which is about 800 times higher than gaseous hydrogen at atmospheric pressure. The evaporation of liquid hydrogen inside the tanks must be very limited, because the consequent pressurisation of the tank could become a danger. So, the storage temperature being very low, it is necessary to have very efficient tank insulation systems, so as to reduce to minimum the heat exchange with the environment. The contemporary necessity for a good insulation and of a high resistance, also at low temperatures, implies a rise in weight and dimensions and therefore a great decrease in energy density: the weight of hydrogen stored is about 18% of the weight of the cryogenic tank and the volume available for hydrogen is 50% of the total of the tank. There are different cycles for the liquefaction of hydrogen. The simplest is the Linde cycle, which is schematically shown in Figure 2 with the schematic of a plant in Figure 3. This cycle is suitable for liquefaction of gases with an inversion temperature (that is the one at which the Joule-Thomson coefficient changes its sign) lower than the ambient temperature: for hydrogen this temperature is 205 K. Figure 2: Linde cycle for hydrogen liquefaction.

Figure 3: Schematic of a Linde plant for hydrogen liquefaction Liquefaction is done by a sequential process of compression-cooling-expansion. Because of the very low liquefaction temperature of hydrogen (20.3 K at atmospheric pressure) it is necessary to use a primary coolant, water, and a secondary coolant with an intermediate liquefaction temperature. This secondary coolant is liquefied with procedures similar to the ones for hydrogen and then it is made to evaporate again taking heat from hydrogen. Usually nitrogen is used as secondary coolant. Making an energy balance of the hydrogen circuit, we have: ( ) m h + m h = m h + m h + m h h 2 r 14 1 r 11 l 6 1 where m is the total mass flow (that is to the compressor) of hydrogen, m l is the mass flow of liquefied hydrogen, m is the mass flow of the coolant, h is the enthalpy. r If x=m r m and y=m l m we have y= ( h1 -h2 ) + x ( h11 -h14 ) ( h1 -h6 ) and the y = h -h h -h where h 3 and h are bound to the evaporation temperature of the secondary coolant. maximum fraction of liquefied hydrogen is ( ) ( ) max 8 3 8 6 8 The minimum specific energy necessary to isobarically liquefy hydrogen from ambient w = T s -s h - h = 14.2 MJ/kg temperature and atmospheric pressure is ( ) ( ) min 1 1 6 1 6 where T is the temperature and s the specific entropy and the subscripts 1 and 6 represent the initial and final conditions respectively. Actually the necessary energy is much higher and this is really one of the more critical aspects of the use on large scale of this form of storage. The Linde cycle underwent in time some changes directed to reduce the compression work and to increase the fraction of liquefied hydrogen. Besides, it is necessary to consider that the molecule of hydrogen can be found in two different energetic states, the ortho-hydrogen form and the para-hydrogen form, depending on the configuration of the nuclear spins. At ambient temperature hydrogen is formed by a mixture defined normal-

hydrogen and composed of 25% of para-hydrogen and 75% of ortho-hydrogen. For higher temperatures this mixture remains unchanged, while as the temperature decreases the para-hydrogen percentage increases progressively (Table 2). The energetic level of the para-hydrogen form is lower than that of the ortho-hydrogen form, so the isothermal change ortho-para-hydrogen is exothermic. Without the use of the catalysts this change occurs spontaneously and proceeds very slowly: then liquefaction of normal-hydrogen takes place and the subsequent heat transfer connected to the above mentioned change causes partial evaporation of hydrogen. To prevent this from occurring, it is necessary to use suitable catalysts and remove during the liquefaction phase the heat connected to the ortho-para-hydrogen change. K 10 20 30 40 50 100 200 300 500 % 100.00 99.82 97.02 88.73 77.05 38.62 25.97 25.07 25.00 J/mol 1417 1417 1417 1417 1416 1295 440 74 0 Table 2: Para-hydrogen contents and conversion heat at equilibrium conditions Figure 4: Schematic of a plant for hydrogen liquefaction

It follows that the plants used for the liquefaction of large quantities of hydrogen employ cycles much more complicated than the above seen one. One example is schematically shown in Figure 4. In such a plant the hydrogen to be liquefied is divided into a primary flow and a secondary flow. The first one is converted from ortho-hydrogen to para-hydrogen in two different stages: the former one (3-4) at the nitrogen boiling point and the second one (5-6) at the hydrogen boiling point. The secondary flow is partially liquefied (6-7) without conversion and is used in the second stage for primary flow conversion. The secondary flow is also used for energy recovery by means of two turbines which produce mechanical energy and at the same time provide for gas cooling. The work of liquefaction increases therefore at the increasing of the final concentration of para-hydrogen: to liquefy 1 kg of hydrogen with a final concentration of parahydrogen equal to 95% it is necessary to supply about 35 MJ, corresponding to about 2.5 times the minimum value and also to about the 25% of the energetic content of hydrogen itself. It is however necessary to notice that, because of the slowness of the spontaneous conversion process, the optimum concentration of para-hydrogen depends on the storage time. 4. Slush Hydrogen Slush hydrogen is a mixture composed of 50% of liquid hydrogen and 50% of solid hydrogen at the temperature and at the pressure of the triple point (table 3). The slush presents a fluid-dynamic behavior similar to the liquid hydrogen one and can be distributed exploiting the same infrastructures, besides it has a higher density (+15%) presenting less leaks for evaporation. liquid (NLC) liquid (TP) slush (TP) solid (TP) Temperature (K) 20.39 13.96 13.96 13.96 Pressure (kpa) 101.33 7.20 7.20 7.20 Density (kg/m 3 ) 70.77 77.22 81.50 86.50 Enthalpy (kj/kg) 271.0 212.8 183.7 154.6 Table 3: Hydrogen characteristics at Normal Liquefaction Condition and at Triple Point To produce the slush the method of scraping is used. This method utilizes liquid helium to produce solid hydrogen on the surface of a heat exchanger. The solid hydrogen produced is then scraped with a sort of auger and mixed with the residual liquid hydrogen, while gaseous helium is recycled to be liquefied again.

- - - TO ACCESS ALL THE 21 PAGES OF THIS CHAPTER, Visit: http://www.eolss.net/eolss-sampleallchapter.aspx Bibliography It is not useful to report a traditional bibliography because we can divide the hydrogen storage systems in two groups: established technologies, like gaseous hydrogen stored in cylinders or liquid hydrogen stored in cryogenic tanks: most of data are probably available in every energy related handbook; technologies under development, which require a further division: metal hydrides: it is possible to write some volumes only reporting the titles of articles written about them and most of the articles concern only one or few specific hydrides so that it is not useful at this concern; carbon structures: there are not many articles, but they are not easy to find; moreover, when this article will be published there will be other article more updated because this field is under continuos development. Interesting and updated articles can be found on: International Journal of Hydrogen Energy; Proceedings of World Hydrogen Energy Conference; Proceedings of HYPOTHESIS (Hydrogen Power Theoretical and Engineering Solutions International Symposium). A good suggestion is also to check on some relevant web sites. There are a lot of them more or less related to hydrogen; here are some of them: http://www.eren.doe.gov/hydrogen/ http://www.hyweb.de/gazette-e http://www.hfcletter.com/ http://www.ttcorp.com/nha/ http://www.rebresearch.com/h2links.html Biographical Sketch Giuseppe Spazzafumo was born in Carrara (Italy), 29 Dec. 1959. Graduated (Mechanical Engineering, thesis on plants of biophotoproduction of hydrogen) in 1984 from University of Pisa. PhD (Energetics, thesis on use of MCFC in powerplants) in 1988 at University of Pisa. He works at University of Cassino since 1989: Lecturer, 1989-1991, Assistant Professor, 1992-2002, and Associate Professor since 2002. Reviewer of International Journal of Hydrogen Energy since 1989. Coordinator of HYPOTHESIS (HYdrogen POwer THeoretical and Engineering Solution International Simposium) series. Secretary of Italian Hydrogen Association (Forum Italiano dell Idrogeno) since its creation. Member of IAHE Board of Directors since 1999.

Prof. Spazzafumo was engaged in researches in Hydrogen Energy Fields since 1984. He is author of more than 50 scientific papers related to hydrogen energy systems or fuel cell systems. His work concerns theoretical studies and simulations of hydrogen energy systems, from hydrogen production to hydrogen final use.