Search NASA⌕ Search

SEARCH · Search NASA

Results for “Highly oscillatory solutions”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Receptivity of Supersonic Boundary Layers Due To Acoustic Disturbances Over Blunt Cones

Receptivity and stability of supersonic boundary layers over a 5-degree straight cone with a blunt tip are numerically investigated at a free stream Mach number of 3.5 and at a high Reynolds number of 106/inch. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The linear stability results showed that bluntness has less stabilizing effects on the stability of boundary layers over cones than on flat plates and wedges. The unsteady simulations of the interaction of plane threedimensional acoustic waves with the cone showed that the modulation of wavelength and the generation of instability waves first occurred near the leading edge in the plane where the constant acoustic phase lines are perpendicular to the cone axis. Further downstream, this instability region spreads in the azimuthal direction from this plane.

Balakumar, P.↗
NASA NTRS2009

High-Order Energy Stable WENO Schemes

A third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme developed by Yamaleev and Carpenter was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables 'energy stable' modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes; ESWENO schemes up to eighth order are presented in the appendix. New weight functions are also developed that provide (1) formal consistency, (2) much faster convergence for smooth solutions with an arbitrary number of vanishing derivatives, and (3) improved resolution near strong discontinuities.

Yamaleev, Nail K.↗
NASA NTRS2009

A Systematic Methodology for Constructing High-Order Energy Stable WENO Schemes

