dbs 03.06.01 The link featuren in this does not work. I am
therefore scrapping the whole feature, because In-Form
can do the job much better
871
PHOTON USE
p;;;
msg Temperature contours
con tem1 z m fi;0.001; con tem1 y m fi;0.001
con tem1 x m fi;0.001
msg
msg Press e to END
ENDUSE
GROUP 1. Run title and other preliminaries
TEXT(3D Steady Heat Conduction In A Cube
TITLE
DISPLAY
This modification of Library case 100 exemplifies the use of
the Linear Temperature Link.
Linear Temperature link allows 'pipe-like' heat transfer by
providing sub-grid-scale links within the computational domain. If
there are thermal leads, thermal bridges or thin metal plates
embedded in a poorly-conducting medium that would prevent use of
coarse grid, this feature can directly link originating and
destination cells.
enddis
REAL(XLENGTH,YLENGTH,ZLENGTH)
XLENGTH=1.0;YLENGTH=1.0;ZLENGTH=1.0
NX=5;NY=5;NZ=5
domain size and grid
NX=10;NY=10;NZ=10;xulast=xlength;yvlast=ylength;zwlast=zlength
#unigrid
GROUP 7. Variables stored, solved & named
solve(tem1);SOLUTN(TEM1,Y,Y,Y,N,N,Y)
store(prps)
GROUP 8. Terms (in differential equations) & devices
**For pure conduction, cut out built-in source and convection
terms
TERMS(TEM1,N,N,Y,N,Y,Y)
GROUP 9. Properties of the medium (or media)
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,TEM1,1.E2,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,TEM1,1.E2,1.0)
L($021)
mesg(insert linear-link patches? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y) then
linear link patches
patch(-a,west,1,1,2,2,2,4,1,1)
coval(-a,tem1,1.0e2,0.0)
patch(-b,east,nx,nx,2,2,3,3,1,1)
coval(-b,tem1,1.0e2,0.0)
patch(-c,south,1,2,ny,ny,1,2,1,1)
coval(-c,tem1,1.0,0.0)
patch(-d,north,4,5,ny,ny,4,5,1,1)
coval(-d,tem1,1.0,0.0)
patch(-e,cell,3,4,3,4,nz,nz,1,1)
coval(-e,tem1,1.e5,0.0)
patch(-f,cell,3,4,ny,ny,3,4,1,1)
coval(-f,tem1,1.0e5,0.0)
patch(-g,cell,1,1,ny,ny,nz,nz,1,1)
coval(-g,tem1,1.0e2,0.0)
patch(-h,cell,1,1,1,1,nz,nz,1,1)
coval(-h,tem1,1.e2,0.0)
patch(heater,north,3,4,ny,ny,3,4,1,1)
coval(heater,tem1,fixval,10.0)
mesg(linear-link patches are active
endif
L($021)
GROUP 15. Termination of sweeps
OUTPUT(PRPS,N,N,N,N,N,N);LSWEEP=100
GROUP 16. Termination of iterations
LITER(TEM1)=100; OVRRLX= 1.7
ISOLX=1;ISOLY=1;ISOLZ=1
GROUP 17. Under-relaxation devices
RELAX(TEM1,LINRLX,0.5)
GROUP 21. Print-out of variables
OUTPUT(TEM1,Y,Y,Y,Y,Y,Y); OUTPUT(PRPS,N,N,N,N,N,N)
GROUP 22. Spot-value print-out
IXMON=NX/2+1;IYMON=NY/2+1;IZMON=NZ/2+1;UWATCH=T
GROUP 23. Field print-out and plot control
NXPRIN=NX/5;NYPRIN=NY/5;NZPRIN=NZ/5;nplt=1;uwatch=t
GROUP 24. Dumps for restarts