**********************************************
TITLE T.001: TESTING SOURCES - Source option
**********************************************
GRID NODEs 22 by 22
DEBUG GEOMETRY OFF
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
/
WALL
OUTL X+
/
LOCA (2, 2) TO (3,22)
BLOC SELECTED
LOCA (4, 2) TO (5, 4)
BLOC SELECTED
LOCA (4,19) TO (5,21)
BLOC SELECTED
LOCA (4,10) TO (5,13)
BLOC SELECTED
/
LOCATE ID=IN1 FROM (4, 5) TO (4, 9)
LOCATE ID=IN2 FROM (4,14) TO (4,18)
/
SOUR FLOW AREA X- ID=IN1 200 U=100 T = 1. K = 30 L = 0.1
SOUR FLOW AREA X- ID=IN2 200 U=100 T = 1. K = 50 L = 0.2
/
DENSity 2.0
VISC 1.E-5
/
DIAGnostic output for U V K E at (11,11) every 20 steps
FLUX DEFAult output OFF
SELECt (1,1) TO (999,999) with frequency of (4,2)
OUTPut U V P K L E T NARROW in SELEcted region
/
CONV GLOB 1.E-10
SOLVE STEADY 500
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.002: TESTING SOURCES - Source option
************************************************************************
GRID NODEs 22 by 22
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMETRY OFF
/
WALL
OUTL X+
/
LOCA (2, 2) TO (3,22)
BLOC SELECTED
LOCA (4, 2) TO (5, 4)
BLOC SELECTED
LOCA (4,19) TO (5,21)
BLOC SELECTED
LOCA (4,10) TO (5,13)
BLOC SELECTED
/
LOCATE ID=IN1 FROM (4, 5) TO (4, 9)
LOCATE ID=IN2 FROM (4,14) TO (4,18)
/
FLOW AREA X- ID=IN1 200 U=100 T = 1. K = 30 L = 0.1
FLOW AREA X- ID=IN2 200 U=100 T = 1. K = 50 L = 0.2
/
DENSity 2.0
VISC 1.E-5
/
DIAGnostic output for U V K E at (11,11) every 20 steps
FLUX DEFAult output OFF
SELECT (1,1) TO (99,99) with frequency of (4,2)
OUTPut U V P K L E T NARROW in SELEcted region
/
CONV GLOB 1.E-10
SOLVE STEADY 500
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.003: TESTING SOURCES - INLEt option
************************************************************************
GRID NODEs 22 by 22
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMETRY OFF
/
WALL
OUTL X+
/
LOCA (2, 2) TO (3,22)
BLOC SELECTED
LOCA (4, 2) TO (5, 4)
BLOC SELECTED
LOCA (4,19) TO (5,21)
BLOC SELECTED
LOCA (4,10) TO (5,13)
BLOC SELECTED
/
LOCATE ID=IN1 FROM (4, 5) TO (4, 9)
LOCATE ID=IN2 FROM (4,14) TO (4,18)
/
INLET X- ID=IN1
INLET X- ID=IN2
BOUN X- ID=IN1 U=100
BOUN X- ID=IN2 U=100
BOUN X- ID=IN1 T=1
BOUN X- ID=IN2 T=1
BOUN X- ID=IN1 K=30
BOUN X- ID=IN2 K=50
BOUN X- ID=IN1 L=0.1
BOUN X- ID=IN2 L=0.2
/
DENSity 2.0
VISC 1.E-5
/
DIAGnostic output for U V K E at (11,11) every 20 steps
FLUX DEFAult output OFF
SELEct (1,1) TO (99,99) with frequency of (4,2)
OUTPut U V P K L E T NARROW in SELEcted region
/
CONV GLOB 1.E-10
SOLVE STEADY 500
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.004: TESTING SOURCES - Source option
************************************************************************
GRID NODEs 22 by 22
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMETRY OFF
/
WALL
OUTL X+

LOCA (2, 2) TO (3,22)
BLOC SELECTED
LOCA (4, 2) TO (5, 4)
BLOC SELECTED
LOCA (4,19) TO (5,21)
BLOC SELECTED
LOCA (4,10) TO (5,13)
BLOC SELECTED
/
LOCATE ID=IN1 FROM (4, 5) TO (4, 9)
LOCATE ID=IN2 FROM (4,14) TO (4,18)
/
SOUR FLOW AREA X- ID=IN1 200 U=100 T = 1. K = 30 L = 0.1
SOUR FLOW AREA X- ID=IN2 200 U=100 T = 1. K = 50 L = 0.2
/
DENSity 2.0
VISC 1.E-5
/
DIAGnostic output for U V K E at (11,11) every 20 steps
FLUX DEFAult output OFF
SELEct (1,1) TO (99,99) with frequency of (4,2)
OUTPut U V P K L E T NARROW in SELEcted region
/
CONV GLOB 1.E-10
SOLVE STEADY 500
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.005: TESTING SOURCES - Source option
************************************************************************
GRID NODEs 22 by 22
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMETRY OFF
/
WALL
OUTL X+
/
LOCA (2, 2) TO (3,22)
BLOC SELECTED
LOCA (4, 2) TO (5, 4)
BLOC SELECTED
LOCA (4,19) TO (5,21)
BLOC SELECTED
LOCA (4,10) TO (5,13)
BLOC SELECTED
/
LOCATE ID=IN1 FROM (4, 5) TO (4, 9)
LOCATE ID=IN2 FROM (4,14) TO (4,18)
/
FLOW AREA X- ID=IN1 50 U=100 T=1. K=30 L=0.1 DENSITY = 4
FLOW AREA X- ID=IN2 100 U=100 T=1. K=50 L=0.2 NORMAL 0.5 0 DENSity 4.0
/
DENSity 2.0
VISC 1.E-5
/
DIAGnostic output for U V K E at (11,11) every 20 steps
FLUX DEFAult output OFF
SELEct (1,1) TO (99,99) with frequency of (4,2)
OUTPut U V P K L E T NARROW in SELEcted region
/
CONV GLOB 1.E-10
SOLVE STEADY 500
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.006: TESTING SOURCES
************************************************************************
GRID NODEs 22 by 22
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMETRY OFF
/
WALL
OUTL X+

