View Issue Details
ID | Project | Category | View Status | Date Submitted | Last Update |
---|---|---|---|---|---|
0001857 | OpenFOAM | Bug | public | 2015-09-28 16:50 | 2016-02-03 21:43 |
Reporter | projectionist | Assigned To | henry | ||
Priority | low | Severity | minor | Reproducibility | always |
Status | closed | Resolution | fixed | ||
Platform | GNU/Linux | OS | Ubuntu | OS Version | 14.04 |
Product Version | dev | ||||
Summary | 0001857: possible bug in mixtureKEpsilon turbulence model | ||||
Description | The 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 Reproduce | The 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 Information | The reported issue is also present in the sources of OF-2.3.x. Other versions were not checked. | ||||
Tags | No tags attached. | ||||
|
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. |
|
Orphaned report |
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 |