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 1,009 records · Page 56

STABILITY PROBLEMS IN MISSILE STRUCTURES

A number of specific buckling problems relating to missile propellant tanks and transition sections are enumerated. In some cases approximate techniques used by missile designers to determine stability are mentioned. The inadequacies of these approximate methods and the unavailability of appropriate usable information is stressed. Furthermore, the mathematical complexities of nonlinear buckling problems are considered. Some possibly fruitful areas of mathematics to simplify these complexities are discussed with the hope that research in these areas will provide a batter means through which the analyst may supply the designer with more realistic and accurate theoretical data.

MISSILE↗

A computer program for the geometrically nonlinear static and dynamic analysis of arbitrarily loaded shells of revolution, theory and users manual

A digital computer program known as SATANS (static and transient analysis, nonlinear, shells) for the geometrically nonlinear static and dynamic response of arbitrarily loaded shells of revolution is presented. Instructions for the preparation of the input data cards and other information necessary for the operation of the program are described in detail and two sample problems are included. The governing partial differential equations are based upon Sanders' nonlinear thin shell theory for the conditions of small strains and moderately small rotations. The governing equations are reduced to uncoupled sets of four linear, second order, partial differential equations in the meridional and time coordinates by expanding the dependent variables in a Fourier sine or cosine series in the circumferential coordinate and treating the nonlinear modal coupling terms as pseudo loads. The derivatives with respect to the meridional coordinate are approximated by central finite differences, and the displacement accelerations are approximated by the implicit Houbolt backward difference scheme with a constant time interval. The boundaries of the shell may be closed, free, fixed, or elastically restrained. The program is coded in the FORTRAN 4 language and is dimensioned to allow a maximum of 10 arbitrary Fourier harmonics and a maximum product of the total number of meridional stations and the total number of Fourier harmonics of 200. The program requires 155,000 bytes of core storage.

Ball, R. E.↗

Electrostatic analysis of charge-coupled structures

An analysis is presented which is based upon a numerical solution of Gauss's law for the multiple dielectric semiconductor and electrode structure. The formulation resulted in a large set of equations, usually nonlinear, which were solved by an iterative scheme based upon the Gauss-Seidel method employing a relaxation parameter. Proper consideration of the nonlinearity of the space-charge in the semiconductor allowed convergent solutions which can give reasonable approximations to the channel potential profiles though not as good for estimating the space charge itself. The results of this work, which are three computer programs, are listed. None of the programs give the electric field explicitly as output data. However, when this information is desired, the programs can be expanded to include field estimations based on a polynominal fit to the potential field, or they may serve as subprograms for main programs dealing with general aspects of charge transfer.

Gassaway, J. D.↗

Hydromagnetic waves and cosmic ray diffusion theory

Pitch angle diffusion of cosmic rays in hydromagnetic wave fields is considered strictly within the quasilinear approximation. It is shown that the popular assumption of an isotropic power spectrum tensor of magnetic fluctuations requires in this case equal forms and magnitudes of Alfven and magnetosonic wave spectra - a situation which is generally unlikely. The relative contributions to the pitch angle diffusion coefficient from the cyclotron resonances and Landau resonance due to the different types of waves are evaluated for a typical situation in the solar wind. Since in this approximation also the Landau resonance does not lead to particle reflections a proper consideration of the nonlinear particle orbits is indeed necessary to overcome the well known difficulties of quasilinear scattering theory for cosmic rays near 90 degrees pitch angle.

Lee, M. A.↗

Iterative methods based upon residual averaging

Iterative methods for solving boundary value problems for systems of nonlinear partial differential equations are discussed. The methods involve subtracting an average of residuals from one approximation in order to arrive at a subsequent approximation. Two abstract methods in Hilbert space are given and application of these methods to quasilinear systems to give numerical schemes for such problems is demonstrated. Potential theoretic matters related to the iteration schemes are discussed.

Neuberger, J. W.↗

Wave-driven winds from cool stars. I - Some effects of magnetic field geometry

