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.
Connecting Functions of Nonlinear Ordinary Differential Equations.
Abstract not provided.
Steyer, Andrew↗
Neural Ordinary Differential Equations for Nonlinear System Identification
Neural ordinary differential equations (NODE) have been recently proposed as a promising approach for nonlinear system identification tasks. In this work, we systematically compare their predictive performance with current state-of-the-art nonlinear and classical linear methods. In particular, we present a quantitative study comparing NODE's performance against neural state-space models and classical linear system identification methods and evaluate their inference speed and prediction performance on open-loop errors across eight different dynamical systems. The experiments show that NODEs can consistently improve the prediction accuracy by order of magnitude compared to benchmark methods. Besides improved accuracy, we also observed that NODEs are less sensitive to hyperparameters compared to neural state-space models by paying the cost of increased computation at the inference time.
machine leaning, system identification, physics in↗
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↗
PETN spark-gap detonators
The well-developed theory of Lorentz plasma that is dominated by electron–ion interactions is used to calculate the PETN arc characteristics. The spark-gap discharge current is a ramp with 10 to 25 ns rise time to peak and remaining constant subsequently. The approximate formulas for the arc channel conductivity, arc temperature, arc radius, and shock pressure from the arc are obtained from a system of nonlinear ordinary differential equations, which is the similarity solution of hydrodynamic equations similar to the Braginskii approximation. These arc parameters are given for the peak current ranging from 100 A to 1000 A and with different rise times. Representative cases are compared to the nonlinear ordinary differential equation code results. The shock pressures at the peak current are comparable to those from a typical commercial EBW bridgewire burst reported in the literature; the arc radius at the peak current is comparable to a typical bridgewire diameter of 0.0375 mm (e.g., RISI detonators, RP-1, and RP-80). The relevant Pop-Plot for low-density PETN is converted into an empirical detonation criterion, which is applicable to explosives subject to shocks of variable pressure. Finally, this criterion is then used to determine the detonation thresholds, which are comparable with test data obtained by Tucker, et al.
42 ENGINEERING↗
Integrable nonlinear evolution of the qubit
Highlights: •We study the nonlinear generalization of the von Neumann equation. •We discuss the nonlinear von Neumann equation for the qubit states. •We analyze solutions to the Riccati system corresponding to evolution of a qubit. -- Abstract: The nonlinear generalization of the von Neumann equation is studied. The global version of this equation preserves convexity of the space of states. The particular case of the evolution of pure states refers to the nonlinear Schrödinger equation discovered by Gisin. The concrete realization of the nonlinear von Neumann equation is investigated in detail for the density matrix representing the qubit states. Such equation reduces to the integrable classical Riccati system of nonlinear ordinary differential equations. An interesting property of the nonlinear dynamics described by this system is the global asymptotical stability of stationary solutions related to evolution from mixed states to pure states.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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.↗
Improved quantum algorithms for linear and nonlinear differential equations
We present substantially generalized and improved quantum algorithms over prior work for inhomogeneous linear and nonlinear ordinary differential equations (ODE). Specifically, we show how the norm of the matrix exponential characterizes the run time of quantum algorithms for linear ODEs opening the door to an application to a wider class of linear and nonlinear ODEs. In [1], a quantum algorithm for a certain class of linear ODEs is given, where the matrix involved needs to be diagonalizable. The quantum algorithm for linear ODEs presented here extends to many classes of non-diagonalizable matrices including singular matrices. The algorithm here is also exponentially faster than the bounds derived in [1] for certain classes of diagonalizable matrices. Our linear ODE algorithm is then applied to nonlinear differential equations using Carleman linearization (an approach taken recently by us in [2]). The improvement over that result is two-fold. First, we obtain an exponentially better dependence on error. This kind of logarithmic dependence on error has also been achieved by [3], but only for homogeneous nonlinear equations. Second, the present algorithm can handle any sparse matrix (that models dissipation) if it has a negative log-norm (including non-diagonalizable matrices), whereas [2] and [3] additionally require normality.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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.↗