Hello all,
I've been looking into including a heat flux term in the Euler equations, along the lines of the usual heat flux that appears in the energy equation in the Navier-Stokes equations (but no viscous stresses).
Is something like this already included in some part of RAMSES? My understanding is that including this in the Godunov solver is very challenging - the Riemann problem becomes much harder to solve I think. So I've been looking into adding this term as another source term, in a similar manner to what is done already in the routine "set_uold" (godunov_fine.f90) with "add_pdv_source_terms". I suspect, however, that ensuring numerical stability in the solver will be difficult, and my preliminary tests support this.
Any help would be much appreciated, thanks!
Cheers,
Graeme