Stability in presence of cubic non-linearity in two dimensions
Initial conditions are chosen as one-dimensional (no x-dependence. Since a one-dimensional soliton does not collapse, the pattern stays stable in y-direction, until instabilities in x-direction break translational symmetry. Since in two dimensions solitons exhibit self-focusing and collapse, the pattern finally diverges.
Equation to solve
+ ∆Ψ(x,y) + (|Ψ(x,y)|2)Ψ(x,y) = 0.
Square 96×96 (a two-dimensional area with coordinates (x,y)).
Mesh grading is 512×512×1 was obtained with use of blockMesh utility.
The function at zero time is Gaussian distribution along y direction, with small modulations along x direction,
Ψ(x,y,z) = e−[(y2)/4] [1+ 0.0001sin(πx/4)].
Neumann (zero gradient normal to the boundary).
The time step is ∆t = 0.005. Calculation time is around 2 hours.
Dynamics of cubic non-linear Schrödinger equation over a 2D square domain.