Search NASA⌕ Search

SEARCH · Search NASA

Results for “background model”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

kynema-sgf [SWR-20-85]

Kynema-SGF (formerly AMR-Wind), wherein SGF stands for structured-grid fluid dynamics, is a massively parallel, block-structured adaptive-mesh, incompressible flow solver. The codebase was initiated in 2019 from incflo. The solver is built on top of the AMReX library. AMReX library provides the mesh data structures, mesh adaptivity, as well as the linear solvers used for solving the governing equations. Kynema-SGF is actively developed and maintained by a dedicated multi-institutional team from Lawrence Berkeley National Laboratory, National Laboratory of the Rockies, and Sandia National Laboratories. The primary applications for Kynema-SGF are: performing large-eddy simulations (LES) of atmospheric boundary layer (ABL) flows, simulating wind farm turbine-wake interactions using actuator disk or actuator line models for turbines, and as a background solver when coupled with a near-body solver (e.g., Kynema-UGF) with overset methodology to perform blade-resolved simulations of multiple wind turbines within a wind farm. For offshore applications, the ability to model the air-sea interaction effects and its impact on the ABL characteristics is another focus for the code development effort. As with other codes in the Kynema ecosystem, Kynema-SGF shares the following objectives: *an open, well-documented implementation of the state-of-the-art computational models for modeling wind farm flow physics at various fidelities that are backed by a comprehensive verification and validation (V&V) process; *be capable of performing the highest-fidelity simulations of flow fields within wind farms; and *be able to leverage the high-performance leadership class computing facilities available at DOE national laboratories.

Ananthan, Shreyas↗
DOE OSTI2020

Little influence of Arctic amplification on mid-latitude climate

Observations and model simulations show enhanced warming in the Arctic under increasing greenhouse gases, a phenomenon known as the Arctic amplification (AA), that is likely caused by sea-ice loss. AA reduces meridional temperature gradients linked to circulation, thus mid-latitude weather and climate changes have been attributed to AA, often on the basis of regression analysis and atmospheric simulations. However, other modelling studies show only a weak link. This inconsistency may result from deficiencies in separating the effects of AA from those of natural variability or background warming. In this work, using coupled model simulations with and without AA, we show that cold-season precipitation, snowfall and circulation changes over northern mid-latitudes come mostly from background warming. AA and sea-ice loss increase precipitation and snowfall above ~60° N and reduce meridional temperature gradients above ~45° N in the lower–mid troposphere. However, minimal impact on the mean climate is seen below ~60° N, with weak reduction in zonal wind over 50°–70° N and 150–700 hPa, mainly over the North Atlantic and northern central Asia. These results suggest that the climatic impacts of AA are probably small outside the high latitudes, thus caution is needed in attributing mid-latitude changes to AA and sea-ice loss on the basis of statistical analyses that cannot distinguish the impact of AA from other correlated changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

HSW-V v1.0: localized injections of interactive volcanic aerosols and their climate impacts in a simple general circulation model

Abstract. A new set of standalone parameterizations is presented for simulating the injection, evolution, and radiative forcing by stratospheric volcanic aerosols against an idealized Held–Suarez–Williamson (HSW) atmospheric background in the Energy Exascale Earth System Model version 2 (E3SMv2). In this model configuration (HSW with enabled volcanism, HSW-V), sulfur dioxide (SO2) and ash are injected into the atmosphere with a specified profile in the vertical, and they proceed to follow a simple exponential decay. The SO2 decay is modeled as a perfect conversion to a long-living sulfate aerosol which persists in the stratosphere. All three species are implemented as tracers in the model framework and are transported by the dynamical core's advection algorithm. The aerosols contribute simultaneously to local heating of the stratosphere and cooling of the surface by a simple plane-parallel Beer–Lambert law applied on two zonally symmetric radiation broadbands in the longwave and shortwave ranges. It is shown that the implementation parameters can be tuned to produce realistic temperature anomaly signatures of large volcanic events. In particular, results are shown for an ensemble of runs that mimic the volcanic eruption of Mt. Pinatubo in 1991. The design requires no coupling to microphysical subgrid-scale parameterizations and thus approaches the computational affordability of prescribed aerosol forcing strategies. The idealized simulations contain a single isolated volcanic event against a statistically uniform climate, where no background aerosols or other sources of externally forced variability are present. HSW-V represents a simpler-to-understand tool for the development of climate source-to-impact attribution methods.

