SEARCH · Engineering Papers

Results for “Parametric Capability”

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 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Code for the manuscript "Lagrangian Large Eddy Simulations vis Physics-informed Machine Learning"

We would like to disclose a code based on Python/Pytorch, in which we implement the data-driven physics-informed machine learning of Lagrangian large eddy simulation that is described in the manuscript "Lagrangian Large Eddy Simulations via Physics-Informed Machine Learning ". High-Reynolds number Homogeneous Isotropic Turbulence is fully described within the Navier-Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as Large Eddy Simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about the effects of sub-grid scales on the resolved scales. Here, we take an alternative approach and design novel LES heuristics stated in terms of Lagrangian particles moving with the flow. Our \emph{Lagrangian LES}, thus L-LES, is described by equations generalizing the weakly compressible Smoothed Particle Hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from Direct Numerical Simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The sub-grid scale contributions are modeled separately with physical constraints to account for the effects from un-resolved scales. We build the resulting model under the Differentiable Programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our Lagrangian LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

Tian, Yifeng↗

NASA NTRS2004

Spitzer Instrument Pointing Frame (IPF) Kalman Filter Algorithm

This paper discusses the Spitzer Instrument Pointing Frame (IPF) Kalman Filter algorithm. The IPF Kalman filter is a high-order square-root iterated linearized Kalman filter, which is parametrized for calibrating the Spitzer Space Telescope focal plane and aligning the science instrument arrays with respect to the telescope boresight. The most stringent calibration requirement specifies knowledge of certain instrument pointing frames to an accuracy of 0.1 arcseconds, per-axis, 1-sigma relative to the Telescope Pointing Frame. In order to achieve this level of accuracy, the filter carries 37 states to estimate desired parameters while also correcting for expected systematic errors due to: (1) optical distortions, (2) scanning mirror scale-factor and misalignment, (3) frame alignment variations due to thermomechanical distortion, and (4) gyro bias and bias-drift in all axes. The resulting estimated pointing frames and calibration parameters are essential for supporting on-board precision pointing capability, in addition to end-to-end 'pixels on the sky' ground pointing reconstruction efforts.

calibration↗

NASA NTRS2012

Parametric 3D Atmospheric Reconstruction in Highly Variable Terrain with Recycled Monte Carlo Paths and an Adapted Bayesian Inference Engine

We describe a method for accelerating a 3D Monte Carlo forward radiative transfer model to the point where it can be used in a new kind of Bayesian retrieval framework. The remote sensing challenge is to detect and quantify a chemical effluent of a known absorbing gas produced by an industrial facility in a deep valley. The available data is a single low resolution noisy image of the scene in the near IR at an absorbing wavelength for the gas of interest. The detected sunlight has been multiply reflected by the variable terrain and/or scattered by an aerosol that is assumed partially known and partially unknown. We thus introduce a new class of remote sensing algorithms best described as "multi-pixel" techniques that call necessarily for a 3D radaitive transfer model (but demonstrated here in 2D); they can be added to conventional ones that exploit typically multi- or hyper-spectral data, sometimes with multi-angle capability, with or without information about polarization. The novel Bayesian inference methodology uses adaptively, with efficiency in mind, the fact that a Monte Carlo forward model has a known and controllable uncertainty depending on the number of sun-to-detector paths used.

Monte Carlo↗

DOE OSTI2024

Recent Advances of PyROS: A Pyomo Solver for Nonconvex Two-Stage Robust Optimization in Process Systems Engineering

This poster highlights uncertainty and technical risk reduction capabilities in CCSI2, with a focus on robust optimization. It presents recent advances of the two-stage robust optimization (RO) solver PyROS and applications to advanced energy systems optimization. To demonstrate the computational performance and reliability of PyROS, a benchmarking study on a library of over 8,500 small-scale RO problems is presented. Further, PyROS is used to obtain robust system designs of a MEA-based CO2 absorber under uncertainty in the thermodynamic property models for a variety of CO2 capture rate threshold requirements. Overall, the results demonstrate that the PyROS solver, including recent extensions to multi-stage RO settings, provides a reliable avenue to optimize the design and operation of advanced energy systems subject to various sources of parametric uncertainty.

Sherman, Jason↗

DOE OSTI2024

