Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “actual”

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 199 records · Page 11

Cesium Removal from 5.5 and 7.0 M Na AP-105 Using Crystalline Silicotitanate

The Tank Side Cesium Removal (TSCR) system, currently operational by Washington River Protection Solutions LLC (WRPS), prepares initial low-activity Hanford waste tank supernate feeds for the Hanford Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Facility. In addition to entrained solids removal from the supernate, the primary goal of TSCR is to remove cesium-137 ( 137 Cs) by ion exchange, allowing contact handling of the liquid effluent product at the WTP as governed by a waste acceptance criterion (WAC). Specific to 137 Cs, this requirement is < 3.18E-5 Ci 137 Cs/mole of Na. Crystalline silicotitanate (CST) ion exchange media, manufactured by Honeywell UOP, LLC (product IONSIVTM R9140-B), has been selected as the ion exchange media for TSCR. CST is a non-elutable inorganic material that has demonstrated robust chemical, physical, and radiation tolerance while maintaining functionality. However, testing to date on actual tank waste samples has been limited to Na concentrations between 5 and 6 M Na (Fiskum et al. 2019a, 2021b and Westesen et al. 2021a, 2021b, 2022) while actual tank conditions can reach upwards of 9 M Na. Testing with feed from Hanford tank AP-105 incorporated testing at both 5.5 and 7 M Na in order to evaluate the impact of Na concentration on volume of waste processed before reaching the WAC. A 6-L volume of 5.5 M Na AP-105 and an 8-volume of 7 M Na AP-105 was processed through the Radioactive Waste Test Platform system, established at Pacific Northwest National Laboratory to support small-scale waste qualification efforts. The columns consisted of 6- to 9-mL CST beds (CST Lot 2002009604, sieved to screen out >30-mesh particles) placed in 1.5-cm-inner-diameter columns. Feed was processed at 1.9 bed volumes (BV) per hour; the flowrate, in terms of contact time with the CST bed, matched the expected flowrate at TSCR. Table S.1 and Figure S.1 summarize the measured AP-105 Cs load performance for each feed condition.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Comparing Occupant Light Exposure in Factory and Office Spaces

There is an increasing desire to provide healthy workplace environments that better serve people while minimizing the impact on the planet. Building recommendations help designers to incorporate the latest research and best practices, but translating research studies conducted in different settings to broader recommendations is challenging. This is particularly true for lighting, with researchers continuing to learn more about circadian, neuroendocrine, and neurobehavioral responses to light. The recommendations are regularly adjusting to incorporate the latest research as well as feedback from practitioners implementing these recommendations. There remains a need to better document the actual light exposure that people experience over the course of a day and across seasons, so that researchers can better understand the influence of light on people and recommendations can be refined. Thankfully, advances in sensor technology are making this more possible than ever. PNNL collaborated with ETC Inc. to understand actual occupant light exposure in factory and office spaces at ETC’s Headquarters in Middleton, WI. The primary building at ETC headquarters opened in 2004 and houses approximately 800 employees that contribute to central business functions, customer support, research and development, as well as the company’s primary 60,000 ft2 manufacturing facility. Vertical light exposure measurements including electric light and daylight were recorded in factory, research and development, and office work areas at occupant eye height. This document will focus on lighting measurements collected in six factory and office locations across work location, season, and time of day. Ultimately, detailed lighting measurements can be coupled with employee surveys to support development of satisfactory work environments.

42 ENGINEERING↗

Evaluation of Naturally Occurring Be-7 Activity on LANL Stack Filter Clumps

