Qualitative Aspects in Nonlocal Dynamics

In this paper, we investigate, through numerical studies, the dynamical evolutions encoded in a linear one-dimensional nonlocal equation arising in peridynamics. The different propagation regimes ranging from the hyperbolic to the dispersive, induced by the nonlocal feature of the equation, are carefully analyzed. The study of an initial value Riemann-like problem suggests the formation of a singularity.


Introduction
In the previous works [3,2] some relevant mathematical features characterizing the dynamical evolution ruled by a peridynamic model introduced by S.A. Silling ([14,15,16,17,18]) have been studied in detail by focusing on the mathematical consistency (see also [6]) of this new formulation of continuum mechanics.The nonlocal effects related to this model represent a very interesting problem in mathematical physics which has captured the interest of many researchers during the past decades (see e.g.[7,8,11,10,12]).
In the present paper we exploit a qualitative analysis of a one-dimensional linear peridynamical model to highlight the behavior of the solutions of an initial value problem in dependence of the characteristics affecting the nonlocal properties of the equation.More precisely, we show through numerical investigations the way nonlocality rules a wide range of modes of wave propagation (see also [5,13]).Indeed, the dynamics predicted by the equation under study exhibits a behavior ranging from the hyperbolic-like to the dispersive-like propagation ( [21]), according to its highly nontrivial dispersive relation.These new phenomena, unpredictable in classical continuum mechanics ( [9]), are ruled by two constitutive parameters both related to the nonlocal character of the equation.
The paper is organized as follows.In Section 2 we introduce the main initial value problem in terms of a linear peridynamic model and the related dispersion relation focusing on the asymptotic limits.Morever, the important notion of group velocity ruling the transport of energy is analyzed.In Section 3 we study the evolution of an initial localized profile in the range of the main parameters characterizing the problem.The shape of the initial datum is contrasted with the parameters affecting nonlocality, thus revealing an entire catalogue of different interesting behaviors.In Section 4 we show how a suitable choice of parameters delivers the case of purely hyperbolic propagation exhibiting traveling wave solution.Eventually, in Section 5 we explore the possible occurrence of singularities ( [4]), allowed in principle for energy finite solutions.To this end, the study of a Riemann-like problem allows room for the formation of a singularity.

One-dimensional peridynamic model
We perform a qualitative study regarding the phenomenology of the linear peridynamic model for α ∈ (0, 1).As shown in [2], this problem admits the unique formal solution where v 0 (ξ) and v 1 (ξ) are respectively the Fourier transforms of the initial conditions v 0 (x) and v 1 (x) and , is the dispersive relation actually governing the dynamics.Eq. (2.3) admits solution for a wide class of initial conditions (see [2] for the rigorous discussion in this regard).However, in this manuscript, we limit our phenomenological investigation to the particular case of initial conditions formally corresponding to a (possibly approximate) traveling wave.This case is described with any loss of generality by the constraint where v is a number denoting the initial velocity of the perturbations (or, in frequency space, one can recast Eq. (2.4) in the form v 1 (ξ) = ivξ v 0 (ξ)).Eq. (2.4) can be easily understood as the limit for t → 0 of the velocity of a traveling waveform f (x − vt): indeed, its velocity is just ∂ t f (x − vt) = −v ∂ q f (q), where q = x − vt.Hence, in the above-mentioned limit, q → x and the constraint (2.4) follows.
The physical interest in the traveling wave initial conditions stands in the comparison with the linear 1-D wave equation.Indeed, since the latter admits non-dispersive propagations, any traveling wave initial conditions will evolve as at most two independent traveling waves, according to the value of v. On the opposite, Eq. (2.3) exhibits a highly non-trivial relation between ω and ξ.More precisely, it has been proven in [2] that Eq. (2.3) admits two characteristic behaviors discriminated by the length scale δ, namely where Γ is the Euler gamma function.Hence, for large scale modes, the dynamics is practically non-dispersive and these modes travel with group velocity v g ≡ ω (ξ) ≈ δ 1−α κ ρ(1−α) .On the other hand, small scale modes experience a sub-linear dispersive dynamics with group velocity v g ≈ ξ α−1 4κ ρ ´∞ 0 1−cos τ τ 1+2α dτ .According to the initial conditions, qualitative differences with respect to the classic case are expected when small scales are switched on along the dynamics.In the following section, we will explore several scenarios to point out this peculiarity.
Let us remark that α ∈ (0, 1) plays a crucial role in the subsequent analysis.Indeed, besides ruling the smoothness of the admissible solutions for the problem (2.1), it is responsible in tuning the effects of the nonlocal feature in the present theory.More precisely, the higher the α, the more hyperbolic the propagation, whereas the lower the α, the stronger the nonlocal effects, leading to dispersive propagation.
We conclude this section by proving for the Cauchy problem in (2.1) that the group velocity given by v g := ω (ξ) represents (see [1]) the velocity of propagation of the energy density of monochromatic waves.To this end, we define the energy density e e := ρ and the energy flux G G := ω∂ ξ W , where W is the potential energy density In this setting, we consider a monochromatic wave Here, for simplicity, up to a phase shift, we assume that A ∈ R. At a formal level, this wave would correspond to picking a ξ ∈ R and choosing as initial conditions v 0 such that v 0 is the Dirac delta at ξ and v 1 such that (2.4) is satisfied.We claim that (2.7) To check this, we observe that the potential energy density of the monochromatic wave is given by On the other hand, the total energy density is given by hence the relation in (2.7) is proven.This result is interesting since it extends the classical result for the wave equation [1,21] to the linear peridynamic model.

