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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Transition Analysis for Isolated Trips on BOLT-II Wind-Tunnel and Flight Configuration

The BOLT-II configuration designed for an upcoming flight experiment exemplifies the complexities of hypersonic boundary layer transition in the presence of multiple and potentially interacting instability mechanisms. The present paper is related to the NASA-led roughness experiment on the secondary side of the BOLT-II configuration. First, computations are performed to aid mean-flow based correlations for the transition measurements acquired during an experimental campaign in the Ludwieg Tube facility at the United States Air Force Academy. To that end, both laminar and turbulent flow computations are performed for selected conditions of the experiment with a 1:3 scale model of the slightly shorter BOLT configuration. The variation in well-known correlation parameters for transition due to isolated roughness elements with respect to element position, height, and flow Reynolds number is delineated. Furthermore, the wake instability theory is used to examine the amplification of unstable disturbances behind the centerline trip at the design conditions of the BOLT-II flight experiment and for its dynamically scaled version on the 1:3 scale model of the wind tunnel experiment. To our knowledge, this study represents the first analysis of the wake instabilities in the presence of a prominent entropy layer outside the boundary layer. At the nominal flight design condition of Re ∞ ) = 5.44 and Re ∞ = 2.5 x 10 6 /m, a diamond planform (“pizza-box”) trip with k/δ = 0.70 and planform-halfwidth-to-height ratio of b/k = 3.0 at X/L = 0.254 induces a strong distortion of the boundary-layer flow in the vicinity of the symmetry plane, such that the wake structure includes a centerline streak produced by the flow retardation behind the trip along with multiple additional, prominent vorticity structures representing the necklace vortex system. Stability results reveal that each of these streaks can support a strong amplification of wake instabilities within a frequency band of 95 kHz–255 kHz, with a peak N-factor of greater than 15 beyond X/L ≈ 0.42. The dynamically scaled wind tunnel configuration with Re ∞ = 20.25 x 10 6 /m is also found to support N > 15 at similar axial stations.

BOLT↗
NASA NTRS1989

Modeling and control of flexible space stations (slew maneuvers)

Large orbiting space structures are expected to experience mechanical vibrations arising from several disturbing forces such as those induced by shuttle takeoff or docking and crew movements. The problem is considered of modeling and control of large space structures subject to these and other disturbing forces. The system consists of a (rigid) massive body, which may play the role of experimental modules located at the center of the space station and flexible configurations, consisting of several beams, forming the space structure. A complete dynamic model of the system was developed using Hamilton's principle. This model consists of radial equations describing the translational motion of the central body, rotational equations describing the attitude motions of the body and several beam equations governing the vibration of the flexible members (platform) including appropriate boundary conditions. In summary, the dynamics of the space structure is governed by a complex system of interconnected partial and ordinary differential equations. Using Lyapunov's approach the asymptotic stability of the space structure is investigated. For asymptotic stability of the rest state (nominal trajectory), feedback controls are suggested. In the investigation, stability of the slewing maneuvers is also considered. Several numerical results are presented for illustration of the impact of coupling and the effectiveness of the stabilizing controls. Some insight is provided into the complexity of modeling, analysis and stabilization of actual space structures.

Ahmed, N. U.↗
NASA NTRS2007

Effects of Nose Bluntness on Stability of Hypersonic Boundary Layers over Blunt Cone

Receptivity and stability of hypersonic boundary layers are numerically investigated for boundary layer flows over a 5-degree straight cone at a free-stream Mach number of 6.0. To compute the shock and the interaction of shock with the instability waves, we solve the Navier-Stokes equations in axisymmetric coordinates. The governing equations are solved 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. After the mean flow field is computed, disturbances are introduced at the upstream end of the computational domain. Generation of instability waves from leading edge region and receptivity of boundary layer to slow acoustic waves are investigated. Computations are performed for a cone with nose radii of 0.001, 0.05 and 0.10 inches that give Reynolds numbers based on the nose radii ranging from 650 to 130,000. The linear stability results showed that the bluntness has a strong stabilizing effect on the stability of axisymmetric boundary layers. The transition Reynolds number for a cone with the nose Reynolds number of 65,000 is increased by a factor of 1.82 compared to that for a sharp cone. The receptivity coefficient for a sharp cone is about 4.23 and it is very small, approx.10(exp -3), for large bluntness.

