Search NASA⌕ Search

SEARCH · Search NASA

Results for “Closure problem”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Harmonizing direct and indirect anthropogenic land carbon fluxes indicates a substantial missing sink in the global carbon budget since the early 20th century

Inconsistencies in the calculation of the two anthropogenic land flux terms of the global carbon cycle are investigated. The two terms—the direct anthropogenic flux (caused by direct human disturbance in anthromes, currently a carbon source to the atmosphere) and the indirect anthropogenic flux (caused indirectly by human activities that lead to global change and affecting all biomes, currently an atmospheric carbon sink)—are typically calculated independently, resulting in inconsistent underlying assumptions. We harmonize the estimation of the two anthropogenic land flux terms by incorporating previous estimates of these inconsistencies. We recalculate the global carbon budget (GCB) and apply change-point analysis to the cumulative budget imbalance. Cumulative over 1850–2018 (1959–2018), harmonization results in a 13% lesser (4% greater) land use source from anthromes and a 20% (23%) lesser land sink. This recalculation yields a greater non-closure of the GCB, indicating a missing carbon sink averaging 0.65 Pg C year -1 since the early 20th century. The imbalance likely results from a combination of method discontinuity and structural errors in the assessment of the direct anthropogenic land use flux, greater ocean carbon uptake, structural errors in land models, and in how these land terms are quantified for the budget. We caution against overconfidence in considering the GCB a solved problem and recommend further study of methodological discontinuities in budget terms. We strongly recommend studies that quantify the direct and indirect anthropogenic land fluxes simultaneously to ensure consistency, with a deeper understanding of human disturbance and legacy effects in anthromes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Turbulent Magnetic Relaxation in Pulsar Wind Nebulae

We present a model for magnetic energy dissipation in a pulsar wind nebula. A better understanding of this process is required to assess the likelihood that certain astrophysical transients may be powered by the spin-down of a “millisecond magnetar.” Examples include superluminous supernovae, gamma-ray bursts, and anticipated electromagnetic counterparts to gravitational wave detections of binary neutron star coalescence. Our model leverages recent progress in the theory of turbulent magnetic relaxation to specify a dissipative closure of the stationary magnetohydrodynamic (MHD) wind equations, yielding predictions of the magnetic energy dissipation rate throughout the nebula. Synchrotron losses are self-consistently treated. To demonstrate the model’s efficacy, we show that it can reproduce many features of the Crab Nebula, including its expansion speed, radiative efficiency, peak photon energy, and mean magnetic field strength. Unlike ideal MHD models of the Crab (which lead to the so-called σ -problem), our model accounts for the transition from ultra to weakly magnetized plasma flow and for the associated heating of relativistic electrons. We discuss how the predicted heating rates may be utilized to improve upon models of particle transport and acceleration in pulsar wind nebulae. We also discuss implications for the Crab Nebula’s γ -ray flares, and point out potential modifications to models of astrophysical transients invoking the spin-down of a millisecond magnetar.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Interpretable and flexible non-intrusive reduced-order models using reproducing kernel Hilbert spaces

This paper develops an interpretable, non-intrusive reduced-order modeling technique using regularized kernel interpolation. Existing non-intrusive approaches approximate the dynamics of a reduced-order model (ROM) by solving a data-driven least-squares regression problem for low-dimensional matrix operators. Our approach instead leverages regularized kernel interpolation, which yields an optimal approximation of the ROM dynamics from a user-defined reproducing kernel Hilbert space. We show that our kernel-based approach can produce interpretable ROMs whose structure mirrors full-order model structure by embedding judiciously chosen feature maps into the kernel. The approach is flexible and allows a combination of informed structure through feature maps and closure terms via more general nonlinear terms in the kernel. We also derive a computable a posteriori error bound that combines standard error estimates for intrusive projection-based ROMs and kernel interpolants. In conclusion, the approach is demonstrated in several numerical experiments that include comparisons to operator inference using both proper orthogonal decomposition and quadratic manifold dimension reduction.

Data-driven model reduction↗
DOE OSTI2013

RECOVERY ACT - Robust Optimization for Connectivity and Flows in Dynamic Complex Networks