Localized initial conditions
In this section we study the dispersive phenomena of the linear peridynamic model in the range of parameters of α and δ.As a first step, we fix our constitutive parameters as ρ = 1 and κ = 1/2 and leave them unchanged for the rest of the discussion.Aiming to study the traveling waves, we then consider Gaussian initial conditions The interesting feature of those initial conditions relies in their localization close to the origin.Indeed, we will provide numerical evidences that these localized initial conditions remain localized as time evolves.This suggests that the evolution experiences almost finite propagation velocity.From now on, the key point to have in mind will be always the relative amplitude between the two parameters δ and group velocity v g characterizing the nonlocality of (2.1) and the two parameters σ and v involved in the initial conditions.The rest of the section is divided into three parts, devoted to the analysis of the interplay between the nonlocal range δ and the amplitude of the initial condition σ, since their ratio plays a crucial role in occurrence of dispersive effects in the wave propagation.More precisely, when σ is much smaller than δ, the solution develops high oscillations typical of the dispersive regime almost immediately, showing the deep difference between the peridynamics and the wave equation.When σ is of the same order of δ we capture the strong transport effect in short times and the occurrence of dispersion on the long range.Finally, when σ is much bigger than δ, the dynamics experiences the (almost) total absence of dispersion and the occurrence only of hyperbolic propagation, typical of the wave equation.In order to support our claim about the immediate occurrence of dispersive behavior, independently of the choice of α, we show different cases regarding three possible choices of α.
In Fig. 1, we consider α = 10 −1 .In Fig. 2 we consider α = 1/2.In Fig. 3 we consider α = 9/10.In order to support our claim about the competitive occurrence of dispersive behavior and hyperbolic propagation, according of the value of α, we show different cases regarding three possible choices of α.
In Fig. 4, we consider α = 10 −1 .In Fig. 5 we consider α = 1/2.In Fig. 6 we consider α = 9/10.In order to support our claim about the pure hyperbolic propagation and (almost) total absence of dispersive behavior, with respect to the value of α, we show different cases regarding three possible choices of α.
In Fig. 7 we consider α = 10 −1 .In Fig. 8 we consider α = 1/2.In Fig. 9 we consider α = 9/10.From the previous numerical analysis, we emphasize an important feature affecting the propagation problem under study.Indeed, all the previous plots clearly put in evidence the existence of a propagation velocity v p bounded by Heuristically, one can appreciate that for every t, the solutions live almost entirely within the region |x| ≤ v max t.The physical reason behind the value of v max stands in the fact that this is the upper bound of the group velocity v g , namely no frequency can travel at a propagation velocity higher than v max itself.In the next section, we will investigate in detail the finite propagation velocity ruling the purely hyperbolic case.

Hyperbolic-like propagation
This section is devoted to the special case when the solution is a (possibly approximated) traveling wave.For this case we choose In particular, we study the case σ = 10 and δ = 1.In Fig. 10 we consider α = 10 −1 .In Fig. 11 we consider α = 1/2.In Fig. 12 we consider α = 9/10.

