TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.6, Date 12/09/05
CPVNAM=VDI;SPPNAM=Core
Echo DISPLAY / USE settings
DISPLAY
Two-dimensional steady-state heat conduction problem
tests symmetrical solution.
|-----------------|
|\\\\\\\\\\\\\ T=1|
| \\\\\\ \\\\\\\ |
| \\\\ k2 \\\\\|
| \\\\ \\\\\\|
| \\\\\\\\\|
| \\\\\\\|
| \\\\\|
| k1 \\\|
Y |T=0 \|
| |-----------------|
+-- X
The solution should be lineary for k1=k2.
The 161 number property is the solid with incompressible
air properties.
ENDDIS
PHOTON USE
p;;;;
vi z;up y
msg Temperature contours
con TEM1 z 1 fil;.001
msg Press e to END
use patgeo
ENDUSE
************************************************************
IRUNN = 1 ;LIBREF = 709
************************************************************
REAL(XLEN,K1,K2)
INTEGER(NCEL)
NCEL=3
mesg(NX number of cells is :ncel: OK? If not, insert uneven desired number,
readvdu(ncel,int,:ncel:)
ncel
XLEN=3.0
K1=1.0;K2=K1 ! heat conductivity
Group 1. Run Title
TEXT(2D PARSOL heat-conduction, xy plane
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,NCEL,NCEL,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(148) =TEM1 ; NAME(150) =PRPS
* Solved variables list
SOLVE(TEM1)
* Stored variables list
STORE(PRPS,KOND,DVAR,DRES,DCOR)
STORE(S_EX,D_EX) ! exact solution at sunny and dark
STORE(SE/S,DE/D) ! exact/newparsol_solution for sunny and dark
* Additional solver options
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
* Domain material is: 0 Air at 20 deg C, 1 atm, treated as i
SETPRPS(1, 0)
ENUT = 0.000000E+00
DVO1DT = 3.410000E-03
PRNDTL(TEM1) = -K1
CSG10='Q1' ! signal use of the following properties line
! which correspond to steel
MATFLG=T;NMAT=1
<161> artificial solid material having same properties as material 0
161 1.189 1.544E-05 1005.0 1. 3.41E-3 0.0
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS) = -1.000000E+00
FIINIT(DVAR) = -1.234E-11
FIINIT(DRES) = -1.234E-11
FIINIT(DCOR) = -1.234E-11
FIINIT(S_EX) = 0.; FIINIT(D_EX) = 0.
FIINIT(SE/S) = 0.; FIINIT(DE/D) = 0.
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
PATCH(HOT,VOLUME,NX,NX,NY,NY,NZ,NZ,1,1)
COVAL(HOT,TEM1,1.E6,1.0)
PATCH(COLD,VOLUME,1,1,1,1,1,1,1,1)
COVAL(COLD,TEM1,1.E6,0.0)
inform7begin
(stored var rst1 is resi(tem1))
(stored var crt1 is corr(tem1))
(stored s_ex at hot is 1.0)
patch(struct1,cell,2,2,1,1,1,nz,1,1)
(stored s_ex at struct1 is .325)
patch(struct2,cell,1,1,2,2,1,nz,1,1)
(stored s_ex at struct2 is .325)
patch(struct3,cell,nx-1,nx-1,ny,ny,1,nz,1,1)
(stored s_ex at struct3 is .675)
patch(struct4,cell,nx-1,nx-1,ny,ny,1,nz,1,1)
(stored s_ex at struct4 is .675)
patch(cut1,cell,nx,nx,1,1,1,nz,1,1)
(stored s_ex at cut1 is .4708)
(stored d_ex at cut1 is .5292)
patch(cut2,cell,1,1,ny,ny,1,nz,1,1)
(stored s_ex at cut2 is .4708)
(stored d_ex at cut2 is .5292)
patch(cut3,cell,nx/2+1,nx/2+1,ny/2+1,ny/2+1,1,nz,1,1)
(stored s_ex at cut3 is .45)
(stored d_ex at cut3 is .55)
(stored var SE/S is s_ex/tem1 with if(tem1!=0.))
(stored var DE/D is d_ex/dvar with if(dvar!=0.))
inform7end
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 100
************************************************************
Group 16. Terminate Iterations
LITER (TEM1) = 100
RESREF(TEM1)=0.
ENDIT(TEM1)=0.
************************************************************
Group 17. Relaxation
relax(tem1,linrlx,1.)
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
PARSOL = T
************************************************************
Group 20. Preliminary Printout
ECHO = T
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL = 1.000000E-02
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 1 ;IZMON = -3
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 1 ; NYPRIN = 1
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
patch(plot,profil,1,nx,1,1,1,1,1,1)
coval(plot,tem1,0,0)
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, XLEN, XLEN, 1.000000E+00
> DOM, MONIT, 5.000000E-01, 5.000000E-01, 5.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> GRID, RSET_X_1, NCEL, 1.000000E+00
> GRID, RSET_Y_1, NCEL, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> OBJ, NAME, BLOCK
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, XLEN, XLEN, 1.000000E+00
> OBJ, GEOMETRY, wedge1
> OBJ, ROTATION24, 8
> OBJ, GRID, NO
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 161, air
> OBJ, INI_TEMP, 0.000000E+00
STOP
isg60=1 ! set active NewParsol