Search NASA⌕ Search

SEARCH · Search NASA

Results for “Modeling & Simulation (M&S)”

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

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
DOE OSTI2025

Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape: Modeling Archive

This modeling archive is in support of the Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) publication "Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape", by Murphy et al. (2025). This archive contains model input files and outputs from landscape-scale simulations conducted using ELM, the land model component of the Department of Energy’s Energy Exascale Earth System Model (E3SM), at the Council NGEE Arctic field site (Council Road mile marker 71) on Alaska’s Seward Peninsula. Input data and model output from two sets of ELM simulations are provided. The first set of simulations were conducted with the two default ELM Arctic plant functional types (PFTs; broadleaf deciduous boreal shrub and a C3 grass) and the second set of simulations were conducted with a set of nine Arctic-specific PFTs including nonvascular mosses and lichens, graminoids, forbs, evergreen dwarf shrubs, three height classes of deciduous shrubs (dwarf, low, and low to tall), and deciduous alder shrubs (Sulman et al., 2021). Parameter names and major parameter changes in the Arctic-specific PFT configuration are described in Sulman et al. (2021) and archived in the Sulman et al. (2021) dataset (see below). Simulations were spatially explicit, covering an approximately 6.4X3.3 km domain at the Council site with a spatial resolution of 100 m for a total of 2,112 simulated grid cells under each ELM PFT configuration. The modeling archive contains meteorological forcing (seven *.nc files and one *.txt file), a domain definition file (one *.nc files), land surface configuration files (two *.nc files), parameter files (two *.nc files), annual ELM output files spanning 1980-2014 (68 *.nc files), and a User’s Guide (*pdf file). Additional information on the provided files is in the “Modeling Archive Contents” section of the User’s Guide. Model outputs are aggregated to the column scale (i.e. PFT-specific outputs are not provided here).

Murphy, Bailey [ORNL] (ORCID:0000000203995221)↗
DOE OSTI2023

Nested solitons in two-field fuzzy dark matter

ABSTRACT Dark matter as scalar particles consisting of multiple species is well motivated in string theory where axion fields are ubiquitous. A two-field fuzzy dark matter (FDM) model features two species of ultralight axion particles with different masses, m1 ≠ m2, which is extended from the standard one-field model with $m_a \sim 10^{-22} \, {\rm eV}$. Here we perform numerical simulations to explore the properties of two-field FDM haloes. We find that the central soliton has a nested structure when m2 ≫ m1, which is distinguishable from the generic flat-core soliton in one-field haloes. However, the formation of this nested soliton is subject to many factors, including the density fraction and mass ratio of the two fields. Finally, we study non-linear structure formation in two-field cosmological simulations with self-consistent initial conditions and find that the small-scale structure in two-field cosmology is also distinct from the one-field model in terms of DM halo counts and soliton formation time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Dark Energy Survey year 6 results: Magnification modeling and its impact on galaxy clustering and galaxy-galaxy lensing cosmology

Gravitational lensing magnification alters the observed spatial distribution of galaxies and must be accounted for to prevent biases in cosmological probes of the large-scale structure. We investigate its effects on the Dark Energy Survey Year 6 galaxy clustering and galaxy-galaxy lensing analyses using the fiducial lens (position tracer) sample M ag L im++. Magnification bias is parameterized by a coefficient that describes the response of the number of selected objects per unlensed area element to a change in the lensing convergence. We quantify this coefficient using the BALROG synthetic source injection catalog to account for the complexity of the selection function, and compare these results with simplified estimates. The resulting values of the magnification coefficients for each redshift bin are [3.16 ± 0.08, 2.76 ± 0.21, 4.09 ± 0.15, 4.42 ± 0.16, 4.90 ± 0.29, 4.83 ± 0.25]. Relative to Year 3, this analysis provides more precise and accurate magnification bias estimates through a larger BALROG area and reweighting to better match the data properties. Here, the cosmological results are robust when tested against various magnification parameter prior choices and also when adding cross-clustering between lens redshift bins. Neglecting magnification, however, introduces significant systematic shifts: relative to the fiducial analysis with Gaussian priors centered on the BALROG -derived estimates, we observe shifts of 1.37σ in S 8 and -0.84σ in Ω m (with cosmic shear included: -0.61σ in S 8 and -0.71σ in Ω m ), in agreement with findings from simulated data, demonstrating that magnification must be modeled to avoid biases. Freeing the magnification bias in lens bin 2 leads to unphysical negative values, further justifying its exclusion from the fiducial Year 6 analysis.