Hollowed, Joseph P. (ORCID:0000000286581672)↗
DOE OSTI2025

Coupled geomechanical investigation of depletion-induced fault reactivation

Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.

Coupled geomechanical modeling↗
DOE OSTI2026

Inference of the linear matter power spectrum at z = 0 using DESI DR1 Full-Shape data

Measurements of galaxy distributions at large cosmic distances capture clustering from the past. In this study, we use a cosmological model to translate these observations into the present-day galaxy distribution. Specifically, we reconstruct the 3D linear matter power spectrum at redshift z = 0 using Dark Energy Spectroscopic Instrument (DESI) Year 1 (DR1) galaxy clustering data and Cosmic Microwave Background (CMB) observations, assuming the ΛCDM model, and compare it to the result assuming the w 0 w a CDM model. Building on previous state-of-the-art methods, we apply Effective Field Theory (EFT) modelling of the galaxy power spectrum to account for small-scale effects in the 2-point statistics of galaxy data. Implementation of the EFT approach improves the modelling of the galaxy power spectrum, providing a more robust consistency test of the assumed cosmological model. By casting both CMB and galaxy clustering observations, spanning distinct redshift regimes, into k-space, we can identify discrepancies between the datasets of different redshifts, which would indicate potential inaccuracies in the assumed expansion history. While previous studies have shown consistency with ΛCDM, this work extends the analysis with higher-quality data to further test the expansion histories of both ΛCDM and w 0 w a CDM. Our findings show that both ΛCDM and w 0 w a CDM provide consistent fits to the linear matter power spectrum recovered from DESI DR1 data.

cosmological parameters from LSS↗
DOE OSTI2025

HDG-2 & VIC Design Updates

HDG-2 and VIC Design updates. Including background, HDG-1 challenges, HDG-2 Mitigation efforts, Effect of updated ring model, Pressure boundary surface finish, and conclusions. Vendor Irradiation Capsule (VIC) to include, background, design overview, testing and neutronics results, and models.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological implications from two decades of spectroscopic surveys at the Apache Point Observatory

We present the cosmological implications from final measurements of clustering using galaxies, quasars, and Lyα forests from the completed Sloan Digital Sky Survey (SDSS) lineage of experiments in large-scale structure. These experiments, composed of data from SDSS, SDSS-II, BOSS, and eBOSS, offer independent measurements of baryon acoustic oscillation (BAO) measurements of angular-diameter distances and Hubble distances relative to the sound horizon, r d , from eight different samples and six measurements of the growth rate parameter, fσ 8 , from redshift-space distortions (RSD). This composite sample is the most constraining of its kind and allows us to perform a comprehensive assessment of the cosmological model after two decades of dedicated spectroscopic observation. We show that the BAO data alone are able to rule out dark-energy-free models at more than eight standard deviations in an extension to the flat, ΛCDM model that allows for curvature. When combined with Planck Cosmic Microwave Background (CMB) measurements of temperature and polarization, under the same model, the BAO data provide nearly an order of magnitude improvement on curvature constraints relative to primary CMB constraints alone. Independent of distance measurements, the SDSS RSD data complement weak lensing measurements from the Dark Energy Survey (DES) in demonstrating a preference for a flat ΛCDM cosmological model when combined with Planck measurements. The combined BAO and RSD measurements indicate σ 8 =0.85±0.03, implying a growth rate that is consistent with predictions from Planck temperature and polarization data and with General Relativity. When combining the results of SDSS BAO and RSD, Planck, Pantheon Type Ia supernovae (SNe Ia), and DES weak lensing and clustering measurements, all multiple-parameter extensions remain consistent with a ΛCDM model. Regardless of cosmological model, the precision on each of the three parameters, Ω Λ , H 0 , and σ 8 , remains at roughly 1%, showing changes of less than 0.6% in the central values between models. In a model that allows for free curvature and a time-evolving equation of state for dark energy, the combined samples produce a constraint Ω K =-0.0022±0.0022. The dark energy constraints lead to w 0 =-0.909±0.081 and w a =$-0.49_{-0.30}^{+0.35}$, corresponding to an equation of state of w p =-1.018±0.032 at a pivot redshift z p =0.29 and a Dark Energy Task Force Figure of Merit of 94. The inverse distance ladder measurement under this model yields H 0 =68.18±0.79 km s -1 Mpc -1 , remaining in tension with several direct determination methods; the BAO data allow Hubble constant estimates that are robust against the assumption of the cosmological model. In addition, the BAO data allow estimates of H 0 that are independent of the CMB data, with similar central values and precision under a ΛCDM model. Our most constraining combination of data gives the upper limit on the sum of neutrino masses at Σm ν <0.115 eV (95% confidence). Finally, we consider the improvements in cosmology constraints over the last decade by comparing our results to a sample representative of the period 2000-2010. We compute the relative gain across the five dimensions spanned by w, Ω k , Σm ν , H 0 , and σ 8 and find that the SDSS BAO and RSD data reduce the total posterior volume by a factor of 40 relative to the previous generation. Finally, adding again the Planck, DES, and Pantheon SN Ia samples leads to an overall contraction in the five-dimensional posterior volume of 3 orders of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Quantifying scalar field dynamics with DESI 2024 Y1 BAO measurements

