dealiasing and pressure Hessian
TurTLE allows for two dealiasing methods now (2/3 and smooth dealiasing). We need to also implement 1/2 dealiasing. Opening issue for easier tracking of development.
See below document prepared by @mcarbone :
There are two points laid out in the above document:
-
We may test the new dealiasing method by simply computing spectra of dealiased fields and comparing spectrum for 2/3 with spectrum for 1/2. There are analytic expressions for wavenumbers where spectra go to 0.
-
Pressure Hessian and squared velocity gradients satisfy known equations that fail when using 2/3 dealiasing, but work for 1/2.
Therefore I see the following steps for the work:
- Implement 1/2 dealiasing (i.e. new template for
kspace). - Write new TurTLE test that generates spectra comparisons for different dealiasing methods, i.e. point 1 from above.
- Write post-processing tool that computes the pressure Hessian and the velocity gradient, and verifies the equations detailed in 2.
Afterwards we see exactly what type of postprocessing tool we write for future data analysis.