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 CARREAU-YUSADA 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(116 3D SPARSOL Pipe Flow-3Carreau-Yusada) ************************************************************ Echo save-block settings for Group 1 save1begin This case concerns the steady laminar flow of a Carreau- Yusada non-Newtonian fluid in a circular pipe. This case is similar to library case J104 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 Carreau-Yusada fluid is given by: emu = Ft*emu0*[1+(Ft*T.G)^n]^{(m-1)/n} where G is the mean strain rate, and Ft is the temperature function defined by: Ft = exp (-beta*(T-Tr)) where beta = 0.01303111 degC^-1 and Tr = 240 degC. The rheology parameters are set to: low-shear viscosity, emu0=1659.677 Pa.s time constant T = 0.1153758 s flow-behaviour index m=0.3113008 flow-behaviour index n=0.9734119 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)=BTAU NAME(141)=TFNN ;NAME(142)=TEM1 NAME(144)=WDIS ;NAME(145)=GEN1 NAME(147)=PRPS ;NAME(148)=ENUL NAME(149) =DEN1 * Solved variables list SOLVE(P1,U1,V1,W1,TEM1) * Stored variables list STORE(DEN1,ENUL,PRPS,GEN1,WDIS,TFNN,BTAU,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 =1659.677002 ;ENULB =0. ;ENULC =0.311309 CP1 =2100. DISWAL ENUT =1.0E-10 DVO1DT =-6.8E-05 PRNDTL(TEM1)=-0.15 ************************************************************ Echo save-block settings for Group 9 save9begin REAL(AX,RAM,DD,CC,BB,TR) The required model for the apparent dynamic viscosity is: ETA = AX*Ft*{1.+(r*Ft**G)^a)}^[(n-1)/a] where Ft = EXP(-b*(T-Tr)) AX =1659.677 ! Pa.s RAM=0.1153758 ! s DD =0.9734119 CC =0.3113088 BB =0.0103111 ! deg C^-1 TR =240.0 ! deg C Carreau-Yusada model - built-in settings The built-in form of the apparent dynamic viscosity of the Carreau-Yusada fluid is given by: ETA=(ETAI+(ETA0-ETAI)*(1.+(T*G)^m)**(n-1.)/m) and so the following settings are made: ETAI = 0.0 = (ENULB) infinite-shear viscosity, ETA0 = AX = (ENULA) = 1659.677 Pa.s - low-shear viscosity n = CC = (ENULC) = 0.3113088 flow-behaviour index T = RR = (ENULD) = 0.1153758 s m = DD = (ENULE) = 0.9734119 and Ft = exp (beta*(T-Tr)) where beta= -b = -0.00202 degC^-1 (=ENULF) and Tr = TR = 240 degC (=ENULG) IENULA=4 ENUL=GRND4 ENULA=AX ENULB=0.0 ENULC=CC ENULD=RAM ENULE=DD Temperature dependence IENULB=2 STORE(TFNN) ENULG=TR ENULF=-BB save9end ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;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-04 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ************************************************************ 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. ) ************************************************************ 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,GXMONI,PLOTALL,L,T) SPEDAT(SET,RLXFAC,REFLEN,R,1.0E-02) SPEDAT(SET,MAXINC,P1,R,1.0E+07) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=3.125E+06 ;EX(U1)=0.8862 EX(V1)=1.395E-03 ;EX(W1)=1.424E-03 EX(SRM1)=772.5 ;EX(BTAU)=1.155E+05 EX(TFNN)=0.6627 ;EX(TEM1)=240. EX(WDIS)=5.499E-04 ;EX(GEN1)=3.656E+13 EX(PRPS)=42.82 ;EX(ENUL)=0.196 EX(DEN1)=807.5 ;EX(LTLS)=6.734E-07 ************************************************************ 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,N,N) 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 DISTIL=T EX(P1 )=3.125E+06 EX(U1 )=8.862E-01 EX(V1 )=1.395E-03 EX(W1 )=1.424E-03 EX(TFNN)=6.627E-01 EX(TEM1)=2.400E+02 EX(LTLS)=6.363E-07 EX(ENUL)=1.960E-01 EX(DEN1)=8.075E+02 EX(WDIS)=5.264E-04 EX(SRM1)=7.725E+02 EX(LTLS)=6.734E-07 EX(WDIS)=5.499E-04 EX(BTAU)=1.155E+05 EX(GEN1)=3.656E+13 save24end GVIEW(P,-0.241541,-0.609825,0.754832) GVIEW(UP,0.157865,0.742808,0.650627) GVIEW(WINDOW,1022,726) GVIEW(DEPTH,22.2903) 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, TYPE, INLET > OBJ, COLOR-MODE, USER > OBJ, COLOR-VAL, 192 > OBJ, WIREFRAME, YES > 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, WIREFRAME, YES > 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