Recent Advances of PyROS: A Pyomo Solver for Nonconvex Two-Stage Robust Optimization in Process Systems Engineering

Sherman, Jason↗

NASA NTRS1978

Investigation of a cryogenic thermal diode

This paper describes a series of parametric investigations to determine the effect of various fluid charges on the performance of a 0.635-cm-diameter spiral-artery, liquid-trap diode in both the forward and reverse modes. Specific parameters such as forward- and reverse-mode conductances, shutdown times and energies, and recovery to forward-mode operation, are evaluated for ethane as a working fluid in the temperature range 170 K to 220 K. Results indicate that the heat pipe will not reliably start up in the forward mode. However, startup can be initiated when preceded by a diode reversal. Also concluded are data which show the susceptibility of the diode to fluid charge and tilt. The optimum fluid charge was found to be 2.67 g, and transport capability at this charge was in excess of 1200 W-cm at 200 K. The diode in the reverse mode exhibited a rapid shutdown (within 9 min) with a shutdown energy of 1150 J (0.32 Wh).

Williams, R. J.↗

NASA NTRSDate not supplied

The Porous Microstructure Analysis (PuMA) software

The open-source Porous Microstructure Analysis (PuMA) software was implemented to offer an efficient framework for determining material characteristics from 3D microstructures. Its development was inspired by progress in X-ray microtomography, an imaging technology that captures the internal structure of materials in 3D, and even in a 4D temporal context. Over recent years, this method has transformed the domain of materials science due to its capability to non-destructively examine material microstructures while presenting digital data about their geometrical details. It has provided insights into materials relevant to several NASA missions, including heatshields, parachute fabrics, meteorites, and other advanced composites. PuMA, in its current version 3, delivers an array of features, spanning from basic geometric insights of a microstructure to intricate anisotropic thermo-elastic and chemical behavior. Specifically, the software evaluates morphological attributes (specific surface area, volume fractions, mean intercept lengths, orientation) and physical characteristics (conductivity, elasticity, permeability, and tortuosity). Additionally, it can model material degradation processes, such as oxidation and surface chemistry interactions. The software can generate synthetic microstructures, from straightforward geometrical designs to intricate woven and non-woven geometries. Coupling material generation and characterization enable parametric studies and sensitivity analysis to optimize the microstructural performance and inform design decisions and reliability assessment based on uncertainty quantification. A recent addition to PuMA includes the TomoSAM plugin, devised to incorporate the cutting-edge Segment Anything Model (SAM). SAM is a promptable deep learning model that can identify objects and create image masks in a zero-shot manner, based only on a few user clicks. The synergy between these tools aids in the segmentation of complex 3D datasets from tomography or other imaging techniques, which would otherwise require a laborious manual segmentation process.

Tomography↗

NASA NTRSDate not supplied

The Porous Microstructure Analysis (PuMA) software

The open-source Porous Microstructure Analysis (PuMA) software was created to offer an efficient framework for determining material properties from 3D microstructures. Its development was inspired by progress in X-ray microtomography, an imaging technology that captures the internal structure of materials in 3D, and even in a 4D temporal context. Over recent years, this method has transformed the domain of materials science due to its capability to non-destructively examine material microstructures while presenting digital data about their geometrical details. It has provided insights into materials relevant to several NASA missions, including heatshields, parachute fabrics, meteorites, and other advanced composites. PuMA, in its current version 3, delivers an array of features, spanning from basic geometric insights of a microstructure to intricate anisotropic thermo-elastic and chemical behavior. Specifically, the software evaluates morphological attributes (specific surface area, volume fractions, mean intercept lengths, orientation) and physical characteristics (conductivity, elasticity, permeability, and tortuosity). Additionally, it can model material degradation processes, such as oxidation and surface chemistry interactions. The software can generate synthetic microstructures, from straightforward geometrical designs to intricate woven and non-woven geometries. Coupling material generation and characterization enables parametric studies and sensitivity analysis to optimize the microstructural performance and inform design decisions and reliability assessment based on uncertainty quantification. A recent addition to PuMA includes the TomoSAM plugin, devised to incorporate the cutting-edge Segment Anything Model (SAM) into our image segmentation workflow. SAM is a promptable deep learning model that can identify objects and create image masks in a zero-shot manner, based only on a few user clicks. The synergy between these tools aids in the segmentation of complex 3D datasets from tomography and other imaging techniques, which would otherwise require a laborious manual segmentation process.

