TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2020, Date 13/01/21 CPVNAM=VDI; SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings PHOTON USE AUTOPLOT file phida 3 clear msg CROSS-FLUID 3d PIPE FLOW msg Pressure (P1) profile pause da 1 p1 z 1 y 1 col3 1 redr pause msg pressto continue pause clear msg Pressure (u1) profile pause da 1 u1 z 1 y 1 colf 1 redr pause msg press to end end END_USE ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(117 3D SPARSOL Pipe Flow-4 Cross ) ************************************************************ Echo save-block settings for Group 1 save1begin This case concerns the steady laminar flow of a Cross non-Newtonian fluid in a circular pipe. This case is similar to library case J102 except the simulation is carried out using SPARSOL on a 3d Cartesian mesh assuming 1/4 symmetry. The flow conditions and geometry are similar, and can be summarised as follows: Pipe diameter = 3mm Pipe length = 150mm Inlet velocity = 1.337m/s Density = 735 kg/m^3 Specific heat = 2.1 kJ/kg.K Thermal conductivity = 0.15 W/m.K The apparent dynamic viscosity of the Cross fluid is given by: emu = Ft*emu0/[1+(Ft*T.G)^n] where G is the mean strain rate, and Ft is the temperature function defined by: Ft = exp (beta/(T+T0)) where beta = 5795.964 K and T0 = TEMP0 = 273.15K. The rheology parameters are set to: low-shear viscosity, emu0=0.02398645 Pa.s time constant T = 1.370996E-6 s flow-behaviour index n=0.7174. The case should be run in parallel mode on 4 processors rather than serial mode for the most rapid convergence. save1end ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,150,13,13,1.0E-05) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(139)=SRM1 ;NAME(140)=STRS NAME(141)=TFNN ;NAME(142)=TEM1 NAME(144)=WDIS ;NAME(145)=GEN1 NAME(147)=PRPS ;NAME(148)=ENUL NAME(149)=DEN1 ;NAME(150)=BTAU * Solved variables list SOLVE(P1,U1,V1,W1,TEM1) * Stored variables list STORE(BTAU,DEN1,ENUL,PRPS,GEN1,WDIS,TFNN,STRS) STORE(SRM1) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(TEM1,Y,Y,Y,N,N,Y) ************************************************************ Group 8. Terms & Devices TERMS (U1 ,Y,Y,Y,N,Y,Y) TERMS (V1 ,Y,Y,Y,N,Y,Y) TERMS (W1 ,Y,Y,Y,N,Y,Y) TERMS (TEM1,N,Y,Y,N,Y,Y) ADDDIF = T NEWENL = T ************************************************************ Group 9. Properties PRESS0 =1.01325E+05 ;TEMP0 =273. * Domain material index is 51 signifying: * Freezing water at 1 atm SETPRPS(1, 51) RHO1 =735. ENUL = GRND4 ENULA =0.023986 ;ENULB =0. ;ENULC =0.717393 CP1 =2100. DISWAL ENUT =1.0E-10 DVO1DT =-6.8E-05 PRNDTL(TEM1)=-0.15 ************************************************************ Echo save-block settings for Group 9 save9begin REAL(ETA0,AA,CC,TB,TAU) The required model for the apparent dynamic viscosity is: ETA = ETA0/(1.+(ETA0*G/TAU)^(1-C)) where ETA0 = A*EXP(Tb/(T+T0)) ETA0=1928.6723 ! Pa.s CC=0.2826074 TB=5795.984-TEMP0 ! degC AA=0.02398645 ! Pa.s TAU=17495.76 ! Pa.s Cross model - built-in settings The built-in form of the apparent dynamic viscosity of the Cross fluid is given by: emu = Ft*(emui+(emu0-emui)/[1+(Ft*T.G)^n]) and so the following settings are made: emui = 0.0 (ENULB) infinite-shear viscosity, emu0 = A = (ENULB) = 0.02398645 Pa.s - low-shear viscosity n = 1.-C = (ENULC) = n =0.7174 flow-behaviour index T = A/TAU = (ENULD) = 1.370996E-6 s and Ft = exp (beta/(T+T0)) where beta = 5795.964 K (=ENULF)and T0 = TEMP0 = 273.15K. ENUL=GRND4; IENULA=2 ENULA=AA ENULB=0.0 ENULC=1.-CC ENULD=AA/TAU STORE(TFNN) Temperature dependence ! IENULB > 0 ENUL is multipled by F(T) < 0 As above, but time constant (ENULD) also multiplied by F(T) IENULB=-4 ENULF=TB+TEMP0 save9end ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;FIINIT(U1)=1.337 FIINIT(TEM1)=240. ;FIINIT(WDIS)=0.1 FIINIT(PRPS)=-1. No PATCHes used for this Group INIADD = F NAMFI ='phi' ************************************************************ Group 12. Convection and diffusion adjustments No PATCHes used for this Group ************************************************************ Group 13. Boundary & Special Sources No PATCHes used for this Group BUOYA =0. ; BUOYB =0. ; BUOYC =-9.81 BUOYE =240. EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 10000 RESFAC =1.0E-05 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(TEM1)=50 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,0.99 ) RELAX(U1 ,LINRLX,0.5 ) RELAX(V1 ,LINRLX,0.5 ) RELAX(W1 ,LINRLX,0.5 ) RELAX(TEM1,LINRLX,0.25 ) RELAX(LTLS,LINRLX,1. ) RELAX(ENUL,LINRLX,0.3 ) ************************************************************ Group 18. Limits VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-204.75 ************************************************************ Group 19. EARTH Calls To GROUND Station PARSOL = F CONWIZ = T ISG62 = 1 SPEDAT(SET,CUTCELL,FRCMN,R,1.0E-03) SPEDAT(SET,CUTCELL,FRCMX,R,0.999) SPEDAT(SET,OUTPUT,NOFIELD,L,T) SPEDAT(SET,RLXFAC,REFLEN,R,1.0E-02) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,MAXINC,P1,R,1.0E+08) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=2.551E+06 ;EX(U1)=0.8861 EX(V1)=8.748E-04 ;EX(W1)=1.002E-03 EX(SRM1)=663.700012 ;EX(STRS)=6.101 EX(TFNN)=5.347E+04 ;EX(TEM1)=240. EX(WDIS)=5.499E-04 ;EX(GEN1)=7.436E+13 EX(LTLS)=6.734E-07 ;EX(PRPS)=42.82 EX(ENUL)=0.2316 ;EX(DEN1)=807.5 EX(BTAU)=1.359 ************************************************************ Group 21. Print-out of Variables OUTPUT(P1 ,N,N,Y,Y,Y,Y) OUTPUT(U1 ,N,N,Y,Y,Y,Y) OUTPUT(V1 ,N,N,Y,Y,Y,Y) OUTPUT(W1 ,N,N,Y,Y,Y,Y) OUTPUT(TEM1,N,N,Y,Y,Y,Y) OUTPUT(GEN1,N,N,Y,N,N,N) OUTPUT(ENUL,N,N,Y,N,Y,Y) OUTPUT(DEN1,N,N,Y,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 75 ;IYMON = 1 ;IZMON = 7 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts ************************************************************ Echo save-block settings for Group 24 save24begin EX(P1 )=2.551E+06 EX(U1 )=8.861E-01 EX(V1 )=8.748E-04 EX(W1 )=1.002E-03 EX(SRM1)=6.637E+02 EX(STRS)=6.101E+00 EX(TFNN)=5.347E+04 EX(TEM1)=2.400E+02 EX(GEN1)=7.436E+13 EX(ENUL)=2.316E-01 EX(DEN1)=8.075E+02 EX(PRPS)=4.282E+01 EX(LTLS)=6.734E-07 EX(WDIS)=5.499E-04 EX(BTAU)=1.359E+00 save24end GVIEW(P,-0.557753,-0.83,-3.369044E-03) GVIEW(UP,-0.012671,4.456025E-03,0.99991) GVIEW(DEPTH,3.0E+04) GVIEW(VDIS,0.072051) GVIEW(CENTRE,0.075,2.0E-03,2.0E-03) > DOM, SIZE, 1.500000E-01, 4.000000E-03, 4.000000E-03 > DOM, MONIT, 7.450000E-02, 1.250000E-04, 1.625000E-03 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, INCREMENT, 1.000000E-01, 1.000000E-01, 1.000000E-01 > GRID, RSET_X_1, 150, 1.000000E+00 > GRID, RSET_Y_1, 12, 1.000000E+00 > GRID, RSET_Y_2, 1, 1.000000E+00 > GRID, RSET_Z_1, 12, 1.000000E+00 > GRID, RSET_Z_2, 1, 1.000000E+00 > DOM, T_AMBIENT, 2.400000E+02 > DOM, INI_AMB, YES > DOM, INI_BUOY, YES > OBJ, NAME, WEST0 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 3.000000E-03, 3.000000E-03 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, INLET > OBJ, COLOR-MODE, USER > OBJ, COLOR-VAL, 192 > OBJ, PRESSURE, P_AMBIENT > OBJ, VELOCITY, 1.337 ,0. ,0. > OBJ, TEMPERATURE, 240. > OBJ, NAME, EAST0 > OBJ, POSITION, 1.500000E-01, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 3.000000E-03, 3.000000E-03 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, OUTLET > OBJ, COLOR-MODE, USER > OBJ, COLOR-VAL, 214 > OBJ, PRESSURE, P_AMBIENT > OBJ, TEMPERATURE, T_AMBIENT > OBJ, COEFFICIENT, 1000. > OBJ, NAME, CBLOCK > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, 3.000000E-03, 3.000000E-03 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, corner > OBJ, ROTATION24, 10 > OBJ, VISIBLE, NO > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 125,Plasterboard > OBJ, FIXED_TMP, 0. ,240. > OBJ, NAME, NBLOCK > OBJ, POSITION, 0.000000E+00, 3.000000E-03, 0.000000E+00 > OBJ, SIZE, TO_END, TO_END, 3.000000E-03 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, ROTATION24, 10 > OBJ, VISIBLE, NO > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 125,Plasterboard > OBJ, FIXED_TMP, 0. ,240. > OBJ, NAME, HBLOCK > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 3.000000E-03 > OBJ, SIZE, TO_END, TO_END, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, ROTATION24, 10 > OBJ, VISIBLE, NO > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 125,Plasterboard > OBJ, FIXED_TMP, 0. ,240. STOP