SEARCH · Engineering Papers

Results for “Nonlinear equations”

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 829 records · Page 46

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Parallel processors and nonlinear structural dynamics algorithms and software

Techniques are discussed for the implementation and improvement of vectorization and concurrency in nonlinear explicit structural finite element codes. In explicit integration methods, the computation of the element internal force vector consumes the bulk of the computer time. The program can be efficiently vectorized by subdividing the elements into blocks and executing all computations in vector mode. The structuring of elements into blocks also provides a convenient way to implement concurrency by creating tasks which can be assigned to available processors for evaluation. The techniques were implemented in a 3-D nonlinear program with one-point quadrature shell elements. Concurrency and vectorization were first implemented in a single time step version of the program. Techniques were developed to minimize processor idle time and to select the optimal vector length. A comparison of run times between the program executed in scalar, serial mode and the fully vectorized code executed concurrently using eight processors shows speed-ups of over 25. Conjugate gradient methods for solving nonlinear algebraic equations are also readily adapted to a parallel environment. A new technique for improving convergence properties of conjugate gradients in nonlinear problems is developed in conjunction with other techniques such as diagonal scaling. A significant reduction in the number of iterations required for convergence is shown for a statically loaded rigid bar suspended by three equally spaced springs.

Belytschko, Ted↗

NASA NTRS1987

A linearized Euler analysis of unsteady flows in turbomachinery

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗

NASA NTRS1989

Calculation of unsteady flows in turbomachinery using the linearized Euler equations

Hall, Kenneth C.↗

NASA NTRS1977

A numerical scheme to solve unstable boundary value problems

The considered scheme makes it possible to determine an unstable steady state solution in cases in which, because of lack of symmetry, such a solution cannot be obtained analytically, and other time integration or relaxation schemes, because of instability, fail to converge. The iterative solution of a single complex equation is discussed and a nonlinear system of equations is considered. Described applications of the scheme are related to a steady state solution with shear instability, an unstable nonlinear Ekman boundary layer, and the steady state solution of a baroclinic atmosphere with asymmetric forcing. The scheme makes use of forward and backward time integrations of the original spatial differential operators and of an approximation of the adjoint operators. Only two computations of the time derivative per iteration are required.

Kalnay-Rivas, E.↗

DOE OSTI2021

IMEX Runge-Kutta Parareal for Non-diffusive Equations

Parareal is a widely studied parallel-in-time method that can achieve meaningful speedup on certain problems. However, it is well known that the method typically performs poorly on non-diffusive equations. This paper analyzes linear stability and convergence for IMEX Runge-Kutta Parareal methods on non-diffusive equations. By combining standard linear stability analysis with a simple convergence analysis, we find that certain Parareal configurations can achieve parallel speedup on non-diffusive equations. These stable configurations all posses low iteration counts, large block sizes, and a large number of processors. Numerical examples using the nonlinear Schrodinger equation demonstrate the analytical conclusions.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1981

Linear and nonlinear aspects of snow albedo feedbacks in atmospheric models

Namias' hypothesis, that anomalous snowcover on the eastern side of the North American continent can generate an anomalous east coast low pressure system and an anomalous inland high pressure system, is consistent with the time-averaged anomalous response from a nonlinear, primitive equation channel model with an idealized, flat land-sea arrangement. An attempt to understand and describe this anomalous response in the nonlinear model as a linear response to anomalous diabatic heating was largely unsuccessful, primarily because the anomalous eddy fluxes were also important. This unsuccessful attempt to describe the nonlinear model's time averages by linear theory then motivated several comparisons between linear and nonlinear severely truncated quasi-geostrophic models. It was also found in these models that the eddy fluxes were extremely important for forcing or dissipating the stationary eddies.

Roads, J. O.↗

NASA NTRS1985

Parametric Stiffness Control of Flexible Structures

An unconventional method for control of flexible space structures using feedback control of certain elements of the stiffness matrix is discussed. The advantage of using this method of configuration control is that it can be accomplished in practical structures by changing the initial stress state in the structure. The initial stress state can be controlled hydraulically or by cables. The method leads, however, to nonlinear control equations. In particular, a long slender truss structure under cable induced initial compression is examined. both analytical and numerical analyses are presented. Nonlinear analysis using center manifold theory and normal form theory is used to determine criteria on the nonlinear control gains for stable or unstable operation. The analysis is made possible by the use of the exact computer algebra system MACSYMA.

Moon, F. C.↗

NASA NTRS1987

Active stability augmentation of large space structures: A stochastic control problem

