Differences between production term in k-epsilon!

TomWaits · May 13, 2024, 14:12

Reading about the production term for k-epsilon here:
https://www.cfd-online.com/Wiki/Stan...-epsilon_model
it is given as
$P_k=2\mu_t S_{ij}S_{ij}$

but on the NASA website here:
https://turbmodels.larc.nasa.gov/ke-chien.html
a much more complicated form is given as:
$P=\tau_{ij}\frac{du_i}{dx_j}$

where
$\tau_{ij}=\mu_t(2S_{ij}-2/3\frac{du_k}{dx_k}\delta_{ij})-2/3\rho k\delta_{ij}$

I have coded both up and get slightly different results for the production term in my k-epsilon model! Where do these two forms come from and why are they different?

Many thanks in advance,
Tom Waits

sbaffini · May 15, 2024, 08:03

Quote:

Originally Posted by TomWaits

Reading about the production term for k-epsilon here:
https://www.cfd-online.com/Wiki/Stan...-epsilon_model
it is given as
$P_k=2\mu_t S_{ij}S_{ij}$

but on the NASA website here:
https://turbmodels.larc.nasa.gov/ke-chien.html
a much more complicated form is given as:
$P=\tau_{ij}\frac{du_i}{dx_j}$

where
$\tau_{ij}=\mu_t(2S_{ij}-2/3\frac{du_k}{dx_k}\delta_{ij})-2/3\rho k\delta_{ij}$

I have coded both up and get slightly different results for the production term in my k-epsilon model! Where do these two forms come from and why are they different?

Many thanks in advance,
Tom Waits

The second one is the only correct form of the production term in k equation. The first one is, nonetheless, a very common approximation (fully correct only for soleinodal velocity fields).

The approximation is certainly more stable, that's one of the reasons it got large acceptance.

FMDenaro · May 15, 2024, 09:26

The isotropic terms could be absorbed by the “pressure”

May 13, 2024, 14:12	Differences between production term in k-epsilon!	#1
TomWaits Member Tom Waits Join Date: Aug 2018 Posts: 42 Rep Power: 8	Reading about the production term for k-epsilon here: https://www.cfd-online.com/Wiki/Stan...-epsilon_model it is given as $P_k=2\mu_t S_{ij}S_{ij}$ but on the NASA website here: https://turbmodels.larc.nasa.gov/ke-chien.html a much more complicated form is given as: $P=\tau_{ij}\frac{du_i}{dx_j}$ where $\tau_{ij}=\mu_t(2S_{ij}-2/3\frac{du_k}{dx_k}\delta_{ij})-2/3\rho k\delta_{ij}$ I have coded both up and get slightly different results for the production term in my k-epsilon model! Where do these two forms come from and why are they different? Many thanks in advance, Tom Waits

May 15, 2024, 09:26		#3
FMDenaro Senior Member Filippo Maria Denaro Join Date: Jul 2010 Posts: 6,897 Rep Power: 73	The isotropic terms could be absorbed by the “pressure” TomWaits likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
How to add Source term (2) for PYROLYSIS - reactingOneDim	Kummi	OpenFOAM	10	April 12, 2021 22:11
Adding a source term to k and epsilon equation for buoyancy	rajeshkunwar	OpenFOAM Running, Solving & CFD	2	April 19, 2019 08:41
Production term in Finite Volumes Schemes	Ivan	Main CFD Forum	10	October 17, 2012 18:39
Extra term in Epsilon equation	Valdemir G. Ferreira	Main CFD Forum	0	April 30, 2000 09:57
bouyancy term in epsilon equation	Michael	Main CFD Forum	1	June 25, 1999 11:20