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

Impact of Vibrational Nonequilibrium on the Supersonic Mode in Hypersonic Boundary Layers

The supersonic mode in hypersonic boundary layers has recently been shown to achieve a larger peak growth rate than the traditional second mode in particular circumstances relevant to free flight. In this environment, high levels of aerodynamic heating are encountered such that incorporation of chemical and vibrational nonequilibrium effects can become critical. The impact of thermal nonequilibrium (TNE) on the supersonic mode in particular has not been thoroughly and systematically investigated. This work uses thermochemical nonequilibrium (TCNE) direct numerical simulation (DNS) and TCNE linear stability theory (LST) to obtain a more complete investigation of the supersonic mode using both nonequilibrium and frozen gas models. The simulation is Mach 5 flow over a 1-mm-nose-radius axisymmetric cone 1 m in length, representative of experimental conditions in ground test facilities. Both LST and DNS results indicate that TNE effects are destabilizing compared with frozen flow and can lead to an N-factor difference of approximately +1 at the end of the cone. LST results also indicate that TNE effects are stabilizing to the second and supersonic modes in comparison with thermal equilibrium. The supersonic mode appeared to be slightly more affected by TNE than the second mode. In conclusion, collectively, the findings in this work suggest that even in ground test facilities where TCNE effects are relatively small compared with free flight, incorporating nonequilibrium models in both the mean flow and stability solvers can have a practical effect on hypersonic boundary-layer transition prediction analyses.

42 ENGINEERING↗
DOE OSTI2022

Stochastic agent-based model for predicting turbine-scale raptor movements during updraft-subsidized directional flights

Rapid expansion of wind energy development across the world has highlighted the need to better understand turbine-caused avian mortality. The risk to golden eagles (Aquila chrysaetos) is of particular concern due to their small population size and conservation status. Golden eagles subsidize their flight in part by soaring in orographic updrafts, which can place them in conflict with wind turbines utilizing the same low-altitude wind resource. Understanding the behavior of soaring raptors in varying atmospheric conditions can therefore be relevant to predicting and mitigating their risk of collision. We present a predictive movement model that simulates individual paths of golden eagles during directional flight (such as migration) that is subsidized by orographic updraft. We modeled eagles in a 50 km by 50 km study area in Wyoming containing three wind power plants with documented golden eagle collisions with turbines. The movement model is applicable to any region where ground elevation is known at turbine scale (50 m) and wind conditions are known at facility scale (3 km). For a given set of atmospheric conditions, the model simulates movements of thousands of orographic soaring eagles to produce a density map quantifying the relative probability of eagle presence. We validated the simulated tracks with GPS telemetry data showing four directional tracks made by golden eagles transiting through the area in 2019 and 2020. For each eagle track, validation was performed using the ratio of the model-simulated eagle presence likelihood with uniform eagle presence and the presence computed using directed random-walk movements. We found that the predictive performance of the model was significantly better (likelihood ratio 1) for low-altitude movements than high-altitude movements that can involve thermal-soaring. We employed the model to produce seasonal presence maps for migrating golden eagles. We found significant turbine-level variations in eagle presence between northerly and southerly migration routes through the study area. Overall, the proposed model offers a generalizable, probabilistic, and predictive tool to assist wind energy developers, ecologists, wildlife managers, and industry consultants in estimating the potential for conflict between soaring birds and wind turbines, thereby reducing the need for site-specific data on golden eagle movements.

17 WIND ENERGY↗
DOE OSTI2021

Three-dimensional Hydrodynamic Simulations of Convective Nuclear Burning in Massive Stars Near Iron Core Collapse

