CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

linearization of k-e source terms

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   May 19, 2006, 05:44
Default linearization of k-e source terms
  #1
khaiching
Guest
 
Posts: n/a
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-
  Reply With Quote

Old   May 19, 2006, 14:50
Default Re: linearization of k-e source terms
  #2
ganesh
Guest
 
Posts: n/a
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

  Reply With Quote

Old   May 20, 2006, 00:08
Default Re: linearization of k-e source terms
  #3
khaiching
Guest
 
Posts: n/a
Hi Ganesh, my email id is: ngkhaiching2000@yahoo.com. Thanks for your help..
  Reply With Quote

Old   July 10, 2009, 06:17
Default
  #4
Senior Member
 
Aram Amouzandeh
Join Date: Mar 2009
Location: Vienna, Vienna, Austria
Posts: 186
Rep Power: 8
mabinty is on a distinguished road
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
mabinty is offline   Reply With Quote

Old   July 10, 2009, 09:22
Default sample of fortran program
  #5
New Member
 
Asghar
Join Date: Apr 2009
Posts: 16
Rep Power: 8
bohluly is on a distinguished road
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
bohluly is offline   Reply With Quote

Old   July 11, 2009, 04:26
Default
  #6
Member
 
Mohammad Reza Hadian
Join Date: Mar 2009
Location: Yazd, Iran
Posts: 50
Rep Power: 8
hadian is on a distinguished road
would you please put the full title and address of the paper here?
hadian is offline   Reply With Quote

Old   May 23, 2013, 02:46
Default
  #7
Member
 
Join Date: Aug 2011
Posts: 74
Rep Power: 5
idefix is on a distinguished road
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
idefix is offline   Reply With Quote

Old   June 6, 2013, 11:00
Default Prof. Dick's paper
  #8
New Member
 
ali
Join Date: Feb 2012
Posts: 8
Rep Power: 5
cfdman_aero is on a distinguished road
Quote:
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
cfdman_aero is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
swak4foam building problem GGerber OpenFOAM Installation 54 April 24, 2015 16:02
pisoFoam compiling error with OF 1.7.1 on MAC OSX Greg Givogue OpenFOAM Programming & Development 3 March 4, 2011 18:18
Using source terms jsm Main CFD Forum 4 August 20, 2009 06:44
DxFoam reader update hjasak OpenFOAM Post-Processing 69 April 24, 2008 01:24
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51


All times are GMT -4. The time now is 17:56.