Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D equilibrium”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Erosion and impurity transport for the edge localized mode suppression window in KSTAR

A 3D analysis of plasma wall interactions and global impurity transport for the edge localized mode suppression window in KSTAR during H-Mode discharges has shown that carbon erosion at the divertor plates is a strong function of the resonant magnetic perturbation (n = 1) coil current and relative phasing. The Generalized Perturbation Equilibrium Code was used to determine a realistic initial perturbed plasma equilibrium, and EMC3-EIRENE was used to calculate the resulting scrape-off layer plasma used in this study as a fixed background for the ERO2.0 plasma–material interaction model. The resulting transport leads to deposition of impurities along the targets positioned at the high-field side of the device. An attempt at calculating the resulting effective charge state has demonstrated a similar dependence on the perturbation coil current and has been able to determine a window for the experimentally observed values of Zeff by including contributions of all ionized carbon charge states and deuterium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Neoclassical transport of impurities in tokamaks with non-axisymmetric perturbations

The effect of resonant magnetic perturbations (RMPs) on neoclassical transport of impurities is calculated in DIII-D and NSTX. Here, the neoclassical fluxes are evaluated using the NEO code with nonlinear one-fluid nonaxisymmetric equilibrium calculated using M3D-C1. Neoclassical fluxes of impurities show significant changes with RMPs if the impurity resides in low-collisionality regime, but are weakly affected by RMPs in the Pfirsch–Schlüter (P–S) regime. Charge number (Z) of the impurity affects the collisionality of impurity species, which then determines the effect of 3D fields on neoclassical impurity transport. This suggests that RMPs can be possibly used for exhaust of low-Z impurities in these facilities, but have little effect for impurities with charge number greater than 10 or so. Additionally, it is shown that it is the change in the convective fluxes which is a main contributor in modifying the neoclassical impurity fluxes in the presence of RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fast ion stabilization of tilt in large radius FRCs

The field reversed configuration (FRC) has been a curious case in plasma physics research in that early MHD analysis suggested FRCs should be grossly unstable, while experimental results contradicted that prediction. Later, this theory was able to resolve this contradiction by understanding that finite Larmor radius effects largely negated the MHD predictions. Similarly, previous theoretical studies of beam driven FRCs predicted that such system would be unstable to beam driven modes while, again, experimental results indicated the contradiction. In this paper, we reconcile the theoretical understanding of beam driven modes with experimental observations of stability in these systems. By self-consistently capturing fast ion generation from neutral beam injection and its impact on the plasma equilibrium, we show that low amplitude perturbations in the magnetic field, driven by betatron particles, modify the precession frequencies of the betatron particles such that the drive for compressional Alfvén waves in the thermal plasma is reduced. Finally, we are able to demonstrate, for the first time, stable beam driven FRC evolution at high S*/E in 3D kinetic simulations.

Betatrons↗
DOE OSTI2025

Forced 3D Reconnection in an Exponentially Separating Magnetic Field

We present a solvable scenario for 3D reconnection in a sheared magnetic field. We consider a localized external force that is applied slowly to a flux tube and then maintained, generating an Alfvénic perturbation that spreads along the field lines. Separation of the sheared field lines reduces the scale of the perturbation across the field, enhancing magnetic diffusion. For a fusion-motivated equilibrium with exponential field-line separation, we find a reconnection timescale proportional to 𝒮/ ln⁡ 𝒮 under magnetohydrodynamics (MHD) and to 𝒮 1/3 for semicollisional electron-only reconnection, where 𝒮 is the Lundquist number of the perturbed flux tube. We generalize these results to arbitrary magnetic geometries, showing that the semicollisional case is geometry independent. Interestingly, we find that slower field-line separation yields an increased reconnection rate in MHD.

Magnetic confinement fusion↗
DOE OSTI2020

A Monte Carlo Approach to Approximating the Effects of Pore Geometry on the Phase Behavior of Soil Freezing

Abstract Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Photoinduced Dirac semimetal in ZrTe5

