CFD Online Discussion Forums - SpectEddyVisc LES model

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Bugs (https://www.cfd-online.com/Forums/openfoam-bugs/)

- - SpectEddyVisc LES model (https://www.cfd-online.com/Forums/openfoam-bugs/62337-specteddyvisc-les-model.html)

Choosing spectEddyVisc as LESm

Choosing spectEddyVisc as LESmodel in channel 395 tutorial gives the following error using V1.3:

--> FOAM FATAL ERROR : : gsf not dimensionless

From function erfc(const tmp<geometricfield<scalar,> >& gsf)
in file /home/maka/OpenFOAM/OpenFOAM-1.3/src/OpenFOAM/lnInclude/GeometricScalarField.C at line 742.

FOAM aborting

Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&)
Foam::error::abort()
Foam::tmp<foam::geometricfield<double,> > Foam::erfc<foam::fvpatchfield,>(Foam::tmp<foam::ge ometricfield<double,> > const&)
Foam::LESmodels::spectEddyVisc::k() const
Foam::LESmodels::GenEddyVisc::B() const
channelOodles [0x41c042]
__libc_start_main
__gxx_personality_v0

Best regards,
Maka.

Try with

tmp<volscalarfield> spectEddyVisc::k() const
{
volScalarField Eps = 2*nuEff()*magSqr(symm(fvc::grad(U())));

return
cK1_*pow(delta(), 2.0/3.0)*pow(Eps, 2.0/3.0)
*exp(-cK2_*pow(delta(), -4.0/3.0)*nu()/pow(Eps, 1.0/3.0))
- cK3_*pow(Eps*nu(), 1.0/2.0)
*erfc(cK4_*pow(delta(), -2.0/3.0)*sqrt(nu())*pow(Eps, -1.0/6.0));
}

H

Thanks Henry. It is running no

Thanks Henry. It is running now.
Best regards,
Maka.