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.

用Velocity Verlet积分器改进HMC抽样方法

Hamilton Monte Carlo (HMC)方法是一种常用的快速抽样方法.在对哈密顿方程进行抽样时,HMC方法使用Leapfrog积分器,这可能造成方程的位置及动量的迭代值在时间上不同步,其产生的误差会降低抽样效率及抽样结果的稳定性.为此,本文提出了IHMC(Improved HMC)方法,该方法用Velocity Verlet积分器替代Leapfrog积分器,每次迭代时都计算两变量在同一时刻的值.为验证方法的效果,本文进行了两个实验,一个是将该方法应用于非对称随机波动率模型(RASV模型)的参数估计,另一个是将方法应用于方差伽马分布的抽样,结果显示:IHMC方法比HMC方法的效率更高、结果更稳定.

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.

Feasibility of the 2-point method to determine the load-velocity relationship variables during the countermovement jump exercise

Purpose:This study aimed to examine the reliability and validity of load-velocity(L-V)relationship variables obtained through the 2-point method using different load combinations and velocity variables.Methods:Twenty men performed 2 identical sessions consisting of 2 countermovement jumps against 4 external loads(20 kg,40 kg,60 kg,and80 kg)and a heavy squat against a load linked to a mean velocity(MV)of 0.55 m/s(load_(0.55)).The L-V relationship variables(load-axis intercept(L_(0)),velocity-axis intercept(v_(0)),and area under the L-V relationship line(A_(line)))were obtained using 3 velocity variables(MV,mean propulsive velocity(MPV),and peak velocity)by the multiple-point method including(20-40-60-80-load_(0.55))and excluding(20-40-60-80)the heavy squat,as well as from their respective 2-point methods(20-load_(0.55)and 20-80).Results:The L-V relationship variables were obtained with an acceptable reliability(coefncient of variation(CV)≤7.30%;intra-class correlation coefficient>0.63).The reliability of L_(0)and v_(0)was comparable for both methods(CV_(ratio)(calculated as higher value/lower value):1.11-1.12),but the multiple-point method provided Al_(ine)with a greater reliability(CV_(ratio)=1.26).The use of a heavy squat provided the L-V relationship variables with a comparable or higher reliability than the use of a heavy countermovement jump load(CV_(ratio):1.06-1.19).The peak velocity provided the load-velocity relationship variables with the greatest reliability(CV_(ratio):1.15-1.86)followed by the MV(CV_(ratio):1.07-1.18),and finally the MPV.The 2-point methods only revealed an acceptable validity for the MV and MPV(effect size≤0.19;Pearson s product-moment correlation coefficient≥0.96;Lin's concordance correlation coefficient≥0.94).Conclusion:The 2-point method obtained from a heavy squat load and MV or MPV is a quick,safe,and reliable procedure to evaluate the lower-body maximal neuromuscular capacities through the L-V relationship.

Energy-Saving and Punctuality Combined Velocity Planning for the Autonomous-Rail Rapid Tram with Enhanced Pseudospectral Method

Autonomous-rail rapid transit(ART)is a new medium-capacity rapid transportation system with punctuality,comfort and convenience,but low-cost construction.Combined velocity planning is a critical approach to meet the requirements of energy-saving and punctuality.An ART velocity pre-planning and re-planning strategy based on the combination of punctuality dynamic programming(PDP)and pseudospectral(PS)method is proposed in this paper.Firstly,the longitudinal dynamics model of ART is established by a multi-particle model.Secondly,the PDP algorithm with global optimal characteristics is adopted as the pre-planning strategy.A model for determining the number of collocation points of the real-time PS method is proposed to improve the energy-saving effect while ensuring computation efficiency.Then the enhanced PS method is utilized to design the velocity re-planning strategy.Finally,simulations are conducted in the typical scenario with sloping roads,traffic lights,and intrusion of the pedestrian.The simulation results indicate that the ART with the proposed velocity trajectory optimization strategy can meet the punctuality requirement,and obtain better economy efficiency compared with the punctuality green light optimal speed advisory(PGLOSA).

Herrmann Method of Analyzing Structure Design Velocity Field

The shape optimization is studied by adopting the domain integrated method which is based on the calculus of variations during the shape design sensitivity analysis. A new method of improving the efficiency of the design velocity field analysis and the quality of the finite element method (FEM) mesh is put forward. The sensitivity analysis which is based on the calculus of variations is used in the shape optimization. The design velocity field is solved by Herrmann method. An example shows that both the quality of the FEM mesh and the efficiency of the computing of the design velocity field are improved by Herrmann method. So the effect and the efficiency of the shape optimization are guaranteed. If using sensitivity analysis which is based on the calculus of variations in the shape optimization, the sensitivity analysis can be a relatively independent module. The efficiency of computing the design velocity field and the quality of mesh will be improved by using Herrmann method.

Low-velocity Impact Damage Analysis of Composite Laminates Using Self-adapting Delamination Element Method

On the basis ofa 2D 4-node Mindlin shell element method, a novel self-adapting delamination finite element method is presented, which is developed to model the delamination damage of composite laminates. In the method, the sublaminate elements are generated automatically when the delamination damage occurs or extends. Thus, the complex process and state of delamination damage can be simulated practically with high efficiency for both analysis and modeling. Based on the self-adapting delamination method, linear dynamic finite element damage analysis is performed to simulate the low-velocity impact damage process of three types of mixed woven composite laminates. Taking the frictional force among sublaminations during delaminating and the transverse normal stress into account, the analytical results are consistent with those of the experimental data.

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

This paper is devoted to investigate the accuracy of the Pseudo spectral scheme with the Chebyshev tau method and Chebyshev collocation method. The computational results of the nonlinear disturbance development in plane Poiseuille flow for both methods are presented and compared in detail. It is acknowledged that the Chebyshev collocation method has higher precision than the other one, especially for near netural situation.

Debris cloud structure and hazardous fragments distribution under hypervelocity yaw impact

This study investigates how the debris cloud structure and hazardous fragment distribution vary with attack angle by simulating a circular cylinder projectile hypervelocity impinging on a thin plate using the finite element-smoothed particle hydrodynamics(FE-SPH)adaptive method.Based on the comparison and analysis of the experimental and simulation results,the FE-SPH adaptive method was applied to address the hypervelocity yaw impact problem,and the variation law of the debris cloud structure with the attack angle was obtained.The screening criterion of the hazardous fragment at yaw impact is given by analyzing the debris formation obtained by the FE-SPH adaptive method,and the distribution characteristics of hazardous fragments and their relationship with the attack angle are given.Moreover,the velocity space was used to evaluate the distribution range and damage capability of asymmetric hazardous fragments.The maximum velocity angle was extended from fully symmetrical working conditions to asymmetrical cases to describe the asymmetrical debris cloud distribution range.In this range,the energy density was calculated to quantitatively analyze how much damage hazardous fragments inflict on the rear plate.The results showed that the number of hazardous fragments generated by the case near the 35°attack angle was the largest,the distribution range was the smallest,and the energy density was the largest.These results suggest that in this case,debris cloud generated by the impact had the strongest damage to the rear plate.

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.

Computing Streamfunction and Velocity Potential in a Limited Domain of Arbitrary Shape.Part II:Numerical Methods and Test Experiments

Built on the integral formulas in Part I,numerical methods are developed for computing velocity potential and streamfunction in a limited domain.When there is no inner boundary（around a data hole） inside the domain,the total solution is the sum of the internally and externally induced parts.For the internally induced part,three numerical schemes（grid-staggering,local-nesting and piecewise continuous integration） are designed to deal with the singularity of the Green＇s function encountered in numerical calculations.For the externally induced part,by setting the velocity potential（or streamfunction） component to zero,the other component of the solution can be computed in two ways：（1） Solve for the density function from its boundary integral equation and then construct the solution from the boundary integral of the density function.（2） Use the Cauchy integral to construct the solution directly.The boundary integral can be discretized on a uniform grid along the boundary.By using local-nesting（or piecewise continuous integration）,the scheme is refined to enhance the discretization accuracy of the boundary integral around each corner point（or along the entire boundary）.When the domain is not free of data holes,the total solution contains a data-hole-induced part,and the Cauchy integral method is extended to construct the externally induced solution with irregular external and internal boundaries.An automated algorithm is designed to facilitate the integrations along the irregular external and internal boundaries.Numerical experiments are performed to evaluate the accuracy and efficiency of each scheme relative to others.

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.

Signal processing method for shear wave velocity measurement

Soil shear wave velocity （SWV） is an important parameter in geotechnical engineering. To measure the soil SWV, three methods are generally used in China, including the single-hole method, cross-hole method and the surface-wave technique. An optimized approach based on a correlation function for single-hole SWV measurement is presented in this paper. In this approach, inherent inconsistencies of the artificial methods such as negative velocities, and too-large and too-small velocities, are eliminated from the single-hole method, and the efficiency of data processing is improved. In addition, verification using the cross-hole method of upper measuring points shows that the proposed optimized approach yields high precision in signal processing.

A method to model the effect of pre-existing cracks on P-wave velocity in rocks

Crack closure is one of the reasons inducing changes of P-wave velocity of rocks under compression.In this context,a method is proposed to investigate the relationships among P-wave velocity,pre-existing cracks,and confining pressure based on the discrete element method(DEM).Pre-existing open cracks inside the rocks are generated by the initial gap of the flat-joint model.The validity of the method is evaluated by comparing the P-wave velocity tested on a sandstone specimen with numerical result.As the crack size is determined by the diameter of particles,the effects of three factors,i.e.number,aspect ratio,and orientation of cracks on the P-wave velocity are discussed.The results show that P-wave velocity is controlled by the(i.e.number) of open micro-cracks,while the closure pressure is determined by the aspect ratio of crack.The reason accounting for the anisotropy of P-wave velocity is the difference in crack number in measurement paths.Both of the number and aspect ratio of cracks can affect the responses of P-wave velocity to the applied confining pressure.Under confining pressure,the number of open cracks inside rocks will dominate the lowest P-wave velocity,and the P-wave velocity of the rock containing narrower cracks is more sensitive to the confining pressure.In this sense,crack density is difficult to be back-calculated merely by P-wave velocity.The proposed method offers a means to analyze the effect of pre-existing cracks on P-wave velocity.

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.

Effect of particle shape on packing fraction and velocity profiles at outlet of a silo

Many studies on how the particle shape affects the discharge flow mainly focus on discharge rates and avalanche statistics. In this study, the effect of the particle shape on the packing fraction and velocities of particles in the silo discharge flow are investigated by using the discrete element method. The time-averaged packing fraction and velocity profiles through the aperture are systematically measured for superelliptical particles with different blockinesses. Increasing the particle blockiness is found to increase resistance to flow and reduce the flow rate. At an identical outlet size, larger particle blockiness leads to lower velocity and packing fraction at the outlet. The packing fraction profiles display evidently the self-similar feature that can be appropriately adjusted by fractional power law. The velocity profiles for particles with different shapes obey a uniform self-similar law that is in accord with previous experimental results, which is compatible with the hypothesis of free fall arch. To further investigate the origin of flow behaviors, the packing fraction and velocity field in the region above the orifice are computed. Based on these observations, the flow rate of superelliptical particles is calculated and in agreement with the simulated data.