SEARCH · Search NASA

Results for “nonlinear eigenvalues”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Application of the reduced basis method to nonlinear transient thermal analysis

An effort to apply the reduced basis method to nonlinear transient thermal analysis is described. The method combines the classical Rayleigh-Ritz and modal superposition techniques with contemporary finite element methods to retain modeling versatility as the degrees of freedom in a problem are reduced. The essence of the method is to use a few thermal modes from eigenvalue analyses as basis vectors to represent the temperature response for a given thermal problem similar to the use of vibration modes to represent displacements in a dynamic response problem. Approximate temperature distributions were obtained using the reduced basis method for a small section of the Shuttle Orbiter lower wing undergoing reentry heating. Good agreement was obtained between the reduced basis method solutions and full system solutions with reductions in the degrees of freedom of up to a factor of four. The good agreement indicates the reduced basis method has the potential for significant reduction in computing effort for thermal analysis; however, considerable work remains to determine techniques for selecting the type and number of basis vectors needed for approximate solutions to more complex transient thermal problems.

Shore, C. P.↗

NASA NTRS1989

Sensitivity of actively damped structures to imperfections and modeling errors

The sensitivity of actively damped response of structures with respect to errors in the structural modeling is studied. Two ways of representing errors are considered. The first approach assumes errors in the form of spatial variations (or imperfections) in the assumed mass and stiffness properties of the structures. The second approach assumes errors due to such factors as unknown joint stiffnesses, discretization errors, and nonlinearities. These errors are represented here as discrepancies between experimental and analytical mode shapes and frequencies. The actively damped system considered here is a direct-rate feedback regulator based on a number of colocated velocity sensors and force actuators. The response of the controlled structure is characterized by the eigenvalues of the closed-loop system. The effects of the modeling errors are thus presented as the sensitivity of the eigenvalues of the closed-loop system. Results are presented for two examples: (1) a three-span simply supported beam controlled by three sensors and actuators, and (2) a laboratory structure consisting of a cruciform beam supported by cables.

Haftka, Raphael T.↗

NASA NTRS1988

The dynamic behaviour of postbuckled composite plates under acoustic excitation

The Rayleigh-Ritz method was used to find the postbuckling static displacement pattern of a composite plane (CFRP) under uniaxial in-plane compression of uniform edge-shortening. The resonance frequencies and mode shapes at various postbuckled states are then evaluated by eigenvalue analysis of the dynamical matrix equation consisting of up-dated tangential stiffness matrix at corresponding static configuration. The theoretical results are compared with experimental results obtained in 16-layered CFRP laminate of aspect ratio 1.5. The resonance frequencies and mode shapes obtained are used to interpret the multimodal and nonlinear strain responses to high level of acoustic excitation. The dominance of second-mode contribution and softening-spring behavior are found in the strain response of postbuckled plates.

Ng, C. F.↗

NASA NTRS1978

Gaussian step-pressure loading of rigid viscoplastic plates

The response of a thin, rigid viscoplastic plate subjected to a spatially axisymmetric Gaussian step pressure impulse loading was studied analytically. A Gaussian pressure distribution in excess of the collapse load was applied to the plate, held constant for a length of time, and then suddenly removed. The plate deforms with monotonically increasing deflections until the dynamic energy is completely dissipated in plastic work. The simply supported plate of uniform thickness obeys the von Mises yield criterion and a generalized constitutive equation for rigid viscoplastic materials. For the small deflection bending response of the plate, the governing system of equations is essentially nonlinear. Transverse shear stress is neglected in the yield condition and rotary inertia in the equations of dynamic equilibrium. A proportional loading technique, known to give excellent approximations of the exact solution for the uniform load case, was used to linearize the problem and to obtain the analytical solutions in the form of eigenvalue expansions. The effects of load concentration, of an order of magnitude change in the viscosity of the plate material, and of load duration were examined while holding the total impulse constant.

Hayduk, R. J.↗

NASA NTRS1987

A combined stochastic feedforward and feedback control design methodology with application to autoland design

