TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2022, Date 05/07/22
CPVNAM=VDI; SPPNAM=Core
Echo DISPLAY / USE settings
DISPLAY
Library case p133
2-phase Rising liquid drop
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(2-Phase Rising Bubble )
************************************************************
Group 2. Transience
STEADY=F
* Set overall time and no. of steps
RSET(U,0.,3.,3000)
* Modify regions
RSET(T,1,3000,1.)
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,320,640,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(146)=SURT ;NAME(147)=PRPS
NAME(148)=DEN1 ;NAME(149)=VFOL
NAME(150) =SURN
* Solved variables list
SOLVE(P1,U1,V1,VFOL,SURN)
* Stored variables list
STORE(DEN1,PRPS,SURT)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
TERMS (VFOL,N,N,N,N,Y,N)
TERMS (SURN,N,N,N,N,Y,N)
GALA = T
RUPLIM =0.999 ;RLOLIM =1.0E-03
************************************************************
Group 9. Properties
PRESS0 =1.01325E+05 ;TEMP0 =273.149994
RHO1 =1000.
ENUL =0.01006
CP1 =4181.799805
ENUT =0.
DVO1DT =1.18E-04
* List of user-defined materials to be read by EARTH
MATFLG=T;IMAT=2
* Name
*Ind. Dens. Viscos. Spec.heat Conduct. Expans. Compr.
*
34 1.0 0.1 1005.0 0.0258 3.41E-3 0.0
*
76 1000.0 1.006E-2 4181.8 0.597 1.18E-4
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=0. ;FIINIT(PRPS)=-1.
FIINIT(VFOL)=1. ;FIINIT(SURN)=1.
No PATCHes used for this Group
INIADD = F
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
PATCH(BUOYANCY, PHASEM, 0, 0, 0, 0, 0, 0, 1, 3000)
COVAL(BUOYANCY, U1, FIXFLU, GRND2)
COVAL(BUOYANCY, V1, FIXFLU, GRND2)
PATCH(STENS, PHASEM, 0, 0, 0, 0, 0, 0, 1, 3000)
COVAL(STENS, U1, FIXFLU, GRND3)
COVAL(STENS, V1, FIXFLU, GRND3)
BUOYA =0. ; BUOYB =-0.98 ; BUOYC =0.
BUOYD =1000.
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 30
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=200 ;LITER(VFOL)=1
LITER(SURN)=1
ENDIT(P1)=1.0E-04
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
RELAX(VFOL,LINRLX,1. )
RELAX(SURN,LINRLX,1. )
************************************************************
Group 18. Limits
VARMAX(SURN)=1. ;VARMIN(SURN)=0.
************************************************************
Group 19. EARTH Calls To GROUND Station
SURF = T
PARSOL = F
IPRPSA = 76 ;IPRPSB = 34
ISG62 = 1
ISURFA = 5
BTATHC =3.5
ISRFTA = 1
SURFTA =1.96 ;SURFTB=0. ;SURFTC=0.
CSG10 ='q1'
SPEDAT(SET,OUTPUT,NOFIELD,L,T)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=612. ;EX(U1)=0.04869
EX(V1)=0.1446 ;EX(SURT)=0.9019
EX(PRPS)=144.399994 ;EX(DEN1)=902.
EX(VFOL)=0.9019 ;EX(SURN)=0.9019
************************************************************
Group 21. Print-out of Variables
OUTPUT(VFOL,Y,Y,Y,N,Y,Y)
OUTPUT(SURN,Y,Y,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 158 ;IYMON = 158 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
IDISPA = 25 ;IDISPB = 0 ;IDISPC = 0
CSG1 ='A'
GVIEW(P,0.,0.,1.)
GVIEW(UP,0.,1.,0.)
GVIEW(DEPTH,6624.743652)
GVIEW(VDIS,1.175755)
GVIEW(CENTRE,0.5,1.,0.5)
> DOM, SIZE, 1.000000E+00, 2.000000E+00, 1.000000E+00
> DOM, MONIT, 4.921875E-01, 4.921875E-01, 5.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID, RSET_X_1, 80, 1.000000E+00
> GRID, RSET_X_2, 160, 1.000000E+00
> GRID, RSET_X_3, 80, 1.000000E+00
> GRID, RSET_Y_1, 80, 1.000000E+00
> GRID, RSET_Y_2, 160, 1.000000E+00
> GRID, RSET_Y_3, 400, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> DOM, INI_AMB, YES
> DOM, INI_BUOY, YES
> OBJ, NAME, BOTTOMWALL
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, TO_END, 0.000000E+00, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, box
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, NAME, TOPWALL
> OBJ, POSITION, 0.000000E+00, AT_END, 0.000000E+00
> OBJ, SIZE, TO_END, 0.000000E+00, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, box
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, NAME, BUBBLE
> OBJ, POS-CEN, 5.000000E-01, 5.000000E-01, 5.000000E-01
> OBJ, SIZE, 5.000000E-01, 5.000000E-01, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cylinder
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 34, AIR_AT_20_DEG_ (LIGHT)
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
************************************************************
STOP