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

Hicks-Henne Geometric Sensitivities Matrix

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

Like Tree1Likes
  • 1 Post By kmandar

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 26, 2016, 19:39
Default Hicks-Henne Geometric Sensitivities Matrix
  #1
New Member
 
Outside US
Join Date: May 2013
Posts: 6
Rep Power: 12
elio_xh is on a distinguished road
Is there any way to extract the geometric sensitivities (design velocities) matrix for a Hicks-Henne 2D airfoil adjoint optimization?

I see in history output of a continuous adjoint run the "Sens_Geo" column but my understanding is this value is the summation of all sensitivities on the surface, but what I'm looking for is the full matrix that measures the influence of the design variables on the position of each grid point on the surface mesh.

I've tried editing the CAdjEulerSolver::Inviscid_Sensitivity module (in solution_adjoint_mean.cpp), to print all the CSensitivity values calculated at every point on the boundary, but my understanding is even this already sums each design variable contribution. Is there any way to extract that matrix?
elio_xh is offline   Reply With Quote

Old   March 28, 2016, 17:25
Default
  #2
hlk
Senior Member
 
Heather Kline
Join Date: Jun 2013
Posts: 309
Rep Power: 13
hlk is on a distinguished road
Quote:
Originally Posted by elio_xh View Post
Is there any way to extract the geometric sensitivities (design velocities) matrix for a Hicks-Henne 2D airfoil adjoint optimization?

I see in history output of a continuous adjoint run the "Sens_Geo" column but my understanding is this value is the summation of all sensitivities on the surface, but what I'm looking for is the full matrix that measures the influence of the design variables on the position of each grid point on the surface mesh.

I've tried editing the CAdjEulerSolver::Inviscid_Sensitivity module (in solution_adjoint_mean.cpp), to print all the CSensitivity values calculated at every point on the boundary, but my understanding is even this already sums each design variable contribution. Is there any way to extract that matrix?
Thanks for your question.
The output file surface_adjoint.dat contains the surface sensitivity on your plotted marker, as well as the values of the adjoint variables.
The output file of_grad_{} contains the sensitivity of the objective with respect to the design variables defined in the config file.
The routines in SU2_DOT project the surface sensitivites onto particular design variables.
hlk is offline   Reply With Quote

Old   April 1, 2016, 11:16
Default
  #3
Member
 
Mandar Kulkarni
Join Date: Nov 2013
Location: Virginia Tech, Blacksburg, VA
Posts: 52
Rep Power: 12
kmandar is on a distinguished road
Hi elio_xh,

One way I think of getting discrete design velocities from SU2 would be using finite difference as follows:

1. Use the Hicks_Henne design variable definition to perturb the mesh with a small step, say DV_VALUE= 0.01. You can output the deformed mesh using SU2_DEF suite.

2. Use finite difference to get the design velocity at each grid point using the two meshes and the value of the perturbation.

Of course, the result will depend on the step size, so you may have to do a step size study to get the correct step size.

I believe there would be a way to get these values without using finite difference, but I am not aware of it.

Another way would be to get the design velocities from software you use to discretize the domain using this geometry.

Hope this helps,
Mandar
mlooti likes this.
kmandar is offline   Reply With Quote

Reply

Tags
adjoint, geometric, sensitivities

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
Using multiple Hicks Henne deformations jack1980 SU2 Shape Design 3 July 30, 2019 15:01
Convection Diffusion 1-D Matrix Problem balrog Main CFD Forum 1 November 3, 2014 13:34
Force can not converge colopolo CFX 13 October 4, 2011 22:03
OpenFOAM version 1.6 details lakeat OpenFOAM Running, Solving & CFD 42 August 26, 2009 21:47
How to update polyPatchbs localPoints liu OpenFOAM Running, Solving & CFD 6 December 30, 2005 17:27


All times are GMT -4. The time now is 20:04.