Review of Collocation Methods and Applications in Solving Science and Engineering Problems

The collocation method is a widely used numerical method for science and engineering problems governed by partial differential equations.This paper provides a comprehensive review of collocation methods and their applications,focused on elasticity,heat conduction,electromagnetic field analysis,and fluid dynamics.The merits of the collocation method can be attributed to the need for element mesh,simple implementation,high computational efficiency,and ease in handling irregular domain problems since the collocation method is a type of node-based numerical method.Beginning with the fundamental principles of the collocation method,the discretization process in the continuous domain is elucidated,and how the collocation method approximation solutions for solving differential equations are explained.Delving into the historical development of the collocation methods,their earliest applications and key milestones are traced,thereby demonstrating their evolution within the realm of numerical computation.The mathematical foundations of collocation methods,encompassing the selection of interpolation functions,definition of weighting functions,and derivation of integration rules,are examined in detail,emphasizing their significance in comprehending the method’s effectiveness and stability.At last,the practical application of the collocation methods in engineering contexts is emphasized,including heat conduction simulations,electromagnetic coupled field analysis,and fluid dynamics simulations.These specific case studies can underscore collocation method’s broad applicability and effectiveness in addressing complex engineering challenges.In conclusion,this paper puts forward the future development trend of the collocation method through rigorous analysis and discussion,thereby facilitating further advancements in research and practical applications within these fields.

A Collocation Technique via Pell-Lucas Polynomials to Solve Fractional Differential EquationModel for HIV/AIDS with Treatment Compartment

In this study,a numerical method based on the Pell-Lucas polynomials(PLPs)is developed to solve the fractional order HIV/AIDS epidemic model with a treatment compartment.The HIV/AIDS mathematical model with a treatment compartment is divided into five classes,namely,susceptible patients(S),HIV-positive individuals(I),individuals with full-blown AIDS but not receiving ARV treatment(A),individuals being treated(T),and individuals who have changed their sexual habits sufficiently(R).According to the method,by utilizing the PLPs and the collocation points,we convert the fractional order HIV/AIDS epidemic model with a treatment compartment into a nonlinear system of the algebraic equations.Also,the error analysis is presented for the Pell-Lucas approximation method.The aim of this study is to observe the behavior of five populations after 200 days when drug treatment is applied to HIV-infectious and full-blown AIDS people.To demonstrate the usefulness of this method,the applications are made on the numerical example with the help of MATLAB.In addition,four cases of the fractional order derivative(p=1,p=0.95,p=0.9,p=0.85)are examined in the range[0,200].Owing to applications,we figured out that the outcomes have quite decent errors.Also,we understand that the errors decrease when the value of N increases.The figures in this study are created in MATLAB.The outcomes indicate that the presented method is reasonably sufficient and correct.

A Numerical Investigation Based on Exponential Collocation Method for Nonlinear SITR Model of COVID-19

In this work,the exponential approximation is used for the numerical simulation of a nonlinear SITR model as a system of differential equations that shows the dynamics of the new coronavirus(COVID-19).The SITR mathematical model is divided into four classes using fractal parameters for COVID-19 dynamics,namely,susceptible(S),infected(I),treatment(T),and recovered(R).The main idea of the presented method is based on the matrix representations of the exponential functions and their derivatives using collocation points.To indicate the usefulness of this method,we employ it in some cases.For error analysis of the method,the residual of the solutions is reviewed.The reported examples show that the method is reasonably efficient and accurate.

THE ACCURACY COMPARISON BETWEEN CHEBYSHEV-τ METHOD AND CHEBYSHEV COLLOCATION METHOD

采用伪谱方法中的Ｃｈｅｂｙｓｈｅｖ－τ方法和配置点方法，分别研究了平面Ｐｏｉｓｅｕｉｌｅ流中扰动的非线性演化问题，并进行了比较．结果表明，在中性情况附近，配置点方法的精度要高于Ｃｈｅｂｙｓｈｅｖ－τ方法的精度．

