SEARCH · Search NASA

Results for “PARABOLIC EQUATION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Simulation of crossflow instability on a supersonic highly swept wing

A direct numerical simulation (DNS) algorithm has been developed and validated for use in the investigation of crossflow instability on supersonic swept wings, an application of potential relevance to the design of the High-Speed Civil Transport (HSCT). The algorithm is applied to the investigation of stationary crossflow instability on an infinitely long 77-degree swept wing in Mach 3.5 flow. The results of the DNS are compared with the predictions of linear parabolized stability equation (PSE) methodology. In-general, the DNS and PSE results agree closely in terms of modal growth rate, structure, and orientation angle. Although further validation is needed for large-amplitude (nonlinear) disturbances, the close agreement between independently derived methods offers preliminary validation of both DNS and PSE approaches.

Pruett, C. David↗

NASA NTRS1996

Nonlinear stability of supersonic jets

The stability calculations made for a shock-free supersonic jet using the model based on parabolized stability equations are presented. In this analysis the large scale structures, which play a dominant role in the mixing as well as the noise radiated, are modeled as instability waves. This model takes into consideration non-parallel flow effects and also nonlinear interaction of the instability waves. The stability calculations have been performed for different frequencies and mode numbers over a range of jet operating temperatures. Comparisons are made, where appropriate, with the solutions to Rayleigh's equation (linear, inviscid analysis with the assumption of parallel flow). The comparison of the solutions obtained using the two approaches show very good agreement.

Tiwari, S. N.↗

NASA NTRS1995

Nonlinear Stability of Supersonic Jets

This paper presents stability calculations made for a shock-free supersonic jet using the model based on parabolized stability equations. In this analysis the large-scale structures, which play a dominant role in the mixing as well as the noise radiated, are modeled as instability waves. This model takes into consideration non-parallel flow effects and also nonlinear interaction of the instability waves. The stability calculations have been performed for different frequencies and mode numbers over a range of jet operating temperatures. Comparisons are made, where appropriate, with the solutions to Rayleigh's equation (linear, inviscid analysis with the assumption of parallel flow). The comparison of the solutions obtained using the two approaches show very good agreement.

Bhat, T. R. S.↗

NASA NTRS1998

Effects of Suction on Swept-Wing Transition

Stability experiments are conducted in the Arizona State University Unsteady Wind Tunnel on a 45 deg swept airfoil. The pressure gradient is designed to provide purely crossflow-dominated transition; that is, the boundary layer is subcritical to Tollmien-Schlichting disturbances. The airfoil surface is hand polished to a 0.25 microns rms finish. Under these conditions, stationary crossflow disturbances grow to nonuniform amplitude due to submicron surface irregularities near the leading edge. Uniform stationary crossflow waves are produced by controlling the initial conditions with spanwise arrays of micron-sized roughness elements near the attachment line. Hot-wire measurements provide detailed maps of the crossflow wave structure, and accurate spectral decompositions isolate individual-mode growth rates for the fundamental and harmonic disturbances. Roughness spacing, roughness height, and Reynolds number are varied to investigate the growth of all amplified wavelengths. The measurements show early nonlinear mode interaction causing amplitude saturation well before transition. Comparisons with nonlinear parabolized stability equations calculations show excellent agreement in both the disturbance amplitude and the mode-shape profiles.

Saric, William S.↗

NASA NTRS1998

Transition in Turbine Flows

