A Chebyshev interval method for nonlinear dynamic systems under uncertainty

This paper proposes a new interval analysis method for the dynamic response of nonlinear systems with uncertain-but-bounded parameters using Chebyshev polynomial series. Interval model can be used to describe nonlinear dynamic systems under uncertainty with low-order Taylor series expansions. However, the Taylor series-based interval method can only suit problems with small uncertain levels. To account for larger uncertain levels, this study introduces Chebyshev series expansions into interval model to develop a new uncertain method for dynamic nonlinear systems. In contrast to the Taylor series, the Chebyshev series can offer a higher numerical accuracy in the approximation of solutions. The Chebyshev inclusion function is developed to control the overestimation in interval computations, based on the truncated Chevbyshev series expansion. The Mehler integral is used to calculate the coefficients of Chebyshev polynomials. With the proposed Chebyshev approximation, the set of ordinary differential equations (ODEs) with interval parameters can be transformed to a new set of ODEs with deterministic parameters, to which many numerical solvers for ODEs can be directly applied. Two numerical examples are applied to demonstrate the effectiveness of the proposed method, in particular its ability to effectively control the overestimation as a non-intrusive method.     

Breaking the limits: The Taylor series method

This paper discusses several examples of ordinary differential equation (ODE) applications that are difficult to solve numerically using conventional techniques, but which can be solved successfully using the Taylor series method. These results are hard to obtain using other methods such as Runge-Kutta or similar schemes; indeed, in some cases these other schemes are not able to solve such systems at all. In particular, we explore the use of the high-precision arithmetic in the Taylor series method for numerically integrating ODEs. We show how to compute the partial derivatives, how to propagate sets of initial conditions, and, finally, how to achieve the Brouwer’s Law limit in the propagation of errors in long-time simulations. The TIDES software that we use for this work is freely available from a website.     

On the construction of high order <Emphasis Type="Italic">A</Emphasis>(<Emphasis Type="Italic">α</Emphasis>)-stable hybrid linear multistep methods for stiff IVPs in ODEs

In this paper, we present a class of A(α)-stable hybrid linear multistep methods for numerical solving stiff initial value problems (IVPs) in ordinary differential equations (ODEs). The method considered uses a second derivative like the Enright’s second derivative linear multistep methods for stiff IVPs in ODEs.     

A New Integrator for Special Third Order Differential Equations With Application to Thin Film Flow Problem

In recent time, Runge-Kutta methods that integrate special third order ordinary differential equations (ODEs) directly are proposed to address efficiency issues associated with classical Runge-Kutta methods. Albeit, the methods require evaluation of three set of equations to proceed with the numerical integration. In this paper, we propose a class of multistep-like Runge-Kutta methods (hybrid methods), which integrates special third order ODEs directly. The method is completely derivative-free. Algebraic order conditions of the method are derived. Using the order conditions, a four-stage method is presented. Numerical experiment is conducted on some test problems. The method is also applied to a practical problem in Physics and engineering to ascertain its validity. Results from the experiment show that the new method is more accurate and efficient than the classical Runge-Kutta methods and a class of direct Runge-Kutta methods recently designed for special third order ODEs.     

Numerical solutions of index-1 differential algebraic equations can be computed in polynomial time

The cost of solving an initial value problem for index-1 differential algebraic equations to accuracy ɛ is polynomial in ln(1/ɛ). This cost is obtained for an algorithm based on the Taylor series method for solving differential algebraic equations developed by Pryce. This result extends a recent result by Corless for solutions of ordinary differential equations. The results of the standard theory of information-based complexity give exponential cost for solving ordinary differential equations, being based on a different model.     

Reduction of hugoniot-maslov chains for trajectories of solitary vortices of the “shallow water” equations to the hill equation

According to Maslov’s idea, many two-dimensional, quasilinear hyperbolic systems of partial differential equations admit only three types of singularities that are in general position and have the property of “structure self-similarity and stability.” Those are: shock waves, “narrow” solitons, and “square-root” point singularities (solitary vortices). Their propagation is described by an infinite chain of ordinary differential equations (ODE) that generalize the well-known Hugoniot conditions for shock waves. After some reasonable closure of the chain for the case of solitary vortices in the “shallow water” equations, we obtain a nonlinear system of sixteen ODE, which is exactly equivalent to the (linear) Hill equation with a periodic potential. This means that, in some approximations, the trajectory of a solitary vortex can be described by the Hill equation. This result can be used to predict the trajectory of the vortex center if we know its observable part. Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 112, No. 1, pp. 47–66.     

