Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model validation and calibration”

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 505 records · Page 28

Acquisition and Processing Considerations for Infrared Images of Rotating Turbine Blades

Abstract As designers aim to increase efficiency in gas turbines for aircraft propulsion and power generation, spatially resolved experimental measurements are needed to validate computational models and compare improvement gains of new cooling designs. Infrared (IR) thermography is one such method for obtaining spatially resolved temperature measurements. As technological advances in thermal detectors enable faster integration times, surface temperature measurements of rotating turbine blades become possible to capture including the smallest features. This paper outlines opportunities enabled by the latest IR detector technologies for capturing spatially resolved rotating blade temperatures, while also addressing some of the challenges of implementing IR for turbine rigs such as the one in the Steady Thermal Aero Research Turbine (START) Laboratory. This paper documents critical steps in achieving accurate measurements including calibration, integration times, spatial noise, and motion blur. From these results, recommendations are provided for achieving accurate IR measurements collected in a rotating turbine facility to study film cooling.

Engineering↗

Going Forward with the Nancy Grace Roman Space Telescope Transient Survey: Validation of Precision Forward-modeling Photometry for Undersampled Imaging

The Nancy Grace Roman Space Telescope (Roman) is an observatory for both wide-field observations and coronagraphy that is scheduled for launch in the mid-2020s. Part of the planned survey is a deep, cadenced field or fields that enable cosmological measurements with type Ia supernovae (SNe Ia). With a pixel scale of 0$_.^"$11, the Wide Field Instrument will be undersampled, presenting a difficulty for precisely subtracting the galaxy light underneath the SNe. We use simulated data to validate the ability of a forward-model code (such codes are frequently also called “scene-modeling” codes) to perform precision supernova photometry for the Roman SN survey. Our simulation includes over 760,000 image cutouts around SNe Ia or host galaxies (~10% of a full-scale survey). To have a realistic 2D distribution of underlying galaxy light, we use the VELA simulated high-resolution images of galaxies. We run each set of cutouts through our forward-modeling code which automatically measures time-dependent SN fluxes. Given our assumed inputs of a perfect model of the instrument point-spread functions and calibration, we find biases at the millimagnitude level from this method in four red filters (Y106, J129, H158, and F184), easily meeting the 0.5% Roman inter-filter calibration requirement for a cutting-edge measurement of cosmological parameters using SNe Ia. Simulated data in the bluer Z087 filter shows larger ~ 2–3 mmg biases, also meeting this requirement, but with more room for improvement. Our forward-model code has been released on Zenodo.

79 ASTRONOMY AND ASTROPHYSICS↗

Solubility Predictions for Hanford Tank Concentrates with HSC - 20459

The commercially available Outotec HSC 9.9.0{sup R} provides excellent predictions for the solubilities of many species for Hanford waste tank concentrates. The HSC solubility predictions with a Harvie activity model agreed well with five boildowns that showed solubility limits for NaNO{sub 3}, NaNO{sub 2}, Na{sub 2}CO{sub 3}, NaF, Na{sub 2}SO{sub 4}, and Na{sub 3}PO{sub 4}. However, the HSC did not predict density and therefore water content very well and therefore we used the measured densities and derived water content from the difference between measured density and assays. Furthermore, we adjusted the HSC enthalpy for the minor species Na{sub 3}PO{sub 4}*NaF double salt in order to fit the five boildowns as well. Previous work has shown that gibbsite becomes and remains supersaturated in all of these boildowns and so there was no alumina precipitation despite supersaturation. Furthermore, an equilibrium model like HSC always shows complete oxidation of any TOC in the presence of nitrate and nitrite. We therefore suppressed the HSC oxidation of acetate and oxalate in order to fit measured acetate and oxalate solubilities as TOC. Space in Hanford waste tanks is very limited and therefore evaporator campaigns continue to recover needed tank space from dilute tank liquids. Prior to each evaporator campaign, a boildown of a sample of tank liquid for each campaign provides a basis for planning by measuring density, water activity, soluble species, and solids precipitation as a function of concentration at 18 deg. C. These boildowns then provide validations and calibrations for the HSC solubility prediction for Hanford wastes. Previous work has shown the solubilities of minor electrolytes in these five boildowns follow rather simple power law solubility products. Hanford's 222-S laboratory has reported each of five boildowns for concentrated liquids from five waste tanks and tank blends. All boildowns eventually showed the precipitation of the dominant electrolyte, NaNO{sub 3}, and overall solution activity largely drives all solubility. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Developing an experimental-computational framework to investigate the deformation mechanisms and mechanical properties of Al-8Ce-10Mg alloys at micro and macroscales

