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.展开更多
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.展开更多
Deep neural networks(DNNs)are effective in solving both forward and inverse problems for nonlinear partial differential equations(PDEs).However,conventional DNNs are not effective in handling problems such as delay di...Deep neural networks(DNNs)are effective in solving both forward and inverse problems for nonlinear partial differential equations(PDEs).However,conventional DNNs are not effective in handling problems such as delay differential equations(DDEs)and delay integrodifferential equations(DIDEs)with constant delays,primarily due to their low regularity at delayinduced breaking points.In this paper,a DNN method that combines multi-task learning(MTL)which is proposed to solve both the forward and inverse problems of DIDEs.The core idea of this approach is to divide the original equation into multiple tasks based on the delay,using auxiliary outputs to represent the integral terms,followed by the use of MTL to seamlessly incorporate the properties at the breaking points into the loss function.Furthermore,given the increased training dificulty associated with multiple tasks and outputs,we employ a sequential training scheme to reduce training complexity and provide reference solutions for subsequent tasks.This approach significantly enhances the approximation accuracy of solving DIDEs with DNNs,as demonstrated by comparisons with traditional DNN methods.We validate the effectiveness of this method through several numerical experiments,test various parameter sharing structures in MTL and compare the testing results of these structures.Finally,this method is implemented to solve the inverse problem of nonlinear DIDE and the results show that the unknown parameters of DIDE can be discovered with sparse or noisy data.展开更多
This paper develops and analyzes a moving mesh finite difference method for solving partial integro-differential equations. First, the time-dependent mapping of the coordinate transformation is approximated by a a pie...This paper develops and analyzes a moving mesh finite difference method for solving partial integro-differential equations. First, the time-dependent mapping of the coordinate transformation is approximated by a a piecewise linear function in time. Then, piecewise quadratic polynomial in space and an efficient method to discretize the memory term of the equation is designed using the moving mesh approach. In each time slice, a simple piecewise constant approximation of the integrand is used, and thus a quadrature is constructed for the memory term. The central finite difference scheme for space and the backward Euler scheme for time are used. The paper proves that the accumulation of the quadrature error is uniformly bounded and that the convergence of the method is second order in space and first order in time. Numerical experiments are carried out to confirm the theoretical predictions.展开更多
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 study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled b...In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.展开更多
This paper is devoted to understanding the stability of perturbations around the hydrostatic equilibrium of the Boussinesq system in order to gain insight into certain atmospheric and oceanographic phenomena.The Bouss...This paper is devoted to understanding the stability of perturbations around the hydrostatic equilibrium of the Boussinesq system in order to gain insight into certain atmospheric and oceanographic phenomena.The Boussinesq system focused on here is anisotropic,and involves only horizontal dissipation and thermal damping.In the 2D case R^(2),due to the lack of vertical dissipation,the stability and large-time behavior problems have remained open in a Sobolev setting.For the spatial domain T×R,this paper solves the stability problem and gives the precise large-time behavior of the perturbation.By decomposing the velocity u and temperatureθinto the horizontal average(ū,θ)and the corresponding oscillation(ū,θ),we can derive the global stability in H~2 and the exponential decay of(ū,θ)to zero in H^(1).Moreover,we also obtain that(ū_(2),θ)decays exponentially to zero in H^(1),and thatū_(1)decays exponentially toū_(1)(∞)in H^(1)as well;this reflects a strongly stratified phenomenon of buoyancy-driven fluids.In addition,we establish the global stability in H^(3)for the 3D case R^(3).展开更多
In this study,we aimto investigate certain triple integral transformand its application to a class of partial differentialequations.We discuss various properties of the new transformincluding inversion, linearity, exi...In this study,we aimto investigate certain triple integral transformand its application to a class of partial differentialequations.We discuss various properties of the new transformincluding inversion, linearity, existence, scaling andshifting, etc. Then,we derive several results enfolding partial derivatives and establish amulti-convolution theorem.Further, we apply the aforementioned transform to some classical functions and many types of partial differentialequations involving heat equations,wave equations, Laplace equations, and Poisson equations aswell.Moreover,wedraw some figures to illustrate 3-D contour plots for exact solutions of some selected examples involving differentvalues in their variables.展开更多
Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational...Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational domains,etc.Typical reduced order modeling techniques accelerate the solution of the parametric PDEs by projecting them onto a linear trial manifold constructed in the ofline stage.These methods often need a predefined mesh as well as a series of precomputed solution snapshots,and may struggle to balance between the efficiency and accuracy due to the limitation of the linear ansatz.Utilizing the nonlinear representation of neural networks(NNs),we propose the Meta-Auto-Decoder(MAD)to construct a nonlinear trial manifold,whose best possible performance is measured theoretically by the decoder width.Based on the meta-learning concept,the trial manifold can be learned in a mesh-free and unsupervised way during the pre-training stage.Fast adaptation to new(possibly heterogeneous)PDE parameters is enabled by searching on this trial manifold,and optionally fine-tuning the trial manifold at the same time.Extensive numerical experiments show that the MAD method exhibits a faster convergence speed without losing the accuracy than other deep learning-based methods.展开更多
Photonic computing has recently become an interesting paradigm for high-speed calculation of computing processes using light-matter interactions.Here,we propose and study an electromagnetic wave-based structure with t...Photonic computing has recently become an interesting paradigm for high-speed calculation of computing processes using light-matter interactions.Here,we propose and study an electromagnetic wave-based structure with the ability to calculate the solution of partial differential equations(PDEs)in the form of the Helmholtz wave equation,∇^(2)(x,y)T+k^(2)(x,y)=0,with k as the wavenumber.To do this,we make use of a network of interconnected waveguides filled with dielectric inserts.In so doing,it is shown how the proposed network can mimic the response of a network of T-circuit elements formed by two series and a parallel impedances,i.e.,the waveguide network effectively behaves as a metatronic network.An in-depth theoretical analysis of the proposed metatronic structure is presented,showing how the governing equation for the currents and impedances of the metatronic network resembles that of the finite difference representation of the Helmholtz wave equation.Different studies are then discussed including the solution of PDEs for Dirichlet and open boundary value problems,demonstrating how the proposed metatronic-based structure has the ability to calculate their solutions.展开更多
We propose and analyze a spectral Jacobi-collocation approximation for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide a rigorou...We propose and analyze a spectral Jacobi-collocation approximation for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide a rigorous error analysis for the collection method, which shows that the errors of the approximate solution decay exponentially in L^∞ norm and weighted L^2-norm. The numerical examples are given to illustrate the theoretical results.展开更多
In this paper we are looking forward to finding the approximate analytical solutions for fractional integro-differential equations by using Sumudu transform method and Hermite spectral collocation method.The fractiona...In this paper we are looking forward to finding the approximate analytical solutions for fractional integro-differential equations by using Sumudu transform method and Hermite spectral collocation method.The fractional derivatives are described in the Caputo sense.The applications related to Sumudu transform method and Hermite spectral collocation method have been developed for differential equations to the extent of access to approximate analytical solutions of fractional integro-differential equations.展开更多
A lumped mass approximation scheme of a low order Crouzeix-Raviart type noncon- forming triangular finite element is proposed to a kind of nonlinear parabolic integro-differential equations. The L2 error estimate is d...A lumped mass approximation scheme of a low order Crouzeix-Raviart type noncon- forming triangular finite element is proposed to a kind of nonlinear parabolic integro-differential equations. The L2 error estimate is derived on anisotropic meshes without referring to the traditional nonclassical elliptic projection.展开更多
A highly efficient H1-Galerkin mixed finite element method (MFEM) is presented with linear triangular element for the parabolic integro-differential equation. Firstly, some new results about the integral estimation ...A highly efficient H1-Galerkin mixed finite element method (MFEM) is presented with linear triangular element for the parabolic integro-differential equation. Firstly, some new results about the integral estimation and asymptotic expansions are studied. Then, the superconvergence of order O(h^2) for both the original variable u in H1 (Ω) norm and the flux p = u in H(div, Ω) norm is derived through the interpolation post processing technique. Furthermore, with the help of the asymptotic expansions and a suitable auxiliary problem, the extrapolation solutions with accuracy O(h^3) are obtained for the above two variables. Finally, some numerical results are provided to confirm validity of the theoretical analysis and excellent performance of the proposed method.展开更多
In this paper, the author discusses the multiple positive solutions for an infinite boundary value problem of first order impulsive singular integro-differential equations on the half line by means of the fixed point ...In this paper, the author discusses the multiple positive solutions for an infinite boundary value problem of first order impulsive singular integro-differential equations on the half line by means of the fixed point theorem of cone expansion and compression with norm type.展开更多
Using the monotone iterative method and Monch Fixed point theorem, the existence of solutions and coupled minimal and maximal quasisolutions of initial value problems for mixed monotone second-order integro-differenti...Using the monotone iterative method and Monch Fixed point theorem, the existence of solutions and coupled minimal and maximal quasisolutions of initial value problems for mixed monotone second-order integro-differential equations in Banach spaces are studied. Some existence theorems of solutions and coupled minimal and maximal quasisolutions are obtained.展开更多
In this paper, we prove an important existence and uniqueness theorem for a fractional order Fredholm – Volterra integro-differential equation with non-local and global boundary conditions by converting it to the cor...In this paper, we prove an important existence and uniqueness theorem for a fractional order Fredholm – Volterra integro-differential equation with non-local and global boundary conditions by converting it to the corresponding well known Fredholm integral equation of second kind. The considered in this paper has been solved already numerically in [1].展开更多
The current work models a weak(soft) interface between two elastic materials as containing a periodic array of micro-crazes. The boundary conditions on the interfacial micro-crazes are formulated in terms of a system ...The current work models a weak(soft) interface between two elastic materials as containing a periodic array of micro-crazes. The boundary conditions on the interfacial micro-crazes are formulated in terms of a system of hypersingular integro-differential equations with unknown functions given by the displacement jumps across opposite faces of the micro-crazes. Once the displacement jumps are obtained by approximately solving the integro-differential equations, the effective stiffness of the micro-crazed interface can be readily computed. The effective stiffness is an important quantity needed for expressing the interfacial conditions in the spring-like macro-model of soft interfaces. Specific case studies are conducted to gain physical insights into how the effective stiffness of the interface may be influenced by the details of the interfacial micro-crazes.展开更多
Let X be a complex Banach space and let B and C be two closed linear operators on X satisfying the condition D(B)?D(C),and let d∈L^(1)(R_(+))and 0≤β<α≤2.We characterize the well-posedness of the fractional int...Let X be a complex Banach space and let B and C be two closed linear operators on X satisfying the condition D(B)?D(C),and let d∈L^(1)(R_(+))and 0≤β<α≤2.We characterize the well-posedness of the fractional integro-differential equations D^(α)u(t)+CD^(β)u(t)=Bu(t)+∫_(-∞)td(t-s)Bu(s)ds+f(t),(0≤t≤2π)on periodic Lebesgue-Bochner spaces L^(p)(T;X)and periodic Besov spaces B_(p,q)^(s)(T;X).展开更多
In this paper, the author discusses the multiple positive solutions for an infinite boundary value problem of first order impulsive superlinear integro-differential equations on the half line by means of the fixed poi...In this paper, the author discusses the multiple positive solutions for an infinite boundary value problem of first order impulsive superlinear integro-differential equations on the half line by means of the fixed point theorem of cone expansion and compression with norm type.展开更多
文摘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.
文摘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.
文摘Deep neural networks(DNNs)are effective in solving both forward and inverse problems for nonlinear partial differential equations(PDEs).However,conventional DNNs are not effective in handling problems such as delay differential equations(DDEs)and delay integrodifferential equations(DIDEs)with constant delays,primarily due to their low regularity at delayinduced breaking points.In this paper,a DNN method that combines multi-task learning(MTL)which is proposed to solve both the forward and inverse problems of DIDEs.The core idea of this approach is to divide the original equation into multiple tasks based on the delay,using auxiliary outputs to represent the integral terms,followed by the use of MTL to seamlessly incorporate the properties at the breaking points into the loss function.Furthermore,given the increased training dificulty associated with multiple tasks and outputs,we employ a sequential training scheme to reduce training complexity and provide reference solutions for subsequent tasks.This approach significantly enhances the approximation accuracy of solving DIDEs with DNNs,as demonstrated by comparisons with traditional DNN methods.We validate the effectiveness of this method through several numerical experiments,test various parameter sharing structures in MTL and compare the testing results of these structures.Finally,this method is implemented to solve the inverse problem of nonlinear DIDE and the results show that the unknown parameters of DIDE can be discovered with sparse or noisy data.
基金partly supported by SRF for ROCS, SEMsupported by a grant from the "project 211 (phase Ⅲ)" of the Southwestern University of Finance and Economics
文摘This paper develops and analyzes a moving mesh finite difference method for solving partial integro-differential equations. First, the time-dependent mapping of the coordinate transformation is approximated by a a piecewise linear function in time. Then, piecewise quadratic polynomial in space and an efficient method to discretize the memory term of the equation is designed using the moving mesh approach. In each time slice, a simple piecewise constant approximation of the integrand is used, and thus a quadrature is constructed for the memory term. The central finite difference scheme for space and the backward Euler scheme for time are used. The paper proves that the accumulation of the quadrature error is uniformly bounded and that the convergence of the method is second order in space and first order in time. Numerical experiments are carried out to confirm the theoretical predictions.
基金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.
文摘In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.
基金supported by National Natural Science Foundation of China(12071391,12231016)the Guangdong Basic and Applied Basic Research Foundation(2022A1515010860)。
文摘This paper is devoted to understanding the stability of perturbations around the hydrostatic equilibrium of the Boussinesq system in order to gain insight into certain atmospheric and oceanographic phenomena.The Boussinesq system focused on here is anisotropic,and involves only horizontal dissipation and thermal damping.In the 2D case R^(2),due to the lack of vertical dissipation,the stability and large-time behavior problems have remained open in a Sobolev setting.For the spatial domain T×R,this paper solves the stability problem and gives the precise large-time behavior of the perturbation.By decomposing the velocity u and temperatureθinto the horizontal average(ū,θ)and the corresponding oscillation(ū,θ),we can derive the global stability in H~2 and the exponential decay of(ū,θ)to zero in H^(1).Moreover,we also obtain that(ū_(2),θ)decays exponentially to zero in H^(1),and thatū_(1)decays exponentially toū_(1)(∞)in H^(1)as well;this reflects a strongly stratified phenomenon of buoyancy-driven fluids.In addition,we establish the global stability in H^(3)for the 3D case R^(3).
文摘In this study,we aimto investigate certain triple integral transformand its application to a class of partial differentialequations.We discuss various properties of the new transformincluding inversion, linearity, existence, scaling andshifting, etc. Then,we derive several results enfolding partial derivatives and establish amulti-convolution theorem.Further, we apply the aforementioned transform to some classical functions and many types of partial differentialequations involving heat equations,wave equations, Laplace equations, and Poisson equations aswell.Moreover,wedraw some figures to illustrate 3-D contour plots for exact solutions of some selected examples involving differentvalues in their variables.
基金supported by the National Key R&D Program of China under Grant No.2021ZD0110400.
文摘Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational domains,etc.Typical reduced order modeling techniques accelerate the solution of the parametric PDEs by projecting them onto a linear trial manifold constructed in the ofline stage.These methods often need a predefined mesh as well as a series of precomputed solution snapshots,and may struggle to balance between the efficiency and accuracy due to the limitation of the linear ansatz.Utilizing the nonlinear representation of neural networks(NNs),we propose the Meta-Auto-Decoder(MAD)to construct a nonlinear trial manifold,whose best possible performance is measured theoretically by the decoder width.Based on the meta-learning concept,the trial manifold can be learned in a mesh-free and unsupervised way during the pre-training stage.Fast adaptation to new(possibly heterogeneous)PDE parameters is enabled by searching on this trial manifold,and optionally fine-tuning the trial manifold at the same time.Extensive numerical experiments show that the MAD method exhibits a faster convergence speed without losing the accuracy than other deep learning-based methods.
基金support of the Leverhulme Trust under the Leverhulme Trust Research Project Grant scheme(Grant No.RPG-2020-316)V.P-P.and R.G.M would like to thank the support from the Engineering and Physical Sciences Research Council(EPSRC)under the EPSRC DTP PhD scheme(Grant No.EP/T517914/1).
文摘Photonic computing has recently become an interesting paradigm for high-speed calculation of computing processes using light-matter interactions.Here,we propose and study an electromagnetic wave-based structure with the ability to calculate the solution of partial differential equations(PDEs)in the form of the Helmholtz wave equation,∇^(2)(x,y)T+k^(2)(x,y)=0,with k as the wavenumber.To do this,we make use of a network of interconnected waveguides filled with dielectric inserts.In so doing,it is shown how the proposed network can mimic the response of a network of T-circuit elements formed by two series and a parallel impedances,i.e.,the waveguide network effectively behaves as a metatronic network.An in-depth theoretical analysis of the proposed metatronic structure is presented,showing how the governing equation for the currents and impedances of the metatronic network resembles that of the finite difference representation of the Helmholtz wave equation.Different studies are then discussed including the solution of PDEs for Dirichlet and open boundary value problems,demonstrating how the proposed metatronic-based structure has the ability to calculate their solutions.
基金supported by NSFC Project(11301446,11271145)China Postdoctoral Science Foundation Grant(2013M531789)+3 种基金Specialized Research Fund for the Doctoral Program of Higher Education(2011440711009)Program for Changjiang Scholars and Innovative Research Team in University(IRT1179)Project of Scientific Research Fund of Hunan Provincial Science and Technology Department(2013RS4057)the Research Foundation of Hunan Provincial Education Department(13B116)
文摘We propose and analyze a spectral Jacobi-collocation approximation for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide a rigorous error analysis for the collection method, which shows that the errors of the approximate solution decay exponentially in L^∞ norm and weighted L^2-norm. The numerical examples are given to illustrate the theoretical results.
文摘In this paper we are looking forward to finding the approximate analytical solutions for fractional integro-differential equations by using Sumudu transform method and Hermite spectral collocation method.The fractional derivatives are described in the Caputo sense.The applications related to Sumudu transform method and Hermite spectral collocation method have been developed for differential equations to the extent of access to approximate analytical solutions of fractional integro-differential equations.
基金Supported by the National Natural Science Foundation of China (10671184)
文摘A lumped mass approximation scheme of a low order Crouzeix-Raviart type noncon- forming triangular finite element is proposed to a kind of nonlinear parabolic integro-differential equations. The L2 error estimate is derived on anisotropic meshes without referring to the traditional nonclassical elliptic projection.
基金Project supported by the National Natural Science Foundation of China(Nos.10971203,11271340,and 11101381)the Specialized Research Fund for the Doctoral Program of Higher Education(No.20094101110006)
文摘A highly efficient H1-Galerkin mixed finite element method (MFEM) is presented with linear triangular element for the parabolic integro-differential equation. Firstly, some new results about the integral estimation and asymptotic expansions are studied. Then, the superconvergence of order O(h^2) for both the original variable u in H1 (Ω) norm and the flux p = u in H(div, Ω) norm is derived through the interpolation post processing technique. Furthermore, with the help of the asymptotic expansions and a suitable auxiliary problem, the extrapolation solutions with accuracy O(h^3) are obtained for the above two variables. Finally, some numerical results are provided to confirm validity of the theoretical analysis and excellent performance of the proposed method.
基金supported by the National Nature Science Foundation of China (10671167)
文摘In this paper, the author discusses the multiple positive solutions for an infinite boundary value problem of first order impulsive singular integro-differential equations on the half line by means of the fixed point theorem of cone expansion and compression with norm type.
文摘Using the monotone iterative method and Monch Fixed point theorem, the existence of solutions and coupled minimal and maximal quasisolutions of initial value problems for mixed monotone second-order integro-differential equations in Banach spaces are studied. Some existence theorems of solutions and coupled minimal and maximal quasisolutions are obtained.
文摘In this paper, we prove an important existence and uniqueness theorem for a fractional order Fredholm – Volterra integro-differential equation with non-local and global boundary conditions by converting it to the corresponding well known Fredholm integral equation of second kind. The considered in this paper has been solved already numerically in [1].
文摘The current work models a weak(soft) interface between two elastic materials as containing a periodic array of micro-crazes. The boundary conditions on the interfacial micro-crazes are formulated in terms of a system of hypersingular integro-differential equations with unknown functions given by the displacement jumps across opposite faces of the micro-crazes. Once the displacement jumps are obtained by approximately solving the integro-differential equations, the effective stiffness of the micro-crazed interface can be readily computed. The effective stiffness is an important quantity needed for expressing the interfacial conditions in the spring-like macro-model of soft interfaces. Specific case studies are conducted to gain physical insights into how the effective stiffness of the interface may be influenced by the details of the interfacial micro-crazes.
基金the NSF of China(12171266,12171062)the NSF of Chongqing(CSTB2022NSCQ-JQX0004)。
文摘Let X be a complex Banach space and let B and C be two closed linear operators on X satisfying the condition D(B)?D(C),and let d∈L^(1)(R_(+))and 0≤β<α≤2.We characterize the well-posedness of the fractional integro-differential equations D^(α)u(t)+CD^(β)u(t)=Bu(t)+∫_(-∞)td(t-s)Bu(s)ds+f(t),(0≤t≤2π)on periodic Lebesgue-Bochner spaces L^(p)(T;X)and periodic Besov spaces B_(p,q)^(s)(T;X).
文摘In this paper, the author discusses the multiple positive solutions for an infinite boundary value problem of first order impulsive superlinear integro-differential equations on the half line by means of the fixed point theorem of cone expansion and compression with norm type.