Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “General relativity equations & solutions”

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.

At least 379 records · Page 21

Extended molecular eigenmodes treatment of dipole–dipole NMR relaxation in real fluids

Traditional models of NMR relaxation fail to account for the complex, multi-exponential behavior of the autocorrelation function in realistic systems characterized by soft-interactions and molecules that are chemically and physically complex. Here, in this study, we describe the relative diffusion of the spin dipoles by means of a Fokker–Planck equation that includes an interaction potential of mean force to account for the response of the physical/chemical environment around the dipoles. By numerically solving the Fokker–Planck equation for the diffusion propagator, we estimate dipole–dipole NMR relaxation for like- and unlike-spin systems via its eigenmode solution. We test the model against molecular simulations of diffusing dipoles with harmonic potentials and also validate using experimental longitudinal relaxation data from real systems, including Gd(III)–aqua and Gd(III)–DO3A–butrol complexes, the latter being an important MRI contrast agent. Using this novel approach, we predict both the inner- and outer-shell contributions to the relaxivity rates with excellent accuracy at frequencies relevant to MRI. We also show that, under the appropriate assumptions, our framework naturally recovers the Bloembergen–Purcell–Pound, the Solomon–Bloembergen–Morgan, and the Hwang–Freed models. Our implementation is general and publicly available for application to a broad range of systems.