Tomography↗

NASA NTRS2020

Analysis of Open Hole Compression Specimens Using the CompDam Continuum Damage Mechanics Model

This report documents a validation study on the standard Open Hole Compression (OHC) laminate test specimen conducted as part of the NASA Advanced Composites Project (ACP). Tests were conducted on OHC specimens with hard, quasi-isotropic, and soft layups using digital image correlation and X-Ray computed tomography to capture the structural response and damage evolution. Progressive damage models were constructed for use with the CompDam continuum damage mechanics code following the best practices established during the ACP. Detailed interrogation of the analysis results and comparison with experimental measurements provide a basis for assessing the capability of the modeling approach for OHC. The structural response is found to be captured well, with strengths predicted within 3\% of the experimental values for hard and quasi-isotropic laminates. In the soft laminate, the model predicts failure to be more brittle than the nonlinear, ductile response that was measured. Damage states extracted from the models at the same load level as test measurements are overlaid to show directly the similarities and differences between test and analysis results. Studying the damage evolution predicted by the analysis reveals that the failure process is a competition between fiber damage and delamination/sub-laminate buckling, with fiber damage dominating the collapse in the hard laminate and sub-laminate buckling governing in the soft laminate. Finally, a series of parametric studies varying numerical solution parameters (mesh size, mass scaling) and physical properties (fiber direction compressive strength and toughness) reveal sensitivities and deficiencies of the model. To the authors' knowledge, this study is the first for OHC specimens to include detailed evaluation of damage mode interactions, direct overlay of predicted and measured damage states, and sensitivity of the predicted results to difficult-to-measure fiber direction material properties.

Continuum damage mechanics↗

DOE OSTI2021

Validation and Demonstration of Control System Functional Capabilities within the IES Plug-and-Play Simulation Environment

The concept of an integrated energy system (IES) is meant to combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework is used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES. The new RAVEN plugin for grid and capacity optimization (HERON) only addresses the limits that affect the production variables and the corresponding rates of variation (explicit constraints). However, other variables are subject to constraints, and the associated limits should be accounted for (implicit constraints). In particular, for the power dispatch problem, the optimization algorithm takes into account the limits on the electrical power output and the corresponding hourly power variations but does not consider other constraints on process variables whose response affects the service life of the IES. This report describes a scheme that allows accounting for implicit constraints without increasing the size of the optimization problem. To obtain a more accurate approximation of the nonlinear dynamic behavior, a parametric version of the dynamic mode decomposition with control (DMDc) algorithm was developed to derive the state-space representation matrices of the IES components at different scheduling parameter. Thanks to this approach, a more accurate approximation of the system response can be obtained, the limits imposed by thermal mechanical implicit constraints can be translated into power dispatch limits, and the feedbacks to HERON power dispatcher can be provided. To assess the developed methodology, a power dispatching test case composed of three power generating and storage units (Balance of Plant, Secondary Energy Source, Thermal Energy Storage) was developed. The power output of each one of the three units was optimized to meet the imposed time-dependent load demand trajectory and to maximize the IES profitability by meeting both the explicit and implicit constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1977

Parametric Performance of a Spiral-Artery, Liquid-Trap-Diode Heat Pipe

A series of parametric investigations are described which determine the effect of various fluid charges on the performance of a 0.635-cm-diam spiral-artery, liquid-trap diode in both the forward and reverse modes. Specific parameters such as forward- and reverse-mode conductances, shutdown times and energies, and recovery to forward-model operation, are evaluated for ethane as a working fluid in the temperature range 170 K to 220 K. Results indicate that the heat pipe will not reliably start up in the forward mode. However, startup can be initiated when preceded by a diode reversal. Also included are data which show the susceptibility of the diode to fluid charge and tilt. The optimum fluid charge was found to be 2.67 g and transport capability at this charge was in excess of 1200 W-cm at 200 K. The diode in the reverse mode exhibited a rapid shutdown (within 9 min) with a shutdown energy of 1150 J (0.32 Wh).

Williams, R. J.↗

NASA NTRS2014

