Parabolic regularity and scaling
Match space-time scales, quantify positive-time derivative control, and identify the regularity that coefficients can support.
On this page
The question for today
What you’ll understand
- Use the quadratic time scale associated with diffusion.
- Interpret an interior derivative estimate without claiming time-zero or boundary regularity.
- Construct a uniformly elliptic rough-coefficient example whose gradient is discontinuous.
Diffusion links a spatial length to the square root of elapsed time. Regularity estimates respect this geometry: one time derivative counts like two spatial derivatives. This scaling helps organize both derivative bounds and Hölder spaces. It also clarifies the limits of smoothing. Constant-coefficient heat flow is exceptionally regular, but discontinuous coefficients can prevent classical spatial smoothness even when diffusion remains uniformly nondegenerate.
The natural cylinders are not symmetric in space and time
For uₜ=Δu, the rescaled function uᵣ(x,t)=u(rx,r²t) satisfies the same equation on the correspondingly rescaled domain. One spatial derivative gains a factor r and one time derivative gains r². Thus an interior estimate naturally uses a backward cylinder with spatial radius r and time depth r². The word backward here describes the neighborhood feeding a forward estimate; it does not mean the ill-posed backward heat evolution is being solved. With a fixed positive diffusivity κ, the corresponding diffusion time scale is proportional to r²/κ.
For a smooth homogeneous heat solution on a cylinder twice as large, an interior gradient estimate bounds the gradient on the smaller cylinder by a constant times r⁻¹ times the solution supremum. Higher spatial and time derivatives receive their matching powers. Constants depend on dimension and the chosen nested-cylinder proportions. The estimate is local and stays away from the initial or lateral boundary. It cannot be extended to t=0 with the same constant for arbitrary rough data, nor can it guarantee smooth attainment of incompatible boundary and initial values.
Coefficient regularity determines which estimates are available
Parabolic Hölder spaces compare space increments with the square root of time increments. For a scalar nondivergence equation uₜ−aᵢⱼuᵢⱼ=f, uniformly positive definite coefficients that are Hölder continuous with exponents α in space and α/2 in time, together with similarly Hölder forcing, support interior Schauder estimates. In their a priori form these bound the second spatial and first time derivatives, including their Hölder behavior, by the solution size and forcing norm on a larger cylinder. Such an estimate assumes an appropriate solution and does not alone construct one.
Uniform ellipticity by itself does not make discontinuous coefficients smooth. A divergence-form equation with a jump in conductivity can have a continuous temperature but a jump in its gradient while its flux remains continuous. The weak formulation is designed to accommodate that behavior. Differentiating the coefficient as if it were a classical smooth function would introduce unjustified terms. Distinguish divergence from nondivergence form when invoking a regularity theorem, and match the data spaces and coefficient assumptions to that theorem rather than borrowing a conclusion from constant-coefficient heat flow.
Smoothing and existence are different stages of an argument
A regularity argument usually begins with a solution obtained by a kernel, Galerkin method, or another existence construction. It then proves additional local properties under suitable assumptions. A global classical solution up to a boundary needs more: boundary smoothness, regularity of prescribed values, and compatibility where the initial slice meets the lateral boundary. An interior derivative estimate does not address those corners. Similarly, a nonlinear equation may be uniformly parabolic only while its solution remains in a bounded range separated from degeneracy.
Fourier modes provide a quantitative view of how high derivatives become controlled after positive time. A mode n has gradient size proportional to n exp(−n²t) for unit diffusivity. Initially this grows with n, but at fixed positive time the exponential dominates. Maximizing the analogous continuous-frequency expression gives a scale proportional to , consistent with the interior estimate. This explains how initial high-frequency roughness can disappear without any uniform derivative bound at the initial instant. The argument describes smoothing, not a stable inversion of the lost frequencies.
Work it through
Find when a finer mode has the smaller gradient
Compare u₁=e⁻ᵗsin x and u₂=e⁻⁴ᵗsin2x, both solving the unit-diffusivity periodic heat equation. When do their gradient suprema agree?
The gradient suprema are e⁻ᵗ and 2e⁻⁴ᵗ. At time zero the finer mode has twice the gradient size.
Equality requires 2e⁻⁴ᵗ=e⁻ᵗ, or e³ᵗ=2, giving t=(log2)/3. After that time the second mode has the smaller gradient.
At t=1/4 the values are exp(−1/4)≈0.7788 and 2e⁻¹≈0.7358, verifying the crossover. The spatial-frequency factor and the quadratic damping rate must both be retained.
The crossover occurs at (log2)/3. Faster decay can outweigh an initially larger derivative without contradicting time-zero roughness.
The living notebook
Explore the model
Interactive lab · heat
Check spatial and temporal derivative scales
Reset, set n=3 and κ=0.5, then Step once to t=0.5. Predict amplitude exp(−2.25), gradient supremum 3exp(−2.25), and time-derivative supremum 4.5exp(−2.25). The readout gives amplitude and decay rate; multiply them to check the temporal prediction. This analytic periodic mode does not test a rough coefficient or a boundary regularity theorem.
- Current amplitude
- 1
- Spatial L² norm
- 1.772
- Decay rate · κn²
- 3.6 /s
The exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forSpatial differentiation contributes one frequency factor, while temporal differentiation contributes a squared-frequency diffusion rate.
Pause & reason
Check your understanding
Your turn
Try it yourself
On (−1,1), let a(x)=1 for x<0 and a(x)=2 for x>0. Set u(x)=x for x≤0 and u(x)=x/2 for x≥0, independent of time, with its fixed endpoint values. Is it a weak solution of uₜ−(auₓ)ₓ=0? Is it C¹ across zero?
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Compute the flux auₓ on each side before differentiating it distributionally.
Compare with the solution
The function is continuous and belongs to H¹. Its flux equals 1 on both sides, so its distributional derivative is zero; with uₜ=0 it is a stationary weak solution and has the stated endpoint traces. Its one-sided derivatives are 1 and 1/2, so it is not C¹ across zero. The coefficient is bounded and uniformly elliptic, but its jump prevents an automatic classical-gradient smoothing conclusion.
Keep the idea
Take it with you
- Diffusive time scales quadratically with spatial length.
- Interior estimates and global boundary regularity require different information.
- Uniform parabolicity does not erase rough coefficient structure.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.