solution convergence problem

anbuselvan

I am simulating the internal flow through a scramjet engine.

Desired analysis type: 2-D steady viscous-laminar

Inlet condition: pressure far field

Pressure: 48Pa

temperature: 67K

mach number: 8

The geometry consists of convergent type inlet followed by constant height combustion chamber and the divergent type nozzle. width is constant throughout. Since the geometry is symmetrical, i tried the simulation with 2-D symmetry condition.

Solver details: coupled, second order upwind discretization implicit scheme.

I meshed with gambit using unstructured (pave scheme) quadrilateral elements. The solution is well converged for inviscid case. the residuals went beyond 10e-12, but the solution didn't converge for viscous case. the X-momentum residuals starts oscillating at 10e-4, then Y-momentum oscillates, then energy and continuity. sometimes it blew up if i continue to simulate.

what could be the possible reasons for it?