Quintessence scalar fields are a natural candidate for evolving dark energy. Unlike the phenomenological w 0 ⁢w a parameterization of the dark energy equation of state, they cannot accommodate the phantom regime of dark energy w⁡(z) <–1, or crossings into the phantom regime. Recent baryon acoustic oscillation (BAO) measurements by the Dark Energy Spectroscopic Instrument (DESI) indicate a preference for evolving dark energy over a cosmological constant, ranging from 2.6⁢σ –3.9⁢σ when fitting to w 0⁢ w a , and combining the DESI BAO measurements with other cosmological probes. In this work, we directly fit three simple scalar field models to the DESI BAO data, combined with cosmic microwave background anisotropy measurements and supernova datasets. We find the best fit model to include a 2–4% kinetic scalar field energy Ω scf,k , for a canonical scalar field with a quadratic or linear potential. However, only the DESY-Y5 supernova dataset combination shows a preference for quintessence over Λ cold dark matter (CDM) at the 95% confidence level. Fitting to the supernova datasets Pantheon, Pantheon+, DES-Y5, and Union3, we show that the mild tension (n σ < 3.4) under Λ ⁢CDM emerges from a BAO preference for smaller values of fractional mass-energy density Ω m < 0.29, while all supernova datasets, except for Pantheon, prefer larger values, Ω m > 0.3. The tension under Λ ⁢CDM remains noticeable (n σ < 2.8), when replacing two of the DESI BAO redshift bins with effective redshifts z eff = 0.51, and z eff = 0.706 with comparable BOSS DR 12 BAO measurements at z eff = 0.51, and z eff = 0.61. Furthermore, canonical scalar fields as dark energy are successful in mitigating that tension.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Cosmological probes of Dark Radiation from Neutrino Mixing

Abstract Models of stepped dark radiation have recently been found to have an important impact on the anisotropies of the cosmic microwave background, aiding in easing the Hubble tension. In this work, we study models with a sector of dark radiation with a step in its abundance, which thermalizes after big bang nucleosynthesis by mixing with the standard model neutrinos. For this, we extend an earlier work which has focused on the background evolution only until the dark sector thermalizes by deriving the full background and perturbation equations of the model and implementing them in an Einstein-Boltzmann solving code. We expound on the behavior of this model, discussing the wide range of parameters that result in interesting and viable cosmologies that dynamically generate dark radiation during a range of epochs. We find that for the strongly self-coupled regime, there is no large cosmological impact for a tight prior on the mass, whereas larger mass ranges allow a smooth interpolation between a behavior close to the ΛCDM cosmological standard model and close to an additional component of strongly self-interacting dark radiation. In the weakly self-coupled regime we find that we can accommodate a parameter space relevant for the neutrino anomalies as well as one relevant to easing the Hubble tension.

Astronomy & Astrophysics↗
DOE OSTI2024

Life-Cycle Assessment Integration into Scalable Open-Source Numerical Models (LiAISON) for Analyzing Emerging Low-Carbon Technologies

