|
[Sponsors] |
Please help: Multiphase flow boundary condition |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 27, 2004, 11:11 |
Please help: Multiphase flow boundary condition
|
#1 |
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 |
|
September 1, 2004, 14:37 |
Re: Please help: Multiphase flow boundary conditio
|
#2 |
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. |
|
September 1, 2004, 18:21 |
Re: Please help: Multiphase flow boundary conditio
|
#3 |
Guest
Posts: n/a
|
P.S. I increased LSWEEP when running the above mentioned case
|
|
September 2, 2004, 03:24 |
Re: Please help: Multiphase flow boundary conditio
|
#4 |
Guest
Posts: n/a
|
Welcome back, Mick!
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 09: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 08:30 |