Search NASA⌕ Search

SEARCH · Search NASA

Results for “continuum”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Tori, Klein bottles, and modulo 8 parity/time-reversal anomalies of 2+1d staggered fermions

We study the symmetries of lattice staggered fermions in 2+1d. Using the symmetries, we can place the system on any sheared torus or Klein bottle. These different backgrounds provide diagnostics of various ’t Hooft anomalies associated with the crystalline symmetries. We then compare the lattice model to its continuum limit. The symmetries of the lattice system are mapped in a nontrivial way to the symmetries of the continuum theories. Using this map, we match the ’t Hooft anomalies on the lattice and the continuum. Along the way, we develop a general formalism to study Hamiltonian lattice models on nontrivial, compact, flat spaces.

Algorithms and Theoretical Developments↗
DOE OSTI2025

LSM and CPT

We study a number of 1+1d lattice models with anti-unitary symmetries that simultaneously reflect space and reverse time. Some of these symmetries are anomalous, leading to Lieb-Schultz-Mattis-type constraints, thus excluding a trivially gapped phase. Examples include a mod 8 anomaly in the Majorana chain and various mod 2 anomalies in the spin chain. In some cases, there is an exact, non-anomalous lattice symmetry that flows in the continuum to CPT. In some other cases, the CPT symmetry of the continuum theory is emergent or absent. Depending on the model, the anomaly of the lattice model is matched in the continuum in different ways. In particular, it can be mapped to an emergent anomaly of an emanant symmetry.

Anomalies in Field and String Theories↗
DOE OSTI2025

An adaptive, data-driven multiscale approach for dense granular flows

The accuracy of coarse-grained continuum models of dense granular flows is limited by the lack of high-fidelity closure models for granular rheology. One approach to addressing this issue, referred to as the hierarchical multiscale method, is to use a high-fidelity fine-grained model to compute the closure terms needed by the coarse-grained model. The difficulty with this approach is that the overall model can become computationally intractable due to the high computational cost of the high-fidelity model. In this work, we describe a multiscale modeling approach for dense granular flows that utilizes neural networks trained using high-fidelity discrete element method (DEM) simulations to approximate the constitutive granular rheology for a continuum incompressible flow model. Our approach leverages an ensemble of neural networks to estimate predictive uncertainty that allows us to determine whether the rheology at a given point is accurately represented by the neural network model. Additional DEM simulations are only performed when needed, minimizing the number of additional DEM simulations required when updating the rheology. This adaptive coupling significantly reduces the overall computational cost of the approach while controlling the error. In addition, the neural networks are customized to learn regularized rheological behavior to ensure well-posedness of the continuum solution. We first validate the approach using two-dimensional steady-state and decelerating inclined flows. We then demonstrate the efficiency of our approach by modeling three-dimensional sub-aerial granular column collapse for varying initial column aspect ratios, where our multiscale method compares well with the computationally expensive computational fluid dynamics (CFD)-DEM simulation.

Dense granular flows↗
DOE OSTI2024

The STAG code: A fully relativistic super transition array calculation using Green’s functions

Calculating opacities for a wide range of plasma conditions (i.e. temperature, density, element) requires detailed knowledge of the plasma configuration space and electronic structure. For plasmas composed of heavier elements, relativistic effects are important in both the electronic structure and the details of opacity spectra. We extend our previously described superconfiguration and super transition array capabilities (Gill et al., 2023) to include a fully relativistic formalism. Further, the use of hybrid bound-continuum supershells in our superconfigurations demonstrates the importance of a consistent treatment of bound and continuum electrons in dense plasma opacities, and we expand the discussion of these consequences to include issues associated with equation of state and electron correlations between bound and continuum electrons.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Modelling and analysis of salt-convection effect on oxide reduction process for uranium oxides using smoothed particle hydrodynamics

