Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear approximation”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

On the design of nonlinear controllers for flight control systems

A method of approximate input-output linearization by dynamic state feedback is presented, with the flight control of VTOL aircraft used as an example. It is shown that the closed loop system has a graceful degradation of performance as the moment-to-force coupling is increased. It is also demonstrated that the approach proposed here leads to an asymptotically stable closed loop system with guaranteed bounds on the tracking error caused by the nonminumum phase character of the system.

Hauser, John↗
NASA NTRS1995

The large discretization step method for time-dependent partial differential equations

A new method for the acceleration of linear and nonlinear time dependent calculations is presented. It is based on the Large Discretization Step (LDS) approximation, defined in this work, which employs an extended system of low accuracy schemes to approximate a high accuracy discrete approximation to a time dependent differential operator. Error bounds on such approximations are derived. These approximations are efficiently implemented in the LDS methods for linear and nonlinear hyperbolic equations, presented here. In these algorithms the high and low accuracy schemes are interpreted as the same discretization of a time dependent operator on fine and coarse grids, respectively. Thus, a system of correction terms and corresponding equations are derived and solved on the coarse grid to yield the fine grid accuracy. These terms are initialized by visiting the fine grid once in many coarse grid time steps. The resulting methods are very general, simple to implement and may be used to accelerate many existing time marching schemes.

Haras, Zigo↗
NASA NTRS1981

DC-magnetic field vector measurement

A magnetometer experiment was designed to determine the local magnetic field by measuring the total of the Earth's magnetic field and that of an unknown spacecraft. The measured field vector components are available to all onboard experiments via the Spacelab command and data management system. The experiment consists of two parts, an electronic box and the magnetic field sensor. The sensor includes three independent measuring flux-gate magnetometers, each measuring one component. The physical background is the nonlinearity of the B-H curve of a ferrite material. Two coils wound around a ferrite rod are necessary. One of them, a tank coil, pumps the ferrite rod at approximately 20 kilohertz. As a consequence of the nonlinearity, many harmonics can be produced. The second coil (i.e., the detection coil) resonates to the first harmonic. If an unknown dc or low-frequency magnetic field exists, the amplitude of the first harmonic is a measure for the unknown magnetic field. The voltages detected by the sensors are to be digitized and transferred to the command and data management system.

Schmidt, R.↗
NASA NTRS1995

Nonlinear structure formation in flat cosmological models

This paper describes the formation of nonlinear structure in flat (zero curvature) Friedmann cosmological models. We consider models with two components: the usual nonrelativistic component that evolves under gravity and eventually forms the large-scale structure of the universe, and a uniform dark matter component that does not clump under gravity, and whose energy density varies with the scale factor a(t) like a(t)(sup -n), where n is a free parameter. Each model is characterized by two parameters: the exponent n and the present density parameter Omega(sub 0) of the nonrelativistic component. The linear perturbation equations are derived and solved for these models, for the three different cases n = 3, n is greater than 3, and n is less than 3. The case n = 3 is relevant to model with massive neutrinos. The presence of the uniform component strongly reduces the growth of the perturbation compared with the Einstein-de Sitter model. We show that the Meszaros effect (suppression of growth at high redshift) holds not only for n = 4, radiation-dominated models, but for all models with n is greater than 3. This essentially rules out any such model. For the case n is less than 3, we numerically integrate the perturbation equations from the big bang to the present, for 620 different models with various values of Omega(sub 0) and n. Using these solutions, we show that the function f(Omega(sub 0), n) = (a/delta(sub +))d(delta)(sub +)/da, which enters in the relationship between the present density contrast delta(sub 0) and peculiar velocity field u(sub 0) is essentially independent of n. We derive approximate solutions for the second-order perturbation equations. These second-order solutions are tested against the exact solutions and the Zel'dovich approximation for spherically symmetric perturbations in the marginally nonlinear regime (the absolute value of delta is less than or approximately 1). The second-order and Zel'dovich solutions have comparable accuracy, significantly higher than the accuracy of the linear solutions. We then investigate the dependence of the delta(sub 0) - u(sub 0) relationship upon the value of n in the nonlinear regime using the second-order solutions for marginally nonlinear, general perturbations, and the exact solutions for strongly nonlinear, spherically symmetric perturbations. In both cases, we find that the delta(sub 0) - u(sub 0) relationship remains independent of n. We speculate that this result extends to strongly nonlinear, general perturbations as well. This eliminates any hope to determine the presence of the uniform component or the value of n using dynamical methods. Finally, we compute the nonlinear evolution of the skewness of the distribution of values of delta, assuming Gaussian initial conditions. We find that the skewness is not only independent of n, but also of Omega(sub 0). Thus the skewness cannot be used to discriminate among various models with Gaussian initial conditions. However, it can be used for testing the Gaussianity of the initial conditions themselves. We conclude that the uniform component leaves no observable signature in the present large-scale structure of the universe. To determine its presence and nature, we must investigate the relationship between the past and present universe, using redshift-dependent tests.

