SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
This document summarizes peer review comments and scores for the fiscal year 2023 U.S. Department of Energy (DOE) Hydrogen Program Annual Merit Review and Peer Evaluation Meeting (AMR), held June 6-9, 2023. This review process provides project- and program-level evaluations of DOE-funded research, development, demonstration, and analysis of hydrogen and fuel cell technologies.
This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.
Hydrogen has been promoted as a key component of global decarbonization efforts, with various past studies estimating carbon emissions for several production pathways, but little past work has considered its water resource needs. This life cycle analysis considers hydrogen production on a per-kilogram basis for 11 pathways, fossil and non-fossil. It also includes impacts of treating water to the required quality for hydrogen production. Greenhouse gas emissions results were in a range of –15 to +31 kg CO 2 e/kg H 2 produced. Water consumption varied more widely, from about 7 to 55 kg water/kg H 2 for fossil-based pathways and 530 to 3400 kg water/kg H 2 for biomass-based pathways. Electrolysis with various renewable electricity scenarios were also modeled. Altogether, there are challenging tradeoffs to be navigated to achieve a low carbon and water footprint in hydrogen economy.
In this report a series of iron (Fe) modified CuO-ZnO-ZrO 2 -Al 2 O 3 (CZZA) catalysts, with various Fe loadings, were prepared using a co-precipitation method. A bifunctional catalyst, consisting of Fe-modified CZZA and HZSM-5, was studied for dimethyl ether (DME) synthesis via CO 2 hydrogenation. The effects of Fe loading, reaction temperature, reaction pressure, space velocity, and concentrations of precursor for the synthesis of the Fe-modified CZZA catalyst on the catalytic activity of DME synthesis were investigated. Long-term stability tests showed that Fe modification of the CZZA catalyst improved the catalyst stability for DME synthesis via CO 2 hydrogenation. The activity loss, in terms of DME yield, was significantly reduced from 4.2% to 1.4% in a 100 h run of reaction, when the Fe loading amount was 0.5 (molar ratio of Fe to Cu). An analysis of hydrogen temperature programmed reduction revealed that the introduction of Fe improved the reducibility of the catalysts, due to assisted adsorption of H 2 on iron oxide. The good stability of Fe-modified CZZA catalysts in the DME formation was most likely attributed to oxygen spillover that was introduced by the addition of iron oxide. This could have inhibited the oxidation of the Cu surface and enhanced the thermal stability of copper during long-term reactions.
Ignition delay times from undiluted mixtures of natural gas (NG)/H2/Air and NG/NH3/Air were measured using a high-pressure shock tube at the University of Central Florida. The combustion temperatures were experimentally tested between 1000-1500 K near a constant pressure of 25 bar. Mixtures were kept undiluted to replicate the same chemistry pathways seen in gas turbine combustion chambers. Recorded combustion pressures exceeded 200 bar due to the large energy release, hence why these were performed at the high-pressure shock tube facility. The data is compared to the predictions of the NUIGMech 1.1 mechanism for chemical kinetic model validation and refinement. An exceptional agreement was shown for stoichiometric conditions in all cases but strayed at lean and rich equivalence ratios, especially in the lower temperature regime of H2 addition and all temperature ranges of the baseline NG mixture. Hydrogen addition also decreased ignition delay times by nearly 90%, while NH3 fuel addition made no noticeable difference in ignition delay time. NG/NH3 exhibited similar chemistry to pure NG under the same conditions, which is shown in a sensitivity analysis, demonstrating hydrogen chemistry to be dominant in NG/H2 mixtures and hydrocarbon chemistry to be dominant in NG/NH3 mixtures. The reaction CH3 + O2 = CH3O + O is identified and suggested as a possible modification target to improve model performance. Increasing the robustness of chemical kinetic models via experimental validation will directly aid in designing next-generation combustion chambers for use in gas turbines, which in turn will greatly lower global emissions and reduce greenhouse effects.
The overarching goal of the reported project was to develop reversible metal hydrides that can be dehydrogenated and recharged at relatively low temperatures and pressures. The project complemented the research on reversible metal hydrides for hydrogen storage conducted by the Hydrogen Materials – Advanced Research Consortium (HyMARC). The students involved were trained at both the University of Texas at El Paso (UTEP) and Sandia National Laboratories (SNL). The research objectives were to fabricate nanoscale magnesium boride powders with different B/Mg atomic ratios and to investigate the effects of B/Mg ratio and particle size on the hydrogenation properties of the obtained materials. Magnesium diboride (MgB 2 ) and tetraboride (MgB 4 ) were synthesized by using combustion synthesis and by heating in a tube furnace. The latter method resulted in materials with significantly reduced impurity levels. An effective method for the removal of magnesium oxide impurities was identified. High-energy ball milling produced submicron MgB 2 and MgB 4 powders. However, hydrogenation experiments at pressures up to 700 bar have not shown any uptake of hydrogen. Thermogravimetric analysis of the decomposition of the obtained MgB 2 and MgB 4 demonstrated higher thermodynamic stability of the latter. Thus, the boron-rich nature and phase stability of magnesium borides with high B/Mg ratios make them poor candidates for direct hydrogenation. However, the synthesized MgB 4 is a promising energetic additive to solid fuels for high-speed air-breathing propulsion as it offers a lower oxidation onset temperature and greater energy density than MgB 2 .
In this work CuFeMnO 4 (copper iron manganese oxide) oxygen carrier was characterized using DSC/TGA, in-situ XRD, and SEM-EDS to gain a clear elucidation of the chemical looping combustion (CLC) reactions with H2 which is a component of synthesis gas. A reaction model which best described the experimental TGA reaction data was selected from the analysis. Intrinsic rate of reaction parameters obtained from the study could be used for designing a reactor for scale-up. A model with combined contributions from nucleation and growth model and dimensional phase boundary model was found to exhibit the best correlation with the experimental data. The analysis provided validation of intrinsic reaction rates and other rate parameters.
This document summarizes peer review comments and scores for the fiscal year 2021 U.S. Department of Energy (DOE) Hydrogen Program Annual Merit Review and Peer Evaluation Meeting (AMR), held virtually June 7-11, 2021. In response to direction from various stakeholders, including the National Academies, this review process provides project- and program-level evaluations of DOE-funded research, development, demonstration, and analysis of hydrogen and fuel cell technologies.
Compound-specific stable isotope analysis of hydrogen (δD) and carbon (δ 13 C) in organic compounds is a valuable tool in biogeochemical research. A key limitation of this method is the relatively large amount of sample required to achieve desirable precision. We developed a large-volume (20 μL) injection method that allows for high throughput analysis of less concentrated samples and tested it for δ 13 C and δD measurements of n-alkanes. We also conducted a comparison of reference standards and assessed several methods to normalize and correct n-alkane δD and δ13C measurements. The mean precision of the δD method based on 233 environmental n-alkane samples (two to three replications per sample) is 4.0‰ (1σ, estimated from the weighted mean of the pooled unbiased standard deviations) and 0.46‰ (1σ) for δ 13 C from 37 environmental samples (two to three replications per sample). The evaluation of reference standards shows that the use of n-alkane standards with large offsets in δD values in adjacent n-alkane chains can lead to biases in measurement correction. The large-volume injection method shows good reproducibility of δ 13 C and δD measurements of n-alkanes and reduces the required sample concentration by about 80%. We propose that for δD measurements, a reference standard set should be used in which each reference standard has a limited range of δD values and no adjacent n-alkane chains, to minimize memory effects.
This gap analysis study included all parts of the hydrogen value chain, from hydrogen generation, to storage, and end use, with the exception of transportation applications. While these research projects are not always commercially viable, they are important in advancing technology, and especially in reaching goals such as improved efficiency and reducing emissions. Without emissions targets enforced by the government, or incentives for industry for emissions reduction, investing in advanced technology and manufacturing system upgrades can be exceedingly difficult to justify. This analysis identifies and prioritizes sensor R&D development areas that are critical to enable the hydrogen economy.
Alkali metals have been recognized as effective promoters in heterogeneous catalysis, capable of enhancing catalytic activity and tuning product distributions. Over the past few decades, significant efforts have been made aiming to reveal the mechanisms underlying the promoting effect of alkalis. However, the roles that alkali metals play in the catalytic process remain elusive due to challenges in capturing their catalytic behaviours upon exposure to reactive environments. This review summarizes recent surface science and theoretical studies of alkali (potassium, cesium)-decorated metal and metal oxide model catalysts, revealing the crucial tuning by alkalis of activity and selectivity for CO 2 hydrogenation. The analysis of electronic structures identifies the selective binding mechanism of the positively charged alkali ions on the surface, being able to reduce the surface work function and lead to strong electron polarization on the surfaces. Depending on the alkali–support interaction, the deposition of alkalis can selectively modify the bindings of reaction intermediates involved in CO 2 hydrogenation via the interplay among the ionic, covalent and electrostatic tunings. As a result, CO 2 can be effectively activated and converted into diverse products at the alkali–support interface, ranging from formic acid to methanol and ethanol. The identified selective bond-tuning advances the application of alkalis in promoting catalytic activity and controlling catalytic selectivity at alkali–support interfaces.
Polymer materials are used widely as sealing materials and liners in hydrogen storage and transportation infrastructure. During operating conditions, polymers are exposed to high pressure hydrogen gas periodically. When in contact, the hydrogen gas diffuses through the polymer material with ease due to its smaller sized molecules compared to the polymer molecules. In the event of rapid decompression of outside hydrogen gas, the hydrogen molecules within the polymer diffuse out. This process is known as the desorption of polymer and can cause blisters or cracking causing permanent damage to the polymer. Here, we used a combined experimental and modeling approach to investigate the desorption of the hydrogen gas in the polymer. Thermal desorption analysis (TDA) experiments were performed on the polymer sample exposed to the high-pressure hydrogen gas to measure the diffusion coefficient and equilibrium hydrogen concentration. Experiments also provided the quantity of hydrogen within the polymer during the desorption process which followed Fick’s law. However, for higher pressurization values, some anomaly was observed in the quantity of hydrogen remaining in the polymer sample. A continuum mechanics based coupled diffusion–deformation-damage model was deployed to simulate the gas desorption process in the material using the finite element method to understand its effect on the material. It was found that the anomaly in TDA results were due to the blister formation on the polymer surface which resulted in longer hydrogen retention. Simulations also showed that blister formation can be affected by multiple parameters such as the addition of filler particles, and pre-existing cavity size.
This paper develops a new approach for analysis of subgrid closures for turbulent combustion as modeled using direct quadrature (finite-rate chemistry) techniques. The approach, termed multiresolution analysis through mesh-sequenced realizations (MRA-MSR), conducts simultaneous, constrained large-eddy simulations on a set of hierarchically coarsened meshes. Furthermore, the availability of underlying fine-mesh (subgrid) data corresponding to coarse-mesh locations allows a clearer assessment of the effects of unresolved fluctuations on apparent reactivity. A key to MRA-MSR is the correlation of eddy structures at coarser mesh levels, which is facilitated by the transfer of filtered fine-mesh velocity information. A seven-mesh MRA-MSR hierarchy using three resolution levels is applied to one of the Sydney bluff-body stabilized methane–hydrogen flames. Analysis of the simultaneously evolved data at different resolution levels reveals several interesting trends. First, at high Damköhler numbers, there is clear evidence of attenuation of apparent reactivity due to the effects of unresolved fluctuations. Secondly, single-point, single-time filtered density functions of a normalized subgrid Damköhler number show a characteristic beta probability density function (PDF) form and display evidence of scale similarity. Interrogation of the MRA-MSR database also shows that the recently-proposed least-squares minimization (LSM) turbulence–chemistry interaction model can account for the observed diminishment in reactivity at high Damköhler numbers but cannot reduce scatter significantly. A new form of the LSM model, which makes use of the normalized subgrid Damköhler number beta PDF distribution, performs slightly better than the original model, illustrating the potential of MRA-MSR both in assessing existing closure concepts and in developing new ones.
This study presents the techno-economic analysis (TEA) results of integrating electrolysis hydrogen (H 2 ) production, compressed air energy storage (CAES), and H 2 -fired combustion turbines in a high variable renewable energy (VRE) market environment. The inconsistent nature of VRE creates challenges for power producers in maintaining a stable electrical grid as it is increasingly utilized. The mission of the National Energy Technology Laboratory (NETL) is driving innovation and delivering energy solutions, and a multiangle approach involving H 2 production, energy storage, and next-level H 2 -fired combustion turbine generator (CTG) technologies could provide one solution for the nation’s growing electrical grid issues with increasing VRE sources. The H 2 production and CAES concept were investigated due to its ability to provide utility-scale H 2 -fueled power generation with large-scale energy storage capabilities. Two facilities with CAES and natural gas-fired CTGs (in McIntosh, Alabama, and Huntorf, Germany) have been operating for decades. In contrast, this study investigates the potential to replace the natural gas fuel with H 2 fuel. This concept has been publicly presented by both Siemens Energy and Bechtel Global Engineering, Construction & Project Management (Bechtel) at two different power generating levels. The CAES and air expansion/combustion turbine power generation in this study are primarily based on the Siemens Energy system (Bailie, Aug. 10-11, 2021) (Scheller, Feb. 21, 2023). The inclusion of a H 2 turboexpander generator is from Bechtel (Gülen, Sep. 6, 2022). A block flow diagram of the proposed process is illustrated in Exhibit ES-1.
In this work, neutron computed tomography (CT) is employed to investigate the dissolution of porous aluminum electrodes during electrocoagulation (EC). Porous electrodes were chosen in efforts to reduce electric power requirements by using larger surface-area electrodes, having both inner and outer surface, for the EC process. Neutron CT allowed 3D reconstruction of the porous electrodes, and image analysis provided the volume of each electrode vs. thickness, which can indicate whether the inner surface is effectively involved in EC reactions. For the anode, the volume decreased uniformly throughout the thickness of the electrode, indicating that both the outer and inner surface participated in electrochemical dissolution, while the volume of the cathode increased uniformly vs. thickness, indicating deposition of material on both the outer and inner surface. The attenuation coefficient vs. thickness, increased for both anode and cathode, indicating surface chemistry changes. For the anode, the attenuation coefficient increased slightly but uniformly, probably due to aluminum oxide formation on the surface of the anode. For the cathode, the attenuation coefficient increased more than for the anode and nonuniformly. The higher increase in the attenuation coefficient for the cathode is due to precipitation of aluminum hydroxide on the electrode surface, which added hydrogen. Image analysis also showed that, although the attenuation coefficient increased throughout the thickness of the electrode, most of the hydroxide deposition occurred on the outer surface. Energy analysis showed that porous electrodes can be used to reduce process energy requirements by as much as 4 times compared to solid electrodes.
Liquid organic hydrogen carriers (LOHCs) are an attractive fuel source due to their compatibility with existing transportation methods and ease of use. However, they suffer from sluggish (de)hydrogenation kinetics. One promising platform for developing next-generation catalysts is metal–organic frameworks (MOFs), which can enable systematic interrogation into the influence of metal identity and spatial arrangement. In this study, the effect of the coordination environment was investigated using Ni- and Co-based azolate MOFs: MFU-4l-OH (M x Zn 5–x (OH) 4 (BTDD) 3 ; x = 4 for M = Co and x = 3 for M = Ni, H 2 BTDD = bis(1H-1,2,3-triazolo[4,5-b][4′,5′-i])dibenzo[1,4]dioxin), composed of single-site nodes, and M(OH) 2 BBTA (M = Ni, Co; H 2 BBTA = 1H,5H-benzo(1,2-d:4,5-d’)bistriazole), composed of extended chain-type nodes. The catalysts were characterized by isotherms, powder X-ray diffraction (PXRD), scanning electron microscopy (SEM), inductively coupled plasma-optical emission spectroscopy (ICP-OES), and X-ray photoelectron spectroscopy (XPS) analysis. Acetylene hydrogenation activity under steady state conditions (150 °C, 1:1 C 2 H 2 :H 2 ) revealed higher turnover frequencies (TOFs) up 1.17 × 10 –2 mol C 2 H 2 mol −1 Ni min −1 and 1.98 × 10 –3 mol C 2 H 2 mol −1 Co min −1 for Ni-MFU-4l-OH and Co-MFU-4l-OH, respectively, compared to their BBTA analogues. However, the Co-based MOFs, particularly Co 2 (OH) 2 -BBTA, exhibited greater selectivity (up to 19%) for the fully hydrogenated ethane product. Isosteric heat of adsorption (Q st ) measurements for ethylene and ethane revealed that the BBTA framework had stronger binding to the products than MFU-4l. Furthermore, these findings demonstrate that metal identity and coordination environment may modulate acetylene hydrogenation performance, leading to design principles for tuning LOHC hydrogenation catalysts.
Rigorous stakeholder-vetted techno-economic analysis was performed to assess the cost of hydrogen (H 2 ) produced using state-of-the-art Liquid Alkaline (LA) electrolysis. Projected high-volume, untaxed levelized cost of hydrogen (LCOH) range from 2020US $\$ 1.84$ to $\$ 2.88$/kg-H 2 depending on technology year, process design, and electrolyzer project scale, assuming an electricity price of $\$ 0.03$/kWh. The total installed capital cost for a LA electrolysis plant was estimated from bottom-up stack and installed cost models that account for purchased equipment, installation costs, site preparation, and general overhead costs. For this study, the LA electrolysis plant is assumed to be a stick-built, greenfield project developed by an EPC firm with electrolysis stacks purchased directly from an electrolysis stack manufacturer. The price of the electrolysis stacks is based on a bottom-up cost assessment with business markup for the electrolysis company fabricator. Methods from the Hydrogen Analysis (H 2 A) production model, a peer-reviewed national laboratory-developed discounted cash flow model, were used to calculate the LCOH production in 2020$/kg-H 2 . The baseline electricity price case ($\$ 0.03$/kWh) corresponds to average wholesale electricity prices currently possible in U.S. markets with plentiful wind. Similar low-cost electricity pricing is possible from solar Power Purchase Agreements (PPA) although these prices are typically limited by renewable energy capacity factors.
Green hydrogen is emerging as one of the most promising alternatives to replace fossil fuels. While hydrogen gas has a good energy density by weight, its poor energy density by volume requires it to be stored under high pressure for commercial use. The hydrogen infrastructure developed to handle this high pressure hydrogen contains multiple components made of polymers as sealing agents. Although polymers do not react chemically with hydrogen gas, they undergo a mechanical failure when high pressure hydrogen gas is suddenly depressurized. This phenomenon known as 'rapid decompression failure' occurs due to the diffusion of hydrogen through polymer and getting trapped inside pre-existing cavities or voids. Here, in this paper, a continuum mechanics-based coupled diffusion-deformation-damage model was developed to predict the hydrogen distribution, stress distribution, and damage propagation inside the polymer while it undergoes rapid decompression failure. The polymer was modeled as a hyperelastic material because it represents the nonlinear material response observed in uniaxial tensile tests perfectly. The effects of hydrogen diffusivity, pre-existing cavity size, cavity location, applied hydrogen pressure, and depressurization rate on damage initiation were studied. It was found that the coefficient of diffusion plays an important role in damage initiation and damage was mostly concentrated in the inside areas rather than near the surface. Experiments were conducted with EPDM polymer which agreed well with the predicted trends using the given model. The effect of adding carbon black and silica filler particles and plasticizer to the pure EPDM polymer was also studied. It was found that damage during RDF decreases with the addition of fillers, but increases with the addition of the plasticizer. Finally, the damage evolution in the presence of two cavities was also studied, and was found that the interaction of stress fields around the cavities alters the damage occurring during RDF.