Oxygen ion transport via molten salt continuum is a key kinetic feature of electrochemical reduction process for uranium oxides in molten salt. During the electrolysis, the evolution of oxygen ion concentration field in molten salt continuum appears from porous metallic layer to open bulk salt. Despite the important role affecting process throughput, its effect coupled with controlled salt movement was not adequately investigated. To capture the reactive oxygen ion transport phenomena, a multidimensional reduction model based on smoothed particle hydrodynamics (SPH) was combined with the advection-diffusion model. A pelletized oxide fuel bed capturing interconnected salt continuum is considered as a reactant material geometry, which allows forced convective delivery of molten salt through the reactant materials with a relatively low pressure head. Here, in this study, reduction time of a pelletized-fuel bed was quantified with respect to key design parameters of the reactant material configuration. The simulation results demonstrated the advantage of salt convection which effectively removes oxygen ion products from the reactant materials and consequently accelerates the reduction process. Moreover, the reduction-time correlation of a single pelletized fuel was developed as a function of relevant dimensionless parameters and was extended to the application for a pelletized-fuel bed. Finally, several insights and implications in the use of a pelletized oxide fuel were presented. The developed SPH-based framework enabled multi-physics analysis of governing dynamics in appreciating and screening candidate reactant material configurations in a cost-effective manner.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods

Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented here in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N -body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Enabling event-by-event precision in γ-ray cascades for neutron-induced reactions

Neutron-induced γ-ray spectra provide key inputs for modern active interrogation applications. A precise modeling of the nuclear reaction and subsequent emission of γ rays is challenging and often impossible due to limitations on evaluated data file formats and nuclear transport simulation codes. We present a framework that addresses these challenges by combining experimental data and reaction-model calculation outputs into an extended candidate version of the Generalized Nuclear Data Structure (GNDS) file, the successor format for the legacy Evaluated Nuclear Data File (ENDF-6). This proposed GNDS hierarchical format contains all the necessary ingredients for inline γ-ray cascade reproduction with event-by-event precision, including continuum–continuum and continuum–discrete transitions following neutron-capture and inelastic neutron scattering reactions. Cascade-event generation based on our approach demonstrates improved energy conservation on an event-by-event basis and permits the use of γ-γ coincidences in applications. This work offers, for the first time, a method to generate neutron-capture and inelastic neutron-scattering γ-ray cascades where energy conservation, correlations, and experimental primaries are fully accounted for.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Universality of a truncated sigma-model

Bosonic quantum field theories, even when regularized using a finite lattice, possess an infinite dimensional Hilbert space and, therefore, cannot be simulated in quantum computers with a finite number of qubits. A truncation of the Hilbert space is then needed and the physical results are obtained after a double limit: one to remove the truncation and another to remove the regulator (the continuum limit). A simpler alternative is to find a model with a finite dimensional Hilbert space belonging to the same universality class as the continuum model (a “qubitization”), so only the space continuum limit is required. A qubitization of the 1 + 1 dimensional asymptotically free O(3) nonlinear σ-model based on ideas of non-commutative geometry was previously proposed [1] and, in this paper, we provide evidence that it reproduces the physics of the σ-model both in the infrared and the ultraviolet regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

How threshold effects in spectroscopic factors influence heavy-ion knockout reactions

A two-decade-old puzzle in heavy-ion one-nucleon knockout reactions is the strong correlation between the reduction factor R S = σ exp /σ t h and the Fermi surface asymmetry ΔS. Theoretical cross sections typically rely on spectroscopic factors (SFs) from shell model (SM) calculations, which neglect continuum coupling effects. Here, we employ the Gamow shell model (GSM), which explicitly incorporates continuum coupling, to compute SFs for p-shell nuclei and predict corresponding theoretical cross sections. Systematic calculations demonstrate that using GSM-derived SFs substantially reduces discrepancies between theoretical and experimental results. This improvement is particularly significant for deeply bound nucleon knockout in nuclei near the dripline, where traditional SM-based calculations fall short. As a result, using GSM SFs, the ratio R s exhibits no pronounced dependence on ΔS. Furthermore, both the ratio of GSM SFs to SM SFs and their corresponding reaction cross sections ratios exhibit a strong ΔS dependence. We have also compared GSM SFs and cross sections with those from the no-core shell model calculations, giving a similar pronounced sensitivity to ΔS. Detailed analysis attributes these correlations to threshold effects for SFs in weakly bound systems. Overall, incorporating continuum coupling via GSM enhances the reliability of SF predictions for exotic, weakly bound nuclei and provides key insights toward resolving the enduring puzzle in heavy-ion knockout reactions from a nuclear structure perspective.

Gamow shell model↗
DOE OSTI2020

Finite scale theory: Predicting nature’s shocks

