Search NASA⌕ Search

SEARCH · Search NASA

Results for “spectral stability analysis”

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 19 records

Comparison of methods for prediction of transition by stability analysis

The paper compares several methods of transition prediction of linear stability analysis. The spectral stability analysis code SALLY is used to analyze flows over laminar flow control wings. It is shown that transition prediction by the envelope method and a new modified wave packet method are comparable in reliability, but that the envelope method is more efficient computationally. This is based on the results which show that the wave packet method provides N factors which are at best as consistent as those of the envelope method; since the wave packet method is at least 3 times as expensive as the envelope method, the latter is recommended for engineering design calculations.

Malik, M. R.↗

Transition prediction using three dimensional stability analysis

Several methods of transition prediction by linear stability analysis are compared. The spectral stability analysis code SALLY is used to analyze flows over laminar flow control wings. It is shown that transition by the envelope method and a new modified wave packet method are comparable in reliability but that the envelope method is more efficient computationally.

Malik, M. R.↗

Stability analysis of spectral methods for hyperbolic initial-boundary value systems

A constant coefficient hyperbolic system in one space variable, with zero initial data is discussed. Dissipative boundary conditions are imposed at the two points x = + or - 1. This problem is discretized by a spectral approximation in space. Sufficient conditions under which the spectral numerical solution is stable are demonstrated - moreover, these conditions have to be checked only for scalar equations. The stability theorems take the form of explicit bounds for the norm of the solution in terms of the boundary data. The dependence of these bounds on N, the number of points in the domain (or equivalently the degree of the polynomials involved), is investigated for a class of standard spectral methods, including Chebyshev and Legendre collocations.

Gottlieb, D.↗

Discrete-Element and Material-Point Method (DEM and MPM) Based Solvers for Sustainable Technologies

We present the use of discrete element method (DEM) and material point method (MPM) in three relevant green technology applications that include biomass feedstock handling, lithium-ion battery manufacturing, and high-pressure reverse osmosis. Our open-source DEM and MPM solvers are developed using performance portable grid and particle management library, AMReX, thus enabling superior performance on NVIDIA and AMD GPUs with > 100 million particles. Our DEM solver resolves the motion of individual particles in a granular system and includes a bonded sphere method for modeling non-spherical particles along with Hertzian and liquid bridge-based contact models. We simulate highly variable biomass feedstock flows in large-scale hoppers for biofuel production and electrode calendering in battery manufacturing using DEM. Our simulations predict flow blockage in large scale biomass hoppers and electrode microstructure variations, thus providing valuable information for biofuel and battery manufacturers, respectively. The second half of the talk will be on MPM and its application towards pore resolved simulations of reverse osmosis membranes under compressive loads. We present a validation study of our MPM simulations with membrane microscopy imaging thus providing useful insights on membrane stability under high pressure conditions. We also present a spectral stability analysis of using linear hat, quadratic and cubic spline basis in MPM indicating regions of numerical stability.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗

Preliminary calibration results for the BATSE instrument on CGRO

Preliminary results pertaining to spectral reconstruction using Burst and Transient Source (BATSE) Large Area Detector measurements of solar flares are presented. The solar flare measurements are currently being used to fine tune the calibration of our data analysis software. The current status of the stability of spectral analysis, given the systematic errors present in burst location, are given. A brief description is given of enhancements to the input data for the atmospheric scattering algorithm that will be implemented in the data analysis software.

Pendleton, G. N.↗

A code for linear stability analysis

A new spectral code, Linear.x, has been written in FORTRAN 77 for the analysis of the linear stability of some basic state. While being generic and unrelated to any particular physical problem, the code provides for various common tasks, including global (eigenvalue spectra), local (single eigenvalues and eigenfunctions), table (one-dimensional and multidimensional tables of eigenvalues), curve (curves in parameter space), and others. A specific problem can be defined as a set of files some of which are included at compile time, while definitions, tasks, and parameters are read during run time. The code is currently used for the stability analysis of compressible flows.

Herbert, Thorwald↗

Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows

The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.

high pressure reverse osmosis↗

Discovering CO Adsorption and Desorption Pathways from Chemical Reaction Neural Network Modeling of Transient Kinetics Spectroscopy

Here, we demonstrate a data-driven approach to interpreting surface reactions by combining time-resolved gas pulsing infrared spectroscopy with chemical reaction neural networks (CRNNs). Using CO adsorption and desorption on Pd(111) at 460–490 K as a model system, we show how transient kinetic data can reveal detailed reaction mechanisms. Starting with a simple one-species model, we systematically evaluate increasingly complex mechanisms involving hollow and bridge site adsorption. Despite the similar goodness of fit to the same experimental absorbance data, our models predict distinct coverage dynamics for different adsorption sites. Through analysis of spectral peak stability and predicted dynamics, we identify a mechanism in which CO primarily adsorbs on bridge sites followed by rapid conversion to hollow sites as being the most physically consistent with experimental observations. This work provides a framework for extracting mechanistic insights from limited experimental data, demonstrating how machine learning can bridge the gap between transient kinetic measurements and a molecular-level understanding of surface reactions.

