Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow dynamics”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Multiscale Approach to Simulate Non‐Isothermal Multiphase Flow in Deformable Porous Materials

Coupled thermal, hydraulic, and mechanical processes in porous materials play important roles in several energy and environmental technologies. The Darcy-Brinkman-Biot (DBB) framework has proven effective in modeling multiphase fluid flow in deformable porous solids across both pore and Darcy scales, including in systems where fractures coexist with a porous matrix. In this study, we extend the DBB framework, originally designed for isothermal conditions, to address non-isothermal problems by incorporating an energy conservation equation. The resulting solver, hybridBiotThermalInterFoam, enables simulations of coupled multiphase fluid flow, heat transfer, and solid deformation in hybrid-scale systems containing both solid-free regions and ductile porous domains. The new solver is validated through comparisons with analytical solutions and, also, against established heat transfer solvers chtMultiRegionFoam and compressibleInterFoam. Further, a series of 2D and 3D case studies, including two-phase heat transfer in solid-free, static, or deformable porous media, highlights the solver's capacity to simulate complex flow dynamics and heat transport in systems involving high mobility ratios, viscous fingering, and fracture propagation. Our results establish the feasibility of incorporating thermal effects in simulations of a wide variety of energy geotechnics and environmental applications, including enhanced hydrocarbon recovery, soil remediation, and enhanced geothermal energy systems.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2020

Aerotaxy: gas-phase epitaxy of quasi 1D nanostructures

Cost- and resource-efficient growth is necessary for many applications of semiconductor nanowires. We here present the design, operational details and theory behind Aerotaxy, a scalable alternative technology for producing quality crystalline nanowires at a remarkably high growth rate and throughput. Using size-controlled Au seed particles and organometallic precursors, Aerotaxy can produce nanowires with perfect crystallinity and controllable dimensions, and the method is suitable to meet industrial production requirements. Here, we explain why Aerotaxy is an efficient method for fabricating semiconductor nanowires and explain the technical aspects of our custom-built Aerotaxy system. Investigations using SEM (Scanning Electron Microscope), TEM (Transmission Electron Microscope)and other characterization methods are used to support the claim that Aerotaxy is indeed a scalable method capable of producing nanowires with reproducible properties. We have investigated both binary and ternary III–V semiconductor material systems like GaAs and GaAsP. In addition, common aspects of Aerotaxy nanowires deduced from experimental observations are used to validate the Aerotaxy growth model, based on a computational flow dynamics (CFD) approach. We compare the experimental results with the model behaviour to better understand Aerotaxy growth.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Effect of Wettability Gradient on the Scale Formation in Falling Film Flow

Mitigation of scale formation and performance degradation remains a vital challenge for falling film evaporators in various industries. In this work, an experimental study of falling film flow on a horizontal tube is conducted to investigate the effects of wettability gradients on thermal, hydraulic, and fouling behavior. It is revealed that certain hydrophobic coating patterns, such as strip, ring, and grid patterns, lead to unwetted heat transfer area, which results in decreased heat transfer compared to fully wetted plain tube. By adjusting the geometry and position of the wettability gradient, the hybrid coating demonstrates improved heat transfer performance. Based on the characteristics of horizontal tube falling film flow, impinging jet, thin film flow, and liquid retention at the tube bottom, a hybrid coating pattern is developed to improve surface wetting and mitigate the scaling coverage. It is revealed that scale deposition is regulated by wettability gradient. Crystals tend to be dense and compact in hydrophilic areas, while they appear scattered or even absent in hydrophobic regions, depending on the dimension of the hydrophobic area. While at the hydrophilic/hydrophobic boundary, a noticeable scale thickness step is observed, which raises the potential for self-cleaning. Furthermore, the balance of minimization of scaling layer coverage and maximization of wetting area requires an optimal design in coating dimensions, for which a systemic study of both flow dynamics and fouling characteristics on the falling film is necessary in the future.

42 ENGINEERING↗
DOE OSTI2020

Evaluating the Influence of Hemorheological Parameters on Circulating Tumor Cell Trajectory and Simulation Time