Kara, K.↗
NASA NTRS1993

An efficient procedure for cascade aeroelastic stability determination using nonlinear, time-marching aerodynamic solvers

A numerical eigenvalue problem formulation and a practical calculation procedure for exact eigenvalues and corresponding eigenvectors are developed and applied to a nonlinear, two-dimensional, time-marching full potential solver for cascade aeroelastic stability analysis. This procedure is based on the Lanczos recursive method and it directly calculates stability information about a nonlinear steady state. It is compared to conventional approaches in the frequency and time domains developed earlier and is found to be 100-10.000 times more computationally efficient. Eigenvalue constellations and the flutter results for flow through a cascade SR5 propfan airfoil are presented.

Mahajan, Aparajit J.↗
NASA NTRS1987

Stability and transition in supersonic boundary layers

The full three-dimensional time-dependent compressible Navier-Stokes equations are numerically solved by a Fourier-Chebyshev collocation algorithm to study the stability of supersonic flows over a flat plate. Several non-linear numerical experiments suggest the existence of a secondary instability which might provide a possible route to transition. The interaction of the modes involved in this secondary instability is possibly amenable to a Floquet theory. Pertinent differences between this instability and the more common incompressible K-type instabilities are pointed out.

Erlebacher, Gordon↗
NASA NTRS1994

The stabilizing effect of compressibility in turbulent shear flow

Direct numerical simulation of turbulent homogeneous shear flow is performed in order to clarify compressibility effects on the turbulence growth in the flow. The two Mach numbers relevant to homogeneous shear flow are the turbulent Mach number M(t) and the gradient Mach number M(g). Two series of simulations are performed where the initial values of M(g) and M(t) are increased separately. The growth rate of turbulent kinetic energy is observed to decrease in both series of simulations. This 'stabilizing' effect of compressibility on the turbulent energy growth rate is observed to be substantially larger in the DNS series where the initial value of M(g) is changed. A systematic companion of the different DNS cues shows that the compressibility effect of reduced turbulent energy growth rate is primarily due to the reduced level of turbulence production and not due to explicit dilatational effects. The reduced turbulence production is not a mean density effect since the mean density remains constant in compressible homogeneous shear flow. The stabilizing effect of compressibility on the turbulence growth is observed to increase with the gradient Mach number M(g) in the homogeneous shear flow DNS. Estimates of M(g) for the mixing and the boundary layer are obtained. These estimates show that the parameter M(g) becomes much larger in the high-speed mixing layer relative to the high-speed boundary layer even though the mean flow Mach numbers are the same in the two flows. Therefore, the inhibition of turbulent energy production and consequent 'stabilizing' effect of compressibility on the turbulence (over and above that due to the mean density variation) is expected to be larger in the mixing layer relative to the boundary layer in agreement with experimental observations.

Sarkar, S.↗
NASA NTRS1994

Coupled Marangoni-Benard/Rayleigh-Benard Instability with Temperature Dependent Viscosity

The onset of convection induced by coupled surface tension gradient and buoyancy forces is investigated with temperature dependent viscosity. Both surface tension and viscosity are assumed to vary linearly with temperature. The limiting case, Ma=0, is the buoyancy driven convection problem typically referred to as the Rayleigh-Benard stability problem. The other limiting case, Ra=0, is the surface tension (gradient) driven flow problem referred to as the Marangoni-Benard problem. The equations and boundary conditions obtained from the linear analysis are solved numerically as a generalized eigenvalue problem. Neutral stability curves for different viscosity slopes have been generated for the Marangoni-Benard and Rayleigh-Benard problems. It is shown that the curves can be collapsed to a single curve by appropriately scaling the results for each of the limiting cases. The critical Marangoni number is determined as a function of the slope of the viscosity temperature variation, epsilon, for different values of the Rayleigh number. When the viscosity decreases linearly with temperature, the coupled buoyancy-surface tension problem, including the limiting cases of Ra=0 and Ma=0, is found to be more stable than the constant viscosity case.

Skarda, J. Raymond Lee↗
NASA NTRS1971

Nutational stability of a dual-spin satellite under the influence of applied reaction torques

