SEARCH · Engineering Papers

Results for “Analytical sensitivity analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Origin of core-to-core x-ray emission spectroscopy sensitivity to structural dynamics

Recently, coherent structural dynamics in the excited state of an iron photosensitizer was observed through oscillations in the intensity of Kα x-ray emission spectroscopy (XES). Understanding the origin of the unexpected sensitivity of core-to-core transitions to structural dynamics is important for further development of femtosecond time-resolved XES methods and, we believe, generally necessary for interpretation of XES signals from highly non-equilibrium structures that are ubiquitous in photophysics and photochemistry. Here, we use multiconfigurational wavefunction calculations combined with atomic theory to analyze the emission process in detail. The sensitivity of core-to-core transitions to structural dynamics is due to a shift of the minimum energy metal–ligand bond distance between 1s and 2p core-hole states. A key effect is the additional contraction of the non-bonding 3s and 3p orbitals in 1s core-hole states, which decreases electron– electron repulsion and increases overlap in the metal–ligand bonds. The effect is believed to be general and especially pronounced for systems with strong bonds. The important role of 3s and 3p orbitals is consistent with the analysis of radial charge and spin densities and can be connected to the negative chemical shift observed for many transition metal complexes. The XES sensitivity to structural dynamics can be optimized by tuning the emission energy spectrometer, with oscillations up to ±4% of the maximum intensity for the current system. The theoretical predictions can be used to design experiments that separate electronic and nuclear degrees of freedom in ultrafast excited state dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

MatLab Package for Whispering Gallery Mode Data

Dye-doped whispering gallery mode resonator (WGMR) microspheres yield highly structured emission spectra that are extremely sensitive to their environment and are of intense interest for use in a variety of sensing applications. Efforts to leverage the unique sensitivities of WGMRs have relied on stringent experimental requirements to correlate specific spectral shifts/changes to an analyte/stimulus such as 1) precise positional knowledge, 2) reference spectra for each microsphere, and 3) high mechanical stability. Consequently, these can hinder adequate mixing or incorporation of analytes and creates challenges for remote sensing. The MATLAB codes provided here are to be used in conjunction with a continuous flow technique for measuring WGM spectra of dye-doped microspheres suspended in solution. One MATLAB script, smooths the data, automatically baseline corrects it to isolate the whispering gallery modes (WGM) from the unwanted bulk emission, and assesses the similarity of each spectrum to aid in selecting a set of unique WGM spectra for further analysis. The next script is designed to analyze WGM spectra to determine the size of the resonator and the refractive index (RI) of its local environment without a priori knowledge of the individual microsphere. The final script allows the user to cluster spheres based on the product of their RI and radius and the contrast ratio of the RI of the sphere material and that of its environment by using a shared nearest neighbor spectral clustering algorithm.

Lilley, Laura↗

DOE OSTI2023

LDRD23-0184: Resilience and Hazard Risk Assessment to Prioritize Security Operations for Decisions and Impacts (RHAPSODI)

Recent examples provide a significant concern for the resilience of the U.S. electric grid and represent a need for enhanced decision-making to address an increasingly wide range of complex system interactions and potential consequences. In response, this LDRD project produced a proof-of-concept evaluation called the Resilience and Hazard Assessment to Prioritize Security Operations for Decisions and Impacts (RHAPSODI) methodology as an agile and flexible analytic framework capable of addressing multiple, diverse threats to desired electric grid performance. After empirically grounding needs for the future of U.S. electric grid resilience, this project employed the systems-theoretic process analysis (STPA) to develop a systems engineering risk model. The results of a completed feasibility study of a notional high voltage transmission system demonstrate an improved ability to incorporate both spatial (e.g., geographically distributed) and temporal (e.g., dynamic and time-dependent) elements of security risk to the gird. The success of this LDRD project provides the foundation for further evolution of the systems engineering risk model for the grid; derivation of quantitative approaches to evaluate risk and resilience performance; facilitation of agile experimenting and grid sensitivity to a range of vulnerabilities; and development of tools to assist decision-makers in enhancing U.S. electrical grid resilience.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

A new multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential processes under process model and parametric uncertainty

