Numerical renormalization group algorithms for self-similar solutions of partial differential equations
For researchers studying self-similar solutions of PDEs, this work validates and extends a numerical method that can handle models resistant to traditional analysis, though the approach is not entirely new.
The paper validates a numerical renormalization group (nRG) algorithm for revealing self-similar dynamics in PDEs, showing it drives solutions to fixed points exponentially fast. It demonstrates the algorithm's applicability across various PDEs, enabling analysis of long-time behavior in models that are difficult to study otherwise.
We systematically study a numerical procedure that reveals the asymptotically self-similar dynamics of solutions of partial differential equations (PDEs). This procedure, based on the renormalization group (RG) theory for PDEs, appeared initially in a conference proceeding by Braga et al. \cite{BFI04}. This numerical version of RG method, dubbed as the numerical RG (nRG) algorithm, numerically rescales the temporal and spatial variables in each iteration and drives the solutions to a fixed point exponentially fast, which corresponds to the self-similar dynamics of the equations. In this paper, we carefully examine and validate this class of algorithms by comparing the numerical solutions with either the exact or the asymptotic solutions of the model equations in literature. The other contribution of the current paper is that we present several examples to demonstrate that this class of nRG algorithms can be applied to a wide range of PDEs to shed lights on longtime self-similar dynamics of certain physical models that are difficult to analyze, both numerically and analytically.