Computing wrong adjoint for multizone problems
 

Hi there,

I am currently working on a multizone FSI shape optimization problem of a 2D NACA 2412 airfoil. So basically, my design problem will run without problem only when 1 adjoint of a objective function is computed but when multiple adjoints have to be computed, it fails. In my case, I am setting an objective function (OPT_OBJECTIVE) as DRAG and one of the constraints (OPT_CONSTRAINT) as LIFT > 0.5.

Running the solution, it fails because of the ADJOINT_LIFT solution. Looking through the log_adjoint.out file for the ADJOINT_LIFT solution, it is shown that for zone 0 (fluid), it is computing the drag adjoint instead, but for zone 1 (structure), it is computing the lift adjoint, which is why when the python scripts are trying to find the fsisteady_adj_res_cl_0.csv file it is there is no file for it because there is only fsisteady_adj_res_cd_0.csv. It is also noted that the ADJOINT_DRAG has already been computed at this point and I confirm that it is correct for both zones.

Attachment 82442

Attachment 82437

Attachment 82441

Any idea why its computing the wrong adjoint? Is there any possible fix for this, or even so, may I find where within the source code this problem could be arising? I have attached my config files, mesh file and the design solution file below:

https://drive.google.com/drive/folde...QB?usp=sharing

Thank you guys so much in advance, and hope to hear from you soon!

Okay, nevermind, I found the solution. I realized I defined OBJECTIVE_FUNCTION = DRAG in the fluid zone which conflicts with the adjoint lift calculations for that zone.