Modeling and optimization of multi-tubular metal hydride beds for efficient hydrogen storage Constantinos A. Krokos a,b , Dragan Nikolic a , Eustathios S. Kikkinides a, *, Michael C. Georgiadis c , Athanasios K. Stubos b a Department of Mechanical Engineering, University of Western Macedonia, Bakola & Salvera Str., 50100 Kozani, Greece b National Center of Scientific Research DEMOKRITOS, Institute of Nuclear Technology and Radiation Protection, 15310 Ag. Paraskevi Attikis, Athens, Greece c Department of Engineering Informatics and Telecommunications, University of Western Macedonia, Karamanli & Ligeris Str., 50100 Kozani, Greece article info Article history: Received 24 July 2009 Received in revised form 9 September 2009 Accepted 9 September 2009 Available online 2 October 2009 Keywords: Metal hydrides Hydrogen storage Multi-tubular reactors Numerical simulations abstract This work presents a novel systematic approach for the optimal design of a multi-tubular metal hydride tank, containing up to nine tubular metal hydride reactors, used for hydrogen storage. The tank is designed to store enough amount of hydrogen for 25 km range 1 , for a fuel cell vehicle. A detailed 3D Cartesian, mathematical model is developed and validated against a 2D cylindrical developed by Kikkinides et al. [1]. The objective is to find the optimal process design so as to increase the overall thermal efficiency, and thus minimize the storage time. Optimization results indicate that almost 90% improvement of the storage time can be achieved, over the case where the tank is not optimized and for a minimum storage capacity of 99% of the maximum value. ª 2009 Professor T. Nejat Veziroglu. Published by Elsevier Ltd. All rights reserved. 1. Introduction Hydrogen is shown to be the future fuel from the point of view of human fuel evolution. The fuel evolution experienced the history from coal through petroleum to natural gas following the direction of increasing the content of hydrogen, therefore, it must finally reach the destination of pure hydrogen [1–3]. Automotive industry is slowly moving towards this direction by developing fuel cell vehicles, which have the benefit of zero emissions. However, before introducing the fuel cell vehicle in the market, several issues including the hydrogen storage need to be addressed. Currently, high pressure tanks and liquid hydrogen tanks are used for road tests, but both technologies do not meet all the requirements of future fuel cell vehicles [3]. Storage of hydrogen in a pressurized cylinder is not likely to be applied in the future due to low density, high pressurization costs and safety considerations [2]. Liquid hydrogen could be applied if the unit cost becomes compa- rable with gasoline, yet inevitable boiling-off liquid might be in concern [2]. Metal hydride materials seem to exhibit good storage capabilities and constitute a promising direction towards the hydrogen economy [4–6]. This can be attributed to the reduced operating pressure compared to compressed gas technology, thus ensuring less weight and increased safety [1]. The current work is motivated by the desire to design a metal hydride tank, with small charge–discharge duration germane * Corresponding author. E-mail address: [email protected](E.S. Kikkinides). 1 Based on data from Honda FCX Clarity. Available at www.sciencedirect.com journal homepage: www.elsevier.com/locate/he 0360-3199/$ – see front matter ª 2009 Professor T. Nejat Veziroglu. Published by Elsevier Ltd. All rights reserved. doi:10.1016/j.ijhydene.2009.09.021 international journal of hydrogen energy 34 (2009) 9128–9140
13
Embed
Modeling and optimization of multi-tubular metal hydride beds for efficient hydrogen storage
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
i n t e r n a t i o n a l j o u r n a l o f h y d r o g e n e n e r g y 3 4 ( 2 0 0 9 ) 9 1 2 8 – 9 1 4 0
Avai lab le at www.sc iencedi rect .com
journa l homepage : www.e lsev ie r . com/ loca te /he
Modeling and optimization of multi-tubular metal hydridebeds for efficient hydrogen storage
Constantinos A. Krokos a,b, Dragan Nikolic a, Eustathios S. Kikkinides a,*,Michael C. Georgiadis c, Athanasios K. Stubos b
a Department of Mechanical Engineering, University of Western Macedonia, Bakola & Salvera Str., 50100 Kozani, Greeceb National Center of Scientific Research DEMOKRITOS, Institute of Nuclear Technology and Radiation Protection,
15310 Ag. Paraskevi Attikis, Athens, Greecec Department of Engineering Informatics and Telecommunications, University of Western Macedonia, Karamanli & Ligeris Str.,
Fig. 8 – Time-space evolution of temperature in the reactor
in the reactor. (a) Reactor with one metal hydride bed, (b)
reactor with five metal hydride beds. Axial z [ 0.5, t [ 50 s.
a
b
Fig. 9 – Number of metal hydride tubes effect on (a) storage
time and (b) cooling time, T0 [ 290 K.
i n t e r n a t i o n a l j o u r n a l o f h y d r o g e n e n e r g y 3 4 ( 2 0 0 9 ) 9 1 2 8 – 9 1 4 0 9135
reason it can be used in more complex geometries. A second
step was to study the temporal evolution of the average
temperature in the metal hydride bed, for two cooling fluid
velocities, 0 m/s and 1 m/s. The simulation results produced by
the two models, as shown in Fig. 6, are almost identical. Indeed
the error of the results produced by the 3D Cartesian model
against the results produced by the 2D cylindrical is less than 1%
thus indicating the accuracy of the 3D Cartesian model.
From the above results it is evident that the developed
Cartesian model can represent accurately the cylindrical
geometry of the MH storage tanks provided that we use a fine
enough resolution to describe the reactor’s cross-sectional area.
If the mesh resolution is not fine enough we may expect some
error in the treatment of the boundary layer at the powder bed/
wall interface. Nevertheless, this error is much smaller when
comparing volume averaged quantities including the average
amount of hydrogen stored in the MHT as a function of time. For
this reason extensive simulations have been performed using
a single MHT reactor to study the effect of resolution on the
average amount of hydrogen stored in the MHT. Our studies
have shown that for 15 cells/diameter, the relative error is
approximately 4.6% at small to medium storage times (70% of
maximumcapacity) dropping to 1.7% at large storage times (99%
of maximum storage). The error is still tolerable even for
12/cells/diameter becoming 7% and 2.9%, at small and large
times, respectively. The above trend remains when using more
than 1 MHT, so the impact on the reduction of storage time due
to thepresence ofmulti-tubularconfigurations isnot affectedby
errors in resolution. Obviously improvement of the resolution
can be easily achieved provided that the necessary resources of
computer memory are available.
3.2. Optimization of the thermal efficiency of the metalhydride reactor
There are two main factors affecting the storage efficiency of
a metal hydride reactor. The first factor is the type of the metal
hydride alloy used in the reactor, which defines the maximum
capacity of the tank. The second is the thermal efficiency,
which defines the time required for complete hydrogen
storage. In this case the metal hydride used in the simulations
is LaNi5 with a maximum theoretical capacity of 0.0138 kg
hydrogen/kg alloy. So, this work focused on the optimal heat
management of the reactor in order to minimize the storage
time. More specifically, first the cooling fluid velocity is
Fig. 10 – Studied configurations and the temperature profiles inside the reactor for each configuration. Axial z [ 0.5, t [ 50 s.
i n t e r n a t i o n a l j o u r n a l o f h y d r o g e n e n e r g y 3 4 ( 2 0 0 9 ) 9 1 2 8 – 9 1 4 09136
a
b
Fig. 11 – Configuration effect on (a) storage time and (b)
cooling time, T0 [ 290 K.
Fig. 12 – Metal hydride tanks. (a) The base-case tank, (b) the
optimized tank.
i n t e r n a t i o n a l j o u r n a l o f h y d r o g e n e n e r g y 3 4 ( 2 0 0 9 ) 9 1 2 8 – 9 1 4 0 9137
optimized followed by the number and positions of the metal
hydride tubes inside the reactor.
3.2.1. Effect of cooling fluid velocitySeveral cooling fluid velocities were studied for the base case
of the modular reactor shown in Fig. 4. The effect of cooling
fluid velocity on the storage time is shown in Fig. 7(a) while
Fig. 7(b) shows the effect of cooling fluid velocity on cooling
time, which is the time required for the temperature to drop to
a certain value. It is clear that the storage time depends
strongly on the amount of hydrogen adsorbed which in this
case is measured as the percentage of the maximum theo-
retical amount. It is clear that a limit cooling fluid velocity
exists, in this case 1 m/s, beyond which neither the storage
time nor the cooling time are affected.
The exothermic absorption process increases the temper-
ature in the metal hydride bed analogous to the formation
enthalpy of the hydride (DH¼�29879 J mol�1 for LaNi5 case).
The heat generated is then transported to the bed’s boundary
by conduction. The heat transfer rate from the metal hydride
alloy to the boundary is defined by the thermal conductivity of
the alloy (lMH¼ 0.66 W m�1 K�1 for LaNi5 case). The cooling
fluid rejects the available amount of heat from the boundary
to the environment. Increasing the cooling fluid velocity
increases the heat transfer rate from the metal hydride
boundary to the environment. The problem is that if the
cooling fluid velocity is increased above 1 m/s there is no
improvement on heat transfer rate because there is no avail-
able heat to be transferred. Conduction thus becomes the
limiting factor.
Despite the existence of a limiting cooling fluid velocity the
time required for 99% hydrogen storage changes from 6910 s
when the cooling fluid velocity is 0 m/s to 4780 s when the
cooling fluid velocity is 1 m/s, indicating an improvement of
30%. Similarly the time required for the temperature inside the
metal hydride tube, to drop to T/T0¼ 1.01 changes from 8200 s
to 5590 s indicating an improvement of 31%. Note that the tank
was designed to store enough amount of hydrogen for 25 km
using the 99% of its theoretical maximum storage capacity.
It should be noted that the assumption of constant cooling
velocity is only justified at high cooling fluid velocities. For
low cooling fluid velocities on the other hand it is better to
include an appropriate heat transfer coefficient to account for
the convective boundary layer effect. This is easy to imple-
ment but for the purposes of the present work, which is to
provide the methodology and tools to perform simulations
and optimizations of hydrogen storage using multi-tube
metal hydride reactor configurations, it will not change any
of the conclusions in a substantial manner. The developed
in-house software is relatively straightforward to build, easy
to implement and maintain and gives accurate results
compared to cylindrical geometries provided that the mesh is
fine enough. Furthermore an even more accurate treatment
would be to employ a commercial software that allows for
i n t e r n a t i o n a l j o u r n a l o f h y d r o g e n e n e r g y 3 4 ( 2 0 0 9 ) 9 1 2 8 – 9 1 4 09138
more suitable meshes and modeling of more complex flows
for the coolant, however, there are currently some compu-
tational difficulties and limitations in interfacing gPROMS
with such a software. Nevertheless this will be systematically
considered in a future work.
3.2.2. Effect of the number of metal hydride tubesIncreasing the cooling fluid velocity is not the panacea for the
reduction of storage time. More complicated, non-concentric
metal hydride reactor geometries should be studied to achieve
a further decrease of the storage time. The same amount of
metal hydride as the concentric reactor shown in Fig. 4, was
divided to 2, 3, 4, 5 and 9 tubes to test multi-tubular metal
hydride reactor geometries. It is essential to keep the amount
of metal hydride the same in order to ensure a constant
overall storage capacity of the tanks under consideration. The
cooling fluid velocity was also kept constant and equal to 1 m/
s for all the cases studied.
Fig. 8 shows that the division of the metal hydride alloy in
several tubes reduces the zones of the temperature gradients
that occur inside the reactor and thus improve the perfor-
mance of the storage time by reducing the time needed for
complete storage and release of hydrogen. Fig. 9 shows the
effect of the number of the metal hydride tubes on the
storage and cooling time. As the number of metal hydride
tubes increases, the time for complete storage decays expo-
nentially. More specifically the time required for 99%
hydrogen storage reduces from 4789 s, when one metal
hydride tube was used, to 2010 s when metal hydride was
divided to three metal hydride tubes, indicating an
improvement of 58%. If the number of tubes is increased to
nine then the time required for 99% storage becomes 1030 s
which represents a total improvement of 78%. Analogously
the time required for the temperature, inside the metal
hydride tubes, to drop to T/T0¼ 1.01 decreased from 5590 s,
when a single metal hydride tube was used, to 2350 s when
three metal hydride tubes were used, indicating an
improvement of 58%. When the number of tubes is increased
to nine then a further improvement of 48.9% on the cooling
time was achieved.
The improvement of cooling and storage time, when the
metal hydride is divided in several tubes is attributed to the
increment of the Surface/Volume ratio as the number of metal
hydride tubes increases. More specifically:
AV¼ 2$p$r$L
p$r2$L¼ 2
r; (23)
Note that the Surface/Volume ratio increases as the metal
hydride tube radius decreases.
3.2.3. Effect of metal hydride tubes position in the reactorIt is important to investigate if there is an ‘‘optimal’’ config-
uration (metal hydride tubes position in the reactor). Hence
we studied five different configurations for the case of the nine
metal hydride tubes reactor, to explore the effect of the posi-
tion of the tubes on the storage and cooling time. The five
studied configurations are shown in Fig. 10.
The effect of process configuration on cooling and
storage time is shown in Fig. 11. It is evident that the posi-
tion of the metal hydride tubes plays an important role
regarding the temperature profile inside the reactor.
Therefore the configuration with the best storage time is
configuration E, where the metal hydride tubes are placed in
such a manner so as the available space of the reactor
(where cooling fluid resides) is more efficiently utilized. So,
if they the tubes are not evenly distributed with respect to
each other and the storage container then the available
cooling fluid is not effectively utilized. For example in
configuration C the right part of the cross-sectional area is
not-well utilized as the cooling fluid in that region flows
without removing enough heat from the reactors. On the
other hand, the upper and lower regions are over-crowded
with MH reactors and as result there is not enough cooling
fluid per unit time to remove the produced heat due to
absorption of hydrogen in the neighboring MH reactors. If
the reactors touch each other we loose interfacial area of
reactor to fluid. More specifically the time required for 99%
hydrogen storage changes from 1030 s, for the configuration
A to 780 s when the configuration E was used, indicating an
improvement of 25%. Similarly the time required for the
temperature to drop to T/T0¼ 1.01 changed from 1200 s to
900 s implying an improvement of 25%.
In summary, it should be emphasized that the aim of the
present work is to provide a systematic methodology and
modeling tools to perform simulations and optimizations of
hydrogen storage using multi-tube metal hydride reactor
configurations. This has been illustrated by several cases. The
main modeling assumptions of this study relate to the use of
a constant cooling fluid velocity, the rather simple represen-
tation of the wall structure and the negligence of the radial
flow within the hydride bed. The latter has indeed been found
negligible in all cases studied. As for the remaining assump-
tions, they do not affect in any significant manner the
conclusions and principal goals of the work. While the quali-
tative tendencies found and the proposed methodology, do
not bear any dependence on these assumptions, it is true that
the actual calculated tank filling times carry an estimated
uncertainty of about 10–15% as a result of the simplifying
assumptions. On the other hand, a worth mentioning limita-
tion of the computational approach has to do with the mesh
resolution employed. As shown by extensive simulations that
were performed specifically to study the effect of resolution,
the errors induced are less than 10% in all cases, i.e. smaller
than the uncertainty imposed by the assumptions.
4. Conclusions
A 3D Cartesian mathematical model has been developed to
study complex metal hydride reactor geometries. The
Cartesian model has been validated and then used to study
the effect of three basic parameters on the storage and
cooling time of the tank. The first parameter is the cooling
fluid velocity whose optimal value found to 1 m/s. A further
step along this study was to explore the effect of the number
of the metal hydride tubes (for the same amount of metal
hydride alloy) on the storage and cooling time. Then the
effect of the number of metal hydride tubes and their posi-
tion was systematically studied. It was shown that the
configuration with nine tubes resulted in the highest
i n t e r n a t i o n a l j o u r n a l o f h y d r o g e n e n e r g y 3 4 ( 2 0 0 9 ) 9 1 2 8 – 9 1 4 0 9139
improvement on the storage and cooling time. The justifica-
tion of this trend is due to the fact that as the number of
metal hydride tubes increases the radius of the tubes
decreases and thus the Surface/Volume ratio increases.
Finally, the optimal position of the metal hydride tubes was
also studied. It was illustrated that the tubes should be
placed in a uniform manner to utilize more effectively the
available space in which the cooling medium resides that
leading to a more effective heat management.
The base-case tank and the optimized tank, both capable
for storing enough amount of hydrogen for 25 km range, are
shown in Fig. 12. The time required for 99% hydrogen storage
changed from 6910 s, for the base-case tank with cooling fluid
velocity 0 m/s, to 780 s, when the cooling fluid velocity, the
number of metal hydride tubes and the configuration of the
tubes in the reactor were optimized, indicating a total
improvement of 88.7%. Evidently the above figures could be
further improved if more efficient metal hydrides are applied
and this is the subject for a future work.
Acknowledgements
Financial support from the European Commission under the
DIAMANTE Marie Curie ToK project of FP6 (contract No: MTKI-
CT-2005-029544) is gratefully acknowledged. A. K. Stubos and
C.A. Krokos wish to acknowledge the EU/NESSHY Integrated
Project (SES6-518271) for partial financial support.
r e f e r e n c e s
[1] Kikkinides E, Georiadis M, Stubos A. On the optimization ofhydrogen storage in metal hydride beds. Int J HydrogenEnergy 2006;31(6):737–51.
[2] Zhou Li. Progress and problems in hydrogen storage methods.J Renewable Sustainable Energy Rev. 2005;9(4):395–408.
[3] Kikkinides ES. Hydrogen-based energy systems: the storagechallenge. In: Georgiadis MC, Kikkinides ES,Pistikopoulos EN, editors. Chapter 3. Weinheim: Wiley-VCHVerlag, ISBN 978-3-527-31694-6; 2008. p. 85–123.
[4] Mohanty KK. The near-term energy challenge. AIChE J 2003;49(10):2454–60.
[5] Kaplan Y, Veziroglou TN. Mathematical modelling ofhydrogen storage in LaNi5 hydride bed. Int J Energy Res 2003;27(11):1027–38.
[6] Sakintuna B, Lamari-Darkrim F, Hirscher M. Metal hydridematerials for solid hydrogen storage: a review. Int J HydrogenEnergy 2007;32(9):1121–40.
[7] Mellouli S, Askri F, Dhaou H, Jemni A, Nasrallah SB. Numericalstudy of heat exchanger effects on charge/discharge times ofmetal-hydrogen storage vessel. Int J Hydrogen Energy 2009;34(7):3005–17.
[8] Muthukumar P, Prakash Maiya M, Srinivasa Murthy S.Experiment on a metal hydride-based hydrogen storagedevice. Int J Hydrogen Energy 2005;30(15):1569–81.
[9] MacDonald B, Rowe A. Experimental and numerical analysisof dynamic metal hydride hydrogen storage systems. J PowerSources 2007;174(1):282–93.
[10] Jemni A, Nasrallah SB. Study of two-dimensional heat andmass transfer during absorption in a metal-hydrogenreactor. Int J Hydrogen Energy 1995;20(1):43–52.
[11] Nasrallah SB, Jemni A. Study of two-dimensional heat andmass transfer during desorption in a metal-hydrogenreactor. Int J Hydrogen Energy 1997;22(1):67–76.
[12] Jemni A, Nasrallah B, Lamloumi J. Experimental andtheoretical study of a metal-hydrogen reactor. Int J HydrogenEnergy 1999;24(7):631–44.
[13] Nakagawa T, Inomata A, Aoki H, Miura T. Numerical analysisof heat and mass transfer characteristics in the metalhydride bed. Int J Hydrogen Energy 2000;25(4):339–50.
[14] Mat MD, Kaplan Y. Numerical study of hydrogen absorptionin an LM–Ni5 hydride reactor. Int J Hydrogen Energy 2001;26(9):957–63.
[15] Demircan A, Demiralp M, Kaplan Y, Mat MD, Veziroglu TN.Experimental and theoretical analysis of hydrogenabsorption in LaNi5H2 reactors. Int J Hydrogen Energy 2005;30(13):1437–46.
[16] Kaplan Y. Effect of design parameters on enhancement ofhydrogen charging in metal hydride reactors. Int J HydrogenEnergy 2009;34(5):2288–94.
[17] Aldas K, Mat MD, Kaplan Y. A three-dimensional model forabsorption in a metal hydride bed. Int J Hydrogen Energy2002;27(10):1049–56.
[18] Mat MD, Kaplan Y, Aldas K. Investigation of three-dimensional heat and mass transfer in a metal hydridereactor. Int J Energy Research 2002;26(11):973–86.
[19] Askri F, Jemni A, Nasrallah SB. Study of two-dimensionaland dynamic heat and mass transfer in a metal-hydrogenreactor. Int J Hydrogen Energy 2003;28(5):537–57.
[20] Lozano AG, Eigen N, Keller C, Dornheim M, Bormann R.Effects of heat transfer on the sorption kinetics of complexhydride reacting systems. Int J Hydrogen Energy 2009;34(4):1896–903.
[21] Chung CA, Ci-Jyun Ho. Thermal-fluid behavior of the hydridingand dehydriding processes ina metal hydride hydrogenstoragecanister. Int J Hydrogen Energy 2009;34(10):4351–64.
[22] Askri F, Salah BM, Jemni A, Nasrallah BS. Optimization ofhydrogen storage in metal-hydride tanks. Int J HydrogenEnergy 2009;34(2):897–905.
[23] Gambini M, Manno M, Vellini M. Numerical analysis andperformance assessment of metal hydride-basedhydrogen storage systems. Int J Hydrogen Energy 2008;33(21):6178–87.
[24] Brown TM, Brouwer J, Samuelsen GS, Holcomb FH, King J.Accurate simplified dynamic model a metal hydride tank.Int J Hydrogen Energy 2008;33(20):5596–605.
[25] Mayer U, Groll M, Super W. Heat and mass transfer in metalhydride reaction beds: experimental and theoretical results.J Less Common Metals 1987;131(1–2):235–44.
[26] Chaise A, Marty P, De Rango P, Fruchart D. A simple criterionfor estimating the effect of pressure gradients duringhydrogen absorption in a hydride reactor. Int J Heat MassTransfer 2009;52(19–20):4564–72.
[27] de Rango P, Chaise A, Charbonnier J, Fruchart D, Jehan M,Marty P, et al. Nanostructured magnesium hydride for pilottank development. J Alloys Compounds 2007;446–447:52–7.
[28] Huot J, Liang G, Boily S, Van Neste A, Schultz R. Structural studyand hydrogen sorption kinetics of ball-milled magnesiumhydride. J Alloys Compounds 1999;293–295:495–500.
[29] Levesque S, Ciureanu M, Roberge R, Motyka T. Hydrogenstorage for fuel cell systems with stationary applications I.Transient measurement technique for packed bedevaluation. Int J Hydrogen Energy 2000;25(11):1095–105.
[30] Nagel M, Komazaki Y, Suda S. Effective thermal conductivityof a metal hydride bed augmented with a copper wire-netmatrix. J Less Common Metals 1987;131(1–2):426.
[31] Kim KJ, Lloyd G, Razani Jr A, Feldman KT. Thermal analysisof the Ca0.4Mm0.6Ni5 metal-hydride reactor. Applied ThermalEngineering 1998;18(12):1325–36.
i n t e r n a t i o n a l j o u r n a l o f h y d r o g e n e n e r g y 3 4 ( 2 0 0 9 ) 9 1 2 8 – 9 1 4 09140
[32] Kim KJ, Feldman Jr KT, Lloyd G, Razani A, Shanahan KL.Performance of high power metal hydride reactors.Int J Hydrogen Energy 1998;23(5):355–62.
[33] Mohan G, Maiya PM, Murthy SS. Performance simulation ofmetal hydride hydrogen storage device with embeddedfilters and heat exchanger tubes. Int J Hydrogen Energy 2007;32(18):4978–87.
[34] Mellouli S, Askri F, Dhaou H, Jemni A, Ben Nasrallah S.A novel design of a heat exchanger for a metal-hydrogenreactor. Int J Hydrogen Energy 2007;32(15):3501–7.