A high-order, purely frequency based harmonic balance formulation for continuation of periodic solutions: The case of non-polynomial nonlinearities

International audienceIn this paper, we extend the method proposed by Cochelin and Vergez [A high order purely frequency-based harmonic balance formulation for continuation of periodic solutions, Journal of Sound and Vibration, 324 (2009) 243-262] to the case of non-polynomial nonlinearities. This extension allows for the computation of branches of periodic solutions of a broader class of nonlinear dynamical systems. The principle remains to transform the original ODE system into an extended polynomial quadratic system for an easy application of the harmonic balance method (HBM). The transformation of non-polynomial terms is based on the differentiation of state variables with respect to the time variable, shifting the nonlinear non-polynomial nonlinearity to a time-independent initial condition equation, not concerned with the HBM. The continuation of the resulting algebraic system is here performed by the asymptotic numerical method (high order Taylor series representation of the solution branch) using a further differentiation of the non-polynomial algebraic equation with respect to the path parameter. A one dof vibro-impact system is used to illustrate how an exponential nonlinearity is handled, showing that the method works at very high order, 1000 in that case. Various kinds of nonlinear functions are also treated, and finally the nonlinear free pendulum is addressed, showing that very accurate periodic solutions can be computed with the proposed method

Computation Sequences for Series and Polynomials

Approximation to the solutions of non-linear differential systems is very useful when the exact solutions are unattainable. Perturbation expansion replaces the system with a sequences of smaller problems, only the first of which is typically nonlinear. This works well by hand for the first few terms, but higher order computations are typically too demanding for all but the most persistent. Symbolic computation is thus attractive; however, symbolic computation of the expansions almost always encounters intermediate expression swell, by which we mean exponential growth in subexpression size or repetitions. A successful management of spatial complexity is vital to compute meaningful results. This thesis contains two parts. In the first part, we investigate a heat transfer problem where two-dimensional buoyancy-induced flow between two concentric cylinders is studied. Series expansion with respect to Rayleigh number is used to compute an approximation of a solution, using a symbolic- numerical algorithm. Computation sequences are used to help reduce the size of intermediate expressions. Up to 30th order solutions are computed. Accuracy, validity and stability of the computed series solution are studied. In the second part, Hilbert’s 16th problem is investigated to find the maximum number of limit cycles of certain systems. Focus values of the systems are computed using perturbation theory, which form multivariate polynomial sys- tems. The real roots of such systems leads to possible limit cycle conditions. A modular regular chains approach is used to triangularize the polynomial systems and help to compute the real roots. A system with 9 limit cycles is constructed using the computed real roots

Order-of-magnitude differences in computational performance of analog Ising machines induced by the choice of nonlinearity

Ising machines based on nonlinear analog systems are a promising method to accelerate computation of NP-hard optimization problems. Yet, their analog nature is also causing amplitude inhomogeneity which can deteriorate the ability to find optimal solutions. Here, we investigate how the system's nonlinear transfer function can mitigate amplitude inhomogeneity and improve computational performance. By simulating Ising machines with polynomial, periodic, sigmoid and clipped transfer functions and benchmarking them with MaxCut optimization problems, we find the choice of transfer function to have a significant influence on the calculation time and solution quality. For periodic, sigmoid and clipped transfer functions, we report order-of-magnitude improvements in the time-to-solution compared to conventional polynomial models, which we link to the suppression of amplitude inhomogeneity induced by saturation of the transfer function. This provides insights into the suitability of systems for building Ising machines and presents an efficient way for overcoming performance limitations

Symbolic computation of exact solutions expressible in hyperbolic and elliptic functions for nonlinear PDEs

Algorithms are presented for the tanh- and sech-methods, which lead to closed-form solutions of nonlinear ordinary and partial differential equations (ODEs and PDEs). New algorithms are given to find exact polynomial solutions of ODEs and PDEs in terms of Jacobi's elliptic functions. For systems with parameters, the algorithms determine the conditions on the parameters so that the differential equations admit polynomial solutions in tanh, sech, combinations thereof, Jacobi's sn or cn functions. Examples illustrate key steps of the algorithms. The new algorithms are implemented in Mathematica. The package DDESpecialSolutions.m can be used to automatically compute new special solutions of nonlinear PDEs. Use of the package, implementation issues, scope, limitations, and future extensions of the software are addressed. A survey is given of related algorithms and symbolic software to compute exact solutions of nonlinear differential equations.Comment: 39 pages. Software available from Willy Hereman's home page at http://www.mines.edu/fs_home/whereman

A Collection of Challenging Optimization Problems in Science, Engineering and Economics

Function optimization and finding simultaneous solutions of a system of nonlinear equations (SNE) are two closely related and important optimization problems. However, unlike in the case of function optimization in which one is required to find the global minimum and sometimes local minima, a database of challenging SNEs where one is required to find stationary points (extrama and saddle points) is not readily available. In this article, we initiate building such a database of important SNE (which also includes related function optimization problems), arising from Science, Engineering and Economics. After providing a short review of the most commonly used mathematical and computational approaches to find solutions of such systems, we provide a preliminary list of challenging problems by writing the Mathematical formulation down, briefly explaning the origin and importance of the problem and giving a short account on the currently known results, for each of the problems. We anticipate that this database will not only help benchmarking novel numerical methods for solving SNEs and function optimization problems but also will help advancing the corresponding research areas.Comment: Accepted as an invited contribution to the special session on Evolutionary Computation for Nonlinear Equation Systems at the 2015 IEEE Congress on Evolutionary Computation (at Sendai International Center, Sendai, Japan, from 25th to 28th May, 2015.

Symbolic computation of hyperbolic tangent solutions for nonlinear differential-difference equations

A new algorithm is presented to find exact traveling wave solutions of differential-difference equations in terms of tanh functions. For systems with parameters, the algorithm determines the conditions on the parameters so that the equations might admit polynomial solutions in tanh. Examples illustrate the key steps of the algorithm. Parallels are drawn through discussion and example to the tanh-method for partial differential equations. The new algorithm is implemented in Mathematica. The package DDESpecialSolutions.m can be used to automatically compute traveling wave solutions of nonlinear polynomial differential-difference equations. Use of the package, implementation issues, scope, and limitations of the software are addressed.Comment: 19 pages submitted to Computer Physics Communications. The software can be downloaded at http://www.mines.edu/fs_home/wherema