SEARCH · Search NASA

Results for “Nonlinear solution”

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 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Hierarchically partitioned nonlinear equation solvers

By partitioning solution space into a number of subspaces, a new multiply constrained partitioned Newton-Raphson nonlinear equation solver is developed. Specifically, for a given iteration, each of the various separate partitions are individually and simultaneously controlled. Due to the generality of the scheme, a hierarchy of partition levels can be employed. For finite-element-type applications, this includes the possibility of degree-of-freedom, nodal, elemental, geometric substructural, material and kinematically nonlinear group controls. It is noted that such partitioning can be continuously updated, depending on solution conditioning. In this context, convergence is ascertained at the individual partition level.

Padovan, Joseph↗

NASA NTRS1991

Higher-order harmonics in electromagnetic scattering from a nonlinear anisotropic cylinder

The general solution of the scattering of harmonic signals, by a medium which has nonlinear constitutive relations, contains higher harmonics. The scattering amplitude consists of the linear component and additional correction terms that depend on the amplitude of the incident wave. In this paper, we examine the effects of nonlinearities on the field components and on the radar cross section at the second and third harmonics.

Hasan, Moh'd A.↗

NASA NTRS1972

A finite element procedure for nonlinear prebuckling and initial postbuckling analysis

A procedure cast in a form appropriate to the finite element method is presented for geometrically nonlinear prebuckling and postbuckling structural analysis, including the identification of snap-through type of buckling. The principal features of this procedure are the use of direct iteration for solution of the nonlinear algebraic equations in the prebuckling range, an interpolation scheme for determination of the initial bifurcation point, a perturbation method in definition of the load-displacement behavior through the postbuckling regime, and extrapolation in determination of the limit point for snap-through buckling. Three numerical examples are presented in illustration of the procedure and in comparison with alternative approaches.

Mau, S. T.↗

NASA NTRS1991

On the nonlinear development of the most unstable Goertler vortex mode

The nonlinear development of the most unstable Gortler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Gortler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Gortler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Gortler equations valid for O(1) Gortler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Denier, James P.↗

NASA NTRS1993

On the nonlinear development of the most unstable Goertler vortex mode TLSP - Final report

The nonlinear development of the most unstable Goertler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Goertler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Goertler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Goertler equations valid for O(1) Goertler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Denier, James P.↗

NASA NTRS1988

A nonlinear analysis of an equilibrium craze. I - Problem formulation and solution. II - Simulations of craze and crack growth

A numerical method is developed for the determinations of the equilibrium shape of a craze in an infinite elastic plane whose fibrils exhibit very general nonlinear force-displacement behavior. The problem formulation is based on the superposition of the relevant elasticity Green's function; the solution of the resulting nonlinear problem is obtained by using Picard's successive approximation scheme. The model is used to investigate the effect of nonlinear fibril behavior on the mechanics of craze and crack growth, and the results are compared with the Dugdale model.

Ungsuwarungsri, T.↗

NASA NTRS1992

Chaotic structures of nonlinear magnetic fields. I - Theory. II - Numerical results

A study of the evolutionary properties of nonlinear magnetic fields in flowing MHD plasmas is presented to illustrate that nonlinear magnetic fields may involve chaotic dynamics. It is shown how a suitable transformation of the coupled equations leads to Duffing's form, suggesting that the behavior of the general solution can also be chaotic. Numerical solutions of the nonlinear magnetic field equations that have been cast in the form of Duffing's equation are presented.

Lee, Nam C.↗

NASA NTRS1999

Bidirectional Reflectance of Flat, Optically Thick Particulate Layers: An Efficient Radiative Transfer Solution and Applications to Snow and Soil Surfaces

We describe a simple and highly efficient and accurate radiative transfer technique for computing bidirectional reflectance of a macroscopically flat scattering layer composed of nonabsorbing or weakly absorbing, arbitrarily shaped, randomly oriented and randomly distributed particles. The layer is assumed to be homogeneous and optically semi-infinite, and the bidirectional reflection function (BRF) is found by a simple iterative solution of the Ambartsumian's nonlinear integral equation. As an exact Solution of the radiative transfer equation, the reflection function thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. Since this technique bypasses the computation of the internal radiation field, it is by far the fastest numerical approach available and can be used as an ideal input for Monte Carlo procedures calculating BRFs of scattering layers with macroscopically rough surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web at http://ww,,v.giss.nasa.gov/-crmim/brf.html and can be applied to a wide range of remote sensing, engineering, and biophysical problems. We also examine the potential effect of ice crystal shape on the bidirectional reflectance of flat snow surfaces and the applicability of the Henyey-Greenstein phase function and the 6-Eddington approximation in calculations for soil surfaces.

