SEARCH · Search NASA

Results for “Dynamic Interaction Test”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Dynamic Interaction of Grid-Forming and Grid-Following Inverters with Synchronous Generators in Hybrid Power Plants

As the integration of renewable energy resources such as wind and solar generation continues to grow, the use of hybrid power plants with both synchronous machines and power electronic inverters operating in parallel becomes more common. This paper studies the impact on small-signal stability for these pairs of devices, first as an isolated two-device system and then as part of a 9-bus test system. Each operational scenario was studied with both a grid-following inverter and a grid-forming inverter. In each scenario, the state variables associated with the synchronous generator's excitation system were found to be responsible for the inception of instability. These results underscore the importance of accurate modeling of excitation systems when modeling power systems with hybrid power plants.

hybrid power plant↗

DOE OSTI2020

Nanobodies against the metal binding domains of ATP7B as tools to study copper transport in the cell

Abstract Nanobodies are genetically engineered single domain antibodies derived from the unusual heavy-chain only antibodies found in llamas and camels. The small size of the nanobodies and flexible selection schemes make them uniquely versatile tools for protein biochemistry and cell biology. We have developed a panel of nanobodies against the metal binding domains of the human copper transporter ATP7B, a multidomain membrane protein with a complex regulation of enzymatic activity and intracellular localization. To enable the use of the nanobodies as tools to investigate copper transport in the cell, we characterized their binding sites and affinity by isothermal titration calorimetry and NMR. We have identified nanobodies against each of the first four metal binding domains of ATP7B, with a wide affinity range, as evidenced by dissociation constants from below 10−9 to 10−6 M. We found both the inhibitory and activating nanobodies among those tested. The diverse properties of the nanobodies make the panel useful for the structural studies of ATP7B, immunoaffinity purification of the protein, modulation of its activity in the cell, protein dynamics studies, and as mimics of copper chaperone ATOX1, the natural interaction partner of ATP7B.

Uhlemann, Eva-Maria E.↗

DOE OSTI2021

Nutrient Inputs Stimulate Mercury Methylation by Syntrophs in a Subarctic Peatland

Climate change dramatically impacts Arctic and subarctic regions, inducing shifts in wetland nutrient regimes as a consequence of thawing permafrost. Altered hydrological regimes may drive changes in the dynamics of microbial mercury (Hg) methylation and bioavailability. Important knowledge gaps remain on the contribution of specific microbial groups to methylmercury (MeHg) production in wetlands of various trophic status. Here, we measured aqueous chemistry, potential methylation rates (k meth ), volatile fatty acid (VFA) dynamics in peat-soil incubations, and genetic potential for Hg methylation across a groundwater-driven nutrient gradient in an interior Alaskan fen. We tested the hypotheses that (1) nutrient inputs will result in increased methylation potentials, and (2) syntrophic interactions contribute to methylation in subarctic wetlands. We observed that concentrations of nutrients, total Hg, and MeHg, abundance of hgcA genes, and rates of methylation in peat incubations (k meth ) were highest near the groundwater input and declined downgradient. hgcA sequences near the input were closely related to those from sulfate-reducing bacteria (SRB), methanogens, and syntrophs. Hg methylation in peat incubations collected near the input source (FPF2) were impacted by the addition of sulfate and some metabolic inhibitors while those down-gradient (FPF5) were not. Sulfate amendment to FPF2 incubations had higher k meth relative to unamended controls despite no effect on k meth from addition of the sulfate reduction inhibitor molybdate. The addition of the methanogenic inhibitor BES (25 mM) led to the accumulation of VFAs, but unlike molybdate, it did not affect Hg methylation rates. Rather, the concurrent additions of BES and molybdate significantly decreased k meth , suggesting a role for interactions between SRB and methanogens in Hg methylation. The reduction in k meth with combined addition of BES and molybdate, and accumulation of VFA in peat incubations containing BES, and a high abundance of syntroph-related hgcAsequences in peat metagenomes provide evidence for MeHg production by microorganisms growing in syntrophy. Collectively the results suggest that wetland nutrient regimes influence the activity of Hg methylating microorganisms and, consequently, Hg methylation rates. Our results provide key information about microbial Hg methylation and methylating communities under nutrient conditions that are expected to become more common as permafrost soils thaw.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Anomalously high elastic modulus of a poly(ethylene oxide)-based composite electrolyte