Los Alamos National Laboratory (LANL) is responsible for providing solutions to challenging issues through multidisciplinary science, engineering, and technology that impact national security. As part of the Radioactive Air Emissions Management (RAEM) Team, it is our responsibility to ensure LANL operates within the requirements of the Code of Federal Regulations 40 CFR 61, Subpart H, National Emission Standards for Emissions of Radionuclides Other Than Radon From Department of Energy Facilities (Rad-NESHAP). The standard states, “Emissions of radionuclides to the ambient air from Department of Energy facilities shall not exceed those amounts that would cause any member of the public to receive in any year an effective dose equivalent of 10 mrem/yr” (EPA, 2002) During a routine assessment of LANL’s Rad-NESHAP compliance program, an opportunity for improvement was identified that addresses the evaluation of routine stack samples which potentially have naturally occurring radioactive materials (NORM) on the filter media. In the past and based on professional judgment, detections of NORM such as K-40 and Be-7 typically were rejected from the analytical results. The assessment finding suggested developing a technical criterion to ensure that such rejections were justified, comparing actual measured values with anticipated NORM concentrations for various stacks located at LANL. The significance of accepting or rejecting samples as “true” or “NORM” and its impact on the RAEM’s mission and compliance to the regulations is explained later in the document. Additionally, this document will describe the process and calculations for how such results are to be accepted as actual samples “true positive” or rejected as NORM.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Hanford Waste Treatment Plant Low Activity Waste Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation (Rev.1)

The Hanford Tank Waste Treatment and Immobilization Plant Low Activity Waste (LAW) facility stack monitor locations were qualified using scale model stacks to mitigate the risk of discovering that sampling locations do not meet the qualification criteria on the full-scale stacks. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling location were geometrically similar to the actual stack and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the LAW stacks were performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing, determines the range of stack flow rates for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. For this analysis, the range of qualified flow rates listed is conservatively based on the average DV through 6 DV for LV-S1, LV-S2, and LV-C2, and 1/3 DV to 3 DV for LV-S3. Table S1 lists the operating flow rates along with the conservative lower and upper qualified stack flow rates for each of the LAW facility stacks. For each stack, the operating flow is below the upper qualified stack.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Predicting seismic amplitudes with machine learning

The accurate estimation of seismic wave amplitude is vital to precisely determine the yield, magnitude, and event discrimination possible for a given network – a critical element in nuclear explosion monitoring. This task is complicated by several factors, including but not limited to radiation pattern, scattering effects, and crustal variations, which can lead to the attenuation or amplification of amplitude along a given raypath. In this report, we explore the novel application of machine learning to the task of seismic amplitude estimation by training a simple Artificial Neural Network (ANN) on an S-wave amplitude dataset from Lai et al. (2019). Attributes from this dataset used as input to the ANN included event-station distances, station locations (latitude, longitude), event locations (latitude, longitude), event depths, event magnitudes, radiation patterns, signal-to noise ratio (SNR) measurements (average-amplitude, peak-to-trough, maximum peak), and signal periods. We find that the trained ANN predicts S-wave amplitudes with a modest tendency toward underestimating the actual values, as indicated by a linear regression between predicted and actual data (slope: 0.892, intercept: -0.651). These results suggest that an ANN can perform this task, with potential for significant improvements through improved datasets, architectures, and parameter tuning.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Integration and commissioning plan of Full Flow Purifier at Muon Campus

The Full Flow Purifier for Fermilab’s Muon Campus uses a charcoal bed surrounded by a liquid nitrogen jacket to purify up to 240 g/s of helium gas. Fabrication by Ability Engineering Technology Inc. has been completed and the purifier delivered to Fermilab. It is the largest purifier to be used at Fermilab based on both capacity and size. A previous paper discussed the design of purifier for various operational conditions and horizontal shipping. The purifier is designed to withstand 5g force in vertical and 2g force in lateral and longitudinal directions. Transportation experience from vendor to Fermilab and within site is discussed. Integration of the purifier involved design and fabrication of a liquid nitrogen transfer line, regeneration system, and helium piping to connect it to Muon Campus cryogenic system. It also involved establishing electrical, instrumentation and controls connections to the system. Integration of the purifier is discussed in detail. The purifier is to be commissioned using up to 4 MYCOM helium compressors to supply helium and liquid nitrogen supplied by a 15,000 -gallon tank. Actual effectiveness of the 3-stream heat exchanger is to be estimated based on measured temperatures and flow rate. Impurity levels shall be monitored at inlet and outlet of the purifier. Theoretical adsorption capacity of the purifier is calculated based on the measured temperatures and flowrates and is compared to actual adsorption capacity over time.

