Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radiative transfer simulations”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Multi-Angle Polarimetry: The Once and Future King of Aerosol Remote Sensing

Although aerosols (and their interactions with clouds) are widely known to be one of the most uncertain components of the climate, they remain largely unconstrained in climate simulations. This is because global observations of all the parameters relevant to such simulations - quantity, size, shape, optical properties and chemical composition - are very difficult to simultaneously retrieve from existing remote sensing instruments. The problem can be addressed by maximizing the scene information gathered by a remote sensing instrument, by the use of (passive) multi-spectral, multi-angle and polarimetrically sensitive sensors. These observations, coupled with a radiative transfer model, can be inverted to solve for aerosol parameters. However, the choices to be made when designing an such an observing system and retrieval algorithm are complex, and a variety of approaches have been undertaken by the scientific community. I will review the various multi-spectral, multi-angle, polarimetric observation systems employed for aerosol remote sensing and their corresponding retrieval algorithms. This includes the French Polarization and Directionality of Earth Reflectance (POLDER) instrument, which has been the only such instrument successfully deployed in orbit thus far (most recently from 2004-2013), the NASA Aerosol Polarimetry Sensor (APS) on the ill-fated NASA Glory Mission (launch failure in 2011), potential or planned polarimeters on the NASA Aerosol-Cloud-Ecosystem (ACE) and Pre-Aerosol, Clouds and ocean Ecosystems (PACE) missions, and airborne prototypes from around the world.

Knobelspiesse, Kirk↗
NASA NTRSDate not supplied

The NASA Ames Aerothermodynamics Branch

The Ames Aerothermodynamics branch provides NASA and the country an integrated modeling, simulation, and test capability in the areas of aerothermodynamics, shock layer radiation, and entry system design. With roots dating back to the founding of Ames in the 1940s, the branch has contributed to every major NASA mission with an atmospheric entry component and actively partners with commercial space companies to assist with entry system development and technology transfer. This presentation provides a brief overview of the branch, its core competencies, and recent accomplishments.

Jeffrey Hill↗
DOE OSTI2024

Automation and optimization of stopping and range of ions in matter simulation runtime

Prior to every ion implantation experiment a simulation of the ion range and other relevant parameters is performed using Monte-Carlo based codes. Although increasing computing power has improved the speed of these calculations, the demands on Monte-Carlo codes are also increasing, requiring evaluation of the optimal number of simulations while ensuring accuracy within threshold bounds. We evaluate the “Stopping and Range of Ions in Matter” (SRIM) code due to its widespread usage. We show how dividing simulations into multiple parallel simulations with different random seeds can lead to calculation speedup and find lower bounds for the required number of ion traces simulated based on an exemplar system of a Ga focused ion beam and a high energy C beam as used in high linear energy transfer testing. Here our results indicate simulations can yield results within the underlying data accuracy of SRIM at 10X and 100X shorter simulation time than the SRIM default values.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Effect of soot concentration on radiative transport in a 2 m JP-8 pool fire

Hubbard et al. (2024) investigated conjugate heat transfer in a 7.9 m JP-8 pool fire, utilizing both heavily sooting and non-sooting fuel models. Their findings indicated that all models produced pool radiative heat fluxes within the measurement uncertainty for large-scale outdoor experiments. This study aims to further investigate the influence of soot concentration on radiative transport. To this end, Large Eddy Simulations (LES) of a 2 m diameter JP-8 pool fire were conducted, employing four turbulent combustion chemistry models: the Eddy Dissipation Concept (EDC) and three variants of flamelet model. Additionally, three variations of the one-equation subgrid scale turbulence kinetic energy model, and multiple mesh resolutions were utilized. The primary validation metrics were pool radiative heat flux and plume radiative fraction, with simulation data compared to experimental results published by Blanchat et al. (2010). The EDC model produced high soot concentrations near the pool surface, resulting in a significant shielding effect. In contrast, the flamelet combustion model has different kinetic limiters for soot production, leading to soot formation higher in the flame, where it can be more readily oxidized. Flamelet models incorporating both acetylene and benzene contributions to soot growth demonstrated better agreement with experimental data. To assess the sensitivity of pool radiative heat flux to soot mass concentration, we used a one-dimensional approximation over a wide range of soot concentrations and compared it to three-dimensional solutions obtained using Discrete Ordinates. Furthermore, this analysis illustrates the impacts of soot shielding and Turbulence Radiation Interactions (TRI), with TRI emission enhancement factors predicted to be in the range of 3–5 near the edges of the plume.

