SEARCH · Engineering Papers

Results for “Dynamics”

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 1,603 records · Page 89

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Overcoming potential energy distortions in constrained internal coordinate molecular dynamics simulations

The Internal Coordinate Molecular Dynamics (ICMD) method is an attractive molecular dynamics (MD) method for studying the dynamics of bonded systems such as proteins and polymers. It offers a simple venue for coarsening the dynamics model of a system at multiple hierarchical levels. For example, large scale protein dynamics can be studied using torsional dynamics, where large domains or helical structures can be treated as rigid bodies and the loops connecting them as flexible torsions. ICMD with such a dynamic model of the protein, combined with enhanced conformational sampling method such as temperature replica exchange, allows the sampling of large scale domain motion involving high energy barrier transitions. Once these large scale conformational transitions are sampled, all-torsion, or even all-atom, MD simulations can be carried out for the low energy conformations sampled via coarse grained ICMD to calculate the energetics of distinct conformations. Such hierarchical MD simulations can be carried out with standard all-atom forcefields without the need for compromising on the accuracy of the forces. Using constraints to treat bond lengths and bond angles as rigid can, however, distort the potential energy landscape of the system and reduce the number of dihedral transitions as well as conformational sampling. We present here a two-part solution to overcome such distortions of the potential energy landscape with ICMD models. To alleviate the intrinsic distortion that stems from the reduced phase space in torsional MD, we use the Fixman compensating potential. To additionally alleviate the extrinsic distortion that arises from the coupling between the dihedral angles and bond angles within a force field, we propose a hybrid ICMD method that allows the selective relaxing of bond angles. This hybrid ICMD method bridges the gap between all-atom MD and torsional MD. We demonstrate with examples that these methods together offer a solution to eliminate the potential energy distortions encountered in constrained ICMD simulations of peptide molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Role of density modulation in the spatially resolved dynamics of strongly confined liquids

Confinement by walls usually produces a strong modulation in the density of dense liquids near the walls. Using molecular dynamics simulations, we examine the effects of the density modulation on the spatially resolved dynamics of a liquid confined between two parallel walls, using a resolution of a fraction of the interparticle distance in the liquid. The local dynamics is quantified by the relaxation time associated with the temporal autocorrelation function of the local density. We find that this local relaxation time varies in phase with the density modulation. The amplitude of the spatial modulation of the relaxation time can be quite large, depending on the characteristics of the wall and thermodynamic parameters of the liquid. To disentangle the effects of confinement and density modulation on the spatially resolved dynamics, we compare the dynamics of a confined liquid with that of an unconfined one in which a similar density modulation is induced by an external potential. We find several differences indicating that density modulation alone cannot account for all the features seen in the spatially resolved dynamics of confined liquids. We also examine how the dynamics near a wall depends on the separation between the two walls and show that the features seen in our simulations persist in the limit of large wall separation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Structure and dynamics of ethane confined in silica nanopores in the presence of CO 2

Fundamental understanding of the subcritical/supercritical behavior of key hydrocarbon species inside nano-porous matrices at elevated pressure and temperature is less developed compared to bulk fluids, but this knowledge is of great importance for chemical and energy engineering industries. This study explores in detail the structure and dynamics of ethane (C 2 H 6 ) fluid confined in silica nanopores, with a focus on the effects of pressure and different ratios of C 2 H 6 and CO 2 at non-ambient temperature. Quasi-elastic neutron scattering (QENS) experiments were carried out for the pure C 2 H 6 , C 2 H 6 :CO 2 = 3:1, and 1:3 mixed fluids confined in 4-nm cylindrical silica pores at three different pressures (30 bars, 65 bars, and 100 bars) at 323 K. Two Lorentzian functions were required to fit the spectra, corresponding to fast and slow translational motions. No localized motions (rotations and vibrations) were detected. Higher pressures resulted in hindrances of the diffusivity of C 2 H 6 molecules in all systems investigated. Pore size was found to be an important factor, i.e., the dynamics of confined C 2 H 6 is more restricted in smaller pores compared to the larger pores used in previous studies. Molecular dynamics simulations were performed to complement the QENS experiment at 65 bars, providing supportive structure information and comparable dynamic information. The simulations indicate that CO 2 molecules are more strongly attracted to the pore surface compared to C 2 H 6 . The C 2 H 6 molecules interacting with or near the pore surface form a dense first layer (L1) close to the pore surface and a second less dense layer (L2) extending into the pore center. Both the experiments and simulations revealed the role that CO 2 molecules play in enhancing C 2 H 6 diffusion (“molecular lubrication”) at high CO 2 :C 2 H 6 ratios. The energy scales of the two dynamic components, fast and slow, quantified by both techniques, are in very good agreement. Herein, the simulations identified the fast component as the main contributor to the dynamics. Molecule motions in the L2 region are mostly responsible for the dynamics (fast and slow) that can be detected by the instrument.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Method of information entropy for convergence assessment of molecular dynamics simulations

