Explore the wave equation
Let’s next look at the wave equation,
\[\pdd{u}{t}=D \nabla^2 u,\]on a rectangular 2D domain with homogeneous Neumann (aka noflux) boundary conditions,
\[\pd{u}{x}(0,y,t) = \pd{u}{x}(L_x,y,t) = \pd{u}{y}(x,0,t) = \pd{u}{y}(x,L_y,t) = 0.\]
Load the interactive simulation, which has been set up for this tutorial.

Click on the screen to visualise a disturbance in the medium which will then propagate in all directions.

Now press
, paint some initial data, and then press to set it in motion. 
Press
to clear the screen.
Playing with the diffusion coefficient, $D$
What does changing the diffusion coefficient, $D$, do?

Change its value by clicking and editing the value of $D$: try increasing it by a factor of 10 or even 100. → Parameters

Now click again on the screen and see how fast the disturbance spreads out throughout the domain.
Explore how the speed depends on the diffusion coefficient.
Numerical notes
The VisualPDE solver only works for systems of firstorder (in time) equations. So in fact what is being simulated is the system
\[\begin{aligned}\pd{u}{t}&=v+CD\nabla^2 u,\\ \pd{v}{t} &= D \nabla^2 u, \end{aligned}\]which becomes the wave equation when $C=0$. Note that the term in the second equation is a kind of crossdiffusion, with $u$ diffusing into $v$.
The parameter $C$ is used to prevent spurious oscillations due to the equation being hyperbolic. Try varying the value of $C$ to observe how it changes the solution structure.
Standing wave solutions
If we take initial conditions of
\[\begin{aligned}u(x,y,0) &= \cos\left(\frac{n \pi x}{L_x}\right)\cos\left(\frac{m \pi y}{L_y}\right),\\ \pd{u}{t}(x,y,0)&=0,\end{aligned}\]with Neumann boundary conditions, we can find a standing wave solution of the form
\[u(x,y,t) = \cos\left(D\pi\sqrt{\frac{n^2}{L_x^2}+\frac{m^2}{L_y^2}}\,t\right)\cos\left(\frac{n \pi x}{L_x}\right)\cos\left(\frac{m \pi y}{L_y}\right),\]which oscillates in time and space.
You can play with such an initial condition in this initialised simulation, changing $n$ and $m$ in and restarting the simulation by pressing → Parameters to see how these parameters influence the solution.
The damping factor $C$ is also set to zero in this case. If you increase its value, the solution amplitude will decay over time.
If you change the boundary conditions to Dirichlet (and set $C=0.01$), the simulation will exhibit some fascinatingly symmetric oscillations.
1D and 2D waves
In a 1D domain, d’Alembert’s solution to the wave equation can be used to show that an initial disturbance in $u$ (and not $\partial u/\partial t$) will cause rightward and leftward moving waves.
 See this in a 1D simulation.
In principle the same concept works for higher dimensions, though reflections from boundary conditions can lead to more complicated behaviour, as seen in this 2D simulation plotted as a surface.