This is a highly efficient and elegant code, but I suspect there might be an issue with the equation. After reviewing the source code, it appears that the code solves the dimensionless Navier-Stokes equations (I noticed that the kinematic viscosity ν is replaced with 1/Re). However, in the dimensionless form, the external force term should also be scaled by 1/Re .
After correcting the equation, I observed that the divergence issues have significantly improved, which indicates the revised formulation now better aligns with physical principles.
Reference Support:
For the correct dimensionless form of the Kolmogorov flow equation, see:
Kolmogorov flow: Linear stability and energy transfers in a minimal low-dimensional model
Chaos 30, 073110 (2020); doi: 10.1063/5.0002751
This is a highly efficient and elegant code, but I suspect there might be an issue with the equation. After reviewing the source code, it appears that the code solves the dimensionless Navier-Stokes equations (I noticed that the kinematic viscosity ν is replaced with 1/Re). However, in the dimensionless form, the external force term should also be scaled by 1/Re .
After correcting the equation, I observed that the divergence issues have significantly improved, which indicates the revised formulation now better aligns with physical principles.
Reference Support:
For the correct dimensionless form of the Kolmogorov flow equation, see:
Kolmogorov flow: Linear stability and energy transfers in a minimal low-dimensional model
Chaos 30, 073110 (2020); doi: 10.1063/5.0002751