Nonlocal Trace Spaces and Extension Results for Nonlocal Calculus
Abstract
For a given Lipschitz domain , it is a classical result that the trace space of is , namely any function has a well-defined trace on its codimension-1 boundary and any function on can be extended to a function. Recently, [26] characterizes the trace space for nonlocal Dirichlet problems involving integrodifferential operators with infinite interaction ranges, where the boundary datum is provided on the whole complement of the given domain . In this work, we study function spaces for nonlocal Dirichlet problems with a finite range of nonlocal interactions, which naturally serves a bridging role between the classical local PDE problem and the nonlocal problem with infinite interaction ranges. For these nonlocal Dirichlet problems, the boundary conditions are normally imposed on a region with finite thickness volume which lies outside of the domain. We introduce a function space on the volumetric boundary region that serves as a trace space for these nonlocal problems and study the related extension results. Moreover, we discuss the consistency of the new nonlocal trace space with the classical space as the size of nonlocal interaction tends to zero. In making this connection, we conduct an investigation on the relations between nonlocal interactions on a larger domain and the induced interactions on its subdomain. The various forms of trace, embedding and extension theorems may then be viewed as consequences in different scaling limits.
keywords:
Trace Theorem; Nonlocal Function Space; Inverse Trace Theorem; Trace Space; Finite Nonlocal Interactions; Extension Operator.MSC:
[2010] 46E35 , 47G10 , 35A23 , 35R111 Introduction
In [30], Gagliardo characterizes the trace space of the Sobolev space () for a given bounded Lipschitz domain . The result consists of the following two parts. First, the trace operator from to is linear and continuous, and conversely, one can define a continuous linear extension operator from to . Our goal in this work is to study the trace spaces of some nonlocal function spaces denoted by related to the Dirichlet energies of a class of nonlocal problems. Here, represents the horizon parameter that characterizes the ranges of nonlocal interactions, and with being viewed as a nonlocal “boundary” set of the given domain . The function space is defined as the completion of with respect to the norm
(1.1) |
with the associated semi-norm given by
(1.2) |
Notice that the space also depends on , and we always assume in this paper. The kernel function in (1.2), for and is taken as
(1.3) |
where for any , denotes the characteristic function on the set , and normalizes the moment of . In particular, we have
where denotes the area of the -sphere since
We note that, for different , can be obtained from a rescaling of a given -independent nonnegative kernel defined on by:
(1.4) |
It is easy to see that the nonlocal function space contains all -integrable functions on with finite norms with respect to . Moreover, for any finite and given , the kernel is integrable for and the corresponding space is equivalent to the space; while is non-integrable for and is equivalent to the standard fractional Sobolev space . Moreover, we have the convergence of the space to the local limit as , see e.g., discussions in [13, 50, 44, 28].
Given any domain , let denote the space of all -functions with the norm defined as
(1.5) |
Here the semi-norm is defined as
(1.6) |
where and . Our main result is to show that is the trace space of , i.e., to establish the existence of trace operator and extension operator that define continuous linear maps in between and .
For the ease of presentations, in the following we denote:
(1.7) |
and we note that when , we have and . With these notations, we first show the main results on half spaces as follows.
Theorem 1.1 (Trace theorem on half spaces).
Let and , then there exists a constant independent of and such that for any ,
and therefore
Theorem 1.2 (Inverse trace theorem on half spaces).
Let for some fixed number , , and , then there exists an extension operator such that
where is a constant independent of , and .
Using partition of unity techniques, the above trace theorems in special domains can then be extended to more general domains. which are stated in the theorems below.
Theorem 1.3 (General trace and inverse trace theorems).
Assume that is a bounded and simply connected Lipschitz domain in and is its nonlocal boundary set. There exists a constant depending on the domain , such that for any and ,
On the other hand, for any , and , there exists an extension operator such that
Here , are constants independent of , .
The paper is organized as follows. In Section 2 we discuss the motivation of this work, together with additional definitions and notation relevant to our main results. To provide some insights on the various nonlocal spaces under consideration, we also investigate their scaling properties and consistency with the classical trace spaces in the local limit, as . Moreover, in making these connections, our study also represents an investigation on the relations between nonlocal interactions on a larger domain and the induced interactions on a subdomain of a smaller size or dimension. This further leads to a new way of viewing the various forms of trace, embedding and extension theorems in different function spaces as consequences in different scaling limits, further illustrating the contributions of our study. Sections 3-5 contain the proofs of the aforementioned trace theorems. In particular, to show these trace theorems, while following the footsteps of the proofs for the local trace theorems, we take into account the effect of nonlocal interactions. In Section 3 we provide the proof of Theorem 1.1 with a special case first, which captures the intrinsic effect of nonlocal interactions defined on a larger domain for subdomains. We then extend the results to cases with general using a scaling argument. For the inverse trace theorem, in Section 4 we present the proof of Theorem 1.2 by constructing an extension operator based on the Whitney decomposition. Then, in Section 5 we prove Theorem 1.3 for the general bounded simply connected Lipschitz domain using partition of unity techniques. Lastly, Section 6 summarizes our findings and discusses future research directions.
2 Motivation and Notation
In this section, we first make a few comments on the motivation of our work. We then investigate the consistency and connection between our nonlocal trace space with the trace space in the classical calculus in Section 2.3. We also provide notations and several useful lemmas for the later proofs, including some scaling properties in Section 2.4 and the Whitney decomposition of in Section 2.6.
2.1 Local, nonlocal and fractional modeling
A major motivation of our work comes from nonlocal modeling that are represented by integro-differential equations, in particular, equations involving nonlocal interactions with a finite interaction length. The latter have drawn much attention recently in modelling certain physical systems where the classical models are not most effective. Comparing with the classical local partial differential equation (PDE) models, these equations have the ability to describe these physical phenomena in a setting with reduced regularity requirements allowing singularities and discontinuities to naturally occur [23, 32, 29]. On the other hand, when comparing with the nonlocal integro-differential equations characterized by an infinite lengthscale, compactly supported nonlocal models are computationally more efficient and therefore a more feasible choice for scientific and engineering applications. These extra flexibility and efficiency allow this framework to be used in many different situations involving physical discontinuity such as dynamic fracture [55, 56, 36, 33, 18, 67, 65, 62], corrosion models [16, 17, 37, 40], and heat conduction [7]. The development in this subject has also produced other applications in image processing [41] and population models [14] among many other different fields which can be further seen in [8]. Particularly, nonlocal problems with boundary constraints have become of recent interest in works such as [39, 9, 47, 42, 25, 43, 49, 48, 59, 63, 64, 66]. In nonlocal models, the boundary conditions are normally not imposed on a sharp interface. Rather, they are imposed on a region with non-zero volume which lies outside of the domain, and treating the nonlocal boundary problem improperly can cause artificial phenomena such as a “surface” or “skin” effect [10, 34, 51, 15]. Differs from the local problems, in some nonlocal problems boundary effects play a major role. For example, in nonlocal minimal surface problems, the “stickiness” effect arises and the boundary datum may not be attained continuously [21, 12]. All the above examples indicate that studying the nonlocal boundary conditions and the associated nonlocal trace spaces are critical for the development of nonlocal models.
In this work, we aim to introduce a function space that serves as a trace space for nonlocal problems with constant finite interaction length (the so-called interaction radius or horizon ), and study related extension results. Extension and trace theorems are well-known in the study of classical local problems with boundary constraints. For the case of Sobolev spaces of integer order, these results are well-established long time ago (see, e.g., [1, 57]). For Sobolev spaces with fractional order of differentiability, which can be seen as one type of nonlocal problems with infinite interaction length, the trace space and extension results are studied in [38, 26, 11, 53]. The latter can be useful in studying nonlocal problems with non-homongeneous boundary data, such as those associated with the nonlocal Laplacian and nonlocal -Laplacian, see for example [2, 3, 4, 11, 20, 52]. In [24, 61, 27], trace theorems are developed for nonlocal problems with varying influence horizon , where as approaches the boundary, in a way that the trace spaces of classical Sobolev spaces are recovered. The trace results are also applied to the study of the coupling of nonlocal and local models [60]. To our best knowledge, the definition of trace space and extension results for nonlocal problems with constant finite horizon have not been dealt with so far. These results would extend the knowledge on the trace space in nonlocal calculus and its connection with the trace space in classical calculus. Moreover, the trace theorem and the inverse trace theorem would also provide important mathematical tools for developing well-posed nonlocal models with volumetric boundary conditions, such as discussed in [63].
2.2 Nonlocal Space , Associated Nonlocal Problems and Their Local Limits
Before discussing their connections in the following sections, in this section we introduce the classical and nonlocal Laplacian operators and their corresponding nonlocal function spaces relevant to this paper. The discussions in this subsection are restricted to the Hilbert space setting where .
Given a scalar function , the classical Laplacian operator is defined as and boundary value problems on the domain related to are often associated with the Sobolev space with its norm defined by
On the other hand, when incorporating long-range interactions into the model such that where every point is interacting with a finite neighborhood of points, a nonlocal Laplacian operator is then given by
where is a kernel function that will be prescribed shortly, and
is the interaction domain of . The nonlocal Laplacian operator is associated with the following nonlocal norm
In this paper we further assume that such neighborhood is a Euclidean ball surrounding , i.e., . Here is the interaction radius or horizon. This fact has implications on the boundary conditions that are prescribed on a collar of thickness outside the domain , that we have the interaction domain . In particular, we can take a popular class of kernels as in (1.3). We note that when the constant , we have the following property
(2.1) |
since the moment of kernel in (1.3) is normalized to . Then it is well-known (see, e.g., [23]) that the nonlocal diffusion operator converge to its local for all counterpart pointwise: for any and ,
Moreover, when , its nonlocal norm converges to the norm:
Naturally, we can extend the above conclusion to more general cases of nonlocal and local -Laplacians corresponding to .
2.3 Nonlocal Space and Connection to Classical Local Trace Spaces
First of all, we may view (1.6) as a nonlocal counterpart of the classical trace semi-norm
and seek a nonlocal analog of the classical trace theorem. The relation between the classical trace space and the new nonlocal trace space can be seen from the limiting process as in the following proposition. In the rest of the paper, we use if for a generic constant independent of and . We also write if and .
Proposition 2.4.
Let and be defined as in (1.7) for , then
for any for some . Here denotes the ball centered at with radius in d-1.
Proof.
In this proof, we denote any point by . Similarly is also denoted by . We first have the estimate
To estimate the part, we first note that implies
(2.2) |
for a constant independent of , , and . Notice also that for some . Therefore
For the part, we first want to show for any and , we have
(2.3) |
To show (2.3), we can first assume and without loss of generality. Then by rewriting as and the fact that
we can estimate by two different cases where or . If , then we must have and therefore
On the other hand if , then we have . Therefore, (2.3) is true and this leads to
Lastly, for we note that
and therefore
With the fact that
where the limits are achieved by the dominated convergence theorem, we then obtain
∎
2.4 Change of Variables and Scaling Identities
To further understand the trace theorem and nonlocal spaces of a given and the connections with existing studies in the literature, we consider some scaling identities here. We recall the notation introduced in (1.7) so that for any , we have , which leads to the following scaling argument:
Lemma 2.5.
Given , for any in or , let , then belongs to or and we have respectively
Moreover, the above results also hold for :
Proof.
The proof is obtained by a change of variables. In particular, denoting and , we have
(2.4) |
Similarly, for the trace norm we have
All other identities can be proved similarly. ∎
2.5 Equivalent Semi-norms
We now introduce a lemma that allows us to compare the nonlocal spaces with different sizes of .
Lemma 2.6.
Let and a convex domain. There exists and such that for any ,
Proof.
First, note that
Choose an with . Notice that for . We split into the following
Since is convex by assumption, for any and , we have for each and so is well defined at these points. Now, applying the inequality , we have
Notice that by the change of variables and we have and the Jacobian matrix
(2.5) |
where is the identity matrix. Thus and then
where we have used and . So the left half of the inequality is true with .
Lastly, the right half of the inequality is true with , since
∎
2.6 Dyadic Cubes and Whitney Type Decomposition
The proof of Theorem 1.2 relies on extension results of Whitney type, the subject of which can be found in [58]. Here we focus on defining Whitney type decompositions for the half space and its subdomain . For any , we define the collection of dyadic cubes in , i.e., the cubes of the form where and is the shifted unit cube for . Let denote the side length of the cube , and the collection of cubes with . For , we now define two types of decomposition of the domain using the dyadic cubes for (), which will be useful in Section 4 to prove the inverse trace result. For any , we define .
-
1.
Type I decomposition. Let for some and , then
(2.6) -
2.
Type II decomposition. Let for some , then we define
(2.7)
Naturally, we will write and . Notice that coincides with the classical Whitney decomposition of the half space, where the length of each cube is proportional to the distance between the cube and the boundary of the domain. This type of decomposition is also used to prove the classical and fractional extension results [26, 38]. The Type I decomposition, however, has a special set which touches the boundary and it is used later to construct extension operator for the case .
3 Nonlocal Trace Theorem
In this section we consider the trace theorem on half spaces and provide the proof for Theorem 1.1. We recall that the result stated corresponds to , and . In particular, with the help of the scaling arguments in Lemma 2.5, we first prove the results for the special case , then extend the results to general . Since forms a dense set in , it suffices to prove the conclusion for , which is the case presented in the proofs.
When , we will prove the following theorem:
Theorem 3.7 (Trace theorem on half spaces when ).
Let , then there exist a generic constants depending only on and , such that for any and any ,
(3.1) | ||||
(3.2) |
To prove Theorem 3.7, for any , we define the (hyper)cube associated with by . Now for any , we write , where
for . Then for all , we have . Notice that has non-trivial intersections with the set . Now we use to denote the union of all (hyper)cubes in that have non-trivial intersections with the diagonal line from the center of to the center of . Then we have the following lemma.
Lemma 3.8.
Let , then for any , we have
where is a constant only dependent on and .
Proof.
Let be the number of (hyper)cubes in the set . Then we know that . We denote these (hyper)cubes by , where for each with and . Moreover, for . Therefore, we can connect and by the points for . Then
Now integrate the above equation with respect to over for each , we get
where denotes the ball of radius in . The lemma is then a result of the above estimate and the fact that . ∎
Lemma 3.8 has the following implication.
Corollary 3.9.
For any and any , we have
(3.3) |
where is a constant that depends only on and .
Proof.
Lastly, we can show the estimate in Theorem 3.7. Note that
(3.4) |
where only depends on and . Therefore it suffices to prove the estimate with and then invoke the above inequality for general .
Lemma 3.10 (A nonlocal embedding lemma).
For any and ,
(3.5) |
where is a constant independent of and .
Proof.
In the following we proceed to show the proof of (3.2). In a similar spirit to Corollary 3.9 but with a more refined consideration than the application of a direct Hölder’s inequality, we first state an intermediate result.
Lemma 3.11.
There exists a positive constant depending only on and such that for any positive integer and for any , the following estimate holds:
Proof.
For , let us consider the covering of the path from to given by for . We take for and set . Then, since
integrating the above inequality over with respect to for each , and taking both sides to the power of yields:
Now, we integrate the above inequality over with respect to :
(3.6) |
For the first term above,
where the last step is obtained by Corollary 3.9. For the second term in (3.6) we use the Minkowski’s integral inequality [54] and Lemma 3.8:
In the above derivation, we have used the fact that or for and therefore, the corresponding sets of (hyper)cube, as defined earlier, satisfy . ∎
Remark 1.
We note that the terms in the inequality of the above lemma are well defined by Hölder’s inequality. Indeed, it is easy to see that
for a positive constant depending only on and . Our goal is, however, to show a much refined bound so that the dependence on in the above inequality can be dropped.
Next, we derive a result that helps us to obtain an estimate related to the second term in the above lemma. For , , we present some lemmas to bound the estimates on the nonlocal differences from to and the nonlocal difference from to , respectively. For any , we let be an integer lattice point associated with such that be a (hyper)cube containing . Note that the association may not be unique if is on the boundary of some open (hyper)cube , integer lattice point. In such a case, we may select any of the neighboring to be the associated (hyper)cube. Naturally, if is an integer lattice point itself, we can use the default choice . For , we let be the collection of (hyper)cubes associated with and defined previously. We denote as the number of (hyper)cubes in . We again use for any integer .
Lemma 3.12.
There exists a positive constant depending only on and such that for any and any , the following holds:
(3.7) |
Proof.
we write in the spherical coordinate to get
where
and we have used the fact that .
With the above lemma, we can then bound the nonlocal difference from to with the trace semi-norm. Such an estimate can be seen as the norm of the nonlocal variations along the normal direction (with respect to the strip domain) being controlled by the nonlocal semin-norm for .
Lemma 3.13.
There exists a positive constant depending only on and such that for any , the following holds:
(3.9) |
Remark 2.
Apply a similar argument to the nonlocal difference from to , we can similarly obtain:
(3.10) |
where is a positive constant depending only on and .
We now proceed to investigate the nonlocal difference from to in the following lemma, which can be seen analogously as the norm of the nonlocal tangential variations being controlled by the nonlocal semin-norm:
Lemma 3.14.
There exists a positive constant depending on and such that for any , the following estimate holds:
(3.11) |
Proof.
For notation simplicity, we drop the dependence of on its argument in the derivation here. First notice that from Lemma 3.8 we have
The last inequality in the above estimate is due to the fact that there are at most (hyper)cubes in the set . Therefore,
where we have also used Lemma 2.6 in the last inequality. ∎
We now have the following lemma for the trace semi-norm:
Lemma 3.15.
There exist a positive constant depending only and such that for any ,
(3.12) |
Proof.
We first note that
Moreover, we note that for , it then suffices to show
since for any where only depends on and . Taking the path from to , and then finally , we have
Denoting and integrating the above over with respect to and taking both hand sides to the power of yields:
Considering a fixed , we integrate the above inequality over with respect to , , respectively, then integrate over d-1 with respect to :
Multiplying the above inequalities with and integrating with respect to over d-1 yield:
where the last inequality follows immediate from Lemmas 3.13-3.14. ∎
Proof of Theorem 3.7 and Theorem 1.1. From Lemma 3.10, we have
for any . By taking for in the above inequality, we obtain (3.1). (3.2) is an immediate result in Lemma 3.15.
The proof of the general nonlocal trace Theorem 1.1 on half spaces then follows from Theorem 3.7 and the scaling argument in Lemma 2.5: for any and let , then and
∎
Remark 3.
We have proved in Lemma 3.10 that is bounded by the norm and the semi-norm on a stripe domain , which might appear as a stronger statement that implies the result on a half space in Theorem 3.7. Though, we note that the result on a stripe domain can also be a consequence of Theorem 3.7. Likewise, a bound on the semi-norm can also be obtained on a stripe domain :
To verify the above conclusions, let us derive the results on the stripe domain from (3.1) and (3.2). We take a smooth cutoff function such that supp and for . Denoting such that for and for , then we note that there exists a generic constant independent of and , . Substituting into (3.2), we have
Notice that
otherwise
since . So we have for all . Similarly, substituting into (3.1) for yields
4 Nonlocal Inverse Trace Theorem
For and (), we now define the extension operator . Notice that the kernel is defined in (1.4) with the two cases and . We have the following two cases for the definition of .
Case 1: . We define a partition of unity for according to the decomposition defined in (2.6). For any , let be a smooth function associated with such that is bounded below uniformly on , Lip and supp is contained in an -neighborhood of . Moreover, on . Notice that should also depend on and here we drop the dependence for simplicity of notations. The extension operator is then defined as
(4.1) |
where
and the map is defined for any as
(4.2) |
Case 2: . We similarly define as a partition of unity for according to defined in (2.7). More specifically, for any , is a smooth function bounded below uniformly on , Lip and supp is contained in an -neighborhood of . Moreover, on . Then the extension operator is given by
(4.3) |
where
and the map for any as
(4.4) |
Notice that for such in , if , then we have and .
Remark 4.
The two types of extensions in (4.1) and (4.3) work for and respectively. Notice that if , then the extension in (4.3) gives a continuously function across the boundary to have necessary regularity. On the other hand, the extended function in (4.1) is discontinuous across . Such an extension is fine in this case, since is equivalent to for , and it accepts discontinuous functions. We also note that the map in (4.4) characterizes two regimes – the “fractional regime”, where any cube in is mapped to its symmetric reflection in , and the “classical regime”, where any cube in is mapped to a (hyper)rectangle in . Related discussions on extension operators for the fractional and classical Sobolev spaces using Whitney decompositions can be found in [26] and [38].
Theorem 4.16.
Proof.
We first take the case , where is defined in (4.1). Notice that by the construction of the support of each overlaps with only a finite number of the supports of other function. For the estimate, we have
Thus (4.5) is true.
Now to estimate , we first note that
So we only need to estimate
(4.7) |
and
(4.8) |
Notice that for any ,
where denotes the collections of all the cubes that have nontrivial overlaps with the -neighborhood of . We then have the estimate
We now estimate . Note that
(4.9) |
we then have
where we have used in the last inequality as a result of and . Notice that by (1.4), so
For any ,
Notice that is the constant that appears in (4.6) and it blows up as so we have to take a fixed . Therefore,
where we have used the fact that for any and . So
Now since
as a result of the fact that and , we have
Together we have shown that in (4.7) is bounded by .
To estimate in (4.8), we first note that
(4.10) |
since all cubes in have length greater than or equal to and if . Now suppose for , and then
On the other hand if and , then we know that both and are in , and therefore
In turn, we have
Taking into account the two cases, we can show
where
and
We first estimate . Since the number of sets in is uniformly bounded by a constant for any , and Lip, we have
where we have also used for any and . Therefore
For I.2.a, we first notice that
where we have used the fact that contains only finite number of sets and Lip. Now from (4.9), we obtain
Notice also that
So we have
for some , and denotes the ball of radius centered at . Now chose , then for all . Therefore
Together we have shown (4.6) for the case .
Now for , we take the extension operator defined in (4.3). The estimate of can be shown similarly as in the first case. For the estimate of , similar to the first case considered earlier, it is not hard to see that we only need to estimate
(4.11) |
and
(4.12) |
Similar to the first case, we can split into two parts.
From a similar equation to (4.9), we have
(4.13) |
where we have used and . So
For any and ,
Notice again that the constant blows up as . Then
where we have used the fact that . Using this estimate, one can show
To estimate , we first define a decomposition of
(4.14) |
Then
(4.15) |
Notice that for any and for , we have
Moreover, we must have for the double integral in (4.15) to be non-zero. We thus have the estimate
Therefore
Now for , we can first write
(4.16) |
Observe that
where the second part in the above equation is only nonzero for for (and therefore in this case because of the nonlocal interaction length). Similarly as before, we have
where
and
It is easy to see that can be estimated similarly as and we have . Now for the estimate of , we have two different cases where and . For the case , the estimate follows similarly to the estimate of which is omitted here. For the case , we proceed by noticing that if and , then we must have so that the double integral in (4.16) becomes zero. Therefore we only need to consider (i.e., for ) in this case. Notice that
for and . For , we use (4.13) and
to get
On the other hand, if , then , and we can use
and
to arrive at
Together we have shown and as a result the estimate (4.6) is proved for . ∎
5 Extension to General Lipschitz Domains
In this section we will extend the trace theorems for the half plane to a general Lipschitz domain. We take the strategy to first generalize to a “special” Lipschitz domain before showing the fully general case.
5.1 Some technical lemmas
We will first list some lemmas which are used to show the transformations from the special Lipschitz domain to the half-space are continuous, with detailed proofs elaborated in A.
Lemma 5.17.
Let be a Lipschitz function with Lipschitz constant , , and define the function where for all . Then is a Lipschitz function with the same Lipschitz constant , which is independent of .
Lemma 5.18.
Let and be as defined in Lemma 5.17. Then for any .
Lemma 5.19.
Let and be defined as in the Lemma 5.17. Then for , if
(5.1) | ||||
and , , then
for some positive constants and , independent of .
Lemma 5.20.
5.2 Extension to special Lipshitz domains
We now first present a nonlocal trace theorem to generalize the earlier result shown on a half space with a flat boundary on one side to the case of any infinite domain whose boundary is defined by a Lipshitz graph.
Theorem 5.22.
Consider a special Lipschitz domain, , where there is a Lipschitz function with Lipschitz constant such that
with a collar boundary
Here is as defined in Lemma 5.17. Then the nonlocal trace theorem holds on this domain, i.e.,
for some constant which is independent of , and .
Proof.
Define the operators and
(5.3) | ||||
(5.4) |
Then the operator is defined as
(5.5) |
We will show that is a bounded operator from to with a bounded inverse . Moreover, is a bounded operator from to .
To show is a bounded operator from to , we let , and the constant as defined in Lemma 5.19. From definition, we have , and . Then,
(5.6) |
and
where we have used Lemma 5.18 for the last estimate as well as (5.6). Together these estimates show
where the constant is independent of and .
To show that exists and it is a bounded operator from to , we note first that one can get (5.2) from (5.1). It is then easy to check that can be defined by for and where is given by (5.2). Using a change of variable estimate given by Lemma 5.19 we can then deduce the continuity of the same way as we have done for with details omitted.
Now to show that is a bounded operator from to , we note that , and by applying Lemma 2.6 with and where is defined in Lemma 5.20, we have
Here the constant is independent of and .
Applying Lemma 5.20 part (a) to and then using a change of variable, we have
Similarly applying Lemma 5.20 to and along with the proper change of variables, we finally have where is independent of and . Thus, we have the continuity of . Using proven properties of and along with the trace theorem for the half-plane, we have
∎
Using the transformation operators and , we can also generalize the inverse nonlocal trace theorem to the special Lipschitz domain.
Theorem 5.23.
Proof.
Notice that the inverse operator can be defined as
Now we wish to show that is a bounded operator from to . Using Lemma 5.21 and the defined there along with the appropriate changes of the variables, we can show similarly as we have done in Theorem 5.22 that
where we have used Lemma 2.6 with and in the last step. The continuity of is thus shown.
Finally, using the continuity properties of , and we have
∎
5.3 Extension to more general Lipshitz domains
The extension to more general Lipshitz domains can be otained by using the partition of unity technique. We first decompose the boundary collar region into finitely many balls so that we can locally view the boundary as multiple special Lipschitz domains. From there we apply Theorem 5.22 to each part and join the estimates together with a partition of unity. The detailed derivation is given as follows.
First, for all discussions given in this subsection, we state some assumptions on the domains and define the necessary spaces and functions.
We consider a general bounded simply connected Lipschitz domain , which naturally makes a compact set for any finite . Since has a Lipschitz boundary (see, e.g., [31] Def 1.2.1.1), there exist local coordinate systems for , a collection of balls , and Lipschitz functions for some , such that ,
Notice that in the above definition, are fixed points on . Letting , we have and also a positive number such that . Defining from as in Lemma 5.17 we have
where and . Notice here the Lipschitz constants of and depend on the domain , and the collection of balls only and are thus independent of and . Then given we have along with functions so that and where
Define functions such that
-
1.
,
-
2.
and on .
Since is covered by , we have for a fixed constant depending only on the maximum number of overlapped balls in the set . We also define
for each . Then since is uniformly bounded above and below, we also have .
Furthermore, we would like to define the extension operator on the general domain. By Theorem 5.23, there exists an extension operator such that for and
For any , we assume . Then we can define the extension operator by
(5.7) |
for any . Before continuing we show some useful estimates which will together swiftly prove Theorem 1.3.
Lemma 5.24.
Let be a simply connected Lipschitz domain with an interaction domain for where , , and , as well as , are defined as above. For any , , where is a constant independent of and .
Proof.
First note for ,
and hence
Notice that we can extend by zero so that can be viewed as a function on . Using the fact that we have
(5.8) |
Notice that for each , and we have
(5.9) | |||
where depends on , and . Therefore, the integral term in (5.8) is bounded by a multiple of . Moreover, using the compact support of ,
(5.10) |
Therefore (5.10) along with the estimates of (5.8) and (5.9)
where the constant is independent of and . ∎
Lemma 5.25.
Let be a simply connected Lipschitz domain with an interaction domain for where , , and are defined as above, along with and for each . For any , we have
where is independent of and .
Proof.
First, note that and also if and , then
Moreover, using the compact support of ,
where the last equality is because since . Then for each ,
(5.11) |
where is independent of and . Also, since ,
The estimates for can be done similarly. ∎
Lemma 5.26.
Let be a simply connected Lipschitz domain with an interaction domain for where , , and , along with , are defined as above. For any , we have
where is independent of and .
Proof.
Since is supported in and , it is obvious that . Now
(5.12) |
For the first term in the last line, we have
Now for any ,
where only depends on . Note that in the last step above, we have adopted the -th local coordinate system defined before to represent the points and without adding labels to them. Define such that for any , with the same -th local coordinate system representation
(5.13) |
We can see from Lemma 5.19 that , so . By the change of variable,
Let and , the last line can be estimated by
where is independent of and . Collecting the above estimates, we have
Now for the second term in (5.12), we notice that if and , then , therefore for each
Note that in the last line, we have again adopted the -th local coordinate system. Let be defined in (5.13) and . Then by the same reasoning as above,
Since , the above quantity can be bounded by
where is a constant that depends on , , and . Combining the estimates, we get
∎
Lemma 5.27.
Let be a simply connected Lipschitz domain with an interaction domain for where , , and are defined as above. Assume each as above and defined by (5.7), then for any , where is independent of and .
Proof.
6 Conclusion and Discussion
This work gives suitable characterizations of the trace spaces of a class of nonlocal function spaces denoted by , where the parameter is the nonlocal interaction length and characterizes the singularity of the nonlocal interaction kernels. Such nonlocal function spaces have been extensively used recently as the energy spaces associated with nonlocal diffusion and nonlocal mechanics models [22, 55, 46, 19, you2020data, FR2]. However, a clear understanding of the trace spaces of has been largely limited [61]. In the current work, we have introduced the function space as the trace space of and demonstrated that the trace map from to is continuous (given by the trace theorem) and conversely, there is a continuous extension operator from to (given by the inverse trace theorem). Moreover, the estimates on the trace and the inverse trace maps are uniform with respect to the horizon parameter , so that one can recover the classical trace and inverse trace theorems in the local limit as the nonlocal interaction length . This is also important since there are many instances of nonlocal models that recover the classical diffusion or elasticity equations as [22, 45, 46].
The investigation of trace spaces of Sobolev spaces has been a classical research area that has important implications in the mathematical and numerical studies of boundary value problems of local PDEs. The results of this work therefore are expected to be helpful in the rigorous studies of nonlocal equations with possible nonlocal boundary constraints similar to their PDEs counterparts. Studies in this direction are currently underway. Moreover, nonlocal functions spaces on vector fields such as those appear in [45, 46] can also be studied in the future. Another interesting direction for the future is to investigate, when , the consistency of suitably defined nonlocal spaces, similar to those discussed in this work, with their fractional limits.
Acknowledgments
Y. Yu is supported by the National Science Foundation under award DMS 1753031. Q. Du is supported in part by the National Science Foundation under award DMS-2012562 and the ARO MURI Grant W911NF-15-1-0562. X. Tian is supported by the National Science Foundation under award DMS-2111608.
References
- Adams and Fournier [2003] Adams, R. A. and Fournier, J. J. (2003). Sobolev spaces. Elsevier.
- Andreu et al. [2008] Andreu, F., Mazón, J. M., Rossi, J. D., and Toledo, J. (2008). A nonlocal p-laplacian evolution equation with neumann boundary conditions. Journal de mathématiques pures et appliquées, 90(2):201–227.
- Andreu et al. [2009] Andreu, F., Mazón, J. M., Rossi, J. D., and Toledo, J. (2009). A nonlocal p-laplacian evolution equation with nonhomogeneous dirichlet boundary conditions. SIAM journal on mathematical analysis, 40(5):1815–1851.
- Andreu-Vaillo et al. [2010] Andreu-Vaillo, F., Mazón, J. M., Rossi, J. D., and Toledo-Melero, J. J. (2010). Nonlocal diffusion problems. Number 165. American Mathematical Soc.
- Avkhadiev [2006] Avkhadiev, F. G. (2006). Hardy-type inequalities on planar and spatial open sets. Proceedings of the Steklov Institute of Mathematics, 255(1):2–12.
- Avkhadiev and Nasibullin [2014] Avkhadiev, F. G. and Nasibullin, R. (2014). Hardy-type inequalities in arbitrary domains with finite inner radius. Siberian Mathematical Journal, 55(2):191–200.
- Bobaru and Duangpanya [2010] Bobaru, F. and Duangpanya, M. (2010). The peridynamic formulation for transient heat conduction. International Journal of Heat and Mass Transfer, 53:4047–4059.
- Bobaru et al. [2015] Bobaru, F., Foster, J., Geubelle, P., and Silling, S. (2015). Handbook of Peridynamic Modeling. Modern Mechanics and Mathematics. Taylor & Francis.
- Bobaru et al. [2016] Bobaru, F., Foster, J. T., Geubelle, P. H., and Silling, S. A. (2016). Handbook of peridynamic modeling. CRC press.
- Bobaru and Ha [2011] Bobaru, F. and Ha, Y. D. (2011). Adaptive refinement and multiscale modeling in 2D peridynamics. International Journal for Multiscale Computational Engineering, 9(6).
- Bogdan et al. [2020] Bogdan, K., Grzywny, T., Pietruska-Pałuba, K., and Rutkowski, A. (2020). Extension and trace for nonlocal operators. Journal de Mathématiques Pures et Appliquées, 137:33–69.
- Borthagaray et al. [2019] Borthagaray, J. P., Li, W., and Nochetto, R. H. (2019). Finite element discretizations of nonlocal minimal graphs: convergence. arXiv preprint arXiv:1905.06395.
- Bourgain et al. [2001] Bourgain, J., Brezis, H., and Mironescu, P. (2001). Another look at Sobolev spaces, pages 439–455. IOS Press.
- Carrillo and Fife [2005] Carrillo, C. and Fife, P. (2005). Spatial effects in discrete generation population models. Journal of Mathematical Biology, 50(2):161–188.
- Chen et al. [2020] Chen, J., Jiao, Y., Jiang, W., and Zhang, Y. (2020). Peridynamics boundary condition treatments via the pseudo-layer enrichment method and variable horizon approach. Mathematics and Mechanics of Solids, page 1081286520961144.
- Chen and Bobaru [2015] Chen, Z. and Bobaru, F. (2015). Peridynamic modeling of pitting corrosion damage. Journal of the Mechanics and Physics of Solids, 78:352 – 381.
- Chen et al. [2015] Chen, Z., Zhang, G., and Bobaru, F. (2015). The influence of passive film damage on pitting corrosion. Journal of the Electrochemical Society, 163:C19–C24.
- Cheng et al. [2015] Cheng, Z., Zhang, G., Wang, Y., and Bobaru, F. (2015). A peridynamic model for dynamic fracture in functionally graded materials. Composite Structures, 133:529–546.
- D’Elia et al. [2020] D’Elia, M., Du, Q., Glusa, C., Gunzburger, M., Tian, X., and Zhou, Z. (2020). Numerical methods for nonlocal and fractional models. Acta Numerica, 29:1–124.
- Dipierro et al. [2017a] Dipierro, S., Ros-Oton, X., and Valdinoci, E. (2017a). Nonlocal problems with neumann boundary conditions. Revista Matemática Iberoamericana, 33(2):377–416.
- Dipierro et al. [2017b] Dipierro, S., Savin, O., and Valdinoci, E. (2017b). Boundary behavior of nonlocal minimal surfaces. Journal of Functional Analysis, 272(5):1791–1851.
- Du et al. [2012] Du, Q., Gunzburger, M., Lehoucq, R. B., and Zhou, K. (2012). Analysis and approximation of nonlocal diffusion problems with volume constraints. SIAM review, 54(4):667–696.
- Du et al. [2013] Du, Q., Gunzburger, M., Lehoucq, R. B., and Zhou, K. (2013). A nonlocal vector calculus, nonlocal volume-constrained problems, and nonlocal balance laws. Mathematical Models and Methods in Applied Sciences, 23(03):493–540.
- Du et al. [2021] Du, Q., Mengesha, T., and Tian, X. (2021). Fractional hardy-type and trace theorems for nonlocal function spaces with heterogeneous localization. arXiv preprint arXiv:2103.09387.
- Du et al. [2017] Du, Q., Tao, Y., and Tian, X. (2017). A peridynamic model of fracture mechanics with bond-breaking. Journal of Elasticity, pages 1–22.
- Dyda and Kassmann [2019] Dyda, B. and Kassmann, M. (2019). Function spaces and extension results for nonlocal dirichlet problems. Journal of Functional Analysis, 277(11):108134.
- Foss [2020] Foss, M. (2020). Traces on general sets in for functions with no differentiability requirements. arXiv preprint arXiv:2007.00863.
- Foss and Radu [2016] Foss, M. and Radu, P. (2016). Differentiability and integrability properties for solutions to nonlocal equations. In New Trends in Differential Equations, Control Theory and Optimization: Proceedings of the 8th Congress of Romanian Mathematicians, pages 105–119. World Scientific.
- Foss et al. [2018] Foss, M. D., Radu, P., Wright, C., et al. (2018). Existence and regularity of minimizers for nonlocal energy functionals. Differential and Integral Equations, 31(11/12):807–832.
- Gagliardo [1957] Gagliardo, E. (1957). Caratterizzazioni delle tracce sulla frontiera relative ad alcune classi di funzioni in variabili. Rendiconti del seminario matematico della universita di Padova, 27:284–305.
- Grisvard [1985] Grisvard, P. (1985). Elliptic Problems in Nonsmooth Domains. Monographs and studies in mathematics. Pitman Advanced Pub. Program.
- Gunzburger and Lehoucq [2010] Gunzburger, M. and Lehoucq, R. B. (2010). A nonlocal vector calculus with application to nonlocal boundary value problems. Multiscale Modeling & Simulation, 8(5):1581–1598.
- Ha and Bobaru [2010] Ha, Y. D. and Bobaru, F. (2010). Studies of dynamic crack propagation and crack branching with peridynamics. International Journal of Fracture, 162(1):229–244.
- Ha and Bobaru [2011] Ha, Y. D. and Bobaru, F. (2011). Characteristics of dynamic brittle fracture captured with peridynamics. Engineering Fracture Mechanics, 78(6):1156–1168.
- Hardy et al. [1952] Hardy, G., Littlewood, J., and Polya, G. (1952). Inequalities cambridge univ. Press, Cambridge, (1988).
- Hu et al. [2012] Hu, W., Ha, Y. D., and Bobaru, F. (2012). Peridynamic model for dynamic fracture in unidirectional fiber-reinforced composites. Computer Methods in Applied Mechanics and Engineering, 217:247–261.
- Jafarzadeh et al. [2017] Jafarzadeh, S., Chen, Z., and Bobaru, F. (2017). Peridynamic modeling of repassivation in pitting corrosion of stainless steel. CORROSION, 74.
- Koskela et al. [2017] Koskela, P., Soto, T., and Wang, Z. (2017). Traces of weighted function spaces: dyadic norms and whitney extensions. Science China Mathematics, 60(11):1981–2010.
- Le and Bobaru [2018] Le, Q. and Bobaru, F. (2018). Surface corrections for peridynamic models in elasticity and fracture. Computational Mechanics, 61(4):499–518.
- Li et al. [2017] Li, S., Chen, Z., Tan, L., and Bobaru, F. (2017). Corrosion-induced embrittlement in zk60a mg alloy. Materials Science and Engineering: A, 713.
- Lou et al. [2010] Lou, Y., Zhang, X., Osher, S., and Bertozzi, A. (2010). Image recovery via nonlocal operators. Journal of Scientific Computing, 42(2):185–197.
- Macek and Silling [2007] Macek, R. W. and Silling, S. A. (2007). Peridynamics via finite element analysis. Finite Elements in Analysis and Design, 43(15):1169–1178.
- Madenci and Oterkus [2014] Madenci, E. and Oterkus, E. (2014). Peridynamic theory. In Peridynamic Theory and Its Applications, pages 19–43. Springer.
- Mengesha [2012] Mengesha, T. (2012). Nonlocal korn-type characterization of sobolev vector fields. Communications in Contemporary Mathematics, 14(04):1250028.
- Mengesha and Du [2014a] Mengesha, T. and Du, Q. (2014a). The bond-based peridynamic system with dirichlet-type volume constraint. Proc. Roy. Soc. Edinburgh Sect. A, 144(1):161–186.
- Mengesha and Du [2014b] Mengesha, T. and Du, Q. (2014b). Nonlocal constrained value problems for a linear peridynamic navier equation. Journal of Elasticity, 116(1):27–51.
- Oterkus [2010] Oterkus, E. (2010). Peridynamic theory for modeling three-dimensional damage growth in metallic and composite structures. PhD thesis, The University of Arizona.
- Oterkus [2015] Oterkus, S. (2015). Peridynamics for the solution of multiphysics problems.
- Oterkus et al. [2014] Oterkus, S., Madenci, E., and Agwai, A. (2014). Peridynamic thermal diffusion. Journal of Computational Physics, 265:71–96.
- Ponce [2004] Ponce, A. C. (2004). An estimate in the spirit of Poincaré’s inequality. J. Eur. Math. Soc, 6(1):1–15.
- Prudhomme and Diehl [2020] Prudhomme, S. and Diehl, P. (2020). On the treatment of boundary conditions for bond-based peridynamic models. Computer Methods in Applied Mechanics and Engineering, 372:113391.
- Ros-Oton [2016] Ros-Oton, X. (2016). Nonlocal elliptic equations in bounded domains: a survey. Publicacions matematiques, pages 3–26.
- Rutkowski [2020] Rutkowski, A. (2020). Function spaces and the Dirichlet problem for nonlocal operators. Ph.D. thesis, Wrocław University of Science and Technology, Wrocław, Poland.
- Schep [1995] Schep, A. R. (1995). Minkowski’s integral inequality for function norms. In Operator theory in function spaces and Banach lattices, pages 299–308. Springer.
- Silling [2000] Silling, S. A. (2000). Reformulation of elasticity theory for discontinuities and long-range forces. Journal of the Mechanics and Physics of Solids, 48(1):175–209.
- Silling et al. [2007] Silling, S. A., Epton, M., Weckner, O., Xu, J., and Askari, E. (2007). Peridynamic states and constitutive modeling. Journal of Elasticity, 88(2):151–184.
- Slobodetskiĭ [1958] Slobodetskiĭ, L. N. (1958). S.l. sobolev’s spaces of fractional order and their application to boundary problems for partial differential equations. In Doklady Akademii Nauk, volume 118, pages 243–246. Russian Academy of Sciences.
- Stein [1970] Stein, E. M. (1970). Singular Integrals and Differentiability Properties of Functions (PMS-30), Volume 30. Princeton University Press.
- Tao et al. [2017] Tao, Y., Tian, X., and Du, Q. (2017). Nonlocal diffusion and peridynamic models with neumann type constraints and their numerical approximations. Applied Mathematics and Computation, 305:282–298.
- Tao et al. [2019] Tao, Y., Tian, X., and Du, Q. (2019). Nonlocal models with heterogeneous localization and their application to seamless local-nonlocal coupling. Multiscale Modeling & Simulation, 17(3):1052–1075.
- Tian and Du [2017] Tian, X. and Du, Q. (2017). Trace theorems for some nonlocal function spaces with heterogeneous localization. SIAM Journal on Mathematical Analysis, 49(2):1621–1644.
- Trask et al. [2019] Trask, N., You, H., Yu, Y., and Parks, M. L. (2019). An asymptotically compatible meshfree quadrature rule for nonlocal problems with applications to peridynamics. Computer Methods in Applied Mechanics and Engineering, 343:151–165.
- You et al. [2020a] You, H., Lu, X., Trask, N., and Yu, Y. (2020a). An asymptotically compatible approach for neumann-type boundary condition on nonlocal problems. ESAIM: Mathematical Modelling and Numerical Analysis, 54(4):1373–1413.
- You et al. [2020b] You, H., Yu, Y., and Kamensky, D. (2020b). An asymptotically compatible formulation for local-to-nonlocal coupling problems without overlapping regions. Computer Methods in Applied Mechanics and Engineering, 366:113038.
- Yu et al. [2018] Yu, Y., Bargos, F. F., You, H., Parks, M. L., Bittencourt, M. L., and Karniadakis, G. E. (2018). A partitioned coupling framework for peridynamics and classical theory: analysis and simulations. Computer Methods in Applied Mechanics and Engineering, 340:905–931.
- Yu et al. [2021] Yu, Y., You, H., and Trask, N. (2021). An asymptotically compatible treatment of traction loading in linearly elastic peridynamic fracture. Computer Methods in Applied Mechanics and Engineering, 377:113691.
- Zhang et al. [2018] Zhang, G., Gazonas, G. A., and Bobaru, F. (2018). Supershear damage propagation and sub-rayleigh crack growth from edge-on impact: A peridynamic analysis. International Journal of Impact Engineering, 113:73 – 87.
Appendix A Proof of Lemmas in Section 5
Proof of Lemma 5.17. First, note that the function is well-defined for each since
there must exist at least one such that by the continuity of . The minimum of such is the value of and clearly by how the were chosen.
Let and without loss of generality assume . By definition of , there is a such that . Now define
Notice that we must have . If this was not the case, then consider function
for . By the assumption, we have
Since is continuous, there must be a such that for some . Then,
which is a contradiction of the definition of . With this inequality we finally have
∎
Proof of Lemma 5.18. Let and define and . By the definition of , there is a such that where . Note that where . If , then we are done since then
Otherwise, consider the triangle made by the points , , and as shown in Figure 1. ***Here we are viewing the triangle by taking the natural isometry to from an orthonormal basis for span. Note that and are linearly independent under the assumption . If this was not the case, then requiring . Since is a single-valued function, , hence , giving a contradiction. Let be the angle associated with the vertex and note that is a Lipschitz function with Lipschitz constant . Then,
So if is the angle associated with the vertex , then by the law of sines
and so . By the definition of we clearly have . Since was arbitrary we finish the proof.
∎
Proof of Lemma 5.19. By Lemma 5.17 is a Lipschitz function with Lipschitz constant . Then in combination with Lemma 5.18,
For left inequality, first notice from (5.1) that
Additionally, implies so then,
∎
Proof of Lemma 5.20. For (a) if , then
Therefore we have
For (b), with and noting we have
By the same reason, we have
We similarly have (c) by noticing that . ∎