Process-based models have been widely used for hydrologic modeling, and it is a common practice to use sensitivity analysis methods for excluding non-influential hydrologic processes from further investigation and/or model improvement. This study develops a new method called multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential system processes and parameters. MMADS is conceptually similar to the Morris method for addressing parametric uncertainty, but has a unique feature to address both process model uncertainty (i.e., a process may be represented by multiple process models) and process model parameter uncertainty (i.e., parameters associated with a process model are random). MMADS first evaluates absolute differences of a quantity of interest (i.e., a system model output) by varying process models and/or process model parameter values, and then calculates the mean and variance of the differences for investigating process influence. The mean measures overall influence of the process on the quantity of interest, and the variance estimates influence of nonlinear effects of the process and/or its interactions with other processes. MMADS is an extension of the Morris method from a parameter space to a joint parameter-model space for explicitly addressing both process model uncertainty and model parameter uncertainty. The performance of MMADS is evaluated by using two numerical experiments. One experiment is based on Sobol’s G*-function with ten product elements, and has analytical solutions of the MMADS mean and variance of absolute differences. The other experiment is for groundwater flow modeling which considers three processes (i.e., recharge, geology, and snowmelt) that interact with each other. Finally, results indicate that MMADS is computationally efficient and can identify non-influential processes of complex hydrological systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Neutron Total Scattering Analysis of Materials for the Nuclear Fuel Cycle

Here, this paper reviews recent experimental efforts at the University of Tennessee and Oak Ridge National Laboratory to comprehensively characterize the structural details of materials relevant for the nuclear fuel cycle by employing advanced neutron scattering techniques. For the study of nuclear ceramics, neutron scattering offers distinct advantages over traditional laboratory or synchrotron X-ray diffraction, including enhanced sensitivity to elements with a low atomic mass, such as oxygen, nitrogen, and carbon. The key to these efforts is the recent advancement in the neutron scattering infrastructure at the high-flux diffractometers at the Spallation Neutron Source. The high neutron flux at these instruments enables neutron total scattering, a nondestructive bulk technique that simultaneously captures both short-range structural effects through pair distribution function analysis and long-range order through diffraction pattern analysis. This approach is particularly important for a comprehensive description of defective, disordered, or amorphous nuclear materials. The case studies presented here include analyses of the local defect structure in hyperstoichiometric uranium oxides and short-range order of ion-irradiated ceramics. This advanced analytical methodology will improve our understanding of the behavior of materials in extreme environments and contribute to the development of more resilient nuclear materials.

Neutron scattering↗

DOE OSTI2022

Comparative Technoeconomic Analysis of Pathways for Electrochemical Reduction of CO 2 with Methanol to Produce Methyl Formate

Electrochemical CO 2 reduction has promise as a technology that could help society reach carbon neutrality while producing valuable fuels and chemicals. Herein, the electrochemical synthesis of methyl formate, a product not observed in aqueous CO 2 electrolysis, has been analyzed by a rigorous technoeconomic model to evaluate its commercial viability. Methyl formate synthesis has been demonstrated with high faradaic efficiency through the electroreduction of CO 2 in methanol. Four competing approaches were analyzed: (1) Electroreduction of captured CO 2 in a dual CH 3 OH/H 2 O electrolyzer, (2) Direct electroreduction of flue gas CO 2 in a dual CH 3 OH/H 2 O electrolyzer, (3) Electroreduction of captured CO 2 in a CH 3 OH/CH 3 OH electrolyzer, and (4) Electroreduction of captured CO 2 in a H 2 O/H 2 O electrolyzer with a downstream CH 3 OH reactor. Sensitivity analyses, cost contour plots, and comparison plots were generated. The dual methanol/water electrolysis approach was the most cost-competitive, with a levelized cost of methyl formate below the present market price. Here, the all-methanol electrolysis route was more expensive due to increased methanol consumption and greater distillation costs. Methyl formate production through aqueous CO 2 electrolysis to formic acid with a secondary esterification reaction was by far the most expensive approach, primarily due to the energy-intensive nature of distilling formic acid from water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Advances in mass spectrometry-enabled multiomics at single-cell resolution

We report biological organisms are multifaceted, intricate systems where slight perturbations can result in extensive changes in gene expression, protein abundance and/or activity, and metabolic flux. These changes occur at different timescales, spatially across cells of heterogeneous origins, and within single-cells. Hence, multimodal measurements at the smallest biological scales are necessary to capture dynamic changes in heterogeneous biological systems. Of the analytical techniques used to measure biomolecules, mass spectrometry (MS) has proven to be a powerful option due to its sensitivity, robustness, and flexibility with regard to the breadth of biomolecules that can be analyzed. Recently, many studies have coupled MS to other analytical techniques with the goal of measuring multiple modalities from the same single-cell. It is with these concepts in mind that we focus this review on MS-enabled multiomic measurements at single-cell or near-single- cell resolution.

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Enabling Cybersecurity, Situational Awareness and Resilience in Distribution Grids with High Penetration of Photovoltaics (CARE-PV) (Final Report)

