Two numerical methods in the family of weighted residual methods; the orthogonal collocation and least squares methods, are used within the spectral framework to solve a linear reaction-diffusion pellet problem with slab and spherical geometries. The node points are in this work taken as the roots of orthogonal polynomials in the Jacobi family. Two Jacobi polynomial parameters, alpha and beta, can be used to tune the distribution of the roots within the domain. Further, the internal points and the boundary points of the boundary-value problem can be given according to: i) Gauss-Lobatto-Jacobi points, or ii) Gauss-Jacobi points plus the boundary points. The objective of this paper is thus to investigate the influence of the distribution of the node points within the domain adopting the orthogonal collocation and least squares methods. Moreover, the results of the two numerical methods are compared to examine whether the methods show the same sensitivity and accuracy to the node point distribution. The notifying findings are as follows: i) The Legendre polynomial, i.e., alpha=beta=0, is a very robust Jacobi polynomial giving the better condition number of the coefficient matrix and the polynomial also give good behavior of the error as a function of polynomial order. This polynomial gives good results for small and large gradients within both slab and spherical pellet geometries. This trend is observed for both of the weighted residual methods applied. ii) Applying the least squares method the error decreases faster with increasing polynomial order than observed with the orthogonal collocation method. However, the orthogonal collocation method is not so sensitive to the choice of Jacobi polynomial and the method also obtains lower error values than the least squares method due to favorable lower condition numbers of the coefficient matrices. Thus, for this particular problem, the orthogonal collocation method is recommended above the least squares method. iii) The orthogonal collocation method show minor differences between Gauss-Lobatto-Jacobi points and Gauss-Jacobi plus boundary points.
Direct Trajectory Optimization and Costate Estimation via an Orthogonal Collocation Method