Pinheiro dos Santos, Thiago J. [Rice Univ., Housto↗

Numerical analysis of turbine blade tip treatments

Three-dimensional solutions of the Navier-Stokes equations for a turbine blade with a turning angle of 180 degrees have been computed, including blade tip treatments involving cavities. The geometry approximates a preliminary design for the GGOT (Generic Gas Oxidizer Turbine). The data presented here will be compared with experimental data to be obtained from a linear cascade using original GGOT blades. Results have been computed for a blade with 1 percent clearance, based on chord, and three different cavity sizes. All tests were conducted at a Reynolds number of 4 x 10 exp 7. The grid contains 39,440 points with 10 spanwise planes in the tip clearance region of 5.008E-04 m. Streamline plots and velocity vectors together with velocity divergence plots reveal the general flow behavior in the clearance region. Blade tip temperature calculations suggest placement of a cavity close to the upstream side of the blade tip for reduction of overall blade tip temperature. The solutions do not account for the relative motion between the endwall and the turbine blade. The solutions obtained are generally consistent with previous work done in this area,

Gopalaswamy, Nath S.↗

A new theory for multistep discretizations of stiff ordinary differential equations: Stability with large step sizes

A large set of variable coefficient linear systems of ordinary differential equations which possess two different time scales, a slow one and a fast one is considered. A small parameter epsilon characterizes the stiffness of these systems. A system of o.d.e.s. in this set is approximated by a general class of multistep discretizations which includes both one-leg and linear multistep methods. Sufficient conditions are determined under which each solution of a multistep method is uniformly bounded, with a bound which is independent of the stiffness of the system of o.d.e.s., when the step size resolves the slow time scale, but not the fast one. This property is called stability with large step sizes. The theory presented lets one compare properties of one-leg methods and linear multistep methods when they approximate variable coefficient systems of stiff o.d.e.s. In particular, it is shown that one-leg methods have better stability properties with large step sizes than their linear multistep counter parts. The theory also allows one to relate the concept of D-stability to the usual notions of stability and stability domains and to the propagation of errors for multistep methods which use large step sizes.

Majda, G.↗

Redshifting the Cosmological Constant in Unimodular Gravity via Nonlinear Quantum Mechanics

The cosmological constant problem represents a profound conflict between quantum field theory and general relativity. Unimodular gravity offers a compelling starting point by de-gravitating the vacuum energy of the Standard Model, but this framework traditionally trades the problem of vacuum energy for a fine-tuning of initial conditions, which manifest as a ``shadow" cosmological constant. In this paper, we resolve this initial conditions problem by proposing a novel modification to gravity based on nonlinear quantum mechanics. We introduce specific state-dependent terms to the Hamiltonian, constructed from expectation values of the metric such as the average Ricci scalar. These terms alter the dynamical equations of gravity such that the shadow energy density associated with unconstrained initial conditions redshifts away with cosmic expansion, rendering it negligible at late times. The resulting cosmology is naturally dominated by matter and radiation without fine-tuning. We demonstrate that this significant infrared modification of gravity is consistent with local and cosmological tests of gravity. We comment on the possibility of testing this solution in cosmological measurements of Newton's constant.

Kaplan, David E. [Johns Hopkins U.; Tokyo U., IPMU↗

Uniform-density Bose-Einstein condensates of the Gross-Pitaevskii equation found by solving the inverse problem for the confining potential

Here, in this work, we consider a “reverse-engineering” approach to construct confining potentials that support exact, constant density kovaton solutions to the classical Gross-Pitaevskii equation (GPE) also known as the nonlinear Schr¨odinger equation (NLSE). In the one-dimensional case, the exact solution is the sum of stationary kink and anti-kink solutions, i.e. a kovaton, and in the overlapping region, the density is constant. In higher dimensions, the exact solutions are generalizations of this wave function. In the absence of self-interactions, the confining potential is similar to a smoothed out finite square well with minima also at the edges. When self-interactions are added, a term proportional to ±gψ*ψ gets added to the confining potential and ±gM, where M is the norm, gets added to the total energy. In the realm of stability analysis, we find (linearly) stable solutions in the case with repulsive self-interactions which also are stable to self-similar deformations. For attractive interactions, however, the minima at the edges of the potential get deeper and a barrier in the center forms as we increase the norm. This leads to instabilities at a critical value of M (related to the number of particles in the BEC). Comparing the stability criteria from Derrick’s theorem and Bogoliubov-de Gennes analysis stability results, we find that both predict stability for repulsive self-interactions and instability at a critical mass M for attractive interactions. However, the numerical analysis gives a much lower critical mass. The numerical analysis shows further that the initial instabilities violate the symmetry x → -x assumed by Derrick’s theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A High-Order Eulerian–Lagrangian Runge–Kutta Finite Volume (EL–RK–FV) Method for Scalar Nonlinear Conservation Laws

Abstract We present a class of high-order Eulerian–Lagrangian Runge–Kutta finite volume methods that can numerically solve Burgers’ equation with shock formations, which could be extended to general scalar conservation laws. Eulerian–Lagrangian (EL) and semi-Lagrangian (SL) methods have recently seen increased development and have become a staple for allowing large time-stepping sizes. Yet, maintaining relatively large time-stepping sizes post shock formation remains quite challenging. Our proposed scheme integrates the partial differential equation on a space-time region partitioned by linear approximations to the characteristics determined by the Rankine–Hugoniot jump condition. We trace the characteristics forward in time and present a merging procedure for the mesh cells to handle intersecting characteristics due to shocks. Following this partitioning, we write the equation in a time-differential form and evolve with Runge–Kutta methods in a method-of-lines fashion. High-resolution methods such as ENO and WENO-AO schemes are used for spatial reconstruction. Extension to higher dimensions is done via dimensional splitting. Numerical experiments demonstrate our scheme’s high-order accuracy and ability to sharply capture post-shock solutions with large time-stepping sizes.

Chen, Jiajie↗

Nonlocal Kernel Network (NKN): a Stable and Resolution-Independent Deep Neural Network.

Neural operators have recently become popular tools for designing solution maps between function spaces in the form of neural networks. Differently from classical scientific machine learning approaches that learn parameters of a known partial differential equation (PDE) for a single instance of the input parameters at a fixed resolution, neural operators approximate the solution map of a family of PDEs [6, 7]. Despite their success, the uses of neural operators are so far restricted to relatively shallow neural networks and confined to learning hidden governing laws. In this work, we propose a novel nonlocal neural operator, which we refer to as nonlocal kernel network (NKN), that is resolution independent, characterized by deep neural networks, and capable of handling a variety of tasks such as learning governing equations and classifying images. Our NKN stems from the interpretation of the neural network as a discrete nonlocal diffusion reaction equation that, in the limit of infinite layers, is equivalent to a parabolic nonlocal equation, whose stability is analyzed via nonlocal vector calculus. The resemblance with integral forms of neural operators allows NKNs to capture long-range dependencies in the feature space, while the continuous treatment of node-to-node interactions makes NKNs resolution independent. The resemblance with neural ODEs, reinterpreted in a nonlocal sense, and the stable network dynamics between layers allow for generalization of NKN’s optimal parameters from shallow to deep networks. This fact enables the use of shallow-to-deep initialization techniques [8]. Our tests show that NKNs outperform baseline methods in both learning governing equations and image classification tasks and generalize well to different resolutions and depths.

97 MATHEMATICS AND COMPUTING↗

Revisiting matrix-based inversion of scanning mobility particle sizer (SMPS) and humidified tandem differential mobility analyzer (HTDMA) data

Abstract. Tikhonov regularization is a tool for reducing noise amplification during data inversion. This work introduces RegularizationTools.jl, a general-purpose software package for applying Tikhonov regularization to data. The package implements well-established numerical algorithms and is suitable for systems of up to ∼ 1000 equations. Included is an abstraction to systematically categorize specific inversion configurations and their associated hyperparameters. A generic interface translates arbitrary linear forward models defined by a computer function into the corresponding design matrix. This obviates the need to explicitly write out and discretize the Fredholm integral equation, thus facilitating fast prototyping of new regularization schemes associated with measurement techniques. Example applications include the inversion involving data from scanning mobility particle sizers (SMPSs) and humidified tandem differential mobility analyzers (HTDMAs). Inversion of SMPS size distributions reported in this work builds upon the freely available software DifferentialMobilityAnalyzers.jl. The speed of inversion is improved by a factor of ∼ 200, now requiring between 2 and 5 ms per SMPS scan when using 120 size bins. Previously reported occasional failure to converge to a valid solution is reduced by switching from the L-curve method to generalized cross-validation as the metric to search for the optimal regularization parameter. Higher-order inversions resulting in smooth, denoised reconstructions of size distributions are now included in DifferentialMobilityAnalyzers.jl. This work also demonstrates that an SMPS-style matrix-based inversion can be applied to find the growth factor frequency distribution from raw HTDMA data while also accounting for multiply charged particles. The outcome of the aerosol-related inversion methods is showcased by inverting multi-week SMPS and HTDMA datasets from ground-based observations, including SMPS data obtained at Bodega Marine Laboratory during the CalWater 2/ACAPEX campaign and co-located SMPS and HTDMA data collected at the US Department of Energy observatory located at the Southern Great Plains site in Oklahoma, USA. Results show that the proposed approaches are suitable for unsupervised, nonparametric inversion of large-scale datasets as well as inversion in real time during data acquisition on low-cost reduced-instruction-set architectures used in single-board computers. The included software implementation of Tikhonov regularization is freely available, general, and domain-independent and thus can be applied to many other inverse problems arising in atmospheric measurement techniques and beyond.

54 ENVIRONMENTAL SCIENCES↗

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗

Quantifying the role of the lattice in metal–insulator phase transitions

Many materials exhibit phase transitions at which both the electronic properties and the crystal structure change. Some authors have argued that the change in electronic order is primary, with the lattice distortion a relatively minor side-effect, and others have argued that the lattice distortions play an essential role in the energetics of the transition. In this paper, we introduce a formalism that resolves this long-standing problem. The methodology works with any electronic structure method that produces solutions of the equation of state determining the electronic order parameter as a function of lattice distortion. We use the formalism to settle the question of the physics of the metal–insulator transitions in the rare-earth perovskite nickelates (RNiO 3 ) and Ruddlesden–Popper calcium ruthenates (Ca 2 RuO 4 ) in bulk, heterostructure, and epitaxially strained thin film forms, finding that electron-lattice coupling is key to stabilizing the insulating state in both classes of materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Approximate Turbulent Boundary-Layer Development in Plane Compressible Flow Along Thermally Insulated Surfaces with Application to Supersonic-Tunnel Contour Correction

Numerical solutions of the differential equation obtained from the momentum theorem for the development of a turbulent boundary layer along a thermally insulated surface in two-dimensional and in radial shock-free flow are presented in tabular form for a range of Mach numbers from 0.100 to 10. The solution can be used in a step-wise procedure with any given distribution of favorable pressure gradients and for zero pressure gradients. Solutions are also given for use with moderate adverse pressure gradients. The mean velocity in the boundary layer is approximated by a power-law profile. In view of the stepwise integration methods to be used, the exponent designated the profile shape can be varied along the surface between the integral fraction limits 1/5 and 1/11 through interpolation. Agreement obtained between theoretical and experimental boundary-layer development in a supersonic nozzle at a nominal Mach number of 2 indicates the general validity of the approximations used in the analysis - in particular, the method of extrapolating low-speed skin-friction relations to high Mach number flows. The extrapolation method used assumes that the skin-friction coefficient depend primarily on Reynolds number, provided that the density and the kinematic viscosity are evaluated at surface conditions.

BOUNDARY LAYER - WIND TUNNEL WALLS↗

Current-conserving relativistic linear response for collisional plasmas

Here we investigate the response of a relativistic plasma to electromagnetic fields in the framework of the Boltzmann equation incorporating a collision term in the relaxation rate approximation selected in a form assuring current conservation. We obtain an explicit solution for the linearized perturbation of the Fermi–Dirac equilibrium distribution in terms of the average relaxation rate k. We study the resulting covariant, gauge invariant, and current conserving form of the polarization tensor in the ultrarelativistic and non-relativistic limits. We evaluate the susceptibility in the ultrarelativistic limit and explore their dependence on k. Finally, we study the dispersion relations for the longitudinal and transverse poles of the propagator. We show that for k >2ω p , where ω p is the plasma frequency, the plasma wave modes are overdamped. In the opposite case, k$\ll$ ω p , the propagating plasma modes are weakly damped.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Guidance law development for aeroassisted transfer vehicles using matched asymptotic expansions

This report addresses and clarifies a number of issues related to the Matched Asymptotic Expansion (MAE) analysis of skip trajectories, or any class of problems that give rise to inner layers that are not associated directly with satisfying boundary conditions. The procedure for matching inner and outer solutions, and using the composite solution to satisfy boundary conditions is developed and rigorously followed to obtain a set of algebraic equations for the problem of inclination change with minimum energy loss. A detailed evaluation of the zeroth order guidance algorithm for aeroassisted orbit transfer is performed. It is shown that by exploiting the structure of the MAE solution procedure, the original problem, which requires the solution of a set of 20 implicit algebraic equations, can be reduced to a problem of 6 implicit equations in 6 unknowns. A solution that is near optimal, requires a minimum of computation, and thus can be implemented in real time and on-board the vehicle, has been obtained. Guidance law implementation entails treating the current state as a new initial state and repetitively solving the zeroth order MAE problem to obtain the feedback controls. Finally, a general procedure is developed for constructing a MAE solution up to first order, of the Hamilton-Jacobi-Bellman equation based on the method of characteristics. The development is valid for a class of perturbation problems whose solution exhibits two-time-scale behavior. A regular expansion for problems of this type is shown to be inappropriate since it is not valid over a narrow range of the independent variable. That is, it is not uniformly valid. Of particular interest here is the manner in which matching and boundary conditions are enforced when the expansion is carried out to first order. Two cases are distinguished-one where the left boundary condition coincides with, or lies to the right of, the singular region, and another one where the left boundary condition lies to the left of the singular region. A simple example is used to illustrate the procedure where the obtained solution is uniformly valid to O(Epsilon(exp 2)). The potential application of this procedure to aeroassisted plane change is also described and partially evaluated.

Calise, Anthony J.↗

Physics Mining of Multi-Source Data Sets

Powerful new parallel data mining algorithms can produce diagnostic and prognostic numerical models and analyses from observational data. These techniques yield higher-resolution measures than ever before of environmental parameters by fusing synoptic imagery and time-series measurements. These techniques are general and relevant to observational data, including raster, vector, and scalar, and can be applied in all Earth- and environmental science domains. Because they can be highly automated and are parallel, they scale to large spatial domains and are well suited to change and gap detection. This makes it possible to analyze spatial and temporal gaps in information, and facilitates within-mission replanning to optimize the allocation of observational resources. The basis of the innovation is the extension of a recently developed set of algorithms packaged into MineTool to multi-variate time-series data. MineTool is unique in that it automates the various steps of the data mining process, thus making it amenable to autonomous analysis of large data sets. Unlike techniques such as Artificial Neural Nets, which yield a blackbox solution, MineTool's outcome is always an analytical model in parametric form that expresses the output in terms of the input variables. This has the advantage that the derived equation can then be used to gain insight into the physical relevance and relative importance of the parameters and coefficients in the model. This is referred to as physics-mining of data. The capabilities of MineTool are extended to include both supervised and unsupervised algorithms, handle multi-type data sets, and parallelize it.

Helly, John↗

New heavenly double copies

The double copy relates scattering amplitudes and classical solutions in Yang-Mills theory, gravity, and related field theories. Previous work has shown that this has an explicit realisation in self-dual YM theory, where the equation of motion can be written in a form that maps directly to Plebański’s heavenly equation for self-dual gravity. The self-dual YM equation involves an area-preserving diffeomorphism algebra, two copies of which appear in the heavenly equation. In this paper, we show that this construction is a special case of a wider family of heavenly-type examples, by (i) performing Moyal deformations, and (ii) replacing the area-preserving diffeomorphisms with a less restricted algebra. As a result, we obtain a double-copy interpretation for hyper-Hermitian manifolds, extending the previously known hyper-Kähler case. We also introduce a double-Moyal deformation of the heavenly equation. The examples where the construction of Lax pairs is possible are manifestly consistent with Ward’s conjecture, and suggest that the classical integrability of the gravity-type theory may be guaranteed in general by the integrability of at least one of two gauge-theory-type single copies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigation of vortex breakdown on a delta wing using Euler and Navier-Stokes equations

A numerical investigation of leading-edge vortex breakdown on a delta wing at high angles of attack is presented. The analysis has been restricted to low speed flows on a flat plate wing with sharp leading edges. Both Euler and Navier-Stokes (assuming fully laminar and fully turbulent flows) equations have been used in this study and the results are compared against experimental data. Predictions of vortex breakdown progression with angle of attack with both Euler and Navier-Stokes equations are shown to be consistent with the experimental data. However, the Navier-Stokes predictions show significant improvements in breakdown location at angles of attack where the vortex breakdown approaches the wing apex. The location of the primary vortex and the level of vorticity in the pre-breakdown regions in these flowfield solutions are affected very little by the viscous effects, even though the Navier-Stokes solutions exhibit viscous phenomena such as secondary and tertiary vortices which are not present in the Euler solutions. In the post-breakdown regions, however, the levels of vorticity in the primary vortex have increased differences between the Euler and Navier-Stokes solutions at comparable locations. Navier-Stokes indicates the presence of a secondary vortex even after the primary vortex is burst. The predicted trajectories of the primary vortex are in very good agreement with the test data, the laminar solutions providing the overall best comparison. The Euler shows a small displacement of the primary vortex, relative to experiment, due to the lack of secondary vortices. The turbulent Navier-Stokes, in general, fall between the Euler and laminar solutions. These findings are based on solutions from meshes that are not usually considered fine enough for resolving vortical flows. To further understand the vortex breakdown phenomenon and the fine details of the vortical flow structure, further calculations with the finer meshes and improved turbulence models are necessary.

S. Agrawal↗

Containerless Ripple Turbulence

One of the longest standing unsolved problems in physics relates to the behavior of fluids that are driven far from equilibrium such as occurs when they become turbulent due to fast flow through a grid or tidal motions. In turbulent flows the distribution of vortex energy as a function of the inverse length scale [or wavenumber 'k'] of motion is proportional to 1/k(sup 5/3) which is the celebrated law of Kolmogorov. Although this law gives a good description of the average motion, fluctuations around the average are huge. This stands in contrast with thermally activated motion where large fluctuations around thermal equilibrium are highly unfavorable. The problem of turbulence is the problem of understanding why large fluctuations are so prevalent which is also called the problem of 'intermittency'. Turbulence is a remarkable problem in that its solution sits simultaneously at the forefront of physics, mathematics, engineering and computer science. A recent conference [March 2002] on 'Statistical Hydrodynamics' organized by the Los Alamos Laboratory Center for Nonlinear Studies brought together researchers in all of these fields. Although turbulence is generally thought to be described by the Navier-Stokes Equations of fluid mechanics the solution as well as its existence has eluded researchers for over 100 years. In fact proof of the existence of such a solution qualifies for a 1 M$ millennium prize. As part of our NASA funded research we have proposed building a bridge between vortex turbulence and wave turbulence. The latter occurs when high amplitude waves of various wavelengths are allowed to mutually interact in a fluid. In particular we have proposed measuring the interaction of ripples [capillary waves] that run around on the surface of a fluid sphere suspended in a microgravity environment. The problem of ripple turbulence poses similar mathematical challenges to the problem of vortex turbulence. The waves can have a high amplitude and a strong nonlinear interaction. Furthermore, the steady state distribution of energy again follows a Kolmogorov scaling law; in this case the ripple energy is distributed according to 1/k (sup 7/4). Again, in parallel with vortex turbulence ripple turbulence exhibits intermittency. The problem of ripple turbulence presents an experimental opportunity to generate data in a controlled, benchmarked system. In particular the surface of a sphere is an ideal environment to study ripple turbulence. Waves run around the sphere and interact with each other, and the effect of walls is eliminated. In microgravity this state can be realized for over 2 decades of frequency. Wave turbulence is a physically relevant problem in its own right. It has been studied on the surface of liquid hydrogen and its application to Alfven waves in space is a source of debate. Of course, application of wave turbulence perspectives to ocean waves has been a major success. The experiment which we plan to run in microgravity is conceptually straightforward. Ripples are excited on the surface of a spherical drop of fluid and then their amplitude is recorded with appropriate photography. A key challenge is posed by the need to stably position a 10cm diameter sphere of water in microgravity. Two methods are being developed. Orbitec is using controlled puffs of air from at least 6 independent directions to provided the positioning force. This approach has actually succeeded to position and stabilize a 4cm sphere during a KC 135 segment. Guigne International is using the radiation pressure of high frequency sound. These transducers have been organized into a device in the shape of a dodecahedron. This apparatus 'SPACE DRUMS' has already been approved for use for combustion synthesis experiments on the International Space Station. A key opportunity presented by the ripple turbulence data is its use in driving the development of codes to simulate its properties.

Putterman, Seth↗

Technology Strategy Assessment: Findings from Storage Innovations 2030 Bidirectional Hydrogen Storage

Hydrogen is the most common element in the universe, comprising nearly 75% of all normal matter, and it has been used by scientists for centuries, but it was not fully recognized as an element until 1766, when it was isolated by Henry Cavendish. Early work focused on the generation of hydrogen through the oxidation of metals in water, which released hydrogen gas. Hydrogen’s lighter-than-air and flammable properties were immediately used in engines, zeppelins, and as feedstock for a wide variety of chemical reactions. Several approaches were developed for the production of hydrogen with the most common being associated with the production and conversion of hydrocarbon-based fuels. Coal gasification, steam methane reforming, and other reformation processes provide the majority of current hydrogen production due to the relatively low cost of hydrogen produced through these processes. More than 95% of hydrogen production is used for industrial processes rather than energy storage. To facilitate affordable decarbonization of these industrial processes and to advance the use of hydrogen as a fuel in transportation, DOE launched the Hydrogen Shot as part of the Energy Earthshots Initiative. The goal of the Hydrogen Shot is to reduce the cost of clean hydrogen by 80% to $1/kg of clean hydrogen production within one decade (known as the “1 1 1” goal). This is distinct from the Long-Duration Storage Shot, which is the primary focus of this report; however, it is intrinsically linked to bidirectional hydrogen storage. Several important chemical synthesis processes are dependent upon hydrogen, and the production and use of hydrogen is generally driven by its connection to one of these markets. For example, ammonia is one of the most highly produced chemicals in the world and it depends chiefly on hydrogen. Ammonia is primarily used for agricultural fertilizer and is considered to be largely responsible for a doubling of agricultural production per unit of land over the last century. Another one of hydrogen’s primary uses is as a catalyst in petroleum refining during the desulfurization process. Beyond chemical production, hydrogen is used as a reductant in the production of steel and has been demonstrated as a substitute for metallurgical coal in the production of raw iron. It is even used in the hydrogenation reaction for food products to create more shelf-stable semi-solid fats. However, while hydrogen is produced on the order of 100 million metric tons/year globally to feed these industries, more than 95% of hydrogen is produced from hydrocarbons that emit CO2 during the process. Conversely, electrolysis is a process by which electricity is used to separate hydrogen and oxygen in water molecules, usually across a membrane. Hydrogen production via electrolysis lowers the carbon intensity of produced hydrogen when coupled with low-carbon electricity. Currently, global electrolysis capacity is on the order of 1 GW, which equates to about 500 metric tons/day of hydrogen production. To support large-scale industrial decarbonization, capacity will likely need to increase by two to three orders of magnitude. Electrolysis technology is broadly separated into groups that are defined by the electrolyte used, with further subdivision based on the operating characteristics. The majority of commercial electrolyzer systems are based around three main technology groups: liquid alkaline, proton exchange membrane, and solid oxide. Liquid Alkaline (LA) electrolysis is the oldest, most mature, least expensive, and most common commercial technology, with 400 plants in operation by 1902. Its hydrogen output is low relative to the size of the system due to a low current density. LA electrolysis utilizes a liquid potassium hydroxide solution as the electrolyte. Proton exchange membrane (PEM) electrolysis (also known as polymer electrolyte membrane electrolysis), described in 1960, relies on an acid-impregnated polymer membrane as the electrolyte and typically offers three to six times higher hydrogen production per unit cell area than LA electrolysis. Solid oxide electrolysis, or high-temperature electrolysis, utilizes a ceramic cell as the electrolyte and operates on steam rather than liquid water, enabling electrical efficiencies of more than 90%, which is up from 60% with PEM. Two pre-commercial electrolyzer technologies to note are alkaline exchange membrane (AEM) and proton-conducting solid oxide electrolysis cell (SOEC). AEM potentially has the advantages of both LA and PEM technologies in that it is able to use low-cost materials like LA but with the ability to operate at higher output pressures with a smaller footprint like PEM. Proton-conducting SOEC is similar to commercial SOEC, which uses an oxide-conducting ceramic; however, it uses a proton-conducting ceramic that has the potential to operate at lower temperatures and has lower capital costs. Each of these technologies is experiencing a rapid improvement in performance and a reduction in installed cost, and each appears to be well suited to specific applications. Besides differences in the type of electrolyzer used, the main difference in the architecture of bidirectional hydrogen systems is how the hydrogen is stored. Currently, the most cost-effective way to store large amounts of hydrogen gas is underground, such as in large salt caverns that have been hollowed out. These salt caverns are geographically concentrated in small portions of the United States and are not generally near large metropolitan areas; however, other subsurface architectures are being investigated to expand this reach. A more widely deployable option is aboveground pressurized tanks. These systems are about 10 times as expensive because of the materials and safety margins required to hold hydrogen at high pressures. A third option is using materials-based storage, such as liquid organic hydrogen carriers. By reversibly attaching the produced hydrogen to other molecules, it can be stored at near atmospheric pressure and room temperature. This has the potential to reduce the material cost of storage but may result in a reduction in the efficiency of the process because there are both hydrogen uptake and release processes. While materials-based storage has not been used extensively for large-scale hydrogen storage in the past, there is currently significant activity regarding developing materials and processes for use in large-scale hydrogen storage applications. Electrolysis-produced hydrogen offers an unusual opportunity for energy storage applications. Unlike more conventional energy storage approaches, such as batteries, which operate entirely within electrical markets, hydrogen is a valuable product beyond the electric market and can be directed to the most lucrative use. Hydrogen also can be directly converted back to electricity using either a fuel cell or turbine, or it can be sold to other markets, such as chemical synthesis, steel production, or even export. In this way, excess electricity can be upgraded to the most valuable product. Finally, its use can be actively managed between multiple off-takers; for example, local hydrogen storage can provide a specific amount of stored electricity and any excess can be exported to ammonia production. This flexibility is amplified by the fact that hydrogen storage has fully decoupled power and energy components, which allows for affordable scaling options. Together, this allows a substantial amount of creativity to enable the economic utilization of variable power resources while supporting decarbonization of the industry.

08 HYDROGEN↗