CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

compressibleTwoPhaseEulerFoam pEqn missing code for particle pressure

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 15, 2012, 11:28
Default compressibleTwoPhaseEulerFoam pEqn missing code for particle pressure
  #1
Member
 
bruce
Join Date: May 2009
Location: Germany
Posts: 42
Rep Power: 8
bruce is on a distinguished road
Hello,

I have found that compressible version of twoPhaseEulerFoam has missing code as shown below, which accounts for particle pressure in solid momentum equation. Which is being accounted in pressure loop.


Code:
    if (g0.value() > 0.0)
    {
        phiDraga -= ppMagf*fvc::snGrad(alpha)*mesh.magSf();
    }

    if (kineticTheory.on())
    {
        phiDraga -= rUaAf*fvc::snGrad(kineticTheory.pa()/rhoa)*mesh.magSf();
    }
The above code is not included in pEqn of compressibleTwoPhaseEulerFoam ! Is this really missing or it is being implemented in different way?
bruce is offline   Reply With Quote

Old   June 23, 2012, 13:41
Default
  #2
Member
 
bruce
Join Date: May 2009
Location: Germany
Posts: 42
Rep Power: 8
bruce is on a distinguished road
Resolved by commit 5e27d9890131ea032a51bf2d0d0d3fe6ae724673
bruce is offline   Reply With Quote

Old   August 14, 2012, 07:14
Default
  #3
New Member
 
Join Date: Feb 2012
Posts: 3
Rep Power: 5
stern is on a distinguished road
Dear Bruce,

from your post I found that compressibleTwoPhaseEulerFoam is not working well.
Can you inform me at which location the missing a "ppDrag" term has to be included in pEqn.H of compressibleTwoPhaseEulerFoam.

A second question is: Is it possible to simulate an isothermal problem with compressibleTwoPhaseEulerFoam? If not, did you know if twoPhaseEulerFoam is now working well?

Best regards,
Nicole
stern is offline   Reply With Quote

Old   August 14, 2012, 08:54
Default
  #4
Member
 
bruce
Join Date: May 2009
Location: Germany
Posts: 42
Rep Power: 8
bruce is on a distinguished road
Hi,

Which version of the solver are you using? There are many recent updates to compressibleTwoPhaseEulerFoam in git version 2.1.x. I have no issues using the latest version.

compressibleTwoPhaseEulerFoam can solve temperature for both phases. But twoPhaseEulerFoam has no such facility. You have two choices while using compressibleTwoPhaseEulerFoam. Set isothermal BC everywhere. Or, comment out temperature solution in the solver and make it cold flow solver.

Best,
Bruce
bruce is offline   Reply With Quote

Old   August 14, 2012, 09:25
Default
  #5
New Member
 
Join Date: Feb 2012
Posts: 3
Rep Power: 5
stern is on a distinguished road
Hi Bruce,

now I use the latest release version 2.1.1

With twoPhaseEulerFoam of previous versions I have bad experiences, because in the freeboard region of a fluidization bed downward high velocities of the solid arise. Because in the release notes of version 2.1.1 nothing about this solver exists, I am unsure if the module twoPhaseEulerFoam was corrected. That is why I decide to work with compressibleTwoPhaseEulerFoam, although my problem is isothermal.

Now, my problem is that after certain time the field "pa" is blowing up in certain cells at the baffle boundary situated in the bed.

Actually, I am looking for a bug-free 2 phase Euler-Euler code, because I am trying since 2 years to use OpenFoam for fluidization problems and after each test I realize that the cod has bugs.

Many thanks for your help,
Nicole
stern is offline   Reply With Quote

Old   April 8, 2013, 04:08
Default
  #6
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 667
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by bruce View Post
Hi,

Which version of the solver are you using? There are many recent updates to compressibleTwoPhaseEulerFoam in git version 2.1.x. I have no issues using the latest version.

compressibleTwoPhaseEulerFoam can solve temperature for both phases. But twoPhaseEulerFoam has no such facility. You have two choices while using compressibleTwoPhaseEulerFoam. Set isothermal BC everywhere. Or, comment out temperature solution in the solver and make it cold flow solver.

