Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Using a Canopy Model Framework to Improve Large-Eddy Simulations of the Neutral Atmospheric Boundary Layer in the Weather Research and Forecasting Model

A canopy model framework is implemented in the Weather Research and Forecasting Model to improve the accuracy of large-eddy simulations (LES) of the atmospheric boundary layer (ABL). The model includes two options that depend on the scale of surface roughness elements. A resolved canopy model, typically used to model flow through vegetation canopies, is employed when roughness elements are resolved by the vertical LES grid. In the case of unresolved roughness, a modified “pseudocanopy model” is developed to distribute drag over a shallow layer above the surface. Both canopy model options are validated against idealized test cases in neutral stability conditions and are shown to improve surface layer velocity profiles relative to simulations employing Monin–Obukhov similarity theory (MOST), which is commonly used as a surface boundary condition in ABL models. Use of the canopy model framework also leads to increased levels of resolved turbulence kinetic energy and turbulent stresses. Because LES of the ABL has a well-known difficulty recovering the expected logarithmic velocity profile (log law) in the surface layer, particular focus is placed on using the pseudocanopy model to alleviate this issue over a range of model configurations. Tests with varying surface roughness values, LES closures, and grid aspect ratios confirm that the pseudocanopy model generally improves log-law agreement relative to simulations that employ a standard MOST boundary condition. The canopy model framework thus represents a low-cost, easy-to-implement method for improving LES of the ABL.

17 WIND ENERGY↗
DOE OSTI2010

Improving Small Signal Stability through Operating Point Adjustment

ModeMeter techniques for real-time small signal stability monitoring continue to mature, and more and more phasor measurements are available in power systems. It has come to the stage to bring modal information into real-time power system operation. This paper proposes to establish a procedure for Modal Analysis for Grid Operations (MANGO). Complementary to PSS’s and other traditional modulation-based control, MANGO aims to provide suggestions such as increasing generation or decreasing load for operators to mitigate low-frequency oscillations. Different from modulation-based control, the MANGO procedure proactively maintains adequate damping for all time, instead of reacting to disturbances when they occur. Effect of operating points on small signal stability is presented in this paper. Implementation with existing operating procedures is discussed. Several approaches for modal sensitivity estimation are investigated to associate modal damping and operating parameters. The effectiveness of the MANGO procedure is confirmed through simulation studies of several test systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Nanocellulose‐Carboxymethylcellulose Electrolyte for Stable, High‐Rate Zinc‐Ion Batteries

Abstract Aqueous Zn ion batteries (ZIBs) are one of the most promising battery chemistries for grid‐scale renewable energy storage. However, their application is limited by issues such as Zn dendrite formation and undesirable side reactions that can occur in the presence of excess free water molecules and ions. In this study, a nanocellulose‐carboxymethylcellulose (CMC) hydrogel electrolyte is demonstrated that features stable cycling performance and high Zn 2+ conductivity (26 mS cm −1 ), which is attributed to the material's strong mechanical strength (≈70 MPa) and water‐bonding ability. With this electrolyte, the Zn‐metal anode shows exceptional cycling stability at an ultra‐high rate, with the ability to sustain a current density as high as 80 mA cm −2 for more than 3500 cycles and a cumulative capacity of 17.6 Ah cm −2 (40 mA cm −2 ). Additionally, side reactions, such as hydrogen evolution and surface passivation, are substantially reduced due to the strong water‐bonding capacity of the CMC. Full Zn||MnO 2 batteries fabricated with this electrolyte demonstrate excellent high‐rate performance and long‐term cycling stability (>500 cycles at 8C). These results suggest the cellulose‐CMC electrolyte as a promising low‐cost, easy‐to‐fabricate, and sustainable aqueous‐based electrolyte for ZIBs with excellent electrochemical performance that can help pave the way toward grid‐scale energy storage for renewable energy sources.

25 ENERGY STORAGE↗
DOE OSTI2022

Anomaly Detection and Mitigation in FACTS-based Wide-Area Voltage Control Systems using Machine Learning

