Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid simulation”

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 739 records · Page 41

Simulations of Cryogenic Line Chilldown with Advanced Sub-Grid Wall Boiling Models

A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.

Computational Fluid Dynamics↗

Simulations of Cryogenic Line Chilldown with Advanced Sub-Grid Wall Boiling Models

A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.

Computational Fluid Dynamics↗

A Parameterization of Sub-grid Topographical Effects on Solar Radiation in the E3SM Land Model (Version 1.0): Implementation and Evaluation Over the Tibetan Plateau

Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth System Models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to evaluate the performance of ELM. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover and surface temperature especially in the high-elevation and snow-cover regions over the TP. Topography has non-negligible effects on surface energy budget, snow cover, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20 W/m 2 , 10 W/m 2 , and 5 W/m 2 , respectively. The differences in land surface albedo, snow cover fraction, and surface temperature between TOP and PP exceed 0.1, 20%, and 1K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations, and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are more significant with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

58 GEOSCIENCES↗

Cause and Cure - Deterioration in Accuracy of CFD Simulations with Use of High-Aspect-Ratio Triangular/Tetrahedral Grids

Traditionally high-aspect ratio triangular/tetrahedral meshes are avoided by CFD researchers in the vicinity of a solid wall, as it is known to reduce the accuracy of gradient computations in those regions. Although for certain complex geometries, the use of high-aspect ratio triangular/tetrahedral elements in the vicinity of a solid wall can be replaced by quadrilateral/prismatic elements, ability to use triangular/tetrahedral elements in such regions without any degradation in accuracy can be beneficial from a mesh generation point of view. The benefits also carry over to numerical frameworks such as the space-time conservation element and solution element (CESE), where simplex elements are the mandatory building blocks. With the requirement of the CESE method in mind, a rigorous mathematical framework that clearly identifies the reason behind the difficulties in use of such high-aspect ratio simplex elements is formulated using two different approaches and presented here. Drawing insights from the analysis, a potential solution to avoid that pitfall is also provided as part of this work. Furthermore, through the use of numerical simulations of practical viscous problems involving high-Reynolds number flows, how the gradient evaluation procedures of the CESE framework can be effectively used to produce accurate and stable results on such high-aspect ratio simplex meshes is also showcased.

CESE Method↗

Range Safety Flight Elevation Limit Calculation Software

This program was developed to fill a need within the Wallops Flight Facility workflow for automation of the development of vertical plan limit lines used by flight safety officers during the conduct of expendable launch vehicle missions. Vertical plane present-position-based destruct lines have been used by range safety organizations at numerous launch ranges to mitigate launch vehicle risks during the early phase of flight. Various ranges have implemented data submittal and processing workflows to develop these destruct lines. As such, there is significant prior art in this field. The ElLimits program was developed at NASA's Wallops Flight Facility to automate the process for developing vertical plane limit lines using current computing technologies. The ElLimits program is used to configure launch-phase range safety flight control lines for guided missiles. The name of the program derives itself from the fundamental quantity that is computed - flight elevation limits. The user specifies the extent and resolution of a grid in the vertical plane oriented along the launch azimuth. At each grid point, the program computes the maximum velocity vector flight elevation that can be permitted without endangering a specified back-range location. Vertical plane x-y limit lines that can be utilized on a present position display are derived from the flight elevation limit data by numerically propagating 'streamlines' through the grid. The failure turn and debris propagation simulation technique used by the application is common to all of its analysis options. A simulation is initialized at a vertical plane grid point chosen by the program. A powered flight failure turn is then propagated in the plane for the duration of the so-called RSO reaction time. At the end of the turn, a delta-velocity is imparted, and a ballistic trajectory is propagated to impact. While the program possesses capability for powered flight failure turn modeling, it does not require extensive user inputs of vehicle characteristics (e.g., thrust and aerodynamic data), nor does it require reams of turn data after the traditional fashion of the Air Force ranges. The program requires a nominal trajectory table (time, altitude, range, velocity, and flight elevation) and makes heavy use of it to initialize and model a failure turn.

Lanzi, Raymond J↗

Numerical Behaviour of a Smooth Local Correlation-based Transition Model in a Newton-Krylov Flow Solver

