SEARCH · Engineering Papers

Results for “nonlinear eigenvalue problem”

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 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Numerical simulation of spatially-evolving instability

A computational study of the spatial stability of plane Poiseuille flow is presented. The numerical scheme employs a time-splitting method to integrate the full Navier-Stokes equations using spectral collocation/finite-difference discretization on a non-staggered mesh. The eigenvalue decomposition procedure is applied for the solution of the Poisson equations using the capacitance matrix technique. The buffer domain method is incorporated for the outflow boundary conditions. The input perturbation velocities are obtained by solving the Orr-Sommerfeld equation for the nonlinear eigenvalue problem employing the companion matrix method. Computational results are compared with the linear theory for two-dimensional disturbances.

Danabasoglu, G.↗

NASA NTRS1992

Buckling of imperfect, anisotropic, ring-stiffened cylinders under combined loads

The objectives of this study is to develop an analysis to predict loads of ring-stiffened anisotropic cylinders subject to axial compression, torsion and internal pressure. This structure is modeled as a branched shell. A nonlinear axisymmetric prebuckling equilibrium state is assumed which is amenable to exact solution within each branch. Axisymmetric geometric imperfections are included. Buckling displacements are represented by a Fourier series in the circumferential coordinate and the finite element method in the axial coordinate. Application of Trefftz criterion to the second variation of the total potential energy leads to a nonlinear eigenvalue problem for the buckling load and mode. Results are presented for both unstiffened and ring-stiffened cylinders in the form of buckling interaction diagrams. Imperfections can cause an unexpected buckling mode in the ring web which would not occur for the perfect structure, and pressurization diminishes the benefit of adding rings to the unstiffened shell to increase the buckling load. The implementation of the analysis methodology into a structural sizing algorithm is discussed.

Ley, Robert P.↗

NASA NTRS1992

Comparison of level-1 and level-2 buckling and postbuckling solutions

The accuracy of the solutions obtained with the Level-1 computational modules of DISDECO (Delft Interactive Shell Design Code) are investigated by comparing them with solutions obtained with ANILISA, a Level-2 computational module. Whereas the Level-1 solutions are based on the use of trigonometric functions to represent the out-of-plane displacement W and involve Galerkin type approximations to reduce the governing partial differential equations to simple algebraic expressions, the Level-2 solutions employ a truncated Fourier expansion in the circumferential direction and solve the resulting ordinary differential equations forming a nonlinear eigenvalue problem by a Stodola like technique. This involves the numerical integration of the set of ordinary differential equations, which in turn makes it possible to satisfy the specified boundary conditions exactly. It is shown that in many applications the Level-1 predictions concerning the critical buckling load of the perfect structure and its imperfection sensitivity are quite accurate and reliable. For those cases where this is not true, it is shown how the pre and postprocessing facilities of DISDECO can be used efficiently to arrive at a reliable prediction. This involves always the use of the Level-2 computational module ANILISA with its rigorous prebuckling solution and its feature which makes it possible to satisfy the specified boundary conditions exactly.

Arbocz, J.↗

NASA NTRS1994

Buckling of imperfect, anisotropic, ring-stiffened cylinders under combined loads

The objective of this study is to develop an anlysis to predict buckling loads of ring-stiffened anisotropic cylinders subject to axial compression, torsion, and internal pressure. This structure is modeled as a branched shell. A nonlinear axisymmetric prebuckling equilibrium state is assumed which is amenable to an exact solution within each branch. Axisymmetric geometric imperfections are included. Buckling displacements are represented by a Fourier series in the circumferential coordination and the finite-element method in the axial coordinate. Application of the Trefftz criterion to the second variation of the total potential energy leads to a nonlinear eigenvalue problem for the buckling load and mode. Results are presented for both unstiffened and ring-stiffened cylinders in the form of buckling interaction diagrams. Imperfections can cause an unexpected buckling mode in the ring web which would not occur for the perfect structure, and pressurization diminishes the benefit of adding rings to the unstiffened shell to increase the buckling load. The implementation of the analysis methodology into a structural sizing algorithm is discussed.

Ley, Robert P.↗

DOE OSTI2019

Two-tier model reduction of viscoelastically damped finite element models