Since legacy distribution systems have very limited visibility beyond the substation, high penetration of PV at the grid edge presents some unique operational challenges. One approach to address these challenges is to use information from advanced metering infrastructure (AMI) and µPMUs. However, exploiting this information is impacted by a number of factors, including multi-timescale measurements, volume of data generated, communication network impairments (e.g., information loss and latency) and susceptibility to cyber-attacks. Therefore, one of the critical tasks involved in the management of a distribution grid is to develop complete situational awareness by integrating cyber-security mechanisms with state estimation strategies and leveraging this situational awareness to assure energy services at strategic locations while exploiting AMI/PV inverter/ µPMU data. This CARE-PV project addresses the fundamental challenges in situational awareness and resilience to cyber and physical vectors by exploiting the synergy between innovative modeling, estimation, data analytics, testing and validation using smart PV inverters designed at K-State and facilities at NREL. Specifically, the project involved the development, testing and validation of the following novel enabling technologies: (Thrust 1) Resilience to cyber vectors that impact data integrity was addressed via a two-level defense strategy that combines cyber intrusion detection using self-learning, cooperative smart PV inverters, and a novel moving target defense framework to combat data integrity attacks. (Thrust 2) Resilience to cyber-physical vectors that impact situational awareness by limiting data availability was addressed via novel centralized and decentralized, sparsity-based static and dynamic state estimation approaches that enhance observability even when the underlying system is unobservable. (Thrust 3) Leveraging a unique probabilistic sensitivity analysis approach accompanied by one-of-a-kind dominant influencer set computation, the vulnerability of critical infrastructure at strategic locations was evaluated so that proactive PV-based control strategies can be used to support operations under normal/outage scenarios. These CARE-PV project innovations were demonstrated on both small-scale IEEE and larger utility-scale testbeds (Thrust 4). Feedback from Industry Advisory Board members was used to formulate a commercialization pathway for a subset of CARE-PV technologies. These CARE-PV technologies will ultimately lead to reliable and secure, large-scale integration of renewable energy and mitigate the risk of energy disruption resulting from cyber incidents and other emerging threats within the energy environment.

14 SOLAR ENERGY↗

DOE OSTI2020

Single Bimodular Sensor for Differentiated Detection of Multiple Oxidative Gases

Semiconductive metal-oxide sensors suffer from cross-sensitivities under mixed chemical condition, specifically upon mixture of multiple oxidative or reductive gases. Herein, a single bimodular sensor is demonstrated for smart differentiation of multiple oxidative analytes by relating the resistance-metric mode to impedance-metric mode. The sensor construct based on ZnO nanorods readily outputs three response datasets upon exposure of oxidative-gas mixture including O 2 , SO 2 , and NO 2 , the resistance, real part impedance, and imaginary part impedance. The differentiative and correlated nature between these response signals allows such a single sensor platform to differentiate these oxidative gases accurately and robustly. Linear and non-linear decision boundaries are established over a large gas-concentration range from 2 ppm to 3% through a combination of principal component analysis and artificial neural network training. A facile user interface is demonstrated for recognition and measurement of unknown gas analytes, with the error of the predicted analyte-concentration as low as 2%.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Solving the structure of “single-atom” catalysts using machine learning – assisted XANES analysis

We show that "single-atom” catalysts (SACs) have demonstrated excellent activity and selectivity in challenging chemical transformations such as photocatalytic CO 2 reduction. For heterogeneous photocatalytic SAC systems, it is essential to obtain sufficient information of their structure at the atomic level in order to understand reaction mechanisms. In this work, a SAC was prepared by grafting a molecular cobalt catalyst on a light-absorbing carbon nitride surface. Due to the sensitivity of the X-ray absorption near edge structure (XANES) spectra to subtle variances in the Co SAC structure in reaction conditions, different machine learning (ML) methods, including principal component analysis, K-means clustering, and neural network (NN), were utilized for in situ Co XANES data analysis. As a result, we obtained quantitative structural information of the SAC nearest atomic environment thereby extending the NN-XANES approach previously demonstrated for nanoparticles and size-selective clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Interpreting and Stabilizing Machine-Learning Parametrizations of Convection

