Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled problems”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

COG11.3 Abstract

COG is a high-resolution code for the Monte Carlo simulation of coupled particle transport in arbitrary3-D geometry. COG will transport neutrons, protons, deuterons, alpha particles with energies up to hundreds of GeV, and photons with energy ranges limited by the available cross section sets and physics models. Electrons can be transported via the EGS5 electron transport kernel, electrons can also be transported. The COG code is a significant upgrade from earlier Monte Carlo transport codes and has been written specifically to makes it more versatile, accurate, and easy to use. COG has provisions for calculating deep penetration (shielding) problems, criticality problems, and neutron activation problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Fermilab's controls development with virtual accelerator

Control Systems development is often the last thing considered when designing and building new equipment, e.g. a new detector or superconducting RF LINAC; however when the new equipment is installed, it is the first thing desired to be operational for testing. Due to frequent delays in building new equipment and project deadlines, control system development and testing is often curtailed. A way to alleviate this problem is to simulate the control system, though this will be challenging for complex systems.The Fermilab PIP-II (proton improvement plan - II) project is being constructed at Fermilab to deliver $800\,MeV$ protons of $>1\,MW$ beam power to replace the present LINAC for the remainder of the existing accelerator complex. The new LINAC consists of a warm front end (WFE), 23 superconducting RF cryomodules (of 5 types), and a beam transfer line (BTL) to the existing complex.The accelerator physics group has a parallel project to create a digital twin (DT) of the PIP-II accelerator. We have coupled the EPICS controls to this DT and are developing both the DT and EPICS software in parallel. This will allow us to develop the EPICS software framework, the HMIs, sequences, high level physics applications, and other services for use in a fully functional control system.This presentation will detail the work that we have performed to date and show demonstrations of controlling and monitoring the status of the accelerator, as well as future plans for this work.

Hanlet, Pierrick [Fermilab]↗
DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗
DOE OSTI2026

A conservative discontinuous Galerkin algorithm for particle kinetics on smooth manifolds

A novel, conservative discontinuous Galerkin algorithm is presented for particle kinetics on manifolds. The motion of particles on the manifold is represented using both canonical and non-canonical Hamiltonian formulations. Our schemes apply to both formulations, but the canonical formulation results in a particularly efficient scheme that also conserves particle density and energy exactly. The collisionless update is coupled to a Bhatnagar-Gross-Krook (BGK) collision operator that provides a simplified model for relaxation to local thermodynamic equilibrium. An iterative scheme is constructed to ensure collisional invariants (density, momentum and energy) are preserved numerically. Rotation of the manifold is incorporated by modifying the Hamiltonian while ensuring a canonical formulation. Several test problems, including a kinetic version of the classical Sod shock problem, Kelvin-Helmholtz instability on the surfaces of a sphere and a hyperboloid, with and without rotations, are presented. A prospectus for further development of this approach to simulation of kinetic theory in general relativity is presented.

Discontinuous Galerkin↗
DOE OSTI2025

SympGNNs: Symplectic Graph Neural Networks for identifying high-dimensional Hamiltonian systems and node classification

Existing neural network models to learn Hamiltonian systems, such as SympNets, although accurate in low-dimensions, struggle to learn the correct dynamics for high-dimensional many-body systems. Herein, we introduce Symplectic Graph Neural Networks (SympGNNs) that can effectively handle system identification in high-dimensional Hamiltonian systems, as well as node classification. SympGNNs combine symplectic maps with permutation equivariance, a property of graph neural networks. Specifically, we propose two variants of SympGNNs: (i) G-SympGNN and (ii) LA-SympGNN, arising from different parameterizations of the kinetic and potential energy. We demonstrate the capabilities of SympGNN on two physical examples: a 40-particle coupled Harmonic oscillator, and a 2000-particle molecular dynamics simulation in a two-dimensional Lennard-Jones potential. Furthermore, we demonstrate the performance of SympGNN in the node classification task, achieving accuracy comparable to the state-of-the-art. Finally, we also empirically show that SympGNN can overcome the oversmoothing and heterophily problems, two key challenges in the field of graph neural networks.