The numerical behaviour of transport-equation-based transition models, including both iterative and grid convergence, is influenced by the source terms. Transition models contain source terms that are large and highly nonlinear, and can be destabilizing in a strong implicit solver. Linearization strategies with varying levels of coupling are evaluated in conjunction with a source-term time step restriction to determine best-practices for solving the SA-sLM2015smooth local correlation-based transition model in an implicit Newton-Krylov flow solver. Achieving deep iterative convergence facilitates a detailed investigation of the grid convergence of these free-transition simulations, which are evaluated relative to fully-turbulent simulations performed using the Spalart-Allmaras turbulence model. Simulations of the NLF0416 general aviation airfoil, VA-2 supercritical airfoil, and NASA CRM-NLF wing-body geometry are performed over a range of grid levels. The results demonstrate that both a fully-coupled linearization strategy and a source-term time step restriction improve nonlinear convergence as the complexity of the free-transition simulations increases. In general, additional grid resolution is required for free-transition simulations relative to fully-turbulent simulations in order to achieve a similar level of accuracy, with the grid convergence of free-transition simulations sensitive to the streamwise grid spacings in the transition regions.

AATT↗

Commercial Building Planning and Retrofitting Strategy for Grid Services

The increasing integration of distributed energy resources (DERs) plays an important role in improving energy consumption efficiency. In September 2020, the Federal Energy Regulatory Commission (FERC) approved Order 2222 which opens wholesale electricity markets to small capacity DERs. The benefit of this new FERC Order 2222 is that DERs, such as rooftop solar panels and batteries, will be able to participate in regional electricity markets and provide grid services. Meanwhile, the planning and operation strategies of DERs are facing new challenges to account for the impact of the wholesale market with numerous uncertainty factors. Therefore, in this paper, we propose a new planning and retrofitting model for long-term commercial buildings that considers both DER investment and market participation. Specifically, we explore the capability of implementing DERs for grid services. The effectiveness of the proposed model is validated using real-world data. Simulation results also validate that participating in grid services can significantly increase revenues through appropriate building energy management and shorten the payback period of DER investments.

building energy management↗

Wave Resource Characterization Using an Unstructured Grid Modeling Approach

This paper presents a modeling study conducted on the central Oregon coast for wave resource characterization using the unstructured-grid SWAN model coupled with a nested-grid WWIII model. The flexibility of models of various spatial resolutions and the effects of open- boundary conditions simulated by a nested-grid WWIII model with different physics packages were evaluated. The model results demonstrate the advantage of the unstructured-grid modeling approach for flexible model resolution and good model skills in simulating the six wave resource parameters recommended by the International Electrotechnical Commission in comparison to the observed data in Year 2009 at National Data Buoy Center Buoy 46050. Notably, spectral analysis indicates that the ST4 physics package improves upon the model skill of the ST2 physics package for predicting wave power density for large waves, which is important for wave resource assessment, device load calculation, and risk management. In addition, bivariate distributions show the simulated sea state of maximum occurrence with the ST4 physics package matched the observed data better than that with the ST2 physics package. This study demonstrated that the unstructured-grid wave modeling approach, driven by the nested-grid regional WWIII outputs with the ST4 physics package, can efficiently provide accurate wave hindcasts to support wave resource characterization. Our study also suggests that wind effects need to be considered if the dimension of the model domain is greater than approximately 100 km, or O (10^2 km).

Wu, Wei-Cheng↗

Stabilizing Inverter-Based Transmission Systems: Power Hardware-in-the-Loop Experiments with a Megawatt-Scale Grid-Forming Inverter

This article presents what the authors believe to be the first experimental verification of the ability of grid-forming (GFM) inverters to stabilize a transmission electric power system that is otherwise unstable. The experiments described here were performed using power hardware-in-the-loop (PHIL) simulation to connect a megawatt-scale battery inverter to a real-time electromagnetic transient (EMT) simulation of the near-future Maui power system. This allows the dynamic interactions between the inverter and the power system to be observed without putting the real power system at risk. The ability to use the actual inverter hardware removes the need to rely on a computer model approximation of the inverter's behavior.

electromagnetic transient simulations↗

Incentive-Based Control and Coordination of Distributed Energy Resources. Final Report