Nonspherical structure in massive stars at the point of iron core collapse can have a qualitative impact on the properties of the ensuing core-collapse supernova explosions and the multimessenger signals they produce. Strong perturbations can aid successful explosions by strengthening turbulence in the postshock region. Here we report on a set of 4π 3D hydrodynamic simulations of O- and Si-shell burning in massive star models of varied initial masses using MESA and the FLASH simulation framework. We evolve four separate 3D models for roughly the final 10 minutes prior to and including iron core collapse. We consider initial 1D MESA models with masses of 14, 20, and 25 M ⊙ to survey a range of O/Si-shell density and compositional configurations. Further, we characterize the convective shells in our 3D models and compare them to the corresponding 1D models. In general, we find that the angle-average convective speeds in our 3D simulations near collapse are three to four times larger than the convective speeds predicted by MESA at the same epoch for our chosen mixing length parameter of α MLT = 1.5. In three of our simulations, we observe significant power in the spherical harmonic decomposition of the radial velocity field at harmonic indices of ℓ = 1–3 near collapse. Our results suggest that large-scale modes are common in massive stars near collapse and should be considered a key aspect of presupernova progenitor models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

High-resolution (30-m) urban land cover projections for Los Angeles California Urban Area: 2010 to 2100 under SSP5

These data represent simulations of future land use and land cover (LULCC) for Los Angeles urban area (U.S. Census Bureau defined area) as raster tiff images at a 30-m pixel resolution and at decadal time steps from 2010 to 2100. LULCC classes in this product follow the National Land Cover Dataset (NLCD) classification. NLCD 21-24 correspond to open developed, low developed, medium developed, and high developed urban land classes, respectively. Only urban land cover classes (NLCD class 21, 22, 23, and 24) are dynamic over time; however, all NLCD classes are included in the final product. Therefore, NLCD classes that do not convert to an urban class will be similar to year 2000. The products were developed using a hybridized statistical and cellular automata approach. Linear mixed models (LMMs) were used to estimate future urban land budgets based on 1-km urban land fraction projections from Gao and Pesaresi (2021), whereas separate generalized linear mixed models (GLMMs) were used to estimate shifts in urban land intensities based on retrospective shifts in NLCD urban class intensities over a 20- year period. Based on urban land allocations from the statistical models, a cellular-automata and downscaling routine was used to simulate dynamic urban land expansion at a 30-m resolution based on suitability criteria. Scenarios of future urban landcover change projections include variant solutions for the Shared Socioeconomic Pathway 5 (SSP5) based on different population assumptions, different land use intensification assumptions, variable land zoning constraints, and iterative adjustments to correct for over allocation of urban expansion across decadal time periods from 2010 to 2100. This results in 320 raster products.

Land↗
DOE OSTI2023

High-resolution (30-m) urban land cover projections for Los Angeles California Urban Area: 2010 to 2100 under SSP3 and SSP5 [Updated simulations based on population-driven urban intensity transitions]

These data (v3) are updated from previous versions (1 and 2) in that they include consider the effects of population on transitions in urban land intensity. This leads to more reasonable differences in urban land projections under variant SSPs. For the present dataset, both SSP3 and SSP5 are provided. These data represent simulations of future land use and land cover (LULCC) for Los Angeles urban area (U.S. Census Bureau defined area) as raster tiff images at a 30-m pixel resolution and at decadal time steps from 2010 to 2100. LULCC classes in this product follow the National Land Cover Dataset (NLCD) classification. NLCD 21-24 correspond to open developed, low developed, medium developed, and high developed urban land classes, respectively. Only urban land cover classes (NLCD class 21, 22, 23, and 24) are dynamic over time; however, all NLCD classes are included in the final product. Therefore, NLCD classes that do not convert to an urban class will be similar to year 2000. The products were developed using a hybridized statistical and cellular automata approach. Linear mixed models (LMMs) were used to estimate future urban land budgets based on 1-km urban land fraction projections from Gao and Pesaresi (2021), whereas separate generalized linear mixed models (GLMMs) were used to estimate shifts in urban land intensities based on retrospective shifts in NLCD urban class intensities over a 20- year period. Based on urban land allocations from the statistical models, a cellular-automata and downscaling routine was used to simulate dynamic urban land expansion at a 30-m resolution based on suitability criteria. Scenarios of future urban landcover change projections include variant solutions for the Shared Socioeconomic Pathway 5 (SSP5) and SSP 3 based on different population assumptions, different land use intensification assumptions, variable land zoning constraints, and iterative adjustments to correct for over allocation of urban expansion across decadal time periods from 2010 to 2100. This results in 320 raster products.