LOCA (1, 1) TO (1,9) ID=INLET
/
LOCA (4,19) TO (5,21)
BLOC SELECTED
LOCA (4,10) TO (7,12)
BLOC SELECTED
/
LOCATE ID=SLOTB1 FROM (4,9) TO (7, 9)
LOCATE ID=SLOTB2 FROM (4,13) TO (7,13)
LOCATE ID=SLOTB3 FROM (8,10) TO (8,12)
/
SOUR FLOW AREA Y+ ID=SLOTB1 20 V=-10 T = 1. K = 30 L = 0.1
FLOW AREA Y- ID=SLOTB2 20 V= 10 T = 1. K = 50 L = 0.1
FLOW AREA X- ID=SLOTB3 20 U= 20 T = 1. K = 50 L = 0.1
/
/
INLET X- ID=INLET
FLOW ONLY AREA X- ID=INLET 200
/
BOUN X- ID=INLET U=100
BOUN X- ID=INLET K=30
BOUN X- ID=INLET L=0.1
/
DENSity 2.0
VISC 1.E-5
/
DIAGnostic output for U V K E at (11,11) every 40 steps
FLUX DEFAult output OFF
SELEct (1,1) TO (99,99) with frequency of (4,2)
OUTPut U V P K L E T NARROW in SELEcted region
/
CONV GLOB 1.E-10
SOLVE STEADY 500
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.007: TESTING SOURCES
************************************************************************
GRID NODEs 22 by 22
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMETRY OFF
WALL
OUTL X+
LOCA (2, 2) TO (3,22)
BLOC SELECTED
LOCA (4, 2) TO (5, 4)
BLOC SELECTED
LOCA (4,10) TO (5,12)
BLOC SELECTED
LOCA (4,20) TO (5,21)
BLOC SELECTED
/
LOCATE ID=IN1 FROM (4, 5) TO (4, 6)
LOCATE ID=IN2 FROM (4, 7) TO (4, 8)
LOCATE ID=IN3 FROM (4, 9) TO (4, 9)
LOCATE ID=IN4 FROM (4,13) TO (4,14)
LOCATE ID=IN5 FROM (4,15) TO (4,16)
LOCATE ID=IN6 FROM (4,17) TO (4,18)
LOCATE ID=IN7 FROM (4,19) TO (4,19)
/
SOUR FLOW AREA X- ID=IN1 200 U=100 T=1. K=10 L=0.04
SOUR FLOW AREA X- ID=IN2 200 T=2. K=20 L=0.08 MOMENTUM
SOUR FLOW AREA X- ID=IN3 50 U=100 T=3. K=30 L=0.12 NORM 4. 0.
SOUR FLOW AREA X- ID=IN4 20 U=100 T=4. K=40 L=0.16 TOTAL ! Area=0.10
SOUR FLOW AREA X- ID=IN5 20 T=5. K=50 L=0.20 TOTAL MOME
SOUR FLOW AREA X- ID=IN6 200 T=6. K=60 L=0.22 NORM 1. 0. MOME
SOUR FLOW AREA X- ID=IN7 20 U=100 T=7. K=70 L=0.24 NORM 2. 0. DENS=5
/
DENSity 2.0
VISC 1.E-5
/
DIAGnostic output for U V K E P at (11,11) every 20 steps
FLUX DEFAult output OFF
SELEct (1,1) TO (99,99) with frequency of (4,2)
OUTPut U V P K L E T NARROW in SELEcted region
/
CONV GLOB 1.E-8
SOLVE STEADY 500
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.008: TESTING SOURCES
************************************************************************
GRID NODEs 22 by 22
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMETRY OFF
WALL
OUTL X+
LOCA (2, 2) TO (3,22)
BLOC SELECTED
LOCA (4, 2) TO (5, 4)
BLOC SELECTED
LOCA (4,10) TO (5,12)
BLOC SELECTED
LOCA (4,20) TO (5,21)
BLOC SELECTED
/
LOCATE ID=IN1 FROM (4, 5) TO (4, 6)
LOCATE ID=IN2 FROM (4, 7) TO (4, 8)
LOCATE ID=IN3 FROM (4, 9) TO (4, 9)
LOCATE ID=IN4 FROM (4,13) TO (4,14)
LOCATE ID=IN5 FROM (4,15) TO (4,16)
LOCATE ID=IN6 FROM (4,17) TO (4,18)
LOCATE ID=IN7 FROM (4,19) TO (4,19)
/
FLOW AREA X- ID=IN1 200 U=100 T=1. K=10 L=0.04
FLOW AREA X- ID=IN2 200 T=2. K=20 L=0.08 MOMENTUM
FLOW AREA X- ID=IN3 50 U=100 T=3. K=30 L=0.12 NORM 4. 0.
FLOW AREA X- ID=IN4 20 U=100 T=4. K=40 L=0.16 TOTAL ! Area=0.10
FLOW AREA X- ID=IN5 20 T=5. K=50 L=0.20 TOTAL MOME
FLOW AREA X- ID=IN6 200 T=6. K=60 L=0.22 NORM 1. 0. MOME
FLOW AREA X- ID=IN7 20 U=100 T=7. K=70 L=0.24 NORM 2. 0. DENS=5
/
DENSity 2.0
VISC 1.E-5
/
DIAGnostic output for U V K E P at (11,11) every 20 steps
FLUX DEFAult output OFF
SELEct (1,1) TO (99,99) with frequency of (4,2)
OUTPut U V P K L E T NARROW in SELEcted region
/
CONV GLOB 1.E-8
SOLVE STEADY 500
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.009: Testing gradients
************************************************************************
Grid 3 by 12
COORdinate X Range 0.5
COORdinate Y Range 1.0
DEBUG GEOMERTY OFF
/
SET U Polynomial of Y 1. 0 1. 0. 0.
ALLOcate DUDY
ALLOcate D2UDY
SET DUDY = Gradient of U with Y ! LINEAR
SET D2UDY = Gradient of DUDY with Y
/
SOLVE OFF
/
ALLOCATE ERR1
SET ERR1 LINEAR 0. +2 * Y
SET ERR1 LINEAR 0. -1. * DUDY ADD
/
ALLOCATE ERR2
SET ERR2 LINEAR 2. -1. * D2UDY
/
OUTPUT U ERR1 ERR2 DUDY D2UDY NARROW
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.010: Test Writing Diagnostic File & WRTe file
************************************************************************
GRID NODEs 22 by 12
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMERTY OFF
/
WALL at undefined outer boundaries
LOCA ID=INLT (1,1) (1,21)
INLET ID=INLT X-
OUTL X+
SYMM Y-
/
SET U at X- boundary of ID=INLT 100.0
SET T at X- boundary of ID=INLT 1.0
/
DENSity 1.0
VISC 1.0
SPEC 1.0
PRAN EFFE 0.7
/ SOLUTION OPTIONS
RELA U=0.8, V=0.8, P=0.5, T=1.0
/
/----------------------------------------------------------------------/
DIAGnostic output U V P T RP at (6,11) every 20 steps 'diag_eqn.tmp'
/
/FLUX DEFAult output OFF
SELEct (1,1) TO (99,99) with frequency of (4,2)
OUTPut U T NARROW in SELEcted region
SAVE U V P T AF AFX AFY AFZ on file 'T10.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-8
SOLVE STEADY 21
/
DIAGnostic NOW
OUTPUT NOW in SELEcted region
SOLVE STEADY 10
WRITE AF FC 'AREAS.TMP'
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.011: Test Variable Diffusion coefficient
************************************************************************
GRID NODEs 22 by 12
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMERTY OFF
/
WALL at undefined outer boundaries
WALLS FIXED FOR C
LOCA ID=INLT (1,1) (1,21)
INLET ID=INLT X-
OUTL X+
SYMM Y-
/
SET U at X- boundary of ID=INLT 100.0
SET C at X- boundary of ID=INLT 1.0
BOUN C at Y+ 0
/
DENSity 1.0
VISC 1.0
SCHMIDT EFFE 1.0
SCHMIDT 1.0
DIFFUSION FOR C XX LINEAR 1. +10. * X
DIFFUSION FOR C YY LINEAR 1. +1 * Y
/ SOLUTION OPTIONS
RELA C=0.9
/
/----------------------------------------------------------------------/
DIAGnostic output U V C at (6,11) every 20 steps
/
SELEct (1,1) TO (99,99) with frequency of (4,2)
OUTPut U C in NARRow mode for SELEcted region
SAVE OFF
/
CONV GLOB 1.E-8
SOLVE STEADY 200
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.012: Test Internal Wall and Mulitple Area Based Outlets
************************************************************************
GRID NODEs 12 by 12
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMERTY OFF
/
WALL at undefined outer boundaries
LOCA ID=WALLIN (1,5) (99,5)
WALL Y+ ID=WALLIN
/
LOCATE ID=INLET1 from ( 1, 1) to ( 1, 5)
LOCATE ID=INLET2 from ( 1, 6) to ( 1,99)
INLET X- ID=INLET1
INLET X- ID=INLET2
/
LOCATE ID=OUTLET1 from (99, 1) to (99, 5)
LOCATE ID=OUTLET2 from (99, 6) to (99,99)
OUTLET X+ ID=OUTLET1 AREA
OUTLET X+ ID=OUTLET2 AREA
/
SYMM Y-
/
SET U at X- boundary of ID=INLET1 1.0
SET U at X- boundary of ID=INLET2 1.0
SET T at X- boundary of ID=INLET1 1.0
SET T at X- boundary of ID=INLET2 1.0
/
DENSity 1.0
VISC 1.0
SPEC 1.0
PRAN EFFE 0.7
/ SOLUTION OPTIONS
RELA U=0.2, T=0.5
/
/----------------------------------------------------------------------/
DIAGnostic output U V P T RP at (6,5) every 10 steps
/
/FLUX DEFAult output OFF
SELEct (1,1) TO (99,99) with frequency of (2,2)
OUTPut U V T P FC in NARROW mode for SELEcted region
/
CONV GLOB 1.E-8
SOLVE STEADY 200
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.013: 2-in-1 split to illustrate OUTPut and FLOW commands
************************************************************************
GRID NODEs 22 by 22
DEBUG GEOMERTY OFF
/
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
WALL at undefined outer boundaries
LOCA ID=INLT (1,1) (1,21)
INLET ID=INLT X-
LOCA (999,1) (999,11)
OUTL X+ SELECTED split in ratio of AREAs
LOCA (999,12) (999,22)
OUTL X+ SELECTED split in ratio of AREAs
/
SET U at X- boundary of ID=INLT 100.0
SET T at X- boundary of ID=INLT 1.0
/
SELECT (1,11) (999,11)
FLOW ONLY INTERNAL 0. AT Y+ OF SELECTED
/
DENSity 1.0
VISC 1.0
SPEC 1.0
PRAN EFFE 0.7
/
RELA T=1.0
/ SOLUTION OPTIONS
/----------------------------------------------------------------------/
DIAGnostic output U V P T RP at (6,11) every 20 steps
/
FLUX DEFAult output OFF
FLUX AVERAGE FOR U to standard output device
FLUX AVERAGE FOR P to 'T13_FLUX_AVG.TMP'
FLUX AVERAGE FOR P to standard output device
SELEct (1,1) TO (99,99) with frequency of (4,2)
OUTPut U V T in NARROW mode for SELEcted region
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-10
/
SOLVE STEADY 1000
/
PRINT FLOW
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.014: Testing of TOTAL option with flow
************************************************************************
GRID NODEs 22 by 22
DEBUG GEOMETRY OFF
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
/
WALL
OUTL X+
/
LOCA (2, 2) TO (3,22)
BLOC SELECTED
LOCA (4, 2) TO (5, 4)
BLOC SELECTED
LOCA (4,19) TO (5,21)
BLOC SELECTED
LOCA (4,10) TO (5,13)
BLOC SELECTED
/
LOCATE ID=IN1 FROM (4, 5) TO (4, 9)
LOCATE ID=IN2 FROM (4,14) TO (4,18)
/
FLOW Y- Q=0.02 U =10, V = 50 K=30 L=.001 T=500 at ID=IN1 TOTAL
FLOW Y+ Q=0.02 U =10, V =-50 K=30 L=.001 T=500 at ID=iN2 TOTAL
/
DENSity 2.0
VISC 1.E-5
/
DIAGnostic output for U V K E at (11,11) every 25 steps
FLUX DEFAult output OFF
SELEct (1,1) TO (99,99) with frequency of (4,2)
OUTPut U V P K L E T in NARROW mode for SELEcted region
/
CONV GLOB 1.E-10
SOLVE STEADY 500
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.015: Testing of TRANSFER option with flow
************************************************************************
GRID NODEs 22 by 22 BY 22
DEBUG GEOMETRY OFF
COOR X RANGE 10.0
COOR Y RANGE 1.0
COOR Z RANGE 1.0
/ Initial and Boundary Conditions
INLET X-
OUTLet X+
SYMMETRY Y-
SYMMETRY Y+
SYMMETRY Z-
SYMMETRY Z+
/
SET U 10.0
SET U at X- 10.0
/
DENSity 1.0
VISC 1.0
/
/OBJECT SWIRLER vane angle 45. 'swirler.loc' ID=SWIRL
/ rin 0.2 rout 0.2 hub_rin 0.01 hub_rout 0.01
/ (1., 0.5, 0.5) (2., 0.5, 0.5)
/
LOCAte ID=BODY (9,10,10) TO (12,13,13)
BLOCK ID=BODY
/
LOCAte MATCH X- of ID=BODY and x+ of ID=BODY ID=SWRL 'BODY_LOC.TMP'
/
TRANSFER FLUX ID=SWRL
SWIRL ID=SWRL VANE ANGLE = 45. CENTER = (5.0 0.5 0.5)
/
RELA U=0.8, p=0.5 T=1.0
/
DIAGnostic output U V W P RP BP at (13,11,11) every 1 steps
/
CONV GLOB 1.E-8
/
LOCAte (12,9,9) TO (12,14,14)
OUTPUT U V W P NARROW YZ PLANE SELECT
/
SOLVE STEADY 10
/
SAVE OFF U V W P T on file 'T15_SAV.TMP'
/
END
/

-Back To Top

************************************************************************
TITLE T.016: Testing of TRANSFER FLUX option with flow
************************************************************************
GRID NODEs 22 by 22 BY 22
DEBUG GEOMETRY OFF
COOR X RANGE 10.0
COOR Y RANGE 1.0
COOR Z RANGE 1.0
/ Initial and Boundary Conditions
INLET X-
OUTLet X+
SYMMETRY Y-
SYMMETRY Y+
SYMMETRY Z-
SYMMETRY Z+
/
SET U 10.0
SET U at X- 10.0
/
DENSity 1.0
VISC 1.0
/
/OBJECT SWIRLER vane angle 45. 'swirler.loc' ID=SWIRL
/ rin 0.2 rout 0.2 hub_rin 0.01 hub_rout 0.01
/ (1., 0.5, 0.5) (2., 0.5, 0.5)
/
LOCAte ID=BODY (9,10,10) TO (12,13,13)
BLOCK ID=BODY
/
LOCAte MATCH X- of ID=BODY and x+ of ID=BODY ID=SWRL 'BODY_LOC.TMP'
/
TRANSFER FLUX ID=SWRL
SWIRL FLUX MODE: ID=SWRL VANE ANGLE = 45. CENTER = (5.0 0.5 0.5)
/
RELA U=0.8, p=0.5 T=1.0
/
DIAGnostic output U V W P RP BP at (13,11,11) every 1 steps
/
CONV GLOB 1.E-8
/
LOCAte (12,9,9) TO (12,14,14)
OUTPUT U V W P NARROW YZ PLANE SELECT
/
SOLVE STEADY 10
/
SAVE OFF U V W P T on file 'T16_SAV.TMP'
/
END
/

-Back To Top

