TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2007, Date 22/06/07 CPVNAM=VDI;SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY This two-dimensional heat conduction problem tests the correctness of NEWPARSOL solver which gives correct result at one sweep. The domain end object materials are air. The object does not fit the grid. Each cut-cell has neighbouring cut-cells. Y|-------!-------!-------|-------| | ! ! | | | 1.5 ! ! | 2 | | ! ! | | |-------!-------!-------|-------| | ! /*!*\ | | | ! /***!***\ | | | ! /*****!*****\ | | |-------!-------!-------|-------| | ! \*****!*****/ | | | ! \***!***/ | | | ! \*!*/ | | |-------!-------!-------|-------| | ! ! | | | 1 ! ! | 1.5 | | ! ! | | |-------!-------!-------|-------|---X ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 710 ************************************************************ Group 1. Run Title TEXT(Ellipsoid simulation ) ************************************************************ Group 2. Transience STEADY=t * Set overall time and no. of steps RSET(U,0.000000E+00,9.000000E+03,75) * Modify regions RSET(T,1,450,1.000000E+00) ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,4,4,1) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T SOLVE(TEM1) * Stored variables list STORE(DVAR,DRES,DCOR,cutc,volu) ************************************************************ Group 8. Terms & Devices ************************************************************ Group 9. Properties PRESS0 = 1.000000E+05 ;TEMP0 = 2.730000E+02 * Domain material index is 0 signifying: * Air at 20 deg C, 1 atm, treated as incompressible SETPRPS(1, 0) DVO1DT = 3.410000E-03 PRNDTL(TEM1) = -2.580000E-02 MATFLG=T; NMAT=1 152 1.189 1.544E-05 1005.0 0.0258 3.41E-3 0.0 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields 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(COLD,CELL,1,1,1,1,1,NZ,1,LSTEP) COVAL(COLD,TEM1,FIXVAL,1.000000E+00) PATCH(HOT,CELL,NX,NX,NY,NY,1,NZ,1,LSTEP) COVAL(HOT,TEM1,FIXVAL, 2.000000E+00) EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 1 RESFAC = 1.000000E-03 ************************************************************ Group 16. Terminate Iterations LITER (TEM1) = 200 ************************************************************ Group 17. Relaxation RELAX(TEM1,linrlx, 1.0) ************************************************************ Group 19. EARTH Calls To GROUND Station USEGRD = T ;USEGRX = T ASAP = T PARSOL = T CSG10 ='Q1' ************************************************************ Group 20. Preliminary Printout ECHO = T ************************************************************ Group 21. Print-out of Variables ************************************************************ Group 22. Monitor Print-Out IXMON = 20 ;IYMON = 1 ;IZMON = 17 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 10000 ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts IDISPA = 1 ;IDISPB = 0 ;IDISPC = 450 CSG1 ='T' GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00) GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00) > DOM, SIZE, 5.000000E-01, 5.000000E-01, 5.000000E-01 > DOM, MONIT, 3.383813E-01, 2.584409E-01, 2.540000E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, SNAPSIZE, 1.000000E-02 > GRID, RSET_Y_2, 2, 1.000000E+00 > OBJ, NAME, Block > OBJ, POSITION, 2.490000E-01, 1.090000E-01, 0.000000E+00 > OBJ, SIZE, 2.000000E-01, 2.000000E-01, 5.000000E-01 > OBJ, ROT-ANGLE, 4.500000E+01, 0.000000E+00, 0.000000E+00 > OBJ, GEOMETRY, cube > OBJ, ROTATION24, 14 > OBJ, TYPE, BLOCKAGE > OBJ, COLOR-MODE, USER > OBJ, COLOR-VAL, 192 > OBJ, OPAQUE, 50 > OBJ, MATERIAL, 152,Egg > OBJ, TIME_LIMITS, ALWAYS_ACTIVE STOP ISG60=1