TY - JOUR

T1 - Revisiting diffusion

T2 - Self-similar solutions and the t-1/2 decay in initial and initial-boundary value problems

AU - Kevrekidis, P. G.

AU - Williams, M. O.

AU - Mantzavinos, D.

AU - Charalampidis, E. G.

AU - Choi, M.

AU - Kevrekidis, I. G.

N1 - Funding Information:
Acknowledgments. The authors are grateful to C. Eugene Wayne for numerous enlightening discussions on the subject and useful comments on the manuscript. The first author gratefully acknowledges the support of NSF-DMS-1312856, NSF-PHY-1602994 and the ERC under FP7, Marie Curie Actions, People, International Research Staff Exchange Scheme (IRSES-605096). The research of the sixth author was partially supported by the U.S. National Science Foundation and by the IAS-TU Munich through a Hans Fischer Senior Fellowship. The third and fourth authors would like to thank the Department of Chemical and Biological Engineering at Princeton University for the kind hospitality.
Publisher Copyright:
© 2017 Brown University.

PY - 2017

Y1 - 2017

N2 - The diffusion equation is a universal and standard textbook model for partial differential equations (PDEs). In this work, we revisit its solutions, seeking, in particular, self-similar profiles. This problem connects to the classical theory of special functions and, more specifically, to the Hermite as well as the Kummer hypergeometric functions. Reconstructing the solution of the original diffusion model from self-similar solutions of the associated self-similar PDE, we infer that the t-1/2 decay law of the diffusion amplitude is not necessary. In particular, it is possible to engineer setups of both the Cauchy problem and the initial-boundary value problem in which the solution decays at a different rate. Nevertheless, we observe that the t-1/2 rate corresponds to the dominant decay mode among integrable initial data, i.e., ones corresponding to finite mass. Hence, unless the projection to such a mode is eliminated, generically this decay will be the slowest one observed. In initial-boundary value problems, an additional issue that arises is whether the boundary data are consonant with the initial data; namely, whether the boundary data agree at all times with the solution of the Cauchy problem associated with the same initial data, when this solution is evaluated at the boundary of the domain. In that case, the power law dictated by the solution of the Cauchy problem will be selected. On the other hand, in the non-consonant cases a decomposition of the problem into a self-similar and a non-self-similar one is seen to be beneficial in obtaining a systematic understanding of the resulting solution.

AB - The diffusion equation is a universal and standard textbook model for partial differential equations (PDEs). In this work, we revisit its solutions, seeking, in particular, self-similar profiles. This problem connects to the classical theory of special functions and, more specifically, to the Hermite as well as the Kummer hypergeometric functions. Reconstructing the solution of the original diffusion model from self-similar solutions of the associated self-similar PDE, we infer that the t-1/2 decay law of the diffusion amplitude is not necessary. In particular, it is possible to engineer setups of both the Cauchy problem and the initial-boundary value problem in which the solution decays at a different rate. Nevertheless, we observe that the t-1/2 rate corresponds to the dominant decay mode among integrable initial data, i.e., ones corresponding to finite mass. Hence, unless the projection to such a mode is eliminated, generically this decay will be the slowest one observed. In initial-boundary value problems, an additional issue that arises is whether the boundary data are consonant with the initial data; namely, whether the boundary data agree at all times with the solution of the Cauchy problem associated with the same initial data, when this solution is evaluated at the boundary of the domain. In that case, the power law dictated by the solution of the Cauchy problem will be selected. On the other hand, in the non-consonant cases a decomposition of the problem into a self-similar and a non-self-similar one is seen to be beneficial in obtaining a systematic understanding of the resulting solution.

UR - http://www.scopus.com/inward/record.url?scp=85028572543&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85028572543&partnerID=8YFLogxK

U2 - 10.1090/qam/1473

DO - 10.1090/qam/1473

M3 - Article

AN - SCOPUS:85028572543

VL - 75

SP - 581

EP - 598

JO - Quarterly of Applied Mathematics

JF - Quarterly of Applied Mathematics

SN - 0033-569X

IS - 4

ER -