************************************************************************
TITLE T.017: Testing of TRANSFER SOURCE option with flow
************************************************************************
GRID NODEs 22 by 22 BY 22
DEBUG GEOMETRY OFF
COOR X RANGE 10.0
COOR Y RANGE 1.0
COOR Z RANGE 1.0
/ Initial and Boundary Conditions
INLET X-
OUTLet X+
SYMMETRY Y-
SYMMETRY Y+
SYMMETRY Z-
SYMMETRY Z+
/
SET U 10.0
SET U at X- 10.0
/
DENSity 1.0
VISC 1.0
/
/OBJECT SWIRLER vane angle 45. 'swirler.loc' ID=SWIRL
/ rin 0.2 rout 0.2 hub_rin 0.01 hub_rout 0.01
/ (1., 0.5, 0.5) (2., 0.5, 0.5)
/
LOCAte ID=BODY (9,10,10) TO (12,13,13)
BLOCK ID=BODY
/
LOCAte MATCH X- of ID=BODY and x+ of ID=BODY ID=SWRL 'BOBY_LOC.TMP'
/
TRANSFER SOURCE 0.394132 TOTAL ID=SWRL
SWIRL SOURCE MODE: ID=SWRL VANE ANGLE = 45. CENTER = (5.0 0.5 0.5)
/
RELA U=0.8, p=0.5 T=1.0
/
DIAGnostic output U V W P RP BP at (13,11,11) every 1 steps
/
CONV GLOB 1.E-8
/
LOCAte (12,9,9) TO (12,14,14)
OUTPUT U V W P NARROW YZ PLANE SELECT
/
SOLVE STEADY 10
/
SAVE OFF U V W P T on file 'T17_SAV.TMP'
/
END
/

-Back To Top

************************************************************************
TITLE T.018: Testing of TRANSFER SOURCE Solid body
************************************************************************
GRID NODEs 22 by 22 BY 22
DEBUG GEOMETRY OFF
COOR X RANGE 10.0
COOR Y RANGE 1.0
COOR Z RANGE 1.0
/ Initial and Boundary Conditions
INLET X-
OUTLet X+
SYMMETRY Y-
SYMMETRY Y+
SYMMETRY Z-
SYMMETRY Z+
/
SET U 10.0
SET U at X- 10.0
/
DENSity 1.0
VISC 1.0
/
/OBJECT SWIRLER vane angle 45. 'swirler.loc' ID=SWIRL
/ rin 0.2 rout 0.2 hub_rin 0.01 hub_rout 0.01
/ (1., 0.5, 0.5) (2., 0.5, 0.5)
/
LOCAte ID=BODY (9,10,10) TO (12,13,13)
BLOCK ID=BODY
/
LOCAte MATCH X- of ID=BODY and x+ of ID=BODY ID=SWRL 'BODY_LOC.TMP'
/
TRANSFER SOURCE 0.394132 TOTAL ID=SWRL
SWIRL SOURCE ID=SWRL SOLID BODY 5. REV/SEC CENTER = (5.0 0.5 0.5)
/
RELA U=0.8, p=0.5 T=1.0
/
DIAGnostic output U V W P RP BP at (13,11,11) every 1 steps
/
CONV GLOB 1.E-8
/
LOCAte (12,9,9) TO (12,14,14)
OUTPUT U V W P NARROW YZ PLANE SELECT
/
SOLVE STEADY 10
/
SAVE OFF U V W P T on file 'T18_SAV.TMP'
/
END
/

-Back To Top

************************************************************************
TITLE T.019: Testing of multivariate SET LINEAR and SUM
************************************************************************
GRID NODEs 12 by 12 BY 12
DEBUG GEOMETRY OFF
/ Initial and Boundary Conditions
INLET X-
OUTLet X+
/
VISC 0.01
SET U X- LINEAR 1. * Y +5
SET V X- LINEAR -1. * Y -5
SET W X- 10.
/
LOCAte ID=BODY (1,1,1) TO (6,6,6)
/
ALLOCATE VAR1
ALLOCATE VAR2
ALLOCATE VAR3
ALLOCATE VAR4
/
SET VAR1 id=body: SUM of X + Y + Z ALWAYS
SET VAR2 id=body: LINEAR of 3 * X + 3 * Y + 3 * Z ALWAYS
SET VAR3 id=body: LINEAR 3. * VAR1 -1. * VAR2 ALWAYS
SET VAR4 id=body: LINEAR 3. * U +3. * V + 1. * W -10. ALWAYS
/
DIAGnostic output U V W P RP BP at (5,5,5) every 1 steps
/
CONV GLOB 1.E-8
/
SELECT (1,1,5) TO (999,999,7) with frequency of (2,2,1)
OUTPUT U V W P VAR1 VAR2 VAR3 VAR4 NARROW XY PLANES SELECT
/
SOLVE STEADY 10
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.020: Testing OPEN Option - Inlet at one end; fixed p at other
************************************************************************
GRID NODEs 22 by 12
DEBUG GEOMETRY OFF
COOR X RANGE 20
COOR Y RANGE 1
DEBUG GEOMETRY OFF
/
WALL BY DEFAULT
INLET X-
/
OPEN X+
BOUND X+ P = 0.
/BOUND X+ U GRAD = 0.
/
SET U X- 1.
DENSity 1.0
VISC 0.01
/
RELAX U= 0.5
/
DIAGnostic output for U V P BP at (21,6) every 30 steps
/
FLUX DEFAult output OFF
/
SELECT (1,1) (999,999) (4,1)
OUTPut U PTOTAL NARROW SELECted
/
CONV GLOB 1.E-10
SOLVE STEADY 1000
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.021: Testing OPEN option - Fixed P at both ends
************************************************************************
GRID NODEs 22 by 12
DEBUG GEOMETRY OFF
COOR X RANGE 20
COOR Y RANGE 1
DEBUG GEOMETRY OFF
/
WALL BY DEFAULT
/
OPEN X-
BOUND X- P = 1. TOTAL
/BOUND X- GRAD U = 0.
/
OPEN X+
BOUND X+ P = 0.
/BOUND X+ GRAD U = 0.
/
DENSity 1.0
VISC 0.01
/
DIAGnostic output for U V P BP at (21,6) every 30 steps
/
FLUX DEFAult output OFF
/
SELECT (1,1) (999,999) (4,1)
OUTPut U PTOTAL NARROW SELECted
/
CONV GLOB 1.E-10
SOLVE STEADY 1000
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.022: PROBlem ATMOsperic Option with a heat source and 2 outlets
************************************************************************
GRID is 22 by 22
DEBUG GEOMETRY OFF
COORdinate X RANGE 50.
COORdinate Y RANGE 50.
/
PROBLEM ATMO
/
ADIABATIC
WALL default
VISCosity is 1.0E-2
/
INIT T=293
/
LOCAte region as: ( 1,1) to ( 1, 5)
INLET X- SELECTED
SET U = 0.1 SELECTED
/
LOCATE (1,16) to (1, 22)
OUTLET X- SELECTED flow fraction 0.8
/
LOCAte region as: (22,16) to (22, 22)
OUTLET X+ SELECTED flow fraction 0.2
/
SELECT (11, 2) (11, 2)
SOURCE FLOW 0.25 U=0.05 V=0.05 T=294 SELECTED
/
GAS PRESSURE 101325
DENSITY GAS LAW
GRAVITY 0. -9.81
PRANDT 0.7
/
SPECific heat 1000.
/
LIMIT T MIN = 293 MAX =294
/
DIAGnostic U V P T RHO output at node (12,21) ; print every 100 step
/
RELAX U=0.2 P=0.8
MATRIX P 16
/
CONVERGENCE P 1.E-10 10
/
SELECT (1,1) TO (99,99) in steps of (4,2)
OUTPUT U V P Narrow in SELEcted region
CONVERGENCE 1.E-12
/
SOLVE in STEADY mode 1000
SAVE OFF
/
END
/

-Back To Top

*******************************************************************
TITLE T.023: Source Angles with Speed option and internal.
*******************************************************************
/
GRID 22 22 22
DEBUG GEOMETRY OFF
COORdinate X RANGe 1.0
COORdinate Y RANGe 1.0
COORdinate Z RANGe 1.0
/
WALL default
OUTLet X+
/
LOCAte ID=SRC1 from (2,22,5) to (4,22,16)
LOCAte ID=SRCA1 from (2,1,5) to (4,1,16)
LOCAte ID=SRC2 from (6,22,5) to (8,22,16)
LOCAte ID=SRCA2 from (6,1,5) to (8,1,16)
LOCAte ID=SRC3 from (10,22,5) to (12,22,16)
LOCAte ID=SRCA3 from (10,1,5) to (12,1,16)
LOCAte ID=SRC4 from (14,22,5) to (16,22,16)
LOCAte ID=SRCA4 from (14,1,5) to (16,1,16)
/
/
// Set 1: Typically used for multihole injection.
/
SOURce ANGLes AREA TOTAl ID=SRC1 Y+ 3.24 kg/s
deflection = 30 degrees rotation = 0 degrees axis = (1 0 0)
SPEEd = 1000 m/s INTErnal
/
/ Almost equivalent source on bottom boundary.
SOURce VOLUme FLOW TOTAl ID=SRCA1 3.24 kg/s U=500.0 V=866.0 W=0.0
/
// Set 2:
/
/ Expect velocity at boundary face: U = 500 m/s, V = -866 m/s,
/ and rapid decrease in magnitude with distance from top boundary.
SOURce ANGLes AREA TOTAl ID=SRC2 Y+ 3.24 kg/s
deflection = 30 degrees rotation = 0 degrees axis = (1 0 0)
SPEEd = 1000 m/s
/
/ Almost equivalent source on bottom boundary.
SOURce AREA FLOW TOTAl ID=SRCA2 Y- 3.24 kg/s U=500.0 V=866.0 W=0.0
/
// Set 3:
/
/ Expect a velocity of U = 17.32, V = -30 m/s at the top boundary.
SOURce ANGLes AREA TOTAl ID=SRC3 Y+ 3.24 kg/s
deflection = 30 degrees rotation = 0 degrees axis = (1 0 0)
/
/ Expect a velocity of U = 15, V = +25.98 m/s at the bot. boundary.
SOURce MOMEntum AREA TOTAl ID=SRCA3 Y- 3.24 kg/s NORMAls = (0.5,0.866,0)
/
/ Set 4: Typically used for dilution holes, Inlets etc.
/
/ Expect a velocity of U = 17.32, V = -30 m/s at the top boundary.
FLOW ANGLes AREA TOTAl ID=SRC4 Y+ 3.24 kg/s
deflection = 30 degrees rotation = 0 degrees axis = (1 0 0)
/
/ Expect a velocity of U = 15, V = +25.98 m/s at the bot. boundary.
FLOW MOMEntum AREA TOTAl ID=SRCA4 Y- 3.24 kg/s NORMAls = (0.5,0.866,0)

DENSity 1.2
VISCosity 1.5E-03

DIAGnostics (6,21,10) print U V W P every 10

RELAxation U 0.4 V 0.4 W 0.4 P 0.7

MATRix SWEEps U 2 V 2 W 2 P 5
CONVergence GLOBal U 1.0E-06 REFErence
CONVergence GLOBal P 1.0E-06 5

LOCAte ID=SLICE1 from (1,1,10) to (22,22,10) with frequency of (2,2,1)
OUTPUT U V W P ID=SLICE1
/
SOLVE U V W P 100 STEPs STEADY mode ! 500 steps for full
/
SAVE OFF
/
END
/

-Back To Top

*******************************************************************
TITLE T.024: Source Momentum and Angle with Speed option.
*******************************************************************
/
GRID 22 22 22
COORdinate X Y Z BLOCk format 1 header line from 'ANST24.GRD'

WALL default
OUTLet X+