The practical use of lithium metal anodes in solid-state batteries requires a polymer membrane with high lithium-ion conductivity, thermal/electrochemical stability, and mechanical strength. The primary challenge is to effectively decouple the ionic conductivity and mechanical strength of the polymer electrolytes. We report a remarkably facile single step synthetic strategy based on in-situ crosslinking of poly(ethylene oxide) (xPEO) in the presence of a woven glass fiber (GF). Such a simple method yields composite polymer electrolytes (CPE) of anomalously high elastic modulus up to 2.5 GPa over a broad temperature range (20 °C – 245 °C) that has never been previously documented. An unsupervised machine learning algorithm, K-mean clustering analysis, was implemented on the hyperspectral Raman mapping at the xPEO/GF interface. Using such a unique means, we show for the first time that the promoted mechanical strength originates from xPEO and GF interactions through dynamic hydrogen and ionic bonding. High ionic conductivity is achieved by the addition plasticizer (e.g. tetraglyme), where trifluoromethanesulfonate anions are tethered to the xPEO matrix and Li + cations are favorably transported through coordination with the plasticizer. Further, stringent galvanostatic cycling tests indicates the CPE can be stably cycled for >3000 h in a Li-metal symmetric cell at a moderate temperature (nearly 1500 Coulombs/cm 2 Li equivalents), outperforming most of the PEO-based electrolytes. The GF reinforced CPE reported here has multifunctional uses, such as solid electrolytes for all solid-state batteries and membranes for redox-flow batteries. Although the focus of this study is on lithium-based batteries, the results are equally promising for other alkali metal based batteries such as sodium and potassium.

25 ENERGY STORAGE↗

DOE OSTI2022

Arbuscular Mycorrhizal Tree Communities Have Greater Soil Fungal Diversity and Relative Abundances of Saprotrophs and Pathogens than Ectomycorrhizal Tree Communities

Trees associating with different mycorrhizas often differ in their effects on litter decomposition, nutrient cycling, soil organic matter (SOM) dynamics, and plant-soil interactions. For example, due to differences between arbuscular mycorrhizal (AM) and ectomycorrhizal (ECM) tree leaf and root traits, ECM-associated soil has lower rates of C and N cycling and lower N availability than AM-associated soil. These observations suggest that many groups of nonmycorrhizal fungi should be affected by the mycorrhizal associations of dominant trees through controls on nutrient availability. To test this overarching hypothesis, we explored the influence of predominant forest mycorrhizal type and mineral N availability on soil fungal communities using next-generation amplicon sequencing. Soils from four temperate hardwood forests in southern Indiana, United States, were studied; three forests formed a natural gradient of mycorrhizal dominance (100% AM tree basal area to 100% ECM basal area), while the fourth forest contained a factorial experiment testing long-term N addition in both dominant mycorrhizal types. We found that overall fungal diversity, as well as the diversity and relative abundance of plant pathogenic and saprotrophic fungi, increased with greater AM tree dominance. Additionally, tree community mycorrhizal associations explained more variation in fungal community composition than abiotic variables, including soil depth, SOM content, nitrification rate, and mineral N availability. Our findings suggest that tree mycorrhizal associations may be good predictors of the diversity, composition, and functional potential of soil fungal communities in temperate hardwood forests. These observations help explain differing biogeochemistry and community dynamics found in forest stands dominated by differing mycorrhizal types.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

A holistic framework integrating plant-microbe-mineral regulation of soil bioavailable nitrogen