Neural networks are a promising technique for parameterizing subgrid-scale physics (e.g., moist atmospheric convection) in coarse-resolution climate models, but their lack of interpretability and reliability prevents widespread adoption. For instance, it is not fully understood why neural network parameterizations often cause dramatic instability when coupled to atmospheric fluid dynamics. This paper introduces tools for interpreting their behavior that are customized to the parameterization task. First, we assess the nonlinear sensitivity of a neural network to lower-tropospheric stability and the midtropospheric moisture, two widely studied controls of moist convection. Second, we couple the linearized response functions of these neural networks to simplified gravity wave dynamics, and analytically diagnose the corresponding phase speeds, growth rates, wavelengths, and spatial structures. To demonstrate their versatility, these techniques are tested on two sets of neural networks, one trained with a super-parameterized version of the Community Atmosphere Model (SPCAM) and the second with a near-global cloud-resolving model (GCRM). Additionally, even though the SPCAM simulation has a warmer climate than the cloud-resolving model, both neural networks predict stronger heating/drying in moist and unstable environments, which is consistent with observations. Moreover, the spectral analysis can predict that instability occurs when GCMs are coupled to networks that support gravity waves that are unstable and have phase speeds larger than 5 m s -1 . In contrast, standing unstable modes do not cause catastrophic instability. Using these tools, differences between the SPCAM-trained versus GCRM-trained neural networks are analyzed, and strategies to incrementally improve both of their coupled online performance unveiled.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Implementation and Verification of Cable Bending Stiffness in MoorDyn: Preprint

Cable bending stiffness capabilities have been added to MoorDyn to enable modeling of dynamic power cables. The relatively large motions experienced by floating wind turbines and wave energy converters pose a challenge for power cables, whose internal components provide significant bending resistance and are sensitive to deformation. The behavior and associated design considerations of power cables in these highly dynamic applications make coupled analysis relevant for design. MoorDyn is a lumped-mass mooring dynamics model that is used for floating wind energy simulation (often coupled with OpenFAST) and floating wave energy converter simulation (often coupled with WEC-Sim). Its existing line elasticity formulation considers axial stiffness only. To properly capture the dynamics of power cables, a cable bending stiffness model has been added that approximates cable curvature based on the difference in tangent vectors of adjacent line elements. The resulting bending moment is realized by applying forces on adjacent nodes, leaving the underlying lumped-mass formulation unchanged. The bending stiffness implementation is verified in static conditions against analytical solutions and then in a dynamic power cable scenario in comparison with OrcaFlex. The dynamic scenario uses prescribed motions and includes wave loadings on the cable. Results indicate correct implementation of bending stiffness and show close agreement with OrcaFlex.