Design Concepts for a Small Space-Based GEO Relay Satellite for Missions Between Low Earth and near Earth Orbits

The main purpose of the Small Space-Based Geosynchronous Earth orbiting (GEO) satellite is to provide a space link to the user mission spacecraft for relaying data through ground networks to user Mission Control Centers. The Small Space Based Satellite (SSBS) will provide services comparable to those of a NASA Tracking Data Relay Satellite (TDRS) for the same type of links. The SSBS services will keep the user burden the same or lower than for TDRS and will support the same or higher data rates than those currently supported by TDRS. At present, TDRSS provides links and coverage below GEO; however, SSBS links and coverage capability to above GEO missions are being considered for the future, especially for Human Space Flight Missions (HSF). There is also a rising need for the capability to support high data rate links (exceeding 1 Gbps) for imaging applications. The communication payload on the SSBS will provide S/Ka-band single access links to the mission and a Ku-band link to the ground, with an optical communication payload as an option. To design the communication payload, various link budgets were analyzed and many possible operational scenarios examined. To reduce user burden, using a larger-sized antenna than is currently in use by TDRS was considered. Because of the SSBS design size, it was found that a SpaceX Falcon 9 rocket could deliver three SSBSs to GEO. This will greatly reduce the launch costs per satellite. Using electric propulsion was also evaluated versus using chemical propulsion; the power system size and time to orbit for various power systems were also considered. This paper will describe how the SSBS will meet future service requirements, concept of operations, and the design to meet NASA users' needs for below and above GEO missions. These users' needs not only address the observational mission requirements but also possible HSF missions to the year 2030. We will provide the trade-off analysis of the communication payload design in terms of the number of links looking above and below GEO; the detailed design of a GEO SSBS spacecraft bus and its accommodation of the communication payload, and a summary of the trade study that resulted in the selection of the Falcon 9 launch vehicle to deploy the SSBS and its impact on cost reductions per satellite. ======================================================================== Several initiatives have taken place within NASA1 and international space agencies2 to create a human exploration strategy for expanding human presence into the solar system; these initiatives have been driven by multiple factors to benefit Earth. Of the many elements in the strategy one stands out: to send robotic and human missions to destinations beyond Low Earth Orbit (LEO), including cis-lunar space, Near-Earth Asteroids (NEAs), the Moon, and Mars and its moons.3, 4 The time frame for human exploration to various destinations, based on the public information available,1,4 is shown in Figure 1. Advance planning is needed to define how future space communications services will be provided in the new budget environment to meet future space communications needs. The spacecraft for these missions can be dispersed anywhere from below LEO to beyond GEO, and to various destinations within the solar system. NASA's Space Communications and Navigation (SCaN) program office provides communication and tracking services to space missions during launch, in-orbit testing, and operation phases. Currently, SCaN's space networking relay satellites mainly provide services to users below GEO, at Near Earth Orbit (NEO), below LEO, and in deep space. The potential exists for using a space-based relay satellite, located in the vicinity of various solar system destinations, to provide communication space links to missions both below and above its orbit. Such relays can meet the needs of human exploration missions for maximum connectivity to Earth locations and for reduced latency. In the past, several studies assessed the ability of satellite-based relays working above GEO in conjunction with Earth ground stations. Many of these focused on the trade between space relay and direct-to-Earth station links5,6,7. Several others focused on top-level architecture based on relays at various destinations8,9,10,11,12. Much has changed in terms of microwave and optical technology since the publication of the referenced papers; Ka-band communication systems are being deployed, optical communication is being demonstrated, and spacecraft buses are becoming increasingly more functional and operational. A design concept study was undertaken to access the potential for deploying a Small Space-Based Satellite (SSBS) relay capable of serving missions between LEO and NEO. The needs of future human exploration missions were analyzed, and a notional relay-based architecture concept was generated as shown in Fig. 1. Relay satellites in Earth through cis-Lunar orbits are normally located in stable orbits requiring low fuel consumption. Relay satellites for Mars orbit are normally selected based on the mission requirement and projected fuel consumption. Relay satellites have extreme commonalities of functions between them, differing only in the redundancy and frequencies used; therefore, the relay satellite in GEO was selected for further analysis since it will be the first step in achieving a relay-based architecture for human exploration missions (see Fig.Figure 2). The mission design methodology developed by the Collaborative Modeling for Parametric Assessment of Space Systems (COMPASS) team13 was used to produce the satellite relay design and to perform various design trades. At the start of the activity, the team was provided with the detailed concept of the notional architecture and the system and communication payload drivers.

