Search NASA⌕ Search

SEARCH · Search NASA

Results for “multi-zone”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Air quality and comfort constrained energy efficient operation of multi-zone buildings

Maintaining indoor air quality (IAQ) through effective ventilation is essential for the well-being and productivity of building occupants. Control strategies aimed at improving the efficiency of heating, ventilation and air conditioning (HVAC) systems must jointly determine ventilation and heating and cooling processes. Here, in this paper, we study the problem of minimizing the energy consumption of the HVAC system in a multi-zone building, while meeting thermal comfort and IAQ requirements. We first perform a steady state analysis of the zonal carbon dioxide (CO 2 ) concentration and the temperature dynamics. The resulting expressions are convex in the zonal mass flow rates and zonal temperatures. Guided by the steady state solutions for meeting the thermal comfort constraints, we develop two control policies for improving the energy efficiency of building HVAC systems while jointly satisfying indoor temperature and IAQ constraints. We compare the performance of our proposed approaches with those of multiple baseline approaches which implement separate regimes for controlling zonal temperature and IAQ for a typical work-day in a multi-zone campus building. We have evaluated the performance of our proposed approaches under varying levels of flexibility in zonal temperatures. We have shown that zonal temperature flexibility can result in energy savings up to 32% (for the same control strategies) as compared to the case where no such flexibility is permitted. Our proposed approaches were seen to offer potential savings of nearly 29% compared to the baseline.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Physics vs structure: A systematic benchmark of learning strategies for multi-zone building thermal dynamics

Recent advances in physics-informed and data-driven machine learning promise improved thermal models for advanced building control, yet there is limited quantitative evidence on when added physics structure and architectural complexity are beneficial. Here, this work presents a systematic benchmark of five representative system identification methods for modeling multi-zone building thermal dynamics: linear state-space models, multi-layer perceptrons, neural state-space models, neural ordinary differential equations, and physically-consistent neural networks. The methods are evaluated across multiple data regimes and zone coupling strategies. Using a high-fidelity multi-zone commercial building emulator, we examine short-term and long-term prediction accuracy, computational efficiency, and ease of development. Our results reveal critical trade-offs between prediction performance, model complexity, and physical consistency. We demonstrate that decoupled, nonlinear black-box models consistently outperform coupled physics-constrained architectures in both predictive accuracy and out-of-distribution robustness in majority of the test cases for the building type considered in the study. Our findings quantify the cost of complexity in building thermal modeling and provide concrete, actionable, scenario-based guidelines for selecting model classes for control-oriented applications.

Building thermal modeling↗
DOE OSTI2025

Online Dynamic Mode Decomposition Based System Identification of Multi-Zone Building HVAC Systems

Many works have recently been conducted to reduce the electricity consumption of smart buildings and allow them to support various grid services. Most of these works require accurate system models for the various appliances in the building including heating, ventilation, and air conditioning (HVAC) units. In this paper, we investigate a recursive data-driven system identification strategy to construct the thermal model for a time-varying building with a multi-zone HVAC unit. The online dynamic mode decomposition (DMD)-based strategy is employed to identify the multi-zone thermal building dynamics, where a simple information update (rank-1) is selected to avoid computational complexity. The DMD-based identification strategy is validated using a real gymnasium building equipped with a 4-zone HVAC unit, and its performance is compared with that of the traditional nuclear-norm subspace identification (N2SID) strategy.

