CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

add a phase diffusion term in k-epsilon turbulence model

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

LinkBack Thread Tools Search this Thread Display Modes
Old   February 12, 2020, 12:23
Question add a phase diffusion term in k-epsilon turbulence model
Join Date: Mar 2019
Posts: 35
Rep Power: 2 is on a distinguished road
Hi guys,

I would like to add a diffusion term into the k-epsilon equations so I wrote the following code:

- fvm::div(fvc::snGrad(alpha)*mesh.magSf()*fvc::inte rpolate(nut*rho), k_)

tmp<fvScalarMatrix> epsEqn
fvm::ddt(alpha, rho, epsilon_)
+ fvm::div(alphaRhoPhi, epsilon_)
- fvm::laplacian(alpha*rho*DepsilonEff(), epsilon_)
- fvm::SuSp(((2.0/3.0)*C1_ - C3_)*alpha()*rho()*divU, epsilon_)
- fvm::Sp(C2_*alpha()*rho()*epsilon_()/k_(), epsilon_)
- fvm::div(fvc::snGrad(alpha)*mesh.magSf()*fvc::inte rpolate(nut*rho), epsilon_)

However it cannot be compiled. Please help me to check where I need to modified. Thanks a lot!
Attached Images
File Type: jpg Picture3.jpg (16.9 KB, 3 views) is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
k and Epsilon is not converging in kEpsilon Turbulence model Raza Javed OpenFOAM Running, Solving & CFD 1 June 21, 2019 03:28
Define different turbulence model for each phase (multiphaseEulerFoam) remidemol OpenFOAM Running, Solving & CFD 2 February 27, 2018 18:57
Error in Two phase (condensation) modeling adilsyyed CFX 15 June 24, 2015 20:42
Adding an energy source term to a phase with reactions lxlxlxt Fluent Multiphase 0 February 10, 2014 21:42
Discussion: Reason of Turbulence!! Wen Long Main CFD Forum 3 May 15, 2009 10:52

All times are GMT -4. The time now is 09:01.