On a class of nonlocal wave equations from applications

Beyer, Horst Reinhard; Aksoylu, Burak; Celiker, Fatih

doi:10.1063/1.4953252

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 ℝⁿ, 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,⁵⁵ 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 dynamics^52,54 and has been demonstrated as a viable multiscale material model for length scales ranging from molecular dynamics to classical elasticity.¹⁰ Also see other related engineering applications,^{13,33,35,46,45} the review, and news articles^14,16,21,57 for a comprehensive discussion, and the book.³⁸

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 studies^{6–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 L¹(ℝⁿ), 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.⁶⁹ 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.³² 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.

FIG. 1.

View large Download slide

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.

FIG. 2.

View large Download slide

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 ∈ ℕ^∗,

ρ \frac{\partial^{2} u}{\partial t^{2}} (x, t) = \int_{R^{n}} C (x^{'} - x) \cdot (u (x^{'}, t) - u (x, t)) d x^{'} + b (x, t) .

Here, “⋅” indicates matrix multiplication, or equivalently by the system

ρ \frac{\partial^{2} u_{j}}{\partial t^{2}} (x, t) = \sum_{k = 1}^{n} \int_{R^{n}} C_{j k} (x^{'} - x) \cdot (u_{k} (x^{'}, t) - u_{k} (x, t)) d x^{'} + b_{j} (x, t),

(1.1)

where x ∈ ℝⁿ, t ∈ ℝ, C : ℝⁿ → 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 : ℝⁿ × ℝ → ℝⁿ is the prescribed body force density, and u : ℝⁿ × ℝ → ℝⁿ is the displacement field.

For comparison, e.g., the corresponding wave equation in classical elasticity in 1-space dimension is given by

ρ \frac{\partial^{2} u}{\partial t^{2}} = E \frac{\partial^{2} u}{\partial x^{2}} + b,

(1.2)

where E > 0 is the so called “Young’s modulus,” and describing compression waves in a rod.

If j, k ∈ {1, …, n} and C_jk ∈ L¹(ℝⁿ), we can rewrite (1.1) as

ρ \frac{\partial^{2} u_{j}}{\partial t^{2}} (x, t) = - \sum_{k = 1}^{n} \{[\int_{R^{n}} C_{j k} (x^{'}) d x^{'}] u_{k} (x, t) - (C_{j k} * u_{k} (\cdot, t)) (x)\} + b_{j} (x, t),

(1.3)

for all x ∈ ℝⁿ, t ∈ ℝ, and j ∈ {1, …, n} where ∗ denotes the convolution product. The system (1.3) is the starting point for a functional analytic interpretation, which leads on a well-posed initial value problem. For this purpose, we use methods from operator theory; see, e.g., Refs. 11 and 47.

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.

Theorem 1

Wave Equations

Let $(X, 〈 | 〉)$ 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 satisfying

$u^{″} (t) = - A u (t)$
(2.1)
for all t ∈ ℝ as well as

$u (0) = ξ, u^{'} (0) = η .$
For this u, the corresponding energy function E_u : ℝ → ℝ, defined by

$E_{u} (t) ≔ \frac{1}{2} () 〈 u^{'} (t) | u^{'} (t) 〉 + 〈 u (t) | A u (t) 〉 ()$

for all t ∈ ℝ, is constant.
Moreover, this u is represented by the following solution operators for all t ∈ ℝ:⁷⁰

$u (t) = cos (t \sqrt{A}) ξ + \frac{sin (t \sqrt{A})}{\sqrt{A}} η .$
(2.2)

Proof.

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.

Corollary 2

Stability of solutions

If A is positive, then for every t ∈ ℝ, we have

‖u (t)‖ ⩽ ‖ξ‖ + |t| \cdot ‖η‖ .

Proof.

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,

u^{″} (t) = - A u (t) + b (t), t \in R .

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].

Theorem 3

Solutions of inhomogeneous wave equations

Let

(X, 〈 | 〉)

, A : D(A) → X, σ(A) be as in Theorem 1 and, in addition, A be positive. Let b : ℝ → X be a continuous map, assuming values in D(A²) such that Ab, A² b are continuous. Then, v : ℝ → X, for every t ∈ ℝ defined by

\begin{align} v (t) ≔ \int_{I_{t}} \frac{sin ((t - τ) \sqrt{A})}{\sqrt{A}} b (τ) d τ, \end{align}

where ∫ denotes weak integration in X,

I_{t} ≔ {\begin{cases} [0, t] & if t ⩾ 0 \\ [t, 0] & if t < 0 \end{cases},

is twice continuously differentiable, assumes values in D(A), is such that

v^{″} (t) + A v (t) = b (t), t \in R,

v (0) = v^{'} (0) = 0 .

Proof.

See Section VI D. □

III. THE GOVERNING OPERATOR AND PROPERTIES

The standard data space for the classical wave equation is a L²-space with constant weight, on a non-empty open subset of ℝⁿ, n ∈ ℕ^∗, for instance, $L_{C}^{2} (R)$ 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 $L_{C}^{2} (R)$ 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.

Remark 4

Operator matrices

If 𝕂 ∈ {ℝ, ℂ}, (X, 〈|〉) a non-trivial 𝕂-Hilbert space, (A_jk)_{j,k∈{1,…,n}} a family of elements of L(X, X).

Then by

$A (ξ_{1}, \dots, ξ_{n}) ≔ (\sum_{k = 1}^{n} A_{1 k} ξ_{k}, \dots, \sum_{k = 1}^{n} A_{n k} ξ_{k})$

for every (ξ₁, …, ξ_n) ∈ Xⁿ, there is defined a bounded linear operator with adjoint A^∗ given by

$A^{*} (ξ_{1}, \dots, ξ_{n}) = (\sum_{k = 1}^{n} A_{k 1}^{*} ξ_{k}, \dots, \sum_{k = 1}^{n} A_{k n}^{*} ξ_{k})$

for every (ξ₁, …, ξ_n) ∈ Xⁿ.
If the members of (A_jk)_{j,k∈{1,…,n}} are pairwise commuting, then A is bijective if and only if det(A) is bijective, where

$\begin{align} det (A) ≔ \sum_{σ \in S_{n}} sign (σ) A_{1 σ (1)} \dots A_{n σ (n)}, \end{align}$

S_n denotes the set of permutations of {1, …, n},

$sign (σ) ≔ \prod_{i, j = 1, i < j}^{n} sign (σ (j) - σ (i))$

for all σ ∈ S_n 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 $L_{C}^{2} (R^{n})$ ⁠. Hence, the boundedness and self-adjointness of A follows from those of A_C. The boundedness of A has appeared in various forms^{5,6,8,18,23,24,67} and sometimes for special class of kernel functions. We give a result using kernel functions that are in L¹(ℝⁿ) by utilizing a well-known criterion for integral operators; see, e.g., corollary to Ref. 64 [Theorem 6.24].

Lemma 5

Boundedness of matrix entries

Let n ∈ ℕ^∗, ρ > 0 and C ∈ L¹(ℝⁿ) be even. Then,

A_{C} f ≔ \frac{1}{ρ} [(\int_{R^{n}} C d v^{n}) . f - C * f],

(3.1)

for every

f \in L_{C}^{2} (R^{n})

, where ∗ denotes the convolution product, there is defined a self-adjoint bounded linear operator on

L_{C}^{2} (R^{n})

with operator norm

‖A_{C}‖

satisfying

‖A_{C}‖ ⩽ \frac{1}{ρ} (|\int_{R} C d v^{n}| + {‖C‖}_{1}) ⩽ \frac{2 {‖C‖}_{1}}{ρ} .

Proof.

By a consequence of a well-known criterion for integral operators on L²-spaces, see, e.g., Corollary to Ref. 64 [Theorem 6.24]. □

For other governing operators related to A_C 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 A_C 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 ℝⁿ to bounded domains as functions of the corresponding classical operator. This is the subject of our companion papers.^2,3

Theorem 6

Main result 1

Let n ∈ ℕ^∗,

L_{n}

be the closure of the positive symmetric, essentially self-adjoint operator in

L_{C}^{2} (R^{n})

, given by

(C_{0}^{\infty} (R^{n}, C) \to L_{C}^{2} (R^{n}), f \mapsto - \frac{E}{ρ} △ f),

where ρ > 0 and E > 0. In addition, for n > 1, let C be spherically symmetric, i.e., such that

C \circ R = C,

for every R ∈ SO(n), where SO(n) denotes the map of group of special orthogonal transformations on ℝⁿ. Then

A_{C} = \{) \frac{1}{ρ} [(F_{1} C) (0) - F_{1} C] \circ ι (\} (L_{n}),

where ι : [0, ∞) → ℝⁿ is defined by

ι (s) ≔ (\sqrt{\frac{ρ}{E} s}) . e_{1},

for every s ⩾ 0 and e₁, …, e_n denotes the canonical basis of ℝⁿ.

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 studies^8,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 convergence^64,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.

Theorem 7

Main result 2

Let (X, 〈|〉) be a non-trivial complex Hilbert space and A : D(A) → X a densely defined, linear, and self-adjoint operator with spectrum σ(A). Furthermore, let f₁, f₂, … be a sequence of real-valued functions in

U_{C}^{s} (σ (A))

⁷¹ that is everywhere on σ(A) pointwise convergent to id_σ(A), and for which there is M > 0 such that

|f_{ν}| ⩽ M [(1 + ||) |_{σ (A)}]

for all ν ∈ ℝ. Then for every g ∈ BC(ℝ, ℂ), where BC(ℝⁿ, ℂ) is the space of complex-valued bounded continuous functions on ℝⁿ,

s - lim_{ν \to \infty} [g |_{σ (f_{ν} (A))}] (f_{ν} (A)) = [g |_{σ (A)}] (A),

where for every ν ∈ ℕ^∗, σ(f_ν(A)) denotes the spectrum of f_ν(A).

Remark 8.

The mechanism to obtain strong convergence of nonlocal solutions to classical ones is as follows. We first choose a sequence of micromoduli C_ν so that

f_{C_{ν}} (λ) ≔ F_{1} (C_{ν}) (0) - F_{1} (C_{ν}) (λ)

(4.1)

is pointwise convergent to

{id}_{σ (L_{1})}

and that satisfies condition (6.1), where

L_{1}

is the classical governing operator in 1 dimension defined in Theorem 6. We choose

A = L_{1}

and f_{C_ν} in (4.1). Then, a well-known result [Ref. 64 , Theorem 9.16] indicates that Lemma 15 implies strong resolvent convergence of

f_{C_{ν}} (L_{1})

to

L_{1}

. Then, Theorem 7 implies that for every g ∈ BC(ℝ, ℂ), we obtain the following strong convergence:

s - lim_{ν \to \infty} [g |_{σ (f_{C_{ν}} (L_{1}))}] (f_{C_{ν}} (L_{1})) = [g |_{σ (L_{1})}] (L_{1}) .

(4.2)

We denote the nonlocal operator corresponding to micromoduli C_ν by

A_{C_{ν}} ≔ f_{C_{ν}} (L_{1})

. Then, we apply (4.2) to solution operators g in (2.2), namely,

cos (t \sqrt{A_{C_{ν}}})

and

sin (t \sqrt{A_{C_{ν}}}) / \sqrt{A_{C_{ν}}}

, both of which are in BC(ℝ, ℂ). Consequently, we obtain strong convergence of nonlocal solution operators

[g |_{σ (f_{C_{ν}} (L_{1})}] (f_{C_{ν}} (L_{1}))

in (3.1) to the classical solution operators

[g |_{σ (L_{1})}] (L_{1}) .

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.

Theorem 9

Main result 3

Let n ∈ ℕ^∗, ρ > 0, C ∈ L¹(ℝⁿ) be even such that

c ≔ \int_{R^{n}} C d v^{n} > 0 .

Then, for A_C in (3.1) and t ∈ ℝ, the solutions have the following representation:

\begin{aligned} cos (t \sqrt{A_{C}}) f & = \sum_{k = 0}^{\infty} \frac{1}{2^{k} k!} {(\sqrt{c t^{2} / ρ})}^{k + 1} j_{k - 1} (\sqrt{c t^{2} / ρ}) c^{- k} . C^{k} * f, \\ \frac{sin (t \sqrt{A_{C}})}{\sqrt{A_{C}}} f & = t \sum_{k = 0}^{\infty} \frac{1}{2^{k} k!} {(\sqrt{c t^{2} / ρ})}^{k} j_{k} (\sqrt{c t^{2} / ρ}) c^{- k} . C^{k} * f, \end{aligned}

(4.3)

for every

f \in L_{C}^{2} (R^{n})

, where the spherical Bessel functions j₀, j₁, … are defined as in Ref. 43,

j_{- 1} (z) ≔ \frac{cos (z)}{z}, z \in C^{*},

and the members of the sums are defined for t = 0 by continuous extension.

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.

Example 10.

For ρ, σ, σ_d, a > 0, we define C_σ ∈ L¹(ℝ) and

f \in L_{C}^{2} (R)

by

\begin{align} C_{σ} ≔ \frac{a}{\sqrt{2 π} σ} e^{- [1 / (2 σ^{2})] . {id}_{R}^{2}}, f ≔ \frac{1}{\sqrt{2 π} σ_{d}} e^{- [1 / (2 σ_{d}^{2})] . {id}_{R}^{2}} . \end{align}

Then for k ∈ ℕ^∗

\begin{align} F_{1} C_{σ} = a e^{- (σ^{2} / 2) . {id}_{R}^{2}}, F_{1} C_{σ}^{k} = {(F_{1} C_{σ})}^{k} = a^{k} e^{- k (σ^{2} / 2) . {id}_{R}^{2}}, \\ F_{2} (C_{σ}^{k} * f) = (F_{1} C_{σ}^{k}) \cdot F_{2} f = \frac{a^{k}}{\sqrt{2 π}} e^{- k (σ^{2} / 2) . {id}_{R}^{2}} \cdot e^{- (σ_{d}^{2} / 2) . {id}_{R}^{2}} \\ = \frac{a^{k}}{\sqrt{2 π}} e^{- [(k σ^{2} + σ_{d}^{2}) / 2] . {id}_{R}^{2}} = \frac{1}{\sqrt{2 π}} F_{1} \frac{a^{k}}{\sqrt{2 π} \sqrt{k σ^{2} + σ_{d}^{2}}} e^{- {1 / [2 (k σ^{2} + σ_{d}^{2})]} . {id}_{R}^{2}} \\ = F_{2} \frac{a^{k}}{\sqrt{2 π} \sqrt{k σ^{2} + σ_{d}^{2}}} e^{- {1 / [2 (k σ^{2} + σ_{d}^{2})]} . {id}_{R}^{2}} \end{align}

and hence

C_{σ}^{k} * f = \frac{a^{k}}{\sqrt{2 π} \sqrt{k σ^{2} + σ_{d}^{2}}} e^{- {1 / [2 (k σ^{2} + σ_{d}^{2})]} . {id}_{R}^{2}} .

Since

c = \int_{R} C_{σ} d v^{1} = a > 0,

we conclude from Theorem 9 that for A_C in (3.1) and t ∈ ℝ,

\begin{aligned} cos (t \sqrt{A_{C}}) f & = \sum_{k = 0}^{\infty} \frac{1}{2^{k} k!} {(\sqrt{a t^{2} / ρ})}^{k + 1} j_{k - 1} (\sqrt{a t^{2} / ρ}) \frac{1}{\sqrt{2 π} \sqrt{k σ^{2} + σ_{d}^{2}}} e^{- {1 / [2 (k σ^{2} + σ_{d}^{2})]} . {id}_{R}^{2}}, \\ \frac{sin (t \sqrt{A_{C}})}{\sqrt{A_{C}}} f & = t \sum_{k = 0}^{\infty} \frac{1}{2^{k} k!} {(\sqrt{a t^{2} / ρ})}^{k} j_{k} (\sqrt{a t^{2} / ρ}) \frac{1}{\sqrt{2 π} \sqrt{k σ^{2} + σ_{d}^{2}}} e^{- {1 / [2 (k σ^{2} + σ_{d}^{2})]} . {id}_{R}^{2}} . \end{aligned}

(5.1)

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 earlier^23,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.

Example 11.

As in the previous example, for ρ, σ, a, b, ε > 0, we define C_σ ∈ L¹(ℝ) by

\begin{align} C_{σ} ≔ \frac{a}{\sqrt{2 π} σ} e^{- [1 / (2 σ^{2})] . {id}_{R}^{2}} . \end{align}

Then

F_{1} C_{σ} = a e^{- (σ^{2} / 2) . {id}_{R}^{2}},

and for k ∈ ℕ^∗,

F_{1} C_{σ}^{k} = {(F_{1} C_{σ})}^{k} = a^{k} e^{- k (σ^{2} / 2) . {id}_{R}^{2}} = F_{1} \frac{a^{k}}{\sqrt{2 π k} σ} e^{- [1 / (2 k σ^{2})] . {id}_{R}^{2}},

and hence

C_{σ}^{k} = \frac{a^{k}}{\sqrt{2 π} \sqrt{k σ^{2}}} e^{- [1 / (2 k σ^{2})] . {id}_{R}^{2}} .

Furthermore, we define

f \in L_{C}^{2} (R)

by

\begin{align} f ≔ b e^{- ε . {id}_{R}} \cdot χ_{_{[0, \infty)}} . \end{align}

Then for x ∈ ℝ,

\begin{align} (C_{σ}^{k} * f) (x) = \frac{a^{k} b}{\sqrt{2 π} \sqrt{k σ^{2}}} \int_{0}^{\infty} e^{- [{(x - y)}^{2} / (2 k σ^{2})]} \cdot e^{- ε y} d y \\ = \frac{2 a^{k} b}{π} e^{{(\frac{ε σ}{2} \sqrt{2 k})}^{2}} e^{- ε x} erfc (\frac{ε σ}{2} \sqrt{2 k} - \frac{x}{σ \sqrt{2 k}}), \end{align}

where erfc denotes the error function defined according to DLMF.⁴³ We note for x ∈ ℝ that

lim_{k \to 0} \frac{a^{k} b}{2} e^{{(\frac{ε σ}{2} \sqrt{2 k})}^{2}} e^{- ε x} erfc (\frac{ε σ}{2} \sqrt{2 k} - \frac{x}{σ \sqrt{2 k}}) = {\begin{cases} 0 & if x < 0 \\ \frac{b}{2} & if x = 0 \\ b e^{- ε x} & if x > 0 \end{cases} .

Since

c = \int_{R} C_{σ} d v^{1} = a > 0,

we conclude from Theorem 9 that for t ∈ ℝ,

\begin{aligned} cos (t \sqrt{A_{C}}) f & = \sum_{k = 0}^{\infty} \frac{1}{2^{k} k!} {(\sqrt{a t^{2} / ρ})}^{k + 1} j_{k - 1} (\sqrt{a t^{2} / ρ}) f_{k}, \\ \frac{sin (t \sqrt{A_{C}})}{\sqrt{A_{C}}} f & = t \sum_{k = 0}^{\infty} \frac{1}{2^{k} k!} {(\sqrt{a t^{2} / ρ})}^{k} j_{k} (\sqrt{a t^{2} / ρ}) f_{k}, \end{aligned}

(5.2)

for A_C in (3.1), where for every x ∈ ℝ and k ∈ ℕ^∗,

\begin{align} f_{0} & ≔ b e^{- ε . {id}_{R}} \cdot χ_{_{[0, \infty)}}, \\ f_{k} (x) & ≔ \frac{2 b}{π} e^{{(\frac{ε σ}{2} \sqrt{2 k})}^{2}} e^{- ε x} erfc (\frac{ε σ}{2} \sqrt{2 k} - \frac{x}{σ \sqrt{2 k}}) \\ = \frac{b}{\sqrt{2 π} \sqrt{k σ^{2}}} e^{- ε x} \int_{- \infty}^{x} e^{- u^{2} / (2 k σ^{2})} \cdot e^{ε u} d u . \end{align}

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 ∈ ℕ^∗. With the help of the unitary Fourier transform F₂, 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.

Assumption 12.

In the following, for n ∈ ℕ^∗, F₂ denotes the unitary Fourier transformation on

L_{C}^{2} (R^{n})

which, for every rapidly decreasing test function f ∈ 𝒮_ℂ(ℝ), is defined by

(F_{2} f) (k) ≔ \frac{1}{{(2 π)}^{n / 2}} \int_{R^{n}} e^{- i k \cdot {id}_{R^{n}}} f d v^{n}, k \in R^{n} .

Also, we denote by F₁ the map from

L_{C}^{1} (R^{n})

to C_∞(ℝⁿ, ℂ), the space of continuous functions on ℝⁿ vanishing at infinity, which for every

f \in L_{C}^{1} (R^{n})

, is defined by

(F_{1} f) (k) ≔ \int_{R^{n}} e^{- i k \cdot {id}_{R^{n}}} f d v^{n}, k \in R^{n} .

Theorem 13

Fourier transforms of the entries

Let

T_{\frac{1}{ρ} [(F_{1} C) (0) - F_{1} C]}

denote the maximal multiplication operator by the bounded continuous function

\frac{1}{ρ} [(F_{1} C) (0) - F_{1} C]

on

L_{C}^{2} (R^{n})

. Then

F_{2} \circ A_{C} \circ F_{2}^{- 1} = T_{\frac{1}{ρ} [(F_{1} C) (0) - F_{1} C]} .

Proof.

First, we define an operator K ≔ C∘(p₁ − p₂), with projections p₁, p₂ : ℝ²ⁿ → ℝⁿ, where

p_{1} (x_{1}, \dots, x_{n}, y_{1}, \dots, y_{n}) ≔ (x_{1}, \dots, x_{n}), p_{2} (x_{1}, \dots, x_{n}, y_{1}, \dots, y_{n}) ≔ (y_{1}, \dots, y_{n}),

for (x₁, …, x_n, y₁, …, y_n) ∈ ℝ²ⁿ. Since C is even, K is symmetric. Furthermore, for every x ∈ ℝⁿ and y ∈ ℝⁿ,

K (x, \cdot) = C (x - \cdot) = C (\cdot - x), K (\cdot, y) = C (\cdot - y) \in L^{1} (R^{n}) .

K induces a self-adjoint bounded linear integral operator Int(K) on

L_{C}^{2} (R^{n})

defined by

[Int (K) f] (x) ≔ \int_{R^{n}} K (x, \cdot) \cdot f d v^{n} = \int_{R^{n}} C (x - \cdot) \cdot f d v^{n} = (C * f) (x) .

We note that for every

L_{C}^{1} (R^{n}) \cap L_{C}^{2} (R^{n})

⁠,

\begin{align} [F_{2} \circ Int (K)] f & = F_{2} (C * f) = \frac{1}{{(2 π)}^{n / 2}} . F_{1} (C * f) = \frac{1}{{(2 π)}^{n / 2}} . (F_{1} C) (F_{1} f) \\ = (F_{1} C) (F_{2} f) = [T_{F_{1} C} \circ F_{2}] f . \end{align}

Since

L_{C}^{1} (R^{n}) \cap L_{C}^{2} (R^{n})

is dense in

L_{C}^{2} (R^{n})

⁠, the bounded linear operators F₂∘Int(K) and T_F₁C∘F₂ coincide on a dense subspace of

L_{C}^{2} (R^{n})

⁠, they coincide on the whole of

L_{C}^{2} (R^{n})

⁠. Consequently, the result follows from

[F_{2} \circ Int (K)] f = [T_{F_{1} C} \circ F_{2}] f, f \in L_{C}^{2} (R^{n}) .

□

We finally arrive at the proof of main result 1.

Proof.

First, we note that

F_{2} \circ L_{n} \circ F_{2}^{- 1} = T_{\frac{E}{ρ} {||}^{2}},

where

T_{\frac{E}{ρ} {||}^{2}}

denotes the maximal multiplication operator in

L_{C}^{2} (R^{n})

by the function

\frac{E}{ρ} {||}^{2}

⁠. In particular, this implies that the spectrum of

L_{n}

⁠,

σ (L_{n})

⁠, is given by [0, ∞) and for every

g \in U_{C}^{s} ([0, \infty))

that

g (L_{n}) = F_{2}^{- 1} \circ T_{g \circ (\frac{E}{ρ} {||}^{2})} \circ F_{2},

where

T_{g \circ (\frac{E}{ρ} {||}^{2})}

denotes the maximal multiplication operator on

L_{C}^{2} (R^{n})

by the function

g \circ (\frac{E}{ρ} {||}^{2}) .

Furthermore, we note that (F₁C)(0) − F₁C ∈ BC(ℝⁿ, ℝ) and that (F₁C)(0) − F₁C is even, since for every k ∈ ℝⁿ

(F_{1} C) (- k) = \int_{R^{n}} e^{i k \cdot {id}_{R^{n}}} C d v^{n} = \int_{R^{n}} e^{- i k \cdot {id}_{R^{n}}} [C \circ (- {id}_{R^{n}})] d v^{n}

= \int_{R^{n}} e^{- i k \cdot {id}_{R^{n}}} C d v^{n} = (F_{1} C) (k),

(F_{1} C) (k) = \frac{1}{2} [\int_{R^{n}} e^{- i k \cdot {id}_{R^{n}}} C d v^{n} + \int_{R^{n}} e^{i k \cdot {id}_{R^{n}}} C d v^{n}]

= \int_{R^{n}} cos (k \cdot {id}_{R^{n}}) C d v^{n},

(F_{1} C) (0) - (F_{1} C) (k) = \int_{R^{n}} [1 - cos (k \cdot {id}_{R^{n}})] C d v^{n} = 2 \int_{R^{n}} {sin}^{2} (\frac{k}{2} \cdot {id}_{R^{n}}) C d v^{n} .

Furthermore for n > 1, we note that

(F_{1} C) (R (k)) = \int_{R^{n}} e^{- i R (k) \cdot {id}_{R^{n}}} C d v^{n} = \int_{R^{n}} e^{- i R (k) \cdot R} (C \circ R) d v^{n}

= \int_{R^{n}} e^{- i k \cdot {id}_{R^{n}}} (C \circ R) d v^{n} = \int_{R^{n}} e^{- i k \cdot {id}_{R^{n}}} C d v^{n} = (F_{1} C) (k)

for every R ∈ SO(n) and k ∈ ℝⁿ and hence that

(F_{1} C) (k) = (F_{1} C) (|k| . e_{1})

for every k ∈ ℝⁿ. In particular,

\begin{align} \frac{1}{ρ} [(F_{1} C) (0) - F_{1} C] \circ ι \in U_{R}^{s} ([0, \infty)) \end{align}

and

\begin{align} \{) \frac{1}{ρ} [(F_{1} C) (0) - F_{1} C] \circ ι (\} (L_{n}) = F_{2}^{- 1} \circ T_{\{) \frac{1}{ρ} [(F_{1} C) (0) - F_{1} C] \circ ι (\} \circ (\frac{E}{ρ} {||}^{2})} \circ F_{2} \\ = F_{2}^{- 1} \circ T_{\frac{1}{ρ} [(F_{1} C) (0) - F_{1} C] \circ (|| . e_{1})} \circ F_{2} = F_{2}^{- 1} \circ T_{\frac{1}{ρ} [(F_{1} C) (0) - F_{1} C]} \circ F_{2} = A_{C} . \end{align}

□

We give the spectrum and point spectrum of A_C. The spectrum is an essential ingredient in constructing functional calculus for A_C; see Lemma 20.

Corollary 14

Spectral Properties of A_C

σ (A_{C}) = \bar{Ran \frac{1}{ρ} . [(F_{1} C) (0) - F_{1} C]},

σ_{p} (A_{C}) = \{λ \in R : \{k \in R : \frac{1}{ρ} . [(F_{1} C) (0) - (F_{1} C) (k)] = λ\} is no Lebesgue null set\},

where the overline denotes the closure in ℝ. Finally, for every λ ∈ σ(A_C), A_C − λ is not surjective.

Proof.

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.

Lemma 15

Convergence of bounded functions of a self-adjoint operator to that operator

Let (X, 〈|〉) be a non-trivial complex Hilbert space and A : D(A) → X a densely defined, linear, and self-adjoint operator with spectrum σ(A). Furthermore, let f₁, f₂, … be a sequence in

U_{C}^{s} (σ (A))

that is everywhere on σ(A) pointwise convergent to id_σ(A), and for which there is M > 0 such that

|f_{ν}| ⩽ M [(1 + ||) |_{σ (A)}]

(6.1)

for all ν ∈ ℝ. Then

lim_{ν \to \infty} f_{ν} (A) ξ = A ξ, ξ \in D (A) .

(6.2)

Proof.

Let ξ ∈ D(A) and ψ_ξ the corresponding spectral measure. According to the spectral theorem for densely defined, self-adjoint linear operators in Hilbert spaces,

{id}_{R}^{2}

is ψ_ξ-summable and

\begin{align} {‖f_{μ} (A) ξ - f_{ν} (A) ξ‖}^{2} = {‖(f_{μ} - f_{ν}) (A) ξ‖}^{2} \\ = 〈 (f_{μ} - f_{ν}) (A) ξ | (f_{μ} - f_{ν}) (A) ξ 〉 = 〈 ξ | {|f_{μ} - f_{ν}|}^{2} (A) ξ 〉 \\ = \int_{σ (A)} {|f_{μ} - f_{ν}|}^{2} d ψ_{ξ} = {‖f_{μ} - f_{ν}‖}_{2, ψ_{ξ}}^{2} = {‖f_{μ} - {id}_{σ (A)} + {id}_{σ (A)} - f_{ν}‖}_{2, ψ_{ξ}}^{2} \\ ⩽ {({‖f_{μ} - {id}_{σ (A)}‖}_{2, ψ_{ξ}} + {‖{id}_{σ (A)} - f_{ν}‖}_{2, ψ_{ξ}})}^{2}, \end{align}

for μ, ν ∈ ℕ^∗. As a consequence of the pointwise convergence of f₁, f₂, … on σ(A) to id_σ(A), and Lebesgue’s dominated convergence theorem, it follows that

lim_{μ \to \infty} {‖f_{μ} - {id}_{σ (A)}‖}_{2, ψ_{ξ}} = 0

and hence that f₁(A) ξ, f₂(A) ξ, … is a Cauchy sequence in X. Since

(X, ‖‖)

is in particular complete, the latter implies that f₁(A) ξ, f₂(A) ξ, … is convergent in

(X, ‖‖)

⁠. Furthermore,

\begin{align} 〈 ξ | lim_{ν \to \infty} f_{ν} (A) ξ 〉 = lim_{ν \to \infty} 〈 ξ | f_{ν} (A) ξ 〉 = lim_{ν \to \infty} \int_{σ (A)} f_{ν} d ψ_{ξ} \\ = \int_{σ (A)} {id}_{σ (A)} d ψ_{ξ} = 〈 ξ | A ξ 〉, \end{align}

where again the pointwise convergence of f₁, f₂, … on σ(A) to id_σ(A), , Lebesgue’s dominated convergence theorem and the spectral theorem for densely defined, self-adjoint linear operators in Hilbert spaces has been applied. From the polarization identity for 〈|〉, it follows that

〈 ξ | lim_{ν \to \infty} f_{ν} (A) η 〉 = 〈 ξ | A η 〉

for all ξ, η ∈ D(A). Since D(A) is dense in X, the latter implies that

〈 ξ | lim_{ν \to \infty} f_{ν} (A) η 〉 = 〈 ξ | A η 〉

for all ξ ∈ X, η ∈ D(A) and hence for every η ∈ D(A) that

lim_{ν \to \infty} f_{ν} (A) η = A η .

□

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 $L_{C}^{2} (R)$ to the corresponding classical solution at time t.

Lemma 16.

For every ν ∈ ℕ^∗, we define C_ν ∈ L¹(ℝ) by

C_{ν} ≔ 3 E ν^{3} χ_{_{[- \frac{1}{ν}, \frac{1}{ν}]}} .

(6.3)

Then, the result (6.2) preparing for strong convergence of solutions is satisfied

lim_{ν \to \infty} \{\frac{1}{ρ} [(F_{1} C_{ν}) (0) - F_{1} C] \circ ι\} (L_{1}) f = L_{1} f, f \in D (L_{1}) = W_{C}^{2} (R) .

Remark 17.

By applying Theorem 7 as in Remark 8, we obtain that the solutions of the nonlocal operator A_{C_ν} strongly converge to that of the local operator for the micromoduli C_ν in (6.3).

Proof.

For ν ∈ ℕ^∗,

F_{1} C_{ν} = 6 E ν^{3} \frac{sin (ν^{- 1} . {id}_{R})}{{id}_{R}} .

Furthermore, for ν ∈ ℕ^∗, λ ⩾ 0,

\begin{align} \frac{1}{ρ} [(F_{1} C_{ν}) (0) - F_{1} C] \circ ι (λ) = \frac{1}{ρ} [(F_{1} C_{ν}) (0) - F_{1} C_{ν}] (\sqrt{\frac{ρ}{E} λ}) \\ = \frac{6 E ν^{2}}{ρ} [1 - \frac{sin (ν^{- 1} . {id}_{R})}{ν^{- 1} . {id}_{R}}] (\sqrt{\frac{ρ}{E} λ}) \end{align}

and k > 0

\begin{align} 1 - \frac{sin (k / ν)}{k / ν} = ν \int_{0}^{1 / ν} [1 - cos (k x)] d x = \int_{0}^{1} [1 - cos (k u / ν)] d u \\ = \int_{0}^{1} [\int_{0}^{k / ν} u sin (u y) d y] d u = \frac{k}{ν} \int_{0}^{1} [\int_{0}^{1} u sin (k u v / ν) d v] d u \\ = \frac{k^{2}}{ν^{2}} \int_{{[0, 1]}^{2}} u^{2} v \frac{sin (k u v / ν)}{k u v / ν} d u d v \end{align}

and hence that

ν^{2} [1 - \frac{sin (k / ν)}{k / ν}] = k^{2} \int_{{[0, 1]}^{2}} u^{2} v \frac{sin (k u v / ν)}{k u v / ν} d u d v .

From the latter, we conclude with the help of Lebesgue’s dominated convergence theorem that

lim_{ν \to \infty} ν^{2} [1 - \frac{sin (k / ν)}{k / ν}] = \frac{k^{2}}{6}

as well as that

|ν^{2} [1 - \frac{sin (k / ν)}{k / ν}]| ⩽ k^{2} \int_{{[0, 1]}^{2}} u^{2} v |\frac{sin (k u v / ν)}{k u v / ν}| d u d v ⩽ k^{2} \int_{{[0, 1]}^{2}} u^{2} v d u d v = \frac{k^{2}}{6} .

In particular, we conclude for λ ⩾ 0 that

lim_{ν \to \infty} \frac{1}{ρ} [(F_{1} C_{ν}) (0) - F_{1} C] \circ ι (λ) = \frac{6 E ν^{2}}{ρ} \cdot \frac{ρ λ}{6 E ν^{2}} = λ

as well as that

\begin{align} |\frac{1}{ρ} [(F_{1} C_{ν}) (0) - F_{1} C] \circ ι (λ)| ⩽ \frac{6 E}{ρ} \cdot \frac{1}{6} \frac{ρ λ}{E} = λ . \end{align}

Finally, we conclude the result from Lemma 15. □

Lemma 18.

For every ν ∈ ℕ^∗, we define C_ν ∈ L¹(ℝ) by

C_{ν} ≔ \frac{2 E ν^{3}}{\sqrt{2 π}} e^{- (ν^{2} / 2) . {id}_{R}^{2}} = 2 E ν^{2} \cdot \frac{ν}{\sqrt{2 π}} e^{- (ν^{2} / 2) . {id}_{R}^{2}} .

(6.4)

Then, the result (6.2) preparing for strong convergence of solutions is satisfied,

lim_{ν \to \infty} \{\frac{1}{ρ} [(F_{1} C_{ν}) (0) - F_{1} C] \circ ι\} (L_{1}) f = L_{1} f, f \in D (L_{1}) = W_{C}^{2} (R) .

Remark 19.

By applying Theorem 7 as in Remark 8, we obtain that the solutions of the nonlocal operator A_{C_ν} strongly converge to that of the local operator for the micromoduli C_ν in (6.4).

Proof.

For ν ∈ ℕ^∗, λ ⩾ 0,

\begin{align} F_{1} C_{ν} = 2 E ν^{2} \cdot e^{- [1 / (2 ν^{2})] . {id}_{R}^{2}}, \\ \frac{1}{ρ} [(F_{1} C_{ν}) (0) - F_{1} C_{ν}] \circ ι (λ) = \frac{1}{ρ} [(F_{1} C_{ν}) (0) - F_{1} C_{ν}] (\sqrt{\frac{ρ}{E} λ}) \\ = \frac{2 E ν^{2}}{ρ} \{1 - e^{- [1 / (2 ν^{2})] . {id}_{R}^{2}}\} (\sqrt{\frac{ρ}{E} λ}) \end{align}

and k ⩾ 0

\begin{align} ν^{2} [1 - e^{- k^{2} / (2 ν^{2})}] = ν^{2} \int_{0}^{k^{2} / (2 ν^{2})} e^{- u} d u = \int_{0}^{k^{2} / 2} e^{- v / ν^{2}} d v . \end{align}

From the latter, we conclude for k ⩾ 0, with the help of Lebesgue’s dominated convergence theorem, that

lim_{ν \to \infty} ν^{2} [1 - e^{- k^{2} / (2 ν^{2})}] = \frac{k^{2}}{2}

as well as that

|ν^{2} [1 - e^{- k^{2} / (2 ν^{2})}]| ⩽ \frac{k^{2}}{2}

and hence for λ ⩾ 0 that

\begin{align} lim_{ν \to \infty} \frac{1}{ρ} [(F_{1} C_{ν}) (0) - F_{1} C_{ν}] \circ ι (λ) = \frac{2 E}{ρ} \frac{ρ}{2 E} λ = λ, \\ |\frac{1}{ρ} [(F_{1} C_{ν}) (0) - F_{1} C_{ν}] \circ ι (λ)| ⩽ \frac{2 E}{ρ} \frac{ρ}{2 E} λ = λ . \end{align}

Finally, we conclude the result from Lemma 15. □

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.

Lemma 20

Holomorphic functional calculus

Let (X, 〈|〉) be a non-trivial complex Hilbert space, A ∈ L(X, X) self-adjoint, and σ(A) ⊂ ℝ the (non-empty, compact) spectrum of A. Furthermore, let

R > ‖A‖

and f : U_R(0) → ℂ be holomorphic. Then, the sequence

{(\frac{f^{(k)} (0)}{k!} . A^{k})}_{k \in N}

is absolutely summable in L(X, X) and

(f |_{σ (A)}) (A) = \sum_{k = 0}^{\infty} \frac{f^{(k)} (0)}{k!} . A^{k} .

Proof.

First, we note that according to Taylor’s theorem, general properties of power series, and the compactness of σ(A) that

{(\frac{f^{(k)} (0)}{k!} . z^{k})}_{k \in N}

is absolutely summable for every z ∈ U_R(0) as well as, since

σ (A) \subset B_{‖A‖} (0) \subset U_{R} (0),

that the sequence of continuous functions

{(\sum_{k = 0}^{n} \frac{f^{(k)} (0)}{k!} . {({id}_{R} |_{σ (A)})}^{n})}_{n \in N}

converges uniformly to the continuous function f|_σ(A). In particular, since

‖A‖ < R

⁠, this implies that the sequence

{(\frac{f^{(k)} (0)}{k!} . A^{k})}_{k \in N}

is absolutely summable in L(X, X), and it follows from the spectral theorem for bounded self-adjoint operators in Hilbert spaces that

(\sum_{k = 0}^{n} \frac{f^{(k)} (0)}{k!} . {({id}_{R} |_{σ (A)})}^{n}) (A) = \sum_{k = 0}^{n} \frac{f^{(k)} (0)}{k!} . A^{k},

as well as that

(f |_{σ (A)}) (A) = \sum_{k = 0}^{\infty} \frac{f^{(k)} (0)}{k!} . A^{k} .

□

Now, we give the exact representation of the terms present in (2.2) in terms of power series in A.

Lemma 21

Power series representation

Let (X, 〈|〉) be a non-trivial complex Hilbert space,

\sqrt{}

the complex square-root function, with domain ℂ∖((−∞, 0] × {0}). A ∈ L(X, X) self-adjoint and σ(A) ⊂ ℝ the (non-empty, compact) spectrum of A. For every t ∈ ℝ, the sequences

{({(- 1)}^{k} \frac{t^{2 k}}{(2 k)!} . A^{k})}_{k \in N}, {({(- 1)}^{k} \frac{t^{2 k + 1}}{(2 k + 1)!} . A^{k})}_{k \in N}

are absolutely summable in L(X, X) and

\begin{align} cos (t \sqrt{A}) = \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{t^{2 k}}{(2 k)!} A^{k}, \\ \frac{sin (t \sqrt{A})}{\sqrt{A}} = \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{t^{2 k + 1}}{(2 k + 1)!} A^{k} . \end{align}

Proof.

We note that for every t ∈ ℝ,

\begin{align} cos (t \sqrt{}) : C ∖ ((- \infty, 0] \times {0}) \to C, \\ cosh (t \sqrt{} \circ (- {id}_{C})) : C ∖ ([0, \infty) \times {0}) \to C \end{align}

are holomorphic functions such that

\begin{align} cos (t \sqrt{z}) = \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{{(t \sqrt{z})}^{2 k}}{(2 k)!} = \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{t^{2 k}}{(2 k)!} z^{k}, \\ cosh (t \sqrt{- z}) = \sum_{k = 0}^{\infty} \frac{{(t \sqrt{- z})}^{2 k}}{(2 k)!} = \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{t^{2 k}}{(2 k)!} z^{k} \end{align}

for every z ∈ ℂ∖((−∞, 0] × {0}) and z ∈ ℂ∖[0, ∞) × {0}), respectively. Furthermore,

\begin{align} \frac{sin (t \sqrt{})}{\sqrt{}} : C ∖ ((- \infty, 0] \times {0}) \to C, \\ \frac{sinh (t \sqrt{})}{\sqrt{}} \circ (- {id}_{C}) : C ∖ ([0, \infty) \times {0}) \to C \end{align}

are holomorphic functions such that

\begin{align} \frac{sin (t \sqrt{z})}{\sqrt{z}} = \frac{1}{\sqrt{z}} \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{{(t \sqrt{z})}^{2 k + 1}}{(2 k + 1)!} = \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{t^{2 k + 1}}{(2 k + 1)!} z^{k}, \\ \frac{sinh (t \sqrt{- z})}{\sqrt{- z}} = \frac{1}{\sqrt{- z}} \sum_{k = 0}^{\infty} \frac{{(t \sqrt{- z})}^{2 k + 1}}{(2 k + 1)!} = \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{t^{2 k + 1}}{(2 k + 1)!} z^{k} \end{align}

for every z ∈ ℂ∖((−∞, 0] × {0}) and z ∈ ℂ∖[0, ∞) × {0}), respectively. Then, it follows from Lemma 20 that the sequences

{({(- 1)}^{k} \frac{t^{2 k}}{(2 k)!} . A^{k})}_{k \in N}, {({(- 1)}^{k} \frac{t^{2 k + 1}}{(2 k + 1)!} . A^{k})}_{k \in N}

are absolutely summable in L(X, X) and that

\begin{align} cos (t \sqrt{A}) & = \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{t^{2 k}}{(2 k)!} A^{k}, \\ \frac{sin (t \sqrt{A})}{\sqrt{A}} & = \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{t^{2 k + 1}}{(2 k + 1)!} A^{k} . □ \end{align}

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.

Theorem 22.

Let (X, 〈|〉) be a non-trivial complex Hilbert space,

\sqrt{}

the complex square-root function, with domain ℂ∖((−∞, 0] × {0}). A, B ∈ L(X, X) self-adjoint such that [A, B] = 0 and σ(A), σ(A + B) ⊂ ℝ the (non-empty, compact) spectra of A and A + B, respectively. Then

\begin{align} cos (t \sqrt{A + B}) & = \sum_{k = 0}^{\infty} {(- 1)}^{k} \cdot \frac{t^{2 k}}{(2 k)!} . \{[_{0} F_{1} (-; k + \frac{1}{2}; - \frac{t^{2}}{4} . {i d}_{σ (A)})] (A)\} B^{k}, \\ \frac{sin (t \sqrt{A + B})}{\sqrt{A + B}} & = \sum_{k = 0}^{\infty} {(- 1)}^{k} \cdot \frac{t^{2 k + 1}}{(2 k + 1)!} . \{[_{0} F_{1} (-; k + \frac{3}{2}; - \frac{t^{2}}{4} . {i d}_{σ (A)})] (A)\} B^{k}, \end{align}

where ₀F₁ denotes the generalized hypergeometric function, defined as in Ref. 43 .

Proof.

In a first step, we note for every t ∈ ℝ that the family

{({(- 1)}^{k + l} (\binom{k + l}{l}) \frac{t^{2 (k + l)}}{[2 (k + l)]!} A^{k} B^{l})}_{(k, l) \in N^{2}}

is absolutely summable in L(X, X), since for (k, l) ∈ ℕ²

\begin{align} ‖{(- 1)}^{k + l} (\binom{k + l}{l}) \frac{t^{2 (k + l)}}{[2 (k + l)]!} A^{k} B^{l}‖ ⩽ \frac{t^{2 (k + l)}}{[2 (k + l)]!} (\binom{k + l}{l}) {‖A‖}^{k} {‖B‖}^{l} \\ = \frac{(k + l)!}{l! k! [2 (k + l)]!} {(t^{2} ‖A‖)}^{k} {(t^{2} ‖B‖)}^{l} ⩽ \frac{1}{k!} {(t^{2} ‖A‖)}^{k} \frac{1}{l!} {(t^{2} ‖B‖)}^{l} \end{align}

and hence for every finite subset S ⊂ ℕ²,

\begin{align} \sum_{(k, l) \in S} ‖{(- 1)}^{k + l} (\binom{k + l}{l}) \frac{t^{2 (k + l)}}{[2 (k + l)]!} A^{k} B^{l}‖ ⩽ exp (t^{2} ‖A‖) exp (t^{2} ‖B‖) . \end{align}

Also, we note that the family

{({(- 1)}^{k + l} (\binom{k + l}{l}) \frac{t^{2 (k + l) + 1}}{[2 (k + l) + 1]!} A^{k} B^{l})}_{(k, l) \in N^{2}}

is absolutely summable in L(X, X), since for (k, l) ∈ ℕ²,

\begin{align} {‖{(- 1)}^{k + l} (\binom{k + l}{l}) \frac{t^{2 (k + l) + 1}}{[2 (k + l) + 1]!} A^{k} B^{l}‖}_{2} ⩽ \frac{{|t|}^{2 (k + l) + 1}}{[2 (k + l) + 1]!} (\binom{k + l}{l}) {‖A‖}^{k} {‖B‖}^{l} \\ = |t| \frac{(k + l)!}{l! k! [2 (k + l) + 1]!} {(t^{2} ‖A‖)}^{k} {(t^{2} ‖B‖)}^{l} ⩽ |t| \frac{1}{k!} {(t^{2} ‖A‖)}^{k} \frac{1}{l!} {(t^{2} ‖B‖)}^{l}, \end{align}

leading to

\begin{align} \sum_{(k, l) \in S} ‖{(- 1)}^{k + l} (\binom{k + l}{l}) \frac{t^{2 (k + l) + 1}}{[2 (k + l) + 1]!} A^{k} B^{l}‖ ⩽ |t| exp (t^{2} ‖A‖) exp (t^{2} ‖B‖), \end{align}

for every finite subset S ⊂ ℕ². Hence, we conclude the following:

\begin{align} \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{t^{2 k}}{(2 k)!} {(A + B)}^{k} = \sum_{k = 0}^{\infty} \sum_{l = 0}^{k} {(- 1)}^{k} \frac{t^{2 k}}{(2 k)!} (\binom{k}{l}) A^{k - l} B^{l} \\ = \sum_{l = 0}^{\infty} \sum_{k = l}^{\infty} {(- 1)}^{k} \frac{t^{2 k}}{(2 k)!} (\binom{k}{l}) A^{k - l} B^{l} = \sum_{l = 0}^{\infty} [\sum_{k = l}^{\infty} {(- 1)}^{k} (\binom{k}{l}) \frac{t^{2 k}}{(2 k)!} A^{k - l}] B^{l} \\ = \sum_{l = 0}^{\infty} [\sum_{k = 0}^{\infty} {(- 1)}^{k + l} (\binom{k + l}{l}) \frac{t^{2 (k + l)}}{[2 (k + l)]!} A^{k}] B^{l} \\ = \sum_{l = 0}^{\infty} {(- 1)}^{l} t^{2 l} [\sum_{k = 0}^{\infty} {(- 1)}^{k} (\binom{k + l}{l}) \frac{t^{2 k}}{[2 (k + l)]!} A^{k}] B^{l} \\ = \sum_{l = 0}^{\infty} {(- 1)}^{l} t^{2 l} [\sum_{k = 0}^{\infty} \frac{(k + l)!}{[2 (k + l)]! \cdot l!} \cdot \frac{1}{k!} {(- t^{2} A)}^{k}] B^{l} . \end{align}

In the following, we show the auxiliary result that for every k, l ∈ ℕ,

\frac{(k + l)!}{[2 (k + l)]! \cdot l!} = 2^{- k} \cdot \frac{1}{\prod_{m = 0}^{k - 1} [2 (l + m) + 1]} \cdot \frac{1}{(2 l)!} .

(6.5)

The proof proceeds by induction on k. First, we note that

\frac{l!}{(2 l)! \cdot l!} = 2^{- 0} \cdot \frac{1}{\prod_{m = 0}^{- 1} [2 (l + m) + 1]} \cdot \frac{1}{(2 l)!} .

In the following, we assume that is true for some k ∈ ℕ. Then

\begin{align} \frac{(k + l + 1)!}{[2 (k + l + 1)]! \cdot l!} = \frac{1}{2} \cdot \frac{1}{2 (k + l) + 1} \cdot \frac{(k + l)!}{[2 (k + l)]! \cdot l!} \\ = \frac{1}{2} \cdot \frac{1}{2 (k + l) + 1} \cdot 2^{- k} \cdot \frac{1}{\prod_{m = 0}^{k - 1} [2 (l + m) + 1]} \cdot \frac{1}{(2 l)!} \\ = 2^{- (k + 1)} \cdot \frac{1}{\prod_{m = 0}^{k} [2 (l + m) + 1]} \cdot \frac{1}{(2 l)!}, \end{align}

and hence is true for k + 1. The equality implies for every k, l ∈ ℕ that

\begin{align} \frac{(k + l)!}{[2 (k + l)]! \cdot l!} = 2^{- k} \cdot \frac{1}{\prod_{m = 0}^{k - 1} [2 (l + m) + 1]} \cdot \frac{1}{(2 l)!} \\ = 4^{- k} \cdot \frac{1}{\prod_{m = 0}^{k - 1} (l + m + \frac{1}{2})} \cdot \frac{1}{(2 l)!} = 4^{- k} \cdot \frac{1}{(2 l)! {(l + \frac{1}{2})}_{k}} . \end{align}

Hence

\begin{align} \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{t^{2 k}}{(2 k)!} {(A + B)}^{k} = \sum_{l = 0}^{\infty} {(- 1)}^{l} t^{2 l} [\sum_{k = 0}^{\infty} \frac{1}{(2 l)! {(l + \frac{1}{2})}_{k}} \cdot \frac{1}{k!} {(- \frac{t^{2}}{4} . A)}^{k}] B^{l} \\ \sum_{l = 0}^{\infty} {(- 1)}^{l} \frac{t^{2 l}}{(2 l)!} [\sum_{k = 0}^{\infty} \frac{1}{{(l + \frac{1}{2})}_{k} \cdot k!} \cdot {(- \frac{t^{2}}{4} . A)}^{k}] B^{l} . \end{align}

By definition of the generalized hypergeometric function ₀F₁, for every l ∈ ℕ, z ∈ ℂ,

\begin{align} _{0} F_{1} (-; l + \frac{1}{2}; z) = \sum_{k = 0}^{\infty} \frac{z^{k}}{{(l + \frac{1}{2})}_{k} \cdot k!} . \end{align}

Hence

\begin{align} \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{t^{2 k}}{(2 k)!} {(A + B)}^{k} = \sum_{l = 0}^{\infty} {(- 1)}^{l} \cdot \frac{t^{2 l}}{(2 l)!} . \{[_{0} F_{1} (-; l + \frac{1}{2}; - \frac{t^{2}}{4} . {i d}_{σ (A)})] (A)\} B^{l} . \end{align}

Furthermore,

\begin{align} \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{t^{2 k + 1}}{(2 k + 1)!} {(A + B)}^{k} = \sum_{k = 0}^{\infty} \sum_{l = 0}^{k} {(- 1)}^{k} \frac{t^{2 k + 1}}{(2 k + 1)!} (\binom{k}{l}) A^{k - l} B^{l} \\ = \sum_{l = 0}^{\infty} \sum_{k = l}^{\infty} {(- 1)}^{k} \frac{t^{2 k + 1}}{(2 k + 1)!} (\binom{k}{l}) A^{k - l} B^{l} = \sum_{l = 0}^{\infty} [\sum_{k = l}^{\infty} {(- 1)}^{k} (\binom{k}{l}) \frac{t^{2 k + 1}}{(2 k + 1)!} A^{k - l}] B^{l} \\ = \sum_{l = 0}^{\infty} [\sum_{k = 0}^{\infty} {(- 1)}^{k + l} (\binom{k + l}{l}) \frac{t^{2 (k + l) + 1}}{[2 (k + l) + 1]!} A^{k}] B^{l} \\ = \sum_{l = 0}^{\infty} {(- 1)}^{l} t^{2 l + 1} [\sum_{k = 0}^{\infty} {(- 1)}^{k} (\binom{k + l}{l}) \frac{t^{2 k}}{[2 (k + l) + 1]!} A^{k}] B^{l} \\ = t \sum_{l = 0}^{\infty} {(- 1)}^{l} \cdot t^{2 l} [\sum_{k = 0}^{\infty} \frac{(k + l)!}{[2 (k + l) + 1]! \cdot l!} \cdot \frac{1}{k!} . {(- t^{2} A)}^{k}] B^{l} . \end{align}

Since for every k, l ∈ ℕ,

\begin{align} \frac{(k + l)!}{[2 (k + l) + 1]! \cdot l!} = \frac{1}{2} \cdot \frac{1}{k + l + \frac{1}{2}} \cdot \frac{(k + l)!}{[2 (k + l)]! \cdot l!} \\ = \frac{1}{2} \cdot \frac{1}{k + l + \frac{1}{2}} \cdot 4^{- k} \cdot \frac{1}{(2 l)! {(l + \frac{1}{2})}_{k}} = \frac{1}{2} \cdot 4^{- k} \cdot \frac{1}{(2 l)! (l + \frac{1}{2}) {(l + \frac{3}{2})}_{k}} \\ = 4^{- k} \cdot \frac{1}{(2 l + 1)! {(l + \frac{3}{2})}_{k}}, \end{align}

we conclude that

\begin{align} \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{t^{2 k + 1}}{(2 k + 1)!} {(A + B)}^{k} \\ = t \sum_{l = 0}^{\infty} {(- 1)}^{l} \cdot t^{2 l} [\sum_{k = 0}^{\infty} \frac{1}{(2 l + 1)! {(l + \frac{3}{2})}_{k}} \cdot \frac{1}{k!} . {(- \frac{t^{2}}{4} . A)}^{k}] B^{l} \\ = \sum_{l = 0}^{\infty} {(- 1)}^{l} \cdot \frac{t^{2 l + 1}}{(2 l + 1)!} [\sum_{k = 0}^{\infty} \frac{1}{{(l + \frac{3}{2})}_{k} \cdot k!} . {(- \frac{t^{2}}{4} . A)}^{k}] B^{l} \\ = \sum_{l = 0}^{\infty} {(- 1)}^{l} \cdot \frac{t^{2 l + 1}}{(2 l + 1)!} . \{[_{0} F_{1} (-; l + \frac{3}{2}; - \frac{t^{2}}{4} . {i d}_{σ (A)})] (A)\} B^{l} . \end{align}

□

We give a connection between generalized hypergeometric and spherical Bessel functions which will lead to main result 3.

Lemma 23.

For every k ∈ ℕ and x > 0,

\begin{align} \frac{x^{2 k}}{(2 k)!} ._{0} F_{1} (-; k + \frac{1}{2}; - \frac{x^{2}}{4}) & = \frac{1}{2^{k} k!} x^{k + 1} j_{k - 1} (x), \\ \frac{x^{2 k + 1}}{(2 k + 1)!} \cdot_{0} F_{1} (-; k + \frac{3}{2}, - \frac{x^{2}}{4}) & = \frac{1}{2^{k} k!} x^{k + 1} j_{k} (|x|), \end{align}

where the spherical Bessel functions j₀, j₁, … are defined as in Ref. 43 and

j_{- 1} (x) ≔ \frac{cos (x)}{x}, x > 0 .

Proof.

We note that for every ν ∈ (0, ∞), k ∈ ℕ, and x > 0,

\begin{align} J_{ν} (x) & ≔ {(\frac{x}{2})}^{ν} \sum_{k = 0}^{\infty} \frac{{(- 1)}^{k}}{k! Γ (ν + k + 1)} {(\frac{x^{2}}{4})}^{k} = \frac{1}{Γ (ν + 1)} \cdot {(\frac{x}{2})}^{ν} \sum_{k = 0}^{\infty} \frac{{(- 1)}^{k}}{k! \frac{Γ (ν + k + 1)}{Γ (ν + 1)}} {(\frac{x^{2}}{4})}^{k} \\ = \frac{1}{Γ (ν + 1)} \cdot {(\frac{x}{2})}^{ν} \sum_{k = 0}^{\infty} \frac{{(- 1)}^{k}}{k! {(ν + 1)}_{k}} {(\frac{x^{2}}{4})}^{k} = \frac{1}{Γ (ν + 1)} \cdot {(\frac{x}{2})}^{ν} \cdot_{0} F_{1} (-; ν + 1, - x^{2} / 4) . \end{align}

\begin{align} j_{k} (x) & ≔ \sqrt{\frac{π}{2 x}} J_{k + \frac{1}{2}} (x) = \sqrt{\frac{π}{2 x}} \frac{1}{Γ (k + \frac{3}{2})} \cdot {(\frac{x}{2})}^{k + \frac{1}{2}} \cdot_{0} F_{1} (-; k + \frac{3}{2}, - x^{2} / 4) \\ = \frac{\sqrt{π}}{2 Γ (k + \frac{3}{2})} \cdot {(\frac{x}{2})}^{k} \cdot_{0} F_{1} (-; k + \frac{3}{2}, - x^{2} / 4) . \end{align}

Hence for every k ∈ ℕ, x > 0

\begin{align} _{0} F_{1} (-; k + \frac{3}{2}, - x^{2} / 4) = \frac{2 Γ (k + \frac{3}{2})}{\sqrt{π}} {(\frac{x}{2})}^{- k} j_{k} (x) = 2^{k + 1} {(\frac{1}{2})}_{k + 1} x^{- k} j_{k} (x) \end{align}

as well as

\begin{align} \frac{x^{2 k + 1}}{(2 k + 1)!}_{0} F_{1} (-; k + \frac{3}{2}, - x^{2} / 4) = \frac{x^{2 k + 1}}{(2 k + 1)!} 2^{k + 1} {(\frac{1}{2})}_{k + 1} x^{- k} j_{k} (x) \\ = \frac{x^{2 k + 1}}{(2 k + 1)!} 2^{k + 1} 2^{- (k + 1)} \frac{(2 k + 2)!}{2^{k + 1} (k + 1)!} x^{- k} j_{k} (x) = \frac{(2 k + 2)}{2^{k + 1} (k + 1)!} x^{k + 1} j_{k} (x) \\ = \frac{1}{2^{k} k!} x^{k + 1} j_{k} (x) . \end{align}

Furthermore, for every k ∈ ℕ^∗, x > 0

\frac{x^{2 k}}{(2 k)!} ._{0} F_{1} (-; k + \frac{1}{2}; - x^{2} / 4) = \frac{x}{2 k} \frac{1}{2^{k - 1} (k - 1)!} x^{k} j_{k - 1} (x) = \frac{1}{2^{k} k!} x^{k + 1} j_{k - 1} (x) .

(6.6)

Since for x > 0

\begin{align} _{0} F_{1} (-; \frac{1}{2}, - x^{2} / 4) = \sum_{k = 0}^{\infty} \frac{1}{{(\frac{1}{2})}_{k} \cdot k!} \cdot {(- \frac{x^{2}}{4})}^{k} = \sum_{k = 0}^{\infty} \frac{1}{2^{- k} \cdot \frac{(2 k)!}{2^{k} \cdot k!} \cdot k!} \cdot {(- \frac{x^{2}}{4})}^{k} \\ = \sum_{k = 0}^{\infty} \frac{4^{k}}{(2 k)!} \cdot {(- \frac{x^{2}}{4})}^{k} = \sum_{k = 0}^{\infty} {(- 1)}^{k} \frac{x^{2 k}}{(2 k)!} = cos (x), \end{align}

the equality is true also for k = 0, if we define

j_{- 1} (x) ≔ \frac{cos (x)}{x} .

□

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.

Theorem 24.

Let (X, 〈|〉) be a non-trivial complex Hilbert space,

\sqrt{}

the complex square-root function, with domain ℂ∖((−∞, 0] × {0}). c > 0, C ∈ L(X, X) self-adjoint and σ(c − C) ⊂ ℝ the (non-empty, compact) spectrum of c − C. Then for every t ∈ ℝ,

\begin{align} cos (t \sqrt{c - C}) & = \sum_{k = 0}^{\infty} \frac{1}{2^{k} k!} {(\sqrt{c t^{2}})}^{k + 1} j_{k - 1} (\sqrt{c t^{2}}) {(\frac{1}{c} . C)}^{k}, \\ \frac{sin (t \sqrt{c - C})}{\sqrt{c - C}} & = t \sum_{k = 0}^{\infty} \frac{1}{2^{k} k!} {(\sqrt{c t^{2}})}^{k} j_{k} (\sqrt{c t^{2}}) {(\frac{1}{c} . C)}^{k}, \end{align}

where the spherical Bessel functions j₀, j₁, … are defined as in Ref. 43 and

j_{- 1} (x) ≔ \frac{cos (x)}{x}, x > 0

and the members of the sums are defined for t = 0 by continuous extension.

Proof.

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⁻¹³.

Corollary 25

Error estimates

Let

(X, 〈 | 〉), \sqrt{}, c, C, σ (c - C), j_{- 1}, j_{0}, j_{1}, \dots

as in Theorem 24 and N ∈ ℕ. Then for every t ∈ ℝ,

\begin{align} ‖cos (t \sqrt{c - C}) - \sum_{k = 0}^{N} \frac{1}{2^{k} k!} {(\sqrt{c t^{2}})}^{k + 1} j_{k - 1} (\sqrt{c t^{2}}) {(\frac{1}{c} . C)}^{k}‖ & ⩽ \frac{π}{N!} min \{1, {(\frac{t^{2} ‖C‖}{4})}^{N + 1}\} e^{t^{2} ‖C‖ / 4}, \\ ‖\frac{sin (t \sqrt{c - C})}{\sqrt{c - C}} - t \sum_{k = 0}^{\infty} \frac{1}{2^{k} k!} {(\sqrt{c t^{2}})}^{k} j_{k} (\sqrt{c t^{2}}) {(\frac{1}{c} . C)}^{k}‖ & ⩽ \frac{π}{2 (N + 1)!} |t| min \{1, {(\frac{t^{2} ‖C‖}{4})}^{N + 1}\} e^{t^{2} ‖C‖ / 4} . \end{align}

Proof.

As a consequence of Theorem 24, for t ∈ ℝ,

‖cos (t \sqrt{c - C}) - \sum_{k = 0}^{N} \frac{1}{2^{k} k!} {(\sqrt{c t^{2}})}^{k + 1} j_{k - 1} (\sqrt{c t^{2}}) {(\frac{1}{c} . C)}^{k}‖ ⩽ \sum_{k = N + 1}^{\infty} \frac{1}{2^{k} k!} {(\sqrt{c t^{2}})}^{k + 1} |j_{k - 1} (\sqrt{c t^{2}})| {(\frac{‖C‖}{c})}^{k} ⩽ \sum_{k = N + 1}^{\infty} \frac{1}{2^{k} k!} {(\sqrt{c t^{2}})}^{k + 1} π \frac{{(\sqrt{c t^{2}})}^{k - 1}}{2^{k} (k - 1)!} {(\frac{‖C‖}{c})}^{k} = π \sum_{k = N + 1}^{\infty} \frac{1}{k! (k - 1)!} {(\frac{t^{2} ‖C‖}{4})}^{k} ⩽ \frac{π}{N!} \sum_{k = N + 1}^{\infty} \frac{1}{k!} {(\frac{t^{2} ‖C‖}{4})}^{k},

where the integral representation DLMF 10.54.1 of Ref. 43 (http://dlmf.nist.gov/10.54) for spherical Bessel functions has been used. Since

\begin{align} \sum_{k = N + 1}^{\infty} \frac{1}{k!} {(\frac{t^{2} ‖C‖}{4})}^{k} = {(\frac{t^{2} ‖C‖}{4})}^{N + 1} \sum_{k = N + 1}^{\infty} \frac{1}{k!} {(\frac{t^{2} ‖C‖}{4})}^{k - N - 1} \\ ⩽ {(\frac{t^{2} ‖C‖}{4})}^{N + 1} \sum_{k = N + 1}^{\infty} \frac{1}{(k - N - 1)!} {(\frac{t^{2} ‖C‖}{4})}^{k - N - 1} ⩽ {(\frac{t^{2} ‖C‖}{4})}^{N + 1} e^{t^{2} ‖C‖ / 4}, \end{align}

this implies that

\begin{align} ‖cos (t \sqrt{c - C}) - \sum_{k = 0}^{N} \frac{1}{2^{k} k!} {(\sqrt{c t^{2}})}^{k + 1} j_{k - 1} (\sqrt{c t^{2}}) {(\frac{1}{c} . C)}^{k}‖ & ⩽ \frac{π}{N!} min \{1, {(\frac{t^{2} ‖C‖}{4})}^{N + 1}\} e^{t^{2} ‖C‖ / 4} . \end{align}

Furthermore,

\begin{align} ‖\frac{sin (t \sqrt{c - C})}{\sqrt{c - C}} - t \sum_{k = 0}^{N} \frac{1}{2^{k} k!} {(\sqrt{c t^{2}})}^{k} j_{k} (\sqrt{c t^{2}}) {(\frac{1}{c} . C)}^{k}‖ \\ ⩽ |t| \sum_{k = N + 1}^{\infty} \frac{1}{2^{k} k!} {(\sqrt{c t^{2}})}^{k} |j_{k} (\sqrt{c t^{2}})| {(\frac{‖C‖}{c})}^{k} \\ ⩽ |t| \sum_{k = N + 1}^{\infty} \frac{1}{2^{k} k!} {(\sqrt{c t^{2}})}^{k} π \frac{{(\sqrt{c t^{2}})}^{k}}{2^{k + 1} k!} {(\frac{‖C‖}{c})}^{k} \\ = \frac{π}{2} |t| \sum_{k = N + 1}^{\infty} \frac{1}{{(k!)}^{2}} {(\frac{t^{2} ‖C‖}{4})}^{k} ⩽ \frac{π}{2 (N + 1)!} |t| \sum_{k = N + 1}^{\infty} \frac{1}{k!} {(\frac{t^{2} ‖C‖}{4})}^{k} \\ ⩽ \frac{π}{2 (N + 1)!} |t| min \{1, {(\frac{t^{2} ‖C‖}{4})}^{N + 1}\} e^{t^{2} ‖C‖ / 4} . \end{align}

□

D. Proof of Theorem 3 (solutions of inhomogeneous wave equations)

We give a proof of Theorem 3.

Proof.

In a first step, we note for λ > 0 that

\begin{align} \frac{sin [(t - τ) \sqrt{λ}]}{\sqrt{λ}} = \frac{sin (t \sqrt{λ})}{\sqrt{λ}} cos (τ \sqrt{λ}) - cos (t \sqrt{λ}) \frac{sin (τ \sqrt{λ})}{\sqrt{λ}} . \end{align}

Then,

\begin{align} \frac{sin ((t - τ) \sqrt{A})}{\sqrt{A}} b (τ) = \{\frac{sin (t \sqrt{A})}{\sqrt{A}}) cos (τ \sqrt{A}) - (cos (t \sqrt{A}) \frac{sin (τ \sqrt{A})}{\sqrt{A}}\} b (τ) \\ = β (t) α (τ) b (τ) - α (t) β (τ) b (τ) \end{align}

for all t, τ ∈ ℝ, where α, β : ℝ → L(X, X) are defined by

\begin{align} α (t) ≔ cos (t \sqrt{A}), β (t) ≔ \frac{sin (t \sqrt{A})}{\sqrt{A}}, \end{align}

for every t ∈ ℝ. In the following, for ξ ∈ D(A), we are going to use that the maps

\begin{align} (R \to X, t \mapsto α (t) ξ) and (R \to X, t \mapsto β (t) ξ) \end{align}

are differentiable with derivatives

\begin{align} (R \to X, t \mapsto - β (t) A ξ) and (R \to X, t \mapsto α (t) ξ), \end{align}

respectively. We note that, as a consequence of the spectral theorem for densely defined, self-adjoint linear operators in Hilbert spaces, that a, b are strongly continuous and that

α (t) D (A) \subset D (A), β (t) D (A) \subset D (A),

for every t ∈ ℝ. Also for every k ∈ U_ℂ(σ(A))^s, k(A) D(A) ⊂ D(A), and for ξ ∈ D(A),

{‖k (A) ξ‖}_{A}^{2} = {‖k (A) ξ‖}^{2} + {‖A k (A) ξ‖}^{2} = {‖k (A) ξ‖}^{2} + {‖k (A) A ξ‖}^{2} [⩽ {‖k (A)‖}_{Op}^{2} \cdot {‖ξ‖}_{A}^{2}] .

Hence a, b induce strongly continuous maps from ℝ to X_A, which we indicate with the same symbols, and where

X_{A} ≔ (D (A), {‖‖}_{A})

⁠. In addition, we note that the inclusion ι of X_A into X is continuous. In the next step, we observe for a strongly continuous c : ℝ → L(X_A, X_A) and a continuous g : ℝ → X_A that

\begin{align} {‖c (t + h) g (t + h) - c (t) g (t)‖}_{A} & = {‖c (t + h) g (t + h) - c (t + h) g (t) + c (t + h) g (t) - c (t) g (t)‖}_{A} \\ = {‖c (t + h) {[g (t + h) - g (t)]}_{A} + [c (t + h) - c (t)] g (t)‖}_{A} \\ ⩽ ‖c (t + h)‖ \cdot {‖g (t + h) - g (t)‖}_{A} + {‖c (t + h) g (t) - c (t) g (t)‖}_{A} \end{align}

and hence that (ℝ → X_A, t ↦ c(t) g(t)) is continuous as well as that

\begin{align} (R \to X_{A}, t \mapsto \int_{I_{t}}^{A} c (τ) g (τ) d τ), \end{align}

where ∫^A denotes weak integration in X_A, is differentiable with derivative

(R \to X_{A}, t \mapsto c (t) g (t)) .

We conclude for every t ∈ ℝ that

\begin{align} β (t) \int_{I_{t}}^{A} α (τ) b (τ) d τ - α (t) \int_{I_{t}}^{A} β (τ) b (τ) d τ \\ = \int_{I_{t}}^{A} [β (t) α (τ) b (τ) - α (t) β (τ) b (τ)] d τ \\ = \int_{I_{t}}^{A} \frac{sin ((t - τ) \sqrt{A})}{\sqrt{A}} b (τ) d τ = v (t) . \end{align}

Furthermore, we observe for c : ℝ → L(X, X), g : ℝ → X such that Ran(g) ⊂ D(A), t ∈ ℝ, and h ∈ ℝ^∗ that

\begin{align} \frac{1}{h} [c (t + h) g (t + h) - c (t) g (t)] \\ = \frac{1}{h} [c (t + h) g (t + h) - c (t + h) g (t) + c (t + h) g (t) - c (t) g (t)] \\ = c (t + h) \frac{1}{h} [g (t + h) - g (t)] + \frac{1}{h} [c (t + h) - c (t)] g (t) \\ = c (t) \frac{1}{h} [g (t + h) - g (t)] + \frac{1}{h} [c (t + h) g (t) - c (t) g (t)] \\ + [c (t + h) - c (t)] \frac{1}{h} [g (t + h) - g (t)] \end{align}

and hence that

\begin{align} \frac{1}{h} [α (t + h) g (t + h) - α (t) g (t)] - α (t) g^{'} (t) + β (t) A g (t) \\ = α (t) \{\frac{1}{h} [g (t + h) - g (t)] - g^{'} (t)\} + \frac{1}{h} [α (t + h) g (t) - α (t) g (t)] + β (t) A g (t) \\ + [α (t + h) - α (t)] \{\frac{1}{h} [g (t + h) - g (t)] - g^{'} (t)\} + [α (t + h) - α (t)] g^{'} (t), \\ \frac{1}{h} [β (t + h) g (t + h) - β (t) g (t)] - β (t) g^{'} (t) - α (t) g (t) \\ = β (t) \{\frac{1}{h} [g (t + h) - g (t)] - g^{'} (t)\} + \frac{1}{h} [β (t + h) g (t) - β (t) g (t)] - α (t) g (t) \\ + [β (t + h) - β (t)] \{\frac{1}{h} [g (t + h) - g (t)] - g^{'} (t)\} + [β (t + h) - β (t)] g^{'} (t) . \end{align}

This implies that

(R \to X, t \mapsto α (t) g (t)), (R \to X, t \mapsto β (t) g (t))

are differentiable with derivatives

(R \to X, t \mapsto α (t) g^{'} (t) - β (t) A g (t)), (R \to X, t \mapsto β (t) g^{'} (t) + α (t) g (t)),

respectively. Application of the latter to v gives for t ∈ ℝ,

\begin{align} v^{'} (t) & = β (t) α (t) b (t) + α (t) \int_{I_{t}}^{A} α (τ) b (τ) d τ - α (t) β (t) b (t) + β (t) A \int_{I_{t}}^{A} β (τ) b (τ) d τ \\ = α (t) \int_{I_{t}}^{A} α (τ) b (τ) d τ + β (t) \int_{I_{t}} β (τ) A b (τ) d τ \\ = α (t) \int_{I_{t}}^{A} α (τ) b (τ) d τ + β (t) \int_{I_{t}}^{A} β (τ) A b (τ) d τ, \end{align}

where ∫ denotes weak integration in X, and that

\begin{align} v^{″} (t) & = α (t) α (t) b (t) - β (t) A \int_{I_{t}}^{A} α (τ) b (τ) d τ + β (t) β (t) A b (t) + α (t) \int_{I_{t}} β (τ) A b (τ) d τ \\ = α (t) α (t) b (t) + β (t) β (t) A b (t) - β (t) A \int_{I_{t}}^{A} α (τ) b (τ) d τ + α (t) A \int_{I_{t}}^{A} β (τ) b (τ) d τ \\ = b (t) - A v (t) . \end{align}

□

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 A_C 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

1.

H. G.

Aksoy

and

E.

Senocak

, “

Discontinuous Galerkin method based on peridynamic theory for linear elasticity

,”

Int. J. Numer. Methods Eng.

88

,

673

–

692

(

2011

).

https://doi.org/10.1002/nme.3196

Google Scholar

Crossref

2.

B.

Aksoylu

,

H. R.

Beyer

, and

F.

Celiker

, “

Application and implementation of incorporating local boundary conditions into nonlocal problems

” (unpublished).

3.

B.

Aksoylu

,

H. R.

Beyer

, and

F.

Celiker

, “

Theoretical foundations of incorporating local boundary conditions into nonlocal problems

” (unpublished).

4.

B.

Aksoylu

and

F.

Celiker

, “

Comparison of nonlocal operators utilizing perturbation analysis

,” in

Contributions from European Conference on Numerical Mathematics and Advanced Applications ENUMATH 2015

,

Lecture Notes in Computational Science and Engineering

, edited by

B.

Karasözen

, et al. (

Springer

,

2016

).

Google Scholar

5.

B.

Aksoylu

and

T.

Mengesha

, “

Results on nonlocal boundary value problems

,”

Numer. Funct. Anal. Optim.

31

,

1301

–

1317

(

2010

).

https://doi.org/10.1080/01630563.2010.519136

Google Scholar

Crossref

6.

B.

Aksoylu

and

M. L.

Parks

, “

Variational theory and domain decomposition for nonlocal problems

,”

Appl. Math. Comput.

217

,

6498

–

6515

(

2011

).

https://doi.org/10.1016/j.amc.2011.01.027

Google Scholar

Crossref

7.

B.

Aksoylu

and

Z.

Unlu

, “

Conditioning analysis of nonlocal integral operators in fractional Sobolev spaces

,”

SIAM J. Numer. Anal.

52

,

653

–

677

(

2014

).

https://doi.org/10.1137/13092407X

Google Scholar

Crossref

8.

B.

Alali

and

R.

Lipton

, “

Multiscale dynamics of heterogeneous media in the peridynamic formulation

,”

J. Elasticity

106

,

71

–

103

(

2012

).

https://doi.org/10.1007/s10659-010-9291-4

Google Scholar

Crossref

9.

F.

Andreu-Vaillo

,

J. M.

Mazon

,

J. D.

Rossi

, and

J.

Toledo-Melero

,

Nonlocal Diffusion Problems

,

Mathematical Surveys and Monographs

Vol.

165

(

American Mathematical Society and Real Socied Matematica

,

Espanola

,

2010

).

Google Scholar

Crossref

10.

E.

Askari

,

F.

Bobaru

,

R. B.

Lehoucq

,

M. L.

Parks

,

S. A.

Silling

, and

O.

Weckner

, [“

Peridynamics for multiscale materials modeling

,”

J. Phys.: Conf. Ser.

125

,

012078

(

2008

)].

https://doi.org/10.1088/1742-6596/125/1/012078

Google Scholar

Crossref

11.

H. R.

Beyer

,

Beyond Partial Differential Equations: A Course on Linear and Quasi-linear Abstract Hyperbolic Evolution Equations

,

Lecture Notes in Mathematics

Vol.

1898

(

Springer

,

Berlin

,

2007

).

Google Scholar

12.

L.

Caffarelli

and

L.

Silvestre

, “

An extension problem related to the fractional Laplacian

,”

Commun. Partial Differ. Equations

32

,

1245

–

1260

(

2007

).

https://doi.org/10.1080/03605300600987306

Google Scholar

Crossref

13.

E.

Celik

,

I.

Guven

, and

E.

Madenci

, “

Simulations of nanowire bend tests for extracting mechanical properties

,”

Theor. Appl. Fract. Mech.

55

,

185

–

191

(

2011

).

https://doi.org/10.1016/j.tafmec.2011.07.002

Google Scholar

Crossref

14.

Q.

Du

,

M.

Gunzburger

,

R. B.

Lehoucq

, and

K.

Zhou

, “

Analysis and approximation of nonlocal diffusion problems with volume constraints

,”

SIAM Rev.

54

,

667

–

696

(

2012

).

https://doi.org/10.1137/110833294

Google Scholar

Crossref

15.

Q.

Du

,

L.

Ju

,

L.

Tian

, and

K.

Zhou

, “

A posteriori error analysis of finite element method for linear nonlocal diffusion and peridynamic models

,”

Math. Comput.

82

,

1889

–

1922

(

2013

).

https://doi.org/10.1090/S0025-5718-2013-02708-1

Google Scholar

Crossref

16.

Q.

Du

and

R.

Lipton

,

Peridynamics, Fracture, and Nonlocal Continuum Models

(

SIAM News

,

2014

), Vol.

47

.

Google Scholar

17.

Q.

Du

,

L.

Tian

, and

X.

Zhao

, “

A convergent adaptive finite element algorithm for nonlocal diffusion and peridynamic models

,”

SIAM J. Numer. Anal.

51

,

1211

–

1234

(

2013

).

https://doi.org/10.1137/120871638

Google Scholar

Crossref

18.

Q.

Du

and

K.

Zhou

, “

Mathematical analysis for the peridynamic nonlocal continuum theory

,”

ESAIM: Math. Modell. Numer. Anal.

45

,

217

–

234

(

2011

).

https://doi.org/10.1051/m2an/2010040

Google Scholar

Crossref

19.

N.

Duruk

,

H. A.

Erbay

, and

A.

Erkip

, “

Global existence and blow-up for a class of nonlocal nonlinear Cauchy problems arising in elasticity

,”

Nonlinearity

23

,

107

–

118

(

2010

).

https://doi.org/10.1088/0951-7715/23/1/006

Google Scholar

Crossref

20.

N.

Duruk

,

H. A.

Erbay

, and

A.

Erkip

, “

Blow-up and global existence for a general class of nonlocal nonlinear coupled wave equations

,”

J. Differ. Equations

250

,

1448

–

1459

(

2011

).

https://doi.org/10.1016/j.jde.2010.09.002

Google Scholar

Crossref

21.

E.

Emmrich

,

R. B.

Lehoucq

, and

D.

Puhst

, “

Peridynamics: A nonlocal continuum theory

,” in

Meshfree Methods for Partial Differential Equations VI

, edited by

M.

Griebel

and

M. A.

Schweitzer

(

Springer

,

2013

), Vol.

89

, pp.

45

–

65

.

Google Scholar

Crossref

22.

E.

Emmrich

and

O.

Weckner

, “

The peridynamic equation of motion in non-local elasticity theory

,” in

III European Conference on Computational Mechanics Solids, Structures and Coupled Problems in Engineering, edited by C. M. Soares et al.

,

Lisbon, Portugal

,

June, 1996

.

Crossref

23.

E.

Emmrich

and

O.

Weckner

, “

Analysis and numerical approximation of integro-differantial equation modeling non-local effects in linear elasticity

,”

Math. Mech. Solids

12

,

363

–

384

(

2007

).

https://doi.org/10.1177/1081286505059748

Google Scholar

Crossref

24.

E.

Emmrich

and

O.

Weckner

, “

On the well-posedness of the linear peridynamic model and its convergence towards the Navier equation of linear elasticity

,”

Commun. Math. Sci.

5

,

851

–

864

(

2007

).

https://doi.org/10.4310/CMS.2007.v5.n4.a6

Google Scholar

Crossref

25.

E.

Emmrich

and

O.

Weckner

, “

The peridynamic equation and its spatial discretization

,”

Math. Model. Anal.

12

,

17

–

27

(

2007

).

https://doi.org/10.3846/1392-6292.2007.12.17-27

Google Scholar

Crossref

26.

H. A.

Erbay

,

A.

Erkip

, and

G. M.

Muslu

, “

The Cauchy problem for a one-dimensional nonlinear elastic peridynamic model

,”

J. Differ. Equations

252

,

4392

–

4409

(

2012

).

https://doi.org/10.1016/j.jde.2012.01.008

Google Scholar

Crossref

27.

G.

Gilboa

and

S.

Osher

, “

Nonlocal operators with applications to image processing

,”

Multiscale Model. Simul.

7

,

1005

–

1028

(

2008

).

https://doi.org/10.1137/070698592

Google Scholar

Crossref

28.

B.

Hinds

and

P.

Radu

, “

Dirichlet’s principle and wellposedness of steady state solutions for a nonlocal peridynamics model

,”

Appl. Math. Comput.

219

,

1411

–

1419

(

2012

).

https://doi.org/10.1016/j.amc.2012.07.045

Google Scholar

Crossref

29.

V.

Hutson

,

J. S.

Pym

, and

M. J.

Cloud

,

Applications of Functional Analysis and Operator Theory

, 2 ed. (

Elsevier

,

2005

).

Google Scholar

30.

K.

Joergens

,

Spectral Theory of Second-order Ordinary Differential Operators

(

University of Aarhus

,

1962

).

Google Scholar

31.

J. F.

Kalthoff

and

S.

Winkler

, “

Failure mode transition at high rates of shear loading

,” in

Impact Loading and Dynamic Behavior of Materials

, edited by

C.

Chiem

,

H.-D.

Kunze

, and

L.

Meyer

(

DGM Informationsgesellschaft Verlag

,

1988

), Vol. 1, pp. 185–195.

Google Scholar

32.

B.

Kilic

, “

Peridynamic theory for progressive failure prediction in homogeneous and heterogeneous materials

,” Ph.D. thesis,

Department of Aerospace and Mechanical Engineering, University of Arizona

, Tucson, AZ,

2008

.

Google Scholar

33.

B.

Kilic

,

A.

Agwai

, and

E.

Madenci

, “

Peridynamic theory for progressive damage prediction in centre-cracked composite laminates

,”

Compos. Struct.

90

,

141

–

151

(

2009

).

https://doi.org/10.1016/j.compstruct.2009.02.015

Google Scholar

Crossref

34.

B.

Kilic

and

E.

Madenci

, “

Prediction of crack paths in a quenched glass plate by using peridynamic theory

,”

Int. J. Fract.

156

,

165

–

177

(

2009

).

https://doi.org/10.1007/s10704-009-9355-2

Google Scholar

Crossref

35.

B.

Kilic

and

E.

Madenci

, “

Coupling of peridynamic theory and finite element method

,”

J. Mech. Mater. Struct.

5

,

707

–

733

(

2010

).

https://doi.org/10.2140/jomms.2010.5.707

Google Scholar

Crossref

36.

S.

Kindermann

,

S.

Osher

, and

P. W.

Jones

, “

Deblurring and denoising of images by nonlocal functionals

,”

Multiscale Model. Simul.

4

,

1091

–

1115

(

2005

).

https://doi.org/10.1137/050622249

Google Scholar

Crossref

37.

R.

Lipton

, “

Dynamic brittle fracture as a small horizon limit of peridynamics

,”

J. Elasticity

117

,

21

(

2014

).

https://doi.org/10.1007/s10659-013-9463-0

Google Scholar

Crossref

38.

E.

Madenci

and

E.

Oterkus

,

Peridynamic Theory and Its Applications

(

Springer

,

2014

).

Google Scholar

Crossref

39.

T.

Mengesha

and

Q.

Du

, “

The bond-based peridynamic system with Dirichlet-type volume constraint

,”

Proc. R. Soc. Edinburgh, Sect. A

144A

,

161

–

186

(

2014

).

https://doi.org/10.1017/S0308210512001436

Google Scholar

Crossref

40.

T.

Mengesha

and

Q.

Du

, “

Nonlocal constrained value problems for a linear peridynamic Navier equation

,”

J. Elasticity

116

,

27

–

51

(

2014

).

https://doi.org/10.1007/s10659-013-9456-z

Google Scholar

Crossref

41.

Y.

Mikata

, “

Analytical solutions of peristatic and peridynamics problems for a 1D infinite rod

,”

Int. J. Solids Struct.

49

,

2887

–

2897

(

2012

).

https://doi.org/10.1016/j.ijsolstr.2012.02.012

Google Scholar

Crossref

42.

S. G.

Mikhlin

,

Mathematical Physics, An Advanced Course

(

North-Holland

,

Amsterdam

,

1970

).

Google Scholar

43.

NIST Handbook of Mathematical Functions

, edited by

F. W. J.

Olver

,

D. W.

Lozier

,

R. F.

Boisvert

, and

C. W.

Clark

(

Cambridge University Press

,

New York

,

2010

).

Google Scholar

44.

E.

Oterkus

and

E.

Madenci

, “

Peridynamic analysis of fiber reinforced composite materials

,”

J. Mech. Mater. Struct.

7

,

45

–

84

(

2012

).

https://doi.org/10.2140/jomms.2012.7.45

Google Scholar

Crossref

45.

E.

Oterkus

and

E.

Madenci

, “

Peridynamic theory for damage initiation and growth in composite laminate

,”

Key Eng. Mater.

488–489

,

355

–

358

(

2012

).

https://doi.org/10.4028/www.scientific.net/kem.488-489.355

Google Scholar

Crossref

46.

E.

Oterkus

,

E.

Madenci

,

O.

Weckner

,

S. S.

Silling

,

P.

Bogert

, and

A.

Tessler

, “

Combined finite element and peridynamic analyses for predicting failure in a stiffened composite curved panel with a central slot

,”

Compos. Struct.

94

,

839

–

850

(

2012

).

https://doi.org/10.1016/j.compstruct.2011.07.019

Google Scholar

Crossref

47.

A.

Pazy

,

Semigroups of Linear Operators and Applications to Partial Differential Equations

(

Springer

,

New York

,

1983

).

Google Scholar

Crossref

48.

M.

Reed

and

B.

Simon

,

Methods of Modern Mathematical Physics

(

Academic

,

New York

,

1980

), Vol. IV; Vol. I, 1975; Vol. II, 1978; Vol. III, 1979.

Google Scholar

49.

W.

Rudin

,

Functional Analysis

(

McGraw-Hill

,

New York

,

1991

).

Google Scholar

50.

P.

Seleson

,

S.

Beneddine

, and

S.

Prudhomme

, “

A force-based coupling scheme for peridynamics and classical elasticity

,”

Comput. Mater. Sci.

66

,

34

–

49

(

2013

).

https://doi.org/10.1016/j.commatsci.2012.05.016

Google Scholar

Crossref

51.

P.

Seleson

,

M.

Gunzburger

, and

M. L.

Parks

,

Interface problems in nonlocal diffusion and sharp transitions between local and nonlocal domains

Comput. Methods Appl. Mech. Eng.

266

,

185

–

204

(

2013

).

https://doi.org/10.1016/j.cma.2013.05.018

Google Scholar

Crossref

52.

P.

Seleson

,

M.

Gunzburger

, and

M. L.

Parks

, “

Peridynamic state-based models and the embedded-atom model

,”

Commun. Comput. Phys.

15

,

179

–

205

(

2014

).

https://doi.org/10.4208/cicp.081211.300413a

Google Scholar

Crossref

53.

P.

Seleson

and

M. L.

Parks

, “

On the role of the influence function in the peridynamic theory

,”

Int. J. Multiscale Comput. Eng.

9

,

689

–

706

(

2011

).

https://doi.org/10.1615/IntJMultCompEng.2011002527

Google Scholar

Crossref

54.

P.

Seleson

,

M. L.

Parks

,

M.

Gunzburger

, and

R. B.

Lehoucq

, “

Peridynamics as an upscaling of molecular dynamics

,”

Multiscale Model. Simul.

8

,

204

–

227

(

2009

).

https://doi.org/10.1137/09074807X

Google Scholar

Crossref

55.

S.

Silling

, “

Reformulation of elasticity theory for discontinuities and long-range forces

,”

J. Mech. Phys. Solids

48

,

175

–

209

(

2000

).

https://doi.org/10.1016/S0022-5096(99)00029-0

Google Scholar

Crossref

56.

S.

Silling

and

R.

Lehoucq

, “

Convergence of peridynamics to classical elasticity

,”

J. Elasticity

93

,

13

–

37

(

2008

).

https://doi.org/10.1007/s10659-008-9163-3

Google Scholar

Crossref

57.

S.

Silling

and

R. B.

Lehoucq

, “

Peridynamic theory of solid mechanics

,”

Adv. Appl. Mech.

44

,

73

–

168

(

2010

).

https://doi.org/10.1016/S0065-2156(10)44002-8

Google Scholar

Crossref

58.

S. A.

Silling

, “

Dynamic fracture modeling with a meshfree peridynamic code

,” in

Computational Fluid and Solid Mechanics

, edited by

K.

Bathe

(

Elsevier

,

2003

), pp.

641

–

644

.

Google Scholar

Crossref

59.

S. A.

Silling

and

F.

Bobaru

, “

Peridynamic 3D problems of nanofiber networks and carbon nanotube-reinforced composites

,”

AIP Conf. Proc.

712

,

1565

–

1570

(

2004

).

https://doi.org/10.1063/1.176675

Google Scholar

60.

S. A.

Silling

and

F.

Bobaru

, “

Peridynamic modeling of membranes and fibers

,”

Int. J. Nonlinear Mech.

40

,

395

–

409

(

2005

).

https://doi.org/10.1016/j.ijnonlinmec.2004.08.004

Google Scholar

Crossref

61.

S. A.

Silling

,

M.

Zimmermann

, and

R.

Abeyaratne

, “

Deformation of a peridynamic bar

,”

J. Elasticity

73

,

173

–

190

(

2003

).

https://doi.org/10.1023/B:ELAS.0000029931.03844.4f

Google Scholar

Crossref

62.

X.

Tian

and

Q.

Du

, “

Analysis and comparison of different approximations to nonlocal diffusion and linear peridynamic equations

,”

SIAM J. Numer. Anal.

51

,

3458

–

3482

(

2013

).

https://doi.org/10.1137/13091631X

Google Scholar

Crossref

63.

O.

Weckner

and

R.

Abeyaratne

, “

The effect of long-range forces on the dynamics of a bar

,”

J. Mech. Phys. Solids

53

,

705

–

728

(

2005

).

https://doi.org/10.1016/j.jmps.2004.08.006

Google Scholar

Crossref

64.

J.

Weidmann

,

Linear Operators in Hilbert Spaces

,

Graduate Texts in Mathematics

Vol.

68

(

Springer-Verlag

,

1980

).

Google Scholar

Crossref

65.

J.

Weidmann

, “

Strong operator convergence and spectral theory of ordinary differential operators

,”

Univ. Iagellon. Acta Math.

34

,

153

–

163

(

1997

).

Google Scholar

66.

K.

Yosida

,

Functional Analysis

(

Springer

,

Berlin

,

1980

).

Google Scholar

67.

K.

Zhou

and

Q.

Du

, “

Mathematical and numerical analysis of linear peridynamic models with nonlocal boundary conditions

,”

SIAM J. Numer. Anal.

48

,

1759

–

1780

(

2010

).

https://doi.org/10.1137/090781267

Google Scholar

Crossref

68.

M.

Zimmermann

, “

A continuum theory with long-range forces for solids

,” Ph.D. thesis,

Massachusetts Institute of Technology, Department of Mechanical Engineering

, Cambridge, MA,

2005

.

Google Scholar

69.

Throughout the paper, we refer to the Laplace operator as the local or the classical operator interchangeably.

70.

Solution operators, $cos (t \sqrt{A})$ and $sin (t \sqrt{A}) / \sqrt{A},$ denote operators corresponding to unique extensions to entire holomorphic functions. A pedantic notation would be $[\bar{cos (t \sqrt{})} |_{σ (A)}] (A)$ and $[\bar{\frac{sin (t \sqrt{})}{\sqrt{}}} |_{σ (A)}] (A)$ ⁠, respectively. For brevity, we prefer to use the former.

71.

$U_{C}^{s} (σ (A), C)$ 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.

2016

Author(s)

On a class of nonlocal wave equations from applications

I. INTRODUCTION AND MOTIVATION

II. OPERATOR-THEORETIC TREATMENT OF SYSTEMS OF WAVE EQUATIONS

III. THE GOVERNING OPERATOR AND PROPERTIES

IV. MAIN RESULTS

V. EXAMPLES AND ILLUSTRATIONS OF MAIN RESULT 3

VI. PROOFS AND RELATED RESULTS

A. Proof of main result 1 (nonlocal operator is a bounded function of the classical operator)

B. Proof of main result 2 (strong convergence of nonlocal solutions to classical ones through strong resolvent convergence)

C. Proof of main result 3 (representation of the solution in terms of spherical Bessel functions)

D. Proof of Theorem 3 (solutions of inhomogeneous wave equations)

VII. CONCLUSION

Acknowledgments

REFERENCES

Citing articles via

Submit your article

Sign up for alerts

Resources

Explore

pubs.aip.org

Connect with AIP Publishing

On a class of nonlocal wave equations from applications Free

I. INTRODUCTION AND MOTIVATION

II. OPERATOR-THEORETIC TREATMENT OF SYSTEMS OF WAVE EQUATIONS

III. THE GOVERNING OPERATOR AND PROPERTIES

IV. MAIN RESULTS

V. EXAMPLES AND ILLUSTRATIONS OF MAIN RESULT 3

VI. PROOFS AND RELATED RESULTS

A. Proof of main result 1 (nonlocal operator is a bounded function of the classical operator)

B. Proof of main result 2 (strong convergence of nonlocal solutions to classical ones through strong resolvent convergence)

C. Proof of main result 3 (representation of the solution in terms of spherical Bessel functions)

D. Proof of Theorem 3 (solutions of inhomogeneous wave equations)

VII. CONCLUSION

Acknowledgments

REFERENCES

Citing articles via

Submit your article

Sign up for alerts

Related Content

Resources

Explore

pubs.aip.org

Connect with AIP Publishing

This Feature Is Available To Subscribers Only

On a class of nonlocal wave equations from applications