************************************************************ Q1 created by VDI menu, Version 3.5, Date 15/03/02 CPVNAM=VDI;SPPNAM=Core ************************************************************ DISPLAY This In-Form case provides a VR object for use in the setting of source into one half of T-junction. This example simulates a flow in a T-Junction of which contains resistance in one part. In-Form sets momentum source within the limits of VR object with HALFCYL3 name. The Q1 contains PHOTON USE commands ENDDIS Echo DISPLAY / USE settings PHOTON USE p;;; rot y ang 180 gr ou z 1; gr ou x 1; gr ou y m vec x 1 sh; vec z 1 sh; red msg Velocity vectors msg Press return for p1 pause; vec off;red msg Pressure contours con p1 x 1 fi;0.001; con p1 z 1 fi;0.001 enduse ************************************************************ IRUNN = 1 ; LIBREF = 774 ************************************************************ Group 1. Run Title TEXT(T-Junction: In-Form Source. ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,12,20,30) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(150) =PRPS * Solved variables list SOLVE(P1 ,U1 ,V1 ,W1 ) * Stored variables list STORE(PRPS,VABS) * Additional solver options SOLUTN(P1 ,Y,Y,Y,N,N,Y) SOLUTN(U1 ,Y,Y,Y,N,N,Y) SOLUTN(V1 ,Y,Y,Y,N,N,Y) SOLUTN(W1 ,Y,Y,Y,N,N,Y) ************************************************************ Group 8. Terms & Devices TERMS (P1 ,Y,Y,Y,Y,Y,Y) ************************************************************ Group 9. Properties PRESS0 = 1.000000E+05 SETPRPS(1,111) ENUT = 0.000000E+00 DRH1DP = 2.000000E+11 DVO1DT = 3.700000E-06 EMISS = 0.000000E+00 ;SCATT = 0.000000E+00 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(PRPS) = 1.110000E+02 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 No PATCHes used for this Group ************************************************************ Echo InForm settings for Group 13 INFORM13BEGIN (SOURCE of U1 at HALFCYL3 is coval(fixflu,-10.*VABS*U1)) (SOURCE of V1 at HALFCYL3 is coval(fixflu,-10.*VABS*V1)) (SOURCE of W1 at HALFCYL3 is coval(fixflu,-10.*VABS*W1)) INFORM13END ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 200 RESFAC = 1.000000E-03 ************************************************************ Group 16. Terminate Iterations LITER (P1 ) = 50 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 7.000000E-01) ************************************************************ Group 18. Limits VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06 VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06 VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06 ************************************************************ Group 19. EARTH Calls To GROUND Station USEGRD = T ;USEGRX = T ASAP = T CONWIZ = T ************************************************************ Group 20. Preliminary Printout ECHO = T DISTIL = T ;NULLPR = F NDST = 0 DSTTOL = 1.000000E-02 EX (P1 ) = 2.350000E-01 ;EX (U1 ) = 2.025000E-02 EX (V1 ) = 1.611000E-01 ;EX (W1 ) = 1.616000E-01 EX (VABS) = 2.849000E-01 ;EX (PRPS) = 5.233000E+01 ************************************************************ Group 21. Print-out of Variables OUTPUT(P1 ,Y,Y,Y,Y,Y,Y) OUTPUT(U1 ,Y,Y,Y,Y,Y,Y) OUTPUT(V1 ,Y,Y,Y,Y,Y,Y) OUTPUT(W1 ,Y,Y,Y,Y,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 4 ;IYMON = 14 ;IZMON = 15 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts NOWIPE = T GVIEW(P,-8.604726E-01,0.000000E+00,-5.094969E-01) GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00) > DOM, SIZE, 6.000000E-01, 1.500000E+00, 3.000000E+00 > DOM, MONIT, 1.750000E-01, 9.642860E-01, 1.450000E+00 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, SNAPSIZE, 1.000000E-02 > OBJ, NAME, HALFCYL1 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.100000E+00 > OBJ, SIZE, 4.000000E-01, 1.000000E+00, 8.000000E-01 > OBJ, CLIPART, hfcyl > OBJ, ROTATION24, 6 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 1 > OBJ, SCAL_FIXF, 1.000000E+00 > OBJ, NAME, INLET > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.100000E+00 > OBJ, SIZE, 4.000000E-01, 0.000000E+00, 8.000000E-01 > OBJ, CLIPART, hfcyl > OBJ, ROTATION24, 6 > OBJ, TYPE, INLET > OBJ, DENSITY, 1.189000E+00 > OBJ, VELOCITY, 0.000000E+00, 1.000000E+00, 0.000000E+00 > OBJ, TEMPERATURE, 0.000000E+00 > OBJ, NAME, HALFCYL2 > OBJ, POSITION, 0.000000E+00, 5.000000E-01, 0.000000E+00 > OBJ, SIZE, 5.000000E-01, 1.000000E+00, 2.000000E+00 > OBJ, CLIPART, hfcyl > OBJ, ROTATION24, 4 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 1 > OBJ, SCAL_FIXF, 1.000000E+00 > OBJ, NAME, OUTLET1 > OBJ, POSITION, 0.000000E+00, 5.000000E-01, 0.000000E+00 > OBJ, SIZE, 5.000000E-01, 1.000000E+00, 0.000000E+00 > OBJ, CLIPART, hfcyl > OBJ, ROTATION24, 4 > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0.000000E+00 > OBJ, TEMPERATURE, 0.000000E+00 > OBJ, COEFFICIENT, 1.000000E+00 > OBJ, NAME, OUTLET2 > OBJ, POSITION, 0.000000E+00, 5.000000E-01, 3.000000E+00 > OBJ, SIZE, 5.000000E-01, 1.000000E+00, 0.000000E+00 > OBJ, CLIPART, hfcyl > OBJ, ROTATION24, 4 > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0.000000E+00 > OBJ, TEMPERATURE, 0.000000E+00 > OBJ, COEFFICIENT, 1.000000E+00 > OBJ, NAME, HALFCYL3 > OBJ, POSITION, 0.000000E+00, 5.000000E-01, 2.000000E+00 > OBJ, SIZE, 5.000000E-01, 1.000000E+00, 1.000000E+00 > OBJ, CLIPART, hfcyl > OBJ, ROTATION24, 4 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 1 > OBJ, SCAL_FIXF, 1.000000E+00 STOP