A third-order Energy Stable Weighted Essentially Non{Oscillatory (ESWENO) finite difference scheme developed by Yamaleev and Carpenter [1] was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables "energy stable" modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes; ESWENO schemes up to eighth order are presented in the appendix. New weight functions are also developed that provide (1) formal consistency, (2) much faster convergence for smooth solutions with an arbitrary number of vanishing derivatives, and (3) improved resolution near strong discontinuities.

Yamaleev, Nail K.↗
NASA NTRS2006

Stability of Supersonic Boundary Layers Over Blunt Wedges

Receptivity and stability of supersonic boundary layers over blunt flat plates and wedges are numerically investigated at a free stream Mach number of 3.5 and at a high Reynolds number of 10(exp 6)/inch. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. Computations are performed for a flat plate with leading edge thicknesses of 0.0001, 0.001, 0.005 and 0.01 inches that give Reynolds numbers based on the leading edge thickness ranging from 1000 to 10000. Calculations are also performed for a wedge of 10 degrees half angle with different leading edge radii 0.001 and 0.01 inches. The linear stability results showed that the bluntness has a strong stabilizing effect on the stability of two-dimensional boundary layers. The transition Reynolds number for a flat plate with a leading edge thickness of 0.01 inches is about 3.5 times larger than it is for the Blasius boundary layer. It was also revealed that boundary layers on blunt wedges are far more stable than on blunt flat plates.

Balakumar, Ponnampalam↗
NASA NTRS2008

Boundary Layer Receptivity due to Roughness and Freestream Sound for Supersonic Flows over Axisymmetric Cones

The receptivity of supersonic boundary layers due to the interaction of a two-dimensional acoustic wave with a three-dimensional roughness is numerically investigated over a 5-degree sharp tipped straight cone at a free stream Mach number of 3.5 and a high Reynolds number of 10(exp 6)/inch. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for spatial discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for temporal integration. The simulations are performed with different roughness heights and roughness locations. The steady computations with the roughness reveal that the perturbations induced by the roughness remain very close to the wall and do not penetrate to the outer part of the boundary layer. The effects of the roughness persist for about 10 boundary layer thicknesses downstream of each element. The unsteady simulation shows the generation of unstable three-dimensional first mode disturbances within a short distance from the roughness. However, the amplitudes of the instability waves are very small near the neutral point, on the order of 10(exp -4) compared to the free stream disturbances. The receptivity coefficients due to the roughness elements located very close to the tip are one order of magnitude smaller than those due to the roughness elements located closer to the neutral point.

Balakumar, P.↗
NASA NTRS1993

Developmental testing resulting in a simplified liquid oxygen check valve for the Space Shuttle Main Propulsion System

The coil spring in a Space Shuttle liquid oxygen check valve failed due to cyclic fatigue in September, 1991. The dual-flapper, swing check valve is used to prevent reverse flow to the Space Shuttle Main Engines. Upon inspection of the failed component, the spring tangs were found to be missing and heavy wear was observed on the inner diameter of the spring coils. The fracture surfaces revealed that the metal had been steadily worn away until a simple overload caused the final fracture. A series of flow tests using water and a water/gas mixture was conducted to identify the flow phenomenon which produced the cyclic wear. A Plexiglas outlet housing was utilized to view the flapper behavior under different flow conditions and to aid in high speed photography. The tests revealed that flow instabilities induced two oscillatory flapper responses: a rocking mode and a chattering mode. Initially, attempts were made to reduce the spring-flapper oscillations. However, the final solution to the problem was a springless configuration which satisfied the valve's design requirements and eliminated the oscillations. The springless design relied on the inherent ability of the reverse flow momentum to close the flappers.

Aber, Gregory S.↗
NASA NTRS1985

Uniformly high-order accurate non-oscillatory schemes, 1

The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws was begun. These schemes share many desirable properties with total variation diminishing schemes (TVD), but TVD schemes have at most first order accuracy, in the sense of truncation error, at extreme of the solution. A uniformly second order approximation was constucted, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell.

Harten, A.↗
NASA NTRS2008

A Systematic Methodology for Constructing High-Order Energy-Stable WENO Schemes

A third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme developed by Yamaleev and Carpenter (AIAA 2008-2876, 2008) was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables \energy stable" modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes; ESWENO schemes up to eighth order are presented in the appendix. New weight functions are also developed that provide (1) formal consistency, (2) much faster convergence for smooth solutions with an arbitrary number of vanishing derivatives, and (3) improved resolution near strong discontinuities.

Yamaleev, Nail K.↗
NASA NTRS1979

The promise of multicyclic control

Several types of rotors which employ multicyclic control are reviewed and compared. Their differences are high-lighted and their potential advantages and disadvantages are discussed. The flow field these rotors must operate in is discussed, and it is shown that simultaneous elimination of vibration and oscillatory blade loads is not an inherent solution to the roughness problem. The use of rotor blades as energy absorbers is proposed. Input-output relations are considered and a gain control for ROMULAN, a multicyclic controlling computer program, is introduced. Implications of the introduction of multicyclic systems into helicopters are discussed.

Mccloud, J. L., III↗
NASA NTRS1985

Recovering pointwise values of discontinuous data within spectral accuracy

The pointwise values of a function, f(x), can be accurately recovered either from its spectral or pseudospectral approximations, so that the accuracy solely depends on the local smoothness of f in the neighborhood of the point x. Most notably, given the equidistant function grid values, its intermediate point values are recovered within spectral accuracy, despite the possible presence of discontinuities scattered in the domain. (Recall that the usual spectral convergence rate decelerates otherwise to first order, throughout). To this end, a highly oscillatory smoothing kernel is employed in contrast to the more standard positive unit-mass mollifiers. In particular, post-processing of a stable Fourier method applied to hyperbolic equations with discontinuous data, recovers the exact solution modulo a spectrally small error. Numerical examples are presented.

Gottlieb, D.↗
NASA NTRS1985

Recovering pointwise values of discontinuous data within spectral accuracy

The pointwise values of a function, f(x), can be accurately recovered either from its spectral or pseudospectral approximations, so that the accuracy solely depends on the local smoothness of f in the neighborhood of the point x. Most notably, given the equidistant function grid values, its intermediate point values are recovered within spectral accuracy, despite the possible presence of discontinuities scattered in the domain. (Recall that the usual spectral convergence rate decelerates otherwise to first order, throughout). To this end, a highly oscillatory smoothing kernel is employed in contrast to the more standard positive unit-mass mollifiers. In particular, post-processing of a stable Fourier method applied to hyperbolic equations with discontinuous data, recovers the exact solution modulo a spectrally small error. Numerical examples are presented.

Gottlieb, D.↗
NASA NTRS1988

Non-oscillatory spectral Fourier methods for shock wave calculations

A non-oscillatory spectral Fourier method is presented for the solution of hyperbolic partial differential equations. The method is based on adding a nonsmooth function to the trigonometric polynomials which are the usual basis functions for the Fourier method. The high accuracy away from the shock is enhanced by using filters. Numerical results confirm that no oscillations develop in the solution. Also, the accuracy of the spectral solution of the inviscid Burgers equation is shown to be higher than a fixed order.

Cai, Wei↗
NASA NTRS1989

Recent progress on essentially non-oscillatory shock capturing schemes

An account is given of the construction of efficient implementations of 'essentially nonoscillatory' (ENO) schemes that approximate systems of hyperbolic conservation laws. ENO schemes use a local adaptive stencil to automatically obtain information from regions of smoothness when the solution develops discontinuities. Approximations employing ENOs can thereby obtain uniformly high accuracy to the very onset of discontinuities, while retaining a sharp and essentially nonoscillatory shock transition. For ease of implementation, ENO schemes applying the adaptive stencil concept to the numerical fluxes and employing a TVD Runge-Kutta-type time discretization are constructed.

Osher, Stanley↗
NASA NTRS1977

Influence of transverse surface waves on turbulent boundary layers

Wavy wall experiments using solid waves and progressive waves have been reported. For this paper, the major effects of waviness of the wall on the flow are identified as due to oscillatory curvature (convex-concavity) and oscillatory acceleration/deceleration of the flow, which imposes a highly nonequilibrium influence upon the turbulence structure. The theoretical analysis in this presentation takes into account proper turbulence modeling (including the nonequilibrium effects) for the wavy wall problem. The analysis proceeds in three stages: (1) inviscid solution for induced pressure due to the physical wall, (2) solution of a turbulent boundary layer with pressure gradients and curvature effects in the modeling from which the profile correction is computed, and (3) induced pressure computations for the corrected profile. The phase shift of pressure perturbations with respect to the physical wall can be predicted, and pressure drag and skin friction drag can be estimated, with nonlinear viscous effects included. Comparison of the theoretical estimates with experimental data are also presented.

Balasubramanian, R.↗
NASA NTRS1990

Flux-difference split scheme for turbulent transport equations

This paper describes a high accuracy upwind method to solve the mean compressible flow equations closed with two-equation turbulence models. A flux-difference splitting scheme is developed to compute the flowfield without introducing oscillations into the solution. A high-order accuracy scheme is used to predict both the mean flow variables and the turbulence variables. Results comparing two different two-equation turbulence models are presented for supersonic flow over a flat plate and a compressible, free shear-layer. The non-oscillatory behavior of the upwind scheme is demonstrated on the free shear-layer.

Morrison, J. H.↗
NASA NTRS1987

Uniformly high order accurate essentially non-oscillatory schemes. III

In the present continuation of development and analysis efforts towards essentially nonoscillatory shock-capturing approximations of hyperbolic conservation laws, the Godunov (1959) scheme and its MUSCL second-order-accurate extension are generalized to an arbitrary order of accuracy by means of a hierarchy of uniformly high-order-accurate schemes. The design method, which involves the solution's essentially nonoscillatory piecewise polynomial reconstruction on the basis of its cell averages, followed by time-evolution through an approximate solution of the resulting initial value problem and its averaging over each cell, employs an adaptive stencil of grid points to yield schemes that are highly nonlinear.

Harten, Ami↗
NASA NTRS1994

Convective flow effects on protein crystal growth

The long-term stability of the interferometric setup for the monitoring of protein morphologies has been improved. Growth or dissolution of a crystal on a 100 A scale can now be clearly distinguished from dimensional changes occurring within the optical path of the interferometer. This capability of simultaneously monitoring the local interfacial displacement at several widely-spaced positions on the crystal surface with high local depth resolution, has already yielded novel results. We found with lysozyme that (1) the normal growth rate is oscillatory, and (2) the mean growth step density is greater at the periphery of a facet than in its center. The repartitioning of Na(+) and Cl(-) ions between lysozyme solutions and crystals was studied for a wide range of crystallization conditions. A nucleation-growth-repartitioning model was developed to interpret the large body of data in a unified way. The results strongly suggests that (1) the ion to lysozyme ratio in the crystal depends mostly on kinetic rather than crystallographic parameters, and (2) lysozyme crystals possess a salt-rich core with a diameter on the order of 10 microns. The computational model for diffusive-convective transport in protein crystallization (see the First Report) has been applied to a realistic growth cell geometry, taking into account the findings of the above repartitioning studies. These results show that some elements of a moving boundary problem must be incorporated into the model in order to obtain a more realistic description. Our experimental setup for light scattering investigations of aggregation and nucleation in protein solutions has been extensively tested. Scattering intensity measurements with a true Rayleigh scatterer produced systematically increased forward scattering, indicating problems with glare. These have been resolved. Preliminary measurements with supersaturated lysozyme solutions revealed that the scatterers grow with time. Work has begun on a computer program for the unified evaluation of simultaneously obtained, multi-angle static and dynamic light scattering data.

Rosenberger, Franz↗
NASA NTRS1996

Interaction of a shock with a longitudinal vortex

In this paper we study the shock/longitudinal vortex interaction problem in axisymmetric geometry. Linearized analysis for small vortex strength is performed, and compared with results from a high order axisymmetric shock-fitted Euler solution obtained for this purpose. It is confirmed that for weak vortices, predictions from linear theory agree well with results from nonlinear numerical simulations at the shock location. To handle very strong longitudinal vortices, which may ultimately break the shock, we use an axisymmetric high order essentially non-oscillatory (ENO) shock capturing scheme. Comparison of shock-captured and shock-fitted results are performed in their regions of common validity. We also study the vortex breakdown as a function of Mach number ranging from 1.3 to 10, thus extending the range of existing results. For vortex strengths above a critical value. a triple point forms on the shock and a secondary shock forms to provide the necessary deceleration so that the fluid velocity can adjust to downstream conditions at the shock.

Erlebacher, Gordon↗