Abstract Soil organic nitrogen (N) is a critical resource for plants and microbes, but the processes that govern its cycle are not well-described. To promote a holistic understanding of soil N dynamics, we need an integrated model that links soil organic matter (SOM) cycling to bioavailable N in both unmanaged and managed landscapes, including agroecosystems. We present a framework that unifies recent conceptual advances in our understanding of three critical steps in bioavailable N cycling: organic N (ON) depolymerization and solubilization; bioavailable N sorption and desorption on mineral surfaces; and microbial ON turnover including assimilation, mineralization, and the recycling of microbial products. Consideration of the balance between these processes provides insight into the sources, sinks, and flux rates of bioavailable N. By accounting for interactions among the biological, physical, and chemical controls over ON and its availability to plants and microbes, our conceptual model unifies complex mechanisms of ON transformation in a concrete conceptual framework that is amenable to experimental testing and translates into ideas for new management practices. This framework will allow researchers and practitioners to use common measurements of particulate organic matter (POM) and mineral-associated organic matter (MAOM) to design strategic organic N-cycle interventions that optimize ecosystem productivity and minimize environmental N loss.

Environmental Sciences & Ecology↗

DOE OSTI2020

Measurement of temperature of a dusty plasma from its configuration

A new concept called “configurational temperature” is introduced in the context of dusty plasma, where the temperature of the dust particles submerged in the plasma can be measured directly from the positional information of the individual dust particles and the interaction potential between the dust grains. This method does not require the velocity information of individual particles, which is a key parameter to measure the dust temperature in the conventional method. The technique is initially tested using two-dimensional (2D) OpenMP parallel molecular dynamics and Monte Carlo simulation and then compared with the temperature evaluated from experimental data. The experiment have been carried out in the Dusty Plasma Experimental (DPEx) device, where a 2D stationary plasma crystal of melamine formaldehyde particles is formed in the cathode sheath of a DC glow discharge argon plasma. The kinetic temperature of the dust is calculated using the standard particle image velocimetry technique at different pressures. An extended simulation result for the three-dimensional case is also presented, which can be employed for the temperature measurement of a three-dimensional dust crystal in laboratory devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Performance evaluation of surface riding wave energy converter with linear electric generator