Computational fluid dynamics↗
DOE OSTI2026

Mesoscale Organization in Cumulus-Coupled Stratocumulus

Marine cloud systems cover a substantial portion of the world’s oceans. Most of these clouds form relatively close to the ocean surface, typically within one to two kilometers, a region referred to by meteorologists as the marine boundary layer. They are composed predominantly of liquid water, although ice particles can occur in mid- and high-latitude marine clouds during winter. In satellite imagery, these clouds appear bright against the darker ocean surface below, reflecting a large fraction of incoming sunlight back into space that would otherwise warm the ocean. Because marine boundary layer clouds cover such an extensive area of the ocean, they exert a significant influence on Earth’s overall transfer of solar energy absorbed by the surface and thermal energy emitted to space, a balance known as the planetary radiation budget. Marine boundary layer clouds are typically thin, and their formation and dissipation depend on a delicate balance between processes acting at the ocean surface below and the warm, dry air above. They are notoriously difficult to simulate accurately in weather forecast models, which often produce too few marine low clouds in midlatitudes and clouds in tropical regions that are excessively bright, meaning they reflect too much solar radiation. The marine boundary layer is frequently characterized by widespread overcast cloud cover that often transitions from a continuous, single-layer deck to more broken cloud fields toward the tropics. These transitions typically proceed through an intermediate stage in which shallow, broken clouds form beneath the overlying stratiform cloud deck. Once broken clouds develop below the overcast, they frequently self-organize into cloud clusters known as marine boundary layer convective complexes (MBLCCs), although the mechanisms governing the formation and organization of MBLCCs remain poorly understood. Accurately representing these transitions in long-range weather forecast models is essential because they influence the properties of air masses advected over the continental United States and Europe, and they become increasingly important for forecasts on seasonal and longer timescales. We employed two complementary approaches to investigate the processes controlling MBLCCs and their impact on marine cloud cover. Long-term observations from the U.S. Department of Energy’s Eastern North Atlantic (ENA) Observatory provided a unique dataset that allowed us to characterize fundamental properties of MBLCCs, including their typical size and frequency of occurrence. These observations were combined with high-resolution numerical simulations performed on supercomputers to examine the evolution of MBLCCs during cold-air outbreaks over the ENA region.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1988

Dependence of warm and cold climate depiction on climate model resolution

Climate change simulations run with two different resolutions (8 deg by 10 deg and 4 deg by 5 deg) of the NASA Goddard Institute for Space Studies model are used to investigate radiative, dynamical, and regional sensitivities to model grid size. The results show that model resolution affected two key processes in the control runs, moist convection and the nonlinear transfer of kinetic energy into the zonal mean flow. It was found that the finer resolution model has more penetrative convection but less convection overall, with the results that its temperature and wind structure were altered with respect to the coarser grid model.

Rind, David↗
NASA NTRS1995

Ocean observations with EOS/MODIS: Algorithm development and post launch studies

An investigation of the influence of stratospheric aerosol on the performance of the atmospheric correction algorithm was carried out. The results indicate how the performance of the algorithm is degraded if the stratospheric aerosol is ignored. Use of the MODIS 1380 nm band to effect a correction for stratospheric aerosols was also studied. The development of a multi-layer Monte Carlo radiative transfer code that includes polarization by molecular and aerosol scattering and wind-induced sea surface roughness has been completed. Comparison tests with an existing two-layer successive order of scattering code suggests that both codes are capable of producing top-of-atmosphere radiances with errors usually less than 0.1 percent. An initial set of simulations to study the effects of ignoring the polarization of the the ocean-atmosphere light field, in both the development of the atmospheric correction algorithm and the generation of the lookup tables used for operation of the algorithm, have been completed. An algorithm was developed that can be used to invert the radiance exiting the top and bottom of the atmosphere to yield the columnar optical properties of the atmospheric aerosol under clear sky conditions over the ocean, for aerosol optical thicknesses as large as 2. The algorithm is capable of retrievals with such large optical thicknesses because all significant orders of multiple scattering are included.