A problem in SCOLE is that of slewing an offset antenna on a long flexible beam-like truss attached to the space shuttle, with rather stringent pointing accuracy requirements. The relevant methodology aspects in robust feedback-control design for stability augmentation of the beam using on-board sensors is examined. It is framed as a stochastic control problem, boundary control of a distributed parameter system described by partial differential equations. While the framework is mathematical, the emphasis is still on an engineering solution. An abstract mathematical formulation is developed as a nonlinear wave equation in a Hilbert space. That the system is controllable is shown and a feedback control law that is robust in the sense that it does not require quantitative knowledge of system parameters is developed. The stochastic control problem that arises in instrumenting this law using appropriate sensors is treated. Using an engineering first approximation which is valid for small damping, formulas for optimal choice of the control gain are developed.

Balakrishnan, A. V.↗

NASA NTRS1991

Model atmospheres and X-ray spectra of bursting neutron stars

Model atmosphere equations are derived which correspond to plane-parallel nongray atmospheres of very hot neutron stars in hydrostatic and radiative equilibrium. The nonlinear transfer equation used implements an exact photon redistribution function which precisely traces even those scattering events having large photon-electron energy exchange. The numerical results thus obtained encompass tables of the surface fluxes for 20 model atmospheres whose T(eff) ranges from 6 to 30 million K. The spectra of high gravity models are virtually identical with the blackbody spectrum shifted toward higher energies, while models approaching the Eddington limit develop a large low-energy hump.

Madej, Jerzy↗

NASA NTRS2000

A Low Order and a High Order Solution for a Converging-Diverging Nozzle Problem

The Computational Aeroacoustics Workshops on Benchmark problems are conducted in order to generate solutions with a variety of methods for problems that are chosen both to be representative of significant and relevant acoustic phenomena. It is generally recognized that CAA applications require the resolution in time and in space of solution details over a significantly broader range of scales than typical computational fluid dynamics applications. These requirements are forcing the development of new methodologies for CAA applications. This paper presents numerical results from a third and a seventh order algorithm for the propagation of an acoustic signal through a converging-diverging nozzle. These two algorithms are members of a new family of high accuracy methods that have the same order accuracy in both space and time, and are an extension of previous work for linearized Euler equations to fully nonlinear time dependent problems. The simulations are all with the fully nonlinear quasi-1D Euler equations for the total solution, which includes both the fluid dynamics and the acoustics. The acoustic solution is obtained from the time dependent nonlinear solution by subtracting the steady solution.

Goodrich, John W.↗

DOE OSTI2024

Three-dimensional continuum point cloud method for large deformation and its verification

This study presents a strong form based meshfree collocation method, which is named Continuum Point Cloud Method, to solve nonlinear field equations derived from classical mechanics for deformed bodies in three-dimensional Euclidean space. The method and its implementation are benchmarked against a nonlinear vector field using manufactured solutions. The analysis of mechanical fields firstly focuses on the study of St. Venant Kirchhoff and compressible neo-Hookean materials. Results for various initial boundary value problems are presented, including benchmark cases involving unidirectional tension and simple shear. Subsequently, the study concludes with an analysis of a displacement-controlled simulation of a compressible neo-Hookean material, specifically a bar that is pulled to 50% of its original length and rotated 90°. The pure tension case yields a 1.5% error in displacement between computed and expected values and a combined tension and torsion loading case provides further insight into material behavior under complex loading conditions. The resulting normal axial and transverse stress-strain curves are also presented. Lastly, the consistency and robustness of the proposed nonlinear numerical schemes are successfully demonstrated through various numerical experiments.

Compressible neo-Hookean materials↗

DOE OSTI2021

Soliton Generation in Negative Thermal Expansion Materials

Strain solitons have been observed statically in several 2D materials and dynamically in substrate materials using ultrafast laser pulses. The latter case relies on lattice relaxation in response to ultrafast heating in a light-absorbing transducer material, a process which is sensitive to the thermal expansion coefficient. Here we consider an unusual case where the sign of the thermal expansion coefficient is negative, a scenario which is experimentally feasible in light of rapid and recent advances in the discovery of negative thermal expansion materials. We present numerical solutions to a nonlinear differential equation which has been repeatedly demonstrated to quantitatively model experimental data and discuss the salient results using realistic parameters for material linear and nonlinear elasticity. The solitons that emerge from the initial value problem with negative and positive thermal expansion are qualitatively different in several ways. The new case of negative thermal expansion gives rise to a nearly-periodic soliton train with chirped profile and free of an isolated shock front. We suggest this unanticipated result may be realized experimentally and assess the potential for certain applications of this generic effect.

36 MATERIALS SCIENCE↗

NASA NTRS1980

Flutter analysis of an airplane with multiple structural nonlinearities in the control system

Experience has shown that the flutter prediction process for airplanes can be greatly affected by strong concentrated nonlinearities which may be localized in the linking elements of the control mechanism, in the pivot joints of variable-sweep-wing systems, and in the connecting points between wing and pylon-mounted external stores. The principle of equivalent linearization offers an efficent possibility for solving the related nonlinear flutter equations in the frequency domain as a complement to the well-known time domain procedures. Taking as an example an airplane with nonlinear control characteristics, it is demonstrated how the equivalent linearization approach can be extended to rather complicated systems with multiple sets of strongly interacting, concentrated nonlinearities.