The goal of this project was to study robust connectivity and flow patterns of complex multi-scale systems modeled as networks. Networks provide effective ways to study global, system level properties, as well as local, multi-scale interactions at a component level. Numerous applications from power systems, telecommunication, transportation, biology, social science, and other areas have benefited from novel network-based models and their analysis. Modeling and optimization techniques that employ appropriate measures of risk for identifying robust clusters and resilient network designs in networks subject to uncertain failures were investigated in this collaborative multi-university project. In many practical situations one has to deal with uncertainties associated with possible failures of network components, thereby affecting the overall efficiency and performance of the system (e.g., every node/connection has a probability of partial or complete failure). Some extreme examples include power grid component failures, airline hub failures due to weather, or freeway closures due to emergencies. These are also situations in which people, materials, or other resources need to be managed efficiently. Important practical examples include rerouting flow through power grids, adjusting flight plans, and identifying routes for emergency services and supplies, in the event network elements fail unexpectedly. Solutions that are robust under uncertainty, in addition to being economically efficient, are needed. This project has led to the development of novel models and methodologies that can tackle the optimization problems arising in such situations. A number of new concepts, which have not been previously applied in this setting, were investigated in the framework of the project. The results can potentially help decision-makers to better control and identify robust or risk-averse decisions in such situations. Formulations and optimal solutions of the considered problems need to capture uncertainty and risk using appropriate probabilistic, statistical and optimization concepts. The main difficulty arising in addressing these issues is the dramatic increase in the computational complexity of the resulting optimization problems. This project studied novel models and methodologies for risk-averse network optimization- specifically, network design, network flows and cluster detection problems under uncertainty. The approach taken was to incorporate a quantitative risk measure known as conditional value-at-risk that is widely used in financial applications. This approach presents a viable alternate modeling and optimization framework to chance-constrained optimization and mean-variance optimization, one that also facilitates the detection of risk-averse solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗
DOE OSTI2022

Simulation and Reynolds-averaged Navier-Stokes modeling of a three-component Rayleigh-Taylor mixing problem with thermonuclear burn

Rayleigh-Taylor mixing in the presence of a third component with intermediate density is investigated through three-dimensional large-eddy simulation (LES) with a high-order compact finite-difference code. Two configurations are considered: (1) a symmetric configuration in which the Atwood number between the heavy and intermediate components matches the Atwood number between the intermediate and light components and (2) an asymmetric configuration in which the Atwood number between the heavy and intermediate components is an order of magnitude greater than the Atwood number between the intermediate and light components. Mass fraction covariances are extracted, and proposed Reynolds-averaged Navier-Stokes (RANS) closures for density-specific-volume and density-mass-fraction covariances are evaluated in an a priori fashion. Additionally, a multicomponent extension of the k - Φ - L - a - V RANS model [Morgan, Phys. Rev. E 104, 015107 (2021)] is presented which includes model equations for the upper-triangular elements of the mass fraction covariance matrix. This model, referred to as the k - Φ - L - a - C model, is compared against results from LES and against other RANS models. Profiles of average mass fraction, mass-fraction covariance, and density-specific-volume covariance obtained with the k - Φ - L - a - C model are found to agree well with LES data. Finally, the impact of three-component turbulent mixing on average reaction rate is investigated in both premixed and nonpremixed cases by heating the mixing layer and allowing it to undergo thermonuclear (TN) burn. A closure model for average reaction rate is proposed for use with the k - Φ - L - a - C model, and when this model is applied, improved agreement is obtained between LES and RANS in total TN neutron production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A Fokker–Planck approach to a moment closure for mixing in variable-density turbulence

We develop a theory for the cascade mixing terms in a moment closure approach to binary active scalar mixing in variable-density turbulence. To address the variable-density complications we apply, as a principle and constraint, the conservation of the probability density function (PDF) through a Fokker–Planck equation with bounded sample space whose attractor is the beta PDF with skewness. Mixing is related to a single-point PDF as a realisability principle to provide mathematically rigorous expressions for the small scale statistics in terms of largescale moments. The problem of the unknown small-scale mixing is replaced with the determination of the drift and diffusion terms of a Fokker–Planck equation in a beta-PDF-convergent stochastic process. We find that realisability of a beta-convergent process requires the mixing time-scale ratio, taken as a constant in passive scalar mixing, to be a function of the mean mass fraction, mean fluid density, the Atwood number, the density-volume correlation and moments of the density field. We develop and compare the new model with direct numerical simulations data of non-stationary homogeneous variable-density turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