50 EE - Wind and Water Power Program - Water (EE-4↗

DOE OSTI2024

An Open-Source Decarbonization Analytics Framework: Designing for Low-Carbon Emission Districts and Communities: Preprint

This paper introduces an open-source analytics framework designed to assist in creating low or net-zero carbon buildings and urban districts. Integrated within URBANopt, an open-source platform for energy analysis in districts and communities, this framework equips researchers, architects, engineers, and other stakeholders with tools to evaluate the carbon footprint implications of their design choices. The framework enables the analysis of various scenarios, incorporating both historical and future emission factors, and can span across different climate zones, each with distinct grid and emissions characteristics. The results showcase the framework's capability to evaluate the impact of design upgrades and control strategies on carbon emissions in districts and communities. An illustrative analysis using a hypothetical district in Denver, Colorado, shows reduced emissions from energy efficiency upgrades and control strategies, highlighting the sensitivity in their effects on emissions and energy use.

buildings energy efficiency↗

DOE OSTI2023

Comparative Techno-economic Analysis and Life Cycle Assessment of Producing High-Value Chemicals and Fuels from Waste Plastic via Conventional Pyrolysis and Thermal Oxo-degradation

The rapid rise in global plastic production in recent decades has resulted in the massive generation of plastic waste. Over 75% of the plastic waste generated in the United States was sent to landfills, with a meager 8.7% recycled. Plastics are valuable feedstocks for platform chemicals and fuels. Chemical upcycling of waste high-density polyethylene (HDPE) is gaining more attention as a potentially feasible and environmentally friendly plastic waste management technology. Conventional pyrolysis (CPY) and thermal oxo-degradation (TOD) are two chemical upcycling technologies actively researched for decomposing waste HDPE into valuable chemicals and fuels. However, there are few studies on the techno-economic analysis (TEA) and life cycle assessment (LCA) of these technologies for converting waste HDPE to valuable products. This study conducts a comparative TEA and LCA study of the thermochemical decomposition of waste HDPE to produce gaseous (ethylene and propylene) and liquid (naphtha, diesel, and wax) products by CPY and TOD. The study elucidates and compares the impact of hydrocracking longer chain hydrocarbons to produce more valuable products on the TEA and LCA. The TEA showed that the fixed capital investment could range from $\$32.5$ million for TOD without hydrocracking to $\$244$ million for CPY with hydrocracking scenarios. Annual revenues range from $\$28.1$ million to $\$71.5$ million in favor of scenarios with hydrocracking. However, the net present value ranges from $\$1.4$ million to $\$265.8$ million in favor of scenarios without hydrocracking. Sensitivity analysis showed that fixed capital cost, facility capacity, and product prices have the biggest impact on the process economics of the facilities, while utilities and waste transportation to refineries have the biggest impact on environmental impacts. Here, the LCA showed that primary products from scenarios without hydrocracking can be more environmentally friendly than virgin products from petroleum processes. However, TOD and CPY with hydrocracking primary products have more emissions than those of virgin products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗

DOE OSTI2026

Parity-odd four-point correlation function from the DESI data release 1 luminous red galaxy sample

The parity-odd four-point function provides a unique probe of fundamental symmetries and potential new physics in the large-scale structure of the Universe. We present measurements of the parity-odd four-point function using the Dark Energy Spectroscopic Instrument (DESI) DR1 luminous red galaxy (LRG) sample and assess its detection significance. Our analysis considers both auto- and cross-correlations, using two complementary approaches to the covariance: (i) the full analytic covariance matrix applied to the uncompressed data vector, and (ii) a compressed data vector combined with a hybrid covariance matrix constructed from simulations and analytic estimates. When using the full analytic covariance matrix without corrections, we observe apparent auto-correlation signals with significance up to 4⁢𝜎. However, this excess is also consistent with a mismatch between the statistical fluctuations estimated from the simulations and those present in the real data. Our findings therefore suggest that the parity-odd signal in the current DESI DR1 LRG sample is consistent with zero. We note, however, that the low completeness of this sample may have a non-negligible impact on the detection sensitivity. Future data releases with improved completeness will be crucial for further investigation.

Hou, Jiamin [Ludwig-Maximilians-Universität; Unive↗

DOE OSTI2022

Study of Radiative Heat Transfer and Flow Physics from Medium-scale Methanol Pool Fire Simulations [Slides]

Analysis of methanol pool fire conducted as part of validation study for SIERRA/Fuego. Radiation model was effectively calibrated by modifying radiation model parameters for methanol. Computing integrated buoyancy flux, entrainment rate, and turbulent kinetic energy allowed for evaluation of less typical quantities in this validation study. Quantities were compared with experimental data or correlations and generally showed agreement. TKE needed fine mesh to be computed accurately. Predicted flame height less sensitive to variations in mixture fraction than temperature. Mixture fraction is a preferable threshold variable for this application.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Probing plasma physics with spectral index maps of accreting black holes on event horizon scales

The Event Horizon Telescope (EHT) collaboration has produced the first resolved images of the supermassive black holes at the centre of our galaxy and at the centre of the elliptical galaxy M87. As both technology and analysis pipelines improve, it will soon become possible to produce spectral index maps of black hole accretion flows on event horizon scales. Here, in this work, we predict spectral index maps of both M87* and Sgr A* by applying the general relativistic radiative transfer (GRRT) code IPOLE to a suite of general relativistic magnetohydrodynamic (GRMHD) simulations. We analytically show that the spectral index increases with increasing magnetic field strength, electron temperature, and optical depth. Consequently, spectral index maps grow more negative with increasing radius in almost all models, since all of these quantities tend to be maximized near the event horizon. Additionally, photon ring geodesics exhibit more positive spectral indices, since they sample the innermost regions of the accretion flow with the most extreme plasma conditions. Spectral index maps are sensitive to highly uncertain plasma heating prescriptions (the electron temperature and distribution function). However, if our understanding of these aspects of plasma physics can be tightened, even the spatially unresolved spectral index around 230 GHz can be used to discriminate between models. In particular, Standard and Normal Evolution (SANE) flows tend to exhibit more negative spectral indices than Magnetically Arrested Disc (MAD) flows due to differences in the characteristic magnetic field strength and temperature of emitting plasma.

(magnetohydrodynamics) MHD↗