[aeroiitkgp]

Hello,

I am new to SU2 and might ask very basic questions. I tried reading as much posts as possible, but didn't find what I am looking for.

I am trying to run a transonic case with SA model on an aircraft. The mesh I am using has a Y+ around 1. The problem runs very well when I choose first order for convective as well as turbulence terms. If any of the two I switch to second order, the problem diverges almost instantaneously.

I ran the case with first order for almost 11200 iterations. By this time the force and moment coefficients had converged well. But if I switch to second order for even the convective terms it diverges from the very first iteration.

The multigrid is also turned off. I tried turning on the multigrid without switching to second order, but that too diverged the solution after some 20 iterations.

I have attached the configuration file and the residual convergence plots for first order with and without multigrid. Please help me solve this issue.

Note: The angle of attack in the config file is 20, but the same thing happens for even very low angles of attack.

Other questions with some relation to the problem:

1. After I restart, the iteration number starts from 0 again and it overwrites earlier history.dat file. (For restart I had copied restart_flow.dat to solution_flow.dat). But I can see that the solution has indeed restarted as the coefficients value is from the last iteration.

2. The values of primary variables I get after non-dimesionalization is not 1.
Pressure is 0.9985....., velocity is 1.13, and so on. Please list down the equations used to derive other quantities from Mach Number, Temperature, angles and Reynolds Number if possible. What are the other possible input parameters, like, directly specifying the pressure, temperature and velocity components?

[fpalacios]

Hi, thanks for using SU2,

Have you tried to reduce the CFL number? sometimes it is also important to play with
LIMITER_COEFF= 0.3

I'm not sure what version of SU2 are you using but please note that we have improved the stability of the code in the latest SU2 releases.

Best Regards,

Francisco
SU2 lead developer

[aeroiitkgp]

Thank you Francisco,

I am using the latest build of SU2.

I will try to change the limiter coeff as suggested by you. For CFL, how much should I start with? I am using an implicit formulation. Also should I use CFL ramping?

[fchan]

Hi,

I'm having similar issues with a second order simulation, although mine is unsteady and inviscid. It's subsonic (M=0.5), so I doubt it has anything to do with the turbulent solver.

Mine seems to work fine with first order (haven't tested it more than a couple of hundred time steps). However, I can't get more than 10 time steps done with 2nd order spatial. I have limited the number of iterations per time steps, so it eventually diverges. Additionnally, I'm not using MG nor CFL slopes (CFL=0.5 for now). I'm also using implicit time discretization.

@fpalacios: There are a couple of parameters for the slope limiter in SU2 and I'm having trouble identifying which parameters they relate to Venkatakrishnan's paper (AIAA 93-0880). Could you please elaborate on these:

• REF_ELEM_LENGTH= 0.1 (default)
• LIMITER_COEFF= 0.3 (default)
• LIMITER_ITER= 999999 (default)
• SHARP_EDGES_COEFF= 3.0 (default)
• REF_SHARP_EDGES= 3.0 (default)
• SENS_REMOVE_SHARP= NO (default)

Best regards,
Frederik

[aeroiitkgp]

Hi,

By just changing the limiter coeff to 0.3 from 10.0 earlier, I was able to run with second_order_limiter selected as discritization for both convective and turbulent terms. However, still multigrid does not work. All other settings are as shown in the previous configuration file. Kindly help.

It would be nice if we have all the possible parameters of configuration file available in one sample file. I am not able to find such a file.

Also the non-dimensionalization is not leading to unity for many variables. Please tell me how to solve this.

[Santiago Padron]

The file config_template.cfg in the main directory of the SU2 code contains all the possible parameters. You can find it here https://github.com/su2code/SU2/blob/...g_template.cfg
Also there is an option for different types of non-dimensionalization.
% Flow non-dimensionalization (DIMENSIONAL, FREESTREAM_PRESS_EQ_ONE,
% FREESTREAM_VEL_EQ_MACH, FREESTREAM_VEL_EQ_ONE)
REF_DIMENSIONALIZATION= DIMENSIONAL