\(\triangleright \)Finally, define the main function with a MPI session (for parallel computation). The 2D experiment corresponds to the use of the function NS_steady while iterating through orders and refinement levels. And, the 3D experiment is running the function NS_3D with order 4. No refinements were done because the original mesh (done by MFEM [1]}) gives a good solution already.
