On Numerical Methods for Aeothermodynamic Applications