SEARCH · Search NASA
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.
This presentation discusses Spectral Proper Orthogonal Decomposition (SPOD) of Unsteady Pressure-Sensitive Paint (uPSP) measurements in the Space Launch System (SLS) Ascent Transient Aerodynamics Test (ATAT) and its application in the structural dynamics analysis.
This paper discusses Spectral Proper Orthogonal Decomposition (SPOD) of the Unsteady Pressure-Sensitive Paint (uPSP) measurements in recent NASA Ames wind tunnel test. The uPSP measurements were collected using Innovative Scientific Solutions, Inc. (ISSI) porous, fast-response pressure-sensitive paint, 40 ISSI four-inch air-cooled Light-Emitting Diodes, and 8 Phantom v2512 high-speed cameras at 10,000 frames per second in the uPSP Launch Vehicle Demonstration Test (LVDT) of the Space Launch System (SLS) vehicle in the 11-by 11-foot transonic test section of the Unitary Plan Wind Tunnel at NASA Ames Research Center in April 2024. SPOD is derived from a space-time proper orthogonal decomposition problem for statistically stationary flows. SPOD modes are determined in the frequency domain. Each SPOD mode oscillates at a single frequency. SPOD can be viewed as an extension of the Discrete Fourier Transform composition and the Dynamic Mode Decomposition. In this paper, the outputs of SPOD of the uPSP measurements in the uPSP LVDT are presented and the effectiveness of SPOD in the identification, diagnosis and analysis of the aerodynamic and aeroacoustic phenomena is demonstrated. The unsteady and dynamic property of the pressure field on the surface of the SLS Block 1B crew vehicle is presented with the visualization of the SPOD modes of the uPSP measurements in the tests of a Mach sweep run of the uPSP LVDT. The SPOD outputs were generated with the execution in parallel of a code in Python, with the library of Message Passing Interface for parallel processing, on the NASA Pleiades supercomputer. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. Funding was provided by the NASA Aerosciences Evaluation and Test Capabilities Portfolio Office.
This paper discusses Spectral Proper Orthogonal Decomposition (SPOD) of the Unsteady Pressure-Sensitive Paint (uPSP) measurements in recent NASA Ames wind tunnel test. The uPSP measurements were collected using Innovative Scientific Solutions, Inc. (ISSI) Porous, Fast-Response Pressure-Sensitive Paint (FP-XXX), 40 ISSI four-inch air-cooled Light-Emitting Diodes (LEDs), and 8 Phantom v2512 high-speed cameras at 10,000 frames per second in the Launch Vehicle Demo Test (LVDT) of the Space Launch System (SLS) Block 1B crew vehicle in the 11-by 11-foot transonic test section of the Unitary Plan Wind Tunnel (UPWT) at NASA Ames Research Center in April 2024. Fig. 1 shows the scale model of the SLS launch vehicle installed in the wind tunnel during the test.
A novel method for three-dimensional flowfield reconstruction is presented that utilizes spectral proper orthogonal decomposition (SPOD). The time-resolved PIV data are acquired in multiple two-dimensional spanwise planes, with adequate streamwise spacing to resolve the out-of-plane wavelength. Knowledge of the out-of-plane phase speed is also required, and for the present case this value was obtained using instantaneous hot-wire measurements. The method is validated by comparing phase-averaged results acquired with a forcing signal to the SPOD reconstructed results using the same data. The method is also applied to data with no forcing applied, and the results are shown to be similar to the forced case. This method is particularly useful for cases in which there are multiple frequencies of interest, or cases in which the phenomena occur over a broad range of frequencies. This method allows reconstruction of the frequencies of interest without the need for a forcing signal, so the flow can be studied in its natural (undisturbed) state.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Time-resolved, high-speed self-aligned focusing schlieren images were acquired in the NASA Glenn Research Center Transonic Turbine Blade Cascade facility to help understand the aerodynamic behavior of high-pressure, thick trailing edge turbine blades. The trailing edge thickness of 9% of axial chord tested represents simulated ceramic matrix composite fabrication constraints, which was verified previously to possess a high-loss flow regime at high inlet turbulence conditions over a narrow range of Reynolds numbers and at a fixed design exit Mach number of 0.74. Our high-speed images, which were acquired at 10 distinct Reynolds numbers, show a significant increase in energy from flow oscillations due to transonic vortex shedding at Reynolds numbers corresponding to the high loss conditions. For those conditions, strong acoustic waves turn into shock waves. Spectral Proper Orthogonal Decomposition of the high-speed images shows acoustic waves from trailing edge vortex shedding at all conditions, with increased spectral energy at the high-loss conditions and slightly increasing frequency (about 6%) as a function of Reynolds number. Analysis of potential feedback timing is performed using velocity fields from a previous LES simulation, considering different feedback mechanisms. Most noteworthy is the acoustic/shock-boundary layer interaction mechanism on the suction surface at the blade geometric throat, which likely plays an important role in realistic curved blade passages.
The second-mode instability on a 7◦ half-angle sharp cone at Mach 6 is analyzed using high-speed calibrated schlieren imagery at a frame rate near the expected fundamental frequency. Experiments were conducted in the NASA Langley 20-Inch Mach 6 facility at unit Reynolds number between 6.56×106 and 9.71×106 m−1. Time-resolved pixel intensity signals throughout the boundary layer are reconstructed using spatially available data in the schlieren images to recover an effective sampling rate of over 10 MHz; these are then converted to quantitative density gradients using a thin-lens-based calibration technique. A global analysis is performed on the schlieren data to investigate the nonlinear growth of the second-mode fundamental and harmonic content. Point-wise measures of the auto-bicoherence are used to identify specific triadic interactions and the locations of their highest levels of quadratic phase coupling. Significant resonance interactions between the second-mode fundamental and harmonic instabilities were found along with interactions between these and the mean flow. Bispectral mode decomposition is employed to educe the flow structures associated with these interactions. A similar analysis is performed for the power spectrum, with power spectral densities computed for each pixel’s time-series and spectral proper orthogonal decomposition employed to derive the modal structure and energy of the flow at specific frequencies. Comparisons between the bispectral quantities and second-mode power show that nonlinear interactions, particularly resonance interactions, are closely correlated with space-time modulation of disturbances during the nonlinear stage of transition.
Step excrescences are a type of surface imperfection encountered on swept wings of commercial aircraft, commonly because accessibility requirements prevent creating the wing’s surface from a single panel. If the step height is too large, super-critical, the flow can undergo an early transition to turbulence, which can be highly detrimental for the performance of the wing. In the case of a forward-facing step on a swept wing in a low-disturbance environment, stationary crossflow vortices can develop a significant amplitude upstream of the step and hence dominate the structure of the boundary-layer flow over the step. The main goal of the present investigation is to illuminate the path to transition supported by the fascinatingly complex flow field in the direct downstream vicinity of a step with a super critical height. The high-resolution, stereographic Particle Image Velocimetry (PIV) measurement dataset presently available for this flow field provides a complete description of the laminar flow for the execution of BiGlobal stability analysis in a plane parallel to the step. Although the notorious sensitivity of stability results to the description of the base flow demands a very careful uncertainty analysis of those results, it is argued that this very fact can be leveraged to produce new insight into the supported perturbation dynamics. In performing the analysis, several unsteady mode families are discovered that display the explosive perturbation expected for early transition to be induced. In considering domain widths equal to an integer-multiple of the incident crossflow-vortex wavelength and analyzing an extent of 5 crossflow-vortex wavelengths parallel to the step, it is found that the stability results converge while increasing the domain width. It is demonstrated, moreover, that the results for the wider domains can be approximated by appropriately averaging the results on neighboring single-crossflow-vortex-wavelength domains covering the same region. Besides being useful for computational purposes, this observed property suggests interpreting the instability mechanism as a distorted primary mechanism rather than a “proper” secondary mechanism. This follows in the context of the secondary in-stability analysis of three-dimensional boundary layers, because the secondary mechanism is usually characterized by being localized in a pocket of strong shear, while the distorted primary mechanism typically has an infinite support in the direction parallel to the step. Even though the growth rates are found to be sensitive to the interrogation-window size inherent to the PIV post-processing procedure, the spatial structure of the eigenfunctions is found to be relatively insensitive. Lastly, the spatial structure of the eigenfunctions corresponding to all velocity components are matched with the shape functions determined by computing the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content.
Understanding the flow of a crossflow-vortex-dominated boundary layer over a forward-facing step excrescence is necessary in order to mitigate an early transition scenario on the wings and tail of commercial aircraft. By performing BiGlobal stability analysis on the flowfield measured with high-resolution, stereographic Particle Image Velocimetry (PIV), previous work has shown that a family of unstable disturbances exists in the direct downstream vicinity of a supercritical, forward-facing step. The goal of the present paper is to improve the previously used stability approach on two fronts: 1) the stability problem is formulated in a non-orthogonal coordinate system and 2) the contribution of the out-of-plane base-flow derivatives is accounted for locally. The resulting eigen solutions feature significantly stronger growth rates and eigenfunctions that are localized above regions of reverse flow. The changes in the problem formulation furthermore establish a significant improvement in the comparison of the stability solutions with the results from the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content. A large destabilizing effect by the out-of-plane base-flow derivatives is determined to be justified, despite the assumption that these derivatives are small, because the effect can be reconstructed upon using an eigenvalue-correction formula that assumes the responsible terms are infinitesimally small. The unstable perturbation mechanisms are demonstrated to have a convective nature by assessing the relation between the frequency and the out-of-plane wavenumber, which indicates that their group speed does not approach zero. The last 2 facts, that 1) the perturbations are convective and 2) that the out-of-plane base-flow-derivative terms are small, remove any qualitative suspicions that the perturbation problem is not conducive to parabolization. This opens the path to analyzing this problem with a plane-marching approach.
Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.
Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.
The Unsteady Pressure-Sensitive Paint (uPSP) Development Team outof NASA Ames Research Center (ARC) has spent the past five yearsimproving the systems and processes to advance the uPSP technology for production-level wind tunnel testing. Already considered turnkey for small-scale and research applications, development in acquisition, calibration, data transfer, and data processing were needed to be useful to customers testing at NASA wind tunnels. This development focused at ARC at the Unitary Plan Wind Tunnel (UPWT) 11-by 11-ft Transonic Wind Tunnel due to the large optical access of the test section and the NASA Advanced Supercomputer(NAS), also located at ARC. A Launch Vehicle Demonstration Test (LVDT) at the UPWT represents a milestone of this initial phase of development where several new improvements were demonstrated in a production wind tunnel environment for the first time. LVDT was conducted in April 2024 and used a 4% forebody Space Launch System (SLS) Block 1B model as the test article. Both a crew and cargo configuration were tested, with varying Mach numbers, pressures, model positions, and camera magnifications. This paper summarizes the details of the test and is part of a collection with four additional papers that provide greater detail on: high-speed lifetime methodology, spectral proper orthogonal decomposition analysis, quality of high-resolution data compared to Corcos model, and data quality, calibration, and uncertainty.
Reliable noise prediction capabilities are essential to enable novel fuel efficient open rotor designs that can meet the community and cabin noise standards. Toward this end, immersed boundary methods have reached a level of maturity so that they are being frequently employed for specific real world applications within NASA. This paper demonstrates that our higher-order immersed boundary method provides the ability for aeroacoustic analysis of wake-dominated flow fields generated by highly complex geometries. This is the first of a kind aeroacoustic simulation of an open rotor propulsion system employing an immersed boundary method. In addition to discussing the peculiarities of applying the immersed boundary method to this moving boundary problem, we will provide a detailed aeroacoustic analysis of the noise generation mechanisms encountered in the open rotor flow. The simulation data is compared to available experimental data and other computational results employing more conventional CFD methods. The noise generation mechanisms are analyzed employing spectral analysis, proper orthogonal decomposition and the causality method.
An experimental study is conducted of the fluid-thermal-structural interaction of a clamped compliant panel exposed to a three dimensional shock-wave/boundary-layer interaction (SWBLI) induced by a Mach-6 compression ramp with a spanwise nonuniform incoming boundary layer. The nonuniform boundary layer was produced by placing trips on one side of the upstream flat plate, resulting in largely turbulent flow on the tripped side and transitional flow on the untripped side. Measurements of the flowfield confirmed that the tripped boundary layer contained elevated levels of unsteadiness, and the SWBLI was observed to vary from attached to fully separated as the ramp angle was increased from 10◦ to 38◦; the separation region on the tripped side of the panel was noticeably smaller, showing the elevated turbulence levels of the tripped-side flow to remain relatively localized rather than diffusing across the whole model. Full-field, time-resolved panel deformations were measured using high-speed photogrammetry and the vibrational response at each compression angle was characterized. Although the measured modes conformed largely to those from classical clamped-plate theory, some skewing of the mode shapes was observed. IR thermography highlighted regions of the compliant region where elevated temperatures were likely to promote thermal softening effects to the transient panel response. The quasi-static deformation and stress field was used to characterize the internal stress factor of each mode and showed a meaningful relationship between transient panel response and stress contained within each mode: modes with antinodes lying in high-stress areas of the plate tended to exhibit increases in vibrational frequency and decreases in vibrational power, whereas the opposite was true for modes with antinodes in low-stress areas.
A method that produces realistic simulations of atmospheric turbulence is developed and analyzed. The procedure makes use of a generalized spectral analysis, often called a proper orthogonal decomposition or the Karhunen-Loeve expansion. A set of criteria, emphasizing a realistic appearance, a correct spectral shape, and non-Gaussian statistics, is selected in order to evaluate the model turbulence. An actual turbulence record is analyzed in detail, providing both a standard for comparison and input statistics for the generalized spectral analysis. The simulation utilizes the eigenfunction-expansion procedure to produce preliminary time histories of the three velocity components simultaneously. It is found that important cross-statistical features are reasonably well-behaved.
A low-dimensional dynamic model for transitional buoyancy-driven flow in a differentially heated tall enclosure is presented. The full governing partial differential equations with the associated boundary conditions are solved by a spectral element method for a cavity of aspect ratio A=20. Proper orthogonal decomposition is applied to the oscillatory solution at Prandtl number Pr=P tau (omega) = 0.71 and Grashof number G tau (omega) = 3.2 x 10 (exp 4) to construct empirical eigenfunctions. Using the four most energetic empirical eigenfunctions for the velocity and temperature as basis functions and applying Galerkin's method, a reduced model consisting of eight nonlinear ordinary differential equations is obtained. Close to the 'design' conditions (P tau(omega) G tau(omega)), the low-order model (LOM) predictions are in excellent agreement with the predictions of the full model. In particular, the critical Grashof number at the onset of the first temporal flow instability (Hopf bifurcation) was well as the frequency and amplitude of oscillations at supercritical conditions are in excellent agreement with the predictions of the full model. Far from the 'design' conditions, the LOM predicts the existence of multiple stable steady solutions at large values of G tau, and a unique stable steady solution at small values of G tau, and exhibits hysteretic behavior that is qualitatively similar to that observed in direct numerical simulations based on the full model.
Galerkin reduced-order models (ROMs) often struggle to accurately capture multiscale fluid physics in challenging flow regimes such as flows experiencing compressibility effects. This in part stems from the global nature of both the basis construction problem and the spectral formulation itself. In this work, a multi-basis ROM is developed in an affine space based on proper orthogonal decomposition (POD) by projecting the full Navier-Stokes equations expressed in terms of the specific volume, velocity, and pressure primitive variables. The model is applied to high-fidelity numerical simulation datasets obtained for a canonical compressible flow configuration: the flow over a backward facing step at different subsonic Mach numbers. It is observed that application of an eigenvalue reassignment (ER) stabilization method is required to avoid early divergence of the ROM predictions for this configuration in the three Mach numbers tested. The sensitivity of the POD-ROM results to the choice of parameters in the stabilization algorithm is discussed.