TALK=T;RUN( 1, 1)
#cls
display
1-D PARSOL HEAT CONDUCTION WITH HEAT FLUX ON OBJECT
1000 [W/M3] BY VR
Ý--------------------Ý\\\\\\Ý---------------------------Ý
AIR AL AIR
x= 0. 0.45 0.65 1.1
BC.T=0.[k] T=0.[K]
Analytical Solution
-----------------------------------------------------------
ix = 1 2 3 4 5 6
T = 193.7 581. 968.4 1356. 1744 1744
-----------------------------------------------------------
ix = 7 8 9 10 11
T = 1744 1356. 968.4 581. 193.7
-----------------------------------------------------------
Note: This case was used to test PARSOL accuracy
TEM is analytical solution
ERR is TEM/TEM1
enddis
#pause
************************************************************
Q1 created by VDI menu, Version 3.4, Date 17/12/01
CPVNAM=VDI;SPPNAM=Core
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(PARSOL 1D HEAT TRANSFER)
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,11,1,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(148) =PRPS ; NAME(149) =TEM1
* Solved variables list
SOLVE(TEM1)
* Stored variables list
STORE(PRPS,TEM,ERR)
* Additional solver options
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
save7begin
(stored tem is 193.7 with if(ix.eq.1))
(stored tem is 581.0 with if(ix.eq.2))
(stored tem is 968.4 with if(ix.eq.3))
(stored tem is 1356. with if(ix.eq.4))
(stored tem is 1744. with if(ix.eq.5))
(stored tem is 1744. with if(ix.eq.6))
(stored tem is 1744. with if(ix.eq.7))
(stored tem is 1356. with if(ix.eq.8))
(stored tem is 968.4 with if(ix.eq.9))
(stored tem is 581.0 with if(ix.eq.10))
(stored tem is 193.7 with if(ix.eq.11))
(stored err is tem/tem1)
save7end
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
SETPRPS(1, 0)
ENUT = 0.000000E+00
DRH1DP = 1.000000E-05
DVO1DT = 3.410000E-03
PRNDTL(TEM1) = -2.580000E-02
EMISS = 0.000000E+00 ;SCATT = 0.000000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS) = -1.000000E+00 ;FIINIT(TEM1) = 0.000000E+00
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
No PATCHes used for this Group
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 350
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
LITER(TEM1)=1000
************************************************************
Group 17. Relaxation
RELAX(TEM1,FALSDT, 1.000000E+07)
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
PARSOL = T
************************************************************
Group 20. Preliminary Printout
distil=t; dsttol=1e-02
EX(PRPS)= 1.000E+01
EX(TEM1)= 1.039E+03
EX(ERR )= 1.000E+00
EX(TEM )= 1.039E+03
ECHO = T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 5 ;IYMON = 1 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 1
IXPRF = 1 ;IXPRL = 14
ISWPRF = 1 ;ISWPRL = 100000
IPROF = 3
PATCH (MYPLOT ,PROFIL,1,14,1,1,1,1,1,1)
PLOT(MYPLOT ,TEM1, 0.000000E+00, 0.000000E+00)
************************************************************
Group 24. Dumps For Restarts
NOWIPE = T
************************************************************
Group 25. DEBUG Group
DEBUG = T
DBEDGE = T
DBCOMP = T
DBCMPE = T
IREGDB = -1
DBSOL2 = T
DBGPHI(TEM1) = T
GVIEW(P,0.000000E+00,0.000000E+00,-1.000000E+00)
GVIEW(UP,1.000000E+00,0.000000E+00,0.000000E+00)
> DOM, SIZE, 1.100000E+00, 1.000000E+00, 1.000000E+00
> DOM, MONIT, 4.500000E-01, 5.000000E-01, 5.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> GRID, RSET_X_1, 11, 1.000000E+00
> GRID, RSET_Y_1, 1, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> OBJ, NAME, B3
> OBJ, POSITION, 0.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.100000E+00, 0.000000E+00, 1.000000E+00
> OBJ, CLIPART, cube11
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, PLATE
> OBJ, NAME, B5
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.100000E+00, 0.000000E+00, 1.000000E+00
> OBJ, CLIPART, cube11
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, PLATE
> OBJ, NAME, B7
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ, CLIPART, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, SURF_TEMP, 0.000000E+00, 0.000000E+00
> OBJ, NAME, B8
> OBJ, POSITION, 1.100000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ, CLIPART, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, SURF_TEMP, 0.000000E+00, 0.000000E+00
> OBJ, NAME, B6
> OBJ, POSITION, 4.500000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 2.000000E-01, 1.000000E+00, 1.000000E+00
> OBJ, CLIPART, cube4
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 100
> OBJ, HEAT_FLVX, 0.000000E+00, 1.000000E+03
> OBJ, INI_TEMP, -1.000000E+01
M