I need to write a DNS code for simulating reacting flows in turbulent conditions. The code has to be highly scalable, because the computational cost of simulations is expected very large.My idea is to use C++ and I was wondering if you have suggestions about a very efficient C++ library to manage vector/matrix operations in distributed memory architectures. As an example, I know that PETSc and Trilinos packages are widely used, but I have no idea about the best choice for my case.Additional details:compressible formulation structured 3D grids finite difference discretization (staggered arrangement) high order FDgood scalability on 400-500 coresonly for use in my group (no need of complex object oriented design)Thanks in advance!
Direct Numerical Simulation of reacting flows
fluid dynamics;hpc;vector