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

Problem setting up constrained optimization

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 5, 2021, 03:22
Default Problem setting up constrained optimization
  #1
New Member
 
Join Date: Dec 2019
Posts: 2
Rep Power: 0
wJhP is on a distinguished road
Hi,

I have trouble defining an optimization problem. The idea is to maximize the heat transfer on a surface while minimizing the pressure at the inlet.

I tried to setup a constraint in this way:

Code:
MARKER_MONITORING = ( INFLOW, TRANSITION )
MARKER_ANALYZE = ( INFLOW )

OBJECTIVE_FUNCTION= TOTAL_HEATFLUX
OBJECTIVE_WEIGHT = -1.0

OPT_OBJECTIVE=  TOTAL_HEATFLUX*-1
OPT_CONSTRAINT= (  SURFACE_STATIC_PRESSURE < 720 ) * 0.015
The optimization runs and the first deformation is correct: I can improve the objective (heat flux increases) and the pressure is still below the constraint. The following iteration, however, the deformation is much less than the 1st one, basically the same geometry as the original mesh, and then the optimizer stops "successfully".

I am clearly doing something wrong, but I don't understand what.

The design tutorial has this to say about the constraints:


Quote:
Note that this is only compatible with aerodynamic functions, although geometry-based functions can also be included in the OPT_CONSTRAINT setting.
Is the surface pressure an aerodynamic function? Or only lift/drag/pitch are supported?

Also, the parameters OBJECTIVE_FUNCTION and OBJECTIVE_WEIGHT are overwritten when doing the optimization via shape_optimization.py, correct?

I am also trying to do the same optimization using penalty function, which is setup like this:


Code:
OPT_OBJECTIVE= (SURFACE_STATIC_PRESSURE < 720)*0.1 ; TOTAL_HEATFLUX*-1
Here I have the opposite problem: the deformation keeps increasing even tough the constraint on the pressure is not met.


Any ideas are welcome!

Last edited by wJhP; July 5, 2021 at 05:04.
wJhP is offline   Reply With Quote

Old   July 9, 2021, 11:33
Default
  #2
New Member
 
Join Date: Dec 2019
Posts: 2
Rep Power: 0
wJhP is on a distinguished road
For future reference: the optimization with the penalty function runs correctly when the penalty scaling is negative:

Code:
OPT_OBJECTIVE= (SURFACE_STATIC_PRESSURE < 720)*-0.1 ; TOTAL_HEATFLUX*-1



wJhP is offline   Reply With Quote

Reply

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
[snappyHexMesh] determining displacement for added points CFDnewbie147 OpenFOAM Meshing & Mesh Conversion 1 October 22, 2013 09:53
Cells with t below lower limit Purushothama Siemens 2 May 31, 2010 21:58
Problem setting C_MU_T with DEFINE_ADJUST sarah_ron FLUENT 1 February 24, 2005 12:11
Help setting up combustor problem Tom FLUENT 0 February 5, 2005 01:20
Warning 097- AB Siemens 6 November 15, 2004 04:41


All times are GMT -4. The time now is 08:03.