Land↗
DOE OSTI2024

Implementation and Validation of a Generalized Actuator Disk Parameterization for Wind Turbine Simulations Within the FastEddy Model

ABSTRACT Fast and accurate large‐eddy simulation (LES) of the atmospheric boundary layer plays a crucial role in advancing wind energy research. Long‐duration wind farm studies at turbine‐resolving scales have become increasingly important to understand the intricate interactions between large wind farms and the atmospheric boundary layer. However, the prohibitive computational cost of these turbulence‐ and turbine‐resolving simulations has precluded such modeling to be exercised on a regular basis. To that end, we implement and validate the generalized actuator disk (GAD) model in the computationally efficient, graphics processing unit (GPU)–resident, LES model FastEddy. We perform single‐turbine simulations under three atmospheric stabilities (neutral, unstable, and stable) and compare them against observations from the Scaled Wind Farm Technology (SWiFT) facility and other LES codes from the recent Wakebench turbine wake model benchmark. Our idealized LES results agree well with observed wake velocity deficit and downstream recovery across stability regimes. Turbine response in terms of rotational speed, generated power, torque, and thrust coefficient are well predicted across stability regimes and are consistent with the LES results from the benchmark. The FastEddy simulations are found to be at least two orders of magnitude more efficient than the traditional CPU‐based LES models, opening the door for realistic LES simulations of full wind plants as a viable standard practice.

17 WIND ENERGY↗
DOE OSTI2019

Microstructure evolution during near- T g annealing and its effect on shear banding in model alloys

By performing extensive molecular dynamics simulations, we investigate herein the deformation behavior in A l 90 S m 10 and C u 64.5 Z r 35.5 alloys after elongated isothermal annealing in the vicinity of the glass-transition temperature ( T g ). Different microstructural response to the annealing process was observed: A l 90 S m 10 maintains the glassy structure with improved energetic stability, enhanced short-range order (SRO), and a more pronounced spatial network that extends beyond the first atomic shell, while C u 64.5 Z r 35.5 forms nanocrystalline Laves C u 2 Zr phases. Shear banding occurs in both annealed systems under shear loading. For A l 90 S m 10 , the spatial network formed by the local clusters characterizing the SRO of the system is significantly weakened but remains appreciable in the shear band. In contrast, the crystalline ordering in the C u 64.5 Z r 35.5 is completely destroyed during shear banding. Consequently, while displaying higher yield strength, the annealed C u 64.5 Z r 35.5 sample appears to be less ductile. By carefully examining the effect of microstructures on the structural ordering in the shear band and the consequent mechanical response, our work contributes to a better understanding of the deformation mechanism of amorphous alloys as compared with that in crystalline materials.

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

Development of a high charge 10 GeV laser electron accelerator

Low-density meter-scale plasma waveguides produced in meter-scale supersonic gas jets have paved the way for recent demonstrations of all-optical multi-gigaelectronvolt laser wakefield acceleration (LWFA). This paper reviews recent advances by the University of Maryland, which have enabled these results, focusing on the development of elongated supersonic gas jets up to ∼1 m in length, experimental and simulation studies of plasma waveguide formation, and a new three-stage model for relativistic pulse propagation dynamics in these waveguides. We also present results from recent LWFA experiments conducted at the Laboratory for Advanced Lasers and Extreme Photonics at Colorado State University demonstrating high charge, low divergence electron bunches to ∼10 GeV, with laser-to-electron beam efficiency of at least ∼30%.

