TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2022, Date 21/03/22
CPVNAM=VDI; SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(J118-1dFDTurb.PipeFlow Modified PowerLaw)
************************************************************
Echo save-block settings for Group 1
save1begin
This case concerns the steady fully-developed turbulent
flow of a pseudo-plastic fluid with the apparent viscosity
computed from the following modified power-law fluid:
emua = K*G^(n-1) for G > Glim
emua = [ a G + b ] for G ≤ Glim
where G is the mean strain rate, n is the power law index
and K is the consistency. These rheology parameters are
set to: K=54.36 n=0.6 and Glim=258.4
The turbulence is simulated by using the low-Reynolds-
number two-equation Lam-Bremhorst k-e model.
The pipe diameter is 0.5m and the fluid density is
1000 kg/m^3. The single-slab solver is used to compute the
mass flow rate for a specified axial pressure gradient
of 15.904 kNs/m^3.
The main results are written to the inforout file, as
follows:
mass flow rate = 451.1 kg/s
bulk velocity = 36.1 m/s
Friction Re = 321
Power-law Re = 2788
Dodge-Metzner ReDM = 5841
Friction factor f = 0.0061
For the same axial pressure drop, the DNS results of
Singh et al (2017) produce the following:
bulk velocity = 34.6 m/s,
ReDM = 5498 and f = 0.00664.
The Dodge-Metzner empirical correlations indicate that:
ReDM=5443 and f = 0.00667.
J.Singh,M. Rudman & H.M.Blackburn
The influence of shear-dependent rheology on
turbulent pipe flow
J.Fluid Mech. (2017), vol. 822, pp. 848–879.
save1end
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,1,400,1,1.0E-12,1.0E-15,1.0E-12)
* Cylindrical-polar grid
CARTES=F
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(136)=BTAU ;NAME(137)=ENUL
NAME(138)=AREH ;NAME(139)=MUL
NAME(140)=GEN1 ;NAME(141)=SRM1
NAME(143)=YPLS ;NAME(144)=STRS
NAME(145)=SKIN ;NAME(146)=GENK
NAME(147)=WDIS ;NAME(148)=DEN1
NAME(149)=EL1 ;NAME(150)=ENUT
* Solved variables list
SOLVE(W1)
* Stored variables list
STORE(ENUT,EL1,DEN1,WDIS,GENK,SKIN,STRS,YPLS)
STORE(SRM1,GEN1,MUL,AREH,ENUL,BTAU)
TURMOD(KEMODL-LOWRE)
************************************************************
Group 8. Terms & Devices
TERMS (W1 ,N,N,Y,Y,Y,Y)
NEWENL = T
USOURC = T
************************************************************
Group 9. Properties
RHO1 =1000.
ENUL = GRND4
ENULA =54.362885 ;ENULB =0.6 ;ENULC =258.400238 ;ENULD =0.
ENULE =0. ;ENULF =0. ;ENULG =0.
IENULA = 10 ;IENULB = 0
CP1 =1.
* Non-default TERMS setting for KE
TERMS (KE ,N,N,Y,Y,Y,N)
* Non-default TERMS setting for EP
TERMS (EP ,N,N,Y,Y,Y,N)
PRT(EP)=1.314
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(W1)=34.34 ;FIINIT(KE)=3.169
FIINIT(EP)=230.482407 ;FIINIT(ENUL)=6.9E-03
FIINIT(LTLS)=1. ;FIINIT(YPLS)=0.05
FIINIT(WDIS)=1.0E-02 ;FIINIT(ENUT)=0.0155
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(PRESDIF, HIGH, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(PRESDIF, W1, FIXFLU, 1.6E+04)
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 10000
RESFAC =1.0E-08
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(W1 ,LINRLX,0.1 )
RELAX(KE ,LINRLX,0.1 )
RELAX(EP ,FALSDT,0.1 )
RELAX(LTLS,LINRLX,1. )
KELIN = 1
************************************************************
Group 18. Limits
VARMAX(W1)=1.0E+10 ;VARMIN(W1)=1.0E-10
************************************************************
Group 19. EARTH Calls To GROUND Station
GENK = T
FDFSOL = T
PARSOL = F
IENUTA = 3
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,GXMONI,CLASSIC,L,T)
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
OUTPUT(GEN1,Y,N,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 6 ;IZMON = 1
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
GVIEW(P,-0.019838,0.252223,-0.967466)
GVIEW(UP,2.140179E-03,0.967665,0.252231)
GVIEW(VDIS,0.311528)
GVIEW(CENTRE,0.024958,0.25,0.5)
> DOM, SIZE, 1.000000E-01, 5.000000E-01, 1.000000E+00
> DOM, MONIT, 5.000000E-02, 5.158764E-02, 5.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E-01
> DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID, RSET_X_1, 1, 1.000000E+00
> GRID, RSET_Y_1, 399,-1.020004E+00,G
> GRID, RSET_Y_2, 1, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> DOM, T_AMBIENT, 0.000000E+00
> OBJ, NAME, NULL1
> OBJ, POSITION, 0.000000E+00, 4.998485E-01, 0.000000E+00
> OBJ, SIZE, TO_END, 0.000000E+00, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polwire
> OBJ, TYPE, NULL
> OBJ, NAME, PLAT2
> OBJ, POSITION, 0.000000E+00, AT_END, 0.000000E+00
> OBJ, SIZE, TO_END, 0.000000E+00, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polbox
> OBJ, TYPE, PLATE
DISTIL=T
EX(W1 )=1.747E+01;EX(KE )=6.294E+00
EX(EP )=3.104E+02
EX(BTAU)=2.612E+03;EX(ENUL)=4.771E-03
EX(AREH)=3.125E-05;EX(MUL )=4.771E+00
EX(GEN1)=9.300E+05;EX(SRM1)=7.799E+02
EX(LTLS)=1.182E-02;EX(YPLS)=5.000E-02
EX(STRS)=9.930E-03;EX(SKIN)=1.060E+00
EX(GENK)=1.551E+02;EX(WDIS)=5.247E-02
EX(DEN1)=1.000E+03;EX(EL1 )=2.580E-02
EX(ENUT)=2.069E-02
STOP