The lack of a reliable method to evaluate the convergence of molecular dynamics simulations has contributed to discrepancies in different areas of molecular dynamics. Here, the method of information entropy is introduced to molecular dynamics for stationarity assessment. The Shannon information entropy formalism is used to monitor the convergence of the atom motion to a steady state in a continuous spatial domain and is also used to assess the stationarity of calculated multidimensional fields such as the temperature field in a discrete spatial domain. It is demonstrated in this work that monitoring the information entropy of the atom position matrix provides a clear indicator of reaching steady state in radiation damage simulations, non-equilibrium molecular dynamics thermal conductivity computations, and simulations of Poiseuille and Couette flow in nanochannels. A main advantage of the present technique is that it is non-local and relies on fundamental quantities available in all molecular dynamics simulations. Unlike monitoring average temperature, the technique is applicable to simulations that conserve total energy such as reverse non-equilibrium molecular dynamics thermal conductivity computations and to simulations where energy dissipates through a boundary as in radiation damage simulations. The method is applied to simulations of iron using the Tersoff/ZBL splined potential, silicon using the Stillinger–Weber potential, and to Lennard–Jones fluid. Its applicability to both solids and fluids shows that the technique has potential for generalization to other areas in molecular dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2020

Dynamic PRA Methods to Evaluate the Impact on Accident Progression of Accident Tolerant Fuels

Accident tolerant fuels (ATFs) are new nuclear fuels developed in response to the accident at the Fukushima power station in March 2011. The goal of ATFs is to withstand accident scenarios through better performance compared to currently employed fuels (e.g., small-scale hydrogen generation). This paper targets a method for evaluating and comparing ATF performance from a probabilistic risk assessment (PRA) perspective by employing a newly developed combination of event trees and dynamic PRA methods. Compared to classical PRA methods based on event trees and fault trees, dynamic PRA can evaluate with higher resolution the safety impacts of physics dynamics and the timing/sequencing of events on the accident progression without the need to introduce overly conservative modeling assumptions and success criteria. In this paper, we analyze the impact on the accident progression of three different cladding configurations for two initiating events [a large break loss-of-coolant accident (LB-LOCA) and a station blackout (SBO)] by employing dynamic PRA methods. The goal is to compare the safety performance of ATFs (FeCrAl and Cr-coated cladding) and the currently employed Zr-based clad fuel. We employ two different strategies. The first focuses on the identification of success criteria discrepancies between the accident sequences generated by the classical PRA model and the set of simulation runs generated by dynamic PRA using ATF. The second one, on the other hand, directly uses dynamic PRA to evaluate the impact of timing of events (e.g., recovery actions) on accident progression. In conclusion, by applying these methods to the LB-LOCA and SBO initiating events, we show how dynamic PRA methods can provide analysts with detailed and quantitative information on the safety impact of ATFs.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Spin dynamics of the antiferromagnetic Heisenberg model on a kagome bilayer

