A domain decomposition solver for acoustic scattering by elastic objects in layered media
Introduction
We develop an efficient numerical method for computing time-harmonic acoustic wave scattering by an elastic object in three-dimensional layered media. It is capable of modeling the scattering of sonar signals by undersea targets located in or near the seabed in littoral environments with plane or rippled interface between water and sediment. One application for such problems is the detection of hazardous or/and lost objects buried in sediment. For this purpose it is essential to have a numerical approximation which can accurately predict the scattered field by such targets. Our model problem in this paper is the scattering by a solid aluminum cylinder buried in sediment. A similar experimental test setup was used in [30] to perform measurements. This is an exterior problem which has to be truncated into a bounded domain for a finite element discretization; see Fig. 1.
We consider scattering problems in a frequency range in which the size of the computational domain is of the order of 100 wavelengths. The discretization has to have sufficiently many nodes per wavelength. In the considered frequency range the phase (pollution) error of the solution with a second-order accurate discretization is large unless tens of nodes are used per wavelength [20]. This would lead to huge systems of linear equations which are computationally intractable. Many ways have been proposed to reduce the phase error; see [20], for example. We employ a generalization of the approach proposed in [13] in which standard bilinear finite elements are used for two-dimensional Helmholtz problems together with a modified quadrature rule. This leads to fourth-order phase accuracy and a higher order approximation on orthogonal uniform meshes for a constant speed of sound. Actually it coincides with the nine-point fourth-order compact finite difference scheme [29] in the homogeneous medium. A related fourth-order finite difference method for the Helmholtz equation was considered in [12]. The resulting matrix has exactly the same structure and complexity as the one obtained by using the standard Gauss quadrature rule. As shown in Section 5.1 it drastically increases the performance and capability of our finite element modeling. The scheme with eight points per wavelength sufficiently well achieves the resolution we need and thus leads to a linear system with much smaller numbers of unknowns.
The resulting system of linear equations is too large to be solved using a direct method. We propose an iterative solution procedure with a decomposition preconditioner based on a domain embedding approach [15], [16], [17], [21], [22], [24], [26]. In our solution procedure the computational domain is decomposed into the near field subdomain which encloses an elastic target (the interior box in Fig. 1) and the far-field subdomain (the rest of the rectangle computational domain Π). In the far-field subdomain the discretization is based on an orthogonal uniform mesh which is locally adapted to the interface between the water and sediment. We use a multiplicative domain decomposition, i.e. we develop a preconditioning method for the resulting linear equation that partitions the matrix into the decoupled blocks. We will use a separable preconditioner based on the perfectly vertically layered media in our solution procedure for the far-field subdomain. The far-field preconditioner is based on a fast direct solver [16], [33]. Since the media is vertically layered with the wavy interface, our preconditioner coincides with the system matrix except the rows corresponding to unknowns near the interfaces. Thus, we can reduce iterations on a small sparse subspace as has been shown in [22], [24], for example, and Section 4.3. This procedure reduces GMRES iterates onto a sparse subspace which leads to much reduced storage and computational requirements. The GMRES iteration is evaluated by the partial solution method based on a cyclic reduction type fast direct solver. This reduction makes our iterative method extremely efficient and enables us to solve systems with billions of unknowns as our numerical examples demonstrate. With multigrid preconditioners like in [1], [9], [10] and with nonorthogonal grids in [28] such reduction cannot be made and the solution of the model problem would require much more memory and computation. With the separation-of-variables preconditioner in [31] the reduction of the iterations onto a sparse subspace can be made with the model problem.
A similar domain decomposition approach to the one considered here was used for two-dimensional problems in [4], [5], [23]. Specifically, the approach in [23] is extended to the three-dimensional case. Such extension is nontrivial and involves the use of higher order discretization in the layered media and efficient direct solver for the 3D separable preconditioner. Another preconditioning technique for scattering problems in layered media without an object has been considered in [1], [10], [28], [31], for example, and with an object in [19].
The outline of the paper is the following. We begin by defining a model problem in Section 2. For the finite element discretization a weak formulation, meshes, a domain decomposition, the forming of linear systems and a modified quadrature rule are described in Section 3. The iteration with a domain decomposition preconditioner is consider in Section 4. Also, the solution procedures for the far-field and near-field preconditioners are given and how iterations can be performed in a sparse subspace is described. Numerical experiments are presented in Section 5. They include an accuracy study for the discretizations and an efficiency study of the iterative solution. The paper ends with the conclusions in Section 6.
Section snippets
Model problem
We truncate the exterior domain into a rectangular parallelepiped domain Π shown in Fig. 1. We pose an absorbing boundary condition on the truncation boundary ∂Π. For the water and sediment fluid region (Π⧹Ω) we use the Helmholtz equation to model the time-harmonic pressure variations p. In the elastic object (Ω) the displacement u is described by the time-harmonic solution for the linear elastic wave equation. Coupling these equations leads to the partial differential equation model:
Finite element discretization
In this section, we describe our finite element discretization of the model problem.
Preconditioned iteration
The domain decomposition and discretization in Section 3.2 leads the system of linear equations (9) to have the block formwhere the first block row corresponds to the interior of the near-field subdomain Ω1 and the second block row corresponds to the far-field subdomain Ω2 and the interface between the subdomains formed by four faces of a rectangular parallelepiped.
Based on the block structure in (10), we define an upper triangular domain decomposition preconditioner
Grid refinement analysis
In this section, we present a grid refinement analysis for the second-order and fourth-order discretization. Our analysis is conducted for an acoustic scattering problem in the layered media with the straight sediment and no target. The test frequency is 3.75 kHz. The computational domain is [0 m,12 m] × [−0.4 m, 0.4 m] × [−0.8 m, 4 m] and the point source is located at [0.4 m, 0 m, 3.6 m]. The plane interface between the water and sediment is at x3 = 0 m. The density of water is 1000 kg/m3 and the speed of
Conclusions
We have developed an efficient numerical method for computing time-harmonic acoustic scattering by an elastic object in three-dimensional layered media. The infinite domain is truncated to a rectangular parallelepiped and a second-order absorbing boundary condition is posed on the truncation boundary which usually leads to sufficient accuracy for practical purposes. Our discretization uses a modified trilinear finite element discretization with fourth-order phase accuracy. We studied the
Acknowledgments
The authors like to thank Dr. David Burnett, Dr. Quyen Huynh, and Dr. Joseph Lopes for helpful discussions on littoral acoustic scattering problems. The research was supported by the Office of Naval Research Grant N00014-06-1-0067, the National Science Foundation Grant DMS-0610661, and the Academy of Finland Grant #207089.
References (36)
- et al.
An algebraic multigrid based shifted-Laplacian preconditioner for the Helmholtz equation
J. Comput. Phys.
(2007) - et al.
Numerical solution of the nonlinear Helmholtz equation using nonorthogonal expansions
J. Comput. Phys.
(2005) - et al.
Modified integration rules for reducing dispersion in finite element methods
Comput. Methods Appl. Mech. Eng.
(2004) - et al.
Fictitious domain methods for the numerical solution of two-dimensional scattering problems
J. Comput. Phys.
(1998) - et al.
Preconditioned iterative methods on sparse subspaces
Appl. Math. Lett.
(2006) - et al.
A fast iterative solver for scattering by elastic objects in layered media
Appl. Numer. Math.
(2007) - et al.
Separation-of-variables as a preconditioner for an iterative Helmholtz solver
Appl. Numer. Math.
(2003) - et al.
Second-order absorbing boundary conditions for the wave equation: a solution for the corner problem
SIAM J. Numer. Anal.
(1990) Fast Poisson solvers for problems with sparsity
Math. Comp.
(1978)- et al.
Material surface design to counter electromagnetic interrogation of targets
SIAM J. Appl. Math.
(2006)
Determination of interrogating frequencies to maximize electromagnetic backscatter from objects with material coatings
Commun. Comput. Phys.
A triangulation algorithm for fast elliptic solvers based on domain imbedding
SIAM J. Numer. Anal.
Algorithm 832: UMFPACK V4.3—an unsymmetric-pattern multifrontal method
ACM Trans. Math. Software
A multigrid method enhanced by Krylov subspace iteration for discrete Helmhotz equations
SIAM J. Sci. Comput.
A novel multigrid based preconditioner for heterogeneous Helmholtz problems
SIAM J. Sci. Comput.
Compact fourth-order finite difference schemes for Helmholtz equation with high wave numbers
J. Comput. Math.
Cited by (24)
A fast preconditioned iterative method for the electromagnetic scattering by multiple cavities with high wave numbers
2019, Journal of Computational PhysicsCitation Excerpt :Li and Qiao [20] employed a fast Krylov subspace iterative method to solve the scattering problem by the single cavity. Other efficient preconditioners, such as domain decomposition preconditioner [21] and incomplete LU preconditioner [22] were also applied for solving scattering problems. This paper is concentrated on the scattering by multiple cavities with high wave numbers.
A composite preconditioner for the electromagnetic scattering from a large cavity
2011, Journal of Computational PhysicsCitation Excerpt :For such case, a large part of rows of the coefficient matrix of the discrete Helmholtz problem and the preconditioner proposed in this paper coincide. Therefore, preconditioned Krylov subspace methods can be carried out on sparse subspaces [11,40–42,44]. In particular, we propose a sparse preconditioned COCG solver combined with the new preconditioner for a model problem.
Numerical solution of electromagnetic scattering from a large partly covered cavity
2011, Journal of Computational and Applied MathematicsCitation Excerpt :For example, a shifted-Laplacian preconditioner [12,13] was used for heterogeneous Helmholtz problems. A domain decomposition method was proposed in [14] for acoustic scattering by elastic objects in layered media. The rest of the paper is organized as follows.
High-performance modeling acoustic and elastic waves using the parallel Dichotomy Algorithm
2011, Journal of Computational PhysicsCitation Excerpt :Although early results on elastic wavefield modeling have been obtained long ago [42–44] however, they were rather qualitative because of a large step of the space mesh h = 1/5λ ÷ 1/2λ. Considerably increased computer performance and also development of multiprocessor computer systems have made it possible to increase the calculation accuracy [45–47]. However, in spite of available theoretical estimates of the dependence of solution accuracy on mesh step [8], the problem of practical choosing a space step of meshes is still urgent.
A damping preconditioner for time-harmonic wave equations in fluid and elastic material
2009, Journal of Computational PhysicsIs Pollution Effect of Finite Difference Schemes Avoidable for Multi-Dimensional Helmholtz Equations with High Wave Numbers?
2017, Communications in Computational Physics