Singularity formation
In this section we explore the possible spontaneous creation of singularities.Following the analysis performed in [19,20] one can try to study the following Riemann-like problem augmenting (2.1) with the following initial conditions (5.1) Due to the fact that the Fourier transform of v 1 exhibits a strong singularity at ξ = 0, rather than using the representation formula (2.2) in this case, we confine our analysis to some numerical experiments.For this, we truncate (5.1) and choose v − = 0 and v + = 1, namely which does not exhibit the Dirac delta at ξ = 0.In Figures 13 and 14 we plot the real and imaginary parts of the Fourier transform of the solution expressed in (2.2) with initial data given in (5.2).We can appreciate that as time evolves the real part develops a peak at ξ = 0 whereas the imaginary one becomes highly oscillating.Inspired by the analysis carried out in [20] we observe numerically a possible development of a stationary singularity at x = 0.

Conclusion
In this paper, we have explored the features of the peridynamic model proposed in [3,2] as a viable generalization of the wave equation for the study of nonlocal phenomena in the mechanics of continuous media.We have focused our phenomenological analysis to the 1-D linear scenario.Despite its simplicity, this case study already exhibits deep differences with respect to the standard wave equation.In particular, dispersive propagation can occur when small scales dynamics is non-negligible.
This competitive behavior between classical hyperbolic propagation and dispersive phenomena is intimately related to the multiscale essence of the peridynamic framework and the crossover from one regime to the other one is dictated by the interaction length δ and the interaction scaling α.Moreover, we have also shown that the energy propagates with the group velocity, regardless of the scale lengths involved in the dynamics.
From these general features, we have then exploited the phenomenology of the case of Gaussian initial condition with non-trivial velocity.In this scenario, our numerical investigation has shown that dispersive phenomena are relevant when the width σ of the Gaussian initial profile is smaller than the peridynamic radius δ.On the contrary, when σ is much greater than δ, our numerical investigation exhibits a hyperbolic-like propagation.This regime can be interpreted as a direct consequence of the multiscale nature of the dispersive relation ω(ξ).This is nearly linear when ξ δ −1 and grows with a sub-linear rate when ξ δ −1 .As a final comparison between the hyperbolic and dispersive behaviors, for the case σ δ, we have shown the existence of a typical initial velocity of the Gaussian initial condition leading to traveling wave-like solutions.Once again, the cornerstone to understand this behavior is the dispersive relation ω(ξ): indeed this fine-tuned choice for the initial velocity is given by the limit for ξ → 0 of the group velocity ω (ξ), which is a constant, as a consequence of the linearity of ω(ξ) on large scale.This limit is also responsible for the almost finite propagation of the solutions over time.
Finally, in the last section, we have discussed the possible spontaneous formation of singularities developed in an initial-value problem with discontinuous velocity datum.

Date:
June 28, 2021.2010 Mathematics Subject Classification.74A70, 74H05, 35Q70, 35R09, 65R20.GMC, SD, FM, MR and EV are members of the Gruppo Nazionale per l'Analisi Matematica, la Probabilità e le loro Applicazioni (GNAMPA) of the Istituto Nazionale di Alta Matematica (INdAM).GF is member of the Gruppo Nazionale per la Fisica Matematica (GNFM) of the Istituto Nazionale di Alta Matematica (INdAM).SD and EV are members of the AustMS.GMC, FM and MR have been partially supported by the Research Project of National Relevance "Multiscale Innovative Materials and Structures" granted by the Italian Ministry of Education, University and Research (MIUR Prin 2017, project code 2017J4EAYB and the Italian Ministry of Education, University and Research under the Programme Department of Excellence Legge 232/2016 (Grant No. CUP -D94I18000260001).GF acknowledges support by FCT under the program Stimulus with the grant no.CEECIND/04399/2017/CP1387/CT0026.EV has been supported by the Australian Laureate Fellowship FL190100081 "Minimal surfaces, free boundaries and partial differential equations".

Figure 13 .Figure 14 .
Figure 13.Real part of the Fourier transform of the solution with initial conditions (5.2) at different times.This case refers to the parameters α = 1/2 and δ = 1.