Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grid Stability”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Numerical Treatment of Degenerate Diffusion Equations via Feller's Boundary Classification, and Applications

A numerical method is devised to solve a class of linear boundary-value problems for one-dimensional parabolic equations degenerate at the boundaries. Feller theory, which classifies the nature of the boundary points, is used to decide whether boundary conditions are needed to ensure uniqueness, and, if so, which ones they are. The algorithm is based on a suitable preconditioned implicit finite-difference scheme, grid, and treatment of the boundary data. Second-order accuracy, unconditional stability, and unconditional convergence of solutions of the finite-difference scheme to a constant as the time-step index tends to infinity are further properties of the method. Several examples, pertaining to financial mathematics, physics, and genetics, are presented for the purpose of illustration.

Cacio, Emanuela↗
NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗
NASA NTRS2014

Comparitive Study of High-Order Positivity-Preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density ρ and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order Weighted Essentially Non-Oscillatory (WENO) schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative ρ may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case as well as for the case consisting of two species and one reaction.. This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Schemes↗
NASA NTRS2010

Analysis of the Effects of Streamwise Lift Distribution on Sonic Boom Signature

The streamwise lift distribution of a wing-canard-stabilator-body configuration was varied to study its effect on the near-field sonic boom signature. The investigation was carried out via solving the three-dimensional Euler equation with the OVERFLOW-2 flow solver. The computational meshes were created using the Chimera overset grid topology. The lift distribution was varied by first deflecting the canard then trimming the aircraft with the wing and the stabilator while maintaining constant lift coefficient of 0.05. A validation study using experimental results was also performed to determine required grid resolution and appropriate numerical scheme. A wide range of streamwise lift distribution was simulated. The result shows that the longitudinal wave propagation speed can be controlled through lift distribution thus controlling the shock coalescence.

Yoo, Seung Yeun (Paul)↗
DOE OSTI2025

Oscillation Risks of Grid-Following and Grid-Forming Inverter-Based Resources in Series-Compensated Networks

Here, this paper investigates the dynamic behavior of a grid-connected inverter-based resource (IBR) when connected radially to a series compensated line. Potential interactions between the series compensation and the IBR have been identified for both types: grid-following (GFL) or grid-forming (GFM). The study begins with electromagnetic transient (EMT) simulations to demonstrate stability issues. Subsequently, nonlinear analytical models are formulated in the dq frame, validated against the EMT simulation, and leveraged to assess eigenvalues and participation factors. Influencing factors of the dominant oscillation modes have been identified. The analysis results show that series compensation may make a mode associated with the synchronization unit unstable. Furthermore, customized feedback systems are built for the synchronizing loop. Series compensation can increase the sensitivity of the voltage phase angle towards the synchronizing angle, and introduce phase lag in the real power response towards the synchronizing angle. These factors may cause interactions with the phase-locked loop in GFL-IBR systems and with power-based synchronization in GFM-IBR systems, potentially leading to instability.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1984

A pseudospectral Legendre method for hyperbolic equations with an improved stability condition

A new pseudospectral method is introduced for solving hyperbolic partial differential equations. This method uses different grid points than previously used pseudospectral methods: in fact the grid are related to the zeroes of the Legendre polynomials. The main advantage of this method is that the allowable time step is proportional to the inverse of the number of grid points 1/N rather than to 1/n(2) (as in the case of other pseudospectral methods applied to mixed initial boundary value problems). A highly accurate time discretization suitable for these spectral methods is discussed.

Tal-Ezer, H.↗
NASA NTRS1986

A pseudospectral Legendre method for hyperbolic equations with an improved stability condition

A new pseudospectral method is introduced for solving hyperbolic partial differential equations. This method uses different grid points than previously used pseudospectral methods: in fact the grid points are related to the zeroes of the Legendre polynomials. The main advantage of this method is that the allowable time step is proportional to the inverse of the number of grid points 1/N rather than to 1/n(2) (as in the case of other pseudospectral methods applied to mixed initial boundary value problems). A highly accurate time discretization suitable for these spectral methods is discussed.

Tal-Ezer, Hillel↗
DOE OSTI2024

Comparative Study of Data-Driven Area Inertia Estimation Approaches on WECC Power Systems

With the increasing integration of inverter-based resources into the power grid, there has been a notable reduction in system inertia, potentially compromising frequency stability. To assess the suitability of existing area inertia estimation techniques for real-world power systems, this paper presents a rigorous comparative analysis of system identification, measurement reconstruction, and electromechanical oscillation-based area inertia estimation methodologies, specifically applied to the large-scale and multi-area WECC 240-bus power system. Comprehensive results show that the system identification-based approach exhibits superior robustness and accuracy relative to its counterparts.