The rapid growth of low-inertia renewable energy resources represents an immense opportunity for the U.S. to minimize its carbon footprint, while presenting a challenge for system operators as traditional “spinning” generation resources are displaced. This transformation requires solutions to robustly and cost-effectively manage dynamic changes on the grid while ensuring quality of service. There is growing recognition that distributed energy resources (DERs – loads, distributed generation, storage, electric vehicles, etc.) represent a great potential to perform this function. This report details a scalable hierarchical control strategy for coordinating thousands of heterogeneous DERs in order to provide three ancillary gird services: 1) frequency response, a rapid (within 2 seconds) service used to maintain the grid at 60 Hz; 2) regulation response, a larger response used to balance on-going differences in supply and demand on the grid; and, 3) ramping reserve, used to bring on capacity to replace lost generation. The developed control architecture consists of 1) a Distribution Reliability Coordinator (DRC), 2) aggregator controllers, and finally, 3) DER controllers. The DRC coordinates all the flexible loads and distributed energy resources on the distribution system and bids the aggregate flexibility into the wholesale market. The DRC also performs resource allocation (deciding how the flexible devices on the distribution system will be dispatched) to meet the requirements for each of the three ancillary services. The aggregators control a collection of devices that are geographically co-located and provide a platform to aggregate their flexibility and mitigate uncertainty associated with a single flexible load by managing a large aggregation of similar loads (so that any stochastic effects average out). Finally, DER controllers (for example, thermostats for air conditioners, chargers for electric vehicles) adjust the power consumption of individual devices to ensure accurate and reliable delivery of the requested grid services. Evaluation and validation of this approach was conducted with a mixture of simulation and hardware testing for a range of DERs in residential (air conditioners and water heaters), commercial HVAC (ventilation fans and central chiller plants), and micro-grids (battery energy storage, solar PV etc.). This culminated in DER hardware and simulation models operating in conjunction with grid simulations as part of a large scale federated test-bed demonstration. Simulations of over 10,000 controllable residential devices (a combination of air conditioners, water heaters, batteries, and electric vehicles) demonstrated that this advanced control approach could simultaneously meet all three grid services. Furthermore, hardware-in-the-loop testing was performed with operational micro-grid and commercial building HVAC equipment. The micro-grid equipment (representing PV, battery storage, and back-up generators) demonstrated the ability to meet all frequency response requirements (except ramping time). The commercial building HVAC ventilation equipment successfully provided the regulating response service. An Industry Advisory Board (IAB) was formed and engaged throughout the project to get feed-back on project risks and potential commercialization barriers. Their input was invaluable in identify key risks stemming from regulatory and market access challenges, financial risks from under-performance, and competitive threats from utility-scale solutions such as integrated storage. A phased commercialization plan was developed that addresses continued risk reduction (value proposition refinement and performance uncertainty assessments), piloting through field demonstrations with key stakeholders, and finally, deployment as an emerging standard for DER control.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The importance of considering sub-grid cloud variability when using satellite observations to evaluate the cloud and precipitation simulations in climate models

Abstract. Satellite cloud observations have become an indispensable tool for evaluating general circulation models (GCMs). To facilitate the satellite and GCM comparisons, the CFMIP (Cloud Feedback Model Inter-comparison Project) Observation Simulator Package (COSP) has been developed and is now increasingly used in GCM evaluations. Real-world clouds and precipitation can have significant sub-grid variations, which, however, are often ignored or oversimplified in the COSP simulation. In this study, we use COSP cloud simulations from the Super-Parameterized Community Atmosphere Model (SPCAM5) and satellite observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and CloudSat to demonstrate the importance of considering the sub-grid variability of cloud and precipitation when using the COSP to evaluate GCM simulations. We carry out two sensitivity tests: SPCAM5 COSP and SPCAM5-Homogeneous COSP. In the SPCAM5 COSP run, the sub-grid cloud and precipitation properties from the embedded cloud-resolving model (CRM) of SPCAM5 are used to drive the COSP simulation, while in the SPCAM5-Homogeneous COSP run only grid-mean cloud and precipitation properties (i.e., no sub-grid variations) are given to the COSP. We find that the warm rain signatures in the SPCAM5 COSP run agree with the MODIS and CloudSat observations quite well. In contrast, the SPCAM5-Homogeneous COSP run which ignores the sub-grid cloud variations substantially overestimates the radar reflectivity and probability of precipitation compared to the satellite observations, as well as the results from the SPCAM5 COSP run. The significant differences between the two COSP runs demonstrate that it is important to take into account the sub-grid variations of cloud and precipitation when using COSP to evaluate the GCM to avoid confusing and misleading results.

54 ENVIRONMENTAL SCIENCES↗

Weak influence of anthropogenic emissions on aerosol, cloud, and rain in the wet season of the Amazon rainforest

