I am recently interested in the properties of Bose-Einstein condensation in two-dimension hydrodynamic turbulence. Without the Ekman term (linear friction), the time to reach a saturated state is proportional to the Reynolds number, which can be very long. Therefore, I want to improve my GPU code to use adaptive time stepping.
I spent couple weeks implementing a simple 2D hydrodynamic code in CUDA C, which is a technology developed by nVidia. The code runs on the GPU instead of CPU. It uses Fourier spectral method with the 2/3 dealiasing rule and 4-th order Runge-Kutta integrator to solve the vorticity equation