Subedi, Jeewan [Fermilab]↗

An Advanced Machine Learning and Artificial Intelligence System for Demonstrating Radiation Regulatory Compliance in DOE Accelerator Facilities

In this Phase II proposal, Applied Research LLC (ARLLC), Thomas Jefferson National Accelerator Facility (Jefferson Lab), and Old Dominion University (ODU) propose the combination of domain knowledge (beam characteristics, fixed structural shielding, earthen burden (the soil and foliage added to the dome of the experimental halls as additional shielding), etc.), machine learning (ML) and/or artificial intelligence (AI) to correlate a variety of multi-modal onsite signals and the radiation fields seen in accessible areas of the accelerator site and the site boundary. The ML/AI will consider the complex influence of environmental parameters affecting the radon contribution of the measurements, focusing on actual data obtained from Jefferson Lab. In Phase I, the coded beam and location data were fed into a deep learning model to predict doses at several designated locations in Jefferson Lab’s facility. Moreover, a dense radiation map was generated using only a sparse collection of the samples in a facility. In Phase II, we will develop a software prototype containing a radiation prediction algorithm, dense radiation map algorithms, and background noise prediction algorithms, with actual data used to evaluate the prototype. This work will provide a framework for evaluation of radiation measurement results around the site based on learned responses. In addition, the proposed approach allows more granular mapping of radiation levels. Better understanding and communication of these levels is related to the overall approach in keeping doses to personnel ALARA.

43 PARTICLE ACCELERATORS↗

Evaluation of the Energy, Hygrothermal, and Thermal Capacity Performance of Cross-Laminated Timber

Cross-laminated timber (CLT) construction is gaining momentum in the US because it offers multiple advantages over traditional construction methods. Benefits that have received the most attention focus on constructability, the environment, and protection (e.g., blast resistance), although CLT construction is likely to offer other benefits, as well. Still, these have not been studied at length because such evaluations are costly, requiring long-term assessments in an actual building and specialized technical knowledge. Among the possible benefits, CLT construction likely provides a higher-performing building envelope. Using CLT panels to enclose a building means fewer joints in the opaque envelope than what is required in traditional stick-framed construction. Fewer joints mean fewer locations where the air- and water-resistive barrier (WRB) could be compromised; thus, a CLT building enclosure may require less maintenance and have a longer lifespan than a traditionally built structure because of fewer air and water leaks. In addition, CLT’s thermal mass moderates indoor temperatures, allowing the heating, ventilation, and air conditioning (HVAC) system to operate more efficiently during peak hours, reducing operational energy consumption throughout the lifetime of the CLT building (Salonvaara et al., 2022). Furthermore, more stable indoor temperatures can increase occupant comfort. The CLT’s thermal mass can also reduce energy costs by adjusting to utility time-of-use pricing without affecting occupant comfort. The ability of CLT buildings to bridge periods without HVAC operation prepares them for future grid interaction and provides a certain level of resilience against power outages. Researchers have attempted to quantify these benefits; however, their work is based on simplified simulations with numerous assumptions. To correctly understand the benefits, an actual building must be monitored. Therefore, information needs to be gathered on indoor and outdoor temperatures, HVAC energy consumption, thermostat setpoints, temperatures, and thermal transport in CLT components to comprehend how these parameters are affected by the CLT’s thermal mass. These data are needed to reduce the number of assumptions and calibrate simulation models to optimize HVAC controls to minimize overall energy consumption, reduce energy use and higher fees during peak demand, and maintain occupant comfort. Additionally, the calibrated simulation model allows the optimization exercise to be repeated in various US climates. Potential benefits can be tailored to buildings in various locations, and decisions can be made on where CLT construction could be most advantageous. Furthermore, monitoring and simulation results are needed to evaluate the durability of the CLT structures in different climates. This project’s researchers gathered information to help understand and quantify the benefits of CLT buildings concerning operational energy, moderated indoor temperatures, and comfort; the dynamic operation to provide grid services; and resilience in times of power outage. Through the corroboration of simulation models with real-world measurements, this study paves the way for extrapolating findings to other climatic zones and building typologies, thereby broadening the understanding of CLT’s multifaceted benefits and reinforcing its position as a material of choice in sustainable construction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Frequency-Nadir-Constrained Unit Commitment for Low-Inertia, High-IBR Island Power Systems [Slides]

