C++ CFD. the initial interest is in hyperbolic conservation laws in 1D, 2D, and 3D, simple structured meshes to begin.
I rolled my own arrays for this. Probably want to update with Eigen, but allow for hybridized approaches as we will go to the GPU.
Here is shocked flow entering a box and reflecting off the sides:
- van Albada limiter
- Roe flux
Here is the Shock Tube solution computed with the 1D partition of this code:
- minmod limiter
- Roe flux