
[Sponsors] 
May 19, 2006, 05:44 
linearization of ke source terms

#1 
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 ke equations? Could anyone please advise on the available references? (2) If I would like to include the buoyant terms into the ke 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 

May 19, 2006, 14:50 
Re: linearization of ke source terms

#2 
Guest
Posts: n/a

Dear khaiching,
There is a paper by Prof. Dick's group on linearisation of source terms for kepsilon and kw equations. I can mail it to you if you give your mail id Regards, Ganesh 

May 20, 2006, 00:08 
Re: linearization of ke source terms

#3 
Guest
Posts: n/a

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


July 10, 2009, 06:17 

#4 
Senior Member
Aram Amouzandeh
Join Date: Mar 2009
Location: Vienna, Vienna, Austria
Posts: 186
Rep Power: 9 
Hi Ganesh!!
I would kindly ask you to send me the paper on linearisation of source terms of kepsilon equations by Erik Dick s group: aram.amouzandeh@tuwien.ac.at Thx in advance, Aram 

July 10, 2009, 09:22 
sample of fortran program

#5 
New Member
Asghar
Join Date: Apr 2009
Posts: 16
Rep Power: 9 
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 

July 11, 2009, 04:26 

#6 
Member
Mohammad Reza Hadian
Join Date: Mar 2009
Location: Yazd, Iran
Posts: 52
Rep Power: 9 
would you please put the full title and address of the paper here?


May 23, 2013, 02:46 

#7 
Member
Join Date: Aug 2011
Posts: 82
Rep Power: 6 
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 kepsilonmodel 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 epsilonequation it looks like: + fvm::SuSp(4.0*k_*Cmu_*magSqr(symm(fvc::grad(U_)))/(epsilon_+const), epsilon_) with const = dimensiond scalar (1e05) to avoid dividing by 0 Is it correct to implement this term like this? Iīve read: http://www.cfdonline.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 

June 6, 2013, 11:00 
Prof. Dick's paper

#8  
New Member
ali
Join Date: Feb 2012
Posts: 8
Rep Power: 6 
Quote:
It would be my pleasure if you send me that paper via email. My email: ali_or_aero@yahoo.com 

Thread Tools  
Display Modes  


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 PostProcessing  69  April 24, 2008 01:24 
DecomposePar links against liblamso0 with OpenMPI  jens_klostermann  OpenFOAM Bugs  11  June 28, 2007 17:51 