The process of energy decarbonization in island power systems is accelerated due to the swift integration of inverter-based renewable energy resources (IBRs). The unique features of such systems, including rapid frequency changes resulting from potential generation outages or imbalances due to the unpredictability of renewable power, pose a significant challenge in maintaining the frequency nadir without external support. This paper presents a unit commitment (UC) model with data-driven frequency nadir constraints, including either frequency nadir or minimum inertia requirements, helping to limit frequency deviations after significant generator outages. The constraints are formulated using a linear regression model that takes advantage of real-world, year-long generation scheduling and dynamic simulation data. The efficacy of the proposed UC model is verified through a year-long simulation in an actual island power system using historical weather data. The alternative minimum inertia constraint, derived from actual system operation assumptions, is also evaluated. Findings demonstrate that the proposed frequency nadir constraint notably improves the system's frequency nadir under high photovoltaic (PV) penetration levels, albeit with a slight increase in generation costs, when compared to the alternative minimum inertia constraint.

14 SOLAR ENERGY↗

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗

5-minute Wind Power Data based on WFIP2 WRF Simulation

The second Wind Forecast Improvement Project (WFIP2) was a public-private partnership funded by the U.S. Department of Energy and NOAA, aimed at enhancing the forecast skill of numerical weather prediction models for turbine-height winds in regions with complex terrain. An 18-month Weather Research and Forecasting (WRF) model simulation was conducted over the Pacific Northwest, with model outputs validated against observational data collected during WFIP2. Simulated wind speeds were used to estimate wind power generation using reV (the Renewable Energy Potential Model developed by NREL) at ten wind project sites. Two sets of results were produced: one using wind speeds extracted from the model grid cell at the project centroid, and another using wind speeds from the actual turbine locations. For each dataset, power output was calculated using both actual turbine-specific power curves and nine generic power curves to convert wind speed into power.

17 WIND ENERGY↗

HydroFish: freshwater fish co-occurrence with hydropower plants and non-powered dams in conterminous United States sub-basins

