Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parabolized Stability Equations”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Compressible stability of growing boundary layers using parabolized stability equations

The parabolized stability equation (PSE) approach is employed to study linear and nonlinear compressible stability with an eye to providing a capability for boundary-layer transition prediction in both 'quiet' and 'disturbed' environments. The governing compressible stability equations are solved by a rational parabolizing approximation in the streamwise direction. Nonparallel flow effects are studied for both the first- and second-mode disturbances. For oblique waves of the first-mode type, the departure from the parallel results is more pronounced as compared to that for the two-dimensional waves. Results for the Mach 4.5 case show that flow nonparallelism has more influence on the first mode than on the second. The disturbance growth rate is shown to be a strong function of the wall-normal distance due to either flow nonparallelism or nonlinear interactions. The subharmonic and fundamental types of breakdown are found to be similar to the ones in incompressible boundary layers.

Chang, Chau-Lyan↗
NASA NTRS1994

Parabolized stability equations

The parabolized stability equations (PSE) are a new approach to analyze the streamwise evolution of single or interacting Fourier modes in weakly nonparallel flows such as boundary layers. The concept rests on the decomposition of every mode into a slowly varying amplitude function and a wave function with slowly varying wave number. The neglect of the small second derivatives of the slowly varying functions with respect to the streamwise variable leads to an initial boundary-value problem that can be solved by numerical marching procedures. The PSE approach is valid in convectively unstable flows. The equations for a single mode are closely related to those of the traditional eigenvalue problems for linear stability analysis. However, the PSE approach does not exploit the homogeneity of the problem and, therefore, can be utilized to analyze forced modes and the nonlinear growth and interaction of an initial disturbance field. In contrast to the traditional patching of local solutions, the PSE provide the spatial evolution of modes with proper account for their history. The PSE approach allows studies of secondary instabilities without the constraints of the Floquet analysis and reproduces the established experimental, theoretical, and computational benchmark results on transition up to the breakdown stage. The method matches or exceeds the demonstrated capabilities of current spatial Navier-Stokes solvers at a small fraction of their computational cost. Recent applications include studies on localized or distributed receptivity and prediction of transition in model environments for realistic engineering problems. This report describes the basis, intricacies, and some applications of the PSE methodology.

Herbert, Thorwald↗
NASA NTRS2016

Transient Growth Analysis of Compressible Boundary Layers with Parabolized Stability Equations

The linear form of parabolized linear stability equations (PSE) is used in a variational approach to extend the previous body of results for the optimal, non-modal disturbance growth in boundary layer flows. This methodology includes the non-parallel effects associated with the spatial development of boundary layer flows. As noted in literature, the optimal initial disturbances correspond to steady counter-rotating stream-wise vortices, which subsequently lead to the formation of stream-wise-elongated structures, i.e., streaks, via a lift-up effect. The parameter space for optimal growth is extended to the hypersonic Mach number regime without any high enthalpy effects, and the effect of wall cooling is studied with particular emphasis on the role of the initial disturbance location and the value of the span-wise wavenumber that leads to the maximum energy growth up to a specified location. Unlike previous predictions that used a basic state obtained from a self-similar solution to the boundary layer equations, mean flow solutions based on the full Navier-Stokes (NS) equations are used in select cases to help account for the viscous-inviscid interaction near the leading edge of the plate and also for the weak shock wave emanating from that region. These differences in the base flow lead to an increasing reduction with Mach number in the magnitude of optimal growth relative to the predictions based on self-similar mean-flow approximation. Finally, the maximum optimal energy gain for the favorable pressure gradient boundary layer near a planar stagnation point is found to be substantially weaker than that in a zero pressure gradient Blasius boundary layer.

Compressible boundary layer↗
NASA NTRS1993

Role of secondary instability theory and parabolized stability equations in transition modeling

In modeling the laminar-turbulent transition region, the designer depends largely on benchmark data from experiments and/or direct numerical simulations that are usually extremely expensive. An understanding of the evolution of the Reynolds stresses, turbulent kinetic energy, and quantifies in the transport equations like the dissipation and production is essential in the modeling process. The secondary instability theory and the parabolized stability equations method are used to calculate these quantities, which are then compared with corresponding quantities calculated from available direct numerical simulation data for the incompressible boundary-layer flow of laminar-turbulent transition conditions. The potential of the secondary instability theory and the parabolized stability equations approach in predicting these quantities is discussed; results indicate that inexpensive data that are useful for transition modeling in the early stages of the transition region can be provided by these tools.

El-Hady, Nabil M.↗
NASA NTRS1995

Weakly nonparallel and curvature effects on stationary crossflow instability: Comparison of results from multiple-scales analysis and parabolized stability equations

A multiple-scales approach is used to approximate the effects of nonparallelism and streamwise surface curvature on the growth of stationary crossflow vortices in incompressible, three-dimesional boundary layers. The results agree with results predicted by solving the parabolized stability equations in regions where the nonparallelism is sufficiently weak. As the nonparallelism increases, the agreement between the two approaches worsens. An attempt has been made to quantify the nonparallelism on flow stability in terms of a nondimensional number that describes the rate of change of the mean flow relative to the disturbance wavelength. We find that the above nondimensional number provides useful information about the adequacy of the multiple-scales approximation for different disturbances for a given flow geometry, but the number does not collapse data for different flow geometries onto a single curve.

