We study equations from the area of peridynamics, which is a nonlocal extension of elasticity. The governing equations form a system of nonlocal wave equations. We take a novel approach by applying operator theory methods in a systematic way. On the unbounded domain ℝn, we present three main results. As main result 1, we find that the governing operator is a bounded function of the governing operator of classical elasticity. As main result 2, a consequence of main result 1, we prove that the peridynamic solutions strongly converge to the classical solutions by utilizing, for the first time, strong resolvent convergence. In addition, main result 1 allows us to incorporate local boundary conditions, in particular, into peridynamics. This avenue of research is developed in companion papers, providing a remedy for boundary effects. As main result 3, employing spherical Bessel functions, we give a new practical series representation of the solution which allows straightforward numerical treatment with symbolic computation.
I. INTRODUCTION AND MOTIVATION
Classical elasticity has been successful in characterizing and measuring the resistance of materials to crack growth. On the other hand, peridynamics (PD), a nonlocal extension of continuum mechanics developed by Silling,55 is capable of quantitatively predicting the dynamics of propagating cracks, including bifurcation. Its effectiveness has been established in sophisticated applications such as Kalthoff-Winkler experiments of the fracture of a steel plate with notches,31,58 fracture and failure of composites, nanofiber networks, and polycrystal fracture.34,44,60,59 Further applications are in the context of multiscale modeling, where PD has been shown to be an upscaling of molecular dynamics52,54 and has been demonstrated as a viable multiscale material model for length scales ranging from molecular dynamics to classical elasticity.10 Also see other related engineering applications,13,33,35,46,45 the review, and news articles14,16,21,57 for a comprehensive discussion, and the book.38
We study a class of nonlocal wave equations. The driving application is PD whose equation of motion corresponds exactly to the nonlocal wave equation under consideration. The same operator is also employed in nonlocal diffusion.9,14,51 Similar classes of operators are used in numerous applications such as population models, particle systems, phase transition, and coagulation. Nonlocal operators have also been used in image processing.27,36 In addition, we witness a strong interest for PD, its applications, and related nonlocal problems addressing, for instance, conditioning analysis, domain decomposition, and variational theory,5–7 discretization,1,7,25,62 numerical methods,15,17,50 nonlinear PD,26,37 nonlinearity in nonlocal wave equations,19,20 well-posedness in various forms,5,6,8,18,22–24,39,40,67 and other aspects.12,28,41,53
It is part of the folklore in physics that the point particle model, which is the root for locality in physics, is the cause of unphysical singular behavior in the description of the underlying phenomena. This fact is a strong indication that, in the long run, the development of nonlocal theories is necessary for description of natural phenomena. Similar set of operator theory tools used for local problems can be transferred to study nonlocal problems because operator theory does not discern the locality or nonlocality of the governing operator. This article adds valuable tools to the arsenal of methods to analyze nonlocal problems.
The equation of interest falls into the class of abstract evolution equations, more precisely, abstract linear wave equations. Methods from operator theory are ideal to treat such equations. One can directly gain access to powerful tools such as functional calculus for bounded self-adjoint operators, spectral theorems (of densely defined, linear, and self-adjoint operators in Hilbert spaces), and strong resolvent convergence of operators. Well-posedness of the initial value problem, conservation of energy, stability of solutions (merely determined by the spectrum) are all immediately available; see Theorem 1 and Corollary 2. Furthermore, representation of solutions can easily be constructed through functional calculus for bounded operators; see our main result 3.
There are many studies6–8,18,24,39,40,56,67 that aim to establish a connection between the nonlocal operator and the classical one. We denote this nonlocal-to-local connection by NL2L. PD formulation utilizes a distance parameter δ, called the horizon, in order to limit interactions to nearby points. The source of nonlocality is the horizon and it is embedded in the support of the micromodulus function. Studies turn to the limit of vanishing nonlocality, i.e., δ → 0 when they aim to establish NL2L. The exact expression of NL2L is lost due to taking limit. We take general micromodulus functions in L1(ℝn), not necessarily parametrized by δ. Without resorting to a δ → 0 argument, we identify the exact expression of the NL2L connection, what we call as main result 1.
Identifying the NL2L connection leads to a fruitful research direction because it helped us to extend the construction for unbounded domains to the bounded domain case by using the same function of the classical operator used in PD. That way, we kept a close proximity to PD in the bounded domain case. More importantly, NL2L connection turned out to be the notable result that the governing operator is a bounded function of the classical (local) operator.69 This has far reaching consequences. It enables the incorporation of local boundary conditions into nonlocal theories; see the companion papers.2,3 It takes a lot of effort to mitigate the boundary effects. See the comprehensive discussion [Ref. 38, Chaps. 4, 5, 7, and 12] and the numerical study.32 Incorporation of local boundary conditions provides a remedy for boundary effects seen in PD.
The rest of the article is structured as follows. We provide mathematical introduction in Section I. In Section II, we set the operator theory framework to treat the nonlocal wave equation. We prove basic properties of the solutions such as well-posedness of the initial value problem and conservation of energy. We provide a representation of the solutions in terms of bounded functions of the governing operator. We study the stability of solutions. We also consider the class of inhomogeneous wave equations and prove well-posedness of the corresponding initial value problem as well as a representation of the solutions in terms of bounded functions of the governing operator.
In Section III, in the vector-valued case, we note that the governing operator becomes an operator matrix. The generality of operator theory allows a simple extension of the results established for the scalar-valued functions to the vector-valued ones. We prove the boundedness of the entries of the governing operator matrix. The proof is natural due to operator theory again, because it relies on a well-known criterion for integral operators. We present a “diagonalization” of the matrix entries. This is accomplished by employing the unitary Fourier transform and connecting the entries to maximal multiplication operators.
In Section IV, we collect the statements of the three main results:
Main Result 1: Nonlocal operator is a bounded function of the classical operator.
Main Result 2: Strong convergence of nonlocal solutions to classical ones through strong resolvent convergence.
Main Result 3: Representation of the solution in terms of spherical Bessel functions.
In Section V, we construct two examples to demonstrate the usefulness of main result 3. We give explicit representation of the solutions in terms of Bessel functions. Since the explicit representation is available, we easily compute the solutions by symbolic computation, depict the resulting solutions of PD wave equation, and compare to the classical solutions.
In Section VI, we collect the proofs of the three main results together with related supporting material. In particular, in Section VI A, utilizing Fourier transforms, we turn the nonlocal governing operator into maximal multiplication operator. This process can be viewed as a form of “diagonalization.” The spectra of maximal multiplication operators are well understood. In addition, the functional calculus associated to maximal multiplication operators is known. Spectra and functional calculus allow the construction of the functional calculi for the governing operator, which in turn is used to prove that the governing operators with spherically symmetric micromoduli are functions of the Laplace operator. In Section VI B, we prove strong convergence of nonlocal solutions to classical ones through the notion of strong resolvent convergence of functions of the classical operator. Different types of resolvent convergences, i.e., norm, strong, and weak, are available; see Ref. 65 for comparison. We utilize the notion of strong resolvent convergence in order to obtain strong convergence of solution operators. This allows us to prove main result 3. We give examples of sequences of micromoduli that are instances of this result. In Section VI C, we consider the calculation of the solution of the wave equation. Since the governing operator is bounded, holomorphic functions of that operator can be represented in the form of power series in the operator. Then, we give a representation of holomorphic functions, present in the solution of the initial value problem, utilizing the fact that the governing operator is a sum of two commuting operators. We establish main result 3 by discovering that the corresponding power series can be given in terms of a series of Bessel functions. In addition, we provide an error estimate for the series representation. We use this estimate in the plots of solutions of the PD wave equation; see Figures 1 and 2. In Section VI D, we give the explicit representation of the solutions of inhomogeneous wave equations and prove well-posedness of the corresponding initial value problem. We conclude in Section VII.
Evolution of the local and nonlocal wave equation solutions with vanishing initial velocity ((a) and (b)) and vanishing initial displacement ((c) and (d)). For (a) and (c), we use ρ = E = 1, b = 0, values in (1.2). For (b) and (d), we use c = a = 1, ρ = 1, σ = 1, σd = 1/2 values in Example 10. Solutions are generated using symbolic computation. The infinite series in (5.1) is truncated after adding 46 terms. No difference has been observed visually if more terms are added.
Evolution of the local and nonlocal wave equation solutions with vanishing initial velocity ((a) and (b)) and vanishing initial displacement ((c) and (d)). For (a) and (c), we use ρ = E = 1, b = 0, values in (1.2). For (b) and (d), we use c = a = 1, ρ = 1, σ = 1, σd = 1/2 values in Example 10. Solutions are generated using symbolic computation. The infinite series in (5.1) is truncated after adding 46 terms. No difference has been observed visually if more terms are added.
Evolution of the local and nonlocal wave equation solutions with discontinuous initial displacement ((a) and (b)) and discontinuous initial velocity ((c) and (d)). For (a) and (c), we use ρ = E = 1, b = 0, values in (1.2). For (b) and (d), we use c = a = 1, ρ = 1, σ = 1, and b = ϵ = 1 values in Example 11. Solutions are generated using symbolic computation. The infinite series in (5.2) is truncated after adding 46 terms. No difference has been observed visually if more terms are added.
Evolution of the local and nonlocal wave equation solutions with discontinuous initial displacement ((a) and (b)) and discontinuous initial velocity ((c) and (d)). For (a) and (c), we use ρ = E = 1, b = 0, values in (1.2). For (b) and (d), we use c = a = 1, ρ = 1, σ = 1, and b = ϵ = 1 values in Example 11. Solutions are generated using symbolic computation. The infinite series in (5.2) is truncated after adding 46 terms. No difference has been observed visually if more terms are added.
Next, we present the formal system of linear PD wave equations in n-space dimensions [Ref. 55, Eq. (54)], n ∈ ℕ∗,
Here, “⋅” indicates matrix multiplication, or equivalently by the system
where x ∈ ℝn, t ∈ ℝ, C : ℝn → M(n × n, ℝ) is the micromodulus tensor, assumed to be even and assuming values inside the subspace of symmetric matrices, ρ > 0 is the mass density, b : ℝn × ℝ → ℝn is the prescribed body force density, and u : ℝn × ℝ → ℝn is the displacement field.
For comparison, e.g., the corresponding wave equation in classical elasticity in 1-space dimension is given by
where E > 0 is the so called “Young’s modulus,” and describing compression waves in a rod.
If j, k ∈ {1, …, n} and Cjk ∈ L1(ℝn), we can rewrite (1.1) as
II. OPERATOR-THEORETIC TREATMENT OF SYSTEMS OF WAVE EQUATIONS
Analogous to the majority of evolution equations from classical and quantum physics, (1.3) can be treated with methods from operator theory, see, e.g., Refs. 11 and 47 for substantiation of this claim and Ref. 29 for applications of operator theory in engineering. More specifically, this system falls into the class of abstract linear wave equations from Theorem 1. For the proof of this theorem see, e.g., Ref. 11 [Theorem 2.2.1 and Corollary 2.2.2]. Special cases of this theorem are proved in Refs. 30 and 42 and Ref. 48 [Vol. II]. Statements and proofs make use of the spectral theorems of (densely defined, linear and) self-adjoint operators in Hilbert spaces, including the concept of functions of such operators, see, e.g., Ref. 48 [Vol. I], or standard books on functional analysis, such as Refs. 49 and 66. These methods are also used throughout the paper.
This section provides the basic properties of the solutions of abstract wave equations of the form (2.1). Theorem 1 gives the well-posedness of the initial value problem for a class of abstract wave equations, conservation of energy, and a representation of the solutions in terms of bounded functions of the governing operator. Theorem 3 provides special solutions of the associated class of inhomogeneous wave equations. Together with Theorem 1, these solutions provide the well-posedness of the initial value problem of the latter equations as well as a representation of the solutions in terms of bounded functions of the governing operator.
Let be some non-trivial complex Hilbert space. Furthermore, let A : D(A) → X be some densely defined, linear, semibounded self-adjoint operator in X with spectrum σ(A). Finally, let ξ, η ∈ D(A).
- Then there is a unique twice continuously differentiable map u : ℝ → X assuming values in D(A) and satisfyingfor all t ∈ ℝ as well as
- For this u, the corresponding energy function Eu : ℝ → ℝ, defined byfor all t ∈ ℝ, is constant.
- Moreover, this u is represented by the following solution operators for all t ∈ ℝ:70
See the proofs of Ref. 11 [Theorem 2.2.1 and Corollary 2.2.2]. □
Moreover, if A is positive, the solutions of (2.1) are stable, i.e., there are no solutions that are growing exponentially in the norm.
The proof is obvious. □
Duhamel’s principle leads to a solution of (2.1) for vanishing data, the proof of the well-posedness and a representation of the solutions of the initial value problem of the inhomogeneous equation,
For simplicity, the corresponding subsequent Theorem 3 assumes that A is in addition positive, which is the most relevant case for applications because otherwise there are exponentially growing solutions, indicating that the system is unstable. The same statement is true if σ(A) is only bounded from below. On the other hand, Theorem 3 can also be obtained by application of the corresponding well-known more general theorem for strongly continuous semigroups; see, e.g., Ref. 11 [Theorem 4.6.2]. We give a direct proof of Theorem 3 in Section VI D, which does not rely on methods from the theory of strongly continuous semigroups. For the definition of weak integration; see, e.g., Ref. 11 [Sec. 3.2].
See Section VI D. □
III. THE GOVERNING OPERATOR AND PROPERTIES
The standard data space for the classical wave equation is a L2-space with constant weight, on a non-empty open subset of ℝn, n ∈ ℕ∗, for instance, for a bar of infinite extension in 1-space dimension. It turns out that the classical data spaces are suitable also as data spaces for peridynamics, for instance, again for a bar of infinite extension in 1-space dimension, composed of a “linear peridynamic material.” This simplifies the discussion of the convergence of peridynamic solutions to classical solutions.
In the following, we represent (1.3) in the form of (2.1), where the governing operator A is an “operator matrix,” consisting of sums of multiples of the identity and convolution operators, as indicated in (1.3). These matrix entries will turn out to be pairwise commuting. The following remark provides some known relevant information on operator matrices of bounded operators. On the other hand, we avoid explicit matrix notation.
If 𝕂 ∈ {ℝ, ℂ}, (X, 〈|〉) a non-trivial 𝕂-Hilbert space, (Ajk)j,k∈{1,…,n} a family of elements of L(X, X).
- Then byfor every (ξ1, …, ξn) ∈ Xn, there is defined a bounded linear operator with adjoint A∗ given byfor every (ξ1, …, ξn) ∈ Xn.
- If the members of (Ajk)j,k∈{1,…,n} are pairwise commuting, then A is bijective if and only if det(A) is bijective, whereSn denotes the set of permutations of {1, …, n},for all σ ∈ Sn and sign denotes the signum function.
The basic properties of the entries of the operator matrix are given in the following lemma. In fact, these operators turn out to be bounded linear operators on . Hence, the boundedness and self-adjointness of A follows from those of AC. The boundedness of A has appeared in various forms5,6,8,18,23,24,67 and sometimes for special class of kernel functions. We give a result using kernel functions that are in L1(ℝn) by utilizing a well-known criterion for integral operators; see, e.g., corollary to Ref. 64 [Theorem 6.24].
By a consequence of a well-known criterion for integral operators on L2-spaces, see, e.g., Corollary to Ref. 64 [Theorem 6.24]. □
For other governing operators related to AC on bounded domains, see Ref. 4.
IV. MAIN RESULTS
The following are the three main results in this paper:
Main Result 1: Nonlocal operator is a bounded function of the classical operator.
Main Result 2: Strong convergence of nonlocal solutions to classical ones through strong resolvent convergence.
Main Result 3: Representation of the solution in terms of spherical Bessel functions.
Our first main result establishes the connection between the nonlocal operator and the classical one. We identify the exact expression of this connection. Namely, the governing operator AC in (3.1) of the peridynamic wave equation is a bounded function of the classical governing operator in (1.2). This observation enables the comparison of peridynamic solutions to those of classical elasticity. In addition, it enables the generalization of peridynamic-type operators from ℝn to bounded domains as functions of the corresponding classical operator. This is the subject of our companion papers.2,3
Employing the notion of strong resolvent convergence, our second main result gives the strong convergence of solutions of the governing equation to that of the classical equation. There are results that establish the convergence of the peridynamic operator to the classical operator with vanishing nonlocality; see Refs. 6, 7, 24, and 56. More important for applications is the corresponding convergence of solutions. This question has also been pointed out in Ref. 24 [p. 862] for the convergence of PD solutions to that of Navier equation. The studies8,18,67 considered convergence of solutions. In particular, Refs. 39 and 40 put an effort to establish strong convergence of solutions. One can find a survey of existing convergence results in Ref. 21 [Chap. 4]. The tool that we develop for this purpose is the notion of strong resolvent convergence64,65 used in Theorem 7. To the best of our knowledge, our study provides the first exploitation of strong resolvent convergence of the functions of the classical operator to obtain strong convergence of solutions.
As an application, we establish strong convergence of solutions generated by standard examples of micromoduli Cν in the literature, given in (6.3) and (6.4); see Lemmata 16 and 18.
Our third main result is the discovery that the solution can be represented in terms of a series of Bessel functions. This property can be exploited in obtaining numerical representation of the solutions. We utilized this for the solutions depicted in Section V.
V. EXAMPLES AND ILLUSTRATIONS OF MAIN RESULT 3
As an application of main result 3, we construct examples of the nonlocal wave equation to illustrate the separation of waves and propagation of discontinuities, Examples 10 and 11, respectively. The depicted solutions in Figures 1 and 2 were computed symbolically using spherical Bessel function representation of the solutions. The infinite series in (5.1) and (5.2), respectively, are truncated after adding 46 terms. Visually, we do not observe a difference in the solutions if more terms are added.
For constructing the case that characterizes separation of waves, we choose an example that involves a micromodulus and input function both of which are normal distributions with mean value zero and standard deviation σ and σd, respectively. Similar micromodulus functions have also been used in Refs. 41, 63, and 68.
We depict and compare the solutions of the classical and nonlocal wave equations in Figures 1 and 2. In the classical case, as expected, we observe the propagation of waves along characteristics; see Figures 1(a) and 1(c) for vanishing initial velocity and displacement, respectively. In the nonlocal case, we observe repeated separation of waves and an oscillation at the center of the initial pulse; see Figures 1(b) and 1(d) for vanishing initial velocity and displacement, respectively. Solution wave patterns were also reported earlier23,25 for the same micromodulus function in Example 11.
Now, we study the propagation of discontinuity in the data for classical and nonlocal wave equations in the following example.
In the classical wave equation, as expected, discontinuities propagate along the characteristics; see Figures 2(a) and 2(c) for vanishing initial velocity and displacement, respectively. On the other hand, in the nonlocal case, the discontinuity remains in the same place for all time; see Figures 2(b) and 2(d) for vanishing initial velocity and displacement, respectively. This confirms the results given in Ref. 63.
VI. PROOFS AND RELATED RESULTS
A. Proof of main result 1 (nonlocal operator is a bounded function of the classical operator)
For the study of the spectral properties of the matrix entries, needed for the application of the results from Section II, we use Fourier transformations. This step parallels the common procedure for constant coefficient differential operators on ℝn, n ∈ ℕ∗. With the help of the unitary Fourier transform F2, Theorem 13 represents the matrix entries as maximal multiplication operators. This process can be viewed as a form of “diagonalization” of the entries. Also, since bounded maximal multiplication operators commute, the entries commute pairwise. The spectra of maximal multiplication operators are well understood, leading to Corollary 14. Also, the functional calculus which is associated to maximal multiplication operators is known and allows the construction of the functional calculi of the entries. The latter is used in the proof of Theorem 6 which proves that matrix entries corresponding to spherically symmetric micromoduli are functions of the Laplace operator.
We finally arrive at the proof of main result 1.
We give the spectrum and point spectrum of AC. The spectrum is an essential ingredient in constructing functional calculus for AC; see Lemma 20.
The result is a consequence of well-known properties of multiplication operators. □
B. Proof of main result 2 (strong convergence of nonlocal solutions to classical ones through strong resolvent convergence)
Lemma 15 gives conditions for the convergence of bounded functions of a self-adjoint operator to converge to that operator, which implies strong resolvent convergence and also the strong convergence of the same bounded continuous function of each member of the sequence against that bounded continuous function of the self-adjoint operator [Ref. 64, Theorem 9.16]; see Theorem 7.
Now, main result 2, Theorem 7, is a consequence of Lemma 15, and, for example, Ref. 48 [Vol. I, Theorems 8.20 and 8.25].
We apply our main result 2 to standard examples treated in the literature. Hence, we easily obtain strong resolvent convergence for these cases, i.e., Lemmas 16 and 18. They provide sequences of micromoduli which satisfy the conditions of Lemma 15. Lemma 16 has also been treated in Refs. 41 and 61 and Lemma 18 has been treated in Ref. 41. Remark 8 applies Theorem 7 to the sequences of micromoduli from Lemmas 16 and 18. As a consequence, for fixed data and t ∈ ℝ, the solutions of the initial value problem at time t corresponding to the members of each sequence of micromoduli converge in to the corresponding classical solution at time t.
By applying Theorem 7 as in Remark 8, we obtain that the solutions of the nonlocal operator ACν strongly converge to that of the local operator for the micromoduli Cν in (6.3).
By applying Theorem 7 as in Remark 8, we obtain that the solutions of the nonlocal operator ACν strongly converge to that of the local operator for the micromoduli Cν in (6.4).
C. Proof of main result 3 (representation of the solution in terms of spherical Bessel functions)
In this section, the goal is to prove that the solutions of the nonlocal wave equation given in (2.2) can be expressed in terms of spherical Bessel functions. This result has a practical implication because it allows symbolic computation of solutions. Indeed, using symbolic computation, we generate solutions of the nonlocal wave equation and illustrate how separation of waves and propagation of discontinuities occur; see Figures 1 and 2. In addition, one can use the explicit representations of solutions for benchmarking numerical results which helps in the development of numerical methods for computing approximate solutions. Such benchmarking could be used, for instance, for verifying and validating numerical solutions.
Main result 1 indicates that the governing peridynamic operator is bounded. Hence, functions of that operator in (2.2) can be represented in the form of power series in the governing operator due to functional calculus of self-adjoint and bounded operators. In Lemma 20, we prove that holomorphic functional calculus is available for power series representation. In Lemma 21, we apply this representation to the functions present in the solution of the initial value problem of the homogeneous wave equation. Lemmas 20 and 21 can be viewed as straightforward applications of the spectral theorems for densely defined, self-adjoint linear operators in Hilbert spaces.
Now, we give the exact representation of the terms present in (2.2) in terms of power series in A.
The goal here was to show that (2.2) can be expressed by spherical Bessel functions. In order to arrive at a Bessel function representation, we first connect (2.2) to generalized hypergeometric functions. In Theorem 22, we give a general result which involves two arbitrary bounded commuting operators A and B. In Theorem 24, we will apply this result to special commuting operators, multiple of the identity operator, and a convolution, i.e., c − C, because the governing operator is of this form. In addition, we expect that the expansion below can be used for large time asymptotic of the solutions of the nonlocal wave equation, a research direction beyond the scope of this paper.
We give a connection between generalized hypergeometric and spherical Bessel functions which will lead to main result 3.
Eventually, we have a representation involving two commuting operators, with one of the operators being a multiple of the identity and a general C which is not necessarily a convolution operator.
Direct consequence of Theorem 22 and Lemma 23. □
Now, main result 3, Theorem 9, is a consequence of Theorem 24. The representations given in Theorems 22, 24, and 9 were used to express the solutions in Examples 10 and 11.
As a side result, we provide an error estimate for the representation in Theorem 24. For the solution plots, Figures 1 and 2, the infinite series is truncated after adding 46 terms, i.e., N = 45. The below error estimate has been used to monitor the error. For N = 45 and t = 20, the error is in the order of 10−13.
D. Proof of Theorem 3 (solutions of inhomogeneous wave equations)
We give a proof of Theorem 3.
VII. CONCLUSION
Our result that the governing operator is a bounded function of the classical local operator for scalar-valued functions should be generalizable to vector-valued case. Our notable result that the governing operator AC of the peridynamic wave equation is a bounded function of the classical governing operator has far reaching consequences. It enables the comparison of peridynamic solutions to those of classical elasticity. The remarkable implication is that it opens the possibly of defining peridynamic-type operators on bounded domains as functions of the corresponding classical operator. Since the classical operator is defined through local boundary conditions, the functions inherit this knowledge. This observation opens a gateway to incorporate local boundary conditions into nonlocal theories, which has vital implications for numerical treatment of nonlocal problems. This is the subject of our companion papers.2,3
We expect that the expansions in Theorems 22 and 24 can be used for obtaining the large time asymptotic of solutions of the nonlocal wave equation. In the classical case, as expected, we observe the propagation of waves along characteristics. In the nonlocal case, we observe oscillatory recurrent wave separation. We think that this phenomenon is worth investigating. On the other hand, we observe that discontinuity remains stationary in the nonlocal case, whereas, it is well-known that discontinuities propagate along characteristics. We hold that this fundamental difference is one of the most distinguishing features of PD. In conclusion, we believe that we added valuable tools to the arsenal of methods to analyze nonlocal problems.
Acknowledgments
Burak Aksoylu was supported in part by National Science Foundation Grant No. DMS 1016190, European Commission Marie Curie Career Integration Grant No. 293978, and Scientific and Technological Research Council of Turkey (TÜBİTAK) Grant Nos. TBAG 112T240 and MAG 112M891. Research visit of Horst R. Beyer was supported in part by No. TÜBİTAK 2221 Fellowship for Visiting Scientist Program. Sabbatical visit of Fatih Celiker was supported in part by No. TÜBİTAK 2221 Fellowship for Scientist on Sabbatical Leave Program. Fatih Celiker was supported in part by National Science Foundation Grant No. DMS 1115280.
REFERENCES
Throughout the paper, we refer to the Laplace operator as the local or the classical operator interchangeably.
Solution operators, and denote operators corresponding to unique extensions to entire holomorphic functions. A pedantic notation would be and , respectively. For brevity, we prefer to use the former.
denotes the space of bounded complex-valued functions on σ(A) that are strongly measurable in the sense that they are everywhere on σ(A) limit of a sequence of step functions.