CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

gammaCoeff in twoPhaseEulerFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 20, 2014, 09:18
Default gammaCoeff in twoPhaseEulerFoam
  #1
New Member
 
Ramon
Join Date: Feb 2014
Location: Eindhoven
Posts: 25
Rep Power: 12
RjwV is on a distinguished road
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))
        );
Kind regards,
Ramon
RjwV is offline   Reply With Quote

Old   October 28, 2014, 06:56
Default
  #2
New Member
 
Ramon
Join Date: Feb 2014
Location: Eindhoven
Posts: 25
Rep Power: 12
RjwV is on a distinguished road
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_)
and:

Code:
volScalarField gammaCoeff
(
    3.0*(1.0 - sqr(e_))*max(sqr(alpha), residualAlpha)*gs0
   *((4.0/da)*(ThetaSqrt/sqrtPi) - divU)
);
I believe in the 2.3.x release gs0 should be replaced with gs0_ .

Kind regards,
Ramon
RjwV is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 21:21.