PHOTON USE
  p
 
 
 
  msg Temperature contours
  con h1 z m fi;0.001
  con h1 y m fi;0.001
  con h1 x m fi;0.001
  msg Press e to END
  ENDUSE
 
    GROUP 1. Run title and other preliminaries
TEXT(3D STEADY HEAT CONDUCTION IN A CUBE:108
  RADIATION AND CONVECTION FROM A HOT BLOCK

  DISPLAY 
  This standard Library Q1 file (case 100) has been modified
  to exemplify the use of PLANT for external heat-transfer laws.
 
  This feature allows inclusion of external heat loss by way of
  radiation (a (Text**4-Ts**4)) and forced convection a(Text-Ts).
  ENDDIS
 
REAL(XLENGTH,YLENGTH,ZLENGTH)
XLENGTH=1.0;YLENGTH=1.0;ZLENGTH=1.0
NX=5;NY=5;NZ=5
    GROUP 3. X-direction grid specification
GRDPWR(X,NX,XLENGTH,1.0)
    GROUP 4. Y-direction grid specification
GRDPWR(Y,NY,YLENGTH,1.0)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,NZ,ZLENGTH,1.0)
    GROUP 7. Variables stored, solved & named
solve(h1)
    GROUP 8. Terms (in differential equations) & devices
   **For pure conduction, cut out built-in source and convection
     terms
TERMS(h1,N,N,Y,N,Y,Y)
    GROUP 9. Properties of the medium (or media)
RHO1=1.0;PRNDTL(h1)=1
    GROUP 12. Local adjustments of convection and diffusion
    GROUP 13. Boundary conditions and special sources
   **Corner at IX=IY=IZ=1
PATCH(COLD,CELL,1,1,1,1,1,1,1,1)
   **Fix temperature to zero
COVAL(COLD,h1,FIXVAL,0.0)
   **Corner at IX=NX, IY=NY, IZ=NZ
PATCH(HOT,CELL,NX,NX,NY,NY,NZ,NZ,1,1)
   **Fix temperature to 1.0
COVAL(HOT,h1,FIXVAL,1.0)
   external radiation patch
real(radco,tmpx)
radco=1.e-6; tmpx=100.
RG(1)=radco;rg(2)=tmpx

PATCH(SHOWRAD,high,1,NX,1,NY,NZ,NZ,1,1)
  PLANTBEGIN
PATCH(RAD,high,1,NX,1,NY,NZ,NZ,1,1)
    CO=RG(1)*(RG(2)**2+H1**2)*(RG(2)+H1)
    VAL=RG(2)
COVAL(RAD,h1,GRND,GRND)
  PLANTEND

   external  convection patch
real(convco)
convco=10.
PATCH(CONV,north,1,NX,ny,NY,1,NZ,1,1)
COVAL(CONV,h1,convco,tmpx/2.)
    GROUP 15. Termination of sweeps
LSWEEP=100
    GROUP 17. Under-relaxation devices
RELAX(h1,LINRLX,0.5)
    GROUP 19.
NAMSAT=MOSG
    GROUP 20. Preliminary print-out
    GROUP 21. Print-out of variables
OUTPUT(h1,Y,Y,Y,Y,Y,Y)
    GROUP 22. Spot-value print-out
IXMON=NX/2+1;IYMON=NY/2+1;IZMON=NZ/2+1
    GROUP 23. Field print-out and plot control
NXPRIN=NX/5;NYPRIN=NY/5;NZPRIN=NZ/5;nplt=1
tstswp=-1
dmpstk=t
DISTIL=T
EX(H1)=6.820E+01
 LIBREF=108
STOP