sumedhss |
August 9, 2022 15:27 |
Problems with implementing sink term in Gamma-Re-Theta turbulence model
2 Attachment(s)
Hello all,
I'm trying to augment the Gamma-Re-Theta model to include crossflow effects as per Langtry (Langtry, Robin. "Extending the gamma-rethetat correlation based transition model for crossflow effects." In 45th AIAA fluid dynamics conference, p. 2474. 2015.) In this approach, an additional sink term is added to the transport equation for , which is described as-
Here, is the stationery crossflow Reynolds number, which is computed implicitly using a correlation. Since (freestream transition momentum thickness Reynolds number) is also computed implicitly, I defined in a similar fashion as
and called it in
Code:
correctReThetatGammaInt()
, but I get the following warning while compiling-
Code:
In file included from ../turbulenceModels/lnInclude/kOmegaSSTLMN.H:320,
from turbulentTransportModels/myTurbulentTransportModels.C:60:
../turbulenceModels/lnInclude/kOmegaSSTLMN.C: In instantiation of ‘Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > Foam::RASModels::kOmegaSSTLMN<BasicTurbulenceModel>::ReSCF(const Internal&, const Internal&, const Internal&) const [with BasicTurbulenceModel = Foam::IncompressibleTurbulenceModel<Foam::transportModel>; Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Internal = Foam::DimensionedField<double, Foam::volMesh>]’:
../turbulenceModels/lnInclude/kOmegaSSTLMN.C:707:21: required from ‘void Foam::RASModels::kOmegaSSTLMN<BasicTurbulenceModel>::correctReThetatGammaInt() [with BasicTurbulenceModel = Foam::IncompressibleTurbulenceModel<Foam::transportModel>]’
../turbulenceModels/lnInclude/kOmegaSSTLMN.C:785:5: required from ‘void Foam::RASModels::kOmegaSSTLMN<BasicTurbulenceModel>::correct() [with BasicTurbulenceModel = Foam::IncompressibleTurbulenceModel<Foam::transportModel>]’
../turbulenceModels/lnInclude/kOmegaSSTLMN.C:768:6: required from here
I tried running simulations with this model, but they would crash instantly with the error log mentioning . When I replaced with a constant value in
Code:
correctReThetatGammaInt()
, the model compiles without issues and simulations also run without issues. It is clear that the problem is related to the instantiation of the function template, but I am not able to pinpoint the exact reason. I have attached the code and header file along with the post. I'm using OpenFOAM v2012. Any help on this topic would be sincerely appreciated.
Thanks a lot!
|