Results for “Coherent Mode Decomposition”

NASA NTRSDate not supplied

Coherent Mode Decomposition for Kolmogorov Optical Turbulence in a Finite Aperture

An analysis of the coherent mode decomposition of an optical field after propagation through atmospheric turbulence is presented. The coherent modes represent an ideal basis by which to decompose the field for design of mode-limited optical systems. Using the rotational symmetry of the Fredholm integral operator for Kolmogorov optical turbulence, it is shown that the coherent modes exhibit separable solutions classified by radial and azimuthal quantum numbers. The study of the coherent modes is then reduced that of the radial functions determined by the aperture-coherence diameter ratio and obscuration ratio. Analysis of the spectrum of eigenvalues yields sharp bounds on the efficiency of receivers using incoherent or coherent combining with mode-limited photonic devices in the presence of Kolmogorov turbulence. The effective number of modes needed to represent Kolmogorov optical turbulence is studied via the von Neumann entropy, purity, and largest eigenvalue, and the differences in the different definitions is discussed. The similarity of the coherent modes to linearly polarized (LP) fiber modes is quantified yielding a precise characterization of the maximum gain that can be achieved in mode-limited systems via mode shaping techniques. As a final application, a mode sorting technique is presented for optimally splitting power from atmospherically degraded light into a finite number of modes simultaneously maximizing total coupling efficiency and minimizing the average and instantaneous power ratio between channels.

free-space optical communication↗

NASA NTRSDate not supplied

Coherent Mode Decomposition for Kolmogorov Optical Turbulence in a Finite Aperture

An analysis of the coherent mode decomposition of an optical field after propagation through atmospheric turbulence is presented. The coherent modes represent an ideal basis by which to decompose the field for design of mode-limited optical systems. Due to rotational symmetry of the Fredholm integral operator for Kolmogorov optical turbulence, the coherent modes exhibit separable solutions classified by a radial and azimuthal mode index. The study of the coherent modes is then reduced to that of the radial functions determined by the coherence ratio and obscuration ratio. Analysis of the spectrum of eigenvalues yields sharp bounds on the efficiency of receivers using incoherent or coherent combining with mode-limited photonic devices in the presence of Kolmogorov turbulence. The effective number of modes needed to represent Kolmogorov optical turbulence is studied via the von Neumann entropy, purity, and largest eigenvalue, and the differences in the different definitions is discussed. The similarity of the coherent modes to linearly polarized (LP) fiber modes is quantified yielding a precise characterization of the maximum gain that can be achieved in mode-limited systems via mode shaping techniques. As a final application, a mode sorting technique is presented for optimally splitting power from atmospherically degraded light into a finite number of modes simultaneously maximizing total coupling efficiency and minimizing the difference in average power between channels.

free-space optical communication↗

NASA NTRSDate not supplied

Implementation of an Unsteady PSP System in the NASA Transonic Dynamics Tunnel

An unsteady pressure-sensitive paint (uPSP) system has been developed to provide time-resolved pressure measurements in the NASA Langley Transonic Dynamics Tunnel (TDT). Obtaining these measurements necessitated the development of environmental enclosures to protect the high-speed camera and ultraviolet lights required for uPSP from the harsh environment present during tunnel operation. Since the facility main drive was non-functioning during the testing window, performance of the uPSP system was demonstrated using an impinging jet with a passive oscillator attachment to provide unsteady flow with a known frequency independent of amplitude. Measurements were obtained for tunnel pressures ranging from 565 to 2116 psf, and model angles of attack between -4 degrees and 4 degrees. Results indicate that the system is capable of measuring surface pressure differentials on the order of 0.01 psi at full scale with a camera frame rate of at least 10 kHz. Spectral analysis shows that the fundamental frequency of the oscillating jet is captured by the uPSP system, as are the second and third harmonics. Dynamic mode decomposition highlights the dominant coherent spatial structures of the surface pressure, along with the associated frequency and growth rate of each mode, allowing for a de-noised reconstruction of the uPSP measurements. The experimental campaign outlined within this report also confirmed compatibility of the uPSP system with the TDT facility data acquisition system, and verified the successful integration with existing processing capabilities within the NASA advanced supercomputing environment.

