SEARCH · Search NASA

Results for “loading simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Comparative Analysis of HEATNETS for Geothermal Network Performance: Preprint

Thermal energy networks (TENs), also known as 5th generation district energy systems, or more specifically geothermal networks when exchanging heat with geothermal boreholes, are an important technology for decarbonization. In these networks an ambient loop connects buildings and thermal sources, such as a borehole field, to exchange energy and maintain a desired loop temperature. Water-source heat pumps are used at the buildings to connect to the ambient or thermal loop to meet to the building heating and cooling loads and maintain comfort. A semi-transient, reduced-order technical model and techno-economic model, called HEATNETS, has been developed at NREL that captures the flow of energy around a TEN. In this work, a comparison of the HEATNETS technical model and a well-known coding platform used for modeling geothermal networks, TRNSYS, has been completed for a proposed geothermal network as a verification process. Hourly data provided from the TRNSYS simulation included building loads, pumping power, heat pump power, temperature entering and leaving the borehole field, and mass flow rates. The hourly borehole temperatures were used to create a linear regression model utilized in HEATNETS to estimate the borehole field heat exchange. The building loads and mass flow rates were direct inputs to HEATNETS while the pumping power, heat pump power, borehole temperatures, and coefficients of performance were all simulated and calculated by HEATNETS, allowing for direct comparison of the thermal energy transfer, rather than also comparing control systems responses. HEATNETS considers the full process from design inputs to economic outputs and can provide modeling options for high-level initial system design and operational optimization. This study shows that HEATNETS, while not intended to replace other modeling tools, can be a unique modeling tool for the performance of a full geothermal network system.

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Comparative Analysis of HEATNETS for Geothermal Network Performance

Thermal energy networks (TENs), also known as 5th generation district energy systems, or more specifically geothermal networks when exchanging heat with geothermal boreholes, are an important technology for decarbonization. In these networks an ambient loop connects buildings and thermal sources, such as a borehole field, to exchange energy and maintain a desired loop temperature. Water-source heat pumps are used at the buildings to connect to the ambient or thermal loop to meet to the building heating and cooling loads and maintain comfort. A semi-transient, reduced-order technical model and techno-economic model, called HEATNETS, has been developed at NREL that captures the flow of energy around a TEN. In this work, a comparison of the HEATNETS technical model and a well-known coding platform used for modeling geothermal networks, TRNSYS, has been completed for a proposed geothermal network as a verification and validation process. Hourly data provided from the TRNSYS simulation included building loads, pumping power, heat pump power, temperature entering and leaving the borehole field, and mass flow rates. The hourly borehole temperatures were used to create a linear regression model utilized in HEATNETS to estimate the borehole field heat exchange. The building loads and mass flow rates were direct inputs to HEATNETS while the pumping power, heat pump power, borehole temperatures, and coefficients of performance were all simulated and calculated by HEATNETS, allowing for direct comparison of the thermal energy transfer HEATNETS considers the full process from design inputs to economic outputs and can provide modeling options for high-level initial system design and operational optimization. This study focuses on a validation of HEATNETS using results from TRNSYS. HEATNETS is not intended to replace other modeling tools, but this work demonstrates, via a comparison with an industry standard code, that HEATNETS can be a unique, high-level and rapid modeling tool for estimating the performance of a full geothermal network system.

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

Simulating dispatchable grid services provided by flexible building loads: State of the art and needed building energy modeling improvements

End-use electrical loads in residential and commercial buildings are evolving into flexible and cost-effective resources to improve electric grid reliability, reduce costs, and support increased hosting of distributed renewable generation. This article reviews the simulation of utility services delivered by buildings for the purpose of electric grid operational modeling. We consider services delivered to (1) the high-voltage bulk power system through the coordinated action of many, distributed building loads working together, and (2) targeted support provided to the operation of low-voltage electric distribution grids. Although an exhaustive exploration is not possible, we emphasize the ancillary services and voltage management buildings can provide and summarize the gaps in our ability to simulate them with traditional building energy modeling (BEM) tools, suggesting pathways for future research and development.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

Heat Integration Optimization and Dynamic Modeling Investigation for Advancing the Coal-Direct Chemical Looping Process

