CHAPTER 3: FORCED CONVECTION FLOW

Shortcuts

Problem 3.01 - SOLID BODY ROTATION : Z-Y COORDINATES

Problem 3.02 - AXISYMMETRIC SWIRLING FLOW

Problem 3.03 - LAMINAR FULLY DEVELOPED COUETTE FLOW

Problem 3.04 - TURBULENT FULLY DEVELOPED COUETTE FLOW

Problem 3.05 - LAMINAR DEVELOPING FLOW IN A CHANNEL (Same as 3.06 without blocked solid slab at top)

Problem 3.06 - LAMINAR DEVELOPING FLOW IN A CHANNEL (Same as 3.05 with blocked solid slab at bottom)

Problem 3.07 - LAMINAR FULLY DEVELOPED FLOW IN A CHANNEL (Same as 3.08 without blocked slab at top)

Problem 3.08 - LAMINAR FULLY DEVELOPED FLOW IN A CHANNEL (Same as 3.07 with blocked slab at bottom)

Problem 3.09 - FULLY DEVELOPED TURBULENT FLOW IN CHANNEL (Same as 3.10 without blocked slab)

Problem 3.10 - FULLY DEVELOPED TURBULENT FLOW IN CHANNEL (Same as 3.09 with blocked slab)

Problem 3.11 - DEVELOPING TURBULENT FLOW IN A PIPE

Problem 3.12 - LAMINAR DEVELOPING FLOW IN A CHANNEL

Problem 3.13 - TURBULENT DEVELOPING FLOW IN A CHANNEL

Problem 3.14 - TEST OF EQUIVALENCE OF INLET & INJECTION

Problem 3.15 - PARALLELOPIPED CAVITY WITH MOVING WALL - RE = 100

Problem 3.16 - FLOW BEHIND A STEP IN A FLAT CHANNEL - H/h=2, Re=50

Problem 3.17 - FORWARD FACING STEP IN A FLAT CHANNEL - H/h=2, Re=50

Problem 3.18 - DOUBLET WITH HEAT SOURCE IN AN OPEN DOMAIN

Problem 3.19 - 3D ROOM VENTILATION PROBLEM

-Back To Top

************************************************************************
TITLE Problem 3.02 - AXISYMMETRIC SWIRLING FLOW
***** Cartesian mode: few cells in the theta direction over
***** a small angle and periodic boundary condition.
************************************************************************
/
GRID NODEs 12 by 12 by 7
COOR CORN (0.,0.,-0.05) (1.,0.,-0.05) (0.,1.,-0.05) (1.,1.,-0.05)
(0.,0.,0.05) (1.,0.,0.05) (0.,1.,0.05) (1.,1.,0.05) CYLI
GEOM CART
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
/
SET W LINE Y 0.0 4.0
SET V LINE Z 0.0 -4.0
SET FF POLYnomial in Z: coefficients: 0 0 80 0 0
SET FF POLYnomial in Y: coefficients: 0 0 80 0 0 ADD to existing
SET FH LINE FF 0.0 1.0
SET FH LINE P ADD 0.0 -1.0
/
PERI Z THETA 5.729578
/INLET Y-
/INLET Y+
SYMM X-
SYMM X+
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 10.0
VISC 1.0
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS
/DIFF SECO SKEW U V W P
/----------------------------------------------------------------------/
DIAGnostic output U V W P RW at (2,5,5) every 100 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEt (1,1,1) to (999,999,999) with interval of (2,2,2)
OUTPut for SELEcted
SAVE OFF U V W P FF FH on file '302.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-10
SOLVE for U V W P max 1000 steps min 1000 in STEAdy mode
/
END
/
/

-Back To Top

