1
2
3
4
5
6
7
8
9
10
CHAPTER 2: CONDUCTION WITH SPECIFIED CONVECTION
Shortcuts
TITLE Problem 2.01 - CONDUCTION + CONVECTION
TITLE Problem 2.02 - CONVECTION ~ DIFFUSION
TITLE Problem 2.03 - CONVECTION ~ DIFFUSION
TITLE Problem 2.04 - CONVECTION ~ DIFFUSION
TITLE Problem 2.05 - CONVECTION ~ DIFFUSION
TITLE Problem 2.06 - CONVECTION ~ DIFFUSION
TITLE Problem 2.07 - CONVECTION ~ DIFFUSION
TITLE Problem 2.08 - CONVECTION ~ DIFFUSION
TITLE Problem 2.09 - UNIFORM FLOW
TITLE Problem 2.10 - UNIFORM FLOW
TITLE Problem 2.11 - UNIFORM FLOW
TITLE Problem 2.12 - UNIFORM FLOW
TITLE Problem 2.13 - SHEAR FLOW
TITLE Problem 2.14 - SHEAR FLOW
TITLE Problem 2.15 - POISEUILLE FLOW
TITLE Problem 2.16 - POISEUILLE FLOW
************************************************************************
TITLE Problem 2.01 - CONDUCTION + CONVECTION
************************************************************************
/ GRID NODEs 12 by 12 COOR X RANGE 1.0 COOR Y RANGE 2.0
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions SET U CONS 6.0 /IN THE FOLLOWING EXPO, SET A=1./(EXP(1.5U)-1.) ; B=1.5U ; C=0.0 ; D=-A / FOR U=4.0 , A=2.48491E-3 ; B=6.0 ; D=-2.48491E-3 / FOR U=6.0 , A=1.23425E-4; B=9.0; D=-1.23425E-4 SET T EXPO X 1.23425E-4 9.0 0.0 -1.23425E-4 SYMM Y- SYMM Y+
/----------------------------------------------------------------------/
/ Fluid Properties and Constants DENSity 6.0 VISC 8.0 SPEC 4.0 PRAN EFFE 2.0
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30
/----------------------------------------------------------------------/
DIAGnostic output T H at (6,6) every 20 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL CONV T REFE GLOB 1.E-10 1 SOLVE T STEADY 200 / SET FF EXPO X 1.23425E-4 9.0 0.0 -1.23425E-4 SET FF linear function 0. -1. * T ADD to existing value SAVE OFF U V T FF on file '201.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END /
-Back To Top
************************************************************************ TITLE Problem 2.02 - CONVECTION ~ DIFFUSION ************************************************************************ / GRID NODEs 12 by 12 COOR X RANGE 10.0 COOR Y RANGE 10.0 /----------------------------------------------------------------------/ / Initial and Boundary Conditions SET U CONS 5.0 SET T 0.0 / OUTLet at X+ SYMM Y- SYMM Y+ /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 6.0 VISC 8.0 SPEC 0.4 PRAN EFFE 1.E30 /----------------------------------------------------------------------/ / Source and Sink Specifications SOUR T VOLU 30.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS INTE T HYBR RELA T 1.0 LIMI T -1.E30 1.E30 / DIAGnostic output at (6,6) every 20 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF SAVE OFF U V T on file '202.SAV' /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-10 1 SOLVE T STEADY 200 SELEct (1,1) to (999,999) with interval (2,2) OUTPut T for SELEcted subdomain in NARRow mode / END /
-Back To Top
************************************************************************
TITLE Problem 2.03 -
CONVECTION ~ DIFFUSION
************************************************************************
/
GRID NODEs 12 by 12
COOR X RANGE 1.0
COOR Y RANGE 2.0
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
SET U CONS 5.0
/
OUTLET X+
SYMM Y-
SYMM Y+
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 6.0
VISC 8.0
SPEC 0.4
PRAN EFFE 1.E30
/----------------------------------------------------------------------/
/ Source and Sink Specifications
LOCA ID=HEAT (5,1) (7,99)
SOUR T VOLU ID=HEAT 30.0
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS
INTE T HYBR
RELA T 1.0
LIMI T -1.E30 1.E30
/
DIAGnostic output for T and H at (6,6) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SAVE OFF U V T on file '203.SAV'
/------------------------ OPERATIONAL CONTROL
CONV T REFE GLOB 1.E-10 1
SOLVE T STEADY 200
SELEct (1,1) to (999,999) with interval (2,2)
OUTPut T for SELEcted subdomain in NARRow mode
/
END
/
-Back To Top
************************************************************************
TITLE Problem 2.04 -
CONVECTION ~ DIFFUSION
************************************************************************
/
GRID NODEs 12 by 12
COOR X RANGE 10.0
COOR Y RANGE 10.0
GEOMetry ROTAte X and Y by 30 degrees
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
/ SET U=UMAG*COS(THETA) ; V=UMAG*SIN(THETA) : HERE UMAG=5.
SET U CONS 4.330127
SET V CONS 2.5
SET T 0.0
/
OUTLET X+
SYMM Y-
SYMM Y+
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 6.0
VISC 8.0
SPEC 0.4
PRAN EFFE 1.E30
/----------------------------------------------------------------------/
/ Source and Sink Specifications
SOUR T VOLU 30.0
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS
RELA T 1.0
LIMI T -1.E30 1.E30
/
DIAGnostic output for T and H at (6,6) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SAVE OFF U V T on file '204.SAV'
SELEct (1,1) to (999,999) with interval (2,2)
OUTPut T for SELEcted subdomain in NARRow mode
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV T REFE GLOB 1.E-10 1
SOLVE T STEADY 200
/
END
/
-Back To Top
************************************************************************
TITLE Problem 2.05 -
CONVECTION ~ DIFFUSION
************************************************************************
/
GRID NODEs 12 by 12
COORdinate corners (0.,0.) (10.,0.) (5.,8.6602535) (15.,8.6602535)
GEOMetry ROTAte X and Y by 60 degrees
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
/ SET U=UMAG*COS(THETA) ; V=UMAG*SIN(THETA) : HERE UMAG=5.
SET V CONS 4.330127
SET U CONS 2.5
SET T 0.0
/
OUTLET X+
SYMM Y-
SYMM Y+
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 6.0
VISC 8.0
SPEC 0.4
PRAN EFFE 1.E30
/----------------------------------------------------------------------/
/ Source and Sink Specifications
SOUR T VOLU 30.0
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS
RELA T 1.0
LIMI T -1.E30 1.E30
/
DIAGnostic output for T and H at (6,6) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SAVE OFF U V T on file '205.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV T REFE GLOB 1.E-10 1
SOLVE T STEADY 200
SELEct (1,1) to (999,999) with interval (2,2)
OUTPut T for SELEcted subdomain in NARRow mode
/
END
/
-Back To Top
************************************************************************
TITLE Problem 2.06 -
CONVECTION ~ DIFFUSION
************************************************************************
/
GRID NODEs 12 by 12
COORdinate corners (0.,0.) (10.,0.) (5.,8.) (20.,10.)
GEOMetry ROTAte X and Y by 30 degrees
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
SET U CONS 5.0
SET V CONS 2.0
SET T LINE X 0.0 3.0
SET T LINE Y ADD 0.0 4.0
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 6.0
VISC 8.0
SPEC 7.0
PRAN EFFE 1.E30
/----------------------------------------------------------------------/
/ Source and Sink Specifications
SOUR T VOLU 996.0
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS
RELA T 1.0
LIMI T -1.E30 1.E30
/
DIAGnostic output for T and H at (6,6) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV T REFE GLOB 1.E-10 1
SOLVE T STEADY 200
/
SET FF LINEar function: 0. 3 * X
SET FF LINEar function: 0. 4 * Y ADD to existing
SET FF LINEar function: 0. -1 * T ADD to existing
SAVE OFF U V T FF on file '206.SAV'
SELEct (1,1) to (999,999) with interval (2,2)
OUTPut FF for SELEcted subdomain in NARRow mode
/
END
/
-Back To Top
************************************************************************
TITLE Problem 2.07 -
CONVECTION ~ DIFFUSION
************************************************************************
/
GRID NODEs 12 by 12
COORdinate corners (0.,0.) (10.,0.) (5.,8.) (20.,10.)
GEOMetry ROTAte X and Y by 30 degrees
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
SET U CONS 5.0
SET V CONS 2.0
SET T POLYnomial in X: coefficients: 0 0 3 0 0
SET T POLYnomial in Y: coefficients: 0 0 4 0 0 ADD to existing
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 6.0
VISC 8.0
SPEC 7.0
PRAN EFFE 1.E30
/
SOUR T LINE X 0.0 1260. VOLU
SOUR T LINE Y 0.0 672. VOLU
/
RELA T 1.0
LIMI T -1.E30 1.E30
/
DIAGnostic output for T and H at (6,6) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
/
CONV T REFE GLOB 1.E-10 1
SOLVE T STEADY 200
/
SET FF POLYnomial in X: coefficients: 0 0 3 0 0
SET FF POLYnomial in Y: coefficients: 0 0 4 0 0 ADD to existing
SET FF LINEar function 0. -1. * T ADD to existing
SAVE OFF U V T FF on file '207.SAV'
SELEct (1,1) to (999,999) with interval (2,2)
OUTPut FF for SELEcted subdomain in NARRow mode
/
END
/
-Back To Top
************************************************************************
TITLE Problem 2.08 -
CONVECTION ~ DIFFUSION
************************************************************************
/
GRID NODEs 12 by 5 by 12
COORdinate corners (0,0,0) (10,0,0) (0,1,0) (10,1,0)
(5,0,8) (20,0,10) (5,1,8) (20,1,10)
GEOMetry ROTAte X and Z by 30 degrees
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
/
SET U CONS 5.0
SET V CONS 0.0
SET W CONS 2.0
SET T POLYnomial in X: coefficients: 0 0 3 0 0
SET T POLYnomial in Z: coefficients: 0 0 4 0 0 ADD to existing
/
SYMM Y-
SYMM Y+
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 6.0
VISC 8.0
SPEC 7.0
PRAN EFFE 1.E30
/----------------------------------------------------------------------/
/ Source and Sink Specifications
SOUR T LINE X 0.0 1260. VOLU
SOUR T LINE Z 0.0 672. VOLU
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS
/
INTE T HYBR
RELA T 1.0
LIMI T -1.E30 1.E30
/
DIAGnostic output for T and H at (6,3,6) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEct (1,1,1) to (999,999,999) with interval (2,1,4)
OUTPut FF for SELEcted subdomain in NARRow mode by XZ planes
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV T REFE GLOB 1.E-10 1
SOLVE T STEADY 200
/
SET FF POLYnomial in X: coefficients: 0 0 3 0 0
SET FF POLYnomial in Z: coefficients: 0 0 4 0 0 ADD to existing
SET FF LINEar function: 0 -1 * T ADD to existing
SAVE OFF U V W T FF on file '208.SAV'
/
END
/
-Back To Top
************************************************************************
TITLE Problem 2.09 -
UNIFORM FLOW
************************************************************************
/
GRID NODEs 12 by 12 by 12
COOR X RANGE = 1.0
COOR Y RANGE = 2.0
COOR Z RANGE = 3.0
GEOMetry ROTAte X and Y by 30 degrees
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
SET U CONS 4.0
SET V CONS 5.0
SET W CONS 6.0
/
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 7.0
VISC 0.0
/----------------------------------------------------------------------/
DIAGnostic output for U V W at (5,5,5) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SAVE OFF U V W P on file '209.SAV'
SELEct (1,1,5) to (999,999,999) with interval (2,4,4)
OUTPut U,V,W,P for SELEcted subdomain in NARRow mode by XZ planes
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-6
SOLVE FOR 1000 STEPS IN STEADY MODE
/
END
/
-Back To Top
************************************************************************
TITLE Problem 2.10 -
UNIFORM FLOW
************************************************************************
/
GRID NODEs 12 by 12 by 12
COORdinate corners (0,0,0) (10,0,0) (0,5,0) (10,7,0)
(5,0,8) (20,0,10) (5,8,8) (20,10,10)
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
SET U CONS 4.0
SET V CONS 5.0
SET W CONS 6.0
/
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 7.0
VISC 0.0
/----------------------------------------------------------------------/
DIAGnostic output U V W at (5,5,5) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEct (1,1,5) to (999,999,999) with interval (2,4,4)
OUTPut U,V,W,P for SELEcted subdomain in NARRow mode by XZ planes
SAVE OFF U V W P on file '210.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-6
SOLVE FOR 1000 STEPS IN STEADY MODE
/
END
/
-Back To Top
************************************************************************
TITLE Problem 2.11 -
UNIFORM FLOW
************************************************************************
/
GRID NODEs 12 by 12
GEOM Circular arc: r=0.5, 180 degress
GEOMetry ROTAte X and Y by 30 degrees
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
SET U CONS 4.0
SET V CONS 5.0
/
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 7.0
VISC 0.0
/----------------------------------------------------------------------/
/----------------------------------------------------------------------/
DIAGnostic output U V at (5,5) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEct (1,1) to (999,999) with interval (2,2)
OUTPut U,V,P for SELEcted subdomain in NARRow mode
SAVE OFF U V P on file '211.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-6
SOLVE FOR 1000 STEPS IN STEADY MODE
/
END
/
-Back To Top
************************************************************************
TITLE Problem 2.12 -
UNIFORM FLOW
************************************************************************
/
GRID NODEs 12 by 12
GEOMetry is ANNUlus with outer rad = 2, inner = 1 offset = 0.5
GEOMetry ROTAte X and Y by 30 degrees
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
SET U CONS 4.0
SET V CONS 5.0
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 7.0
VISC 0.0
/----------------------------------------------------------------------/
DIAGnostic output U V and P at (5,5) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEct (1,1) to (999,999) with interval (2,2)
OUTPut U,V,P for SELEcted subdomain in NARRow mode
SAVE OFF U V P on file '212.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-6
SOLVE FOR 1000 STEPS IN STEADY MODE
/
END
/
-Back To Top
************************************************************************
TITLE Problem 2.13 -
SHEAR FLOW
************************************************************************
/
GRID NODEs 12 by 12
COORdinate corners (0.,0.) (10.,0.) (5.,8.) (20.,10.)
GEOMetry ROTAte X and Y by 30 degrees
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
SET U LINE Y 0.0 4.0
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 7.0
VISC 0.0
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS
RELA P 0.2
/----------------------------------------------------------------------/
DIAGnostic output of U V P RP at (5,5) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEct (1,1) to (999,999) with interval (2,2)
OUTPut U,V,P for SELEcted subdomain in NARRow mode
SAVE OFF U V P on file '213.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-6
SOLVE FOR 1000 STEPS IN STEADY MODE
/
END
/
-Back To Top
************************************************************************
TITLE Problem 2.14 -
SHEAR FLOW
************************************************************************
/
GRID NODEs 12 by 5 by 12
COORdinate corners (0,0,0) (10,0,0) (0,1,0) (10,1,0)
(5,0,8) (20,0,10) (5,1,8) (20,1,10)
GEOMetry ROTAte X and Z by 30 degrees
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
SET U LINE Z 0.0 4.0
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 7.0
VISC 0.0
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS
RELA P 0.2
/----------------------------------------------------------------------/
DIAGnostic output U V W P and RP at (5,3,5) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEct (1,1,1) to (999,999,999) with interval (2,4,2)
OUTPut U and P for SELEcted subdomain in NARRow mode by XZ planes
SAVE OFF U V W P on file '214.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV GLOB 1.E-6
SOLVE FOR 1000 STEPS IN STEADY MODE
/
END
/
-Back To Top
************************************************************************
TITLE Problem 2.15 -
POISEUILLE FLOW
************************************************************************
/
GRID NODEs 42 by 42
COORdinate corners (0.,0.) (10.,0.) (5.,8.) (20.,10.)
GEOMetry ROTAte X and Y by 30 degrees
INLET X-
INLET X+
INLET Y-
INLET Y+
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
SET U POLYnomial in Y: coefficients: 0 0 3 0 0
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 7.0
VISC 5.0
/----------------------------------------------------------------------/
/ Source and Sink Specifications
SOUR U VOLU -30.0
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS
RELA P 0.2
DIFF SECO SKEW U V W
/----------------------------------------------------------------------/
DIAGnostic output U V P RU RP at (5,5) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEct (1,1) to (999,999) with interval (4,4)
OUTPut for SELEcted subdomain in NARRow mode
SAVE OFF U V P on file '215.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV REFErence for U tolerance 1.E-3
SOLVE FOR 1000 STEPS IN STEADY MODE
/
END
/
-Back To Top
************************************************************************
TITLE Problem 2.16 -
POISEUILLE FLOW
************************************************************************
/
GRID NODEs 32 by 5 by 32
COORdinate corners (0,0,0) (10,0,0) (0,1,0) (10,1,0)
(5,0,8) (20,0,10) (5,1,8) (20,1,10)
GEOMetry ROTAte X and Z by 30 degrees
INLET X-
INLET X+
INLET Y-
INLET Y+
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
/
SET U POLYnomial in Z: coefficients: 0 0 3 0 0
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 7.0
VISC 5.0
/----------------------------------------------------------------------/
/ Source and Sink Specifications
SOUR U VOLU -30.0
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS
RELA P 0.2
DIFF SECO SKEW U V W
/----------------------------------------------------------------------/
DIAGnostic output U V W P RW at (5,3,5) every 20 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEct (1,1,1) to (999,3,999) with interval (6,2,3)
OUTPut for SELEcted subdomain in NARRow mode by XZ planes
SAVE OFF U W P on file '216.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
CONV REFErence for U tolerance 1.E-3
SOLVE FOR 1000 STEPS IN STEADY MODE
/
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
|