Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 37 records · Page 2

Improving Frequency Stability and Minimizing Load Shedding Events by Adopting Grid-Scale Energy Storage with Grid Forming Inverters

The upward adoption trend of renewable generation not only means cleaner energy integrated into modern power grids, but also that most new generation sources are based on front-end inverter bridges, used as interfaces to most wind generation and all the solar PV. It is well known that due to their power electronics-based construction rather than rotational shafts, these sources do not provide inertia inherently, nor substantial amounts of short-circuit currents. However, stable energy such as what can be stored in energy storage systems, although interfaced via inverters, can be controlled to respond to system disturbances in a manner that emulates inertial behavior. This paper focuses on the application of such energy storage systems to augment inertia in the island of Puerto Rico. To do so, a user defined inverter model that contains grid forming capabilities and fast frequency response is modeled and integrated into the real transmission system in power flow and dynamics software. Energy storage is then connected to two selected areas so that it not only provides frequency regulation to avoid widespread load shedding events, but also other tangible benefits. The simulated cases suggest that even relatively small energy storage systems can avert load shedding events if adequately placed in the transmission network.

Grid-forming inverters, IBR, Inertia↗

Creation of Simulated Test Cases for the Oscillation Source Location Contest

This paper describes in detail the design process and creation of 13 oscillatory test cases for the IEEE and NASPI co-hosted Oscillation Source Locating Contest conducted in 2021. The challenges behind the 13 cases are fully demystified. Based on the philosophy, implementation considerations and techniques used for the case design presented in this paper, interested readers can create additional interesting cases for testing the efficiency of their oscillation source location methods.

ENERGY PLANNING, POLICY, AND ECONOMY↗

Creation of Simulated Test Cases for the Oscillation Source Location Contest: Preprint

This paper describes in detail the design process and creation of 13 oscillatory test cases for the IEEE and NASPI co-hosted Oscillation Source Locating Contest conducted in 2021. The challenges behind the 13 cases are fully demystified. Based on the philosophy, implementation considerations and techniques used for the case design presented in this paper, interested readers can create additional interesting cases for testing the efficiency of their oscillation source location methods.

ENERGY PLANNING, POLICY, AND ECONOMY↗

Detection of DoS Attacks Using ARFIMA Modeling of GOOSE Communication in IEC 61850 Substations

Integration of Information and Communication Technology (ICT) in modern smart grids (SGs) offers many advantages including the use of renewables and an effective way to protect, control and monitor the energy transmission and distribution. To reach an optimal operation of future energy systems, availability, integrity and confidentiality of data should be guaranteed. Research on the cyber-physical security of electrical substations based on IEC 61850 is still at an early stage. In the present work, we first model the network traffic data in electrical substations, then, we present a statistical Anomaly Detection (AD) method to detect Denial of Service (DoS) attacks against the Generic Object Oriented Substation Event (GOOSE) network communication. According to interpretations on the self-similarity and the Long-Range Dependency (LRD) of the data, an Auto-Regressive Fractionally Integrated Moving Average (ARFIMA) model was shown to describe well the GOOSE communication in the substation process network. Based on this ARFIMA-model and in view of cyber-physical security, an effective model-based AD method is developed and analyzed. Two variants of the statistical AD considering statistical hypothesis testing based on the Generalized Likelihood Ratio Test (GLRT) and the cumulative sum (CUSUM) are presented to detect flooding attacks that might affect the availability of the data. Our work presents a novel AD method, with two different variants, tailored to the specific features of the GOOSE traffic in IEC 61850 substations. The statistical AD is capable of detecting anomalies at unknown change times under the realistic assumption of unknown model parameters. The performance of both variants of the AD method is validated and assessed using data collected from a simulation case study. We perform several Monte-Carlo simulations under different noise variances. The detection delay is provided for each detector and it represents the number of discrete time samples after which an anomaly is detected. In fact, our statistical AD method with both variants (CUSUM and GLRT) has around half the false positive rate and a smaller detection delay when compared with two of the closest works found in the literature. Our AD approach based on the GLRT detector has the smallest false positive rate among all considered approaches. Whereas, our AD approach based on the CUSUM test has the lowest false negative rate thus the best detection rate. Depending on the requirements as well as the costs of false alarms or missed anomalies, both variants of our statistical detection method can be used and are further analyzed using composite detection metrics.

IEC 61850 electrical substations↗

Effects of grain size and small-scale bedform architecture on CO 2 saturation from buoyancy-driven flow

