TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2009, Date 14/05/09 CPVNAM=VDI;SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY This case modeling a heat exchange by heat conductivity about a filament. The filament is simulated by the cylinder shape with properties of steel. It is surrounded by air. The air temperature has constant significance 18 degree at all domain boundaries. The filament is heated up by a electrical current which occurs in the conductor at presence the filament votage. The emitted heat is calculated as the filament votage divided on thermal resistance which is function of the conductor temperature. Source = Voltage^2 / Resistance where Resistance = 0.1 * Temperature The heating source for the HEATING object is set by In-Form as per unit volume (SOURCE of TEM1 at HEATING is 220^2/(0.1*TEM1) with VOLU) The Q1 contains PHOTON USE commands ENDDIS PHOTON USE p;;; con tem1 z 1 fi;0.0001 enduse ************************************************************ IRUNN = 1 ;LIBREF = 148 ************************************************************ Group 1. Run Title TEXT(Heat exchange about a steel conductor ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,55,55,1) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(149)=PRPS ;NAME(150)=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 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 PRNDTL(TEM1)=-0.0258 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(PRPS)=-1. 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 EGWF = T ************************************************************ Echo save-block settings for Group 13 SAVE13BEGIN (SOURCE of TEM1 at HEATING is 220^2/(0.1*TEM1) with VOLU) SAVE13END ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 500 RESFAC =1.0E-03 ************************************************************ Group 16. Terminate Iterations LITER(TEM1)=500 ENDIT(TEM1)=1.0E-20 ************************************************************ Group 17. Relaxation ************************************************************ Group 18. Limits VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-204.75 ************************************************************ Group 19. EARTH Calls To GROUND Station ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(PRPS)=1.582; EX(TEM1)=20.93 ************************************************************ Group 21. Print-out of Variables ************************************************************ Group 22. Monitor Print-Out IXMON = 41 ;IYMON = 30 ;IZMON = 1 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 NXPRIN = 1 IXPRF = 25 ;IXPRL = 32 NYPRIN = 1 IYPRF = 25 ;IYPRL = 32 ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts GVIEW(P,0.,0.,-1.) GVIEW(UP,0.,1.,0.) > DOM, SIZE, 1.000000E-01, 1.000000E-01, 1.000000E-01 > DOM, MONIT, 7.281251E-02, 5.250000E-02, 5.000000E-02 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, RSET_X_1, 25, 1.000000E+00 > GRID, RSET_X_2, 6, 1.000000E+00 > GRID, RSET_X_3, 24, 1.000000E+00 > GRID, RSET_Y_1, 25, 1.000000E+00 > GRID, RSET_Y_2, 6, 1.000000E+00 > GRID, RSET_Y_3, 24, 1.000000E+00 > GRID, RSET_Z_1, 1, 1.000000E+00 > OBJ, NAME, HEATING > OBJ, POSITION, 4.500000E-02, 4.500000E-02, 0.000000E+00 > OBJ, SIZE, 1.000000E-02, 1.000000E-02, 1.000000E-01 > OBJ, GEOMETRY, cylinder > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 111,STEEL at 27 deg c (C = 1%) > OBJ, INI_TEMP, 18. > OBJ, NAME, W_WALL > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 1.000000E-01, 1.000000E-01 > OBJ, GEOMETRY, cube13 > OBJ, TYPE, PLATE > OBJ, SURF_TEMP, 0. ,18. > OBJ, NAME, E_WALL > OBJ, POSITION, 1.000000E-01, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 1.000000E-01, 1.000000E-01 > OBJ, GEOMETRY, cube13 > OBJ, TYPE, PLATE > OBJ, SURF_TEMP, 0. ,18. > OBJ, NAME, S_WALL > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E-01, 0.000000E+00, 1.000000E-01 > OBJ, GEOMETRY, cube13 > OBJ, TYPE, PLATE > OBJ, SURF_TEMP, 0. ,18. > OBJ, NAME, N_WALL > OBJ, POSITION, 0.000000E+00, 1.000000E-01, 0.000000E+00 > OBJ, SIZE, 1.000000E-01, 0.000000E+00, 1.000000E-01 > OBJ, GEOMETRY, cube13 > OBJ, TYPE, PLATE > OBJ, SURF_TEMP, 0. ,18. STOP