36 MATERIALS SCIENCE↗

Research in computational fluid dynamics and analysis of algorithms

Recently, higher-order compact schemes have seen increasing use in the DNS (Direct Numerical Simulations) of the Navier-Stokes equations. Although they do not have the spatial resolution of spectral methods, they offer significant increases in accuracy over conventional second order methods. They can be used on any smooth grid, and do not have an overly restrictive CFL dependence as compared with the O(N(exp -2)) CFL dependence observed in Chebyshev spectral methods on finite domains. In addition, they are generally more robust and less costly than spectral methods. The issue of the relative cost of higher-order schemes (accuracy weighted against physical and numerical cost) is a far more complex issue, depending ultimately on what features of the solution are sought and how accurately they must be resolved. In any event, the further development of the underlying stability theory of these schemes is important. The approach of devising suitable boundary clusters and then testing them with various stability techniques (such as finding the norm) is entirely the wrong approach when dealing with high-order methods. Very seldom are high-order boundary closures stable, making them difficult to isolate. An alternative approach is to begin with a norm which satisfies all the stability criteria for the hyperbolic system, and look for the boundary closure forms which will match the norm exactly. This method was used recently by Strand to isolate stable boundary closure schemes for the explicit central fourth- and sixth-order schemes. The norm used was an energy norm mimicking the norm for the differential equations. Further research should be devoted to BC for high order schemes in order to make sure that the results obtained are reliable. The compact fourth order and sixth order finite difference scheme had been incorporated into a code to simulate flow past circular cylinders. This code will serve as a verification of the full spectral codes. A detailed stability analysis by Carpenter (from the fluid Mechanics Division) and Gottlieb gave analytic conditions for stability as well as asymptotic stability. This had been incorporated in the code in form of stable boundary conditions. Effects of the cylinder rotations had been studied. The results differ from the known theoretical results. We are in the middle of analyzing the results. A detailed analysis of the effects of the heating of the cylinder on the shedding frequency had been studied using the above schemes. It has been found that the shedding frequency decreases when the wire was heated. Experimental work is being carried out to affirm this result.

Gottlieb, David↗

A spectral collocation solution to the compressible stability eigenvalue problem

A newly developed spectral compressible linear stability code (SPECLS) (staggered pressure mesh) is presented for analysis of shear flow stability, and applied to high speed boundary layers and free shear flows. The formulation utilizes the first application of a staggered mesh for a compressible flow analysis by a spectral technique. An order of magnitude less number of points is needed for equivalent accuracy of growth rates compared to those calculated by a finite difference formulation. Supersonic disturbances which are found to have oscillatory structures were resolved by a spectral multi-domain discretization, which requires a factor of three fewer points than the single domain spectral stability code. It is indicated, as expected, that stability of mixing layers is enhanced by viscosity and increasing Mach number. The mean flow involves a jet being injected into a quiescent gas. Higher temperatures of the injected gas is also found to enhance stability characteristics of the free shear layer.

Macaraeg, Michele G.↗

Suppression of pressure oscillations in an open cavity by passive pneumatic control

A computational investigation has been conducted to determine the effectiveness of a passive control technique in suppressing the oscillations in an open cavity exposed to a supersonic flow. Time-accurate solutions of the unsteady, Reynolds-averaged, Navier-Stokes equations were obtained with an explicit predictor-corrector algorithm; the passive control was implemented through the use of a simple linear pressure-velocity law along the porous cavity floor. The computational code was validated by comparisons with experimental data for the cavity flow without control. The computational results with control demonstrate that the oscillations are suppressed; the fluid dynamic mechanism of the control is seen to be a stabilization of the free-shear layer. Spectral analysis of the unsteady data shows that the resonant frequencies are essentially unchanged by the control. Beneficial reductions in the cavity resonant drag are also observed.

Chokani, N.↗

Spectral Clustering-Based Partitioning of Large-Scale Power Electronics-Based Power Systems for Small-Signal Stability Analysis

