TALK=F;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2009, Date 17/02/09
CPVNAM=VDI;SPPNAM=Core
************************************************************
Group 1. Run Title
TEXT(USP labirinth with conjugate heat transfer. )
************************************************************
title
DISPLAY
This case solves a two-dimensional steady hydrodynamics
problem with conjugate heat transfer in labirinth channel.
The analytical solution is absent. The solution can be compared
with solution of Structured PHOENICS.
Cases changes inlet velocity as Uin=0.1*caseno
Domain material is air.
Upper material is epoxy. Lower material is gold.
|------------------------------------|
| |
|------| EPOXY |
In Flow -> |----| | |
| | | |
| | -------| |
| |-------| | |
| | | |
| | | |
| | ---------| |
| |--------| | |
| | | |
| | | |
| GOLD | ---------| |
| |--------| | |
| | | |
| | | |
| | | |
|------------------------------------|
ENDDIS
REAL(Uin)
INTEGER(caseno)
caseno=1
Uin=0.1*caseno
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
boolean(AGG)
AGG=F
mesg(Do you want to use AGG (y) or CGG (n)? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
AGG=T
endif
if(AGG)then
RSET(M,10,10,1)
else
INTEGER(NCUBE,NCHAN)
NCUBE = 80
NCHAN = NCUBE/20
RSET(M,NCUBE,NCUBE,1)
endif
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(148) = PRPS; NAME(147) = KOND
NAME(149) = TEM1
* Solved variables list
SOLVE(P1 ,U1 ,V1, TEM1)
* Stored variables list
STORE(PRPS,KOND)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
SOLUTN(TEM1 ,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
USP = T ; USPDBG= F
UTCPLT = F ; USPIMB= F
UAUTO = AGG ; UGTGRD= F
USPVTK = T ; USPTER= F
CELLST = 10 ;FACEST = 10
if(AGG)then
USCANX = T ; USCANY= T
USCANZ = T
MXLEV = 0 ;MYLEV = 0 ;MZLEV = 0
DOMAT = -1
MINPRP = -1 ; MAXPRP= 100
UMXLEV = 4 ; UMXCLS= 10000
UNOLAY = 2 ; UIOLEV= 0
UWALEV = 0 ; UGRITR= 2
UMXCRS = 16
else
MXLEV = 4 ;MYLEV = 4 ;MZLEV = 4
DOMAT = 0
MINPRP = -1 ;MAXPRP = 250
endif
************************************************************
Group 9. Properties
PRESS0 =1.0E+05 ;TEMP0 =273.
* Domain material index is 0 signifying:
* Air at 20 deg C, 1 atm, treated as incompressible
SETPRPS(1, 0)
ENUT =0.
DVO1DT =3.41E-03
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS)=-1.
No PATCHes used for this Group
if(AGG)then
else
PATCH(B1,INIVAL, 4*NX/5+NCHAN+1,NX,1,NY,1,NZ,1,1)
COVAL(B1, PRPS, 0.0, 104)
PATCH(B2,INIVAL, 3*NX/5+NCHAN+1,4*NX/5+NCHAN,NY/5+NCHAN+1,NY,1,NZ,1,1)
COVAL(B2, PRPS, 0.0, 104)
PATCH(B3,INIVAL, 2*NX/5+NCHAN+1,3*NX/5+NCHAN,2*NY/5+NCHAN+1,NY,1,NZ,1,1)
COVAL(B3, PRPS, 0.0, 104)
PATCH(B4,INIVAL, NX/5+NCHAN+1,2*NX/5+NCHAN,3*NY/5+NCHAN+1,NY,1,NZ,1,1)
COVAL(B4, PRPS, 0.0, 104)
PATCH(B5,INIVAL, 1,NX/5+NCHAN,4*NY/5+NCHAN+1,NY,1,NZ,1,1)
COVAL(B5, PRPS, 0.0, 104)
PATCH(B6,INIVAL, 1,4*NX/5,1,NY/5,1,NZ,1,1)
COVAL(B6, PRPS, 0.0, 107)
PATCH(B7,INIVAL, 1,3*NX/5,NY/5+1,2*NY/5,1,NZ,1,1)
COVAL(B7, PRPS, 0.0, 107)
PATCH(B8,INIVAL, 1,2*NX/5,2*NY/5+1,3*NY/5,1,NZ,1,1)
COVAL(B8, PRPS, 0.0, 107)
PATCH(B9,INIVAL, 1,NX/5,3*NY/5+1,4*NY/5,1,NZ,1,1)
COVAL(B9, PRPS, 0.0, 107)
endif
INIADD = F
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
No PATCHes used for this Group
if(AGG)then
PATCH(AA11,CELL, 1,8,1,1,1,1,1,1)
COVAL(AA11,TEM1, FIXVAL, 100.)
PATCH(AA12,CELL, 1,NX,NY,NY,1,1,1,1)
COVAL(AA12,TEM1, FIXVAL, 0.)
else
PATCH(AA10,NORTH, 1,NX,NY,NY,1,NZ,1,1)
COVAL(AA10,TEM1, FIXVAL, 0.)
PATCH(AA11,SOUTH, 1,4*NX/5,1,1,1,NZ,1,1)
COVAL(AA11,TEM1, FIXVAL, 100.)
endif
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 3000
RESFAC = 1.000000E-07
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,0.7 )
RELAX(U1 ,FALSDT, 1.000000E-01)
RELAX(V1 ,FALSDT, 1.000000E-01)
mesg(Do you want to use collocated arrangement (y) or staggered one (n)? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
SPEDAT(SET,USP,METHOD,I,1)
RELAX(P1 ,LINRLX, 3.000000E-01)
RELAX(U1 ,FALSDT, 5.000000E-02)
RELAX(V1 ,FALSDT, 5.000000E-02)
mesg(Do you want to use SIMPLEST (y) or SIMPLE (n)? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
SPEDAT(SET,USP,SIMPLEST,L,T)
endif
endif
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
PARSOL = F
ISG62 = 0
ISG60 = 1
SPEDAT(SET,USP,METHOD,I,1)
SPEDAT(SET,USP,SIMPLEST,L,T)
SPEDAT(SET,USPIO,VERTCENT,L,F)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,SOLPHI,P1,I,24)
SPEDAT(SET,SOLPHI,U1,I,0)
SPEDAT(SET,PREPHI,U1,I,0)
SPEDAT(SET,SOLPHI,V1,I,0)
SPEDAT(SET,PREPHI,V1,I,0)
SPEDAT(SET,SOLPHI,TEM1,I,24)
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
if(AGG)then
IXMON = 5 ;IYMON = 4 ;IZMON = 1
else
IXMON = 50 ;IYMON = 42 ;IZMON = 1
endif
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
mesg(Do you want to view results in the centres of cells? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
SPEDAT(SET,USPIO,VERTCENT,L,F)
endif
************************************************************
Group 24. Dumps For Restarts
if(AGG)then
GVIEW(P,-0.042817,0.054706,0.997584)
GVIEW(UP,0.017703,0.998385,-0.05399)
> DOM, SIZE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, MONIT, 4.750000E-01, 3.750000E-01, 5.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID, RSET_X_1, 10, 1.000000E+00
> GRID, RSET_Y_1, 10, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> OBJ, NAME, IN
> OBJ, POSITION, 0.000000E+00, 8.000000E-01, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 5.000000E-02, 1.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, TYPE, INLET
> OBJ, GRID, NO
> OBJ, PRESSURE, 0.
> OBJ, VELOCITY, 0.1 ,0. ,0.
> OBJ, NAME, OUT
> OBJ, POSITION, 8.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 5.000000E-02, 0.000000E+00, 1.000000E+00
> OBJ, GEOMETRY, cube12t
> OBJ, TYPE, OUTLET
> OBJ, GRID, NO
> OBJ, PRESSURE, 0.
> OBJ, COEFFICIENT, 1000.
> OBJ, NAME, B3
> OBJ, POSITION, 0.000000E+00, 8.500000E-01, 0.000000E+00
> OBJ, SIZE, TO_END, TO_END, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cube14
> OBJ, GRID, NO
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 104,Epoxy
> OBJ, NAME, B4
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 2.000000E-01, 8.000000E-01, 1.000000E+00
> OBJ, GEOMETRY, cube14
> OBJ, GRID, NO
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 107, Gold
> OBJ, NAME, B5
> OBJ, POSITION, 2.500000E-01, 6.000000E-01, 0.000000E+00
> OBJ, SIZE, TO_END, 2.500000E-01, 1.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cube14
> OBJ, GRID, NO
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 104,Epoxy
> OBJ, NAME, B6
> OBJ, POSITION, 2.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 2.500000E-01, 5.500000E-01, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cube14
> OBJ, GRID, NO
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 107, Gold
> OBJ, NAME, B7
> OBJ, POSITION, 5.000000E-01, 4.000000E-01, 0.000000E+00
> OBJ, SIZE, TO_END, 2.000000E-01, 1.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cube14
> OBJ, GRID, NO
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 104,Epoxy
> OBJ, NAME, B8
> OBJ, POSITION, 4.500000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 2.000000E-01, 3.500000E-01, 1.000000E+00
> OBJ, GEOMETRY, cube14
> OBJ, GRID, NO
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 107, Gold
> OBJ, NAME, B9
> OBJ, POSITION, 7.000000E-01, 2.000000E-01, 0.000000E+00
> OBJ, SIZE, TO_END, 2.000000E-01, 1.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cube14
> OBJ, GRID, NO
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 104,Epoxy
> OBJ, NAME, B10
> OBJ, POSITION, 6.500000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.500000E-01, 1.500000E-01, 1.000000E+00
> OBJ, GEOMETRY, cube14
> OBJ, GRID, NO
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 107, Gold
> OBJ, NAME, B11
> OBJ, POSITION, 8.500000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, TO_END, 2.000000E-01, 1.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cube14
> OBJ, GRID, NO
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 104,Epoxy
else
GVIEW(P,0.000000E+00,0.000000E+00,-1.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, MONIT, 4.950000E-01, 4.150000E-01, 5.000000E-02
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> OBJ, NAME, IN
> OBJ, POSITION, 0.000000E+00, 8.000000E-01, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 5.000000E-02, 1.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, 1.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, NAME, OUT
> OBJ, POSITION, 8.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 5.000000E-02, 0.000000E+00, 1.000000E+00
> OBJ, GEOMETRY, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.000000E+03
endif
STOP