2025 Advances in NekRS: Supporting improved performance for nuclear applications

This report presents several 2025 advancements in NekRS, a high-fidelity spectral element CFD code developed at Argonne National Laboratory to support the NEAMS thermal-hydraulics program. The forthcoming v25 release consolidates several of these advances, adding new features for portability across heterogeneous GPU architectures, real-time in situ visualization, improved turbulence modeling, and conjugate heat transfer coupling. Over the past year, NekRS has demonstrated strong scalability and performance on DOE’s leading exascale platforms, including Aurora and Frontier, confirming its readiness for some of the largest and most complex simulations attempted to date. These achievements provide a powerful new platform for high-fidelity data generation, which in turn supports the development and validation of advanced closure models critical for reactor safety and design. Significant algorithmic innovations have also been introduced. A new global runtime h-refinement capability simplifies workflows by reducing mesh preparation burdens and enabling coarse-to-fine restarts. Building on this, a novel multigrid strategy was implemented to accelerate pressure and transport solves at scale, addressing long-standing bottlenecks in exascale CFD. Together, these developments improve both the efficiency and accessibility of high-fidelity simulations for reactor-relevant problems. Collectively, these enhancements represent a major step forward in simulation technology, positioning NekRS as a cornerstone of NEAMS efforts to enable accurate, efficient, and scalable high-fidelity analysis of advanced nuclear systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Impulse position control algorithms for nonlinear systems

The article is devoted to the formalization and description of impulse-sliding regime in nonlinear dynamical systems that arise in the application of impulse position controls of a special kind. The concept of trajectory impulse-sliding regime formalized as some limiting network element Euler polygons generated by a discrete approximation of the impulse position control This paper differs from the previously published papers in that it uses a definition of solutions of systems with impulse controls, it based on the closure of the set of smooth solutions in the space of functions of bounded variation. The need for the study of such regimes is the fact that they often arise when parry disturbances acting on technical or economic control system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Public involvement on closure of Asse II radioactive waste repository in Germany

From 1967 to 1978, about 125,800 barrels of low- and intermediate level waste were disposed of - nominally for research purposes - in the former 'Asse' salt mine which had before been used for the production of potash for many years. Since 1988 an inflow of brine is being observed which will cause dangers of flooding and of a collapse due to salt weakening and dissolution if it should increase. Since several years the closure of the Asse repository is planned with the objective to prevent the flooding and collapse of the mine and the release of radioactive substances to the biosphere. The first concept that was presented by the former operator, however, seemed completely unacceptable to regional representatives from politics and NGOs. Their activities against these plans made the project a top issue on the political agenda from the federal to the local level. The paper traces the main reasons which lead to the severe safety problems in the past as well as relevant changes in the governance system today. A focus is put on the process for public involvement in which the Citizens' Advisory Group 'A2B' forms the core measure. Its structure and framework, experience and results, expectations from inside and outside perspectives are presented. Furthermore the question is tackled how far this process can serve as an example for a participatory approach in a siting process for a geological repository for high active waste which can be expected to be highly contested in the affected regions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Initial use of Nek5000/Cardinal to improve closure models in Pronghorn

Heat transfer coefficient closure models for pebble bed reactors are built using a data-driven approach by leveraging the capabilities of an Evolutionary Algorithm entitled Particle Swarm Optimization (PSO). In the present work, the Computational Fluid Dynamics code nekRS was used in order to collect the high-fidelity flow data for a core with 1,568 pebbles. To characterize the heat transfer, multiple concentric regions were considered to extract the physical quantities of interest, e.g./ the Reynolds number. The PSO algorithm is employed as part of an inverse problem targeting determine what are the coefficients for a Nusselt number correlation to match the collected data. Such correlation should follow any given format that is defined a priori. Finally, two correlations are proposed, one with an implicit dependence on the pebbles’ wall temperatures and another expressed as a fully explicit correlation depending on the flow conditions and the position within the core. Anyway, given the generic nature of the proposed approach, correlations following different formats could be tested. Preliminary results for the high-fidelity simulation of a fast MSR core are presented. The target Reynolds number is currently 20K, with the expectation that this will increase, pending the availability of further computational resources. These simulations will be used to inform lower fidelity models, including a coarse CFD turbulence model in Pronghorn. Additionally, they will serve as a reference for the RANS models in Nek5000/NekRS.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Improvement in the Modeled Representation of North American Monsoon Precipitation Using a Modified Kain–Fritsch Convective Parameterization Scheme

