CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > SU2 > SU2 Shape Design

Discrete adjoint cannot be run in parallel if variable load actuator disc model used

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

LinkBack Thread Tools Search this Thread Display Modes
Old   June 15, 2022, 05:45
Default Discrete adjoint cannot be run in parallel if variable load actuator disc model used
New Member
Fabian BŁeler
Join Date: Jun 2022
Posts: 1
Rep Power: 0
FabianBueeler is on a distinguished road
Dear community,

I am trying to create a 3D shape optimisation of a wing with a propeller. I use the discrete adjoint and the variable load actuator disc model from Ettore Saetta for the propeller.

To create the shape optimisation, I first created a direct flow simulation of the wing with the propeller. Then I created a shape optimization of the wing without propeller. This worked well, both the direct flow and the adjoint solvers converged, and the optimizer found a reasonable result. For these simulations I used a parallel compilation of SU2 7.3.1 on linux and between 8 and 32 cores.

Next I added the actuator disc to the shape optimization and ran into a problem. The direct flow simulation worked, but the discrete adjoint simulation did not. The error is:

Error in function AMPI_Request_free (../externals/medi/include/medi/ampi/async.hpp:249)
The message is: Freeing a handle that is not finish with wait, waitall, etc..

I also wanted to upload an image of the error, but CFD online seems to have problems with that. Have a look at log_output.txt in the link below to see the error in the log file.

The error occurred before the discrete adjoint solver started. I had a look at the file async.hpp mentioned in the error message. It seemed to me as if it belongs to MeDiPack and has something to do with running the MPI. So I ran the same simulation with a serial compilation of SU2 7.3.1. And the optimization with the actuator disc worked. The discrete adjoint worked, created reasonable results, and the optimization converged and created reasonable results too.

Moreover, the error does not appear if either the actuator disc model is not used or if the continuous adjoint is used. So my interpretation is, that the problem does not lie in the discrete adjoint or in the actuator disc model itself, but somewhere in the MPI.

I have no knowledge on the MPI, the discrete adjoint or the variable load actuator disc model on code basis. Furthermore, it is important for my project to do shape optimization with the actuator disc model in parallel, since serial optimizations take far too long. So I would be very happy about any inputs or help from anybody of you. Does anybody have an idea, what the problem might be or how it could be solved?

You can find all necessary files to recreate the error here:

log_output.txt contains the output of the erroneous simu.
log_output_serial.txt contains the output of the serial simu (without error). flow.su2 is the mesh with the FFD box, flow_old.su2 is the mesh without FFD box.

Thanks a lot for your help

FabianBueeler is offline   Reply With Quote

Old   June 23, 2023, 15:56
New Member
Join Date: Aug 2022
Posts: 8
Rep Power: 3
L C is on a distinguished road
Hi Fabian, I tried to open your link but it seems to be broken. I've experimented with parallel SU2 and sometimes it is finicky indeed. Not something straightforward to get working right.
L C is offline   Reply With Quote


actuator disc, adjoint, mpi

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
Discrete Element Method Collision Model sfn Fluent Multiphase 5 May 27, 2018 09:25
Explicitly filtered LES saeedi Main CFD Forum 16 October 14, 2015 11:58
simpleFoam in parallel issue plucas OpenFOAM Running, Solving & CFD 3 July 17, 2013 11:30
different amounts of a variable while parallel run with clusters Alish1984 OpenFOAM Programming & Development 0 September 30, 2011 17:15
Replace periodic by inlet-outlet pair lego CFX 3 November 5, 2002 20:09

All times are GMT -4. The time now is 04:56.