Breitbach, E. J.↗

NASA NTRS1984

Optical systolic solutions of linear algebraic equations

The philosophy and data encoding possible in systolic array optical processor (SAOP) were reviewed. The multitude of linear algebraic operations achievable on this architecture is examined. These operations include such linear algebraic algorithms as: matrix-decomposition, direct and indirect solutions, implicit and explicit methods for partial differential equations, eigenvalue and eigenvector calculations, and singular value decomposition. This architecture can be utilized to realize general techniques for solving matrix linear and nonlinear algebraic equations, least mean square error solutions, FIR filters, and nested-loop algorithms for control engineering applications. The data flow and pipelining of operations, design of parallel algorithms and flexible architectures, application of these architectures to computationally intensive physical problems, error source modeling of optical processors, and matching of the computational needs of practical engineering problems to the capabilities of optical processors are emphasized.

Neuman, C. P.↗

NASA NTRS2006

Rapid Aeroelastic Analysis of Blade Flutter in Turbomachines

The LINFLUX-AE computer code predicts flutter and forced responses of blades and vanes in turbomachines under subsonic, transonic, and supersonic flow conditions. The code solves the Euler equations of unsteady flow in a blade passage under the assumption that the blades vibrate harmonically at small amplitudes. The steady-state nonlinear Euler equations are solved by a separate program, then equations for unsteady flow components are obtained through linearization around the steady-state solution. A structural-dynamics analysis (see figure) is performed to determine the frequencies and mode shapes of blade vibrations, a preprocessor interpolates mode shapes from the structural-dynamics mesh onto the LINFLUX computational-fluid-dynamics mesh, and an interface code is used to convert the steady-state flow solution to a form required by LINFLUX. Then LINFLUX solves the linearized equations in the frequency domain to calculate the unsteady aerodynamic pressure distribution for a given vibration mode, frequency, and interblade phase angle. A post-processor uses the unsteady pressures to calculate generalized aerodynamic forces, response amplitudes, and eigenvalues (which determine the flutter frequency and damping). In comparison with the TURBO-AE aeroelastic-analysis code, which solves the equations in the time domain, LINFLUX-AE is 6 to 7 times faster.

Trudell, J. J.↗

NASA NTRS1975

Prediction of axial-flow instabilities in a turbojet engine by use of a multistage compressor simulation on the digital computer

A method of estimating the undistorted stall line for an axial-flow compressor by using the digital computer is presented. The method involves linearization of nonlinear dynamic equations about an operating point on a speed line, and then application of the first method of Lyapunov to determine the stability of the nonlinear system from the stability of the linear system. The method is applied to a simulation of the J85 compressor, which utilizes stage stacking and lumped volume techniques for the interstage regions to simulate steady-state and dynamic compressor performance. The stability boundary predicted by the digital simulation compares quite well with the stall line predicted by a dynamic simulation of the J85 compressor programmed on the analog computer. Since previous studies have shown that the analog-predicted stall line agrees well with the stall line of the compressor, the digital method presented is also a good means of estimating the stall line.

Daniele, C. J.↗

NASA NTRS1989

A stable numerical solution method in-plane loading of nonlinear viscoelastic laminated orthotropic materials

In response to the tremendous growth in the development of advanced materials, such as fiber-reinforced plastic (FRP) composite materials, a new numerical method is developed to analyze and predict the time-dependent properties of these materials. Basic concepts in viscoelasticity, laminated composites, and previous viscoelastic numerical methods are presented. A stable numerical method, called the nonlinear differential equation method (NDEM), is developed to calculate the in-plane stresses and strains over any time period for a general laminate constructed from nonlinear viscoelastic orthotropic plies. The method is implemented in an in-plane stress analysis computer program, called VCAP, to demonstrate its usefulness and to verify its accuracy. A number of actual experimental test results performed on Kevlar/epoxy composite laminates are compared to predictions calculated from the numerical method.

Gramoll, K. C.↗

NASA NTRS1996

Numerical Study of Multigrid Methods with Various Smoothers for the Elliptical Grid Generation Equations

A robust solver for the elliptic grid generation equations is sought via a numerical study. The system of PDEs is discretized with finite differences, and multigrid methods are applied to the resulting nonlinear algebraic equations. Multigrid iterations are compared with respect to the robustness and efficiency. Different smoothers are tried to improve the convergence of iterations. The methods are applied to four 2D grid generation problems over a wide range of grid distortions. The results of the study help to select smoothing schemes and the overall multigrid procedures for elliptic grid generation.

Golik, W. L.↗