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

Discrete adjoint - Derivatives of the target function w.r.t. boundary conditions

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 14, 2020, 14:05
Default Discrete adjoint - Derivatives of the target function w.r.t. boundary conditions
  #1
New Member
 
Aldo Serafino
Join Date: Feb 2020
Location: Paris
Posts: 2
Rep Power: 0
CFDEng1234 is on a distinguished road
Hi everyone.

I'm definitely a new user of SU2 v7.0.1 and I'm considering the idea to use it for my PhD thesis.

For the moment, I've been working on a linear cascade of NACA 0012 airfoils.

I model my domain with an INLET, an OUTLET, two periodic surfaces (low and top) and a FIXED WALL (which is the airfoil).

As boundary conditions, I impose the inlet total pressure, the inlet total temperature and the outlet static pressure.

I have two doubts about the discrete adjoint problem.

1) Can I use the discrete adjoint to obtain the gradient of my target function (which is the lift) with respect to my three boundary conditions? In my case, can I get the derivative of the lift with respect to the inlet total pressure, to the inlet total temperature and to the outlet static pressure?

2) For the discrete adjoint problem, instead of using the built-in target functions (i.e. lift or draft), can I define some custom target functions (for instance, the difference of entropy between the outlet and the inlet sections)?

Thank you for your time and have a nice day. Best regards.
CFDEng1234 is offline   Reply With Quote

Old   February 17, 2020, 12:35
Default
  #2
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 466
Rep Power: 13
pcg is on a distinguished road
Yes, you can. But it will require some coding, for boundary conditions something similar is already done for angle of attack and other variables, see:
CDiscAdjSolver::RegisterVariables(...)
and
CDiscAdjSolver::ExtractAdjoint_Variables(...)
You can register/extract other variables in those methods.

To define new objectives functions start at
CDiscAdjSinglezoneDriver::SetObjFunction()
And follow the methods called from there (highly recommended to setup an IDE for that...).

By the way I think at Dirichlet boundaries the adjoint variables are the sensitivities.
pcg is offline   Reply With Quote

Old   February 19, 2020, 11:39
Default
  #3
New Member
 
Aldo Serafino
Join Date: Feb 2020
Location: Paris
Posts: 2
Rep Power: 0
CFDEng1234 is on a distinguished road
Quote:
Originally Posted by pcg View Post
Yes, you can. But it will require some coding, for boundary conditions something similar is already done for angle of attack and other variables, see:
CDiscAdjSolver::RegisterVariables(...)
and
CDiscAdjSolver::ExtractAdjoint_Variables(...)
You can register/extract other variables in those methods.

To define new objectives functions start at
CDiscAdjSinglezoneDriver::SetObjFunction()
And follow the methods called from there (highly recommended to setup an IDE for that...).

By the way I think at Dirichlet boundaries the adjoint variables are the sensitivities.
Thank you for your answer.I'll try to do what you suggest.



Regards.
CFDEng1234 is offline   Reply With Quote

Reply

Tags
derivatives, discrete adjoint, target function


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
Question about adaptive timestepping Guille1811 CFX 25 November 12, 2017 18:38
Boundary Conditions Of Discrete Adjoint Method For 2D NS Equations Ashkaari Main CFD Forum 4 January 27, 2016 21:05
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 07:20
Compile problem ivanyao OpenFOAM Running, Solving & CFD 1 October 12, 2012 10:31
Error with Wmake skabilan OpenFOAM Installation 3 July 28, 2009 01:35


All times are GMT -4. The time now is 17:50.