Gordon, Howard R.↗
NASA NTRS2011

SMOS/SMAP Synergy for SMAP Level 2 Soil Moisture Algorithm Evaluation

Soil Moisture Active Passive (SMAP) satellite has been proposed to provide global measurements of soil moisture and land freeze/thaw state at 10 km and 3 km resolutions, respectively. SMAP would also provide a radiometer-only soil moisture product at 40-km spatial resolution. This product and the supporting brightness temperature observations are common to both SMAP and European Space Agency's Soil Moisture and Ocean Salinity (SMOS) mission. As a result, there are opportunities for synergies between the two missions. These include exploiting the data for calibration and validation and establishing longer term L-band brightness temperature and derived soil moisture products. In this investigation we will be using SMOS brightness temperature, ancillary data, and soil moisture products to develop and evaluate a candidate SMAP L2 passive soil moisture retrieval algorithm. This work will begin with evaluations based on the SMOS product grids and ancillary data sets and transition to those that will be used by SMAP. An important step in this analysis is reprocessing the multiple incidence angle observations provided by SMOS to a global brightness temperature product that simulates the constant 40 degree incidence angle observations that SMAP will provide. The reprocessed brightness temperature data provide a basis for evaluating different SMAP algorithm alternatives. Several algorithms are being considered for the SMAP radiometer-only soil moisture retrieval. In this first phase, we utilized only the Single Channel Algorithm (SCA), which is based on the radiative transfer equation and uses the channel that is most sensitive to soil moisture (H-pol). Brightness temperature is corrected sequentially for the effects of temperature, vegetation, roughness (dynamic ancillary data sets) and soil texture (static ancillary data set). European Centre for Medium-Range Weather Forecasts (ECMWF) estimates of soil temperature for the top layer (as provided as part of the SMOS ancillary data) were used to correct for surface temperature effects and to derive microwave emissivity. ECMWF data were also used for precipitation forecasts, presence of snow, and frozen ground. Vegetation options are described below. One year of soil moisture observations from a set of four watersheds in the U.S. were used to evaluate four different retrieval methodologies: (1) SMOS soil moisture estimates (version 400), (2) SeA soil moisture estimates using the SMOS/SMAP data with SMOS estimated vegetation optical depth, which is part of the SMOS level 2 product, (3) SeA soil moisture estimates using the SMOS/SMAP data and the MODIS-based vegetation climatology data, and (4) SeA soil moisture estimates using the SMOS/SMAP data and actual MODIS observations. The use of SMOS real-world global microwave observations and the analyses described here will help in the development and selection of different land surface parameters and ancillary observations needed for the SMAP soil moisture algorithms. These investigations will greatly improve the quality and reliability of this SMAP product at launch.

Bindlish, Rajat↗
NASA NTRS1986

Simulation of jet cooling effects on Czochralski crystal growth

The effects of cooling the crystal side surface by blowing a jet of an inert gas are examined in detail for Czochralski crystal growth. A combined model of the crystal + melt, which incorporates the detailed radiation calculations, the shape of the melt-gas meniscus, predicts the growth rate and the crystal-melt interface shape, is used for this study. The convective heat transfer coefficient for the jet is estimated from the correlation available in the literature. The effect of the jet cooling on the interface shape and the pulling rate is significant. The crystal diameter as well as the interface shape tend to be more stable in the environment of the rapid cooling of the crystal by the jet. The crystal diameter or the interface shape can be easily controlled by adjusting the gas flow rate through the jet. This gives the Czochralski pulling an additional degree of freedom facilitating the control of crystal diameter and interface shape.

Srivastava, R. K.↗
DOE OSTI2024

Numerical investigation of mixing and heat transfer in a 7.9 m JP-8 pool fire