We research the spin dynamics of a classical Heisenberg antiferromagnet with nearest-neighbor interactions on a quasi-two-dimensional kagome bilayer. This geometrically frustrated lattice consists of two kagome layers connected by a triangular-lattice linking layer. By combining Monte Carlo with precessional spin dynamics simulations, we compute the dynamical structure factor of the classical spin liquid in the kagome bilayer and investigate the thermal and dilution effects. While the low-frequency and long-wavelength dynamics of the cooperative paramagnetic phase is dominated by spin diffusion, weak magnon excitations persist at higher energies, giving rise to the half-moon pattern in the dynamical structure factor. In the presence of spin vacancies, the dynamical properties of the diluted system can be understood within the two-population picture. The spin diffusion of the “correlated” spin clusters is mainly driven by the zero-energy weather-vane modes, giving rise to an autocorrelation function that decays exponentially with time. On the other hand, the diffusive dynamics of the quasifree “orphan” spins leads to a distinctive longer-time power-law tail in the autocorrelation function. We discuss the implications of our work for the glassy behaviors observed in the archetypal frustrated magnet SrCr 9p Ga 12-9 pO 19 (SCGO).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Probing Ground-State Phase Transitions through Quench Dynamics

The study of quantum phase transitions requires the preparation of a many-body system near its ground state, a challenging task for many experimental systems. The measurement of quench dynamics, on the other hand, is now a routine practice in most cold atom platforms. Here we show that quintessential ingredients of quantum phase transitions can be probed directly with quench dynamics in integrable and nearly integrable systems. As a paradigmatic example, we study global quench dynamics in a transverse-field Ising model with either short-range or long-range interactions. When the model is integrable, we discover a new dynamical critical point with a nonanalytic signature in the short-range correlators. The location of the dynamical critical point matches that of the quantum critical point and can be identified using a finite-time scaling method. We extend this scaling picture to systems near integrability and demonstrate the continued existence of a dynamical critical point detectable at prethermal timescales. We quantify the difference in the locations of the dynamical and quantum critical points away from (but near) integrability. Thus, we demonstrate that this method can be used to approximately locate the quantum critical point near integrability. Furthermore, the scaling method is also relevant to experiments with finite time and system size, and our predictions are testable in near-term experiments with trapped ions and Rydberg atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Inverse transformation: unleashing spatially heterogeneous dynamics with an alternative approach to XPCS data analysis

X-ray photon correlation spectroscopy (XPCS), an extension of dynamic light scattering (DLS) in the X-ray regime, detects temporal intensity fluctuations of coherent speckles and provides scattering-vector-dependent sample dynamics at length scales smaller than DLS. The penetrating power of X-rays enables XPCS to probe the dynamics in a broad array of materials, including polymers, glasses and metal alloys, where attempts to describe the dynamics with a simple exponential fit usually fail. In these cases, the prevailing XPCS data analysis approach employs stretched or compressed exponential decay functions (Kohlrausch functions), which implicitly assume homogeneous dynamics. This paper proposes an alternative analysis scheme based upon inverse Laplace or Gaussian transformation for elucidating heterogeneous distributions of dynamic time scales in XPCS, an approach analogous to the CONTIN algorithm widely accepted in the analysis of DLS from polydisperse and multimodal systems. In conclusion, using XPCS data measured from colloidal gels, it is demonstrated that the inverse transform approach reveals hidden multimodal dynamics in materials, unleashing the full potential of XPCS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Dynamic Parameter Estimation with Physics-based Neural Ordinary Differential Equations

Accurate estimation of dynamic parameters of gen-erators is crucial to building a reliable model for dynamical studies and reliable operation of the power system. This paper develops a physics-based neural ordinary differential equations (ODE) approach to learn the parameters of generator dynamic model using phasor measurement units (PMU) data. We design a physics-based neural network to represent the swing equations of the power system dynamics. A loss function is defined as the difference between dynamic simulation results from the physics-based neural networks and pseudo PMU measurements. The parameters of generator dynamic model are iteratively updated using the neural ODEs and the adjoint method. By exploiting the mini-batch scheme in neural ODE training, the parameter estimation performance is significantly improved. Numerical study results on a 3-machine 9-bus system show that the proposed algorithm outperforms state-of-the-art baseline method in both computation time and dynamic parameter estimation accuracy.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2014

Dynamic-Feature Extraction, Attribution and Reconstruction (DEAR) Method for Power System Model Reduction