Deep learning↗
DOE OSTI2025

Moment method and continued fraction expansion in Floquet operator Krylov space

Recursion methods such as Krylov techniques map complex dynamics to an effective noninteracting problem in one dimension. For example, the operator Krylov space for Floquet dynamics can be mapped to the dynamics of an edge operator of the one-dimensional Floquet inhomogeneous transverse field Ising model (ITFIM), where the latter, after a Jordan-Wigner transformation, is a Floquet model of noninteracting Majorana fermions and the couplings correspond to Krylov angles. We present an application of this showing that a moment method exists where given an autocorrelation function, one can construct the corresponding Krylov angles and from that the corresponding Floquet ITFIM. Consequently, when no solutions for the Krylov angles are obtained, it indicates that the autocorrelation is not generated by unitary dynamics. We highlight this by studying certain special cases: stable m-period dynamics derived using the method of continued fractions, exponentially decaying, and power-law decaying stroboscopic dynamics. Remarkably, our examples of stable m-period dynamics correspond to m-period edge modes for the Floquet ITFIM where, deep in the chain, the couplings correspond to a critical phase. Furthermore our results pave the way to engineer Floquet systems with desired properties of edge modes and also provide examples of persistent edge modes in gapless Floquet systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Higgs Inflation and the Electroweak Gauge Sector

Abstract We introduce a method that allows the Higgs to be the inflaton. The Higgs is considered as a pseudo‐Nambu‐Goldstone (pNG) boson of a global coset symmetry , which is spontaneously breaks at an energy scale . A suitable Chern−Simons (CS) interaction is given to it, with representing the dimensionless CS coupling strength and an decay constant. As a result, slow‐roll inflation occurs via ‐induced friction down a steep sinusoidal potential. To obey electroweak symmetry, the lowest‐order CS interaction is required to be quadratic in the Higgs, with the coupling strength . Higher‐order interaction terms keep the full Lagrangian nearly invariant under the approximate pNG shift symmetry. Employing the simplest symmetry coset , ‐folds of inflation occur when . Successfully explaining inflation necessitates small values of the decay constant, ; this in turn requires large , which is ruled out by electric dipole measurements. Although the electroweak hierarchy problem while achieving successful inflation, the real benefit is found in providing a different path to identifying the Higgs as the inflaton, outside the standard modified‐gravity framework.

