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 415 records · Page 23

Slow Wake Recovery and Low Turbulence Behind Wind Farms Parameterized in Mesoscale Simulations

Numerical weather prediction (NWP) and climate models equipped with wind-farm parameterizations (WFPs) can simulate cluster wake effects affecting downstream wind farms in both onshore and offshore environments. This study evaluates wake recovery behind a wind farm represented by the NWP-WFP approach in the Weather Research and Forecasting (WRF) model using either the Fitch et al. (2012) or Ma et al. (2022a, b) WFPs. Results are benchmarked against large-eddy simulations (LES) of an idealized offshore wind farm with aligned and staggered layouts under neutral atmospheric stability. Near-farm wake recovery is underestimated in NWP-WFP simulations due to its representation on a coarse mesoscale grid. This limitation leads to slow wake recovery through two interconnected mechanisms: (i) spatial gradients in the wind velocity field are weaker compared to LES and (ii) turbulence kinetic energy (TKE) remains low not because of excessive dissipation but due to insufficient shear production caused by these weakened gradients. For the scenario considered here, a wind-speed bias develops in the near-farm wake and persists into the far wake. Differences between the NWP-WFP simulations and LES emerge within a short distance downstream of the farm exit, where the mesoscale simulations recover too slowly. This reduced recovery contributes approximately 0.15-0.50 m s-1 to the near-farm wind-speed bias. The bias established in this region is not subsequently compensated for downstream but instead propagates into the far wake, where wind-speed differences of approximately 0.4-0.6 m s-1 remain up to 50 km downstream. Higher-resolution mesoscale simulations partially reduce this bias. Increasing turbine-added TKE or including subgrid wake effects provides additional improvement, but neither fully addresses the underlying cause. The slow wake recovery is not caused by limitations of the WFPs themselves, as it also occurs outside their region of influence, and adding subgrid wake effects does not significantly impact recovery. Rather, the slow wake recovery is a consequence of mesoscale flow representation. This behavior is not limited to regions downstream of the wind farm but is less visible within the farm, where wake recovery occurs simultaneously with turbine-induced momentum extraction. These results highlight the need for improved representations of wake recovery both within and downstream of wind farms. While enhanced subgrid modeling, shear-driven TKE production, and refined WFP formulations may improve intra-farm dynamics, accurately capturing near-farm wake recovery downstream remains challenging, as WFPs do not act in this region.

17 WIND ENERGY↗

Does the resistive tearing instability nonlinearly evolve to a fast reconnection mode

A fundamental problem in applying linear tearing instability theory to the rapid processes (particle acceleration, heating) in flares was the characteristically slow rate of reconnection. This problem can be at least partially overcome if the tearin mode nonlinearly evolves to a regime in which the reconnection rate is substantially enhanced, such as that for the Petschek configuration. This possibility was often suggested, and some numerical simulations appear to provide support for such a view. Numerical simulation are used to study the nonlinear evolution of the tearing stability and show that a fast Petschek-like regime may not be achieved. This conclusion follows when there are sufficient grid points within the diffusion region to completely resolve the nonlinear dynamic interactions in the diffusion layer. When the numerical resolution is not adequate, the solution does appear to approach a Petschek configuration. The resolved solution contains reverse flow vortices and current sheets, terminated with a current reversal, similar to those obtained by Syrovatsky (JEPT, 33, 933, 1971).

Steinolfson, R. S.↗

Interannual Variability of OLR as Observed by AIRS and CERES

The paper examines spatial anomaly time series of Outgoing Longwave Radiation (OLR) and Clear Sky OLR (OLR(sub CLR)) as determined using observations from CERES Terra and AIRS over the time period September 2002 through June 2011. We find excellent agreement of the two OLR data sets in almost every detail down to the x11deg spatial grid point level. The extremely close agreement of OLR anomaly time series derived from observations by two different instruments implies high stability of both sets of results. Anomalies of global mean, and especially tropical mean, OLR are shown to be strongly correlated with an El Nino index. These correlations explain that the recent global and tropical mean decreases in OLR over the time period studied are primarily the result of a transition from an El Nino condition at the beginning of the data record to La Nina conditions toward the end of the data period. We show that the close correlation of mean OLR anomalies with the El Nino Index can be well accounted for by temporal changes of OLR within two spatial regions, one to the east of, and one to the west of, the NOAA Nino-4 region. Anomalies of OLR in these two spatial regions are both strongly correlated with the El Nino Index as a result of the strong anti-correlation of anomalies of cloud cover and mid-tropospheric water vapor in these two regions with the El Nino Index.

Susskind, Joel↗

Overview of the NASA Glenn Flux Reconstruction Based High-Order Unstructured Grid Code

