Discretely Self-Similar Solutions in Nonlinear PDEs

Updated 22 January 2026

Discretely self-similar solutions are exact invariants of nonlinear PDEs under discrete scaling, recovering their shape only at specific scales.
They model critical phenomena like singularity formation, gravitational collapse, and shock dynamics by exhibiting periodic behavior in rescaled variables.
Analytical and numerical techniques such as similarity variables, energy estimates, and spectral analysis underpin their construction and stability investigations.

A discretely self-similar (DSS) solution is an exact solution to a nonlinear partial differential equation (PDE) that is invariant under a discrete one-parameter scaling group, but not under the full continuous scaling group. DSS solutions appear in a range of contexts, including fluid dynamics (Navier–Stokes, Euler, MHD), geometric PDEs (wave maps, Einstein equations), and critical phenomena in gravitational collapse. They serve as key models for singularity formation, critical blow-up, and as attractors or non-attractors for dynamics near singularities.

1. Definition and Structural Properties

Let $\lambda > 1$ be a fixed scale factor. Given a system with scaling symmetry (e.g., Navier–Stokes, Euler), a solution $u(x,t)$ is called $\lambda$ -discretely self-similar (DSS) if

$u(\lambda x, \lambda^k t) = \lambda^{\alpha k} u(x,t)$

for a specific $\alpha$ determined by the scaling of the equation (often, $\alpha = -1$ in critical cases) and all (suitable) integers $k$ . In practice,

For the 3D Navier–Stokes equations, DSS with factor $\lambda$ means $u(\lambda x, \lambda^2 t) = \lambda u(x,t)$ , with the pressure scaling as $p(\lambda x, \lambda^2 t) = \lambda^2 p(x,t)$ .
For the Euler equations, the scaling exponents involve the blow-up rate; e.g., for exponent $\alpha$ , $v(\lambda x, \lambda^{\alpha+1} t) = \lambda^\alpha v(x,t)$ .

A DSS solution is not fixed under all positive dilations (not continuous self-similar), but recovers its shape at the specific scales defined by $\lambda$ . In self-similar (or similarity) variables, DSS corresponds to time-periodic solutions in a rescaled dynamic, with period $T$ in the logarithmic time variable $s$ (e.g., $T = 2\log \lambda$ for parabolic scaling) (Bradshaw et al., 2015, Tsai, 2012).

2. DSS Solutions in Fluid Dynamics

In the context of the 3D incompressible Navier–Stokes equations, there is an extensive theory on the existence and qualitative properties of DSS solutions:

Existence: For arbitrary large DSS initial data in the critical weak $L^{3,\infty}$ space, global-in-time DSS local Leray solutions exist; see (Bradshaw et al., 2015, Tsai, 2012). The framework adapts to data in critical Besov spaces $\dot{B}^{-1+3/p}_{p,\infty}$ ($3 < p < 6$) (Bradshaw et al., 2017), and $L^2_{\mathrm{loc}}$ (Chae et al., 2016, Bradshaw et al., 2018).
Uniqueness: Uniqueness is known in the small-data regime (e.g., small $L^{3,\infty}$ ), but is open for large data or for weak Leray-type solutions (Bradshaw et al., 2015, Tsai, 2012).
Regularity and Decay: DSS solutions are smooth away from the origin in space and regular for $|x| > R_0\sqrt{t}$ . Sharp pointwise decay rates have recently been established, with the nonlinear part decaying strictly faster than the linear evolution (Bradshaw et al., 2022, Bradshaw et al., 2024). For instance, for rough initial data in $L^q_{\mathrm{loc}}$ (with $q > 3$ ), the solution decomposes into an $O((|x|+\sqrt t)^{-1})$ term and a part that is $o((|x|+\sqrt t)^{-1})$ (Bradshaw et al., 2024).
Besov Data: The class of initial data supporting DSS solutions can be extended to critical Besov spaces, where nontrivial solutions exist for arbitrarily large data (Bradshaw et al., 2017).
Physical Interpretation: DSS solutions describe flows that "echo" themselves at discrete time intervals (e.g., $t \to \lambda^2 t$ ), capturing intermittent or intermittently singular dynamics not visible in ordinary self-similar solutions (Bradshaw et al., 2015, Bradshaw et al., 2017).

2.2 Euler and Magnetohydrodynamics