In interconnected power systems, dynamic model reduction can be applied on generators outside the area of interest to mitigate the computational cost with transient stability studies. This paper presents an approach of deriving the reduced dynamic model of the external area based on dynamic response measurements, which comprises of three steps, dynamic-feature extraction, attribution and reconstruction (DEAR). In the DEAR approach, a feature extraction technique, such as singular value decomposition (SVD), is applied to the measured generator dynamics after a disturbance. Characteristic generators are then identified in the feature attribution step for matching the extracted dynamic features with the highest similarity, forming a suboptimal ‘basis’ of system dynamics. In the reconstruction step, generator state variables such as rotor angles and voltage magnitudes are approximated with a linear combination of the characteristic generators, resulting in a quasi-nonlinear reduced model of the original external system. Network model is un-changed in the DEAR method. Tests on several IEEE standard systems show that the proposed method gets better reduction ratio and response errors than the traditional coherency aggregation methods.

Dynamic response, model reduction, orthogonal deco↗

DOE OSTI2015

Coarse-Grained Clustering Dynamics of Heterogeneously Coupled Neurons

The formation of oscillating phase clusters in a network of identical Hodgkin–Huxley neurons is studied, along with their dynamic behavior. The neurons are synaptically coupled in an all-to-all manner, yet the synaptic coupling characteristic time is heterogeneous across the connections. In a network of N neurons where this heterogeneity is characterized by a prescribed random variable, the oscillatory single-cluster state can transition—through (possibly perturbed) period-doubling and subsequent bifurcations—to a variety of multiple-cluster states. The clustering dynamic behavior is computationally studied both at the detailed and the coarse-grained levels, and a numerical approach that can enable studying the coarse-grained dynamics in a network of arbitrarily large size is suggested. Among a number of cluster states formed, double clusters, composed of nearly equal sub-network sizes are seen to be stable; interestingly, the heterogeneity parameter in each of the double-cluster components tends to be consistent with the random variable over the entire network: Given a double-cluster state, permuting the dynamical variables of the neurons can lead to a combinatorially large number of different, yet similar “fine” states that appear practically identical at the coarse-grained level. For weak heterogeneity we find that correlations rapidly develop, within each cluster, between the neuron’s “identity” (its own value of the heterogeneity parameter) and its dynamical state. For single- and double-cluster states we demonstrate an effective coarse-graining approach that uses the Polynomial Chaos expansion to succinctly describe the dynamics by these quickly established “identity-state” correlations. This coarse-graining approach is utilized, within the equation-free framework, to perform efficient computations of the neuron ensemble dynamics.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Evaluating crown scorch predictions from a computational fluid dynamics wildland fire simulator

Abstract Background Crown scorch—the heating of live leaves, needles, and buds in the vegetative canopy to lethal temperatures without widespread combustion—is one of the most common fire effects shaping post-fire canopies. Despite the ability of computational fluid dynamic models to finely resolve fire activity and buoyant plume dynamics including heterogenous 3D distributions of forest canopy heating, these models have had only limited use in simulating fire effects and have not been used to evaluate crown scorch. Here, we demonstrate a method of evaluating crown scorch using a computational fluid dynamics model, FIRETEC, and validate this approach by simulating the experiments that were used to develop Van Wagner’s 1973 crown scorch model. Results The average scorch height prediction from FIRETEC compares well with the empirical model derived by Van Wagner, which is the most widely used empirical model for crown scorch. We further find that the 3D buoyant plume dynamics from a steady and homogeneous idealized heat source on the ground results in a spatially heterogenous crown scorch pattern reflecting complex heating dynamics that are best represented by percent scorch rather than height of scorch. Conclusions The ability of the computational fluid dynamics model to capture variation in crown scorch due to 3D buoyant plume dynamics provides direct links between forest structure, fire behavior, and fire effects that can be used by forest managers and researchers to better understand how fires result in crown damage under various environmental and management scenarios.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Fluctuations in the Dynamics of Glasses