************************************************************************
TITLE Problem 3.03 - LAMINAR FULLY DEVELOPED COUETTE FLOW
***
*** Couette flow: set top wall moving & source of u to zero
*** Poiseuelle flow: set top wall vel=0 and give finite u source
*** For this case must get ucl=(source of u)*(h**2)/(8*visc)
*** where ucl=center line velocity and h is plate spacing.
************************************************************************
/
GRID NODEs 3 by 42
COORdinate corners (0.,0.) (1.,0.) (0,1.) (1.,1.)
PERIodic in X direction
WALL at Y-
WALL at Y+
/
/ Initial and Boundary Conditions
BOUN U VALU CONS Y+ 100.0
BOUN V VALU CONS Y+ 0.0
/
/ Fluid Properties and Constants
DENSity 1.0
VISC 1.0
SPEC 1.0
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS
/----------------------------------------------------------------------/
DIAGnostic output U at (2,21) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEt (1,1) to (999,999) with interval of (1,2)
OUTPut U for SELEcted
SAVE OFF U V P on file '303.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
RELA U=0.9, V=0.9
CONV U REFE 1.E-3
SOLVE U V P STEADY 1000
/
END
/
/

-Back To Top

************************************************************************
TITLE Problem 3.04 - TURBULENT FULLY DEVELOPED COUETTE FLOW
*****
***** Couette flow: set top wall moving & set source of u to zero
***** Poiseuelle flow: set top wall vel=0 and give finite u source
************************************************************************
/
GRID NODEs 3 by 42
COORdinate corners (0.,0.) (1.,0.) (0,1.) (1.,1.)
PERIodic in X direction
WALL at Y-
WALL at Y+
/
BOUN U VALU CONS Y+ 100.0
/
DENSity 1.0
VISC .001
/
RELA U=0.9, V=0.9, K=0.8, E=0.8
/
DIAGnostic output U P at (2,21) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEt (1,1) to (999,999) with interval of (1,2)
OUTPut U for SELEcted
SAVE OFF U V P K E L MU on file '304.SAV'
/
CONV U REFE 1.E-6
SOLVE U,V,P STEADY 1000
/
END
/

-Back To Top

************************************************************************
TITLE Problem 3.05 - LAMINAR DEVELOPING FLOW IN A CHANNEL
***** Same as 3.06 without blocked solid slab at top
************************************************************************
/
GRID NODEs 42 by 22
COOR X RANGE 10.0 1.05
COOR Y RANGE 1.0 1.0
//SELEct (1,22) (42,42)
//BLOCKAGE at SELEcted zone
/ Initial and Boundary Conditions
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
DIFF SECO SKEW U V P T
/----------------------------------------------------------------------/
DIAGnostic output U V P T RP at (6,21) every 50 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEt (1,1) to (999,999) with interval of (4,2)
OUTPut U T for selected domain in NARRow mode
SAVE OFF U V P T on file '305.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-8
SOLVE U V P T STEADY 1000
/
END
/
/

-Back To Top

************************************************************************
TITLE Problem 3.06 - LAMINAR DEVELOPING FLOW IN A CHANNEL
***** Same as 3.05 with blocked solid slab at bottom
************************************************************************
/
GRID NODEs 42 by 42
COOR X RANGE 10.0 1.05
COOR Y RANGE 2.0 1.0
/
WALL at undefined outer boundaries
LOCA ID=INLT (1,1) (1,21)
INLET ID=INLT X-
LOCA ID=OUTLET (99,1) (99,21)
OUTL X+ ID=OUTLET
SYMM Y-
SELEct (1,22) (42,42)
BLOCKAGE at SELEcted zone
/ Initial and Boundary Conditions
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
DIFF SECO SKEW U V P T
/----------------------------------------------------------------------/
DIAGnostic output U V P T RP at (6,21) every 50 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEt (1,1) to (999,999) with interval of (4,4)
OUTPut U T for selected domain in NARRow mode
SAVE OFF U V P T on file '306.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-8
SOLVE U V P T STEADY 1000
/
END
/
/

-Back To Top