Martel, Hugo↗
NASA NTRS1995

A hybrid neurocomputing/numerical strategy for nonlinear structural analysis

A hybrid neurocomputing/numerical strategy is presented for geometrically nonlinear analysis of structures. The strategy combines model-free data processing capabilities of computational neural networks with a Pade approximants-based perturbation technique to predict partial information about the nonlinear response of structures. In the hybrid strategy, multilayer feedforward neural networks are used to extend the validity of solutions by using training samples produced by Pade approximations to the Taylor series expansion of the response function. The range of validity of the training samples is taken to be the radius of convergence of Pade approximants and is estimated by setting a tolerance on the diverging approximants. The norm of residual vector of unbalanced forces in a given element is used as a measure to assess the quality of network predictions. To further increase the accuracy and the range of network predictions, additional training data are generated by either applying linear regression to weight matrices or expanding the training data by using predicted coefficients in a Taylor series. The effectiveness of the hybrid strategy is assessed by performing large-deflection analysis of a doubly-curved composite panel with a circular cutout, and postbuckling analyses of stiffened composite panels subjected to an in-plane edge shear load. In all the problems considered, the hybrid strategy is used to predict selective information about the structural response, namely the total strain energy and the maximum displacement components only.

Szewczyk, Z. Peter↗
NASA NTRS2006

Initial On-Orbit Spatial Resolution Characterization of OrbView-3 Panchromatic Images

Characterization was conducted under the Memorandum of Understanding among Orbital Sciences Corp., ORBIMAGE, Inc., and NASA Applied Sciences Directorate. Acquired five OrbView-3 panchromatic images of the permanent Stennis Space Center edge targets painted on a concrete surface. Each image is available at two processing levels: Georaw and Basic. Georaw is an intermediate image in which individual pixels are aligned by a nominal shift in the along-scan direction to adjust for the staggered layout of the panchromatic detectors along the focal plane array. Georaw images are engineering data and are not delivered to customers. The Basic product includes a cubic interpolation to align the pixels better along the focal plane and to correct for sensor artifacts, such as smile and attitude smoothing. This product retains satellite geometry - no rectification is performed. Processing of the characterized images did not include image sharpening, which is applied by default to OrbView-3 image products delivered by ORBIMAGE to customers. Edge responses were extracted from images of tilted edges in two directions: along-scan and cross-scan. Each edge response was approximated with a superposition of three sigmoidal functions through a nonlinear least-squares curve-fitting. Line Spread Functions (LSF) were derived by differentiation of the analytical approximation. Modulation Transfer Functions (MTF) were obtained after applying the discrete Fourier transform to the LSF.

Blonski, Slawomir↗
NASA NTRS1993

A new flux conserving Newton's method scheme for the two-dimensional, steady Navier-Stokes equations

