TALK=F;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.5, Date 15/05/03
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
TEXT(Temperature-controlled fan
TITLE
DISPLAY
This case simulates a FAN-like device located inside to the domain
area and turned revolved under some angle concerning the horizon.
The fan is described by a box objects.
The fan is switched on automatically when air temperature value
in some cell will be more installed limits. Otherwise it is
switched off.
The In-Form is used for the setting of a momentum source for
whole object as
(SOURCE of U1 at FAN is :USORC: with WHOLOB!IF(SWTH.GT.0))
(SOURCE of V1 at FAN is :VSORC: with WHOLOB!IF(SWTH.GT.0))
The fan will work when the SWTH single real variable will be
more zero. Otherwise the fan will be switched off.
The significance of SWTH single variable is defined by
temperature value in cell with indexes IX=IXP and IY=IYP
and TLIM limiting temperature
(STORE1 SWTH is TEM1[:IXP:,:IYP:]-:TLIM: with TSTFIN)
The alternate way is the direct indication of x and y coordinates
{0.05,0.05} of sensor position inside domain
(STORE1 SWTH is TEM1{0.05,0.05}-:TLIM: with TSTFIN)
The Q1 contains PHOTON USE commands
ENDDIS
PHOTON USE
p
p1;;;
con tem1 z 1 fi;.0001
ve z 1
use patgeo
upause 2
p
p2;;;
con tem1 z 1 fi;.0001
ve z 1
use patgeo
upause 2
p
p3;;;
con tem1 z 1 fi;.0001
ve z 1
use patgeo
upause 2
p
p4;;;
con tem1 z 1 fi;.0001
ve z 1
use patgeo
upause 2
p
p5;;;
con tem1 z 1 fi;.0001
ve z 1
use patgeo
upause 2
p
p6;;;
con tem1 z 1 fi;.0001
ve z 1
use patgeo
enduse
************************************************************
IRUNN = 1 ;LIBREF = 152
************************************************************
Group 1. Run Title
TEXT(Temperature-controlled fan
Group 2. Transience
STEADY=F
* Set overall time and no. of steps
RSET(U,0.000000E+00,3.600000E+03,6)
Groups 3, 4, 5 Grid Information
RSET(M,40,40,1)
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(148) =SPH1 ; NAME(149) =TEM1
NAME(150) =PRPS
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,TEM1)
* Stored variables list
STORE(PRPS,SPH1)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
Group 9. Properties
PRESS0 = 1.000000E+05 ;TEMP0 = 2.730000E+02
SETPRPS(1, 0)
PRNDTL(TEM1) = 1.000000E+00
Group 11.Initialise Var/Porosity Fields
FIINIT(TEM1) = 1.800000E+01 ;FIINIT(PRPS) = -1.000000E+00
Group 13. Boundary & Special Sources
Echo InForm settings for Group 13
INFORM13BEGIN
INTEGER(IXP,IYP); IXP=NX/2; IYP=NY/2
REAL(ANGL); ANGL=3.14159/4.
REAL(INVEL); INVEL=0.4
REAL(AREA); AREA=0.05*0.1
REAL(PSORC); PSORC=RHO1*INVEL*AREA
REAL(USORC); USORC=INVEL*COS(ANGL)*PSORC
REAL(VSORC); VSORC=INVEL*SIN(ANGL)*PSORC
(SOURCE of U1 at FAN is :USORC: with WHOLOB!IF(SWTH.GT.0))
(SOURCE of V1 at FAN is :VSORC: with WHOLOB!IF(SWTH.GT.0))
REAL(TLIM); TLIM=40.
(MAKE of SWTH is 0)
InForm settings the indexes of sensor
(STORE1 SWTH is TEM1[:IXP:,:IYP:]-:TLIM: with TSTFIN)
InForm settings the x and y coordinates of sensor
(STORE1 SWTH is TEM1{0.05,0.05}-:TLIM: with TSTFIN)
INFORM13END
Group 15. Terminate Sweeps
LSWEEP = 500
Group 16. Terminate Iterations
LITER (TEM1) = 100
ENDIT (TEM1) = 1.000000E-06
Group 17. Relaxation
RELAX(P1 ,LINRLX, 8.000000E-01)
RELAX(U1 ,FALSDT, 1.000000E+00)
RELAX(V1 ,FALSDT, 1.000000E+00)
RELAX(TEM1,FALSDT, 1.000000E+01)
Group 18. Limits
VARMIN(TEM1) = 1.800000E+01
Group 20. Preliminary Printout
DISTIL = T
EX (P1 ) = 5.564000E-08 ;EX (U1 ) = 1.288000E-05
EX (V1 ) = 1.242000E-05 ;EX (SPH1) = 9.997000E+02
EX (TEM1) = 4.519000E+01 ;EX (PRPS) = 3.722000E+00
Group 22. Monitor Print-Out
IXMON = 20 ;IYMON = 20 ;IZMON = 1
TSTSWP = -1
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
Group 24. Dumps For Restarts
IDISPA = 1
CSG1 ='P'
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-01, 1.000000E-01, 1.000000E-01
> DOM, MONIT, 5.000000E-02, 5.000000E-02, 5.000000E-02
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> OBJ, NAME, FAN
> OBJ, POSITION, 5.000000E-02, 1.500000E-02, 0.000000E+00
> OBJ, SIZE, 1.000000E-02, 5.000000E-02, 1.000000E-01
> OBJ, ROT-ANGLE, 4.500000E+01, 0.000000E+00, 0.000000E+00
> OBJ, CLIPART, cubet1
> OBJ, ROTATION24, 1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, -1
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, HEAT_FLUX, 0.000000E+00, 1.000000E-10
> OBJ, INI_PRESS, 0.000000E+00
> OBJ, SCAL_FIXF, 0.000000E+00
> OBJ, NAME, E_OUTLET
> OBJ, POSITION, 1.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 1.000000E-01, 1.000000E-01
> OBJ, CLIPART, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, SURF_TEMP, 0.000000E+00, 1.800000E+01
> OBJ, NAME, N_OUTPUT
> OBJ, POSITION, 0.000000E+00, 1.000000E-01, 0.000000E+00
> OBJ, SIZE, 1.000000E-01, 0.000000E+00, 1.000000E-01
> OBJ, CLIPART, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, SURF_TEMP, 0.000000E+00, 1.800000E+01
> OBJ, NAME, W_OUTLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 1.000000E-01, 1.000000E-01
> OBJ, CLIPART, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, SURF_TEMP, 0.000000E+00, 1.800000E+01
> OBJ, NAME, S_OUTLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E-01, 0.000000E+00, 1.000000E-01
> OBJ, CLIPART, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, SURF_TEMP, 0.000000E+00, 1.800000E+01
> OBJ, NAME, HEATBLOC
> OBJ, POSITION, 8.400000E-02, 4.900000E-02, 0.000000E+00
> OBJ, SIZE, 1.000000E-02, 5.000000E-02, 1.000000E-01
> OBJ, ROT-ANGLE, 4.500000E+01, 0.000000E+00, 0.000000E+00
> OBJ, CLIPART, cube4
> OBJ, ROTATION24, 1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 100
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, HEAT_FLUX, 0.000000E+00, 5.000000E+00
> OBJ, NAME, PR_REL
> OBJ, POSITION, 2.440777E-03, 2.500000E-03, 0.000000E+00
> OBJ, SIZE, 2.440777E-03, 2.500000E-03, 1.000000E-01
> OBJ, CLIPART, cubet
> OBJ, ROTATION24, 1
> OBJ, TYPE, PRESSURE_RELIEF
> OBJ, PRES_RELIEF, 1.000000E-01, 0.000000E+00
STOP