We have further developed our capabilities to analyze transition in turbine boundary layers from first principles by integrating the nonlinear parabolized stability equations (PSE) with improved initial and boundary conditions. With modified iteration schemes, we are able to proceed deeper into the transition region where skin friction coefficient and heat transfer coefficient significantly increase. Initial and boundary conditions at elevated turbulence levels can be derived by receptivity analysis. Test runs for ERCOFTAC test case T3A at 2.4\% turbulence level provide results in good agreement with the experimental data. The sharper minimum of the skin coefficient also shown by DNS results is likely due to the missing intermittency. The method has been applied to various experimentally studied turbine blades (UTRC, VKI, Zierke, Langston, Hippensteele, and others). The PSE results, though physically reasonable, do not agree as well as expected with the experimental findings. We have, therefore, performed an extensive search for the reasons of the seemingly systematic deviations. A first source of uncertainty has been found in the often insufficient documentation of the experiments (e.g. on blockage by end-wall boundary layers). However, variation of the relevant parameters does not lead to more satisfactory agreement. A second reason has been found in the "standard procedure" which considers a 2D flow at midspan and uses a panel code and subsequent boundary-layer code to obtain the laminar basic flow for the transition analysis. Comparison with the pressure distribution obtained with a 3D design code (RVC3D) shows significant three-dimensionality of the flow (e.g. in the UTRC experiments). The spanwise variation has been neglected in our original PSE code. To overcome this problem, we have developed the PSE/3D for fully 3D boundary layers to account for streamwise and spanwise variations. Since the design code does not provide the boundary-layer flow with sufficient resolution, we have generated the Euler solution and employed a 3D boundary-layer code to obtain the viscous basic flow. Although only the linear stability level of PSE/3D has been implemented so far, the discrepancies with the experiments change but do not disappear. We still find deviations between the computed and experimental variations of C(sub f), and St along the blade for laminar flow. The main reason can be seen by comparing the solution of the boundary-layer code with the viscous results of the design code. The conventional boundary-layer solution exhibits an asymptotic behavior appropriate in external aerodynamics but does not match the steep gradients of the inviscid flow through the passage and consequently provides biased results for C(sub f), and St. An attempt is currently being made to correct this deficiency. Before attempting to perform the transition analysis for the viscous flow provided by the design code, we have analyzed the implementation and "best possible" results. Code and results exhibit flaws that may negatively affect the design and are intolerable for transition analysis. Therefore, we have decided to develop a new code to obtain a reliable basis for stability and transition studies. We expect to report improved results by the time of the meeting.

Herbert, Thorwald↗

NASA NTRS1999

Subsonic and Supersonic Jet Noise Calculations Using PSE and DNS

Noise radiated from a supersonic jet is computed using the Parabolized Stability Equations (PSE) method. The evolution of the instability waves inside the jet is computed using the PSE method and the noise radiated to the far field from these waves is calculated by solving the wave equation using the Fourier transform method. We performed the computations for a cold supersonic jet of Mach number 2.1 which is excited by disturbances with Strouhal numbers St=.2 and .4 and the azimuthal wavenumber m=l. Good agreement in the sound pressure level are observed between the computed and the measured (Troutt and McLaughlin 1980) results.

Balakumar, P.↗

NASA NTRS1999

On the Stability of Three-Dimensional Boundary Layers: Linear and Nonlinear Stability - Part 1

The primary stability of incompressible three-dimensional boundary layers is investigated using the Parabolized Stability Equations (PSE). We compute the evolution of stationary and traveling disturbances in the linear and nonlinear region prior to transition. As model problems, we choose Swept Hiemenz Flow and the DLR Transition Experiment. The primary stability results for Swept Hiemenz Flow agree very well with computations by Malik et al. For the DLR Experiment, the mean flow profiles are obtained by solving the boundary layer equations for the measured pressure distribution. Both linear and nonlinear results show very good agreement with the experiment.

Janke, Erik↗

NASA NTRS1999

Unified CFD Methods Via Flowfield-Dependent Variation Theory

This paper addresses the flowfield-dependent variation (FDV) methods in which complex physical phenomena are taken into account in the final form of partial differential equations to be solved so that finite difference methods (FDM) or finite element methods (FEM) themselves will not dictate the physics, but rather are no more than simply the options how to discretize between adjacent nodal points or within an element. The variation parameters introduced in the formulation are calculated from the current flowfield based on changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkohler numbers between adjacent nodal points, which play many significant roles such as adjusting the governing equations (hyperbolic, parabolic, and/or e!liptic), resolving various physical phenomena, and controlling the accuracy and stability of the numerical solution. The theory is verified by a number of example problems addressing the physical implications of the variation parameters which resemble the flowfield itself, shock capturing mechanism, transitions and interactions between inviscid/viscous, compressibility/incompressibility, and laminar/turbulent flows.

Chung, T. J.↗

NASA NTRS1998

Angular Broadening: Effects of Nonzero, Spatially Varying Plasma Frequency Between the Source and Observer