************************************************************************
TITLE Problem 3.07 - LAMINAR FULLY DEVELOPED FLOW IN A CHANNEL
***** Same as 3.08 without blocked slab at top
************************************************************************
/
//// Case A: No blockage
GRID NODEs 3 by 22
COOR Y RANGE 1.0
WALL at undefined outer boundaries
//// Case B: Block top half
/GRID NODEs 3 by 42
/COOR Y RANGE 2.0
/SELEct (1,22) (999,999) <<<< Activate for Blockage
/BLOCKAGE at SELEcted zone
/
COOR X RANGE 1.0
/
SYMM Y-
PERIodic in X direction
/ Fluid Properties and Constants
DENSity 1.0
VISC 1.0
/ Source and Sink Specifications
LOCA ID=REG1 (1,1) (3,21)
SOUR U CONS VOLU 200.0 for ID=REG1
/----------------------------------------------------------------------/
DIAGnostic output U P at (2,21) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SAVE OFF U V P on file '307.SAV'
SELEt (1,1) to (999,999) with interval of (1,2)
OUTPut U for SELEcted subdomain
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
RELAX U=1.0, V=1.0
CONV U REFE 1.E-3
SOLVE U V P STEADY 1000
/
END
/

-Back To Top

************************************************************************
TITLE Problem 3.08 - LAMINAR FULLY DEVELOPED FLOW IN A CHANNEL
***** Same as 3.07 with blocked slab at bottom
************************************************************************
/
//// Case A: No blockage
/GRID NODEs 3 by 22
/COOR Y RANGE 1.0
WALL at undefined outer boundaries
//// Case B: Block top half
GRID NODEs 3 by 42
COOR Y RANGE 2.0
SELEct (1,22) (999,999) <<<< Activate for Blockage
BLOCKAGE at SELEcted zone
/
COOR X RANGE 1.0
/
SYMM Y-
PERIodic in X direction
/ Fluid Properties and Constants
DENSity 1.0
VISC 1.0
/ Source and Sink Specifications
LOCA ID=REG1 (1,1) (3,21)
SOUR U CONS VOLU 200.0 for ID=REG1
/----------------------------------------------------------------------/
DIAGnostic output U P at (2,21) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SAVE OFF U V P on file '308.SAV'
SELEt (1,1) to (999,999) with interval of (1,2)
OUTPut U for SELEcted subdomain
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
RELAX U=1.0, V=1.0
CONV U REFE 1.E-3
SOLVE U V P STEADY 1000
/
END
/
/

-Back To Top

************************************************************************
TITLE Problem 3.09 - FULLY DEVELOPED TURBULENT FLOW IN CHANNEL
***** Same as 3.10 without blocked slab
************************************************************************
/
GRID NODEs 3 by 22
COOR Y RANGE 1.0
WALL at undefined outer boundaries
/
/GRID NODEs 3 by 42
/COOR Y RANGE 2.0
/SELEct (1,22) (3,42)
/BLOCKAGE at SELEcted zone
/
COOR X RANGE 1.0
/ Initial and Boundary Conditions
SET U=1.0E3
SET K=1.0E3
SET E=1.0E6
PERIodic in X direction
SYMM Y-
/ Fluid Properties and Constants
DENSity 1.0
VISC .01
/ Source and Sink Specifications
LOCA ID=REG1 (1,2) (3,21)
SOUR U CONS VOLU 100.0 for ID=REG1
/ SOLUTION OPTIONS
RELA U=1.0, V=1.0, K=0.8, E=0.8
/
DIAGnostic output for U P K E L at (2,11) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEt (1,1) to (999,999) with interval of (1,2)
OUTPut K L for SELEcted subdomain NARRow mode
SAVE OFF U V P K E L MU on file '309.SAV'
/ OPERATIONAL CONTROL
CONV U REFE 1.E-3
SOLVE STEADY 1000
/
END
/
/

-Back To Top

