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

twoPhaseEulerFoam variables

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 31, 2011, 09:51
Default twoPhaseEulerFoam variables
  #1
Member
 
Costas K
Join Date: Nov 2009
Location: UK
Posts: 51
Rep Power: 7
thekay is on a distinguished road
Send a message via ICQ to thekay Send a message via MSN to thekay
Hi all

I was going through the code of twoPhaseEulerFoam and I have some questions about the variables in a part of alphaEqn.H and I'd be thankful if anyone could help me identify their use.

What is phic used for? It is initialised as
Code:
surfaceScalarField phic = phi;
What about alphaf? It is defined as
Code:
surfaceScalarField alphaf = fvc::interpolate(alpha);
And finally I cannot understand what phipp represents and how it is used.

Code:
surfaceScalarField phipp = ppMagf*fvc::snGrad(alpha)*mesh.magSf();
Cheers all
thekay is offline   Reply With Quote

Old   April 1, 2011, 02:45
Default
  #2
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Hi,

phipp is a term used only if the ppMagf model is used. In such case, the momentum equation for phase a has an additional pressure gradient grad(pa). If you proceed as in the derivation of phase a "pressure equation", but keep that equation written in terms of alpha, understand the origin of phipp.

Hint: grad(pa) = G(alpha) * grad(alpha), where G(alpha) = ppMagf = d(pa)/d(alpha).

Best,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   November 26, 2011, 00:48
Default
  #3
New Member
 
Kshitij kunte
Join Date: Jun 2011
Posts: 16
Rep Power: 6
k.kshitij is on a distinguished road
Hi Alberto,

I have one question about the ppmagf's usage, according to my understanding the ppmagf model is implemented when the kinetic theory is off, ie when the value of g0 > 0.
Then why does its definition involve the kinetic theory particle pressure "pa". Or I'm wrong in getting this.
I'm a newbee to C++ and have not understood the code yet would you please clarify.

Thanks in advance.
k.kshitij is offline   Reply With Quote

Old   November 26, 2011, 22:00
Default
  #4
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
The currently implemented ppMagf does not involve the kinetic theory model. I used "pa" in comment #2 to suggest how to extend the current treatment of the particle pressure in the alpha equation only implemented when the ppMagf model is used also to the case when the kinetic theory is used.

Best,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto 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
ADDITIONAL VARIABLES LIMITERS N.P CFX 5 December 1, 2011 08:47
$FOAM_USER_APPBIN & LIBBIN variables not set rassilon OpenFOAM Installation 2 February 16, 2010 21:54
Thread Variables in Case File Bernhardt FLUENT 1 June 15, 2007 01:30
PHI file structure Eugene Phoenics 9 November 2, 2001 23:00
Saving variables on a user patch with physical coordinates Zoltan Turzo CFX 2 April 20, 2000 15:05


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