|
[Sponsors] |
August 17, 2015, 12:51 |
Errors in twoPhaseEulerFoam's Kinetic Theory
|
#1 |
New Member
Ramon
Join Date: Feb 2014
Location: Eindhoven
Posts: 25
Rep Power: 12 |
Hello everyone,
I believe I have found two errors in the Kinetic Theory of Granular Flow as implemented in OpenFOAM. I am currently working with v. 2.3.0, but counts for all versions up to 2.4.x. Number one; the frictional pressure and viscosity have not been added to the solids pressure and solids shear viscosity before solving ThetaEqn. The code snippets as posted below here define nut and PsCoeff, which are directly used in ThetaEqn. In literature, usually the frictional pressure and viscosity are added before solving the granular energy equation. Code:
nut_ = viscosityModel_->nu(alpha, Theta_, gs0_, rho, da, e_); Code:
volScalarField PsCoeff ( granularPressureModel_->granularPressureCoeff ( alpha, gs0_, rho, e_ ) ); Code:
if (alpha1[celli] > alphaMax.value() - 5e-2) { nuf[celli] = 0.5*pf[celli]*sin(phi_.value()) /( sqrt(1.0/6.0*(sqr(D[celli].xx() - D[celli].yy()) + sqr(D[celli].yy() - D[celli].zz()) + sqr(D[celli].zz() - D[celli].xx())) + sqr(D[celli].xy()) + sqr(D[celli].xz()) + sqr(D[celli].yz())) + I2Dsmall ); } Kind regards, Ramon |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Building OpenFOAM1.7.0 from source | ata | OpenFOAM Installation | 46 | March 6, 2022 13:21 |
How to write k and epsilon before the abnormal end | xiuying | OpenFOAM Running, Solving & CFD | 8 | August 27, 2013 15:33 |
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 | bookie56 | OpenFOAM Installation | 8 | August 13, 2011 04:03 |
IcoFoam parallel woes | msrinath80 | OpenFOAM Running, Solving & CFD | 9 | July 22, 2007 02:58 |
Could anybody help me see this error and give help | liugx212 | OpenFOAM Running, Solving & CFD | 3 | January 4, 2006 18:07 |