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

 March 5, 2014, 19:14 About the computation of the gradient of Volume grid node movement #1 Member   Tommy Chen Join Date: Mar 2011 Location: University of Michigan Posts: 92 Rep Power: 6 Hi Guys: I want to know some details about the computation of the gradient of Volume grid node movement respect to the surface grid node in SU2. Could anyone point out that part of code related to this functionality ? Thanks !

March 10, 2014, 17:17
#2
Super Moderator

Francisco Palacios
Join Date: Jan 2013
Location: Stanford, CA
Posts: 301
Rep Power: 5
Quote:
 Originally Posted by Tommy Chen Hi Guys: I want to know some details about the computation of the gradient of Volume grid node movement respect to the surface grid node in SU2. Could anyone point out that part of code related to this functionality ? Thanks !
Mmm, maybe I have misunderstood the question. But, if you are asking for the volumetric sensitivity (how changes the volumetric mesh when we perturb the surface geometry), the continuous adjoint formulation does not require that step (everything is computed only using surface information.

Best,
Francisco

March 10, 2014, 18:40
#3
Member

Tommy Chen
Join Date: Mar 2011
Location: University of Michigan
Posts: 92
Rep Power: 6
Quote:
 Originally Posted by fpalacios Mmm, maybe I have misunderstood the question. But, if you are asking for the volumetric sensitivity (how changes the volumetric mesh when we perturb the surface geometry), the continuous adjoint formulation does not require that step (everything is computed only using surface information. Best, Francisco
Hi Francisco:

Is this partially due to the assumption that the perturbation of the surface grid node is normal to the local wall surface when you formulate the sensitivity equation?

Since I am much more familiar with discrete adjoint, in which the chain rule has to be used , including the derivatives of volume mesh node displacements respect to every surface mesh node displacement, in order to compute the final gradient of dv. I do not understand why on surface information is necessary.

So , the outputed surface adjoint file after running the SU2_CFD adjoint solver is the only file needed to compute the projected gradient, is it right ?

many thanks

 March 11, 2014, 01:56 #4 Super Moderator   Thomas D. Economon Join Date: Jan 2013 Location: Stanford, CA Posts: 260 Rep Power: 5 Hi Tommy, Indeed, this is the case for the type of continuous adjoint based on shape calculus that we have implemented in SU2. The resulting variation of the objective function is simply an integral over the design surface (no volume mesh sensitivities), and the integrand is what we call the surface sensitivity, which is computed at each mesh node. The surface sensitivity is a measure of the change in the objective function given a small perturbation in the local normal direction for each surface mesh node. Hope this helps clear things up, T

March 12, 2014, 17:50
#5
Member

Tommy Chen
Join Date: Mar 2011
Location: University of Michigan
Posts: 92
Rep Power: 6
Quote:
 Originally Posted by economon Hi Tommy, Indeed, this is the case for the type of continuous adjoint based on shape calculus that we have implemented in SU2. The resulting variation of the objective function is simply an integral over the design surface (no volume mesh sensitivities), and the integrand is what we call the surface sensitivity, which is computed at each mesh node. The surface sensitivity is a measure of the change in the objective function given a small perturbation in the local normal direction for each surface mesh node. Hope this helps clear things up, T
Hi Tom

Could you recommend some related material to me about the adjoint method used in SU2 ? Paper, books or whatever that can interpret the theory

Thank you

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Presly ParaView 0 June 12, 2013 20:34 Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00 Marcus Main CFD Forum 3 December 5, 2000 01:25 Chuck Leakeas Main CFD Forum 2 May 26, 2000 11:18 Erich F. Main CFD Forum 8 October 23, 1999 14:43

All times are GMT -4. The time now is 01:51.