This research program addressed the question of the presence of dynamical heterogeneities – strong spatial fluctuations in the dynamics – in glass forming liquids and jamming systems near dynamical arrest, with particular emphasis on aging systems, i.e. systems that are in the (non-equilibrium) glass regime. The main goals proposed for this research were: to perform numerical simulations of atomistic structural glass models to characterize dynamical heterogeneities in out of equilibrium (aging) glassy systems; to test the hypothesis that a certain symmetry (“time reparametrization symmetry”) is present in microscopic models of glassy systems and that this symmetry can explain the main features of dynamical heterogeneities; and to test to what extent these phenomena are universal across different glassy/jamming systems. It was found that: most of the important features of dynamical heterogeneities in the aging regime could be described in terms of simple scaling behaviors; that some of the most basic theoretical models of glassy systems indeed have time reparametrization symmetry; that all tests performed in numerical simulation data were consistent with the predictions from time reparametrization symmetry; and that to a large degree, the main features of dynamical heterogeneities were universal across different glassy systems. Most of the findings that came out of this research have been reported in detail in eight papers in high quality journals, two unpublished but publicly accessible manuscripts, and 27 invited and contributed talks.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2015

Modeling dynamic stall on wind turbine blades under rotationally augmented flow fields

It is well known that airfoils under unsteady flow conditions with a periodically varying angle of attack exhibit aerodynamic characteristics different from those under steady flow conditions, a phenomenon commonly known as dynamic stall. It is also well known that the steady aerodynamic characteristics of airfoils in the inboard region of a rotating blade differ from those under steady two-dimensional (2D) flow conditions, a phenomenon commonly known as rotational augmentation. This paper presents an investigation of these two phenomena together in the inboard parts of wind turbine blades. This analysis is carried out using data from three sources: (1) the National Renewable Energy Laboratory’s Unsteady Aerodynamics Experiment Phase VI experimental data, including constant as well as continuously pitching blade conditions during axial operation, (2) data from unsteady Delayed Detached Eddy Simulations (DDES) carried out using the Technical University of Denmark’s in-house flow solver Ellipsys3D, and (3) data from a simplified model based on the blade element momentum method with a dynamic stall subroutine that uses rotationally augmented steady-state polars obtained from steady Phase VI experimental sequences, instead of the traditional 2D nonrotating data. The aim of this work is twofold. First, the blade loads estimated by the DDES simulations are compared to three select cases of the N sequence experimental data, which serves as a validation of the DDES method. Results show reasonable agreement between the two data in two out of three cases studied. Second, the dynamic time series of the lift and the moment polars obtained from the experiments are compared to those from the dynamic stall subroutine that uses the rotationally augmented steady polars. This allowed the differences between the stall phenomenon on the inboard parts of harmonically pitching blades on a rotating wind turbine and the classic dynamic stall representation in 2D flow to be investigated. Results from the dynamic stall subroutine indicated a good qualitative agreement between the model and the experimental data in many cases, which suggests that the current 2D dynamic stall model as used in BEM-based aeroelastic codes may provide a reasonably accurate representation of three-dimensional rotor aerodynamics when used in combination with a robust rotational augmentation model.

17 WIND ENERGY↗

NASA NTRS2013

Session 6: Dynamic Modeling and Systems Analysis

These presentations cover some of the ongoing work in dynamic modeling and dynamic systems analysis. The first presentation discusses dynamic systems analysis and how to integrate dynamic performance information into the systems analysis. The ability to evaluate the dynamic performance of an engine design may allow tradeoffs between the dynamic performance and operability of a design resulting in a more efficient engine design. The second presentation discusses the Toolbox for Modeling and Analysis of Thermodynamic Systems (T-MATS). T-MATS is a Simulation system with a library containing the basic building blocks that can be used to create dynamic Thermodynamic Systems. Some of the key features include Turbo machinery components, such as turbines, compressors, etc., and basic control system blocks. T-MAT is written in the Matlab-Simulink environment and is open source software. The third presentation focuses on getting additional performance from the engine by allowing the limit regulators only to be active when a limit is danger of being violated. Typical aircraft engine control architecture is based on MINMAX scheme, which is designed to keep engine operating within prescribed mechanical/operational safety limits. Using a conditionally active min-max limit regulator scheme, additional performance can be gained by disabling non-relevant limit regulators

control↗

NASA NTRS2015

Dynamic Characterization of Galfenol

