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 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 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.

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.

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.

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.

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.

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.

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.

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.

A Jacobi Spectral Collocation Method for Solving Fractional Integro-Differential Equations

The aim of this paper is to obtain the numerical solutions of fractional Volterra integrodifferential equations by the Jacobi spectral collocation method using the Jacobi-Gauss collocation points.We convert the fractional order integro-differential equation into integral equation by fractional order integral,and transfer the integro equations into a system of linear equations by the Gausssian quadrature.We furthermore perform the convergence analysis and prove the spectral accuracy of the proposed method in L∞norm.Two numerical examples demonstrate the high accuracy and fast convergence of the method at last.

Jacobi Collocation Methods for Solving Generalized Space-Fractional Burgers Equations

The aim of this paper is to obtain the numerical solutions of generalized space-fractional Burgers' equations with initial-boundary conditions by the Jacobi spectral collocation method using the shifted Jacobi-Gauss-Lobatto collocation points. By means of the simplifed Jacobi operational matrix, we produce the diferentiation matrix and transfer the space-fractional Burgers' equation into a system of ordinary diferential equations that can be solved by the fourth-order Runge-Kutta method. The numerical simulations indicate that the Jacobi spectral collocation method is highly accurate and fast convergent for the generalized space-fractional Burgers' equation.

A numerical approach based on the meshless collocation method in elastodynamics

In this paper, a collocation technique with the modified equilibrium on line method （ELM） for imposition of Neumann （natural） boundary conditions is presented for solving the two-dimensional problems of linear elastic body vibrations. In the modified ELM, equilibrium over the lines on the natural boundary is satisfied as Neumann boundary condition equations. In other words, the natural boundary conditions are satisfied naturally by using the weak formulation. The performance of the modified version of the ELM is studied for collocation methods based on two different ways to construct meshless shape functions： moving least squares approximation and radial basis point interpolation. Numerical examples of two-dimensional free and forced vibration analyses show that by using the modified ELM, more stable and accurate results would be obtained in comparison with the direct collocation method.

The Adaptive Wavelet Collocation Method and Its Application in Front Simulation

The adaptive wavelet collocation method （AWCM） is a variable grid technology for solving partial differential equations （PDEs） with high singularities. Based on interpolating wavelets, the AWCM adapts the grid so that a higher resolution is automatically attributed to domain regions with high singularities. Accuracy problems with the AWCM have been reported in the literature, and in this paper problems of efficiency with the AWCM are discussed in detail through a simple one-dimensional （1D） nonlinear advection equation whose analytic solution is easily obtained. A simple and efficient implementation of the AWCM is investigated. Through studying the maximum errors at the moment of frontogenesis of the 1D nonlinear advection equation with different initial values and a comparison with the finite difference method （FDM） on a uniform grid, the AWCM shows good potential for modeling the front efficiently. The AWCM is also applied to a two-dimensional （2D） unbalanced frontogenesis model in its first attempt at numerical simulation of a meteorological front. Some important characteristics about the model are revealed by the new scheme.

Stability and dispersion analysis of reproducing kernel collocation method for transient dynamics

A reproducing kernel collocation method based on strong formulation is introduced for transient dynamics. To study the stability property of this method, an algorithm based on the von Neumann hypothesis is proposed to predict the critical time step. A numerical test is conducted to validate the algorithm. The numerical critical time step and the predicted critical time step are in good agreement. The results are compared with those obtained based on the radial basis collocation method, and they axe in good agreement. Several important conclusions for choosing a proper support size of the reproducing kernel shape function are given to improve the stability condition.

A Meshless Collocation Method with Barycentric Lagrange Interpolation for Solving the Helmholtz Equation

In this paper,Chebyshev interpolation nodes and barycentric Lagrange interpolation basis function are used to deduce the scheme for solving the Helmholtz equation.First of all,the interpolation basis function is applied to treat the spatial variables and their partial derivatives,and the collocation method for solving the second order differential equations is established.Secondly,the differential matrix is used to simplify the given differential equations on a given test node.Finally,based on three kinds of test nodes,numerical experiments show that the present scheme can not only calculate the high wave numbers problems,but also calculate the variable wave numbers problems.In addition,the algorithm has the advantages of high calculation accuracy,good numerical stability and less time consuming.

Higher Order Collocation Methods for Nonlocal Problems and Their Asymptotic Compatibility

We study the convergence and asymptotic compatibility of higher order collocation methods for nonlocal operators inspired by peridynamics,a nonlocal formulation of continuum mechanics.We prove that the methods are optimally convergent with respect to the polynomial degree of the approximation.A numerical method is said to be asymptotically compatible if the sequence of approximate solutions of the nonlocal problem converges to the solution of the corresponding local problem as the horizon and the grid sizes simultaneously approach zero.We carry out a calibration process via Taylor series expansions and a scaling of the nonlocal operator via a strain energy density argument to ensure that the resulting collocation methods are asymptotically compatible.We fnd that,for polynomial degrees greater than or equal to two,there exists a calibration constant independent of the horizon size and the grid size such that the resulting collocation methods for the nonlocal difusion are asymptotically compatible.We verify these fndings through extensive numerical experiments.

THE SOLUTION OF A CRACK EMANATING FROM AN ARBITRARY HOLE BY BOUNDARY COLLOCATION METHOD

In this paper a group of stress functions has been proposed for the calculation of a crack emanating from a hole with different shape (including circular, elliptical, rectangular, or rhombic hole) by boundary collocation method. The calculation results show that they coincide very well with the existing solutions by other methods for a circular or elliptical hole with a crack in an infinite plate. At the smae time, a series of results for different holes in a finite plate has also been obtained in this paper. The proposed functions and calculation procedure can be used for a plate of a crack emanating from an arbitrary hole.