Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stein equation”

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.

18 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Numerical solution of singular Lyapunov equations

We consider the numerical solution of large scale singular (continuous-time) Lyapunov equations of the form AX + XA T + BB T = 0, where A is semistable, that is, its spectrum is contained in the left half plane, with the exception of a few semisimple eigenvalues at zero. We also consider the case of a few semisimple eigenvalues on the imaginary axis. We assume that we know these few eigenvalues (zero or imaginary), and that we have or can compute the corresponding invariant subspaces. We use this information to build an appropriate newly proposed subspace on which to project the Lyapunov equations, and then compute a low-rank approximation to the least squares solution. Selected illustrative numerical examples are provided.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Uncertainty quantification of graph convolution neural network models of evolving processes

The application of neural network models to scientific machine learning tasks has proliferated in recent years. In particular, neural networks have proved to be adept at modeling processes with spatial–temporal complexity. Nevertheless, these highly parameterized models have garnered skepticism in their ability to produce outputs with quantified error bounds over the regimes of interest. Hence there is a need to find uncertainty quantification methods that are suitable for neural networks. In this work we present comparisons of the parametric uncertainty quantification of neural networks modeling complex spatial–temporal processes with Hamiltonian Monte Carlo and Stein variational gradient descent and its projected variant. Specifically we apply these methods to graph convolutional neural network models of evolving systems modeled with recurrent neural network and neural ordinary differential equations architectures. We show that Stein variational inference is a viable alternative to Monte Carlo methods with some clear advantages for complex neural network models. For our exemplars, Stein variational interference gave similar pushed forward uncertainty profiles through time compared to Hamiltonian Monte Carlo, albeit with generally more generous variance. As a result, projected Stein variational gradient descent also produced similar uncertainty profiles to the non-projected counterpart, but large reductions in the active weight space were confounded by the stability of the neural network predictions and the convoluted likelihood landscape.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The shock physics of giant impacts: Key requirements for the equations of state

We discuss major challenges in modeling giant impacts between planetary bodies, focusing on the equations of state (EOS). During the giant impact stage of planet formation, rocky planets are melted and partially vaporized. However, most EOS models fail to reproduce experimental constraints on the thermodynamic properties of the major minerals over the required phase space. Here, we present an updated version of the widely-used ANEOS model that includes a user-defined heat capacity limit in the thermal free energy term. Our revised model for forsterite (Mg 2 SiO 4 ), a common proxy for the mantles of rocky planets, provides a better fit to material data over most of the phase space of giant impacts. We discuss the limitations of this model and the Tillotson equation of state, a commonly used alternative model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Quantum Fluctuations and Multifractally enhanced Superconductivity in Disordered Thin Films

The interplay between electron–electron interactions and weak localization (or anti-localization) phenomena in two-dimensional systems can significantly enhance the superconducting transition temperature. For this work, we develop the theory of quantum fluctuations within such multifractally enhanced superconducting states in thin films. In conditions of weak disorder, we employ the Finkel’stein nonlinear sigma model to derive an effective action for the superconducting order parameter and the quasiclassical Green’s function, meticulously accounting for the influence of quantum fluctuations. This effective action, applicable for interactions of any strength, reveals the critical role of well-known collective modes in a dirty superconductor, and its saddle-point analysis leads to modified Usadel and gap equations. These equations comprehensively incorporate the renormalizations stemming from the interplay between interactions and disorder, resulting in the non-trivial energy dependence of the gap function. Notably, our analysis establishes a direct relation between the self-consistent gap equation at the superconducting transition temperature and the known renormalization group equations for interaction parameters in the normal state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Nonlinearities in Black Hole Ringdowns

The gravitational wave strain emitted by a perturbed black hole (BH) ringing down is typically modeled analytically using first-order BH perturbation theory. In this Letter, we show that second-order effects are necessary for modeling ringdowns from BH merger simulations. Focusing on the strain’s (ℓ,m)=(4,4) angular harmonic, here we show the presence of a quadratic effect across a range of binary BH mass ratios that agrees with theoretical expectations. We find that the quadratic (4,4) mode’s amplitude exhibits quadratic scaling with the fundamental (2,2) mode—its parent mode. The nonlinear mode’s amplitude is comparable to or even larger than that of the linear (4,4) mode. Therefore, correctly modeling the ringdown of higher harmonics—improving mode mismatches by up to 2 orders of magnitude—requires the inclusion of nonlinear effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

