|
[Sponsors] |
How could I modify the pEqn.H by adding ajustPhi() in twoPhaseEulerFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 2, 2013, 11:32 |
How could I modify the pEqn.H by adding ajustPhi() in twoPhaseEulerFoam
|
#1 |
Member
Charlie
Join Date: Dec 2010
Location: USA
Posts: 85
Rep Power: 15 |
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 |
|
April 2, 2013, 15:36 |
|
#2 |
Member
Charlie
Join Date: Dec 2010
Location: USA
Posts: 85
Rep Power: 15 |
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 |
|
April 3, 2013, 11:05 |
|
#3 | |
Senior Member
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 621
Rep Power: 0 |
Greetings Zhen Cheng,
Just curious about your findings, Quote:
Dan |
||
March 3, 2014, 09:13 |
|
#4 |
New Member
Josebus van Zyl
Join Date: Jan 2013
Location: Stellenbosch
Posts: 1
Rep Power: 0 |
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 |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
how to set periodic boundary conditions | Ganesh | FLUENT | 15 | November 18, 2020 06:09 |
alphaEqn.H in twoPhaseEulerFoam | cheng1988sjtu | OpenFOAM Bugs | 15 | May 1, 2016 16:12 |
Help with chtMultiRegionFoam | jbvw96 | OpenFOAM Running, Solving & CFD | 2 | December 26, 2010 17:16 |
fluent add additional zones for the mesh file | SSL | FLUENT | 2 | January 26, 2008 11:55 |