Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical stability”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Improved 3-D turbomachinery CFD algorithm

The building blocks of a computer algorithm developed for the time-accurate flow analysis of rotating machines are described. The flow model is a finite volume method utilizing a high resolution approximate Riemann solver for interface flux definitions. This block LU implicit numerical scheme possesses apparent unconditional stability. Multi-block composite gridding is used to orderly partition the field into a specified arrangement. Block interfaces, including dynamic interfaces, are treated such as to mimic interior block communication. Special attention is given to the reduction of in-core memory requirements by placing the burden on secondary storage media. Broad applicability is implied, although the results presented are restricted to that of an even blade count configuration. Several other configurations are presently under investigation, the results of which will appear in subsequent publications.

Janus, J. Mark↗
NASA NTRS1988

Transonic aeroelastic analysis of launch vehicle configurations

A numerical study of the aeroelastic stability of typical launch vehicle configurations in transonic flight is performed. Recent computational fluid dynamics techniques are used to simulate the transonic aerodynamic flow fields, as opposed to relying on experimental data for the unsteady aerodynamic pressures. The flow solver is coupled to an appropriate structural representation of the vehicle. The aerodynamic formulation is based on the thin layer approximation to the Reynolds-Averaged Navier-Stokes equations, where the account for turbulent mixing is done by the two-layer Baldwin and Lomax algebraic eddy viscosity model. The structural-dynamic equations are developed considering free-free flexural vibration of an elongated beam with variable properties and are cast in modal form. Aeroelastic analyses are performed by integrating simultaneously in the two sets of equations. By tracing the growth or decay of a perturbed oscillation, the aeroelastic stability of a given constant configuration can be ascertained. The method is described in detail, and results that indicate its application are presented. Applications include some validation cases for the algorithm developed, as well as the study of configurations known to have presented flutter programs in the past.

Filgueirasdeazevedo, Joao Luiz↗
NASA NTRS1989

Application of stochastic robustness to aircraft control systems

Stochastic robustness, a simple numerical procedure for estimating the stability robustness of linear, time-invariant systems, is applied to a forward-swept-wing aircraft control system. Based on Monte Carlo evaluation of the system's closed-loop eignevalues, this analysis approach introduces the probability of instability as a scalar stability robustness measure. The related stochastic root locus provides insight into robustness characteristics of the closed-loop system. Three Linear Quadratic controllers of decreasing robustness are chosen to demonstrate the use of stochastic robustness to analyze and compare control designs. Examples are presented illustrating the use of stochastic robustness analysis to address the effects of actuator dynamics and unmodeled dynamics on the stability robustness of the forward-swept-wing aircraft.

Stengel, Robert F.↗
NASA NTRS1990

Numerical investigation of laminar-turbulent transition in a flat plate wake

Lamina-turbulent transition of high-deficit flat plate wakes is investigated by direct numerical simulations using the complete Navier Stokes equations. The simulations are based on a spatial model so that both the base flow and the disturbance flow can develop in the downstream direction. The Navier Stokes equations are used in a vorticity-velocity form and are solved using a combination of finite difference and spectral approximations. Fourier series are used in the spanwise direction. Second-order finite-differences are used to approximate the spatial derivatives in the streamwise and transverse directions. For the temporal discretion, a combination of ADI, Crank-Nicolson, and Adams-Bashforth methods is employed. The discretized velocity equations are solved using fast Helmholtz solvers. Code validation is accomplished by comparison of the numerical results to both linear stability and to experiments. Calculations of two- and/or three-dimensional sinuous and mode disturbances in the wake of flat plate are undertaken. For calculations of two-dimensional disturbances, the wake is forced at an amplitude level so that nonlinear disturbance development may be observed. In addition, the forcing amplitude is varied in order to determine its effect on the disturbance behavior. To investigate the onset of three-dimensionality, the wake is forced with a small-amplitude three-dimensional disturbance and a larger amplitude two-dimensional disturbance. The two-dimensional forcing amplitude is varied in order to determine its influence on the three-dimensional flow field.

Fasel, Hermann F.↗
NASA NTRS1990

The inviscid stability of supersonic flow past heated or cooled axisymmetric bodies