Anthropogenic emissions have been shown to affect new particle formation, aerosol concentrations, and clouds. Such effects vary with region, environmental conditions, and cloud types. In the wet season of Amazonia, anthropogenic emissions emitted from Manaus, Brazil, can significantly increase the cloud condensation nuclei (CCN) concentrations compared to the background of mainly natural aerosols. However, the regional response of cloud and rain to anthropogenic emissions in Amazonia remains very uncertain. Here, we aim to quantify how aerosol concentration, cloud, and rain respond to changes in anthropogenic emissions through parameterized new particle formation and primary aerosol emission in the Manaus region and to understand the underlying mechanisms. We ran the atmosphere-only configuration of the HadGEM3 climate model with a nested regional domain that covers most of the rainforest region (720 km by 1200 km with 3 km resolution) under scaled regional emissions. The 7 d simulations show that, in the areas that are affected by anthropogenic emissions, when aerosol and precursor gas emissions are doubled from the baseline emission inventories, aerosol number concentrations increase by 13 %. The nucleation rate that involves sulfuric acid and biogenic compounds generally increases with pollution levels. However, nucleation is suppressed very close to the pollution source, resulting in lower nucleation and soluble Aitken mode aerosol number concentrations. We also found that doubling the anthropogenic emission can increase the cloud droplet number concentrations ( N d ) by 9 %, but cloud water and rain mass mixing ratios do not change significantly. Even very strong reductions in aerosol number concentrations by a factor of 4, which is an unrealistic condition, cause only a 4 % increase in rain over the domain. If we assume our simulation has a fine enough grid resolution and an accurate representation of the relevant atmospheric processes, the simulated weak and non-linear response of cloud and rain properties to linearly scaled anthropogenic emissions suggests that the interactions among aerosol, cloud, and precipitation in the Amazonian convective environment are buffered by microphysical processes. It also implies that the convective environment is resilient to the changes in Nd that occur in response to localized anthropogenic aerosol perturbations.

