From 50dc19899e392fcd4603e21d32cffbe1af2bf523 Mon Sep 17 00:00:00 2001 From: Hrvoje Jasak Date: Thu, 4 Jun 2015 18:03:34 +0100 Subject: [PATCH] Alternative form of k-e coupling --- .../RAS/coupledKEpsilon/coupledKEpsilon.C | 10 ++++++++++ 1 file changed, 10 insertions(+) diff --git a/src/turbulenceModels/incompressible/RAS/coupledKEpsilon/coupledKEpsilon.C b/src/turbulenceModels/incompressible/RAS/coupledKEpsilon/coupledKEpsilon.C index c5842af24..613ed4ea5 100644 --- a/src/turbulenceModels/incompressible/RAS/coupledKEpsilon/coupledKEpsilon.C +++ b/src/turbulenceModels/incompressible/RAS/coupledKEpsilon/coupledKEpsilon.C @@ -302,6 +302,16 @@ void coupledKEpsilon::correct() keEqn.insertEquation(0, kEqn); } + // Add coupling term: C1*Cmu*(symm(grad(U))) k but with wall function + // corrections: must be calculated from G. HJ, 27/Apr/2015 + + // Add coupling term: epsilon source depends on k + // k, e sink terms cannot be changed because of boundedness + keEqn.insertEquationCoupling + ( + 1, 0, -C1_*G*epsilon_/sqr(k_) + ); + // Update source coupling: coupling terms eliminated from source keEqn.updateSourceCoupling();