Stochastic wave equations and regularity
Calculate displacement and velocity covariances and separate oscillatory wave propagation from heat-semigroup smoothing.
On this page
The question for today
What you’ll understand
- Write the spectral stochastic wave convolution in displacement and velocity.
- Check their distinct state-space regularities.
- Calculate expected mode energy without confusing oscillation with damping.
Replacing heat diffusion by wave propagation changes the response operator even when the driving noise looks identical. A wave mode oscillates rather than decaying, and its displacement and velocity occupy different spatial regularity spaces. The displacement response includes an inverse frequency, while the velocity response does not. Computing both covariances prevents the common mistake of treating a well-defined displacement field as evidence of finite total wave energy. The linear Dirichlet model also shows how persistent white forcing injects energy into every mode.
The state contains two coupled components
On (0,pi), write the wave equation with Dirichlet Laplacian as dU=Vdt and , with zero initial displacement and velocity and unit cylindrical forcing in the velocity equation. In normalized sine coordinates, the kth displacement is the Brownian integral with response , while its velocity uses . These formulas define scalar coordinates first; a state-space solution requires summing their covariance contributions.
The deterministic wave group acts naturally on paired Sobolev spaces with one derivative of difference between displacement and velocity. It conserves an appropriate quadratic energy rather than regularizing arbitrary initial data as an analytic heat semigroup does. An inverse frequency in the forced displacement response is useful, but it is not exponential high-frequency damping. Nonlinear wave equations, boundary forcing and spatially correlated noise can require other spaces or methods; the present coordinate calculation is a controlled linear case.
Displacement and velocity have different covariance tails
The Itô isometry gives and . At any fixed positive time, displacement variance has order k inverse squared, whereas velocity variance approaches a nonzero constant. Thus the displacement is L²-valued in one dimension, but the velocity is not. A finite collection of coordinate velocities cannot reveal this infinite tail obstruction.
Using the spectral convention from the heat lesson, the paired state belongs in mean square to H_A^s times H_A^(s−1) exactly for s less than one-half on the interval. The expected ordinary energy in H0¹ times L² diverges for unit white forcing. On a smooth bounded d-dimensional Dirichlet domain, Weyl counting gives the analogous fixed-time threshold s less than one minus d/2. These are Sobolev trace calculations, not a claim about pointwise path regularity or a nonlinear random-field existence theorem in every dimension.
Noise injection and regularity are not equilibrium statements
For one unit-forced mode, expected velocity squared plus k squared times displacement squared equals t; the sine corrections cancel. The physical energy with an extra factor one-half therefore has mean t/2. An undamped wave mode has no stationary covariance balancing this growth. Adding damping changes the stochastic dynamics and can permit equilibrium, but it is a new equation. A deterministic animation with bounded amplitude does not contradict energy growth under repeated random forcing.
Regularity conclusions also depend on where the noise enters. Spatial coloring reduces the covariance of high modes, while forcing the displacement equation instead changes the response operators. A numerical method that evolves only a few smooth modes will always show finite velocity energy at that cutoff. To claim convergence, examine covariance tails in the intended paired norm and control the time approximation separately. Oscillatory cancellation inside an ordinary integral cannot be assumed to cancel the nonnegative squared responses appearing in stochastic variance.
Work it through
A phase where covariance corrections vanish
In the unit-forced mode k=2, evaluate displacement variance, velocity variance and unhalved quadratic energy at t=pi/2.
The oscillatory correction contains sin(2kt)=sin(2pi)=0. The displacement variance is t/(2k²)=(pi/2)/8=pi/16.
The velocity variance is t/2=pi/4. The covariance components are different even though they come from the same noise-driven oscillator.
The expected unhalved energy is pi/4+4×pi/16=pi/2, exactly t. Including the conventional energy factor one-half would give pi/4 instead, so the normalization must be stated.
Each unit-forced mode absorbs expected unhalved energy at rate one, regardless of its frequency.
The living notebook
Explore the model
Interactive lab · oscillator
Undamped propagation preserves deterministic amplitude
Pause and Reset, then set both controls to damping ratio ζ=0 and natural frequency ω0=1. Two Steps give time 1; predict position cos(1), about 0.540302. Reset and set ζ=0 and ω0=2; after two Steps predict cos(2), about −0.416147. The instrument starts at displacement one and velocity zero. It shows deterministic oscillation, not stochastic forcing or the growing covariance of a forced wave mode.
- Displacement x(t)
- 1
- Velocity v(t)
- 0
- Regime
- underdamped
Unit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. 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 forWave oscillation and heat damping have different deterministic responses, and noise energy must be calculated separately.
Pause & reason
Check your understanding
Your turn
Try it yourself
At time pi, compare modes one and two under unit velocity forcing. Find both displacement variances, both velocity variances and their unhalved expected energies.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Both sine corrections vanish at this time. Keep the frequency weight in the displacement contribution to energy.
Compare with the solution
For k=1, displacement and velocity variances are both pi/2, giving energy pi. For k=2, displacement variance is pi/8 while velocity variance is still pi/2; weighting the former by four again gives energy pi. Higher frequency reduces displacement size but not the energy injection per unit-forced mode. Summing that equal positive energy over infinitely many modes diverges, even though the displacement variance sum alone converges in one dimension.
Keep the idea
Take it with you
- A wave state pairs displacement with a rougher velocity component.
- Oscillation does not replace a stochastic covariance summability test.
- Finite displacement norm does not establish finite wave energy or an invariant law.
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.