Extravasation of circulating tumor cells (CTCs) occurs primarily in the microvasculature, where flow and cell interactions significantly affect the blood rheology. Capturing cell trajectory at this scale requires the coupling of several interaction models, leading to increased computational cost that scales as more cells are added or the domain size is increased. In this work, we focus on micro-scale vessels and study the influence of certain hemorheological factors, including the presence of red blood cell aggregation, hematocrit level, microvessel size, and shear rate, on the trajectory of a circulating tumor cell. We determine which of the aforementioned factors significantly affect CTC motion and identify those which can potentially be disregarded, thus reducing simulation time. We measure the effect of these elements by studying the radial CTC movement and runtime at various combinations of these hemorheological parameters. To accurately capture blood flow dynamics and single cell movement, we perform high-fidelity hemodynamic simulations at a sub-micron resolution using our in-house fluid dynamics solver, HARVEY. We find that increasing hematocrit increases the likelihood of tumor cell margination, which is exacerbated by the presence of red blood cell aggregation. As microvessel diameter increases, there is no major CTC movement towards the wall; however, including aggregation causes the CTC to marginate quicker as the vessel size increases. Finally, as the shear rate is increased, the presence of aggregation has a diminished effect on tumor cell margination.

Roychowdhury, Sayan↗
DOE OSTI2020

Application of a Heat Integrated Post-combustion CO 2 Capture System with Hitachi Advanced Solvent into Existing Coal-Fired Power Plant (Final Technical Report)

The goal of this final project report is to comprehensively summarize the work conducted on project DE-FE0007395. In accordance with the Project Management Plan (PMP), Revision F dated 5/10/2019, and Statement of Project Objectives (SOPO) within, the University of Kentucky (UK) Center for Applied Energy Research (CAER) (Recipient) has successfully demonstrated a unique, versatile CO 2 capture system (CCS) using a heat integrated process combined with two-stage stripping for process intensification, heat recovery and demineralized (DM) water generation. This project involved the design, fabrication, installation, testing of and data analysis from the UK CAER 0.7 MWe small pilot scale CO2 capture process installed at Kentucky Utilities (KU) E.W. Brown Generating Station in Harrodsburg, KY. During each of the four project Budget Periods (BPs), UK CAER met all project deliverables, all project milestones, with National Energy Technology Laboratory (NETL) approved adjustments made to the campaign long-term hours during BP4. The CCS was constructed in modular skids. Two solvent campaigns were initially conducted; the first with a 30 wt% monoethanolamine (MEA) as a baseline, and the second with the Hitachi H3-1 advanced solvent. Additional tests were performed with two advanced solvents including CAER and Proprietary Solvent C. Short-period testing with a higher concentration of 40 wt% MEA was conducted to evaluate the potential saving with high alkalinity. From the various solvent campaigns, unique aspects of the UK CAER CCS technology, as well as its flexibility and versatility were experimentally validated and demonstrated. With respect to solvent evaluation efforts in identifying candidates with significant operational and capital cost savings potential, performance of solvents were evaluated to determine the energy requirements for regeneration; environmental impacts from secondary emissions and degradation products; degradation rates, solvent make-up rates and stability. The assessments were done from parametric tests that determined optimum operating conditions for the individual solvents to maximize process efficiency and minimize the parasitic load of the power plant, and from long term campaigns (1000 hours for 30 wt% MEA and 1000 hours for H3-1) which collectively informed the techno-economic analyses (TEA) of the process. The long term campaigns included corrosion studies which used three types of metal coupons in different sections of the process: (absorber, primary stripper, lean carbon-loaded and rich carbon loaded flow streams in process) to mimic heat and flow dynamics process equipment were exposed to. The estimated corrosion rates were used to elucidate corrosion mechanisms and to further guide process material selection for potential capital cost savings.

gas separation↗
DOE OSTI2023

Laser heating and evaporation of a single droplet