LOCAte ID=SRC1 from (3,22,5) to (5,22,16)
LOCAte ID=SRC2 from (7,22,5) to (9,22,16)
LOCAte ID=SRC3 from (11,22,5) to (13,22,16)
/
SOURce ANGLes AREA TOTAl ID=SRC1 Y+ 3.24 kg/s
deflection = 30 degrees rotation = 0 degrees axis = (1 0 0)
SPEEd = 1000 m/s

SOURce ANGLes AREA TOTAl ID=SRC2 Y+ 3.24 kg/s
deflection = 30 degrees rotation = 0 degrees axis = (1 0 0)

FLOW ANGLes AREA TOTAl ID=SRC3 Y+ 3.24 kg/s
deflection = 30 degrees rotation = 0 degrees axis = (1 0 0)

DENSity 1.2
VISCosity 1.5E-03

DIAGnostics (6,21,10) print U V W P every 10

RELAxation U 0.3 V 0.3 W 0.3 P 0.6

MATRix SWEEps U 2 V 2 W 2 P 2
CONVergence GLOBal U 1.0E-06 REFErence
CONVergence GLOBal P 1.0E-06 5

LOCAte ID=SLICE1 from (1,1,10) to (22,22,10) with frequency of (2,2,1)
OUTPUT U V W P ID=SLICE1
/
SOLVE U V W P 100 STEPs STEADY mode ! 500 steps for full
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.025: Test Output of DIffusive, Convective and TOTAL fluxes
************************************************************************
/
GRID NODEs 22 by 22
DEBUG GEOMETRY OFF
/
COOR X RANGE = 10.
WALL at undefined outer boundaries
/
SYMM Y-
INLET X-
OUTLET X+
/
SET U 1
SET T X- 100
/
DENSity 1.0
VISC 0.07
PRANDTL NUMEBR EFFECTIVE 1.
/
DIAGnostics at (12,12) every 20 Steps
SELECT (1,1) TO (999,999) (4,4)
OUTPut U TOTAL FLUX ! at step 250
OUTPut T U DIFFUSIVE FLUX ! at step 250
OUTPut U V T CONVECTIVE FLUX ! at step 250
SAVE T U DIFFUSIVE FLUX ! at step 250
OUTPut U T for SELEcted subdomain
/
CONV U REFE 1.E-5
SOLVE U V P T STEADY 1000
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.026: Test FIX OFF
************************************************************************
/
GRID NODEs 22 by 22
DEBUG GEOMETRY OFF
/
COOR X RANGE = 10.
WALL at undefined outer boundaries
/
SYMM Y-
INLET X-
OUTLET X+
/
/ Fluid Properties and Constants
SET U 1
SET T X- 100
/
Select ID=MIDDLE (12,12) TO (13,13)
FIX SELECT U
FIX SELECT T
SET U 0.5 SELECT
SET T 49. SELECT
/
DENSity 1.0
VISC 0.07
PRANDTL NUMEBR EFFECTIVE 1.
/
DIAGnostics at (12,12) every 20 Steps
/
SELECT (1,1) (999,999) (5,5)
OUTPut U T NARROW SELECT
/
CONV U REFE 1.E-5
SOLVE U V P T STEADY 50
FIX T ID=MIDDLE OFF
FIX U ID=MIDDLE OFF
/
SOLVE U V P T STEADY 50
FIX T ID=MIDDLE
FIX U ID=MIDDLE
/
SOLVE U V P T STEADY 50
FIX T ID=MIDDLE OFF
/
SOLVE U V P T STEADY 500
SAVE OFF
/
END
/

-Back To Top

***********************************************************************
TITLE T.027: COMPRESSIBLE FLOW: CLOSED BOX TEST
***********************************************************************
GRID 12 by 3
DEBUG GEOMETRY OFF
COORDINATE X RANGE = 10
WALLS BY DEFAULT
/
COMPRESSIBLE FLOW !ISENTROPIC
LAMINAR FLOW
ADIABATIC WALLS
SET T = 300.
SET H = 300000.
INITIAL P = 9.E4
/
SELECT (1,1) TO (6,99)
SET P = 9.E5 SELECT
/
DENSITY GAS LAW COMPRESSIBLE MODE
GAS molecular weight = 27.71333333 ! to get R/mwt = 300
GAS BASE PRESSURE 0.
SPECIFIC HEAT 1000
VISCOSITY 1.E-10
/
DIAGNOSTIC NODE TIME U T RHO P AT 7,2 PRINT 20 STEPS 'DIAG_EQN.TMP'
CONVERGENCE FLOW 1.E-12 5
CONVERGENCE PP 1.E-12 5
CONVERGENCE T 1.E-12 2
CONVERGENCE U 1.E-12 2
/
OUTPut U P T RHO MACH in NARRow mode
/
SOLVE U V P T for 0.01 in STEPS OF 1.E-4
SAVE OFF
/
END
/

-Back To Top

***********************************************************************
TITLE T.028: Compressible Transient: I/O bc's high delta_P at t0
***********************************************************************
GRID is 12 by 3
COORdinate X RANGE 10
DEBUG GEOMETRY OFF
/
COMPRESSIBLE !ISENTROPIC
VISCosity is 1.5E-6
SPECIFIC HEAT UNIFORM = 1000
/
INLET X-
OUTLET X+
/
INIT U = 100
SET P = 1.0E+5
INITIAL H = 1.0E+6
INITIAL T = 1.0E+3
/
BOUNDARY for U X- 1000
BOUNDARY for P X- 1.0E+7
BOUNDARY for T X- 1000
/
DENSITY GAS LAW COMPRESSIBLE MODE
GAS BASE PRESSURE 0.
GAS CONSTANT = 1000
GAS molecular weight = 10 ! To get rho=1 at p=10^5, T=10^3
/
DIAGnostic TIME U V P T RHO output at node (12,2); print every 20 step
/
LIMIT T MIN = 200 MAX 2000
CONVERGENCE FLOW 1.E-8 10
CONVERGENCE PP 1.E-8 10
/
OUTPUT U P PTOTAL T RHO FC MACH NARROW EVERY ! 50 STEPS
/
SOLVE U V P T for 0.1 in STEPS OF 0.001
SAVE OFF
/
END
/

-Back To Top

***********************************************************************
TITLE T.029: Compressible Steady: I/O bc's high delta_P at t0
***********************************************************************
GRID is 12 by 3
COORdinate X RANGE 10
DEBUG GEOMETRY OFF
/
COMPRESSIBLE !ISENTROPIC
VISCosity is 1.5E-6
SPECIFIC HEAT UNIFORM = 1000
disable global
/
INLET X-
OUTLET X+
/
INIT U = 1000
SET P = 1.0E+6
INITIAL H = 1.0E+6
INITIAL T = 1.0E+3
/
BOUNDARY for U X- 1000
BOUNDARY for P X- 1.0E+7
BOUNDARY for T X- 1000
/
DENSITY GAS LAW COMPRESSIBLE MODE
GAS BASE PRESSURE 0.
GAS CONSTANT = 1000
GAS molecular weight = 10 ! To get rho=1 at p=10^5, T=10^3
/
DIAGnostic TIME U V P T RHO output at node (2,2); print every 20 step
/
CONVERGENCE FLOW 1.E-10 1
CONVERGENCE PP 1.E-10 10
/
OUTPUT U P PTOTAL T RHO FC MACH NARROW EVERY ! 50 STEPS
/
LIMIT T MIN 500 MAX 1500
LIMIT RHO MIN 1 MAX 200
SOLVE U V P T STEADY 200
SAVE OFF
/
END
/

-Back To Top

***********************************************************************
TITLE T.030: Compressible Steady: OPEN/OPEN bc's (Transient/Steady)
***********************************************************************
GRID is 12 by 3
COORdinate X RANGE 10
DEBUG GEOMETRY OFF
/
COMPRESSIBLE ! ISENTROPIC
VISCosity is 1.0E-16
SPECIFIC HEAT UNIFORM = 1000
/
OPEN X-
OPEN X+
/
SET U = 200
SET P = 1.0E+7
SET H = 1.0E+6
SET T = 1.0E+3
/
BOUNDARY for P X- 1.0E+7 TOTAL
BOUNDARY for T X- 1000
BOUNDARY for P X+ 0.5E+7
/BOUNDARY for P X+ 1.0E+5 ! For supersonic
BOUNDARY for H X+ GRAD = 0.
/
DENSITY GAS LAW COMPRESSIBLE MODE
GAS BASE PRESSURE 0.
GAS CONSTANT = 1000
GAS molecular weight = 10 ! To get rho=1 at p=10^5, T=10^3
/
DIAGnostic TIME U V P T RHO output at node (1,2); print every 200 step
/
CONVERGENCE FLOW 1.E-12 10
CONVERGENCE PP 1.E-12 2
/
OUTPUT U P PTOTAL T H RHO MACH NARROW EVERY ! 50 STEPS
/
LIMIT T MIN 500 MAX 1500
LIMIT RHO MIN 1 MAX 200
/SOLVE U V P STEADY 2000 MIN 2000
SOLVE U V P T 20 0.01
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.031: Compressible Transient (Ma=0.675) Flow Over a 4% Bump
************************************************************************
/
GRID NODEs 32 by 12
COOR X Y 'ANST31.GRD' BLOCK FORMAT header 2 LINES
/
/ Initial and Boundary Conditions
COMPRESSIBLE !ISENTHALPIC FLOW
SPECIFIC HEAT UNIFORM = 1000
/
WALL at undefined outer boundaries
OPEN X-
OPEN X+
/OUTLET X+ ! Also works but time step needs to be 1.E-4
/
SET P = 2.00E5
SET T = 300
/
ADIABATIC WALLS
/
BOUNDARY P at X- 2.00E+5 TOTAL
BOUNDARY V at X- 0.
BOUNDARY P at X+ 1.509E+5 ! Mach # = 0.675
/
/BOUNDARY P at X+ 1.798E+5 ! Mach # = 0.400
/BOUNDARY P at X+ 1.750E+5 ! Mach # = 0.450
/BOUNDARY P at X+ 1.705E+5 ! Mach # = 0.496
/BOUNDARY P at X+ 1.340E+5 ! Mach # = 0.800
/BOUNDARY P at X+ 0.800E+5 ! Mach # = 0.850
/
DENSITY GAS LAW COMPRESSIBLE MODE
GAS BASE PRESSURE 0.
GAS molecular weight = 29
/
VISCOSITY 1.E-10 ! Essentially inviscid
/
DIAGnostic output U V P T BP at (26,2) every 50 steps !just after bump
/
LIMIT T MINImum 10 max 2000
LIMIT RHO MINImum 0.01 max 10
LIMIT P MINImum 1.E2 max 1.E6
MATRIX P 20
/CONV GLOB 1.E-16 20 ! Improves residue for U,W,T equations
/CONV PP 1.E-10 10 ! Improves residue for mass (P) equation
/
SELECt (1,1) to (999,999) with frequency of (6,2)
OUTPut U V P PTOTAL RHO MACH in NARROW mode in SELEcted region
SAVE U V P T RHO PTOTAL MACH
/
SOLVE U V P T 0.5 1.E-3
SAVE OFF
/
END
/
/

-Back To Top

