CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Phoenics

Please help: Multiphase flow boundary condition

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 27, 2004, 10:11
Default Please help: Multiphase flow boundary condition
  #1
Alex Costa
Guest
 
Posts: n/a
Dear group, I am trying to understand the following PIL statement: PATCH(TOP,CELL,NX,NX,1,1,1,1,1,LSTEP) COVAL(TOP,P1,1.0,0.0);COVAL(TOP,P2,RHO2,0.0)

Let me describe my understanding: It is applied in just one cell (WHY?) During all the time steps. Applies for mass flow boundary conditions for phase 1 and phase 2. For phase 1, it sets an extra term in the pressure equation equal to T 1 (0 - P1) For phase 2, it sets an extra term equal to T RHO2 (0 - P2) All right, But what is the physical meaning of these two bondary conditions?

Futhermore, After I compile the q1 the above statement is replaced by: PATCH (TOP ,CELL ,1,1,0,0,0,0,1,10) COVAL (TOP ,R1 , 1.000000E+00, 0.000000E+00) COVAL (TOP ,R2 , 1.000000E+03, 0.000000E+00) Does it mean R2 is set to 0 and adds -T 1 R1 to the R1 equation? What is the relationship with the previous statement for pressure and again what is the physical meaning?

Any interpretation about these statements are very welcome.

Alex Costa
  Reply With Quote

Old   September 1, 2004, 13:37
Default Re: Please help: Multiphase flow boundary conditio
  #2
Mick Hughes
Guest
 
Posts: n/a
Hi Alex, I'm not an expert on two phase flows, I may be able to shed some light on your questions. Firstly library case W771. I just tried this case but I made changes as it did not converge straight out of the box.

Firstly, shorten the run time or increase the number of run time discrtisations.

Secondly, This is an old library case so I used the photon text viewer to look at the results, this needs the case to be at least 2d, so I made 2 cells in the Y-direction and stored(V1,V2), FIINIT(v1,v2)=0.

Dont forget to change the "gravity" and "top" patches from ...1,1,... to ...1,ny,....

thirdly, add CSG1='a' and IDISPA=1 to drop out intermediate PHI files.

Hopefully, the case should now run and give some physically realistic result.

About the following boundary condition : PATCH(TOP,CELL,NX,NX,1,ny,1,1,1,LSTEP) COVAL(TOP,P1,1.0,0.0);COVAL(TOP,P2,RHO2,0.0)

Density is used to multiply the coefficient so as in cell mass fractions are used C:\phoenics\D_POLIS\D_LECS\IPSA\IPSA.HTM

However, this cell is not an outlet BC, I think what is happening here is that pressure is being fixed so that a unique solution can be found. Just like flow in a closed box needs a pressure relief cell, this example is essentailly a closed box, so we need to fix the pressure somewhere to zero. This we do for both phases. The "CELL" multiplier is therefore fine for this situation.

Just some extra words, the pressure correction equation is essentially solving a Laplacian type equation. To guarantee a unique solution to this type of equation the value has to be fixed in one cell at least. For flow simulations it is not the absolute value of pressure that drives the flow field, but the pressure differences, so this is generally OK.

Hope this helps some.

  Reply With Quote

Old   September 1, 2004, 17:21
Default Re: Please help: Multiphase flow boundary conditio
  #3
mick
Guest
 
Posts: n/a
P.S. I increased LSWEEP when running the above mentioned case
  Reply With Quote

Old   September 2, 2004, 02:24
Default Re: Please help: Multiphase flow boundary conditio
  #4
Rami
Guest
 
Posts: n/a
Welcome back, Mick!
  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
rotating wall boundary condition in multiphase sdp FLUENT 0 February 16, 2009 22:12
Change wall boundary condition in multiphase? Jacques FLUENT 0 November 27, 2008 07:26
outflow boundary condition for multiphase flow deepak FLUENT 1 June 21, 2007 08:23
How to handle peoridic boundary in multiphase flow Harry Main CFD Forum 0 November 19, 2005 16:10
Mass flow boundary condition Síle FLUENT 0 June 12, 2003 07:30


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