In this work, we devised a new WEC (wave energy converter) called SR-WEC (Surface Riding WEC). The SR-WEC consists of two bodies: the outer cylinder with an armature assembly (body #1) and a magnet assembly (body #2) sliding inside the armature. For the SR-WEC, the relative sliding displacement and velocity are caused by gravity acceleration and the outer cylinder's motions, and they lead to electrical power generation. To evaluate its performance, a numerical simulation tool was developed, which solves the fully-coupled floater-mooring-generator dynamics. During the developing stage, the appropriate hydrodynamics model, sliding mechanics model, mooring dynamics model, and LEG (linear electric generator) electro-magnetic model were independently developed and then fully coupled in time domain to account for the cross-coupling interactions among them. Then, the developed simulation tool was verified component by component against various laboratory tests. Subsequently, systematic parametric studies were conducted with several important design parameters under various wave conditions to enhance power generation. After that, the average output power was evaluated in enlarged operational wave conditions. The present SR-WEC is particularly designed to be efficient at low sea states, which is good since they cover the majority of typical annual sea states.

42 ENGINEERING↗

DOE OSTI2021

Model predictive control for integrated control of air-conditioning and mechanical ventilation, lighting and shading systems

Modern buildings are increasingly automated and often equipped with multiple building services (e.g., air-conditioning and mechanical ventilation (ACMV), dynamic shading, dimmable lighting). These systems are conventionally controlled individually without considering their interactions, affecting the building’s overall energy inefficiency and occupant comfort. A model predictive control (MPC) system that features a multi-objective MPC scheme to enable coordinated control of multiple building services for overall optimized energy efficiency, indoor thermal and visual comfort, as well as a hybrid model for predicting indoor visual comfort and lighting power is proposed. The MPC system was implemented in a test facility having two identical, side-by-side experimental cells to facilitate comparison with a building management system (BMS) employing conventional reactive feedback control. The MPC system coordinated the control of the ACMV, dynamic façade and automated dimmable lighting systems in one cell while the BMS controlled the building services in the other cell in a conventional manner. The MPC side achieved 15.1–20.7% electricity consumption reduction, as compared to the BMS side. Simultaneously, the MPC system improved indoor thermal comfort by maintaining the room within the comfortable range (-0.5 < predicted mean vote < 0.5) for 98.3% of the time, up from 91.8% of the time on the BMS side. Visual comfort, measured by indoor daylight glare probability (DGP) and horizontal illuminance level at work plane height, was maintained for the entire test period on the MPC side, improving from having visual comfort for 94.5% and 85.7% of the time, respectively, on the BMS side.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Going beyond BEM with BEM: an insight into dynamic inflow effects on floating wind turbines

Blade element momentum (BEM) theory is the backbone of many industry-standard wind turbine aerodynamic models. To be applied to a broader set of engineering problems, BEM models have been extended since their inception and now include several empirical corrections. These models have benefitted from decades of development and refinement and have been extensively used and validated, proving their adequacy in predicting aerodynamic forces of horizontal-axis wind turbine rotors in most scenarios. However, the analysis of floating offshore wind turbines (FOWTs) introduces new sets of challenges, especially if new-generation large and flexible machines are considered. In fact, due to the combined action of wind and waves and their interaction with the turbine structure and control system, these machines are subject to unsteady motion and thus unsteady inflow on the wind turbine's blades, which could put BEM models to the test. Consensus has not been reached on the accuracy limits of BEM in these conditions. This study contributes to the ongoing research on the topic by systematically comparing four different aerodynamic models, ranging from BEM to computational fluid dynamics, in an attempt to shed light on the unsteady aerodynamic phenomena that are at stake in FOWTs and whether BEM is able to model them appropriately. Simulations are performed on the UNAFLOW 1:75 scale rotor during imposed harmonic surge and pitch motion. Experimental results are available for these conditions and are used for baseline validation. The rotor is analyzed in both rated operating conditions and low wind speeds, where unsteady aerodynamic effects are expected to be more pronounced. Results show that BEM, despite its simplicity, can adequately model the aerodynamics of FOWTs in most conditions if augmented with a dynamic inflow model.

17 WIND ENERGY↗

DOE OSTI2023

Fluid structural interaction simulation of the MP-1 plate performance irradiated in the advanced test reactor

Due to the high neutron fluxes they generate and increased heat transfer performance, plate type fuels are used in the U.S. high-performance research reactors. During irradiation, a significant amount of fission energy (i.e., ~200 MeV per atom) is released by the U-235 chain reaction that is carried away by the coolant. Through thermal hydraulic analysis, the coolant’s heat transfer capability is investigated to ensure that the desired fuel temperature can be maintained. At high temperatures, the plate undergoes elastic/plastic deformation, creep, and swelling as a result of both the temperature gradients and the fission gas production within the fuel. These effects are studied via fuel performance analysis. When the plate deformation is small enough that the coolant flow remains relatively unchanged, conducting these two types of analyses independently will suffice. But at high fuel burnups, the swollen plates may encroach into the coolant channels that separate the fuel plates from each other and cause these channels to narrow. Large reductions in channel gap size imperil cooling performance, causing fuel temperatures to rise. Further, if the plate deformation is asymmetric, the fuel centerline will shift toward one side of the channel, causing an uneven reduction in coolability. In addition, a boehmite (oxide) layer will, over time, grow on the plate surface, further obstructing the heat transfer process. To precisely predict fuel temperatures/deformation, a complete coupled analysis that considers coolant flow, heat transfer, oxide growth, elastic/plastic deformation, creep, and swelling is needed; however, this type of analysis method is not available in the literature. To fill the gap, this research developed a fluid structure interaction (FSI) approach to the fuel plate analysis, then successfully applied it to the Mini-Plate (MP)-1 experiment, which was irradiated in the Advanced Test Reactor (ATR) for both one and two cycles. The complete analysis coupled STAR-CCM+, a computational fluid dynamic (CFD) software for calculating flow, with Abaqus, a finite element analysis code for calculating plate deformation. Improvements in the results were found when comparing the fully coupled analysis to the independently conducted analyses but they were not significant due to the miniature size of the plates and the relatively short irradiation time. In the future, the fully coupled approach presented herein will be applied to full-size fuel plates with longer irradiation cycles once additional experiments become available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Multi-scale, Multi-disciplinary, and Multi-agent Explainable AI with Koopman-Undergirded Learning, Prediction, and Analysis (M3EA KULPA) (Project Closeout Report)

The goal of this project was to develop and use domain-aware machine learning formulations, based on the Koopman Operator (KO), for modelling multi-scale, multi-disciplinary (e.g., multi-physics), and/or multi-agent systems. The project developed these formulations for the following cases: • Systems with dynamics at two separate time scales, • Systems with a bi-level hierarchical control structure, • Systems with bi-level hierarchical control and dynamics at two separate time scales (the lower level controls operating at the faster time scale), and • Systems with n separate but interacting agents/disciplines (with/without control, respectively); the controls for each agent could include bi-level hierarchical control and dynamics at two separate time scales as described above. The project then defined a set of dynamical systems consisting of different nonlinear oscillators that could be used to test these different formulations and then subsequently learned the KO models for those systems. With the KO models, we were able to do the following: • Quantify system stability, including both long-term and transient behavior, • Quantify the effects of feedbacks between the different time scales and agents/disciplines in terms of those feedbacks’ effects on system stability, • Replace a standard Proportional-Integral (PI) control in the hierarchical control structure with a KO-based Linear-Quadratic Regular (LQR), a form of optimal control, • Calculate optimal supervisory control policies a) with and without time scale separated dynamics at the lower level control levels and b) with both PI and KO-based LQR lower level control policies, and • Calculate dynamic Nash equilibria for multi-agent systems where each agent makes its own control decisions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Deuteron off-resonance rotating frame relaxation for the characterization of slow motions in rotating and static solid-state proteins

