|
[Sponsors] |
March 26, 2016, 20:39 |
Hicks-Henne Geometric Sensitivities Matrix
|
#1 |
New Member
Outside US
Join Date: May 2013
Posts: 6
Rep Power: 13 |
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? |
|
March 28, 2016, 18:25 |
|
#2 | |
Senior Member
Heather Kline
Join Date: Jun 2013
Posts: 309
Rep Power: 14 |
Quote:
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. |
||
April 1, 2016, 12:16 |
|
#3 |
Member
Mandar Kulkarni
Join Date: Nov 2013
Location: Virginia Tech, Blacksburg, VA
Posts: 52
Rep Power: 13 |
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 |
|
Tags |
adjoint, geometric, sensitivities |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Using multiple Hicks Henne deformations | jack1980 | SU2 Shape Design | 3 | July 30, 2019 16:01 |
Convection Diffusion 1-D Matrix Problem | balrog | Main CFD Forum | 1 | November 3, 2014 14:34 |
Force can not converge | colopolo | CFX | 13 | October 4, 2011 23:03 |
OpenFOAM version 1.6 details | lakeat | OpenFOAM Running, Solving & CFD | 42 | August 26, 2009 22:47 |
How to update polyPatchbs localPoints | liu | OpenFOAM Running, Solving & CFD | 6 | December 30, 2005 18:27 |