Relay Satellites↗

DOE OSTI2022

Microstructure-Sensitive modeling of surface roughness and notch effects on extreme value fatigue response

Fatigue life knockdown due to surface roughness induces local stress concentrations. Physical fatigue experiments are crucial to qualify materials for use but are limited in their capability to examine the mechanisms underlying the fatigue crack formation driving force due to these stress concentrations. Microstructure-sensitive computational models can aid in this regard and are employed in this work to examine the effects of surface roughness in Al 7075-T6. Digital microstructure models are subjected to crystal plasticity finite element method (CPFEM) simulations to examine the combined effects of microstructure and grain scale asperities on the driving force for fatigue crack formation. Following elastic–plastic shakedown, mesoscale volume-averaged fatigue indicator parameters (FIPs) are computed within fatigue damage process zones of grains. An increase in the intensity of realistic surface roughness profiles corresponds accordingly to larger FIPs but it is difficult to precisely interpret the mechanisms and effects of this intensification. We thus parametrically investigate surface asperities that couple with microstructure and lack of constraint on slip at the free surface to describe microstructure-sensitive surface roughness knockdown effects on fatigue resistance. The depth of a single valley (i.e., notch) due to surface roughness is found to be more detrimental to fatigue resistance than the radius (i.e., notch width, not to be confused with notch acuity). Finally, we also find that different crystallographic textures, grain morphologies, and number of sampled grains have a limited effect on the depth of influence of these notches.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Dynamics of McMillan mappings I. McMillan multipoles

In this article, we consider two dynamical systems: the McMillan sextupole and octupole integrable mappings, originally proposed by Edwin McMillan. Both represent the simplest symmetric McMillan maps, characterized by a single intrinsic parameter. While these systems find numerous applications across various domains of mathematics and physics, some of their dynamical properties remain unexplored. We aim to bridge this gap by providing a comprehensive description of all stable trajectories, including the parametrization of invariant curves, Poincaré rotation numbers, and canonical action–angle variables. In the second part, we establish connections between these maps and general chaotic maps in standard form. Our investigation reveals that the McMillan sextupole and octupole serve as first-order approximations of the dynamics around the fixed point, akin to the linear map and quadratic invariant (known as the Courant–Snyder invariant in accelerator physics), which represents zeroth-order approximations (referred to as linearization). Furthermore, we propose a novel formalism for nonlinear Twiss parameters, which accounts for the dependence of rotation number on amplitude. This stands in contrast to conventional betatron phase advance used in accelerator physics, which remains independent of amplitude. Notably, in the context of accelerator physics, this new formalism demonstrates its capability in predicting dynamical aperture around low-order resonances for flat beams, a critical aspect in beam injection/extraction scenarios.

43 PARTICLE ACCELERATORS↗

NASA NTRS2021

Modeling of Guided Waves for Aerospace Applications

Advancements in computer hardware has led to new possibilities for rapid modeling and simulation capabilities across many scientific fields. Nondestructive evaluation (NDE) can benefit from increased use of simulation tools to guide optimization of inspection and health monitoring methods, enhance understanding of data, aid in development of defect characterization methods, and generate data sets for use with machine learning and model-assisted probability of detection. Recent work at NASA has entailed development and benchmarking of both custom simulation codes and commercial simulation tools for ultrasonic wave propagation. This paper describes recent work at NASA in modeling of guided waves in composites and other aerospace materials. Results and computational speeds for a composite benchmark case are reported for a custom finite difference Rotated Staggered Grid code and for the commercial finite element software package, Pogo. Recent progress in linking NDE models to parametric analysis tools is also discussed.

Nondestructive evaluation↗

DOE OSTI2025

Mineralized sclerites in the gorgonian coral Leptogorgia chilensis as a natural jamming system