The inviscid, linear, nonaxisymmetric, temporal stability of the boundary layer associated with the supersonic flow past axisymmetric bodies (with particular emphasis on long thin, straight circular cylinders), subject to heated or cooled wall conditions is investigated. The eigenvalue problem is computed in some detail for a particular Mach number or 3.8, revealing that the effect of curvature and the choice of wall conditions both have a significant effect on the stability of the flow. Both the asymptotic, large azimuthal wavenumber solution and the asymptotic, far downstream solution are obtained for the stability analysis and compared with numerical results. Additionally, asymptotic analyses valid for large radii of curvature with cooled/heated wall conditions, are presented. In general, important differences were found to exist between the wall temperature conditions imposed and the adiabatic wall conditions considered previously.

Shaw, Stephen J.↗
NASA NTRS1990

Instability during directional solidification - Gravitational effects

After an elementary introduction to solidification, recent research that attempts to understand the interaction of fluid flow with a crystal-melt interface under gravitational influence using linear stability analysis is reviewed. Recent numerical calculations of fluid flow and interface morphology during alloy solidification are briefly discussed.

Coriell, S. R.↗
NASA NTRS1992

The inviscid stability of supersonic flow past heated or cooled axisymmetric bodies

The inviscid, linear, nonaxisymmetric, temporal stability of the boundary layer associated with the supersonic flow past axisymmetric bodies (with particular emphasis on long thin, straight circular cylinders), subject to heated or cooled wall conditions is investigated. The eigenvalue problem is computed in some detail for a particular Mach number or 3.8, revealing that the effect of curvature and the choice of wall conditions both have a significant effect on the stability of the flow. Both the asymptotic, large azimuthal wavenumber solution and the asymptotic, far downstream solution are obtained for the stability analysis and compared with numerical results. Additionally, asymptotic analyses valid for large radii of curvature with cooled/heated wall conditions are presented. In general, important differences were found to exist between the wall temperature conditions imposed and the adiabatic wall conditions considered previously.

Shaw, Stephen J.↗
NASA NTRS1995

Thermally Stable, Piezoelectric and Pyroelectric Polymeric Substrates and Method Relating Thereto

Production of an electric voltage in response to mechanical excitation (piezoelectricity) or thermal excitation (pyroelectricity) requires a material to have a preferred dipole orientation in its structure. This preferred orientation or polarization occurs naturally in some crystals such as quartz and can be induced into some ceramic and polymeric materials by application of strong electric or mechanical fields. For some materials, a combination of mechanical and electrical orientation is necessary to completely polarize the material. The only commercially available piezoelectric polymer is poly(vinylidene fluoride) (PVF2). However, this polymer has material and process limitations which prohibit its use in numerous device applications where thermal stability is a requirement. By the present invention, thermally stable, piezoelectric and pyroelectric polymeric substrates were prepared from polymers having a softening temperature greater than 1000C. A metal electrode material is deposited onto the polymer substrate and several electrical leads are attached to it. The polymer substrate is heated in a low dielectric medium to enhance molecular mobility of the polymer chains. A voltage is then applied to the polymer substrate inducing polarization. The voltage is then maintained while the polymer substrate is cooled 'freezing in' the molecular orientation. The novelty of the invention resides in the process of preparing the piezoelectric and pyroelectric polymeric substrate. The nonobviousness of the invention is found in heating the polymeric substrate in a low dielectric medium while applying a voltage.

Simpson, Joycelyn O.↗
NASA NTRS2001

Panel Flutter and Sonic Fatigue Analysis for RLV

A methodology is presented for the flutter analysis of the seal of thermal protection system (TPS) panel of X-33 Advanced Technology Demonstrator test vehicle. The seal is simulated as a two-dimensional cantilevered panel with an elastic stopper, which is modeled as an equivalent spring. This cantilever beam-spring model under the aerodynamic pressure at supersonic speeds turns out to be an impact nonlinear dynamic system. The flutter analysis of the seal is thus carried out using, time domain numerical simulation with a displacement stability criterion. The flutter boundary of the seal is further verified with a family of three traditional and one nontraditional panel flutter models. The frequency domain method that applies eigenanalysis on the traditional panel flutter problem was used. The results showed that the critical dynamic pressure could be more than doubled with properly chosen material for the base stopper. The proposed methodology can be easily extended to three-dimensional panel seals with flow angularity.

