Abstract
By means of a particle model that includes interactions only via the local particle concentration, we show that hyperballistic diffusion may result. This is done by findng the exact solution of the corresponding non-linear diffusion equation, as well as by particle simulations. The connection between these levels of description is provided by the Fokker-Planck equation describing the particle dynamics. PACS numbers:
I Introduction
Superdiffusion is characterized by the fact that the root mean square displacement of some kind of particles, increases with time t as with the exponent , the normal diffusion value being . This behavior may arise in physical, biological or geological systems; examples include Levy flights [1, 2], particle motion in random potentials or the seemingly random paths of objects moving in turbulent flows [3, 4].
Biological examples may be found in the foraging movement of spider monkeys [5] and the flight paths of albatrosses [6, 7]; in both cases . These movements are Levy walks, which are random walks of uncorrelated steps of length , that take their value from a distribution . They result in superdiffusive behavior with when [1].
However, the mere observation that the step length distribution has a fat tail, does not by itself provide any physical model to explain the superdiffusive behavior. The simplest physical example of superdiffusion is perhaps provided by the undamped Langevin equation which describes a random walk in momentum space and a corresponding real space displacement with [8]. This kind of behavior is termed hyperballistic as . Quantum- or classical particles in random potentials behave much like those described by the undamped Langevin equation, and yield hyperballistic diffusion with [8] too, though Golubovic et al. [9] studied a case where . In optical experiments [10, 11] where the spatial coordinate in the direction of the light plays the role of the time coordinate, hyperballistic spreading has been observed as well. This effect is linked to Anderson localization [12], and comes from a transition where the light modifies its mean free path as it passes through the medium.
Anomalous diffusion of the subdiffusive kind has been studied in a wide range of contexts: It may be observed in compressible gases flowing through porous media [13, 14], a pulse of energy propagating in vacuum [15], or in filtration processes [16]. Another example is heat diffusion at high temperature [17, 18]. Population dynamics gives rise to this kind of behavior [19–21], as does water ingress in zeolites as studied by Azevedo et al. [22, 23] and Fischer et al. [24]. The diffusion of grains in granular media considered by Christov and Stone [25] is yet another example. Pritchard et al., [26] studied gravity-driven fluid flow in layered porous media finding that the fluid motion could be described by a concentration-dependent diffusivity as did Hansen et al. [27] for the spreading of wetting films in wedges. Anomalous diffusion in random geometries, fractals and tree-like structures has been studied for decades [28–32]. Common to all of these examples is subdiffusion,
Hyperballistic diffusion seems almost a contradiction in terms, for how could a random walker move faster than a directed walker that never changes direction? The explanation lies in the fact that the velocity, and thus the step length, keeps increasing with time without limits. This behavior is of course unphysical in the context of the Langevin equation as there will always be dissipative forces that match the fluctuations, but has a physical basis in random potentials. On the other hand, in a hydrodynamic shear-flow that increases without bounds, a random walker will achieve step-lengths that are umlimited too [33, 34], an effect that may give rise to hyper-ballistic diffusion. Without diverging velocities or step lengths, long range time-correlations are required for superdiffusion, an example being the elephant random walk, so named because both the walkers and elephants have long memories, which in the model give rise to (sub-ballistic) superdiffusion [35].
Generally, superdiffusion has been modeled by independent agents interacting with an environment, or possessing a long term memory [36]. The main question of the present article is if superdiffusion, including the hyperballistic case, could result directly from a Markovian description of particle interactions. Such interactive systems could include crowds of people, bacteria swimmers competing for food [37, 38] or the evolution of the porosity in a granular packing. For the purpose of addressing this question we investigate the potentially simplest description of particle interactions, namely, that where a conserved concentration C of particles is governed by Ficks law . Here the C-dependence in D reflects interactions between the particles; in many cases of interest these interactions are well captured by this type of mean field description.
II Solution to the Non-Linear Diffusion Equation
Already in 1959 did Pattle [39] solve the diffusion equationwhere is the concentration and D is given by the power law where is a constant reference concentration, is the diffusivity at that reference value, and the exponent . Pattle found the root mean square displacement withwhere d is the dimension. For negative γ this will always lead to sub-diffusion. We have recently shown that in there are exact solutions with positive γ as well [40], which still satisfy Eq. 2, thus yielding superdiffusion with as always. In the present article we take this result further by deriving the solution for and for in any dimension. When the corresponding exponent t will then take on any value, including those of the hyperballistic regime, implying that hyperballistic diffusion is a higher-dimensional effect. We coin the term “explosive” for the corresponding time dependence of because the decay of an intially localized C-profile is qualitatively faster than normal diffusive, or even superdiffusive, decay.
To validate the mean field description and provide it with a physical basis, we introduce a particle model that is described by Eq. 1. The step lengths in this model , and therefore correspond to velocities that diverge as . This would correspond to an unlimited access to thermal energy. However, unlike the Langevin equation where [8], this model can produce any τ-value.
Following the same lines as in [40] we rewrite Eq. 1 asHence, we see that we need for the equation to be defined when . The initial condition at is a point source pulse containing particles, . This means that there is no intrinsic length- or time scale in the problem, and the particle number inside a radius r should satisfy the scale invariance condition for some . Differentiating this equation with respect to r, using the fact that leads to the scaling relationWe are free to chose λ such that , that is by requiring that where for simplicity, we have introduced the function , andwhere we have introduced and the reduced variable . Inserting Eq. 5 in Eq. 3 yieldsfor some dimensionless constant c, which can be absorbed in the definition of . The point-like initial condition, implies , and the left hand side of Eq. 6 can be easily integrated to giveNote That This Form Immediately Giveswith τ given by Eq. 2.
From Eq. 6, we also have an expression for ,which can be integrated to give,
For Fick’s law to be valid throughout the domain, , and therefore, , must be differentiable everywhere when . To avoid a spike at the origin we must have and also a finite , which implies that . So, Eq. 10 may be integrated to yieldwhere k is an integration constant. This expression is independent of the dimension d. The value of the constant k can be determined through the normalization, , which givesand yields the concentration field by means of Eq. 5; 11. The mean square displacement is given bywhich is limited to the range of γ-values where the integrals in Eq. 8 converge. Since for large r this range is . However, in any particle simulation there will always be a largest particle position that will act as a cut-off. This means that the factor in survives, but that its prefactor will fluctuate with the value. The behavior with different d and γ is summarized in Table 1.
TABLE 1
| rrms Prediction converges | C Is normalizable | ||
|---|---|---|---|
| d = 1 | Never | All | |
| d = 2 | All | ||
| d = 3 |
Behavior with γ in various dimensions d as predicted by Eq. 2.
Interestingly, there exists an alternative route to the solution given in Eq. 5; 11: Working in Plastino and Plastino [41] showed that by adding a drift term to Eq. 1 that corresponds to the force from a harmonic potential, a stationary solution could be found. Using the anzats that the full time-dependent solution has that same structure as the stationary one, only with time-dependent coefficients, the structure of Eq. 5; 11 is established. Later, Tsallis and Bukman [42] established the full analytic solution to this problem, which structurally reduces to Eq. 5; 11 as the strength of the potential is taken to zero.
III Particle Model that Realizes the Non-Linear Diffusion Equation
We will employ two simulation models, both in with random walkers, labeled i, that have positions . The particles interact only via the value of C, which is the local population density. The steps are chosen isotropically at each time step; to find their length we need to derive the appropriate Fokker-Planck equation and match it to Eq. 3. For every time step the walkers movewhere α is a Cartesian index and the function is to be determined. This defines a Wiener process with η as a random variable with and . Now, following the same steps as in [40, 43] we use the standard Chapman-Kolmogorov, or master equation, to derive the following Fokker-Planck equation for the particle concentration Here is the mean squared jump length per time,where is the probability per unit time that a walker jumps a distance from . Setting givesand requiring equivalence with Eq. 3 thus implies that . This leads to the stepwhere the random variable η is given above.This defines the particle model that is described by Eq. 3.
In the finite interaction range model C is calculated by assuming a maximum interaction range between particles. This is done by calculating C onto a lattice with lattice constant : The local value at the discrete site is simply times the number of particles at positions that satisfy . The step length for a particle that is located at depends on the C-value at the nearest lattice site. The finite interaction range of this model has a discretization effect: Once C is so small that there is only one- or zero particles in each cell, the step length will always be the same, and as a result, there will be a cross-over to normal diffusion, an effect that is observed in the curve of Figure 1A.
FIGURE 1
The other, infinite interaction range model employs no lattice at all, but evaluates C at any particle position as where is a fixed particle number and is the volume of the sphere that contains nearest neighbors, as illustrated by Figure 2 in the case with .
FIGURE 2
There is no upper limit to the size of , and it is in this sense that the model has a potentially infinite interaction range. When this model will never cross over to normal diffusive behavior. A somewhat similar particle, but one-dimensional, particle model was introduced by Borland [44], who in stead of calculating the local C-value from the neighboring particles, employed the analytical solution for C.
In Figure 3 the analytic solution of Eq. 11 is plotted for different γ-values. The term “explosive” seems an appropriate label for the behavior of the concentration for two reasons: First, as close to the critical value of 2/3, the initial concentration drops by more than 10 orders of magnitude in the same time that the negative γ solutions (taken from Pattle [39]), drop by less than two orders. Second, the divergence of the integral in Eq. 8 defining signals a cross-over to a regime where the break-away particles dominate the behavior at ever increasing step lengths.
FIGURE 3
Figure 4 show simulations using dimensionless spatial and time coordinates. If units were assigned to them the background diffusivity would have dimension length2/time as usual. The time step is chosen in order to avoid significant changes in the local concentration from time-step to time-step. Note the increasing presence of particles that separate from the main crowd as γ is increased.
FIGURE 4
In Figure 5 the data collapse anticipated in Eq. 5 is seen to be satisfied. Figures 1A,B demonstrate that the particle displacement is in fact characterized by Eq. 13, the difference between Figures 1A,B, being that the first figure compares simulations and the full analytic prediction of Eq. 13, while the hyperballistic transport shown in Figures 1B, only confirms the prediction of the τ exponent, Eq. 2. Note that in Figures 1A the convergence to the prediction of Eq. 13, happens over a time that increases with γ, signaling the end of the regime where has an exact analytical expression.
FIGURE 5
Figure 6 summarizes this comparison for the full range of relevant γ-values, using the finite-range model for the smaller- and the infinite range model for the larger γ-values.
FIGURE 6
FIGURE 7
IV Conclusion
In conclusion, we have shown that particle interactions described entirely in terms of their local concentration may yield superdiffusion, and even hyperballistic diffusion. This was done by solving the diffusion equation with the diffusivity exactly. The particle interactions were described in terms of this concentration dependence alone. Unlike earlier solutions [42, 44] in the present solutions yield hyper-ballistic diffusion. In () this happens when . The 3-dimensional particle model that was introduced as a realization of this diffusion equation was found to reproduce the exact solution for a range of γ-values, and also the predicted root-mean-square displacement in the range of γ values where this prediction is finite. Notably, also outside this range () did the particle simulation confirm the predicted diffusion exponent.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
All authors contributed to the analytic work and discussions. EGF did the simulations and the writing of the paper.
Funding
This work was partly supported by the Research Council of Norway through its Centers of Excellence funding scheme, project number 262644.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1.
BouchaudJGeorgesA. Anomalous diffusion in disordered media: statisticalmechanisms, models, and physical applications. Phys Rep (1990) 195:125. 10.1016/0370-1573(90)90099-n
2.
GoshSKCherstvyAGGrebenkovDMetzlerR. Anomalous non-gaussiantracer diffusion in crowded two-dimensional environments. New J Phys (2016) 18:013027. 10.1088/1367-2630/18/1/013027
3.
RichardsonL. Atmospheric diffusion shown on a distance-neighbour graph. Proc Roy Soc London A (1926) 100:709. 10.1098/rspa.1926.0043
4.
SchlesingerMWestBKlafterJ. Levy dynamics of enhanced diffusion: application to turbulence. Rev Geophys (1987) 58:1101. 10.1103/PhysRevLett.58.1100
5.
Ramos-FernundezGMateosJMiramontesOCochoGLarraldeHAyala-OrozcoB. Levy walk patterns in the foraging movements of spidermonkeys (Ateles geoffroyi). Behav Ecol Sociobiol (2004) 55:223. 10.1007/s00265-003-0700-6
6.
ViswanathanGMAfanasyevVBuldyrevSVMurphyEJPrincePAStanleyHE. Lévy flight search patterns of wandering albatrosses. Nature (1996) 381:413. 10.1038/381413a0
7.
ViswanathanGMBuldyrevSVHavlinSda LuzMGERaposoEPStanleyHE. Optimizing the success of random searches. Nature (1999) 401:911914. 10.1038/44831
8.
JayannavarAMKumarN. Nondiffusive quantum transport in a dynamically disordered medium. Phys Rev Lett (1982) 48:553. 10.1103/physrevlett.48.553
9.
GolubovicLFengSZengFA. Classical and quantum superdiffusion in a time-dependent random potential. Phys Rev Lett (1991) 67:2115. 10.1103/PhysRevLett.67.2115
10.
LeviLKrivolapovYFishmanSSegevM. Hyper-transport of light and stochastic acceleration by evolving disorder. Nat Phys (2012) 8:912. 10.1038/nphys2463
11.
SagiYBrookMAlmogIDavidsonN. Observation of anomalous diffusion and fractional self-similarity in one dimension. Phys Rev Lett (2012) 108:093002. 10.1103/physrevlett.108.093002
12.
AndersonPW. Absence of diffusion in certain random lattices. Phys Rev (1958) 109:1492. 10.1103/physrev.109.1492
13.
MuskatM. The flow of fluids through porous media. J Appl Physics (1937) 8:274. 10.1063/1.1710292
14.
BarenblattGIEntovVMRyzhikVM. Theory of fluid flows through natural rocks, Vol. ber90. Dordrecht, Netherland: Kluwer (1990). 396 p.
15.
HurtadoPIKrapivskyPL. Compact waves in microscopic nonlinear diffusion. Phys Rev E (2012) 85:060103. 10.1103/physreve.85.060103
16.
KamenomostskayaSK. Similar solutions and the asymptotics of filtration equations. Arch Rational Mech Anal (1976) 60, 171. 10.1007/bf00250678
17.
ZeldovichIBKompaneezAS. “On the theory of heat propagation with heatconduction depending on temperature,” in Lectures dedicated on the 70th anniversary of A. F. Joffe. Moscow, Russia: Akad. Nauk SSSR (1950). p. 61–71.
18.
ZeldovichYRaizerYP. Physics of shock waves and high-temperature hydrodynamic phenomena, Vol. 2. New York, NY: Academic Press (1968). p. 826–7.
19.
GurtinMEMacCamyRC. On the diffusion of biological populations. Math Biosci (1977) 33:35. 10.1016/0025-5564(77)90062-1
20.
NewmanWI. Some exact solutions to a non-linear diffusion problem in population genetics and combustion. J Theor Biol (1980) 85:325. 10.1016/0022-5193(80)90024-7
21.
MurrayJD. Mathematical biology, m89, biomathematics. Berlin, Germany: Springer (1989). 770 p.
22.
de AzevedoENde SousaPLde SouzaREEngelsbergMMirandaMDNDNSilvaMA. Concentration-dependent diffusivity and anomalous diffusion: a magnetic resonance imaging study of water ingress in porous zeolite. Phys Rev E (2006a) 73:011204. 10.1103/physreve.73.011204
23.
de AzevedoENda Silvaandde SouzaandDVEngelsbergREEngelsbergM. Water ingress in y-type zeolite: anomalous moisture-dependent transport diffusivity. Phys Rev E(2006b) 74:041108. 10.1103/physreve.74.041108
24.
FischerDFingerTAngensteinFRSR. Diffusive and subdiffusive axial transport of granular material in rotating mixers. Phys Rev E (2009) 80:061302. 10.1103/physreve.80.061302
25.
StannariusICHasHA. Resolving a paradox of anomalous scalings in the diffusion of granular materials. Proc Natl Acad Sci (2012) 109:16012. 10.1073/pnas.1211110109
26.
StoneDWoodsAWHoggAJ. On the slow draining of a gravity current moving through a layered permeable medium. J Fluid Mech (2001) 444:23. 10.1017/s002211200100516x
27.
HansenASkagerstamB-STøråG. Anomalous scaling and solitary waves in systems with nonlinear diffusion. Phys Rev E (2011) 83:056314. 10.1103/physreve.83.056314
28.
ShlesingerMF. Asymptotic solutions of continuous-time random walks. J Stat Phys (1974) 10:421. 10.1007/bf01008803
29.
HavlinSKieferJEWeissGH. Anomalous diffusion on a random comblike structure. Phys Rev A (1987) 36:1403. 10.1103/physreva.36.1403
30.
HavlinSDjordjevicZVMajidIStanleyHEWeissGH. Relation between dynamic transport properties and static topological structure for the lattice-animal model of branched polymers. Phys Rev Lett (1984) 53:178. 10.1103/physrevlett.53.178
31.
de GennesPG. La percolation un concept unificateur, La Recherche. Singapore: World Scientific Publishing (1976). 919 p.
32.
ben-AvrahamDHavlinS. Diffusion and reactions in fractals and disordered systems. London, United Kingdom: Cambridge University Press (2000). 316 p.
33.
ben-AvrahamDLeyvrazFRednerS. Superballistic motion in a “random-walk” shear flow. Phys Rev A (1992) 45:2315. 10.1103/physreva.45.2315
34.
ben-NaimERednerSben-AvrahamD. Bimodal diffusion in power-law shear flows. Phys Rev A (1992) 45:7207. 10.1103/physreva.45.7207
35.
SchützGMTrimperS. Elephants can always remember: exact long-range memory effects in a non-markovian random walk. Phys Rev E (2004) 70:045101. 10.1103/physreve.70.045101
36.
MorgadoROliveiraFABatrouniGGHansenA. Relation between anomalous and normal diffusion in systems with memory. Phys Rev Lett (2002) 89:100601. 10.1103/physrevlett.89.100601
37.
WeiXBauerW. Starvation-induced changes in motility, chemotaxis, and flagellation of rhizobium melilot. Appl Abd Environ Microbiol (1998) 65:1708–14. 10.1128/AEM.64.5.1708-1714.1998
38.
WuX-LLibchaberA. Particle diffusion in a quasi-two-dimensional bacterial bath. Phys Rev Lett (2000) 84:3017. 10.1103/physrevlett.84.3017
39.
PattleRE. Diffusion from an instantaneous point source with a concentration-dependent coefficient. Q J Mech Appl Math (1959) 12:407. 10.1093/qjmam/12.4.407
40.
HansenAFlekkøyEGBaldelliB. Anomalous diffusion in systems with concentration-dependent diffusivity: exact solutions and particle simulations. Front Phys (2020) 8:519624. 10.3389/fphy.2020.519624
41.
PlastinoARPlastinoA. Non-extensive statistical mechanics and generalized Fokker-Planck equation. Physica A (1995) 222:347. 10.1016/0378-4371(95)00211-1
42.
TsallisCBukmanDJ. Anomalous diffusion in the presence of external forces: exact time-dependent solutions and their thermostatistical basis. Phys Rev E (1996) 54:R2197. 10.1103/physreve.54.r2197
43.
van KampenN. Stochastic processes in physics and chemistry. 3rd ed, Vol. k07. Amsterdam: North Holland (2007). 464 p.
44.
BorlandL. Microscopic dynamics of the nonlinear Fokker-Planck equation: a phenomenological model. Phys Rev E (1998) 57:6634. 10.1103/physreve.57.6634
Appendix
In Pattles classical 1959 paper [39] the solution of Eq. 1 is not actually derived, but only written down. So, for completeness we derive it here along the same lines as those leading up to Eq. 11. In the solutions thus derived has a final support outside which it is strictly zero. For any γ the normalizationwhere we have used the isotropic nature of the problem to perform the angular integration and thus introduced the geometric factor when .
We see from Eq. 11 that, for , the domain of the probability density, , is limited to , so that the normalization condition isyielding the normalization constantfor .
We now combine results, using Eqs. 5, 7; 11. to find the concentration field, whereBy comparison, for we havewith k given by Eq. 12 now. Finally we find that for , withIn Figure 7 this behavior is confirmed by simulations using the finite-range model.
Summary
Keywords
anomalous diffusion, concentration-dependent diffusivity, non-linear diffusion equation, brownian motion (wiener process), random walks
Citation
Flekkøy EG, Hansen A and Baldelli B (2021) Hyperballistic Superdiffusion and Explosive Solutions to the Non-Linear Diffusion Equation. Front. Phys. 9:640560. doi: 10.3389/fphy.2021.640560
Received
11 December 2020
Accepted
18 January 2021
Published
17 March 2021
Volume
9 - 2021
Edited by
Fernando A. Oliveira, University of Brasilia, Brazil
Reviewed by
Marie-Christine Firpo, Center National de la Recherche Scientifique (CNRS), France
Haroldo V. Ribeiro, State University of Maringá, Brazil
Updates
Copyright
© 2021 Flekkøy, Hansen and Baldelli.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Eirik G. Flekkøy, flekkoy@fys.uio.no
This article was submitted to Interdisciplinary Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.