Here, we demonstrate the feasibility of deuterium solid-state NMR off-resonance rotating frame relaxation measurements for studies of slow motions in biomolecular solids. The pulse sequence, which includes adiabatic pulses for magnetization alignment, is illustrated for static and magic-angle spinning conditions away from rotary resonances. We apply the measurements for three systems with selective deuterium labels at methyl groups: a) a model compound, Fluorenylmethyloxycarbonyl methionine-D 3 amino acid, for which the principles of the measurements and corresponding motional modeling based on rotameric interconversions are demonstrated; b) amyloid-β 1-40 fibrils labeled at a single alanine methyl group located in the disordered N-terminal domain. This system has been extensively studied in prior work and here serves as a test of the method for complex biological systems. The essential features of the dynamics consist of large-scale rearrangements of the disordered N-terminal domain and the conformational exchange between the free and bound forms of the domain, the latter one due to transient interactions with the structured core of the fibrils. and c) a 15-residue helical peptide which belongs to the predicted α-helical domain near the N-terminus of apolipoprotein B. The peptide is solvated with triolein and incorporates a selectively labeled leucine methyl groups. The method permits model refinement, indicating rotameric interconversions with a distribution of rate constants.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Aero-thermal numerical characterization of blunt fin-induced shock wave–boundary layer interaction and its control through leading-edge cooling injection

This study represents a novel evaluation of active flow control to alleviate the aerothermal penalties created by the blunt fin-induced shock wave–boundary layer interaction. The manuscript analyzes the effect of flow injection on a blunt fin-induced shock wave–boundary layer interaction via computational fluid dynamics simulations with various degrees of resolution. The impact on the mean flow topology and wall variables was investigated utilizing Reynolds-averaged Navier–Stokes simulations. Detached-eddy simulations revealed the low-frequency shock motion, shock wave–boundary layer, and horseshoe vortex interaction. The test article was exposed to two different incoming boundary layer thicknesses; the thicker boundary layer led to the appearance of larger turbulent scales. The Detached-eddy simulations revealed the time history of the shock wave–boundary layer interaction, focusing on the inception and development of the recirculated flow regions. Ultimately, spectral proper orthogonal decomposition was employed to identify the structures associated with the low-frequency shock motion caused by the shock wave–boundary layer interaction.

Lozano, Francisco (ORCID:0000000277871125)↗

DOE OSTI2022

