{"title":"体积积分方程中奇异积分的有效数值计算","authors":"Cedric Münger;Kristof Cools","doi":"10.1109/JMMCT.2022.3185531","DOIUrl":null,"url":null,"abstract":"We present a method for the numerical evaluation of 6D and 5D singular integrals appearing in Volume Integral Equations. It is an extension of the Sauter-Schwab/Taylor-Duffy strategy for singular triangle-triangle interaction integrals to singular tetrahedron-tetrahedron and triangle-tetrahedron interaction integrals. The general advantages of these kind of quadrature strategy is that they allow the use of different kinds of kernel and basis functions. They also work on curvilinear domains. They are all based on relative coordinates tranformation and splitting the integration domain into subdomains for which quadrature rules can be constructed. We show how to build these tensor-product quadrature rules in 6D and 5D and further show how to improve their efficiency by using quadrature rules defined over 2D, 3D and 4D simplices. Compared to the existing approach, which computes the integral over the subdomains as a sequence of 1D integrations, significant speedup can be achieved. The accuracy and convergence properties of the method are demonstrated by numerical experiments for 5D and 6D singular integrals. Additionally, we applied the new quadrature approach to the triangle-triangle interaction integrals appearing in Surface Integral Equations.","PeriodicalId":52176,"journal":{"name":"IEEE Journal on Multiscale and Multiphysics Computational Techniques","volume":null,"pages":null},"PeriodicalIF":1.8000,"publicationDate":"2022-06-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Efficient Numerical Evaluation of Singular Integrals in Volume Integral Equations\",\"authors\":\"Cedric Münger;Kristof Cools\",\"doi\":\"10.1109/JMMCT.2022.3185531\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"We present a method for the numerical evaluation of 6D and 5D singular integrals appearing in Volume Integral Equations. It is an extension of the Sauter-Schwab/Taylor-Duffy strategy for singular triangle-triangle interaction integrals to singular tetrahedron-tetrahedron and triangle-tetrahedron interaction integrals. The general advantages of these kind of quadrature strategy is that they allow the use of different kinds of kernel and basis functions. They also work on curvilinear domains. They are all based on relative coordinates tranformation and splitting the integration domain into subdomains for which quadrature rules can be constructed. We show how to build these tensor-product quadrature rules in 6D and 5D and further show how to improve their efficiency by using quadrature rules defined over 2D, 3D and 4D simplices. Compared to the existing approach, which computes the integral over the subdomains as a sequence of 1D integrations, significant speedup can be achieved. The accuracy and convergence properties of the method are demonstrated by numerical experiments for 5D and 6D singular integrals. Additionally, we applied the new quadrature approach to the triangle-triangle interaction integrals appearing in Surface Integral Equations.\",\"PeriodicalId\":52176,\"journal\":{\"name\":\"IEEE Journal on Multiscale and Multiphysics Computational Techniques\",\"volume\":null,\"pages\":null},\"PeriodicalIF\":1.8000,\"publicationDate\":\"2022-06-23\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"IEEE Journal on Multiscale and Multiphysics Computational Techniques\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://ieeexplore.ieee.org/document/9804769/\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q3\",\"JCRName\":\"ENGINEERING, ELECTRICAL & ELECTRONIC\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"IEEE Journal on Multiscale and Multiphysics Computational Techniques","FirstCategoryId":"1085","ListUrlMain":"https://ieeexplore.ieee.org/document/9804769/","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"ENGINEERING, ELECTRICAL & ELECTRONIC","Score":null,"Total":0}
Efficient Numerical Evaluation of Singular Integrals in Volume Integral Equations
We present a method for the numerical evaluation of 6D and 5D singular integrals appearing in Volume Integral Equations. It is an extension of the Sauter-Schwab/Taylor-Duffy strategy for singular triangle-triangle interaction integrals to singular tetrahedron-tetrahedron and triangle-tetrahedron interaction integrals. The general advantages of these kind of quadrature strategy is that they allow the use of different kinds of kernel and basis functions. They also work on curvilinear domains. They are all based on relative coordinates tranformation and splitting the integration domain into subdomains for which quadrature rules can be constructed. We show how to build these tensor-product quadrature rules in 6D and 5D and further show how to improve their efficiency by using quadrature rules defined over 2D, 3D and 4D simplices. Compared to the existing approach, which computes the integral over the subdomains as a sequence of 1D integrations, significant speedup can be achieved. The accuracy and convergence properties of the method are demonstrated by numerical experiments for 5D and 6D singular integrals. Additionally, we applied the new quadrature approach to the triangle-triangle interaction integrals appearing in Surface Integral Equations.