A combined stochastic feedforward and feedback control design methodology was developed. The objective of the feedforward control law is to track the commanded trajectory, whereas the feedback control law tries to maintain the plant state near the desired trajectory in the presence of disturbances and uncertainties about the plant. The feedforward control law design is formulated as a stochastic optimization problem and is embedded into the stochastic output feedback problem where the plant contains unstable and uncontrollable modes. An algorithm to compute the optimal feedforward is developed. In this approach, the use of error integral feedback, dynamic compensation, control rate command structures are an integral part of the methodology. An incremental implementation is recommended. Results on the eigenvalues of the implemented versus designed control laws are presented. The stochastic feedforward/feedback control methodology is used to design a digital automatic landing system for the ATOPS Research Vehicle, a Boeing 737-100 aircraft. The system control modes include localizer and glideslope capture and track, and flare to touchdown. Results of a detailed nonlinear simulation of the digital control laws, actuator systems, and aircraft aerodynamics are presented.

Halyo, Nesim↗

NASA NTRS2017

Stabilization of Hypersonic Boundary Layers by Linear and Nonlinear Optimal Perturbations

The effect of stationary, finite-amplitude, linear and nonlinear optimal perturbations on the modal disturbance growth in a Mach 6 axisymmetric flow over a 7 deg. half-angle cone with 0:126 mm nose radius and 0:305 m length is investigated. The freestream parameters (M = 6, Re(exp 1) = 18 x 10(exp. 6) /m) are selected to match the flow conditions of a previous experiment in the VKI H3 hypersonic tunnel. Plane-marching parabolized stability equations are used in conjunction with a partial-differential equation based planar eigenvalue analysis to characterize the boundary layer instability in the presence of azimuthally periodic streaks. The streaks are observed to stabilize nominally planar Mack mode instabilities, although oblique Mack mode and first-mode disturbances are destabilized. Experimentally measured transition onset in the absence of any streaks correlates with an amplification factor of N = 6 for the planar Mack modes. For high enough streak amplitudes, the transition threshold of N = 6 is not reached by the Mack mode instabilities within the length of the cone; however, subharmonic first-mode instabilities, which are destabilized by the presence of the streaks, do reach N = 6 near the end of the cone. The highest stabilization is observed at streak amplitudes of approximately 20 percent of the freestream velocity. Because the use of initial disturbance profiles based on linear optimal growth theory may yield suboptimal control in the context of nonlinear streaks, the computational predictions are extended to nonlinear optimal growth theory. Results show that by using nonlinearly optimal perturbation leads to slightly enhanced stabilization of plane Mack mode disturbances as well as reduced destabilization of subharmonic first-mode disturbances.

Incompressible boundary layer↗

NASA NTRS1994

Benchmark problems in computational aeroacoustics

A recent directive at NASA Langley is aimed at numerically predicting principal noise sources. During my summer stay, I worked with high-order ENO code, developed by Dr. Harold Atkins, for solving the unsteady compressible Navier-Stokes equations, as it applies to computational aeroacoustics (CAA). A CAA workshop, composed of six categories of benchmark problems, has been organized to test various numerical properties of code. My task was to determine the robustness of Atkins' code for these test problems. In one category, we tested the nonlinear wave propagation of the code for the one-dimensional Euler equations, with initial pressure, density, and velocity conditions. Using freestream boundary conditions, our results were plausible. In another category, we solved the linearized two-dimensional Euler equations to test the effectiveness of radiation boundary conditions. Here we utilized MAPLE to compute eigenvalues and eigenvectors of the Jacobian given variable and flux vectors. We experienced a minor problem with inflow and outflow boundary conditions. Next, we solved the quasi one dimensional unsteady flow equations with an incoming acoustic wave of amplitude 10(exp -6). The small amplitude sound wave was incident on a convergent-divergent nozzle. After finding a steady-state solution and then marching forward, our solution indicated that after 30 periods the acoustic wave had dissipated (a period is time required for sound wave to traverse one end of nozzle to other end).

Porter-Locklear, Freda↗

NASA NTRS1990

Efficient numerical simulation of electron states in quantum wires

