1. Introduction
In this work, we use Clifford analysis and, for the case
, quaternionic analysis in the study of inhomogeneous Moisil-Teodorescu systems defined on bounded or unbounded domains of
,
. For higher dimensions, our main goal is to provide an explicit general solution to the
n-dimensional generalization of the equation
under certain conditions imposed on the data
over bounded Lipschitz domains. Thereafter, restricting ourselves to the three-dimensional case, we will analyze the div-curl system without boundary conditions on different classes of exterior domains. When the normal component of the vector field is also known, one speaks of a Neumann problem; one of the first papers to address this boundary-value problem was [
1] using the fundamental theorem of vector calculus. Later, the results were extended to exterior domains in [
2] under the condition that
and
decay faster that
as
.
Historically, one of the first works on the application of quaternionic analysis for elliptic systems in unbounded domains was the article [
3]. In that work, weighted Banach spaces
,
were employed in order to guarantee a good behavior at infinity of the Teodorescu transform and the Cauchy operator, which are some classical operators in quaternionic analysis. Unfortunately, the main disadvantage of considering the Teodorescu transform on unbounded domains in the usual form (that is, using the same Cauchy kernel
as in the bounded case) is that boundedness is lost in the classical integrable spaces
. A different approach was presented in the articles [
4,
5,
6]. This approach relies on the use of a perturbation of the Teodorescu transform through the addition of a monogenic term whose singularity lies outside of the unbounded domain under consideration. More precisely, the studies in those papers employed modified Cauchy kernels. Unlike those works, we will employ the usual kernel in the present manuscript. Moreover, in spite of the fact that the integrability of the integral operators is not improved, the operators exhibit a good asymptotic behavior at infinity for our purposes. To that end, we will employ results from [
4,
5,
6] with modified Cauchy kernels. Indeed, they are still valid for the present case, only the integrability ranges will be different.
One of the main results presented in this work (Theorem 5) establishes that a weak solution of the div-curl system in exterior domains star-shaped with respect to (w.r.t.) infinity has the form
where
is harmonic in
. In fact, it is this property together with the asymptotic condition
which allow us to construct a general solution with no boundary conditions. The non-uniqueness of the solution is clear. Indeed, if we add the gradient of a harmonic function in
to Equation (
1), then it will have the same divergence and rotational over all
. Moreover, the solution Equation (
1) admits the following Helmholtz-type decomposition in exterior domains (Corollary 3):
Here,
with
defined in Equation (
2). In addition, we obtained a second Helmholtz-type decomposition for the class of strong local Lipschitz exterior domains in terms of the layer potentials (Corollary 2):
where
is as in Equation (
4),
and
satisfies
where
defined in Equation (
2). Moreover,
Note that in the solenoidal part of both the Helmholtz-type decompositions Equations (
3) and (
6), the following operator appears:
Moreover, it coincides with the Biot–Savart operator defined over
. In particular, when
represents the current density, we obtain the Biot–Savart law in electromagnetism. This important Biot–Savart operator is also part of a strategic decomposition of the Teodorescu transform Equations (
52) and (
53), in which much of our Clifford analysis is based. More precisely, the operator Equation (
9) coincides with the vector part of the Teodorescu transform
applied to
, that is,
Interestingly, the non-uniqueness of these general solutions Equations (
3) and (
6) allows us to transform the associated Neumann BVP for the div-curl system into a Neuman BVP for the Laplacian, with results available in the literature on its existence, regularity and uniqueness. This gives rise to the main result of
Section 6, Theorem 7.
It is worth pointing out that each of the operators appearing in the expression of the general solution Equation (
1) is important an operator from quaternionic analysis. In the present work, we will frequently employ an important decomposition of the Teodorescu operator used in [
7,
8,
9] for bounded domains of
. In turn, the radial operator appearing in the last term of Equation (
1) was defined in [
10,
11] as a generalization to exterior domains of an important family of radial operators. We will talk briefly about this operator at the beginning of
Section 5. In this work, we will express the general solution following a quaternionic approach. More precisely, we will embed the div-curl system in the algebra of quaternions and then project this quaternionic solution to a purely vector one without affecting the system.
The outline of the paper is as follows. In
Section 2, we present the notation with basic theory of Clifford analysis as well as some facts about the regularity of the domain. In
Section 3, a general weak solution of the generalization to
n-dimensions of the equation
is provided using an appropriate embedding argument to the Clifford algebra
. In
Section 4, we present the construction of hyper-conjugate harmonic pairs in unbounded domains in terms of certain layer potentials and give explicit formulas for a solution of the div-curl system without boundary conditions for exterior domains satisfying the strong local Lipschitz condition. In
Section 5, we derive another explicit solution of the div-curl system in exterior domains now under the geometric condition that
is star-shaped w.r.t. infinity. The construction of this second solution relies on the properties of a family of radial integral operators restricted to a family of harmonic functions with good behavior at infinity. In
Section 6, we analyze in detail the regularity and asymptotic behavior of the solution of the div-curl system Equation (
1). Thereafter, we adjust this solution to construct a weak solution of the Neumann BVP of the div-curl system in exterior domains. Due to the easy handling of the radial operator that appears in the construction of the general solution (Equation (
1)) of the exterior div-curl system, we will thoroughly analyze and adjust this general solution instead of the other general solution Equation (
3) found in this work (see Theorem 3 for more details).
In
Section 7, we found an equivalence between the solutions of the inhomogeneous Lamé-Navier Equation (
113) in elasticity and the solutions of a inhomogeneous div-curl system (Lemma 1). Later, we apply the results obtained in the previous sections and provide a weak general solution of the inhomogeneous Lamé-Navier Equation (
113). Moreover, we give explicit solutions in appropriate interior or exterior domain in
. We close the section showing that these weak solutions are in fact strong solutions through an embedding argument.
3. An n-Dimensional Generalization of the Div-Curl System
In this section, we are interested in the analysis of an
n-dimensional inhomogeneous Moisil–Teodorescu system, whose component equations give rise to an
n-dimensional generalization of the div-curl system
where
is a paravector-valued function with a vanishing scalar part and
is a
-valued function. More precisely, we will assume that
and
when
and
is paravector-valued when
. By applying the operator
D in both sides of Equation (
32), we can check that the condition
is necessary for Equation (
32) to have a solution. On the other hand, observe that
solves Equation (
32) if and only if
is a solution of
where
and
is the coefficient of
in the expression of
g. Note that the left-hand side of Equation (
33) is the additive inverse of the divergence of
. Meanwhile, the left-hand side of Equation (34) is the generalization of the curl operator to
n-dimensions, which was previously studied in [
19]. It is worth mentioning that the main goal in that article was to find necessary and sufficient conditions to obtain a unique solution
, which depends continuously on
.
In the present work, we will follow a different approach. More precisely, our approach will hinge on solving the embedding of Equations (
33) and (34) into the Clifford structure provided by the equivalent Equation (
32). The steps of our construction method are the following:
Step 1. Find a
-valued solution to Equation (
32) using the Teodorescu transform, that is,
.
Step 2. Proceed then to describe the kernel of the non-paravector component operator
and restrict the right-hand side of Equation (
32) to this class.
Step 3. In turn, our -valued solution becomes a paravector-valued solution after the restriction in Step 2. Afterwards, we use the theory of hyper-conjugate harmonic pairs in order to construct a paravector and monogenic function whose scalar part coincides with the scalar part of w given in Step 1.
Proposition 2. If , then
- (i)
is harmonic if and only if .
- (ii)
is harmonic if and only if
- (ii)
is harmonic if and only if
Proof. Note that, beforehand, . Since in and , we obtain . Finally, taking the scalar part, the paravector part and the non-paravector part of the above-mentioned equation, we reach the conclusions. □
We say that
and
are
hyper-conjugate harmonic pairs in
if
is monogenic in
. We will illustrate now a way to generate monogenic paravector-valued functions when only the scalar part is known. The idea is the same as that used for the three-dimensional case in [
8] [Cor. A.2]. It is worth pointing out here that this is not the only procedure. colorredSee, for instance, the radial integral operator used in reference [
7] [Prop. 2.3] for star-shaped bounded domains in
.
The three-dimensional
singular Cauchy integral operator
satisfies both the identity
and the Plemelj–Sokhotski formulas:
Here,
,
and
. The notation n.t.- means that the limit must be taken non-tangential. The scalar component operator of
acting over scalar valued functions will be denoted by
, and it is of particular interest for the scope of this work. More precisely,
This scalar operator is well known in harmonic analysis, and it is fundamental in the classical Dirichlet problem. Moreover, whenever
be a bounded Lipschitz domain, then
is invertible in
for
[
20], where the value of
depends only of the Lipschitz character of
. Returning to the construction of hyper-conjugate harmonic pairs, a natural way to construct them is through the Cauchy operator Equation (
19), which generates monogenic functions. More precisely,
.
Proposition 3. Let Ω be a bounded Lipschitz domain in with . If is a scalar harmonic function, then , where is monogenic in Ω and . In other words, and are hyper-conjugate harmonic pairs.
Proof. Observe that the Plemelj–Sokhotski formula Equation (
36) describes the trace of the Cauchy operator. As a consequence
By the maximum principle for harmonic functions, we conclude that
in
, as desired. □
As mentioned at the beginning of the present section, the necessary conditions for the equivalent system Equation (
32) to have a solution coincide with the first two hypotheses on
in the following result. Meanwhile, the third boundary hypothesis imposed on
is used to ensure that the solution has vanishing non-paravector part.
Theorem 2. Let Ω be a bounded Lipschitz domain in with , and let . If and in Ω, and on , then a weak solution of the n-dimensional div-curl system in Equations (33) and (34) is given bywhereThis solution is unique up to the gradient of a scalar harmonic function in Ω. Moreover, in the case , we have . Proof. In the proof, we will follow Steps 1, 2 and 3 described above. Using Gauss’ theorem on
(see [
14] [Rmk. A.2.23]), we obtain
By hypothesis,
in
and
on
. It follows that
in
, as desired. On the other hand, Proposition 2 guarantees that
is a scalar harmonic function in
. In turn, Proposition 3 implies that
is monogenic and its scalar part is
, where Equation (
40) is satisfied. Let
, and note that it is purely vectorial by virtue of
. As a consequence, Equation (
39) is reached. Moreover,
This means that
satisfies the equivalent system, Equation (
32). The fact that
belongs to the Sobolev space
is a direct consequence of the properties of the Teodorescu and Cauchy operators. Finally, if
, then
. In turn, this identity implies that
, which means that
vanishes in
. □
Corollary 1. Let Ω and be as in Theorem 2. Then, a right inverse of the n-dimensional generalized curl operator iswith . Moreover, is divergence-free in Ω. Furthermore, if , then . Proof. Taking
in Equation (
39), we readily obtain the expression for the right inverse of the generalized curl operator. To see if
is divergence-free, we will use the alternative expression
Since
is the right inverse of
D in
and
is always monogenic in
, then
which is what we wanted to prove. □
Before closing this section, we must point out that Theorem 2 and Corollary 1 generalize to
n dimensions those results recorded as [
8] [Th. A.1] and [
8] [Cor. A.3], respectively, which were valid for bounded Lipschitz domains in
. On the other hand, as illustrated in the last part of Theorem 2, this construction is mathematically much more interesting in the three-dimensional case, needless to say that this is the physically most relevant. In fact, if
, then the assumptions on
in Theorem 2 and Corollary 1 that involve the non-paravector part disappear, while the hypothesis
becomes the irrevocable condition that
has zero divergence. Finally, we point out that the present work is not the first to make use of Clifford analysis and the construction of hyper-conjugate harmonic pairs to address inhomogeneous Moisil–Teodorescu systems. A recent work in which these tools were employed is [
21].
4. Unbounded Domains
From now on, we will restrict our attention to the case , and analyze the classical div-curl system in unbounded domains . To that end, we require some hypothesis on to guarantee that the operator is invertible in .
In
Section 3, we use the fact that the operator
is invertible in
for all
when
is Lipschitz. We are interested now in the inversion of the operator
in
. To analyze in more detail the range of
p for which this operator is invertible, let us define the
boundary averaging operator A as
choosing
in such way that
. This operator induces a natural mapping
from
to
, where
is the subspace of functions in
with a mean equal to 0. Using the Banach closed range Theorem, it was established in [
8] [Prop. 3.3] that
on
, for all
, when
is Lipschitz. Note also that
does not interfere with the averaging process:
because
This and [
8] [Prop. 3.3]) show that the operator
sends
into itself. Moreover, this operator has an inverse
with
when
is Lipschitz, and
when
is
Lipschitz and
.
The following result provides an alternative form to complete a scalar-valued harmonic function to a paravector-valued monogenic function in the exterior domain .
Proposition 4. Let satisfy the strong local Lipschitz condition and let . If is a scalar harmonic function in such that as , then , where is monogenic in and . In other words, and are hyper-conjugate harmonic pairs in .
Proof. Mimicking the analysis at the end of
Section 2.2, we can assure
, for
. By using the Plemelj–Sokhotski formula in
Equation (
36), we obtain
that is,
. Moreover,
where the last inequality comes from the continuity of the operator
in
. Using the asymptotic behavior of
as
and Equation (
50), it follows that
as
. Due to the the uniqueness of the Dirichlet problem in exterior domains, we readily conclude that
.
in
comes from Equation (
50) and Cauchy’s integral formula is used for exterior domains [
14] [Th. 7.14]. □
4.1. Teodorescu Transform over Unbounded Domains
To start with, note that the Teodorescu transform defined previously in Equation (
17) reduces in the three-dimensional case to
Furthermore, its decomposition is simplified to the expression
In turn, the operators
, defined previously in the
n-dimensional case as Equations (
25)–(27), respectively, are reduced to the following expressions, respectively:
Note that the last identities in the above expressions of
,
and
are derived from Proposition 1. Observe that
are the divergence, gradient and curl of the
Newton potential , respectively. This potential operator (also known as
volume potential) has been extensively studied in various works, such as [
22] [Sec. 2.2] and [
23]. In addition to the role that these component operators have in our construction of solutions for Equation (
57), they also provide a lot of analytical information. For instance,
and
are the Biot–Savart operators for bounded and unbounded domains, respectively.
We recall next some properties of the Teodorescu transform in classical
spaces. As mentioned above, one of the disadvantages of using a kernel without any modification is that its integrability range is reduced. Indeed, if
and
, then
for
(see [
6] [Lem. 2]) in that
Integrating over
yields
which is finite for
. Let
. By utilizing the Fubini–Tonelli theorem, we obtain
From the fact that the kernel and w belong to for , we readily obtain .
4.2. The Div-Curl System over Unbounded Domains
In this stage, we will give an explicit solution for the Equation (
57) on unbounded domains of
satisfying the strong Lipschitz condition with weaker topological constraints. To that end, we will recall some auxiliary results reported in [
7,
8,
24]. Fortunately, the operator theory needed for the quaternionic integral operators over unbounded domains is already well developed [
4,
5,
6]. The novelty now lies in the use of the monogenic completion method discussed in Proposition 4 via the single layer operators.
Let us consider the div-curl system without boundary conditions
where
,
and
in
. Note that the equivalence of the systems in Equations (
32)–(34) is readily verified when
. Moreover, the system in Equation (
57) is equivalent to Equation (
58)
Due to the action of the operator
D to a vector-valued function,
is rewritten in quaternionic notation as
. In the same way as for the bounded case, the mentioned equivalence will be the key in the analysis of the exterior div-curl system.
Theorem 3. Suppose that satisfies the strong local Lipschitz condition and . Let and in . Then, a weak solution of the div-curl system Equation (57) in is given bywhere . This general solution is unique up to the gradient of a harmonic function in . Proof. Proposition 2(i) states that
is harmonic if and only if
. Let us check first that
as
. Let
and define
. The restriction of the scalar part of the Teodorescu transform to the bounded domain
can be estimated by
This means that
as
. Due to
as
, we readily obtain
as
. Now, we will examine the regularity of
. To start with, note that
Moreover, for
.
By Equations (
61) and (
62), we can conclude that
. As a consequence, by Proposition 4,
is monogenic in
and
, where
satisfies
. It follows that
Using this and the decomposition Equation (
52), we obtain
is a purely vector solution of Equation (
58), whence the conclusion of this result follows. □
Note that Equation (
59) can be rewritten as
where
. Define the
single layer potential [
25] as
It is worth pointing out that the Cauchy operator evaluated in scalar functions possesses a decomposition in terms of the operators div and curl [
8]. More precisely,
Using the last equation with
and
as above, and replacing the second and third expressions of Equation (
53) in Equation (
59), we observe that the solution of the div-curl system can be rewritten in a way similar to the classic
Helmholtz decomposition theorem. More precisely, we have the following result.
Corollary 2. Under the same hypothesis of Theorem 3, the solution, Equation (59), admits a Helmholtz-type decomposition of the formwhereand . Moreover, Comparing the decomposition Equation (
69) with the classical Helmholtz decomposition on the entire three-dimensional space [
26] [p. 166] [
27] [Lem. 3.1, 3.2], we readily observe that they adopt similar forms. Note that the vector field
is divergence-free. This follows from
which is a consequence of the first equation of Equation (
53) and the proof of Theorem 3.
5. Div-Curl System in Exterior Domains
In this section, we derive another explicit solution to the div-curl system Equation (
57), this time using another method to generate hyper-conjugate harmonic pairs. The cornerstone now is a radial integral operator defined on an infinite ray instead of the integral equation method provided by the layer potentials.
For the remainder of this manuscript and for the sake of convenience, we will suppose that
is star-shaped w.r.t. the origin. It is worth pointing out that if
is star-shaped w.r.t. any other point, then a simple translation would make it star-shaped w.r.t. the origin. The radial integral operator mentioned above was recently proposed and firstly analyzed in [
10,
11]. There exists an important family of radial integral operators in star-shaped domains, which takes on the form
where usually
. Using standard relations such as
, one may readily verify the following relations:
;
;
;
;
;
;
and
. This family of operators plays an important role in the theory of special functions as well as in mathematical physics. Another interesting application appears in quaternionic analysis when
. Indeed, this radial integral operator generates harmonic functions
for each
which is a harmonic function in
[
7] [Prop. 2.3]. Moreover,
is a quaternion-valued monogenic function in
. This means that the radial operator
provides an explicit way to generate hyper-conjugate harmonic pairs in the star-shaped domain with respect to the origin when the scalar part is known. For convenience, we recall next the main result of [
7].
Theorem 4. (Delgado and Porter [
7] [Th. 4.4])
. Let Ω be a bounded star-shaped domain. If in Ω, and if for , then a general weak solution of the div-curl system is given by Moreover, this solution is unique up to the gradient of a harmonic function in Ω. We now turn our attention to unbounded domains. To start with, note that a similar radial integral operator
acting on functions defined on
was defined in [
10,
11], for star-shaped domains
. Equivalently,
will be
star-shaped w.r.t. infinity, which means that any infinite ray emanating from
is entirely contained in
. In other words,
is
star-shaped w.r.t. infinity if
for all
and
. More precisely, the following operator preserves the above-mentioned properties of the operator
when it is restricted to a class of functions with a suitable behavior at infinity:
Note that this class of scalar-valued harmonic functions is non-empty in that
belongs to
component-wise. The harmonicity in
is straightforward in that it is monogenic in
. Meanwhile, the radiation condition at infinity readily follows. Indeed, observe that
Let us define
by
We will call the exterior monogenic completion operator in light of the next result.
Proposition 5. Let be star-shaped w.r.t. infinity. Then, the operator in Equation (79) sends to . Moreover, for every real-valued harmonic function , In other words, there is a monogenic function w in such that .
Proof. Beforehand, note that
satisfies Equation (
80) if and only if
and
satisfy the div-curl system
Using the fact that
is harmonic in
and some identities from the vector calculus, it follows that the following is satisfied for
:
The action of the rotational operator to the integrand of
is given by
However,
by that hypothesis. We conclude finally that Equation (
84) reduces to
, as desired. □
We must mention that the operator
played a fundamental role in [
10] [Th. 2] to obtain the general solution of the biharmonic equation. It was also crucial to obtain the general solution of the div-curl system in exterior domains when the known data
and
belong to the class of functions
component-wise [
10] [Th. 3]. Our next result is more general in that we consider arbitrary integrable functions in
and not only harmonic functions in the class
.
Theorem 5. Suppose that Ω is a bounded domain and is star-shaped w.r.t. infinity. Let , and in . Then, a weak solution of the div-curl system Equation (57) in is given by This general solution is unique up to the gradient of a harmonic function in .
Proof. The proof is similar to that of Theorem 3; only the generation of the monogenic function whose scalar part coincides with the operator
changes. Consider
By Proposition 2,
is harmonic in
. To check that
is monogenic in
, it only remains to verify that the hypothesis of Proposition 5 holds. In other words, we will show that
belongs to the family of functions in
: let
for some
. By Equation (
61), then
Using the estimation Equation (
87) and letting
, we obtain
as
and
which means that the harmonic function
belongs to the class
. Proposition 5 establishes now that
is monogenic in
. Finally, since this monogenic function has the same scalar part as the quaternionic solution
, then we obtain
satisfies the equivalent system Equation (
58), which is what we wanted to establish. □
In our derivation of the solution to the exterior div-curl problem Equation (
57), we followed a path different from the classical works by Girault and Raviart [
28]. The present solution hinges on the exterior monogenic completion operator
defined in Equation (
79) (which was firstly introduced in [
10,
11]), and on the properties derived in the present work for the component operators of
Corollary 3. Under the hypothesis of Theorem 5, the solution Equation (85) admits a Helmholtz-type decompositionwhere Moreover, is harmonic in .
Proof. By Equation (
53), we obtain
Meanwhile, a simple computation shows that
which yields Equation (
91). On the other hand,
due to
being harmonic in
. It only remains to prove that the second term of Equation (
97) is also harmonic in the exterior domain, but this fact is derived from the fact that
is harmonic in
. □
Unfortunately, unlike the Helmholtz-type decomposition given in Equation (
69), the new decomposition, Equation (
69), is not divergence-free in the exterior domain
. Later, in Theorem 7, the regularity of the solution Equation (
85) will be analyzed as well as its asymptotic behavior.
6. Neumann Boundary-Value Problems
In this stage, we will analyze an exterior Neumann boundary-value problem associated with the div-curl system Equation (
57). More precisely, we will check that there exists a Helmholtz-type solution of the boundary-value problem which preserves the optimal behavior at infinity whenever
belongs to
.
Firstly, note that the normal trace of the solution Equation (
85) is well defined. As a consequence, solving for the following
exterior Neumann boundary-value problem gives:
which is equivalent to solving the Neumann boundary-value problem for the Laplace equation in exterior domains
Here,
and
is the general solution provided by Theorem 5. More precisely,
solves Equation (
98) if and only if
solves Equation (
99). It is the non-uniqueness of the solution of the div-curl system without boundary conditions and the fact that the normal trace of
is well defined that allow us to formulate the equivalent Neumann problem Equation (
99).
Theorem 6. (Neudert and von Wahl [
27] [Th. 2.1])
. Let be a bounded domain with a smooth boundary, let and assume . Then, the Neumann boundary-value problem has a unique solution . Before introducing the main theorem of this section, we will establish some crucial results. To start with, we will prove next that the composition of the exterior monogenic completion operator with preserves the regularity and asymptotic behavior of the Teodorescu transform .
Proposition 6. Let Ω be a bounded domain and be star-shaped w.r.t. infinity. Let and . Then for each , and as .
Proof. That
as
follows from Equation (
87). This implies that
, for each
. Using Equation (
61), we obtain
As a consequence,
as
. By virtue of the fact that
is the radial integral of
in the variable
t over the interval
, we readily obtain
as
. This implies that
, for each
. The
i-th component of
is represented by
. Without loss of generality, we will only calculate the asymptotic behavior of
. Taking the
component of Equation (
102), note that
Thus, as , for all . Since is the radial integral of in the variable t over the interval , we readily obtain as , for all . We conclude that for each . □
Theorem 7. Let Ω be a bounded domain with smooth boundary, let be star-shaped w.r.t. infinity, and suppose that satisfies the strong local Lipschitz condition. Let , , in and . Then, the exterior Neumann boundary-value problem Equation (98) has a unique solutionwhere is the general solution described by Theorem 5 and, in turn, is the unique solution of the Neumann boundary-value problem Equation (99). The asymptotic behavior of the solution is as . Moreover, if , then as Proof. The first part of the proof is reduced to verifying that
in light of Theorem 6 and the equivalence between the systems of Equations (
98) and (
99). Meanwhile, the uniqueness of solutions of Equation (
98) is derived from the uniqueness of solutions of Equation (
99). From Proposition 6,
. In turn, the Sobolev Imbedding Theorem [
17] [Th. 4.12, Part II] assures that
, for
. Consequently,
, for
. The fact that the remaining terms of
in the identity Equation (
100) belong to
follows from known properties of the Teodorescu transform—namely, that
.
Note that as follows from Theorem 6. The fact that satisfies the same decay condition at infinity results from Proposition 6 (when we obtained that as ) and from the asymptotic behavior of the Teodorescu transform which reads as . Therefore, as as we desired.
To establish the last part of the proof, note that
is divergence-free in
in that the Teodorescu transform
is well defined over all
and it is monogenic in
. From Equation (
82), we have
Therefore,
is divergence-free in
, and
By [
27] [Lem. 2.2], it follows that
as
, as needed. □
Corollary 4. Under the same hypothesis of Theorem 7, if and in , and , then the solution Equation (104) satisfies as .□ In [
27], [Th. 3.2] was given an exhaustively classification of the asymptotic behavior of the solutions of the Neumann BVP Equation (
98) under an appropriate functional setting. In that work, the authors used the solutions of the div-curl system in the entire three-dimensional space and correct the boundary values by harmonic vector fields; the second part is similar to the equivalent BVP Equation (
99) considered in this work.
Regularity of the solution: We can apply the Sobolev embedding theorem so that, if we require a higher regularity for the function
, then the range in which the embedding is achieved is improved. If
, then
and
(analogous to the proof of Proposition 6). Consequently,
. Observe that up to this step, we have still not used the extra geometric condition of the domain. If
satisfies the strong local Lipschitz condition, then
for
follows from [
17] [Th. 4.12, Part II]. Therefore, the range of integrability in the hypotheses of Theorem 7 can be improved from
to
.
Corollary 5. Let be as in Theorem 7. If , , in and , then the exterior Neumann boundary-value problem Equation (98) has a unique solutionwhere and . Moreover, the asymptotic behavior of the solution is as . Obviously, we can modify the functional framework of our Neumann boundary-value problem for the div-curl system in the context of weighted Sobolev spaces which, as shown in [
29], gives a correct functional setting to the exterior Neumann problem for the Laplace equation. This approach was also used to analyzed the regularity of the Teodorescu transform in exterior domains [
3].
Right Inverse of the Curl and Double Curl Operator
Let us define the subspace of divergence-free
- functions
as the set
. Taking
in the general solution Equation (
85), we readily obtain a right inverse operator for the curl operator in exterior domains whose complement is a star-shaped domain, namely,
Meanwhile, when
satisfies the strong local Lipschitz condition, the operator
is also a right inverse for the curl operator in
when
. Moreover, both operators are divergence-free invariant in
. Due to the Helmholtz-type decomposition Equation (
91), a right inverse operator of the
operator in exterior domains of star-shaped domains is given by
Similarly, we can obtain a right inverse operator for
by taking
in the Helmholtz-type decomposition Equation (
69). Indeed, if
, then
defined by
is also a right inverse of the
operator in exterior domains satisfying the strong local Lipschitz condition. As a consequence of this discussion, given
, there exists
with the property that
7. Lamé–Navier Equation
In this section, we will apply the results obtained in the previous sections in order to provide an explicit solution to the well known Lamé–Navier problem in elasticity [
30]. Let us consider the inhomogeneous Lamé–Navier equation
where
and
are known as the
first and
second Lamé parameters, respectively. The parameters on the right-hand side of Equation (
113) have physical significance. For example,
denotes the temperature field,
represents the body forces, and the residual strain
defines the vector field
as follows:
It is worth recalling that this system with the right-hand side of Equation (
113) equal to zero was originally introduced by G. Lamé while he was studying the method of separation of variables for solving the wave equation in elliptic coordinates [
31]. Recently, several works have addressed the homogeneous Lamé–Navier equation using the tools of quaternionic analysis. For instance, [
32,
33].
Quaternion algebra was also used in [
34] to give an extension of the classical Kolosov–Muskhelishvili formulas from elasticity to three dimension. This approach is also based on the classical harmonic potential representation due to Papkovich and Neuber, as well as on a monogenic representation. In the latter technique, the main tool is the decomposition of harmonic functions as the sum of a monogenic with an antimonogenic function in the quaternion setting. For the complete details of this decomposition, see [
35]. Here, we will proceed following a completely different path. In fact, we will show that the solutions of Equation (
113) can be constructed by solving a specific div-curl system whose solutions are readily at hand with the theory developed herein.
Lemma 1. Let be a bounded or unbounded domain in . Let satisfy the systemwith being a quaternionic solution of the inhomogeneous Moisil–Teodorescu system in , and Then, is a solution of the inhomogeneous Lamé–Navier Equation (113) in , respectively. Proof. Since
, the identity Equation (
113) can be rewritten as
Let
be a quaternionic solution of the Moisil–Teodorescu system
in
, respectively. Moreover,
, where
is an arbitrary monogenic function in
. Equating the scalar parts of
, we readily obtain
in
, respectively. Now, equating the vector parts of
, we have
The conclusion of the result follows now by Equation (
115). □
If
, then
is monogenic in
, respectively. Moreover, the relation between the homogeneous Lamé–Navier equation, Equation (
113), and the monogenic functions taking values in the quaternions was first observed by G. Moisil [
36]. Using [
7] [Th. 4.4] for star-shaped domains and Theorem 5 for star-shaped domains w.r.t. infinity, we obtain an explicit solution of the following Lamé–Navier equation, Equation (
113).
Theorem 8. Let be a bounded star-shaped domain in , and let satisfy the strong local Lipschitz condition. Let for in the bounded case, and in the unbounded case. Then a general solution of the Lamé–Navier equation, Equation (113), is given bywhere and are quaternionic solutions of the inhomogeneous Moisil–Teodorescu systems in Ω and in , and F and H are arbitrary monogenic functions in , respectively. Moreover, , and these general solutions are unique up to the gradient of a harmonic function in respectively. Proof. By construction,
and
in
, with
and
being the vector parts of the arbitrary monogenic functions in
and
, respectively. It only remains to verify that
f and
h belong to
and
, for
and
, respectively. Since
for
, then
and
for
and
respectively. The result readily follows from Lemma 1 and Theorem 4 for the star-shaped domain
, and from Lemma 1 and Theorem 5 for the exterior domain
. For the unbounded case scenario, the regularity of the solution comes from Corollary 5, due to the fact that
. In the bounded case, the proof of Proposition 6 yields that
defined in Equation (
74) belongs to
, as desired. □
Theorem 9. Let be a bounded Lipschitz domain, and let satisfy the strong local Lipschitz condition. Let for in the bounded case, and in the unbounded case. Then, a general solution of the Lamé–Navier Equation (113) is given bywhere Here, and are quaternionic solutions of the inhomogeneous Moisil–Teodorescu systems in Ω and in , and F and H are arbitrary monogenic functions in , respectively. These general solutions are unique up to the gradient of a harmonic function in , respectively.
Proof. The result readily follows from Lemmas 1 and [
8] [Appendix, Th. A.1] or, alternatively, from the particular case
of Theorem 2 for the interior domain
, Lemma 1 and Theorem 3 for the exterior domain
. □
Note that we can guarantee that the weak solutions in Theorem 8 are continuous solutions in , respectively. Using the Sobolev embedding theorem, for . Alternatively, for , we reach :
Corollary 6. Let as in Theorem 8. Let and , assume that Ω has Lipschitz boundary or satisfies the strong local Lipschitz condition, respectively. Then the weak solutions Equation (118) and () are strong solutions of the Lamé–Navier equation, Equation (113), in , respectively, that is, . Before concluding this work and motivated by one of the reviewers of this work, it is worth noting that the exterior problems considered here are posed in the complement of a star-like region. As a consequence, the cohomological aspects that are handled by Hodge theory are trivial [
37]. On the other hand, traditional Hodge theory is usually in the context of compact manifolds. If one sticks to complements of star-like regions, then operators such as the
operator introduced above seem to be related to the Poincare homotopy operator acting on differential forms.