Mei, Chuh↗
NASA NTRS2004

On the Total Variation of High-Order Semi-Discrete Central Schemes for Conservation Laws

We discuss a new fifth-order, semi-discrete, central-upwind scheme for solving one-dimensional systems of conservation laws. This scheme combines a fifth-order WENO reconstruction, a semi-discrete central-upwind numerical flux, and a strong stability preserving Runge-Kutta method. We test our method with various examples, and give particular attention to the evolution of the total variation of the approximations.

Bryson, Steve↗
NASA NTRS2010

Detailed Comparison of DNS to PSE for Oblique Breakdown at Mach 3

A pair of oblique waves at low amplitudes is introduced in a supersonic flat-plate boundary layer. Their downstream development and the concomitant process of laminar to turbulent transition is then investigated numerically using Direct Numerical Simulations (DNS) and Parabolized Stability Equations (PSE). This abstract is the last part of an extensive study of the complete transition process initiated by oblique breakdown at Mach 3. In contrast to the previous simulations, the symmetry condition in the spanwise direction is removed for the simulation presented in this abstract. By removing the symmetry condition, we are able to confirm that the flow is indeed symmetric over the entire computational domain. Asymmetric modes grow in the streamwise direction but reach only small amplitude values at the outflow. Furthermore, this abstract discusses new time-averaged data from our previous simulation CASE 3 and compares PSE data obtained from NASA's LASTRAC code to DNS results.

Mayer, Christian S. J.↗
NASA NTRS2012

Transition Within Leeward Plane of Axisymmetric Bodies at Incidence in Supersonic Flow

Boundary layer transition along the leeward symmetry plane of axisymmetric bodies at nonzero angle of attack in supersonic flow was investigated experimentally and numerically as part of joint research between the Japan Aerospace Exploration Agency (JAXA) and National Aeronautics and Space Administration (NASA). Transition over four axisymmetric bodies (namely, Sears-Haack body, semi-Sears-Haack body, straight cone and flared cone) with different axial pressure gradients was measured in two different facilities with different unit Reynolds numbers. The semi-Sears-Haack body and flared cone were designed at JAXA to broaden the range of axial pressure distributions. For a body shape with an adverse pressure gradient (i.e., flared cone), the experimentally measured transition patterns show an earlier transition location along the leeward symmetry plane in comparison with the neighboring azimuthal locations. For nearly zero pressure gradient (i.e.,straight cone), this feature is only observed at the larger unit Reynolds number. Later transition along the leeward plane was observed for the remaining two body shapes with a favorable pressure gradient. The observed transition patterns are only partially consistent with the numerical predictions based on linear stability analysis. Additional measurements are used in conjunction with the stability computations to explore the phenomenon of leeward line transition and the underlying transition mechanism in further detail.

Tokugawa, Naoko↗
NASA NTRS2015

Observations and Numerical Modeling of the Jovian Ribbon

Multiple wavelength observations made by the Hubble Space Telescope in early 2007 show the presence of a wavy, high-contrast feature in Jupiter's atmosphere near 30 degrees North. The "Jovian Ribbon," best seen at 410 nanometers, irregularly undulates in latitude and is time-variable in appearance. A meridional intensity gradient algorithm was applied to the observations to track the Ribbon's contour. Spectral analysis of the contour revealed that the Ribbon's structure is a combination of several wavenumbers ranging from k equals 8-40. The Ribbon is a dynamic structure that has been observed to have spectral power for dominant wavenumbers which vary over a time period of one month. The presence of the Ribbon correlates with periods when the velocity of the westward jet at the same location is highest. We conducted numerical simulations to investigate the stability of westward jets of varying speed, vertical shear, and background static stability to different perturbations. A Ribbon-like morphology was best reproduced with a 35 per millisecond westward jet that decreases in amplitude for pressures greater than 700 hectopascals and a background static stability of N equals 0.005 per second perturbed by heat pulses constrained to latitudes south of 30 degrees North. Additionally, the simulated feature had wavenumbers that qualitatively matched observations and evolved throughout the simulation reproducing the Jovian Ribbon's dynamic structure.

Saturnian Ribbon↗
NASA NTRSDate not supplied

Laminar to Turbulence Transition in Boundary Layers due to Tripping Devices