A new algorithm is presented for the numerical simulation of electrons in a quantum wire as described by a two-dimensional eigenvalue problem for Schroedinger's equation coupled with Poisson's equation. Initially, the algorithm employs an underrelaxed fixed point iteration to generate an approximation which is reasonably close to the solution. Subsequently, this approximate solution is employed as an initial guess for a Jacobian-free implementation of an approximate Newton method. In this manner the nonlinearity in the model is dealt with effectively. The effectiveness of this approach is demonstrated in a set of numerical experiments which study the electron states on the cross section of a quantum wire structure based on III-V semiconductors at 4.2 and 77 K.

Kerkhoven, Thomas↗

NASA NTRS2001

Discontinuous Galerkin Methods for NonLinear Differential Systems

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 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 PDE (partial differential equation) system. Central to the development of the simplified DG methods is the Eigenvalue Scaling Theorem which characterizes right symmetrizers of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobian matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler equations of gas dynamics and extended conservation law systems derivable as moments of the Boltzmann equation. Using results from kinetic Boltzmann moment closure theory, we then derive and prove energy stability for several approximate DG fluxes which have practical and theoretical merit.

Barth, Timothy↗

NASA NTRS1990

On a class of unsteady three-dimensional Navier Stokes solutions relevant to rotating disc flows: Threshold amplitudes and finite time singularities

A class of exact steady and unsteady solutions of the Navier Stokes equations in cylindrical polar coordinates is given. The flows correspond to the motion induced by an infinite disc rotating with constant angular velocity about the z-axis in a fluid occupying a semi-infinite region which, at large distances from the disc, has velocity field proportional to (x,-y,O) with respect to a Cartesian coordinate system. It is shown that when the rate of rotation is large, Karman's exact solution for a disc rotating in an otherwise motionless fluid is recovered. In the limit of zero rotation rate a particular form of Howarth's exact solution for three-dimensional stagnation point flow is obtained. The unsteady form of the partial differential system describing this class of flow may be generalized to time-periodic equilibrium flows. In addition the unsteady equations are shown to describe a strongly nonlinear instability of Karman's rotating disc flow. It is shown that sufficiently large perturbations lead to a finite time breakdown of that flow whilst smaller disturbances decay to zero. If the stagnation point flow at infinity is sufficiently strong, the steady basic states become linearly unstable. In fact there is then a continuous spectrum of unstable eigenvalues of the stability equations but, if the initial value problem is considered, it is found that, at large values of time, the continuous spectrum leads to a velocity field growing exponentially in time with an amplitude decaying algebraically in time.

Hall, Philip↗

NASA NTRS1975

Sustained small oscillations in nonlinear control systems

Some results of bifurcation theory were used to study the existence of small-amplitude periodic behavior in launch vehicle dynamics, assuming that nonlinearity exists as a cubic term in the rudder response. The analysis follows closely Sattinger's (1973) approach to the theory of periodic bifurcations. The conditions under which a bifurcating branch of orbitally stable periodic solutions will exist are determined. It is shown that in more complicated cases, the conditions under which the system matrix has a pair of simple purely imaginary eigenvalues can be determined with the aid of linear stability techniques.

George, J. H.↗

NASA NTRS1985

Three-dimensional baroclinic instability of a Hadley cell for small Richardson number

A three-dimensional, linear stability analysis of a baroclinic flow for Richardson number, Ri, of order unity is presented. The model considered is a thin horizontal, rotating fluid layer which is subjected to horizontal and vertical temperature gradients. The basic state is a Hadley cell which is a solution of the complete set of governing, nonlinear equations and contains both Ekman and thermal boundary layers adjacent to the rigid boundaries; it is given in a closed form. The stability analysis is also based on the complete set of equations; and perturbation possessing zonal, meridional, and vertical structures were considered. Numerical methods were developed for the stability problem which results in a stiff, eighth-order, ordinary differential eigenvalue problem. The previous work on three-dimensional baroclinic instability for small Ri was extended to a more realistic model involving the Prandtl number, sigma, and the Ekman number, E, and to finite growth rates and a wider range of the zonal wavenumber.