unsteady pressure-sensitive paint↗

NASA NTRSDate not supplied

Implementation of an Unsteady PSP System in the NASA Transonic Dynamics Tunnel

An unsteady pressure-sensitive paint (uPSP) system has been developed to provide time-resolved pressure measurements in the NASA Langley Transonic Dynamics Tunnel (TDT). Obtaining these measurements necessitated the development of environmental enclosures to protect the high-speed camera and ultraviolet lights required for uPSP from the harsh environment present during tunnel operation. Since the facility main drive was non-functioning during the testing window, performance of the uPSP system was demonstrated using an impinging jet with a passive oscillator attachment to provide unsteady flow with a known frequency independent of amplitude. Measurements were obtained for tunnel pressures ranging from 565 to 2116 psf, and model angles of attack between -4 degrees and 4 degrees. Results indicate that the system is capable of measuring surface pressure differentials on the order of 0.01 psi at full scale with a camera frame rate of at least 10 kHz. Spectral analysis shows that the fundamental frequency of the oscillating jet is captured by the uPSP system, as are the second and third harmonics. Dynamic mode decomposition highlights the dominant coherent spatial structures of the surface pressure, along with the associated frequency and growth rate of each mode, allowing for a de-noised reconstruction of the uPSP measurements. The experimental campaign outlined within this report also confirmed compatibility of the uPSP system with the TDT facility data acquisition system, and verified the successful integration with existing processing capabilities within the NASA advanced supercomputing environment.

unsteady pressure-sensitive paint↗

NASA NTRSDate not supplied

Implementation of an Unsteady PSP System in the NASA TDT

An unsteady pressure-sensitive paint (uPSP) system has been developed to provide time-resolved pressure measurements in the NASA Langley Transonic Dynamics Tunnel (TDT). Obtaining these measurements necessitated the development of environmental enclosures to protect the high-speed camera and ultraviolet lights required for uPSP from the harsh environment present during tunnel operation. Since the facility main drive was non-functioning during the testing window, performance of the uPSP system was demonstrated using an impinging jet with a passive oscillator attachment to provide unsteady flow with a known frequency independent of amplitude. Measurements were obtained for tunnel pressures ranging from 565 to 2116 psf, and model angles of attack between -4 degrees and 4 degrees. Results indicate that the system is capable of measuring surface pressure differentials on the order of 0.01 psi at full scale with a camera frame rate of at least 10 kHz. Spectral analysis shows that the fundamental frequency of the oscillating jet is captured by the uPSP system, as are the second and third harmonics. Dynamic mode decomposition highlights the dominant coherent spatial structures of the surface pressure, along with the associated frequency and growth rate of each mode, allowing for a de-noised reconstruction of the uPSP measurements. The experimental campaign outlined within this report also confirmed compatibility of the uPSP system with the TDT facility data acquisition system, and verified the successful integration with existing processing capabilities within the NASA advanced supercomputing environment. Note: this presentation is an MP4 video with sound, color with a run time of 10 minutes 37 seconds.

unsteady pressure-sensitive paint↗

NASA NTRS1994

Global-scale modes of surface temperature variability on interannual to century timescales

Using 100 years of global temperature anomaly data, we have performed a singluar value decomposition of temperature variations in narrow frequency bands to isolate coherent spatio-temporal modes of global climate variability. Statistical significance is determined from confidence limits obtained by Monte Carlo simulations. Secular variance is dominated by a globally coherent trend; with nearly all grid points warming in phase at varying amplitude. A smaller, but significant, share of the secular variance corresponds to a pattern dominated by warming and subsequent cooling in the high latitude North Atlantic with a roughly centennial timescale. Spatial patterns associated with significant peaks in variance within a broad period range from 2.8 to 5.7 years exhibit characteristic El Nino-Southern Oscillation (ENSO) patterns. A recent transition to a regime of higher ENSO frequency is suggested by our analysis. An interdecadal mode in the 15-to-18 years period and a mode centered at 7-to-8 years period both exhibit predominantly a North Atlantic Oscillation (NAO) temperature pattern. A potentially significant decadal mode centered on 11-to-12 years period also exhibits an NAO temperature pattern and may be modulated by the century-scale North Atlantic variability.