Direct numerical simulations (DNS)and linear stability analysis were performed to identify the mechanisms by which two-and three-dimensional trips promote early transition in subsonic boundary layers over a flat plate. Investigations were conducted for two boundary layer transition trip configurations: one is a two-dimensional rectangular-shaped rod and the other is a three-dimensional zigzag tape. Reynolds numbers based on the trip heights are 683 and 300 for the two-and three-dimensional trips, respectively. Simulations with the two-dimensional trip showed that the flow separates at the top of the trip and forms a long recirculation zone downstream. However, the flow remains steady in the simulation without any external disturbances. Linear stability analysis showed that the inflectional profile in the recirculation zone is strongly unstable to frequencies in the range of 100-300 Hz. Simulations with the external disturbances forced at the frequency of 200 Hz showed that disturbances grow strongly inside the recirculation zone as observed in the Klebanoff1 experiment.

Transition↗
NASA NTRSDate not supplied

Laminar to Turbulence Transition in Boundary Layers due to Tripping Devices

Direct numerical simulations (DNS)and linear stability analyses were performed to identify the mechanisms by whichtwo-and three-dimensional trips promote early transition in subsonic boundary layers over a flat plate. Investigations were conducted for threeboundary-layer transition trip configurations:a two-dimensional smooth hump,a two-dimensional rectangular-shapedrod and a three-dimensional zig-zagtape. The smooth hump was used by Park et al.1in the linear and nonlinear parabolized stability (PSE) computations. Reynolds numbersbased on the tripheights are 710, 683,and 300 for the two-and three-dimensional trips, respectively. Simulations withthe smooth hump showed that the flow becomes stable upstream and strongly unstable downstream of the hump. Simulations with the two-dimensional rectangular trip showed that the flow separates at the top of the trip and forms a long recirculation zonedownstream. However, the flow remains steady in the simulation without any external disturbances. Linear stability analysis showed that the inflectional profile in the recirculation zone is strongly unstable to frequencies in the range of 100-300 Hz. Simulations with the external disturbances forced at the frequency of 200 Hz showed that disturbances grow strongly inside the recirculation zone as observed in the Klebanoff2 experiment. Three-dimensional simulations with secondary disturbances of fundamental, subharmonic, and oblique types produced turbulent flows shortlydownstream. The three-dimensional zig-zag trip with a small roughness Reynolds number lead to early breakdown to a turbulent flow.

Stability↗
NASA NTRS2020

Entropy Stable Nonconforming Discretizations with the Summation-By-Parts Property for Curvilinear Coordinates

The entropy conservative/stable algorithm of Friedrichet al.(2018) for hyperbolic conservation laws on nonconforming p-refined/coarsened Cartesian grids, is extended to curvilinear grids for the compressible Euler equations. The primary focus is on constructing appropriate coupling procedures across the curvilinear nonconforming interfaces. A simple and flexible approach is proposed that uses interpolation operators from one element to the other. On the element faces,the analytic metrics are used to construct coupling terms, while metric terms in the volume are approximated to satisfy a discretization of the geometric conservation laws. The resulting scheme is entropy conservative/stable, elementwise conservative, and freestream preserving. The accuracy and stability properties of the resulting numerical algorithm are shown to be comparable to those ofthe original conforming scheme (∼p+ 1 convergence) in the context of the isentropic Euler vortex and the inviscid Taylor–Green vortex problems on manufactured high order grids.

David Del Rey Fernandez↗
NASA NTRS1999

Practical Aspects of Stabilized FEM Discretizations of Nonlinear Conservation Law Systems with Convex Extension

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the Galerkin least-squares (GLS) and the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the POE system. Central to the development of the simplified GLS and DG methods is the Degenerative Scaling Theorem which characterizes right symmetrizes of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobean matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler, Navier-Stokes, and magnetohydrodynamic (MHD) equations. Linear and nonlinear energy stability is proven for the simplified GLS and DG methods. Spatial convergence properties of the simplified GLS and DO methods are numerical evaluated via the computation of Ringleb flow on a sequence of successively refined triangulations. Finally, we consider a posteriori error estimates for the GLS and DG demoralization assuming error functionals related to the integrated lift and drag of a body. Sample calculations in 20 are shown to validate the theory and implementation.

Barth, Timothy↗