[yashar.afarin]

Hi Alberto,

thanks for sharing your corrected dynamicSmagorinsky. I have complied it without any problem. but when I want to decompose my mesh, the below warning wolud be appear( I am using OF-1.6.x and Fedora-12):
--> FOAM Warning :
From function dlLibraryTable:pen(const fileName& functionLibName)
in file db/dlLibraryTable/dlLibraryTable.C at line 79
could not load /home/Yashar/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libdynamicSmagorinskyModel.so: undefined symbol: _ZNK4Foam14incompressible9LESModels

I would appreciate if you could help me to solve it.

[alberto]

Hi Yashfar,

nothing to worry about. You can safely ignore the warning. It happens sometime when using libraries through controlDict.

Best,

[yashar.afarin]

Originally Posted by alberto

Hi Yashfar,

nothing to worry about. You can safely ignore the warning. It happens sometime when using libraries through controlDict.

Best,

Hi Alberto,

thanks for your reply and your sharing.

sincerely yours.

[yashar.afarin]

Hi Alberto,

I would appreciate if you could answer me that " how could this model be used in compressible solver?"

Best regard,

[alberto]

Hi Yashar,

it cannot. The implementation I released works for incompressible solvers. For the compressible version you need to implement a compressible dynamic Smagorinsky model (see for example Piomelli's papers to have an idea).

Best,

[BernhardGrieser]

Alberto,

thank you very much for your efforts on implementing a truly dynamic Smagorinsky model! I realized the flaw of the old implementation a few weeks ago and I am more than happy that someone already took care of the issue.

I'm running my simulations as a non-root user on my university's HP XC3000, where OpenFOAM 1.7 is pre-installed.
Logically it produces an error when compiling with wmake libso:

...
/usr/bin/ld: cannot open output file /software/all/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/libdynamicSmagorinskyModel.so: Permission denied
collect2: ld returned 1 exit status
make: *** [/software/all/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/libdynamicSmagorinskyModel.so] Error 1

Do you know a work-around, so I can compile it as non-root?

Best regards,
Bernhard

[BernhardGrieser]

ok, problem solved. One only has to change the content of ./Make/files:
I replaced FOAM_LIBBIN with FOAM_USER_LIBBIN

[alberto]

Good the code is useful to someone :-)

[JamesJCFD]

Dear Alberto,

Many thanks for your contribution! This is just what I have been looking for...

kind regards,

James

[alberto]

Happy to share!

[PerryLJohnson]

Alberto,

Just trying to understand your code (and OpenFOAM coding in general), not being picky.

If I have understood your code correctly, by changing the average(LLMM/MMMM) to fvc::average(LLMM/MMMM), you have averaged over the faces belonging to the cells, which are interpolated from their respective cells.

So essentially, you have applied the box filter (or simple filter as OpenFOAM calls it) and could accomplish the same thing using filter_(LLMM/MMMM) assuming the user chooses the simple filter. This would also allow use of the other filters to perform this averaging at the choice of the user. Please correct me if I have misunderstood anything.

FYI, the fluent manual states that it averages locally using the test filter, which appears to be what you have done as well (fluent also uses box filter).

Regards,
Perry

[alberto]

Originally Posted by PerryLJohnson

Alberto,

Just trying to understand your code (and OpenFOAM coding in general), not being picky.

If I have understood your code correctly, by changing the average(LLMM/MMMM) to fvc::average(LLMM/MMMM), you have averaged over the faces belonging to the cells, which are interpolated from their respective cells.

So essentially, you have applied the box filter (or simple filter as OpenFOAM calls it) and could accomplish the same thing using filter_(LLMM/MMMM) assuming the user chooses the simple filter. This would also allow use of the other filters to perform this averaging at the choice of the user. Please correct me if I have misunderstood anything.

Yes, correct.

FYI, the fluent manual states that it averages locally using the test filter, which appears to be what you have done as well (fluent also uses box filter).