Singer, Bart A.↗
NASA NTRS1992

Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

Joslin, Ronald D.↗
NASA NTRS1993

Stability of mixing layers

The research program for the first year of this project (see the original research proposal) consists of developing an explicit marching scheme for solving the parabolized stability equations (PSE). Performing mathematical analysis of the computational algorithm including numerical stability analysis and the determination of the proper boundary conditions needed at the boundary of the computation domain are implicit in the task. Before one can solve the parabolized stability equations for high-speed mixing layers, the mean flow must first be found. In the past, instability analysis of high-speed mixing layer has mostly been performed on mean flow profiles calculated by the boundary layer equations. In carrying out this project, it is believed that the boundary layer equations might not give an accurate enough nonparallel, nonlinear mean flow needed for parabolized stability analysis. A more accurate mean flow can, however, be found by solving the parabolized Navier-Stokes equations. The advantage of the parabolized Navier-Stokes equations is that its accuracy is consistent with the PSE method. Furthermore, the method of solution is similar. Hence, the major part of the effort of the work of this year has been devoted to the development of an explicit numerical marching scheme for the solution of the Parabolized Navier-Stokes equation as applied to the high-seed mixing layer problem.

Tam, Christopher↗
NASA NTRS2004

LASTRAC.3d: Transition Prediction in 3D Boundary Layers

Langley Stability and Transition Analysis Code (LASTRAC) is a general-purpose, physics-based transition prediction code released by NASA for laminar flow control studies and transition research. This paper describes the LASTRAC extension to general three-dimensional (3D) boundary layers such as finite swept wings, cones, or bodies at an angle of attack. The stability problem is formulated by using a body-fitted nonorthogonal curvilinear coordinate system constructed on the body surface. The nonorthogonal coordinate system offers a variety of marching paths and spanwise waveforms. In the extreme case of an infinite swept wing boundary layer, marching with a nonorthogonal coordinate produces identical solutions to those obtained with an orthogonal coordinate system using the earlier release of LASTRAC. Several methods to formulate the 3D parabolized stability equations (PSE) are discussed. A surface-marching procedure akin to that for 3D boundary layer equations may be used to solve the 3D parabolized disturbance equations. On the other hand, the local line-marching PSE method, formulated as an easy extension from its 2D counterpart and capable of handling the spanwise mean flow and disturbance variation, offers an alternative. A linear stability theory or parabolized stability equations based N-factor analysis carried out along the streamline direction with a fixed wavelength and downstream-varying spanwise direction constitutes an efficient engineering approach to study instability wave evolution in a 3D boundary layer. The surface-marching PSE method enables a consistent treatment of the disturbance evolution along both streamwise and spanwise directions but requires more stringent initial conditions. Both PSE methods and the traditional LST approach are implemented in the LASTRAC.3d code. Several test cases for tapered or finite swept wings and cones at an angle of attack are discussed.

Chang, Chau-Lyan↗
NASA NTRS2016

Direct Numerical Simulation of Transition Due to Traveling Crossflow Vortices

Previous simulations of laminar breakdown mechanisms associated with stationary crossflow instability over a realistic swept-wing configuration are extended to investigate the alternate scenario of transition due to secondary instability of traveling crossflow modes. Earlier analyses based on secondary instability theory and parabolized stability equations have shown that this alternate scenario is viable when the initial amplitude of the most amplified mode of the traveling crossflow instability is greater than approximately 0.03 times the initial amplitude of the most amplified stationary mode. The linear growth predictions based on the secondary instability theory and parabolized stability equations agree well with the direct numerical simulation. Nonlinear effects are initially stabilizing but subsequently lead to a rapid growth followed by the onset of transition when the amplitude of the secondary disturbance exceeds a threshold value. Similar to the breakdown of stationary vortices, the transition zone is rather short and the boundary layer becomes completely turbulent across a distance of less than 15 times the boundary layer thickness at the completion of transition.

Li, Fei↗
NASA NTRS1993

Linear and nonlinear PSE for compressible boundary layers

Compressible stability of growing boundary layers is studied by numerically solving the partial differential equations under a parabolizing approximation. The resulting parabolized stability equations (PSE) account for nonparallel as well as nonlinear effects. Evolution of disturbances in compressible flat-plate boundary layers are studied for freestream Mach numbers ranging from 0 to 4.5. Results indicate that the effect of boundary-layer growth is important for linear disturbances. Nonlinear calculations are performed for various Mach numbers. Two-dimensional nonlinear results using the PSE approach agree well with those from direct numerical simulations using the full Navier-Stokes equations while the required computational time is less by an order of magnitude. Spatial simulation using PSE were carried out for both the fundamental and subharmonic type breakdown for a Mach 1.6 boundary layer. The promising results obtained show that the PSE method is a powerful tool for studying boundary-layer instabilities and for predicting transition over a wide range of Mach numbers.

