首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
This work is concerned with the numerical solution of a nonlinear weakly singular Volterra integral equation. Owing to the singular behavior of the solution near the origin, the global convergence order of product integration and collocation methods is not optimal. In order to recover the optimal orders a hybrid collocation method is used which combines a non-polynomial approximation on the first subinterval followed by piecewise polynomial collocation on a graded mesh. Some numerical examples are presented which illustrate the theoretical results and the performance of the method. A comparison is made with the standard graded collocation method.  相似文献   

2.
Summary. To solve 1D linear integral equations on bounded intervals with nonsmooth input functions and solutions, we have recently proposed a quite general procedure, that is essentially based on the introduction of a nonlinear smoothing change of variable into the integral equation and on the approximation of the transformed solution by global algebraic polynomials. In particular, the new procedure has been applied to weakly singular equations of the second kind and to solve the generalized air foil equation for an airfoil with a flap. In these cases we have obtained arbitrarily high orders of convergence through the solution of very-well conditioned linear systems. In this paper, to enlarge the domain of applicability of our technique, we show how the above procedure can be successfully used also to solve the classical Symm's equation on a piecewise smooth curve. The collocation method we propose, applied to the transformed equation and based on Chebyshev polynomials of the first kind, has shown to be stable and convergent. A comparison with some recent numerical methods using splines or trigonometric polynomials shows that our method is highly competitive. Received October 1, 1998 / Revised version received September 27, 1999 / Published online June 21, 2000  相似文献   

3.
In this paper, a multi-parameter error resolution technique is introduced and applied to the collocation method for Volterra integral equations. By using this technique, an approximation of higher accuracy is obtained by using a multi-processor in parallel. Additionally, a correction scheme for approximation of higher accuracy and a global superconvergence result are presented.  相似文献   

4.
We consider the problem of scattering of a time-harmonic acoustic incident plane wave by a sound soft convex polygon. For standard boundary or finite element methods, with a piecewise polynomial approximation space, the computational cost required to achieve a prescribed level of accuracy grows linearly with respect to the frequency of the incident wave. Recently Chandler–Wilde and Langdon proposed a novel Galerkin boundary element method for this problem for which, by incorporating the products of plane wave basis functions with piecewise polynomials supported on a graded mesh into the approximation space, they were able to demonstrate that the number of degrees of freedom required to achieve a prescribed level of accuracy grows only logarithmically with respect to the frequency. Here we propose a related collocation method, using the same approximation space, for which we demonstrate via numerical experiments a convergence rate identical to that achieved with the Galerkin scheme, but with a substantially reduced computational cost.  相似文献   

5.
Summary. A residual-based a posteriori error estimate for boundary integral equations on surfaces is derived in this paper. A localisation argument involves a Lipschitz partition of unity such as nodal basis functions known from finite element methods. The abstract estimate does not use any property of the discrete solution, but simplifies for the Galerkin discretisation of Symm's integral equation if piecewise constants belong to the test space. The estimate suggests an isotropic adaptive algorithm for automatic mesh-refinement. An alternative motivation from a two-level error estimate is possible but then requires a saturation assumption. The efficiency of an anisotropic version is discussed and supported by numerical experiments. Received November 29, 1999 / Revised version received August 10, 2000 / Published online May 30, 2001  相似文献   

6.
The modified regularized long wave (MRLW) equation is solved numerically by collocation method using cubic B-splines finite element. A linear stability analysis of the scheme is shown to be marginally stable. Three invariants of motion are evaluated to determine the conservation properties of the algorithm, also the numerical scheme leads to accurate and efficient results. Moreover, interaction of two and three solitary waves are studied through computer simulation and the development of the Maxwellian initial condition into solitary waves is also shown.  相似文献   

7.
This paper is concerned with obtaining the approximate solution for VolterraHammerstein integral equation with a regular kernel. We choose the Gauss points associated with the Legendre weight function ω(x) = 1 as the collocation points. The Legendre collocation discretization is proposed for Volterra-Hammerstein integral equation. We provide an error analysis which justifies that the errors of approximate solution decay exponentially in L~2 norm and L~∞ norm. We give two numerical examples in order to illustrate the validity of the proposed Legendre spectral collocation method.  相似文献   