The wave-driven wind theory of Hartmann and MacGregor (1980) is extended to include effects due to non-radial divergence of the flow. Specifically, isothermal expansion within a flow tube whose cross-sectional area increases outward faster than the square of the radius near the stellar surface is considered. It is found that the qualitative conclusions of Hartmann and MacGregor concerning the physical properties of Alfven wave-driven winds are largely unaffected. In particular, mass fluxes of similar magnitude are obtained, and wave dissipation is still necessary to produce acceptably small terminal velocities. Increasingly divergent flow geometries generally lead to higher initial wind speeds and slightly lower terminal velocities. For some cases of extremely rapid flow tube divergence, steady supersonic wind solutions which extend to infinity with vanishing gas pressure cannot be obtained. In addition, departures from spherical symmetry can cause the relative Alfven wave amplitude delta-B/B to become approximately greater than 1 within several stellar radii of the base of the wind, suggesting that nonlinear processes may contribute to the wave dissipation required by the theory.

Hartmann, L.↗

Formulation and solution of rotary-wing aeroelastic stability and response problems

The state of the art in the formulation and solution of rotary-wing aeroelastic stability and response problems is reviewed in detail. The approximations used in the structural, inertia and aerodynamic operators are discussed. The important role of geometric nonlinearities, due to moderate deflections, and aerodynamic stall in the aeroelastic stability and response problem are identified. It is also shown that geometric nonlinearities are of primary importance in aeroelastic stability calculations, and have a more limited, though important, role in response calculations. Next, formulation of coupled rotor/fuselage problems is described, for both air and ground resonance type problems. Both topics, the isolated blade problem and the coupled rotor/fuselage problem, are treated for both hover and forward flight. Solution of aeroelastic stability and response problems proceeds in two stages. First, the spatial dependence is eliminated by using Galerkin's method, or by using the finite element method. Next the nonlinear, or linear, ordinary differential equation with periodic coefficients have to be solved for stability or response. Efficient numerical methods for accomplishing these objectives are presented in a comprehensive manner. The paper contains a number of illustrative numerical results which are intended to clarify various aspects of the modeling process and serve as representative results for both aeroelastic stability and response calculations for a variety of blade and rotor configurations.

Friedmann, P. P.↗

Statistical theories of Langmuir turbulence. II - Subsonic to sonic transition

The subsonic limit of the quadratic direct interaction approximation (DIA) applied to the Zakharov equations is compared with the cubic DIA applied to the nonlinear Schroedinger equation, which is the subsonic limit of the Zakharov equations. Comparisons with Monte Carlo simulations of a truncated system show that the first theory more accurately describes the regime of stationary turbulence, while the second theory more accurately describes the subsonic evolution of the modulational instability. The weak turbulence limits of the two theories describe the sonic and subsonic regimes, respectively. The addition of vertex corrections to the DIA leads to a hybrid weak turbulence theory that smoothly interpolates between the sonic and subsonic regimes.

Dubois, D. F.↗

Postbuckling of laminated anisotropic panels

A two-part study of the buckling and postbuckling of laminated anisotropic plates with bending-extensional coupling is presented. The first part involves the development and application of a modified Rayleigh-Ritz analysis technique. Modifications made to the classical technique can be grouped into three areas. First, known symmetries of anisotropic panels are exploited in the selection of approximation functions. Second, a reduced basis technique based on these same symmetries is applied in the linear range. Finally, geometric boundary conditions are enforced via an exterior penalty function approach, rather than relying on choice of approximation functions to satisfy these boundary conditions. Numerical results are presented for both the linear and nonlinear range, with additional studies made to determine the effect of variation in penalty parameter and number of basis vectors. In the second part, six panels possessing anisotropy and bending-extensional coupling are tested. Detailed comparisons are made between experiment and finite element results in order to gain insight into the postbuckling and failure characteristics of such panels. The panels are constructed using two different lamination sequences, and panels with three different aspect ratios were constructed for each lamination sequence.

Jeffrey, Glenda L.↗

