photon use
p
up y
msg the grid. press return for temperature contours
gr z 1;pause; gr off;red;gr ou z 1
msg temperature contours. press return for velocity vectors
con tem1 z 1 fi;0.0002;pause; con off;red;se re ve 2;vec z 1 sh
msg velocity vectors. press return for reduced-pressure contours
pause;vec off;red;con p1 z 1 fi;0.0002
enduse
************************************************************
Group 1. Run Title
TEXT(NON-ISOTHERMAL 2D LAMINAR MIXING LAYER
************************************************************
Group 2. Transience
STEADY = T
DISPLAY
This case is a illustratiom of the Colocated-Covariant
Method (CCM). The geometry is with two inlets at differents
temperatures. Both sides walls are symetry planes. The
velocities difference and the temperatures entail a
recirculation zone behind the solid edge.
The solid conductivity is temperature dependent:
lenda=14.17+6.25E-3 T
(T in K)
The kinematic viscosity of air follows the Sutherland law.
ENDDIS
L(PAUSE
DISPLAY
Geometry: outlet
--------
P1=P2=1.013E5 Pa
| | ^ H2=0.1 m
| | |
| | |
| | |
| | |
| | |
| |\ | - H1=0.1 m
| |//\ | |
Symetry | |///| | Symetry |
plane | |///| | plane |
|^^^^|///| | |
| V1 |///|^^^^^^| |
| |///| V2 | |
air air
T1=600K T2=300K
V1=2m/s V2=0.5m/s
ENDDIS
L(PAUSE
**********************|**************************************
Groups 3, 4, 5 Grid Information
INTEGER(NX1,NX2,NX3,NXTOT)
NX1=20;NX2=10;NX3=25;NXTOT=NX1+NX2+NX3
INTEGER(NY1,NY2,NY3,NYTOT)
NY1=25;NY2=20;NY3=30;NYTOT=NY1+NY2+NY3
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,NXTOT,NYTOT,1)
* Set overall domain extent:
* xulast yvlast zwlast name
XSI= 8.000E-03;YSI= 2.000E-01;ZSI= 1.000E-03;RSET(D,CHAM )
************************************************************
Group 6. Body-Fitted coordinates
BFC=T
* Set points
XPO= 0.0000E+00;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P1 )
XPO= 0.0000E+00;YPO= 9.2000E-02;ZPO= 0.0000E+00;GSET(P,P2 )
XPO= 0.0000E+00;YPO= 1.0000E-01;ZPO= 0.0000E+00;GSET(P,P3 )
XPO= 0.0000E+00;YPO= 2.0000E-01;ZPO= 0.0000E+00;GSET(P,P4 )
XPO= 1.5000E-03;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P5 )
XPO= 1.5000E-03;YPO= 9.2000E-02;ZPO= 0.0000E+00;GSET(P,P6 )
XPO= 1.5000E-03;YPO= 1.0000E-01;ZPO= 0.0000E+00;GSET(P,P7 )
XPO= 3.0000E-03;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P9 )
XPO= 3.0000E-03;YPO= 9.2000E-02;ZPO= 0.0000E+00;GSET(P,P10 )
XPO= 2.0000E-03;YPO= 1.0000E-01;ZPO= 0.0000E+00;GSET(P,P11 )
XPO= 8.0000E-03;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P13 )
XPO= 8.0000E-03;YPO= 9.2000E-02;ZPO= 0.0000E+00;GSET(P,P14 )
XPO= 8.0000E-03;YPO= 1.0000E-01;ZPO= 0.0000E+00;GSET(P,P15 )
XPO= 8.0000E-03;YPO= 2.0000E-01;ZPO= 0.0000E+00;GSET(P,P16 )
* Set lines/arcs
GSET(L,L1,P6,P7,NY2,1.0)
GSET(L,L2,P10,P11,NY2,1.0)
GSET(L,L3,P14,P15,NY2,1.0)
GSET(L,L4,P2,P3,NY2,1.0)
GSET(L,L5,P7,P11,NX2,1.0)
GSET(L,L6,P1,P5,NX1,-1.3)
GSET(L,L7,P5,P9,NX2,1.0)
GSET(L,L8,P9,P13,NX3,1.3)
GSET(L,L9,P13,P14,NY1,S1.5)
GSET(L,L10,P9,P10,NY1,S1.5)
GSET(L,L11,P5,P6,NY1,S1.5)
GSET(L,L12,P1,P2,NY1,S1.5)
GSET(L,L13,P15,P16,NY3,1.8)
GSET(L,L14,P3,P4,NY3,1.8)
GSET(L,L18,P3,P7,NX1,-1.3)
GSET(L,L19,P11,P15,NX3,1.3)
GSET(L,L20,P10,P14,NX3,1.3)
GSET(L,L21,P4,P16,NX1+NX2+NX3,1.0)
* Set frames
GSET(F,F1,P1,-,P5,P6,P7,-,P3,P2)
GSET(F,F2,P5,-,P9,P10,P11,-,P7,P6)
GSET(F,F3,P9,-,P13,-,P14,-,P10,-)
GSET(F,F4,P10,-,P14,-,P15,-,P11,-)
GSET(F,F5,P3,P7.P11,P15,-,P16,-,P4,-)
* Match a grid mesh
GSET(M,F1,+I+J,1,1,1,TRANS)
GSET(M,F2,+I+J,NX1+1,1,1,TRANS)
GSET(M,F3,+I+J,NX1+NX2+1,1,1,TRANS)
GSET(M,F4,+I+J,NX1+NX2+1,NY1+1,1,TRANS)
GSET(M,F5,+I+J,1,NY1+NY2+1,1,TRANS)
* Copy/Transfer/Block grid planes
GSET(C,K2,F,K1,1,NXTOT,1,NYTOT,+,0,0,1.0000E-03,INC,1)
**********
NONORT = T
* X-cyclic boundaries switched
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME( 42) =ENUL ; NAME( 43) =KOND
NAME( 44) =BLOK ; NAME( 45) =TEM1
NAME( 46) =WCRT ; NAME( 47) =VCRT
NAME( 48) =UCRT ; NAME( 49) =DEN1
NAME( 50) =PRPS
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,TEM1)
* Stored variables list
STORE(PRPS,DEN1,UCRT,VCRT,WCRT,BLOK,KOND,ENUL)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
CCM=T
IVARBK = -1 ;ISOLBK = 1
************************************************************
Group 8. Terms & Devices
NEWRH1 = T
NEWENL = T
ISOLX = 0 ;ISOLY = 0 ;ISOLZ = 0
************************************************************
Group 9. Properties
PRESS0 = 1.013E+05
ENUT = 0.000E+00
* List of user-defined materials to be read by EARTH
MATFLG=T;IMAT=2
* Name
*Ind. Dens. Viscos. Spec.heat Conduct. Expans. Compr.
*
152 8000 0 500 GRND4 0 0
* constants for GRND option no 1
1.4170E+01 6.2500E-03
*
31 GRND5 GRND6 1004.0 0.0263 3.33E-3 GRND5
* constants for GRND option no 1 (Ideal gas law)
0 3.4843E-03
* constants for GRND option no 2 (Sutherland law)
1.4600E-06 110
* constants for GRND option no 3
0 0 7.1429E-01
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(BLOK) = 1.000E+00 ;FIINIT(TEM1) = 4.500E+02
FIINIT(PRPS) = 3.100E+01
CONPOR(STEEL , -1.00,CELL ,-#2,-#2,-#1,-#2,-#1,-#1)
INIT(STEEL ,BLOK, 0.000E+00, 2.000E+00)
INIT(STEEL ,PRPS, 0.000E+00, 1.520E+02)
INIADD = F
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
INLET (BFCIN1 ,SOUTH ,#1,#1,#1,#1,#1,#1,#1,#1)
VALUE (BFCIN1 ,P1 , GRND1 )
VALUE (BFCIN1 ,U1 , GRND1 )
VALUE (BFCIN1 ,V1 , GRND1 )
VALUE (BFCIN1 ,TEM1, 6.000E+02)
VALUE (BFCIN1 ,VCRT, 2.000E+00)
VALUE (BFCIN1 ,DEN1, 5.883E-01)
INLET (BFCIN2 ,SOUTH ,#3,#3,#1,#1,#1,#1,#1,#1)
VALUE (BFCIN2 ,P1 , GRND1 )
VALUE (BFCIN2 ,U1 , GRND1 )
VALUE (BFCIN2 ,V1 , GRND1 )
VALUE (BFCIN2 ,TEM1, 3.000E+02)
VALUE (BFCIN2 ,VCRT, 5.000E-01)
VALUE (BFCIN2 ,DEN1, 1.176E+00)
OUTLET(OUT ,NORTH ,#1,#3,#3,#3,#1,#1,#1,#1)
VALUE (OUT ,P1 , 0.000E+00)
VALUE (OUT ,U1 , 0.000E+00)
VALUE (OUT ,V1 , 0.000E+00)
BFCA = 1.189E+00
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 600
SELREF = T
RESFAC = 1.000E-02
************************************************************
Group 16. Terminate Iterations
LITER (U1 ) = 2 ;LITER (V1 ) = 2
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 9.000E-01)
RELAX(U1 ,FALSDT, 5.062E-02)
RELAX(V1 ,FALSDT, 5.062E-02)
RELAX(TEM1,FALSDT, 1.062E+03)
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
CSG10 =Q1
************************************************************
Group 22. Monitor Print-Out
IXMON = NX1+NX2-1 ;IYMON = NY1+NY2+2 ;IZMON = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NXPRIN = 1; NYPRIN = 1
IPROF = 3; ITABL = 3