Decarbonizing the industrial sector is a significant challenge in achieving a net-zero greenhouse gas (GHG) emissions economy by 2050 and the Paris Agreement, i.e., a global climate change mitigation target of achieving a maximum average temperature change potential of 1.5 degrees Celsius or less by 2100 with respect to pre-industrial levels. In the United States (US), the industrial sector accounts for 23% of total GHG emissions and is home to a number of hard-to-electrify activities. The chemicals subsector has the single largest subsector emissions profile after direct emissions from fossil fuel combustion and leakage from fossil fuel distribution systems. Within the chemicals subsector, many processes depend on hydrogen or ammonia precursors. Decarbonizing these two commodities would contribute significantly to decarbonizing the industrial sector as hydrogen could also be used for low carbon steel production (e.g., hydrogen-based direct reduction of iron) and other industrial applications. Emerging technologies require the application of prospective life cycle assessment (LCA), which can account for technology (foreground) scaling and process improvements via learning-by-doing, among others. In many cases, the future system context (background) in which the technologies are assumed to operate in is equally relevant. Background scenarios generated by integrated assessment models (IAM) can coherently incorporate potential future dynamics of the energy-climate-human-land system. Further, IAM scenarios are harmonized across socioeconomic and climate change mitigation pathways, which facilitates the comparability of prospective LCAs using different IAMs. We introduce an open source prospective LCA framework, the Life-cycle Assessment Integration into Scalable Open-source Numerical models (LiAISON), to analyze the non-linear relationships between technology foreground and the future energy system background across a series of midpoint and resource use metrics. The integration of LCA and IAM data is achieved using prospective environmental Impact assessment (PREMISE). We showcase it by assessing two Power-to-Hydrogen (PtH2) processes, namely Solid Oxide Electrolysis (SOE) and Polymer Electrolyte Membrane Electrolysis (PEME). We compare the technologies to a baseline of hydrogen production via natural gas-based Steam Methane Reforming (SMR) in a US context of multiple energy system and climate change mitigation futures. Besides providing an analysis that specifies the LCA results ranges with temporal and geospatial explicitness across the two technologies, metrics, and impact assessment methods, this research also aims to establish a base framework that can be expanded to use other IAM generated scenarios and US open-source life cycle inventory (LCI) databases. We find that the temporal environmental performance of either technology or their difference to SMR is directly influenced by the underlying background dynamics. Under baseline projections (i.e., no decarbonization goals), neither process reaches parity with the incumbent technology across several environmental metrics. Under the decarbonization scenarios, the underlying sectoral shifts result in declining impacts over time, compared to 2020 levels, except for metal depletion levels, which increase. The background shifts postulate a heavily decarbonized economy and energy system, which help technologies reach parity with SMR between 2040-2050 (RCP2.6) and 2030-2040 (RCP1.9) for global warming. Despite declines across several other metrics over time, neither PtH2 technology break even with SMR by 2100 besides for global warming.

decarbonizing↗
DOE OSTI2024

Life-Cycle Assessment Integration into Scalable Open-Source Numerical Models (LiAISON) for Analyzing Emerging Low-Carbon Technologies

Decarbonizing the industrial sector is a significant challenge in achieving a net-zero greenhouse gas (GHG) emissions economy by 2050 and the Paris Agreement, i.e., a global climate change mitigation target of achieving a maximum average temperature change potential of 1.5 degrees C or less by 2100 with respect to pre-industrial levels. In the United States (US), the industrial sector accounts for 23% of total GHG emissions and is home to a number of hard-to-electrify activities. The chemicals subsector has the single largest subsector emissions profile after direct emissions from fossil fuel combustion and leakage from fossil fuel distribution systems. Within the chemicals subsector, many processes depend on hydrogen or ammonia precursors. Decarbonizing these two commodities would contribute significantly to decarbonizing the industrial sector as hydrogen could also be used for low carbon steel production (e.g., hydrogen-based direct reduction of iron) and other industrial applications. Emerging technologies require the application of prospective life cycle assessment (LCA), which can account for technology (foreground) scaling and process improvements via learning-by-doing, among others. In many cases, the future system context (background) in which the technologies are assumed to operate in is equally relevant. Background scenarios generated by integrated assessment models (IAM) can coherently incorporate potential future dynamics of the energy-climate-human-land system. Further, IAM scenarios are harmonized across socioeconomic and climate change mitigation pathways, which facilitates the comparability of prospective LCAs using different IAMs. We introduce an open source prospective LCA framework, the Life-cycle Assessment Integration into Scalable Open-source Numerical models (LiAISON), to analyze the non-linear relationships between technology foreground and the future energy system background across a series of midpoint and resource use metrics The integration of LCA and IAM data is achieved using prospective environmental Impact assessment (PREMISE). We showcase it by assessing two Power-to-Hydrogen (PtH2) processes, namely Solid Oxide Electrolysis (SOE) and Polymer Electrolyte Membrane Electrolysis (PEME). We compare the technologies to a baseline of hydrogen production via natural gas-based Steam Methane Reforming (SMR) in a US context of multiple energy system and climate change mitigation futures. Besides providing an analysis that specifies the LCA results ranges with temporal and geospatial explicitness across the two technologies, metrics, and impact assessment methods, this research also aims to establish a base framework that can be expanded to use other IAM generated scenarios and US open-source life cycle inventory (LCI) databases. We find that the temporal environmental performance of either technology or their difference to SMR is directly influenced by the underlying background dynamics. Under baseline projections (i.e., no decarbonization goals), neither process reaches parity with the incumbent technology across several environmental metrics. Under the decarbonization scenarios, the underlying sectoral shifts result in declining impacts over time, compared to 2020 levels, except for metal depletion levels, which increase. The background shifts postulate a heavily decarbonized economy and energy system, which help technologies reach parity with SMR between 2040-2050 (RCP2.6) and 2030-2040 (RCP1.9) for global warming. Despite declines across several other metrics over time, neither PtH2 technology break even with SMR by 2100 besides for global warming.