A new numerical method is developed for the solution of the two-dimensional, steady Navier-Stokes equations. The method that is presented differs in significant ways from the established numerical methods for solving the Navier-Stokes equations. The major differences are described. First, the focus of the present method is on satisfying flux conservation in an integral formulation, rather than on simulating conservation laws in their differential form. Second, the present approach provides a unified treatment of the dependent variables and their unknown derivatives. All are treated as unknowns together to be solved for through simulating local and global flux conservation. Third, fluxes are balanced at cell interfaces without the use of interpolation or flux limiters. Fourth, flux conservation is achieved through the use of discrete regions known as conservation elements and solution elements. These elements are not the same as the standard control volumes used in the finite volume method. Fifth, the discrete approximation obtained on each solution element is a functional solution of both the integral and differential form of the Navier-Stokes equations. Finally, the method that is presented is a highly localized approach in which the coupling to nearby cells is only in one direction for each spatial coordinate, and involves only the immediately adjacent cells. A general third-order formulation for the steady, compressible Navier-Stokes equations is presented, and then a Newton's method scheme is developed for the solution of incompressible, low Reynolds number channel flow. It is shown that the Jacobian matrix is nearly block diagonal if the nonlinear system of discrete equations is arranged approximately and a proper pivoting strategy is used. Numerical results are presented for Reynolds numbers of 100, 1000, and 2000. Finally, it is shown that the present scheme can resolve the developing channel flow boundary layer using as few as six to ten cells per channel width, depending on the Reynolds number.

Scott, James R.↗
NASA NTRS1994

The Local Supercluster as a test of cosmological models

The Local Supercluster kinematic properties (the Local Group infall toward the Virgo Cluster and galaxy density distribution about the Virgo Cluster) in various cosmological models are examined utilizing large-scale N-body (PM) simulations 500(exp 3) cells, 250(exp 3) particles, and box size of 400 h(exp -1) Mpc) and are compared to observations. Five models are investigated: (1) the standard, Cosmic Background Explorer Satellite (COBE)-normalized cold dark matter (CDM) model with omega = 1, h = 0.5, and sigma(sub 8) = 1.05; (2) the standard Hot Dark Matter (HDM) model with omega = 1, h = 0.75, and sigma(sub 8) = 1; (3) the tilted CDM model with omega = 1, h = 0.5, n = 0.7, and sigma(sub 8) = 0.5; (4) a CDM + lambda model with omega = 0.3, lambda = 0.7, h = 2/3, and sigma(sub 8) = 2/3; (5) the PBI model with omega = 0.2, h = 0.8, x = 0.1, m = -0.5, and sigma(sub 8) = 0.9. Comparison of the five models with the presently available observational measurements v(sub LG) = 85 - 305 km/s (with mean of 250 km/s), delta(n(sub g))/(n(sub g)-bar) = 1.40 + or - 0.35) suggests that an open universe with omega approximately 0.5 (with or without lambda) and sigma(sub 8) approximately 0.8 is preferred, with omega = 0.3-1.0 (with or without lambda) and sigma(sub 8) = 0.7-1.0 being the acceptable range. At variance with some previous claims based on either direct N-body or spherical nonlinear approaches, we find that a flat model with sigma(sub 8) approximately 0.7-1.0 seems to be reasonably consistent with observations. However, if one favors the low limit of v(sub LG) = 85 km/s, then an omega approximately 0.2-0.3 universe seems to provide a better fit, and flat (omega = 1) models are ruled out at approximately 95% confidence level. On the other hand, if the high limit of v(sub LG) = 350 km/s is closer to the truth, then it appears that omega approximately 0.7-0.8 is more consistent. This test is insensitive to the shape of the power spectrum, but rather sensitive to the normalization of the perturbation amplitude on the relevant scale (e.g., sigma(sub 8)) and omega. We find that neither linear nor nonlinear relations (with spherical symmetry) are good approximations for the relation between radial peculiar velocity and density perturbation, i.e., nonspherical effects and gravitational tidal field are important. The derived omega using either of the two relations is underestimated. In some cases, this error is as large as a factor of 2-4.

Cen, Renyue↗
NASA NTRS1984

A finite element method for nonlinear forced vibrations of beams