In this paper, we describe a continuum model that accurately reproduces the experimentally measured structure of physical shocks in a perfect gas. In this work, we begin by presenting a history of shock structure research, theoretical, experimental and numerical, to quantify the significant discrepancies between Navier–Stokes predictions and laboratory measurements. In our first main result, we discuss modifications that generalize the Chapman–Enskog approximation and lead to a new continuum model termed the Finite Scale Equations. In our second main result, we use this continuum model to calculate shock structure with results comparable in accuracy to numerical simulations based on the Boltzmann equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Field control of quasiparticle decay in a quantum antiferromagnet

Dynamics in a quantum material is described by quantized collective motion: a quasiparticle. The single-quasiparticle description is useful for a basic understanding of the system, whereas a phenomenon beyond the simple description such as quasiparticle decay which affects the current carried by the quasiparticle is an intriguing topic. The instability of the quasiparticle is phenomenologically determined by the magnitude of the repulsive interaction between a single quasiparticle and the two-quasiparticle continuum. Although the phenomenon has been studied in several materials, thermodynamic tuning of the quasiparticle decay in a single material has not yet been investigated. Here we show, by using neutron scattering, magnetic field control of the magnon decay in a quantum antiferromagnet RbFeCl 3 , where the interaction between the magnon and continuum is tuned by the field. At low fields where the interaction is small, the single magnon decay process is observed. In contrast, at high fields where the interaction exceeds a critical magnitude, the magnon is pushed downwards in energy and its lifetime increases. Our study demonstrates that field control of quasiparticle decay is possible in the system where the two-quasiparticle continuum covers wide momentum-energy space, and the phenomenon of the magnon avoiding decay is ubiquitous.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Modeling the formation of Sedan Crater using the FLAG and HOSS codes

Numerical modeling of explosion crater formation requires accounting for complex physical processes. Numerical validation of explosion cratering is an important step in modeling and requires experimental data for comparison. Models using discrete elements and continuum models have both benefits and drawbacks to their approaches. In this work, we consider both an arbitrary Lagrangian–Eulerian (ALE) hydrocode and a finite discrete element method (FDEM) approach to modeling the formation of the Sedan crater, the largest human-made crater in the United States. The Sedan crater formed from an underground nuclear detonation in the Nevada desert as part of Project Plowshare. Our models show that the continuum approach of the hydrocode matched well compared to early test time prior to the mound rupture and subsequent fireball venting, when most of the alluvium exhibited fluid behavior. Our FDEM approach matched the final crater dimensions well, after material had settled back into the crater, when material strength and solid mechanics play key roles. Our work shows how leveraging the benefits of multiple numerical approaches can lead to better understanding of complex physical problems, especially problems with limited experimental data. By using a continuum approach to early-time hydrodynamics and an FDEM approach to later-time solid mechanics, we can better understand the different physical regimes of explosion crater formation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Organizing bacterial vortex lattices by periodic obstacle arrays

Recent experiments have shown that the complex spatio-temporal vortex structures emerging in active fluids are susceptible to weak geometrical constraints. This observation poses the fundamental question of how boundary effects stabilize a highly ordered pattern from seemingly turbulent motion. Here we show, by a combination of continuum theory and experiments on a bacterial suspension, how artificial obstacles guide the flow profile and reorganize topological defects, which enables the design of bacterial vortex lattices with tunable properties. To this end, the continuum model is extended by appropriate boundary conditions. Beyond the stabilization of square and hexagonal lattices, we also provide a striking example of a chiral, antiferromagnetic lattice exhibiting a net rotational flow, which is induced by arranging the obstacles in a Kagome-like array. Recent experiments showed that weak geometrical constraints can organize topological defects in turbulent bacterial suspensions. Here, the authors use a continuum model to study the connection between symmetry and stability of emergent vortex patterns and the geometry of constraining pillar arrays.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

The outburst of the changing-look AGN IRAS 23226-3843 in 2019

