Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Unsteady optimization in V5.0 #377

Closed
mghaffari opened this issue Feb 19, 2017 · 2 comments
Closed

Unsteady optimization in V5.0 #377

mghaffari opened this issue Feb 19, 2017 · 2 comments
Assignees
Labels

Comments

@mghaffari
Copy link

mghaffari commented Feb 19, 2017

Unsteady optimization in V5.0 has some issues. The problem is that after direct iterations SU2 cannot calculate adjoint for each time step properly.
It just keep calculating adjoint equation for one time step. For example a part of the log_adjoint file for NACA64A010 optimization is as follow:

Loading flow solution from direct iteration 250.
Begin direct solver to store flow data (single iteration).
Compute residuals to check the convergence of the direct problem.
log10[Maximum residual]: -2.77153, located at point 1018.
End direct solver, begin adjoint problem.

Min DT: 0.000323558.Max DT: 1.20067.Dual Time step: 0.677451.
log10[Maximum residual]: -3.34769.
Maximum residual point 134, located at (0.63649, 0.0371117).

IntIter ExtIter Res[Psi_Rho] Res[Psi_E] Sens_Geo Sens_AoA
0 0 -4.444881 -4.159977 4.4996e+00 -6.5737e-04
10 0 -4.488648 -4.764869 4.4053e+00 -6.5737e-04
20 0 -4.433964 -4.864120 2.7194e+00 -6.5737e-04
30 0 -4.449170 -4.851557 1.7148e+00 -6.5737e-04
40 0 -4.669063 -4.954710 1.1593e+00 -6.5737e-04
50 0 -4.689637 -5.116023 9.2103e-01 -6.5737e-04
60 0 -4.797857 -5.289007 8.5846e-01 -6.5737e-04
70 0 -4.922113 -5.428807 8.9570e-01 -6.5737e-04
80 0 -5.021316 -5.534515 9.5224e-01 -6.5737e-04
90 0 -5.102257 -5.604468 9.9700e-01 -6.5737e-04
100 0 -5.169358 -5.654203 1.0424e+00 -6.5737e-04

Min DT: 0.0003236.Max DT: 1.201.Dual Time step: 0.6775.
log10[Maximum residual]: -4.172.
Maximum residual point 693, located at (0.03989, -0.0254).

IntIter ExtIter Res[Psi_Rho] Res[Psi_E] Sens_Geo Sens_AoA
0 0 -5.221491 -5.704096 1.0838e+00 -6.5737e-04
10 0 -5.262729 -5.764520 1.1144e+00 -6.5737e-04
20 0 -5.306548 -5.835208 1.1305e+00 -6.5737e-04
30 0 -5.363881 -5.911928 1.1291e+00 -6.5737e-04
40 0 -5.436887 -5.991739 1.1172e+00 -6.5737e-04
50 0 -5.521833 -6.072162 1.1056e+00 -6.5737e-04
60 0 -5.611948 -6.152996 1.0993e+00 -6.5737e-04
70 0 -5.701126 -6.235540 1.0974e+00 -6.5737e-04
...........

the ExtIter is locked to 0

@hlkline hlkline added the bug label Mar 14, 2017
@talbring talbring self-assigned this Jun 14, 2017
@talbring
Copy link
Member

There is indeed a problem there, I will look into it.

@talbring
Copy link
Member

Solved in develop.

@talbring talbring removed their assignment Jul 13, 2017
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

4 participants