Best,
Bruce
Hi,

Bruce,can compressibletwophaseeulerfoam works on a tet mesh?
I have tried twophaseeulerfoam but failed.

regards,
sharonyue is offline   Reply With Quote

Old   July 31, 2013, 09:20
Default is there any tutorials or solved examples on compressibleTwoPhaseEulerFoam
  #7
Member
 
Anil Kunwar
Join Date: Jun 2013
Posts: 52
Rep Power: 3
Annier is an unknown quantity at this point
i need to learn about compressibleTwoPhaseEulerFoam. are there any tutorials or explanation of examples related to this solver in OpenFoam?
Annier is offline   Reply With Quote

Old   July 31, 2013, 09:45
Default
  #8
Member
 
bruce
Join Date: May 2009
Location: Germany
Posts: 42
Rep Power: 8
bruce is on a distinguished road
annier,

there are 3 tutorials being shipped with OF. see for example in,

https://github.com/OpenFOAM/OpenFOAM...haseEulerFoam/

but no explanations of problem setup.
bruce is offline   Reply With Quote

Old   August 1, 2013, 08:30
Default thank you
  #9
Member
 
Anil Kunwar
Join Date: Jun 2013
Posts: 52
Rep Power: 3
Annier is an unknown quantity at this point
Dear Bruce,
thank you very much for your kindness.
Annier is offline   Reply With Quote

Old   August 1, 2013, 08:47
Default my query on the tutorials
  #10
Member
 
Anil Kunwar
Join Date: Jun 2013
Posts: 52
Rep Power: 3
Annier is an unknown quantity at this point
dear bruce,
the tutorials bubbleColumn and mixerVessel2D are also present in twoPhaseEulerFoam. what codes or files are added in those tutorials of compressibleTwoPhaseEulerFoam inorder to account for the compressibility? it would be easier for me to learn about these tutorials if you would explain me the difference between these tutorials with the same names. thanking you in advance.
Annier is offline   Reply With Quote

Old   August 1, 2013, 09:32
Default
  #11
Member
 
bruce
Join Date: May 2009
Location: Germany
Posts: 42
Rep Power: 8
bruce is on a distinguished road
Hi, Sorry I ve been away using these solvers for long time. but first you could carefully observe the differences in the problem setup for both cases and come up with general understanding. and later on dig into the differences in the code itself. may be it is easier to start with "fluidisedBed" problem.

Regards,
bruce is offline   Reply With Quote

Old   August 1, 2013, 22:53
Default thank you
  #12
Member
 
Anil Kunwar
Join Date: Jun 2013
Posts: 52
Rep Power: 3
Annier is an unknown quantity at this point
dear bruce,
thank you very much.
Annier is offline   Reply With Quote

Old   August 11, 2013, 16:55
Default energy equation in compressibleTwoPhaseEulerFoam
  #13
New Member
 
Josefine Wilms
Join Date: Sep 2009
Posts: 9
Rep Power: 7
Fine is on a distinguished road
Hi

Does anybody know what the following two lines in the EEqn.H file of compressibleTwoPhaseEulerFoam mean?

Code:
+ heatTransferCoeff*he1/Cpv1/rho1
      - fvm::Sp(heatTransferCoeff/Cpv1/rho1, he1)
It appears as though the same term is added and then subtracted?

Regards
Josefine
Fine is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Pressure Inlet VS velocity Inlet difference Mohsin FLUENT 6 January 15, 2015 18:15
Pressure BC for combustion chamber Giuki FLUENT 1 July 19, 2011 11:35
Setup/monitor points of pressure and force coefficients siw CFX 3 October 22, 2010 06:07
FLOW AROUND A PLATE_NEGATIVE ABSOLUTE PRESSURE???? tania FLUENT 11 March 23, 2004 09:51
Hydrostatic pressure in 2-phase flow modeling (long) DS & HB Main CFD Forum 0 January 8, 2000 16:00


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