Sensor-actuator placement for flexible structures with actuator dynamics

A novel approach for placement of sensors and actuators in control of flexible space structures is developed. Using an approximation of the control forces and output measurements by spatially continuous functions, the approach follows a nonlinear programming technique to determine optimal locations for sensors and actuators. Two different criteria are considered for the placement of sensors and actuators. The first criterion optimizes the location of the sensors and actuators in order to move the transmission zeros of the system farther to the left of the imaginary axis. The second criterion, however, places the sensors and actuators to optimize a function of the singular values of the Hankel matrix, which includes both measures of controllability and observability. Moreover, the effect of actuator dynamics in the placement of sensors and actuators is investigated.

Maghami, P. G.↗

Transformation of two and three-dimensional regions by elliptic systems

A reliable linear system is presented for grid generation in 2-D and 3-D. The method is robust in the sense that convergence is guaranteed but is not as reliable as other nonlinear elliptic methods in generating nonfolding grids. The construction of nonfolding grids depends on having reasonable approximations of cell aspect ratios and an appropriate distribution of grid points on the boundary of the region. Some guidelines are included on approximating the aspect ratios, but little help is offered on setting up the boundary grid other than to say that in 2-D the boundary correspondence should be close to that generated by a conformal mapping. It is assumed that the functions which control the grid distribution depend only on the computational variables and not on the physical variables. Whether this is actually the case depends on how the grid is constructed. In a dynamic adaptive procedure where the grid is constructed in the process of solving a fluid flow problem, the grid is usually updated at fixed iteration counts using the current value of the control function. Since the control function is not being updated during the iteration of the grid equations, the grid construction is a linear procedure. However, in the case of a static adaptive procedure where a trial solution is computed and used to construct an adaptive grid, the control functions may be recomputed at every step of the grid iteration.

Mastin, C. Wayne↗

Role of cross-field current instability in substorm onsets and intensifications

A cross field current instability is investigated as a potential mechanism for current reduction/disruption during substorm onsets and intensifications. Linear stability analysis shows that, for sufficiently strong current density, the instability can occur in the plasma sheet with a growth time comparable to the substorm onset time and excite waves with a significant electromagnetic component. Nonlinear analysis shows that the wave growth reaches the nonlinear stage in less than one ion gyroperiod and can reduce the cross tail current by approximately 15 to 28 of its initial value at saturation. The resulting anomalous resistivity is 11 to 12 orders of magnitude above the classical value. For a typical current reduction, the plasma in the disruption region is subjected to an earthward force. The substorm development scenario constructed based on this instability can readily account for a large number of substorm features.

Lui, Anthony T. Y.↗

Time Varying Compensator Design for Reconfigurable Structures Using Non-Collocated Feedback

Analysis and synthesis tools are developed to improved the dynamic performance of reconfigurable nonminimum phase, nonstrictly positive real-time variant systems. A novel Spline Varying Optimal (SVO) controller is developed for the kinematic nonlinear system. There are several advantages to using the SVO controller, in which the spline function approximates the system model, observer, and controller gain. They are: The spline function approximation is simply connected, thus the SVO controller is more continuous than traditional gain scheduled controllers when implemented on a time varying plant; ft is easier for real-time implementations in storage and computational effort; where system identification is required, the spline function requires fewer experiments, namely four experiments; and initial startup estimator transients are eliminated. The SVO compensator was evaluated on a high fidelity simulation of the Shuttle Remote Manipulator System. The SVO controller demonstrated significant improvement over the present arm performance: (1) Damping level was improved by a factor of 3; and (2) Peak joint torque was reduced by a factor of 2 following Shuttle thruster firings.

Scott, Michael A.↗

Nonlinear Ultrasonic Characterization of Fatigue Microstructures