Wu, Tumin [University of Tennessee, Knoxville (UTK↗
DOE OSTI2023

Performance analysis and comparison of data-driven models for predicting indoor temperature in multi-zone commercial buildings

Building thermal models, which characterize the properties of a building’s envelope and thermal mass, are essential for accurate indoor temperature and cooling/heating demand prediction. Because of their flexibility and ease of use, data-driven models are increasingly used. Here, this study compared and analyzed the performance of gray-box (resistance-capacitance) and black-box (recurrent neural network) models for predicting indoor air temperature in a real multi-zone commercial building. The developed resistance-capacitance model served as a benchmark model for which full sets of temporal data and building information were used as inputs. The recurrent neural network models were trained and tested assuming various available types and amounts of temporal data and known building physical information to investigate the effects of data and information availability. Feature importance analysis was conducted to select the key variables for different prediction targets under different scenarios. This research provides guidance in selecting an appropriate building thermal response modeling method based on the measured data availability, building physical information, and application.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Datasets of Faults in Variable Air Volume Terminal Units in a Multi-Zone Commercial Building

Faults in HVAC systems can decrease system efficiency and equipment lifespan, leading to 5%–30% of energy consumption being wasted in commercial buildings. We identified two common faults in HVAC variable air volume systems: a stuck damper fault in the variable air volume terminal unit and a discharge airflow sensor fault. We conducted three sets of damper stuck tests and two sets of airflow sensor tests, each including a fault-free scenario and scenarios with varying levels of faults, over one day. The faults were implemented in Oak Ridge National Laboratory’s two-story Flexible Research Platform building to generate a high-quality, well-controlled dataset covering fault-induced and fault-free scenarios. The test building, fault test scenarios, and data validation are described here. The open-source dataset includes 1 min intervals of weather and building data on the presence and absence of building faults. This dataset can be used to analyze the effects of HVAC system faults on system operation and indoor building conditions, and to develop or evaluate a fault detection and diagnosis algorithm.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Non-Rare Earth Magnesium Bumper Beams

Magna will team with the Pacific Northwest National Laboratory (PNNL) to develop Shear Assisted Processing and Extrusion (ShAPE™) of rectangular and multi-zoned extrusion profiles using non-RE Mg. It is envisioned that the ShAPE™ process, developed during this project for non-circular profiles, will serve as the foundation for fabrication of full-scale bumper beams of the future where near net-shapes are extruded directly from castings in a single step, without the need for costly rare earth additives. In this project, the ShAPE™ process will be developed for a rectangular profile measuring 1” x 2” with a 0.06” wall thickness. Successful extrusion of rectangular cross sections is a critical step along the path toward profiles with multiple zones. With the process established for a rectangular profile, the project will then make advancements toward multi-zoned cross sections. A magnesium alloy bumper beam offers a 35% weight savings compared to aluminum alloys and 60% compared to high strength steels. The bumper beam was chosen as the first component to target because the geometry can have a fairly simple 0.6-0.12 inch wall/web thickness, generally compatible with PNNL’s knowledge base and equipment capability surrounding ShAPE™. The bumper component offers commercial value by providing mass reduction in the front of the vehicle along with ride and handling benefit. PNNL has demonstrated the ability to ShAPE™ process non-RE Mg tubing with high strength, ductility and energy absorption. It has been shown that ShAPE™ extruded ZK60 and Mg2Si tubing with 0.30 inch diameter and 0.03 inch wall thickness have similar energy absorption to AA6061-T6. This is due to significant grain refinement, breakdown and dispersion of second phases, and the ability to eliminate anisotropy in compressive/tensile strength by aligning texture 45 degrees to the extrusion axis. Scalability of the ShAPE™ process has been demonstrated in ZK60 by fabricating 2.0 inch diameter tubing with 0.06 and 0.12 inch wall thicknesses having concentricity as tight as 1.3%. These tubes exhibit an ultimate tensile strength of 260-300 MPa and room temperature ductility of >20% parallel, perpendicular and 45 degrees to the extrusion axis. Aluminum bumpers represent approximately 54% of the NA light vehicle market. If the research effort is successful, providing a 30% mass reduction at equal cost relative to aluminum bumpers, it is anticipated that 25% market penetration would be realized over the next 5 year and 50% market share by 2030. 100% of the aluminum automotive bumpers are sourced from the tier network, providing resourcing opportunity based on equal performance, reduced mass and cost equivalent lower than aluminum bumpers. Due to the relatively low vehicle content, magnesium components are currently shredded and 100% recycled with aluminum content in the US, presenting a positive life cycle impact.

42 ENGINEERING↗
DOE OSTI2023

Quantifying spatiotemporal variability in occupant exposure to an indoor airborne contaminant with an uncertain source location

Well-mixed zone models are often employed to compute indoor air quality and occupant exposures. While effective, a potential downside to assuming instantaneous, perfect mixing is underpredicting exposures to high intermittent concentrations within a room. When such cases are of concern, more spatially resolved models, like computational-fluid dynamics methods, are used for some or all of the zones. But, these models have higher computational costs and require more input information. A preferred compromise would be to continue with a multi-zone modeling approach for all rooms, but with a better assessment of the spatial variability within a room. Here to do so, we present a quantitative method for estimating a room's spatiotemporal variability, based on influential room parameters. Our proposed method disaggregates variability into the variability in a room's average concentration, and the spatial variability within the room relative to that average. This enables a detailed assessment of how variability in particular room parameters impacts the uncertain occupant exposures. To demonstrate the utility of this method, we simulate contaminant dispersion for a variety of possible source locations. We compute breathing-zone exposure during the releasing (source is active) and decaying (source is removed) periods. Using CFD methods, we found after a 30 minutes release the average standard deviation in the spatial distribution of exposure was approximately 28% of the source average exposure, whereas variability in the different average exposures was lower, only 10% of the total average. We also find that although uncertainty in the source location leads to variability in the average magnitude of transient exposure, it does not have a particularly large influence on the spatial distribution during the decaying period, or on the average contaminant removal rate. By systematically characterizing a room's average concentration, its variability, and the spatial variability within the room important insights can be gained as to how much uncertainty is introduced into occupant exposure predictions by assuming a uniform in-room contaminant concentration. We discuss how the results of these characterizations can improve our understanding of the uncertainty in occupant exposures relative to well-mixed models.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

A bi-level data-driven framework for fault-detection and diagnosis of HVAC systems

Long-term operation of heating, ventilation, and air conditioning (HVAC) systems will eventually lead to a range of HVAC system failures, resulting in excessive energy consumption and maintenance costs. Here, to avoid HVAC malfunctioning, fault detection diagnostic (FDD) is utilized as a common practice. Machine learning methods have lately received considerable interest for FDD analysis of HVAC systems due to their high detection accuracy. Meanwhile, HVAC malfunctions are regarded as rare occurrences, hence normal operating data samples are much more accessible than data samples in faulty and malfunctioning conditions. The dominating frequency of normal operation in HVAC datasets has also led to heavily biased classification algorithms within the literature. Moreover, the focus of previous literature has been on increasing the accuracy of the models which leads to a high number of false positives (misleading alarms) in the system. In order to enhance the performance of diagnostic procedures and fill the mentioned gaps, this study proposes a novel data-driven framework. A bi-level machine learning framework is developed for diagnosing faults in air handling units (AHUs) and rooftop units (RTUs) based on principal component analysis (PCA), time series anomaly detection, and random forest (RF). It is shown that PCA can reduce the dataset dimension with one principal component accounting for 95% of data variance. Also, the random forest could classify the faults with 89% precision for single-zone AHU, 85% precision for RTU, and 79% for multi-zone AHU. By proposing this framework, three persistent challenges are addressed: (I) minimizing false positives; (II) accounting for data imbalance; and (III) normal condition monitoring of equipment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Data reduction considerations for the burning velocity of spherical constant volume flames of R32 (CH 2 F 2 ) with air

Here, the present work explores data reduction techniques for the measurement of the laminar burning velocities of R32(CH 2 F 2 )-air mixtures using a constant volume combustion device, in which the pressure-time history is the only measured parameter. To allow clear assessment of the accuracy of the data reduction methods, the pressure-time histories used for analysis are synthetically generated via a detailed numerical simulation employing full kinetics and with and without an optically-thin radiation model. Various data reduction models are employed, including a two-zone model and two multi-zone models, and these are compared with the results from the burning velocity obtained from the output of the numerical simulation. The data reduction schemes are shown to be accurate if the same radiation model is employed in the data reduction as was used in the flame simulation to generate the pressure trace used for post-processing. If the incorrect radiation model is employed, however, the errors can be quite large. The effects of stretch, radiation, and different data post-processing methodologies are explored and the errors quantified. Stretch is shown to be important for the early stages and the selected data range that is used for extrapolation has a significant effect on the extrapolated burning velocity. However, with an appropriate choice of data considered for extrapolation, the prediction of the unstretched burning velocity can be quite accurate.

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

Hybrid combustion modeling approach for turbulent jet ignition in natural-gas pre-chamber spark-ignition engines at high EGR

Here, this study presented a hybrid modeling approach for simulating turbulent jet ignition and combustion processes in a natural-gas pre-chamber spark-ignition engine operating under exhaust gas recirculation (EGR) diluted conditions. In-depth analyses of experimental data and simulation results from previous work [Chinnathambi et al., ICEF2021-67836; Kim et al., Fuel 409: 137815, 2026] revealed two key findings: (i) the magnitude of pressure difference between the pre-chamber and main chamber ($∆P_{PC-MC}$) was positively correlated with the combustion duration from the moment of $∆P_{PC-MC}=0$ to the point of 5% mass fraction burned, with larger $∆P_{PC-MC}$ associated with longer duration; and (ii) the turbulent combustion regime in the main chamber transitioned from the broken reaction zone to the corrugated flamelet regime, with the Karlovitz number exceeding 100 immediately after turbulent hot jets were ejected from nozzles, coinciding with observed local extinction events. To accurately simulate the entire combustion process, a hybrid approach was developed under Reynolds-Averaged Navier Stokes framework, combining the G-equation model for pre-chamber combustion with the multi-zone well-stirred reactor approach and a turbulence-chemistry interaction (TCI) submodel for main chamber combustion. The TCI submodel accounted for the attenuation of reaction rates due to turbulent strain and modeled local extinction by suppressing reaction rates under certain flow and flame conditions. When applied to three EGR rate conditions toward the dilution limit, the hybrid modeling approach accurately reproduced experimental data in terms of cylinder pressure, apparent heat release rate, and the observed positive correlation, including the delayed onset of main chamber combustion—a feature not captured by existing combustion models.

computational fluid dynamics simulation↗
DOE OSTI2023

Evaluation of thermostat location for multizone commercial building performance

In multi-zone buildings, it is often found that a single shared thermostat controls more than one conditioned zones. Although these shared zones are supposed to have similar thermal needs (e.g., cooling and heating load), in reality, they are not mainly due to different orientations, sizes of windows, occupancy, space types, etc. This can cause unnecessary energy waste or thermal discomfort for the occupants. How to quantify this impact in multizone buildings remains a research gap. Therefore, this study aims to evaluate the impact of different sensor (i.e., thermostat) locations for multizone commercial buildings through a comprehensive modeling study. Here, two different scenarios for the sensor locations were selected to evaluate the impact in terms of energy and thermal comfort. The scenario (1) is that one to five sensors distributed among the five zones, but the sensor readings from selected zones will be used for no-sensor zones, which is no-mean sensor scenario. The scenario (2) is that one to five sensors distributed among the five zones, but the average temperature from the shared zones will be used for each of the shared zones, which is a mean sensor scenario. The uncertainty analysis was performed for different sensor location scenarios.(a)The major findings from an energy perspective, for scenario (1), the differences of cooling energy go as high as 17% more or 12% less, compared with the baseline. For heating energy consumption, the discrepancies go as high as 51% more or 52% less, compared with baseline. For site energy consumption, the discrepancies go as high as 3.2% more or 3.2% less, compared with baseline. For fan energy consumption, the discrepancies go as high as 3.2% more or as low as 1.0% less, compared with baseline. For scenario (2), the discrepancies of cooling energy go as high as 3% more, or 0.5% less, compared with the baseline. For heating energy consumption, the discrepancies, are go as high as 10.1% more or as low as 3.0% less, compared with baseline. For site energy consumption, the discrepancies go as high as 1.3% more or 0.3% less, compared with baseline. For fan energy consumption, the discrepancies go as high as 3.1% more or as low as 1.0% less, compared with baseline.(b) In terms of the indoor thermal comfort, for the no-mean-sensor scenarios, the discrepancies of unmet hours for cooling mode can be as high as 1,200 h, compared with the baseline. The discrepancies of unmet hours for heating mode can be as high as 740 h, compared with the baseline. For the mean-sensor scenarios, the discrepancies of unmet hours for cooling mode can be as high as 750 h, compared with the baseline. The discrepancies of unmet hours for heating mode can be as high as 50 h, compared with the baseline.

42 ENGINEERING↗
DOE OSTI2025

EnergyPlus-MCP: A model-context-protocol server for ai-driven building energy modeling

Traditional building energy modeling with the EnergyPlus building performance simulation engine requires domain expertise, programming skills, and intensive manual efforts limiting its effective adoption. This paper introduces EnergyPlus-MCP, the first open-source Model Context Protocol (MCP) server specifically designed for EnergyPlus simulation workflows, establishing a new foundational infrastructure for AI-driven building energy modeling. The MCP server implements a layered architecture with 35 specialized tools spanning model management, editing and analysis, HVAC and other systems configuration inspection, and simulation execution, enabling Large Language Models to interact with EnergyPlus through conversational interfaces. The server addresses critical workflow barriers by automating model validation, streamlining energy efficiency measures modification, and providing intelligent output management with interactive visualization. Through practical demonstrations using a multi-zone building retrofit analysis, we show how the EnergyPlus-MCP server significantly reduces manual efforts while maintaining full simulation rigor. By providing accessible natural language interfaces to sophisticated building energy analysis, this approach enables scalable deployment of simulation expertise across public and private organizations, educational institutions, and research teams, fundamentally transforming traditional building energy modeling practices.

AI↗
DOE OSTI2024

Insights into the broadband emission of the TeV blazar Mrk 501 during the first X-ray polarization measurements

Aims.We present the first multiwavelength study of Mrk 501 that contains simultaneous very-high-energy (VHE)γ-ray observations and X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). Methods.We used radio-to-VHE data from a multiwavelength campaign carried out between March 1, 2022, and July 19, 2022 (MJD 59639 to MJD 59779). The observations were performed by MAGIC,Fermi-LAT,NuSTAR,Swift(XRT and UVOT), and several other instruments that cover the optical and radio bands to complement the IXPE pointings. We characterized the dynamics of the broadband emission around the X-ray polarization measurements through its multiband fractional variability and correlations, and compared changes observed in the polarization degree to changes seen in the broadband emission using a multi-zone leptonic scenario. Results.During the IXPE pointings, the VHE state is close to the average behavior, with a 0.2–1 TeV flux of 20%–50% of the emission of the Crab Nebula. Additionally, it shows low variability and a hint of correlation between VHEγ-rays and X-rays. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings, taken in March 2022 (MJD 59646 to 59648 and MJD 59665 to 59667), with a synchrotron peak frequency > 1 keV. For the third IXPE pointing, in July 2022 (MJD 59769 to 59772), the synchrotron peak shifts toward lower energies and the optical/X-ray polarization degrees drop. All three IXPE epochs show an atypically low Compton dominance in theγ-rays. The X-ray polarization is systematically higher than at lower energies, suggesting an energy stratification of the jet. While during the IXPE epochs the polarization angles in the X-ray, optical, and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. Such results further support the hypothesis of an energy-stratified jet. We modeled broadband spectra taken simultaneous to the IXPE pointings, assuming a compact zone that dominates in the X-rays and the VHE band, and an extended zone stretching farther downstream in the jet that dominates the emission at lower energies.NuSTARdata allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change in the magnetization and/or the emission region size, which directly connects the shift in the synchrotron peak to lower energies with the drop in the polarization degree.

Astronomy & Astrophysics↗
DOE OSTI2024

Controlling Air Conditioners for Frequency Regulation: A Real-World Example

Even though thermostatically controlled loads like air conditioners present a great potential for providing ancillary services to the electric power grid, the practical challenges associated with their real-time coordination have not received the necessary attention. In this work, we present a nondisruptive load control application, specifically, we demonstrate how real residential air conditioners can provide frequency regulation. Aggregate power adjustment is achieved by modifying the ON/OFF modes of the air conditioners. To account for both single and multi-zone houses, we extend the currently available techniques and develop an approach that can be used for controlling aggregations that include both types of houses. A discussion of the practical challenges encountered in our field experiments is provided, along with the hardware and software approaches we developed to circumvent them. Here, we argue that limitations of current thermostat APIs introduce significant challenges and are an impediment to widespread adoption of fast load control applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

End-Use Savings Shapes Measure Documentation: Variable Refrigerant Flow with Heat Recovery and Dedicated Outdoor Air System

Building on the successfully completed effort to calibrate and validate the U.S. Department of Energy's ResStock™ and ComStock™ models over the past 3 years, the objective of this work is to produce national data sets that empower analysts working for federal, state, utility, city, and manufacturer stakeholders to answer a broad range of analysis questions. The goal of this work is to develop energy efficiency, electrification, and demand flexibility end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) measures. "Measures" refers to energy efficiency variables that can be applied to buildings during modeling. An end-use savings shape is the difference in energy consumption between a baseline building and a building with an energy efficiency, electrification, or demand flexibility measure applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ComStock is a highly granular, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the commercial building stock across the United States. The baseline model intends to represent the U.S. commercial building stock as it existed in 2018. The methodology and results of the baseline model are discussed in the final technical report of the End-Use Load Profiles project. This documentation focuses on a single heating, ventilation, and air-conditioning (HVAC) end-use savings shape measure - a variable refrigerant flow with heat recovery (VRF HR) heating and cooling system coupled with a dedicated outdoor air system (DOAS) for ventilation. This measure replaces existing multi-zone variable air volume (VAV) systems or single-zone rooftop units (RTU) with a VRF HR system coupled with a DOAS that includes an energy/heat recovery ventilator (E/HRV). The measure covers 53% of exisiting building stock's floor area and is not applicable to HVAC system types using district heating or cooling or buildings/spaces that include high-ventilation spaces such as kitchens where the amount of exhaust air is large. A DOAS with E/HRV is used to provide required outdoor ventilation air to spaces since ventilation air is generally not supplied by a VRF HR system. An exhaust air energy recovery ventilator (ERV ) with sensible and latent heat exchange is added to humid climate zones while a heat recovery ventilator (HRV ) with sensible only exchange is added to drier climate zones. The ERV is modeled as a fixed membrane plate counterflow heat exchanger, while the HRV is modeled as a sensible-only fixed aluminum plate counterflow heat exchanger. Both systems include a bypass (for temperature control and economizer lockout) and minimum exhaust temperature control for frost prevention.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

First-principles-integrated Study of Blazar Synchrotron Radiation and Polarization Signatures from Magnetic Turbulence

Blazar emission is dominated by nonthermal radiation processes that are highly variable across the entire electromagnetic spectrum. Turbulence, which can be a major source of nonthermal particle acceleration, can widely exist in the blazar emission region. The Turbulent Extreme Multi-Zone (TEMZ) model has been used to describe turbulent radiation signatures. Recent particle-in-cell (PIC) simulations have also revealed the stochastic nature of the turbulent emission region and particle acceleration therein. However, radiation signatures have not been systematically studied via first-principles-integrated simulations. In this paper, we perform combined PIC and polarized radiative transfer simulations to study synchrotron emission from magnetic turbulence in the blazar emission region. We find that the multiwavelength flux and polarization are generally characterized by stochastic patterns. Specifically, the variability timescale and average polarization degree (PD) are governed by the correlation length of the turbulence. Interestingly, magnetic turbulence can result in polarization angle swings with arbitrary amplitudes and duration, in either direction, that are not associated with changes in flux or PD. Surprisingly, these swings, which are stochastic in nature, can appear either bumpy or smooth, although large-amplitude swings (>180°) are very rare, as expected. Our radiation and polarization signatures from first-principles-integrated simulations are consistent with the TEMZ model, except that in the latter, there is a weak correlation, with zero lag, between flux and degree of polarization.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
DOE OSTI2025

Evaluating Chemical Kinetics Predictions for Propane Using 3-D and 0-D Models in a Boosted Spark-Ignited Engine

Propane has been shown to be a promising alternative fuel to reduce emissions while simultaneously achieving high efficiencies in medium- and heavy-duty engines. These high-power density applications require boosted engines which, combined with high compression ratio, can lead to auto-ignition and knock. While three-dimensional (3-D) computational fluid dynamics (CFD) models are often used for resolving the complex fluid flow in engines, these models can become computationally expensive when simulating detailed chemical kinetics. Likewise, zero-dimensional (0-D) models are computationally concise enough for kinetics development, but lack any flow-field information which governs the flame propagation processes in spark ignition (SI) engines. This work presents a comprehensive comparison between 3-D and 0-D closed cycle simulations at knocking conditions in a high compression ratio high stroke-to-bore ratio propane engine. In order to initialize the flow-field for the 3-D closed cycle (intake valve closing, (IVC) to exhaust valve opening, (EVO)) simulation, a motored multi-cycle 3-D model was run using Converge to create a map at IVC, reducing the computational time. The map allowed a non-homogeneous 3-D closed cycle simulation to be satisfactorily validated against experiments, while a homogeneous case using only the turbulence field mapping was also simulated, mimicking 0-D modeling. The 3-D simulations were used to prescribe the initial conditions (e.g., IVC thermodynamics, speciation, burn-rate profile) for a 2-zone 0-D SI engine model in Chemkin Pro for both cases. It was found that 2-zone 0-D modeling underpredicted the knock onset timing, likely due to the lack of thermal stratification in the unburned gas region. Future work will carry multi-zone 0-D modeling to capture the fuel auto-ignition in the unburned region.

Douvry-Rabjeau, Julien [Oakland University, Roches↗