************************************************************************
TITLE T.032: Compressible Steady (Ma=1.65) Flow over a 4% Bump in Channel
************************************************************************
GRID NODEs 32 by 12
COOR X Y 'ANST31.GRD' BLOCK FORMAT header 2 LINES
/
COMPRESSIBLE !ISENTHALPIC FLOW
SPECIFIC HEAT UNIFORM = 1000
/
WALL at undefined outer boundaries
INLET X-
OUTLET X+
/
SET P = 2.00E5
SET T = 300
/
ADIABATIC WALLS
BOUNDARY P at X- 2.00E+5
BOUNDARY U at X- 573
/
DENSITY GAS LAW COMPRESSIBLE MODE
GAS BASE PRESSURE 0.
GAS molecular weight = 29
/
VISCOSITY 1.E-8
DIAGnostic output U V P T BP at (22,11) every 50 steps
/
LIMIT T MINImum 10 max 2000
LIMIT RHO MINImum 0.01 max 10
LIMIT P MINImum 1.E2 max 1.E6
MATRIX U 2 V 2 T 2 P 16
CONV GLOB 1.E-12 1
CONV PP 1.E-10 20
/
SELECt (1,1) to (999,999) with frequency of (6,2)
OUTPut U V P PTOTAL RHO MACH in NARROW mode in SELEcted region
SAVE U V P T RHO PTOTAL MACH
/
SOLVE STEADY U V P T 1000 500
SAVE OFF
/
END
/
/

-Back To Top

************************************************************************
TITLE T.033: Testing SET DISTANCE, SET SQUARE & SET ROOT SQUARE Options
************************************************************************
/
GRID NODEs 22 by 22
DEBUG GEOMETRY OFF
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
/
WALLs by default
INLET X-
OUTLET X+
SYMMETRY Y-
/
LOCATE ID=Subdmn01 from (1, 1) to (999, 5)
LOCATE ID=Subdmn02 from (1, 6) to (999, 7)
LOCATE ID=Subdmn03 from (1, 8) to (999, 9)
LOCATE ID=Subdmn04 from (1, 10) to (999, 11)
LOCATE ID=Subdmn05 from (1, 12) to (999, 13)
LOCATE ID=Subdmn06 from (1, 14) to (999, 15)
LOCATE ID=Subdmn07 from (1, 16) to (999, 99)
/
SET U = 1 for X- boundary of domain
/
SET DISTANCE based 3 sets ID=Subdmn01 ALWAYS
X Y T C
(0.2, 0.1) 100 0.1
(2.0, 0.5) 200 0.2
(8.0, 0.9) 500 1.0
/
SET DISTANCE 3 sets ID=Subdmn02 ALWAYS (X,Y,T,C) file='ANST33.FIL'
/
SET LINEar function T = X + Y ID=Subdmn03 ALWAYS
SET LINEar function T = 10. * X + 100. * Y ID=Subdmn04 ALWAYS
/
SET Sum of SQUAres function T = X Y ID=Subdmn05 ALWAYS
SET ROOT of Sum of SQUAres function T = X Y ID=Subdmn06 ALWAYS
/
SET Sum of SQUAres function T = 10 X 10 Y +100 ID=Subdmn07 ALWAYS
SET ROOT SQUAres function T = 10 X +10 Y ID=Subdmn07 ADD ALWAYS
/
DENSity 1.0
VISC 1.E-3
/
DIAGnostic output for U V P T C at (11,2) every 1 steps
/
FLUX DEFAult output OFF
SELECt (1,1) to (999,999) with frequency of (4,2)
OUTPut T C in NARROW mode in SELEcted region
/
SOLVE U V W P STEADY 10
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.034: Testing PRINT AVERAGE Commands and Options
// These commands can be used to generate output of the integrated
// value of a variable at a number of locations. If a 3D subdomain is
// specified, then integration is performed plane by plane at a number
// of locations. If a plane is specified, then the integration is
// performed at a number of strips in one direction.
************************************************************************
/
GRID NODEs 12 by 9 BY 7
DEBUG GEOMETRY OFF
COOR X RANGE 10. 1.05
WALLs by default
INLET X-
OUTLET X+
SYMMETRY Y-
/
SET U = 10 for X- boundary of domain
SET LINEar function T = Z + Y
SET LINEar function T = 10. * Z + 100. * Y add
SET T Y+ 0.
SET T Y- 100
SET T Z+ FLUX = 0.
SET T Z- FLUX = 0.
/
DENSity 1.0
VISC 1.E-3
DIAGnostic output for U V P T C at (11,7,6) every 10 steps
/
/======================================================================/
**** Output# 1
***>>> Pattern & Profile factor (T vs. R) at outlet
* Specify T3=10, T4=arithmetic mean
* Print to standard output

PRINT PROFile T at 5 locations at X+ boundary to standard output
at R = 0.1 0.3 0.5 0.7 0.9 TOLErance 0.1 BASE value=10.
/======================================================================/
/
**** Output# 2
***>>> Same as #1 but T4 as mass-avearged value

PRINT PROFile T at 5 locations at X+ boundary to standard output
at R = 0.1 0.3 0.5 0.7 0.9 TOLErance 0.1, BASE=10. MASS averaged
/
/======================================================================/
/
**** Output# 3
***>>> Same as #2; but send output to a file called T34_AVG.TMP

PRINT PROFile T at 5 locations at X+ boundary to file 'T34_AVG.TMP'
at R = 0.1 0.3 0.5 0.7 0.9 TOLErance 0.1, BASE=10. MASS averaged
/
/======================================================================/
/
**** Output# 4
***>>> Write T (Averaged Y & Z) vs X at 5 axial locations
* Print to standard output every 100 steps

PRINT AVERage in X of T at 5 locations every 100 steps standard output
/
/======================================================================/
/
**** Output# 5
***>>> Same as #4 but output to acrAVERAGE.TMP every 100 steps

PRINT AVERage in X of T at 5 locations to file 'acr_AVERAGE.TMP'
/
/======================================================================/
/
**** Output# 6
***>>> Compute T_mean, T_min, T_max, T_s.d. at 5 axial locations
***>>> Axial planes (yz) will be averaged
* Print to standard output

PRINT AVERage STATISTICS X of T at 5 locations to standard output device
/
/======================================================================/
**** Output# 7
***>>> Compute T vs. Theta at outlet; large tolerance
* Print to previous file

PRINT PROFile T at 5 locations at X+ boundary of domain
at THETA = 0.1 0.3 0.5 0.7 0.9 TOLErance 0.5 to 'acr_AVERAGE.TMP'
/
/======================================================================/
**** Output# 8
***>>> Non-dimensionalized T vs. Y (averged in Z) at outlet
* Print to previous file every 200 steps

PRINT PROFile T vs Y at 5 locations NORMalize X+
every 200 steps to 'acr_AVERAGE.TMP'
/
/======================================================================/
**** Output# 9
***>>> Same as #8
* Output to new file

PRINT AVERage T vs Y at 5 locs NORMalize X+ every 200 steps 'T_NORM.TMP'
/
/======================================================================/
**** Output# 10
***>>> U and T vs. X (averged in Y & Z) at 30 axial locations
* Any missing values computed by LINEar interpolation
* Output to previous file

PRINT PROFile for U and T at 15 X locations LINEAr 'acr_AVERAGE.TMP'
/
/======================================================================/
**** Output# 11
***>>> U and T vs. X (averged in Y & Z) at 30 axial locations
* Any missing values computed by 1/r^2 interpolation
* Output to previous file
PRINT PROFile for U and T at 15 X locations SQUAre 'acr_AVERAGE.TMP'

/
/======================================================================/
**** Output# 12
***>>> Same as 11 but very narrow tolerance to get "thin slice" planes
* Any missing values computed by SQUARE interpolation
* Output to previous file

PRINT PROFile U, T at 15 X TOLE=0.005 SQUAre interpolation 'acr_AVERAGE.TMP'
/
/======================================================================/
/
/
CONVERGENCE 1.E-4
FLUX DEFAult output OFF
SELECT (1,1,2) TO (99,99,999) at frequency (3,3,5)
OUTPut U T NARROW SELECTED
/
SOLVE U V W P T STEADY 1000
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.035: Testing HISTory SOURce and STORage options
// These commands can be used to generate output of the history of
// the SOURce or STORage (accumulation) term in the transport equation
// as a function of time. The output is directed to a file; also
// tables for each variable are generated at the end of simulation.
************************************************************************
/
DEBUG GEOMETRY OFF
GRID NODEs 12 by 12
COOR X RANGE 10. 1.05
WALLs by default
INLET X-
OUTLET X+
SYMMETRY Y-
/
LOCATE ID=MIDDLE from (7,3) to (9,7)
/
SET U = 10 for X- boundary of domain
SET LINEar function T = 100 + 100 * Y
SET T Y+ 0.
SET T Y- 100
/
DENSity 1.0
VISC 1.E-3
DIAGnostic output for U V P T C at (11,7) every 20 steps
/
/======================================================================/
**** Output# 1
***>>> History of changes in accumulative term (storage) of U

HISTory STORage for U t EVERY 30 steps on file 'T35_STOR.TMP'
/======================================================================/
/
**** Output# 2
***>>> History of source term of U for subregion
/
HISTory of SOURCE for U EVERY 50 steps for ID=MIDDLE on file 'T35_SOUR.TMP'
/
/======================================================================/
/
**** Output# 3
***>>> History of changes in accumulative term (storage) of T

HISTory STORage for T EVERY 25 steps ID=MIDDLE region
/======================================================================/
/
**** Output# 4
***>>> Get History of source term of T for subregion
/
HISTory of SOURCE for T every 100 step
/======================================================================/
/
CONVERGENCE 1.E-8
FLUX DEFAult output OFF
OUTPut off
/
SOLVE U V P T STEADY 1000
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.036: Testing HISTory COORdinate Command options
// These commands can be used to generate output of the history of
// the values of variable at arbitrary locations which are specified
// by their (x,y,z) coordinates. Values are interpolated to these
// locations from the computed grid nodes by 1/r or 1/r^2 algorithm.
************************************************************************
/
DEBUG GEOMETRY OFF
GRID NODEs 12 by 12
COOR X RANGE 10. 1.05
WALLs by default
INLET X-
OUTLET X+
SYMMETRY Y-
/
LOCATE ID=MIDDLE from (7,3) to (9,7)
/
SET U = 10 for X- boundary of domain
SET LINEar function T = 100 + 100 * Y
SET T Y+ 0.
SET T Y- 100
/
DENSity 1.0
VISC 1.E-3
DIAGnostic output for U V P T C at (11,7) every 20 steps
/
/======================================================================/
**** Output# 1
***>>> History at fixed location Linear interpolation

