2.6 Kuramoto-Sivashinsky - Using PDE class

This example implements a scalar PDE using the PDE. We here consider the Kuramoto–Sivashinsky equation, which for instance describes the dynamics of flame fronts:

\[\partial_t u = -\frac12 |\nabla u|^2 - \nabla^2 u - \nabla^4 u\]
pde custom expression
  0%|          | 0/10.0 [00:00<?, ?it/s]
Initializing:   0%|          | 0/10.0 [00:00<?, ?it/s]
  0%|          | 0/10.0 [00:13<?, ?it/s]
  0%|          | 0.01/10.0 [00:24<6:41:01, 2408.55s/it]
  0%|          | 0.02/10.0 [00:24<3:20:18, 1204.29s/it]
  1%|          | 0.12/10.0 [00:24<33:03, 200.72s/it]
 54%|█████▍    | 5.41/10.0 [00:24<00:20,  4.45s/it]
 54%|█████▍    | 5.41/10.0 [00:24<00:20,  4.46s/it]
100%|██████████| 10.0/10.0 [00:24<00:00,  2.41s/it]
100%|██████████| 10.0/10.0 [00:24<00:00,  2.41s/it]

from pde import PDE, ScalarField, UnitGrid

grid = UnitGrid([32, 32])  # generate grid
state = ScalarField.random_uniform(grid)  # generate initial condition

eq = PDE({"u": "-gradient_squared(u) / 2 - laplace(u + laplace(u))"})  # define the pde
result = eq.solve(state, t_range=10, dt=0.01)
result.plot()

Total running time of the script: (0 minutes 24.276 seconds)