Abstract Novel phases of matter with unique properties that emerge from quantum and topological protection present an important thrust of modern research. Of particular interest is to engineer these phases on demand using ultrafast external stimuli, such as photoexcitation, which offers prospects of their integration into future devices compatible with optical communication and information technology. Here, we use MeV Ultrafast Electron Diffraction (UED) to show how a transient three-dimensional (3D) Dirac semimetal state can be induced by a femtosecond laser pulse in a topological insulator ZrTe 5 . We observe marked changes in Bragg diffraction, which are characteristic of bond distortions in the photoinduced state. Using the atomic positions refined from the UED, we perform density functional theory (DFT) analysis of the electronic band structure. Our results reveal that the equilibrium state of ZrTe 5 is a topological insulator with a small band gap of ~ 25 meV, consistent with angle-resolved photoemission (ARPES) experiments. However, the gap is closed in the presence of strong spin-orbit coupling (SOC) in the photoinduced transient state, where massless Dirac fermions emerge in the chiral band structure. The time scale of the relaxation dynamics to the transient Dirac semimetal state is remarkably long, τ ~ 160 ps, which is two orders of magnitude longer than the conventional phonon-driven structural relaxation. The long relaxation is consistent with the vanishing density of states in Dirac spectrum and slow spin-repolarization of the SOC-controlled band structure accompanying the emergence of Dirac fermions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Modeling of HTTF test PG-26 using RELAP5-3D and SAM

The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled nuclear reactors (HTGR) with prismatic fuel and reflector blocks [1]. Several tests have been completed at the HTTF including depressurized conduction cooldown (DCC) and pressurized conduction cooldown (PCC) transients. This summary reports on the analysis of test PG-26 using the INL system code REALP5-3D [2] as well as the ANL system code SAM [3]. Test PG-26 is a progression of the Double Ended Inlet-Outlet Crossover Duct Break transient that is referred to as a DCC [4]. Core initial conditions (i.e., before the DCC started) have been met using low power (<100 kW) and two of ten available electric heater banks. The DCC transient was initiated during the 50th hour of the test. The break valves were opened, and hot helium from the core and cold helium from the reactor cavity simulation tank (RCST) started mixing. The gases flowed in a countercurrent fashion, where the top half of the hot duct contained hot helium that flowed in one direction and cold helium that flowed in the other direction in the bottom half of the duct. After the pressure and density reached equilibrium, the event entered a diffusion mode. The onset of a reverse natural circulation was not observed during the DCC period of the test.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Structure and Flow-Viscosity of Filled-Polymer-Based 3D Printing Ink: Exploration through Coarse-Grained Molecular Dynamics

The addition of nanofiller particles to a polymer matrix has long been known to enhance or modify the composite’s mechanical and rheological properties. However, quantitatively capturing such changes with molecular level simulations remains computationally challenging. Toward that goal, we performed coarse-grained molecular dynamics of a nanocomposite system at a fixed (25 vol %) filler loading under nonspecific, weak polymer–filler interactions representative of a broad class of technologically important materials. We report several interesting results, including: (1) the equilibrium chain-configuration remains Gaussian-like as in an unfilled melt; (2) smaller filler particles display a stronger tendency to cluster; (3) larger fillers act as plasticizers by reducing the entanglement density and accelerating the chain mobility; and (4) fillers enhance the tensile response modulus, with the effect being stronger for larger particles. We also simulate cluster breakup, yielding, and elongational flow under an applied time-linear tensile strain and study the flow viscosity as a function of filler-size and chain-length.

Materials science↗
DOE OSTI2026

Validation of an integrated modeling framework for investigating 3D plasma responses in tokamak plasmas

As contemporary experimental tokamaks are pushed toward reactor-relevant operation, they provide essential testbeds for demonstrating and exploring ELM control strategies for deployment in future fusion pilot plants. Accurate predictions of full plasma responses are essential to guide and optimize these demonstrations. This paper introduces and validates an integrated modeling framework over a historical range of DIII-D operational space. The integrated modeling framework uses only scalar plasma parameters and optional reference boundary to self-consistently and flexibly scan through tokamak operational space and estimate the corresponding ELM-suppression relevant plasma response over 3D coil phase space. The framework generates tightly converged equilibria that satisfy a target set of plasma parameters (I p , β N , l i ), with kinetic profiles constrained by an EPED(NN)-computed pedestal and empirical core model. The plasma response of these modeled equilibria is calculated with GPEC. Validated against 55 distinct DIII-D equilibria, the framework consistently reproduces experimentally constrained equilibrium pressure, q, and other representative profiles using the time-varying information from only evolving scalar plasma parameters; plasma response validation is performed for the n=3 perturbation. In particular, the plasma response of a DIII-D discharge scanning q 95 in search of ELM suppression windows is quantitatively reproduced. By accurately reproducing experimental equilibria and plasma responses across wide parameter variation, the framework supports potential for synthetic parameter scans in key operational, stability, and plasma response dimensions to investigate RMP ELM-suppression experiments and inform predictive RMP scenario optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Verification of a quasi-linear model for gyrokinetic turbulent transport