It is assumed that since the solar paddle attachments to the hub of the spacecraft are not rigidly locked, the effect of the solar panels can be replaced by a constant reaction torque acting on the hub of the spacecraft. This could result in the satellite having an equilibrium motion about an axis displaced from the nominal axis of symmetry. The variational equations of motion are developed about such an equilibrium position using the SAS-A spacecraft as a model. Energy dissipation on the rotor as well as the main body is included. This nonautonomous set of differential equations are linearized and transformed to an autonomous set using the Liapunov reducibility theorem. The stability of the kinematically similar system is examined numerically using representative SAS-A parameters for the case when either pair of solar panels is assumed to be loosely attached. Stability is verified for small system nutation angles (0.1 degree).

Bainum, P. M.↗
NASA NTRS1972

Nutational stability of a dual-spin satellite under the influence of applied reaction torques.

The nutational stability of a dual-spin satellite under the influence of an applied reaction torque is analyzed. It is assumed that since the solar paddle attachments to the hub of the spacecraft are not rigidly locked, the effect of the solar panels can be replaced by a constant reaction torque acting on the hub of the spacecraft. This could result in the satellite having an equilibrium motion about an axis displaced from the nominal axis of symmetry. The variational equations of motion are developed about such an equilibrium position using the SAS-A spacecraft as a model. Energy dissipation on the rotor as well as the main body is included. This nonautonomous set of differential equations are linearized and transformed to an autonomous set using the Lyapunov Reducibility Theorem. The stability of the kinematically similar system is examined numerically using representative SAS-A parameters for the case when either pair of solar panels is assumed to be loosely attached. Stability is verified for small system nutation angles (0.1 degree) although the time constant associated with the least damped mode is approximately one order of magnitude larger than that for the SAS-A system without the influence of the reaction torques.

Bainum, P. M.↗
NASA NTRS1994

A Numerical Comparison of Barrier and Modified Barrier Methods for Large-Scale Bound-Constrained Optimization

When a classical barrier method is applied to the solution of a nonlinear programming problem with inequality constraints, the Hessian matrix of the barrier function becomes increasingly ill-conditioned as the solution is approached. As a result, it may be desirable to consider alternative numerical algorithms. We compare the performance of two methods motivated by barrier functions. The first is a stabilized form of the classical barrier method, where a numerically stable approximation to the Newton direction is used when the barrier parameter is small. The second is a modified barrier method where a barrier function is applied to a shifted form of the problem, and the resulting barrier terms are scaled by estimates of the optimal Lagrange multipliers. The condition number of the Hessian matrix of the resulting modified barrier function remains bounded as the solution to the constrained optimization problem is approached. Both of these techniques can be used in the context of a truncated-Newton method, and hence can be applied to large problems, as well as on parallel computers. In this paper, both techniques are applied to problems with bound constraints and we compare their practical behavior.

Nash, Stephen G.↗
NASA NTRS2021

Wake Instability Behind Isolated Trip Near the Leading Edge of the BOLT-II Configuration

The BOLT-II (Holden Mission) configuration is an extended version of the BOLT flight article and depicts hypersonic boundary-layer transition in the presence of multiple and potentially interacting instability mechanisms. Several numerical studies of the boundary-layer instabilities over these configurations have been reported in the recent literature, including our previous studies of the modal instability characteristics of boundary-layer streaks adjacent to the minor-axis symmetry plane of the BOLT configuration and the wake instabilities behind a diamond planform (“pizza-box”) trip along the symmetry plane on the secondary side of the BOLT-II configuration. The present work extends the latter study to a scaled version of the same trip that is located in the region of nonzero crossflow in the vicinity of the leading edge at X/L = 0.5. Collectively, the trips along the symmetry plane and near the leading edge are the focus of the NASA roughness experiment on the secondary side of the BOLT-II configuration. The laminar basic state computation at the nominal flight design condition of Re ∞ = 5.44 x 10 6 /m and Re ∞ = 2.5 x 10 6 /m shows that the leading-edge trip with k/δ ≈ 0.70 and planform-side-length-to-height ratio of b/k = 3.0 induces multiple asymmetric, longitudinal streaks within the trip wake. The most prominent streak among these resembles a finite amplitude crossflow vortex, and it supports the amplification of multiple families of unstable modes. The application of multi-dimensional instability analysis to the wake flow reveals that the most amplified unstable mode can achieve a peak N-factor of up to 20 by X/L = 0.80, indicating that the onset of transition is more than likely to occur within eighty percent of the model length. To the best of our knowledge, the present study represents the first analysis including nonparallel and curvature effects on hypersonic tripwake instabilities in the presence of boundary-layer crossflow over a three-dimensional configuration.