The laser technology is being abundantly studied for controlled energy deposition for a range of applications in aerodynamic flow control, material processing, ignition, and combustion. The absorption of laser radiation by liquid droplets affects further propagation of laser in the atmosphere and causes bleaching of suspended droplets while the ignition and combustion characteristics in combustors are influenced by the evaporation rate of the sprayed fuel. In this work, we present a multi-dimensional mathematical model built on OpenFOAM for laser heating and evaporation of a single droplet in the diffusion dominated regime taking into account absorption of the laser radiation, evaporation process and vapor flow dynamics. The developed solver is validated against available experimental and numerical data for the ethanol and water droplet heating and evaporation. For continuous heating the peak temperature is established by the balance of cooling, evaporation and heating and results in high temperature for larger droplets. It has been shown that for heating by a single laser pulse the maximum temperature of droplets depends only on the peak intensity of the laser radiation. Furthermore, for the peak irradiance close to the transition to the boiling regime, temporal dynamics of the droplet temperature is independent of the droplet size. With proper normalization of time, the dynamics of the droplet shrinkage and cooling is shown to be independent of droplet sizes and peak laser intensities. The influence of cooling and evaporation processes on droplet heating was found to be controlled by the pulse repetition rate for repeated pulse operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Scaling hyporheic nitrogen cycling in large river alluvial aquifers

Research Element 1: Denitrification is a difficult flux to measure when fluxes are low and when the ecosystem is difficult to sample. These attributes describe Nyack floodplain, a large gravel-bedded aquifer adjacent to the Middle Fork Flathead River. Given that denitrification is simply respiration with nitrate as a terminal electron acceptor, we reasoned that N 2 production should negatively covary with oxygen respiration so long that oxygen is available. We sampled 16 wells 8 times throughout a year and analyzed these water samples on a membrane-inlet mass spectrometer for N 2 :Ar and O 2 . We fit multilevel models to N 2 :Ar as a function of dissolved oxygen. Within wells, there was no relationship between O 2 and N 2 :Ar. Among wells N 2 :Ar increased as O 2 declined. Based on the slope of this relationship, ~ -0.07, and estimated oxygen respiration in Nyack derived from prior studies, a possible range for denitrification is 0-8 mmol m -2 d -1 corresponding to a turnover rate for nitrate of 0-0.2/d. This estimated turnover rate is in the lower 20% of rates measured in the LINX 2 study, yet shows that even in a low nitrate (0.005 mmol/L) aquifer denitrification can be a substantial sink for dissolved N. Research Element 3 and 4: We developed a new theoretical approach to describe the relationship between water inflow, water outflow, and water age in porous media control volumes under dynamic flow conditions. Our approach had three distinct steps. First we described the application of well-established equations borrowed from Chemical Engineering to describe water age distributions in the hyporheic zone under steady state conditions. We published the associated results in the scientific journal PLOS ONE. Second, we developed a new theory ("flow spans") to describe how predicted steady-state water age distributions will respond to time varying imbalances in hyporheic discharge and recharge. Finally, we developed a proof-of-concept simulation model based on the new theory and used the model to examine the effects of flow regulation on a simulated river with characteristics representative of western U.S. montane rivers with expansive hyporheic zones.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Journal paper

Laser technology is being widely studied for controlled energy deposition for a range of applications, including flow control, ignition, combustion, and diagnostics. The absorption and scattering of laser radiation by liquid droplets in aerosols affects propagation of the laser beam in the atmosphere, while the ignition and combustion characteristics in combustion chambers are influenced by the evaporation rate of the sprayed fuel. In this work, we present a mathematical model built on OpenFOAM for laser heating and evaporation of a single droplet in the diffusion-dominated regime taking into account absorption of the laser radiation, evaporation process, and vapor flow dynamics. The developed solver is validated against available experimental and numerical data for heating and evaporation of ethanol and water droplets. The two main regimes—continuous and pulsed laser heating—are explored. For continuous laser heating, the peak temperature is higher for larger droplets. For pulsed laser heating, when the peak irradiance is close to transition to the boiling regime, the temporal dynamics of the droplet temperature does not depend on the droplet size. With the empirical normalization of time, the dynamics of the droplet shrinkage and cooling are found to be independent of droplet sizes and peak laser intensities.

42 ENGINEERING↗
DOE OSTI2024

Testing of Helium-Cooled Metal Molybdenum Disk Target