Mishchenko, Michael I.↗

NASA NTRS1975

Thermomechanical coupling in fatigue fracture of viscoelastic materials

Subcritical crack growth in a linear viscoelastic material subjected to cyclic loading is investigated starting with the thermodynamic power balance. Physically, it is supposed that the subcritical slow crack growth is due to local weakening of the material in the neighborhood of the crack tip. Under fatigue loading it is assumed that all energy dissipation goes into heat and that this heat build-up is the dominant irreversible process governing crack growth. The cycle averaged temperature distribution around the crack tip is obtained from local application of the first law of thermodynamics (i.e., conservation of energy). The analysis leads to the solution of a nonlinear integro-differential equation for crack length as a function of time which is coupled with the local energy equation. A regular perturbation technique is used to obtain an analytical solution which compares very well with experimental results.

Hufferd, W. L.↗

NASA NTRS1995

Development of a linearized unsteady Euler analysis for turbomachinery blade rows

A linearized unsteady aerodynamic analysis for axial-flow turbomachinery blading is described in this report. The linearization is based on the Euler equations of fluid motion and is motivated by the need for an efficient aerodynamic analysis that can be used in predicting the aeroelastic and aeroacoustic responses of blade rows. The field equations and surface conditions required for inviscid, nonlinear and linearized, unsteady aerodynamic analyses of three-dimensional flow through a single, blade row operating within a cylindrical duct, are derived. An existing numerical algorithm for determining time-accurate solutions of the nonlinear unsteady flow problem is described, and a numerical model, based upon this nonlinear flow solver, is formulated for the first-harmonic linear unsteady problem. The linearized aerodynamic and numerical models have been implemented into a first-harmonic unsteady flow code, called LINFLUX. At present this code applies only to two-dimensional flows, but an extension to three-dimensions is planned as future work. The three-dimensional aerodynamic and numerical formulations are described in this report. Numerical results for two-dimensional unsteady cascade flows, excited by prescribed blade motions and prescribed aerodynamic disturbances at inlet and exit, are also provided to illustrate the present capabilities of the LINFLUX analysis.

Verdon, Joseph M.↗

NASA NTRS1974

Stability analysis of nonlinear systems by multiple time scaling

The asymptotic solution for the transient analysis of a general nonlinear system in the neighborhood of the stability boundary was obtained by using the multiple-time-scaling asymptotic-expansion method. The nonlinearities are assumed to be of algebraic nature. Terms of order epsilon to the 3rd power (where epsilon is the order of amplitude of the unknown) are included in the solution. The solution indicates that there is always a limit cycle which is stable (unstable) and exists above (below) the stability boundary if the nonlinear terms are stabilizing (destabilizing). Extension of the solution to include fifth order nonlinear terms is also presented. Comparisons with harmonic balance and with multiple-time-scaling solution of panel flutter equations are also included.

Morino, L.↗

NASA NTRS1975

K/S Lambert problem

The Lambert problem in orbital mechanics is formulated in Kustaanheimo/Stiefel variables. The problem is to determine the required energy and the value of the generalized eccentric anomaly such that a particle at the initial position vector will transfer to the final position vector in a physical time interval. The fictitious time solution results in two nonlinear equations in the two unknowns, energy and fictitious time. The generalized eccentric anomaly solution, however, results in only one nonlinear equation in the one unknown, the eccentric anomaly. This simplification is possible because the energy equation is separable in the eccentric anomaly formulation.

Jezewski, D. J.↗

NASA NTRS1989

Multi-grid for structures analysis

In structural analysis the amount of computational time necessary for a solution is proportional to the number of degrees of freedom times the bandwidth squared. In implicit time analysis, this must be done at each discrete point in time. If, in addition, the problem is nonlinear, then this solution must be iterated at each point in time. If the bandwidth is large, the size of the problem that can be analyzed is severely limited. The multi-grid method is a possible algorithm which can make this solution much more computationally efficient. This method has been used for years in computational fluid mechanics. It works on the fact that relaxation is very efficient on the high frequency components of the solution (nearest neighbor interactions) and not very good on low frequency components of the solution (far interactions). The multi-grid method is then to relax the solution on a particular model until the residual stops changing. This indicates that the solution contains the higher frequency components. A coarse model is then generated for the lower frequency components to the solution. The model is then relaxed for the lower frequency components of the solution. These lower frequency components are then interpolated to the fine model. In computational fluid mechanics the equations are usually expressed as finite differences.

Kascak, Albert F.↗

NASA NTRS2003

On the Hilbert-Huang Transform Data Processing System Development

