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 POWER-LAW FLUID 3d PIPE FLOW
   msg Pressure (P1) profile
   pause
   da 1 p1 z 1 y 1
   col3 1
   redr
   pause
   msg press  to 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(114 3D SPARSOL Pipe Flow-1Power-Law     )
 ************************************************************
  Echo save-block settings for Group  1
  save1begin
  This case concerns the steady laminar flow of a Power-
  Law non-Newtonian fluid in a circular pipe. This case
  is similar to library case J100 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 Power-Law fluid
  is given by:
 
    emu = Ft*A*G^(n-1)
 
  where G is the mean strain rate, and Ft is the temperature
  function defined by:
 
    Ft = exp (-beta*(T-Tr))
 
  where beta = 0.01037222 degC^-1 and Tr = 240 degC.
 
  The rheology parameters are set to:
 
   consistency index,   A=2855.322 Pa.s
   flow-behaviour index n=0.4413787
 
  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)=WDIS ;NAME(142)=TEM1
 NAME(144)=TFNN ;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,TFNN,WDIS,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 =2855.322021 ;ENULB =0.441379 ;ENULC =1600.
 CP1 =2100.
 DISWAL
 ENUT =1.0E-10
 DVO1DT =-6.8E-05
 PRNDTL(TEM1)=-0.15
 ************************************************************
  Echo save-block settings for Group  9
  save9begin
real(AA,BB,CC,TR,E0TR)
E0TR=1600.
AA=2855.322 ; CC=0.01037222 ; BB=0.4413787
TR=240.
 
     Power-law model - built-in settings
IENULA=0
ENUL=GRND4
ENULA=AA      ! Consistency index
ENULB=BB      ! Flow-behaviour index
ENULC=E0TR    ! Zero-shear dynamic viscosity
     Temperature dependence
IENULB=2
STORE(TFNN)
ENULG=TR      ! Reference temperature
ENULF=-CC      !
  save9end
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(P1)=0. ;FIINIT(U1)=1.337
 FIINIT(WDIS)=0.1 ;FIINIT(TEM1)=240.
 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,REFVEL,R,1.)
 SPEDAT(SET,RLXFAC,REFLEN,R,1.0E-02)
 SPEDAT(SET,MAXINC,U1,R,100.)
 SPEDAT(SET,MAXINC,V1,R,100.)
 SPEDAT(SET,MAXINC,W1,R,100.)
 SPEDAT(SET,MAXINC,TEM1,R,1000.)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=3.118E+06 ;EX(U1)=0.8865
 EX(V1)=1.447E-03 ;EX(W1)=1.474E-03
 EX(SRM1)=763.099976 ;EX(BTAU)=3.118E+09
 EX(WDIS)=5.499E-04 ;EX(TEM1)=240.
 EX(TFNN)=0.6627 ;EX(GEN1)=4.57E+13
 EX(PRPS)=42.82 ;EX(ENUL)=0.1684
 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.118E+06
EX(U1  )= 8.865E-01
EX(V1  )= 1.447E-03
EX(W1  )= 1.474E-03
EX(TEM1)= 2.400E+02
EX(ENUL)= 1.684E-01
EX(DEN1)= 8.075E+02
EX(LTLS)= 6.363E-07
EX(TFNN)= 6.627E-01
EX(WDIS)= 5.499E-04
EX(LTLS)= 6.734E-07
EX(SRM1)= 7.631E+02
EX(BTAU)= 3.118E+09
EX(GEN1)= 4.570E+13
  save24end
 
 GVIEW(P,-0.328126,-0.943846,0.038572)
 GVIEW(UP,0.048123,0.024078,0.998551)
 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