View Issue Details

IDProjectCategoryView StatusLast Update
0001857OpenFOAMBugpublic2016-02-03 21:43
Reporterprojectionist Assigned Tohenry  
PrioritylowSeverityminorReproducibilityalways
Status closedResolutionfixed 
PlatformGNU/LinuxOSUbuntuOS Version14.04
Product Versiondev 
Summary0001857: possible bug in mixtureKEpsilon turbulence model
DescriptionThe mixtureKEpsilon turbulence model refers in its header file to the publications of A. Behzadi, R.I. Issa and H. Rusche as well as R.T. Lahey.
Unfortunately, I have only access to the publication of A. Behzadi, R.I. Issa and H. Rusche, so this report might be unfounded.
However, if I follow the derivation of the factor beta from the formulas present in the paper of A. Behzadi, R.I. Issa and H. Rusche, I find that the implemented formula in mixtureKEpsilon.C is off by a factor of 2.

However, there might be additional information in the paper of Lahey to which I have no access.
Steps To ReproduceThe following is implemented:

volScalarField beta
    (
        (6*this->Cmu_/(4*sqrt(3.0/2.0)))
       *fluid.drag(gas).K()/liquid.rho()
       *(liquidTurbulence.k_/liquidTurbulence.epsilon_)
    );

whereas, I obtain the following formula:

volScalarField beta
    (
        (6*this->Cmu_/(2.0*sqrt(3.0/2.0)))
       *fluid.drag(gas).K()/liquid.rho()
       *(liquidTurbulence.k_/liquidTurbulence.epsilon_)
    );

Notice, the factor of 2.0 in the denominator of the first term instead of 4 as in the upper code part.

Furthermore, the formulation of the factor 6/(2.0*sqrt(3.0/2.0)) - or 6/(4*sqrt(3.0/2.0)) as in the original code - seems strangely complicated.
Additional InformationThe reported issue is also present in the sources of OF-2.3.x. Other versions were not checked.
TagsNo tags attached.

Relationships

related to 0001916 closedhenry epsilonm seems not to converge in $FOAM_TUTORIAL/multiphase/twoPhaseEulerFoam/RAS/bubbleColumn 

Activities

henry

2015-09-28 16:57

manager   ~0005377

Many of Lahey's papers and reports are available on-line and there is reasonable correspondence between the formulation of the standard parts of the models in them. It would help a lot if you could check all the sources rather than have me spend time digging up the papers and studying the code again.

henry

2016-02-03 21:43

manager   ~0005909

Orphaned report

Issue History

Date Modified Username Field Change
2015-09-28 16:50 projectionist New Issue
2015-09-28 16:57 henry Note Added: 0005377
2015-11-29 19:29 wyldckat Relationship added related to 0001916
2016-02-03 21:43 henry Note Added: 0005909
2016-02-03 21:43 henry Status new => closed
2016-02-03 21:43 henry Assigned To => henry
2016-02-03 21:43 henry Resolution open => fixed