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

How could I modify the pEqn.H by adding ajustPhi() in twoPhaseEulerFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 2, 2013, 11:32
Default How could I modify the pEqn.H by adding ajustPhi() in twoPhaseEulerFoam
  #1
Member
 
Charlie
Join Date: Dec 2010
Location: 415 Kinross Dr. Newark, DE 19711
Posts: 78
Rep Power: 6
cheng1988sjtu is on a distinguished road
Hi Foamers,

When I use zeroGradient pressure condition at the wall, I could see pressure sigularity at the pReff(), and I searched in the forum, and found that that's due to that fact that there is no adjustPhi after the buoyancy term.

So, I added the following lines before phi = alphaf*phia + betaf*phib;

adjustPhi(phia, Ua, p);
adjustPhi(phib, Ub, p);

however, after I added this lines, I could see that my Dispersed volume fraction is decreasing though slowly.

Could anyone give me a hint on how adjustPhi() should work?

Thanks

Zhen Cheng
cheng1988sjtu is offline   Reply With Quote

Old   April 2, 2013, 15:36
Default
  #2
Member
 
Charlie
Join Date: Dec 2010
Location: 415 Kinross Dr. Newark, DE 19711
Posts: 78
Rep Power: 6
cheng1988sjtu is on a distinguished road
Never mind, I figured out that the mass reduction is not due to the adjustPhi, and FYI, this approach do eliminate the pressure singularity at the zeroGradient boundary (p).

Charlie
cheng1988sjtu is offline   Reply With Quote

Old   April 3, 2013, 11:05
Default
  #3
Senior Member
 
chegdan's Avatar
 
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 546
Rep Power: 18
chegdan will become famous soon enough
Greetings Zhen Cheng,

Just curious about your findings,

Quote:
Never mind, I figured out that the mass reduction is not due to the adjustPh
What was the source of mass reduction in your simulations? Just in case others run into the same issue.

Dan
__________________
Dan

Find me on twitter @dancombest and LinkedIn
chegdan is offline   Reply With Quote

Old   March 3, 2014, 10:13
Default
  #4
New Member
 
Josebus van Zyl
Join Date: Jan 2013
Location: Stellenbosch
Posts: 1
Rep Power: 0
Popmech is on a distinguished road
Hi, I am using a modified version of the twoPhaseEulerFoam solver to model a mixing vessel with settling particles. I have a zero-gradient pressure condition on all surfaces including the free-surface which is represented by a symmertryPlane type boundary condition. I suffer from a similar problem with a singularity occurring at the pReffCell in the simulation which causes the volume fractions to be solved incorrectly. I have tried implementing the adjustPhi call as described in this thread, though this causes an error stating there is no outflow boundary to adjust. Does someone know how to work around this problem?

Kind Regards, JVZ
Popmech 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
how to set periodic boundary conditions Ganesh FLUENT 13 January 22, 2014 05:11
alphaEqn.H in twoPhaseEulerFoam cheng1988sjtu OpenFOAM Bugs 14 September 18, 2011 13:46
Help with chtMultiRegionFoam jbvw96 OpenFOAM Running, Solving & CFD 2 December 26, 2010 18:16
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55


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