Cosmological parameters↗
DOE OSTI2025

Glacial geology of the Hudson Mountains, Amundsen Sea sector, West Antarctica

The Hudson Mountains are situated in the eastern Amundsen Sea sector of the West Antarctic Ice Sheet, adjacent to Pine Island Glacier. They form a volcanic field of 17 stratovolcanoes and parasitic vents, preserved as nunataks. Two former tributaries of Pine Island Glacier (Larter and Lucchitta glaciers) flow through the mountains. Here we present a detailed study of the glacial geology of the area. We describe field observations and measurements of geomorphological features from 15 of the nunataks, meltwater ponds found on the surface of three nunataks and supraglacial features (ice dolines) from two sites near the present grounding line. Together these provide constraints on the past ice sheet extent, flow pathways and thermal regime, and enhance our understanding of the present hydrological regime – all of which are important as context for the observed modern ice sheet behaviour. We find evidence suggesting that all nunataks in the Hudson Mountains were covered by ice during the Last Glacial Maximum (defined here as 26.5-19 ka) and have since deglaciated. Faceted and polished erratic cobbles and boulders of exotic lithologies (syenites, alkali granites, granites, granodiorites, tonalites and gabbros) are numerous and perched on nunatak surfaces. A marked difference between the dominant erratic lithologies on nunataks adjacent to Pine Island Glacier (granite) and Lucchitta Glacier (granodiorite-tonalite) indicates that the ice sheet was transporting clasts from at least two distinct upstream source regions. The similarity in degree of weathering suggests, however, that all the erratics were transported by one phase of (warm-based) glaciation; their presence on or close to the summits of all except one nunatak indicates that the ice sheet during that time was at least 700 m thicker than present. These results are consistent with ice sheet model simulations which suggest that all nunataks in the Hudson Mountains were completely submerged by the Last Glacial Maximum ice sheet.

58 GEOSCIENCES↗
DOE OSTI2021

Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star

A realistic parametrization of convection and convective boundary mixing in conventional stellar evolution codes is still the subject of ongoing research. Furthermore, to improve the current situation, multidimensional hydrodynamic simulations are used to study convection in stellar interiors. Such simulations are numerically challenging, especially for flows at low Mach numbers which are typical for convection during early evolutionary stages. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. Simulations of convection for a realistic stellar profile are analyzed regarding the properties of convective boundary mixing. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25 M ⊙ star model. The results obtained with the low-Mach AUSM + -up solver were compared to results when using its non low-Mach variant AUSM B + -up. We applied well-balancing of the gravitational source term to maintain the initial hydrostatic background stratification. The computational grids have resolutions ranging from 180 × 90 2 to 810 × 540 2 cells and the nuclear energy release was boosted by factors of 3 × 10 3 , 1 × 10 4 , and 3 × 10 4 to study the dependence of the results on these parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Wake Statistics of Different-Scale Wind Turbines under Turbulent Boundary Layer Inflow

Subscale wind turbines can be installed in the field for the development of wind technologies, for which the blade aerodynamics can be designed in a way similar to that of a full-scale wind turbine. However, it is not clear whether the wake of a subscale turbine, which is located closer to the ground and faces different incoming turbulence, is also similar to that of a full-scale wind turbine. In this work we investigate the wakes from a full-scale wind turbine of rotor diameter 80 m and a subscale wind turbine of rotor diameter of 27 m using large-eddy simulation with the turbine blades and nacelle modeled using actuator surface models. The blade aerodynamics of the two turbines are the same. In the simulations, the two turbines also face the same turbulent boundary inflows. The computed results show differences between the two turbines for both velocity deficits and turbine-added turbulence kinetic energy. Such differences are further analyzed by examining the mean kinetic energy equation.

17 WIND ENERGY↗
DOE OSTI2021

Thermal-Hydrological-Mechanical Modelling of Stockton University Reservoir Cooling System, Fine Scale Stress Test Modelling