Angular broadening of radiation due to scattering by density irregularities is usually described using geometric optics (GO) or the parabolic wave equation (PWE) with the assumptions that the radiation frequency f greatly exceeds the local plasma frequency f(sub p0) or that f(sub p0)/f is constant along the path. These assumptions are inappropriate for many solar system radio phenomena. Here the PWE and GO formalisms are extended to treat angular broadening in plasmas with nonzero, spatially varying ratios, f(sub p0)(z)/f < 1. The new PWE results show that the correlation function, scattered angular spectrum, and other quantities are modified by inclusion of a denominator factor [1 - f(sup 2, sub p0)(z prime)/f prime] inside the path integral over z prime, while the mean-square scattering angle depends on both f(sub p0)(z)/f at the observer and the foregoing factor inside the path integral. The PWE and GO predictions for are identical and involve equivalent assumptions. Previous GO and PWE results are recovered in the limits that f(sub p0)(z prime)/f is constant or zero. The new PWE and GO results will permit more accurate calculation of angular broadening for solar system and astrophysical sources. Moreover and importantly, due to the PWE and GO results for being identical, previous GO analyses of in solar system contexts are essentially correct, except for the neglect of or minor deficiencies in the treatment of nonzero, spatially varying f(sub p0)/f effects. The identical GO and PWE results for and the form of the PWE equation for the correlation function raise questions as to whether diffraction is unimportant for angular broadening (under the usual PWE conditions). Future direct comparisons of the PWE predictions with angular spectra calculated using existing GO ray-tracing codes should answer these questions. Diffraction effects are probably important when the medium and turbulence are not sufficiently homogeneous transverse to the central ray path. The implications are also discussed for studies of the 2-3.5 kHz radiation observed in the outer heliosphere.

Cairns, Iver H.↗

NASA NTRS2003

The Langley Stability and Transition Analysis Code (LASTRAC) : LST, Linear and Nonlinear PSE for 2-D, Axisymmetric, and Infinite Swept Wing Boundary Layers

During the past two decades, our understanding of laminar-turbulent transition flow physics has advanced significantly owing to, in a large part, the NASA program support such as the National Aerospace Plane (NASP), High-speed Civil Transport (HSCT), and Advanced Subsonic Technology (AST). Experimental, theoretical, as well as computational efforts on various issues such as receptivity and linear and nonlinear evolution of instability waves take part in broadening our knowledge base for this intricate flow phenomenon. Despite all these advances, transition prediction remains a nontrivial task for engineers due to the lack of a widely available, robust, and efficient prediction tool. The design and development of the LASTRAC code is aimed at providing one such engineering tool that is easy to use and yet capable of dealing with a broad range of transition related issues. LASTRAC was written from scratch based on the state-of-the-art numerical methods for stability analysis and modem software technologies. At low fidelity, it allows users to perform linear stability analysis and N-factor transition correlation for a broad range of flow regimes and configurations by using either the linear stability theory (LST) or linear parabolized stability equations (LPSE) method. At high fidelity, users may use nonlinear PSE to track finite-amplitude disturbances until the skin friction rise. Coupled with the built-in receptivity model that is currently under development, the nonlinear PSE method offers a synergistic approach to predict transition onset for a given disturbance environment based on first principles. This paper describes the governing equations, numerical methods, code development, and case studies for the current release of LASTRAC. Practical applications of LASTRAC are demonstrated for linear stability calculations, N-factor transition correlation, non-linear breakdown simulations, and controls of stationary crossflow instability in supersonic swept wing boundary layers.

Chang, Chau-Lyan↗

NASA NTRS2004

Langley Stability and Transition Analysis Code (LASTRAC) Version 1.2 User Manual

LASTRAC is a general-purposed, physics-based transition prediction code released by NASA for Laminar Flow Control studies and transition research. The design and development of the LASTRAC code is aimed at providing an engineering tool that is easy to use and yet capable of dealing with a broad range of transition related issues. It was written from scratch based on the state-of-the-art numerical methods for stability analysis and modern software technologies. At low fidelity, it allows users to perform linear stability analysis and N-factor transition correlation for a broad range of flow regimes and configurations by using either the linear stability theory or linear parabolized stability equations method. At high fidelity, users may use nonlinear PSE to track finite-amplitude disturbances until the skin friction rise. This document describes the governing equations, numerical methods, code development, detailed description of input/output parameters, and case studies for the current release of LASTRAC.

Chang, Chau-Lyan↗

NASA NTRS1996

Simulation of Three-Dimensional Symmetric and Asymmetric Instabilities in Attachment-Line Boundary Layers

On a swept wing, contamination along the leading edge, Tollmien-Schlichting waves, stationary or traveling crossflow vortices, and/or Taylor-Gortler vortices can cause the catastrophic breakdown of laminar to turbulent flow, which leads to increased skin-friction drag for the aircraft. The discussion in this Note will be limited to disturbances which evolve along the attachment line (leading edge of swept wing). If the Reynolds number of the attachment-line boundary layer is greater than some critical value, then the complete wing is inevitably engulfed in turbulent flow. Essentially, there are two critical Reynolds number points that must be considered. The first is for small-amplitude disturbances, and the second is for bypass transition. The present study will use direct numerical simulations to validate a linear 2D-eigenvalue prediction method based on parabolized stability equations by Lin and Malik. This method is considered because it suggests that a number of symmetric and asymmetric modes exist and are stable or unstable on the attachment line depending on the Reynolds number. If validated, the approach would predict a number of modes which are linearly damped in the Reynolds number regime 100 to 245; however, these modes may grow nonlinearly and provide an explanation to this region.