Wang, Xuemei [University of Leeds (United Kingdom)↗

A hierarchical stochastic modeling approach for representing point bar geometries and petrophysical property variations

The flow of fluids in point bars is affected by the existence of heterogeneities like shale drapes that are found on the surfaces of inclined heterolithic stratifications. In fact, these shale drapes can act as fluid flow baffles; therefore, developing a framework for modeling point bars and their associated heterogeneities is vital. In this study, a stochastic process-based modeling approach is presented for capturing the main point bar heterogeneities: accretion surfaces (i.e., the aerial heterogeneity) and inclined heterolithic stratifications (i.e., the vertical heterogeneity). The former was modeled using a sine-generation function and the latter, with a sigmoidal function after which they were combined into a 3D point bar model. To ensure proper modeling of petrophysical properties, we developed a more representative gridding scheme which generates curvilinear grids representative of the point bar geometry. This grid was then transformed into a rectilinear grid to allow for geostatistical simulation after which all petrophysical properties were mapped back into the original curvilinear grid. An essential element of this modeling approach is the stochastic representation of shale drapes at the interface between successive accretion surfaces. The workflow was tested using a real field dataset for the Cranfield field, Mississippi. The constructed model was then subjected to a flow simulation study mimicking a CO 2 storage scenario. Various sensitivities were simulated to evaluate the effect of heterogeneities on CO 2 flow within the point-bar. Results demonstrate the importance of representing point-bar related heterogeneity and the spatial distribution of shale drapes on CO 2 plume migration and storage.

58 GEOSCIENCES↗

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within ≲3σ of the mean of all simulations on a given grid (128 3 and 256 3 grid cells), where σ describes the statistical variation due to the flow’s time dependence. They also agree well with a 512 3 reference run. The 128 3 and 256 3 simulations agree within 9% and 4%, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

79 ASTRONOMY AND ASTROPHYSICS↗

Loci-STREAM Version 0.9

Loci-STREAM is an evolving computational fluid dynamics (CFD) software tool for simulating possibly chemically reacting, possibly unsteady flows in diverse settings, including rocket engines, turbomachines, oil refineries, etc. Loci-STREAM implements a pressure- based flow-solving algorithm that utilizes unstructured grids. (The benefit of low memory usage by pressure-based algorithms is well recognized by experts in the field.) The algorithm is robust for flows at all speeds from zero to hypersonic. The flexibility of arbitrary polyhedral grids enables accurate, efficient simulation of flows in complex geometries, including those of plume-impingement problems. The present version - Loci-STREAM version 0.9 - includes an interface with the Portable, Extensible Toolkit for Scientific Computation (PETSc) library for access to enhanced linear-equation-solving programs therein that accelerate convergence toward a solution. The name "Loci" reflects the creation of this software within the Loci computational framework, which was developed at Mississippi State University for the primary purpose of simplifying the writing of complex multidisciplinary application programs to run in distributed-memory computing environments including clusters of personal computers. Loci has been designed to relieve application programmers of the details of programming for distributed-memory computers.

Wright, Jeffrey↗

Simulation of Hawaiian Electric Companies Feeder Operations with Advanced Inverters and Analysis of Annual Photovoltaic Energy Curtailment

The Hawaiian Electric Companies achieved a consolidated Renewable Portfolio Standard (RPS) of approximately 26% at the end of 2016. This significant RPS performance was achieved using various renewable energy sources - biomass, geothermal, solar photovoltaic (PV) systems, hydro, wind, and biofuels - and customer-sited, grid-connected technologies (primarily private rooftop solar PV systems). The Hawaiian Electric Companies are preparing grid-modernization plans for the island grids. The plans outline specific near-term actions to accelerate the achievement of Hawai'i's 100% RPS by 2045. A key element of the Companies' grid-modernization strategy is to utilize new technologies - including storage and PV systems with grid-supportive inverters - that will help to more than triple the amount of private rooftop solar PV systems. The Hawaiian Electric Companies collaborated with the Smart Inverter Technical Working Group Hawai'i (SITWG) to partner with the U.S. Department of Energy's National Renewable Energy Laboratory (NREL) to research the implementation of advanced inverter grid support functions (GSF). Together with the technical guidance from the Companies's planning engineers and stakeholder input from the SITWG members, NREL proposed a scope of work that explored different modes of voltage-regulation GSF to better understand the trade-offs of the grid benefits and curtailment impacts from the activation of selected advanced inverter grid support functions. The simulation results presented in this report examine the effectiveness in regulating voltage as well as the impact to the utility and the customers of various inverter-based grid support functions on two Hawaiian Electric distribution substations.

14 SOLAR ENERGY↗

Benefits of Dual Fuel Heat Pump Grid-responsive Control: A Model-based Control Optimization Approach Using Building and Equipment Co-simulation

Conventional dual fuel heat pumps lack the intelligent control mechanisms to efficiently manage the switch between heat pump and furnace, leading to sub-optimal energy usage and, in some cases, increased operating costs. To resolve this gap, this study applies optimized control on hybrid heat pumps. With a focus on equipment control strategies, we compare the performances of five spacing heating equipment, including a conventional heat pump (HP), a conventional furnace, a dual fuel heat pump (DFHP) with conventional control, a dual fuel heat pump with smart control, and a novel seamlessly fuel flexible heat pump (SFFHP). While DFHP runs on either gas or electricity at any given moment, SFFHP concurrently consumes gas and electricity by continuously optimizing the proportion of each. In this research, a co-simulation framework is developed by integrating a building envelope model with a physics-based heat pump simulation model to analyze the benefits of grid-responsive controls of DFHP and SFFHP. The model-based optimal controls adjust the operation of the heat pump and gas furnace based on utility price signals and marginal grid emission to minimize utility cost and CO 2 emissions for multiple climate zones, different utility tariffs, and marginal grid emission scenarios. Case studies in Chicago and Los Angeles demonstrate that SFFHP and DFHP, with model-based optimal control, can deliver significant reductions in peak demand, utility cost, and CO 2 emission. In Chicago, SFFHP and smart controlled DFHP yield up to 64.7% and 61.7% utility cost reduction and up to 15.7% and 8.5% CO 2 emission reduction compared to the gas furnace. In Los Angeles, SFFHP and smart controlled DFHP achieve up to 43.6% and 40.1% utility cost reduction and up to 13.8% and 14.1% CO2 emission reduction compared to conventional heat pumps. In conclusion, by leveraging the fuel flexibility nature of dual fuel heat pumps, the model-based control optimization approach makes dual fuel heat pump an attractive option for demand response programs.

Control↗

Novel Real- Time EMT- TS Modeling Architecture for Feeder Blackstart Simulations

This paper presents the development and benchmarking of a novel real-time electromagnetic-transient and transient-stability (EMT-TS) modeling architecture for distribution feeder restoration studies. The work presents for the first time in the literature a real-time EMT-TS testbed in which the grid-forming unit is simulated in EMT domain, operating as the slack bus of the phasor domain while including unbalanced voltage conditions. To evaluate the performance and limitations of the proposed model, an equivalent EMT testbed is developed to serve as a benchmark. First, the co-simulation framework and the domain coupling methodology are introduced. Next, the steady state operation of the EMT-TS and EMT models are matched. Then, tests are conducted to analyze the transient performance of the proposed EMT-TS model when compared to its EMT benchmark. Results reveal that when utilizing a battery energy storage system (BESS) as the grid-forming unit, the EMT-TS testbed can maintain high accuracy for typical load steps.

Paduani, Victor↗