With the increasing deployment of Flexible AC Transmission System (FACTS) devices in wide-area voltage control systems (WAVCS) for achieving improved voltage stability of bulk power systems, the possibility for cyber attacks on these systems is also increasing. Successful stealthy cyber attacks that are difficult to detect by traditional informational technology (IT)-based cybersecurity solutions or threshold-based bad data detectors can lead to a voltage collapse in power grid. This paper presents the testbed-based attacks implementation and real-time evaluation of machine learning (ML) algorithm for detecting and mitigating stealthy cyber attacks on FACTS-based WAVCS on a hardware-in-the-loop (HIL) testbed. Initially, we discuss the implementation of a fuzzy logic controller (FLC) that controls a Static VAR Compensator (SVC) device deployed in a two-area four-machine Kundur power system for improving transient voltage stability. Later, the ML-based Anomaly Detection and Mitigation (ADM) system is implemented on the cyber-physical HIL testbed to detect and mitigate various stealthy cyber attacks, which are injected in real-time over the wide-area network (WAN). The experimental results show accurate and effective performance of ADM system in detecting and mitigating anomalies while keeping the grid stable and within the system operating limits, as defined by the North America Electric Reliability Corporation (NERC).

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2010

Constraint-preserving boundary conditions in the 3+1 first-order approach

A set of energy-momentum constraint-preserving boundary conditions is proposed for the first-order Z4 case. The stability of a simple numerical implementation is tested in the linear regime (robust stability test), both with the standard corner and vertex treatment and with a modified finite-differences stencil for boundary points which avoids corners and vertices even in Cartesian-like grids. Moreover, the proposed boundary conditions are tested in a strong-field scenario, the Gowdy waves metric, showing the expected rate of convergence. The accumulated amount of energy-momentum constraint violations is similar or even smaller than the one generated by either periodic or reflection conditions, which are exact in the Gowdy waves case. As a side theoretical result, a new symmetrizer is explicitly given, which extends the parametric domain of symmetric hyperbolicity for the Z4 formalism. The application of these results to first-order Baumgarte-Shapiro-Shibata-Nakamura-like formalisms is also considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1997

Hypersonic Flow Computations on Unstructured Meshes

