CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > SU2 > SU2 Shape Design

Adjoin Output Options

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree2Likes
  • 1 Post By pcg
  • 1 Post By pcg

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 10, 2020, 17:38
Default Adjoin Output Options
  #1
New Member
 
Eduardo de Oliveira Carvalho
Join Date: Mar 2018
Location: São José dos Campos - SP. Brazil
Posts: 15
Rep Power: 5
Eduardo Carvalho is on a distinguished road
Hi all,

What are the options of the adjoint solver's output?

I failed in executing a dryrun to verify the options of output available for the adjoint solver. I also didn't find this in the user's guide or in the template.

Thanks
Eduardo Carvalho is offline   Reply With Quote

Old   August 10, 2020, 19:16
Default
  #2
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 301
Rep Power: 7
pcg is on a distinguished road
Hello,

There are no options, you get the coordinates, the solution (the adjoints of the primal solution variables), and the sensitivities (volume and surface).
In the case of the continuous adjoint only the surface sensitivities matter, for discrete adjoints only the volume sensitivities matter (the surface values are just for post processing and they do not reflect what would happen if the mesh was deformed).

There is not much else we could output for adjoint solvers.

Cheers.
Eduardo Carvalho likes this.
pcg is offline   Reply With Quote

Old   August 10, 2020, 19:42
Default
  #3
New Member
 
Eduardo de Oliveira Carvalho
Join Date: Mar 2018
Location: São José dos Campos - SP. Brazil
Posts: 15
Rep Power: 5
Eduardo Carvalho is on a distinguished road
Thank you for your answer, Pedro.

I'm trying to reproduce the plot from Figure 2 of the unsteady optimization's tutorial (https://su2code.github.io/tutorials/..._Opt_NACA0012/).

I'm trying to find the "drag sensitivity" output, but the original tutorial's adjoin history file doesn't have any variable with this name. I've thought that maybe was some variable not listed in the history output command (maybe the SENS_GEO for instance).

Do you know how to get this output?
Eduardo Carvalho is offline   Reply With Quote

Old   August 11, 2020, 04:35
Default
  #4
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 301
Rep Power: 7
pcg is on a distinguished road
Ah ok, that's different (I thought you were talking about volume outputs).
History outputs are specified in groups, those sensitivities are in group SENSITIVITY (and they should be part of the defaults).
pcg is offline   Reply With Quote

Old   August 11, 2020, 14:03
Default
  #5
New Member
 
Eduardo de Oliveira Carvalho
Join Date: Mar 2018
Location: São José dos Campos - SP. Brazil
Posts: 15
Rep Power: 5
Eduardo Carvalho is on a distinguished road
Ok, so just to clarify: which output variable is being plotted as the drag sensitivity in the Figure 2 of the Unstady Shape Optimization Tutorial???
(https://su2code.github.io/tutorials/..._Opt_NACA0012/)
Eduardo Carvalho is offline   Reply With Quote

Old   August 11, 2020, 17:42
Default
  #6
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 301
Rep Power: 7
pcg is on a distinguished road
I'm afraid that is yet a different thing, the author is looking at drag sensitivity w.r.t. the design variable setup in the config file, this is not a direct output of SU2 as some post processing is required.
Derivatives w.r.t. shape parameters are output to a file named according to GRAD_OBJFUNC_FILENAME (default is of_grad.dat) when running SU2_DOT_AD.
When you run the tutorial do you see that file?
The author might have post processed all time steps manually to build that plot, in practice (i.e. for optimization) you do not need the history of the sensitivity, just the averaged value.
Eduardo Carvalho likes this.
pcg is offline   Reply With Quote

Old   August 12, 2020, 13:57
Default
  #7
New Member
 
Eduardo de Oliveira Carvalho
Join Date: Mar 2018
Location: São José dos Campos - SP. Brazil
Posts: 15
Rep Power: 5
Eduardo Carvalho is on a distinguished road
Yes, I've found the file. It cointains a value of "GRADIENT" for each "VARIABLE", as you mentioned.

So, to summarize:
1 - The SENS_GEO described in the history file is a vector containing a projection of the volumetric sensitivity into the surface for each time step.
2 - The gradient described in the of_grad_cd.csv file is calculated once per design and is a "windowed mean" of the volumetric sensitivities multiplied by the sensitivities w.r.t. your deformation method.
Is that correct?

Thanks again for your help.
Eduardo Carvalho is offline   Reply With Quote

Old   August 13, 2020, 05:44
Default
  #8
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 301
Rep Power: 7
pcg is on a distinguished road
SENS_GEO is the sum of the sensitivities at the surface, it does not include the projection step required to compute the values that go into of_grad_cd.csv.
The purpose of SENS_GEO is just to monitor convergence.
The projection step is, like you say, the sensitivities w.r.t. the deformation method, or more formally the partial derivatives of the mesh coordinates w.r.t. the design parameters.
pcg is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
OpenFOAM - Control and Editing of Lagrangian Output Variables gmag OpenFOAM Programming & Development 1 June 25, 2019 12:29
[ICEM] Appropriate bocos for POPINDA output FlorisvdBeek ANSYS Meshing & Geometry 1 November 19, 2017 04:22
Projected area to output parameters Daniil85 Fluent UDF and Scheme Programming 1 March 16, 2016 15:44
Generating output data on AQWA Naut afitriadhy ANSYS 0 July 2, 2015 02:37
lift and drag on ship superstructures vaina74 OpenFOAM Running, Solving & CFD 3 June 8, 2010 12:30


All times are GMT -4. The time now is 13:53.