area inertia estimation↗
NASA NTRS1978

On implicit finite-difference simulations of three-dimensional flow

An implicit finite-difference procedure for unsteady three-dimensional flow capable of handling arbitrary geometry through the use of general coordinate transformations is described. Viscous effects are optionally incorporated with a 'thin layer' approximation of the Navier-Stokes equations. An implicit approximate factorization technique is employed so that the small grid sizes required for spatial accuracy and viscous resolution do not impose stringent stability limitations. Results obtained from the program include transonic inviscid and laminar-turbulent solutions about simple body configurations. Comparisons with existing theories and experiments are made. Numerical accuracy and the effect of three-dimensional coordinate singularities are also discussed.

Pulliam, T. H.↗
NASA NTRS1980

Component-adaptive grid embedding

One of the major problems related to transonic flow prediction about realistic aircraft configuration is the generation of a suitable grid which encompasses such configurations. In general, each aircraft component (wing, fuselage, nacelle) requires a grid system that is usually incompatible with the grid systems of the other components; thus, the implementation of finite-difference methods for such geometrically-complex configurations is a difficult task. An approach was developed to treat such a problem. The basic idea is to generate different grid systems, each suited for a particular component. Thus, the flow field domain is divided into overlapping subdomains of different topology. These grid systems are then interfaced with each other in such a way that stability, convergence speed and accuracy are maintained.

Atta, E. H.↗
NASA NTRS1980

Implicit Finite-Difference Simulations of Three-Dimensional Compressible Flow

An implicit finite-difference procedure for unsteady three-dimensional flow capable of handling arbitrary geometry through the use of general coordinate transformations is described. Viscous effects are optionally incorporated with a "thin-layer" approximation of the Navier-Stokes equations. An implicit approximate factorization technique is employed so that the small grid sizes required for spatial accuracy and viscous resolution do not impose stringent stability limitations. Results obtained from the program include transonic inviscid or viscous solutions about simple body configurations. Comparisons with existing theories and experiments are made. Numerical accuracy and the effect of three-dimensional coordinate singularities are also discussed.

Pulliam, Thomas H.↗
NASA NTRS1999

Comparison of AGE and Spectral Methods for the Simulation of Far-Wakes

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This report demonstrates that a good compromise is possible with the Advected Grid Explicit (AGE) method. AGE has proven to be both efficient and accurate for simulating turbulent free-shear flows, including planar mixing layers and planar jets. Its efficiency results from its localized fully explicit finite difference formulation (Bisset 1998a,b) that is very straightforward to compute, outweighing the need for a fairly small timestep. Also, most of the successful simulations were slightly under-resolved, and therefore they were, in effect, large-eddy simulations (LES) without a sub-grid-scale (SGS) model, rather than direct numerical simulations (DNS). The principle is that the role of the smallest scales of turbulent motion (when the Reynolds number is not too low) is to dissipate turbulent energy, and therefore they do not have to be simulated when the numerical method is inherently dissipative at its resolution limits. Such simulations are termed 'auto-LES' (LES with automatic SGS modeling) in this report.

Bisset, D. K.↗
NASA NTRS2004

Astrophysics of Reference Frame Tie Objects

The Astrophysics of Reference Frame Tie Objects Key Science program will investigate the underlying physics of SIM grid objects. Extragalactic objects in the SIM grid will be used to tie the SIM reference frame to the quasi-inertial reference frame defined by extragalactic objects and to remove any residual frame rotation with respect to the extragalactic frame. The current realization of the extragalactic frame is the International Celestial Reference Frame (ICRF). The ICRF is defined by the radio positions of 212 extragalactic objects and is the IAU sanctioned fundamental astronomical reference frame. This key project will advance our knowledge of the physics of the objects which will make up the SIM grid, such as quasars and chromospherically active stars, and relates directly to the stability of the SIM reference frame. The following questions concerning the physics of reference frame tie objects will be investigated.

Johnston, Kenneth J.↗
NASA NTRS2017

A Comparative Study of Simulated and Measured Main Landing Gear Noise for Large Civil Transports

Computational results for the NASA 26%-scale model of a six-wheel main landing gear with and without a toboggan-shaped noise reduction fairing are presented. The model is a high-fidelity representation of a Boeing 777-200 aircraft main landing gear. A lattice Boltzmann method was used to simulate the unsteady flow around the model in isolation. The computations were conducted in free-air at a Mach number of 0.17, matching a recent acoustic test of the same gear model in the Virginia Tech Stability Wind Tunnel in its anechoic configuration. Results obtained on a set of grids with successively finer spatial resolution demonstrate the challenge in resolving/capturing the flow field for the smaller components of the gear and their associated interactions, and the resulting effects on the high-frequency segment of the farfield noise spectrum. Farfield noise spectra were computed based on an FWH integral approach, with simulated pressures on the model solid surfaces or flow-field data extracted on a set of permeable surfaces enclosing the model as input. Comparison of these spectra with microphone array measurements obtained in the tunnel indicated that, for the present complex gear model, the permeable surfaces provide a more accurate representation of farfield noise, suggesting that volumetric effects are not negligible. The present study also demonstrates that good agreement between simulated and measured farfield noise can be achieved if consistent post-processing is applied to both physical and synthetic pressure records at array microphone locations.

