About SmagorinskyEffectiveOmega::ForcedWithStress implementation

[aaron]

Hello,

Sourced from the code in src/dynamics/collisionLEC.cse.h (OLB version 1.7, around line 986).

I’m particularly interested in the Smagorinsky LES model used for the collision process in the following class:

struct SmagorinskyEffectiveOmega<BGK, descriptors::D3Q27<FIELDS…>, momenta::ForcedWithStress<momenta::BulkTuple>, equilibria::SecondOrder>

Since this model combines turbulent viscosity and an external force field, I would appreciate a more structured and formatted equation that clarifies the details of this collision process.

(I am familiar with the LES equation and Guo’s force field equation, and just curious about the implementation here).

The numerous xNumber variables seem optimized for performance but make readability challenging.

Moreover, if you have any relevant references, I would greatly appreciate them.

Thank you!

Best regards,
Aaron

[shota]

Dear Aaron,

the “.cse.h” files are automatically generated and optimized regarding common subexpressions.
As you said, they are not designed to be understandable but to run efficiently.

The code is generated from the implementations in collisionLES.h so there you will find readable implementations to have insights into the model.

BR
Shota

[aaron]

Thank you, Shota.

I followed your suggestions, disabled the CSE flag in the Makefile, and GDB successfully captured what you described.

Problem solved!

Best regards,
Aaron

[Author]

Posts