[MaxGains]

Hello,

I'm looking to use SU2 to manually compute derivatives of an objective function
with respect to flight conditions (Mach, AoA) and some shape design variables.
I'm not looking to perform shape optimization here; I only need to evaluate
gradients at a range of input points.
I've parameterized my geometry with a 3D FFD box, which I am manually
adjusting by setting the DV_PARAM and DV_VALUE entries in my config file.
The following steps are running correctly:
+ mesh deformation using SU2_DEF
+ evaluation of primal solution using SU2_CFD
+ evaluation of adjoint solution using SU2_CFD.
Using the SENS_MACH and SENS_AOA outputs, I've validated my adjoint solution
against finite difference gradients.

When I run SU2_DOT <config_file>.cfg, I get the following error message:
Surface sensitivity not found in file
Here, my simulation is 3D, so my surface adjoint results are output as
`surface_adjoint.vtu'; could this be the cause of the problem?

If I run the same steps on the QuickStart example, everything works and I get a
gradient in the `of_grad.dat' file.

Thanks for any help.

[Tongtong1919]

Hi MaxGains,

Could you please describe the problem in more detail or provide a screenshot of the error?

From my experience, the SU2_DOT module uses the result file from SU2_CFD (restart_adj.dat, which should be renamed to solution_adj.dat) to obtain the projection results for the design variables.

Tongtong

[MaxGains]

Thanks Tongtong.


As a reference point, the QuickStart examples works fine for me.
After running the primal and adjoint simulations and renaming the restart files, my directory contains the following files:

adjoint.vtu  flow.vtu  history.csv  inv_NACA0012.cfg  mesh_NACA0012_inv.su2  solution_adj_cd.dat  solution_flow.dat  surface_adjoint.csv  surface_flow.csv

I then run

SU2_DOT inv_NACA0012.cfg

and successfully generate a file called

of_grad.dat

I'm now trying to replicate this workflow for my own example.
I first converge my primal and adjoint simulations, then rename restart_adj_cl.dat to solution_adj_cl.dat.
I run my job file

SU2_DOT run_adjoint.cfg

and get the following error output:

----------------- Physical Case Definition ( Zone 0 ) -------------------
Input mesh file name: mesh_deformed.su2
Input sensitivity file name: surface_adjoint.

-------------- Surface deformation parameters ( Zone 0 ) ----------------
Design variables definition (markers <-> value <-> param):
FFD (control point) <-> WING3D <-> 0 <-> ( WING, 0, 0, 0, 1, 0, 0 )
FFD (control point) <-> WING3D <-> 0 <-> ( WING, 1, 0, 0, 1, 0, 0 )
FFD (control point) <-> WING3D <-> 0.0406 <-> ( WING, 0, 1, 0, 1, 0, 0 )
FFD (control point) <-> WING3D <-> 0.0406 <-> ( WING, 1, 1, 0, 1, 0, 0 )
FFD (control point) <-> WING3D <-> 0 <-> ( WING, 0, 0, 1, 1, 0, 0 )
FFD (control point) <-> WING3D <-> 0 <-> ( WING, 1, 0, 1, 1, 0, 0 )
FFD (control point) <-> WING3D <-> 0.0406 <-> ( WING, 0, 1, 1, 1, 0, 0 )
FFD (control point) <-> WING3D <-> 0.0406 <-> ( WING, 1, 1, 1, 1, 0, 0 )

---------------- Design problem definition  ( Zone 0 ) ------------------
CL objective function.

-------------------- Output Information ( Zone 0 ) ----------------------
Output gradient file name: of_grad.dat.

------------- Config File Boundary Information ( Zone 0 ) ---------------
+-----------------------------------------------------------------------+
|                        Marker Type|                        Marker Name|
+-----------------------------------------------------------------------+
|                          Far-field|                           FARFIELD|
+-----------------------------------------------------------------------+
|                     Symmetry plane|                           SYMMETRY|
+-----------------------------------------------------------------------+
|                     Heat flux wall|                             WING3D|
+-----------------------------------------------------------------------+
Three dimensional problem.
2057 grid points.
1920 volume elements.
3 surface markers.
120 boundary elements in index 0 (Marker = WING3D).
192 boundary elements in index 1 (Marker = SYMMETRY).
120 boundary elements in index 2 (Marker = FARFIELD).
1728 hexahedra.
192 prisms.
Identify vertices.
Storing a mapping from global to local point index.

----------------------- Preprocessing computations ----------------------
Setting local point connectivity.
Checking the numerical grid orientation of the elements.
All volume elements are correctly orientend.
All surface elements are correctly orientend.
Identify edges and vertices.
Setting the bound control volume structure.
Storing a mapping from global to local point index.
Reading surface sensitivities at each node from file.


Error in "virtual void CPhysicalGeometry::SetBoundSensitivity(CConfig*)":
-------------------------------------------------------------------------
Surface sensitivity not found in file.
------------------------------ Error Exit -------------------------------

Here, since my case is 3D, my surface adjoint file is surface_adjoint.vtu (rather than surface_adjoint.dat).
Could this be causing the problem?


Regards,
Viv

[bigfootedrockmidget]

I don't think the reader can read paraview vtu files. You need SURFACE_CSV in OUTPUT_FILES I guess.
If this is indeed the problem, maybe we can change the error message.

[MaxGains]

Yep -- that was the problem!


Very easy to convert between the two with vtk library though