Nonexistence of DSS Blow-up: In the Euler case, as well as for inviscid MHD, robust nonexistence theorems show that under suitable decay or integrability of the profile, there are no nontrivial DSS blow-up solutions; any such profile must be identically zero (Chae, 2013, Xue, 2014, Chae, 2013). If a DSS solution vanishes in an open subset for a full period in similarity variables, it must vanish identically—a strong unique continuation property (Chae, 2013).
Infinite Energy and Weighted Framework: In MHD, infinite-energy DSS solutions exist globally in time for weighted $L^2_{w_\gamma}$ initial data with $\gamma < 2$ (Fernández-Dalgo et al., 2019). These solutions are constructed in weighted spaces and can display nontrivial DSS behavior even when global energy is infinite.

2.3 Oberbeck–Boussinesq

The Oberbeck–Boussinesq system with Newtonian gravity also admits global-in-time large data DSS solutions, with precise a priori energy bounds in similarity variables and no smallness assumption on the norms of the initial data (Tsai, 2024).

3. DSS Solutions in Geometric and Relativistic PDEs

3.1 Wave Maps and Critical Blowup

For wave maps (e.g., from $\mathbb{R}^{1+d}$ to $\mathbb{S}^1$ ), there exists a countable family of discretely self-similar blow-up profiles, constructed explicitly via similarity variables and hypergeometric functions (Glogić et al., 18 Dec 2025). The nonlinear stability of all such DSS blow-up profiles (for all dimensions and for both even and radial data) is established via resolvent construction, spectral analysis, and Lyapunov–Perron methods. The codimension of instability is explicitly computed for each profile.

3.2 General Relativity and Critical Collapse

Einstein–Klein–Gordon at Large $D$ : Analytic closed-form DSS solutions of the Einstein–Klein–Gordon system at large spacetime dimension $D$ have been constructed using systematic $1/D$ expansions, revealing universal and non-universal features, echoing periods, mass scaling laws, and connections to numerically observed Choptuik critical solutions (Ecker et al., 20 Jan 2026).
Exterior-Naked Singularities: In spherically symmetric Einstein–scalar field systems, DSS "exterior-naked singularity" regions have been constructed that are smooth up to the past light cone of the singularity. These display bounded scalar fields but infinite oscillations in both the scalar field and the mass-aspect ratio; they differ from previously known Christodoulou CSS constructions by allowing $C^\infty$ regularity and capturing infinite-frequency behavior in the singular limit. Open questions remain regarding their connection to global DSS spacetimes arising from smooth Cauchy data (Cicortas et al., 2024).

4. DSS in Shock Dynamics and Explosions

Strong Explosion Problem: For perturbations of classical Sedov–Taylor blast waves with a log-periodic density profile, discrete self-similarity arises naturally, and the perturbations repeat themselves after discrete rescalings in time and space. These "echoing" DSS modes have direct relevance for astrophysical situations involving resonant interaction of blast waves with log-periodic density inhomogeneities (Yalinewich et al., 2014).

5. Nonexistence, Uniqueness Theorems, and Structural Rigidity

Several nonexistence results and unique continuation properties have been proven for DSS solutions:

Euler and Inviscid MHD: No nontrivial DSS blow-up profile exists under physically reasonable (decay or integrability) assumptions on the profile function, both for decaying and mild algebraic growth at infinity. Profiles satisfying periodicity and certain energy inequalities must be trivial (Chae, 2013, Xue, 2014, Chae, 2013).
Unique continuation: If a DSS solution (in self-similar or similarity variables) vanishes in a spatial neighborhood for one period, it vanishes globally (Chae, 2013).
Absence of DSS Bifurcation: In stationary Navier–Stokes, no bifurcation of swirling DSS solutions branches off the family of continuously self-similar Landau solutions. Numerical and analytical evidence confirms no small DSS branches emanate from the Landau family (Kwon et al., 2020).
Wave Maps and Stability: All DSS blow-up profiles for the wave maps-type system are codimensionally unstable except for the fundamental mode, with precise stability indices derived from spectral theory (Glogić et al., 18 Dec 2025).

6. Analytical and Numerical Construction Techniques