The response of objects engulfed in, and adjacent to, large-scale pool fires is of interest in accident and safety assessments. In this study, Fuego, a low-Mach turbulent reacting flow code was used to study conjugate heat transfer in a 7.9 m diameter JP-8 pool fire. Simulations were designed to replicate past experimental measurements (Blanchat et al., 2006) of incident heat flux to three cylindrical calorimeters in and around the pool fire. Two turbulent combustion models were compared directly - the eddy dissipation concept and a more recently developed unsteady flamelet model. First and second order spatial and temporal discretization schemes were also compared to assess the performance of low-dissipation numerical operators. Heat flux predictions to the transportation size calorimeter outside the fire were within experimental uncertainties. Inside the fire, experimental measurements were higher than predicted values and may have been a consequence of soot deposition and augmented participating media radiation from soot and fuel vapor. Simulation predictions improved in cases where turbulent kinetic energy and mixing were more resolved. This work, and others referenced herein, suggest that spatial resolution on the order of 0.5–1.0 cm may be required to fully resolve fluid instabilities, vortex production between fire plumes and crosswind, soot production, and fuel-air mixing. This presents a substantial computational challenge for safety assessments of engulfed objects in fully turbulent pool fires.

Conjugate heat transfer↗
DOE OSTI2024

Impact of mid- Z gas fill on dynamics and performance of shock-driven implosions at the OMEGA laser

Shock-driven implosions with 100% deuterium (D 2 ) gas fill compared to implosions with 50:50 nitrogen-deuterium (N 2 ⁢D 2 ) gas fill have been performed at the OMEGA laser facility to test the impact of the added mid-Ζ fill gas on implosion performance. Ion temperature (Τ ion ) as inferred from the width of measured DD-neutron spectra is seen to be 34%±6% higher for the N 2⁢ D 2 implosions than for the D 2 -only case, while the DD-neutron yield from the D 2 -only implosion is 7.2±0.5 times higher than from the N 2⁢ D 2 gas fill. The T ion enhancement for N 2 ⁢D 2 is observed in spite of the higher Z, which might be expected to lead to higher radiative loss, and higher shock strength for the D 2 -only versus N 2 ⁢D 2 implosions due to lower mass, and is understood in terms of increased shock heating of N compared to D, heat transfer from N to D prior to burn, and limited amount of ion-electron-equilibration-mediated additional radiative loss due to the added higher-Z material. Further, this picture is supported by interspecies equilibration timescales for these implosions, constrained by experimental observables. The one-dimensional (1D) kinetic Vlasov-Fokker-Planck code ifp and the radiation hydrodynamic simulation codes hyades (1D) and xrage [1D, two-dimensional (2D)] are brought to bear to understand the observed yield ratio. Comparing measurements and simulations, the yield loss in the N 2 ⁢D 2 implosions relative to the pure D 2 -fill implosion is determined to result from the reduced amount of D 2 in the fill (fourfold effect on yield) combined with a lower fraction of the D 2 fuel being hot enough to burn in the N 2 ⁢D 2 case. The experimental yield and T ion ratio observations are relatively well matched by the kinetic simulations, which suggest interspecies diffusion is responsible for the lower fraction of hot D 2 in the N 2 ⁢D 2 relative to the D 2 -only case. The simulated absolute yields are higher than measured; a comparison of 1D versus 2D XRAGE simulations suggest that this can be explained by dimensional effects. The hydrodynamic simulations suggest that radiative losses primarily impact the implosion edges, with ion-electron equilibration times being too long in the implosion cores. The observations of increased T ion and limited additional yield loss (on top of the fourfold expected from the difference in D content) for the N 2 ⁢D 2 versus D 2 -only fill suggest it is feasible to develop the platform for studying CNO-cycle-relevant nuclear reactions in a plasma environment.

47 OTHER INSTRUMENTATION↗
NASA NTRS1972

Thermal properties of granulated materials.

Review of the thermophysical properties of granular materials or silicates believed to simulate the lunar surface layer. Emphasis is placed on thermal conductivity data and the effects of material and environmental variables on the thermal conductivity. There are three basic mechanisms of heat transfer in particulate materials: conduction by the gas contained in the void spaces between the particles; conduction within the solid particles and across the interparticle contacts; and thermal radiation within the particles, across the void spaces between particle surfaces, and between void spaces themselves. Gas and solid conduction, thermal radiation, and the interaction between conduction and radiation are considered.

Wechsler, A. E.↗
DOE OSTI2025

Ab Initio Molecular Dynamics Study of the Reduction of Acetone by the Hydrated Electron