Mann, Michael E.↗

NASA NTRSDate not supplied

Implicit Large-Eddy Simulation of Axisymmetric and Asymmetric Turbulent Shock Wave Boundary Layer Interaction at Mach 2.5

High-speed inlet flows feature turbulent shock wave boundary layer interactions at supersonic Mach numbers. A canonical Mach 2.5 shock wave boundary layer interaction experiment at the NASA Glenn Research Center was designed for axisymmetric and asymmetric interactions depending on the placement of the shock generator. Four different wall-resolved implicit large-eddy simulations of the experiment were carried out at one-tenth of the Reynolds number of the experiment: (1) Two axisymmetric interactions with shock generator radii of 9.2% and 14.7% of the test section diameter and (2) two asymmetric interactions with 14.7% shock generator radius and shock generator centerline offsets of one-third and one-fourth of the test section radius. For the two axisymmetric interactions, various sub-domains upstream of the interaction, within the interaction, and downstream of the interaction were analyzed with the proper orthogonal decomposition. For all three regions, the dominant modes are streamwise coherent structures with Strouhal numbers based on separation length of 0.15 and below. Fourier spectra of the wall-pressure coefficient for the interaction region have peaks at Strouhal numbers based on separation length of 0.15 and 0.04. For the two asymmetric interactions, the streamwise extent of the separation varies considerably in the azimuthal direction. For both shock generator offsets, the interaction is strongest for the azimuthal location that is closest to the shock generator. Compared to the axisymmetric interaction, the streamwise extent of the separated flow region is noticeably reduced and a substantial crossflow is observed. Fourier analyses of the unsteady flow indicate low-frequency content f or the upstream part of the interaction which is more two-dimensional. Away from the upstream region, with increasing sweep angle and crossflow, the low-frequency content is diminished.

Computational Fluid Dynamics↗

NASA NTRSDate not supplied

Self-Aligned Focusing Schlieren and OH Planar Laser-Induced Fluorescence Flow Visualization in a Dual-Mode Scramjet

Ahigh-speed self-aligned focusing schlieren (SAFS) system was used to visualize density gradients in and around the cavity flameholder of the combustor section of the University of Virginia Supersonic Combustion Facility (UVASCF). Images with this system were acquired at a framing rate of 110 kHz with no fuel injection, with fuel injection but no flame, and for fuel injection with combustion corresponding to a global equivalence ratio of 𝜙 = 0.18. Images with an air throttle in operation to modify the shock train location with fuel injection and with flame were also acquired. Simultaneous OH planar laser-induced fluorescence (PLIF) images were also acquired at a framing rate of 20 Hz. Results obtained with both visualization techniques are compared to one another to highlight how SAFS can complement more advanced flow visualization techniques and resolve dynamic behavior that may not otherwise be captured. Both proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) analysis techniques are applied to the SAFS image sequences to identify coherent periodic structures for the runs with fuel injection and combustion.

Brett F Bathel↗

NASA NTRS2023

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↗

NASA NTRS2023

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↗

NASA NTRS2000

CVS Decomposition of 3D Homogeneous Turbulence Using Orthogonal Wavelets

This paper compares the filtering used in Coherent Vortex Simulation (CVS) decomposition with an orthogonal wavelet basis, with the Proper Orthogonal Decomposition (POD) or Fourier filtering. Both methods are applied to a field of Direct Numerical Simulation (DNS) data of 3D forced homogeneous isotropic turbulence at microscale Reynolds number R(sub lambda) = 168. We show that, with only 3%N retained modes, CVS filtering separates the coherent vortex tubes from the incoherent background flow. The latter is structureless, has an equipartition energy spectrum, and has a Gaussian velocity probability distribution function (PDF) and an exponential vorticity PDF. On the other hand, the Fourier basis does not extract the coherent vortex tubes cleanly and leaves organized structures in the residual high wavenumber modes whose PDFs are stretched exponentials for both the velocity and the vorticity.

Farge, Marie↗

NASA NTRS2014