The purpose of the project is to address the optimization and startup operation of a modular coal direct chemical looping (CDCL) combustion system integrated with a steam cycle for power generation to reduce the risks involved in further scale-up of the technology. The modular reactor design of the CDCL process provides flexibility in the fabrication of the reactor and in its operating capacity (i.e. turndown ratio) at the cost of a more complex heat exchange network (HEN) design and integration. To address the technology gaps and advance the efficiency and economic feasibility of the CDCL technology, the project will perform a detailed and comprehensive analysis of the integration of a modular CDCL reactor system and a steam cycle system under both static and transient conditions via HEN process performance simulations and system dynamic modeling, respectively. The scope of work consists of 1) Experimental and computational studies of the CDCL combustor reactor 2) Comprehensive static (i.e. steady-state) system HEN design analysis in CDCL 550 MWe commercial unit for power generation and 3) Dynamic modeling of site specific design of 10MWe CDCL large pilot plant. The project team has successfully developed and validated a kinetic model for the oxidation of oxygen carriers in the combustor using the unreacted shrinking core model (UCSM). The model is capable of capturing the oxidation kinetics of fully or partially reduced oxygen carrier particles. A computational fluid dynamics (CFD) model is developed to simulate the hydrodynamics, heat transfer, and chemical reaction occurring in the CDCL combustor. The model is developed in MFIX and ANSYS Fluent. Key aspects of CDCL combustor operation, including heat transfer, oxygen carrier oxidation, and the transport of oxygen carrier particles, are simulated using this CFD model. The HEN for a commercial scale 550 MWe CDCL power plant is simulated and optimized using ASPEN Plus. Practical design considerations are incorporated based on industrial experiences. The performance and cost for the commercial CDCL plant is updated based on these analyses. A dynamic model for the 10 MWe CDCL pilot plant is developed in ProTRAX simulation software. The model is based on the pilot plant design developed in project DE-FE0027654 “10 MWe CDCL Large Pilot Plang – Pre-FEED Study” and the steam cycle data obtained from Dover Light & Power plant. The transient behaviors during pilot plant load variation are simulated using the dynamic model.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2024

Simulation-Guided Decision-Making for Enhancing Energy Resilience in a Remote Alaskan Community

The increasing frequency and severity of extreme weather events underscore the need to bolster the resilience of energy infrastructure in remote coastal communities exposed to climate hazards. In recent years, considerable effort has been made to harden the grid infrastructure of remote communities through investment in energy storage, advanced metering, renewable generation and energy-efficient loads. However, simulation-based studies are still needed to identify appropriate locations for investment and determine the adequacy of existing infrastructure in supporting new resources. Preparing the required models may often be challenging due to insufficient metering, disparate data sources and workforce limitations. This paper describes modeling efforts undertaken to represent, within a unified co-simulation platform: (a) the electric distribution network and (b) the thermal behavior of a community medical center, in the remote community of Cordova, Alaska. With the help of case-studies, it is shown how the developed simulation platform can help the local utility in making decisions regarding capital investment aimed at enhancing community resilience.

Microgrid, energy storage, resilience↗

DOE OSTI2022

The Effect of Inverter Loading Ratio on Energy Estimate Bias: Preprint

Subhourly effects, particularly variability in solar irradiance, can lead to underestimation of inverter clipping losses and overestimation of energy in hourly photovoltaic system performance models, particularly for systems with high inverter loading ratios. Direct simulation of this error can be complicated by factors such as the representation of spatial and temporal variability in hourly weather data and transient system conditions. In this work we take an alternative approach using real system power measurements to show that energy predictions from typical industry models suffer from a bias that increases with inverter loading ratio. We also show that this loading ratio-dependent bias is strongly correlated with an empirical subhourly inverter clipping bias derived from real power plant data. Finally, we show that this bias is not necessarily specific to any one model or weather dataset by recreating similar biases with alternatives of each.

clipping↗

DOE OSTI2025

Spatial description of dislocation nucleation in the shock response of single-crystal aluminum