Transient simulations of linear viscoelastically damped structures necessitate critical computational resources to directly integrate the full-order finite element model with time-stepping algorithms. Traditional modal reduction techniques are not directly applicable to these systems since viscoelastic materials depend on time and frequency. A more appropriate reduction basis is obtained from the nonlinear, complex eigenvalue problem, whose eigenvectors capture the appropriate kinematics and enable frequency-based mode selection; unfortunately, the computational cost is prohibitive for computing these modes from large-scale engineering models. To address this shortcoming, this work proposes a novel two-tier reduction procedure to reduce the upfront cost of solving the complex, nonlinear eigenvalue problem. The first reduction step reduces the full-order model with real mode shapes linearized about various centering frequencies to capture the kinematics over a full range of viscoelastic material behavior (glassy, rubbery, and glass-transition zones). This tier-one reduction preserves time-temperature superposition and allows the equations to depend parametrically on operating temperature. The second-level reduction then solves the complex, nonlinear eigenmode solutions in the tier-one reduced space about a fixed temperature to further reduce the equations-of-motion. The method is demonstrated on a cantilevered sandwich plate to display its accuracy and efficiency in comparison to full-order model predictions.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1984

Differential eigenvalue problems in which the parameter appears nonlinearly

Several methods are examined for determining the eigenvalues of a system of equations in which the parameter appears nonlinearly. The equations are the result of the discretization of differential eigenvalue problems using a finite Chebyshev series. Two global methods are considered which determine the spectrum of eigenvalues without an initial estimate. A local iteration scheme with cubic convergence is presented. Calculations are performed for a model second order differential problem and the Orr-Sommerfeld problem for plane Poiseuille flow.

Bridges, T. J.↗

NASA NTRS1987

Boundary layer stability calculations

In this paper numerical calculation of the spatial stability of disturbances in the parallel and nonparallel Blasius boundary layers is considered. Chebyshev polynomials are used for discretization. The problem with the boundary condition at infinity is overcome, and the resulting nonlinear matrix eigenvalue problem is attacked directly. The secondary eigenvalue problem for three-dimensional disturbances is shown to be uniformly stable, and particular solutions of this problem generated by the Orr-Sommerfeld equation are shown. A numerical solution of the nonparallel problem is considered using Chebyshev polynomials. The matrix equations are analyzed directly and the problem of uniqueness of the nonparallel correction is settled by careful application of the Fredholm alternative. Nonparallel corrections to the streamwise eigenfunction are shown.

Bridges, Thomas J.↗

DOE OSTI2018

On nonlinear inexact Uzawa algorithms for stabilized saddle point problems

In this paper, the nonlinear inexact Uzawa (NIU) methods for saddle point problems are studied. A further improved convergence result of the NIU algorithm for stabilized saddle point problems is presented, which results in much smaller convergence factor. This result can also be viewed as a generalization of some previous work for classical saddle point problems. Based on this result, we generalize some adaptive versions of the NIU algorithm with variable iteration parameters for classical saddle point problems to solve stabilized saddle point problems. Convergence analyses of these algorithms are presented. It shows that they converge under similar conditions as those in classical cases, which are more practical and without estimates on the extreme eigenvalues of the involved preconditioned systems. Numerical experiments are given to demonstrate the effectiveness of these algorithms over some original NIU algorithms.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1981

Multigrid solvers on parallel computers

Massively parallel computers, as considered in this investigation, are not yet available. However, a large-scale parallel computer cannot usefully be designed before the hypothetical algorithms which will employ it are studied. Most of the studies of parallel partial differential equations (PDE) solvers are based on solution techniques much slower (on sequential machines) than multigrid methods. Multigrid methods are highly parallelizable. Each of their processes can simultaneously be performed at all grid points. The present investigation is concerned with a preliminary exploration of the potential of multigrid, or, more generally, Multi-Level Adaptive Techniques (MLAT) on computers with many processors. Basic processes are considered, taking into account coarse-grid approximation, relaxation, coarse-grid corrections, full multigrid algorithms, nonlinear problems and eigenvalue problems, fine-to-coarse correction, and chains of problems. Details of parallel multigrid processing are also examined.

Brandt, A.↗

NASA NTRSDate not supplied

Adjoint-based Sensitivities of Flutter Predictions based on the Linearized Frequency-domain Approach

