[vachan]

Hello,

I am using SU2 to simulate a Mach 2 steady laminar shock-boundary layer interaction problem. To generate the shock, I have imposed pre- and post-shock conditions (using MARKER_SUPERSONIC_INLET) on two parts of the top boundary. To prevent the reflection of leading edge shock, for the third part of the top boundary, I have given the MARKER_SUPERSONIC_OUTLET condition. I am listing the config file here.

SOLVER = NAVIER_STOKES

VISCOSITY_MODEL= SUTHERLAND

REYNOLDS_NUMBER = 1E5
REYNOLDS_LENGTH = 80E-3
MACH_NUMBER = 2.15
FREESTREAM_OPTION = TEMPERATURE_FS
FREESTREAM_PRESSURE = 1117.3
FREESTREAM_TEMPERATURE = 157
FREESTREAM_VELOCITY = (540, 0, 0)

MESH_FILENAME = mesh2_62x53.su2
MESH_FORMAT = SU2
MARKER_EULER = (front_back)
MARKER_HEATFLUX = (wall, 0.0)
MARKER_SUPERSONIC_INLET = (preshock_inflow, 157, 1117.3, 540, 0, 0, postshock_inflow, 167.284, 1393.6, 520.562, -34.6182, 0)
MARKER_SUPERSONIC_OUTLET = (outflow_top)
MARKER_OUTLET = (outflow_side, 1718.81)

CONV_NUM_METHOD_FLOW = SLAU2

INIT_OPTION = REYNOLDS
TIME_DOMAIN = YES
TIME_MARCHING = TIME_STEPPING
TIME_DISCRE_FLOW = RUNGE-KUTTA_EXPLICIT
UNST_CFL_NUMBER = 0.5
MAX_TIME = 8e-3
TIME_ITER = 1000000000

OUTPUT_WRT_FREQ = 10000
SCREEN_WRT_FREQ_TIME = 1000
HISTORY_WRT_FREQ_TIME = 100
SCREEN_OUTPUT = (CUR_TIME, TIME_STEP, TIME_ITER, RMS_DENSITY)
HISTORY_OUTPUT = (CUR_TIME, TIME_ITER, RMS_RES)

However, I am still seeing the reflection of the leading edge shock from the top boundary (attached figure temp.png). It reflects onto the bottom boundary (a wall) and therefore disturbs the flow there, thereby hampering the wall pressure comparison with reference data.

Is there anything I can do to prevent the reflection of leading edge shock? Please let me know in case any additional information is required.

Thanking in anticipation,
Vachan

[pcg]

If I understand correctly, I would try to use farfield for the second and third parts combined

[vachan]

Dear Pedro,

Thank you very much for responding. If I understand correctly, the farfield BC treats the boundary to have free stream conditions? So wouldn't that be wrong here since the flow conditions are not free stream (pre-shock), but post-shock values?

Please correct me if I am wrong.

Vachan

[pcg]

Yes but you can change the freestream conditions to be what you want post shock.

[vachan]

Never thought of that! Thank you very much. I will get back with an update.