One of the main heritage tools used in scientific and engineering data spectrum analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). The Fourier view of nonlinear mechanics that had existed for a long time, and the associated FFT (fairly recent development), carry strong a-priori assumptions about the source data, such as linearity and of being stationary. Natural phenomena measurements are essentially nonlinear and nonstationary. A very recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT) proposes a novel approach to the solution for the nonlinear class of spectrum analysis problems. Using the Empirical Mode Decomposition (EMD) followed by the Hilbert Transform of the empirical decomposition data (HT), the HHT allows spectrum analysis of nonlinear and nonstationary data by using an engineering a-posteriori data processing, based on the EMD algorithm. This results in a non-constrained decomposition of a source real value data vector into a finite set of Intrinsic Mode Functions (IMF) that can be further analyzed for spectrum interpretation by the classical Hilbert Transform. This paper describes phase one of the development of a new engineering tool, the HHT Data Processing System (HHTDPS). The HHTDPS allows applying the "T to a data vector in a fashion similar to the heritage FFT. It is a generic, low cost, high performance personal computer (PC) based system that implements the HHT computational algorithms in a user friendly, file driven environment. This paper also presents a quantitative analysis for a complex waveform data sample, a summary of technology commercialization efforts and the lessons learned from this new technology development.

Kizhner, Semion↗

NASA NTRS1984

Forced stationary solutions in a barotropic channel: Multiple equilibria and theory of nonlinear resonance

A barotropic flow over sinusoidal topography in a beta channel in the presence of forcing and dissipation was studied. The work is divided into two parts: (1) computation of the stationary solutions of the model, by Newton's method together with a finite difference approximation, and comparison of solutions with the solutions of the linear theory; (2) development of a nonlinear theory that reproduce the numerical solutions.

Rambaldi, S.↗

DOE OSTI2024

Validation of a Tapered Impact Damper for Traffic Signal Structure Vibrations Using the Method of Harmonic Balance

Significant vibration amplitudes and cycles can be produced when traffic signal structures with low inherent damping are excited near one of their natural frequencies. For the mitigation of wind-induced vibrations, dynamic vibration absorbers coupled to the structure are often used. Here, this research investigates the performance of a tapered impact damper, consisting of a hanging spring-mass oscillator inside a housing capable of reducing vibration amplitude over a broader frequency range than the conventional tuned mass damper. A nonlinear, two degree-of-freedom model is developed with coordinates representing the traffic structure and the tapered impact damper. This research focuses on the application of the harmonic balance method to approximate the periodic solutions of the nonlinear equations to compute the nonlinear dynamics of the damped traffic signal structure. After designing and manufacturing a tapered impact damper, the traffic signal structure is tested with and without the damper using free vibration snapback tests. The experimental frequency and damping backbone curves are used to validate the analytical model, and the effectiveness of the damper is discussed.

42 ENGINEERING↗

NASA NTRS1985

Three-dimensional vibration analysis of a uniform beam with offset inertial masses at the ends

Analysis of a flexible beam with displaced end-located inertial masses is presented. The resulting three-dimensional mode shape is shown to consist of two one-plane bending modes and one torsional mode. These three components of the mode shapes are shown to be linear combinations of trigonometric and hyperbolic sine and cosine functions. Boundary conditions are derived to obtain nonlinear algebraic equations through kinematic coupling of the general solutions of the three governing partial differential equations. A method of solution which takes these boundary conditions into account is also presented. A computer program has been written to obtain unique solutions to the resulting nonlinear algebraic equations. This program, which calculates natural frequencies and three-dimensional mode shapes for any number of modes, is presented and discussed.

Robertson, D. K.↗

NASA NTRS1997

Nonlinear Rescaling and Proximal-Like Methods in Convex Optimization

The nonlinear rescaling principle (NRP) consists of transforming the objective function and/or the constraints of a given constrained optimization problem into another problem which is equivalent to the original one in the sense that their optimal set of solutions coincides. A nonlinear transformation parameterized by a positive scalar parameter and based on a smooth scaling function is used to transform the constraints. The methods based on NRP consist of sequential unconstrained minimization of the classical Lagrangian for the equivalent problem, followed by an explicit formula updating the Lagrange multipliers. We first show that the NRP leads naturally to proximal methods with an entropy-like kernel, which is defined by the conjugate of the scaling function, and establish that the two methods are dually equivalent for convex constrained minimization problems. We then study the convergence properties of the nonlinear rescaling algorithm and the corresponding entropy-like proximal methods for convex constrained optimization problems. Special cases of the nonlinear resealing algorithm are presented. In particular a new class of exponential penalty-modified barrier functions methods is introduced.

Polyak, Roman↗