X‑ray Coherent Diffractive Imaging of Large Helium Nanodroplets Doped with Small Molecules

We report the first X-ray coherent diffractive imaging experiment on molecule-doped helium nanodroplets. It complements previous work, where we reported single-shot X-ray coherent diffractive imaging studies of Xe dopant clusters formed in 4He and 3He droplets. These noble gas clusters were used to visualize the impact of rotational excitation of the droplets on the spatial distribution of atomic dopants within the droplets, and to study the differences and connections between quantum and classical droplet rotational motion. Here, we expand our studies to the molecular dopants CF4, CHF3, CH3CN, and SF6, imaged with 1.5 keV photons. We find multiple Bragg spots in the diffraction patterns of molecule-doped droplets with radii of approximately 600 nm, which provide evidence that molecules form elongated clusters with preferential alignment along the angular momentum axis of the 4He droplets, in agreement with our previous results on the aggregation of Xe clusters on quantum vortices. Real-space reconstructions of molecular dopant cluster density profiles are obtained for droplets with smaller radii of approximately 300 nm. The diffuse images suggest the formation of low-density, potentially porous, molecular clusters upon aggregation at T = 0.4 K in 4He droplets. In the normal fluid 3He droplets, molecules aggregate into loose clusters on the droplets' equator, similar to previous observations for Xe atoms. Time-of-flight mass spectra reveal that the doped helium nanodroplet moieties fragment extensively into constituent atomic ions, producing only a small fraction of molecular fragment ions. The findings are discussed in the context of previously proposed schemes to use He droplets as potential tamper materials for ultrafast X-ray imaging experiments.

Feinberg, AlexandraJ↗
DOE OSTI2025

Discrete fracture network model benchmarks developed and applied in a DECOVALEX-2023 repository performance assessment study

This study presents newly developed benchmarks for modeling flow and transport within discrete fracture networks (DFNs) and useful methods for analyzing the results. The new benchmarks are designed to test modeling approaches for use in probabilistic performance assessment models of deep geologic repositories in fractured rock. The benchmarks simulate flow and transport through a 1 km 3 block of fractured rock. The first simulates migration of a short pulse of tracer through a simple network of four intersecting fractures. The second adds 1089 stochastically generated fractures. The third changes the pulse to a continuous point source. Evaluation of model performance relies on moment analysis and comparison of the results of different models. The expected nondimensional first moment of the conservative tracer for each benchmark is 1. The benchmarks were simulated by teams from Canada, Czechia, Germany, Korea, Sweden, Taiwan, and the United States as part of a DECOVALEX-2023 study (decovalex.org). The teams used various approaches, including explicit DFN modeling, DFN upscaling to an equivalent continuous porous medium (ECPM), and a combination of both methods. Transport mechanisms are modeled using either the advection-dispersion equation or particle tracking. Results demonstrate strong agreement among the models in breakthrough behavior up to the 75th percentile. Significant deviations in first moments and well-clustered outputs led to the identification of inaccuracies in several models. Such findings exemplify the benefit of exercising these benchmarks and using the presented methods to test DFN flow and transport models.

Benchmark↗
DOE OSTI2022

Improving Common PV Module Temperature Models by Incorporating Radiative Losses to the Sky

PV module operating temperature is the second-most important factor influencing PV system yield–after irradiance–and a substantial contributor to uncertainty in energy system yield predictions. Models commonly used to predict operating temperature in system simulations are based on a simplified energy balance that lumps together different heat loss mechanisms–including radiation–and assumes an overall linear behavior. Radiative heat loss to the sky is usually substantial, but modeling it accurately requires additional information about down-welling long-wave radiation or sky temperature and increases the complexity of temperature model equations. In this work we show how radiative losses to the sky can be separated into two parts to improve the accuracy of modeling without additional complexity. We also predict and demonstrate the variation of these losses at different tilt angles and show that the effective view factor is reduced by the non- isotropic distribution of down-welling long-wave radiation. Finally, we demonstrate substantial reduction in bias (MBE) and scatter (RMSE) when the new radiative loss term is added to the Faiman model using one year of measurements at Sandia National Labs.