The soft corals (Cnidaria, Octocorallia), a diverse group of colonial marine invertebrates, can reversibly tune their body stiffness in response to external stimuli. This capability is attributed to their dynamic skeletal systems, which consist of thousands of mineralized skeletal elements, called sclerites, embedded within a gel-like matrix that swells/deswells and unjams/jams the sclerites, thus modulating skeletal stiffness. While sclerite morphology is widely used for species identification, its role in the mechanical performance of a soft coral’s skeletal system is largely unknown. Here, we investigated structure-jamming relationships in sclerite-based skeletal architectures using the red gorgonian octocoral Leptogorgia chilensis as a model system. The sclerites of L. chilensis exhibit a shaft-like geometry with two axial branches and two sets of triradiate side branches, which are aligned with the crystallographic symmetry of the constituent magnesium-containing calcite. By combining multiscale three-dimensional (3D) structural characterization, parametric geometrical modeling, 3D printing, mechanical testing, and discrete element simulations, we demonstrate how sclerite geometry achieves a balanced jamming performance in terms of stiffness, weight, strength, and fracture resistance in comparison to alternative geometries parametrically modified from the native sclerites (e.g., changes in the length and number of side branches). Here, we also found that these performance metrics are achieved through the effective interlocking among side and axial branches, which is further enhanced by the fractal-like microscopic spikes on the branch tips. The findings in this natural jamming system offer insights for designing synthetic mechanotunable material architectures for a wide range of applications, from soft robotics to mechanical dampeners.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Deep nonparametric estimation of intrinsic data structures by chart autoencoders: Generalization error and robustness

Autoencoders have demonstrated remarkable success in learning low-dimensional latent features of high-dimensional data across various applications. Assuming that data are sampled near a low-dimensional manifold, we employ chart autoencoders, which encode data into low-dimensional latent features on a collection of charts, preserving the topology and geometry of the data manifold. Our paper establishes statistical guarantees on the generalization error of chart autoencoders, and we demonstrate their denoising capabilities by considering n noisy training samples, along with their noise-free counterparts, on a d-dimensional manifold. By training autoencoders, we show that chart autoencoders can effectively denoise the input data with normal noise. We prove that, under proper network architectures, chart autoencoders achieve a squared generalization error in the order of n–$\frac{2}{d+2}$log 4 n, which depends on the intrinsic dimension of the manifold and only weakly depends on the ambient dimension and noise level. We further extend our theory on data with noise containing both normal and tangential components, where chart autoencoders still exhibit a denoising effect for the normal component. As a special case, our theory also applies to classical autoencoders, as long as the data manifold has a global parametrization. Furthermore, our results provide a solid theoretical foundation for the effectiveness of autoencoders, which is further validated through several numerical experiments.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2001

Three-Dimensional Aerodynamic Instabilities In Multi-Stage Axial Compressors

This thesis presents the conceptualization and development of a computational model for describing three-dimensional non-linear disturbances associated with instability and inlet distortion in multistage compressors. Specifically, the model is aimed at simulating the non-linear aspects of short wavelength stall inception, part span stall cells, and compressor response to three-dimensional inlet distortions. The computed results demonstrated the first-of-a-kind capability for simulating short wavelength stall inception in multistage compressors. The adequacy of the model is demonstrated by its application to reproduce the following phenomena: (1) response of a compressor to a square-wave total pressure inlet distortion; (2) behavior of long wavelength small amplitude disturbances in compressors; (3) short wavelength stall inception in a multistage compressor and the occurrence of rotating stall inception on the negatively sloped portion of the compressor characteristic; (4) progressive stalling behavior in the first stage in a mismatched multistage compressor; (5) change of stall inception type (from modal to spike and vice versa) due to IGV stagger angle variation, and "unique rotor tip incidence" at these points where the compressor stalls through short wavelength disturbances. The model has been applied to determine the parametric dependence of instability inception behavior in terms of amplitude and spatial distribution of initial disturbance, and intra-blade-row gaps. It is found that reducing the inter-blade row gaps suppresses the growth of short wavelength disturbances. It is also concluded from these parametric investigations that each local component group (rotor and its two adjacent stators) has its own instability point (i.e. conditions at which disturbances are sustained) for short wavelength disturbances, with the instability point for the compressor set by the most unstable component group. For completeness, the methodology has been extended to describe finite amplitude disturbances in high-speed compressors. Results are presented for the response of a transonic compressor subjected to inlet distortions.

Tan, Choon S.↗