Isogeometric locally-conformal perfectly matched layer for time-harmonic acoustics

Published on Oct 1, 2021in Computer Methods in Applied Mechanics and Engineering5.763
· DOI :10.1016/J.CMA.2021.113925
Yongzhen Mi (Institute of High Performance Computing Singapore), Xiang Yu (Institute of High Performance Computing Singapore)
Abstract null null null This paper presents an isogeometric formulation of the perfectly matched layer (PML) for time-harmonic acoustic simulations. The new formulation is a null generalization null null of the conventional locally-conformal PML, in which the null NURBS null patch supporting the PML domain is transformed from real space to complex space. This is achieved by complex coordinate stretching, based on a stretching vector field indicating the directions in which incident sound waves are absorbed. To solve the difficulty that stretching vectors cannot be defined for control points not lying on the patch, an interpolation method is proposed in which the control-point stretching vectors are calculated by interpolating those at a set of integration points. An automatically matched layer (AML) is also formulated in the isogeometric context. Due to its form-invariant nature, the null null null Helmholtz equation null null null which governs time-harmonic acoustics is kept unchanged in the complex space, except that the null Laplacian operator null null null and the sound pressure should be transformed to their complex forms. The performance of the isogeometric PML formulation is discussed through several numerical examples, spanning one to three dimensions and involving both sound radiation and sound scattering problems. It is found that the proposed method presents superior null computational accuracy , high geometric adaptivity, and good robustness against challenging null null geometric features . The geometry-preserving ability inherent in the isogeometric framework brings extra benefits by eliminating geometric errors that are unavoidable in the conventional PML. Meanwhile, these null null null properties null are not sensitive to the location of the sound source or the depth of the PML domain. The implementation of the isogeometric PML only requires the replacement of real control-point coordinates by complex ones, making it possible to incorporate the proposed method to various existing numerical routines.
#1Hadrien Beriot (Siemens)H-Index: 9
#2Axel Modave (IRIA: French Institute for Research in Computer Science and Automation)H-Index: 8
We address the efficient finite element solution of exterior acoustic problems with truncated computational domains surrounded by perfectly matched layers (PMLs). The PML is a popular non-reflecting technique, which combines accuracy, computational efficiency and geometric flexibility. Unfortunately, the effective implementation of the PML for generally-shaped convex domains is tricky because of the geometric parameters that are required to define the PML medium. In this work, a comprehensive im...
3 CitationsSource
#1Xiangyu You (HUST: Huazhong University of Science and Technology)H-Index: 6
#2Yingbin Chai (HUST: Huazhong University of Science and Technology)H-Index: 14
Last. Wei Li (HUST: Huazhong University of Science and Technology)H-Index: 99
view all 3 authors...
Abstract It is known that for exterior acoustics the domain-based finite element method always encounters the accuracy loss for high wavenumbers and the difficulty in dealing with the Sommerfeld radiation condition. In this work, by using the polynomial and radial bases with specific node-selection schemes, the point interpolation method with the generalized gradient smoothing technique is employed to improve the accuracy for acoustic computations. The appropriately softened acoustical stiffness...
3 CitationsSource
#1Antoine Vermeil de Conchard (École centrale de Nantes)H-Index: 1
#2Huina Mao (KTH: Royal Institute of Technology)H-Index: 5
Last. Romain Rumpler (KTH: Royal Institute of Technology)H-Index: 8
view all 3 authors...
Abstract Effective treatment of unbounded domains using artificial truncating boundaries are essential in numerical simulation, e.g. using the Finite Element Method (FEM). Among these, Perfectly Matched Layers (PML) have proved to be particularly efficient and flexible. However, an efficient handling of frequency sweeps is not trivial with such absorbing layers since the formulation inherently contains coupled space- and frequency-dependent terms. Using the FEM, this may imply generating system ...
3 CitationsSource
#1Shaowei Wu (Chongqing Jiaotong University)H-Index: 1
#2Yang Xiang (WUT: Wuhan University of Technology)H-Index: 1
Abstract To improve the computational accuracy and eliminate the element matching in the traditional method based on the finite and infinite elements, the coupling of the meshfree weak-form method with the variable order infinite acoustic wave envelope element (WEE) is studied in this paper and achieved by forming continuous acoustic shape function on an artificial interface, by which the infinite domain is divided into two subdomains where the meshfree and WEE methods are applied respectively. ...
9 CitationsSource
#1A. El Kacimi (Cadi Ayyad University)H-Index: 6
#2Omar Laghrouche (Heriot-Watt University)H-Index: 19
Last. Jon Trevelyan (Durham University)H-Index: 18
view all 6 authors...
Abstract The aim of this paper is to accurately solve short wave scattering problems governed by the Helmholtz equation using the Bernstein–Bezier Finite Element method (BBFEM), combined with a conformal perfectly matched layer (PML). Enhanced PMLs, where curved geometries are represented by means of the blending map method of Gordon and Hall, are numerically investigated. In particular, the performance of radial and elliptical shaped PMLs, with a parabolic absorption function, are assessed and ...
6 CitationsSource
#1A. F. M. ter Elst (University of Auckland)H-Index: 19
#2El Maati Ouhabaz (University of Bordeaux)H-Index: 20
Abstract We prove Poisson upper bounds for the heat kernel of the Dirichlet-to-Neumann operator with variable Holder coefficients when the underlying domain is bounded and has a C 1 + κ -boundary for some κ > 0 . We also prove a number of other results such as gradient estimates for heat kernels and Green functions G of elliptic operators with possibly complex-valued coefficients. We establish Holder continuity of ∇ x ∇ y G up to the boundary. These results are used to prove L p -estimates for c...
8 CitationsSource
#1Tahsin Khajah (UT Tyler: University of Texas at Tyler)H-Index: 5
#2Vianey Villamizar (BYU: Brigham Young University)H-Index: 8
Abstract This work is concerned with a unique combination of high order local absorbing boundary conditions (ABC) with a general curvilinear Finite Element Method (FEM) and its implementation in Isogeometric Analysis (IGA) for time-harmonic acoustic waves. The ABC employed were recently devised by Villamizar et al. (2017) . They are derived from exact Farfield Expansions representations of the outgoing waves in the exterior of the regions enclosed by the artificial boundary. As a consequence, th...
11 CitationsSource
#1Jie Yang (HEU: Harbin Engineering University)H-Index: 1
#2Xinyu Zhang (HEU: Harbin Engineering University)H-Index: 1
Last. Wenping Zhang (HEU: Harbin Engineering University)H-Index: 1
view all 4 authors...
This paper focuses on the development of absorbing boundary conditions and their implementation in the smoothed particle hydrodynamics (SPH) method for wave propagation problems. A compact perfectly matched layer (C-PML) approach is formulated for transient acoustic problems in an infinite space using SPH models with a computational domain of finite dimensions. The proposed approach is based on the concept of fictitious physical damping acting within the additional perfectly matched layers (PMLs...
3 CitationsSource
#1Yongzhen Mi (SJTU: Shanghai Jiao Tong University)H-Index: 4
#2Hui Zheng (SJTU: Shanghai Jiao Tong University)H-Index: 14
Abstract This paper presents an interpolation method for coupling non-conforming Non-Uniform Rational B-Splines (NURBS) patches in isogeometric analysis of vibro-acoustic systems. This coupling method is based on a master–slave formulation, in which a smooth interpolating function is first constructed to represent the sound pressure field of the master side and then is sampled to determine the discrete sound pressure of the slave side. The coupling constraints are studied by the interpolation ma...
7 CitationsSource
#1Yongzhen Mi (SJTU: Shanghai Jiao Tong University)H-Index: 4
#2Hui Zheng (SJTU: Shanghai Jiao Tong University)H-Index: 14
Last. Yuanyi Huang (SAIC: Science Applications International Corporation)H-Index: 1
view all 4 authors...
This paper develops a weak formulation for steady-state dynamic analysis of vibro-acoustic systems in the isogeometric framework. The Non-Uniform Rational B-Splines (NURBS) polynomials are used as ...
4 CitationsSource
Cited By0