The HydroFish dataset lists all existing hydropower plants (EHAs) and non-powered dams (NPDs; ≥ 0.001 MW potential nominal capacity), delineates the hydrologic sub-basins in which they are situated, and then lists all freshwater fish species reported to occur in those sub-basins. This dataset was compiled using the HydroBio dataset (https://hydrosource.ornl.gov/data/datasets/hydrobio/) and contains 24 total variables that describe hydrologic sub-basins, each unique EHA (plant ID and name, geographic coordinates, permit type, capacity, etc.) and NPD (ID value, known names, geographic coordinates, and estimated potential nominal capacity), and freshwater fish species in the sub-basin (common and scientific name, origin, and migratory and threat status). The HydroBio dataset was built using Oak Ridge National Laboratory’s Existing Hydropower Assets Dataset (2024 version) and Non-powered Dam Technical Potential Dataset (2024 version), and NatureServe’s fish species distribution dataset (2023 version). The dataset also contains summary variables that report the unique number of EHAs, NPDs, and freshwater fish species per sub-basin. The HydroFish dataset contains two unique data files: 1) a .csv metadata file describing the dataset variables, and 2) a .csv data file containing the actual dataset. Note that there may be many rows per unique existing hydropower plant or non-powered dam given that distinct species are listed per existing plant or NPD per sub-basin. The dataset is downloadable as a zip file containing the metadata and dataset files.

Bozeman, Bryan [Oak Ridge National Laboratory (ORN↗

Reinforcement Learning-based Output Structured Feedback for Distributed Multi-Area Power System Frequency Control

Load frequency control (LFC) is a key factor to maintain the stable frequency in multi-area power systems. As the modern power systems evolve from centralized to decentralized paradigm, LFC needs to consider the decentralized scheme that considers limited information from the information-exchange graph for the generator control of each interconnected area. This paper aims to solve a data-driven constrained LQR problem with mean-variance risk constraints and output structured feedback, and applies this framework to solve the LFC problem in multi-area power systems. By reformulating the constrained optimization problem into a minimax problem, the stochastic gradient descent max-oracle (SGDmax) algorithm with zero-order policy gradient (ZOPG) is adopted to find the optimal feedback gain from the learning, while guaranteeing the convergence. In addition, to improve the adaptation of the proposed learning method to new or varying models, we construct an emulator grid that approximates the dynamics of a physical grid and performs training based on this model. Once the feedback gain is obtained from the emulator grid, it is applied to the physical grid with a robustness test to check whether the controller from the approximated emulator applies to the actual system. Numerical tests show that the obtained feedback controller can successfully control the frequency of each area, while mitigating the uncertainty from the loads, with reliable robustness that ensures the adaptability of the obtained feedback gain to the actual physical grid.

Kwon, Kyung-bin↗

Uncontrolled Electric Vehicle Charging Impacts on Distribution Electric Power Systems with Primarily Residential, Commercial or Industrial Loads

An increase in Electric Vehicles (EV) will result in higher demands on the distribution electric power systems (EPS) which may result in thermal line overloading and low voltage violations. To understand the impact, this work simulates two EV charging scenarios (home- and work-dominant) under potential 2030 EV adoption levels on 10 actual distribution feeders that support residential, commercial, and industrial loads. The simulations include actual driving patterns of existing (non-EV) vehicles taken from global positioning system (GPS) data. The GPS driving behaviors, which explain the spatial and temporal EV charging demands, provide information on each vehicles travel distance, dwell locations, and dwell durations. Then, the EPS simulations incorporate the EV charging demands to calculate the power flow across the feeder. Simulation results show that voltage impacts are modest (less than 0.01 p.u.), likely due to robust feeder designs and the models only represent the high-voltage (“primary”) system components. Line loading impacts are more noticeable, with a maximum increase of about 15%. Additionally, the feeder peak load times experience a slight shift for residential and mixed feeders (≈1 h), not at all for the industrial, and 8 h for the commercial feeder.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Factors Impacting Range Hood Use in California Houses and Low-Income Apartments

Venting range hoods can control indoor air pollutants emitted during residential cooktop and oven cooking. To quantify their potential benefits, it is important to know how frequently and under what conditions range hoods are operated during cooking. We analyzed data from 54 single family houses and 17 low-income apartments in California in which cooking activities, range hood use, and fine particulate matter (PM 2.5 ) were monitored for one week per home. Range hoods were used for 36% of cooking events in houses and 28% in apartments. The frequency of hood use increased with cooking frequency across homes. In both houses and apartments, the likelihood of hood use during a cooking event increased with the duration of cooktop burner use, but not with the duration of oven use. Actual hood use rates were higher in the homes of participants who self-reported more frequent use in a pre-study survey, but actual use was far lower than self-reported frequency. Residents in single family houses used range hoods more often when cooking caused a discernible increase in PM 2.5 . In apartments, residents used their range hood more often only when high concentrations of PM 2.5 were generated during cooking.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Adaptively Learned Modeling for a Digital Twin of Hydropower Turbines with Application to a Pilot Testing System

In the development of a digital twin (DT) for hydropower turbines, dynamic modeling of the system (e.g., penstock, turbine, speed control) is crucial, along with all the necessary data interface, virtualization, and dashboard designs. Since the DT must mimic the actual dynamics of the hydropower turbine accurately, adaptive learning is required to train these dynamic models online so that the models in the DT can effectively follow the representation of the actual hydropower turbine dynamics accurately and reliably. This study presents an adaptive learning method for obtaining the hydropower turbine models for DT development of hydropower systems using the recursive least squares algorithm. To simplify the formulation, the hydropower turbine under consideration was assumed to operate near a fixed operating point, where the system dynamics can be well represented by a set of linear differential equations with constant parameters. In this context, the well-known six-coefficient model for the Francis turbine was formulated as the starting point to obtain input and output models for the turbine. Then, an adaptive learning mechanism was developed to learn model parameters using real-time data from a hydropower turbine testing system. This led to semi-physical modeling, in which first principles and data-driven modeling are integrated to produce dynamic models for DT development. Applications to a pilot system at the Norwegian University of Science and Technology (NTNU) were made, and the models learned adaptively using the data collected from the university’s pilot system. Desired modeling and validation results were obtained.

13 HYDRO ENERGY↗

Inverse Compton Cooling in the Coronae of Simulated Black Hole Accretion Flows

Here we present a formulation for a local cooling function to be employed in the diffuse, hot corona region of 3D GRMHD simulations of accreting black holes. This new cooling function calculates the cooling rate due to inverse Compton scattering by considering the relevant microphysics in each cell in the corona and approximating the radiation energy density and Compton temperature thereby integrating over the thermal seed photon flux from the disk surface. The method either assumes the ion and electron temperatures are equal (1T) or calculates them separately (2T) using an instantaneous equilibrium approach predicated on the actual relevant rate equations (Coulomb and Compton). The method is shown to be consistent with a more detailed ray-tracing calculation where the bulk of the cooling occurs, but is substantially less costly to perform. As an example, we apply these methods to a harm3d simulation of a 10M⊙, non-spinning black hole, accreting at nominally 1% the Eddington value. This new approach leads to radiative efficiency values sime65% above Novikov–Thorne, with a larger fraction of total cooling in the corona as compared to simulations performed using the original target-temperature cooling function. Time-averaged post-processing reveals that the continuum spectral observations predicted from these simulations are qualitatively similar to actual X-ray binary data, especially so for the 1T approach, which yields a harder power-law component (Γ = 2.25) compared to the 2T version (Γ = 2.53).

79 ASTRONOMY AND ASTROPHYSICS↗

Can Cooling and Heating Functions Be Modeled with Homogeneous Radiation Fields?

Abstract Cooling and heating functions describe how radiative processes impact the thermal state of a gas as a function of its temperature and other physical properties. In a most general case the functions depend on the detailed distributions of ionic species and on the radiation spectrum. Hence, these functions may vary on a very wide range of spatial and temporal scales. In this paper, we explore cooling and heating functions between 5 ≤ z ≤ 10 in simulated galaxies from the Cosmic Reionization On Computers project. We compare three functions: (1) the actual cooling and heating rates of hydrodynamic cells as a function of cell temperature, (2) the median cooling and heating functions computed using median interstellar medium (ISM) properties (median ISM), and (3) the median of the cooling and heating functions of all gas cells (instantaneous). We find that the median ISM and instantaneous approaches to finding a median cooling and heating function give identical results within the spread due to cell-to-cell variation. However, the actual cooling (heating) rates experienced by the gas at different temperatures in the simulations do not correspond to either summarized cooling (heating) functions. In other words, the thermodynamics of the gas in the simulations cannot be described by a single set of a cooling plus a heating function with a spatially constant radiation field that could be computed with common tools, such as CLOUDY.

79 ASTRONOMY AND ASTROPHYSICS↗