A computational fluid dynamics code based on the flux reconstruction (FR) method is currently being developed at NASA Glenn Research Center to ultimately provide a large- eddy simulation capability that is both accurate and efficient for complex aeropropulsion flows. The FR approach offers a simple and efficient method that is easy to implement and accurate to an arbitrary order on common grid cell geometries. The governing compressible Navier-Stokes equations are discretized in time using various explicit Runge-Kutta schemes, with the default being the 3-stage/3rd-order strong stability preserving scheme. The code is written in modern Fortran (i.e., Fortran 2008) and parallelization is attained through MPI for execution on distributed-memory high-performance computing systems. An h- refinement study of the isentropic Euler vortex problem is able to empirically demonstrate the capability of the FR method to achieve super-accuracy for inviscid flows. Additionally, the code is applied to the Taylor-Green vortex problem, performing numerous implicit large-eddy simulations across a range of grid resolutions and solution orders. The solution found by a pseudo-spectral code is commonly used as a reference solution to this problem, and the FR code is able to reproduce this solution using approximately the same grid resolution. Finally, an examination of the code's performance demonstrates good parallel scaling, as well as an implementation of the FR method with a computational cost/degree- of-freedom/time-step that is essentially independent of the solution order of accuracy for structured geometries.

High-Order Methods↗

Prediction of compliant wall drag reduction, part 1

Computer codes developed to test Bushnell's compliant wall drag reduction model are discussed. One code computes the evolution of mean velocity profiles during the period between bursts as forced by an imposed large-scale pressure pulse due to earlier bursts. Another code computes the local stability characteristics of these computed profiles. The programs use Chebyshev polynomials to resolve the normal boundary layer (y) direction and a staggered grid of mesh points to resolve the x direction. Typically, 257 grid points and 33 Chebyshev polynomials are used in the computations.

Orszag, S. A.↗

Calculations of rotational flows using stream function

The stream function equation is solved for steady two-dimensional (and axisymmetric) rotational flows. Both finite differences and finite volumes discretization techniques are studied, using generalized body fitted coordinates and unstructured staggered grids, respectively. For inviscid transonic flows, a new artificial viscosity scheme which does not produce any artificial vorticity is introduced, for the stability of the mixed flow calculations and for capturing shocks. The solution of Euler equations, in primitive variables, are also considered. The effects of the artificial viscosity and numerical boundary conditions on the total enthalpy and the vorticity distributions are demonstrated.

Hafez, M.↗

Development of an Aeroelastic Code Based on an Euler/Navier-Stokes Aerodynamic Solver

This paper describes the development of an aeroelastic code (TURBO-AE) based on an Euler/Navier-Stokes unsteady aerodynamic analysis. A brief review of the relevant research in the area of propulsion aeroelasticity is presented. The paper briefly describes the original Euler/Navier-Stokes code (TURBO) and then details the development of the aeroelastic extensions. The aeroelastic formulation is described. The modeling of the dynamics of the blade using a modal approach is detailed, along with the grid deformation approach used to model the elastic deformation of the blade. The work-per-cycle approach used to evaluate aeroelastic stability is described. Representative results used to verify the code are presented. The paper concludes with an evaluation of the development thus far, and some plans for further development and validation of the TURBO-AE code.

Bakhle, Milind A.↗

Recent Enhancements to USM3D Unstructured Flow Solver for Unsteady Flows

The NASA USM3D unstructured flow solver is undergoing extensions to address dynamic flow problems in support of NASA and NAVAIR efforts to study the applicability of Computational Fluid Dynamics tools for the prediction of aircraft stability and control characteristics. The initial extensions reported herein include two second-order time stepping schemes, Detached-Eddy Simulation, and grid motion. This paper reports the initial code verification and validation assessment of the dynamic flow capabilities of USM3D. The cases considered are the classic inviscid shock-tube problem, low Reynolds number wake shedding from a NACA 0012 airfoil, high Reynolds number DES-based wake shedding from a 4-to-1 length-to-diameter cylinder, and forced pitch oscillation of a NACA 0012 airfoil with inviscid and turbulent flow.

Pandya, Mohagna J.↗

Composite isogrid structures for parabolic surfaces

The invention relates to high stiffness parabolic structures utilizing integral reinforced grids. The parabolic structures implement the use of isogrid structures which incorporate unique and efficient orthotropic patterns for efficient stiffness and structural stability.

Silverman, Edward M.↗

A computational examination of directional stability for smooth and chined forebodies at high-alpha