Application of Spectral Analysis Techniques in the Intercomparison of Aerosol Data. Part II: Using Maximum Covariance Analysis to Effectively Compare Spatiotemporal Variability of Satellite and AERONET Measured Aerosol Optical Depth

Moderate Resolution Imaging SpectroRadiometer (MODIS) and Multi-angle Imaging Spectroradiomater (MISR) provide regular aerosol observations with global coverage. It is essential to examine the coherency between space- and ground-measured aerosol parameters in representing aerosol spatial and temporal variability, especially in the climate forcing and model validation context. In this paper, we introduce Maximum Covariance Analysis (MCA), also known as Singular Value Decomposition analysis as an effective way to compare correlated aerosol spatial and temporal patterns between satellite measurements and AERONET data. This technique not only successfully extracts the variability of major aerosol regimes but also allows the simultaneous examination of the aerosol variability both spatially and temporally. More importantly, it well accommodates the sparsely distributed AERONET data, for which other spectral decomposition methods, such as Principal Component Analysis, do not yield satisfactory results. The comparison shows overall good agreement between MODIS/MISR and AERONET AOD variability. The correlations between the first three modes of MCA results for both MODIS/AERONET and MISR/ AERONET are above 0.8 for the full data set and above 0.75 for the AOD anomaly data. The correlations between MODIS and MISR modes are also quite high (greater than 0.9). We also examine the extent of spatial agreement between satellite and AERONET AOD data at the selected stations. Some sites with disagreements in the MCA results, such as Kanpur, also have low spatial coherency. This should be associated partly with high AOD spatial variability and partly with uncertainties in satellite retrievals due to the seasonally varying aerosol types and surface properties.

aerosols↗

NASA NTRS1987

An orthogonal decomposition of the axisymmetric jet mixing layer utilizing cross-wire velocity measurements

The structure of the mixing layer of the axisymmetric turbulent jet was investigated using cross-wire velocity measurements, which generated radial and streamwise cross-spectral data at eight radial positions and 25 azimuthal positions. The results show the existence of an azimuthal coherent ringlike structure dominated by axisymmetric mode near the potential core; however, from the center of the mixing layer out toward its low-speed side, the fourth, fifth, and sixth modes are dominating. A mechanism for the life cycle of coherent structures and for turbulence production in the axisymmetric jet mixing layer is proposed, which accounts for the observed phenomena and is amenable to dynamical analysis using the equations of motion.

Glauser, Mark N.↗

NASA NTRS1996

Structure Identification Within a Transitioning Swept-Wing Boundary Layer

Extensive measurements are made in a transitioning swept-wing boundary layer using hot-film, hot-wire and cross-wire anemometry. The crossflow-dominated flow contains stationary vortices that breakdown near mid-chord. The most amplified vortex wavelength is forced by the use of artificial roughness elements near the leading edge. Two-component velocity and spanwise surface shear-stress correlation measurements are made at two constant chord locations, before and after transition. Streamwise surface shear stresses are also measured through the entire transition region. Correlation techniques are used to identify stationary structures in the laminar regime and coherent structures in the turbulent regime. Basic techniques include observation of the spatial correlations and the spatially distributed auto-spectra. The primary and secondary instability mechanisms are identified in the spectra in all measured fields. The primary mechanism is seen to grow, cause transition and produce large-scale turbulence. The secondary mechanism grows through the entire transition region and produces the small-scale turbulence. Advanced techniques use Linear Stochastic Estimation (LSE) and Proper Orthogonal Decomposition (POD) to identify the spatio-temporal evolutions of structures in the boundary layer. LSE is used to estimate the instantaneous velocity fields using temporal data from just two spatial locations and the spatial correlations. Reference locations are selected using maximum RMS values to provide the best available estimates. POD is used to objectively determine modes characteristic of the measured flow based on energy. The stationary vortices are identified in the first laminar modes of each velocity component and shear component. Experimental evidence suggests that neighboring vortices interact and produce large coherent structures with spanwise periodicity at double the stationary vortex wavelength. An objective transition region detection method is developed using streamwise spatial POD solutions which isolate the growth of the primary and secondary instability mechanisms in the first and second modes, respectively. Temporal evolutions of dominant POD modes in all measured fields are calculated. These scalar POD coefficients contain the integrated characteristics of the entire field, greatly reducing the amount of data to characterize the instantaneous field. These modes may then be used to train future flow control algorithms based on neural networks.