decarbonizing↗
DOE OSTI2024

Anomaly detection with flow-based fast calorimeter simulators

Recently, several normalizing flow-based deep generative models have been proposed to accelerate the simulation of calorimeter showers. Using caloflow as an example, we show that these models can simultaneously perform unsupervised anomaly detection with no additional training cost. As a demonstration, we consider electromagnetic showers initiated by one (background) or multiple (signal) photons. The caloflow model is designed to generate single-photon showers, but it also provides access to the shower likelihood. We use this likelihood as an anomaly score and study the showers tagged as being unlikely. As expected, the tagger struggles when the signal photons are nearly collinear but is otherwise effective. This approach is complementary to a supervised classifier trained on only specific signal models using the same low-level calorimeter inputs. While the supervised classifier is also highly effective at unseen signal models, the unsupervised method is more sensitive in certain regions, and, thus, we expect that the ultimate performance will require a combination of these approaches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Modeling 233 Pa Generation in Thorium-fueled Reactors for Safeguards

Thorium has been considered as a possible alternative to uranium for nuclear fuel for many decades. It is three to four times more abundant in the earth than uranium and produces significantly less long-lived transuranic nuclear waste. Some claim thorium poses fewer proliferation concerns than other fuel types largely due to 232 U buildup (and associated high energy gamma-emitting decay products) in the irradiated thorium fuel. However, to fully explore potential proliferation concerns, generation and subsequent decay of 233 Pa produced in the reactor core still must be studied. With its half-life of 27 days, 233 Pa decays to 233 U, which is an International Atomic Energy Agency (IAEA) defined special fissionable material that can be used for nuclear weapons production. With more research being dedicated to thorium-fueled reactors, and several of these reactor designs possessing online fuel processing (allowing for on-site protactinium separation), it is important to understand this potential proliferation pathway. In particular, it is theoretically possible to extract protactinium from the irradiated fuel salt before it decays into 233 U. This hypothetical potential diversion can become an even greater proliferation concern if the extracted protactinium is purified through a second separation of protactinium approximately ten days later to remove the short half-life decay products of 232 Pa and 234 Pa, thus resulting in a higher concentration of the 233 Pa isotope, which decays into weapons usable 233 U with hardly any 232 U or 234 U in it. To estimate the concern of this potential proliferation challenge of thorium, different nuclear material accountancy techniques were reviewed for their viability to quantify 233 Pa if extracted from used thorium fuel. Characteristics of interest included technology maturity, cost, precision, and time taken to acquire results. Some technologies, like hybrid K-edge densitometry and passive gamma spectroscopy, appear to be viable techniques based on current literature. Due to the limited scope of this project, only passive gamma spectroscopy was further investigated. Three different reactor types (PWR, CANDU, MSR) were modeled with mixed thorium-uranium oxide fuels that were burned until the fuel was spent. The protactinium in the used fuel was extracted at the time of shutdown and the change in isotopic content of the protactinium quantified. Gamma spectroscopy simulations were performed for the protactinium isotopes and their decay products at various decay times. Given the simplicity of the models and large assumptions made (e.g. no background, no shielding, no self-attenuation), the initial results indicate that though 233 Pa is detectible for all the reactor types modeled at all decay times (0 to 300 days), more work should be done with higher fidelity models.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2023