Computational Fluid Dynamics (CFD) has been used to study aircraft forebody flowfields at low-speed, angle-of-attack conditions with sideslip. The purpose is to define forebody geometries which provide good directional stability characteristics under these conditions. The flows over the experimentally investigated F-5A forebody and chine type configuration, previously computed by the authors, were recomputed with better grid topology and resolution. The results were obtained using a modified version of CFL3D (developed at NASA Langley) to solve either the Euler equations or the Reynolds equations employing the Baldwin-Lomax turbulence model with the Degani-Schiff modification to account for massive crossflow separation. Based on the results, it is concluded that current CFD methods can be used to investigate the aerodynamic characteristics of forebodies to achieve desirable high angle-of-attack characteristics. An analytically defined generic forebody model is described, and a parametric study of various forebody shapes was then conducted to determine which shapes promote a positive contribution to directional stability at high angle-of-attack. An unconventional approach for presenting the results is used to illustrate how the positive contribution arises. Based on the results of this initial parametric study, some guidelines for aerodynamic design to promote positive directional stability are presented.

Ravi, Ramakrishnan↗

Assessing Climate Change Impacts on the Stability of Small Tidal Inlets: Part 2- Data Rich Environments

Climate change (CC) is likely to affect the thousands of bar-built or barrier estuaries (here referred to as Small tidal inlets - STIs) around the world. Any such CC impacts on the stability of STIs, which governs the dynamics of STIs as well as that of the inlet-adjacent coastline, can result in significant socio-economic consequences due to the heavy human utilisation of these systems and their surrounds. This article demonstrates the application of a process based snap-shot modelling approach, using the coastal morphodynamic model Delft3D, to 3 case study sites representing the 3 main STI types; Permanently open, locationally stable inlets (Type 1), Permanently open, alongshore migrating inlets (Type 2) and Seasonally/Intermittently open, locationally stable inlets (Type 3). The 3 case study sites (Negombo lagoon - Type 1, Kalutara lagoon - Type 2, and Maha Oya river - Type 3) are all located along the southwest coast of Sri Lanka. After successful hydrodynamic and morphodynamic model validation at the 3 case study sites, CC impact assessment are undertaken for a high end greenhouse gas emission scenario. Future CC modified wave and riverflow conditions are derived from a regional scale application of spectral wave models (WaveWatch III and SWAN) and catchment scale applications of a hydrologic model (CLSM) respectively, both of which are forced with IPCC Global Climate Model output dynamically downscaled to approximately 50 km resolution over the study area with the stretched grid Conformal Cubic Atmospheric Model CCAM. Results show that while all 3 case study STIs will experience significant CC driven variations in their level of stability, none of them will change Type by the year 2100. Specifically, the level of stability of the Type 1 inlet will decrease from 'Good' to 'Fair to poor' by 2100, while the level of (locational) stability of the Type 2 inlet will also decrease with a doubling of the annual migration distance. Conversely, the stability of the Type 3 inlet will increase, with the time till inlet closure increasing by approximately 75%. The main contributor to the overall CC effect on the stability of all 3 STIs is CC driven variations in wave conditions and resulting changes in longshore sediment transport, not Sea level rise as commonly believed.

IPC↗

Time-marching solution of incompressible Navier-Stokes equations for internal flow

Primitive variables with central differencing on a staggered grid are used in the present, factored ADI finite-difference scheme for artificial compressibility method solution of the incompressible Navier-Stokes equations, leading to a close coupling between velocity and pressure that both enhances stability and eliminates the need for artificial damping. Computational efficiency is enhanced through the use of a spatially variable, fixed Courant number-based time-step. The numerical results obtained for a driven cavity at Re of 10,000, with local cell Re as high as 100, exhibits no flow variable spatial oscillations on a 40 x 40 stretched grid solution.

Soh, W. Y.↗

Stability of interfaces with mesh refinement

A study is conducted of the stability of mesh refinement in space and time for several different interface equations and finite-difference approximations. First, a root condition which implies stability for the initial-boundary value problem for this type of interface is derived. From the root condition, the stability of several interface equations is proved, using the maximum principle. In some cases, the final verification steps can be done analytically; in other cases, a simple computer program has been written to check the condition for values of a parameter along the boundary of the unit circle. Using this method, stability for Lax-Wendroff with all the interface conditions considered, and for Leapfrog with interpolation interface conditions when the fine and coarse grids overlap is proved.

Berger, M. J.↗

NASA Tech Briefs, September 2005

