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 - 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
SN - 0033-569X
VL - 75
SP - 581
EP - 598
JO - Quarterly of Applied Mathematics
JF - Quarterly of Applied Mathematics
IS - 4
ER -