Mesh, properties, initial conditions, injection/withdrawal rates for modelling thermal, hydrological, and mechanical effects of fluid injection to and withdrawal from ground for Stockton University reservoir cooling system (aquifer storage cooling system), Galloway, New Jersey, for unscheduled two hour injection at 133 % designed capacity, on fine scale grid, with some results. Second simulation of J.T. Smith, E. Sonnenthal, P. Dobson, P. Nico, and M. Worthington, 2021. Thermal-hydrological-mechanical modeling of Stockton University reservoir cooling system, Proceedings of the 46th Workshop on Geothermal Reservoir Engineering, Stanford University, SGP-TR-218, from which Figures 6-9, pertain.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Modulation of the Eastern Equatorial Pacific Seasonal Cycle by Tropical Instability Waves

Feedbacks from tropical instability waves (TIWs) on the seasonal cycle of the eastern Pacific Ocean are studied using two eddy-rich ocean simulations, with and without TIWs. By warming the equatorial waters by up to 0.4°C through nonlinear advection in boreal summer and fall, TIWs reduce the amplitude of the seasonal cycle in upper ocean temperatures. In addition, TIWs stabilize the upper part of the Equatorial Undercurrent (EUC) through enhanced barotropic energy conversion, leading to a year-round weakening by -0.15 m s -1 and preventing an unrealistic re-intensification in boreal fall usually found in non-eddy resolving models. A coarser simulation at 1-degree horizontal resolution fails to reproduce the TIW-induced nonlinear warming of equatorial waters, but succeeds in inhibiting the EUC re-intensification. This suggests a threshold effect in TIW strength, associated with the model's ability to simulate eddies, which may be responsible for long-standing biases displayed by global climate models in this region.

58 GEOSCIENCES↗
DOE OSTI2021

ABACUSSUMMIT : a massive set of high-accuracy, high-resolution N -body simulations

We present the public data release of the AbacusSummit cosmological N-body simulation suite, produced with the Abacus N-body code on the Summit supercomputer of the Oak Ridge Leadership Computing Facility. Abacus achieves $\mathcal {O}(10^{-5})$ median fractional force error at superlative speeds, calculating 70M particle updates per second per node at early times, and 45M particle updates per second per node at late times. The simulation suite totals roughly 60 trillion particles, the core of which is a set of 139 simulations with particle mass $2\times 10^{9}\, h^{-1}\, \mathrm{M}_\odot$ in box size $2\, h^{-1}\, \mathrm{Gpc}$. The suite spans 97 cosmological models, including Planck 2018, previous flagship simulation cosmologies, and a linear derivative and cosmic emulator grid. A subsuite of 1883 boxes of size $500\, h^{-1}\, \mathrm{Mpc}$ is available for covariance estimation. AbacusSummit data products span 33 epochs from z = 8 to 0.1 and include light cones, full particle snapshots, halo catalogues, and particle subsets sampled consistently across redshift. AbacusSummit is the largest high-accuracy cosmological N-body data set produced to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Black hole-neutron star binaries near neutron star disruption limit in the mass regime of event GW230529

In May 2023, the LIGO-Virgo-KAGRA (LVK) Collaboration reported the likely black hole-neutron star (BHNS) merger GW230529_181500. The signal was observed with high significance in only one detector, limiting constraints on the black hole spin and motivating our study of disruption in this mass regime. That event is expected to be the merger of a 2.5–4.5 M ⊙ primary with a secondary compact object of mass between 1.2–2.0 M ⊙ . This makes it the first BHNS merger with a significant potential for the production of electromagnetic (EM) counterparts, and provides further evidence for compact objects existing within the suspected lower mass gap. To produce post-merger EM transients, the component of the black hole spin aligned with the orbital angular momentum must be sufficiently high, allowing the neutron star to be tidally disrupted. The disrupting BHNS binary may then eject a few percent of a solar mass of matter, leading to an observable kilonova driven by radioactive decays in ejecta, and/or a compact-binary gamma-ray burst (cbGRB) resulting from the formation of an accretion disk and relativistic jet. Determining which mergers lead to disruption of the neutron star is necessary to predict the prevalence of EM signals from BHNS mergers, yet most BHNS simulations so far have been performed far from the minimum spin required for tidal disruption. Here, we use the Spectral Einstein Code to explore the behavior of BHNS mergers in a mass range consistent with GW230529_181500 close to that critical spin, and compare our results against the mass remnant model currently used by the LVK Collaboration to predict the probability of tidal disruption. Furthermore, our numerical results reveal the emergence of non-zero accretion disks even below the predicted NS disruption limit, of low mass but capable of powering cbGRBs. Our results also demonstrate that the remnant mass model underpredicts the disk mass for the DD2 equation of state, while they are within expected modeling errors for SFHo. The disruption limit itself, however, is not found to significantly differ from the predictions of the analytical model, unless remnant masses M rem ≲ 0.001 M ⊙ prove interesting observationally. In all of our simulations, any kilonova signal would be dim and most likely dominated by post-merger disk outflows.

