[godfatherBond]

Hi FOAMers,
If anybody missed the compressible dynamic Smagorinsky model in the newer versions of OF (v4, v5,...), I am attaching a working version known to work with v4 and v5 (and above). The code is very well commented, hence it is easier to understand.
Also, some definitions (like filtered strain rate,..etc) are Favre Filtered to be consistent with the literature.It seems that OpenFOAM's dynamic models are not compressible per se. This may be not be useful for compressible flows as well as low Mach reacting flows with a variable density.
You may use the attached version it if you want.

Please let me know you comments.

[Paul_AERO]

Hi Maximus,

I sent you a personal message to ask the compilation of your source codes.
Could you check the inbox?

Thank you.

Paul

[godfatherBond]

Originally Posted by Paul_AERO

Hi Maximus,

I sent you a personal message to ask the compilation of your source codes.
Could you check the inbox?

Thank you.

Paul

For your questions on the personal message.
1. Check the .H file. The reference is mentioned there.
2. For the turbulent properties, you just need to specify the "secondary" filter type. The constants are evaluated dynamically.

Cheers..

[Paul_AERO]

Originally Posted by godfatherBond

For your questions on the personal message.
1. Check the .H file. The reference is mentioned there.
2. For the turbulent properties, you just need to specify the "secondary" filter type. The constants are evaluated dynamically.

Cheers..

Thanks Maximus.

1. For your reply for (2), secondary filter means:
dynSmagorinskyCompCoeffs
{
filter simple;
}
like above? Thanks for your help.

2. In the source code, dynSmagorinsky.C, line 132-172,
The coefficient, CI ( equ(17) in the M. Pino Martin's article) seems to be calculated. I was wondering why both CI and Ck are used. What does Ck mean? The same applies to Cd/Cs, line 74-128.
+[2021.9.10]: I figured out Cs and Ck to be used for the output of Cd and CI, respectively. Cs and Ck are calculated with strain rate and filter.

3. Even though the compilation is made without critical errors, the run fails saying, rhoPimpleFoam:
symbol lookup error: /opt/OpenFOAM/OpenFOAM-v2106/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so: undefined symbol: _ZN4Foam9LESModels18dynSmagorinskyCompINS_15EddyDi ffusivityINS_18ThermalDiffusivityINS_27Compressibl eTurbulenceModelINS_11fluidThermoEEEEEEEE6calcCfEv

I included turbulenceModels/transportModels in folders of 'compressible' and 'turbulenceModels.' Could you give me a hint for what library should be included in the sources?
+[2021.9.10]: I successfully compiled codes with changing one line in the source code. The line, "virtual void calcCf(); //MA 19.JAN.2020", in dynSmagorinskyComp.H is not used in the computing process, which gives arise the problem in compilation.

Thank you.

Paul

[godfatherBond]

Originally Posted by Paul_AERO

Thanks Maximus.

1. For your reply for (2), secondary filter means:
dynSmagorinskyCompCoeffs
{
filter simple;
}
like above? Thanks for your help.

2. In the source code, dynSmagorinsky.C, line 132-172,
The coefficient, CI ( equ(17) in the M. Pino Martin's article) seems to be calculated. I was wondering why both CI and Ck are used. What does Ck mean? The same applies to Cd/Cs, line 74-128.
+[2021.9.10]: I figured out Cs and Ck to be used for the output of Cd and CI, respectively. Cs and Ck are calculated with strain rate and filter.

3. Even though the compilation is made without critical errors, the run fails saying, rhoPimpleFoam:
symbol lookup error: /opt/OpenFOAM/OpenFOAM-v2106/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so: undefined symbol: _ZN4Foam9LESModels18dynSmagorinskyCompINS_15EddyDi ffusivityINS_18ThermalDiffusivityINS_27Compressibl eTurbulenceModelINS_11fluidThermoEEEEEEEE6calcCfEv

I included turbulenceModels/transportModels in folders of 'compressible' and 'turbulenceModels.' Could you give me a hint for what library should be included in the sources?
+[2021.9.10]: I successfully compiled codes with changing one line in the source code. The line, "virtual void calcCf(); //MA 19.JAN.2020", in dynSmagorinskyComp.H is not used in the computing process, which gives arise the problem in compilation.

Thank you.

Paul

Hi Paul,
I am a bit occupied at the moment. Pardon me if my replies are a bit late.

2. You may find the answers in the Ref. mentioned in the code or in any standard article on dynamic modelling. Hence, I will skip this.

3. I think the error is related to the Makefiles. If you do not add appropriate dependent files in "options". This does not have anything to do with the source code of the turbulence model. It is a minor error and I am sure you will be able to find the solution in a jiffy.

Cheers..

[godfatherBond]

Hi,
If you have compilation issues, try adding the following headers explicitly:

#include "CompressibleTurbulenceModel.H"
#include "compressibleTransportModel.H"
#include "fluidThermo.H"
#include "addToRunTimeSelectionTable.H"
#include "makeTurbulenceModel.H"

#include "ThermalDiffusivity.H"
#include "EddyDiffusivity.H"

#include "LESModel.H"

[kcm2102]

Thanks godfatherBond
The turbulent model is working fine.

Regards
kcm2102

[youdongyao]

Originally Posted by kcm2102

Thanks godfatherBond
The turbulent model is working fine.

Regards
kcm2102

Hello, could you please show me your file for reference? Can you provide a calculation example file if it is convenient? I successfully compiled the file, but it did not work properly. Thank you.

[youdongyao]

Originally Posted by godfatherBond

Hi FOAMers,
If anybody missed the compressible dynamic Smagorinsky model in the newer versions of OF (v4, v5,...), I am attaching a working version known to work with v4 and v5 (and above). The code is very well commented, hence it is easier to understand.
Also, some definitions (like filtered strain rate,..etc) are Favre Filtered to be consistent with the literature.It seems that OpenFOAM's dynamic models are not compressible per se. This may be not be useful for compressible flows as well as low Mach reacting flows with a variable density.
You may use the attached version it if you want.

Please let me know you comments.

Hello, I see that you have mentioned on the forum that dyndynSmagorinskyComp has been successfully tested. I tried to compile dynSmagorinskyComp, it was successful, but when using it it showed an unknown turbulence model, I don't know why, I have added libso to ControlDict. Could you please share your model and test examples with me? Thank you very much!! By the way, I'm using openfoam7.

[youdongyao]

Originally Posted by godfatherBond

Hi FOAMers,
If anybody missed the compressible dynamic Smagorinsky model in the newer versions of OF (v4, v5,...), I am attaching a working version known to work with v4 and v5 (and above). The code is very well commented, hence it is easier to understand.
Also, some definitions (like filtered strain rate,..etc) are Favre Filtered to be consistent with the literature.It seems that OpenFOAM's dynamic models are not compressible per se. This may be not be useful for compressible flows as well as low Mach reacting flows with a variable density.
You may use the attached version it if you want.

Please let me know you comments.

Originally Posted by godfatherBond

Hi Paul,
I am a bit occupied at the moment. Pardon me if my replies are a bit late.

2. You may find the answers in the Ref. mentioned in the code or in any standard article on dynamic modelling. Hence, I will skip this.

3. I think the error is related to the Makefiles. If you do not add appropriate dependent files in "options". This does not have anything to do with the source code of the turbulence model. It is a minor error and I am sure you will be able to find the solution in a jiffy.

Cheers..

Could you please provide the makeTurbulence document? Thank you very much