There is a promising future for the use of aluminum-cerium-magnesium alloys in a broad range of applications, including devices that operate at high temperatures. With cerium currently considered a waste product of rare earth mining validation studies of possible applications are essential to reduce the environmental waste. Here, a computational-experimental framework is developed to investigate the role both the intermetallic and matrix have on the mechanical properties of these alloys. A set of experiments, including SEM/EBSD imaging, nanoindentation, in-situ SEM tensile testing, and in-situ SEM-DIC tests are performed to characterize the microstructure and mechanical properties of these alloys. Furthermore, the elastic, plastic, and failure deformation mechanisms of the microstructure, and their correlation with the bulk scale mechanical properties are investigated. Experimental results are also used to calibrate parameters for a crystal plasticity finite element model, by performing a computational framework that minimizes the error between the computational and experimental results. This model is then utilized to investigate how the area percentage of intermetallics and the crystallographic texture control the mechanical properties of the alloy. Simulation results show that an increase in the percentage of intermetallics increase the strength but decrease the ductility of the alloy. Also, a change in material texture improves strength and reduces damage that leads to material failure. The development of the crystal plasticity model, as discussed in this work, opens opportunities for future investigations of similar aluminum-cerium-magnesium alloys.

36 MATERIALS SCIENCE↗

An analytic and complete equation of state for condensed phase materials

Analytic equations of state (EOS) are intended to reproduce theoretical and experimental data in a single phase portion of the thermodynamic space. We devise a complete and thermodynamically consistent model with four distinct features: (1) a reference isotherm that remains thermodynamically stable, (2) a flexible specific heat model based on a fourth-order rational polynomial, (3) a Grüneisen parameter that depends on specific volume and temperature, and (4) pressure and internal energy functions that can be inverted analytically in temperature. The model aims to improve the accuracy of existing equations of state while remaining computationally efficient. To demonstrate its features, we include calibrations for single-crystal pentaerythritol tetranitrate (PETN), liquid nitromethane (NM), and hexagonal close-packed beryllium (Be) metal. The parameter optimization uses the specific heat capacity, Grüneisen parameter, and static compression curves obtained from density functional theory for the crystalline solids and molecular dynamics simulations for liquid NM. We also present a velocity autocorrelation function that yields accurate phonon densities of states for the EOS calibration from the molecular dynamics trajectories. Each of the three calibrations is constrained to enforce the ambient state from experimental measurements and validated against experimental Hugoniot data from multiple sources. We also include one-dimensional hydrodynamic simulations of the isentropic compression experiments for beryllium conducted at the Z facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Data-Driven Constitutive Model for the Inelastic Response of Metals: Application to 316H Steel

Here, predictions of the mechanical response of structural elements are conditioned by the accuracy of constitutive models used at the engineering length-scale. In this regard, a prospect of mechanistic crystal-plasticity-based constitutive models is that they could be used for extrapolation beyond regimes in which they are calibrated. However, their use for assessing the performance of a component is computationally onerous. To address this limitation, a new approach is proposed whereby a surrogate constitutive model (SM) of the inelastic response of 316H steel is derived from a mechanistic crystal plasticity-based polycrystal model tracking the evolution of dislocation densities on all slip systems. The latter is used to generate a database of the expected plastic response and dislocation content evolution associated with several instances of creep loading. From the database, a SM is developed. It relies on the use of orthogonal polynomial regression to describe the evolution of the dislocation content. The SM is then validated against predictions of the dead load creep response given by the polycrystal model across a range of temperatures and stresses. When the SM is used to predict the response of 316H during complex non monotonic loading, extrapolating to new loading conditions, it is found that predictions compare particularly well against those from the physics-based polycrystal model.

36 MATERIALS SCIENCE↗

Long-term calibration and validation of stability of the Auger Engineering Radio Array using the diffuse Galactic radio emission