The nodal admittance matrix (NAM)-based approach is well-suited for small-signal stability analysis of large-scale power electronics-based power systems (PEPSs), as it preserves the system structure through its admittance matrix. Previous studies have explored partitioning such systems into subareas and interconnections to reduce computational burden; however, they lacked a formal algorithmic procedure for determining feasible partitions. While several grid partitioning methods, such as those based on graph theory or machine learning, exist in the literature, they cannot be directly applied to NAM-based analysis due to differing objectives and constraints. Here, this paper addresses this gap by presenting a systematic, step-by-step procedure for applying a spectral partitioning algorithm that yields a division of the system into subareas suitable for NAM-based analysis. The computational complexity of the proposed method is also derived to demonstrate its efficiency and justify the practicality of the resulting subarea decomposition. The performance of the partitioning method is evaluated by applying the spectral clustering-derived subareas and interconnections to the NAM-based partitioning approach on a 140-bus system. Computational times for the full-system and partitioned NAM analyses are compared using MATLAB. Additionally, PSCAD simulations of the complete system and partitioned subareas are carried out to verify the effectiveness of the proposed method.

Nupur [Univ. of Tennessee, Knoxville, TN (United S↗

GX 3+1: The Stability of Spectral Index as a Function of Mass Accretion Rate

We present an analysis of the spectral and timing properties observed in X-rays from neutron star (NS) binary GX 3+1 (4U 1744-26) during long-term transitions between the faint and bright phases superimposed on short-term transitions between lower banana (LB) and upper banana (UB) branches in terms of its color-color diagram, We analyze all observations of this source obtained with the Rossi X-ray Timing Explorer and BeppoSAX satellites, We find that the X-ray broadband energy spectra during these spectral transitions can be adequately reproduced by a composition of a low-temperature blackbody component, a Comptonized component (COMPTB), and Gaussian component We argue that the electron temperature kTe of the Compton cloud monotonically increases from 2.3 keY to 4.5 keY, when GX 3+1 makes a transition from UB to LB. We also detect an evolution of noise components (a very low frequency noise and a high-frequency noise) during these LB-UB transitions. Using a disk seed photon normalization of COMPTB, which is proportional to the mass accretion rate, we find that the photon power-law index Gamma is almost constant (Gamma = 2.00 +/- 0.02) when mass accretion rate changes by factor four. In addition, we find that the emergent spectrum is dominated by the strong Comptonized component We interpret this quasi-stability of the index Gamma and a particular form of the spectrum in the framework of a model in which the energy release in the transition layer located between the accretion disk and NS surface dominates that in the disk. Moreover, this index stability effect now established for GX 3+ I was previously found in the atoll source 4U 1728-34 and suggested for a number of other low-mass X-ray NS binaries. This intrinsic behavior of NSs, in particular for atoll sources, is fundamentally different from that seen in black hole binary sources where the index monotonically increases during spectral transition from the low state to the high state and then finally saturates at high values of mass accretion rate.

Seifana, Elena↗

Bi-stability and period-doubling cascade of frequency combs in exceptional-point lasers

Abstract Recent studies have demonstrated that a laser can self-generate frequency combs when tuned near an exceptional point (EP), where two cavity modes coalesce. These EP combs induce periodic modulation of the population inversion in the gain medium, and their repetition rate is independent of the laser cavity’s free spectral range. In this work, we perform a stability analysis that reveals two notable properties of EP combs, bi-stability and a period-doubling cascade. The period-doubling cascade enables halving of the repetition rate while maintaining the comb’s total bandwidth, presenting opportunities for the design of highly compact frequency comb generators.

exceptional points↗

Early Oscillation Detection for Hybrid DC/DC Converter Fault Diagnosis

This paper describes a novel fault detection technique for hybrid DC/DC converter oscillation diagnosis. The technique is based on principles of feedback control loop oscillation and RF signal modulations, and Is realized by using signal spectral analysis. Real-circuit simulation and analytical study reveal critical factors of the oscillation and indicate significant correlations between the spectral analysis method and the gain/phase margin method. A stability diagnosis index (SDI) is developed as a quantitative measure to accurately assign a degree of stability to the DC/DC converter. This technique Is capable of detecting oscillation at an early stage without interfering with DC/DC converter's normal operation and without limitations of probing to the converter.

Wang, Bright L.↗

Eigenvalue calculation procedure for an Euler/Navier-Stokes solver with application to flows over airfoils

A Lanczos procedure is presently applied to a Navier-Stokes (N-S) solver for eigenvalues and eigenvectors associated with the small-perturbation analysis of the N-S equations' finite-difference representation for airfoil flows; the matrix used is very large, sparse, real, and nonsymmetric. The Lanczos procedure is shown to furnish complete spectral information for the eigenvalues, as required for transient-stability analysis of N-S solvers.

Mahajan, Aparajit J.↗

High-resolution spectral analysis

Analyzer extends the range and resolution of a digital spectrum analyzer without placing stringent stability requirements on the sampling rate. It compares an unknown signal with a stable frequency standard.

Gillmore, W. F., Jr.↗