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 p131 Liquid drop between fluids with different surface tensions ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Liquid lens ) ************************************************************ Group 2. Transience STEADY=F * Set overall time and no. of steps RSET(U,0.,2.5,500) * Modify regions ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,90,90,1) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(144)=PRPS ;NAME(145)=DEN1 NAME(146)=ENUL ;NAME(147)=SURN NAME(148)=VFOL ;NAME(149)=SRN2 NAME(150) =VFL2 * Solved variables list SOLVE(P1,U1,V1,SURN,VFOL,SRN2,VFL2) * Stored variables list STORE(ENUL,DEN1,PRPS) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) ************************************************************ Group 8. Terms & Devices TERMS (SURN,N,N,N,N,Y,N) TERMS (VFOL,N,N,N,N,Y,N) TERMS (SRN2,N,N,N,N,Y,N) TERMS (VFL2,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 =1. ENUL =1. CP1 =1005. ENUT =0. DVO1DT =3.41E-03 * List of user-defined materials to be read by EARTH MATFLG=T;IMAT=3 * Name *Ind. Dens. Viscos. Spec.heat Conduct. Expans. Compr. *67 1.0 1.0 4181.8 0.597 1.18E-4 * 68 1.0 1.0 4181.8 0.597 1.18E-4 * 69 1.0 1.0 1005.0 0.0258 3.41E-3 0.0 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;FIINIT(PRPS)=-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(STENS, PHASEM, 0, 0, 0, 0, 0, 0, 1, 500) COVAL(STENS, U1, FIXFLU, GRND3) COVAL(STENS, V1, FIXFLU, GRND3) EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 30 RESFAC =1.0E-05 ************************************************************ Group 16. Terminate Iterations LITER(SURN)=1 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,1. ) RELAX(SURN,LINRLX,1. ) RELAX(VFOL,LINRLX,1. ) RELAX(SRN2,LINRLX,1. ) RELAX(VFL2,LINRLX,1. ) ************************************************************ Group 18. Limits VARMAX(SURN)=1. ;VARMIN(SURN)=0. VARMAX(SRN2)=1. ;VARMIN(SRN2)=0. ************************************************************ Group 19. EARTH Calls To GROUND Station SURF = T PARSOL = F IPRPSA = 68 ;IPRPSB = 67 ;IPRPSC = 69 ISG62 = 1 ISURFA = 5 BTATHC =3.5 ISRFTA = 1 SURFTA =1. ;SURFTB=0. ;SURFTC=0. ISRFTB= 1 SRFTA2 =1.2 ;SRFTB2=0. ;SRFTC2=0. ISRFTC= 1 SRFTA3 =1.2 ;SRFTB3=0. ;SRFTC3=0. CSG10 ='q1' SPEDAT(SET,OUTPUT,NOFIELD,L,T) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,MAXINC,SURN,R,1.) SPEDAT(SET,MAXINC,VFOL,R,1.) SPEDAT(SET,MAXINC,SRN2,R,1.) SPEDAT(SET,MAXINC,VFL2,R,1.) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=3.056 ;EX(U1)=5.385E-04 EX(V1)=4.354E-04 ;EX(PRPS)=1.185E+05 EX(DEN1)=1. ;EX(ENUL)=1. EX(SURN)=0.358 ;EX(VFOL)=0.358 EX(SRN2)=0.2828 ;EX(VFL2)=0.2828 ************************************************************ Group 21. Print-out of Variables OUTPUT(SURN,Y,Y,Y,Y,Y,Y) OUTPUT(VFOL,Y,Y,Y,N,Y,Y) OUTPUT(SRN2,Y,Y,Y,Y,Y,Y) OUTPUT(VFL2,Y,Y,Y,N,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 55 ;IYMON = 37 ;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 = 5 ;IDISPB = 0 ;IDISPC = 0 CSG1 ='A' GVIEW(P,0.,0.,1.) GVIEW(UP,0.,1.,0.) GVIEW(DEPTH,9000.) GVIEW(VDIS,0.315715) GVIEW(CENTRE,0.249291,0.249291,0.5) > DOM, SIZE, 4.985818E-01, 4.985820E-01, 1.000000E+00 > DOM, MONIT, 3.021237E-01, 2.020194E-01, 5.000000E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02 > GRID, RSET_X_1, 18, 1.000000E+00 > GRID, RSET_X_2, 54, 1.000000E+00 > GRID, RSET_X_3, 18, 1.000000E+00 > GRID, RSET_Y_1, 18, 1.000000E+00 > GRID, RSET_Y_2, 27, 1.000000E+00 > GRID, RSET_Y_3, 27, 1.000000E+00 > GRID, RSET_Y_4, 18, 1.000000E+00 > GRID, RSET_Z_1, 1, 1.000000E+00 > DOM, INI_AMB, YES > DOM, INI_BUOY, YES > OBJ, NAME, PLAT1 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, PLATE > OBJ, WIREFRAME, YES > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, PLAT1_2 > OBJ, POSITION, AT_END, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, PLATE > OBJ, WIREFRAME, YES > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, PLAT3 > 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, WIREFRAME, YES > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, PLAT3_4 > 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, WIREFRAME, YES > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, HEAVY > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, 2.492909E-01, 1.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, WIREFRAME, YES > OBJ, MATERIAL, 68, WATER_ADIM (HEAVY) > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, DROP > OBJ, POS-CEN, 2.492909E-01, 2.492909E-01, 5.000000E-01 > OBJ, SIZE, 3.003130E-01, 3.003130E-01, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cylinder > OBJ, TYPE, BLOCKAGE > OBJ, COLOR-MODE, USER > OBJ, COLOR-VAL, 142 > OBJ, OPAQUE, 50 > OBJ, WIREFRAME, YES > OBJ, MATERIAL, 69, AIRAT20DEGC1AT (THIRD) > OBJ, TIME_LIMITS, ALWAYS_ACTIVE ************************************************************ STOP