The Auger Engineering Radio Array (AERA) measures radio emission from high-energy extensive air showers. Consisting of 153 autonomous radio-detector stations spread over 17 km^2, it detects radio waves in the frequency range of 30 to 80 MHz. Accurate characterization of the detector response is crucial for proper interpretation of the collected data. Previously, this was achieved through laboratory measurements of the analog chain and simulations and measurements of the antenna's directional response. In this paper, we perform an absolute calibration using the continuously monitored sidereal modulation of the diffuse Galactic radio emission. Calibration is done by comparing the average frequency spectra recorded by the stations with predictions from seven different models of the full radio sky, accounting for the system response, which includes the antenna, filters, and amplifiers. The analysis of the calibration constants over a period of seven years shows no relevant and no significant ageing effect in the AERA antennas. This result confirms the long-term stability of the detector stations and demonstrates the possibility for a radio detector to effectively monitor ageing effects of other detectors operating over extended periods.

Antennas↗

Prognostic importance of the preoperative New‐Naples prognostic score for patients with gastric cancer

Abstract Background The wide applicability of the Naples prognostic score (NPS) is still worthy of further study in gastric cancer (GC). This study aimed to construct a New‐NPS based on the differences in immunity and nutrition in patients with upper and lower gastrointestinal tumors to help obtain an individualized prediction of prognosis. Methods This study retrospectively analyzed patients who underwent radical gastrectomy from April 2014 to September 2016. The cutoff values of the preoperative neutrophil‐to‐lymphocyte ratio (NLR), lymphocyte‐to‐monocyte ratio (LMR), serum albumin (Alb), and total cholesterol (TC) were calculated by ROC curve analysis. ROC and t‐ROC were used to evaluate the accuracy of the prognostic markers. The Kaplan–Meier method and log‐rank test were used to analyze the overall survival probability. Univariate and multivariate analyses based on Cox risk regression were used to show the independent predictors. The nomogram was made by R studio. The predictive accuracy of nomogram was assessed using a calibration plot, concordance index (C‐index), and decision curve. Results A total of 737 patients were included in training cohort, 411 patients were included in validation cohort. ROC showed that the New‐NPS was more suitable for predicting the prognosis of GC patients. NPS = 2 indicated a poor prognosis. Multivariate analysis showed that CEA ( P = 0.026), Borrmann type ( P = 0.001), pTNM ( P < 0.001), New‐NPS ( P < 0.001), and nerve infiltration ( P = 0.035) were independent risk factors for prognosis. Conclusion The New‐NPS based on the cutoff values of NLR, LMR, Alb, and TC is not only suitable for predicting prognosis but can also be combined with clinicopathological characteristics to construct a nomogram model for GC patients.

Wang, Hao↗

Energy Impact of Connected and Automated Vehicle Technologies. Final report

The overarching goal of this project is to understand the potential impact of connected and automated vehicles. The goal was achieved through data collection, model development, algorithm designs, simulations, and limited field tests. The main outcomes from this project include: (1) we collected energy consumption and GPS data from 500 vehicles over one year, with a total mileage of 8 million miles; (2) Based on the collected data and other datasets collected at the University of Michigan, we developed a calibrated Ann Arbor model in Polaris (model developed by ANL), and the fuel economy accuracy was found to be around 3.9% by comparison with field collected data; (3) An open-source SUMO model of Ann Arbor was developed; (4) Eco-Routing algorithms in Ann Arbor using the SUMO model shows 6% fuel saving potential; (5) Experiments conducted at the Mcity test facility shows that human drivers roughly follow the Eco-driving suggestions roughly 70% of the time; (6) Based in the Ann Arbor travel patterns, we found that each shared automated vehicle can replace around 4 individually owned vehicles; and (7) Adaptive Traffic Signal Control Algorithm developed through this project has been validated both in simulations and preliminary test results. For connected and automated vehicles, on average the performance is 13% delay reduction, and 10% fuel reduction. While connected and automated vehicles are in their early stage of deployment, the results from this project confirm that there is significant potential for energy saving if the technologies are developed and used properly. The three main technologies studied in this project include eco-routing, shared autonomous rides, and adaptive traffic signal controls. The data collected and model developed through the project can be used to study many other connected and automated vehicle technologies.

02 PETROLEUM↗

APOLLO: a facility-scale differentiable virtual accelerator for Fermilab

