|
[Sponsors] |
October 20, 2014, 09:18 |
gammaCoeff in twoPhaseEulerFoam
|
#1 |
New Member
Ramon
Join Date: Feb 2014
Location: Eindhoven
Posts: 25
Rep Power: 12 |
Dear Foamers,
I would like to adjust the kineticTheoryModel.C file, more specifically the dissipation rate gammaCoeff. Currently the simplified equation 3.24 on page 50 of Van Wachem's thesis has been implemented. However, I would like to implement equation 3.23, where the divergence term of the solids velocity has not been neglected. Could I get some hints on how I could program this into the code properly? Currently I have this, but it is not working. Code:
// Dissipation (Eq. 3.24, p.50) /* volScalarField gammaCoeff ( 12.0*(1.0 - sqr(e_)) *max(sqr(alpha), residualAlpha_) *gs0*(1.0/da)*ThetaSqrt/sqrtPi );*/ // Added eqn. 3.23 from Van Wachem (RJWV) volScalarField gammaCoeff ( 3.0*(1.0 - sqr(e_))*max(sqr(alpha), residualAlpha_)*gs0 *((4.0/da)*(ThetaSqrt/sqrtPi) - fvm::div(U)) ); Ramon |
|
October 28, 2014, 06:56 |
|
#2 |
New Member
Ramon
Join Date: Feb 2014
Location: Eindhoven
Posts: 25
Rep Power: 12 |
Just for the sake of completion, it seems I have managed to tackle the problem myself. It was not that difficult after all, but hopefully it can aid someone in the future.
in kineticTheoryModel.C the following changes can be made: Code:
volScalarField divU = fvc::div(this->U_) Code:
volScalarField gammaCoeff ( 3.0*(1.0 - sqr(e_))*max(sqr(alpha), residualAlpha)*gs0 *((4.0/da)*(ThetaSqrt/sqrtPi) - divU) ); Kind regards, Ramon |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Sliding mesh for twoPhaseEulerFoam | Lapo | OpenFOAM Programming & Development | 4 | November 25, 2019 15:23 |
Is twoPhaseEulerFoam applicable to 3D cases / delivering erroneous results? | ThomasV | OpenFOAM | 0 | November 11, 2013 09:10 |
Something wrong in UEqns.H within twoPhaseEulerFoam | cheng1988sjtu | OpenFOAM | 2 | June 24, 2011 11:48 |
twoPhaseEulerFoam | freemankofi | OpenFOAM | 0 | May 23, 2011 17:24 |
problems in Two Phase flow using twoPhaseEulerFoam with OpenFoam 1.6 | raagh77 | OpenFOAM Running, Solving & CFD | 0 | March 6, 2010 06:11 |