Alexander, Stephon [Brown Center for Theoretical P↗
DOE OSTI2025

Toward scalable bound-to-resonance extrapolations for few- and many-body systems

In open quantum many-body systems, the theoretical description of resonant states of many particles strongly coupled to the continuum can be challenging. Such states are commonplace in, for example, exotic nuclei and hadrons, and can reveal important information about the underlying forces at play in these systems. In this work, we demonstrate that the complex-augmented eigenvector continuation (CA-EC) method, originally formulated for the two-body problem with uniform complex scaling, can reliably perform bound-to-resonance extrapolations for genuine three-body resonances having no bound subsystems. Here, we first establish that three-body bound-to-resonance extrapolations are possible by benchmarking different few-body approaches, and we provide arguments to explain how the extrapolation works in the many-body case. We furthermore pave the way towards scalable resonance extrapolations in many-body systems by showing that the CA-EC method also works in the Berggren basis, studying a realistic application using the Gamow shell model.

Ab initio calculations↗
DOE OSTI2024

Review of recent activities with MOOSE, an open-source finite element & finite volume multi-fidelity simulation framework

Modeling and simulation are an increasing part of engineering. This is undoubtedly driven by the high costs of constructing experimental facilities, but also enabled by the exponential increase in computing powers over the last decades, which allows computational models to be closer than ever to reality. One of the main drivers for the development of MOOSE is supporting advanced nuclear reactor simulations. A challenging aspect of modeling advanced nuclear reactors is the plurality of physics involved, including neutronics, thermal hydraulics and fuel performance. These physics are all coupled to some extent and are generally solved in a sequential but iterative fashion. The United States (U.S.) national laboratories have been developing MOOSE, an open source multiphysics framework since its inception at the Idaho National Laboratory (INL) in 2008. This framework enables seamless coupling of multiphysics simulations and facilitates the implementation of new physics and material governing laws. It is continuously expanded with novel numerical methods and new pre-implemented physics module. Numerous applications, developed within the Department of Energy (DOE) laboratories, academia, and industry, including outside of nuclear engineering, have been developed to study specialized physics problems. International collaborations are welcome on this open-source modeling and simulation project.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Exploiting a Shortcoming of Coupled-Cluster Theory: The Extent of Non-Hermiticity as a Diagnostic Indicator of Computational Accuracy

The fundamental non-Hermitian nature of the forms of the coupled-cluster (CC) theory widely used in quantum chemistry has usually been viewed as a negative, but the present paper shows how this can be used to an advantage. Specifically, the non-symmetric nature of the reduced one-particle density matrix (in the molecular orbital basis) is advocated as a diagnostic indicator of computational quality. In the limit of the full coupled-cluster theory [which is equivalent to full configuration interaction (FCI)], the electronic wave function and correlation energy are exact within a given one-particle basis set, and the symmetric character of the exact density matrix is recovered. The extent of the density matrix asymmetry is shown to provide a measure of “how difficult the problem is” (like the well-known T 1 diagnostic), but its variation with the level of theory also gives information about “how well this particular method works”, irrespective of the difficulty of the problem at hand. The proposed diagnostic is described and applied to a select group of small molecules, and an example of its overall utility for the practicing quantum chemist is illustrated through its application to the beryllium dimer (Be 2 ). Future application of this idea to excited states, open-shell systems, and symmetry-breaking problems and an extension of the method to the two-particle density are then proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Covariance-Free Bifidelity Control Variates Importance Sampling for Rare Event Reliability Analysis

Multifidelity modeling has been steadily gaining attention as a tool to address the problem of exorbitant model evaluation costs that makes the estimation of failure probabilities a significant computational challenge for complex real-world problems, particularly when failure is a rare event. To implement multifidelity modeling, estimators that efficiently combine information from multiple models/sources are necessary. In past works, the variance reduction techniques of control variates (CV) and importance sampling (IS) have been leveraged for this task. In this paper, we present the CVIS framework—a creative take on a coupled CV and IS estimator for bifidelity reliability analysis. The framework addresses some of the practical challenges of the CV method by using an estimator for the control variate mean and sidestepping the need to estimate the covariance between the original estimator and the control variate through a clever choice for the tuning constant. Furthermore, the task of selecting an efficient IS distribution is also considered, with a view towards maximally leveraging the bifidelity structure and maintaining expressivity. Additionally, a diagnostic is provided that indicates both the efficiency of the algorithm as well as the relative predictive quality of the models utilized. Finally, the behavior and performance of the framework is explored through analytical and numerical examples.

Markov chain Monte Carlo↗
DOE OSTI2025

Cosmological moduli and non-perturbative production of axions

Cosmological moduli generically come to dominate the energy density of the early universe, and thereby trigger an early matter dominated era. Such non-standard cosmological histories are expected to have profound effects on the evolution and production of axion cold dark matter and dark radiation, as well as their prospects for detection. We consider moduli-axion couplings and investigate the early history of the coupled system, considering closely the evolution of the homogeneous modulus field, the back-reaction from the axion, and the energy densities of the two fields. A particular point of interest is the enhancement of axion production from modulus decay, due to tachyonic and parametric resonant instabilities, and the implications of such production on the cosmological moduli problem, axion dark radiation, and the available parameter space for axion dark matter. Using an effective field theory approach, WKB-based semi-analytical analysis, and detailed numerical estimates of the co-evolution of the system, we evaluate the expected decay efficiency of the modulus to axions. The effects of higher-order operators are studied and implications for UV-complete frameworks such as the Large Volume Scenarios in Type IIB string theory are considered in detail.

Axions and ALPs↗
DOE OSTI2025

Pulling Order Back from the Brink of Disorder: Observation of a Nodal-Line Spin Liquid and Fluctuation Stabilized Order in K 2 ⁢IrCl 6

Competing interactions in frustrated magnets can give rise to highly degenerate ground states from which correlated liquidlike states of matter often emerge. The scaling of this degeneracy influences the ultimate ground state, with extensive degeneracies potentially yielding quantum spin liquids, while subextensive or smaller degeneracies yield static orders. A long-standing problem is to understand how ordered states precipitate from this degenerate manifold and what echoes of the degeneracy survive ordering. Here, we use neutron scattering to experimentally demonstrate a new “nodal-line” spin liquid, where spins collectively fluctuate within a subextensive manifold spanning one-dimensional lines in reciprocal space. Realized in the spin-orbit-coupled, face-centered-cubic iridate K 2 ⁢IrCl 6 , we show that the subextensive degeneracy is robust, but remains susceptible to fluctuations or longer-range interactions which cooperate to select a magnetic order at low temperatures. Proximity to the nodal-line spin liquid influences the ordered state, enhancing the effects of quantum fluctuations that in turn act to stabilize the sublattice magnetization through the self-consistent opening of a large spin-wave gap. Our results demonstrate how quantum fluctuations can act counterintuitively in frustrated materials: Even in a case where fluctuations are ineffective at selecting an ordered state from a degenerate manifold, at the brink of the nodal spin liquid, they can act to protect the ordered state and dictate its low-energy physics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Parallel-in-Time Solution of Hyperbolic PDE Systems via Characteristic-Variable Block Preconditioning

We consider the parallel-in-time solution of both linear and nonlinear hyperbolic partial differential equation (PDE) systems in one spatial dimension. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner applied in the characteristic variables of the underlying linear(ized) hyperbolic PDE. This change of variables is motivated by the observation that intervariable coupling between characteristic variables is weak, at least locally where spatio-temporal variations in the eigenvectors of the associated flux Jacobian are sufficiently small, while that between the original variables is not. For an ℓ-dimensional system of PDEs, applying the preconditioner consists of solving a sequence of ℓ scalar linear(ized)-advection-like problems, each associated with a different characteristic wave-speed in the underlying linear(ized) PDE. Furthermore, we approximately solve these linear advection problems using multigrid reduction-in-time (MGRIT); however, any other suitable parallel-in-time method could be used. Numerical examples are shown for the (linear) acoustics equations in heterogeneous media and for the (nonlinear) shallow water equations and Euler equations of gas dynamics with shocks and rarefactions. For many test problems, the solver converges in just a handful of iterations and with mesh-independent convergence rates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Projecting U.S. coastal storm surge risks and impacts with deep learning

Storm surge is one of the deadliest hazards posed by tropical cyclones (TCs), yet assessing its current and future risk is difficult due to the phenomenon’s rarity and physical complexity. Recent advances in artificial intelligence applications to natural hazard modeling suggest a new avenue for addressing this problem. We develop a deep learning storm surge model to efficiently estimate coastal surge risk in the United States from 900 000 synthetic TC events, accounting for projected changes in TC behavior and sea levels. The derived historical 100 year surge (the event with a 1% yearly exceedance probability) agrees well with historical observations and other modeling techniques. When coupled with an inundation model, we find that heightened TC intensities and sea levels by the end of the century result in a 50% increase in population at risk. Key findings include markedly heightened risk in Florida, and critical thresholds identified in Georgia and South Carolina.

RAFT↗
DOE OSTI2024

New nonrenormalization theorem from UV/IR mixing

In this paper, we prove a new nonrenormalization theorem which arises from UV/IR mixing. This theorem and its corollaries are relevant for all four-dimensional perturbative tachyon-free closed string theories which can be realized from higher-dimensional theories via geometric compactifications. As such, our theorem therefore holds regardless of the presence or absence of spacetime supersymmetry and regardless of the gauge symmetries or matter content involved. This theorem resolves a hidden clash between modular invariance and the process of decompactification, and enables us to uncover a number of surprising phenomenological properties of these theories. Chief among these is the fact that certain physical quantities within such theories cannot exhibit logarithmic or power-law running and instead enter an effective fixed-point regime above the compactification scale. This cessation of running occurs as the result of the UV/IR mixing inherent in the theory. These effects apply not only for gauge couplings but also for the Higgs mass and other quantities of phenomenological interest, thereby eliminating the logarithmic and/or power-law running that might have otherwise appeared for such quantities. These results illustrate the power of UV/IR mixing to tame divergences—even without supersymmetry—and reinforce the notion that UV/IR mixing may play a vital role in resolving hierarchy problems without supersymmetry. Published by the American Physical Society 2024

Abel, Steven (ORCID:000000031213907X)↗
DOE OSTI2026

Eddy Covariance Theory: A Review

Eddy covariance (EC), the gold standard for measuring ecosystem scale gas and heat exchanges, has transformed our understanding of the breathing of the biosphere, and thus global change biology. Despite numerous methodological improvements and insights gained from the technique, the community faces persistent challenges that have been present since the first EC measurements. Here, we review the theoretical developments underpinning EC. We present theoretical developments in four important areas that have relevance to EC measurements of the net ecosystem exchanges (NEE) of gases and heat from a single tower: (i) measuring the total vertical flux density, (ii) flux attenuation, (iii) coordinate rotations, and (iv) energy balance closure. Persistent problems with EC measurements, such as the inability to close the energy budget, led us to identify two priorities for revisiting the theory underlying: (i) sensible heat flux calculations, and (ii) constraining the mean vertical wind velocity. We present a framework for improved calculation of sensible heat flux derived from first principles of fluid mechanics and thermodynamics that considers coupled heat and mass transfer so that conservation of both is obeyed. These refinements are motivated by the need for unbiased measurements of energy and mass transfer between the land surface and atmosphere for ecosystem research and to validate satellite observations and land surface models.

ecosystem fluxes↗
DOE OSTI2026

Reimagining How Flood Warnings Can Inform Decision‐Making and Community Actions

Society faces increasingly severe flood hazards, intensifying demand for flood early warning systems (FEWS) that deliver accurate and actionable information. However, most existing FEWS remain prediction‐centric, treating decision‐making as a downstream consumer of hazard forecasts while offering limited support for uncertainty interpretation, risk communication, and real‐world response. This Perspective presents a vision and blueprint for a novel inland FEWS‐decision‐making (FEWS‐DM) framework that repositions decision‐making as an equal partner in the forecasting process—not a passive recipient of its outputs. The framework is built on three tightly coupled, co‐evolving thrusts: Physical Science (T1), which advances flood prediction with quantified uncertainty informed by decision relevance; Human Science (T2), which incorporates psychology, behavior, and cultural and institutional context; and Decision Science (T3), which unifies physical predictions and human factors through principled, utility‐based decision support with end‐to‐end uncertainty management. Rather than treating T1 as a solved problem, FEWS‐DM recognizes that forecast development itself must be shaped by decision needs through continuous bidirectional feedback. We identify key scientific, behavioral, and operational challenges limiting such integration and discuss the enabling role of AI, while emphasizing human‐centered design and community feedback as essential for building trust and improving flood risk management.

54 ENVIRONMENTAL SCIENCES↗