The Symmetric Galerkin Boundary Element Method is advantageous for the linear elastic fracture and crackgrowth analysis of solid structures,because only boundary and crack-surface elements are needed.However,for engin...The Symmetric Galerkin Boundary Element Method is advantageous for the linear elastic fracture and crackgrowth analysis of solid structures,because only boundary and crack-surface elements are needed.However,for engineering structures subjected to body forces such as rotational inertia and gravitational loads,additional domain integral terms in the Galerkin boundary integral equation will necessitate meshing of the interior of the domain.In this study,weakly-singular SGBEM for fracture analysis of three-dimensional structures considering rotational inertia and gravitational forces are developed.By using divergence theorem or alternatively the radial integration method,the domain integral terms caused by body forces are transformed into boundary integrals.And due to the weak singularity of the formulated boundary integral equations,a simple Gauss-Legendre quadrature with a few integral points is sufficient for numerically evaluating the SGBEM equations.Some numerical examples are presented to verify this approach and results are compared with benchmark solutions.展开更多
Differential quadrature method is employed by numerous researchers due to its numerical accuracy and computational efficiency,and is mentioned as potential alternative of conventional numerical methods.In this paper,a...Differential quadrature method is employed by numerous researchers due to its numerical accuracy and computational efficiency,and is mentioned as potential alternative of conventional numerical methods.In this paper,a differential quadrature based numerical scheme is developed for solving volterra partial integro-differential equation of second order having a weakly singular kernel.The scheme uses cubic trigonometric B-spline functions to determine the weighting coefficients in the differential quadrature approximation of the second order spatial derivative.The advantage of this approximation is that it reduces the problem to a first order time dependent integro-differential equation(IDE).The proposed scheme is obtained in the form of an algebraic system by reducing the time dependent IDE through unconditionally stable Euler backward method as time integrator.The scheme is validated using a homogeneous and two nonhomogeneous test problems.Conditioning of the system matrix and numerical convergence of the method are analyzed for spatial and temporal domain discretization parameters.Comparison of results of the present approach with Sinc collocation method and quasi-wavelet method are also made.展开更多
In this paper,a Jacobi-collocation spectral method is developed for a Volterraintegro-differential equation with delay,which contains a weakly singular kernel.We use a function transformation and a variable transforma...In this paper,a Jacobi-collocation spectral method is developed for a Volterraintegro-differential equation with delay,which contains a weakly singular kernel.We use a function transformation and a variable transformation to change the equation into a new Volterra integral equation defined on the standard interval[-1,1],so that the Jacobi orthogonal polynomial theory can be applied conveniently.In order to obtain high order accuracy for the approximation,the integral term in the resulting equation is approximated by Jacobi spectral quadrature rules.In the end,we provide a rigorous error analysis for the proposed method.The spectral rate of convergence for the proposed method is established in both the L^(∞)-norm and the weighted L^(2)-norm.展开更多
In this study, a revised version of some numerical methods for a class of hybrid integro-differential equations with weakly singular kernels (Abel types) is presented. These equations were developed from a class of in...In this study, a revised version of some numerical methods for a class of hybrid integro-differential equations with weakly singular kernels (Abel types) is presented. These equations were developed from a class of integro-differential equations of first kind originating from an aeroelasticity problem. By manipulating the bounds of initial conditions with random variations, this study numerically demonstrated the well-posedness properties of the equations. Finally, an assumption of separating variables, allowed for linear splines to be chosen as a basis and for the differentiation and integration of the integro-differential part to be interchanged;hence, a numerical scheme was constructed.展开更多
For fractional Volterra integro-differential equations(FVIDEs)with weakly singular kernels,this paper proposes a generalized Jacobi spectral Galerkin method.The basis functions for the provided method are selected gen...For fractional Volterra integro-differential equations(FVIDEs)with weakly singular kernels,this paper proposes a generalized Jacobi spectral Galerkin method.The basis functions for the provided method are selected generalized Jacobi functions(GJFs),which can be utilized as natural basis functions of spectral methods for weakly singular FVIDEs when appropriately constructed.The developed method's spectral rate of convergence is determined using the L^(∞)-norm and the weighted L^(2)-norm.Numerical results indicate the usefulness of the proposed method.展开更多
We establish the global existence of small-amplitude solutions near a global Maxwellian to the Cauchy problem of the Vlasov-Maxwell-Boltzmann system for non-cutoff soft potentials with weak angular singularity. This e...We establish the global existence of small-amplitude solutions near a global Maxwellian to the Cauchy problem of the Vlasov-Maxwell-Boltzmann system for non-cutoff soft potentials with weak angular singularity. This extends the work of Duan et al.(2013), in which the case of strong angular singularity is considered, to the case of weak angular singularity.展开更多
This article studies the development of two numerical techniques for solving convection-diffusion type partial integro-differential equation(PIDE)with a weakly singular kernel.Cubic trigonometric B-spline(CTBS)functio...This article studies the development of two numerical techniques for solving convection-diffusion type partial integro-differential equation(PIDE)with a weakly singular kernel.Cubic trigonometric B-spline(CTBS)functions are used for interpolation in both methods.The first method is CTBS based collocation method which reduces the PIDE to an algebraic tridiagonal system of linear equations.The other method is CTBS based differential quadrature method which converts the PIDE to a system of ODEs by computing spatial derivatives as weighted sum of function values.An efficient tridiagonal solver is used for the solution of the linear system obtained in the first method as well as for determination of weighting coefficients in the second method.An explicit scheme is employed as time integrator to solve the system of ODEs obtained in the second method.The methods are tested with three nonhomogeneous problems for their validation.Stability,computational efficiency and numerical convergence of the methods are analyzed.Comparison of errors in approximations produced by the present methods versus different values of discretization parameters and convection-diffusion coefficients are made.Convection and diffusion dominant cases are discussed in terms of Peclet number.The results are also compared with cubic B-spline collocation method.展开更多
We target here to solve numerically a class of nonlinear fractional two-point boundary value problems involving left-and right-sided fractional derivatives.The main ingredient of the proposed method is to recast the p...We target here to solve numerically a class of nonlinear fractional two-point boundary value problems involving left-and right-sided fractional derivatives.The main ingredient of the proposed method is to recast the problem into an equivalent system of weakly singular integral equations.Then,a Legendre-based spectral collocation method is developed for solving the transformed system.Therefore,we can make good use of the advantages of the Gauss quadrature rule.We present the construction and analysis of the collocation method.These results can be indirectly applied to solve fractional optimal control problems by considering the corresponding Euler–Lagrange equations.Two numerical examples are given to confirm the convergence analysis and robustness of the scheme.展开更多
In this paper,we propose a hybrid spectral method for a type of nonlocal problems,nonlinear Volterra integral equations(VIEs)of the second kind.The main idea is to use the shifted generalized Log orthogonal functions(...In this paper,we propose a hybrid spectral method for a type of nonlocal problems,nonlinear Volterra integral equations(VIEs)of the second kind.The main idea is to use the shifted generalized Log orthogonal functions(GLOFs)as the basis for the first interval and employ the classical shifted Legendre polynomials for other subintervals.This method is robust for VIEs with weakly singular kernel due to the GLOFs can efficiently approximate one-point singular functions as well as smooth functions.The well-posedness and the related error estimates will be provided.Abundant numerical experiments will verify the theoretical results and show the high-efficiency of the new hybrid spectral method.展开更多
This paper is concerned with a compact difference scheme with the truncation error of order 3/2 for time and order 4 for space to an evolution equation with a weakly singular kernel.The integral term is treated by mea...This paper is concerned with a compact difference scheme with the truncation error of order 3/2 for time and order 4 for space to an evolution equation with a weakly singular kernel.The integral term is treated by means of the second order convolution quadrature suggested by Lubich.The stability and convergence are proved by the energy method.A numerical experiment is reported to verify the theoretical predictions.展开更多
The theory of a class of spectral methods is extended to Volterra integrodifferential equations which contain a weakly singular kernel(t−s)^(−μ) with 0<μ<1.In this work,we consider the case when the underlying...The theory of a class of spectral methods is extended to Volterra integrodifferential equations which contain a weakly singular kernel(t−s)^(−μ) with 0<μ<1.In this work,we consider the case when the underlying solutions of weakly singular Volterra integro-differential equations are sufficiently smooth.We provide a rigorous error analysis for the spectral methods,which shows that both the errors of approximate solutions and the errors of approximate derivatives of the solutions decay exponentially in L^(∞)-norm and weighted L^(2)-norm.The numerical examples are given to illustrate the theoretical results.展开更多
This paper is concerned with obtaining an approximate solution and an approximate derivative of the solution for neutral Volterra integro-differential equation with a weakly singular kernel.The solution of this equati...This paper is concerned with obtaining an approximate solution and an approximate derivative of the solution for neutral Volterra integro-differential equation with a weakly singular kernel.The solution of this equation,even for analytic data,is not smooth on the entire interval of integration.The Jacobi collocation discretization is proposed for the given equation.A rigorous analysis of error bound is also provided which theoretically justifies that both the error of approximate solution and the error of approximate derivative of the solution decay exponentially in L∞norm and weighted L2 norm.Numerical results are presented to demonstrate the effectiveness of the spectral method.展开更多
The z-transform is introduced to analyze a full discretization method fora partial integro-differential equation (PIDE) with a weakly singular kernel. In thismethod, spectral collocation is used for the spatial discre...The z-transform is introduced to analyze a full discretization method fora partial integro-differential equation (PIDE) with a weakly singular kernel. In thismethod, spectral collocation is used for the spatial discretization, and, for the time stepping, the finite difference method combined with the convolution quadrature rule isconsidered. The global stability and convergence properties of complete discretizationare derived and numerical experiments are reported.展开更多
In this paper, a first order semi-discrete method of a partial integro-differential equation with a weakly singular kernel is considered. We apply Galerkin spectral method in one direction, and the inversion technique...In this paper, a first order semi-discrete method of a partial integro-differential equation with a weakly singular kernel is considered. We apply Galerkin spectral method in one direction, and the inversion technique for the Laplace transform in another direction, the result of the numerical experiment proves the accuracy of this method.展开更多
In this paper,the discontinuous Galerkin method is applied to solve the multi-pantograph delay differential equations.We analyze the optimal global convergence and local superconvergence for smooth solutions under uni...In this paper,the discontinuous Galerkin method is applied to solve the multi-pantograph delay differential equations.We analyze the optimal global convergence and local superconvergence for smooth solutions under uniform meshes.Due to the initial singularity of the forcing term f,solutions of multi-pantograph delay differential equations are singular.We obtain the relevant global convergence and local superconvergence for weakly singular solutions under graded meshes.The numerical examples are provided to illustrate our theoretical results.展开更多
We present a new numerical method for solving two-dimensional Stokes flow with deformable interfaces such as dynamics of suspended drops or bubbles.The method is based on a boundary integral formulation for the interf...We present a new numerical method for solving two-dimensional Stokes flow with deformable interfaces such as dynamics of suspended drops or bubbles.The method is based on a boundary integral formulation for the interfacial velocity and is spectrally accurate in space.We analyze the singular behavior of the integrals(single-layer and double-layer integrals)appearing in the equations.The interfaces are formulated in the tangent angle and arc-length coordinates and,to reduce the stiffness of the evolution equation,the marker points are evenly distributed in arc-length by choosing a proper tangential velocity along the interfaces.Examples of Stokes flow with bubbles are provided to demonstrate the accuracy and effectiveness of the numerical method.展开更多
基金support of the National Natural Science Foundation of China(12072011).
文摘The Symmetric Galerkin Boundary Element Method is advantageous for the linear elastic fracture and crackgrowth analysis of solid structures,because only boundary and crack-surface elements are needed.However,for engineering structures subjected to body forces such as rotational inertia and gravitational loads,additional domain integral terms in the Galerkin boundary integral equation will necessitate meshing of the interior of the domain.In this study,weakly-singular SGBEM for fracture analysis of three-dimensional structures considering rotational inertia and gravitational forces are developed.By using divergence theorem or alternatively the radial integration method,the domain integral terms caused by body forces are transformed into boundary integrals.And due to the weak singularity of the formulated boundary integral equations,a simple Gauss-Legendre quadrature with a few integral points is sufficient for numerically evaluating the SGBEM equations.Some numerical examples are presented to verify this approach and results are compared with benchmark solutions.
文摘Differential quadrature method is employed by numerous researchers due to its numerical accuracy and computational efficiency,and is mentioned as potential alternative of conventional numerical methods.In this paper,a differential quadrature based numerical scheme is developed for solving volterra partial integro-differential equation of second order having a weakly singular kernel.The scheme uses cubic trigonometric B-spline functions to determine the weighting coefficients in the differential quadrature approximation of the second order spatial derivative.The advantage of this approximation is that it reduces the problem to a first order time dependent integro-differential equation(IDE).The proposed scheme is obtained in the form of an algebraic system by reducing the time dependent IDE through unconditionally stable Euler backward method as time integrator.The scheme is validated using a homogeneous and two nonhomogeneous test problems.Conditioning of the system matrix and numerical convergence of the method are analyzed for spatial and temporal domain discretization parameters.Comparison of results of the present approach with Sinc collocation method and quasi-wavelet method are also made.
基金supported by the State Key Program of National Natural Science Foundation of China(11931003)the National Natural Science Foundation of China(41974133,11671157)。
文摘In this paper,a Jacobi-collocation spectral method is developed for a Volterraintegro-differential equation with delay,which contains a weakly singular kernel.We use a function transformation and a variable transformation to change the equation into a new Volterra integral equation defined on the standard interval[-1,1],so that the Jacobi orthogonal polynomial theory can be applied conveniently.In order to obtain high order accuracy for the approximation,the integral term in the resulting equation is approximated by Jacobi spectral quadrature rules.In the end,we provide a rigorous error analysis for the proposed method.The spectral rate of convergence for the proposed method is established in both the L^(∞)-norm and the weighted L^(2)-norm.
文摘In this study, a revised version of some numerical methods for a class of hybrid integro-differential equations with weakly singular kernels (Abel types) is presented. These equations were developed from a class of integro-differential equations of first kind originating from an aeroelasticity problem. By manipulating the bounds of initial conditions with random variations, this study numerically demonstrated the well-posedness properties of the equations. Finally, an assumption of separating variables, allowed for linear splines to be chosen as a basis and for the differentiation and integration of the integro-differential part to be interchanged;hence, a numerical scheme was constructed.
基金supported by the State Key Program of National Natural Science Foundation of China(Grant No.11931003)by the National Natural Science Foundation of China(Grant Nos.41974133,12126325)by the Postgraduate Scientific Research Innovation Project of Hunan Province(Grant No.CX20200620).
文摘For fractional Volterra integro-differential equations(FVIDEs)with weakly singular kernels,this paper proposes a generalized Jacobi spectral Galerkin method.The basis functions for the provided method are selected generalized Jacobi functions(GJFs),which can be utilized as natural basis functions of spectral methods for weakly singular FVIDEs when appropriately constructed.The developed method's spectral rate of convergence is determined using the L^(∞)-norm and the weighted L^(2)-norm.Numerical results indicate the usefulness of the proposed method.
基金supported by the Fundamental Research Funds for the Central UniversitiesNational Natural Science Foundation of China(Grant Nos.11601169,11471142,11271160,11571063,11731008 and 11671309)
文摘We establish the global existence of small-amplitude solutions near a global Maxwellian to the Cauchy problem of the Vlasov-Maxwell-Boltzmann system for non-cutoff soft potentials with weak angular singularity. This extends the work of Duan et al.(2013), in which the case of strong angular singularity is considered, to the case of weak angular singularity.
文摘This article studies the development of two numerical techniques for solving convection-diffusion type partial integro-differential equation(PIDE)with a weakly singular kernel.Cubic trigonometric B-spline(CTBS)functions are used for interpolation in both methods.The first method is CTBS based collocation method which reduces the PIDE to an algebraic tridiagonal system of linear equations.The other method is CTBS based differential quadrature method which converts the PIDE to a system of ODEs by computing spatial derivatives as weighted sum of function values.An efficient tridiagonal solver is used for the solution of the linear system obtained in the first method as well as for determination of weighting coefficients in the second method.An explicit scheme is employed as time integrator to solve the system of ODEs obtained in the second method.The methods are tested with three nonhomogeneous problems for their validation.Stability,computational efficiency and numerical convergence of the methods are analyzed.Comparison of errors in approximations produced by the present methods versus different values of discretization parameters and convection-diffusion coefficients are made.Convection and diffusion dominant cases are discussed in terms of Peclet number.The results are also compared with cubic B-spline collocation method.
基金The Russian Foundation for Basic Research(RFBR)Grant No.19-01-00019.
文摘We target here to solve numerically a class of nonlinear fractional two-point boundary value problems involving left-and right-sided fractional derivatives.The main ingredient of the proposed method is to recast the problem into an equivalent system of weakly singular integral equations.Then,a Legendre-based spectral collocation method is developed for solving the transformed system.Therefore,we can make good use of the advantages of the Gauss quadrature rule.We present the construction and analysis of the collocation method.These results can be indirectly applied to solve fractional optimal control problems by considering the corresponding Euler–Lagrange equations.Two numerical examples are given to confirm the convergence analysis and robustness of the scheme.
基金The research of C.Zhang is partially supported by NSFC(Grant Nos.11971207,12071172)the Natural Science Foundation of the Jiangsu Higher Education Institutions of China(Grant No.20KJA11002)The research of S.Chen is partially supported by NSFC(Grant No.11801235).
文摘In this paper,we propose a hybrid spectral method for a type of nonlocal problems,nonlinear Volterra integral equations(VIEs)of the second kind.The main idea is to use the shifted generalized Log orthogonal functions(GLOFs)as the basis for the first interval and employ the classical shifted Legendre polynomials for other subintervals.This method is robust for VIEs with weakly singular kernel due to the GLOFs can efficiently approximate one-point singular functions as well as smooth functions.The well-posedness and the related error estimates will be provided.Abundant numerical experiments will verify the theoretical results and show the high-efficiency of the new hybrid spectral method.
基金supported by the National Natural Science Foundation of China(10971062)the Scientific Research Foundation of Central South University of Forestry and Technology.
文摘This paper is concerned with a compact difference scheme with the truncation error of order 3/2 for time and order 4 for space to an evolution equation with a weakly singular kernel.The integral term is treated by means of the second order convolution quadrature suggested by Lubich.The stability and convergence are proved by the energy method.A numerical experiment is reported to verify the theoretical predictions.
基金This work is supported by the Foundation for Talent Introduction of Guangdong Provincial University,Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme(2008)National Science Foundation of China(10971074).
文摘The theory of a class of spectral methods is extended to Volterra integrodifferential equations which contain a weakly singular kernel(t−s)^(−μ) with 0<μ<1.In this work,we consider the case when the underlying solutions of weakly singular Volterra integro-differential equations are sufficiently smooth.We provide a rigorous error analysis for the spectral methods,which shows that both the errors of approximate solutions and the errors of approximate derivatives of the solutions decay exponentially in L^(∞)-norm and weighted L^(2)-norm.The numerical examples are given to illustrate the theoretical results.
基金supported by Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme(2008)National Science Foundation of China(10971074)Specialized Research Fund for the Doctoral Program of Higher Education(20114407110009).
文摘This paper is concerned with obtaining an approximate solution and an approximate derivative of the solution for neutral Volterra integro-differential equation with a weakly singular kernel.The solution of this equation,even for analytic data,is not smooth on the entire interval of integration.The Jacobi collocation discretization is proposed for the given equation.A rigorous analysis of error bound is also provided which theoretically justifies that both the error of approximate solution and the error of approximate derivative of the solution decay exponentially in L∞norm and weighted L2 norm.Numerical results are presented to demonstrate the effectiveness of the spectral method.
基金The first author was supported in part by Scientific Research Fund of Hunan Provincial Education Department of China(10C0654)the NSF of China(10971059,11101136)+3 种基金the NSF of Hunan Province,China(10JJ6003)the Grant of Science and Technology Commission of Hunan Province,China(2012FJ4116)the NSF of Hunan University of Technology(2011HZX17)The second author was supported in part by NSF of China(10271046,10971062).
文摘The z-transform is introduced to analyze a full discretization method fora partial integro-differential equation (PIDE) with a weakly singular kernel. In thismethod, spectral collocation is used for the spatial discretization, and, for the time stepping, the finite difference method combined with the convolution quadrature rule isconsidered. The global stability and convergence properties of complete discretizationare derived and numerical experiments are reported.
基金Supported by the Natural Science Foundation of Fujian Province (2001J009, Z0511015).
文摘In this paper, a first order semi-discrete method of a partial integro-differential equation with a weakly singular kernel is considered. We apply Galerkin spectral method in one direction, and the inversion technique for the Laplace transform in another direction, the result of the numerical experiment proves the accuracy of this method.
基金supported by the Natural Science Foundation of China(No.11571027),the International Research Cooperation Seed of Beijing University of Technology(No.2018B32)Science and Technology Projects of Beijing Education Commission Foundatio(No.KM201510005032),and the 16th graduate science and technology fund of Beijing university of technology(No.ykj-2017-00127).
文摘In this paper,the discontinuous Galerkin method is applied to solve the multi-pantograph delay differential equations.We analyze the optimal global convergence and local superconvergence for smooth solutions under uniform meshes.Due to the initial singularity of the forcing term f,solutions of multi-pantograph delay differential equations are singular.We obtain the relevant global convergence and local superconvergence for weakly singular solutions under graded meshes.The numerical examples are provided to illustrate our theoretical results.
基金supported by the grants NSF-DMS 0511411,0914923 and 0923111.
文摘We present a new numerical method for solving two-dimensional Stokes flow with deformable interfaces such as dynamics of suspended drops or bubbles.The method is based on a boundary integral formulation for the interfacial velocity and is spectrally accurate in space.We analyze the singular behavior of the integrals(single-layer and double-layer integrals)appearing in the equations.The interfaces are formulated in the tangent angle and arc-length coordinates and,to reduce the stiffness of the evolution equation,the marker points are evenly distributed in arc-length by choosing a proper tangential velocity along the interfaces.Examples of Stokes flow with bubbles are provided to demonstrate the accuracy and effectiveness of the numerical method.