Open
Description
Hi,
I'm trying to solve the Swift-Hohenbeg model's PDE as such:
l = 3
k_c = 0.5
grid = pde.CartesianGrid([[0, 100], [0, 100]], [100, 100])
field = pde.ScalarField.random_normal(grid, mean=10)
storage = pde.MemoryStorage()
eq = pde.PDE({"u": f"u -{l}*u**3 -laplace(laplace(u)) -{k_c**4}*u - 2*{k_c**2}*laplace(u)"})
result = eq.solve(field, t_range=100, dt=1e-2, tracker=["progress", storage.tracker(1)])
result.plot()
But for many values of l and k_c, The final result becomes a scalar field filled with NaN.
The current values causes this error.
Thanks
Metadata
Metadata
Assignees
Labels
No labels