CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Phoenics (https://www.cfd-online.com/Forums/phoenics/)

 9527 March 27, 2009 22:34

where is the superman ---about the IPSA-full

dear all:
i want to simulate the two-flow fluid with IPSA-full model .but
(1)how can i get the value of P2,and how can i modify the P2 value .for example make p2=p2+constant
(2)in the IPSA-full model ,how does it trace the interface of two-flow.could you tell me where the associated subroutine is

 9527 April 9, 2009 21:14

no one?？？/？/???/

 phoenics_cfd April 27, 2009 07:57

P2 in IPSA

If you enter STORE(P2) in the Q1 input file , then PHOENICS will compute the pressure acting on the 2nd-phase momentum momentum to be P1 + P2, where P1 is the pressure acting on phase 1. Therefore, P2 is the EXTRA pressure experienced by phase 2.

It is the responsibility of the user to specify the formula for P2 in Group 19 Section of GROUND, or alternatively by use of the INFORM facility within the Q1 input file.

A GROUND-coding example of setting a P2-formula is given in Group 19 Section 4 of the file GREX3.FOR, as follows:

C
IF(STORE(P2)) THEN
IF(STORE(VPOR)) THEN
C.... Extra second-phase pressure proportional to volume fraction
C divided by porosity.
CALL FN15(P2,R2,VPOR,0.0,PHS2PA)
ELSE
IF(PHS2P) THEN
C.... Extra second-phase pressure proportional to:
C a*(r2-r2crit)**b if r2 gt r2crit.
L0R2=L0F(R2)
L0P2=L0F(P2)
RLX=ABS(DTFALS(P2))
DO 1941 I=1,NXNY
IF(F(L0R2+I).GT.PHS2PA) THEN
GP2=PHS2PB*(F(L0R2+I)-PHS2PA)**PHS2PC
ELSE
GP2=0.0
ENDIF
F(L0P2+I)=RLX*GP2+(1.0-RLX)*F(L0P2+I)
1941 CONTINUE
ELSE
C.... Extra second-phase pressure proportional to volume fraction
CALL FN2(P2,R2,0.0,PHS2PA)
ENDIF
ENDIF
ENDIF
ENDIF

 9527 April 28, 2009 03:04