A commonly noted problem in the simulation of warm season convection in the North American monsoon region has been the inability of atmospheric models at the meso-β scales (10 s to 100 s of kilometers) to simulate organized convection, principally mesoscale convective systems. With the use of convective parameterization, high precipitation biases in model simulations are typically observed over the peaks of mountain ranges. To address this issue, the Kain–Fritsch (KF) cumulus parameterization scheme has been modified with new diagnostic equations to compute the updraft velocity, the convective available potential energy closure assumption, and the convective trigger function. The scheme has been adapted for use in the Weather Research and Forecasting (WRF). A numerical weather prediction-type simulation is conducted for the North American Monsoon Experiment Intensive Observing Period 2 and a regional climate simulation is performed, by dynamically downscaling. In both of these applications, there are notable improvements in the WRF model-simulated precipitation due to the better representation of organized, propagating convection. The use of the modified KF scheme for atmospheric model simulations may provide a more computationally economical alternative to improve the representation of organized convection, as compared to convective-permitting simulations at the kilometer scale or a super-parameterization approach.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Interpreting mixing relationships in energetic melts to estimate vapor contribution and composition

Heterogeneous, glassy debris can be produced by energetic events such as near surface nuclear explosions and asteroid impacts. In these events, the resultant debris comprises a rapidly-cooled mixture of the projectile or bomb material with proximate materials such as soil and anthropogenic structures. Glasses formed in these events preserve inter and/or intra sample compositional variation reflecting complex mixtures of the precursor materials or their chemically fractionated products. Understanding formation mechanisms and chemical processes relating to energetic glass formation based on study of the end products requires an objective method for calculating the precursor components that mixed together to produce the final product. Here, we compare a set of multivariate linear-least squares approaches to address the unmixing problem of heterogeneous glasses including classical least squares (CLS), principal component analysis (PCA) and two multivariate curve resolution – alternating least squares (MCR-ALS) approaches. We generate a synthetic data set representing spatially resolved, compositional data in a set of heterogeneously mixed glasses to understand model sensitivity to precursor abundance and composition. We explore the impact of incomplete knowledge of the system prior to the formation of the melts on the resulting confidence and accuracy of the solution (i.e., the number, composition, and abundances of the precursor materials). We conclude that a closure and equality constrained MCR-ALS approach succeeds in identifying unrecognized precursor components and is a suitable alternative to CLS approaches when the investigator has incomplete knowledge of the starting precursor chemical compositions. Finally, we demonstrate application of this MCR-ALS approach to measurements of Trinitite.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2015

BHR equations re-derived with immiscible particle effects

Compressible and variable density turbulent flows with dispersed phase effects are found in many applications ranging from combustion to cloud formation. These types of flows are among the most challenging to simulate. While the exact equations governing a system of particles and fluid are known, computational resources limit the scale and detail that can be simulated in this type of problem. Therefore, a common method is to simulate averaged versions of the flow equations, which still capture salient physics and is relatively less computationally expensive. Besnard et al. [1] developed such a model for variable density miscible turbulence, where ensemble-averaging was applied to the flow equations to yield a set of filtered equations. Besnard et al. further derived transport equations for the Reynolds stresses, the turbulent mass flux, and the density-specific volume covariance, to help close the filtered momentum and continuity equations. We re-derive the exact BHR closure equations which include integral terms owing to immiscible effects. Physical interpretations of the additional terms are proposed along with simple models. The goal of this work is to extend the BHR model to allow for the simulation of turbulent flows where an immiscible dispersed phase is non-trivially coupled with the carrier phase.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2013

Eulerian adaptive finite-difference method for high-velocity impact and penetration problems

