GROUP 1. Run title and other preliminaries
TEXT(In-Form Represents Hot Moving Spheres
TITLE
libref=765
DISPLAY
The two spheres move from left to right within the upper and
lower halves of the domain, specified by PATCHes FIRST and
SECOND.
Their motion is specified by INFOB commands, wherein the
formula takes the form:
SPHERE(x-centre, y-centre, z-centre, radius) ,
each of the arguments is itself a formula.
SOURCE commands specify how much heat the spheres emit,
also by way of formulae.
Several formulae are provided from which the user can choose.
WARNING:
If the case is to be run by way of the VR-Editor, the patches
should be extended to cover the whole domain; for otherwise
the Editor will convert them to objects with different names.
Then a single patch will suffice, e.g.
PATCH(WHOLEDOM,1,NX,1,NY,1,LSTEP)
to which the INFOB statements should refer.
The Q1 contains PHOTON USE commands
ENDDIS
GROUP 2. Transience; time-step specification
STEADY=F;GRDPWR(T,70,700.,1.0)
GROUP 3. X-direction grid specification
GRDPWR(X,80,4.,1.0)
GROUP 4. Y-direction grid specification
GRDPWR(Y,40,2.,1.0)
GROUP 5. Z-direction grid specification
GRDPWR(Z,1,0.1,1.0)
GROUP 7. Variables stored, solved & named
NAME(150)=TEM1
SOLVE(TEM1)
PRNDTL(TEM1)=0.1
GROUP 9.Properties used
CP1=1000.
GROUP 11. Initialization of variable or porosity fields
FIINIT(TEM1)=20.
GROUP 13. Boundary conditions and special sources
PATCH(WALLE,EWALL,NX,NX,1,NY,1,NZ,1,LSTEP)
COVAL(WALLE,TEM1,1.0,20.0)
PATCH(WALLW,WWALL,1,1,1,NY,1,NZ,1,LSTEP)
COVAL(WALLW,TEM1,1.0,20.0)
PATCH(WALLN,NWALL,1,NX,NY,NY,1,NZ,1,LSTEP)
COVAL(WALLN,TEM1,1.0,20.0)
PATCH(WALLS,SWALL,1,NX,1,1,1,NZ,1,LSTEP)
COVAL(WALLS,TEM1,1.0,20.0)
INFORM13BEGIN
*** Hot sphere - 1
PATCH(FIRST,VOLUME,1,NX,1,NY/2,1,NZ,1,LSTEP)
*** Hot sphere - 2
PATCH(SECOND,VOLUME,1,NX,NY/2+1,NY,1,NZ,1,LSTEP)
Note: The FIRST and SECOND patches cover the parts of the domain
in which the spheres move, namely the upper and the lower
halves.
Where the spheres actually are at any time is conveyed by
the INFOB statements, which also assign the (obligatory)
names INFOB_1 and INFOB_2.
Note: xce,yce and zce are the x,z and z coordinates of the sphere
centre, and radius has its obvious meaning.
They are character variables which are evaluated by
Satellite because they are enclosed within colons.
char(xce,yce,zce,radius,sour)
xce=.5+.5*(tim/100.-1); yce=.5; zce=.05; radius= .25
The above setting of yce causes the first sphere to follow a
horizontal path.
Move the following statement one ot two spaces to the left if the
path is to have a gradient of 0.1/0.5; and also increase the
PATCH(FIRST statement to NY, so as to ensure that the patch is
large enough.
yce=0.1 + 01*(tim/100.-1)
(INFOB at FIRST is SPHERE(:xce:,:yce:,:zce:,:radius:$
) with INFOB_1)
xce=.5+.25*(tim/100.-1); yce=1.5; zce=.05; radius= .25
(INFOB at SECOND is SPHERE(:xce:,:yce:,:zce:,:radius:$
) with INFOB_2)
*** Setting of TEM1 sources into the two spheres
sour=100.
(SOURCE of TEM1 at FIRST is :sour: with INFOB_1)
The following statement, moved to the left so as to become active,
would cause the magnitude of the soulce to increase with
temperature, a potentially explosive situation.
sour=10.0 + 2.0*tem1
Otherwise, the sour=100. setting will be used
(SOURCE of TEM1 at SECOND is :sour: with INFOB_2)
INFORM13END
GROUP 15. Termination of sweeps
LSWEEP=30
GROUP 22. Spot-value print-out
NPRMON=LSWEEP;IXMON=NX/4+1;IYMON=NY/4+1;IZMON=NZ/2+1
GROUP 23. Print-out & plot control
TSTSWP=-1;IDISPA=10; CSG1=P
SELREF=T; RESFAC=1.E-2
NXPRIN=1;NYPRIN=1
DISTIL=T;EX(TEM1)=2.295E+01 ;EX(SPH1)= 1.000E+03
PHOTON USE
p
p10;
gr ou z 1
con tem1 z 1 fil;.00001
msg at 100 sec.
pause
p
p20;
gr ou z 1
con tem1 z 1 fil;.00001
msg at 200 sec.
pause
p
p30;
gr ou z 1
con tem1 z 1 fil;.00001
msg at 300 sec.
pause
p
p40;
gr ou z 1
con tem1 z 1 fil;.00001
msg at 400 sec.
pause
p
p50;
gr ou z 1
con tem1 z 1 fil;.00001
msg at 500 sec.
pause
p
p60;
gr ou z 1
con tem1 z 1 fil;.00001
msg at 600 sec.
pause
p
p70;
gr ou z 1
con tem1 z 1 fil;.00001
msg at 700 sec.
ENDUSE