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.
Trajectories of constant chemical potential in the P-T plane serve as an integral formulation of London's equation. The trajectories are useful for analysis and synthesis of fountain effect pump performance. A family of trajectories is generated from available numerical codes.
In November 2017, the Transient Reactor Test Facility returned to operation. Since that time, many transient test series have been completed, such as the Transient Heatsink Overpower Response capsule (THOR), the Transient Water Irradiation System for TREAT (TWIST), and Sirius. Each has provided valuable data for materials performance and reactor safety that can be applied in future designs. During each experimental series, detector count rates provided important information on the core behavior during transients. However, a limitation of these data is that variations in the neutron distribution during experiments can cause errors when attempting to infer reactivity evolution from detector signals. Neutron physics codes can be used to compute the flux shape variations. However, this is a poor solution when the experimental data is used for code verification, validation and uncertainty quantification. Indeed, if the output of the code is used both as a reference and to correct what the reference is compared to, the circular dependency limits the quality of the verification, validation and uncertainty quantification approach. To overcome this problem, the autocorrected reactimeter algorithm (ACRA) has been developed. This approach infers a time-dependent reactivity evolution by testing different spatial corrections and selecting the one that minimizes reactivity variations when the core is in a frozen configuration (i.e., when there is no variation in parameters affecting reactivity). However, the scope of this method was limited to transients where there were negligible thermal feedback. Indeed, the core is never in a frozen configuration when the fuel temperature varies during the whole transient. This is our motivation for developing an improved version of the ACRA that does not require frozen configurations. To develop this new algorithm, we need a precise and unbiased implementation of the inverse point kinetic equations (IPKEs) as any error in the reactivity evaluation will be propagated into the choice of the optimal spatial correction. Indeed, the previous reactimeter algorithm would use approximations, such as a negligible flux amplitude derivative, to focus on rapidity. For the numerical validation of ACRA, we aim at absolute error under for reactivity derived from signals similar to the one of this study. In this summary, we test eight different IPKE implementations. Each will process a mockup signal built for this study, similar to those that the future ACRA will process. Each reactivity output will be compared to the reference reactivity that has been used to generate the mockup signal. The implementation minimizing the difference with the reference reactivity will be used in the development of a new ACRA formulation.
Discretized solution algorithms, which find solutions of field equations in a two or three dimensional field, generally use meshes which are fitted to the field boundary to allow convenient formulation of boundary conditions there. A mesh is defined to be the image of a rectangular grid in computational space under a mesh mapping which maps computational space into physical space. It is not necessary that all of computational space be mapped onto the region of interest in physical space. Parts of it can be excised to give a better fit to the boundary. Many different excisions can be made to fit a single boundary; the choice depends on the mesh arrangement desired in the field.
A methodology for surface grid generation in three dimensional space is described. The method solves a Poisson equation for each coordinate on arbitrary surfaces using successive line over-relaxation. The complete surface curvature terms were discretized and retained within the nonhomogeneous term in order to preserve surface definition; there is no need for conventional surface splines. Control functions were formulated to permit control of grid orthogonality and spacing. A method for interpolation of control functions into the domain was devised which permits their specification not only at the surface boundaries but within the interior as well. An interactive surface generation code which makes use of this methodology is currently under development.
The determination of time-varying lift force and pitch moment generated by a purely pitching Earth-entry capsule is investigated. Experiments were conducted in the 12-foot Low-Speed Tunnel at NASA Langley Research Center, testing a range of oscillation frequencies at a zero-mean angle of attack with a pitching amplitude of 10 degrees. In light of these measurements, a closed-form set of analytically derived equations for lift and moment was used to develop a semi-empirical formulation, incorporating empirically determined values from the experimental runs. The equations are grounded in potential flow theory, Theodorsen's classical theory of unsteady aerodynamics, and the Joukowski theorem of conformal mapping. The unsteady aerodynamics generated by the oscillating body are then modeled by constructing frequency response functions, with quasi-steady forces and moments serving as inputs and unsteady forces and moments as outputs. The experimentally determined gain and phase variations characterize the unsteady nature of the flow and the system's response and flow time-lag to input flow parameters for a blunt-body entry vehicle. The final semi-empirical model is validated with a set of parameters beyond the initial test matrix.
The determination of time-varying lift force and pitch moment generated by a purely pitching Earth-entry capsule is investigated. Experiments were conducted in the 12-foot Low-Speed Tunnel at NASA Langley Research Center, testing a range of oscillation frequencies at a zero-mean angle of attack with a pitching amplitude of 10 degrees. In light of these measurements, a closed-form set of analytically derived equations for lift and moment was used to develop a semi-empirical formulation, incorporating empirically determined values from the experimental runs. The equations are grounded in potential flow theory, Theodorsen's classical theory of unsteady aerodynamics, and the Joukowski theorem of conformal mapping. The unsteady aerodynamics generated by the oscillating body are then modeled by constructing frequency response functions, with quasi-steady forces and moments serving as inputs and unsteady forces and moments as outputs. The experimentally determined gain and phase variations characterize the unsteady nature of the flow and the system's response and flow time-lag to input flow parameters for a blunt-body entry vehicle. The final semi-empirical model is validated with a set of parameters beyond the initial test matrix.
Large Language Models (LLMs) show great promise as a powerful tool for scientific literature exploration. However, their effectiveness in providing scientifically accurate and comprehensive answers to complex questions within specialized domains remains an active area of research. Using the field of high-temperature cuprates as an exemplar, we evaluate the ability of LLM systems to understand the literature at the level of an expert. We construct an expert-curated database of 1,726 scientific papers that covers the history of the field, and a set of 67 expert-formulated questions that probe deep understanding of the literature. We then evaluate six different LLM-based systems for answering these questions, including both commercially available closed models and a custom retrieval-augmented generation (RAG) system capable of retrieving images alongside text. Experts then evaluate the answers of these systems against a rubric that assesses balanced perspectives, factual comprehensiveness, succinctness, and evidentiary support. Among the six systems, two using RAG on curated literature outperformed existing closed models across key metrics, particularly in providing comprehensive and well-supported answers. We discuss promising aspects of LLM performances as well as critical short-comings of all the models. The set of expert-formulated questions and the rubric will be valuable for assessing expert level performance of LLM based reasoning systems.
Scaling laws relating storm electrical generator power (and hence lightning flash rate) to charge transport velocity and storm geometry were originally posed by Vonnegut (1963). These laws were later simplified to yield simple parameterizations for lightning based upon cloud top height, with separate parameterizations derived over land and ocean. It is demonstrated that the most recent ocean parameterization: (1) yields predictions of storm updraft velocity which appear inconsistent with observation, and (2) is formally inconsistent with the theory from which it purports to derive. Revised formulations consistent with Vonnegut's original framework are presented. These demonstrate that Vonnegut's theory is, to first order, consistent with observation. The implications of assuming that flash rate is set by the electrical generator power, rather than the electrical generator current, are examined. The two approaches yield significantly different predictions about the dependence of charge transfer per flash on storm dimensions, which should be empirically testable. The two approaches also differ significantly in their explanation of regional variability in lightning observations.
Chirality-induced spin selectivity (CISS) describes how chiral molecules and materials generate spin polarization even at thermal equilibrium. This observation has challenged established principles of microscopic reversibility and Onsager reciprocity. We resolve this paradox by formulating a pseudo-Hermitian quantum framework that separates thermodynamic equilibrium from time-reversal symmetry. Within this approach, structural chirality and electron correlations, irrespective of their microscopic origin, are sufficient to produce CISS observables. Chirality enters through a non-local metric η that couples spin and spatial motion, leading to real spectra, unitary evolution, and thermodynamic consistency. The framework predicts a chirality-induced spin magnetic ordering characterized by a spin-displacement order, which reconciles equilibrium spin polarization with detailed balance and explains the persistence of CISS in materials composed of light elements. We derive generalized Onsager-Casimir relations that respect the observed CISS symmetry, i.e., parity $(\mathscr{P})$-odd and time-reversal $(\mathscr{T})$-odd, but exhibiting $\mathscr{PT}$-even symmetry. This approach establishes a coherent foundation for equilibrium CISS and provides a route to link chemical chirality with measurable spin-to-charge conversion effects.
Damage in the microstructures of energetic materials (EMs), such as propellants and plastic bonded explosives (PBXs), can significantly alter their response to external loads. Both sensitization and desensitization can occur, causing concerns with safety and performance in the field; predictive models that connect damage and the sensitivity of EMs can enable design and provide confidence in their robustness and reliability. However, modeling of damage evolution is challenging for real microstructures of EMs; samples of damaged EMs are difficult to obtain, thereby hindering experiments and direct numerical simulations to determine the sensitivity of EMs at various stages of damage. Here, we develop an approach to generate synthetic, i.e., in silico produced, damaged microstructures for use in simulations to connect damage levels to sensitivity. The development of the present workflow to generate and impose varying levels of damage in microstructures, known as HEDS (Heterogeneous Energetic Material Damage Simulator), begins with a small set of images of damaged PBXs and combines a collection of deep neural network techniques to generate microstructures with varying levels of damage. By making the synthetic microstructures conform closely to those observed in available real, imaged microstructures, we develop an ensemble of damaged microstructures that can be used for in silico shock experiments. HEDS develops these microstructure ensembles as level set fields, which are directly employed in a sharp interface Eulerian hydrocode where shock simulations are performed to quantify the energy release rate from hotspot fields generated in the microstructure. These capabilities can be useful for the analysis and assessment of changes in the sensitivity of EMs and to design formulations that are less susceptible to damage-induced changes in sensitivity and performance.
The streaming operator, which generates a displacement of a particle on a straight line at a constant speed in transport theory, is derived algebraically from a spherical coordinate formulation of Newton’s second law. This derivation leads to an operator that has more partial derivatives than a Cartesian coordinate formulation of the operator. The additional partial derivatives, which are with respect to the normalized velocity variables of a particle, take into account the intrinsic curvature of a ball. Moreover, these partial derivatives mitigate ray effects, which arise when a finite number of normalized velocities (also called directions or discrete ordinates) are used to simulate a continuous S 2 sphere of directions, by rotating the polar axis of the S 2 sphere into the radial direction of the coordinate system. As a consequence of this rotation, the number of actual discrete ordinates is greatly amplified to an enormous number of effective discrete ordinates by a multiplier that is equal to the number of patches that partitions a spherical surface. In addition to the derivation of the streaming operator, we provide in closed form a solution to the system of characteristic equations that is equivalent to the streaming operator. Furthermore, the solution to the system of characteristic equations enables the construction of an integral operator that is the inverse to the streaming operator. Examples in which ray effects are immensely mitigated by spherical coordinates are presented.
This study investigated the effects of aluminum concentration on simulated high level waste glass properties to eventually establish an aluminum limit (as single-component or multiple-component constraints) for glass formulations for high-aluminum Hanford wastes. A test matrix of 25 high-aluminum glasses (20 = Al 2 O 3 = 28.57 wt%) was generated, and the chemical compositions were measured. The following properties were measured: crystal formation after centerline canister cooling, crystallinity as a function of temperature, density, viscosity, electrical conductivity, toxic leaching characteristics using the Toxicity Characteristic Leaching Procedure, product consistency using the Product Consistency Test, and SO 3 solubility. These results are reported here.
Geologic hydrogen (geoH2) is an emerging technology with massive current market interest and distinct potential to change the paradigm of hydrogen production. Two major subsurface processes influence the amount of geoH2 that are available for energy extraction: 1) geochemical reactions of H2O and Fe2+-bearing rocks which can produce hydrogen in the subsurface environment, where 2) various active microbial communities consume hydrogen as an energy source before the hydrogen reaches the surface. The net gain of hydrogen from these two competing processes dictates the production rate of geoH2. A recent study (Templeton et al., 2024) suggested that for most natural geoH2 systems, five orders of magnitude of production rate enhancement are needed to make geoH2 production economical in the near term. Effective enhancement of the production rate requires an in-depth understanding of the two geoH2 processes, in order to promote the H2-generating geochemical processes and suppress the H2-consuming microbial processes. However, current significant knowledge gaps in these two processes hinders the efforts to formulate strategies to enhance geoH2 production. The National Renewable Energy Laboratory (NREL) is a U.S. Department of Energy National Laboratory with the core mission of leading research, innovation, and strategic partnership to deliver solutions for a clean energy based economy. NREL's extensive research portfolio in hydrogen, bioenergy, geothermal, industrial decarbonization, and energy analysis makes us well positioned to conduct interdisciplinary research and facilitate technology deployment in the geoH2 space. In this presentation, we will discuss ongoing geoH2 research and engagement efforts at NREL, including: 1) geochemical investigation to understand the reaction mechanisms and production rate and potential of different source minerals and rocks, 2) microbiological investigation to understand methanogenesis and acetogenesis in the subsurface geoH2 environment, and identify effective inhibitors for these microbial processes, and 3) preliminary analysis for geoH2 production in the State of Minnesota, where abundant Fe-rich rocks for stimulated geoH2 production and ample opportunity to utilize geoH2 in transforming iron and steel industries are currently available.
As part of the Cyclotron Road program, Anthro Energy investigated the formulation, synthesis, and application of advanced supramolecular compounds for lithium-ion batteries. Anthro Energy developed and tested formulations and chemistries to create polymeric materials with improved mechanical strength, safety, and electrochemical performance. These materials were explored for use in next-generation batteries with a variety of form factors. Initially, Anthro Energy tested these materials to create flexible and deformable batteries, which are in high demand for the wearable/flexible electronics industry. In future work, Anthro Energy’s materials may be applied to create rigid structural batteries, or to enable novel chemistries such as silicon anodes or lithium metal anodes. In the process of testing the company’s unique lithium-ion cells, Anthro Energy explored the development and utilization of novel cell components, cell architectures, and testing protocols. These next-generation flexible and structural batteries may eventually enable new categories of electronic devices.
The stability of spinning flexible satellites in a force-free environment was analyzed. The satellite was modeled as a rigid core having attached to it a flexible appendage idealized as a collection of particles (point masses) interconnected by springs. Both Liapunov and Routh-Hurwitz stability procedures are used. In the former, the Hamiltonian of the system, constrained through the angular momentum integral so as to admit complete damping, is used as a testing function. Equations of motion are written using the hybrid coordinate formulation, which readily accepts a modal coordinate transformation ultimately allowing truncation to a level amenable to literal stability analysis. Closed form stability criteria are generated for the first mode of a restricted appendage model lying in a plane containing the system center of mass and orthogonal to the spin axis. The effects of spin on flexible bodies are discussed by considering a very elementary particle model. Control of passively unstable spacecraft is briefly considered.
Using low frequency satellite observations, a direct comparison of the different theories is possible of the structure of the interplanetary medium. Measurements were made by satellite electron detectors in the same region of the interplanetary media where rapid drift type 3 solar radio bursts are generated. Results show, on the basis of this comparison, the general inadequacies of current type 3 burst theories and indicate a needed emphasis on formulating a new theory that is consistent with the low frequency observations. Representative observations from the RAE satellite are shown. The brightness temperature, which is a measure of the energy flux, versus frequency was plotted. Two events, a large burst and a large storm are also shown, each recorded at frequencies of 2.8 MHz and 700 KHz. High frequency curves obtained from ground-based observations are given.
The Pontryagin Maximum Principle is used to formulate the problem of finding optimum atmospheric vehicular reentry trajectories. The optimization problem is that of minimizing an integral which is a function of the state and control variables. The vehicle's motion is assumed to be influenced by a gravitational force and an aerodynamic force. The problem is formulated and the necessary equations are developed simultaneously for three sets of Euler angles. Computational procedures are suggested so that numerical trajectories may be generated.
An algorithm is developed for producing numerical simulations of large motions of a nonuniform cantilever beam in orbit. Special emphasis is given to the effective formulation of equations of motion and to the use of the finite element method to construct modal functions. Finite element methods are used to generate modal functions in such a way as to permit a particularly harmonious relationship to be established between the disciplines of rigid body dynamics and structural analysis. Sufficient information is provided to enable a reader to create, with relatively little effort, his own simulation program; simulation results are reported, both to provide check cases for other investigators and to illustrate certain important facets of the behavior of flexible spacecraft.