GW230529_181500↗
DOE OSTI2023

Measurements of dielectron production in Au + Au collisions at s N N = 27 , 39, and 62.4 GeV from the STAR experiment

We report systematic measurements of dielectron (e + e – ) invariant-mass M ee spectra at midrapidity in Au + Au collisions at s N N = 27 = 27, 39, and 62.4 GeV taken with the STAR detector at the Relativistic Heavy Ion Collider. For all energies studied, a significant excess yield of dielectrons is observed in the low-mass region ( 0.40 < M ee < 0.75 MeV / c 2 ) compared to hadronic cocktail simulations at freeze-out. Models that include an in-medium broadening of the ρ -meson spectral function consistently describe the observed excess. In addition, we report acceptance-corrected dielectron-excess spectra for Au + Au collisions at midrapidity ( |y ee | < 1) in the 0–80% centrality bin for each collision energy. The integrated excess yields for 0.4 < M ee < 0.75 GeV / c 2 , normalized by the charged particle multiplicity at midrapidity, are compared with previously published measurements for Au + Au at s N N = 27 = 19.6 and 200 GeV. Models that include an in-medium broadening of the ρ -meson spectral function consistently describe the observed excess. Furthermore, the normalized excess yields in the low-mass region show no significant collision energy dependence. The data, however, are consistent with model calculations that demonstrate a modest energy dependence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

X-ray morphology of cluster-mass haloes in self-interacting dark matter

ABSTRACT We perform cosmological zoom-in simulations of 19 relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with $M_{\rm 200} \simeq (1\!-\!3)\times 10^{15}{\, \rm M_\odot }$. Both the dark matter and the intracluster gas distributions in SIDM appear more spherical than their CDM counterparts. Mock X-ray images are generated based on the simulations and are compared to the real X-ray images of 84 relaxed clusters selected from the Chandra and ROSAT archives. We perform ellipse fitting for the isophotes of mock and real X-ray images and obtain the ellipticities at cluster-centric radii of $r\simeq 0.1\!-\!0.2R_{\rm 200}$. The X-ray isophotes in SIDM models with increasing cross-sections are rounder than their CDM counterparts, which manifests as a systematic shift in the distribution function of ellipticities. Unexpectedly, the X-ray morphology of the observed non-cool-core clusters agrees better with SIDM models with cross-section $(\sigma /m)= 0.5\!-\!1\, {\rm cm}^2\, {\rm g}^{-1}$ than CDM and SIDM with $(\sigma /m)=0.1\, {\rm cm}^2\, {\rm g}^{-1}$. Our statistical analysis indicates that the latter two models are disfavoured at the $68{{\ \rm per\ cent}}$ confidence level (as conservative estimates). This conclusion is not altered by shifting the radial range of measurements or applying a temperature selection criterion. However, the primary uncertainty originates from the lack of baryonic physics in the adiabatic model, such as cooling, star formation and feedback effects, which still have the potential to reconcile CDM simulations with observations.

Shen, Xuejian (ORCID:000000026196823X)↗
DOE OSTI2020

Early-Universe Simulations of the Cosmological Axion

Ultracompact dark matter (DM) minihalos at masses at and below 10 –12 $M_⊙$ arise in axion DM models where the Peccei-Quinn (PQ) symmetry is broken after inflation. The minihalos arise from density perturbations that are generated from the nontrivial axion self-interactions during and shortly after the collapse of the axion-string and domain-wall network. We perform high-resolution simulations of this scenario starting at the epoch before the PQ phase transition and ending at matter-radiation equality. We characterize the spectrum of primordial perturbations that are generated and comment on implications for efforts to detect axion DM. In this work, we also measure the DM density at different simulated masses and argue that the correct DM density is obtained for $m_a$ = 25.2 ± 11.0 μ eV .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Field-scale dynamics of planting dates in the US Corn Belt from 2000 to 2020

