Applying L_inf Norm to Surface Static Pressure Objective for Adjoint Simulations #2582
Replies: 2 comments · 3 replies
-
|
CAdjEulerSolver is only for continuous adjoint. |
Beta Was this translation helpful? Give feedback.
-
|
You can find the changes in this file. I only changed 3 lines of code to implement norm. Derivatives are not zero. When I check the history output it should show higher values of pressure at DSN_001, but it shows the same value where average is applied. At the end it doesn't even reach the same pressure that is obtained with average pressure calculations. |
Beta Was this translation helpful? Give feedback.
-
|
Is it not normal that the maximum is larger than the average? Maybe you can attend a developer meeting and explain the problem. The schedule is in our Slack channel. |
Beta Was this translation helpful? Give feedback.
-
|
Well, it could be normal if it would start with higher pressures. I don't understand why the first pressure calculation is the same in both averaged and norm-applied surface static pressure. I will try to attend the meeting tomorrow. |
Beta Was this translation helpful? Give feedback.
-
|
I realized when I put SURFACE_STATIC_PRESSURE in the history output (without any modification) for unsteady flows, outputs are shown as 0. Using a custom output for this actually works. I guess this affects the unsteady adjoint solver as well. output is like this: Is this a bug or am I missing something? |
Beta Was this translation helpful? Give feedback.
Uh oh!
There was an error while loading. Please reload this page.
Uh oh!
There was an error while loading. Please reload this page.
-
Hello everybody,
I am using unsteady discrete adjoint for an optimization case where I want to decrease the maximum pressure over a wall. I was using SURFACE_STATIC_PRESSURE until now and I want to alter this to trick the adjoint code so it can optimize the maximum pressure over this wall. I applied L_infinity norm to the SURFACE_STATIC_PRESSURE part in the CFlowOutput::SetAnalyzeSurface and compared the results for some test cases to see if I am getting the max pressure. It works up to a degree that I can use. You can see the equations from the picture below.
However, when I started a shape optimization simulation I don't think it works. Should I also update some parts in the CAdjEulerSolver.cpp?
Can anybody help me with this problem?
Beta Was this translation helpful? Give feedback.
All reactions