14 SOLAR ENERGY↗
DOE OSTI2021

MoorPy (Quasi-Static Mooring Analysis in Python)

MoorPy is a quasi-static mooring model and a suite of associated functions for mooring system analysis. The core model supports quasi-static analysis of moored floating systems including any arrangement of mooring lines and floating platforms. It solves the distributed position and tension of each mooring line segment using standard catenary equations. Floating platforms can be represented with linear hydrostatic characteristics. MoorPy automatically computes a floating system's equilibrium state and can be queried to identify a mooring system's nonlinear force-displacement relationships. Linearized stiffness matrices are efficiently computed using semi-analytic Jacobians. MoorPy also includes plotting functions and a library of mooring component property and cost coefficients. MoorPy can be used directly from Python scripts to perform mooring design and analysis tasks, or it can be coupled with other tools to compute quasi-static mooring reactions as part of a larger simulation.

Hall, Mathew↗
DOE OSTI2024

Efficient Modeling of Floating Wind Arrays Including Current Loads and Seabed Bathymetry

Capabilities for modeling the effects of seabed bathymetry and current drag loads on a floating wind farm are now available in an open-source model for quasi-static analysis. In this model, mooring lines and dynamic cables are represented by a quasi-static solver that can quickly represent complex mooring/cabling arrangements and arrays of floating bodies. To account for seabed bathymetry, we expand the model to include a surface mesh that captures changes in water depth over a rectangular grid. We formulate modifications to the catenary equations that capture a mooring line's profile and tensions when contacting a slope seabed. To account for current drag loads on mooring lines and dynamic cables, we formulate a novel technique that rotates the reference frame so that the vector sum of the the weight and the current force are used in the catenary equations, while accounting for the seabed orientation. To complete the system, current drag loads on floating substructures are handled by inclusion of strip-theory drag calculations. These new capabilities are verified by comparing with results from the established offshore dynamics models MoorDyn and OrcaFlex in equivalent steady-state scenarios. The results show very good agreement for both sloped seabeds and current loads. With computation times of the quasi-static model typically under one second, the model additions are a useful capability toward rapidly evaluating a floating wind array's response to environmental loads under realistic site conditions.

bathymetry↗
DOE OSTI2024

Efficient Modeling of Floating Wind Arrays Including Current Loads and Seabed Bathymetry: Preprint

Capabilities for modeling the effects of seabed bathymetry and current drag loads on a floating wind farm are now available in an open-source model for quasi-static analysis. In this model, mooring lines and dynamic cables are represented by a quasistatic solver that can quickly represent complex mooring/cabling arrangements and arrays of floating bodies. To account for seabed bathymetry, we expand the model to include a surface mesh that captures changes in water depth over a rectangular grid. We formulate modifications to the catenary equations that capture a mooring line's profile and tensions when contacting a slope seabed. To account for current drag loads on mooring lines and dynamic cables, we formulate a novel technique that rotates the reference frame so that the vector sum of the the weight and the current force are used in the catenary equations, while accounting for the seabed orientation. To complete the system, current drag loads on floating substructures are handled by inclusion of strip-theory drag calculations. These new capabilities are verified by comparing with results from the established offshore dynamics models MoorDyn and OrcaFlex in equivalent steady-state scenarios. The results show very good agreement for both sloped seabeds and current loads. With computation times of the quasi-static model typically under one second, the model additions are a useful capability toward rapidly evaluating a floating wind array's response to environmental loads under realistic site conditions.

bathymetry↗
DOE OSTI2021

New Perspectives on the Exoplanet Radius Gap from a Mathematica Tool and Visualized Water Equation of State

