1 Introduction
As one of the most fundamental numerical methods, the boundary integral equation (BIE) method [25]
has been extensively developed for numerical solutions of partial differential equations problems with various structures including bounded closed-surface
[8, 4, 12], open screens [9, 10, 11, 27, 35], period or non-period infinite surface [13], and so on. The BIE method has a feature of discretization of domains with lower dimensionality, and it is also a feasible method for the numerics of high frequency scattering problems. For large-scale problems with high-frequencies or three-dimensional complicated geometries, such iterative algorithms [19] as the Krylov-subspace linear algebra solver GMRES, together with adequate acceleration techniques [5, 8, 28], are generally required for fast solving the resulting linear system whose coefficient matrix is dense. The efficiency of the GMRES iteration is highly related to the spectral features of the coefficient matrix of the linear system [34] and therefore, appropriate preconditioning, such as the analytical preconditioning based on the Calderón formulas [7, 14] and the algebraic preconditioning strategies [6], are usually employed. However, only a few theoretical properties of the Calderón formulas (also called the Calderón relation in this work) for acoustic/electromagnetic closed-surface problems [7, 14], two-dimensional acoustic open-surface problems [9, 27] and elastic closed-surface problems with the standard traction operator [11, 13], have been studied in open literatures. We also refer to [21, 22, 23, 24] and the references therein for the study of inverses of integral operators on disks and the corresponding preconditioning associated with boundary element Galerkin discretizations.This work is devoted to studying the theoretical properties of the Calderón formulas related to the two-dimensional problems of elastic scattering by closed- or open-surfaces which have many significant applications in science and engineering [31, 36], including geophysics, non-destructive testing of solids materials, mining and energy production, etc. A fundamental purpose of utilizing the Calderón formulas is to construct BIEs, for example, the second-kind Fredholm integral equations, with the highly favorable spectral properties that the eigenvalues of the BIEs are bounded away from zero and infinity. One can refer to the methodologies discussed in [7] for the acoustic case and those in [15] for the electromagnetic case. Although for the acoustic and elastodynamic problems, the Calderón formulas are indeed the composition of the single-layer integral operator and the hyper-singular integral operator , the extension of the theoretical analysis on the Calderón formulas in acoustics to that on the elastodynamic cases, however, encounters additional challenges. More precisely, for the smooth closed-surface case, the acoustic Calderón formula reads where represents the transpose of the double-layer boundary integral operator, and is compact in appropriate Sobolev spaces. This fact ensures that the acoustic Calderón formula is of the second-kind Fredholm type. However, the corresponding operator in the elastic case is not compact, see for example [1, 2]. In addition, the highly singular character of the associated integral kernel in elastodynamic hyper-singular operator is much more complicated than that in the acoustic case.
For the closed-surface case, by applying the polynomial compactness of the statistic elastic Neumann-Poincaré double-layer operator and its transpose , it has been shown in [11, 12] that the elastodynamic Calderón formula involving the standard traction operator (2.1) is exactly a second-kind Fredholm operator (see [38] for two dimensional poroelastic case) whose eigenvalues are bounded away from zero and infinity with accumulation points being dependent on the elastic Lamé parameters. In addition, on the basis of another special choice of the traction operator (see Lemma 3.1(ii)), the corresponding elastodynamic Calderón formulas are Fredholm operators of second-kind in both two and three dimensions as well. In this paper, a generalized traction operator [20, 19] related to the generalized Betti’s formula [3, 26] will be considered, and the general results to be presented in Theorem 3.4 indicate that the generalized elastic Calderón formula is a Fredholm operator of second-kind except for the above two special forms of traction operator.
Unfortunately, the properties of the closed-surface Calderón formula become invalid in the open-surface case. As being verified in [27, 32], the two-dimensional acoustic composite operator , which is not a second-kind Fredholm operator, takes a local singularity like where denotes the distance between the node being considered and the nearby endpoint of the open-arc. The composite operator in the elastic case suffers from analogous character and therefore, it can not be treated in order to obtain favorable features in the classical Sobolev spaces. Instead of discussing , in light of the singular character of the solutions of the single-layer and the hyper-singular BIEs for solving the corresponding acoustic open-surface scattering problems [18], a novel acoustic version weighted Calderón formula is proposed in [9, 10], and it can be written into a sum of an invertible operator and a compact operator [9, 27]. According to what we have known, the similar theoretical analysis for elastic problems, including acoustic and electromagnetic problems in three dimensions, still remains unavailable in existing literatures. As being numerically demonstrated in [11, 12] for elastic open-surface scattering problems, the elastic version formula leads to a significant decreasing on the GMRES iterations compared to the un-preconditioned one for a given residual tolerance. As a significant complement to the above numerical observation, we present in the current work a rigorous theory on the two-dimensional elastic Calderón formula which actually can be viewed as a Fredholm integral operator of second kind and a compact perturbation of a bounded and invertible operator. In addition, the accumulation point of the spectrum of the invertible operator is the same as that of the eigenvalues of the elastic Calderón formula in the two-dimensional closed-surface case, see Remark 4.8.
The remainder of this paper is organized as follows. Section 2 introduces the generalized traction operator together with the elastic single-layer and hyper-singular boundary integral operators. Section 3 investigates the spectral properties of the elastic Calderón formula in the closed-surface case and presents some regularized formulations of the elastic hyper-singular boundary integral operators. The elastic Calderón formula in the open-surface case is studied in Section 4: instead of the original elastic boundary integral operators, weighted single-layer and hyper-singular boundary integral operators under certain edge singularity circumstance of potentials are introduced in Section 4.1; in terms of the analysis results of the elastic Calderón formula on a special straight open-arc and the singularity decompositions of the integration kernels, the spectral properties of the generalized elastic Calderón formula in the universal open-surface case is analyzed in Section 4.2 and 4.3. A conclusion is finally given in Section 5.
2 Preliminaries
Let be a smooth closed-surface or open-surface in . Denote by the Lamé parameters and let be the mass density of a linear isotropic and homogeneous elastic medium. Denote by the frequency. For elastic problems, the standard traction operator on the boundary is defined as
(2.1) |
in which is the unit outer normal to the boundary ,
is the corresponding tangential vector,
denotes the normal derivative and . To produce the generalized elastic Caldrón relations, we consider a modified traction operator [20] defined as follows:(2.2) |
where . Obviously, holds if .
In this work, we are interested in the theoretical properties of the Calderón formulas, i.e., the composite operator , for elastic closed-surface and open-surface scattering problems, where denote the elastodynamic single-layer and hyper-singular boundary integral operators, respectively, in the form of
(2.3) |
and
(2.4) |
Here, we denote by
the fundamental displacement tensor of the time-harmonic Navier equation in
. That is(2.5) |
where denotes the Lamé operator given by
and is the identity matrix. It is known that admits the form [26]
where represents the fundamental solution of the Helmholtz equation in :
(2.6) |
with wave numbers , being the Hankel function of the first kind of order zero, and . The wave numbers are called the wave number of the compressional and shear waves, respectively, where
Remark 2.1.
The modified traction operator can be formulated alternatively as
(2.7) |
where the Günter derivative operator is given by
This modified traction operator is equivalent to another generalized form discussed in [19]
It follows that when .
3 Calderón relation: closed-surface
Let be a smooth closed boundary. As shown in [19], the Calderón identities
(3.1) |
hold where the single-layer boundary integral operator and hyper-singular boundary integral operator are defined in (2.3) and (2.4), respectively, and the operators defined as
are called the double-layer boundary integral operator and transpose of double-layer boundary integral operator, respectively. It is known [25] that and are linear bounded operators.
Lemma 3.1.
(i). If , then and are compact. Here, , is the identity operator and is a constant given by
(3.2) |
(ii). If , then the operators themselves are compact.
The property of the integral operators for general are given in the following lemma.
Lemma 3.2.
For all , and are compact. Here, and is a constant that depends on the Lamé parameters:
(3.3) |
Proof.
Let be the corresponding elastic double-layer operator in the zero-frequency case which is given by [25]
where
with
Denote by the elements of a matrix. Then direct calculation using (2.2) yields that
(3.4) | |||||
and
(3.5) | |||||
Then gives
where
Define the integral operators in the sense of Cauchy principle value as follows:
Due to the property that [16], we conclude that the operator is compact. For the operator , it has been proved in [1, Proposition 3.1] that is compact. Then it follows that is compact. Then the compactness of results immediately from the fact that is compact, since has a at-most weakly singular kernel, and
The compactness of can be proved analogously. ∎
Remark 3.3.
As stated in the following theorem, the generalized elastic Calderón relation for closed-surface problem is a direct corollary of Lemma 3.2. Note that under the special values (3.6) of , it holds that .
Theorem 3.4.
If
(3.6) |
the compositions and are compact. Otherwise, the compositions and are Fredholm operators of second-kind both of which consists of a non-empty sequence of eigenvalues converging to (Fig. 1).
![]() |
![]() |
Unfortunately, the generalized elastic Calderón relation for the closed-surface problem does not hold for the open-surface case any more, see Section 4.1. Before discussing the corresponding Calderón relation for the open-surface problem, we provide with some useful formulations in the closed-surface case which can be extended to the open-surface case for certain weighted integral operators introduced in Section 4.1. Let be the corresponding elastic hyper-singular operator in the zero-frequency case which is given by
(3.7) |
Lemma 3.5.
The hyper-singular operator can be reformulated as
(3.8) |
where the integral operator is defined as
where
Proof.
The proof of this lemma follows from the same steps as [25, Lemma 2.2.3] and is omitted here. ∎
Remark 3.6.
Lemma 3.7.
The hyper-singular boundary integral operator can be expressed alternatively as
(3.9) | |||||
Proof.
This result is a generalized form of [37, Theorem 6.2] and the proof is also omitted. ∎
4 Calderón relation: open-surface
In this section, we study the elastic Calderón relation in the open-surface case. However, unlike the closed-surface case, an appropriate functional setting for the investigation of the composite operator defined on an open-arc seems to be non-existing. As demonstrated in [27, Appendix B], even for the simplest open-surface–straight arc , the acoustic composite operator maps the constant function with value into a function possessing edge singularity which does not belong to , where denotes the distance between and the corresponding end point for any in a neighbourhood of each end point. This difficulty also appears in the elastic open-surface case since both the acoustic and elastic integral operators defined on the straight arc contain similar singular kernels.
4.1 Weighted integral operators
Instead of studying the operators , we discuss the weighted forms resulting from certain regularity of potentials. Let be a smooth open arc in and the unbounded domain is fulfilled with a linear isotropic and homogeneous elastic medium. Then the time-harmonic problem of elastic scattering by an open-surface can be modeled by the Navier equation
(4.1) |
together with the boundary conditions
(4.2) |
and the Kupradze radiation condition at infinity [26]. Here, denotes the displacement field.
Proof.
It is known that [33, 35] the solutions of the Dirichlet and Neumann elastic problems of scattering by an open-arc admit the representations in forms of single- and double-layer potentials, respectively, i.e.,
(4.3) |
and
(4.4) |
respectively. Then the Dirichlet and Neumann problems reduce to the boundary integral equations
(4.5) |
Definition 4.2.
An operator between two Sobolev spaces is called bicontinuous if it is continuous and invertible. As a corollary, the inverse is also continuous. Assume that where is a smooth boundary of a bounded domain in . We denote by the space of all satisfying .
Lemma 4.3.
The operators and are bicontinuous under the assumption that
(4.6) |
Proof.
The bicontinuity of the operator holds analogously to [33]. The assumption (4.6) of gives . By the integration kernel of together with the singularity decomposition (A.2), it can be deduced that the coefficient of the weakly-singular part in the term
is non-zero. Following the proof of the solvability of the hyper-singular operator [35, (3.6)] results into the bicontinuity of the operator . ∎
Assumption 4.1.
In the rest of this work, we always assume the value of to be such that (4.6) and additionally, hold.
Denote by a non-negative smooth function for to represent the distance between and the corresponding end point for any in a neighbourhood of each end point. Assuming that the right-hand sides in (4.5) are both infinitely differentiable, it is known [18] that the density functions in (4.5) can be expressed in the forms
(4.7) |
where denotes a smooth function that reproduces the asymptotic as . It implies that is infinitely differentiable up to the endpoints, and the new solutions are smooth up to the end points of . Taking into account the solution singularities (4.7), we obtain the new boundary integral equations
(4.8) |
where the weighted integral operators are defined as
A similar regularized formulation of can be obtained from (3.9) since the weight function in is smooth boundary-vanishing, see (4.19).
Without loss of generality, suppose that the boundary can be parameterized by means of a smooth vector function satisfying . Here the prime denotes the derivative with respect to . Choosing the smooth weighting function as yields
Then utilizing the changes of variables leads us to
where the operator is given by
The parameterized form corresponding to the integral operator can be deduced in a similar manner.
Note that each smooth dependence function can be extended to be a -periodic and even function. To study the properties of the parameterized operators , we define the following Sobolev spaces:
Definition 4.4.
For , the Sobolev space is defined as the completion of space of infinitely differentiable -periodic and even functions defined in the real line with respect to the norm
where denotes the coefficients in the cosine expansion of :
4.2 Calderón relation: straight arc
Let be, specially, the straight arc which means that for and . In this subsection, we consider the operators on the straight arc at zero frequency that are denoted by , i.e.,
(4.9) |
and
(4.10) |
where
The expression of follows immediately from Lemma 3.5 thanks to the smooth boundary-vanishing weight .
Comments
There are no comments yet.