A Deep Collocation Method for the Bending Analysis of Kirchhoff Plate

In this paper,a deep collocation method(DCM)for thin plate bending problems is proposed.This method takes advantage of computational graphs and backpropagation algorithms involved in deep learning.Besides,the proposed DCM is based on a feedforward deep neural network(DNN)and differs from most previous applications of deep learning for mechanical problems.First,batches of randomly distributed collocation points are initially generated inside the domain and along the boundaries.A loss function is built with the aim that the governing partial differential equations(PDEs)of Kirchhoff plate bending problems,and the boundary/initial conditions are minimised at those collocation points.A combination of optimizers is adopted in the backpropagation process to minimize the loss function so as to obtain the optimal hyperparameters.In Kirchhoff plate bending problems,the C^1 continuity requirement poses significant difficulties in traditional mesh-based methods.This can be solved by the proposed DCM,which uses a deep neural network to approximate the continuous transversal deflection,and is proved to be suitable to the bending analysis of Kirchhoff plate of various geometries.

THE COLLOCATION METHODS FOR SINGULAR INTEGRAL EQUATIONS WITH CAUCHY KERNELS

This paper applies the singular integral operators, singular quadrature operators and discretization matrices associated with singular integral equations with Cauchy kernels, which are established in [1], to give a unified framework for various collocation methods of numerical solutions of singular integral equations with Cauchy kernels. Under the framework, the coincidence of the direct quadrature method and the indirect quadrature method is very simple and obvious.

Probabilistic small signal stability analysis of power system with wind power and photovoltaic power based on probability collocation method

Recently, with increasing improvements in the penetration of wind power and photovoltaic power in the world, probabilistic small signal stability analysis(PSSSA) of a power system consisting of multiple types of renewable energy has become a key problem. To address this problem, this study proposes a probabilistic collocation method(PCM)-based PSSSA for a power system consisting of wind farms and photovoltaic farms. Compared with the conventional Monte Carlo method, the proposed method meets the accuracy and precision requirements and greatly reduces the computation; therefore, it is suitable for the PSSSA of this power system. Case studies are conducted based on a 4-machine 2-area and New England systems, respectively. The simulation results show that, by reducing synchronous generator output to improve the penetration of renewable energy, the probabilistic small signal stability(PSSS) of the system is enhanced. Conversely, by removing part of the synchronous generators to improve the penetration of renewable energy, the PSSS of the system may be either enhanced or deteriorated.

ELASTO-PLASTICITY ANALYSIS BASED ON COLLOCATION WITH THE MOVING LEAST SQUARE METHOD

A meshless approach based on the moving least square method is developed for elasto-plasticity analysis,in which the incremental formulation is used.In this approach,the dis- placement shape functions are constructed by using the moving least square approximation,and the discrete governing equations for elasto-plastic material are constructed with the direct collo- cation method.The boundary conditions are also imposed by collocation.The method established is a truly meshless one,as it does not need any mesh,either for the purpose of interpolation of the solution variables,or for the purpose of construction of the discrete equations.It is simply formu- lated and very efficient,and no post-processing procedure is required to compute the derivatives of the unknown variables,since the solution from this method based on the moving least square approximation is already smooth enough.Numerical examples are given to verify the accuracy of the meshless method proposed for elasto-plasticity analysis.

FRACTURE CALCULATION OF BENDING PLATES BY BOUNDARY COLLOCATION METHOD

Fracture of Kirchhoff plates is analyzed by the theory of complex variables and boundary collocation method. The deflections, moments and shearing forces of the plates are assumed to be the functions of complex variables. The functions can satisfy a series of basic equations and governing conditions, such as the equilibrium equations in the domain, the boundary conditions on the crack surfaces and stress singularity at the crack tips. Thus, it is only necessary to consider the boundary conditions on the external boundaries of the plate, which can be approximately satisfied by the collocation method and least square technique. Different boundary conditions and loading cases of the cracked plates are analyzed and calculated. Compared to other methods, the numerical examples show that the present method has many advantages such as good accuracy and less computer time. This is an effective semi_analytical and semi_numerical method.