8.
This paper is devoted to the approximate solution of the classical first-kind boundary integral equation with logarithmic kernel (Symm's equation) on a closed polygonal boundary in ℝ2. We propose a fully discrete method with a trial space of trigonometric polynomials, combined with a trapezoidal rule approximation of the integrals. Before discretization the equation is transformed using a nonlinear (mesh grading) parametrization of the boundary curve which has the effect of smoothing out the singularities at the corners and yields fast convergence of the approximate solutions. The convergence results are illustrated with some numerical examples.  相似文献   

9.
We discuss the application of spline collocation methods to a certain class of weakly singular Volterra integral equations. It will be shown that, by a special choice of the collocation parameters, superconvergence properties can be obtained if the exact solution satisfies certain conditions. This is in contrast with the theory of collocation methods for Abel type equations. Several numerical examples are given which illustrate the theoretical results.  相似文献   

10.
We discuss the convergence properties of spline collocation and iterated collocation methods for a weakly singular Volterra integral equation associated with certain heat conduction problems. This work completes the previous studies of numerical methods for this type of equations with noncompact kernel. In particular, a global convergence result is obtained and it is shown that discrete superconvergence can be achieved with the iterated collocation if the exact solution belongs to some appropriate spaces. Some numerical examples illustrate the theoretical results.  相似文献   

11.
Summary The collocation method is a popular method for the approximate solution of boundary integral equations, but typically does not achieve the high order of convergence reached by the Galerkin method in appropriate negative norms. In this paper a quadrature-based method for improving upon the collocation method is proposed, and developed in detail for a particular case. That case involves operators with even symbol (such as the logarithmic potential) operating on 1-periodic functions; a smooth-spline trial space of odd degree, with constant mesh spacingh=1/n; and a quadrature rule with 2n points (where ann-point quadrature rule would be equivalent to the collocation method). In this setting it is shown that a special quadrature rule (which depends on the degree of the splines and the order of the operator) can yield a maximum order of convergence two powers ofh higher than the collocation method.  相似文献   

12.
Midpoint collocation for Cauchy singular integral equations   总被引:1,自引:0,他引:1  
Summary A Cauchy singular integral equation on a smooth closed curve may be solved numerically using continuous piecewise linear functions and collocation at the midpoints of the underlying grid. Even if the grid is non-uniform, suboptimal rates of convergence are proved using a discrete maximum principle for a modified form of the collocation equations. The same techniques prove negative norm estimates when midpoint collocation is used to determine piecewise constant approximations to the solution of first kind equations with the logarithmic potential.This work was supported by the Australian Research Council through the program grant Numerical analysis for integrals, integral equations and boundary value problems  相似文献   

13.
We establish the uniform convergence of a collocation method for solving a class of singular integral equations. This method uses the Jacobi polynomials {P n (, ) } as basis elements and the zeros of a Chebyshev polynomial of the first kind as collocation points. Uniform convergence is shown to hold under the weak assumption that the kernel and the right-hand side are Hölder-continous functions. Convergence rates are also given.  相似文献   

14.
15.
Summary A method which combines quadrature with trigonometric interpolation is proposed for singular integral equations on closed curves. For the case of the circle, the present method is shown to be equivalent to the trigonometric -collocation method together with numerical quadrature for the compact term, and is shown to be stable inL 2 provided the operatorA is invertible inL 2. The results are extended to arbitraryC curves, to give a complete error analysis in the scale of Sobolev spacesH s . In the final section the case of a non-invertible operatorA is considered.  相似文献   

16.
17.
In this paper, we investigate the convergence behavior of a Runge–Kutta type modified Landweber method for nonlinear ill-posed operator equations. In order to improve the stability and convergence of the Landweber iteration, a 2-stage Gauss-type Runge–Kutta method is applied to the continuous analogy of the modified Landweber method, to give a new modified Landweber method, called R–K type modified Landweber method. Under some appropriate conditions, we prove the convergence of the proposed method. We conclude with a numerical example confirming the theoretical results, including comparisons to the modified Landweber iteration.  相似文献   

18.
This paper is concerned with the applicability of maximum defect polynomial (Galerkin) spline approximation methods with graded meshes to Wiener-Hopf operators with matrix-valued piecewise continuous generating function defined on R. For this, an algebra of sequences is introduced, which contains the approximating sequences we are interested in. There is a direct relationship between the stability of the approximation method for a given operator and invertibility of the corresponding sequence in this algebra. Exploring this relationship, the methods of essentialization, localization and identification of the local algebras are used in order to derive stability criteria for the approximation sequences.Supported by grant Praxis XXI/BD/4501/94 from FCT.Partly supported by FCT/BMFT grant 423.  相似文献   

19.
20.
A boundary integral representation for the conformal mappingfrom a simply connected domain onto a given canonical regionis presented. The conformal map can be constructed via the solutionof a first-kind integral equation with logarithmic kernel: Symm'sintegral equation. A novel method for solving this equationis presented and a partial error analysis, including convergencerates, is given. The analysis is also relevant to the solutionof Symm's equation on a smooth, open contour. A superconvergenceresult for the conformal mapping is proved. Now in the Department of Mathematics and Computer Science, Universtyof Leicester, Leicester LE1 7RH, UK.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号