research40 min studyLesson 847 in Mathematics

Stochastic wave equations and regularity

Calculate displacement and velocity covariances and separate oscillatory wave propagation from heat-semigroup smoothing.

On this page

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 dV=ΔDUdt+dWdV=\Delta_DUdt+dW, 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 sin⁡(k(t−r))/k\sin(k(t-r))/k, while its velocity uses cos⁡(k(t−r))\cos(k(t-r)). 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.

Uk(t)=∫0tsin⁡(k(t−r))k dβk(r),Vk(t)=∫0tcos⁡(k(t−r)) dβk(r)U_k(t)=\int_0^t\frac{\sin(k(t-r))}{k}\,d\beta_k(r),\qquad V_k(t)=\int_0^t\cos(k(t-r))\,d\beta_k(r)
The same Brownian motion drives each mode's displacement and velocity, while different modes use independent Brownian motions.↔ Scroll wide equations horizontally.

Displacement and velocity have different covariance tails

The Itô isometry gives Var⁡(Uk)=t/(2k2)−sin⁡(2kt)/(4k3)\operatorname{Var}(U_k)=t/(2k^2)-\sin(2kt)/(4k^3) and Var⁡(Vk)=t/2+sin⁡(2kt)/(4k)\operatorname{Var}(V_k)=t/2+\sin(2kt)/(4k). 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.

Var⁡(Uk(t))=t2k2−sin⁡(2kt)4k3,Var⁡(Vk(t))=t2+sin⁡(2kt)4k\operatorname{Var}(U_k(t))=\frac{t}{2k^2}-\frac{\sin(2kt)}{4k^3},\qquad\operatorname{Var}(V_k(t))=\frac t2+\frac{\sin(2kt)}{4k}
The formulas are exact at fixed t for unit velocity noise. Their sum after weighting displacement by k² gives expected unhalved mode energy t.↔ Scroll wide equations horizontally.

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.

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.

  1. The oscillatory correction contains sin(2kt)=sin(2pi)=0. The displacement variance is t/(2k²)=(pi/2)/8=pi/16.

  2. The velocity variance is t/2=pi/4. The covariance components are different even though they come from the same noise-driven oscillator.

  3. 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.

The result

Each unit-forced mode absorbs expected unhalved energy at rate one, regardless of its frequency.

Explore the model

Interactive lab

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.

Undamped propagation preserves deterministic amplitudeUnit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. Current readings: Displacement x(t): 1; Velocity v(t): 0; Regime: underdamped. All parameters are adjustable with the labeled controls below.02.557.510−1.15−0.5700.571.15Time · sDisplacementunderdamped
Displacement x(t)Current time
Displacement x(t)
1
Velocity v(t)
0
Regime
underdamped
t = 0 s

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.

Check your understanding

Why can the one-dimensional white-forced wave displacement be L²-valued while its velocity is not?

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.

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.