IRAS 23226-3843 has previously been classified as a changing-look active galactic nucleus (AGN) based on observations taken in the 1990s in comparison to X-ray data (Swift, XMM-Newton, and NuSTAR) and optical spectra taken after a very strong X-ray decline in 2017. In 2019, Swift observations revealed a strong rebrightening in X-ray and UV fluxes. We aimed to study this outburst in greater detail. We took follow-up Swift, XMM-Newton, and NuSTAR observations of IRAS 23226-3843 together with optical spectra (SALT and SAAO 1.9 m telescope) from 2019 until 2021. IRAS 23226-3843 showed a strong X-ray and optical outburst in 2019. It varied in the X-ray continuum by a factor of 5 and in the optical continuum by a factor of 1.6 within two months. This corresponds to a factor of 3 after correction for the host galaxy contribution. The Balmer and Fe II emission-line intensities showed comparable variability amplitudes during the outburst in 2019. The Hα emission-line profiles of IRAS 23226-3843 changed from a blue-peaked profile in the years 1997 and 1999 to a broad double-peaked profile in 2017 and 2019. However, there were no major profile variations in the extremely broad double-peaked profiles despite the strong intensity variations in 2019. One year after the outburst, IRAS 23226-3843 changed its optical spectral type and became a Seyfert type 2 object in 2020. Blue outflow components are present in the optical Balmer lines and in the Fe band in the X-rays. A deep broadband XMM-Newton/NuSTAR spectrum was taken during IRAS 23226-3843’s maximum state in 2019. This spectrum is qualitatively very similar to a spectrum taken in 2017, but by a factor of 10 higher. The soft X-ray band appears featureless. The soft excess is well modeled with a Comptonization model. A broadband fit with a power-law continuum, Comptonized soft excess, and Galactic absorption gives a good fit to the combined EPIC-pn and NuSTAR spectrum. In addition, we see a complex and broadened Fe K emission-line profile in the X-rays. The changing-look character in IRAS 23226-3843 is most probably caused by changes in the accretion rate – based on the short-term variations on timescales of weeks to months.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The contribution of Compton ionization to ultrafast x-ray scattering

We investigate the role of Compton ionization in ultrafast non-resonant x-ray scattering using a molecular model system, which includes the ionization continuum via an orthonormalized plane wave ansatz. Elastic and inelastic components of the scattering signal, as well as coherent-mixed scattering that arises from electron dynamics, are calculated. By virtue of a near-quantitative distinction between scattering related to electronic transitions into bound and continuum states, we demonstrate how Compton ionization contributes to the coherent-mixed component. Analogous to inelastic scattering, the contribution to the coherent-mixed signal is significant and particularly manifests at intermediate and high-momentum transfers. Strikingly, for molecules with inversion symmetry, the exclusion of bound or continuum transitions may lead to the prediction of spurious coherent-mixed signals. We conclude that qualitative and quantitative accuracies of predicted scattering signals on detectors without energy resolution require that elements of the two-electron density operator are used. This approach inherently accounts for all accessible electronic transitions, including ionization.

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

MHD stability trends and improved performance of LHD inward-shifted configurations: The role of the neutral beam current drive and thermal plasma density

The aim of the present study is to analyze the effect of the neutral beam current drive (NBCD), thermal plasma density, and NBI operational regime on the stability of pressure gradient-driven modes (PGDM) and Alfvén eigenmodes (AE) in LHD inward-shifted configurations. The stabilization of n/m = 1/2 PGDM (n toroidal mode and m poloidal mode) is observed in the discharge 167 800 during the co-NBCD phase. The iota profile evolution measured by motional stark effect diagnostic may indicate the iota profile up-shift caused by the co-NBCD can induce a non-resonant transition of the rational surface 1/2 before the mode stabilization. The evolution of the iota profile and continuum gaps in the discharge 167 805 during the ctr-NBCD phase leads to the stabilization of the AE, caused by the narrowing of the continuum gap as the iota profile down-shift. Opposite stability trends are identified for PGDM and AE stability with respect to the thermal plasma density. A larger thermal plasma density (larger thermal β) further enhances PGDM although the continuum gaps are narrower leading to configurations with stable AEs. The linear stability of AEs is analyzed using the gyro-fluid FAR3d code to reproduce the AE stability trends observed in the experiments with respect to the NBCD and thermal plasma density. The analysis of hypothetical scenarios dedicated to study different NBI operational regimes with respect to EP energy, and β and radial density profiles indicate off-axis NBI operation shows a higher EP β threshold to destabilize AEs compared to on-axis configuration. This is explained by the presence of a TAE gap in the inner plasma region, easily destabilized by an on-axis NBI injection. The control of the NBCD and thermal plasma in the discharge 167 800 shows a transitory stabilization of PGDM and AEs, as well as an improved discharge performance identified by an increment of the neutron fluxes.

Varela, Jacobo↗
DOE OSTI2025