Small-scale (mm-dm scale) heterogeneity has been shown to significantly impact CO 2 migration and trapping. To investigate how and why different aspects of small-scale heterogeneity affect the amount of capillary trapping during buoyancy-driven upward migration of CO 2 , we conducted modified invasion percolation simulations on heterogeneous domains. Realistic simulation domains are constructed by varying two important aspects of small-scale geologic heterogeneity: sedimentary bedform architecture and grain size contrast between the matrix and the laminae facies. Buoyancy-driven flow simulation runs cover 59 bedform architecture and 40 grain size contrast cases. Simulation results show that the domain effective CO 2 saturation is strongly affected by both grain size and bedform architecture. At high grain size contrasts, bedforms with continuous ripple lamination at the cm scale tend to retain higher CO 2 saturation than bedforms with discontinuous or cross lamination. In addition, the “extremely well sorted” grain sorting cases tend to have lower CO 2 saturation than expected for cross-laminated domains. Finally, both a denser CO 2 phase and greater interfacial tension increase CO 2 saturation. Again, variation in fluid properties seems to have a greater effect on CO 2 saturation for cross-laminated domains. This result suggests that differences in bedform architecture can impact how CO 2 saturation values respond to other variables such as grain sorting and fluid properties.

58 GEOSCIENCES↗

Numerical simulation of asteroid geometry variance on airburst threat

For an atmospheric airburst the primary source of concern when assessing uncertainty is the size and velocity. Determining these properties provides the basis for threat assessment, as the total energy of the asteroid may then be estimated, and the threat investigated thoroughly. Even with clarity as to how much energy an asteroid may deposit, a great deal of uncertainty still exists for the actual energy deposition process. One such source of uncertainty is the geometry of the incoming asteroid. The geometry of an asteroid will alter the stress distribution during entry, which adds uncertainty to when fracture will occur. Here, in this study, we use Smoothed Particle Hydrodynamics to model the atmospheric airburst of Tunguska-scale asteroids with varying geometric profiles, including a sphere, ellipsoid, binary and superellipsoid. Each asteroid is modeled as a homogenous structure with strength. We assess uncertainty through a series of planar 2D simulation cases for each geometry, comparing the source of stochasticity across geometries. A single 3D airburst simulation for each geometry is also analyzed. Additionally, the 3D cases are compared to the highly uncertain Tunguska event, predicting variance in burst height across geometries, but all bounded by theoretical burst heights proposed for Tunguska.

Airburst↗

On the applicability of various levels of detail for occupant behavior representation and modeling in building performance simulation

Occupant behavior (OB) is one of the significant sources of uncertainty in building performance simulation. While OB modeling has received increased attention in the past decade, research on the degree of granularity or level of detail (LoD) required for representing occupants is still in the nascent stages. This paper analyzes the modeling and applicability of three LoDs to represent occupants in building performance assessment. A medium-sized prototype office building located in Chicago, Illinois is used as the simulation case study. Ten occupant-centric attributes are adopted to develop the LoDs for OB representation. We first demonstrate the different modeling approaches required for simulating the three fidelity levels. Later, we illustrate the suitability of the developed LoDs in supporting six building performance use cases across different lifecycle stages. Furthermore, this study intends to provide guidance for the building simulation community on appropriate OB representation to support various use cases.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Progress in the validation of rotor aerodynamic codes using field data

Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.

17 WIND ENERGY↗

Variational fast forwarding for quantum simulation beyond the coherence time

Abstract Trotterization-based, iterative approaches to quantum simulation (QS) are restricted to simulation times less than the coherence time of the quantum computer (QC), which limits their utility in the near term. Here, we present a hybrid quantum-classical algorithm, called variational fast forwarding (VFF), for decreasing the quantum circuit depth of QSs. VFF seeks an approximate diagonalization of a short-time simulation to enable longer-time simulations using a constant number of gates. Our error analysis provides two results: (1) the simulation error of VFF scales at worst linearly in the fast-forwarded simulation time, and (2) our cost function’s operational meaning as an upper bound on average-case simulation error provides a natural termination condition for VFF. We implement VFF for the Hubbard, Ising, and Heisenberg models on a simulator. In addition, we implement VFF on Rigetti’s QC to demonstrate simulation beyond the coherence time. Finally, we show how to estimate energy eigenvalues using VFF.

97 MATHEMATICS AND COMPUTING↗

risingMicrobubbleLattice

This code is a simulation case to be run with OpenFOAM, an open-source computational fluid dynamics software. A gmsh mesh file is also included. Specifically, this simulation demonstrates the transport of a single microbubble rising through an ordered lattice due to an applied flow. The bubble deforms as it squeezes through the pores of the lattice. Is