Techniques for defining a finite element model (FEM) for analysis of nonlinear vibrations in beam structures subjected to harmonic excitation are presented. The resulting model covers longitudinal deformation and inertial effects. The nonlinear oscillations of a beam element under forced excitation are modeled by a harmonic force matrix based on first order approximations of the Jacobian elliptic forcing function. Harmonic force and nonlinear stiffness matrices are derived and the nonlinear forced responses of beams are calculated under various boundary conditions. The results of FEM computations for simply-supported and clamped beams show that midplane stretching caused by large deflections increases the nonlinearity. Axially-restrained beams experience only hardening nonlinearity, while axially-free beams have reduced nonlinearity in deformation and inertia and an increase in linearity due to large deflection.

Mei, C.↗
NASA NTRS1984

The application of satellite data to study the effects of latent heat release on cyclones

Generalized energetics were studied for nonlinear inviscid symmetric instability (SI). It was found that the linear theory fails to predict the stability in certain cases where the basic state is transitional between stability and instability. The initial growth of the SI perturbations can be fairly well approximated by linear theory, but the long time nonlinear evaluations will be bonded energetically if the SI region is finite. However, a further extension of the energetics to conditional symmetric instability (CSI) shows that the nonlinear evolution of circulation will energetically depend much more on the precipitation in a complicated way. By treating the latent heat as a source which is implicitly related to the motion field, the existence, uniqueness and stability of steady viscous (CSI) circulations are studied. Viscous CSI circulations are proved to be unique and asymptotically stable when the heat sources are weak and less sensitive to the motion perturbations. By considering the fact that moist updrafts are narrow and using eddy viscosity of 0(1,000 m squared/s) the stability criterion suggests that some frontal rainbands were probably dominated by the CSI mechanism even in their mature quasi-steady stage.

Clark, J. H. E.↗
NASA NTRS1985

An Algorithm for Efficient Maximum Likelihood Estimation and Confidence Interval Determination in Nonlinear Estimation Problems

An algorithm for maximum likelihood (ML) estimation is developed with an efficient method for approximating the sensitivities. The algorithm was developed for airplane parameter estimation problems but is well suited for most nonlinear, multivariable, dynamic systems. The ML algorithm relies on a new optimization method referred to as a modified Newton-Raphson with estimated sensitivities (MNRES). MNRES determines sensitivities by using slope information from local surface approximations of each output variable in parameter space. The fitted surface allows sensitivity information to be updated at each iteration with a significant reduction in computational effort. MNRES determines the sensitivities with less computational effort than using either a finite-difference method or integrating the analytically determined sensitivity equations. MNRES eliminates the need to derive sensitivity equations for each new model, thus eliminating algorithm reformulation with each new model and providing flexibility to use model equations in any format that is convenient. A random search technique for determining the confidence limits of ML parameter estimates is applied to nonlinear estimation problems for airplanes. The confidence intervals obtained by the search are compared with Cramer-Rao (CR) bounds at the same confidence level. It is observed that the degree of nonlinearity in the estimation problem is an important factor in the relationship between CR bounds and the error bounds determined by the search technique. The CR bounds were found to be close to the bounds determined by the search when the degree of nonlinearity was small. Beale's measure of nonlinearity is developed in this study for airplane identification problems; it is used to empirically correct confidence levels for the parameter confidence limits. The primary utility of the measure, however, was found to be in predicting the degree of agreement between Cramer-Rao bounds and search estimates.

Murphy, Patrick Charles↗
NASA NTRS2022

Alternative regularizations for Outer-Approximation algorithms for convex MINLP

In this work, we extend the regularization framework from Kronqvist et al. (Math Program 180(1):285–310, 2020) by incorporating several new regularization functions and develop a regularized single-tree search method for solving convex mixed-integer nonlinear programming (MINLP) problems. We propose a set of regularization functions based on distance metrics and Lagrangean approximations, used in the projection problem for finding new integer combinations to be used within the Outer-Approximation (OA) method. The new approach, called Regularized Outer-Approximation (ROA), has been implemented as part of the open-source Mixed-integer nonlinear decomposition toolbox for Pyomo—MindtPy. We compare the OA method with seven regularization function alternatives for ROA. Moreover, we extend the LP/NLP Branch and Bound method proposed by Quesada and Grossmann (Comput Chem Eng 16(10–11):937–947, 1992) to include regularization in an algorithm denoted RLP/NLP. We provide convergence guarantees for both ROA and RLP/NLP. Finally, we perform an extensive computational experiment considering all convex MINLP problems in the benchmark library MINLPLib. The computational results show clear advantages of using regularization combined with the OA method.