************************************************************************
TITLE Problem 3.10 - FULLY DEVELOPED TURBULENT FLOW IN CHANNEL
***** Same as 3.09 with blocked slab
************************************************************************
/
/GRID NODEs 3 by 22
/COOR Y RANGE 1.0
/
GRID NODEs 3 by 42
COOR Y RANGE 2.0
PERIodic in X direction
SYMM at Y-
WALL at Y+
SELEct (1,22) (3,42)
BLOCKAGE at SELEcted zone
/
COOR X RANGE 1.0
/ Initial and Boundary Conditions
SET U=1.0E3
SET K=1.0E3
SET E=1.0E6
/ Fluid Properties and Constants
DENSity 1.0
VISC .01
/ Source and Sink Specifications
LOCA ID=REG1 (1,2) (3,21)
SOUR U CONS VOLU 100.0 for ID=REG1
/ SOLUTION OPTIONS
RELA U=1.0, V=1.0, K=0.8, E=0.8
/
DIAGnostic output for U P K E L at (2,11) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEt (1,1) to (999,999) with interval of (1,2)
OUTPut K L for SELEcted subdomain NARRow mode
SAVE OFF U V P K E L MU on file '310.SAV'
/ OPERATIONAL CONTROL
CONV U REFE 1.E-3
SOLVE STEADY 1000
/
END
/

-Back To Top

************************************************************************
TITLE Problem 3.11 - DEVELOPING TURBULENT FLOW IN A PIPE
************************************************************************
/
GRID NODEs 22 by 22
COOR X RANG 1.0 1.05
COOR R RANG 0.1
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
INLE at X-
OUTL at X+
SYMM at Y-
WALL at Y+
SET U 1.0
SET K .003
SET L .025
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 1.0
VISC 1.E-5
/
DIAGnostic output for U V P K L at (2,21) every 100 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEt (1,1) to (999,999) with interval of (4,2)
OUTPut U V K L in SELEcted subdomain NARRow mode
SAVE OFF U V P K E L MU on file '312.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-10 1
SOLVE STEADY 1000
/
END
/
/

-Back To Top

************************************************************************
TITLE Problem 3.12 - LAMINAR DEVELOPING FLOW IN A CHANNEL
**** Testing internal mass injection capability
**** Case A --> inlet at X- boundary
**** Case B --> mass injection just inside the X- boundary
************************************************************************
/
GRID NODEs 22 by 22
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,22) !! Case A
/SET U CONS ID=INLT 100.0 !! Case A
/SET T CONS ID=INLT 1.0 !! Case A
/INLET ID=INLT X- !! Case A
/
SOUR AREA X- FLOW ID=INLT 200 U=100 V=0 T=1.0 !! Case B
BOUN U VALU X- 100.0 !! Case B
BOUN T VALU X- 1.0 !! Case B
/
OUTL X+
SYMM Y+
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 2.0
VISC 1.0
SPEC 1.0
PRAN EFFE 0.7
/
DIAGnostic output for U V P T at (11,11) every 50 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SAVE OFF U V P T on file '313.SAV'
SELEct (1,1) to (999,999) with frequency (4,2)
OUTPut U V P for SELEcted subdomain
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-10 1
SOLVE U V P T STEADY 1000
/
END
/

-Back To Top

************************************************************************
TITLE Problem 3.13 - TURBULENT DEVELOPING FLOW IN A CHANNEL
**** Testing internal mass injection capability
**** Case A --> inlet at X- boundary
**** Case B --> mass injection just inside the X- boundary
************************************************************************
/
GRID NODEs 22 by 22
COOR X RANGE 2.0 1.05
COOR Y RANGE 0.1
WALL at undefined outer boundaries
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
/
/ InitialISE U, K, L AT ALL NODES EXCEPT THOSE AT INLET (1,2)-(1,21)
/
LOCA ID=PAT1 (2,1) (22,22)
SET U CONS ID=PAT1 1.0
SET K CONS ID=PAT1 .003
SET L CONS ID=PAT1 .0025
LOCA ID=PAT2 (1,1) (1,1)
SET U CONS ID=PAT2 1.0
SET K CONS ID=PAT2 .003
SET L CONS ID=PAT2 .0025
LOCA ID=PAT3 (1,22) (1,22)
SET U CONS ID=PAT3 1.0
SET K CONS ID=PAT3 .003
SET L CONS ID=PAT3 .0025
/
/ CASE A ---> USE FOLLOWING FIVE LINES TO GIVE INLET BC
/
LOCA ID=INLT (1,2) (1,21)
SET U CONS ID=INLT 1.0
SET K CONS ID=INLT .003
SET L CONS ID=INLT .0025
INLET ID=INLT X-
/
/ CASE B --> USE FOLLOWING TWO LINES
/
/SYMM X-
/LOCA ID=MFLO (2,2) (2,21)
/SOUR FLOW ID=MFLO .005 U=1 V=0 K=.003 L=.0025
/
OUTL X+
SYMM Y+
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 1.0
VISC 1.E-5
/
DIAGnostic U V P K E output at (11,11) every 100 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SAVE OFF U V P K E L MU on file '314.SAV'
SELEt (1,1) to (999,999) with interval of (2,2)
OUTPut U V P K E in SELEcted subdomain
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-10 1
SOLVE STEADY 1000
/
END
/