Application of Boundary Collocation Method to Two-Dimensional Wave Propagation

Boundary Collocation Method （BCM） based on Eigenfunction Expansion Method （EEM）, a new numerical method for solving two-dimensional wave problems, is developed. To verify the method, wave problems on a series of beaches with different geometries are solved, and the errors of the method are analyzed. The calculation firmly confirms that the results will be more precise if we choose more rational points on the beach. The application of BCM, available for the problems with irregular domains and arbitrary boundary conditions, can effectively avoid complex calculation and programming. It can be widely used in ocean engineering.

Simulation of electrically driven jet using Chebyshev collocation method

The model of electrically driven jet is governed by a series of quasi 1D dimensionless partial differential equations(PDEs).Following the method of lines,the Chebyshev collocation method is employed to discretize the PDEs and obtain a system of differential-algebraic equations(DAEs).By differentiating constrains in DAEs twice,the system is transformed into a set of ordinary differential equations(ODEs) with invariants.Then the implicit differential equations solver 'ddaskr' is used to solve the ODEs and post-stabilization is executed at the end of each step.Results show the distributions of radius,linear charge density,stretching ratio and also the horizontal velocity at a time point.Meanwhile,the spiral and expanding projections to X-Y plane of the jet centerline suggest the occurring of bending instability.

MHD Falkner-Skan flow of Maxwell fluid by rational Chebyshev collocation method

The magnetohydrodynamics （MHD） Falkner-Skan flow of the Maxwell fluid is studied. Suitable transform reduces the partial differential equation into a nonlinear three order boundary value problem over a semi-infinite interval. An efficient approach based on the rational Chebyshev collocation method is performed to find the solution to the proposed boundary value problem. The rational Chebyshev collocation method is equipped with the orthogonal rational Chebyshev function which solves the problem on the semi-infinite domain without truncating it to a finite domain. The obtained results are presented through the illustrative graphs and tables which demonstrate the affectivity, stability, and convergence of the rational Chebyshev collocation method. To check the accuracy of the obtained results, a numerical method is applied for solving the problem. The variations of various embedded parameters into the problem are examined.

Laguerre-Gauss collocation method for initial value problems of second order ODEs

This paper proposes a new collocation method for initial value problems of second order ODEs based on the Laguerre-Gauss interpolation. It provides the global numerical solutions and possesses the spectral accuracy. Numerical results demonstrate its high efficiency.

Simulation of Seawater Intrusion in Coastal Confined Aquifer Using a Point Collocation Method Based Meshfree Model

Seawater intrusion caused by groundwater over-exploitation from coastal aquifers poses a severe problem in many regions. Formulation of proper pumping strategy using a simulation model can assure sustainable supply of fresh water from the coastal aquifers. The focus of the present study is on the development of a numerical model based on Meshfree (MFree) method to study the seawater intrusion problem. For the simulation of seawater intrusion problem, widely used models are based on Finite Difference (FDM) and Finite Element (FEM) Methods, which demand well defined grids/meshes and considerable pre-processing efforts. Here, MFree Point Collocation Method (PCM) based on the Radial Basis Function (RBF) is proposed for the simulation. Diffusive interface approach with density-dependent dispersion and solution of flow and solute transport is adopted. These equations are solved using PCM with appropriate boundary conditions. The developed model has been verified with Henry’s problem, and found to be satisfactory. Further the model has been applied to another established problem and an attempt is made to examine the influence of important system parameters including pumping and recharge on the seawater intrusion. The PCM based MFree model is found computationally efficient as preprocessing is avoided when compared to other numerical methods.