As the design complexity of modern accelerators grows, there is more interest in using advanced simulations that have fast execution time or yield additional insights like gradients. The FAST/IOTA facility has been working on implementing and experimentally validating an end-to-end digital twin that is both fast and gradient-aware, allowing for rapid prototyping of new software and experiments with minimal beam time costs. Our framework integrates physics and ML codes for linac and ring simulation through a set of generic interfaces between surrogate and physics-based sections. To reproduce device inputs and outputs, system state is exposed as a deterministic discrete event simulator. Because Fermilab is undergoing control system transition, both EPICS and ACNET frontends are supported. Recently, we have begun transitioning to a new community lattice standard, PALS, as well as developing standardized infrastructure for data ingest and normalization to prepare for model calibration during FAST proton injector commissioning. We discuss implementation details as well as challenges, and future plans to extend modelling to main complex proton accelerators like PIPII and Booster.

Kuklev, Nikita [Fermilab]↗

Transit Rider/Travel Behavior Inventory Survey - Minneapolis-St. Paul Metro - 1990

The 1990 transit on-board survey aimed to update the 1988 survey, which was conducted as part of the Preliminary Engineering Study for the Hennepin County Light Rail Transit System. The 1988 study received Regional Transit Board funding, and survey results could be applied to a mode split model for projecting ridership on the proposed Light Rail Transit System. However, this survey was designed with the 1990 Travel Behavior Inventory in mind. The intention had been to update the 1988 survey in 1990 to be compatible with the 1990 Travel Behavior Inventory data. The results of the 1990 survey were used primarily to create a table of observed transit trips between each of the 1,200 traffic analysis zones in the region. This trip table was used to calibrate a new mode split model, which estimated future year travel by mode. The 1990 update survey focused on new routes and routes that had changed significantly since 1988. To preserve compatibility with the 1988 survey, the same survey questionnaire card was used, together with the same survey procedures for data collection. The procedures randomly sampled bus patrons during the transit trip, asking key questions about the patron and the transit trip. The survey card was intended for patrons to fill out quickly so it could be completed during the transit trip. The questions focused on conditions that have proven over time to significantly influence ridership. In all, a total of 20,126 valid survey records were processed. Adding surrogate trips, the survey data file is composed of a total of 27,159 trip records . About 10% of the records were filled out by persons who had answered more than one questionnaire.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A better understanding of the mechanics of borehole breakout utilizing a finite strain gradient-enhanced micropolar continuum model

Borehole breakout denotes the failure in rock mass subjected to drilling, caused by stress concentrations exceeding the material strength. Depending on the material properties, the preexisting in situ stress state, and the borehole dimensions, different types of borehole breakout, such as spiral-shaped breakout or v-shaped breakout, are distinguished in the literature. In the present work, we address the influence of the material properties on the predicted borehole breakout mode in a comprehensive finite element study. To this end, herein we employ a gradient-enhanced micropolar damage-plasticity model based on the Mohr–Coulomb strength criterion, formulated in the finite strain regime, which is calibrated based on experimental results from plane strain compression tests on Red Wildmoor sandstone. In the numerical study, the influence of the in situ stress state, the material friction angle, plastic dilation, post peak residual strength, and the inherent material length scale parameters are investigated. Thereby, we demonstrate that depending on the material parameters, substantially different failure modes, characterized by strongly localized shear bands or diffuse failure zones, are predicted. It is shown that in particular the brittleness of the material in the post peak regime has a major influence on the predicted breakout type. Moreover, a statistical validation of the results is obtained by considering different random field distributions of the initial material strength.

58 GEOSCIENCES↗

Achieving 50% Energy Savings in Chicago Homes: A Case Study for Advancing Equity and Climate Goals

