I have this question for quite a long time. I want to know more about how the code deal with non-dimensional values. The non-dimensional compressible NS should have a Reynolds number appear in the viscous term. But the compressible solver only asks for viscosity as an input. So how do you calculate other reference values that Reynolds number needs. And viscosity should be a function of temperature as well. So my guessing is everything is based on initial condition I set?

Thanks for reply. I still have question about it. So PyFR doesn’t solve non-dimensional equations, it treats any input as dimensional? As long as I can write done my relationship to get desired Mach number and Reynolds number, it will take those input variables rho,u,v,w,p as reference from initial condition? Am I right in this way?

PyFR doesn’t care about dimensions at all. So it neither uses reference density, temperature, and pressure to non-dimensionalise things nor does it output results in set dimensions. It is up to the user to interpret the output as either m/s, inches/hours or light years/sol.