Numerical quadrature for high-dimensional singular integrals over parallelotopes
نویسندگان
چکیده
We introduce and analyze a family of algorithms for an efficient numerical approximation of integrals of the form I = ∫ C(1) ∫ C(2) F (x, y, y−x)dydx where C, C are d-dimensional parallelotopes (i.e. affine images of d-hypercubes) and F has a singularity at y − x = 0. Such integrals appear in Galerkin discretization of integral operators in R. We construct a family of quadrature rules QN with N function evaluations for a class of integrands F which may have algebraic singularities at y − x = 0 and are Gevrey-δ regular for y − x 6= 0. The main tool is a regularizing coordinate transformation, simultaneously simplifying the singular support and the domain of integration. For the full tensor product variant of the suggested quadrature family we prove that QN achieves the exponential convergence rate O(exp(−rN )) with the exponent γ = 1/(2dδ + 1). In the special case of a singularity of the form ‖y − x‖ with real α we prove that the improved convergence rate of γ = 1/(2dδ) is achieved if a certain modified one-dimensional Gauss-Jacobi quadrature rule is used in the singular direction. We give numerical results for various types of the quadrature rules, in particular based on tensor product rules, standard (Smolyak), optimized and adaptive sparse grid quadratures and Sobol’ sequences.
منابع مشابه
A General New Algorithm for Regulaization of Singular Integrals in Three-Dimensional Boundary Elemnts
In this paper an algorithm is presented for the regularization of singular integrals with any degrees of singularity, which may be employed in all three-dimensional problems analyzed by Boundary Elements. The integrals in Boundary Integrals Equations are inherently singular. For example, one can mention the integrals confronted in potential problems to evaluate the flow or the gradient of the f...
متن کاملA General New Algorithm for Regulaization of Singular Integrals in Three-Dimensional Boundary Elemnts
In this paper an algorithm is presented for the regularization of singular integrals with any degrees of singularity, which may be employed in all three-dimensional problems analyzed by Boundary Elements. The integrals in Boundary Integrals Equations are inherently singular. For example, one can mention the integrals confronted in potential problems to evaluate the flow or the gradient of the f...
متن کاملExponential convergence of Gauss-Jacobi quadratures for singular integrals over high dimensional simplices
Galerkin discretizations of integral operators in R d require the evaluation of integrals R S (1) R S (2) f (x, y) dydx where S (1) , S (2) are d-dimensional simplices and f has a singularity at x = y. In [3] we constructed a family of hp-quadrature rules Q N with N function evaluations for a class of integrands f allowing for algebraic singularities at x = y, possibly non-integrable with respe...
متن کاملNearly Singular Integrals in 3D Stokes Flow
A straightforward method is presented for computing three-dimensional Stokes flow, due to forces on a surface, with high accuracy at points near the surface. The flow quantities arewritten as boundary integrals using the free-spaceGreen’s function. To evaluate the integrals near the boundary, the singular kernels are regularized and a simple quadrature is applied in coordinate charts. High orde...
متن کاملGauss-Jacobi-type quadrature rules for fractional directional integrals
Fractional directional integrals are the extensions of the Riemann-Liouville fractional integrals from oneto multi-dimensional spaces and play an important role in extending the fractional differentiation to diverse applications. In numerical evaluation of these integrals, the weakly singular kernels often fail the conventional quadrature rules such as Newton-Cotes and Gauss-Legendre rules. It ...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Computers & Mathematics with Applications
دوره 66 شماره
صفحات -
تاریخ انتشار 2013