Dislocation dipole substructures formed during metal fatigue are shown to produce a substantial distortion of ultrasonic waves propagating through the fatigued material. A model of ultrasonic wave-dislocation dipole interactions is developed that quantifies the wave distortion by means of a material nonlinearity parameter (beta). Application of the model to AA2024-T4 predicts a value of p approximately 300% larger in material cyclically loaded for 100 kcycles in stress-control at 276 MPa and R=0 than that measured for virgin material. Experimental measurements show a monotonic increase in p as a function of the number of fatigue cycles that closely approaches the predicted increase. The experiments also suggest that the relevant dislocation substructures are localized in the material.

Cantrell, John H.↗

A Learn-to-Fly Approach for Adaptively Tuning Flight Control Systems

A method is presented for adaptively tuning feedback control gains in a flight control system to achieve desired closed-loop performance. The method combines efficient aircraft parameter estimation, for identifying closed-loop dynamics models, with online nonlinear optimization, for sequentially perturbing and updating control gains to improve performance. Access to flight measurements and the control gains is required, but no prior information about the aircraft dynamics or the flight control architecture is needed. After the procedure, the optimized control gains (with uncertainties), the open-loop dynamics model, and the closed-loop dynamics model are available. The method is demonstrated for tuning a longitudinal stability augmentation system using a realistic nonlinear flight dynamics simulation of a subscale airplane. Convergence was attained using five maneuvers that spanned approximately one minute of flight test time. Although shown for a relatively simple case, the method is general and can be applied to other aircraft, axes of motion, performance metrics, and control system designs.

Adaptive tuning↗

Hybrid perturbation/Bubnov-Galerkin technique for nonlinear thermal analysis

A two step hybrid analysis technique to predict the nonlinear steady state temperature distribution in structures and solids is presented. The technique is based on the regular perturbation expansion and the classical Bubnov-Galerkin approximation. The functions are obtained by using the regular perturbation method. These functions are selected as coordinate functions and the classical Bubnov-Galerkin technique is used to compute their amplitudes. The potential of the proposed hybrid technique for the solution of nonlinear thermal problems is discussed. The effectiveness of this technique is demonstrated by the effects of conduction, convection, and radiation modes of heat transfer. It is indicated that the hybrid technique overcomes the two major drawbacks of the classical techniques: (1) the requirement of using a small parameter in the regular perturbation method; and (2) the arbitrariness in the choice of the coordinate functions in the Bubnov-Galerkin technique. The proposed technique extends the range of applicability of the regular perturbation method and enhances the effectiveness of the Bubnov-Galerkin technique.

Noor, A. K.↗

A Fast Temporal Decomposition Procedure for Long-Horizon Nonlinear Dynamic Programming

We propose a fast temporal decomposition procedure for solving long-horizon nonlinear dynamic programs. The core of the procedure is sequential quadratic programming (SQP) that utilizes a differentiable exact augmented Lagrangian as the merit function. Within each SQP iteration, we approximately solve the Newton system using an overlapping temporal decomposition strategy. We show that the approximate search direction is still a descent direction of the augmented Lagrangian provided the overlap size and penalty parameters are suitably chosen, which allows us to establish the global convergence. Moreover, we show that a unit step size is accepted locally for the approximate search direction and further establish a uniform, local linear convergence over stages. This local convergence rate matches the rate of the recent Schwarz scheme (Na et al. 2022). However, the Schwarz scheme has to solve nonlinear subproblems to optimality in each iteration, whereas we only perform a single Newton step instead. Numerical experiments validate our theories and demonstrate the superiority of our method.

97 MATHEMATICS AND COMPUTING↗

Reconstructing the primordial spectrum of fluctuations of the universe from the observed nonlinear clustering of galaxies

It was discovered that the nonlinear evolution of the two point correlation function in N-body experiments of galaxy clustering with Omega = 1 appears to be described to good approximation by a simple general formula. The underlying form of the formula is physically motivated, but its detailed representation is obtained empirically by fitting to N-body experiments. In this paper, the formula is presented along with an inverse formula which converts a final, nonlinear correlation function into the initial linear correlation function. The inverse formula is applied to observational data from the CfA, IRAs, and APM galaxy surveys, and the initial spectrum of fluctuations of the universe, if Omega = 1.

Hamilton, A. J. S.↗