Topivs include: Diamond-Coated Carbon Nanotubes for Efficient Field Emission; Improved Anode Coatings for Direct Methanol Fuel Cells; Advanced Ablative Insulators and Methods of Making Them; PETIs as High-Temperature Resin-Transfer-Molding Materials; Stable Polyimides for Terrestrial and Space Uses; Low-Density, Aerogel-Filled Thermal-Insulation Tiles; High-Performance Polymers Having Low Melt Viscosities; Nonflammable, Hydrophobic Aerogel Composites for Insulation; Front-Side Microstrip Line Feeding a Raised Antenna Patch; Medium-Frequency Pseudonoise Georadar; Facilitating Navigation Through Large Archives; Program for Weibull Analysis of Fatigue Data; Comprehensive Micromechanics-Analysis Code - Version 4.0; Component-Based Visualization System; Software for Engineering Simulations of a Spacecraft; LabVIEW Interface for PCI-SpaceWire Interface Card; Path Following with Slip Compensation for a Mars Rover; International Space Station Electric Power System Performance Code-SPACE; Software for Automation of Real-Time Agents, Version 2; Software for Optimizing Plans Involving Interdependent Goals; Computing Gravitational Fields of Finite-Sized Bodies; Custom Sky-Image Mosaics from NASA's Information Power Grid; ANTLR Tree Grammar Generator and Extensions; Generic Kalman Filter Software; Alignment Stage for a Cryogenic Dilatometer; Rugged Iris Mechanism; Treatments To Produce Stabilized Aluminum Mirrors for Cryogenic Uses; Making AlNx Tunnel Barriers Using a Low-Energy Nitrogen-Ion Beam; Making Wide-IF SIS Mixers with Suspended Metal-Beam Leads; Sol-Gel Glass Holographic Light-Shaping Diffusers; Automated Counting of Particles To Quantify Cleanliness; Phase Correction for GPS Antenna with Nonunique Phase Center; Compact Infrasonic Windscreen; Broadband External-Cavity Diode Laser; High-Efficiency Solar Cells Using Photonic-Bandgap Materials; Generating Solid Models from Topographical Data; Computationally Lightweight Air-Traffic-Control Simulation; Spool Valve for Switching Air Flows Between Two Beds; Partial Model of Insulator/ Insulator Contact Charging; Asymmetric Electrostatic Radiation Shielding for Spacecraft; and Reusable Hybrid Propellant Modules for Outer-Space Transport.

Source record↗

Distribution Grid Incentive Design with Unknown Agent Behavior

Motivation: During extreme events, traditional grid regulation methods (e.g., energy prices, net power injection limits) may be insufficient. While system operators typically lack control over end-user grid interactions, (e.g., energy demand), incentives can influence behavior - for example, a user that receives a grid-driven incentive may adjust their consumption or expose relevant control variables in response. Problem: Optimize for the best incentive subject to system stability constraints. However, user behavior is unknown to the SO - i.e., for a given incentive, the amount of curtailed load or control variables exposed is unknown.

feedback based control↗

Vorticity imbalance and stability in relation to convection

A complete synoptic-scale vorticity budget was related to convection storm development in the eastern two-thirds of the United States. The 3-h sounding interval permitted a study of time changes of the vorticity budget in areas of convective storms. Results of analyses revealed significant changes in values of terms in the vorticity equation at different stages of squall line development. Average budgets for all areas of convection indicate systematic imbalance in the terms in the vorticity equation. This imbalance resulted primarily from sub-grid scale processes. Potential instability in the lower troposphere was analyzed in relation to the development of convective activity. Instability was related to areas of convection; however, instability alone was inadequate for forecast purposes. Combinations of stability and terms in the vorticity equation in the form of indices succeeded in depicting areas of convection better than any one item separately.

Read, W. L.↗

Demonstration of Grid Services Using Mixed Grid-Forming and Grid-Following Technologies at the Wheatridge Renewable Energy Facility (Final Technical Report)

This Final Technical Report summarizes the analytical, modeling, and engineering work performed to evaluate grid-forming (GFM) inverter capabilities within a large hybrid renewable energy facility. The project focused on assessing the ability of advanced inverter-based resources to provide essential reliability services through coordinated operation of grid-forming (GFM) and grid-following (GFL) technologies. During Budget Period 1, the project developed a comprehensive framework of GFM performance metrics, including angle support, voltage regulation, frequency response, damping behavior, and current-limiting performance. Extensive electromagnetic transient (EMT) studies and hardware-in-the-loop (HIL) testing were conducted to validate GFM behavior under a range of operating conditions, including weak-grid scenarios, voltage disturbances, and multi-resource interactions. The project also established validated modeling approaches, plant-level control integration strategies, commissioning frameworks, and high-speed measurement infrastructure to support future field demonstration. Although the project concluded prior to field demonstration, the results provide a utility-scale foundation for evaluating, modeling, and deploying grid-forming technologies. The methodologies and tools developed contribute to industry understanding of inverter-based resource behavior and support future power system reliability under increasing renewable penetration.

14 SOLAR ENERGY↗

MOD silver metallization for photovoltaics

Photovoltaic cells require back side metallization and a collector grid system on the front surface. Both front and back surface metallizations should have good adhesion, low contact resistance, low sheet resistance, long term stability, and their deposition methods should not degrade the n-p junction. Advantages and disadvantages of different deposition methods are discussed.

Vest, G. M.↗