On the Diffusivity of Moist Static Energy and Implications for the Polar Amplification Response to Climate Warming

Energy balance models (EBMs) have been widely used in a range of climate problems, but the assumption of constant diffusivity in the parameterization of the moist static energy (MSE) flux can be hardly justified. We demonstrate in this study that the diffusive MSE flux can be derived from the basic energy balance equation with a few tolerable assumptions. The estimated diffusivity is both spatially and seasonally dependent, and its midlatitude average is then tested against several scaling theories for the midlatitude eddy diffusivity. The result supports the diffusivity theory of Held and Larichev (1996) modified for the moist atmosphere, affording a dynamics-based parameterization of MSE diffusivity. The implementation of the parameterization in an EBM leads to an interactive MSE diffusivity that accounts for the midlatitude eddy response to climate forcing perturbations. Under a uniform radiative forcing, the EBM with a diffusivity so parameterized produces a weakening of the midlatitude diffusivity and a modestly polar-amplified surface temperature response as an inevitable outcome under the dual constraints of the nonlinear Clausius-Clapeyron relation and the temperature gradient-dependent diffusivity, even in the absence of any poleward amplifying radiative feedbacks. As the consequence of more isothermal temperature and reduced diffusivity, the variance of the midlatitude surface temperature also decreases with warming.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗

DOE OSTI2020

Hardening in Au-Ag nanoboxes from stacking fault-dislocation interactions

Porous, nano-architected metals with dimensions down to ~10 nm are predicted to have extraordinarily high strength and stiffness per weight, but have been challenging to fabricate and test experimentally. Here, we use colloidal synthesis to make ~140 nm length and ~15 nm wall thickness hollow Au-Ag nanoboxes with smooth and rough surfaces. In situ scanning electron microscope and transmission electron microscope testing of the smooth and rough nanoboxes show them to yield at 130 ± 45 MPa and 96 ± 31 MPa respectively, with significant strain hardening. A higher strain hardening rate is seen in rough nanoboxes than smooth nanoboxes. Finite element modeling is used to show that the structure of the nanoboxes is not responsible for the hardening behavior suggesting that material mechanisms are the source of observed hardening. Molecular dynamics simulations indicate that hardening is a result of interactions between dislocations and the associated increase in dislocation density.

42 ENGINEERING↗

DOE OSTI2025

Integration of a grey-box refrigerated case model in EnergyPlus via Python plugin

Commercial buildings, in particular grocery stores (due mainly to their large refrigeration load), provide opportunities for energy cost reductions. Grocery stores could offer substantial load flexibility to the power grid through participation in demand response programs because of their usage patterns and relatively high energy intensity. This load flexibility could come from modifying the control of heating, ventilation, and air conditioning (HVAC) systems, refrigeration systems, or both. Although estimation of the HVAC system’s load flexibility potential is relatively targeted in the literature, estimating load flexibility of refrigeration systems is nascent and has been a challenge, in part because of the lack of proper simulation tools that capture the dynamics in the refrigeration cases. The existing refrigerated case model within EnergyPlus, a whole building energy simulation program, assumes a constant case temperature throughout the simulation period and does not explicitly model the cycling of the compressor serving the refrigerated case. In addition, it does not encompass modeling of temperatures of the product inside the refrigerated case. This difference between modeled and actual operation can be a barrier to the development of demand control algorithm and accurate analysis of load flexibility potential. In this paper, we present a grey-box model for modeling refrigerated cases in grocery stores, which include medium temperature and low temperature. Four cases are modeled; two are low-temperature closed cases and two are medium-temperature cases with one closed and one open. Data from an experimental facility are used to train and test the models. Results demonstrate the efficacy of the grey-box models in predicting the temperatures. This model is integrated into EnergyPlus to capture the dynamic effects of case temperature on the environment and enhance the calculation of sensible and latent heat exchange with the environment (case credits). These enhancements can be leveraged more broadly to model advanced refrigeration controls such as defrost, develop and test unique algorithms that could affect refrigeration interactions with HVAC, and refine store design for any commercial building with refrigeration.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