[WIP] Implementation of Amplification Factor Transport(AFT) 2019b tranistion model by sun5k · Pull Request #2422 · su2code/SU2

sun5k · 2025-01-15T13:05:57Z

Proposed Changes

Hi, all.

This PR is for implementing well known subsonic transition model (SA-AFT2019b).

Coder's SA-AFT2019b model is composed of the amplification factor and the logarithmic intermittency transport equation as shown below.

Current State

The primary implementation of the model has been completed and is now in the validation. The validation problems are refer : https://doi.org/10.2514/6.2023-3530 and https://doi.org/10.2514/6.2019-0039.

The transition onset location and the amplification factor field are slightly different. I'm still debugging and checking various things.

Validation Problem of Current State

Firstly, I'm trying the validation problem for 2-D problems (flat plate and NLF-0416 airfoil).
The current results are shown below.

S&K Flat plate from https://doi.org/10.2514/6.2023-3530

NLF-0416 airfoil from https://doi.org/10.2514/6.2023-3530

PR Checklist

Put an X by all that apply. You can fill this out after submitting the PR. If you have any questions, don't hesitate to ask! We want to help. These are a guide for you to know what the reviewers will be looking for in your contribution.

I am submitting my contribution to the develop branch.
My contribution generates no new compiler warnings (try with --warnlevel=3 when using meson).
My contribution is commented and consistent with SU2 style (https://su2code.github.io/docs_v7/Style-Guide/).
I used the pre-commit hook to prevent dirty commits and used pre-commit run --all to format old commits.
I have added a test case that demonstrates my contribution, if necessary.
I have updated appropriate documentation (Tutorials, Docs Page, config_template.cpp), if necessary.

Primary back up

Test for SAAFT2019b

Renew_Min_Max

HL

Common/include/option_structure.hpp

SU2_CFD/include/numerics/turbulent/transition/trans_sources.hpp

SU2_CFD/src/output/CFlowOutput.cpp

      break;

+    case TURB_TRANS_MODEL::AFT:
+      //nodes -> SetAFT_Wonder_Func(iPoint, HL, H12, dNdRet, Ret0, D_H12, l_H12, m_H12, kv, Rev, Rev0, F_growth, F_crit, PAF, Pg, Dg);


SU2_CFD/src/solvers/CTransAFTSolver.cpp

AFT2017b

SU2_CFD/include/numerics/turbulent/transition/trans_sources.hpp

Test_Strain

pcarruscag

Thanks 👍 Looks good but we can re-use a lot of code from the LM solver and see the question about the aux var gradient which looks suspicious.

Common/include/option_structure.hpp

Common/src/CConfig.cpp

SU2_CFD/include/numerics/turbulent/transition/trans_correlations.hpp

SU2_CFD/include/numerics/turbulent/transition/trans_diffusion.hpp

SU2_CFD/src/solvers/CTransAFTSolver.cpp

pcarruscag · 2025-01-19T23:13:52Z

SU2_CFD/src/solvers/CTransAFTSolver.cpp

+  switch (config->GetKind_Gradient_Method()) {
+    case GREEN_GAUSS:
+      SetAuxVar_Gradient_GG(geometry, config);
+      break;
+    case WEIGHTED_LEAST_SQUARES:
+      SetAuxVar_Gradient_LS(geometry, config);
+    default:
+      break;
+    }


Why are you computing the AuxVarGradient after you use it above instead of before?

I wanted to discuss this with someone.

On the SA-AFT2019b model, HL is defined as follows :

where, d is the distance from the wall.

The AuxVar[0] is for the d and AuxVarGradient[0] is for the . So, the d and are the constant value when the geometry and mesh are not deformed.

And, the Auxvar[1] is for the and AuxVarGradient[1] is for the .

The velocity values change with each iteration. So, the Auxvar[1] and it's gradient must be redefined and calculated each iteration. Therefore, I decided to do the gradient calculation in the Postprocessing.

On second thought, the Auxvar assign in "trans_sources.hpp", or would you have any other ideas?

In the paper it is mentioned that grad d is the wall-normal unit vector, don't we have that stored already? I remember seeing it in another PR, maybe it was not merged yet.

In the paper it is mentioned that grad d is the wall-normal unit vector, don't we have that stored already? I remember seeing it in another PR, maybe it was not merged yet.

OK. I'll try to find it. :D

I think it is in the simplified 2015 PR

If I understood correctly then yes, it is in the simplified transition model PR #1901. The branch is still under development (currently working on the validation on T3 series flat plates and Eppler airfoil) but that part has already been implemented and verified.

@rois1995

Hi, if I understood correctly, the modifications in #1901 were introduced to obtain the unit normal vector at nodes. Is that right?

Commit:
b1edf61

Yes, it sets at each point in the mesh the wall normal of the associated point on a viscous surface.

SU2_CFD/src/solvers/CTransAFTSolver.cpp

pcarruscag · 2025-01-19T23:17:29Z

SU2_CFD/src/solvers/CTransAFTSolver.cpp

+    su2double Temp3 = flowNodes->GetVelocity(iPoint, 0) * nodes->GetAuxVarGradient(iPoint, 0, 0);
+    su2double Temp4 = flowNodes->GetVelocity(iPoint, 1) * nodes->GetAuxVarGradient(iPoint, 0, 1);


Add some documentation please, if the initial aux var grad is 0 then you set the aux var to 0, then its gradient is 0 and this remains 0 forever? Looks like a bug.

pcarruscag · 2025-01-19T23:18:50Z

SU2_CFD/src/solvers/CTransAFTSolver.cpp

+
+    /*-- destruction term of Intermeittency(Gamma) --*/
+    const su2double Dg = c_2 * Density_i * VorticityMag * F_turb * (c_3 * exp(lnIntermittency) - 1.0);
+    nodes -> SetAFT_Wonder_Func(iPoint, HL, H12, dNdRet, Ret0, D_H12, l_H12, m_H12, kv, Rev, Rev0, F_growth, F_crit, PAF, Pg, Dg);


This is a lot of memory and these variables are only used for post processing (output) right?
If so these can be computed in the output class when setting the output values.
That saves doing this work every iteration and storing all these variables.

Thank you for your advise.
Current class, this function is for debugging. Later, I'll change to compute in the output class.

sun5k · 2025-01-31T07:07:04Z

Thanks 👍 Looks good but we can re-use a lot of code from the LM solver and see the question about the aux var gradient which looks suspicious.

Thank you for comment @pcarruscag. I'll try to address each of your comments one by one.

For now, the SA-AFT2017b model is for testing purposes, so I'll delete the relevant parts.

Common/src/CConfig.cpp

+                case TURB_MODEL::NONE: SU2_MPI::Error("No turbulence model has been selected but AFT transition model is active.", CURRENT_FUNCTION); break;
+                case TURB_MODEL::SST: SU2_MPI::Error("k-w SST turbulence model has been selected but AFT transition model is active.", CURRENT_FUNCTION); break;
+              }
+              cout << "-2017b" << endl;  break;


