TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2019, Date 28/03/19 CPVNAM=VDI; SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings VRV USE VARIABLE VFOL CONTOUR ON MSG Press OK to animate the contour plot of volume fraction PAUSE ANIMATE MSG Press OK to animate the contour plot of temperature PAUSE VARIABLE TEM1 VARIABLE RANGE 20 60 ANIMATE ENDUSE ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(XY Water-Paraffin VOF-CICSAM x-y hot ************************************************************ Group 2. Transience STEADY=F * Set overall time and no. of steps RSET(U,0.,1.,500) * Modify regions RSET(T,1,500,1.) ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,70,70,1,1.0E-04) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(145)=SPH1 ;NAME(146)=TEM1 NAME(147)=PRPS ;NAME(148)=DEN1 NAME(149)=SURN ;NAME(150)=VFOL * Solved variables list SOLVE(P1,U1,V1,TEM1,SURN,VFOL) * Stored variables list STORE(DEN1,PRPS,SPH1) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(TEM1,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) GALA = T RUPLIM =0.99 ;RLOLIM =1.0E-02 ************************************************************ Group 9. Properties PRESS0 =1.01325E+05 ;TEMP0 =273. RHO1 =998.22998 ENUL =1.006E-06 CP1 =4181.799805 ENUT =0. DVO1DT =1.18E-04 PRNDTL(TEM1)=-0.597 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;FIINIT(TEM1)=20. 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(BUOYANCY, PHASEM, 0, 0, 0, 0, 0, 0, 1, 500) COVAL(BUOYANCY, U1, FIXFLU, GRND2) COVAL(BUOYANCY, V1, FIXFLU, GRND2) PATCH(BUOYVOF, PHASEM, 0, 0, 0, 0, 0, 0, 1, 500) COVAL(BUOYVOF, U1, FIXFLU, GRND3) COVAL(BUOYVOF, V1, FIXFLU, GRND3) BUOYA =0. ; BUOYB =-9.81 ; BUOYC =0. BUOYD =861.150024 BUOYE =20. EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 60 RESFAC =1.0E-04 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(SURN)=1 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,1. ) RELAX(SURN,LINRLX,1. ) RELAX(VFOL,LINRLX,1. ) ************************************************************ Group 18. Limits VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-204.75 VARMAX(SURN)=1. ;VARMIN(SURN)=0. ************************************************************ Group 19. EARTH Calls To GROUND Station SURF = T IPRPSA = 67 ;IPRPSB = 59 ISG62 = 0 ISURFA = 1 SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout ************************************************************ Group 21. Print-out of Variables OUTPUT(SPH1,N,N,Y,Y,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 35 ;IYMON = 22 ;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,3000.) GVIEW(VDIS,5.831162E-03) GVIEW(CENTRE,5.0E-03,5.0E-03,5.0E-03) > DOM, SIZE, 1.000000E-02, 1.000000E-02, 1.000000E-02 > DOM, MONIT, 4.928572E-03, 3.071429E-03, 5.000000E-03 > 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, 28, 1.000000E+00 > GRID, RSET_X_2, 7, 1.000000E+00 > GRID, RSET_X_3, 7, 1.000000E+00 > GRID, RSET_X_4, 28, 1.000000E+00 > GRID, RSET_Y_1, 28, 1.000000E+00 > GRID, RSET_Y_2, 14, 1.000000E+00 > GRID, RSET_Y_3, 14, 1.000000E+00 > GRID, RSET_Y_4, 14, 1.000000E+00 > GRID, RSET_Z_1, 1, 1.000000E+00 > DOM, INI_AMB, YES > DOM, INI_BUOY, YES > OBJ, NAME, SIDE1 > OBJ, POSITION, AT_END, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 1.000000E-02, 1.000000E-02 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube11 > OBJ, TYPE, PLATE > OBJ, WIREFRAME, YES > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, SIDE2 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 1.000000E-02, 1.000000E-02 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube11 > OBJ, TYPE, PLATE > OBJ, WIREFRAME, YES > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, SIDE3 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E-02, 0.000000E+00, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube11 > OBJ, TYPE, PLATE > OBJ, WIREFRAME, YES > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, OUTLH > OBJ, POSITION, 0.000000E+00, AT_END, 0.000000E+00 > OBJ, SIZE, 1.000000E-02, 0.000000E+00, 1.000000E-02 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube12t > OBJ, TYPE, OUTLET > OBJ, COLOR-MODE, DEFAULT2 > OBJ, COLOR-VAL, 214 > OBJ, OPAQUE, 100 > OBJ, WIREFRAME, YES > OBJ, PRESSURE, P_AMBIENT > OBJ, TEMPERATURE, T_AMBIENT > OBJ, COEFFICIENT, 1000. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, WATER > OBJ, POSITION, 4.000000E-03, 6.000000E-03, 0.000000E+00 > OBJ, SIZE, 2.000000E-03, 2.000000E-03, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cylinder > OBJ, ROTATION24, 13 > OBJ, TYPE, BLOCKAGE > OBJ, WIREFRAME, YES > OBJ, MATERIAL, 67, WATER at 20. deg C > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, INI_TEMP, 60. > OBJ, NAME, WEDGE1 > OBJ, POSITION, 5.000000E-03, 4.000000E-03, 0.000000E+00 > OBJ, SIZE, 1.000000E-03, 2.000000E-03, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, wedge > OBJ, ROTATION24, 9 > OBJ, TYPE, BLOCKAGE > OBJ, WIREFRAME, YES > OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, INI_TEMP, 30. > OBJ, NAME, WEDGE2 > OBJ, POSITION, 4.000000E-03, 4.000000E-03, 0.000000E+00 > OBJ, SIZE, 1.000000E-03, 2.000000E-03, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, wedge > OBJ, ROTATION24, 18 > OBJ, TYPE, BLOCKAGE > OBJ, WIREFRAME, YES > OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, INI_TEMP, 40. STOP