Joslin, Ronald D.↗

NASA NTRS1995

Effect of Far-Field Boundary Conditions on Boundary-Layer Transition

The effect of far-field boundary conditions on the evolution of a finite-amplitude two-dimensional wave in the Blasius boundary layer is assessed. With the use of the parabolized stability equations (PSE) theory for the numerical computations, either asymptotic, Dirichlet, Neumann or mixed boundary conditions are imposed at various distances from the wall. The results indicate that asymptotic and mixed boundary conditions yield the most accurate mean-flow distortion and unsteady instability modes in comparison with the results obtained with either Dirichlet or Neumann conditions.

Bertolotti, Fabio P.↗

NASA NTRS1994

Effect of Far-Field Boundary Conditions on Boundary-Layer Transition

Bertolotti, Fabio P.↗

NASA NTRS2000

Nonlinear Stability and Structure of Compressible Reacting Mixing Layers

The parabolized stability equations (PSE) are used to investigate issues of nonlinear flow development and mixing in compressible reacting shear layers. Particular interest is placed on investigating the change in flow structure that occurs when compressibility and heat release are added to the flow. These conditions allow the 'outer' instability modes- one associated with each of the fast and slow streams-to dominate over the 'central', Kelvin-Helmholtz mode that unaccompanied in incompressible nonreacting mixing layers. Analysis of scalar probability density functions in flows with dominant outer modes demonstrates the ineffective, one-sided nature of mixing that accompany these flow structures. Colayer conditions, where two modes have equal growth rate and the mixing layer is formed by two sets of vortices, offer some opportunity for mixing enhancement. Their extent, however, is found to be limited in the mixing layer's parameter space. Extensive validation of the PSE technique also provides a unique perspective on central- mode vortex pairing, further supporting the view that pairing is primarily governed perspective sheds insight on how linear stability theory is able to provide such an accurate prediction of experimentally-observed, fully nonlinear flow phenomenon.

Day, M. J.↗

NASA NTRS2007

Simulations of Boundary-Layer Transition

For incompressible benchmark flows, we have demonstrated the capability of the parabolized stability equations (PSE) to simulate the transition process in excellent agreement with microscopic experiments and direct Navier-Stokes simulations at modest computational cost. Encouraged by these results, we have developed the PSE methodology of three-dimensional boundary-layers in general curvilinear coordinates for the range from low to hypersonic speeds, and for both linear and nonlinear problems. For given initial and boundary conditions, the approach permits simulations from receptivity through linear and secondary instabilities into the late stages of transition where significant changes in skin friction and heat transfer coefficients occur. We have performed transition simulations for a variety of two- and three-dimensional similarity solutions and for realistic flows over swept wings at subsonic and supersonic speeds, the pressure ans suction side of turbine blades at low and medium turbulence levels, and over a blunt cone at Mach number Ma = 8. We present selected results for different transition mechanisms with emphasis on the late stage of transition and the evolution of wall-shear stress and heat transfer.

Herbert, Thorwald↗

NASA NTRS2005

Distributed Blowing and Suction for the Purpose of Streak Control in a Boundary Layer Subjected to a Favorable Pressure Gradient

An analysis of the optimal control by blowing and suction in order to generate stream- wise velocity streaks is presented. The problem is examined using an iterative process that employs the Parabolized Stability Equations for an incompressible uid along with its adjoint equations. In particular, distributions of blowing and suction are computed for both the normal and tangential velocity perturbations for various choices of parameters.

Forgoston, Eric↗

NASA NTRS2009

Roughness Based Crossflow Transition Control: A Computational Assessment

A combination of parabolized stability equations and secondary instability theory has been applied to a low-speed swept airfoil model with a chord Reynolds number of 7.15 million, with the goals of (i) evaluating this methodology in the context of transition prediction for a known configuration for which roughness based crossflow transition control has been demonstrated under flight conditions and (ii) of analyzing the mechanism of transition delay via the introduction of discrete roughness elements (DRE). Roughness based transition control involves controlled seeding of suitable, subdominant crossflow modes, so as to weaken the growth of naturally occurring, linearly more unstable crossflow modes. Therefore, a synthesis of receptivity, linear and nonlinear growth of stationary crossflow disturbances, and the ensuing development of high frequency secondary instabilities is desirable to understand the experimentally observed transition behavior. With further validation, such higher fidelity prediction methodology could be utilized to assess the potential for crossflow transition control at even higher Reynolds numbers, where experimental data is currently unavailable.

Li, Fei↗