CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   SU2 (https://www.cfd-online.com/Forums/su2/)
-   -   Possible bug in SU2_DOT with DISCRETE_ADJOINT (https://www.cfd-online.com/Forums/su2/169587-possible-bug-su2_dot-discrete_adjoint.html)

kmandar April 13, 2016 12:25

Possible bug in SU2_DOT with DISCRETE_ADJOINT
 
Hi SU2 developers,

Here is what I have

1. Run SU2_CFD with the quick start NACA0012 example/
2. Run SU2_CFD_AD with MATH_PROBLEM= DISCRETE_ADJOINT
3. Run SU2_DOT for Hicks-Henne design variable -- I get the following error message in the log file (or at the command prompt):
Code:

There is no adjoint restart file!! solution_adj_cl.dat.
BUT if I repeat the last step:

3. Run SU2_DOT with configuration file option MATH_PROBLEM= CONTINUOUS_ADJOINT, the job runs correctly and gives me the the result for the derivative of CL with the design variable.

It seems there is some bug with DISCRETE_ADJOINT option for MATH_PROBLEM if run with SU2_DOT.

-- Mandar

talbring April 13, 2016 13:31

Hi kmandar,

you have to use SU2_DOT_AD to get the correct projection. But this requires the volume adjoint solution, this is why you need to copy restart_adj_cl.dat to solution_adj_cl.dat first.

Or just run the script discrete_adjoint.py, this performs all necessay steps automatically.

Tim

kmandar April 13, 2016 17:32

Hi talbring,

That worked well.

The script runs well too.

Thanks,
Mandar

kmandar April 13, 2016 17:52

After running the discrete_adjoint.py script, the gradient for lift are stored in the file of_grad_cl.dat and it looks like this:


Code:

VARIABLES="VARIABLE"      , "GRADIENT"      , "FINDIFF_STEP" 
          0              , 2.83854        , 0.001

What is the meaning of the last column (FinDiff_Step) in the context of discrete adjoint method?

Thanks,
Mandar

talbring April 15, 2016 04:17

Hi Mandar,

actually this is a leftover of an earlier version where the projection was also in the case of the discrete adjoint done using Finite Differences.

However, its not done anymore so you can ignore this column.

Tim

kmandar April 15, 2016 16:27

Hi Tim,

Ok, thanks,
Mandar


All times are GMT -4. The time now is 03:21.