High-order adaptive Gegenbauer integral spectral element method for solving non-linear optimal control problems

Abstract

In this work, we propose an adaptive spectral element algorithm for solving non-linear optimal control problems. The method employs orthogonal collocation at the shifted Gegenbauer–Gauss points combined with very accurate and stable numerical quadratures to fully discretize the multiple-phase integral form of the optimal control problem. The proposed algorithm relies on exploiting the underlying smoothness properties of the solutions for computing approximate solutions efficiently. In particular, the method brackets discontinuities and ‘points of nonsmoothness’ through a novel local adaptive algorithm, which achieves a desired accuracy on the discrete dynamical system equations by adjusting both the mesh size and the degree of the approximating polynomials. A rigorous error analysis of the developed numerical quadratures is presented. Finally, the efficiency of the proposed method is demonstrated on three test examples from the open literature.

Keywords:

• Adaptive strategy
• Gegenbauer polynomials
• integration matrix
• optimal control problems
• spectral element methods

Notes

No potential conflict of interest was reported by the author.

1 The kth element shifted Gegenbauer–Gauss (KESGG) nodes in , for all .

2 The points at which a certain order derivative of a real-valued function is discontinuous.