#cls display 1-D PARSOL HEAT CONDUCTION WITH LINEAR HEAT SOURCE ON OBJECT: SU = 0.5*(1000-PHI) BY INFORM Ý--------------------Ý\\\\\\Ý---------------------------Ý AIR AL AIR x= 0. 0.45 0.65 1.1 BC.T=0.[k] T=0.[K] Correct Solution Heat Flux (out) = 53.4 [W] ----------------------------------------------------------- ix = 1 2 3 4 5 6 T = 51.8 155.3 258.9 362.5 466. 466. ----------------------------------------------------------- ix = 7 8 9 10 11 T = 466. 362.5 258.9 155.3 51.8 ----------------------------------------------------------- Note: This case was used to test PARSOL accuracy enddis #pause TALK=T;RUN( 1, 1) ************************************************************ 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-X 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) * Additional solver options SOLUTN(TEM1,Y,Y,Y,N,N,Y) ************************************************************ 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 = 50 RESFAC = 1.000000E-03 ************************************************************ Group 16. Terminate Iterations ************************************************************ 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)= 2.777E+02 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) INFORM13BEGIN (SOURCE OF TEM1 AT b6 IS 0.5*(1000-tem1) with line!volu) INFORM13END > 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 M