CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Bugs (http://www.cfd-online.com/Forums/openfoam-bugs/)
-   -   SpectEddyVisc LES model (http://www.cfd-online.com/Forums/openfoam-bugs/62337-specteddyvisc-les-model.html)

maka January 23, 2009 14:12

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.

henry January 23, 2009 14:23

Try with tmp
 
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

maka January 24, 2009 18:43

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


All times are GMT -4. The time now is 08:53.