HISTory U, P, T at COORdinates (0.035, 0.577) every 30 steps on 'T36_XYZ.TMP'
/======================================================================/
/
**** Output# 2
***>>> History at fixed location inverse Square interpolation
/
HISTory U, P, T at COORdinates (0.035, 0.577) SQUAre method 30 steps
/
/======================================================================/
/
CONVERGENCE 1.E-8
FLUX DEFAult output OFF
OUTPut OFF
/
SOLVE U V P T STEADY 1000
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.037: Testing SET NODE multiple variable Options; file read
// Input both direct and from file read
************************************************************************
/
DEBUG GEOMETRY OFF
GRID NODEs 12 by 12
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
/
WALLs by default
INLET X-
OUTLET X+
SYMMETRY Y-
/
LOCATE ID=Subdmn01 from (1, 1) to ( 5, 5)
LOCATE ID=Subdmn02 from (1, 1) to ( 5, 5) FIELD ONLY
LOCATE ID=Subdmn03 LIST IJK (2, 8) (2,9) (2,10) (2,11)
LOCATE ID=Subdmn07 from (7, 7) to (10, 10)
/
SET NODE ID=Subdmn01 X-
U SKIP T C
10 -1 100 0.1
20 -1 200 0.2
30 -1 300 0.3
40 -1 400 0.4
/
SET NODE ID=Subdmn02
U SKIP T C
15 -1 110 0.1
16 -1 120 0.2
17 -1 130 0.3
18 -1 140 0.4
15 -1 110 0.1
16 -1 120 0.2
17 -1 130 0.3
18 -1 140 0.4
15 -1 110 0.1
16 -1 120 0.2
17 -1 130 0.3
18 -1 140 0.4
15 -1 110 0.1
16 -1 120 0.2
17 -1 130 0.3
18 -1 140 0.4
/
SET NODE ID=Subdmn03 X-
U SKIP T C
10 -1 100 0.1
20 -1 200 0.2
30 -1 300 0.3
40 -1 400 0.4
/
SET NODE ID=Subdmn07 V T SKIP SKIP C on file 'ANST37.FIL'
/
DENSity 1.0
VISC 1.E-3
/
DIAGnostic output for U V P T C at (11,2) every 1 steps
/
SELECT from (1,1) to (999,999) with frequency of (2,1)
OUTPut U V T C NARROW SELECTED
/
FLUX DEFAult output OFF
/
SOLVE OFF U V W P STEADY 10
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.038: Testing PLOT3D File (Formatted) Input
************************************************************************
/
GRID PLOT3D 'ANST38.GRD' !! GRID is 39 20
/
WALL BY DEFAULT
INLET X-
OUTL X+
SYMM Y-
/
SET U at X- boundary 1
SET T at X- boundary 100
/
DENSity 1.0
VISC 1.0
SPEC 1.0
PRAN EFFE 0.7
/ SOLUTION OPTIONS
MATRIX P = 128
DIAGnostic output U V P T RP at (4,11) every 10 steps
CONVERGENCE GLOBAL 1.E-8
/
FLUX DEFAult output OFF
/
SELECT from (1,1) to (999,999) with frequency of (6,2)
OUTPUT in Narrow mode for SELEcted region
SOLVE U V P T STEADY 500
SAVE OFF
/
END
/

-Back To Top

*********************************************************************
TITLE T.039: Test 1 for swirler profile option.
*********************************************************************
/
DEBUG GEOMETRY OFF
GRID 22 x 22 x 22
COORdinate X RANGe = 2.0
COORdinate Y RANGe = 1.0
COORdinate Z RANGe = 1.0
/
LOCAte PAIR ID=RGN1
1983, 4 1984, 4 1985, 4 1986, 4 1987, 4 1988, 4 2383, 4 2384, 4 2385, 4
2386, 4 2387, 4 2388, 4 2783, 4 2784, 4 2785, 4 2786, 4 2787, 4 2788, 4
3183, 4 3184, 4 3185, 4 3186, 4 3187, 4 3188, 4 3583, 4 3584, 4 3585, 4
3586, 4 3587, 4 3588, 4 3983, 4 3984, 4 3985, 4 3986, 4 3987, 4 3988, 4
4383, 4 4384, 4 4385, 4 4386, 4 4387, 4 4388, 4 4783, 4 4784, 4 4785, 4
4786, 4 4787, 4 4788, 4 5183, 4 5184, 4 5185, 4 5186, 4 5187, 4 5188, 4
5583, 4 5584, 4 5585, 4 5586, 4 5587, 4 5588, 4 5983, 4 5984, 4 5985, 4
5986, 4 5987, 4 5988, 4 6383, 4 6384, 4 6385, 4 6386, 4 6387, 4 6388, 4
! BBox: (0.2,1,0.2) to (0.8,1,0.8)
/
LOCAte PAIR ID=RGN2
3361, 1 3381, 1 3401, 1 3421, 1 3761, 1 3781, 1 3801, 1 3821, 1 4161, 1
4181, 1 4201, 1 4221, 1 4561, 1 4581, 1 4601, 1 4621, 1
! BBox: (0,0.4,0.4) to (0,0.6,0.6)
/
LOCAte ID=RGN3 from (3,8,22) to (8,15,22)
! BBox: (0.1,0.3,1) to (0.7,0.7,1)
/
SET U 1.0E-06
SET V 0.0
SET W 0.0
SET T 0.0
SET K 1.0
SET L 0.1
/
WALLs by default
INLET ID=RGN1
INLET ID=RGN2
INLET ID=RGN3 Z+
OUTLET X+
/
SWIRler PROFile RADIal mdot=50.0 kg/s center=(0.5,0.5,0.5)
AXIS=(1,0,0) file='ANST39R.FIL'
additional variables in file: T K L
at ID=RGN1
detailed DIAGnostics on.
/
SWIRler PROFile mdot=10.0 kg/s center=(1.0,0.5,0.5)
SCALing (3.53553E-02,0.1060660) file='ANST39.FIL'
additional variables in file: T K L
at ID=RGN2
detailed DIAGnostics on.
/
SWIRler PROFile mdot=50.0 kg/s center=(0.4,0.5,1.0)
file='ANST39.FIL'
additional variables in file: T K L
at ID=RGN3 Z+
detailed DIAGnostics on.
/
SCREen WIDE
DIAGnostic NODE (1,11,11) print U V W P T K L every 1 iteration

CONVergence GLOBal 1.0E-06

SELEct (1,1,22) to (22,22,22) in frequency of (4,2,1)
OUTPUT U V W T K L in NARRow mode for SELEcted region

SOLVe U V W P K E H STEADy 5
/
WRITE U V W T K L ID=RGN1
WRITE U V W T K L ID=RGN2
PRINT FLOW and STATs
/
SAVE 'T39_SAV.TMP'
END
/

-Back To Top

*********************************************************************
TITLE T.040: Test 2 for swirler profile option.
*********************************************************************
/
DEBUG GEOMETRY OFF
GRID 22 x 22 x 22
COORdinate X RANGe = 2.0
COORdinate Y RANGe = 1.0
COORdinate Z RANGe = 1.0
/
LOCAte PAIR ID=RGN1
1983, 4 1984, 4 1985, 4 1986, 4 1987, 4 1988, 4 2383, 4 2384, 4 2385, 4
2386, 4 2387, 4 2388, 4 2783, 4 2784, 4 2785, 4 2786, 4 2787, 4 2788, 4
3183, 4 3184, 4 3185, 4 3186, 4 3187, 4 3188, 4 3583, 4 3584, 4 3585, 4
3586, 4 3587, 4 3588, 4 3983, 4 3984, 4 3985, 4 3986, 4 3987, 4 3988, 4
4383, 4 4384, 4 4385, 4 4386, 4 4387, 4 4388, 4 4783, 4 4784, 4 4785, 4
4786, 4 4787, 4 4788, 4 5183, 4 5184, 4 5185, 4 5186, 4 5187, 4 5188, 4
5583, 4 5584, 4 5585, 4 5586, 4 5587, 4 5588, 4 5983, 4 5984, 4 5985, 4
5986, 4 5987, 4 5988, 4 6383, 4 6384, 4 6385, 4 6386, 4 6387, 4 6388, 4
! BBox: (0.2,1,0.2) to (0.8,1,0.8)
/
LOCAte PAIR ID=RGN2
3361, 1 3381, 1 3401, 1 3421, 1 3761, 1 3781, 1 3801, 1 3821, 1 4161, 1
4181, 1 4201, 1 4221, 1 4561, 1 4581, 1 4601, 1 4621, 1
! BBox: (0,0.4,0.4) to (0,0.6,0.6)
/
LOCAte ID=RGN3 from (3,8,22) to (8,15,22)
! BBox: (0.1,0.3,1) to (0.7,0.7,1)
/
SET U 1.0E-06
SET V 0.0
SET W 0.0
SET T 0.0
SET K 1.0
SET L 0.1
/
WALLs by default
INLET ID=RGN1
INLET ID=RGN2
INLET ID=RGN3 Z+
OUTLET X+
/
SWIRler PROFile RADIal mdot=50.0 kg/s center=(0.5,0.5,0.5)
AXIS=(1,0,0) file='ANST39R.FIL'
additional variables in file: T K L
at ID=RGN1
ALWAys detailed DIAGnostics on.
/
SWIRler PROFile mdot=10.0 kg/s center=(1.0,0.5,0.5)
SCALing (3.53553E-02,0.1060660) file='ANST39.FIL'
additional variables in file: T K L
at ID=RGN2
ALWAys detailed DIAGnostics on.
/
SWIRler PROFile mdot=50.0 kg/s center=(0.4,0.5,1.0)
file='ANST39.FIL'
additional variables in file: T K L
at ID=RGN3 Z+
ALWAys detailed DIAGnostics on.
/

SCREen WIDE
DIAGnostic NODE (1,11,11) print U V W P T K L every 1 iteration

CONVergence GLOBal 1.0E-06

SELEct (1,1,22) to (22,22,22) in frequency of (4,2,1)
OUTPUT U V W T K L in NARRow mode for SELEcted region

SOLVe U V W P K E H STEADy 5
/
WRITE U V W T K L ID=RGN1
WRITE U V W T K L ID=RGN2
PRINT FLOW and STATs
/
SAVE 'T40_SAV.TMP'
END
/

-Back To Top

*********************************************************************
TITLE T.041: Test 3 for swirler profile option with Fuel Injection.
*********************************************************************
/
DEBUG GEOMETRY OFF
GRID 22 x 22 x 22
COORdinate X RANGe = 2.0
COORdinate Y RANGe = 1.0
COORdinate Z RANGe = 1.0
/
LOCAte PAIR ID=RGN1
1983, 4 1984, 4 1985, 4 1986, 4 1987, 4 1988, 4 2383, 4 2384, 4 2385, 4
2386, 4 2387, 4 2388, 4 2783, 4 2784, 4 2785, 4 2786, 4 2787, 4 2788, 4
3183, 4 3184, 4 3185, 4 3186, 4 3187, 4 3188, 4 3583, 4 3584, 4 3585, 4
3586, 4 3587, 4 3588, 4 3983, 4 3984, 4 3985, 4 3986, 4 3987, 4 3988, 4
4383, 4 4384, 4 4385, 4 4386, 4 4387, 4 4388, 4 4783, 4 4784, 4 4785, 4
4786, 4 4787, 4 4788, 4 5183, 4 5184, 4 5185, 4 5186, 4 5187, 4 5188, 4
5583, 4 5584, 4 5585, 4 5586, 4 5587, 4 5588, 4 5983, 4 5984, 4 5985, 4
5986, 4 5987, 4 5988, 4 6383, 4 6384, 4 6385, 4 6386, 4 6387, 4 6388, 4
! BBox: (0.2,1,0.2) to (0.8,1,0.8)
/
LOCAte PAIR ID=RGN2
3361, 1 3381, 1 3401, 1 3421, 1 3761, 1 3781, 1 3801, 1 3821, 1 4161, 1
4181, 1 4201, 1 4221, 1 4561, 1 4581, 1 4601, 1 4621, 1
! BBox: (0,0.4,0.4) to (0,0.6,0.6)
/
LOCAte ID=RGN3 from (3,8,22) to (8,15,22)
! BBox: (0.1,0.3,1) to (0.7,0.7,1)
/
SET U 1.0E-06
SET V 0.0
SET W 0.0
SET T 0.0
SET K 1.0
SET L 0.1
SET FU 1.0E-06
/
WALLs by default
INLET ID=RGN1
INLET ID=RGN2
INLET ID=RGN3 Z+
OUTLET X+
/
/
SWIRler PROFile RADIal mdot=50.0 kg/s center=(0.5,0.5,0.5) AXIS=(1,0,0)
FUEL mdotf=5 kg/s file='ANST41R.FIL'
additional variables in file: T K L FU
at ID=RGN1
ALWAys
detailed DIAGnostics on.
/
SWIRler PROFile mdot=10.0 kg/s center=(1.0,0.5,0.5)
FUEL mdotf=1 kg/s file='ANST41A.FIL'
additional variables in file: T K L FU
at ID=RGN2
ALWAys
detailed DIAGnostics on.
/
SWIRler PROFile mdot=50.0 kg/s center=(0.4,0.5,1.0)
FUEL mdotf=5 kg/s file='ANST41.FIL'
additional variables in file: T K L FU
at ID=RGN3 Z+
ALWAys
detailed DIAGnostics on.
/