Owing to the complex processes involved, faithful prediction of high-velocity impact events demands a simulation method delivering efficient calculations based on comprehensively formulated constitutive models. Such an approach is presented herein, employing a weighted essentially non-oscillatory (WENO) method within an adaptive mesh refinement (AMR) framework for the numerical solution of hyperbolic partial differential equations. Applied widely in computational fluid dynamics, these methods are well suited to the involved locally non-smooth finite deformations, circumventing any requirement for artificial viscosity functions for shock capturing. Application of the methods is facilitated through using a model of solid dynamics based upon hyper-elastic theory comprising kinematic evolution equations for the elastic distortion tensor. The model for finite inelastic deformations is phenomenologically equivalent to Maxwell s model of tangential stress relaxation. Closure relations tailored to the expected high-pressure states are proposed and calibrated for the materials of interest. Sharp interface resolution is achieved by employing level-set functions to track boundary motion, along with a ghost material method to capture the necessary internal boundary conditions for material interactions and stress-free surfaces. The approach is demonstrated for the simulation of high velocity impacts of steel projectiles on aluminium target plates in two and three dimensions.

Barton, Philip T.↗
DOE OSTI2014

Convexity properties of images under nonlinear integral operators

Conditions are obtained for the image of a given set under a general completely continuous nonlinear integral operator to have convex closure. These results are used to establish the uniqueness of quasi-solutions of nonlinear integral equations of the first kind and to prove the solvability of equations of the first kind on a dense subset of the right-hand sides. Bibliography: 11 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Parallel expansion of a fuel pellet plasmoid

The problem of the assimilation of a cryogenic fuel pellet injected into a hot plasma is considered. Due to the transparency to ambient particles of the plasmoid, the localised region of high-density plasma created by ionisation of the ablated pellet material, electrons reach a ‘quasiequilibrium’ (QE) state which is characterised by a steady-state on the fastest collisional time scale. The simplified electron kinetic equation of the QE state is solved. Taking a velocity moment of the higher-order electron kinetic equation, which is valid on the expansion time scale, permits a fluid closure, yielding an evolution equation for the macroscopic parameters describing the QE distribution function. In contrast to the Braginskii equations, the closure does not require that electrons have a short mean free path compared with the size of density perturbations, and permits an anisotropic and highly non-Maxwellian distribution function. As the QE distribution function accounts for both trapped and passing electrons, the self-consistent electric potential that causes the expansion can be properly described, in contrast to earlier models of pellet plasmoid expansion with an unbounded potential. The plasmoid expansion is simulated using both a Vlasov model and a cold-fluid model for the ions. During the expansion plasmoid ions and electrons obtain nearly equal amounts of energy; as hot ambient electrons provide this energy in the form of collisional heating of plasmoid electrons, the expansion of a pellet plasmoid is expected to be a potent mechanism for the transfer of energy from electrons to ions on a time scale shorter than that of ion–electron thermalisation.

Physics↗
DOE OSTI2012

AN APPLICATION OF THE WIENER HERMITE EXPANSION TO THE NONLINEAR EVOLUTION OF DARK MATTER

We apply the Wiener Hermite (WH) expansion to the nonlinear evolution of the large-scale structure and obtain an approximate expression for the matter power spectrum in the full order of the expansion. This method allows us to expand any random function in terms of an orthonormal basis in the space of random functions in such a way that the first order of the expansion expresses the Gaussian distribution, and others are the deviations from Gaussianity. It is proved that the WH expansion is mathematically equivalent to the {Gamma}-expansion approach in the renormalized perturbation theory (RPT). While exponential behavior in the high-k limit has been proved for the mass density and velocity fluctuations of dark matter in the RPT, we prove the behavior again in the context of the WH expansion using the result of the standard perturbation theory (SPT). We propose a new approximate expression for the matter power spectrum which interpolates the low-k expression corresponding to the 1-loop level in SPT and the high-k expression obtained by taking a high-k limit of the WH expansion. The validity of our prescription is specifically verified by comparing with the 2-loop solutions of the SPT. The proposed power spectrum agrees with the result of the N-body simulation with accuracy better than 1% or 2% in a range of baryon acoustic oscillation scales, where the wave number is about k = 0.2-0.4 h Mpc{sup -1} at z = 0.5-3.0. This accuracy is comparable to or slightly less than the ones in the closure theory, the fractional difference of which from the N-body result is within 1%. One merit of our method is that the computational time is very short because only single and double integrals are involved in our solution.

79 ASTRONOMY AND ASTROPHYSICS↗