Boundary layer transition↗
NASA NTRS1980

On the dynamics of large orbiting flexible beams and platforms oriented along the local horizontal

The dynamics and stability of large orbiting flexible beams and platforms oriented along the local horizontal are treated both analytically and numerically. It is assumed that such structures could be gravitationally stabilized by attaching a rigid lightweight dumbbell at the center of mass by a spring loaded hinge which also could provide viscous damping. For the beam it is seen that the small amplitude inplane pitch motion, dumbbell librational motion, and the anti-symmetric elastic modes are all coupled. The three dimensional equations of motion for a circular flat plate in orbit with a two-degree-of-freedom gimballed dumbbell are also developed and show that only those elastic modes described by a single nodal diameter line are influenced by the dumbbell motion. Stability criteria are developed for both examples and a parametric study of the least damped mode characteristics together with numerically simulated transient responses are carried out.

Bainum, P. M.↗
NASA NTRS1978

Investigation of ODE integrators using interactive graphics

Two FORTRAN programs using an interactive graphic terminal to generate accuracy and stability plots for given multistep ordinary differential equation (ODE) integrators are described. The first treats the fixed stepsize linear case with complex variable solutions, and generates plots to show accuracy and error response to step driving function of a numerical solution, as well as the linear stability region. The second generates an analog to the stability region for classes of non-linear ODE's as well as accuracy plots. Both systems can compute method coefficients from a simple specification of the method. Example plots are given.

Brown, R. L.↗
NASA NTRS1992

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems. Task 4: Advanced fan section aerodynamic analysis

The purpose of this study is the development of a three-dimensional Euler/Navier-Stokes flow analysis for fan section/engine geometries containing multiple blade rows and multiple spanwise flow splitters. An existing procedure developed by Dr. J. J. Adamczyk and associates and the NASA Lewis Research Center was modified to accept multiple spanwise splitter geometries and simulate engine core conditions. The procedure was also modified to allow coarse parallelization of the solution algorithm. This document is a final report outlining the development and techniques used in the procedure. The numerical solution is based upon a finite volume technique with a four stage Runge-Kutta time marching procedure. Numerical dissipation is used to gain solution stability but is reduced in viscous dominated flow regions. Local time stepping and implicit residual smoothing are used to increase the rate of convergence. Multiple blade row solutions are based upon the average-passage system of equations. The numerical solutions are performed on an H-type grid system, with meshes being generated by the system (TIGG3D) developed earlier under this contract. The grid generation scheme meets the average-passage requirement of maintaining a common axisymmetric mesh for each blade row grid. The analysis was run on several geometry configurations ranging from one to five blade rows and from one to four radial flow splitters. Pure internal flow solutions were obtained as well as solutions with flow about the cowl/nacelle and various engine core flow conditions. The efficiency of the solution procedure was shown to be the same as the original analysis.

Crook, Andrew J.↗
NASA NTRS1978

Investigation of ODE integrators using interactive graphics

In the present paper, two FORTRAN programs are described which use an interactive graphic terminal to generate accuracy and stability plots for a given multistep integrator of ordinary differential equations (ODE). One program treats the fixed-stepsize linear case with complex variable solutions, and generates plots to show accuracy and error response to step driving function of a numerical solution. It also generates the linear stability region. The other program generates an analog with an accuracy plot. Both systems are capable of computing method coefficients from a simple specification of the method. Example plots are given.

Brown, R. L.↗
NASA NTRS1978

Modification of the Douglas Neumann program to improve the efficiency of predicting component interference and high lift characteristics

The Douglas Neumann method for low-speed potential flow on arbitrary three-dimensional lifting bodies was modified by substituting the combined source and doublet surface paneling based on Green's identity for the original source panels. Numerical studies show improved accuracy and stability for thin lifting surfaces, permitting reduced panel number for high-lift devices and supercritical airfoil sections. The accuracy of flow in concave corners is improved. A method of airfoil section design for a given pressure distribution, based on Green's identity, was demonstrated. The program uses panels on the body surface with constant source strength and parabolic distribution of doublet strength, and a doublet sheet on the wake. The program is written for the CDC CYBER 175 computer. Results of calculations are presented for isolated bodies, wings, wing-body combinations, and internal flow.

Bristow, D. R.↗