Note
Click here to download the full example code
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\]
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:23<6:29:45, 2340.84s/it]
0%| | 0.02/10.0 [00:23<3:14:41, 1170.45s/it]
1%|1 | 0.1/10.0 [00:23<38:37, 234.09s/it]
42%|####2 | 4.2/10.0 [00:23<00:32, 5.58s/it]
42%|####2 | 4.2/10.0 [00:23<00:32, 5.58s/it]
100%|##########| 10.0/10.0 [00:23<00:00, 2.34s/it]
100%|##########| 10.0/10.0 [00:23<00:00, 2.34s/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 23.648 seconds)