
[Sponsors] 
Sensitivity in surface_sens.vtk V.S. Surface_Sensitivity in surface_adjoint.vtk 

LinkBack  Thread Tools  Search this Thread  Display Modes 
December 19, 2018, 04:18 
Sensitivity in surface_sens.vtk V.S. Surface_Sensitivity in surface_adjoint.vtk

#1 
Member
Join Date: May 2017
Posts: 39
Rep Power: 8 
Hi everyone,
I am using the discrete adjoint solver and I am interested in surface sensitivities in the normal direction. However, I found that the "Surface_Sensitivity" inside "surface_adjoint.vtk" is not the same as the "Sensitivity" in "surface_sens.vtk". The Surface_Sensitivity is actually by 5 orders of magnitude larger than the Sensitivity. I thought they're both computed by taking the inner product of the surface mesh sensitivities and the cell normal vector, but apparently I was wrong. I would be glad if someone could explain how they are calculated. Many thanks in advance 

December 29, 2018, 11:33 

#2 
Member
na
Join Date: Jul 2018
Posts: 78
Rep Power: 6 
Hi Ryan,
Upfront: I use a recent develop version of SU2 in which 'Surface_Sensitivity' is used for the variable in surface_sens.* (* is vtk for your paraview output) and surface_adjoint.*, I try to be clear with what I mean. I'll try to explain briefly what I found in the code and what is happening here in general: 1. The surface_adjoint.* is written by SU2_CFD_AD (just as adjoint.*) and surface_sens.* by SU2_DOT_AD. 2. The sensitivity information in both files "mean different things". In adjoint.* and surface_adjoint.* you'll find the "Volume Sensitivities", i.e. the sensitivity with respect to every(!) mesh node coordinate alone. In surface_sens.* the information of the mesh deformation algorithm is incorporated and they are only available on the surface nodes (which control the geometry/shape). 3. The computation of Surface_Sensitivity is, as you said, the dot product between normal and sensitivity vector. But I am not sure whether that computation is correct for surface_adjoint.* and therefore also adjoint.*. I have to take a closer look to be sure though, and one has to be careful whether one has a unit normal or Area normal. For now I think a factor of 1/Area is the difference to the computation in surface_sens (based on a given sensitivity vector). Summing up: The sensitivities in both files are different in what they represent and shouldn't be the same (in general). However, the computation of Surface_Sensitivity looks inconsistent to me for both files which might be an error. I hope what I wrote was somewhat understandable. Let me know if an explanation was bad. I'll take a closer look at the computation in the next days and report here. If you want to look in the code yourself and try out things I can give you the lines of code where to computation takes place. Cheers, Tobi 

January 4, 2019, 16:05 

#3 
Member
Join Date: May 2017
Posts: 39
Rep Power: 8 
Hi Tobi,
thanks a bunch for the thorough explanation. I believe you're right, just by comparing the result in Paraview it seems the difference between them comes from the fact that one of them uses the area normal vector and the other uses the unit normal vector. Cheers and happy new year, Ryan 

July 27, 2019, 02:25 
Sensitivity analysis failure

#4 
New Member
Angu Praveen
Join Date: May 2019
Posts: 14
Rep Power: 6 
Dear all,
I'm trying to do a 3D SE2A Aircraft wing optimization using SU2. Unfortunately, I am unable to get the following things working. 1. The solution doesn't converge even after 5000 iterations both in ADJOINT and DIRECT solvers. I mean it keeps computing until the max iter value. 2. After the iterations are done, I tried to compute the sensitivity using SU2_DOT_AD and it gives of_grad.dat file with '0' as the value with the following message before "Design variable (FFD_SETTING) number 0." Any help would be highly appreciated Thanks in advance Angu 

September 6, 2020, 07:03 

#5  
New Member
Haoge Li
Join Date: Oct 2017
Posts: 2
Rep Power: 0 
Hi TKatt,
Thanks for your thorough explanation, as I computed the naca0012 airfoil with the continuous adjoint solver, and i'm confused about the Surface_sensitivity in the surface_adjoint.dat. Does the surface area normal vector point to the flow field, or the interior of the body? Any help and explanation will be appreciated! Thank you a lot Haoge Li Quote:


March 2, 2022, 18:29 
Is there a good example for surface sensitivity using SU2

#6  
New Member
California
Join Date: Jan 2019
Posts: 13
Rep Power: 6 
Hi Ryan,
I want to compute the surface sensitivity for my configuration as well. But it seems to me that SU2 surface sensitivity requires some deformation by FFD or by user defined deformation function. What I want is small perturbance on surface and I want to calculate the current surface sensitivity. I wonder how to do that? Quote:


March 5, 2022, 06:36 

#7  
Member
Join Date: Jun 2020
Posts: 30
Rep Power: 4 
Quote:


Tags 
discrete adjoint, su2, su2 ad, su2 adjoint 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
How to implement a sensitivity analysis for certain items  Mason liu  CFX  10  August 17, 2017 08:46 
Estimating lift change using sensitivity  arrow201  SU2  0  January 31, 2017 13:28 
About Surface Sensitivity computation  Tommy Chen  SU2  1  March 10, 2014 16:14 
Adjoint sensitivity on specified surface  robyTKD  OpenFOAM Programming & Development  0  April 29, 2013 05:05 
What exactly is a sensitivity analysis?  Josh  CFX  3  August 19, 2009 09:18 