3D isogeometric indirect BEM solution based on virtual surface sources on the boundaries of Helmholtz acoustic problems

Ahmed Mostafa Shaaban,Jon Trevelyan,Timon Rabczuk

DOI: https://doi.org/10.1007/s00366-023-01933-5

IF: 8.083

2024-02-08

Engineering With Computers

Abstract:A solution for 3D Helmholtz acoustic problems is introduced based on an indirect boundary element method (indirect BEM) coupled with isogeometric analysis (IGA). The novelty of this work arises from using virtual surface sources placed directly on the scatterer boundaries, producing robust results. These virtual surface sources are discretized by the same Non-Uniform Rational B-Splines (NURBS) approximating the scatterer CAD model. This allows modeling of general irregular geometries. The proposed solution has the same features of BEM approaches, which do not need any domain discretization or truncation boundaries at the far-field. It shows an additional merit by arranging the linear system of equations directly depending on a single coefficient matrix, consuming less computational time compared to other BEM methods. A Greville abscissae collocation scheme is proposed with offsets at -continuities. This collocation scheme allows for easy evaluation for both free-terms and normals at the collocation points. The performance of the proposed solution is discussed on 3D numerical exterior problems and compared against other BEM methods. Then, the practical interior muffler problem with internal extended thin tubes is studied and the obtained results are compared against other numerical methods in addition to the available experimental data, showing the capability of the proposed solution in handling thin-walled geometries.

engineering, mechanical,computer science, interdisciplinary applications

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to use the indirect boundary element method (indirect BEM) combined with isogeometric analysis (IGA) to solve Helmholtz acoustic problems in three - dimensional space. Specifically, the paper proposes a new method. By placing virtual surface sources on the boundary of the scatterer and using the same non - uniform rational B - spline (NURBS) as the scatterer CAD model for discretization, complex geometric shapes can be modeled. This method not only avoids the singularity problem in traditional methods, but also improves computational efficiency and reduces computing time.

### Main contributions of the paper:
1. **Use of virtual surface sources**: The paper proposes to directly place virtual surface sources on the boundary of the scatterer. These sources are discretized by NURBS to adapt to complex and irregular geometric shapes.
2. **Combination of isogeometric analysis (IGA)**: By using the same NURBS basis functions as the scatterer CAD model, the combination of the boundary element method and isogeometric analysis is achieved, which improves the accuracy and computational efficiency of the model.
3. **Improved collocation scheme**: The Greville collocation scheme with offset is introduced. This scheme is offset at C0 continuity, which simplifies the calculation of free terms and normal vectors and at the same time avoids the singularity problem.
4. **Direct construction of linear equations**: By directly relying on a single coefficient matrix to construct linear equations, the computing time is reduced and it is more efficient compared to other BEM methods.

### Specific problems solved:
- **Modeling of complex geometric shapes**: Traditional numerical methods have difficulties in dealing with complex geometric shapes, while the IGAinBEM method can effectively handle these complex shapes.
- **Singularity problem**: By placing virtual surface sources on the boundary, the singularity problem in traditional methods is avoided and the stability of the calculation is improved.
- **Computational efficiency**: By directly constructing linear equations, the computing time is reduced and the computational efficiency is improved.

### Application scenarios:
- **Communication systems**: Used for optimizing antenna design and signal propagation analysis.
- **Ultrasonic devices**: Used for ultrasonic imaging and the design of ultrasonic treatment devices.
- **Noise control**: Used for noise control and noise reduction techniques in architectural design.
- **Biomedical applications**: Used for the study of sound wave propagation in biological tissues.
- **Seismology**: Used for the simulation and analysis of seismic wave propagation.
- **Damage control**: Used for structural health monitoring and damage assessment.

In conclusion, the paper proposes an innovative method. By combining the indirect boundary element method and isogeometric analysis, it solves the problems of complex geometric modeling and singularity in three - dimensional Helmholtz acoustic problems and has broad application prospects.