NorthStar Medical Radioisotopes LLC is planning to produce an important medical radioisotope, molybdenum-99 (Mo-99), through photonuclear reaction on molybdenum-100 (Mo-100). In this approach, molybdenum metal will be bombarded with a 40-MeV electron beam. Because enriched Mo-100 is expensive, it is desired to use as much beam power as possible to achieve maximum production yield and minimize target mass. This objective leads to very high beam power density (heat deposition in the target), which sets challenging requirements for cooling. Together with scientists at Los Alamos National Laboratory, a team at Argonne National Laboratory has developed and demonstrated a cooling approach using pressurized helium, which allows for efficient heat removal. One of the main challenges in this approach is the management of the heat load on the target window. The target window separates the high-pressure helium inside the target from the vacuum in the beamline, so it is constantly under stress from differential pressure. Also, the window is cooled only by the helium gas flowing on one side, making the window cooling even more challenging. High heat deposition in the target disks also imposes a strict requirement on performance of the helium cooling system and thickness of the target disks. The target disks are produced from metal powder via a press-and-sinter process. The resulting disks do not possess tensile strength as high as solid molybdenum and might not survive the vibration from the high-velocity helium coolant and high thermal stress from beam heating. An Argonne team of scientists performed a series of tests at Argonne’s Low Energy Accelerator Facility (LEAF) [3-7]. This report describes two series of tests for scale down production target designs that utilize full-scale 29 millimeters diameter, 0.75 mm thick press-and-sintered disks. We performed two thermal tests with different beam parameters and configurations of the disk laminations. We compared the results of the window temperature measurements and cooling system parameters obtained in the experiments with those predicted by analytical calculations and Computation Flow Dynamic (CFD) simulations. Results of the experiments and calculations are presented below.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Filament Formation and Melt Spinning of Coal-Based Mesophase Pitch for Carbon Fiber Production

High-performance carbon fibers excel as high specific strength and modulus materials, and are utilized in composites applications ranging across aerospace, automotive, energy, infrastructure, and sports equipment sectors. Compared to PAN-based carbon fiber, mesophase pitch-based carbon fiber has lower tensile strength but provides a higher modulus, higher thermal conductivity, and a lower cost potential. However, challenges in stable, continuous melt spinning processing are a serious limiting factor. Unlike typical melt spinning of long, linear chain polymers, mesophase pitch is a comprised of relatively shorter polycyclic aromatic hydrocarbons which form a liquid crystal. In its ‘green’ or ‘as-spun’ state, the fibers are very fragile. Moreover, its temperature of processing is quite high, often approaching 400ºC. Complex flow dynamics combined with short length and time scale heat transfer render its melt spinning a formidable processing challenge. Improved understanding of the root causes for nascent filament breakage and spinning instabilities are needed. This work aims to determine fundamental phenomena that govern mesophase pitch filament formation by comparison of required draw force and uninterrupted spinning minutes with applied draw down ratio. The effect of perturbations on spinning stability are analyzed using capillary rheometry, microscopy of elongating filaments, and filtration.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Iodine Capture Studies of Copper- and Bismuth-Based Sorbents

The release of radioiodine, one of several radionuclides of concern when recycling used nuclear fuel (UNF), is an important consideration in the fuel cycle. In this study, two sorbent materials, Cu 0 -polyacrylonitrile (PAN) and Bi 0 -PAN, were tested as solid sorbent candidates for iodine capture. Experiments using a thin bed of sorbent material were exposed to vaporized iodine for over 300 hours (~2 weeks) under varied conditions in a dynamic flow environment. The overall performance was monitored in real-time using thermogravimetric analysis, and the materials were subsequently characterized for surface and bulk analysis using scanning electron microscopy – energy-dispersive spectroscopy (SEM-EDS) and powder x-ray diffraction (pXRD), respectively. Iodine (in the form of I 2 ) is expected to be released primarily in the dissolver off-gas (DOG) stream; therefore, this study demonstrates the effects of elemental iodine (I 2 ), water vapor (H 2 O), and nitrogen dioxide (NO 2 ). When exposed to ‘ideal’ conditions (in which I 2 is carried by dry air), Cu 0 -PAN and Bi 0 -PAN behave differently, with TGA analysis indicating that Cu 0 sorption performance is higher than that of Bi 0 , as evidenced by a larger mass change. Under ‘harsh’ conditions—such as a gaseous feed containing I 2 , H 2 O, and NO 2 vapors carried by air,—iodine capture performance for both Cu 0 - and Bi 0 -PAN are affected. SEM-EDS and pXRD analysis of these materials is discussed herein.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

