Hi: APpreciate if anyone could enlighten me on this:?

(1) What are the standard techniques designed to linearize the source terms for the k-e equations? Could anyone please advise on the available references?

(2) If I would like to include the buoyant terms into the k-e source term, what would be the correct implementation? Any idea on how to make the buoyant term more "implicit", i.e. transfer that to the main diagonal of the coefficient matrix?

TQ

-khai ching-

Dear khaiching,

There is a paper by Prof. Dick's group on linearisation of source terms for k-epsilon and k-w equations. I can mail it to you if you give your mail id

Regards,

Ganesh

Hi Ganesh, my email id is: ngkhaiching2000@yahoo.com. Thanks for your help..

[mabinty]

Hi Ganesh!!

I would kindly ask you to send me the paper on linearisation of source terms of k-epsilon equations by Erik Dick s group:

aram.amouzandeh@tuwien.ac.at

Thx in advance,
Aram

[bohluly]

source_kEx(ij)= gen ! explicit source of k
source_kIm(ij)= -Cmu(ij)*k1/ (nut(ij)+.000000001) ! implicit source of k

cturb1 = 1.44d0 ; cturb2 = 1.92d0

k_t0 = k_t(ij)+.0000001
ep_0 = ep_(ij)

source_eEx(ij)= cturb1 * ep_0 / k_t0 * gen ! explicit source of e
source_eIm(ij)= -cturb2 * ep_0 / k_t0 ! implicit source of e

.
.
.
! source term k
do ij=1,no_nodes
k_t(ij) = ( k_t(ij) + source_kEx(ij) *dt) / (1.-dt*source_kIm(ij))
enddo
! source term e
do ij=1,no_nodes
ep_(ij) = ( ep_(ij) + source_eEx(ij) *dt) / (1.-dt*source_eIm(ij))
enddo

[hadian]

would you please put the full title and address of the paper here?

[idefix]

Hello together,

even the thread is old I hope somebody is reading it

I want to add a source term in the epsilon equation in the k-epsilon-model

The term is:

4 * k * C_mu * HLP

with:
HLP = (du/dx)^2 + (dv/dy)^2 + (dw/dz)^2 + 1/2 * (du/dy + dv/dx)^2 + 1/2 * (du/dz + dw/dx)^2 + 1/2 * (dv/dz + dw/dy)^2

After adding this term to the epsilon-equation it looks like:
+ fvm::SuSp(4.0*k_*Cmu_*magSqr(symm(fvc::grad(U_)))/(epsilon_+const), epsilon_)

with const = dimensiond scalar (1e-05) to avoid dividing by 0

Is it correct to implement this term like this?
I´ve read: http://www.cfd-online.com/Wiki/Sourc..._linearization
but I don´t know how to divide my term in a constant and a linear part.

Can anybody help?
thanks a lot

greetings
Idefix

[cfdman_aero]

Originally Posted by ganesh
;42758

Dear khaiching,

There is a paper by Prof. Dick's group on linearisation of source terms for k-epsilon and k-w equations. I can mail it to you if you give your mail id

Regards,

Ganesh

Dear ganesh,

It would be my pleasure if you send me that paper via e-mail. My e-mail: ali_or_aero@yahoo.com