We have investigated the reaction dynamics of the reduction of acetone by the solvated electron in water using ab initio molecular dynamics simulations at 298 and 373 K. The rate constants derived from the simulations are consistent with experimental observations that the reaction has a low activation energy. Detailed analyses carried out to shed light on the mechanism of the reaction show that solvent reorganization plays a key role as a reaction coordinate, as expected from the Marcus electron transfer theory. Furthermore, constrained density functional theory calculations indicate that the electronic coupling is large, placing the reaction in the adiabatic limit. Indeed, the activation energies and rate constants of the simulations are in accordance with the predictions of adiabatic Marcus theory.

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

Exploring the Unusual Reactivity of the Hydrated Electron with CO 2

Many questions remain about the reactions of the hydrated electron despite decades of study. Of particular note is the fact that they do not appear to follow the Marcus theory of electron transfer reactions, a feature that has yet to be explained. To investigate these issues, we use ab initio molecular dynamics (AIMD) simulations to investigate a one of the better studied reactions, the hydrated electron reduction of CO 2 . The rate constant for the hydrated electron-CO 2 reaction complex to react to form CO 2 - is, for the first time, estimated from AIMD simulations. Results at 298 and 373 K show the rate constant is insensitive to temperature, consistent with the low measured activation energy for the reaction, and the implications of this behavior are examined. The sampling provided by the simulations yields insight into the reaction mechanism. The reaction is found to involve both solvent reorganization and changes in the carbon dioxide structure. The latter lead to significant vibrational excitation of the bending and symmetric stretch vibrations in the CO 2 - product, indicating the reaction is vibrationally nonadiabatic. The former is estimated from calculation of an approximate collective solvent coordinate and the free energy in this coordinate is determined. Furthermore, these results indicate that AIMD simulations can reasonably estimate hydrated electron reaction activation energies and provide new insight into the mechanism that can help illuminate the features of this unusual chemistry.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS2021

The Thermal Infrared Optical Depth of Mineral Dust Retrieved from Integrated CALIOP and IIR Observations

Recent studies reveal that a higher fraction of coarse mineral dust particles than that estimated by climate model simulations has been observed in the atmosphere, leading to a more significant positive (i.e., warming) longwave (LW) thermal infrared (TIR) direct radiative effect (DRE). However, the magnitude of this DRE remains highly uncertain because our understanding of the radiative properties, quantitatively represented by the optical depth of dust, especially information on the TIR, remains limited. This study presents a simple approach to retrieve the thermal infrared dust aerosol optical depth (DAODTIR) over oceans during nighttime using the observations from the Infrared Imaging Radiometer (IIR) and the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) onboard CALIPSO. For each cloud-free dust-laden profile identified by the IIR-CALIOP observation, a Lookup-Table (LUT) of the 10.6 μm IIR band brightness temperatures difference (dBT) under different DAODTIR with respect to their dust-free BTs is constructed based on the CALIOP retrieved dust vertical profile and pre-assumed dust scattering properties using a fast radiative transfer model. Then the DAODTIR is retrieved by projecting the IIR-observed dBT on the LUT. Sensitivity studies show that the DAODTIR retrieval at 10.6 μm is more susceptible to the dust particle size distribution (PSD) assumption than dust refractive indices. To estimate the uncertainty caused by PSD assumption, two DAODTIR retrieval products, one based on the dust PSD from the AERONET at Cape Verde and the other on an in situ measured PSD from the recent Fennec campaign, are provided. The retrieval uncertainty is mainly contributed by the BT difference between the observation and simulation using auxiliary atmospheric data. The climatology of the retrieval from 2013 to 2019 shows confident spatiotemporal variations of DAODTIR with the global-averaged value of 0.006 and 0.008 based on different pre-assumed dust PSDs. Climatological results agree reasonably well with two independent DAODTIR retrieval products based on the Infrared Atmospheric Sounding Interferometer (IASI) over the active dust transport regions, such as North and Tropical Atlantic (r = 0.904 and 0.819) and Indian Ocean (r = 0.832). The seasonal and interannual variation is also well-compared (r = 0.758) with AERONET coarse-mode AOD at 97 selected sites. The synergic CALIOP observation allows the retrieved DAODTIR to directly compare with the extrapolated DAODTIR from DAOD in visible (i.e., 532 nm), which helps evaluate the observational constraints on DAODTIR. This study offers a unique prospect of collocating active lidar and passive IR observations for retrieving dust DAODTIR.