Approximated solutions in rational form for systems of differential equations

In this paper we present a technique to study the existence of rational solutions for systems of differential equations — for an ordinary differential equation, in particular. The method is relatively straightforward; it is based on a rationality characterisation that involves matrix Padé approximants. It is important to note that, when the solution is rational, we use formal power series “without taking into account” their circle of convergence; at the end of this paper we justify this. We expound the theory for systems of linear first-order ordinary differential equations in the general case. However, the main ideas are applied in numerical resolution of partial differential equations. This revised version was published online in June 2006 with corrections to the Cover Date.     

Computational bounds on polynomial differential equations

In this paper we study from a computational perspective some properties of the solutions of polynomial ordinary differential equations.We consider elementary (in the sense of Analysis) discrete-time dynamical systems satisfying certain criteria of robustness. We show that those systems can be simulated with elementary and robust continuous-time dynamical systems which can be expanded into fully polynomial ordinary differential equations in Q[π]. This sets a computational lower bound on polynomial ODEs since the former class is large enough to include the dynamics of arbitrary Turing machines.We also apply the previous methods to show that the problem of determining whether the maximal interval of definition of an initial-value problem defined with polynomial ODEs is bounded or not is in general undecidable, even if the parameters of the system are computable and comparable and if the degree of the corresponding polynomial is at most 56.Combined with earlier results on the computability of solutions of polynomial ODEs, one can conclude that there is from a computational point of view a close connection between these systems and Turing machines.     

Averaging and hopf bifurcation in degenerate systems of differential equations with delays’

This article considers bifurcation of families of periodic solutions from equilibrium states for systems of differential equations with constant delays that, upon linearization at zero, become systems of ordinary differential equations. An averaging transformation is constructed to simplify solution of the bifurcation problem. There is an error in the title of Belan’s article in Dinamicheskie Sistemy, No. 11. 1992. It should read “On quasiperiodic solutions of semilinear parabolic equations”. Translated fromDinamicheskie Sistemy. Vol. 12, pp. 85–91, 1993.     

On Galilean connections and the first jet bundle

We see how the first jet bundle of curves into affine space can be realized as a homogeneous space of the Galilean group. Cartan connections with this model are precisely the geometric structure of second-order ordinary differential equations under time-preserving transformations — sometimes called KCC-theory. With certain regularity conditions, we show that any such Cartan connection induces “laboratory” coordinate systems, and the geodesic equations in this coordinates form a system of second-order ordinary differential equations. We then show the converse — the “fundamental theorem” — that given such a coordinate system, and a system of second order ordinary differential equations, there exists regular Cartan connections yielding these, and such connections are completely determined by their torsion.     

Half-explicit Runge-Kutta methods with explicit stages for differential-algebraic systems of index 2

Usually the straightforward generalization of explicit Runge-Kutta methods for ordinary differential equations to half-explicit methods for differential-algebraic systems of index 2 results in methods of orderq≤2. The construction of higher order methods is simplified substantially by a slight modification of the method combined with an improved strategy for the computation of the algebraic solution components. We give order conditions up to orderq=5 and study the convergence of these methods. Based on the fifth order method of Dormand and Prince the fifth order half-explicit Runge-Kutta method HEDOP5 is constructed that requires the solution of 6 systems of nonlinear equations per step of integration.     

On the nonresonance property of linear differential-algebraic systems

We consider a system of linear ordinary differential equations in which the coefficient matrix multiplying the derivative of the unknown vector function is identically singular. For systems with constant and variable coefficients, we obtain nonresonance criteria (criteria for bounded-input bounded-output stability). For single-input control systems, we consider the problem of synthesizing a nonresonant system in the stationary and nonstationary cases. An arbitrarily high unsolvability index is admitted. The analysis is carried out under assumptions providing the existence of a so-called “equivalent form” with separated “algebraic” and “differential” components.     

Piecewise finite series solution of nonlinear initial value differential problem

In this paper, we apply a piecewise finite series as a hybrid analytical-numerical technique for solving some nonlinear systems of ordinary differential equations. The finite series is generated by using the Adomian decomposition method, which is an analytical method that gives the solution based on a power series and has been successfully used in a wide range of problems in applied mathematics. We study the influence of the step size and the truncation order of the piecewise finite series Adomian (PFSA) method on the accuracy of the solutions when applied to nonlinear ODEs. Numerical comparisons between the PFSA method with different time steps and truncation orders against Runge-Kutta type methods are presented. Based on the numerical results we propose a low value truncation order approach with small time step size. The numerical results show that the PFSA method is accurate and easy to implement with the proposed approach.     