Since 2019, the National Renewable Energy Laboratory (NREL) and Elevate have collaborated to identify pathways to deep energy retrofits in Chicago's housing stock, document equity implications and co-benefits of this transition, and validate the findings by implementing retrofits in real Chicago homes. This document summarizes our analysis process to model advanced retrofit packages that lead to greater than 50% energy savings in Chicago homes. Based on these findings, we have also developed a roadmap with the City to guide implementation, and are deploying the recommended retrofit packages in real Chicago homes to realize these energy savings. This work was developed in collaboration with two key stakeholders - the City of Chicago and Commonwealth Edison (ComEd) - and funded by the U.S. Department of Energy (DOE). NREL's Residential Buildings team maintains the best-in-class ResStock™ energy model of the U.S. residential building stock. For this work, we calibrated ResStock to Chicago's unique local housing stock to accurately simulate energy use in Chicago homes both for current conditions and with various retrofit scenarios. We simulated a wide range of potential building retrofits covering all aspects of residential energy use and then grouped these into packages based on energy and utility bill savings. ResStock can model diverse building types and housing characteristics, so we're able to observe the range of outcomes that might occur when these upgrades are deployed across the entire housing stock. We can then estimate potential energy savings from an advanced retrofit program on Chicago's housing stock by comparing the modeled energy use before versus after a retrofit. This novel version of ResStock, calibrated to Chicago with data from Elevate, can help City officials, ComEd, and other partners plan for community-scale decarbonization via residential retrofits. Specifically, this work contributes the following project goals: Develop a building retrofit prioritization strategy for Chicago single-family and 2- to 4-unit buildings; Identify neighborhoods and home types that have the highest potential for savings from electrification; and Assess the impact of advanced building retrofits on energy use, utility bills, and CO2 emissions at the city and building level. Although this study is specific to Chicago, its methods and learnings are applicable across the United States. These findings are especially notable for heat pumps and electrification retrofits in cold climates.

building energy modeling↗

Coupling Remote Sensing With a Process Model for the Simulation of Rangeland Carbon Dynamics

Rangelands provide significant environmental benefits through many ecosystem services, which may include soil organic carbon (SOC) sequestration. However, quantifying SOC stocks and monitoring carbon (C) fluxes in rangelands are challenging due to the considerable spatial and temporal variability tied to rangeland C dynamics as well as limited data availability. We developed the Rangeland Carbon Tracking and Management (RCTM) system to track long-term changes in SOC and ecosystem C fluxes by leveraging remote sensing inputs and environmental variable data sets with algorithms representing terrestrial C-cycle processes. Bayesian calibration was conducted using quality-controlled C flux data sets obtained from 61 Ameriflux and NEON flux tower sites from Western and Midwestern US rangelands to parameterize the model according to dominant vegetation classes (perennial and/or annual grass, grass-shrub mixture, and grass-tree mixture). The resulting RCTM system produced higher model accuracy for estimating annual cumulative gross primary productivity (GPP) (R 2 > 0.6, RMSE <390 g C m -2 ) relative to net ecosystem exchange of CO 2 (NEE) (R 2 > 0.4, RMSE <180 g C m -2 ). Model performance in estimating rangeland C fluxes varied by season and vegetation type. The RCTM captured the spatial variability of SOC stocks with R 2 = 0.6 when validated against SOC measurements across 13 NEON sites. Model simulations indicated slightly enhanced SOC stocks for the flux tower sites during the past decade, which is mainly driven by an increase in precipitation. Future efforts to refine the RCTM system will benefit from long-term network-based monitoring of vegetation biomass, C fluxes, and SOC stocks.

54 ENVIRONMENTAL SCIENCES↗

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs: Preprint

Subsurface data analysis, reservoir modeling, and machine learning (ML) techniques have been applied to the Brady Hot Springs (BHS) geothermal field in Nevada, USA to further characterize the subsurface and assist with optimizing reservoir management. Hundreds of reservoir simulations have been conducted in TETRAD-G and CMG STARS to explore different injection and production fluid flow rates and allocations and to develop a training data set for ML. This process included simulating the historical injection and production since 1979 and prediction of future performance through 2040. ML networks were created and trained using TensorFlow based on multilayer perceptron (MLP), long short-term memory (LSTM), and convolutional neural network (CNN) architectures. These networks took as input selected flow rates, injection temperatures, and historical field operation data and produced estimates of future production temperatures. This approach was first successfully tested on a simplified single fracture doublet system, followed by the application to the BHS reservoir. Using an initial BHS dataset with 37 simulated scenarios, the trained and validated network predicted the production temperature for 6 production wells with the mean absolute percentage error of less than 8%. In a complementary analysis effort, the principal component analysis applied to 13 BHS geological parameters revealed that vertical fracture permeability shows the strongest correlation with fault density and fault intersection density. A new BHS reservoir model was developed considering the fault intersection density as proxy for permeability. This new reservoir model helps to explore under-exploited zones in the reservoir. A data gathering plan to obtain additional subsurface data was developed; it includes temperature surveying for three idle injection wells, at which the reservoir simulations indicate high bottom-hole temperatures. The collected data assist with calibrating the reservoir model and may lead to converting these wells to producers to access under-exploited zones in the reservoir. Data gathering activities are planned for the first quarter of 2021.

