Chebyshev series solution of nonlinear ordinary differential equations - Initial-value problems
Chebyshev series solution of nonlinear ordinary differential equations - initial value problems
Byrd, P. F.↗
SEARCH · Engineering Papers
Results for “nonlinear ordinary differential 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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Asymptotic integration algorithms for nonhomogeneous, nonlinear, first order, ordinary differential equations
New methods for integrating systems of stiff, nonlinear, first order, ordinary differential equations are developed by casting the differential equations into integral form. Nonlinear recursive relations are obtained that allow the solution to a system of equations at time t plus delta t to be obtained in terms of the solution at time t in explicit and implicit forms. Examples of accuracy obtained with the new technique are given by considering systems of nonlinear, first order equations which arise in the study of unified models of viscoplastic behaviors, the spread of the AIDS virus, and predator-prey populations. In general, the new implicit algorithm is unconditionally stable, and has a Jacobian of smaller dimension than that which is acquired by current implicit methods, such as the Euler backward difference algorithm; yet, it gives superior accuracy. The asymptotic explicit and implicit algorithms are suitable for solutions that are of the growing and decaying exponential kinds, respectively, whilst the implicit Euler-Maclaurin algorithm is superior when the solution oscillates, i.e., when there are regions in which both growing and decaying exponential solutions exist.
Walker, K. P.↗
An indirect optimization method with improved convergence characteristics
A method is presented for obtaining numerical solutions to a certain class of two-point boundary value problems which arise in optimal control theory. The problems are characterized by systems of nonlinear ordinary differential equations with nonlinear boundary conditions.
Doiron, H. H.↗
A general conclusion regarding the large amplitude flexural vibration of beams and plates.
Use is made of the two coupled partial differential equations that govern the large amplitude flexural vibration of flat plates. On the basis of an assumed vibration mode taken to be separable in the space coordinates and time, the modal equation is obtained. In the case of a one-term solution, it is shown that the modal equation is the same as the dynamic equation of a mass-spring combination where the restoring force of the nonlinear spring is a cubic and an odd function of the displacement. Since the modal equation is a nonlinear, ordinary differential equation with positive coefficients, it is shown that the nonlinearity associated with the large amplitude flexural vibration of beams and plates is of the hardening type, that is, the frequency increases with amplitude.
Pandalai, K. A. V.↗
On the construction of highly stable, explicit, numerical methods for integrating coupled ordinary differential equations with parasitic eigenvalues
Explicit numerical methods for integrating coupled nonlinear ordinary differential equations with parasitic eignevalues
Lomax, H.↗
Light Diffraction by Large Amplitude Ultrasonic Waves in Liquids
Light diffraction from ultrasound, which can be used to investigate nonlinear acoustic phenomena in liquids, is reported for wave amplitudes larger than that typically reported in the literature. Large amplitude waves result in waveform distortion due to the nonlinearity of the medium that generates harmonics and produces asymmetries in the light diffraction pattern. For standing waves with amplitudes above a threshold value, subharmonics are generated in addition to the harmonics and produce additional diffraction orders of the incident light. With increasing drive amplitude above the threshold a cascade of period-doubling subharmonics are generated, terminating in a region characterized by a random, incoherent (chaotic) diffraction pattern. To explain the experimental results a toy model is introduced, which is derived from traveling wave solutions of the nonlinear wave equation corresponding to the fundamental and second harmonic standing waves. The toy model reduces the nonlinear partial differential equation to a mathematically more tractable nonlinear ordinary differential equation. The model predicts the experimentally observed cascade of period-doubling subharmonics terminating in chaos that occurs with increasing drive amplitudes above the threshold value. The calculated threshold amplitude is consistent with the value estimated from the experimental data.
Adler, Laszlo↗
A recursive approach to the equations of motion for the maneuvering and control of flexible multi-body systems
Interest lies in a mathematical formulation capable of accommodating the problem of maneuvering a space structure consisting of a chain of articulated flexible substructures. Simultaneously, any perturbations from the 'rigid body' maneuvering and any elastic vibration must be suppressed. The equations of motion for flexible bodies undergoing rigid body motions and elastic vibrations can be obtained conveniently by means of Lagrange's equations in terms of quasi-coordinates. The advantage of this approach is that it yields equations in terms of body axes, which are the same axes that are used to express the control forces and torques. The equations of motion are nonlinear hybrid differential quations. The partial differential equations can be discretized (in space) by means of the finite element method or the classical Rayleigh-Ritz method. The result is a set of nonlinear ordinary differential equations of high order. The nonlinearity can be traced to the rigid body motions and the high order to the elastic vibration. Elastic motions tend to be small when compared with rigid body motions.
Kwak, Moon K.↗
Postbuckling of orthotropic composite plates loaded in compression
The nonlinear large deflection equations of von Karman are written for 'specially' orthotropic plates. The equations are then manipulated to determine the parameters required to establish postbuckling behavior. It is found that only two new parameters are needed beyond those required for buckling. By assuming trigonometric functions in one direction, the plate equations are converted into ordinary nonlinear differential equations which are solved numerically using a two point boundary problem solver that makes use of Newton's method. The postbuckling behavior is obtained for simply supported and clamped, long, rectangular, orthotropic plates covering the complete range of dimensions and material properties.
Stein, M.↗
Applying integrals of motion to the numerical solution of differential equations
A method is developed for using the integrals of systems of nonlinear, ordinary differential equations in a numerical integration process to control the local errors in these integrals and reduce the global errors of the solution. The method is general and can be applied to either scaler or vector integrals. A number of example problems, with accompanying numerical results, are used to verify the analysis and support the conjecture of global error reduction.
Jezewski, D. J.↗
Applying integrals of motion to the numerical solution of differential equations
A method is developed for using the integrals of systems of nonlinear, ordinary, differential equations in a numerical integration process to control the local errors in these integrals and reduce the global errors of the solution. The method is general and can be applied to either scalar or vector integrals. A number of example problems, with accompanying numerical results, are used to verify the analysis and support the conjecture of global error reduction.
Vezewski, D. J.↗
Analytical results for postbuckling behavior of plates in compression and shear
The postbuckling behavior of long rectangular isotropic and orthotropic plates is determined. By assuming trigonometric functions in one direction, the nonlinear partial differential equations of von Karman large deflection plate theory are converted into nonlinear ordinary differential equations. The ordinary differential equations are solved numerically using an available boundary value problem solver which makes use of Newton's method. Results for longitudinal compression show different postbuckling behavior between isotropic and orthotropic plates. Results for shear show that change in inplane edge constraints can cause large change in postbuckling stiffness.
Stein, M.↗
Analytical results for post-buckling behaviour of plates in compression and in shear
The postbuckling behavior of long rectangular isotropic and orthotropic plates is determined. By assuming trigonometric functions in one direction, the nonlinear partial differential equations of von Karman large deflection plate theory are converted into nonlinear ordinary differential equations. The ordinary differential equations are solved numerically using an available boundary value problem solver which makes use of Newton's method. Results for longitudinal compression show different postbuckling behavior between isotropic and orthotropic plates. Results for shear show that change in inplane edge constraints can cause large change in postbuckling stiffness.
Stein, M.↗
Nonlinear evolution of magnetic flux ropes. 2: Finite beta plasma
In this second paper on the evolution of magnetic flux ropes we study the effects of gas pressure. We assume that the energy transport is described by a polytropic relationship and reduce the set of ideal MHD equations to a single, second-order, nonlinear, ordinary differential equation for the evolution function. For this conservative system we obtain a first integral of motion. To analyze the possible motions, we use a mechanical analogue -- a one-dimensional, nonlinear oscillator. We find that the effective potential for such an oscillator depends on two parameters: the polytropic index gamma and a dimensionless quantity kappa the latter being a function of the plasma beta, the strength of the azimuthal magnetic field relative to the axial field of the flux rope, and gamma. Through a study of this effective potential we classify all possible modes of evolution of the system. In the main body of the paper, we focus on magnetic flux ropes whose field and gas pressure increase steadily towards the symmetry axis. In this case, for gamma greater than 1 and all values of kappa, only oscillations are possible. For gamma less than 1, however, both oscillations and expansion are allowed. For gamma less than 1 and kappa below a critical value, the energy of the nonlinear oscillator determines whether the flux rope will oscillate or expand to infinity. For gamma less than 1 and kappa above critical, however, only expansion occurs. Thus by increasing kappa while keeping gamma fixed (less than 1), a phase transition occurs at kappa = kappa(sub critical) and the oscillatory mode disappears. We illustrate the above theoretical considerations by the example of a flux rope of constant field line twist evolving self-similarly. For this example, we present the full numerical MHD solution. In an appendix to the paper we catalogue all possible evolutions when (1) either the magnetic field or (2) the gas pressure decreases monotonically toward the axis. We find that in these cases critical conditions can occur for gamma greater than 1. While in most cases the flux rope collapses, there are notable exceptions when, for certain ranges of kappa and gamma, collapse may be averted.
Osherovich, V. A.↗
Nonlinear Equations of Motion for Cantilever Rotor Blades in Hover with Pitch Link Flexibility, Twist, Precone, Droop, Sweep, Torque Offset, and Blade Root Offset
Nonlinear equations of motion for a cantilever rotor blade are derived for the hovering flight condition. The blade is assumed to have twist, precone, droop, sweep, torque offset and blade root offset, and the elastic axis and the axes of center of mass, tension, and aerodynamic center coincident at the quarter chord. The blade is cantilevered in bending, but has a torsional root spring to simulate pitch link flexibility. Aerodynamic forces acting on the blade are derived from strip theory based on quasi-steady two-dimensional airfoil theory. The equations are hybrid, consisting of one integro-differential equation for root torsion and three integro-partial differential equations for flatwise and chordwise bending and elastic torsion. The equations are specialized for a uniform blade and reduced to nonlinear ordinary differential equations by Galerkin's method. They are linearized for small perturbation motions about the equilibrium operating condition. Modal analysis leads to formulation of a standard eigenvalue problem where the elements of the stability matrix depend on the solution of the equilibrium equations. Two different forms of the root torsion equation are derived that yield virtually identical numerical results. This provides a reasonable check for the accuracy of the equations.
Hodges, D. H.↗
Nonlinear stability and control study of highly maneuverable high performance aircraft, phase 2
This research should lead to the development of new nonlinear methodologies for the adaptive control and stability analysis of high angle-of-attack aircraft such as the F18 (HARV). The emphasis has been on nonlinear adaptive control, but associated model development, system identification, stability analysis and simulation is performed in some detail as well. Various models under investigation for different purposes are summarized in tabular form. Models and simulation for the longitudinal dynamics have been developed for all types except the nonlinear ordinary differential equation model. Briefly, studies completed indicate that nonlinear adaptive control can outperform linear adaptive control for rapid maneuvers with large changes in alpha. The transient responses are compared where the desired alpha varies from 5 degrees to 60 degrees to 30 degrees and back to 5 degrees in all about 16 sec. Here, the horizontal stabilator is the only control used with an assumed first-order linear actuator with a 1/30 sec time constant.
Mohler, R. R.↗
Dynamic Nonlinear Elastic Stability of Helicopter Rotor Blades in Hover and in Forward Flight
Equations for large coupled flap-lag motion of hingeless elastic helicopter blades are consistently derived. Only torsionally-rigid blades excited by quasi-steady aerodynamic loads are considered. The nonlinear equations of motion in the time and space variables are reduced to a system of coupled nonlinear ordinary differential equations with periodic coefficients, using Galerkin's method for the space variables. The nonlinearities present in the equations are those arising from the inclusion of moderately large deflections in the inertia and aerodynamic loading terms. The resulting system of nonlinear equations has been solved, using an asymptotic expansion procedure in multiple time scales. The stability boundaries, amplitudes of nonlinear response, and conditions for existence of limit cycles are obtained analytically. Thus, the different roles played by the forcing function, parametric excitation, and nonlinear coupling in affecting the solution can be easily identified, and the basic physical mechanism of coupled flap-lag response becomes clear. The effect of forward flight is obtained with the requirement of trimmed flight at fixed values of the thrust coefficient.
Friedmann, P.↗
On numerical nonlinear analysis of highly flexible spinning cantilevers
The general nonlinear discretized equations of motion of spinning elastic solids and structures are derived as a set of nonlinear ordinary differential equations for the case when the strain-displacement and velocity-displacement relations are nonlinear up to the second order. It is shown that the cost of generation of such equations is proportional to the fourth power of the number of degrees of freedom. A computer program is written to automatically generate the equations for the case of spinning cantilevers with initial imperfections. The types and the number of the coordinate functions used in the trial solution are parameters of the program.
Utku, S.↗
Initial values for the integration scheme to compute the eigenvalues for propagation in ducts
A scheme for the calculation of eigenvalues in the problem of acoustic propagation in a two-dimensional duct is described. The computation method involves changing the coupled transcendental nonlinear algebraic equations into an initial value problem involving a nonlinear ordinary differential equation. The simplest approach is to use as initial values the hardwall eigenvalues and to integrate away from these values as the admittance varies from zero to its actual value with a linear variation. The approach leads to a powerful root finding routine capable of computing the transverse and axial wave numbers for two-dimensional ducts for any frequency, lining, admittance and Mach number without requiring initial guesses or starting points.
Eversman, W.↗