SCREen WIDE
DIAGnostic NODE (1,11,11) print U V W P T K L every 1 iteration

CONVergence GLOBal 1.0E-06

SELEct (1,1,22) to (22,22,22) in frequency of (4,2,1)
OUTPUT U V W T K L FU in NARRow mode for SELEcted region

SOLVe U V W P K E H FU STEADy 5
/
WRITE U V W T K L FU ID=RGN1
WRITE U V W T K L FU ID=RGN2
PRINT FLOW and STATs
/
SAVE 'T41_SAV.TMP'
END
/

-Back To Top

*********************************************************************
TITLE T.042: Conjugate Heat Transfer: Variable CP: Pure Conduction
*********************************************************************
/
DEBUG GEOMETRY OFF
GRID NODEs 12 by 22
/
COOR X RANGE = 10.0
COOR Y RANGE = 1.0
/
WALL by default
SYMM Y-
/
LOCATE ID=BLOCK1 COORDINATES (0,0.7) to (1000.,0.8)
BLOCK ID=BLOCK1
LOCATE ID=BLOCK2 COORDINATES (0,0.8) to (1000.,0.9)
BLOCK ID=BLOCK2
LOCATE ID=BLOCK3 COORDINATES (0,0.9) to (1000.,1.1)
BLOCK ID=BLOCK3
/
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
SET T VALU LINEAR 100.0 -10. * X
/
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENS 1.0
VISC 0.01
Prandtl number for fluid effective value = 1
/
SOLID RHO=2 CP=2 KT=0.01
/
SPECIFIC HEAT FOR ID=BLOCK1 is POLYNOMIAL function in T ALWAYS
Coefficients are: 2., 0.1, -0.001, 0.000025
/
SPECIFIC HEAT for ID=BLOCK2 is uniform at 2
/
SPECIFIC HEAT for ID=BLOCK3 is POWER law 0.2 * ( X + 0. ) ^ 2.5 +2.
/
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS
RELA T 0.5
LIMI T -1.E30 1.E30
/
DIAGnostic output at (3,3) every 10 steps
SAVE OFF U V P T H CP on file 'T42_SAV.TMP'
/
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV REFErence value is T: criterion 1.E-3 innner iterations 1
SOLVE T STEADY 500
/
SELECT from (1,1) to (999,999) with frequency of (2,2)
OUTPut T, H, CP NARRow mode for SELEcted domain
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.043: Flow in Square Cavity with User-Defined Transport & Reaction
************************************************************************
// This data set for ACRi Software ANSWER illustrates the use of
// the FREEFORM command language to dynamically define:
// >> New Field Variables,
// >> New Transport Equations, and
// >> Complex Reactions and Sources
//
************************************************************************
/
// The problem is set in the context of the classic Square Cavity
// with a moving lid at the top.
//
// A new species call DOGS is defined on the fly which is governed by
// a transport equation. The transport properties and boundary conditions
// for DOGS are defined by simple commands.
//
// The distribution of DOGS is effected by the availability of MEAT
// as a sinusoidal function of Time, and a 1st-order Arrhenius reaction
// called KILLed. These are defined on-the-fly.
//
// MEAT = 2.0 * SIN ( 0.06886 * Time + 1. ) +2.0
//
// KILLed = -19.6 * exp(-937/T) DOGS^0.5
/
// In addition in one part of the domain, the DOGS are regenerated as
// as a function of U velocity if it is positive and eliminated if it is
// negative. The function is defined as a arbitrary tabulated function
************************************************************************
/
DEBUG GEOMETRY OFF ! Cartesian geom; no need to debug
GRID NODEs BY 22 BY 22 ! Use structured mode of ANSWER
/
COORDINATE X RANGE 1 ! Auto gridding; range specified
COORDINATE Y RANGE 1 ! Auto gridding; range specified
/
*****==========>>>>> Flow conditions
WALL at all outer boundaries
BOUNdary U: at boundary Y+, VALUE=1. ! Specify moving top wall
LAMINAR flow
VISCOSITY is 0.01
SET T (temperature) = 293.15
/
*****==========>>>>> Define new variables on the fly
ALLOCATe "DOGS- Members of Canine Family"
ALLOCATe "MEAT- Food for Dogs"
/
//////////////////////////////////////////////////////////////////////
***===>> Define Transport equation for DOGS and its attributes
//////////////////////////////////////////////////////////////////////
INITial value of DOGS = 10
PRANDtl Number for DOGS is 0.1

! Boundary conditions for DOGS including a linear function at top

BOUNdary X- (left) for DOGS: Value = 2 ! Fixed value
BOUNdary X+ (right) for DOGS: FLUX = 0 ! No DOGS leave or enter
BOUNdary Y- (bottom) for DOGS: FLUX = 1.E-4 ! dogs enter from bottom
BOUNdary Y+ (top) for DOGS: VALUe = Linear function 2. +8. * X

/
//////////////////////////////////////////////////////////////////////
// Define MEAT as sinusoidal function of TIME; make it source for DOGS
//////////////////////////////////////////////////////////////////////
SET MEAT = 2. * SIN ( 0.06886 * Time + 1. ) +2. ALWAys

SOURce for DOGS is LINEAR function: 0. +1. * MEAT per unit VOLUme

! Assume dogs die as a non-linear Arrhenius function of population

ALLOCATE 'KILL- Killed from various causes'
REACtion KILL DOGS Z0 = 19.6 TEMPrature_activation = 937 power = 0.5
SOURce REACTION type for DOGS LINEar -1. * KILL
/

! Assume in one subregion dogs die when U <0 & are born when U >0

SELECT subregion 5,5 to 7,7
SOURCE for DOGS for SELECted region is a TABLE of 3 sets in terms of U
// U SOURCE
-1 -0.01234
0. 0
1. 2.1234
//////////////////////////////////////////////////////////////////////
**** END of specification for DOGS' transport equation
//////////////////////////////////////////////////////////////////////
/
REACtion default hydrocarbon OFF
DIAGNOSTIC NODE print Time, U, V and DOGS AT (6,6) every 100 step
CONVERGENCE REFErence U in LOCAL mode 1.E-6
/
SELEct from (1,1) to (999,999) every (4,2) grid node
OUTPUT U DOGS, MEAT, KILL NARRow mode in SELEcted subdomain
/
HISTORY of DOGS, MEAT, KILL every 16 steps PLOT File="DOGS.HIS"
HISTORY at (2,2) (5,5) (7,7) and (10,10)
/
SOLVE U V P & DOGS for 365 days in step of 0.1
/
END
/

-Back To Top

************************************************************************
TITLE T.044: Testing of Droplet Evaporation Algorithm
************************************************************************
/
GRID NODEs 42 by 37 BY 12
/
COORDINATE X RANGE = 2
COORDINATE Y RANGE = 0.7
COORDINATE Z MINIMUM -0.5 MAX 0.5
/
////////////////////////////////////////////////////////////////////////
* GEOMETRY SPECIFICATION
////////////////////////////////////////////////////////////////////////
/
SYMMetry at Y- boundary
SYMMetry at Y+ boundary
SYMMetry at Z- boundary
SYMMetry at Z+ boundary
/PERIODIC AT Z boundary
INLET X-
OUTLET X+
/
DEFINE GASTEMP = HIGH ! 1500K; Change HIGH to LOW for 373K T_Gas
DEFINE GASPRES = HIGH ! 15Atm; Change HIGH to LOW for 1Atm
/
IF( GASTEMP = LOW ) THEN
SET T = 373. ! Low Temperature S.K. Agarwal et al. AIAA 98-0157
ELSE
SET T = 1500. ! High Temperature S.K. Agarwal et al. AIAA 98-0157
ENDIF
/
IF( GASPRES = LOW ) THEN
GAS PRESSURE = 1.01325E5 ! 1 Atm case
ELSE
GAS PRESSURE = 1.5199E6 ! 15 Atm case
ENDIF
/
SET U = 0.01
SET FU = 0.0
SET O2 = 0.233
SET N2 = 0.767
SET K = 7.5
SET L = 0.014
VISCOSITY 1.8E-5 ! Should change also but smaller influence
/
DENSITY BY GAS LAW
/
////////////////////////////////////////////////////////////////////////
* Initial and Boundary Conditions
////////////////////////////////////////////////////////////////////////
/
GAS FU = 139.26293 CO = 28.01055 H2O= 18.01534
CO2= 44.00995 O2 = 31.99880 N2= 28.01340
/
NOZZLE frequency: 1 2
NOZZLE TRAJECTORY 1 2 50
/NOZZLE TIME STEP = 1.E-5 ITERATIONS 5000 MINIMum T = 200
/
IF( GASTEMP=HIGH ) THEN
NOZZLE DROPLETS Q=0.00001 only 1 SET
DIA X Y Z T U V W
100.0 0.02 0.35 0.0 300. 0.0000001 0.0 0.0
ELSE
NOZZLE DROPLETS Q=0.00001 only 1 SET
DIA X Y Z T U V W
100.0 0.02 0.35 0.0 233. 0.0000001 0.0 0.0
ENDIF
/
SPECific Heat by GORDon-McBride relations on 'CPHTDFLT.ACR'
/
DIAGNOSTIC node (08,08,08) every 1 steps
FLUX DEFault output OFF
SELEct from (1,1,5) to (999,999,7) interval (1,1,1)
OUTPut OFF
/
////////////////////////////////////////////////////////////////////////
* SOLUTION OPTIONS
////////////////////////////////////////////////////////////////////////
/
MATRIX Sweeps in X & Y directions only
CONVERGENCE MAX MASS SOURCE; GLOBAL SUM = 1.0E-8
/
SOLVE 4 STEPS in STEADY MODE
/
END
/

-Back To Top

