Search NASA⌕ Search

SEARCH · Search NASA

Results for “spectral proper orthogonal decomposition”

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

Spectral Proper Orthogonal Decomposition of uPSP Measurements in Recent NASA Ames Wind Tunnel Test

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.

Aeroacoustics↗

Spectral Proper Orthogonal Decomposition of uPSP Measurements in Recent NASA Ames Wind Tunnel Test

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.

Aeroacoustics↗

Three-Dimensional Instantaneous Flow-field Reconstruction Using Planar Spectral Proper Orthogonal Decomposition

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.

boundary-layer transition↗

Spectral proper orthogonal decomposition of active wake mixing dynamics in a stable atmospheric boundary layer

Recent advancements in the use of active wake mixing (AWM) to reduce wake effects on downstream turbines open new avenues for increasing power generation in wind farms. However, a better understanding of the fluid dynamics underlying AWM is still needed to make wake mixing a reliable strategy for wind farm flow control. In this work, a spectral proper orthogonal decomposition (SPOD) is used to analyze the dynamics of coherent flow structures that are induced in the wake through blade pitch actuation. The data are generated using the ExaWind software suite to perform large eddy simulations of an NREL 2.8 MW turbine operating in a stable atmospheric boundary layer. SPOD tracks the modal behavior of flow structures from their generation in the turbine induction field through their growth in the near-wake region and to their subsequent evolution and energy transfers in the far wake. SPOD is shown to be a useful tool in the context of AWM because it translates the wavenumber and frequency inputs to the turbine controller to structures in the wake. A decomposition of the radial shear stress flux in the wake is also developed using SPOD to measure the contribution of coherent flow structures to mean flow turbulent entrainment and wake recovery. The effectiveness of AWM is connected to its ability to excite inherent structures in the wake of the turbine that arise using baseline controls. The effects of AWM on blade loading are also analyzed by connecting the axial force along the blade to the SPOD analysis of the turbine induction field. Lastly, the performance of different AWM strategies is demonstrated in a two-turbine array.

17 WIND ENERGY↗

Surface Properties of Double-Fin Generated Shock-Wave/Boundary-Layer Interactions

Shock-Wave/Boundary-Layer Interaction (SBLI) is ubiquitous in high-speed air vehicles. SBLI surface flowfield properties of crossing shocks induced by symmetric Double-Fins (DFs) at fin angles of 8° and 10° at Mach 2 are experimentally investigated; a 10° Single-Fin (SF) is also explored as a comparison. Surface visualization is used to capture distinct flowfield structures such as separation and upstream influence. Steady and unsteady Pressure-Sensitive Paints (PSPs) are used to obtain global surface pressure fields. Characteristics of DF-SBLI surface topology are discussed by comparing the corresponding SF-SBLI. Here, the effect of fin angle on DF-SBLI is characterized, and unsteady dynamics of the pressure field are examined. Spectral Proper Orthogonal Decomposition (SPOD) is performed on unsteady PSP, and traveling surface pressure waves and their frequency-dependent behavior are discussed. Dispersion relation of surface pressure waves is revealed by performing two-dimensional space–time Fourier transform on unsteady PSP, and the results are compared with the phase velocity distributions obtained from SPOD. Accompanying group velocity behavior is examined for distinct SBLI regions wherein wave propagation is found dependent on its frequency. Steady surface pressures along the centerline are compared with previous experimental data. Fin-tip spacing is found to be an essential parameter for scaling of upstream influence length.

33 ADVANCED PROPULSION SYSTEMS↗

Dynamic Analysis of Reynolds Number Effects on Trailing Edge Transonic Vortex Shedding and Its Impact on Turbine Blade Aerodynamic Performance

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.

Trailing Edge↗

Global Analysis of Nonlinear Second-Mode Development in a Mach-6 Boundary Layer From High-Speed Schlieren Data

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.

boundary layer transition↗

Neural-network learning of SPOD latent dynamics

Here, we aim to reconstruct the latent space dynamics of high dimensional, quasi-stationary systems using model order reduction via the spectral proper orthogonal decomposition (SPOD). The proposed method is based on three fundamental steps: in the first, once that the mean flow field has been subtracted from the realizations (also referred to as snapshots), we compress the data from a high-dimensional representation to a lower dimensional one by constructing the SPOD latent space; in the second, we build the time-dependent coefficients by projecting the snapshots containing the fluctuations onto the SPOD basis and we learn their evolution in time with the aid of recurrent neural networks; in the third, we reconstruct the high-dimensional data from the learnt lower -dimensional representation. The proposed method is demonstrated on two different test cases, namely, a compressible jet flow, and a geophysical problem known as the Madden-Julian Oscillation. An extensive comparison between SPOD and the equivalent POD-based counterpart is provided and differences between the two approaches are highlighted. The numerical results suggest that the proposed model is able to provide low rank predictions of complex statistically stationary data and to provide insights into the evolution of phenomena characterized by specific range of frequencies. The comparison between POD and SPOD surrogate strategies highlights the need for further work on the characterization of the interplay of error between data reduction techniques and neural network forecasts.

97 MATHEMATICS AND COMPUTING↗

Aero-thermal numerical characterization of blunt fin-induced shock wave–boundary layer interaction and its control through leading-edge cooling injection

This study represents a novel evaluation of active flow control to alleviate the aerothermal penalties created by the blunt fin-induced shock wave–boundary layer interaction. The manuscript analyzes the effect of flow injection on a blunt fin-induced shock wave–boundary layer interaction via computational fluid dynamics simulations with various degrees of resolution. The impact on the mean flow topology and wall variables was investigated utilizing Reynolds-averaged Navier–Stokes simulations. Detached-eddy simulations revealed the low-frequency shock motion, shock wave–boundary layer, and horseshoe vortex interaction. The test article was exposed to two different incoming boundary layer thicknesses; the thicker boundary layer led to the appearance of larger turbulent scales. The Detached-eddy simulations revealed the time history of the shock wave–boundary layer interaction, focusing on the inception and development of the recirculated flow regions. Ultimately, spectral proper orthogonal decomposition was employed to identify the structures associated with the low-frequency shock motion caused by the shock wave–boundary layer interaction.

Lozano, Francisco (ORCID:0000000277871125)↗

Stability Analysis of the Flow Over a Swept Forward-Facing Step Using PIV Base Flows

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.

forward-facing step↗

Stability Analysis of the Flow over a Swept Forward-Facing Step using PIV Base Flowsin a Non-Orthogonal Coordinate System

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.

Boundary- Layer Transition↗

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

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.

ARMD↗

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

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.

ARMD↗

Invited: uPSP Launch Vehicle Demonstration Test at NASA Ames Research Center

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.

SLS↗