Equations for the Missing Boundary Values in the Hamiltonian Formulation of Optimal Control Problems

Partial differential equations for the unknown final state and initial costate arising in the Hamiltonian formulation of regular optimal control problems with a quadratic final penalty are found. It is shown that the missing boundary conditions for Hamilton’s canonical ordinary differential equations satisfy a system of first-order quasilinear vector partial differential equations (PDEs), when the functional dependence of the H-optimal control in phase-space variables is explicitly known. Their solutions are computed in the context of nonlinear systems with ℝ ⁿ -valued states. No special restrictions are imposed on the form of the Lagrangian cost term. Having calculated the initial values of the costates, the optimal control can then be constructed from on-line integration of the corresponding 2n-dimensional Hamilton ordinary differential equations (ODEs). The off-line procedure requires finding two auxiliary n×n matrices that generalize those appearing in the solution of the differential Riccati equation (DRE) associated with the linear-quadratic regulator (LQR) problem. In all equations, the independent variables are the finite time-horizon duration T and the final-penalty matrix coefficient S, so their solutions give information on a whole two-parameter family of control problems, which can be used for design purposes. The mathematical treatment takes advantage from the symplectic structure of the Hamiltonian formalism, which allows one to reformulate Bellman’s conjectures concerning the “invariant-embedding” methodology for two-point boundary-value problems. Results for LQR problems are tested against solutions of the associated differential Riccati equation, and the attributes of the two approaches are illustrated and discussed. Also, nonlinear problems are numerically solved and compared against those obtained by using shooting techniques.     

B-Theory of Runge-Kutta methods for stiff volterra functional differential equations

B-stability andB-convergence theories of Runge-Kutta methods for nonlinear stiff Volterra functional differential equations (VFDEs) are established which provide unified theoretical foundation for the study of Runge-Kutta methods when applied to nonlinear stiff initial value problems (IVPs) in ordinary differential equations (ODEs), delay differential equations (DDEs), integro-differential equations (IDEs) and VFDEs of other type which appear in practice.     

B-Theory of Runge-Kutta methods for stiff Volterra functional differential equations

B-stability and B-convergence theories of Runge-Kutta methods for nonlinear stiff Volterra func-tional differential equations(VFDEs)are established which provide unified theoretical foundation for the studyof Runge-Kutta methods when applied to nonlinear stiff initial value problems(IVPs)in ordinary differentialequations(ODEs),delay differential equations(DDEs),integro-differential equatioons(IDEs)and VFDEs of     

An algebraic foundation for factoring linear boundary problems

Motivated by boundary problems for linear differential equations, we define an abstract boundary problem as a pair consisting of a surjective linear map (“differential operator”) and an orthogonally closed subspace of the dual space (“boundary conditions”). Defining the composition of boundary problems corresponding to their Green’s operators in reverse order, we characterize and construct all factorizations of a boundary problem from a given factorization of the defining operator. For the case of ordinary differential equations, the main results can be made algorithmic. We conclude with a factorization of a boundary problem for the wave equation. This work was supported by the Austrian Science Fund (FWF) under the SFB grant F1322.     

Asymptotic behavior of multistep runge-kutta methods for systems of delay differential equations

1. IntroductionIn order to assess the asymptotic behavior of numerical methods for DDEs, much attention has been given in the literature to the scalar case (cL [1-6]). UP to now) only partialresults (of. [7-10]) have dealt with the delay systemswhere y(t) = (yi(t), so(t),' ) yp(t))" E Cd, which is unknown for t > 0, L and M areconstat complex p x Hmatrices, T > 0 is a constat delay and W(t) 6 CP is a specifiedinitial function.In [111, C.J. Zhang and S.Z. Zhou made an investigation on…     

Dissipativity of Runge-Kutta methods for Volterra functional differential equations

This paper is concerned with the numerical dissipativity of nonlinear Volterra functional differential equations (VFDEs). We give some dissipativity results of Runge-Kutta methods when they are applied to VFDEs. These results provide unified theoretical foundation for the numerical dissipativity analysis of systems in ordinary differential equations (ODEs), delay differential equations (DDEs), integro-differential equations (IDEs), Volterra delay integro-differential equations (VDIDEs) and VFDEs of other type which appear in practice. Numerical examples are given to confirm our theoretical results.