DISPLAY
TWO TUBES IN A SEMI-RECTANGULAR CHAMBER
The natural convection around two differently heated tubes
placed in an adiabatic rectangular chamber with circular roof
is solved here by PARSOL.
The case demonstrates the PARSOL 'cut-cell' technique of
representing curvilinear shapes in a Cartesian grid.
The case considers the natural convection arising
in a 2D adiabatic rectangular enclosure. It is covered by the
roof of cylindrical shape. Two cylindrical tubes are placed at
the bottom part of the chamber. They have different temperatures
provoking the gravity-induced circulation in the surrounding air.
ENDDIS
************************************************************
Q1 created by VDI menu, Version 3.3, Date 24/07/00
CPVNAM=VDI;SPPNAM=Core
************************************************************
IRUNN = 1 ;LIBREF = 400
************************************************************
Group 1. Run Title
TEXT(Two tubes in a semi-rectangular chamber )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,25,25,1)
* Set overall domain extent:
* xulast yvlast zwlast
name
XSI= 3.200000E-02; YSI= 3.200000E-02; ZSI= 1.400000E+00
RSET(D,CHAM )
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(150) =TEM1
SOLVE(P1 ,U1 ,V1 ,TEM1)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
SETPRPS(1, 0)
RHO1 = 1.189000E+00
RHO1A = 0.000000E+00 ;RHO1B = 3.496503E-03
RHO1C = 7.142857E-01
PRESS0 = 1.000000E+05
TEMP0 = 2.730000E+02
CP1 = 1.005000E+03
ENUL = 1.544000E-05 ;ENUT = 0.000000E+00
ENULA = 1.460000E-06 ;ENULB = 1.100000E+02
ENULC = 0.000000E+00
DVO1DT = 3.410000E-03
PRNDTL(TEM1) = -2.580000E-02
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
No PATCHes used for this Group
FIINIT(TEM1)=0.5 * (20.0 + 150.0)
INIADD = F
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
PATCH (BUOYANCY,PHASEM,0,0,0,0,0,0,1,1)
COVAL (BUOYANCY,U1 , FIXFLU , GRND3 )
COVAL (BUOYANCY,V1 , FIXFLU , GRND3 )
BUOYA = 0.000000E+00 ; BUOYB =-9.810000E+00
BUOYC = 0.000000E+00
BUOYD = 3.410000E-03
BUOYE = 2.000000E+01
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 350
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
RELAX(U1 ,FALSDT, 1.000000E-01)
RELAX(V1 ,FALSDT, 1.000000E-01)
RELAX(TEM1,FALSDT, 1.000000E+01)
************************************************************
Group 18. Limits
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
VARMAX(TEM1) = 1.000000E+10 ;VARMIN(TEM1) = 2.000000E+01
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
LSG3 = T
CSG3 ='FGEM'
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
OUTPUT(P1 ,Y,Y,Y,N,Y,Y)
OUTPUT(U1 ,Y,Y,Y,N,Y,Y)
OUTPUT(V1 ,Y,Y,Y,N,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 13 ;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
NOWIPE = T
GVIEW(P,0.000000E+00,0.000000E+00,-1.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, 3.200000E-02, 3.200000E-02, 1.400000E+00
> DOM, MONIT, 1.550000E-02, 7.500000E-03, 7.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 3.000000E-02
> DOM, SNAPSIZE, 1.000000E-03
> DOM, RELAX, 5.000000E-01
> OBJ1, NAME, B2
> OBJ1, POSITION, 1.600000E-02, 1.600000E-02, 0.000000E+00
> OBJ1, SIZE, 1.600000E-02, 1.600000E-02, 1.400000E+00
> OBJ1, CLIPART, hh1
> OBJ1, ROTATION24, 1
> OBJ1, GRID, 2.000000E+00
> OBJ1, TYPE, BLOCKAGE
> OBJ1, MATERIAL, 198
> OBJ2, NAME, B3
> OBJ2, POSITION, 0.000000E+00, 1.600000E-02, 0.000000E+00
> OBJ2, SIZE, 1.600000E-02, 1.600000E-02, 1.400000E+00
> OBJ2, CLIPART, hh1
> OBJ2, ROTATION24, 2
> OBJ2, GRID, 2.000000E+00
> OBJ2, TYPE, BLOCKAGE
> OBJ2, MATERIAL, 198
> OBJ3, NAME, B6
> OBJ3, POSITION, 4.999997E-03, 7.300000E-03, 0.000000E+00
> OBJ3, SIZE, 8.000000E-03, 8.000000E-03, 1.400000E+00
> OBJ3, CLIPART, cylinder
> OBJ3, ROTATION24, 1
> OBJ3, GRID, 2.000000E+00
> OBJ3, TYPE, BLOCKAGE
> OBJ3, MATERIAL, 100
> OBJ3, FIXED_TMP, 0.000000E+00, 2.000000E+01
> OBJ4, NAME, B7
> OBJ4, POSITION, 3.200001E-03, 1.440001E-02, 0.000000E+00
> OBJ4, SIZE, 1.600000E-03, 1.600000E-03, 1.400000E+00
> OBJ4, CLIPART, cubet
> OBJ4, ROTATION24, 1
> OBJ4, VISIBLE, -1.000000E+00
> OBJ4, TYPE, PRESSURE_RELIEF
> OBJ4, PRES_RELIEF, 1.000000E+03, 0.000000E+00
> OBJ5, NAME, B8
> OBJ5, POSITION, 1.900001E-02, 7.300000E-03, 0.000000E+00
> OBJ5, SIZE, 8.000000E-03, 8.000000E-03, 1.400000E+00
> OBJ5, CLIPART, cylinder
> OBJ5, ROTATION24, 1
> OBJ5, GRID, 2.000000E+00
> OBJ5, TYPE, BLOCKAGE
> OBJ5, MATERIAL, 100
> OBJ5, FIXED_TMP, 0.000000E+00, 1.500000E+02
STOP
distil=t
EX(P1 )= 2.133E-02
EX(U1 )= 1.225E-02
EX(V1 )= 1.384E-02
EX(PRPS)= 2.152E+01
EX(TEM1)= 8.837E+01