Rockafellow, E. [Univ. of Maryland, College Park, ↗
DOE OSTI2024

Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins

This repository of global hydrogeologic datasets contains aquifer properties on 0.5° scale, including depth to groundwater (Fan et al., 2013), aquifer thickness (de Graaf et al., 2015), WHYMap aquifer classes (Richts et al., 2011), recharge (Döll and Fiedler, 2008; Gleeson et al., 2016), lakes (Messager et al., 2016), porosity and permeability (Gleeson et al., 2014), digitized and geo-processed from their respective sources. Globally gridded aquifer properties could be used independently to estimate global groundwater availability or used as critical inputs to the superwell model to simulate groundwater extraction and provide estimates of pumped volumes and unit costs under user-specific scenarios. Key resources related to this data are: Niazi, H., Ferencz, S. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 superwell model repository which uses this data to simulate groundwater extraction and provides estimates of the global extractable volumes and unit-costs ($/km3) of accessible groundwater production under user-specified extraction scenarios. Repository Overview Main output: aquifer_properties_rec.csv contains all processed outputs, including aquifer properties like porosity, permeability, recharge, lake areas, aquifer thickness, and depth to groundwater. shapefiles.zip: contains all digitized GIS databases and shapefile for all aquifer properties prep_inputs.R and prep_inputs_recharge_lakes.R: R scripts that process the shapefiles to produce the aquifer_properties_rec.csv file plot_inputs.R: R script for plotting the maps and conducting preliminary analysis on the available groundwater volume basin_to_country_mapping.csv, basin_country_region_mapping.csv and continent_county_mapping.csv provide the mapping between continents, 32 energy-economic macro regions, countries, and water basins for post-processing aquifer_properties_rec.csv Maps: Each map visualizes the spatial distribution of one of the aquifer properties across the globe map_in_Porosity.png map_in_Permeability.png map_in_Aquifer_thickness.png map_in_Depth_to_water.png map_in_Recharge.png map_in_Grid_area_km.png map_in_Lake_area_km.png map_in_WHYClass.png Sample inputs sample_inputs.py: this script samples inputs from the aquifer_properties_rec dataset, ensuring the sampled and original inputs maintain the same distributions sampled_data_100.csv contains 100 sampled data points and sampled_data_100.png compares their distributions Dataset Overview The main outputs are consolidated in a comprehensive aquifer_properties_rec.csv file and include the following fields: GridCellID: Unique identifier for each (roughly 0.5°) grid cell Continent: Continent name Country: Country name GCAM_basin_ID: Identifier for GCAM hydrologic basin Basin_long_name: Full name of the basin WHYClass: Hydrogeologic classification based on WHYMap aquifer classes (Richts et al., 2011) Porosity: Soil porosity (%) (Gleeson et al., 2014) Permeability: Soil permeability (in square meters; Gleeson et al., 2014) Aquifer_thickness: Thickness of the aquifer (in meters; de Graaf et al., 2015) Depth_to_water: Depth to groundwater (in meters; Fan et al., 2013) Recharge: long-term annual averaged recharge rates (in m/yr; Döll and Fiedler, 2008; Gleeson et al., 2016) Grid_area: Area of the grid cell (in square meters) Lakes_area: Area of inland lakes (in square meters; Messager et al., 2016) Key References The datasets are digitized versions of global hydrogeologic properties from the following key literature sources: Depth to Groundwater: Fan, Y., Li, H., & Miguez-Macho, G. (2013). Global Patterns of Groundwater Table Depth. Science, 339(6122), 940-943. https://doi.org/10.1126/science.1229881 Aquifer Thickness: de Graaf, I. E. M., Sutanudjaja, E. H., van Beek, L. P. H., & Bierkens, M. F. P. (2015). A high-resolution global-scale groundwater model. Hydrol. Earth Syst. Sci., 19(2), 823-837. https://doi.org/10.5194/hess-19-823-2015 Porosity and Permeability: Gleeson, T., Moosdorf, N., Hartmann, J., & van Beek, L. P. H. (2014). A glimpse beneath earth's surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability and porosity. Geophysical Research Letters, 41(11), 3891-3898. https://doi.org/10.1002/2014GL059856 Aquifer classes: Richts, A., Struckmeier, W. F., & Zaepke, M. (2011). WHYMAP and the Groundwater Resources Map of the World 1:25,000,000. In J. A. A. Jones (Ed.), Sustaining Groundwater Resources: A Critical Element in the Global Water Crisis (pp. 159-173). Springer Netherlands. https://doi.org/10.1007/978-90-481-3426-7_10 Recharge: Döll, P., & Fiedler, K. (2008). Global-scale modeling of groundwater recharge. Hydrol. Earth Syst. Sci., 12(3), 863-885. https://doi.org/10.5194/hess-12-863-2008; Gleeson, T., Befus, K. M., Jasechko, S., Luijendijk, E., & Cardenas, M. B. (2016). The global volume and distribution of modern groundwater. Nature Geoscience, 9(2), 161-167. https://doi.org/10.1038/ngeo2590 Inland Lakes: Messager, M. L., Lehner, B., Grill, G., Nedeva, I., & Schmitt, O. (2016). Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nature Communications, 7(1), 13603. https://doi.org/10.1038/ncomms13603 Cite as Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MultiSector Dynamics-Living, Intuitive, Value-adding, Environment. https://doi.org/10.57931/2484226 Contact Reach out to Hassan Niazi or open an issue in the superwell repository for questions or suggestions.

aquifer thickness↗
DOE OSTI2026

Data for Impact of Vertical and Seasonal Variation in Leaf Traits on Simulating Soybean Canopy Photosynthesis via 1D and 3D Modeling

Accurate modeling of photosynthesis is crucial for predicting crop productivity and quantifying the carbon cycle in agroecosystems. Leaf traits are essential inputs for modeling canopy photosynthesis. Yet, many existing models still use fixed plant functional type (PTF)-based values to parameterize leaf traits under a big-leaf or two-big-leaf assumption, neglecting their vertical profiles and seasonal changes. This simplification may introduce significant uncertainties in estimating gross primary productivity (GPP). In this study, we simulated soybean GPP and tested the effects of vertical and seasonal variation in three key leaf photosynthetic traits: the maximum carboxylation rate at 25 °C (Vcmax25), leaf chlorophyll content (LCC), and leaf mass per area (LMA) in the 1D-SCOPE and 3D-Helios models. Weekly field measurements were conducted during the growing season of 2024 to support the simulation. We designed ten leaf trait parameterization schemes by incorporating different combinations of vertical profiles and seasonal changes, while assuming homogeneous canopy architecture in both models. Our results revealed that Vcmax25 vertical and seasonal variation had the strongest influence on simulated GPP in both 1D and 3D models, while LCC and LMA effects were minimal. Particularly, the scheme with an empirically parameterized Vcmax25 profile achieved comparable performance to the scheme with the measured Vcmax25 profile. Both 1D-SCOPE and 3D-Helios accurately modeled GPP (SCOPE: R2 = 0.87, Bias = 0.55 µmol m⁻² s⁻¹; Helios: R2 = 0.9, Bias = 0.22 µmol m⁻² s⁻¹) under the most complex scheme, and their responses to vertical and seasonal variation in leaf traits were consistent, demonstrating the robustness of our findings. Based on our findings, we propose a scalable framework for parameterizing leaf traits to improve GPP simulations. This study contributes to improving the representation of leaf trait dynamics in canopy-level photosynthesis models, potentially enhancing our ability to predict crop productivity and understand agroecosystem carbon dynamics.

Photosynthesis↗
DOE OSTI2025

Reduced Model of Ionization Lag in Intense Laser-Produced Plasmas

A physics-based empirical formula is derived to predict the ionization lag in underdense plasmas generated by an intense laser. Time-dependent nonlocal thermodynamic equilibrium calculations demonstrate significantly delayed ionization, due to rapid changes in plasma conditions, which critically impacts plasma properties such as opacity, emissivity, and heat transport. The reduced model, based on these calculations, enables the estimation of ionization lag without requiring in-depth knowledge of nonlocal thermodynamic equilibrium modeling. Furthermore, modeling reveals that the two-step ionization process—collisional excitation followed by photoionization—plays a crucial role in this ionization delay, with collisional excitation setting the timescale for ionization. Simulations across a range of elements, from beryllium to germanium, demonstrate that ionization lag is a widespread phenomenon, underscoring the importance of incorporating such processes into ionization modeling in radiation hydrodynamic simulations for various laser-plasma experiments.

direct drive↗
DOE OSTI2022

Evaluation of sky spectra and sky models in daylighting simulations

Sky models in daylight simulations represent the luminance variation across the sky-dome for different locations, dates, times and weather conditions, but skies are typically modelled as colourless. Recent studies explore techniques for incorporating the spectral content of daylighting in simulations. This paper provides an evaluation of the existing spectral sky models in lighting simulation software. The comparisons are made between the available mathematical sky models and naturally occurring skies that were recorded using high dynamic range photography and spectrophotometric measurements. The results show that recently developed sky models present progress compared to colourless sky models, but further research is needed to accurately simulate daylight spectra.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Using MARCUS, MICRE, and COMBLE data to improve understanding and modeling of cloud, aerosol, and boundary layer processes at high-latitudes

Because it is believed general circulation (GCM) and numerical weather prediction (NWP) models underestimate shortwave radiation over the Southern Ocean (SO) due to inadequate representations of boundary layer (BL) and cloud processes, it is critical to improve our understanding of key aerosol, cloud, precipitation, and BL processes. Over the north Atlantic Ocean (NA), cold-air outbreaks are common, yet few studies of the aerosol and environmental controls of the associated convective BL clouds exist as needed to develop and evaluate GCM representations. Although processes cannot be observed, cloud and aerosol properties can be measured in-situ or remotely retrieved, which when combined with numerical simulations enable process level understanding required to improve model representations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Conditional H I Mass Functions and the H I -to-halo Mass Relation in the Local Universe

We present a new H I mass estimator that relates log 10 (M HI /M * ) to a linear combination of four galaxy properties: stellar surface mass density, color index u - r, stellar mass, and concentration index, with the scatter of individual galaxies around the mean H I mass modeled with a Gaussian distribution function. We calibrate the estimator using the xGASS sample, including both H I detection and nondetection, and constrain the model parameters through Bayesian inferences. Tests with mock catalogs demonstrate that our estimator provides unbiased H I masses for optical samples like SDSS. We apply our estimator to the SDSS spectroscopic sample to estimate the H I mass function (HIMF) of local galaxies, as well as the conditional H I mass function in galaxy groups and the H I –halo mass relation. Our HIMF agrees with the ALFALFA measurements at M HI ≳ 5 × 10 9 M ⊙ , but with higher amplitude and a steeper slope at lower masses. We show that this discrepancy is caused primarily by the cosmic variance, which is corrected for the SDSS sample but not for ALFALFA. The total CHIMFs for all halo masses can be described by a single Schechter function, while those of central galaxies show a double-Gaussian profile. The total H I mass in a group increases monotonically with halo mass, but for central galaxies, the H I mass shows weak dependence on halo mass when M h ≳ 10 12 M ⊙ . The observed H I –halo mass relation is not reproduced by current hydrodynamic simulations and semianalytic models of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Baryonic imprints on DM haloes: the concentration–mass relation in the C amels simulations

ABSTRACT The physics of baryons in haloes, and their subsequent influence on the total matter phase space, has a rich phenomenology and must be well understood in order to pursue a vast set of questions in both cosmology and astrophysics. We use the Cosmology and Astrophysics with MachinE Learning Simulation (Camels) suite to quantify the impact of four different galaxy formation parameters/processes (as well as two cosmological parameters) on the concentration–mass relation, cvir−Mvir. We construct a simulation-informed non-linear model for concentration as a function of halo mass, redshift, and six cosmological/astrophysical parameters. This is done for two galaxy formation models, IllustrisTNG and Simba, using 1000 simulations of each. We extract the imprints of galaxy formation across a wide range in mass $M_{\rm vir}\in [10^{11}, 10^{14.5}] \, {\rm M}_\odot \, h^{-1}$ and in redshift z ∈ [0, 6] finding many strong mass- and redshift-dependent features. Comparisons between the IllustrisTNG and Simba results show the astrophysical model choices cause significant differences in the mass and redshift dependence of these baryon imprints. Finally, we use existing observational measurements of cvir−Mvir to provide rough limits on the four astrophysical parameters. Our non-linear model is made publicly available and can be used to include Camels-based baryon imprints in any halo model-based analysis.

Astronomy & Astrophysics↗