************************************************************************
TITLE T.045: Testing Locate POLYgon option in 2D
************************************************************************
/
DEBUG GEOMETRY OFF
COOR X RANGE 10.0 1.0
COOR Y RANGE 1.0 1.0
/
DEFINE DIMENSION = 2
***** ==================================================
IF( DIMENSION .EQ. 2 ) THEN
***** ==================================================
/
GRID NODEs 22 by 22
LOCATE POLYGON ID=DIAMOND output elements on 'DIAMOND_LOC.TMP'
(4, 0.5) (5, 0.25), (6, 0.5), (5, 0.75)
DIAGNOSTIC (14,12)
/
***** ==================================================
ELSE this is 3D case
***** ==================================================
/
GRID NODEs 22 by 22 BY 12
LOCATE POLYGON ID=DIAMOND output elements on 'DIAMOND_LOC.TMP'
(4, 0.5, 0.) (5, 0.25, 0.), (6, 0.5, 0.), (5, 0.75, 0.)
(4, 0.5, 1.) (5, 0.25, 1.), (6, 0.5, 1.), (5, 0.75, 1.)
DIAGNOSTIC (14,12,7)
/
***** ==================================================
ENDIF
***** ==================================================
/
BLOCK ID=DIAMOND
/
SET T X- 100.
/
OUTPut T NARROW
/
CONVergence T 1.E-3 5
SOLVE T STEADY 10 10
/
OUTPut OFF
/
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.046: Testing PRINT FLUX and PRINT SHEAR Commands
************************************************************************
GRID NODEs 22 by 22
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMETRY OFF
/
WALL by default
INLET X-
OUTLET X+
/
LOCA (6,6) TO (8,8)
BLOC SELECTED
/
SET at X- U = 1
SET at X- T = 100
SET at X- K = 0.003
/
DENSity 1.0
VISC 1.E-05
/
DIAGnostic output for U V K E at (11,11) every 20 steps
FLUX DEFAult output OFF
/
PRINT FLUX of U and T INLET, OUTLET and WALLS at the end
PRINT FLUX of U, V,W,T at ALL with DETAIL 'T46_FLUX.TMP' every 70 STEPS
PRINT SHEAR at WALLS at the end of simulations
PRINT SHEAR at WALLS DETAIL 'T46_SHEAR.TMP' every 100 STEP
/
OUTPut OFF
/
CONV GLOB 1.E-7
SOLVE STEADY 200
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.047: Testing CORRELATION Commands
************************************************************************
GRID NODEs 22 by 22
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMETRY OFF
/
WALL by default
INLET X-
OUTLET X+
/
LOCA (6,6) TO (8,8)
BLOC SELECTED
/
SET at X- U = 1
SET at X- T = 100
SET at X- K = 0.003
/
DENSity 1.0
VISC 1.E-05
/
DIAGnostic output for U V K E at (11,11) every 20 steps
FLUX DEFAult output OFF
/
SELECT ID=CORRU (10,10) TO (12,12)
CORRELATION OF U 'T47_CORR_U.TMP' ID=CORRU
/
SELECT ID=CORRUV (11,11) TO (14,15)
CORRELATION OF U V 'T47_CORR_UV.TMP' ID=CORRUV
/
LOCATE LIST IJ (2,2) (3,2) (4,2) (5,2) ID=PAIR1
LOCATE LIST IJ (2,3) (3,3) (4,3) (5,3) ID=PAIR2
LOCATE CORRELATED ELEMENTS ID=PAIR1 AND ID=PAIR2 AS ID=PAIR
CORRELATION TWO POINT U V 'T47_CORR_UV2.TMP' for ID=PAIR
/
OUTPut OFF
/
CONV GLOB 1.E-7
SOLVE STEADY 200
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.048: Testing PRINT FORCE Commands
************************************************************************
GRID NODEs 22 by 22
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMETRY OFF
/
WALL by default
INLET X-
OUTLET X+
/
LOCA (6,6) TO (8,8) ID=BODY
BLOC ID=BODY

/
SET at X- U = 1
SET at X- T = 100
SET at X- K = 0.003
/
DENSity 1.0
VISC 1.E-05
/
DIAGnostic output for U V P at (11,11) every 40 steps
/
FLUX DEFAult output OFF
/
PRINT P FORCE for ID=BODY pivot at (0.,0.) reference_P=0. every 10 steps
output on file 'FORCE.TMP'
/
OUTPut OFF
/
CONV GLOB 1.E-7
SOLVE STEADY 200
/
SAVE OFF
/
END
/

-Back To Top

************************************************************************
TITLE T.049: Testing WRITE commands for boundary output
************************************************************************
GRID NODEs 22 by 22
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
DEBUG GEOMETRY OFF
/
WALL by default
INLET X-
OUTLET X+
/
LOCA (6,6) TO (8,8) ID=BODY
BLOC ID=BODY

/
SET at X- U = 1
SET at X- T = 100
SET at X- K = 0.003
/
DENSity 1.0
VISC 1.E-05
/
FLUX DEFAult output OFF
/
DIAGnostic output for U V P at (11,11) every 40 steps
/
OUTPut OFF
/
CONV GLOB 1.E-7
SOLVE STEADY 200
/
SAVE OFF
/
WRITE U P for IO bndries to file 'IO.TMP' ! All walls of the domain
WRITE U P at BOUNdaries file 'EXTERNAL.TMP' ! External domain bndries
WRITE U P for WALLS to file 'WALLS.TMP' ! All walls of the domain
WRITE P X- ID=BODY to file 'BODYX.TMP' ! Specific boundary of subdomain
WRITE P BOUNdaries of ID=BODY file'BODY.TMP' ! All boudnaries of subdomain
/
END
/

-Back To Top

************************************************************************
TITLE T.050: Testing Locate Cylinder command
************************************************************************
GRID NODEs 22 by 22
COOR X Minimum = -0.5 Maximaum = 0.5
COOR Y Minimum = -0.5 Maximaum = 0.5
DEBUG GEOMETRY OFF
/
WALLs by default
/
LOCAte CIRCLE ID=CYLINDER dia=0.2 center (0.,0.) 'CIRCLE.TMP'
/
BLOCk ID=CYLINDER

/
BOUNDARY Y+ U = 1. ! Move top boundary
/
DIAGNOSTIC (20,3)
/
DENSity 1.0
VISC 1.E-02
/
DIAGnostic output for U V P at (20,03) every 40 steps
/
PRINT FORCE ID=CYLINDER pivot (0.,0.) pref=0., every 10 steps
FLUX DEFAult output OFF
/
OUTPut OFF
/
CONV GLOB 1.E-7
SOLVE STEADY 500
/
SAVE OFF
/
END
/
/

-Back To Top

************************************************************************
TITLE T.051: USING BOUNDARY or SOURCE alternative for input of flux
****
**** This problem illustrates the alternative use of either a BOUNDARY
**** or a SOURCE command to specify the input of a flux in to the
**** computational domain.
****
************************************************************************
/
GRID NODEs 22 by 22
COOR X RANGE 10.0 grid expansion ratio of 1.05
COOR Y RANGE 1.0 with uniform grid
/
WALLs by default
INLET at X- boundary
OUTLET at X+ boundary
SYMMETRY at Y- boundary
/
****** The subdomain must consist of elements just inside domain
LOCAte ID=TOP_WALL from (2, 21) TO (21,21) ! must be elements just inside
/
SET U = 1. ! Uniform plug flow to start with
SET T = 0. ! Uniform initial temperature
/
****** Alternative and equivalent use of BOUNdary or SOURce command
****** Only one of the next two options to be active for any run
/
******
//// OPTION A: Explicit flux boundary conditions at top wall
/
/BOUNDARY for T at Y+ of ID=TOP_WALL FLUX = 1 ! per unit area by default
/
******
///// OPTION B: Insulated top wall and source boundary conditions at top wall
/
SOURCE for T at Y+ of ID=TOP_WALL = 1 per unit AREA
BOUNDARY for T at Y+ of ID=TOP_WALL FLUX = 0 ! insulate the boundary
/
******
/
DENSity 1.0
VISC 0.01
PRANDTL number 0.7
SPECIFIC HEAT = 1
/
DIAGnostic output for U P T at (11,11) every 20 steps
/
SELECT (1,1) TO (99,99) freq (4,2)
OUTPut U V P T NARROW for SELEcted region
/
FLUX DEFAult output OFF
CONVergence REFERENCE T Local 1.E-2
/
SOLVE STEADY 1000
/
END
/
/

-Back To Top

************************************************************************
TITLE T.052: Specifying Variable Thermal Conductivity
****
**** This problem illustrates the specification of variable
**** thermal conductivity in a subregion of the computational domain.
****
**** First a temprary variable, KTVARY, is defined as:
****
**** KTVARY = 2. -0.01 * T ; KTVARY = 0, if it is negative
****
**** Then the thermal conductivity in a selected subregion is defined as
****
**** KT = 0.001 +1. * KTVARY
****
**** This ensures that the minimum value of KT= 0.001.
****
**** With the built-in functions, KT can also be defined directly as:
****
**** KT = 2. -0.01 * T ; KT = 0, if it is negative
****
**** But in this case, the minimum value will be zero. The use of a
**** a temporary holding variable for the function allows more
**** flexibility in defining complex or multi-variate relations.
****
************************************************************************
/
GRID NODEs 22 by 22
COOR X RANGE 1.0 with uniform grid
COOR Y RANGE 1.0 with uniform grid
/
SET T = 0. ! Uniform initial temperature
/
BOUNDARY for T at X- VALUE = 100
BOUNDARY for T at X+ VALUE = 0
BOUNDARY for T at Y- FLUX = 1000
BOUNDARY for T at Y+ VALUE = 0
/
/
******
/
DENSity 1.0
SPECIFIC HEAT = 1
/
VISCosity 1.0
/
****
******>> Define conductivity from PRANtl Number everywhere
/
PRANDTL number effective = 0.5 ! Thermal K = viscosity/Prandtl = 2
/
/
******>> Define subregion for variable thermal conductivity
/
LOCATE ID=VARY_KT from (4,4) to (15,15)
/
******>> Define temporary holding variable for the KT function
ALLOCATE KTVARY
SET KTVARY ALWAYS as LINEAr function: 2. -0.01 * T POSItive values
/
/
******>> Define variable conductivity for the subregion
CONDUCTIVITY for T equation in subdomain ID=VARY_KT is a
LINEAr function: 0.001 +1 * KTVARY
/
/
DIAGnostic output for T H KTVARY at (11,11) every 20 steps
/
/
/
SELECT (1,1) TO (99,99) freq (4,2)
OUTPut DIFFUSION coefficient for T in NARROW mode for SELEcted region
/
CONV REFERENCE T Local 1.E-7
SOLVE T STEADY 1000
/
/
SAVE OFF KTVARY T 'T52.SAV'
/
OUTPut T KTVARY NOW in NARROW mode for SELEcted region
/
END
/
/

-Back To Top

1 2 3 4 5 6 7 8 9 10

 

# Back to ANSWER EXAMPLES DATABASE TOC

 

Related Links:

# ANSWER Applications
# ANSWER Examples
# ANSWER Express
# ANSWER Publications
# ANSWER Users
# ANSWER Price List
# Request CFDStudio/ANSWER Demo
# CFDStudio/ANSWER Tutorials
# ANSWER Manual