Sinc-Collocation Method for Solving Linear and Nonlinear System of Second-Order Boundary Value Problems

Sinc methods are now recognized as an efficient numerical method for problems whose solutions may have singularities, or infinite domains, or boundary layers. This work deals with the sinc-collocation method for solving linear and nonlinear system of second order differential equation. The method is then tested on linear and nonlinear examples and a comparison with B-spline method is made. It is shown that the sinc-collocation method yields better results.

A New Technique for Solving Fractional Order Systems: Hermite Collocation Method

In this study, we establish an approximate method which produces an approximate Hermite polynomial solution to a system of fractional order differential equations with variable coefficients. At collocation points, this method converts the mentioned system into a matrix equation which corresponds to a system of linear equations with unknown Hermite polynomial coefficients. Construction of the method on the aforementioned type of equations has been presented and tested on some numerical examples. Results related to the effectiveness and reliability of the method have been illustrated.

Finite Element Orthogonal Collocation Approach for Time Fractional Telegraph Equation with Mamadu-Njoseh Polynomials

Finite element method (FEM) is an efficient numerical tool for the solution of partial differential equations (PDEs). It is one of the most general methods when compared to other numerical techniques. PDEs posed in a variational form over a given space, say a Hilbert space, are better numerically handled with the FEM. The FEM algorithm is used in various applications which includes fluid flow, heat transfer, acoustics, structural mechanics and dynamics, electric and magnetic field, etc. Thus, in this paper, the Finite Element Orthogonal Collocation Approach (FEOCA) is established for the approximate solution of Time Fractional Telegraph Equation (TFTE) with Mamadu-Njoseh polynomials as grid points corresponding to new basis functions constructed in the finite element space. The FEOCA is an elegant mixture of the Finite Element Method (FEM) and the Orthogonal Collocation Method (OCM). Two numerical examples are experimented on to verify the accuracy and rate of convergence of the method as compared with the theoretical results, and other methods in literature.

Numerical Solution of Troesch’s Problem by Sinc-Collocation Method

A new algorithm is presented for solving Troesch’s problem. The numerical scheme based on the sinc-collocation technique is deduced. The equation is reduced to systems of nonlinear algebraic equations. Some numerical experiments are made. Compared with the modified homotopy perturbation technique (MHP), the variational iteration method and the Adomian decomposition method. It is shown that the sinc-collocation method yields better results.

A Novel Spacetime Collocation Meshless Method for Solving Two- Dimensional Backward Heat Conduction Problems

In this article,a meshless method using the spacetime collocation for solving the two-dimensional backward heat conduction problem(BHCP)is proposed.The spacetime collocation meshless method(SCMM)is to derive the general solutions as the basis functions for the two-dimensional transient heat equation using the separation of variables.Numerical solutions of the heat conduction problem are expressed as a series using the addition theorem.Because the basis functions are the general solutions of the governing equation,the boundary points may be collocated on the spacetime boundary of the domain.The proposed method is verified by conducting several heat conduction problems.We also carry out numerical applications to compare the SCMM with other meshless methods.The results show that the SCMM is accurate and efficient.Furthermore,it is found that the recovered boundary data on inaccessible boundary can be obtained with high accuracy even though the over specified data are provided only at a 1/6 portion of the spacetime boundary.

Error Control Strategies for Numerical Integrations in Fast Collocation Methods

We propose two error control techniques for numerical integrations in fast multiscale collocation methods for solving Fredholm integral equations of the second kind with weakly singular kernels. Both techniques utilize quadratures for singular integrals using graded points. One has a polynomial order of accuracy if the integrand has a polynomial order of smoothness except at the singular point and the other has exponential order of accuracy if the integrand has an infinite order of smoothness except at the singular point. We estimate the order of convergence and computational complexity of the corresponding approximate solutions of the equation. We prove that the second technique preserves the order of convergence and computational complexity of the original collocation method. Numerical experiments are presented to illustrate the theoretical estimates.