TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 3.6, Date 08/09/04 CPVNAM=VDI;SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY This In-Form case shows how heat and momentum sources can be specified for VR objects, in this case the HEATSINK object placed in the upper part of the domain. Hot fluid enters through INFLOW and cooled fluid exits through OUTFLOW. The non-dimensional temperature is represented by enthalpy, H1, which equals 1.0 for the INFLOW fluid, and is 0.0 for the HEATSINK object. The heat-transfer coefficient is proportional to the relative velocity squared(calculated without allowance for staggering) to the power 0.4 (i.e. velocity to the power 0.8); and the friction coefficient is proportion to the relative velocity squared to the power 0.9 (i.e. velocity to the power 1.8). Multiplying coefficients are arbitrary. The Q1 contains PHOTON USE commands. ENDDIS PHOTON USE p msg flow in a duct containing a heat sink vec z 1 sh msg VELOCITY VECTORS msg Pressto continue pause vec off;red msg PRESSURE CONTOURS con p1 z 1 fil;0.01 msg Press to continue pause con off;red msg TEMPERATURE CONTOURS con h1 z 1 fil;0.01 pause con off;red msg HEAT-FLUX CONTOURS con hflx z 1 x 1,m y 15,m fil;0.01 msg Press e to END ENDUSE ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Flow In A Model Furnace: In-Form Source.) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,20,21,1) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(150) =PRPS * Solved variables list SOLVE(P1 ,U1 ,V1 ,H1 ) * Stored variables list STORE(PRPS) * Additional solver options SOLUTN(P1 ,Y,Y,Y,N,N,Y) SOLUTN(H1 ,Y,Y,Y,N,N,Y) ************************************************************ Echo InForm settings for Group 7 inform7begin REAL(HEATCO) HEATCO=5.0 (STORED VAR HFLX IS -:HEATCO:*(U1^2+V1^2)^0.4*H1*VOL) inform7end ************************************************************ Group 8. Terms & Devices TERMS (P1 ,Y,Y,Y,Y,Y,Y) TERMS (H1 ,N,Y,Y,Y,Y,Y) ************************************************************ Group 9. Properties RHO1 = 1.000000E+00 ENUL = 1.000000E-02 CP1 = 1.005000E+03 ENUT = 0.000000E+00 DVO1DT = 3.410000E-03 PRNDTL(H1 ) = 7.000000E-01 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(H1 ) = 0.000000E+00 ;FIINIT(PRPS) = -1.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 EGWF = T ************************************************************ Echo InForm settings for Group 13 INFORM13BEGIN REAL(FRICCO) FRICCO=10.0 (SOURCE of H1 at HEATSINK is -:HEATCO:*(U1^2+V1^2)^0.4*H1*VOL with $ LINE) (SOURCE of U1 at HEATSINK is -:FRICCO:*(U1^2+V1^2)^0.9*U1*VOL with $ LINE) (SOURCE of V1 at HEATSINK is -:FRICCO:*(U1^2+V1^2)^0.9*V1*VOL with $ LINE) INFORM13END ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 400 RESFAC = 1.000000E-04 ************************************************************ Group 16. Terminate Iterations ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 1.000000E+00) ************************************************************ Group 18. Limits VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06 VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06 ************************************************************ Group 19. EARTH Calls To GROUND Station USEGRD = T ;USEGRX = T ASAP = T PARSOL = T LSG57 = T ISG50 = 1 ISG52 = 2 SPEDAT(SET,FACETDAT,NUMOBJ,I,5) ************************************************************ Group 20. Preliminary Printout ECHO = T DISTIL = T ;NULLPR = F NDST = 0 DSTTOL = 1.000000E-02 EX (P1 ) = 2.504000E+02 ;EX (U1 ) = 2.890000E+00 EX (V1 ) = 6.489000E+00 ;EX (H1 ) = 5.909000E-01 EX (PRPS) = 7.905000E-01 ************************************************************ 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 = 18 ;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 LIBREF=775 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, 2.000000E+00, 1.000000E+00 > DOM, MONIT, 6.250000E-01, 1.666500E+00, 5.000000E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, SNAPSIZE, 1.000000E-02 > OBJ, NAME, BLOCK3 > OBJ, POSITION, 5.000000E-01, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 5.000000E-02, 1.333000E+00, 1.000000E+00 > OBJ, GEOMETRY, cube14 > OBJ, ROTATION24, 1 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 102,BRICK at 20 deg c > OBJ, NAME, INFLOW > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 5.714000E-01, 1.000000E+00 > OBJ, GEOMETRY, cube3t > OBJ, ROTATION24, 1 > OBJ, TYPE, INLET > OBJ, PRESSURE, 0.000000E+00 > OBJ, VELOCITY, 8.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, TEMPERATURE, 1.000000E+00 > OBJ, NAME, OUTPUT > OBJ, POSITION, 1.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 5.714000E-01, 1.000000E+00 > OBJ, GEOMETRY, cube12t > OBJ, ROTATION24, 1 > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0.000000E+00 > OBJ, TEMPERATURE, SAME > OBJ, COEFFICIENT, 1.000000E+00 > OBJ, NAME, HEATSINK > OBJ, POSITION, 0.000000E+00, 1.333300E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E+00, 6.667000E-01, 1.000000E+00 > OBJ, GEOMETRY, cubet1 > OBJ, ROTATION24, 1 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, DOMAIN > OBJ, USER_HEAT, 0.000000E+00, 0.000000E+00 > OBJ, SCAL_FIXF, 0.000000E+00 STOP