Case 7: (Convergence Study) B(z) = 0, A = 2 * π * α for
α = 1.d3. Here we study the convergence of the incompressibility condition
by plotting |det(∇ x) - 1| with a range of grid resolutions.
- Part I : nx = ny = 129
- Part II : nx = ny = 257
- Part III : nx = ny = 513
- Part IV : nx = ny = 1025
- Part V : nx = ny = 2049
- Convergence Study: Here we use the loglog plot of
the maximum error = |det(∇ x) - 1| vs the spatial step size h =
400/(nx-1), for nx=129, 257, 513, 1025, 2049 in order to estimate the
order of the convergence for each integration time T = 3, 4.5 and 6
days.