[Improvement] Extend flux solver to higher order

[jaykalinani]

The implemented HLLE flux solver is only 3rd order accurate. We should extend it to higher order, as mentioned, for instance, in the Appendix A.2 of the second Spritz paper: https://arxiv.org/pdf/2012.10174.pdf