mixing_layer_2D

===================== 2D Mixing Layer Test ===================== This test case simulates the Kelvin-Helmholtz instabilty where an initial hyperbolic tangent velocity profile imposed in a fully periodic 2D square box is slightly perturbed to initiate rolling of the shear layers.

The dimension of the side length of the bounding box and the non-dimensional run time are specified through the --grid-size and the --run-time command line arguments respectively. If unspecified default values of 256 and 1.0 are used.

The test conducts a validation run and computes the normalized kinetic energy as a validation metric. Further, the performance of the employed stream-collide algorithm is evaluated and reported in Mega lattice updates per seconds (MLUPS).

The test can be executed run in 2 modes, namely, serial/openmp parallelized.

Serial runs: python mixing_layer_2D.py

OpenMP parallel runs: OMP_NUM_THREADS=4 python mixing_layer_2D.py --openmp