Crop planting dates are a dynamic feature of agricultural systems that respond to short- and long-term climate signals, crop and cultivar selection, and technology changes. Planting date records are essential for yield gap analyses, accurate crop modeling, and tracking farmer adaptations to weather and climate change. Although planting dates have high variation at local scales due to heterogeneity in farm resources and decision-making, available long-term data on planting dates is largely restricted to aggregated regional statistics or, at best, satellite-derived datasets with limited spatiotemporal extent and at resolutions unable to distinguish individual fields (> 250 m). Here, we generated retrospective annual field-scale (30 m) planting date maps for both maize and soybeans spanning 2000-2020 across a 12 state region in the United States Corn Belt based on Landsat satellite data and a large ground sample of over 28,000 maize and soybean fields. Using training data from 2015-2020 for model selection, we found that planting date predictions improved with harmonic regression of Landsat data and additional annual weather covariates. The preferred random forests model approximately doubled performance compared to a null model based on state median planting dates, capturing 47% of field-level variation for maize (mean absolute error, MAE = 7.4 days) and 44% for soybeans (MAE = 7.5 days) against held-out ground truth test data for 2008-2014. We also evaluated the full 2000-2020 dataset with state agricultural statistics, finding strong agreement with median planting dates for maize (R 2 = 0.76, MAE = 4.4 days) and slightly lower agreement for soybeans (R 2 = 0.65, MAE = 5.4 days) when aggregated to the state level. We then used this new dataset to analyze environmental determinants of planting dates at a finer-scale than previously possible, controlling for unobserved variation at the sub-state district level. We found that during 2000-2020, each standard deviation increase in rainfall delayed planting by ~ 2.5 days, and fields with higher soil productivity ratings tended to be planted earlier. We did not find meaningful trends over the last two decades in planting dates for maize or soybeans, in contrast to trends towards earlier planting dates late last century and predicted for this period in climate adaptation studies. We hypothesize increases in early season rainfall may have inhibited these shifts towards earlier planting. Remotely sensed planting dates will be a useful tool for yield gap analyses, crop simulation modeling, and ongoing assessment of climate adaptation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Core collapse supernova gravitational wave emission for progenitors of 9.6, 15, and 25 M ⊙

Here, we present gravitational wave emission predictions based on three core collapse supernova simulations corresponding to three different progenitor masses. The masses span a large range, between 9.6 and 25M ⊙ , are all initially nonrotating, and are of two metallicities: zero and solar. We compute both the temporal evolution of the gravitational wave strains for both the plus and the cross polarizations, as well as their spectral decomposition and characteristic strains. The temporal evolution of our zero metallicity 9.6M ⊙ progenitor model is distinct from the temporal evolution of our solar metallicity 15M ⊙ progenitor model and our zero metallicity 25M ⊙ progenitor model. In the former case, the high-frequency gravitational wave emission is largely confined to a brief time period ~75 ms after bounce, whereas in the latter two cases high-frequency emission does not commence until ~125 ms after bounce or later. The excitation mechanisms of the high-frequency emission in all three cases correspond to proto-neutron star convection and accretion onto the proto-neutron star from the convective gain layer above it, with the former playing the dominant role for most of the evolution. The low-frequency emission in all three models exhibits very similar behavior. At frequencies below ~250 Hz, gravitational waves are emitted by neutrino-driven convection and the standing accretion shock instability (SASI). This emission extends throughout the simulations when a gain region is present. In all three models, explosion is observed at ~125, ~500, and ~250 ms after bounce in the 9.6, 15, and 25M ⊙ progenitor models, respectively. At these times, the low-frequency gravitational wave emission is joined by very low-frequency emission, below ~10 Hz. These very low-frequency episodes are the result of explosion and begin at the above designated explosion times in each of our models. Our characteristic strains tell us that, in principle, all three gravitational wave signals would be detectable by current-generation detectors for a supernova at a distance of 10 kpc. However, our 9.6M ⊙ progenitor model is a significantly weaker source of gravitational waves, with strain amplitudes approximately 5–10 times less than in our other two models. The characteristic strain for this model tells us that such a supernova would be detectable only within a much more narrow frequency range around the maximum sensitivity of today’s detectors. Finally, in our 9.6M ⊙ progenitor model, we see very high-frequency gravitational radiation, extending up to ~2000 Hz. This feature results from the interaction of shock- and deleptonization-induced convection with perturbations introduced in the progenitor by nuclear burning during core collapse. While unique to the 9.6M ⊙ progenitor model analyzed here, this very high-frequency emission may, in fact, be a generic feature of the predictions for the gravitational wave emission from all core collapse supernova models when simulations are performed with three-dimensional progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