The Hydropower Game: An Interactive Learning Tool for the Future Hydropower Workforce

This report documents a Seedling project supported by the U.S. Department of Energy (DOE) Water Power Technologies Office (WPTO), conducted by Argonne National Laboratory (ANL) and the National Laboratory of the Rockies (NLR). The Seedling and Sapling Program provides small grants over short duration to enable early-stage research at national laboratories. This funding mechanism is intended to cultivate innovative ideas and expand research road maps in hydropower and marine energy technologies. Under this Seedling award, the project team developed an educational outreach tool, or “serious game”, built in Python and powered by the Pygame library, aimed at teaching fundamental hydropower concepts through an engaging, interactive learning experience. The game features multiple lessons covering several hydropower topics while allowing for expansion and customization in the instance of future funding availability: • Hydropower Plant Types: Players explore and compare the mechanics and applications of run-of-river, dam hydropower, and pumped-storage hydropower. • Flow Dynamics and Power Generation: Interactive tasks demonstrate how release rate and hydraulic head combine to determine power output. • Grid Operations and Load Following: Simulations illustrate how hydropower facilities respond to fluctuating electricity demand to balance the grid. • Market Integration: Levels scaffold understanding of how hydropower interfaces with the broader energy market, including operational and economic considerations. The tool was showcased at several science, technology, engineering, and mathematics (STEM) outreach events, where it was demoed to students, educators, and the general public. These events provided valuable opportunities to collect feedback on the gameplay, storyline, and educational objectives from a wide audience. The game was presented at STEMapalooza, Introduce a Girl to Engineering Day, STEMCON, and Clean Currents 2025. Insights gathered during these demonstrations informed refinements to the game’s user experience and strengthened its effectiveness as an educational tool for teaching hydropower concepts. This report outlines the game’s design philosophy, educational objectives, technical implementation, user experience insights, and potential for broader deployment within educational and workforce development contexts. It emphasizes how gamified learning can demystify complex hydropower science and inspire interest in water-power technologies. The project’s success demonstrates the value of Seedling funding in fostering creative, lowcost educational tools that support DOE’s mission to advance energy literacy and innovation. The report concludes with recommendations for expanding the tool by adding new levels, integrating assessment metrics, and exploring commercialization or deployment pathways through future Sapling funding. The official webpage of the Hydropower Game, which includes a link to the educational tool, can be accessed at www.anl.gov/hydropower/hydropower-game.

13 HYDRO ENERGY↗
DOE OSTI2021

Multi-Fidelity Active Subspaces for Wind Farm Uncertainty Quantification

Wind plants operate in stochastic environments characterized by complex turbulent flow dynamics and high-dimensional random variables. A key step in uncertainty quantification studies is sensitivity analysis and dimension reduction that can facilitate the development of surrogate models to be used for forward and inverse propagation or optimization under uncertainty. Prior work has shown active subspaces are an effective tool for identifying important directions in the space of stochastic inputs; however, they have only been applied to single-fidelity wind plant models. In this study, we investigate the efficacy of a multi-fidelity active subspace method for analyzing the uncertainty in wind plant power output. The multi-fidelity active subspace estimator offers the promise of increased accuracy in identifying active subspaces as compared to a single-fidelity estimator for the same computational cost, or a reduction in cost for the same accuracy. This makes the study of uncertainty in larger wind plants and with higher fidelity physics tractable. The multi-fidelity active subspace method is applied to gridded and existing wind plant layouts with single and multiple inflow conditions and its performance for surrogate modeling and uncertainty propagation is compared against a single-fidelity active subspace method. This multi-fidelity approach yields substantial computational speedups of 2x - 3.4x across the test cases along with acceptable accuracy in surrogate modeling and computing statistical moments.

active subspace↗
DOE OSTI2022

Numerical Modeling of Laser Heating and Evaporation of a Single Droplet