Entanglement masquerading in the CMB

The simplest single-field inflation models capture all the relevant contributions to the patterns in the Cosmic Microwave Background (CMB) observed today. A key assumption in these models is that the quantum inflationary fluctuations that source such patterns are generated by a particular quantum state — the Bunch-Davies (BD) state. While this is a well-motivated choice from a theoretical perspective, the question arises of whether current data can rule out other, also well motivated, choices of states. In particular, as we previously demonstrated in [1], entanglement is naturally and inevitably dynamically generated during inflation given the presence of a "rolling" spectator scalar field — and the resulting entangled state will yield a primordial power spectrum with potentially measurable deviations compared to the canonical BD result. For this work we developed a perturbative framework to allow a systematic exploration of constraints on (or detection of) entangled states with Planck CMB data using Monte Carlo techniques. We have found that most entangled states accessible with our framework are consistent with the data. One would have to expand the framework to allow a greater variety of entangled states in order to saturate the Planck constraints and more systematically explore any preferences the data may have among the different possibilities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Multimessenger cosmology: Correlating cosmic microwave background and stochastic gravitational wave background measurements

Characterizing the physical properties of the stochastic gravitational wave background (SGWB) is a key step towards identifying the nature of its possible origin. We focus our analysis on SGWB anisotropies. The existence of a nontrivial primordial scalar-tensor-tensor (STT) correlation in the squeezed configuration may be inferred from the effect that a long wavelength scalar mode has on the gravitational wave power spectrum: an anisotropic contribution. Crucially, such a contribution is correlated with temperature anisotropies in the cosmic microwave background (CMB). This paper shows that, for inflationary models that generate suitably large STT non-Gaussianities, cross correlating the CMB with the stochastic background of gravitational waves is a very effective probe of early universe physics. The resulting signal can be a smoking gun for primordial SGWB anisotropies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Improved Characterization of the Astrophysical Muon–neutrino Flux with 9.5 Years of IceCube Data

We present a measurement of the high-energy astrophysical muon–neutrino flux with the IceCube Neutrino Observatory. The measurement uses a high-purity selection of 650k neutrino-induced muon tracks from the northern celestial hemisphere, corresponding to 9.5 yr of experimental data. With respect to previous publications, the measurement is improved by the increased size of the event sample and the extended model testing beyond simple power-law hypotheses. An updated treatment of systematic uncertainties and atmospheric background fluxes has been implemented based on recent models. The best-fit single power-law parameterization for the astrophysical energy spectrum results in a normalization of ${\phi }_{@100\mathrm{TeV}}^{{\nu }_{\mu }+{\bar{\nu }}_{\mu }}={1.44}_{-0.26}^{+0.25}\times {10}^{-18}\,{\mathrm{GeV}}^{-1}{\mathrm{cm}}^{-2}{{\rm{s}}}^{-1}{\mathrm{sr}}^{-1}$ and a spectral index ${\gamma }_{\mathrm{SPL}}={2.37}_{-0.09}^{+0.09}$, constrained in the energy range from 15 TeV to 5 PeV. The model tests include a single power law with a spectral cutoff at high energies, a log-parabola model, several source-class-specific flux predictions from the literature, and a model-independent spectral unfolding. The data are consistent with a single power-law hypothesis, however, spectra with softening above one PeV are statistically more favorable at a two-sigma level.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Multiphysics Degradation Modeling of Energy Storage Materials via RKPM with a Neural Network-Enhancement

In energy storage materials, strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking during charge/discharge cycling, resulting in reduced performance and service life. A coupled electro-chemo-mechanical reproducing kernel particle method (RKPM) formulation is developed, and a patch-test is formulated to certify optimal convergence of the proposed RKPM method for the coupled physics system. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), pixel-based model construction by RKPM is then used to represent the complex material microstructures for modeling the coupled physics of these systems. Further, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1, 2] is introduced to effectively model damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022. [2] Baek, J., Chen, J. S., "A Neural Network-Based Enrichment of Reproducing Kernel Approximation for Modeling Brittle Fracture", Computer Methods in Applied Mechanics and Engineering Vol. 410, 116590, 2024.

electro-chemo-mechanical coupling↗