Chapman, Keith↗

NASA NTRS1986

Coherent fluctuations of extratropical geopotential height and tropical convection in intraseasonal time scales

Northern-Hemisphere twice-daily 500-mb-geopotential-height data for November-March 1975-1980 and November-December 1981 and NOAA satellite measurements of outgoing longwave radiation for the same period are processed to remove seasonal cycles and Fourier-decomposition bandpassed to study variations with 20-70-d periods. The results of analysis using correlation, complex EOF, and composite techniques are presented in extensive maps and graphs and characterized. Extratropical wavetrains are found to evolve systematically from Eurasia eastward to North America and the North Atlantic on 5-6-d time scales, while the intraseasonal variation in tropical convection is dominated by a dipolelike east-west feature propagating from the western Indian Ocean to the dateline with a quasi-period of 40-50 d. The possibility that normal modes coupled between the tropics and midlatitudes may be responsible for these phenomena is considered.

Lau, K.-M.↗

NASA NTRS2001

Impacts of ENSO on the South American Summer Monsoon During 1997-1999

Using the National Center for Environmental Prediction (NCEP) Reanalysis, and CPC Merged Analysis Product (CMAP) rainfall, we have compared and contrasted the anomalies of the South American Summer Monsoon (SASM) during two extreme years of 1997/98 (EI Nino) and 1998/99 (La Nina). The results are assessed against a "canonical" ENSO response (CER) pattern for the SASM obtained from empirical mode decomposition based on a previous period (1980-1995). Overall, the SASM anomalies compare well with CER, but with some important differences. Anomalies occurring in the warm phase of the 1997-98 El Nino are very significant and robust, while those occurring in 1998/99 La Nina, appear to be reversed from 1997/98, but are relatively weak and less well-defined. The most pronounced signal in DJF 1997/98 is the development of drought conditions in northern Brazil, excessive rainfall over northern Peru and Ecuador, and over Uruguay and southern Brazil. The tropical rainfall anomalies are associated with the eastward shift of the Walker circulation, which is represented by pronounced low-level anomalous westerlies over the equatorial eastern Pacific and easterlies over northern Brazil. The easterlies are deflected sharply southeastward by the steep topography of the Andes, enhancing the low-level jet (LLJ) along the eastern foothills of the Andes near 15-20 S. The LLJ penetrates deep into the extratropics, yielding rainfall anomalies further poleward compared to CER. During DJF 1997/98, the eastward expansion of the warm tropospheric temperature over the Nino-3 region causes anomalous geopotential height to develop in the upper troposphere above the Altiplano, leading to a strengthened Bolivian High. An upper-tropospheric jet anomaly maximum is found over the subtropical continent near 30 S, due to increasing meridional gradient of tropospheric temperature, as well as teleconnection patterns linking the South Pacific and the South Atlantic. Consistent with the CER, the South Pacific High is weakened, and the South Atlantic High is strengthened in DJF 1997/98. However, rather than appearing as a coherent large scale signal over the entire Atlantic as in CER, the South Atlantic High anomalies in DJF 1997/98 split into two separate high pressure centers; one located south of the Amazon basin, and another over the southeastern South Atlantic. The former is induced as a Rossby-wave response by ENSO-Induced sinking motion over nor-them and eastern Brazil, and the latter may be associated with extratropical teleconnection signals. Anomalous convection observed over equatorial western Africa may also have contributed to the different circulation response over the Atlantic in DFJ 1997/99 compared to the CER. Comparison of the evolution of the SASM anomalies in 1997-99 with their corresponding annual cycles, suggest that the anomalies are phase-locked to and represent either enhancement or reduction of the annual variations. Results suggest that, from a system perspective the SASM was strengthened during the 1997/98 El Nino and weakened, albeit to a lesser extent, during the 1998/99 La Nina.

Lau, K.-M.↗