Flutter is a critical factor in designing and certifying aircraft. The linearized frequency-domain method offers a lower cost alternative to time-marching computational fluid dynamics for high-fidelity flutter analysis. In this work, adjoint-based sensitivities are added to a flutter analysis based on the linearized frequency-domain method to efficiently compute derivatives of flutter cost functions with respect to design variables or uncertain parameters. The derivation of the adjoint equations, which involve complications such as derivatives of a nonlinear generalized eigenvalue problem with complex-valued inputs and derivatives of the linearized Navier-Stokes equations, is provided. The implemented adjoint terms and derivatives are verified before demonstrating the approach for derivatives of flutter dynamic pressure with respect to Mach number for the AGARD 445.6 wing.

Aeroelasticity↗

DOE OSTI2014

A comparison of acceleration methods for solving the neutron transport k-eigenvalue problem

Over the past several years a number of papers have been written describing modern techniques for numerically computing the dominant eigenvalue of the neutron transport criticality problem. These methods fall into two distinct categories. The first category of methods rewrite the multi-group k-eigenvalue problem as a nonlinear system of equations and solve the resulting system using either a Jacobian-Free Newton–Krylov (JFNK) method or Nonlinear Krylov Acceleration (NKA), a variant of Anderson Acceleration. These methods are generally successful in significantly reducing the number of transport sweeps required to compute the dominant eigenvalue. The second category of methods utilize Moment-Based Acceleration (or High-Order/Low-Order (HOLO) Acceleration). These methods solve a sequence of modified diffusion eigenvalue problems whose solutions converge to the solution of the original transport eigenvalue problem. This second class of methods is, in our experience, always superior to the first, as most of the computational work is eliminated by the acceleration from the LO diffusion system. In this paper, we review each of these methods. Our computational results support our claim that the choice of which nonlinear solver to use, JFNK or NKA, should be secondary. The primary computational savings result from the implementation of a HOLO algorithm. We display computational results for a series of challenging multi-dimensional test problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Applying nonlinear diffusion acceleration to the neutron transport k-Eigenvalue problem with anisotropic scattering

High-order/low-order (or moment-based acceleration) algorithms have been used to significantly accelerate the solution to the neutron transport k-eigenvalue problem over the past several years. Recently, the nonlinear diffusion acceleration algorithm has been extended to solve fixed-source problems with anisotropic scattering sources. In this paper, we demonstrate that we can extend this algorithm to k-eigenvalue problems in which the scattering source is anisotropic and a significant acceleration can be achieved. Lastly, we demonstrate that the low-order, diffusion-like eigenvalue problem can be solved efficiently using a technique known as nonlinear elimination.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2015

Convergence analysis of two-node CMFD method for two-group neutron diffusion eigenvalue problem

In this paper, the nonlinear coarse-mesh finite difference method with two-node local problem (CMFD2N) is proven to be unconditionally stable for neutron diffusion eigenvalue problems. The explicit current correction factor (CCF) is derived based on the two-node analytic nodal method (ANM2N), and a Fourier stability analysis is applied to the linearized algorithm. It is shown that the analytic convergence rate obtained by the Fourier analysis compares very well with the numerically measured convergence rate. It is also shown that the theoretical convergence rate is only governed by the converged second harmonic buckling and the mesh size. It is also noted that the convergence rate of the CCF of the CMFD2N algorithm is dependent on the mesh size, but not on the total problem size. This is contrary to expectation for eigenvalue problem. The novel points of this paper are the analytical derivation of the convergence rate of the CMFD2N algorithm for eigenvalue problem, and the convergence analysis based on the analytic derivations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1989

A hybrid-perturbation-Galerkin technique which combines multiple expansions

A two-step hybrid perturbation-Galerkin method for the solution of a variety of differential equations type problems is found to give better results when multiple perturbation expansions are employed. The method assumes that there is parameter in the problem formulation and that a perturbation method can be sued to construct one or more expansions in this perturbation coefficient functions multiplied by computed amplitudes. In step one, regular and/or singular perturbation methods are used to determine the perturbation coefficient functions. The results of step one are in the form of one or more expansions each expressed as a sum of perturbation coefficient functions multiplied by a priori known gauge functions. In step two the classical Bubnov-Galerkin method uses the perturbation coefficient functions computed in step one to determine a set of amplitudes which replace and improve upon the gauge functions. The hybrid method has the potential of overcoming some of the drawbacks of the perturbation and Galerkin methods as applied separately, while combining some of their better features. The proposed method is applied, with two perturbation expansions in each case, to a variety of model ordinary differential equations problems including: a family of linear two-boundary-value problems, a nonlinear two-point boundary-value problem, a quantum mechanical eigenvalue problem and a nonlinear free oscillation problem. The results obtained from the hybrid methods are compared with approximate solutions obtained by other methods, and the applicability of the hybrid method to broader problem areas is discussed.