Konig, Benedikt↗
NASA NTRS2011

Unstructured CFD Aerodynamic Analysis of a Generic UCAV Configuration

Three independent studies from the United States (NASA), Sweden (FOI), and Australia (DSTO) are analyzed to assess the state of current unstructured-grid computational fluid dynamic tools and practices for predicting the complex static and dynamic aerodynamic and stability characteristics of a generic 53-degree swept, round-leading-edge uninhabited combat air vehicle configuration, called SACCON. NASA exercised the USM3D tetrahedral cell-centered flow solver, while FOI and DSTO applied the FOI/EDGE general-cell vertex-based solver. The authors primarily employ the Reynolds Averaged Navier-Stokes (RANS) assumption, with a limited assessment of the EDGE Detached Eddy Simulation (DES) extension, to explore sensitivities to grids and turbulence models. Correlations with experimental data are provided for force and moments, surface pressure, and off-body flow measurements. The vortical flow field over SACCON proved extremely difficult to model adequately. As a general rule, the prospect of obtaining reasonable correlations of SACCON pitching moment characteristics with the RANS formulation is not promising, even for static cases. Yet, dynamic pitch oscillation results seem to produce a promising characterization of shapes for the lift and pitching moment hysteresis curves. Future studies of this configuration should include more investigation with higher-fidelity turbulence models, such as DES.

Frink, Neal T.↗
DOE OSTI2024

Frequency Support From Electric Vehicles for Advancing Renewable Energy Integration

The integration of renewable energy resources (RERs) in the modern power grid is increasing rapidly because of aggressive decarbonization goals, lower costs, and increased government investment. However, higher penetrations of inverter-based generation can lead to frequency stability issues because of reduced system inertia. This paper develops a framework for quantifying the contribution of electric vehicles (EVs) toward providing frequency support to the grid and thus increasing the penetration limit of renewable energy resources (RERs). EVs are considered to provide both inertial response and primary frequency response support to the grid. A stochastic approach incorporating the uncertainties associated with the behavior of EVs is developed to derive the discharge limit of EV aggregators. A multi-machine system frequency response (MM-SFR) model is developed, which incorporates the dynamic virtual inertia and droop coefficients of EV aggregators derived from the EV control modules. Frequency security constraints are developed from this MM-SFR model, which, along with the converter voltage security and low voltage ride-through constraints, are integrated within a nonlinear optimization framework to determine the RER integration limit. Here, the efficacy of the proposed approach is validated using the RTS-GMLC test system.

25 ENERGY STORAGE↗
DOE OSTI2025

Probability Density Function Control of Frequency Fluctuations in Renewable-Rich Power Systems

The stochastic nature of renewable energy sources (RESs) necessitates treating power system frequency response as a random process with a nonstationary probability density function (PDF). Here, based upon the stochastic distribution control theory originated by the second author, this paper proposes a novel stochastic controller to improve the frequency PDF in power grids when integrating a large amount of RESs, thereby minimizing the effects of uncertainties and enhancing overall system stability. The key idea is to manipulate the controllable power generation resources so that the frequency PDF is make to follow a target PDF by using the stochastic distribution control theory originated by the second author. The proposed method can easily be plugged into existing automatic generation controls for multi-area transmission grids. The proposed method is validated via a modified Kundar's two area system and 240-bus Western Electricity Coordinating Council systems. The simulation results show that the proposed control shapes the frequency PDF narrower and sharper, leading to a notable improvement toward minimizing the effects of randomness and uncertainty during grid operation.

frequency↗
DOE OSTI2024

Voltage stability smart meter for analyzing voltage data and controlling an electrical power source and/or an electric appliance

Systems and methods for voltage stability monitoring and active/reactive power support are disclosed herein. In some embodiments, a smart electric meter of an end user in a grid power system can measure the voltage supplied to the end user via the grid power system, and can analyze the voltage data to detect critical voltage characteristics. The critical voltage characteristics may indicate that a voltage collapse event is likely. The smart electric meter can further estimate a voltage stability margin based on the voltage data. If necessary, the smart electric meter can control an electrical power source and/or an electric appliance positioned at or near the end user to increase the voltage stability margin.

Min, Liang↗