The verification and calibration of a new quasi-linear transport model with a large database of gyrokinetic turbulence simulations is presented in this paper. In a previous paper (Staebler et al 2020 Plasma Phys. Control. Fusion 63 015013), a model for the saturated spectrum of electric potential fluctuations was developed based on the properties of the non-linear 3D spectrum. In this paper, a modification to the overall multiplicative factor in this model is found to be necessary to improve the fit to scans of the temperature and density gradients. The error in the fit of the quasi-linear fluxes of electron and ion energy fluxes is significantly better than for previous saturation models. The spectral shift model for the impact of equilibrium E × B velocity shear (Staebler et al 2013 Phys. Rev. Lett. 110 055003) and the zonal flow mixing model for electron-scale turbulence (Staebler et al 2016 Phys. Plasmas 23 062518) are both revised to be compatible with this new model. Furthermore, the models for the loss of bounce averaging and electron collisions in the TGLF reduced linear equations (Staebler et al 2005 Phys. Plasmas 12 102508) are also changed to improve the linear eigenmodes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A platform for non-local thermodynamic equilibrium atomic physics research using the buried layer target approach at a kJ level laser facility

Here we present a design and first use of a kJ level laser facility for research of non-local thermodynamic equilibrium atomic physics using the buried layer target method. The target design included a metal layer buried inside a plastic tamper with thicknesses tailored to the expected laser intensities. The target was illuminated from each side by two laser beams with intensities of 0.5–5 × 10 14 W/cm 2 . The advanced diagnostic suite included static and time-resolved imagers and spectrometers with various spectral resolutions. A 3D printed dual elliptically curved spectrometer is presented, and its results are compared to a traditional crystal spectrometer. Experimental results and radiation hydrodynamic simulations demonstrate that the target achieved the desired thermodynamic conditions of n e ≈ 10 21 –10 22 cm -3 and T e ≈ 1–2 keV.

3D printing↗
DOE OSTI2023

Status of Multiphysics Modeling of Critical DPCs in a Geological Repository

A summary of previous work focused on the first two scenarios is presented as an established workflow. These canisters could reach an equilibrium at a low quasi–steady-state power level that could be sustained for thousands of years. The majority of this report describes an effort to expand the initial methodology to include the scenario in which a DPC has enough excess reactivity to boil the water inside the canister. This scenario involves a variety of tools and modeling approaches to capture the two-phase behavior, including RELAP5-3D, STAR-CCM+, and TRACE. Efforts using the high-fidelity computational fluid dynamics results to verify the accuracy of RELAP5-3D and TRACE for boiling DPCs are also presented. A new modeling strategy using TRACE appears to be the most promising approach for expanding the analysis of critical DPCs beyond the subcooled scenarios. Finally, this report describes ongoing work investigating the fourth scenario, in which a critical DPC is sealed in saturated bentonite. The primary goal in this research is to understand the power output and associated pressure increase to inform geomechanics simulations that could predict the impact on a geological repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Gradient-based optimization of 3D MHD equilibria

Using recently developed adjoint methods for computing the shape derivatives of functions that depend on magnetohydrodynamic (MHD) equilibria (Antonsen et al. , J. Plasma Phys. , vol. 85, issue 2, 2019; Paul et al. , J. Plasma Phys. , vol. 86, issue 1, 2020), we present the first example of analytic gradient-based optimization of fixed-boundary stellarator equilibria. We take advantage of gradient information to optimize figures of merit of relevance for stellarator design, including the rotational transform, magnetic well and quasi-symmetry near the axis. With the application of the adjoint method, we reduce the number of equilibrium evaluations by the dimension of the optimization space ( ${\sim }50\text {--}500$ ) in comparison with a finite-difference gradient-based method. We discuss regularization objectives of relevance for fixed-boundary optimization, including a novel method that prevents self-intersection of the plasma boundary. We present several optimized equilibria, including a vacuum field with very low magnetic shear throughout the volume.

Physics↗
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Simultaneous x-ray and neutron 4D tomographic study of drying-driven hydro-mechanical behavior of cement-based materials at moderate temperatures

