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.展开更多
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.展开更多
For domains composed by balls in C^n, this paper studies the boundary behaviour of Cauchy type integrals with discrete holomorphic kernels and the corresponding linear singular integral equation on each piece of smoot...For domains composed by balls in C^n, this paper studies the boundary behaviour of Cauchy type integrals with discrete holomorphic kernels and the corresponding linear singular integral equation on each piece of smooth lower dimensional edges on the boundary of the domain.展开更多
The solutions of the nonlinear singular integral equation , t 6 L, are considered, where L is a closed contour in the complex plane, b ≠ 0 is a constant and f(t) is a polynomial. It is an extension of the results obt...The solutions of the nonlinear singular integral equation , t 6 L, are considered, where L is a closed contour in the complex plane, b ≠ 0 is a constant and f(t) is a polynomial. It is an extension of the results obtained in [1] when f(t) is a constant. Certain special cases are illustrated.展开更多
In this paper, we develop a fractional cyclic integral and a Routh equation for fractional Lagrange system defined in terms of fractional Caputo derivatives. The fractional Hamilton principle and the fractional Lagran...In this paper, we develop a fractional cyclic integral and a Routh equation for fractional Lagrange system defined in terms of fractional Caputo derivatives. The fractional Hamilton principle and the fractional Lagrange equations of the system are obtained under a combined Caputo derivative. Furthermore, the fractional cyclic integrals based on the Lagrange equations are studied and the associated Routh equations of the system are presented. Finally, two examples are given to show the applications of the results.展开更多
To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of...To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of the memorial dynamics. Stability criteria of the scheme for an advection equation in certain conditions are derived mathematically. The computations for the advection equation have been conducted with its RT scheme. It is shown that the accuracy of the scheme is much higher than that of the leapfrog (LF) difference scheme.展开更多
The elastodynamic problems of piezoelectric hollow cylinders and spheres under radial deformation can be transformed into a second kind Volterra integral equation about a function with respect to time, which greatly s...The elastodynamic problems of piezoelectric hollow cylinders and spheres under radial deformation can be transformed into a second kind Volterra integral equation about a function with respect to time, which greatly simplifies the solving procedure for such elastodynamic problems. Meanwhile, it becomes very important to find a way to solve the second kind Volterra integral equation effectively and quickly. By using an interpolation function to approximate the unknown function, two new recursive formulae were derived, based on which numerical solution can be obtained step by step. The present method can provide accurate numerical results efficiently. It is also very stable for long time calculating.展开更多
In this paper, we present a brief survey on the updated theory of backward stochas-tic Volterra integral equations (BSVIEs, for short). BSVIEs are a natural generalization of backward stochastic diff erential equati...In this paper, we present a brief survey on the updated theory of backward stochas-tic Volterra integral equations (BSVIEs, for short). BSVIEs are a natural generalization of backward stochastic diff erential equations (BSDEs, for short). Some interesting motivations of studying BSVIEs are recalled. With proper solution concepts, it is possible to establish the corresponding well-posedness for BSVIEs. We also survey various comparison theorems for solutions to BSVIEs.展开更多
文摘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.
文摘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.
基金Project supported by the National Science Foundation of China (10271097)
文摘For domains composed by balls in C^n, this paper studies the boundary behaviour of Cauchy type integrals with discrete holomorphic kernels and the corresponding linear singular integral equation on each piece of smooth lower dimensional edges on the boundary of the domain.
文摘The solutions of the nonlinear singular integral equation , t 6 L, are considered, where L is a closed contour in the complex plane, b ≠ 0 is a constant and f(t) is a polynomial. It is an extension of the results obtained in [1] when f(t) is a constant. Certain special cases are illustrated.
基金Project supported by the National Natural Science Foundations of China(Grant Nos.11272287 and 11472247)the Program for Changjiang Scholars and Innovative Research Team in University(PCSIRT)(Grant No.IRT13097)
文摘In this paper, we develop a fractional cyclic integral and a Routh equation for fractional Lagrange system defined in terms of fractional Caputo derivatives. The fractional Hamilton principle and the fractional Lagrange equations of the system are obtained under a combined Caputo derivative. Furthermore, the fractional cyclic integrals based on the Lagrange equations are studied and the associated Routh equations of the system are presented. Finally, two examples are given to show the applications of the results.
基金The project supported by the National Key Program for Developing Basic Sciences (G1999043408 and G1998040901-1)the National Natural Sciences Foundation of China (40175024 and 40035010)
文摘To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of the memorial dynamics. Stability criteria of the scheme for an advection equation in certain conditions are derived mathematically. The computations for the advection equation have been conducted with its RT scheme. It is shown that the accuracy of the scheme is much higher than that of the leapfrog (LF) difference scheme.
文摘The elastodynamic problems of piezoelectric hollow cylinders and spheres under radial deformation can be transformed into a second kind Volterra integral equation about a function with respect to time, which greatly simplifies the solving procedure for such elastodynamic problems. Meanwhile, it becomes very important to find a way to solve the second kind Volterra integral equation effectively and quickly. By using an interpolation function to approximate the unknown function, two new recursive formulae were derived, based on which numerical solution can be obtained step by step. The present method can provide accurate numerical results efficiently. It is also very stable for long time calculating.
文摘In this paper, we present a brief survey on the updated theory of backward stochas-tic Volterra integral equations (BSVIEs, for short). BSVIEs are a natural generalization of backward stochastic diff erential equations (BSDEs, for short). Some interesting motivations of studying BSVIEs are recalled. With proper solution concepts, it is possible to establish the corresponding well-posedness for BSVIEs. We also survey various comparison theorems for solutions to BSVIEs.