-Back To Top

************************************************************************
TITLE Problem 3.14: TEST OF EQUIVALENCE OF INLET & INJECTION
************************************************************************
/
GRID NODEs 22 by 12
COOR X RANGE 1.0
COOR Y RANGE 0.1
/
/INLET X- !! active for inlet
LOCATE 2,2 TO 2,999 !! active for injection
SOURce: X- FLOW=.01 U=1 K=.003 L=.0025 SELEcted !! active for injection
/
WALL at undefined outer boundaries
OUTL X+
SYMM Y+
/
SET U CONS 1.0
SET K CONS .003
SET L CONS .0025
/
DENS 1.0
VISC 1.E-5
/
DIAG U V P K L (11,11) every 20 steps
DEBUG GEOMERTY OFF
SELEt (1,1) to (999,999) with interval of (2,2)
OUTPut U V P K L in SELEcted subdomain
/
CONV GLOB 1.E-8
SOLVE STEADY 200 minumum 150
SAVE OFF '315.SAV'
/
END
/

-Back To Top

***********************************************************************
TITLE Problem 3.15: PARALLELOPIPED CAVITY WITH MOVING WALL - RE = 100
***********************************************************************
/
GRID NODEs BY 22 BY 22
/
COORDINATE CORNers (0,0), (1,0) (1,1) (2,1)
/
WALL at undefined outer boundaries
BOUNdary U: at boundary Y+, VALUE=1.
/
VISCOSITY 0.01
/
DIAGNOSTIC NODE U V P RP AT (6,6) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
CONVERGENCE REFERENCE for U in LOCAL mode 1.E-6
SOLVE FOR 1000 STEPS IN STEADY MODE
SAVE OFF U,V,P, ON '316.SAV'
/
SELEt (1,1) to (999,999) with interval of (2,2)
OUTPut in SELEcted subdomain
/
END
/

-Back To Top

***********************************************************************
TITLE Problem 3.16: FLOW BEHIND A STEP IN A FLAT CHANNEL - H/h=2, Re=50
***********************************************************************
/
GRID NODEs 20 by 12
/
COORDINATE X
0.00, 0.25, 0.75, 1.25, 1.75, 2.25, 2.75, 3.25, 4.00, 5.00
6.00 7.00, 8.00, 9.00, 11.0, 13.0, 15.0, 17.0, 19.0, 20.0
/
WALL at undefined outer boundaries
SELEct from ( 1,2) to ( 1,6)
INLEt X- at SELEcted boundary segment
SELEct from (20,2) to (20,11)
OUTLet X+ at SELEcted boundary segment
SELEct from (1,7) TO (5,11)
BLOCKAGE at SELEcted zone
/
LOCAte from (1,1) to (1,99)
SET U = 1.0 at X- boundary of SELEcted subdomain
/
LAMINAR FLOW
VISCOSITY VALUE =2.E-2
/
DIAGNOSTIC NODE U V P RP (7,9) print every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
/
RELAX U = 0.25, V=0.25
MATRIX P=10 SWEEPS Y direction only
CONVERGENCE by GLOBAL mode to epsilon = 1.E-6
/
SOLVE for 1000 STEPS in STEADY MODE
/
SELect window from (1,1) to (20,12) skip 2, 1
OUTPut U, P for SELEcted window
SAVE OFF U,V,P ON '317.SAV'
/
END
/