Nonequilibrium molecular dynamics simulations of shock loaded single-crystal Al in the $\langle$100$\rangle$, $\langle$110$\rangle$, $\langle$111$\rangle$, and $\langle$123$\rangle$ orientations are conducted to study elastic and plastic shockwave formation and details associated with dislocation activity. A computer vision-based approach is implemented to capture the presence of dislocations and describe their spatial characteristics in the zone of nucleation behind the propagating shockwave. The methodology developed relies on the sequences of images extracted during shock loading that show dislocation activity within a cross section of the sample. Results reveal that the spacing between activated slip systems is orientation dependent and exhibits a modest reduction for the $\langle$100$\rangle$ and $\langle$111$\rangle$ orientations as shock pressure increases. Comparisons are made to existing theoretical models. Such relationships between shock pressure and dislocation activity, extracted from molecular dynamics simulations, can be used to inform higher length scale simulations or modeling of dislocation-based plasticity during shock.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Wind farm response to mesoscale-driven coastal low level jets: a multiscale large eddy simulation study

Realistic atmospheric turbulence–wind farm interactions during coastal low-level jet (LLJ) events are captured using high-fidelity, mesoscale-driven large eddy simulations (LES) to understand wind turbine loads, wakes and overall performance. The simulation has been carried out using the ExaWind aeroelastic solver, AMR-Wind. The simulations have been compared against a baseline unstable case matching the wind speed, wind direction and TI at hub-height location. Results indicate that the LLJ has negative impacts on the turbine hub and tower loads, and opens up potential avenues for design load mitigation strategies.

17 WIND ENERGY↗

DOE OSTI2022

MEITNER Resource Team Modeling and Simulation Support to Holos-Quad Reactor Development (Final CRADA Report)

The main objective of this CRADA is to provide Argonne National Laboratory (ANL)’s modeling and simulation capabilities via ARPA-E's MEITNER (Modeling-Enhanced Innovations Trailblazing Nuclear Energy Reinvigoration) Resource Team (RT) arrangement to support the demonstration of the viability of HolosGen’s Holos-Quad reactor design. The Holos-Quad reactor design is an advanced reactor concept that incorporates many new design features, such as the neutron-coupled Subcritical Power Modules (SPMs) in its core design, elimination of balance of plant (BOP) by direct integration of a helium Brayton cycle power conversion system with each SPM, among many other innovative features. The demonstration of the viability of such an innovative reactor design warrants iterations of modeling and simulation and testing. The purpose of this project is to utilize ANL’s modeling and simulation capabilities in nuclear reactor analysis and power conversion system analysis to inform HolosGen and the Design Team (DT) in the design and optimization of the Holos-Quad concept. The major work scopes of this project include: to investigate conceptual designs and materials for radiation shielding to protect personnel and internal components such as turbomachinery; Assess the the helium Brayton cycle power conversion system performance in both nominal and load following conditions; Investigate power conversion components and overall system performance; Identify control strategies to enable load following; Perform simulations of HolosGen’s subscale simulator and using available test data for code validation/benchmark purpose; Perform core thermal-hydraulics, safety analysis, and structural analysis.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

Achieving 5,000-h and 8,000-h Low-PGM Electrode Durability on Automotive Drive Cycles

Whereas total Pt loading in anode and cathode catalysts below 0.125 mg cm −2 is required to meet the stringent cost target for automotive fuel cell systems (FCS) for light duty vehicles, low-loaded cathode catalysts are susceptible to unacceptable aging-related performance losses at high current densities. A framework model, validated by accelerated stress test data, has identified cell voltage, relative humidity (RH) and temperature as the key operating variables that affect degradation of a high-activity d-PtCo/C cathode catalyst with 0.1 mg cm −2 Pt loading. Drive cycle simulations indicate that these can be controlled by properly selecting the minimum FCS power, compressor-expander module (CEM) turndown, and stack coolant temperature. The optimum system parameters are 4-kW e minimum power for an 80-kW e FCS, CEM turndown of 12.5, and 66 °C average coolant exit temperature that combine to limit the maximum cell voltage to 850 mV and outlet RH to 90%–100%. Depending on Pt loading, the mismatch between actual and allowable degradation for 10% power loss over 5,000-h lifetime requires the stack to be oversized by 2.4%–5%, resulting in 8.4%–41% lower Pt utilization and 7.1%–20.5% penalty in stack cost. The corresponding results for 8,000-h lifetime are 10.3%-14% stack oversizing, 23%–51.8% lower Pt utilization, and 24.1%–35.4% stack cost penalty.

