This paper describes an approximating solution, based on Lagrange interpolation and spline functions, to treat functional integral equations of Fredholm type and Volterra type. This method extended to functional integ...This paper describes an approximating solution, based on Lagrange interpolation and spline functions, to treat functional integral equations of Fredholm type and Volterra type. This method extended to functional integral and integro-differential equations. For showing efficiency of the method we give some numerical examples.展开更多
It is well-known that interpolation by rational functions results in a more accurate approximation than the polynomials interpolation.However,classical rational interpolation has some deficiencies such as uncontrollab...It is well-known that interpolation by rational functions results in a more accurate approximation than the polynomials interpolation.However,classical rational interpolation has some deficiencies such as uncontrollable poles and low convergence order.In contrast with the classical rational interpolants,the generalized barycentric rational interpolants which depend linearly on the interpolated values,yield infinite smooth approximation with no poles in real numbers.In this paper,a numerical collocation approach,based on the generalized barycentric rational interpolation and Gaussian quadrature formula,was introduced to approximate the solution of Volterra-Fredholm integral equations.Three types of points in the solution domain are used as interpolation nodes.The obtained numerical results confirm that the barycentric rational interpolants are efficient tools for solving Volterra-Fredholm integral equations.Moreover,integral equations with Runge’s function as an exact solution,no oscillation occurrs in the obtained approximate solutions so that the Runge’s phenomenon is avoided.展开更多
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.展开更多
In this study, we prove the of existence of solutions of a convolution Volterra integral equation in the space of the Lebesgue integrable function on the set of positive real numbers and with the standard norm defined...In this study, we prove the of existence of solutions of a convolution Volterra integral equation in the space of the Lebesgue integrable function on the set of positive real numbers and with the standard norm defined on it. An operator P was assigned to the convolution integral operator which was later expressed in terms of the superposition operator and the nonlinear operator. Given a ball B<sub>r</sub> belonging to the space L it was established that the operator P maps the ball into itself. The Hausdorff measure of noncompactness was then applied by first proving that given a set M∈ B r the set is bounded, closed, convex and nondecreasing. Finally, the Darbo fixed point theorem was applied on the measure obtained from the set E belonging to M. From this application, it was observed that the conditions for the Darbo fixed point theorem was satisfied. This indicated the presence of at least a fixed point for the integral equation which thereby implying the existence of solutions for the integral equation.展开更多
The issues of solvability and construction of a solution of the Fredholm integral equation of the first kind are considered. It is done by immersing the original problem into solving an extremal problem in Hilbert spa...The issues of solvability and construction of a solution of the Fredholm integral equation of the first kind are considered. It is done by immersing the original problem into solving an extremal problem in Hilbert space. Necessary and sufficient conditions for the existence of a solution are obtained. A method of constructing a solution of the Fredholm integral equation of the first kind is developed. A constructive theory of solvability and construction of a solution to a boundary value problem of a linear integrodifferential equation with a distributed delay in control, generated by the Fredholm integral equation of the first kind, has been created.展开更多
By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. I...By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. In the numerical analysis, unknown displacement discontinuities are approximated with the products of the fundamental density functions and power series. The fundamental functions are chosen to express a two-dimensional interface crack rigorously. As illustrative examples, the stress intensity factors for two rectangular interface cracks are calculated for various spacing, crack shape and elastic constants. It is shown that the stress intensity factors decrease with the crack spacing.展开更多
In this work,a numerical scheme is constructed for solving nonlinear parabolictype partial-integro differential equations.The proposed numerical scheme is based on radial basis functions which are local in nature like...In this work,a numerical scheme is constructed for solving nonlinear parabolictype partial-integro differential equations.The proposed numerical scheme is based on radial basis functions which are local in nature like finite difference numerical schemes.The radial basis functions are used to approximate the derivatives involved and the integral is approximated by equal width integration rule.The resultant differentiation matrices are sparse in nature.After spatial approximation using RBF the partial integro-differential equations reduce to the system of ODEs.Then ODEs system can be solved by various types of ODE solvers.The proposed numerical scheme is tested and compared with other methods available in literature for different test problems.The stability and convergence of the present numerical scheme are discussed.展开更多
We present the existence of solution for a coupled system of fractional integro-differential equations. The differential operator is taken in the Caputo fractional sense. We combine the diagonalization method with Arz...We present the existence of solution for a coupled system of fractional integro-differential equations. The differential operator is taken in the Caputo fractional sense. We combine the diagonalization method with Arzela-Ascoli theorem to show a fixed point theorem of Schauder.展开更多
In this work, we consider the second order nonlinear integro-differential Equation (IDEs) of the Volterra-Fredholm type. One of the popular methods for solving Volterra or Fredholm type IDEs is the method of quadratur...In this work, we consider the second order nonlinear integro-differential Equation (IDEs) of the Volterra-Fredholm type. One of the popular methods for solving Volterra or Fredholm type IDEs is the method of quadrature while the problem of consideration is a linear problem. If IDEs are nonlinear or integral kernel is complicated, then quadrature rule is not most suitable;therefore, other types of methods are needed to develop. One of the suitable and effective method is homotopy analysis method (HAM) developed by Liao in 1992. To apply HAM, we firstly reduced the IDEs into nonlinear integral Equation (IEs) of Volterra-Fredholm type;then the standard HAM was applied. Gauss-Legendre quadrature formula was used for kernel integrations. Obtained system of algebraic equations was solved numerically. Moreover, numerical examples demonstrate the high accuracy of the proposed method. Comparisons with other methods are also provided. The results show that the proposed method is simple, effective and dominated other methods.展开更多
An entirely new framework is established for developing various single- and multi-step formulations for the numerical integration of ordinary differential equations. Besides polynomials, unconventional base-functions ...An entirely new framework is established for developing various single- and multi-step formulations for the numerical integration of ordinary differential equations. Besides polynomials, unconventional base-functions with trigonometric and exponential terms satisfying different conditions are employed to generate a number of formulations. Performances of the new schemes are tested against well-known numerical integrators for selected test cases with quite satisfactory results. Convergence and stability issues of the new formulations are not addressed as the treatment of these aspects requires a separate work. The general approach introduced herein opens a wide vista for producing virtually unlimited number of formulations.展开更多
It is proved that, for the nondivergence elliptic equations Σi,jn=1aijuxixj=f, if f belongs to the generalized Morrey spaces Lp, (w), then uxixj ∈ Lp, (w), where u is the W2,p-solution of the equations. In order to ...It is proved that, for the nondivergence elliptic equations Σi,jn=1aijuxixj=f, if f belongs to the generalized Morrey spaces Lp, (w), then uxixj ∈ Lp, (w), where u is the W2,p-solution of the equations. In order to obtain this, the author first establish the weighted boundedness for the commutators of some singular integral operators on Lp, (w).展开更多
Using integration by parts and Stokes' formula, the authors give a new definition of Hadamard principal value of higher order singular integrals with Bochner-Martinelli kernel on smooth closed orientable manifolds...Using integration by parts and Stokes' formula, the authors give a new definition of Hadamard principal value of higher order singular integrals with Bochner-Martinelli kernel on smooth closed orientable manifolds in C-n. The Plemelj formula and composite formula of higher order singular integral are obtained. Differential integral equations on smooth closed orientable manifolds are treated by using the composite formula.展开更多
In this paper, the difficulties on calculation in solving singular integral equations are overcome when the restriction of curve of integration to be a closed contour is cancelled. When the curve is an open arc and th...In this paper, the difficulties on calculation in solving singular integral equations are overcome when the restriction of curve of integration to be a closed contour is cancelled. When the curve is an open arc and the solutions for singular integral equations possess singularities of higher order, the solution and the solvable condition for characteristic equations as well as the generalized Noether theorem for complete equations are given.展开更多
A perturbation finite volume(PFV)method for the convective-diffusion integral equa- tion is developed in this paper.The PFV scheme is an upwind and mixed scheme using any higher-order interpolation and second-order in...A perturbation finite volume(PFV)method for the convective-diffusion integral equa- tion is developed in this paper.The PFV scheme is an upwind and mixed scheme using any higher-order interpolation and second-order integration approximations,with the least nodes similar to the standard three-point schemes,that is,the number of the nodes needed is equal to unity plus the face-number of the control volume.For instance,in the two-dimensional(2-D)case,only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized,respectively.The PFV scheme is applied on a number of 1-D linear and nonlinear problems,2-D and 3-D flow model equations.Comparing with other standard three-point schemes,the PFV scheme has much smaller numerical diffusion than the first-order upwind scheme(UDS).Its numerical accuracies are also higher than the second-order central scheme(CDS),the power-law scheme(PLS)and QUICK scheme.展开更多
In this paper, we study a boundary value problem of nonlinear fractional dif- ferential equations of order q (1 〈 q 〈 2) with non-separated integral boundary conditions. Some new existence and uniqueness results a...In this paper, we study a boundary value problem of nonlinear fractional dif- ferential equations of order q (1 〈 q 〈 2) with non-separated integral boundary conditions. Some new existence and uniqueness results are obtained by using some standard fixed point theorems and Leray-Schauder degree theory. Some illustrative examples are also presented. We extend previous results even in the integer case q = 2.展开更多
We consider a kind of scattering problem by a crack F that is buried in a bounded domain D, and we put a point source inside the domain D. This leads to a mixed boundary value problem to the Helmholtz equation in the ...We consider a kind of scattering problem by a crack F that is buried in a bounded domain D, and we put a point source inside the domain D. This leads to a mixed boundary value problem to the Helmholtz equation in the domain D with a crack Г. Both sides of the crack F are given Dirichlet-impedance boundary conditions, and different boundary condition (Dirichlet, Neumann or Impedance boundary condition) is set on the boundary of D. Applying potential theory, the problem can be reformulated as a system of boundary integral equations. We establish the existence and uniqueness of the solution to the system by using the Fredholm theory.展开更多
In this paper, we propose and analyze some schemes of the integral collocation formulation based on Legendre polynomials. We implement these formulae to solve numerically Riccati, Logistic and delay differential equat...In this paper, we propose and analyze some schemes of the integral collocation formulation based on Legendre polynomials. We implement these formulae to solve numerically Riccati, Logistic and delay differential equations with variable coefficients. The properties of the Legendre polynomials are used to reduce the proposed problems to the solution of non-linear system of algebraic equations using Newton iteration method. We give numerical results to satisfy the accuracy and the applicability of the proposed schemes.展开更多
In this article, we study the Volterra integral equations with two kinds of delay that are proportional delay and nonproportional delay. We mainly use Chebyshev spectral collocation method to analyze them. First, we u...In this article, we study the Volterra integral equations with two kinds of delay that are proportional delay and nonproportional delay. We mainly use Chebyshev spectral collocation method to analyze them. First, we use variable transformation to transform the equation into an new equation which is defined in [-1,1]. Then, with the help of Gronwall inequality and some other lemmas, we provide a rigorous error analysis for the proposed method, which shows that the numerical error decay exponentially in L~∞ and L_(ω~c)~2-norm. In the end, we give numerical test to confirm the conclusion.展开更多
We present a fast algorithm based on polynomial interpolation to approximate matrices arising from the discretization of second-kind integral equations where the kernel function is either smooth, non-oscillatory and p...We present a fast algorithm based on polynomial interpolation to approximate matrices arising from the discretization of second-kind integral equations where the kernel function is either smooth, non-oscillatory and possessing only a finite number of singularities or a product of such function with a highly oscillatory coefficient function. Contrast to wavelet-like approximations, ourapproximation matrix is not sparse. However, the approximation can be construced in O(n) operations and requires O(n) storage, where n is the number of quadrature points used in the discretization. Moreover, the matrix-vector multiplication cost is of order O(nlogn). Thus our scheme is well suitable for conjugate gradient type methods. Our numerical results indicate that the algorithm is very accurate and stable for high degree polynomial interpolation.展开更多
The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the bas...The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.展开更多
文摘This paper describes an approximating solution, based on Lagrange interpolation and spline functions, to treat functional integral equations of Fredholm type and Volterra type. This method extended to functional integral and integro-differential equations. For showing efficiency of the method we give some numerical examples.
文摘It is well-known that interpolation by rational functions results in a more accurate approximation than the polynomials interpolation.However,classical rational interpolation has some deficiencies such as uncontrollable poles and low convergence order.In contrast with the classical rational interpolants,the generalized barycentric rational interpolants which depend linearly on the interpolated values,yield infinite smooth approximation with no poles in real numbers.In this paper,a numerical collocation approach,based on the generalized barycentric rational interpolation and Gaussian quadrature formula,was introduced to approximate the solution of Volterra-Fredholm integral equations.Three types of points in the solution domain are used as interpolation nodes.The obtained numerical results confirm that the barycentric rational interpolants are efficient tools for solving Volterra-Fredholm integral equations.Moreover,integral equations with Runge’s function as an exact solution,no oscillation occurrs in the obtained approximate solutions so that the Runge’s phenomenon is avoided.
文摘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.
文摘In this study, we prove the of existence of solutions of a convolution Volterra integral equation in the space of the Lebesgue integrable function on the set of positive real numbers and with the standard norm defined on it. An operator P was assigned to the convolution integral operator which was later expressed in terms of the superposition operator and the nonlinear operator. Given a ball B<sub>r</sub> belonging to the space L it was established that the operator P maps the ball into itself. The Hausdorff measure of noncompactness was then applied by first proving that given a set M∈ B r the set is bounded, closed, convex and nondecreasing. Finally, the Darbo fixed point theorem was applied on the measure obtained from the set E belonging to M. From this application, it was observed that the conditions for the Darbo fixed point theorem was satisfied. This indicated the presence of at least a fixed point for the integral equation which thereby implying the existence of solutions for the integral equation.
文摘The issues of solvability and construction of a solution of the Fredholm integral equation of the first kind are considered. It is done by immersing the original problem into solving an extremal problem in Hilbert space. Necessary and sufficient conditions for the existence of a solution are obtained. A method of constructing a solution of the Fredholm integral equation of the first kind is developed. A constructive theory of solvability and construction of a solution to a boundary value problem of a linear integrodifferential equation with a distributed delay in control, generated by the Fredholm integral equation of the first kind, has been created.
基金supported by the National Natural Science Foundation of China (No. 10872213)
文摘By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. In the numerical analysis, unknown displacement discontinuities are approximated with the products of the fundamental density functions and power series. The fundamental functions are chosen to express a two-dimensional interface crack rigorously. As illustrative examples, the stress intensity factors for two rectangular interface cracks are calculated for various spacing, crack shape and elastic constants. It is shown that the stress intensity factors decrease with the crack spacing.
文摘In this work,a numerical scheme is constructed for solving nonlinear parabolictype partial-integro differential equations.The proposed numerical scheme is based on radial basis functions which are local in nature like finite difference numerical schemes.The radial basis functions are used to approximate the derivatives involved and the integral is approximated by equal width integration rule.The resultant differentiation matrices are sparse in nature.After spatial approximation using RBF the partial integro-differential equations reduce to the system of ODEs.Then ODEs system can be solved by various types of ODE solvers.The proposed numerical scheme is tested and compared with other methods available in literature for different test problems.The stability and convergence of the present numerical scheme are discussed.
文摘We present the existence of solution for a coupled system of fractional integro-differential equations. The differential operator is taken in the Caputo fractional sense. We combine the diagonalization method with Arzela-Ascoli theorem to show a fixed point theorem of Schauder.
文摘In this work, we consider the second order nonlinear integro-differential Equation (IDEs) of the Volterra-Fredholm type. One of the popular methods for solving Volterra or Fredholm type IDEs is the method of quadrature while the problem of consideration is a linear problem. If IDEs are nonlinear or integral kernel is complicated, then quadrature rule is not most suitable;therefore, other types of methods are needed to develop. One of the suitable and effective method is homotopy analysis method (HAM) developed by Liao in 1992. To apply HAM, we firstly reduced the IDEs into nonlinear integral Equation (IEs) of Volterra-Fredholm type;then the standard HAM was applied. Gauss-Legendre quadrature formula was used for kernel integrations. Obtained system of algebraic equations was solved numerically. Moreover, numerical examples demonstrate the high accuracy of the proposed method. Comparisons with other methods are also provided. The results show that the proposed method is simple, effective and dominated other methods.
文摘An entirely new framework is established for developing various single- and multi-step formulations for the numerical integration of ordinary differential equations. Besides polynomials, unconventional base-functions with trigonometric and exponential terms satisfying different conditions are employed to generate a number of formulations. Performances of the new schemes are tested against well-known numerical integrators for selected test cases with quite satisfactory results. Convergence and stability issues of the new formulations are not addressed as the treatment of these aspects requires a separate work. The general approach introduced herein opens a wide vista for producing virtually unlimited number of formulations.
文摘It is proved that, for the nondivergence elliptic equations Σi,jn=1aijuxixj=f, if f belongs to the generalized Morrey spaces Lp, (w), then uxixj ∈ Lp, (w), where u is the W2,p-solution of the equations. In order to obtain this, the author first establish the weighted boundedness for the commutators of some singular integral operators on Lp, (w).
基金the Bilateral Science and Technology Collaboration Program of Australia 1998 the Natural Science Foundation of China (No. 1
文摘Using integration by parts and Stokes' formula, the authors give a new definition of Hadamard principal value of higher order singular integrals with Bochner-Martinelli kernel on smooth closed orientable manifolds in C-n. The Plemelj formula and composite formula of higher order singular integral are obtained. Differential integral equations on smooth closed orientable manifolds are treated by using the composite formula.
基金Foundation item is supported by the NNSF of China(19971064)
文摘In this paper, the difficulties on calculation in solving singular integral equations are overcome when the restriction of curve of integration to be a closed contour is cancelled. When the curve is an open arc and the solutions for singular integral equations possess singularities of higher order, the solution and the solvable condition for characteristic equations as well as the generalized Noether theorem for complete equations are given.
基金The project supported by the National Natural Science Foundation of China(10272106,10372106)
文摘A perturbation finite volume(PFV)method for the convective-diffusion integral equa- tion is developed in this paper.The PFV scheme is an upwind and mixed scheme using any higher-order interpolation and second-order integration approximations,with the least nodes similar to the standard three-point schemes,that is,the number of the nodes needed is equal to unity plus the face-number of the control volume.For instance,in the two-dimensional(2-D)case,only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized,respectively.The PFV scheme is applied on a number of 1-D linear and nonlinear problems,2-D and 3-D flow model equations.Comparing with other standard three-point schemes,the PFV scheme has much smaller numerical diffusion than the first-order upwind scheme(UDS).Its numerical accuracies are also higher than the second-order central scheme(CDS),the power-law scheme(PLS)and QUICK scheme.
文摘In this paper, we study a boundary value problem of nonlinear fractional dif- ferential equations of order q (1 〈 q 〈 2) with non-separated integral boundary conditions. Some new existence and uniqueness results are obtained by using some standard fixed point theorems and Leray-Schauder degree theory. Some illustrative examples are also presented. We extend previous results even in the integer case q = 2.
基金supported by the grant from the National Natural Science Foundation of China(11301405)supported by the grants from the National Natural Science Foundation of China(11171127 and 10871080)
文摘We consider a kind of scattering problem by a crack F that is buried in a bounded domain D, and we put a point source inside the domain D. This leads to a mixed boundary value problem to the Helmholtz equation in the domain D with a crack Г. Both sides of the crack F are given Dirichlet-impedance boundary conditions, and different boundary condition (Dirichlet, Neumann or Impedance boundary condition) is set on the boundary of D. Applying potential theory, the problem can be reformulated as a system of boundary integral equations. We establish the existence and uniqueness of the solution to the system by using the Fredholm theory.
文摘In this paper, we propose and analyze some schemes of the integral collocation formulation based on Legendre polynomials. We implement these formulae to solve numerically Riccati, Logistic and delay differential equations with variable coefficients. The properties of the Legendre polynomials are used to reduce the proposed problems to the solution of non-linear system of algebraic equations using Newton iteration method. We give numerical results to satisfy the accuracy and the applicability of the proposed schemes.
基金supported by National Science Foundation of China(11671157,11626074)Hanshan Normal Uninversity projects(LF201404,Z16027)
文摘In this article, we study the Volterra integral equations with two kinds of delay that are proportional delay and nonproportional delay. We mainly use Chebyshev spectral collocation method to analyze them. First, we use variable transformation to transform the equation into an new equation which is defined in [-1,1]. Then, with the help of Gronwall inequality and some other lemmas, we provide a rigorous error analysis for the proposed method, which shows that the numerical error decay exponentially in L~∞ and L_(ω~c)~2-norm. In the end, we give numerical test to confirm the conclusion.
基金Research supported in part by Hong Kong Research Grant Council grats no.CUHK178/83E
文摘We present a fast algorithm based on polynomial interpolation to approximate matrices arising from the discretization of second-kind integral equations where the kernel function is either smooth, non-oscillatory and possessing only a finite number of singularities or a product of such function with a highly oscillatory coefficient function. Contrast to wavelet-like approximations, ourapproximation matrix is not sparse. However, the approximation can be construced in O(n) operations and requires O(n) storage, where n is the number of quadrature points used in the discretization. Moreover, the matrix-vector multiplication cost is of order O(nlogn). Thus our scheme is well suitable for conjugate gradient type methods. Our numerical results indicate that the algorithm is very accurate and stable for high degree polynomial interpolation.
文摘The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.