Common/src/CConfig.cpp

+                case TURB_MODEL::NONE: SU2_MPI::Error("No turbulence model has been selected but AFT transition model is active.", CURRENT_FUNCTION); break;
+                case TURB_MODEL::SST: SU2_MPI::Error("k-w SST turbulence model has been selected but AFT transition model is active.", CURRENT_FUNCTION); break;
+              }
+              cout << "-2017b" << endl;  break;


sun5k · 2025-05-20T00:08:30Z

The original plan was to finish it before the new version (v8.2.0) was released, but I was busy and couldn't manage it. Now that I have some spare time, I'll start working on it again next week!

Decided to exclude AFT2017b

…ture_AFT_2019_Model

…into feature_AFT_2019_Model

sun5k and others added 5 commits January 6, 2025 23:50

Primary SA_AFT_2019b push

e652a53

Primary back up

Merge branch 'su2code:master' into feature_AFT_2019_Model

e1e12ca

Test_SA_AFT2019b

058ee6e

Test for SAAFT2019b

Renew_Min_Max

448f5b9

Renew_Min_Max

HL update

a583d7c

HL

sun5k added the changelog:feature label Jan 15, 2025

github-advanced-security bot found potential problems Jan 15, 2025

View reviewed changes

Updat_AFT2017b

5349d97

AFT2017b

github-advanced-security bot found potential problems Jan 16, 2025

View reviewed changes

SU2_CFD/include/numerics/turbulent/transition/trans_sources.hpp Fixed Show fixed Hide fixed

SA_AFT2019b

0de4e65

Test_Strain

pcarruscag reviewed Jan 19, 2025

View reviewed changes

Merge branch 'develop' into feature_AFT_2019_Model

5686134

github-advanced-security bot found potential problems Feb 10, 2025

View reviewed changes

sun5k and others added 5 commits May 20, 2025 09:11

Merge branch 'develop' into feature_AFT_2019_Model

c466940

Update Feedback

5e81dd0

Decided to exclude AFT2017b

Merge remote-tracking branch 'origin/feature_AFT_2019_Model' into fea…

c947b3e

…ture_AFT_2019_Model

Merge branch 'develop' into feature_AFT_2019_Model

2960b1a

Merge branch 'feature_AFT_2019_Model' of https://github.com/sun5k/SU2 …

d73a77d

…into feature_AFT_2019_Model

		su2double Temp3 = flowNodes->GetVelocity(iPoint, 0) * nodes->GetAuxVarGradient(iPoint, 0, 0);
		su2double Temp4 = flowNodes->GetVelocity(iPoint, 1) * nodes->GetAuxVarGradient(iPoint, 0, 1);

Conversation

sun5k commented Jan 15, 2025

Proposed Changes

Current State

Validation Problem of Current State

PR Checklist

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Check notice

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Uh oh!

pcarruscag left a comment

Choose a reason for hiding this comment

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

rois1995 Mar 4, 2025 • edited Loading Uh oh! There was an error while loading. Please reload this page.

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

sun5k commented Jan 31, 2025

Uh oh!

Check warning

Check warning

sun5k commented May 20, 2025

Uh oh!

Reviewers

Assignees

Labels

Projects

Milestone

Development

Uh oh!

4 participants

rois1995 Mar 4, 2025 •

edited

Loading