Similarity Variables: Rewriting the PDEs in rescaled ( $y,s$ ) variables is universal, yielding time-periodic PDEs for the profile functions. The DSS property is translated into periodicity in the similarity-time coordinate.
Energy Estimates: Explicit a priori energy inequalities, weighted by cut-off functions and robust under limiting and approximation procedures, yield existence and regularity (Bradshaw et al., 2015, Tsai, 2024).
Galerkin–Fixed Point Methods: Construction frequently proceeds via finite-dimensional Galerkin approximations and Brouwer or Schauder fixed-point theorems in function spaces of time-periodic divergence-free fields, both for fluid and geometric equations (Bradshaw et al., 2015, Bradshaw et al., 2017, Glogić et al., 18 Dec 2025).
Splitting and Decomposition: Initial data are split into parts with desirable decay or integrability, and solutions are decomposed into leading-order profiles and rapidly decaying remainders; this enables refined control over spatial asymptotics (Bradshaw et al., 2022, Bradshaw et al., 2024).
Spectral Theory: Nonlinear stability and instability (codimension) in the wave maps and gravitational collapse contexts are established through the construction of explicit resolvents, Volterra-type solutions, and calculation of unstable modes (Glogić et al., 18 Dec 2025, Ecker et al., 20 Jan 2026).
Numerical Matching: In cases where analytic solutions are not available, matching conditions between boundary layers and interior/near-singularity regions yield quantization of echoing periods, as in the Primakoff limit and critical gravitational collapse (Yalinewich et al., 2014, Ecker et al., 20 Jan 2026).

7. Applications, Physical Interpretation, and Open Problems

Role in Critical Phenomena: DSS solutions function as attractors (or non-attractors) at the threshold of blow-up or black hole formation, with criticality characterized by the echoing period and scaling exponents (Ecker et al., 20 Jan 2026, Cicortas et al., 2024).
Separation of Non-Unique Solutions: For certain initial data, the maximal separation rate between DSS solutions (reflecting possible non-uniqueness) is sharply constrained by spatial decay estimates (Bradshaw et al., 2022, Bradshaw et al., 2024).
Connection to Singularities and Regularity: Incompressible Euler and related systems admit strong rigidity: DSS profiles with sufficient integrability are forced to be trivial, sharply constraining possible singularity scenarios (Chae, 2013, Xue, 2014, Chae, 2013).
Boundary and Interior Matching in GR: For exterior DSS naked singularity solutions, the construction of an interior fill-in from regular data on an asymptotically flat slice remains unresolved, as does the global realization of numerically observed universal DSS attractors (Cicortas et al., 2024).
Generic Global Attractors and Nonuniqueness: For NLS, Euler, MHD, and geometric flows, understanding which DSS solutions serve as universal blow-up or critical attractors remains a central challenge; uniqueness questions for large-data DSS flows in supercritical function spaces remain wide open (Bradshaw et al., 2015, Bradshaw et al., 2017, Glogić et al., 18 Dec 2025).

References:

(Bradshaw et al., 2015) Forward discretely self-similar solutions of the Navier-Stokes equations II
(Bradshaw et al., 2017) Discretely self-similar solutions to the Navier-Stokes equations with Besov space data
(Tsai, 2012) Forward Discretely Self-Similar Solutions of the Navier-Stokes Equations
(Chae et al., 2016) Existence of discretely self-similar solutions to the Navier-Stokes equations for initial value in $L^2_{\mathrm{loc}}$
(Bradshaw et al., 2018) Discretely self-similar solutions to the Navier-Stokes equations with data in $L^2_{\mathrm{loc}}$
(Bradshaw et al., 2022) Spatial decay of discretely self-similar solutions to the Navier-Stokes equations
(Bradshaw et al., 2024) Asymptotic properties of discretely self-similar Navier-Stokes solutions with rough data
(Chae, 2013) Remarks on the asymptotically discretely self-similar solutions of the Navier-Stokes and the Euler equations
(Chae, 2013) Continuation of the zero set for discretely self-similar solutions to the Euler equations
(Xue, 2014) Discretely self-similar singular solutions for the incompressible Euler equations
(Fernández-Dalgo et al., 2019) Existence of infinite-energy and discretely self-similar global weak solutions for 3D MHD equations
(Kwon et al., 2020) On bifurcation of self-similar solutions of the stationary Navier-Stokes equations
(Glogić et al., 18 Dec 2025) Existence and stability of discretely self-similar blowup for a wave maps type equation
(Ecker et al., 20 Jan 2026) Analytic discrete self-similar solutions of Einstein-Klein-Gordon at large D
(Cicortas et al., 2024) Discretely self-similar exterior-naked singularities for the Einstein-scalar field system
(Yalinewich et al., 2014) Discrete Self Similarity in Filled Type I Strong Explosions
(Tsai, 2024) Large discretely self-similar solutions to Oberbeck-Boussinesq system with Newtonian gravitational field
(Lai, 2019) Forward Discretely Self-Similar Solutions of the MHD Equations and the Viscoelastic Navier-Stokes Equations with Damping