-Back To Top

************************************************************************
TITLE Problem 3.17: FORWARD FACING STEP IN A FLAT CHANNEL - H/h=2, Re=50
************************************************************************
/
GRID NODEs 20 by 12
/
COORDINATE X RANGE=20, RATIO=1.
/
WALL at undefined outer boundaries
INLEt at X- boundary
SELEct from (20,7) to (20,12)
OUTLet X+ at SELEcted boundary segment
SELEct from (11,1) TO (20,6)
BLOCKAGE for SELEcted zone
/
LOCAte from (1,1) to (1,99)
SET U = 1.0 at X- boundary of SELEcted subdomain
/
LAMINAR FLOW
VISCOSITY 2.E-2
/
DIAGNOSTIC NODE U V P RP (7,3) print every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
/
RELAx U=0.35 V=0.35
MATRIX SWEEPS in Y direction only
CONVERGENCE by GLOBAL mode to epsilon = 1.E-6
/
SOLVE for 250 STEPS in STEADY MODE
/
SELect window from (1,1) to (20,12) skip 2, 1
OUTPut U, P for SELEcted window
SAVE OFF U,V,P ON '318.SAV'
/
END
/

-Back To Top

************************************************************************
TITLE Problem 3.18: DOUBLET WITH HEAT SOURCE IN AN OPEN DOMAIN
************************************************************************
/ Fluid with finite viscosity.
/ Note that there is an inherent asymmetry since the incoming fluid
/ has zero momentum but the outgoing fluid carries finite momentum
/
GRID NODEs 22 by 21
/
COORDINATE X RANGE 20
COORDINATE Y RANGE 10
/
SYMMetry at Y+ boundary
SYMMetry at Y- boundary
/
SELEct COORdinate window from ( 4.0,4.75) to (5.0, 5.25)
SOURCE for FLOW: Q = 0.001; T = 1.0 at SELEcted location
SELEct COORdinate window from (15.0,4.75) to (16.0, 5.25)
SOURCE for FLOW: Q = -0.001 at SELEcted location
/
VISCOSITY VALUE = 5.E-6
/
DIAGNOSTIC NODE U V P T RP (11,11) print every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
/
CONVERGENCE REFERENCE U 1.E-11 for in LOCAL mode
SOLVE for 500 STEPS in STEADY MODE
/
SAVE OFF U, V, P ON '319.SAV'
SELEt (1,1) to (999,999) with interval of (2,1)
OUTPut in SELEcted subdomain
/
END
/

-Back To Top

***********************************************************************
TITLE Problem 3.19:  3D ROOM VENTILATION PROBLEM
***********************************************************************
/
GRID NODEs 18 by 12 by 8
COORdinate X RANGe 8 m
COORdinate Y RANGe 5 m
COORdinate Z RANGe 3 m
/
SET U = 0
SET V = 0
SET W = 0
SET K = 10
SET L = 0.1
/
WALL at undefined outer boundaries
LOCAte ( 4,7,99)
INLET at Z+ at SELECted location
SET W = -2 at Z+ boundary of SELECted subdomain
LOCAte (15,7,99)
INLEt AT Z+ at SELECted location
SET W = -4 at Z+ boundary of SELECted subdomain
LOCAte ( 9,7,99)
OUTLet at Z+ at SELEcted location
/
DIAGnostics U V W P K E at 9,7,7 every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEct (1,3,2) to (99,99,99) interval of (2,4,2)
OUTPut U, V, W, K L XZ planes in SELECted zone
/
SAVE OFF DATA ON FILE '320.SAV'
/
SOLVe to STEAdy state 500 with a minimum 300 iterations
/
END
/
/

-Back To Top

1 2 3 4 5 6 7 8 9 10

# Back to ANSWER EXAMPLES DATABASE TOC