Guo, Jack [Lawrence Livermore National Laboratory ↗

Control coordination in inverter‐based microgrids using AoI‐based 5G schedulers

Abstract A coordinated set point automatic adjustment with correction enabled (C‐SPAACE) framework that uses 5G communication for real‐time control coordination between inverter‐based resources (IBR) in microgrids is proposed. Utilising slicing capability, 5G offers low‐latency communication to C‐SPAACE under normal conditions. However, given the multitude of power grid use cases, a certain 5G slice for C‐SPAACE may have access only to limited radio spectrum resources, which if not managed well, greatly undermines the communication needs of C‐SPAACE framework. Thus, optimally scheduling the available spectrum resources among IBRs in a sliced 5G network‐based C‐SPAACE framework becomes a critical problem. To address this issue, the authors utilise a novel age of information (AoI) metric and designs an AoI‐based 5G scheduler to provide low‐latency communication to C‐SPAACE. Following this, a co‐simulation environment is designed using PSCAD/EMTDC and Python to simulate a microgrid supported by 5G communication. Time‐domain simulation case studies are performed using the proposed co‐simulation environment to evaluate the performance of C‐SPAACE using 5G with both AoI‐based and other baseline (non‐AoI) schedulers.

Choudhury, Biplav↗

Use of an Alternative Conceptual Model of Vadose Zone Heterogeneity to Evaluate Past Tank Leaks and Other Unplanned Releases within a Tank Farm at the Hanford Site - 20061

Washington State Department of Ecology (Ecology) requested that U.S. Department of Energy (DOE)-Office of River Protection (ORP) consider an evaluation of effects of fine-grained thin sediment layers on transport with a separate alternative conceptual model in its evaluation of the potential impact of tank leaks and other unplanned releases within the tank farm at the Hanford Site. For this alternative model, Ecology recommended that the model be developed based on the general framework of fine-grained units identified by a stakeholder group in their interpretations of variability in moisture content data collected in the vicinity of Waste Management Area (WMA) C. In discussions with DOE-ORP, Ecology acknowledged that the underlying data and interpretations of the occurrence and lateral continuity of the fine-grained thin layers identified by a stakeholder group are uncertain. However, DOE-ORP agreed to Ecology's recommendation and has provided support for the requested evaluation that involved development of an alternative model based on the general framework of the unpublished report by the stakeholder. DOE-ORP considered this evaluation to be a hypothetical evaluation of vadose zone heterogeneities at WMA C. General observations from the range of simulation cases examined in the evaluation of the effects of hypothetical vadose zone heterogeneities at WMA C are as follows. - The movement of the center of mass of the simulated plumes was generally vertically downward below the source for all simulations, including those that incorporated the hypothetical fine-grained units. - All simulations that incorporated hypothetical heterogeneity produced additional plume spreading over what was produced in simulations using Equivalent Homogeneous Media (EHM) model(s)a. The spreading resulted in a broadening of the fringes of the plume, resulting in a wider region of low concentration, but lower peak concentrations associated with the center of mass of the plume. - Simulations that used the silty-sand hydraulic properties, suggested by Ecology for the hypothetical fine-grained units, generally produced similar spreading and slightly lower peak mass flux at the water table when compared to the EHM modeling results. a The EHM-based models used at WMA C do not explicitly include small-scale fine-grained heterogeneities used in the model advocated by Ecology. - Simulations that used the silty-sand hydraulic properties, suggested by Ecology for the hypothetical fine-grained units, generally produced less spreading and an earlier arrival of mass flux at the water table when compared to the use of another set of hydraulic properties from a silty-sand sample collected at a nearby disposal facility. - The EHM representation of the vadose zone generally produced higher peak mass flux and an earlier occurrence of peak fluxes at the water table compared to all analyses incorporating additional hypothetical heterogeneity. Results of this alternative model evaluation of past leaks provided some insight into the transport effects of vadose zone heterogeneities at WMA C that helped resolve Ecology's comments and issues related to the effects of vadose zone heterogeneities on past leaks and losses from the WMA C tank farms area. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Machine learning-enhanced model-based scenario optimization for DIII-D

Abstract Scenario development in tokamaks is an open area of investigation that can be approached in a variety of different ways. Experimental trial and error has been the traditional method, but this required a massive amount of experimental time and resources. As high fidelity predictive models have become available, offline development and testing of proposed scenarios has become an option to reduce the required experimental resources. The use of predictive models also offers the possibility of using a numerical optimization process to find the controllable inputs that most closely achieve the desired plasma state. However, this type of optimization can require as many as hundreds or thousands of predictive simulation cases to converge to a solution; many of the commonly used high fidelity models have high computational burdens, so it is only reasonable to run a handful of predictive simulations. In order to make use of numerical optimization approaches, a compromise needs to be found between model fidelity and computational burden. This compromise can be achieved using neural networks surrogates of high fidelity models that retain nearly the same level of accuracy as the models they are trained to replicate while reducing the computation time by orders of magnitude. In this work, a model-based numerical optimization tool for scenario development is described. The predictive model used by the optimizer includes neural network surrogate models integrated into the fast Control-Oriented Transport simulation framework. This optimization scheme is able to converge to the optimal values of the controllable inputs that produce the target plasma scenario by running thousands of predictive simulations in under an hour without sacrificing too much prediction accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of the Transformer Magnetizing Inductance on the Performance of the Dual-Active Bridge Converter

The magnetizing inductance of the high-frequency transformer in the dual-active bridge (DAB) converter is usually assumed to be much larger than its leakage inductance, thereby the magnetizing current can be negligible. However, ignoring the magnetizing current will lead to an error in the DAB converter design process. This paper analyzes the operation and performance of the DAB converter for a general case in which the transformer magnetizing inductance is taken into consideration. The analysis shows that the magnetizing inductance impacts the soft-switching region, the power transfer capability, and the leakage inductor current. A simulation case study is provided to verify the analysis. A DAB converter prototype is built to experimentally validate the analysis as well.

Yin, Shiyuan↗

A subgrid model for electron-scale turbulent transport in global ion-scale gyrokinetic simulations of tokamak plasmas

A subgrid electron–temperature–gradient (ETG) model is demonstrated here which averages local electron-scale turbulence over intermediate scales in space and time to include in global ion–temperature–gradient (ITG) simulations. This approach results in ion-scale equations, which incorporate the electron heat transport from ETG turbulence and the effects of electron-scale turbulence on the ion scale. Flux-tube ETG Cyclone Base Case simulations are performed at different radial locations and a kinetic form of the flux is added to global ion-scale (IS) simulations as a source term. Analytic radial profiles of ETG heat flux are constructed and compared to flux-tube results at multiple radial locations. Different ratios of ITG to ETG heat flux levels are considered and the results of capturing ETG heat transport in global ion-scale simulations are discussed. The potential coupling of the ETG streamer potential and intermediate-scale zonal flows to the ion scale is further addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Surrogate Model of Flexible Research Platform EnergyPlus Models to Enable Sensitivity Analysis

This letter report describes the surrogate models developed from the EnergyPlus model of Oak Ridge National Laboratory’s Flexible Research Platform. Two data-driven black-box models were developed, and the outputs of the surrogate models were compared with the EnergyPlus model. The two models developed are a multilayer perceptron deep learning model, and a long short-term memory (LSTM) neural network model. The three factors for selecting the black-box models are scalability, computation time, and accuracy. A total of 107 input variables were the dominant variables in determining the outputs of building energy consumptions and thermal comfort. A total of 54 output variables were identified as the prediction targets, including the system- and zone-level outputs. The large set of the simulation cases were generated by integrating sensor errors into an emulator based on EnergyPlus and Python EMS, which includes advanced control sequences from ASHRAE Guideline 36-2018: High-Performance Sequences of Operation. The surrogate models were developed based on a set of large-scale simulation runs (i.e., 4,000 runs) on a cloud platform. The comparison analysis shows that the two black-box models had good accuracy for predicting new outputs for sensitivity analysis using the root mean square error metric. As a next step, the developed surrogate models will be used to perform sensitivity analysis for different sensor impacts (e.g., sensor types, sensor locations).

42 ENGINEERING↗

Flexible Service Contracting for Risk Management within Integrated Transmission and Distribution Systems

The general objective of our project has been to investigate the ability of Independent Distribution System Operators (IDSOs), functioning as linkage agents for Integrated Transmission and Distribution (ITD) systems, to facilitate the flexible availability and usage of reserve in support of ITD system operations. This general objective is in accordance with Order 2222 of the U.S. Federal Energy Regulatory Commission (FERC), titled “Participation of Distributed Energy Resource Aggregations in Markets Operated by Regional Transmission Organizations and Independent System Operators”. The Final Rule for FERC Order 2222 was issued on September 17, 2020. The primary contribution of our project is that we have formulated an innovative energy management Transactive Energy Design (TES) approach for ITD systems that provides promising support for our project objective. Specifically, we have developed a new type of IDSO-managed TES design for distribution system operations, as well as new types of swing contracts permitting IDSOs to participate in transmission system operations as providers of reserve harnessed from distribution system resources in return for appropriate compensation. Together, these design elements constitute a scalable market-based approach facilitating efficient reserve procurement for ITD system operations from a fuller range of power resources. The efficacy of our approach has been demonstrated by means of detailed conceptual analyses as well as test case simulations. To implement the latter, we have developed the ITD TES Platform V2.0, a computational platform that permits the modeling and software implementation of ITD system operations over time. The Electric Reliability Council of Texas (ERCOT) energy region has been used as the empirical anchor for the development of this platform. Our conceptual and test-case work has been reported in a Wiley/IEEE Press book, two refereed book chapters, and seven refereed journal articles. The key components comprising the ITD TES Platform V2.0 have been released as documented open-source software at online GitHub repositories.

24 POWER TRANSMISSION AND DISTRIBUTION↗