Yes. Notice however that FLUENT forces nu_sgs > 0, which is not the case, and I do not do.

Best,

[PerryLJohnson]

Originally Posted by alberto

Yes, correct.



Yes. Notice however that FLUENT forces nu_sgs > 0, which is not the case, and I do not do.

Best,

I did notice that. I know you want to allow for back-scattering, but is there a numerical danger in letting nuEff get *close* to zero? I guess not since the numerical diffusion would take over at this point (provided you were using a dissipative scheme). They also clip nuSgs at 0.23 which seems kinda arbitrary.

Perry

[alberto]

Originally Posted by PerryLJohnson

I did notice that. I know you want to allow for back-scattering, but is there a numerical danger in letting nuEff get *close* to zero? I guess not since the numerical diffusion would take over at this point (provided you were using a dissipative scheme).

Numerically I do not see any problem. It would just be like solving for an inviscid flow. Physically it could not be meaningful however.

They also clip nuSgs at 0.23 which seems kinda arbitrary.

Yes, I agree.

Best,

[kpsl]

Hi,

I am having some trouble compiling this on OpenFOAM 2.0.x.

I get the following error:

dynLocalAverageSmagorinsky.C:40:   instantiated from here
/home/kris/OpenFOAM/OpenFOAM-2.0.x/src/turbulenceModels/incompressible/LES/lnInclude/LESModel.H:115: error: no matching function for call to ‘Foam::incompressible::LESModels::dynLocalAverageSmagorinsky::dynLocalAverageSmagorinsky(const Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&, const Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::transportModel&, const Foam::word&)’
dynLocalAverageSmagorinsky.H:130: note: candidates are: Foam::incompressible::LESModels::dynLocalAverageSmagorinsky::dynLocalAverageSmagorinsky(const Foam::volVectorField&, const Foam::surfaceScalarField&, Foam::transportModel&)
dynLocalAverageSmagorinsky.H:118: note:                 Foam::incompressible::LESModels::dynLocalAverageSmagorinsky::dynLocalAverageSmagorinsky(const Foam::incompressible::LESModels::dynLocalAverageSmagorinsky&)
make: *** [Make/linux64GccDPOpt/dynLocalAverageSmagorinsky.o] Error 1

Is there an updated version for 2.0? Any help would be much appreciated

Kris

[Bernhard]

This was released by Alberto this morning:
https://github.com/AlbertoPa/dynamicSmagorinsky

[grandgo]

hi alberto,

i'm using dynamicSmagorinsky for a LES in parallel processing.
i have this message. the calculation is running but i'm concerned because of this warning.

FOAM Warning :
    From function dlLibraryTable::open(const fileName& functionLibName)
    in file db/dlLibraryTable/dlLibraryTable.C at line 78
    could not load /home/stss8/OpenFOAM/stss8-1.7.x/lib/linux64GccDPOpt/libdynamicSmagorinskyModel.so: undefined symbol: _ZNK4Foam14incompressible9LESModels11GenEddyVisc1BEv

do you know anything about this?

best regards
grandgo

[PerryLJohnson]

grandgo,

I think your question is posed and answered in posts 21 and 22 respectively.

Perry

[grandgo]

Originally Posted by PerryLJohnson

grandgo,

I think your question is posed and answered in posts 21 and 22 respectively.

Perry

right. sorry.

(couldn't find the keywords via search function ...)

[alberto]

Originally Posted by kpsl

Is there an updated version for 2.0? Any help would be much appreciated

I released it this morning. Note that the version for OpenFOAM 2.0.x *breaks* the compatibility with 1.7.x.
If you need both, take a look at the README. However I am not going to update the version for 1.7.x in case bugs are found.

P.S. Notice that OF 2.0.x comes with the dynLagrangian model, which uses Lagrangian averaging to find the dynamic coefficient. It is a more refined model than the traditional dynamic Smagorinsky, so it might be worth considering it too (the reference is in the header file).

Best,