Calculation of Epsilon in LES model

Yuri Almeida · January 15, 2013, 23:05

Dears,

I'm performing flame simulations using the fireFoam solver. The default set for combustion is the eddy dissipation model, which furnishes the equation below for the source term in the species transport eqn:

wFuelNorm_ = C_*rho_*turbulence_.epsilon()
/max(turbulence_.k(),dimensionedScalar("SMALL",dimV elocity*dimVelocity,SMALL))
*min(YFuel , YO2/s.value());

However, the turbulence model utilized is LES OneEqEddy, which does not directly calculate the epsilon. So, where the epsilon calculation is performed? What expression is used for it?

Thanks in advance,

Yuri

gregor · January 17, 2013, 12:55

Hi Yuri,

the oneEqEddy model is derived from the GenEddyVisc model. In GenEddyVisc.H you find that epsilon is computed as

Code:

 ce_*k()*sqrt(k())/delta();

Yuri Almeida · January 18, 2013, 08:43

Hi Gregor,

Now I got it!

Thank you very much!

January 15, 2013, 23:05	Calculation of Epsilon in LES model	#1
Yuri Almeida New Member Yuri Almeida Join Date: Jan 2012 Location: Rio de Janeiro, Brazil Posts: 21 Rep Power: 14	Dears, I'm performing flame simulations using the fireFoam solver. The default set for combustion is the eddy dissipation model, which furnishes the equation below for the source term in the species transport eqn: wFuelNorm_ = C_rho_turbulence_.epsilon() /max(turbulence_.k(),dimensionedScalar("SMALL",dimV elocitydimVelocity,SMALL)) min(YFuel , YO2/s.value()); However, the turbulence model utilized is LES OneEqEddy, which does not directly calculate the epsilon. So, where the epsilon calculation is performed? What expression is used for it? Thanks in advance, Yuri

January 17, 2013, 12:55		#2
gregor Member Gregor Olenik Join Date: Jun 2009 Location: http://greole.github.io/ Posts: 89 Rep Power: 17	Hi Yuri, the oneEqEddy model is derived from the GenEddyVisc model. In GenEddyVisc.H you find that epsilon is computed as Code: ce_k()sqrt(k())/delta();

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
LES Filter in Smagorinsky model on inhomogenous grids	Ivan	Main CFD Forum	1	October 25, 2012 12:30
Implementation of wall model to LES	saeedi	Main CFD Forum	3	August 24, 2012 15:51
Bug in Dynamic 1Eq Eddy Viscosity LES Model Incompressible flow	harishg	OpenFOAM Bugs	5	December 2, 2008 19:18
Smagorinsky model	sylvain91	OpenFOAM Running, Solving & CFD	0	June 22, 2006 07:29
2-equation model of LES and source code	M.R.Hadian	Main CFD Forum	0	February 3, 2002 06:00

January 18, 2013, 08:43		#3
Yuri Almeida New Member Yuri Almeida Join Date: Jan 2012 Location: Rio de Janeiro, Brazil Posts: 21 Rep Power: 14	Hi Gregor, Now I got it! Thank you very much!