Geer, James F.↗

NASA NTRS1990

A hybrid perturbation-Galerkin technique that combines multiple expansions

A two-step hybrid perturbation-Galerkin method for the solution of a variety of differential equations type problems is found to give better results when multiple perturbation expansions are employed. The method assumes that there is parameter in the problem formulation and that a perturbation method can be used to construct one or more expansions in this perturbation coefficient functions multiplied by computed amplitudes. In step one, regular and/or singular perturbation methods are used to determine the perturbation coefficient functions. The results of step one are in the form of one or more expansions each expressed as a sum of perturbation coefficient functions multiplied by a priori known gauge functions. In step two the classical Bubnov-Galerkin method uses the perturbation coefficient functions computed in step one to determine a set of amplitudes which replace and improve upon the gauge functions. The hybrid method has the potential of overcoming some of the drawbacks of the perturbation and Galerkin methods as applied separately, while combining some of their better features. The proposed method is applied, with two perturbation expansions in each case, to a variety of model ordinary differential equations problems including: a family of linear two-boundary-value problems, a nonlinear two-point boundary-value problem, a quantum mechanical eigenvalue problem and a nonlinear free oscillation problem. The results obtained from the hybrid methods are compared with approximate solutions obtained by other methods, and the applicability of the hybrid method to broader problem areas is discussed.

Geer, James F.↗

DOE OSTI2022

Multi-harmonic Rutherford island theory

Rutherford island theory, which governs the nonlinear evolution of tearing modes in tokamak plasmas, is generalized to take into account situations in which the conventional one-harmonic approximation is not valid. The analysis incorporates non-inductive currents driven by radio frequency (RF) electromagnetic waves injected into the plasma. A multi-harmonic tearing mode dispersion relation is derived that takes the form of a nonlinear inhomogeneous matrix eigenvalue problem. The dispersion relation is solved in the so-called two-harmonic approximation, in which only the principal Fourier harmonic of the perturbed magnetic flux and its first overtone are included in the calculation. In the absence of RF current drive, the nonlinear behavior of a tearing mode predicted in the two-harmonic approximation does not differ substantially from that predicted in the one-harmonic approximation. On the other hand, RF current drive that is sufficiently localized in the vicinity of the O-points of the mode's magnetic island chain is capable of triggering bifurcations of the O-points (which is impossible in the one-harmonic approximation). However, the current drive is incapable of triggering bifurcations of the island X-points. This finding is significant because Bardóczi and Evans recently observed bifurcations of magnetic island chain O-points in the presence of RF current drive in the DIII-D tokamak but did not observe bifurcations of the X-points. Finally, the changes in the topology of the magnetic island flux-surfaces induced by RF current drive are found to facilitate the stabilization of the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1977

Aerodynamic loads near cranks, apexes, and tips of thin, lifting wings in incompressible flow

The calculation of the incompressible and irrotational flow in the vicinity of tips and corners of thin, lifting wings is considered. It is shown that the important characteristics of the flow are governed by an eigenvalue problem, which is nonlinear at the trailing edge because of the shed wake (assumed to be in the wing plane). A new solution method was devised because either the existing methods were not valid for the trailing edge case or they would have required excessive amounts of computer time. The new method, which is fundamentally different than the previous ones, was used to calculate solutions for a number of cases, including some for which correct answers had not previously been obtained. Two of these solutions were used to determine the validity of drag and leading-edge-suction distributions near the tips of a delta wing and a swept wing as calculated by using both the vortex lattice method and a kernel function method. The calculations for the swept wing resolved the question of whether or not the induced drag should be zero at the wing tip.

Medan, R. T.↗

DOE OSTI2016

Estimates of the eigenvalues of operator arising in swelling pressure model

Swelling pressures from materials confined by structures can cause structural deformations and instability. Due to the complexity of interactions between expansive solid and solid-liquid equilibrium, the forces exerting on retaining structures from swelling are highly nonlinear. This work is our initial attempt to study a simplistic spectral problem based on the Euler-elastic beam theory and some simplistic swelling pressure model. In this work estimates of the eigenvalues of some initial/boundary value problem for nonlinear Euler-elastic beam equation are obtained.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