Recent astronomical observations obtained with the Kepler and TESS missions and their related ground-based follow-ups revealed an abundance of exoplanets with a size intermediate between Earth and Neptune (1 R ⊕ ≤ R ≤ 4 R ⊕ ). A low occurrence rate of planets has been identified at around twice the size of Earth (2 × R ⊕ ), known as the exoplanet radius gap or radius valley. We explore the geometry of this gap in the mass–radius diagram, with the help of a Mathematica plotting tool developed with the capability of manipulating exoplanet data in multidimensional parameter space, and with the help of visualized water equations of state in the temperature–density (T–ρ) graph and the entropy–pressure (s–P) graph. We show that the radius valley can be explained by a compositional difference between smaller, predominantly rocky planets (<2 × R ⊕ ) and larger planets (>2 × R ⊕ ) that exhibit greater compositional diversity including cosmic ices (water, ammonia, methane, etc.) and gaseous envelopes. In particular, among the larger planets (>2 × R ⊕ ), when viewed from the perspective of planet equilibrium temperature (T eq ), the hot ones (T eq ≳ 900 K) are consistent with ice-dominated composition without significant gaseous envelopes, while the cold ones (T eq ≲ 900 K) have more diverse compositions, including various amounts of gaseous envelopes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Long-term impact of light- and elevated temperature-induced degradation on photovoltaic arrays

Low levelized cost of electricity (LCOE) has been identified as critical for widespread adoption of photovoltaics (PV) without subsidies. Maintaining decades-long high-energy production is often an under-recognized opportunity in meeting cost goals because component lifetimes are not fully quantified at the time of manufacture. Whereas certain standardized tests minimize risk of early failure, there is little guidance to quantitatively predict degradation (or lack thereof) over decades, based on accelerated tests. In this article, we move toward bridging the understanding between indoor accelerated tests and outdoor performance data, with the goal of predicting energy yield with enough accuracy to inform financial decisions. Light- and elevated temperature-induced degradation (LETID) in p-type Si modules is analyzed in terms of impact on long-term module performance and thus LCOE. A method to predict the progression of LETID, using fixed kinetic constants and a numerical solution to the basic reaction rate equations, is detailed. Predictions are compared against both published data and that new to this study. These data include both indoor accelerated tests and fielded modules. We use the results in financial models to derive LCOE of modules in different climates with varying amounts of LETID, including uncertainty. Cost models based on the predictions indicate that LETID has a significant and climate-dependent impact on LCOE. Finally, we show that - even given the uncertainties identified in the study - these financial calculations can provide useful guidance to quantify risk based on accelerated test results. The analysis serves as an example of developing a predictive approach to PV reliability using physics of failure.

14 SOLAR ENERGY↗
DOE OSTI2024

Benchmarking Optimizers for Qumode State Preparation with Variational Quantum Algorithms

Quantum state preparation involves preparing a target state from an initial system, a process integral to applications such as quantum machine learning and solving systems of linear equations. Recently, there has been a growing interest in qumodes due to advancements in the field and their potential applications. However there is a notable gap in the literature specifically addressing this area. This paper aims to bridge this gap by providing performance benchmarks of various optimizers used in state preparation with Variational Quantum Algorithms. We conducted extensive testing across multiple scenarios, including different target states, both ideal and sampling simulations, and varying numbers of basis gate layers. Our evaluations offer insights into the complexity of learning each type of target state and demonstrate that some optimizers perform better than others in this context. Notably, the Powell optimizer was found to be exceptionally robust against sampling errors, making it a preferred choice in scenarios prone to such inaccuracies. Additionally, the Simultaneous Perturbation Stochastic Approximation optimizer was distinguished for its efficiency and ability to handle increased parameter dimensionality effectively.

Kan, Shuwen [Fordham University]↗
DOE OSTI2020

Optical follow-up of the neutron star–black hole mergers S200105ae and S200115j

LIGO and Virgo’s third observing run revealed the first neutron star–black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements creating optical/near-infrared ‘kilonova’ emission. The joint gravitational wave and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter, and independently measure the local expansion rate of the Universe. Here, we present the optical follow-up and analysis of two of the only three high-significance NSBH merger candidates detected to date, S200105ae and S200115j, with the Zwicky Transient Facility. The Zwicky Transient Facility observed ~48% of S200105ae and ~22% of S200115j’s localization probabilities, with observations sensitive to kilonovae brighter than –17.5 mag fading at 0.5 mag d –1 in the g- and r-bands; extensive searches and systematic follow-up of candidates did not yield a viable counterpart. We present state-of-the-art kilonova models tailored to NSBH systems that place constraints on the ejecta properties of these NSBH mergers. Here, we show that with observed depths of apparent magnitude ~22 mag, attainable in metre-class, wide-field-of-view survey instruments, strong constraints on ejecta mass are possible, with the potential to rule out low mass ratios, high black hole spins and large neutron star radii.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