SpectEddyVisc LES model

maka · January 23, 2009, 13:12

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:

stream&)
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, 13:23

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, 17:43

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

January 23, 2009, 13:12	Choosing spectEddyVisc as LESm	#1
maka Senior Member Maka Mohu Join Date: Mar 2009 Posts: 305 Rep Power: 18	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:stream&) 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.

January 23, 2009, 13:23	Try with tmp	#2
henry Senior Member Join Date: Mar 2009 Posts: 854 Rep Power: 22	Try with tmp<volscalarfield> spectEddyVisc::k() const { volScalarField Eps = 2nuEff()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(Epsnu(), 1.0/2.0) erfc(cK4_pow(delta(), -2.0/3.0)sqrt(nu())*pow(Eps, -1.0/6.0)); } H

January 24, 2009, 17:43	Thanks Henry. It is running no	#3
maka Senior Member Maka Mohu Join Date: Mar 2009 Posts: 305 Rep Power: 18	Thanks Henry. It is running now. Best regards, Maka.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Radiation P1 Model with Spectral Model: Multigray	A. Thellmann	CFX	0	October 25, 2008 12:44
Surface Reaction Model or Heterogeneous Rxn Model?	satheeskannan	FLUENT	0	September 2, 2008 18:55
changing model constants in k-e turbulence model	Sunil	CFX	3	October 3, 2006 12:12
Symmetry model half or total geometrical model in LES	vvqf	OpenFOAM Pre-Processing	3	March 22, 2006 09:06
Copying boundary conditions from model to model	Lior	FLUENT	3	September 23, 2004 16:10