A novel and precise characterization of the constitutive behavior of solid and laminated research-grade, polycrystalline Galfenol (Fe81:6Ga18:4) under under quasi-static (1 Hz) and dynamic (4 to 1000 Hz) stress loadings was recently conducted by the authors. This paper summarizes the characterization by focusing on the experimental design and the dynamic sensing response of the solid Galfenol specimen. Mechanical loads are applied using a high frequency load frame. The dynamic stress amplitude for minor and major loops is 2.88 and 31.4 MPa, respectively. Dynamic minor and major loops are measured for the bias condition resulting in maximum, quasi-static sensitivity. Three key sources of error in the dynamic measurements are accounted for: (1) electromagnetic noise in strain signals due to Galfenol's magnetic response, (2) error in load signals due to the inertial force of fixturing, and (3) time delays imposed by conditioning electronics. For dynamic characterization, strain error is kept below 1.2 % of full scale by wiring two collocated gauges in series (noise cancellation) and through lead wire weaving. Inertial force error is kept below 0.41 % by measuring the dynamic force in the specimen using a nearly collocated piezoelectric load washer. The phase response of all conditioning electronics is explicitly measured and corrected for. In general, as frequency increases, the sensing response becomes more linear due to an increase in eddy currents. The location of positive and negative saturation is the same at all frequencies. As frequency increases above about 100 Hz, the elbow in the strain versus stress response disappears as the active (soft) regime stiffens toward the passive (hard) regime.

Dynamic sensing↗

NASA NTRSDate not supplied

Supersonic Free-Flight Dynamics Testing in the Stratosphere

A new technique for obtaining stratospheric free-flight dynamics data for atmospheric entry capsules is described. The Stratospheric Projectile Experiment of Entry Dynamics (SPEED) represents a new approach to characterizing the free-flight dynamics of vehicles in the supersonic and transonic regime of flight. The SPEED test architecture leverages a stratospheric balloon and 3D-printed flight system in a novel two-stage configuration to deliver test articles to supersonic conditions in the atmosphere which achieve dynamic similitude with a full-scale vehicle. Designed to address the limitations of existing test facilities, SPEED captures the complete time evolution of a vehicle’s dynamic state while producing a statistically significant number of observed flight trajectories to address the stochastic nature of the wake-driven dynamic stability phenomenon. SPEED was developed over two years at NASA Ames Research Center. The inaugural flight of the SPEED test platform was conducted in the summer of 2024 where scaled capsules of the Mars Sample Return Earth Entry System and Dragonfly entry vehicle were tested. This paper describes the motivation for a new test architecture, operational constraints and achievable flight test envelopes, the design and development of the SPEED concept, and results from the demonstration flight. Techniques for estimating the dynamic aerodynamic characteristics and atmospheric conditions for each test article individually and as a collective are also discussed.

Entry Vehicle↗

DOE OSTI2020

Multi-Timescale Integrated Dynamic and Scheduling Model (MIDAS-Solar)

Solar photovoltaic (PV) installations have experienced unprecedented growth in the United States. PV will become not only an energy producer but also a necessary provider of ancillary services at multiple timescales. Conventional methods to simulate power system operations - such as long-term production simulation (which typically considers schedules from hours to minutes by using an optimization framework) and short-term transient studies (which simulate dynamics from seconds to sub-seconds using state variables and differential equations) are not sufficient for studying the multiple-timescale variation of solar generation and its impact on system reliability. Long-term system economics and short-term system dynamics are highly coupled, particularly when the penetration level of renewable generation is extremely high, because the uncertainty and variability of solar generation will impact both power systems steady-state and dynamic performance. This project helps meet and exceed the Solar Energy Technologies Office goal of systems integration by directly addressing this stability and reliability challenge for electric grid planning and operation. This will be accomplished by developing temporally comprehensive, closed-loop simulation models that seamlessly simulate power systems operations from economic scheduling (day-ahead to hours) to dynamic response analysis (seconds to sub-seconds). Both a multi-timescale grid model and an integrated PV model will be developed in this project to accurately study the impacts of PV variability and uncertainty on system reliability at multiple timescales. Using quasi-dynamic simulation methods and data-driven security assessment (DSA) criteria will allow the dynamic characteristics of PV to be fed forward into longer-timescale scheduling models for a complete understanding of the effect of short-term PV dynamics on bulk systems operations (e.g., reserve scheduling and deployment). Upon completion of the proposed model, this project will help operators accurately assess system reliability by deploying energy and reserve scheduling under critical contingency conditions and studying interactions among all types of essential reliability services provided by modern PV power plants.

14 SOLAR ENERGY↗