Convex Mixed-integer nonlinear programming↗
NASA NTRS1988

Nonlinearities in spacecraft structural dynamics

In considering nonlinearities in spacecraft structural dynamics, the following are examined: (1) SCOLE Configuration-Equations of Motion; (2) Modeling Error Sources; (3) Approximate Solutions; (4) Comparison of Model Accuracy; (5) Linear and Nonlinear Damping; (6) Experimental Results; and, (7) Future Work.

Taylor, Larry↗
NASA NTRS1992

Collapse of composite tubes under end moments

Cylindrical tubes of moderate wall thickness such as those proposed for the original space station truss, may fail due to the gradual collapse of the tube cross section as it distorts under load. Sometimes referred to as the Brazier instability, it is a nonlinear phenomenon. This paper presents an extension of an approximate closed form solution of the collapse of isotropic tubes subject to end moments developed by Reissner in 1959 to include specially orthotropic material. The closed form solution was verified by an extensive nonlinear finite element analysis of the collapse of long tubes under applied end moments for radius to thickness ratios and composite layups in the range proposed for recent space station truss framework designs. The finite element analysis validated the assumption of inextensional deformation of the cylindrical cross section and the approximation of the material as specially orthotropic.

Stockwell, Alan E.↗
NASA NTRS1994

Effect of nonlinear interactions on p-mode frequencies and line widths

We calculate the effect of nonlinear interactions among solar acoustic modes upon the modal frequencies and energy loss rates (or line widths). The frequency shift for a radial p-mode of frequency 3 mHz is found to be about -0.5 microHz. The magnitude of nonlinear frequency shift increases more rapidly with frequency than the inverse mode mass (mode mass is defined as the ratio of energy in the mode to its surface velocity amplitude squared). This frequency shift is primarily due to nonresonant three-mode interactions and is dominated by high l surface gravity waves (f-modes) and p-modes. The line width of a radial p-mode of frequency 3 mHz, due to resonant nonlinear interactions, is about 0.3 microHz. This result is consistent with that of Kumar and Goldreich (1989). We also find, in agreement with these authors, that the most important nonlinear interactions of trapped p-modes involve f-modes and high-frequency p-modes (frequency greater than about 5 mHz) which propagate in the solar photosphere. Thus, using the arguments advanced by Kumar & Goldreich (1989), we conclude that nonlinear couplings cannot saturate the overstable solar p-modes at their small observed amplitudes. Both the nonlinear frequency shifts and line widths, at a fixed frequency, are proportional to the inverse of mode mass which for modes of degree greater than about 100 is approximately l(exp 0.8). Therefore, the frequency of an f-mode of l = 1000, due to nonlinear interactions, is decreased by approximately 0.4%.

Kumar, Pawan↗
NASA NTRS1976

Approximation concepts for efficient structural synthesis

It is shown that efficient structural synthesis capabilities can be created by using approximation concepts to mesh finite element structural analysis methods with nonlinear mathematical programming techniques. The history of the application of mathematical programming techniques to structural design optimization problems is reviewed. Several rather general approximation concepts are described along with the technical foundations of the ACCESS 1 computer program, which implements several approximation concepts. A substantial collection of structural design problems involving truss and idealized wing structures is presented. It is concluded that since the basic ideas employed in creating the ACCESS 1 program are rather general, its successful development supports the contention that the introduction of approximation concepts will lead to the emergence of a new generation of practical and efficient, large scale, structural synthesis capabilities in which finite element analysis methods and mathematical programming algorithms will play a central role.

Schmit, L. A., Jr.↗