linearization of ke source terms
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 
Re: linearization of ke source terms
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 
Re: linearization of ke source terms
Hi Ganesh, my email id is: ngkhaiching2000@yahoo.com. Thanks for your help..:)

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 
sample of fortran program
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 
would you please put the full title and address of the paper here?

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 
Prof. Dick's paper
Quote:
It would be my pleasure if you send me that paper via email. My email: ali_or_aero@yahoo.com 
All times are GMT 4. The time now is 18:28. 