08 HYDROGEN↗

DOE OSTI2025

Real-time Simulation Model of a Vanadium Redox Flow Battery Energy Storage System with Power Electronics Integration

This paper presents a real-time (RT) model of a 20 kW multi-stack vanadium redox flow battery (VRFB) system including power electronics converters in a single real-time framework. Implemented on a Typhoon HIL 604 device, this model is designed to simulate the operation of the flow battery with system-level controls. To provide a realistic overview of a VRFB system, the model incorporates a closed-loop control system for managing the speeds of two centrifugal pumps and an electro-thermal model of the battery stacks, enabling realtime temperature estimation. Furthermore, a two-stage power electronics topology and associated control solution are presented, operating with distinct strategies for interfacing the battery model with the main grid and a local load. Results from simulations demonstrate that the VRFB model is well suited for RT environment, and system can manage power during both charging and discharging cycles, while a SCADA monitor is used to display critical variables for safe, efficient, and reliable operation of the VRFB systems.

Rezende Da Costa Reis Kimpara, Renata [ORNL] (ORCI↗

DOE OSTI2022

Physics Guided Simulation of Electrostatic Discharge: Technical Report

Triboelectrically-charged objects may create threshold sparks, electrostatic discharge (ESD) events, to equilibrate charge between themselves and other relatively charged objects. ESD events exhibit many complex physical phenomena. They are a nexus of several fields of physics with disparate characteristic scales: plasma physics, chemical kinetics, hydrodynamics, circuit models, etc. These scales can span many orders of magnitude from the varied collisions thermalizing information within a plasma on the $\mathcal{O}(fs/ps)$ to the physical size of the plasma channel on the $\mathcal{O}(100µm)$, to the speed of a nonlinear hydrodynamic wave propagating at $\mathcal{O}(µm, ns)$. These threshold ESD events may occur in situations of programmatic importance, delivering energy and power profiles to a “victim load” generating deleterious consequences. To predict and mitigate these consequences we must answer questions about the spark’s energy budget: how much energy goes into producing the spark channel; how much gets radiated away; how much energy is advected away into the hydrodynamics; and how much energy is delivered to a victim load. An ESD simulation toolset has been created and evolved in order to answer these questions. An appropriate, physics guided implementation for simulation can be done by gaining insight into its constituent physics and leveraging that intuition to choose a suitable numerical operator. We examine in detail the chemical kinetics, circuit discharge, and hydrodynamics to deter mine dominant regimes, values, timescales, and interactions to uncover the underlying physical dynamics. We also examine and propose model reduction schemes for high-dimensional chemical kinetics. We use past and current work with experimentally validated and theoretically-verified hydrodynamics to calculate applicability limits of the non-ionizing strong shock limit. We quantify the energy budget from a hydrodynamic perspective and demonstrate that a significant fraction of the stored energy is “earmarked” for hydrodynamic advection as an energy terminus. Lastly, we combine the constituent physics of an ESD event (chemical kinetics, circuit model, and hydrodynamics) into a cohesive, actionable toolset and obtain promising results from an isothermal test case. We then propose a viable, modular evolution of the ESD toolset based upon the performed examination of the physics uncovering dominant physical scales and the stiffness of the compositional differential system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Evaluation of Distributed Power Apportioning with Net Load Management Engine in Microgrids Using Power Hardware-in-the-Loop Simulation

This article presents the performance evaluation of ratio consensus-based distributed power apportioning engine along with centralized net load management (NLM) engine that ensures viable and stable operation of an islanded microgrid. Managing net load variability in a microgrid with high penetrations of uncertain renewable generation and ever-changing load demands is a crucial need in order to ensure viable and stable operation of the microgrid. Centralized “dispatch-rule”-based and/or multi-agent-based distributed control of distributed energy resources (DERs) in microgrid are well accepted for microgrid by adapting ANSI/ISA-95-based hierarchical control architecture. In the application where microgrid network has large geographical span with multiple DERs dispersed in the network, high penetration of uncertain renewable energy resources, and ever-changing load demands, a judicious selection of techniques/solutions for managing net-load resources for maintaining viability and stability is required. With this motivation, this article proposes a novel solution to mitigate the challenges by incorporating a mixed centralized NLM engine and distributed power apportioning control of DERs and loads. A power-hardware-in-the-loop (PHIL) -based experiment is conducted with the centralized NLM engine and the distributed power apportioning engine along with two commercial inverters. The experimental results validates the efficacy of the proposed method in ensuring viability and stability of a microgrid.

hardware-in-the-loop simulation↗

DOE OSTI2025

Understanding Pore Filling Processes and Adsorption/Desorption Hysteresis in Nanoporous Metal–Organic Frameworks: Insights from Grand Canonical Monte Carlo Simulations and Free Energy Calculations

Grand canonical Monte Carlo (GCMC) simulations were used to investigate pore filling and hysteresis in nanoporous metal-organic frameworks (MOFs). Adsorption and desorption isotherms were calculated for argon at 87 K in 1866 MOFs from the CoRE MOF database and for short n-alkanes in selected MOFs, keeping the adsorbent structure rigid. Analysis of the molecular configurations showed two different mechanisms and origins of hysteresis: one involving a transition of the adsorbate arrangement in the pores similar to a gas-to-liquid transition associated with a large change in the loading and one more similar to a liquid-to-solid transition associated with a relatively small change in the loading. Our GCMC simulations in MOFs with diverse pore topologies indicate exceptions to an empirical relationship for the minimum diameter of a cylindical pore required for hysteresis as a function of the adsorbate diameter and reduced temperature. The simulations reveal some structures where isotherms exhibit two steps in the adsorption branch and only one step in the desorption branch. Hysteresis loops with a different number of adsorption and desorption steps are not common. Here, to better understand why hysteresis is observed in the GCMC simulations, the concept of the transition probability for observing a step in the adsorption isotherm at a given pressure in a GCMC simulation is introduced. We used two different methods to calculate the transition probabilities and find that these yield comparable results. Furthermore, the transition probability provides a measure for the length of GCMC simulations to yield reliable results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Modeling the role of local crystallographic correlations in microstructures of Ti-6Al-4V using a correlated structure visco-plastic self-consistent polycrystal plasticity formulation

This paper presents a multi-level crystal plasticity-based simulation framework for modeling mechanical response and microstructure evolution of Ti-6Al-4V with α-lath/lamellar microstructures. Here, the model is a correlated structure visco-plastic self-consistent (CS-VPSC) formulation linking three scales: a single crystals micro-scale, a lath/lamellar colony meso-scale, and a lath/lamellar aggregate macro-scale. A selected hardening law for the evolution of critical resolved shear stress per slip system used in CS-VPSC is phenomenological. However, it adjusts the resistances of basal and prismatic slip systems based on the geometry of slip transfer between adjacent lamellae. Consistent with experimental evidences, the resolved shear stress on the pyramidal slip planes is dependent not only on the stress in the direction of slip but also on the two orthogonal shear stress components and the three normal stress components (non-Schmid effects). Electron backscatter diffraction (EBSD) data in conjunction with a procedure relying on α→β phase transformation is used to construct paired variants of α-laths/lamellae satisfying their local crystallographic correlations. The procedure fits volume fractions of individual laths/lamellae with the experimental EBSD data and selects a distribution of habit planes between adjacent variants with respect to the loading direction. The simulation framework is applied to interpret the deformation behavior in tension and compression along two sample directions of Ti-6Al-4V fabricated via laser powder bed fusion. Moreover, the model is used to simulate texture evolution during rolling of the material to large strains. It is demonstrated that the model is capable of predicting plastic anisotropy/asymmetry and the concomitant texture evolution. While the model reveals a significant effect of habit plane inclination with respect to the loading direction on yield stress, the comparison of the data and model predictions shows that a random distribution of habit planes fits the flow response. It is further inferred that the tension-compression asymmetry arises from the non-Schmid effects.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Electromagnetic turbulence simulation of tokamak edge plasma dynamics and divertor heat load during thermal quench

Abstract The edge plasma turbulence and transport dynamics, as well as the divertor power loads during the thermal-quench phase of tokamak disruptions, are numerically investigated with BOUT++’s flux-driven six-field electromagnetic turbulence model. Here, transient yet intense particle and energy sources are applied at the pedestal top to mimic the plasma power drive at the edge induced by a core thermal collapse, which flattens the core temperature profile. Interesting features, such as surging of divertor heat load (up to 50 times) and broadening of heat-flux width (up to four times) on the outer-divertor target plate, are observed in the simulation, in qualitative agreement with experimental observations. The dramatic changes in divertor heat load and width are due to the enhanced plasma turbulence activities inside the separatrix. Two cross-field transport mechanisms, namely, the E × B turbulent convection and the stochastic parallel advection/conduction, are identified to play important roles in this process. First, an elevated edge pressure gradient drives instabilities and subsequent turbulence in the entire pedestal region. The enhanced turbulence not only transports particles and energy radially across the separatrix via the E × B convection, which causes the initial divertor heat-load burst, but it also induces amplified magnetic fluctuation B ˜ . Once themagnetic fluctuation is large enough to break the magnetic flux surface, magnetic flutter effect provides an additional radial transport channel. In the late stage of our simulation, | B ˜ r / B 0 | reaches to 10 −4 level that completely breaks magnetic flux surfaces such that stochastic field lines are directly connecting pedestal top plasma to the divertor target plates or first wall, further contributing to the divertor heat-flux width broadening.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Idealized Offshore Low‐Level Jets for Turbine Structural Impact Considerations

Low-level jets (LLJs) describe conditions in which the wind speed reaches a local maximum with respect to altitude near the surface; they have been observed intermittently in the US mid-Atlantic offshore environment. LLJs pose unique operating conditions for future wind turbines operating in the region, such as negative shear and locally strong veer, but they are not typically considered in existing turbine standards. This work builds upon recent research that explains the formation and evolution of the US mid-Atlantic LLJs through a simple analytical governing equation. We generate several LLJ inflow conditions with varying jet characteristics based on this analytical model and create monotonically sheared (MS) analogues with matched veer to assess the impacts of the LLJ on turbine performance and loading. Using aeroelastic simulations with these inflow conditions on the International Energy Agency 15-MW reference turbine, we find that the LLJ leads to a greater range of tower-top pitching and yawing moments, which could contribute to larger accumulated structural fatigue in components compared to MS inflow. These preliminary results demonstrate a path toward a unified set of test cases for low-level wind maxima that can inform the International Electrotechnical Commission standards related to offshore wind turbine design.

17 WIND ENERGY↗

DOE OSTI2023

A Simulation Framework for Analyzing the Impact of Stochastic Occupant Behaviors on Demand Flexibility in Typical Commercial Buildings

As one of the primary users of the electric grid, buildings and building equipment, including heating, ventilation, and air conditioning (HVAC) systems, can be leveraged to provide the flexible demand needed to balance the grid. Typical strategies to achieve demand flexibility are to reduce electricity use during peak or critical periods by shutting down equipment or relaxing system setpoints, which will inevitably impact the occupants’ comfort. When occupants feel uncomfortable, they may take actions to regain their comfort, and some of those actions (such as turning on a personal fan) may have a negative impact on meeting the demand response goal. Therefore, it is important to incorporate occupant behaviors into the assessment ofthe building demand flexibility potential. In this study, a simulation framework that includes simulation of zone thermal loads, an HVAC system, and occupant behaviors, was developed to investigate the impact of occupant behaviors on demand flexibility. A case study was conducted using a small office model from the U.S. Department of Energy (DOE) Commercial Prototype Building Models to simulate the building envelope and zone loads. An agent-based occupant thermal behavior model was adapted to forecast occupants’ thermal comfort and their resulting thermal behaviors. An artificial neural network (ANN) based airflow model trained from a computational fluid dynamics (CFD) model of the zone was adopted to better predict the ambient environment of each occupant. An air-source heat pump simulation model that was calibrated from a real two-stage air-source heat pump system was used as the HVAC system. A typical load shedding event during peak hours was studied. Repeated simulations were conducted to capture the stochastic effects of occupant behaviors. The interplay between the demand flexibility, occupant comfort and behavior were analyzed by evaluating key performance indicators, including the energy use, occupant discomfort duration, and occupant behavior duration during the peak period. The results suggest that this framework can be used to analyze typical commercial buildings and their HVAC systems in terms ofdemand flexibility potential under the impact of occupant behaviors.

Chen, Zhelun↗