2.2.7 Klein-Gordon equation

This example solves the Klein-Gordon equation in one dimension, showing the effect of the mass term on wave propagation. With mass=0 the equation reduces to the standard wave equation; increasing the mass introduces dispersion.

  0%|          | 0/50.0 [00:00<?, ?it/s]
Initializing:   0%|          | 0/50.0 [00:00<?, ?it/s]
  0%|          | 0/50.0 [00:02<?, ?it/s]
  0%|          | 0.01/50.0 [00:02<4:01:07, 289.40s/it]
  1%|          | 0.38/50.0 [00:02<06:17,  7.62s/it]
 59%|█████▉    | 29.65/50.0 [00:02<00:01, 10.24it/s]
 59%|█████▉    | 29.65/50.0 [00:02<00:01, 10.23it/s]
100%|██████████| 50.0/50.0 [00:02<00:00, 17.25it/s]
100%|██████████| 50.0/50.0 [00:02<00:00, 17.25it/s]

from pde import KleinGordonPDE, ScalarField, UnitGrid

grid = UnitGrid([128])  # generate grid
u = ScalarField.from_expression(grid, "exp(-((x - 32) / 5) ** 2)")

eq = KleinGordonPDE(speed=1, mass=0.5)  # define the pde
state = eq.get_initial_condition(u)
result = eq.solve(state, t_range=50, dt=0.01)
result[0].plot()