40 EE - Geothermal Technologies Office (EE-4G)↗

Analysis of Element Size and Element Type on SS304L Material Model Performance

Finite element solutions are only valid if the information inputted into the models is a valid assumption. This is why it is important to determine if the model behaves as expected especially in the material response. Material models that consider deformation in the elastic and plastic regions are desirable. As a material deforms plastically, hardening will occur until the ultimate strength is reached due to the reduction of the cross sectional area in a uniaxial tensile test environment. A phenomenological approximation of this behavior which is known as power law hardening (plh) can be expressed as σ = σ y + K $ϵ^{n}_{p}$ where K and n are empirically determined constants, and σ y is the proportional limit of the material. The first objective of this study is to verify the behavior for a calibrated plh material model for selective deviatoric hexahedral elements when applied to a different material stock. The Sierra SM input for this calibrated model can be found in the appendix. The second objective of this study is to determine the differences in failure characteristics for considerably large shell elements versus that of fine hexahedral elements.

36 MATERIALS SCIENCE↗

MnEdgeNet for accurate decomposition of mixed oxidation states for Mn XAS and EELS L2,3 edges without reference and calibration

Accurate decomposition of the mixed Mn oxidation states is highly important for characterizing the electronic structures, charge transfer and redox centers for electronic, and electrocatalytic and energy storage materials that contain Mn. Electron energy loss spectroscopy (EELS) and soft X-ray absorption spectroscopy (XAS) measurements of the Mn L2,3 edges are widely used for this purpose. To date, although the measurements of the Mn L2,3 edges are straightforward given the sample is prepared properly, an accurate decomposition of the mix valence states of Mn remains non-trivial. For both EELS and XAS, 2+, 3+, and 4+ reference spectra need to be taken on the same instrument/beamline and preferably in the same experimental session because the instrumental resolution and the energy axis offset could vary from one session to another. To circumvent this hurdle, in this study, we adopted a deep learning approach and developed a calibration-free and reference-free method to decompose the oxidation state of Mn L2,3 edges for both EELS and XAS. A deep learning regression model is trained to accurately predict the composition of the mix valence state of Mn. To synthesize physics-informed and ground-truth labeled training datasets, we created a forward model that takes into account plural scattering, instrumentation broadening, noise, and energy axis offset. With that, we created a 1.2 million-spectrum database with 1-by-3 oxidation state composition ground truth vectors. The library includes a sufficient variety of data including both EELS and XAS spectra. By training on this large database, our convolutional neural network achieves 85% accuracy on the validation dataset. We tested the model and found it is robust against noise (down to PSNR of 10) and plural scattering (up to t/λ = 1). We further validated the model against spectral data that were not used in training. In particular, the model shows high accuracy and high sensitivity for the decomposition of Mn 3 O 4 , MnO, Mn 2 O 3 , and MnO 2 . The accurate decomposition of Mn 3 O 4 experimental data shows the model is quantitatively correct and can be deployed for real experimental data. Our model will not only be a valuable tool to researchers and material scientists but also can assist experienced electron microscopists and synchrotron scientists in the automated analysis of Mn L edge data.

25 ENERGY STORAGE↗

Tensorized Interior Radiative Heat Transfer for a Scalable and Calibrated Building Energy Simulator

Building energy simulation is a critical tool for developing and testing advanced control strategies, such as Reinforcement Learning (RL), to provide demand flexibility and affordable energy costs. The recently introduced Smart Buildings Control Suite (sbsim) provides a lightweight, scalable, and data-calibrated simulation environment based on a 2D finite-difference model. However, the initial model primarily focused on conductive and convective heat transfer, neglecting the significant impact of long-wave radiative heat exchange between interior surfaces. This paper presents a significant extension to the sbsim framework by incorporating a physically-grounded model for interior radiative heat transfer. Our primary contribution is the development and integration of a fully tensorized radiative heat transfer module, which preserves the computational efficiency and scalability of the original simulator. This was achieved by developing a pipeline for view factor calculation, including an algorithm to identify directly seeing surfaces within complex floor plans, and formulating the net radiation equations for efficient execution on modern hardware accelerators. We validate the numerical accuracy of our tensorized implementation by comparing its results against a traditional iterative approach, demonstrating identical outcomes. This enhancement increases the physical fidelity of sbsim, enabling more accurate training of RL agents for building energy optimization.

Ham, Sang woo↗