Laser technology is being widely studied for controlled energy deposition for a range of applications, including flow control, ignition, combustion, and diagnostics. The absorption and scattering of laser radiation by liquid droplets in aerosols affects propagation of the laser beam in the atmosphere, while the ignition and combustion characteristics in combustion chambers are influenced by the evaporation rate of the sprayed fuel. In this work, we present a mathematical model built on OpenFOAM for laser heating and evaporation of a single droplet in the diffusion-dominated regime taking into account absorption of the laser radiation, evaporation process, and vapor flow dynamics. The developed solver is validated against available experimental and numerical data for heating and evaporation of ethanol and water droplets. The two main regimes—continuous and pulsed laser heating—are explored. For continuous laser heating, the peak temperature is higher for larger droplets. For pulsed laser heating, when the peak irradiance is close to transition to the boiling regime, the temporal dynamics of the droplet temperature does not depend on the droplet size. With the empirical normalization of time, the dynamics of the droplet shrinkage and cooling are found to be independent of droplet sizes and peak laser intensities.

42 ENGINEERING↗
DOE OSTI2023

Sizing and Selection of Pressure Relief Valves for High-Pressure Thermal–Hydraulic Systems

This study covers the critical concerns related to the sizing, selection, installation, maintenance, and testing of pressure safety valves (PSVs). The aim is to ensure the safety of pressurized systems, hydrostatic transmission systems, and hydraulic plants, including process plants, thermal power plants, and nuclear reactor systems. PSVs are devices that ensure the safety and reliability of pressurized vessels, lines, and systems during overpressure events. The task of selecting which PSV features are of greatest value for a specific purpose is complex—especially in the design of a high-pressure experimental thermal–hydraulic facility for hydrostatic and transient testing of the reactor system—when the systems are in the design and development phases and require qualification and demonstration to prove that they have reached a given level of technological readiness. The present study highlights the required steps for users to follow the associated rules, guidelines, and recommendations. As a part of this research, case studies are presented to help readers better understand the applicable strategy and standards. A discussion and a review of PSV performance degradation and failure are summarized to provide a better understanding of varied process applications and conditions, including fluid flow dynamics, boundary-layer formation and pressure drops, gas bubble formation and collapse, geometric configurations, inlet/outlet piping, abrupt pressure fluctuations, and acoustic resonance. Moreover, this study discusses the servicing and testing of PSVs in a multiphase pressurized system. Overall, it provides a basic overview of how PSVs ensure the safety of pressurized systems, supported by case studies and industrial practices.

02 PETROLEUM↗
DOE OSTI2021

Identification of Preferential Recharge Zones in Karst Systems Based on the Correlation between the Spring Level and Precipitation: A Case Study from Jinan Spring Basin

The Jinan spring basin is located in the karst area of northern China, where springs serve as important sources of water supply. Several studies on spring protection and water supply have been carried out, and scholars have developed some laws on local groundwater flow dynamic and characteristics of aquifer structures. Unfortunately, there is a lack of detailed research on preferential recharge zones, which are the main recharge pathways of springs. Therefore, this research focuses on identifying preferential recharge zones based on the correlation between the spring level and precipitation. The results show that when precipitation is more intense or lasts longer, there is a stronger correlation between spring level and precipitation. It has been established that the precipitation at Donghongmiao station has the closest relationship with the dynamic of Baotu spring, which is found to be the most significant contribution to spring preservation. Two potential preferential recharge zones in the Jinan spring basin are detected through correlation analysis and geological exploration data. These findings support spring protection and water supply projects in karst regions.

58 GEOSCIENCES↗
DOE OSTI2023

Bridging Scales in Black Hole Accretion and Feedback: Magnetized Bondi Accretion in 3D GRMHD

Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magnetohydrodynamic (GRMHD) code KHARMA that allows us to span 7 orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with a Bondi radius of R B ≈ 2 × 10 5 GM•/c 2 , where M• is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as ~r -1 rather than r -3/2 and the accretion rate $\dot{M}$ is consequently suppressed by over 2 orders of magnitude below the Bondi rate $\dot{M}$ B . Here, we find continuous energy feedback from the accretion flow to the external medium at a level of $\thicksim$ 10 -2 $\dot{M}$c 2 $\thicksim$ 5 X 10 -5 $\dot{M}$ B c 2 . Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback.

79 ASTRONOMY AND ASTROPHYSICS↗