A method for computing inviscid hypersonic flow over complex configurations using unstructured meshes is presented. The unstructured grid solver uses an edge{based finite{volume formulation. Fluxes are computed using a flux vector splitting scheme that is capable of representing constant enthalpy solutions. Second{order accuracy in smooth flow regions is obtained by linearly reconstructing the solution, and stability near discontinuities is maintained by locally forcing the scheme to reduce to first-order accuracy. The implementation of the algorithm to parallel computers is described. Computations using the proposed method are presented for a sphere-cone configuration at Mach numbers of 5.25 and 10.6, and a complex hypersonic re-entry vehicle at Mach numbers of 4.5 and 9.8. Results are compared to experimental data and computations made with established structured grid methods. The use of the solver as a screening tool for rapid aerodynamic assessment of proposed vehicles is described.

Bibb, K. L.↗
DOE OSTI2014

Development of Diamond Membrane-Based X-Ray Optics: Fresnel Zone Plates, Beam Stops, and Grids

This project was proposed to develop the design and fabrication of zone plates for hard X-ray focusing up to 25 keV photon energies with 20 nm resolution (with zone doubling), based on diamond membranes. The benefits of using diamond are: improved mechanical robustness which among other advantages prevents collapse of high-aspect-ratio ring structures, a known high-aspect-ratio etch method, and excellent compatibility with other processes, excellent radiation hardness, extremely low X-ray absorption, and significantly improved thermal/dimensional stability as compared to alternative materials. The Phase I findings enable Phase II goals as well as significant progress in other diamond-membrane devices, such as diamond-supported X-ray optics (beam stops, apertures, phase shift disks for phase contrast X-ray microscopy, grids, resolution targets, coded apertures for coherence measurements, etc), X-ray windows and membrane–based detectors. Zone plates with outermost zones of 60 nm were written in 500 nm electron resist and transferred to a depth of up to 1.8 mm (exceeding the expectations) into a stack of engineered multilayer ultra-nano-crystalline diamond alternating non-conductive and conductive layers. The active material for focusing was added on top (Au and Ir). Several processes were investigated and validated to enable in Phase II the enhancement of the thickness of the structures, including polishing the over-plated zone plates and growing diamond on top of Ir-coated zone plates without the destruction of the nanostructures underneath. The zone plates were geometrically characterized.

07 ISOTOPE AND RADIATION SOURCES↗
NASA NTRS2016

Comparison of NACA 0012 Laminar Flow Solutions: Structured and Unstructured Grid Methods

In this paper we consider the solution of the compressible Navier-Stokes equations for a class of laminar airfoil flows. The principal objective of this paper is to demonstrate that members of this class of laminar flows have steady-state solutions. These laminar airfoil flow cases are often used to evaluate accuracy, stability and convergence of numerical solution algorithms for the Navier-Stokes equations. In recent years, such flows have also been used as test cases for high-order numerical schemes. While generally consistent steady-state solutions have been obtained for these flows using higher order schemes, a number of results have been published with various solutions, including unsteady ones. We demonstrate with two different numerical methods and a range of meshes with a maximum density that exceeds 8 × 106 grid points that steady-state solutions are obtained. Furthermore, numerical evidence is presented that even when solving the equations with an unsteady algorithm, one obtains steady-state solutions.

Swanson, R. C.↗
DOE OSTI2017

Importance of turbulence-chemistry interactions at low temperature engine conditions

The role of turbulence-chemistry interaction in autoignition and flame stabilization is investigated for spray flames at low temperature combustion (LTC) conditions by performing high-fidelity three-dimensional computational fluid dynamics (CFD) simulations. A recently developed Tabulated Flamelet Model (TFM) is coupled with a large eddy simulation (LES) framework and validated across a range of Engine Combustion Network (ECN) ambient temperature conditions for n-dodecane fuel. High resolution grids with 0.0625 mm minimum cell size and 25 million total cell count are implemented using adaptive mesh refinement over the spray and combustion regions. Simulations with these grids and multiple LES realizations, with a 103 species n-dodecane mechanism show good agreement with experimental data for all the ambient conditions investigated. This modeling approach with the computational cost advantage of tabulated chemistry is then extended towards understanding the auto-ignition and flame stabilization at an ambient temperature of 750 K. These low temperature conditions lead to substantially higher ignition delays and flame liftoff lengths, and significantly leaner combustion compared to conventional high temperature diesel combustion. These conditions also require the simulations to span significantly larger temporal and spatial dimensions thereby increasing the computational cost. The TFM approach is able to capture autoignition and flame liftoff length at the low temperature conditions. Significant differences with respect to mixing, species formation and flame stabilization are observed under low temperature compared to conventional diesel combustion. At higher ambient temperatures, formation of formaldehyde is observed in the rich region (phi > 1) followed by the formation of OH in the stoichiometric regions. Under low temperature conditions, formaldehyde is observed to form at leaner regions followed by the onset of OH formation in significantly lean regions of the flame. Qualitative differences between species formation and transient flame development for the high and low temperature conditions are presented. The two stage ignition process is further investigated by studying the species formation in mixture fraction space by solving 1D flamelet equations for different scalar dissipation rates and homogeneous reactor assumption. Results show that scalar dissipation causes these radicals to diffuse within the mixture fraction space. As a result, this significantly enhances ignition and plays a dominant role at such low temperature conditions which cannot be captured by the homogeneous reaction assumption based model.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1985

High-Latitude Filtering in a Global Grid-Point Model Using Model Normal Modes

The aim of high-latitude filtering in the vicinity of the poles is to avoid the excessively short time steps imposed on an explicit time-differencing scheme by linear stability due to fast moving inertia-gravity waves near the poles. The model normal mode expansion toward the problem of high-latitude filtering in a global shallow water model using the same philosophy as that used by Daley for the problem for large timesteps in P.E. models with explicit time integration schemes was applied.

Takacs, L. L.↗
NASA NTRS2020

Entropy Stable Nonconforming Discretizations with the Summation-By-Parts Property for Curvilinear Coordinates

The entropy conservative/stable algorithm of Friedrichet al.(2018) for hyperbolic conservation laws on nonconforming p-refined/coarsened Cartesian grids, is extended to curvilinear grids for the compressible Euler equations. The primary focus is on constructing appropriate coupling procedures across the curvilinear nonconforming interfaces. A simple and flexible approach is proposed that uses interpolation operators from one element to the other. On the element faces,the analytic metrics are used to construct coupling terms, while metric terms in the volume are approximated to satisfy a discretization of the geometric conservation laws. The resulting scheme is entropy conservative/stable, elementwise conservative, and freestream preserving. The accuracy and stability properties of the resulting numerical algorithm are shown to be comparable to those ofthe original conforming scheme (∼p+ 1 convergence) in the context of the isentropic Euler vortex and the inviscid Taylor–Green vortex problems on manufactured high order grids.

David Del Rey Fernandez↗
DOE OSTI2015

Simulation of the elementary evolution operator with the motional states of an ion in an anharmonic trap

Following a recent proposal of L. Wang and D. Babikov [J. Chem. Phys. 137, 064301 (2012)], we theoretically illustrate the possibility of using the motional states of a Cd{sup +} ion trapped in a slightly anharmonic potential to simulate the single-particle time-dependent Schrödinger equation. The simulated wave packet is discretized on a spatial grid and the grid points are mapped on the ion motional states which define the qubit network. The localization probability at each grid point is obtained from the population in the corresponding motional state. The quantum gate is the elementary evolution operator corresponding to the time-dependent Schrödinger equation of the simulated system. The corresponding matrix can be estimated by any numerical algorithm. The radio-frequency field which is able to drive this unitary transformation among the qubit states of the ion is obtained by multi-target optimal control theory. The ion is assumed to be cooled in the ground motional state, and the preliminary step consists in initializing the qubits with the amplitudes of the initial simulated wave packet. The time evolution of the localization probability at the grids points is then obtained by successive applications of the gate and reading out the motional state population. The gate field is always identical for a given simulated potential, only the field preparing the initial wave packet has to be optimized for different simulations. We check the stability of the simulation against decoherence due to fluctuating electric fields in the trap electrodes by applying dissipative Lindblad dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Geometric multigrid for an implicit-time immersed boundary method

The immersed boundary (IB) method is an approach to fluid-structure interaction that uses Lagrangian variables to describe the deformations and resulting forces of the structure and Eulerian variables to describe the motion and forces of the fluid. Explicit time stepping schemes for the IB method require solvers only for Eulerian equations, for which fast Cartesian grid solution methods are available. Such methods are relatively straightforward to develop and are widely used in practice but often require very small time steps to maintain stability. Implicit-time IB methods permit the stable use of large time steps, but efficient implementations of such methods require significantly more complex solvers that effectively treat both Lagrangian and Eulerian variables simultaneously. Moreover, several different approaches to solving the coupled Lagrangian-Eulerian equations have been proposed, but a complete understanding of this problem is still emerging. This paper presents a geometric multigrid method for an implicit-time discretization of the IB equations. This multigrid scheme uses a generalization of box relaxation that is shown to handle problems in which the physical stiffness of the structure is very large. Numerical examples are provided to illustrate the effectiveness and efficiency of the algorithms described herein. Finally, these tests show that using multigrid as a preconditioner for a Krylov method yields improvements in both robustness and efficiency as compared to using multigrid as a solver. They also demonstrate that with a time step 100–1000 times larger than that permitted by an explicit IB method, the multigrid-preconditioned implicit IB method is approximately 50–200 times more efficient than the explicit method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Underwater thermomagnetic generator for remote marine thermal energy harvesting and sensing

Thermomagnetic generators offer a promising approach for sustainable power generation in remote marine environments. Here, this study presents the design, thermal modeling, and experimental validation of a passively driven underwater thermomagnetic generator developed for powering ocean observation and monitoring sensors. The generator was evaluated under varying working fluids, thermal boundary conditions, and extended operation to assess real-world applicability. Two fluids, deionized water and silicone-based Thermal C5, were tested under simulated shallow- and deep-ocean conditions. Deionized water outperformed Thermal C5, especially in colder environments (~5°C), achieving a peak output of 2.7 mW due to larger temperature gradients and enhanced convective-evaporative heat transfer. Long-duration tests revealed a transient evaporation-condensation cycle that temporarily reduced rotor immersion and performance before stabilizing. The generator powered commercial marine sensors for over 6 h without external batteries, demonstrating the viability of compact, passively cooled thermomagnetic systems for autonomous, off-grid marine sensing.

13 HYDRO ENERGY↗
DOE OSTI2018

Aqueous Polysulfide–Air Battery with a Mediator-Ion Solid Electrolyte and a Copper Sulfide Catalyst for Polysulfide Redox

Aqueous sodium polysulfide–air batteries are a promising system for enabling low-cost, high energy density, environmentally benign batteries for grid-scale deployment. However, the sluggish redox kinetics of aqueous polysulfide requires the use of a catalyst. Cobalt sulfide is a traditional catalyst with good activity and stability, but it uses expensive and toxic cobalt. This study demonstrates for the first time the long-term cycling performance of a sodium polysulfide–air battery with a low-cost copper sulfide catalyst for polysulfide redox. Alkaline sodium polysulfide anolyte is paired with an acidic sodium phosphate buffered air catholyte by means of a Na+-ion solid-state electrolyte, which eliminates chemical crossover between the two electrodes while serving as a mediator-ion solid electrolyte. Long-term cycling is achieved by means of a decoupled charge and discharge air electrode. Furthermore, the effect of CuS redox activity in the polysulfide–air voltage window is tested by pre-discharging CuS in LiOH and in NaOH. It was found that the sodium polysulfide–air batteries with CuS catalyst exhibit good cycling performance, with the overpotential only showing a modest increase of 0.09 V after 100 cycles. Batteries assembled with CuS pre-discharged in LiOH exhibit a small boost in performance, while batteries assembled with CuS pre-discharged in NaOH show a small performance decline.

25 ENERGY STORAGE↗
NASA NTRS1984

Transonic analysis of canted winglets

A computational method developed to provide a transonic analysis for upper/lower surface wing-tip mounted winglets is described. Winglets with arbitrary planform, cant and toe angle, and airfoil section can be modeled. The embedded grid approach provides high flow field resolution and the required geometric flexibility. In particular, coupled Cartesian/cylindrical grid systems are used to model the complex geometry presented by canted upper/lower surface winglets. A new rotated difference scheme is introduced in order to maintain the stability of the small-disturbance formulation in the presence of large spanwise velocities. Wing and winglet viscous effects are modeled using a two-dimensional 'strip' boundary layer analysis. Correlations with wind tunnel and flight test data for three transport configurations are included.

Rosen, B. S.↗
NASA NTRS1994

Direct Replacement of Arbitrary Grid-Overlapping by Non-Structured Grid

A new approach that uses nonstructured mesh to replace the arbitrarily overlapped structured regions of embedded grids is presented. The present methodology uses the Chimera composite overlapping mesh system so that the physical domain of the flowfield is subdivided into regions which can accommodate easily-generated grid for complex configuration. In addition, a Delaunay triangulation technique generates nonstructured triangular mesh which wraps over the interconnecting region of embedded grids. It is designed that the present approach, termed DRAGON grid, has three important advantages: eliminating some difficulties of the Chimera scheme, such as the orphan points and/or bad quality of interpolation stencils; making grid communication in a fully conservative way; and implementation into three dimensions is straightforward. A computer code based on a time accurate, finite volume, high resolution scheme for solving the compressible Navier-Stokes equations has been further developed to include both the Chimera overset grid and the nonstructured mesh schemes. For steady state problems, the local time stepping accelerates convergence based on a Courant - Friedrichs - Leury (CFL) number near the local stability limit. Numerical tests on representative steady and unsteady supersonic inviscid flows with strong shock waves are demonstrated.

Kao, Kai-Hsiung↗