Stability analysis of Alfvén eigenmodes excited by ion cyclotron resonance heating on EAST

Alfvén wave instabilities driven by energetic particles are common in fusion devices. Understanding their behavior is essential for the good confinement of fusion plasma. For this purpose, a series of experiments are performed on EAST to investigate the excitation of toroidal Alfvén eigenmodes (TAEs). Experimentally, it was found that AEs with frequencies around 80, 134, 157 kHz are excited by ion cyclotron resonance heating (ICRH) hydrogen minority heating scenario. Moreover, the excitation of AEs is independent of the wall-coating materials, but determined with the generation and confinement of the fast ions. Statistical analysis suggests that the TAEs can only be excited when the ICRH power is larger than 2 MW and the H 98 factor is larger than 1.15. In line with the experiments, a set of simulations using TORIC, ASCOT, and NOVA-C is performed. The simulation results show that with the experimental density profile and safety factor, the measured TAE with frequency of 134 kHz is well reproduced by taking into account the plasma toroidal rotation frequency. However, the appearance of the 80 kHz mode cannot be captured in the simulations. The ICRH-generated fast hydrogen ions mainly have banana orbits and that their perpendicular velocity is much larger than that of the parallel ions. The stability analysis with NOVA-C suggests that the radial gradient of the fast H ion distribution is the only driving source of TAEs. For TAEs locate in the center of Alfvén continuum, the D ion Landau damping is the dominant damping term; For most other TAEs which have intersections with the continuum boundary, the continuum damping becomes the dominant damping mechanism.

Alfvén eigenmodes↗
DOE OSTI2022

Disentangling Alzheimer’s disease neurodegeneration from typical brain ageing using machine learning

Abstract Neuroimaging biomarkers that distinguish between changes due to typical brain ageing and Alzheimer’s disease are valuable for determining how much each contributes to cognitive decline. Supervised machine learning models can derive multivariate patterns of brain change related to the two processes, including the Spatial Patterns of Atrophy for Recognition of Alzheimer’s Disease (SPARE-AD) and of Brain Aging (SPARE-BA) scores investigated herein. However, the substantial overlap between brain regions affected in the two processes confounds measuring them independently. We present a methodology, and associated results, towards disentangling the two. T1-weighted MRI scans of 4054 participants (48–95 years) with Alzheimer’s disease, mild cognitive impairment (MCI), or cognitively normal (CN) diagnoses from the Imaging-based coordinate SysTem for AGIng and NeurodeGenerative diseases (iSTAGING) consortium were analysed. Multiple sets of SPARE scores were investigated, in order to probe imaging signatures of certain clinically or molecularly defined sub-cohorts. First, a subset of clinical Alzheimer’s disease patients (n = 718) and age- and sex-matched CN adults (n = 718) were selected based purely on clinical diagnoses to train SPARE-BA1 (regression of age using CN individuals) and SPARE-AD1 (classification of CN versus Alzheimer’s disease) models. Second, analogous groups were selected based on clinical and molecular markers to train SPARE-BA2 and SPARE-AD2 models: amyloid-positive Alzheimer’s disease continuum group (n = 718; consisting of amyloid-positive Alzheimer’s disease, amyloid-positive MCI, amyloid- and tau-positive CN individuals) and amyloid-negative CN group (n = 718). Finally, the combined group of the Alzheimer’s disease continuum and amyloid-negative CN individuals was used to train SPARE-BA3 model, with the intention to estimate brain age regardless of Alzheimer’s disease-related brain changes. The disentangled SPARE models, SPARE-AD2 and SPARE-BA3, derived brain patterns that were more specific to the two types of brain changes. The correlation between the SPARE-BA Gap (SPARE-BA minus chronological age) and SPARE-AD was significantly reduced after the decoupling (r = 0.56–0.06). The correlation of disentangled SPARE-AD was non-inferior to amyloid- and tau-related measurements and to the number of APOE ε4 alleles but was lower to Alzheimer’s disease-related psychometric test scores, suggesting the contribution of advanced brain ageing to the latter. The disentangled SPARE-BA was consistently less correlated with Alzheimer’s disease-related clinical, molecular and genetic variables. By employing conservative molecular diagnoses and introducing Alzheimer’s disease continuum cases to the SPARE-BA model training, we achieved more dissociable neuroanatomical biomarkers of typical brain ageing and Alzheimer’s disease.

Hwang, Gyujoon↗