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

Unsteady Adjoint Shape Design

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 11, 2013, 18:46
Default Unsteady Adjoint Shape Design
  #1
New Member
 
Eduardo Molina
Join Date: Sep 2010
Location: Brazil
Posts: 27
Rep Power: 6
EMolina is on a distinguished road
Hi SU2 Developers,

I wanna know that if it is possible to perform an unsteady Adjoint Shape Optimization, basically a flutter two degrees of freedom case.
If not I would like to know some guidelines to perform such modification on the code. SU2 has everything I need an adjoint implementation, mesh deformation and a direct solver.
I appreciate any help from anyone that is studying unsteady optimization.

Thanks.
EMolina is offline   Reply With Quote

Old   July 25, 2013, 02:16
Default
  #2
Super Moderator
 
Thomas D. Economon
Join Date: Jan 2013
Location: Stanford, CA
Posts: 260
Rep Power: 5
economon is on a distinguished road
Hi Eduardo,

Most of what you need for this type of problem is already in place in the code, as unsteady shape design using the Euler equations for prescribed motion cases has been verified and validated (see http://adl.stanford.edu/papers/AIAA-2013-0632.pdf, for example). We're currently working to extend this to the general RANS case.

Do you also mean that you would like to include the aeroelastic response in the formulation using two degrees of freedom (pitch/plunge, for instance)? If so, we are also working on this capability in the code, but an adjoint for the aeroelastic equations would still need to be implemented.

Please let us know what you had in mind, and perhaps we can point you towards the relevant locations in the code. Thanks for your interest in SU2,
Tom
economon is offline   Reply With Quote

Old   July 29, 2013, 17:59
Default
  #3
New Member
 
Eduardo Molina
Join Date: Sep 2010
Location: Brazil
Posts: 27
Rep Power: 6
EMolina is on a distinguished road
Hi Tom,

Thanks for your reply.

Yes you are right, I would like to include the aeroelastic response using a two degrees of freedom (pitch and plunge) and after implement an adjoint for the aeroelastic equations. At this moment, I would like to start implementing the 2 degrees aeroelastic equations, studying how can I integrate this simple 2 equations on the code, translate the CFD surface pressure to the structural equations and perform the mesh deformation and so on.
I have implemented this system using a Theodorsen function for the aerodynamic part, now I would like to extend to SU2.

I appreciate any guidelines in the code to perform such modifications.

Thanks

Eduardo
EMolina is offline   Reply With Quote

Old   August 22, 2013, 16:19
Default
  #4
New Member
 
Santiago Padron
Join Date: May 2013
Posts: 16
Rep Power: 4
Santiago Padron is on a distinguished road
Hi Eduardo,

As Tom mentioned there is a two degree of freedom aeroelastic capability in the code, but without an adjoint. The code solves the aeroelastic equations as described in this paper http://aero-comlab.stanford.edu/Papers/reno94.pdf

At this moment the options that you will need to include in the configuration file are:

% Mesh motion for unsteady simulations (NO, YES)
GRID_MOVEMENT= YES
%
% Type of mesh motion
GRID_MOVEMENT_KIND= AEROELASTIC
%
% Type of aeroelastic grid movement (RIGID, DEFORM)
% This option either rigidly rotates/translates the entire mesh or the interior of the mesh deforms to accommodate the pitch and plunge.
TYPE_AEROELASTIC_GRID_MOVEMENT = DEFORM
% Type of grid velocities for aeroelastic motion (FD, ANALYTIC)
% Use the grid velocity either from the analytic aeroelastic equations or computed it using finite difference based on movement of the grid points.
TYPE_AEROELASTIC_GRID_VELOCITY = FD

% where to plunge and pitch about.
MOTION_ORIGIN_X= -0.5
MOTION_ORIGIN_Y= 0.0

% The springs natural frequencies (basically the spring constants)
FREQ_PLUNGE_AEROELASTIC= 100.0
FREQ_PITCH_AEROELASTIC= 100.0

Other parameters needed to simulate your particular case would have to be hardcoded in the code.

For the aeroelastic implementation, you can start looking at the AEROELASTIC case in the SetGrid_Movement routine in iteration_structure.cpp.

Santiago
Santiago Padron is offline   Reply With Quote

Old   August 26, 2013, 18:19
Default
  #5
New Member
 
Eduardo Molina
Join Date: Sep 2010
Location: Brazil
Posts: 27
Rep Power: 6
EMolina is on a distinguished road
Hi Tom and Santiago...

Many thanks for your reply.

That was exactly what I was looking for. Any doubts/improvements that I will make in the code, I will post here.

Congratilations for the outstanding job you are doing with SU2.

Keep on rockin' in the free world.

Eduardo
EMolina is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Error in Tutorial-Optimal Shape Design of a Fixed Wing cfdbooks2008 SU2 Shape Design 10 January 29, 2013 10:10
SQP method for adjoint based aerodynamic shape optimization cfdbooks Main CFD Forum 2 May 26, 2009 10:40
Collaboration on continuous adjoint design optimiz Carlos Lozano Main CFD Forum 0 April 4, 2008 04:13
What's the status of Unsteady design of Turbomachi Ji Lucheng & Chen Jiang Main CFD Forum 0 December 24, 2002 04:44
Info: Short Course On Thermal Design of Electronic Equipment Arnold Free Main CFD Forum 0 August 10, 1999 10:18


All times are GMT -4. The time now is 21:23.