Chang, Chau-Lyan↗
NASA NTRS2021

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗
NASA NTRS2025

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗
NASA NTRS2021

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗
NASA NTRS2022

Stability Analysis of Streaks Induced by Optimized Vortex Generators

Numerical computations are performed to investigate the potential for transition control in an axisymmetric boundary layer via fully realizable, streamwise stationary streaks induced by an azimuthally periodic array of surface mounted vortex generators (VGs). Previous work has shown that suitable streaks of this type can significantly reduce the growth of Mack’s second mode instabilities, but large streak amplitudes can make the flow susceptible to previously absent streak instabilities that can become the leading cause of transition. Here, we use the adjoint capabilities of the SU2 flow solver to optimize the VG shape to maximize the reduction in the growth of second-mode disturbances while also preventing the streak amplitudes from reaching large enough values to precipitate an earlier onset of transition via streak instabilities. The geometry and the freestream flow conditions are selected to match a relevant trajectory lo-cation from the HIFiRE-1 flight experiment. Results show that the optimized VGs can increase the mean streak amplitude by 117% with respect to a manually developed baseline design. The stability of this optimized basic state is analyzed via the plane-marching parabolized stability equations, predicting a fully laminar flow over the entire cone, or equivalently, yielding transition delay of 130% versus the 17% for the baseline VGs.

Boundary layer transition↗
NASA NTRS2002

Stochastic Modeling of Laminar-Turbulent Transition

Stochastic versions of stability equations are developed in order to develop integrated models of transition and turbulence and to understand the effects of uncertain initial conditions on disturbance growth. Stochastic forms of the resonant triad equations, a high Reynolds number asymptotic theory, and the parabolized stability equations are developed.

Rubinstein, Robert↗
NASA NTRS2022

Combined Bluntness and Roughness Effects on Cones at Hypersonic Speeds

This computational study investigates the effects of discrete roughness elements on a blunt cone at zero degrees angle of attack in a Mach 6 flow. Motivation was provided by experiments conducted in the Air Force Research Laboratory Mach 6 High Reynolds Number facility on a 7-degree half-angle cone with a roughness array located at 45 degrees from the apex on two nose tips of different blutness but equivalent roughness Reynolds number. Transition was only affected on the blunter cone, indicating that the transition onset is associated with the combined effects of bluntness and roughness. The present study investigates the 15.24 mm nose radius, 420 azimuthal wavenumber case via Navier-Stokes computations of the laminar base flow and instability analysis. Plane-marching parabolized stability equations (PSE) and inflow-resolvent analysis based on the three-dimensional, harmonic linearized Navier-Stokes equations (HLNSE) are used to calculate the amplification of disturbances along the roughness wake as well as over the roughness nearfield. Results show that the roughness shape can have a great impact on the characteristics of the most amplified wake instabilities. For the experimental configuration with cubic roughness elements of 15 μ m height, the flow is marginally unstable. For prismatic elements of 20 μ m height, the PSE predicts a logarithmic disturbance amplification ratio of N = 5.6 along its wake, but this ratio increases to N = 9.6 when the amplification over the roughness and separation regions is included in the inflow-resolvent analysis.

Boundary-layer transition↗
NASA NTRS2021

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper reports the results of a comprehensive linear stability analysis of the boundary layer flow over the common research model with natural laminar flow (CRM-NLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien-Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. However, a significant scatter ( N ∈ (4,9)) was observed in the values of the logarithmic amplification factor along the measured transition front. This previous analysis was based on an approximate basic state (based on a boundary layer code with conical flow approximation) and parallel stability computations without surface curvature effects. Here, we examine the effects of these various approximations with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, both linear stability theory (LST) and the parabolized stability equations (PSE)are used in conjunction with an accurate definition of the laminar boundary layer flow as computed with a Navier-Stokes solver with a Reynolds-Averaged-Navier-Stokes (RANS) based turbulence model within the turbulent parts of the flow. Furthermore, the effects of instability wave propagation within a fully three-dimensional boundary layer are also evaluated by integrating the disturbance growth rates along suitably chosen, curvilinear (i.e., nonplanar) propagation trajectories. The results of this analysis are also used in an accompanying paper by Venkatachari et al. to develop improved, physics based transition predictions for the same CRM-NLF configuration.

Boundary layer transition↗
NASA NTRS1995

Linear Wave Motion from Concentrated Harmonic Sources in Blasius Flow

The motion of individual linear instability waves in shear flows is well described by existing theoretical and numerical methods. However, naturally occuring sources produce coherent wave motions with broadband spanwise wavenumber and frequency spectra, and the different spectral components interact both linearly and nonlinearily. This paper describes a series of calculations for the parameters of three different experiments using locally parallel linear stability theory (LST), the Parabolized Stability Equations (PSE), and Direct Numerical Simulation (DNS). The calculations illustrate the strengths and weaknesses of the different methods, the extent to which the methods agree or disagree, and, finally the extent to which agreement with the measurements can be attained, given that the experiments also have their own difficulties.

linear↗