Shape Aware Quadrature

JOURNAL OF COMPUTATIONAL PHYSICS

卷:374

Article

Thiagarajan, Vaidyanathan¹  Shapiro, Vadim¹ 

[1] Univ Wisconsin, 1513 Univ Ave, Madison, WI 53706 USA

关键词: Numerical integration;    Moment fitting equations;    Shape sensitivity analysis;    Topological sensitivity analysis;    Divergence theorem;    Octree/quadtree;   

DOI  :  10.1016/j.jcp.2018.05.024

来源: Elsevier

【 摘 要 】

With Shape Aware Quadratures (SAQ), for a given set of quadrature nodes, order, and domain of integration, the quadrature weights are obtained by solving a system of suitable moment fitting equations in least square sense. The moments in the moment equations are approximated over a simplified domain (Omega(0)) that is a reasonable approximation of the original domain (Omega) that are then corrected for the deviation of the shape of Omega(0) from Omega via shape correction factors. This idea was already successfully utilized in the Adaptively Weighted Numerical Integration (AW) method [1 3], where moments were computed over simplified but homeomorphic domain and then corrected using first-order shape sensitivity, allowing efficient and accurate integration of integrable functions over arbitrary domains. With SAQ more general shape correction factors can be derived based on a variety of sensitivity analysis techniques such as shape sensitivity (SSA) and topological sensitivity (TSA). Using the right kind/order of shape correction factors for moment approximation enables the resulting quadrature (from the moment fitting equations) to efficiently adapt to the shape of the original domain even in the presence of thousands of small features. We derive shape correction factors based on SSA and TSA of moment integrals. We also demonstrate the efficacy of SAQ in integrating bivariate/trivariate polynomials over several 2D/3D domains in the presence of small features. (C) 2018 Elsevier Inc. All rights reserved.

【 授权许可】

Free