[Blooper]

Hello!

I went through the tutorials, and now I decided to try out an example of my own. I'm simulating a flow inside a tubular channel with a smooth expansion. The precisa data is as follows (I'm using a 3d mesh generated with Blender+enGrid):

Inlet Diameter: 0.2 m
Outlet Diameter: 0.3 m
Expansion Length: 0.3 m (and height = 0.1 m, outlet - inlet diameters)

Pretty simple geometry. Now all I want is to make some air enter this channel at a given speed (22.3 m/s, to be precise), and analyze the pressure loss on the expansion part. So I decided to set the inlet and outlet walls with an euler marker (MARKER_EULER= ( inletwall, outletwall ) and the expansion wall with a heaflux marker, with 0 heat exchange (MARKER_HEATFLUX= ( expansionwall, 0.0 ). I also set the inlet type to be mass-flow, so:
INLET_TYPE= ( MASS_FLOW )
and then
MARKER_INLET= ( inlet, 1.185, 22.3, 1.0, 0.0, 0.0 )

Now, for the outlet condition, i'm obliged to set a back pressure, and so I leave it at the standard 101325. And adding the other necessary definitions, the computation diverges. I've been trying to play around with a couple changes, but no success so far.

I'll post the first part of the program (with constants and boundary definitions, which is where I think the error lies), and if you can spot anything weird, please feel free to point out, it would be of great help! As a complete starter in the CFD business, I'm pretty sure there could be some dumb mistakes in there, so don't feel bad hammering my mistakes, I'm willing to take any constructive criticism.

PHYSICAL_PROBLEM= NAVIER_STOKES

KIND_TURB_MODEL= SA

MATH_PROBLEM= DIRECT

RESTART_SOL= NO

MACH_NUMBER= 0.06408

AoA= 0.0

SIDESLIP_ANGLE= 0.0

FREESTREAM_TEMPERATURE= 300.0

REYNOLDS_NUMBER= 285000.0

REYNOLDS_LENGTH= 0.2

GAMMA_VALUE= 1.4

GAS_CONSTANT= 287.87

PRANDTL_LAM= 0.72

PRANDTL_TURB= 0.9

CONVERT_TO_METER= 1.0

WRITE_CONVERTED_MESH = NO

REF_ORIGIN_MOMENT_X = 0.25

REF_ORIGIN_MOMENT_Y = 0.00

REF_ORIGIN_MOMENT_Z = 0.00

REF_LENGTH_MOMENT= 1.0

REF_AREA= 1.0

REF_PRESSURE= 1.0

REF_TEMPERATURE= 1.0

REF_DENSITY= 1.0

REF_ELEM_LENGTH= 0.1

MARKER_EULER= ( inletwall, outletwall )

MARKER_HEATFLUX= ( expansionwall, 0.0 )

INLET_TYPE= ( MASS_FLOW )

MARKER_INLET= ( inlet, 1.185, 22.3, 1.0, 0.0, 0.0 )

MARKER_OUTLET= ( outlet, 101325 )

MARKER_PLOTTING= ( outlet )

MARKER_MONITORING= ( expansionwall )