I am trying to re-do the simulation of Supersonic Flow in a 2D COnvergent-Divergent Channel provided as validation case in Flow Simulation 2013 Technical Reference.
I used the geometry and the settings provided in the technical reference (inlet boundary condition: M=3, Tstatic=293.2K, Pstatic=1atm, ideal wall, laminar only). Since the boundary condition for the outlet is not given, I simply assumed it to have Pstatic=10Pa. When I run my simulation, I found that my solutions is fluacting for a long time and the mass flow rate never get converse. Can anyone tell me what's going wrong? Is it because the assumption of Pstatic=10Pa is invalid? If yes, what is the appropriate setting for outlet boundary?