One-Dimensional Waves

The classical wave equation in one spatial dimension reads

${\partial }_2^{}u=c^2{\partial }_2^{}u$

where $u\left(x,,,t\right)$ represents the displacement of a string, $c$ is the wave speed, and boundary conditions often fix both ends: $u\left(0,t\right)=0$ and $u\left(1,t\right)=0$. The initial shape $f\left(x\right)$ and initial velocity $g\left(x\right)$ determine the vibration. Analytical solutions exist using separation of variables, but numerical approaches shed light on more complicated setups.

Finite Difference Scheme

This calculator discretizes the string into $N$ equally spaced points with spacing $\mathrm{\Delta x}=\frac{1}{N}$. Time is stepped forward with size $\mathrm{\Delta t}$. The central difference approximation for the second time derivative and second spatial derivative yields the update formula

$u_{j,n+1}^{}=2u_{j,n}^{}-u_{j,n-1}^{}+{\frac{c\mathrm{\Delta t}}{\mathrm{\Delta x}}}^2\left(u_{j+1}^n,-,2,u_j^n,+,u_{j-1}^n\right)$

The ratio $\lambda =\frac{\mathrm{c\; \Delta t}}{\mathrm{\Delta x}}$ controls stability. For the explicit scheme here, we must choose $\lambda \le 1$ to avoid numerical blowup. The simulation initializes displacement using $f$ at $t=0$ and sets the first time step using both $f$ and $g$ to approximate the derivative.

Interpreting Results

The calculator outputs the displacement values after the specified number of time steps. Because the endpoints remain fixed at zero, standing wave patterns emerge when the initial shape aligns with sine modes. For more general initial data, waves reflect from the boundaries and interfere. Plotting the results (not shown here) would reveal how the shape evolves over time.

In physical terms, this scheme models small oscillations of a string under tension. The wave speed $c$ might depend on the string's tension and linear density. The simplest case $c=1$ corresponds to choosing units so that the tension and density factors cancel out.

Broader Context

The wave equation appears in countless applications: sound waves in air, electromagnetic waves in a cavity, vibrations of mechanical structures, and even quantum wave functions in certain limits. Numerical schemes like the one implemented here form the basis for acoustic simulators, seismology models, and physics engines in computer graphics. By experimenting with different initial shapes and velocities, you can glimpse how complex behaviors arise from simple principles.

While this tool focuses on a simple one-dimensional case, higher-dimensional wave equations follow similar discretizations. In two or three dimensions, we approximate Laplacians with finite differences or finite elements. Stability constraints become more intricate, but the underlying ideas remain the same: each time step depends on neighboring spatial points and the preceding time levels.

Beyond explicit schemes, one can use implicit or semi-implicit methods that relax stability restrictions at the cost of solving linear systems each step. More sophisticated algorithms also incorporate damping, variable coefficients, or nonlinear terms. Understanding the basic explicit method, however, is an important stepping stone.

Limitations

The present implementation prioritizes clarity over raw performance. For large grids or many time steps, more efficient loops or even GPU-based computation would be preferable. Additionally, the simple text-based output may not convey the full richness of the wave dynamics. Visualizing the displacement as a moving curve would bring the phenomenon to life. Nonetheless, the code demonstrates the core finite difference technique in a compact form.