Antar, B. N.↗

NASA NTRS1988

Recent developments in structural sensitivity analysis

Recent developments are reviewed in two major areas of structural sensitivity analysis: sensitivity of static and transient response; and sensitivity of vibration and buckling eigenproblems. Recent developments from the standpoint of computational cost, accuracy, and ease of implementation are presented. In the area of static response, current interest is focused on sensitivity to shape variation and sensitivity of nonlinear response. Two general approaches are used for computing sensitivities: differentiation of the continuum equations followed by discretization, and the reverse approach of discretization followed by differentiation. It is shown that the choice of methods has important accuracy and implementation implications. In the area of eigenproblem sensitivity, there is a great deal of interest and significant progress in sensitivity of problems with repeated eigenvalues. In addition to reviewing recent contributions in this area, the paper raises the issue of differentiability and continuity associated with the occurrence of repeated eigenvalues.

Haftka, Raphael T.↗

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 NTRS1984

The collisionless macroscopic Kelvin-Helmholtz instabiity. I Transverse electrostatic mode

The electrostatic limit of the Kelvin-Helmholtz instability in a collisionless plasma is investigated for a magnetic field perpendicular to the velocity flow. A two-fluid eigenvalue analysis of the transverse electrostatic Kelvin-Helmholtz instability is performed using the Chew-Goldberger-Low hydromagnetic equations with first-order finite Larmer radius corrections, establishing that the transverse electrostatic mode reduces to the incompressible transverse MHD mode in the limit of zero ion gyroradius. The linear growth rates are determined and a strong stabilization of the instability for (k-perpendicular) (rho sub i) of 0.2 or greater is observed. The nonlinear stage of the instability displays large vortices whose size is many times the initial velocity shear length. At saturation the simulation is dominated by the longest-wavelength mode permitted in the system. The simulation results can be understood in terms of a simple model which assumes that saturation occurs when the angular frequency of the vortex motion reaches the linear growth rate of the instability.

Prichett, P. L.↗

NASA NTRS1994

Nonlinear equilibration of localized instabilities on a baroclinic jet

Dynamical mechanisms underlying the equilibration of absolute instability are examined in a nonlinear, quasigeostrophic, two-layer model. The key to understanding the nonlinear equilibration is in recognizing that linear absolute instabilities can be stabilized both by a reduction of the vertical shear and by enhancement of the mean barotropic velocity. In a localized domain, the equilibration process proceeds with the creation of locally convectively unstable regions downstream, which encroach onto the locally absolutely unstable region until the local instability is suppressed. That local instabilities exist only if absolutely unstable regions span a minimum size is verified by eigenvalue calculations of three-dimensional flows. Numerical examples suggest that this critical size is at least 9000 km for a wide range of parameter values chosen to investigate the midlatitude storm tracks. Fluctuations arising from local absolute instability obtain maximum amplitude in the downstream convectively unstable regions rather than in the absolutely unstable regions themselves. Together, these results suggest that if an equilibrated absolute instability were to occur in midlatitudes, a zonal band of surface easterlies exceeding 900 km would be required and the associated enhanced variances would not be found coincident with the regions of absolute instability.

Delsole, Timothy M.↗

NASA NTRS1983

Accelerating an iterative process by explicit annihilation

A slowly convergent stationary iterative process can be accelerated by explicitly annihilating (i.e., eliminating) the dominant eigenvector component of the error. The dominant eigenvalue or complex pair of eigenvalues can be estimated from the solution during the iteration. The corresponding eigenvector or complex pair of eigenvectors can then be annihilated by applying an explicit Richardson process over the basic iterative method. This can be done entirely in real arithmetic by analytically combining the complex conjugate annihilation steps. The technique is applied to an implicit algorithm for the calculation of two dimensional steady transonic flow over a circular cylinder using the equations of compressible inviscid gas dynamics. This demonstrates the use of explicit annihilation on a nonlinear problem.

Jespersen, D. C.↗

NASA NTRS1985

Accelerating an iterative process by explicit annihilation

Jespersen, D. C.↗