The drying of cement-based materials is intimately related to their durability, which has significant economic, social and environmental repercussions. The evolution of the saturation of the pore network and the associated drying shrinkage are in fact leading causes of cracking and of the ensuing penetration of aggressive chemicals. This process is highly heterogeneous, due to the thermo-hydric spatial gradients developing in the material from the exposed surfaces to its core and because of local effects, driven by the intrinsically heterogeneous micro-structure (e.g., by the distribution of pores and aggregates). It follows that macroscopic, sample-scale measurements cannot fully disclose the complexity of the underlying processes. In the last few decades, significant advances in full-field techniques have allowed an unprecedented insight into these local processes. For cement-based materials, x-ray and neutron tomography lend themselves as ideal, and highly complementary, tools for the study of their thermo-hydro-mechanical behavior. Notably, the high sensitivity to density variations of x-ray imaging gives access to the developments of fractures, in 4D (3D + time). On the other hand, neutron tomography allows the study of the evolution of the moisture field in 4D, thanks to its high hydrogen sensitivity. The combination of these two techniques provides a unique insight in thermo-hydro-mechanical couplings, e.g., the effect of cracks on the water content field. This contribution presents novel 5D datasets (3D tomographies along time, plus truly simultaneous x-ray and Neutron rapid acquisitions) in-operando of a cement paste and of a concrete sample heated at moderate temperatures (up to 140{sup ∘}C). The analysis of this 5D data-set (once aligned in time and across modalities) allows for example a quantification of the 4D moisture profiles which were found to predict an overall water loss at hydric equilibrium coherent with the corresponding analytical analysis. In the cement paste sample, the x-ray dataset captures the evolution of an extensive cracking network, opening and propagation toward the core of the sample. A novel analysis procedure is here proposed which allows the extraction of these fractures and the analysis of their interplay with local drying as captured through neutron imaging. This for example reveals the depth of penetration of drying in the vicinity of the fractures along time, which is essential for the assessment and calibration of hydro-mechanical coupled models.

36 MATERIALS SCIENCE↗
DOE OSTI2024

MHD stability of spherical tokamak equilibria with non-monotonic q-profiles

We use the 3D magnetohydrodynamic (MHD) code M3D-C1 [Jardin et al., Comput. Sci. Discovery 5, 014002 (2012)] to examine the MHD stability and subsequent evolution of NSTX shot 129169. This discharge had a period with a non-monotonic safety factor profile, q (reversed shear), which was terminated by a MHD event that abruptly lowered the central safety factor, q0, and greatly reduced the peakedness of the pressure profile. We show that the equilibrium just before the MHD event occurred was linearly unstable to many pressure-driven infernal modes. Modes with toroidal mode number n≥3 all had rational surfaces very close to the minimum value of q. However, a non-resonant pressure-driven (1, 1) mode was also present, and this dominated the nonlinear evolution. The final state in the simulation, after the MHD activity subsided, had a reduced and flattened pressure profile and a nearly monotonic q-profile, in qualitative agreement with experimental results. The initial state was also unstable to the resistive interchange criteria in the reversed-shear region, but the final state was stable everywhere. The “double tearing mode” (DTM) does not appear to play a role in the MHD activity of this discharge. In Appendix A, we show that in a torus, the DTM is strongly stabilized by pressure, but it is destabilized in cylindrical geometry (which has been the most extensively analyzed in the literature).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Deconvoluting sources of variability in aerosol jet printing using light scattering measurements

Aerosol jet printing (AJP) is a digital additive manufacturing technique for hybrid and conformal electronics, where its contactless deposition readily enables patterning over 3D surface topography. However, the complexity of aerosol transport physics makes deposition rate sensitive to process variability, hindering widespread industry adoption. Light scattering measurements have recently been demonstrated as a capable tool for measuring deposition rate in real time, enabling closed-loop control and in-situ qualification frameworks. These inline optical measurements present further opportunity as a diagnostic tool to understand the effects of secondary process parameters affecting vapor–liquid equilibrium and heat and mass transport during deposition. Downstream of the optical measurement cell, changing the sheath gas flow rate, increasing temperature via an inline heater, and adding solvent vapor via a sheath gas bubbler were observed to alter drying physics during sheath-collimation, and thereby impaction efficiency. Further upstream, the introduction of solvent vapor via a carrier gas bubbler and liquid build-up in the printhead were seen to affect the quantitative relationship between the light scattering data and the true deposition rate. Using this information, tightened controls of meaningful secondary parameters were implemented to improve the batch-to-batch consistency for printing a dielectric ink. In addition to advancing AJP process reliability for production, this work demonstrates the capability of inline optical measurements to provide insight into process physics and deconvolute competing mechanisms that underly variability.

additive manufacturing↗
DOE OSTI2020

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↗