aerosol↗
NASA NTRS2012

Distributed Assimilation of Satellite-based Snow Extent for Improving Simulated Streamflow in Mountainous, Dense Forests: An Example Over the DMIP2 Western Basins

Snow cover area affects snowmelt, soil moisture, evapotranspiration, and ultimately streamflow. For the Distributed Model Intercomparison Project - Phase 2 Western basins, we assimilate satellite-based fractional snow cover area (fSCA) from the Moderate Resolution Imaging Spectroradiometer, or MODIS, into the National Weather Service (NWS) SNOW-17 model. This model is coupled with the NWS Sacramento Heat Transfer (SAC-HT) model inside the National Aeronautics and Space Administration's (NASA) Land Information System. SNOW-17 computes fSCA from snow water equivalent (SWE) values using an areal depletion curve. Using a direct insertion, we assimilate fSCAs in two fully distributed ways: 1) we update the curve by attempting SWE preservation, and 2) we reconstruct SWEs using the curve. The preceding are refinements of an existing simple, conceptually-guided NWS algorithm. Satellite fSCA over dense forests inadequately accounts for below-canopy snow, degrading simulated streamflow upon assimilation during snowmelt. Accordingly, we implement a below-canopy allowance during assimilation. This simplistic allowance and direct insertion are found to be inadequate for improving calibrated results, still degrading them as mentioned above. However, for streamflow volume for the uncalibrated runs, we obtain: (1) substantial to major improvements (64-81 %) as a percentage of the control run residuals (or distance from observations), and (2) minor improvements (16-22 %) as a percentage of observed values. We highlight the need for detailed representations of canopy-snow optical radiative transfer processes in mountainous, dense forest regions if assimilation-based improvements are to be seen in calibrated runs over these areas.

Yatheendradas, Soni↗
DOE OSTI2020

Collinear dual-pulse laser optical breakdown and energy deposition

A non-equilibrium model for laser-generated plasmas is used to represent collinear dual nano-second-pulse interactions. In the specific case considered, the breakdown is initiated with an ultraviolet (UV) laser pulse at 266 nm, which pre-ionizes the gas, and it is followed by a near-infrared (NIR) pulse at 1064 nm, which deposits significant energy into the ionized mixture. The model is validated against corresponding experiments, and simulation results are interrogated to understand key features of the plasma-kernel dynamics and the post-discharge hydrodynamics. The hydrodynamics of the non-equilibrium plasma is governed by the two-temperature Navier–Stokes equations accounting for both multiphoton ionization and inverse bremsstrahlung. The interaction between the laser beam and the plasma is modeled based on the radiative transfer equation. The temporal and spatial offsets of the two pulses can generate ionization kernels with different topology and dynamics. It is shown that the UV pre-ionization pulse can tailor the plasma region, leading to a larger ionized volume than would occur for a single-pulse breakdown, and it increases the efficiency of the energy deposition for the following NIR discharge. Vorticity in the early post-discharge phase is generated via baroclinic torque from the misalignment of the radial gradient of density (sudden gas expansion) and the strong pressure gradient (initiated by energy deposition from the laser beam). It is predicted that the collinear dual pulse can be designed to specify the plasma kernel size, shape, and maximum temperature in the breakdown phase, and the initial post-breakdown vorticity and plasma-core decay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

iharm3D: Vectorized General Relativistic Magnetohydrodynamics

iharm3D is an open-source C code for simulating black hole accretion systems in arbitrary stationary spacetimes using ideal general-relativistic magnetohydrodynamics (GRMHD). It is an implementation of the HARM (“High Accuracy Relativistic Magnetohydrodynamics”) algorithm outlined in Gammie et al. (2003) with updates as outlined in McKinney & Gammie (2004) and Noble et al. (2006). The code is most directly derived from Ryan et al. (2015) but with radiative transfer portions removed. HARM is a conservative finite-volume scheme for solving the equations of ideal GRMHD, a hyperbolic system of partial differential equations, on a logically Cartesian mesh in arbitrary coordinates.

79 ASTRONOMY AND ASTROPHYSICS↗