//////////////////////////////////////////////////////////////////////// / -- ANSWER EXAMPLES BY ACRi -- / /*********************** IMPORTANT -- PLEASE READ / / THIS FILE IS A COLLECTION OF ANSWER EXAMPLES IN ONE SINGLE FILE. / / ITS USES ARE TWO FOLD: / / 1. IF YOU NEED HELP WITH SETTING A TEST CASE WITH SAY NATURAL / CONVECTION, YOU CAN GO TO CHAPTER 4 (see Table of Contents) / AND FIND A PROBLEM SIMILAR TO YOURS, COPY IT FROM "TITLE" TO / "END" COMMANDS, AND PASTE INTO A NEW FILE. IF USING THE / GRAPHICAL USER INTERFACE (GUI), GOTO "FILE > IMPORT LEGACY FILE" TO / IMPORT THIS NEW FILE INTO THE GUI ENVIRONMENT. / / 2. JUST AS YOU SEE HERE, YOU CAN PUT MORE THAN ONE CASE / (CASE = From one "TITLE" to the first "END" following it) / IN THE SAME FILE. NOTE THAT YOU NEED TO RUN THIS KIND OF / "MULTI-CASE" FILES IN COMMAND MODE. GUI ENVIRONMENT IS DESIGNED / TO WORK WITH "SINGLE-CASE" FILES. / / TO RUN ANSWER IN COMMAND (DOS) MODE, ALL YOU NEED TO DO IS ENTER THE / FOLLOWING SYNTAX AT A DOS SHELL (COMMAND PROMPT): / / answer / / NOTE THAT IF YOU DO NOT WANT TO ENTER THE WHOLE PATH TO ANSWER.EXE, / YOU WILL NEED TO ADD C:\ACRI\BIN TO YOUR SYSTEM PATH. / //////////////////////////////////////////////////////////////////////// / TABLE OF CONTENTS: ***** ***** CHAPTER 1: PURE CONDUCTION OR DIFFUSION ***** CHAPTER 2: CONDUCTION WITH SPECIFIED CONVECTION ***** CHAPTER 3: FORCED CONVECTION FLOW ***** CHAPTER 4: NATURAL CONVECTION ***** See also Benchmark Problems B.02 to B.07 ***** CHAPTER 5: PROBLEMS WITH PRACTICAL PROJECT APPLICATIONS ***** CHAPTER B: BENCHMARK PROBLEMS; KNOWN NUMERICAL SOLUTIONS OR DATA ***** CHAPTER H: HYBRID GRID ELEMENTS - TEST CASES ***** CHAPTER T: TEST PROBLEMS FOR VARIOUS OPTIONS ***** CHAPTER U: UNSTRUCTURED HEXAGONAL GRID CASES ***** CHAPTER V: VERIFICATION PROBLEMS; KNOWN ANALYTIC SOLUTIONS ***** ***** For further questions of clarifications please contact: ***** ACRi, 1931 Stradella Road, Los Angeles, CA 90077 ***** Tel: 310-471-3023; Fax: 310-471-0797 ***** e-mail: acri@ACRiCFD.com //////////////////////////////////////////////////////////////////////// / / *11111111111111111111111111111111111111111111111111111111111111111111111 ***** CHAPTER 1: PURE CONDUCTION AND DIFFUSION *11111111111111111111111111111111111111111111111111111111111111111111111 / ************************************************************************ TITLE Problem 1.01 - Conduction in an Infinite Strip (2D Simulation) ************************************************************************ / GRID NODEs 12 by 12 COOR X RANGE = 10.0 COOR Y RANGE = 1.0 /----------------------------------------------------------------------/ / Initial and Boundary Conditions BOUND T VALU X- 100.0 BOUND T VALU X+ 0.0 SYMM Y+ SYMM Y- /----------------------------------------------------------------------/ / Fluid Properties and Constants DENS 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 / DIAGnostic output at (6,6) every 20 steps DEBUG GEOMERTY OFF SAVE OFF T on file '101.SAV' FLUX DEFAult output OFF /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-6 1 SOLVE T STEADY 200 / SELEct (1,1) to (999,999) with interval (2,1) OUTPut T for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.02 - Conduction in an Infinite Strip (3D Simulation) ************************************************************************ / GRID NODEs 12 by 12 by 12 COOR X RANGE = 1.0 COOR Y RANGE = 1.0 COOR Z RANGE = 10.0 / / Initial and Boundary Conditions BOUND T VALU Z- 100.0 BOUND T VALU Z+ 0.0 SYMM X+ SYMM X- SYMM Y+ SYMM Y- / / Fluid Properties and Constants DENS 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.0 / / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 / / ** NOTE ** IF YOU USE CONJ GRAD SOLVER, CHANGE T TO H IN CONV COMMAND /MATR NSPC T CHOL CONJ DIAG / DIAGnostic output at (6,6,6) every 20 steps DEBUG GEOMERTY OFF SAVE OFF T on file '102.SAV' FLUX DEFAult output OFF / CONV T REFE GLOB 1.E-6 200 / SELEct (1,1,1) to (999,999,999) with interval (2,3,3) OUTPut T for SELEcted subdomain in NARRow mode / SOLVE T STEADY 20 / END / ************************************************************************ TITLE Problem 1.03 - Conduction in a Rotated Infinite Strip (2D) ************************************************************************ / GRID NODEs 12 by 12 COOR X RANGE = 10.0 COOR Y RANGE = 1.0 GEOMetry ROTAte X and Y by 30 degrees /----------------------------------------------------------------------/ / Initial and Boundary Conditions BOUND T VALU X- 100.0 BOUND T VALU X+ 0.0 SYMM Y+ SYMM Y- /----------------------------------------------------------------------/ / Fluid Properties and Constants DENS 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 / DIAGnostic output at (6,6) every 20 steps DEBUG GEOMERTY OFF SAVE OFF T on file '103.SAV' FLUX DEFAult output OFF /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-6 1 SOLVE T STEADY 200 / SELEct (1,1) to (999,999) with interval (2,1) OUTPut T for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.04 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 by 12 COOR X RANGE = 1.0 COOR Y RANGE = 1.0 COOR Z RANGE = 10.0 GEOMetry ROTAte X and Z by 30 degrees /----------------------------------------------------------------------/ / Initial and Boundary Conditions / BOUND T VALU Z- 100.0 BOUND T VALU Z+ 0.0 SYMM X+ SYMM X- SYMM Y+ SYMM Y- /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 / / ** NOTE ** IF YOU USE CONJ GRAD SOLVER, CHANGE T TO H IN CONV COMMAND /MATR NSPC T CHOL CONJ DIAG / DIAGnostic output at (6,6,6) every 20 steps DEBUG GEOMERTY OFF SAVE OFF T on file '104.SAV' FLUX DEFAult output OFF /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-6 200 / SELEct (1,1,1) to (999,999,999) with interval (2,3,3) OUTPut T for SELEcted subdomain in NARRow mode / SOLVE T STEADY 20 / END / ************************************************************************ TITLE Problem 1.05 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 COORdinate corners (0,0) (1,1) (0,1) (1,2) /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T LINE X 100.0 -100.0 / BOUND T VALU X- 100.0 BOUND T VALU X+ 0.0 BOUND T VALU Y- LINE X 100.0 -100.0 BOUND T VALU Y+ LINE X 100.0 -100.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 DIFF SECO SKEW T / DIAGnostic output for T and H at (6,6) every 20 steps DEBUG GEOMERTY OFF / CONV T REFE GLOB 1.E-8 1 FLUX DEFAult output OFF / SOLVE T STEADY 200 / SET FF LINE X 100.0 -100.0 SET FF LINE T 0.0 -1.0 ADD to existing SAVE OFF T FF on file '105.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.06 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 COORdinate corners (0,0) (1,1) (0,1) (1,2) /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T LINE X 100.0 -100.0 / BOUND T FLUX X- 100.0 BOUND T FLUX X+ -100.0 BOUND T VALU Y- LINE X 100.0 -100.0 BOUND T VALU Y+ LINE X 100.0 -100.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS / RELA T 1.0 LIMI T -1.E30 1.E30 DIFF SECO SKEW H T / 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-8 1 SOLVE T STEADY 200 / SET FF LINE X 100.0 -100.0 SET FF LINE T ADD 0.0 -1.0 SAVE OFF T FF on file '106.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.07 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 COORdinate corners (0,0) (1,1) (0,1) (1,2) /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T LINE X 100.0 -100.0 / BOUND T VALU X- 100.0 BOUND T VALU X+ 0.0 BOUND T FLUX Y- -70.710678 BOUND T FLUX Y+ 70.710678 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 DIFF SECO SKEW / 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-8 1 SOLVE T STEADY 200 / SET FF LINE X 100.0 -100.0 SET FF LINE T ADD 0.0 -1.0 SAVE OFF T FF on file '107.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.08 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 COORdinate corners (0,0) (10,1) (0,1) (10,2) /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T LINE X 100.0 -10.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 DIFF SECO SKEW H T / 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-8 1 SOLVE T STEADY 200 / SET FF LINE X 100.0 -10.0 SET FF LINE T ADD 0.0 -1.0 SAVE OFF T FF on file '108.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.09 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 COORdinate corners (0,0) (1,3) (0,1) (1,4) /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T LINE X 100.0 -100.0 / BOUND T VALU X- 100.0 BOUND T VALU X+ 0.0 BOUND T VALU Y- LINE X 100.0 -100.0 BOUND T VALU Y+ LINE X 100.0 -100.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 DIFF SECO SKEW H T / 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-8 1 SOLVE T STEADY 200 / SET FF LINE X 100.0 -100.0 SET FF LINE T ADD 0.0 -1.0 SAVE OFF T FF on file '109.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.10 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 COORdinate corners (0,0) (1,3) (0,1) (1,4) /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T LINE Y 100.0 -100.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 DIFF SECO SKEW H T / 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-8 1 SOLVE T STEADY 200 / SET FF LINE Y 100.0 -100.0 SET FF LINE T ADD 0.0 -1.0 SAVE OFF T FF on file '110.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.11 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 COORdinate corners (0,0) (1,0) (0,1) (1,2) /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T LINE X 0.0 5.0 SET T LINE Y ADD 0.0 6.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 DIFF SECO SKEW H T / 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 LINE X 0.0 5.0 SET FF LINE Y ADD 0.0 6.0 SET FF LINE T ADD 0.0 -1.0 SAVE OFF T FF on file '111.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.12 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 GEOM is a CYLINDER with radius of 0.5 /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T LINE Y 100.0 -100.0 BOUND T VALU X- LINE Y 100.0 -100.0 BOUND T VALU X+ LINE Y 100.0 -100.0 BOUND T VALU Y- LINE Y 100.0 -100.0 BOUND T VALU Y+ LINE Y 100.0 -100.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.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 / CONV T REFE GLOB 1.E-8 1 SOLVE T STEADY 200 / SET FF LINE Y 100.0 -100.0 SET FF LINEar function: 0. -1 * T ADD to existing SAVE OFF T FF on file '112.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.13 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 GEOM is a CYLINDER with radius of 0.5 /----------------------------------------------------------------------/ / Initial and Boundary Conditions / BOUND T VALU X- 100.0 BOUND T VALU X+ 0.0 BOUND T FLUX Y- 0.0 BOUND T FLUX Y+ 0.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS 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 T on file '113.SAV' /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-8 1 SOLVE T STEADY 200 SELEct (1,1) to (999,999) with interval (2,2) OUTPut T for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.14 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 GEOM is a CYLINDER with radius of 0.5 /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T LINE X 100.0 -100.0 BOUND T VALU X- LINE X 100.0 -100.0 BOUND T VALU X+ LINE X 100.0 -100.0 BOUND T VALU Y- LINE X 100.0 -100.0 BOUND T VALU Y+ LINE X 100.0 -100.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.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-8 1 SOLVE T STEADY 200 / SET FF LINE X 100.0 -100.0 SET FF LINEar function: 0. -1 * T ADD to existing SAVE OFF T FF on file '114.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.15 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 GEOMetry is an ANNULUS with outer rad=1. inner rad = 0.5 / /----------------------------------------------------------------------/ / Initial and Boundary Conditions SET T LINE Y 100.0 -100.0 BOUND T VALU X- LINE Y 100.0 -100.0 BOUND T VALU X+ LINE Y 100.0 -100.0 BOUND T VALU Y- LINE Y 100.0 -100.0 BOUND T VALU Y+ LINE Y 100.0 -100.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.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 / CONV T REFE GLOB 1.E-8 1 SOLVE T STEADY 200 / SET FF LINE Y 100.0 -100.0 SET FF LINEar function: 0. -1 * T ADD to existing SAVE OFF T FF on file '115.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.16 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 GEOMetry is an ANNULUS with outer rad=1. inner rad = 0.5 /----------------------------------------------------------------------/ / Initial and Boundary Conditions / BOUND T VALU X- 100.0 BOUND T VALU X+ 0.0 BOUND T FLUX Y- 0.0 BOUND T FLUX Y+ 0.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS 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 T on file '116.SAV' /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-8 1 SOLVE T STEADY 200 SELEct (1,1) to (999,999) with interval (2,2) OUTPut T for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.17 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 GEOMetry is an ANNULUS with outer rad=1. inner rad = 0.5 /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T LINE X 100.0 -100.0 BOUND T VALU X- LINE X 100.0 -100.0 BOUND T VALU X+ LINE X 100.0 -100.0 BOUND T VALU Y- LINE X 100.0 -100.0 BOUND T VALU Y+ LINE X 100.0 -100.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.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-8 1 SOLVE T STEADY 200 / SET FF LINE X 100.0 -100.0 SET FF LINEar function: 0. -1 * T ADD to existing SAVE OFF T FF on file '117.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.18 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 GEOMetry is an ANNULUS with outer rad=1. inner rad = 0.5 /----------------------------------------------------------------------/ / Initial and Boundary Conditions / BOUND T VALU Y- 0.0 BOUND T VALU Y+ 100.0 BOUND T FLUX X- 0.0 BOUND T FLUX X+ 0.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 1.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 T on file '118.SAV' /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-8 1 SOLVE T STEADY 200 SELEct (1,1) to (999,999) with interval (2,2) OUTPut T for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.19 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 GEOMetry is ANNUlus with outer rad = 2, inner = 1 offset = 0.5 /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T LINE X 0.0 5.0 SET T LINE Y ADD 0.0 6.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 DIFF SECO SKEW T / 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 LINE X 0.0 5.0 SET FF LINE Y ADD 0.0 6.0 SET FF LINE T ADD 0.0 -1.0 SAVE OFF T FF on file '119.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.20 - CONDUCTION ************************************************************************ / GRID NODEs 22 by 22 GEOMetry is ANNUlus with outer rad = 2, inner = 1 offset = 0.5 /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T LINE X 0.0 5.0 SET T LINE Y ADD 0.0 6.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 DIFF SECO SKEW H T / 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 LINE X 0.0 5.0 SET FF LINE Y ADD 0.0 6.0 SET FF LINE T ADD 0.0 -1.0 SAVE OFF T FF on file '120.SAV' SELEct (1,1) to (999,999) with interval (4,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.21 - CONDUCTION ************************************************************************ / GRID NODEs 42 by 42 GEOMetry is ANNUlus with outer rad = 2, inner = 1 offset = 0.5 /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T LINE X 0.0 5.0 SET T LINE Y ADD 0.0 6.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 DIFF SECO SKEW H T / DIAGnostic output 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 LINE X 0.0 5.0 SET FF LINE Y ADD 0.0 6.0 SET FF LINE T ADD 0.0 -1.0 SAVE OFF T FF on file '121.SAV' SELEct (1,1) to (999,999) with interval (4,4) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.22 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 COOR X RANGE = 4.0 COOR Y RANGE = 2.0 /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T POWE X 5.0 0.0 2.0 0.0 SET T POWE Y ADD 6.0 0.0 2.0 0.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 /----------------------------------------------------------------------/ / Source and Sink Specifications SOUR T VOLU -132.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 POWE X 5.0 0.0 2.0 0.0 SET FF POWE Y ADD 6.0 0.0 2.0 0.0 SET FF LINE T ADD 0.0 -1.0 SAVE OFF T FF on file '122.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.23 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 COOR X RANGE = 4.0 COOR Y RANGE = 2.0 GEOMetry ROTAte X and Y by 30 degrees /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T POLYnomial in X 0. 0. 5.0 SET T POLYnomial in Y 0. 0. 6.0 ADD /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 /----------------------------------------------------------------------/ / Source and Sink Specifications SOUR T VOLU -132.0 / 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 POLYnomial in X 0. 0. 5.0 SET FF POLYnomial in Y 0. 0. 6.0 ADD SET FF LINE T ADD 0.0 -1.0 SAVE OFF T FF on file '123.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.24 - CONDUCTION ************************************************************************ / GRID NODEs 82 by 82 COORdinate corners (0,0) (4,2) (3,6) (8,8) /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET T POLYnomial in X 0. 0. 5.0 SET T POLYnomial in Y 0. 0. 6.0 ADD /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 /----------------------------------------------------------------------/ / Source and Sink Specifications SOUR T VOLU -132.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS / RELA T 1.0 LIMI T -1.E30 1.E30 DIFF SECO SKEW H T / DIAGnostic output at (6,6) every 20 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-10 1 / SELEct (1,1) to (999,999) with interval (8,8) OUTPut FF for SELEcted subdomain in NARRow mode SAVE OFF T FF on file '124.SAV' / SOLVE T STEADY 200 / SET FF POLYnomial in X 0. 0. 5.0 SET FF POLYnomial in Y 0. 0. 6.0 ADD SET FF LINE T ADD 0.0 -1.0 / END / ************************************************************************ TITLE Problem 1.25 - CONDUCTION ************************************************************************ GRID NODEs 12 by 12 GEOM CIRCular arc: r = 1, angle=10, starting at 85 degrees /----------------------------------------------------------------------/ / Initial and Boundary Conditions SYMM X+ SYMM X- BOUND T VALU Y+ 5.0 SYMM Y- /----------------------------------------------------------------------/ / Fluid Properties and Constants DENS 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 / SOUR T VOLU -120.0 / RELA T 1.0 LIMI T -1.E30 1.E30 / DIAGnostic output at (6,6) every 50 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-10 1 SOLVE T STEADY 200 / SET FF POWEr law in R 5. 0. 2. SET LINEar function: FF = 0. -1. * T ADD to existing value SAVE OFF T FF on file '125.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.26 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 GEOM CIRCular arc: r = 1, angle=10, starting at 85 degrees /----------------------------------------------------------------------/ / Initial and Boundary Conditions PERI X BOUND T VALU Y+ 5.0 SYMM Y- /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 /----------------------------------------------------------------------/ / Source and Sink Specifications SOUR T VOLU -120.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS / RELA T 1.0 LIMI T -1.E30 1.E30 / DIAGnostic output at (6,6) every 50 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-10 1 SOLVE T STEADY 200 / SET FF POWEr law in R 5. 0. 2. 0. SET LINEar function: FF = 0. -1. * T ADD to existing value SAVE OFF T FF on file '126.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.27 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 GEOM CIRCular arc: r = 1, angle=10, starting at 85 degrees /----------------------------------------------------------------------/ / Initial and Boundary Conditions SYMM X+ SYMM X- BOUND T VALU Y+ 5.0 SYMM Y- /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 / DIAGnostic output at (6,6) every 50 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-7 SOLVE T STEADY 200 SET LINEar function: FF = 5. -1. * T SAVE OFF T FF on file '127.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.28 - CONDUCTION ************************************************************************ / GRID NODEs 12 by 12 GEOM CIRCular arc: r = 1, angle=10, starting at 85 degrees /----------------------------------------------------------------------/ / Initial and Boundary Conditions PERI X BOUND T VALU Y+ 5.0 SYMM Y- /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA T 1.0 LIMI T -1.E30 1.E30 / DIAGnostic output at (6,6) every 50 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-7 SOLVE T STEADY 200 / SET LINEar function FF = 5. -1. * T SAVE OFF T FF on file '128.SAV' SELEct (1,1) to (999,999) with interval (2,2) OUTPut FF for SELEcted subdomain in NARRow mode / END / ************************************************************************ TITLE Problem 1.29 - CONDUCTION ************************************************************************ / GRID NODEs 42 by 42 GEOM CIRCular arc: r = 1, angle=10, starting at 85 degrees /----------------------------------------------------------------------/ / Initial and Boundary Conditions PERI X BOUND T VALU Y+ 5.0 SYMM Y- /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 / SOUR T VOLU -120.0 / RELA T 1.0 LIMI T -1.E30 1.E30 / DIAGnostic output at (6,6) every 50 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF SELEct (1,1) to (999,999) with interval (4,4) OUTPut FF for SELEcted subdomain in NARRow mode / CONV T REFE GLOB 1.E-7 / SOLVE T STEADY 200 / SET FF POWEr law in R 5. 0. 2. 0. SET LINEar function FF = 0. -1. * T ADD to existing value SAVE OFF T FF on file '129.SAV' / END / ************************************************************************ TITLE Problem 1.30 - CONDUCTION ************************************************************************ / GRID NODEs 42 by 42 GEOM CIRCular arc: r = 1, angle=10, starting at 85 degrees /----------------------------------------------------------------------/ / Initial and Boundary Conditions PERI X BOUND T VALU Y+ 5.0 SYMM Y- /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 /----------------------------------------------------------------------/ / SOLUTION OPTIONS 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 SELEct (1,1) to (999,999) with interval (4,4) OUTPut FF for SELEcted subdomain in NARRow mode /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-10 200 SOLVE T STEADY 200 / SET LINEr function FF = 5.0 -1. * T SAVE OFF T FF on file '130.SAV' / END / / ************************************************************************ TITLE Problem 1.31 - CONDUCTION WITH INTERNAL OBSTACLES ************************************************************************ / GRID NODEs 22 by 22 COORdinate corners (0.,0.) (1.,0.) (0,1.) (1.,1.) /----------------------------------------------------------------------/ SELEct (7,1) (16,22) BLOCKAGE at SELEcted zone /----------------------------------------------------------------------/ / Initial and Boundary Conditions / BOUND T VALU CONS X- 10.0 BOUND T VALU CONS X+ 20.0 BOUND T FLUX CONS Y- 0.0 BOUND T FLUX CONS Y+ 0.0 ADIA OBST /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 3.0 SPEC 1.0 PRAN EFFE 0.2 /----------------------------------------------------------------------/ / SOLUTION OPTIONS 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 SELEct (1,1) to (999,999) with interval (4,2) OUTPut T for SELEcted subdomain in NARRow mode SAVE OFF T on file '131.SAV' /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV T REFE GLOB 1.E-6 1 SOLVE T STEADY 200 / END / / ************************************************************************ TITLE Problem 1.32 - CONDUCTION IN R-X POLAR COORDINATES ***** T = 100*LN(R) + 100 ************************************************************************ / GRID NODEs 22 by 22 COOR CORN (1,1) (4,1) (1,8) (4,8) CYLI / SET log function: T = 100 * LN ( 1. * Y + 0. ) + 100. / DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 / RELA T 1.0 LIMI T -1.E30 1.E30 DIFF SECO SKEW H T / DIAGnostic output at (6,6) every 20 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF SELEct (1,1) to (999,999) with interval (4,2) OUTPut FF for SELEcted subdomain in NARRow mode / CONV T REFE GLOB 1.E-10 1 SOLVE T STEADY 200 / SET log function: FF = 100 * LN ( 1. * Y + 0. ) + 100. SET LINEar function: FF = 0 -1 * T ADD to existing SAVE OFF T FF on file '132.SAV' / END / ************************************************************************ TITLE Problem 1.33 - CONDUCTION IN R-THETA POLAR COORDINATES ***** T = 100*LN(R) + 100 ************************************************************************ / GRID NODEs 5 by 12 by 12 COOR CORN (0.,1.,0.) (1.,1.,0.) (0.,2.,0.) (1.,2.,0.) (0.,1.,.7) (1.,1.,.7) (0.,2.,.7) (1.,2.,.7) CYLI / SET log function: T = 100 * LN ( 1. * Y + 0. ) + 100. SYMM X- SYMM X+ / DENSity 5.0 VISC 3.0 SPEC 4.0 PRAN EFFE 2.0 / RELA T 1.0 LIMI T -1.E30 1.E30 DIFF SECO SKEW H T / DIAGnostic output T H at (3,6,6) every 20 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF SELEct (1,1,3) to (999,999,999) with interval (2,2,2) OUTPut FF for SELEcted subdomain in NARRow mode SAVE OFF T FF on file '133.SAV' / CONV T REFE GLOB 1.E-10 1 / SOLVE T STEADY 1000 / SET log function: FF = 100 * LN ( 1. * Y + 0. ) + 100. SET LINEar function: FF = 0 -1 * T ADD to existing / END / / *22222222222222222222222222222222222222222222222222222222222222222222222 ***** CHAPTER 2: CONDUCTION WITH SPECIFIED CONVECTION *22222222222222222222222222222222222222222222222222222222222222222222222 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / ************************************************************************ 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 / / *33333333333333333333333333333333333333333333333333333333333333333333333 ***** CHAPTER 3: FORCED CONVECTION FLOW *33333333333333333333333333333333333333333333333333333333333333333333333 / ************************************************************************ TITLE Problem 3.01 - SOLID BODY ROTATION : Z-Y COORDINATES ************************************************************************ / GRID NODEs 3 by 22 by 22 COOR X RANGE 1.0 COOR R RANGE 1.0 COOR Z RANGE 90 DEGREES GEOM convert to CARTesian /----------------------------------------------------------------------/ / Initial and Boundary Conditions / SET W LINE Y 0.0 4.0 SET V LINE Z 0.0 -4.0 PERI Z THETA 90 /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 20 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF SELEt (2,1,1) to (2,999,999) with interval of (2,2,2) OUTPut FF for selected YZ planes SAVE OFF V W P FF on file '301.SAV' / /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV GLOB 1.E-6 SOLVE for 2000 STEPS IN STEADY MODE / / W = OMEGA * R; P = PREF + 0.5*RHO * OMEGA**2 * (R**2 -RREF**2) / In this case: PREF = P(2,2,2) set to 0.; RREF = 0.025 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 FF LINEar -0.05 -1.0 * P ADD to existing / END / / ************************************************************************ 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 / / ************************************************************************ 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 / / ************************************************************************ 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 / ************************************************************************ 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 / / ************************************************************************ 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 / / ************************************************************************ 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 / ************************************************************************ 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 / / ************************************************************************ 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 / / ************************************************************************ 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 / ************************************************************************ TITLE Problem 3.12 - 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 / / ************************************************************************ TITLE Problem 3.13 - 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 / ************************************************************************ TITLE Problem 3.14 - 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 / ************************************************************************ TITLE Problem 3.15: 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 / *********************************************************************** TITLE Problem 3.16: 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 / *********************************************************************** TITLE Problem 3.17: 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 / ************************************************************************ TITLE Problem 3.18: 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 / ************************************************************************ TITLE Problem 3.19: 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 / *********************************************************************** TITLE Problem 3.20: 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 / / *44444444444444444444444444444444444444444444444444444444444444444444444 ***** CHAPTER 4: NATURAL CONVECTION *44444444444444444444444444444444444444444444444444444444444444444444444 / ************************************************************************ TITLE Problem 4.01 - NATURAL CONVECTION IN TWO CONCENTRIC CYLS ************************************************************************ / GRID NODEs 42 by 42 GEOM ANNULus: R_out = 1.625, R_in=0.625 GEOMetry ROTAte X and Y by 270 degrees /----------------------------------------------------------------------/ / Initial and Boundary Conditions WALL at undefined outer boundaries SYMM X- SYMM X+ / BOUN U VALU CONS Y- 0.0 BOUN V VALU CONS Y- 0.0 BOUN T VALU CONS Y- 1.0 / BOUN U VALU CONS Y+ 0.0 BOUN V VALU CONS Y+ 0.0 BOUN T VALU CONS Y+ 0.0 / /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 0.7 /----------------------------------------------------------------------/ / Source and Sink Specifications SOUR V LINE T VOLU 0.0 6.8571429E4 / RELA U=0.3, V=0.3, P=0.1, T=1. DIFF SECO SKEW U V P T / DIAGnostic U V P T RP output at (6,6) every 20 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF SELEct (1,1) to (999,999) interval of (4,2) OUTPut U V T in SELECted zone SAVE OFF U V P T FF on file '403.SAV' /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV GLOB 1.E-5 1 SOLVE U V P T STEADY 1000 / /LOCAte (1,2) (42,2) /SET LINEar FF = 47.775572 -47.775572 * T in SELEcted subdomain / /LOCAte (1,41) (42,41) /SET LINEar FF = 0.0 +124.21649 * T in SELEcted subdomain / / END / * *55555555555555555555555555555555555555555555555555555555555555555555555 ***** CHAPTER 5: PROBLEMS WITH PRACTICAL PROJECT APPLICATIONS *55555555555555555555555555555555555555555555555555555555555555555555555 * / ************************************************************************ TITLE Problem 5.01 - AXISYMMETRIC SWIRLER FLOW >>> C60A ************************************************************************ / GRID NODEs 62 by 18 COOR X R 'ANSD501.GRD' /----------------------------------------------------------------------/ SELEct (1,13) (21,18) BLOCKAGE at SELEcted zone /----------------------------------------------------------------------/ / Initial and Boundary Conditions WALL at undefined outer boundaries SELEct (1,1) (1,12) INLE X- SELECted subdomain OUTL X+ SYMM Y- / SET U 37.0 SET V 0.0 SET W 44.0 SET K 24.3215 SET L 2.54E-3 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 3.48 VISC 2.679E-5 / DIAGnostic U V W P K E output at (30,10) every 50 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF SELEct (1,1) to (999,999) interval of (6,1) OUTPut U V W P K L in SELECted zone SAVE OFF U V W P K E L MU on file '501.SAV' /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV GLOB 1.E-6 1 SOLVE U V W P K E STEADY 1000 / END / ************************************************************************ TITLE Problem 5.02 - AXISYMMETRIC SWIRLER FLOW >>> C60B ************************************************************************ / GRID NODEs 62 by 18 by 7 COOR X Y Z 'ANSD502.GRD' /----------------------------------------------------------------------/ SELEct (1,13,1) (21,999,999) BLOCKAGE at SELEcted zone /----------------------------------------------------------------------/ / Initial and Boundary Conditions / WALL at undefined outer boundaries SELEct (1,1,1) (1,12,999) INLE X- at SELEcted subdomain OUTL X+ SYMM Y- SET U 37.0 SET V SINE THETA -44.0 1.0 0.0 0.0 SET W COSI THETA 44.0 1.0 0.0 0.0 SET K 24.3215 SET L 2.54E-3 PERIodic in Z THETA 5.729578 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 3.48 VISC 2.679E-5 /----------------------------------------------------------------------/ / SOLUTION OPTIONS /DIFF SECO SKEW U V W P / DIAGnostic output U V W P K E at (30,10,4) every 50 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF SELEct from (1,1,4) to (99,99,4) interval (6,2,1) OUTPut U V W P K L in XY planes SELEcted zone SAVE OFF U V W P K E L MU on file '502.SAV' /----------------------------------------------------------------------/ / OPERATIONAL CONTROL // ::::::> akr_note: Only partial convergence with input below CONV GLOB 1.E-06 1 !!! Change to 1.E-10 for full run SOLVE STEADY 200 min 200 !!! Change to 2000 steps for full run / END / *********************************************************************** TITLE CASE 5.04: Ramjet Combustion: Default 4 Step 8 Species Reaction *********************************************************************** / GRID NODEs 32 by 22 / COORDINATE X and R from 'ANSD504.GRD' CORNers / //////////////////////////////////////////////////////////////////////// * GEOMETRY SPECIFICATION //////////////////////////////////////////////////////////////////////// / WALL at undefined outer boundaries SYMMetry at Y- boundary / LOCATE (01,01) to (01,09) with ID=INLEt INLET X- at SELEcted boundary OUTLET X+ / / //////////////////////////////////////////////////////////////////////// * Initial and Boundary Conditions //////////////////////////////////////////////////////////////////////// / REACTION / SET U = 190. at ID=INLEt at X- boundary only SET T = 600 at ID=INLEt at X- boundary only SET FU = 0.05 at ID=INLEt at X- boundary only SET value of K = 1000 everywhere SET value of L = 0.06 everywhere LOCAte (02,02) to (15,11) FIELd nodes only SET T 2000 SELEcted subdomian !!! Start combustion / ADIABATIC WALLS LIMIT T to min = 300, max = 2500 / //////////////////////////////////////////////////////////////////////// * OUTPUT CONTROL //////////////////////////////////////////////////////////////////////// / / DIAGNOSTIC node (08,08) every 50 steps DEBUG GEOMERTY OFF FLUX DEFault output OFF SELEct from (1,1) to (999,999) interval (3,2) OUTPut U,V,P,T,K,L,FF,FU,H2O SELEcted zone / //////////////////////////////////////////////////////////////////////// * SOLUTION OPTIONS //////////////////////////////////////////////////////////////////////// / RELAx U = 0.3 V=0.3 W=0.3 RHO=0.3 P = 0.7 FF=1.0 CONVERGENCE Max mass source = 1.0E-8 CONVERGENCE PP 1.0E-9 max of 10 Iterations on P-Prime Eqn / SOLVE 1000 STEPS in STEADY MODE SAVE OFF ON '504.SAV' / END / *********************************************************************** TITLE CASE 5.05: Ramjet Combustion: Default 2 Step 8 Species Reaction *********************************************************************** / GRID NODEs 32 by 22 / COORDINATE X and R from 'ANSD504.GRD' CORNers / //////////////////////////////////////////////////////////////////////// * GEOMETRY SPECIFICATION //////////////////////////////////////////////////////////////////////// / WALL at undefined outer boundaries SYMMetry at Y- boundary / LOCATE (01,01) to (01,09) with ID=INLEt INLET X- at SELEcted boundary OUTLET X+ / SET U = 190. at ID=INLEt at X- boundary only SET T = 600 at ID=INLEt at X- boundary only SET FU = 0.05 at ID=INLEt at X- boundary only SET value of K = 1000 everywhere SET value of L = 0.06 everywhere / //////////////////////////////////////////////////////////////////////// * Initial and Boundary Conditions //////////////////////////////////////////////////////////////////////// / REACTION TWO STEP / LOCAte (02,02) to (15,11) FIELd nodes only SET T 2000 SELEcted subdomain !!! Start combustion / ADIABATIC WALLS LIMIT T to min = 300, max = 2500 / //////////////////////////////////////////////////////////////////////// * OUTPUT CONTROL //////////////////////////////////////////////////////////////////////// / DIAGNOSTIC node (08,08) every 50 steps DEBUG GEOMERTY OFF FLUX DEFault output OFF SELEct from (1,1) to (999,999) interval (3,2) OUTPut U,V,P,T,K,L,FF,FU,H2O SELEcted zone / //////////////////////////////////////////////////////////////////////// * SOLUTION OPTIONS //////////////////////////////////////////////////////////////////////// / RELAx U = 0.3 V=0.3 W=0.3 RHO=0.3 P = 0.7 FF=1.0 CONVERGENCE Max mass source = 1.0E-8 CONVERGENCE PP 1.0E-9 max of 10 Iterations on P-Prime Eqn / SOLVE 1000 STEPS in STEADY MODE SAVE OFF ON '505.SAV' / END / / *********************************************************************** TITLE CASE 5.06: Ramjet Combustion: User Specified 2 Step 6 Species *********************************************************************** / GRID NODEs 32 by 22 / COORDINATE X and R from 'ANSD504.GRD' CORNers / //////////////////////////////////////////////////////////////////////// * GEOMETRY SPECIFICATION //////////////////////////////////////////////////////////////////////// / WALL at undefined outer boundaries SYMMetry at Y- boundary / LOCATE (01,01) to (01,09) with ID=INLEt INLET X- at SELEcted boundary OUTLET X+ / SET U = 190. at ID=INLEt at X- boundary only SET T = 600 at ID=INLEt at X- boundary only SET FU = 0.05 at ID=INLEt at X- boundary only SET FF = 0.05 at ID=INLEt at X- boundary only SET value of K = 1000 everywhere SET value of L = 0.06 everywhere / //////////////////////////////////////////////////////////////////////// * Initial and Boundary Conditions //////////////////////////////////////////////////////////////////////// / / Set up combustion from scratch REACTION is for a HYDROCARBON fuel / LOCAte (02,02) to (15,11) FIELd nodes only SET T 2000 SELEcted subdomian !!! Start combustion / ADIABATIC WALLS LIMIT T to min = 300, max = 2500 / SPECIFIC heat Gordon-McBride relations: 'cphtDFLT.acr' / / Define molecular weights of the species (FU = C12H22.5) GAS FU =139.26293 CO = 28.01055 H2O= 18.01534 CO2= 44.00995 O2 = 31.99880 N2 = 28.01340 / / Allocate space for reaction rates, called R1, R2 ALLOCATE R1 ALLOCATE R2 / / Define reaction rates / / FU + (10/2 + 19/4) O2 => 10 CO + 19/2 H2O REACtion R1 FU O2 3.3E+14 2.7E+04 0.5 1.0 0.0 MASS TEMP EBU / / CO + 0.5 O2 => CO2 REACtion R2 CO O2 6.0E+08 1.25E+04 1.0 1.0 0.0 MASS TEMP EBU / / EBU coefficients for each reaction / EBU 3.0 10 !! start control after 10 steps EBU R1 FU 3.0 EBU R2 CO 3.0 !! Assume only CO control; default option /EBU R2 CO O2 4.0 1.75 !! 2 * MW_CO /MW_O2; CO & O2 control / / Source terms for FU and CO transport equations / / FU source term SOURce REACtion type FU LINEar -1.0 * R1 / / CO source term ( 2.018 = 12*CO/FU ) <<< C12H22.5 / CO source term ( 2.01134 = 10*CO/FU ) <<< C10H19 SOURce REACtion type CO LINEar +2.01134 * R1 -1.0 * R2 / / Algebraic equations for total H-atom (FH), total O-atom (FO), / CO2, H2O, O2, and N2 by atom balance. FF is total C-atom. / SET ALWAys LINEar FH = 0.0 +1.0 * FF SET ALWAys LINEar FO = 0.2331 -0.2331 * FF SET ALWAys LINEar CO2 = 3.1602 * FF -3.1602 * FU -1.5712 * CO + 0.0 SET ALWAys LINEar H2O = 1.2289 * FH -1.2289 * FU +0.0 SET ALWAys LINEar O2 = 1.0 * FO -0.5714 * CO -0.7273 * CO2 -0.8889 * H2O + 0.0 SET ALWAys LINEar N2 = -1.0 * FU -1.0 * CO -1.0 * CO2 -1.0 * H2O -1.0 * O2 + 1.0 / //////////////////////////////////////////////////////////////////////// * OUTPUT CONTROL //////////////////////////////////////////////////////////////////////// / / DIAGNOSTIC node (08,08) every 50 steps DEBUG GEOMERTY OFF FLUX DEFault output OFF SELEct from (1,1) to (999,999) interval (3,2) OUTPut U,V,P,T,K,L,FF,FU,H2O SELEcted zone / //////////////////////////////////////////////////////////////////////// * SOLUTION OPTIONS //////////////////////////////////////////////////////////////////////// / RELAx U = 0.3 V=0.3 W=0.3 RHO=0.3 P = 0.7 FF=1.0 CONVERGENCE Max mass source = 1.0E-8 CONVERGENCE PP 1.0E-9 max of 10 Iterations on P-Prime Eqn / SOLVE 1000 STEPS in STEADY MODE SAVE OFF ON '506.SAV' / END / *BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB ***** CHAPTER B: BENCHMARK PROBLEMS; KNOWN NUMERICAL SOLUTIONS OR DATA *BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB / ************************************************************************ TITLE Problem B.01: SQUARE CAVITY WITH MOVING WALL ************************************************************************ / GRID NODEs BY 22 BY 22 / COORDINATE X RANGE 1 COORDINATE Y RANGE 1 / WALL at all outer boundaries BOUNdary U: at boundary Y+, VALUE=1. / LAMINAR flow VISCOSITY 0.01 / DIAGNOSTIC NODE U V P RU RV AT (6,6) Every 20 DEBUG GEOMERTY OFF FLUX DEFAult output OFF SELEct from (1,1) to (999,999) interval (2,2) OUTPut for SELEcted window SAVE OFF U,V,P, ON 'B01.SAV' / CONVERGENCE REFERENCE for U in LOCAL mode 1.E-6 SOLVE FOR 1000 STEPS IN STEADY MODE / END / *********************************************************************** TITLE Problem B.02: ATURAL CONVECTION IN A RECTANGULAR ENCLOSURE ***** Height/Length Ratio of 10 ***** Heat Transfer by Laminar Natural Convection Within Rectangular ***** Enclosures: by M.E. Newell and F.W. Schmidt ***** ///// OPTION 1 ***** Natural Convection Simulated by BOUSsinesq approximation *********************************************************************** / GRID NODEs 22 by 32 / WALL at all outer boundaries / COORDINATE X RANGE = 0.1 COORDINATE Y RANGE = 1 / SET T=300 EVERYWHERE / BOUND T X- VALUE 300 BOUND T X+ VALUE 302 BOUND T Y+ FLUX = 0. BOUND T Y- FLUX = 0. / ////////////////////////////////////////////////////////////////////////// ***** Commands to simulate Natural Convection by BOUSsinesq approximation / DENSITY 1.178 everywhere !Reference temperature 300 / ***** From gas law for density: dRHO/dT = -RHO/T ***** Over the range of convern; dRHO/dT_avg = -1.178/301 = -0.0039136 / PROBlem BOUSINESQ T dRHO/dTemperature = -0.0039136, Tavg=301 GRAV 0. -10. / ////////////////////////////////////////////////////////////////////////// / VISCOSITY VALUE 3.086E-5 / DIAGNOSTIC NODE U V P RU RV 5,5 PRINT EVERY 15 STEPS / OUTPut in NARROW mode / MATRIX SWEEPS P=7 times / CONVERGENCE REFERENCE for V in LOCAL mode epsilon = 1.E-9 / SOLVE STEADY MAX 1000 STEPS / END / / / *********************************************************************** TITLE Problem B.03: ATURAL CONVECTION IN A RECTANGULAR ENCLOSURE ***** Height/Length Ratio of 10 ***** Heat Transfer by Laminar Natural Convection Within Rectangular ***** Enclosures: by M.E. Newell and F.W. Schmidt / ///// OPTION 2 ***** Convection Simulated as source term in V equation *********************************************************************** / GRID NODEs 22 by 32 / WALL at all outer boundaries / COORDINATE X RANGE = 0.1 COORDINATE Y RANGE = 1 / SET T=300 EVERYWHERE / BOUND T X- VALUE 300 BOUND T X+ VALUE 302 BOUND T Y+ FLUX = 0. BOUND T Y- FLUX = 0. / ////////////////////////////////////////////////////////////////////////// ***** Natural Convection by source term for Momentum equation / DENSITY 1.178 everywhere !Reference temperature 300 / ***** The value of the source coefficient is -drho/dt * g ***** No explicit specification of GRAVity; it is included in source term / SOURce for V LINEar function of 0. +0.039136 * T per unit VOLUme / ***** In this option a term equal to RHO_avg * G * Y is added to P ***** Hence the following to provide initial estimates of P and wall P / SET LINEAR P = 0.01 -11.78 * Y ! near hydrostatic for smooth start WALL P EXTRapolated Linearly at all walls / ////////////////////////////////////////////////////////////////////////// / VISCOSITY VALUE 3.086E-5 / DIAGNOSTIC NODE U V P RU RV 5,5 PRINT EVERY 15 STEPS / OUTPut in NARROW mode / MATRIX SWEEPS P=7 times / CONVERGENCE REFERENCE for V in LOCAL mode epsilon = 1.E-9 / SOLVE STEADY MAX 1000 STEPS / END / / / *********************************************************************** TITLE Problem B.04: NATURAL CONVECTION IN A RECTANGULAR ENCLOSURE ***** Height/Length Ratio of 10 ***** Heat Transfer by Laminar Natural Convection Within Rectangular ***** Enclosures: by M.E. Newell and F.W. Schmidt / ///// OPTION 3 ***** Convection Simulated by coupling of gravity and variable density *********************************************************************** / GRID NODEs 22 by 32 / WALL at all outer boundaries / COORDINATE X RANGE = 0.1 COORDINATE Y RANGE = 1 / SET T=300 EVERYWHERE / BOUND T X- VALUE 300 BOUND T X+ VALUE 302 BOUND T Y+ FLUX = 0. BOUND T Y- FLUX = 0. / ////////////////////////////////////////////////////////////////////////// ***** Natural Convection by coupling of gravity and variable density / GAS molecular weight 29 DENSITY GAS LAW / / GRAV 0. -10. / ***** In this option a term equal to RHO_avg * G * Y is added to P ***** Hence the following to provide initial estimates of P and wall P / SET LINEAR P = 0.01 -11.78 * Y ! near hydrostatic for smooth start WALL P EXTRapolated Linearly at all walls / ////////////////////////////////////////////////////////////////////////// / VISCOSITY VALUE 3.086E-5 / DIAGNOSTIC NODE U V P RU RV 5,5 PRINT EVERY 15 STEPS / OUTPut in NARROW mode / MATRIX SWEEPS P=7 times / CONVERGENCE REFERENCE for V in LOCAL mode epsilon = 1.E-9 / SOLVE STEADY MAX 1000 STEPS / END / / ************************************************************************ TITLE Problem B.05 - NATURAL CONVECTION IN A TRAPEZOIDAL BOX >>> C41 ***** K C Karki, 1986, Ph.D. Thesis, p.183 ************************************************************************ / GRID NODEs 42 by 42 COORdinate corners (0.,0.) (1.,-.2679491) (0,.3307291) (1.,.5986782) WALL at all outer boundaries /----------------------------------------------------------------------/ / Initial and Boundary Conditions / BOUN U VALU CONS X- 0.0 BOUN V VALU CONS X- 0.0 BOUN T VALU CONS X- 1.0 / BOUN U VALU CONS X+ 0.0 BOUN V VALU CONS X+ 0.0 BOUN T VALU CONS X+ 0.0 / BOUN U VALU CONS Y- 0.0 BOUN V VALU CONS Y- 0.0 BOUN T FLUX CONS Y- 0.0 / BOUN U VALU CONS Y+ 0.0 BOUN V VALU CONS Y+ 0.0 BOUN T FLUX CONS Y+ 0.0 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 0.7 /----------------------------------------------------------------------/ / Source and Sink Specifications SOUR V LINE T VOLU 0.0 1.4285714E5 /----------------------------------------------------------------------/ / SOLUTION OPTIONS / RELA U=0.3, V=0.3, P=0.1, T=1. DIFF SECO SKEW U V P T / DIAGnostic output U V P T RP at (6,6) every 100 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF SELEct from (1,1) to (999,999) interval (2,2) OUTPut for SELEcted window SAVE OFF U V P T on file 'B04.SAV' /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV GLOB 1.E-7 SOLVE U V P T STEADY 1000 / END / ************************************************************************ TITLE Problem B.06: NATURAL CONVECTION IN TWO ECCENTRIC CYLINDERS ***** (Eccentricity, E = 0.623) ***** K C Karki, 1986, Ph.D. Thesis, p.169 ************************************************************************ / GRID NODEs 42 by 42 GEOM ANNULus: R_out = 1.625, R_in=0.625 ecc = -0.623 GEOMetry ROTAte X and Y by 270 degrees /----------------------------------------------------------------------/ / Initial and Boundary Conditions SYMM X- SYMM X+ WALL Y- WALL Y+ / BOUN U VALU CONS Y- 0.0 BOUN V VALU CONS Y- 0.0 BOUN T VALU CONS Y- 1.0 / BOUN U VALU CONS Y+ 0.0 BOUN V VALU CONS Y+ 0.0 BOUN T VALU CONS Y+ 0.0 / LOCA ID=REG1 (1,2) (42,2) SET FF LINE T ID=REG1 0.5971946 -0.5971946 / LOCA ID=REG2 (1,41) (42,41) SET FF LINE T ID=REG2 0.0 1.5527061 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 0.7 /----------------------------------------------------------------------/ / Source and Sink Specifications SOUR V LINE T VOLU 0.0 7.E4 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA U=0.3, V=0.3, P=0.1, T=1. DIFF SECO SKEW U V P T / DIAGnostic output U V P T RP at (6,6) every 50 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF SELEct from (1,1) to (999,999) interval (2,2) OUTPut U, V, T for SELEcted window SAVE OFF U V P T FF on file 'B05.SAV' /----------------------------------------------------------------------/ / OPERATIONAL CONTROL CONV GLOB 1.E-5 1 SOLVE U V P T STEADY 1000 / END / ************************************************************************ TITLE Problem B.07: NATURAL CONVECTION IN TWO ECCENTRIC CYLINDERS ***** (Eccentricity, E = -0.652) ***** K C Karki, 1986, Ph.D. Thesis, p.169 ************************************************************************ / GRID NODEs 42 by 42 GEOM ANNULus: R_out = 1.625, R_in=0.625 ecc = 0.652 GEOMetry ROTAte X and Y by 270 degrees /----------------------------------------------------------------------/ / Initial and Boundary Conditions SYMM at X- SYMM at X+ WALL at Y- WALL at Y+ / BOUN U VALU CONS Y- 0.0 BOUN V VALU CONS Y- 0.0 BOUN T VALU CONS Y- 1.0 / BOUN U VALU CONS Y+ 0.0 BOUN V VALU CONS Y+ 0.0 BOUN T VALU CONS Y+ 0.0 / LOCA ID=REG1 (1,2) (42,2) SET FF LINE T ID=REG1 0.5971946 -0.5971946 / LOCA ID=REG2 (1,41) (42,41) SET FF LINE T ID=REG2 0.0 1.5527061 /----------------------------------------------------------------------/ / Fluid Properties and Constants DENSity 1.0 VISC 1.0 SPEC 1.0 PRAN EFFE 0.7 /----------------------------------------------------------------------/ / Source and Sink Specifications SOUR V LINE T VOLU 0.0 6.8571429E4 /----------------------------------------------------------------------/ / SOLUTION OPTIONS RELA U=0.3, V=0.3, P=0.1, T=1. DIFF SECO SKEW U V P T / DIAGnostic output U V P T RP at (6,6) every 50 steps DEBUG GEOMERTY OFF FLUX DEFAult output OFF SELEct from (1,1) to (999,999) interval (4,2) OUTPut U, V, T for SELEcted window SAVE OFF U V P T FF on file 'B06.SAV' / CONV GLOB 1.E-5 1 / SOLVE U V P T STEADY 2000 / END / ************************************************************************ TITLE Problem B.08 - Coaxial Swirling Flow **** McDonell & Samuelsen (570 swirler) ************************************************************************ / GRID 48 by 44 / WALLS by DEFAult / COORdinate X CORNers 0.0000E+00 1.5500E-03 3.2500E-03 5.6500E-03 8.2500E-03 1.1050E-02 1.4125E-02 1.7675E-02 2.1900E-02 2.6850E-02 3.2400E-02 3.8400E-02 4.4875E-02 5.1825E-02 5.9000E-02 6.6225E-02 7.3500E-02 8.0800E-02 8.8100E-02 9.5400E-02 1.0267E-01 1.0995E-01 1.1725E-01 1.2455E-01 1.3185E-01 1.3915E-01 1.4645E-01 1.5375E-01 1.6105E-01 1.6835E-01 1.7565E-01 1.8295E-01 1.9025E-01 1.9755E-01 2.0485E-01 2.1215E-01 2.1945E-01 2.2675E-01 2.3385E-01 2.4115E-01 2.4865E-01 2.5595E-01 2.6330E-01 2.7060E-01 2.7785E-01 2.8515E-01 2.9245E-01 / COORdinate R CORNers CYLIndrical mode 0.0000E+00 1.0000E-03 3.0000E-03 5.0000E-03 7.0000E-03 9.0000E-03 1.1000E-02 1.3000E-02 1.5000E-02 1.6500E-02 1.7500E-02 1.8500E-02 1.9500E-02 2.0500E-02 2.1500E-02 2.3000E-02 2.5000E-02 2.7000E-02 2.9000E-02 3.1000E-02 3.3000E-02 3.5000E-02 3.7000E-02 3.9000E-02 4.1000E-02 4.3000E-02 4.5000E-02 4.7000E-02 4.9000E-02 5.1000E-02 5.3000E-02 5.5000E-02 5.7000E-02 5.9000E-02 6.1000E-02 6.3000E-02 6.5000E-02 6.7000E-02 6.9000E-02 7.1000E-02 7.3000E-02 7.5000E-02 7.6000E-02 / INLET X- OUTLET X+ SYMMETRY Y- / DENSity = 1.17 VISCOSITY 1.85E-5 / /TURBulence model KE /TURBulence model RNG TURBulence model CUBIc K-E in PASSIVE mode without Vel source terms / /INITIAL GUESSES SET U = 1.0 SET V = 0.0 SET W = 0.0 SET K = 1.0E-05 SET E = 1.0E-05 / /INLET VALUES from experimental data / SET U ON NODE X- -1.1310E+00 -1.1610E+00 -1.2400E+00 -1.2560E+00 -1.3250E+00 -1.3930E+00 -1.4620E+00 -5.3600E-01 2.7500E-01 2.6350E+00 4.8800E+00 6.4630E+00 6.5780E+00 4.8180E+00 2.7950E+00 9.5500E-01 1.0310E+00 1.1070E+00 1.2090E+00 1.3110E+00 1.4060E+00 1.5020E+00 1.6150E+00 1.7280E+00 1.7730E+00 1.8190E+00 1.8560E+00 1.8920E+00 1.9130E+00 1.9340E+00 1.9870E+00 2.0400E+00 2.0730E+00 2.1070E+00 2.1400E+00 2.1740E+00 2.2090E+00 2.2450E+00 2.2480E+00 2.2510E+00 1.7310E+00 1.3400E+00 ! -1.1210E+00 0.0000E+00 / SET V ON NODE X- 0.0000E+00 1.4100E-01 4.3300E-01 7.4300E-01 9.9900E-01 1.2010E+00 1.3680E+00 1.4340E+00 1.5550E+00 1.8240E+00 2.7710E+00 3.4440E+00 3.2030E+00 3.1900E+00 2.2110E+00 5.3800E-01 -4.2700E-01 -9.2400E-01 -9.5200E-01 -7.8800E-01 -7.0300E-01 -6.7800E-01 -5.9200E-01 -4.4400E-01 -2.9600E-01 -1.4800E-01 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 ! 0.0000E+00 0.0000E+00 / SET W ON NODE X- 2.7000E-02 1.0700E-01 2.0900E-01 3.2100E-01 3.9200E-01 4.6400E-01 5.3500E-01 1.1650E+00 1.7170E+00 4.4780E+00 7.1600E+00 9.0060E+00 7.8230E+00 5.6260E+00 3.2160E+00 1.7300E+00 9.0800E-01 8.6000E-02 1.4000E-01 1.9300E-01 1.4700E-01 1.0100E-01 7.9000E-02 5.8000E-02 1.0000E-03 -5.6000E-02 -5.9000E-02 -6.2000E-02 -7.8000E-02 -9.4000E-02 -1.0200E-01 -1.1100E-01 -1.5800E-01 -2.0400E-01 -2.5100E-01 -2.9800E-01 -3.7400E-01 -4.5100E-01 -5.5800E-01 -6.6600E-01 -5.1500E-01 -4.0200E-01 ! 0.0000E+00 0.0000E+00 / SET K ON NODE X- 1.4086E+00 1.3624E+00 1.3438E+00 1.2239E+00 1.2277E+00 1.2236E+00 1.2164E+00 3.2902E+00 6.4943E+00 1.2866E+01 2.1621E+01 2.3127E+01 2.3170E+01 2.2589E+01 1.7133E+01 1.2051E+01 6.0249E+00 1.9246E+00 8.7584E-01 4.2535E-01 3.9442E-01 3.6650E-01 4.3018E-01 4.0374E-01 3.9276E-01 3.8420E-01 3.9235E-01 4.0113E-01 4.0848E-01 4.1789E-01 4.2513E-01 4.3347E-01 4.5188E-01 4.7138E-01 4.9075E-01 5.1122E-01 5.4419E-01 5.8011E-01 6.3156E-01 6.8457E-01 7.5951E-01 8.2400E-01 ! 1.4253E+00 1.3479E-02 / SET E ON NODE X- 2.5307E+02 2.4072E+02 2.3581E+02 2.0496E+02 2.0593E+02 2.0489E+02 2.0310E+02 9.0345E+02 2.5053E+03 6.9859E+03 1.5219E+04 1.6836E+04 1.6883E+04 1.6252E+04 1.0735E+04 6.3327E+03 2.2387E+03 4.0419E+02 1.2408E+02 4.1994E+01 3.7498E+01 3.3587E+01 4.2710E+01 3.8834E+01 3.7261E+01 3.6049E+01 3.7203E+01 3.8458E+01 3.9520E+01 4.0894E+01 4.1961E+01 4.3201E+01 4.5984E+01 4.8991E+01 5.2042E+01 5.5331E+01 6.0770E+01 6.6885E+01 7.5977E+01 8.5741E+01 1.0020E+02 1.1323E+02 ! 2.5760E+02 2.3690E-01 / / DIAGNOSTIC NODE (15,5) print U,V,P,K,E every 20 steps DEBUG GEOMERTY OFF / CONVERGENCE epsilon = 1.0E-8 SOLVE for 1000 steps in STEADY mode !! CONVERGENCE epsilon = 1.0E-9 ! for full convergence !! SOLVE max 2000 steps in STEADY mode ! Fully converged in about 1200 / SELECT (1,1) to (999,999) frequency (4,4) OUTPut U,V,W,P,K,E NARRow mode SAVE OFF U,V,W,P,K,E,L,MU,RHO on 'B07.SAV' / END / ************************************************************************ TITLE Problem B.09: Split flow in a bifurcated channel with open bndries // S.R. Mathur & J.Y. Murthy 1997; Num. Heat Transfer, B, 32, p283-298 ************************************************************************ / /GRID NODEs 122 by 82 <<< used for a finer grid GRID NODEs 62 by 42 COOR X RANGE 6 COOR Y RANGE 4 DEBUG GEOMETRY OFF / WALL BY DEFAULT LOCATE COORDINATE (0.00, 0.00) TO (2.00, 3.00) BLOC SELECTED LOCATE COORDINATE (3.01, 0.00) TO (6.01, 3.00) BLOC SELECTED / LOCATE ID=INLET COORDINATE (0.0, 3.01) TO (0.06, 6.00) INLET X- OF ID=INLET LOCATE ID=OUT1 COORDINATE (5.90, 3.01) TO (6.00, 4.00) LOCATE ID=OUT2 COORDINATE (2.01, 0.00) TO (2.98, 0.06) / OPEN X+ OF ID=OUT1 BOUND X+ OF ID=OUT1 P = 0. /BOUND X+ OF ID=OUT1 GRAD U = 0. / OPEN Y- OF ID=OUT2 BOUND Y- OF ID=OUT2 P = 0. /BOUND Y- OF ID=OUT2 GRAD V = 0. / SET U X- OF ID=INLET POLYNOMIAL in Y -48., +28., -4. DENSity 1.0 VISC 0.01 / MATRIX P=32 / DIAGnostic output for U V P BP at (32,41) every 10 steps / FLUX DEFAult output OFF / SELECT (5,1) (999,999) (5,5) OUTPut U V P NARROW selected / CONV GLOB 1.E-9 SOLVE STEADY 2000 / END / ************************************************************************ TITLE Problem B.10 - TURBULENT DEVELOPING FLOW IN A CHANNEL /// P.L. Stephenson: Int J. H&M Trans. Vol 19, p.413-423, 1976 /// A Theoretical Study of Heat Transfer in Two-Dimensional Turbulent /// Flow in a Cicular Pipe and between Parallel & Diverging Plates ************************************************************************ /// This ANSWER input below explores the fully developed channel flow /// at Re=192000 described by Stephenson. Re number is defined on the /// basis of 2 times plate separation (H). /// /// Three options are explored. All are based on half channel. /// / 1. A normalized input with fixed end pressures / H=2, Mu=4/Re, dpdx=1.80330E-3 (leads to U_bulk=1) /// / 2. A peridic flow with source of momentum = expected dp/dx /// H=5, Mu=0.1593 and dp/dx= 6.7400E3 (leads to U_bulk=3059 cm/s) /// / 3. A developing flow with prescribed inflow at inlet; zero grad outlet /// H=5, U_bulk=3059 and Mu=0.1593 ************************************************************************ / /// Please replace TYPE specification by one of the following 3: / 1. NORMALIZED 2. PERIODIC 3. DEVELOPING / DEFINE TYPE = NORMALIZED <<< Put in the type of flow / **************************************************************** IF( TYPE = NORMALIZED ) THEN **************************************************************** COORdinate Y RANGE 1 in normalized units GRID NODEs 3 by 42 OPEN boundary at X- OPEN boundary at X+ / SET U 1 SET K 0.003 SET L 0.14 / BOUND P 1.80330E-3 at X- ! leads to U_bulk=1. BOUND P 0 at X+ BOUND K GRAD = 0 at X- BOUND K GRAD = 0 at X+ BOUND E GRAD = 0 at X- BOUND E GRAD = 0 at X+ / VISCosity 2.083333E-5 !1/48000 to get Re=192000 of Stephenson DIAGnostics for U K E at node (2,41) every 500 steps / ENDIF / **************************************************************** IF( TYPE = PERIODIC ) THEN **************************************************************** / GRID NODEs 3 by 42 PERIODIC X over unity length COORdinate Y RANGE 2.5 units are cm for this one SET U 3059 velocity in cm/s SET K 3.0E4 turbulent kinetic enegry SET L 0.35 length scale SOURCE for U = 6.74000E3 per unit VOLUME !=dpdx; leads to U_bulk=3059 VISCosity 0.1593 !cm2/s / DIAGnostics for U K E at node (2,41) every 500 steps / ENDIF / **************************************************************** IF( TYPE = DEVELOPING ) THEN **************************************************************** / GRID NODEs 102 by 42 COORdinate X RANGE 100.0 1.05 units are cm this one COORdinate Y RANGE 2.5 units are cm for this one INLET X- OUTL X+ SET U 3059 cm/s SET K 3.0E4 SET L 0.35 DIAGnostic U V P K E output at (61,41) every 100 steps VISCosity 0.1593 !cm2/s / ENDIF / **************************************************************** ////// COMMON SECTION FOR ALL OPTIONS **************************************************************** / / SYMMetry boundary at Y- surface WALL at otherwise undefined boundaries / DENSity of fluid is 1.0 RELAXation factor for U = 1 / DEBUG GEOMERTY option is OFF / PRINT SHEAR STRESS at WALLS with details FLUX AVERAGE OF P, U at X+ boundary / FLUX DEFAult output OFF SAVE U V P K E L MU on file '314.SAV' OUTPut U V P K E in SELEcted subdomain NARROW / CONVergence U LOCAL 1.E-12 / //////////////////////////////////////////////////////////////// ******> Get non-dimensionalized output at outlet //////////////////////////////////////////////////////////////// // **************************************************************** IF( TYPE = NORMALIZED ) THEN **************************************************************** / SOLVE STEADY 30000 30000 / WRITE U at X+ boundary; it is already normalized / /// dpdx = tau = 1.80330E-3 (Channel 1/2-height=1.) STACK scale WRITE K variable with 1.1090778E+3 ! = 2/tau WRITE K at X+ boundary and scale with STACK value / ENDIF / **************************************************************** IF( TYPE = PERIODIC ) THEN **************************************************************** SOLVE STEADY 30000 30000 / STACK scale WRITE U variable with 3.269E-4 ! = 1/3059 (=1/U_b) WRITE U at X+ boundary scale with STACK value / /// dpdx=6.74000E3 ==> tau=dpdx*2.5 = 1.685000E4 (Channel 1/2-height=2.5) STACK scale WRITE K variable with 1.1186944E-4 ! = 2/tau WRITE K at X+ boundary and scale with STACK value / ENDIF / **************************************************************** IF( TYPE = DEVELOPING ) THEN **************************************************************** / SOLVE STEADY 20000 2000 STACK scale WRITE U variable with 3.269E-4 ! = 1/3059 WRITE U at X+ boundary scale with STACK value / /// Here expected tau = 136.8^2 STACK scale WRITE K variable with 1.068E-4 ! = 2/tau WRITE K at X+ boundary and scale with STACK value / ENDIF / END / *HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH ***** CHAPTER H: HYBRID GRID ELEMENTS - TEST CASES *HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH ************************************************************************ TITLE T.101 Triangular Grid - Uni-direction Linear Conduction ************************************************************************ // Illustrates the use of grids with triangular elements // GRID UNSTructured 12 elements COORdinate x y from 'ANST101.XYZ' CONNectivity HYBRId 'ANST101.CNC' / LOCATE ID=XMINUS PAIR (1,1) (7,1) LOCATE ID=XPLUS PAIR (6,2) (12,3) LOCATE ID=YMINUS PAIR (1,2) (4,2) (5,2) LOCATE ID=YPLUS PAIR (7,3) (10,3) (11,3) DENSity 1.0 VISCosity 0.1 ****> For x-direction conduction these are active / BOUNdary T ID=XMINUS X- VALUE = 0. BOUNdary T ID=XPLUS X- VALUE = 3. BOUNdary T ID=YMINUS X- flux = 0. BOUNdary T ID=YPLUS X- flux = 0. ****> For y-direction conduction these are active / /BOUNdary T ID=XMINUS X- flux = 0. /BOUNdary T ID=XPLUS X- flux = 0. /BOUNdary T ID=YMINUS X- value = 0. /BOUNdary T ID=YPLUS X- value = 2. / DIAGnostic element 12 print T RT every 10 steps OUTPUT TABLES X Y T / MATRIX T COOR redblack / SOLVE T STEADY 1000 steps / SAVE OFF T 'T101_SAV.TMP' / END / / ************************************************************************ TITLE T.102 Triangular Grid - Uni-dir Conduction + Specified Convection ************************************************************************ / GRID UNSTructured 12 elements / COORdinate x y from 'ANST101.xyz' CONNectivity HYBRId 'ANST101.cnc' WRITE CORNERS / LOCATE ID=XMINUS PAIR (1,1) (7,1) LOCATE ID=XPLUS PAIR (6,2) (12,3) LOCATE ID=YMINUS PAIR (1,2) (4,2) (5,2) LOCATE ID=YPLUS PAIR (7,3) (10,3) (11,3) DENSity 1.0 VISCosity 0.1 ****> For x-direction conduction these are active SET U = 0.1 BOUNdary T ID=XMINUS X- VALUE = 0. BOUNdary T ID=XPLUS X- VALUE = 3. BOUNdary T ID=YMINUS X- flux = 0. BOUNdary T ID=YPLUS X- flux = 0. ****> For y-direction conduction these are active /SET V = 0.1 /BOUNdary T ID=XMINUS X- flux = 0. /BOUNdary T ID=XPLUS X- flux = 0. /BOUNdary T ID=YMINUS X- value = 0. /BOUNdary T ID=YPLUS X- value = 2. / DIAGnostic element 12 print T RT every 10 steps OUTPUT OFF OUTPUT TABLES X Y T / MATRIX T COOR redblack / SOLVE T STEADY 1000 steps / SAVE OFF T 'T102_SAV.TMP' / END / / ************************************************************************ TITLE T.103: Mixed Quad & Triangular Grid - X-dir Linear Conduction ************************************************************************ / GRID UNSTructured 5 elements COORdinates from 'ANST103.XYZ' CONNectivity HYBRid 'ANST103.CNC' / LOCATE ID=XMINUS1 PAIR (5,3) LOCATE ID=XMINUS2 PAIR (2,2) LOCATE ID=XPLUS1 PAIR (3,3) LOCATE ID=XPLUS2 PAIR (4,2) LOCATE ID=YMINUS1 PAIR (2,3) LOCATE ID=YMINUS2 PAIR (3,2) LOCATE ID=YPLUS1 PAIR (4,3) LOCATE ID=YPLUS2 PAIR (5,2) / LIMIT T MINIMUM -1.E30 MAX 1.E30 / **** ===> AKR: Note **** ===> bc's can be set explicitly or SET used or both; same results / /BOUNdary T LOCATE ID=XMINUS1 VALUE =-50. /BOUNdary T LOCATE ID=XMINUS2 VALUE = 0.0 /BOUNdary T LOCATE ID=XPLUS1 VALUE = 150. /BOUNdary T LOCATE ID=XPLUS2 VALUE = 100. /BOUNdary T LOCATE ID=YMINUS1 VALUE = 50. /BOUNdary T LOCATE ID=YMINUS2 VALUE = 150 /BOUNdary T LOCATE ID=YPLUS1 VALUE = 50. /BOUNdary T LOCATE ID=YPLUS2 VALUE =-50. / SET T LINEAR 0. +100 * X / **** ===> AKR: Note **** ===> If Diff Skew not used error a little higher; still pretty good /DIFFUSTION SKEW SECOND ORDER / /INLET ID=XMINUS X- /OUTLET ID=XPLUS X- /SYMMetry ID=YPLUS X- /WALL ID=YMINUS X- DENSity 1.0 VISCosity 0.1 /BOUNdary T ID=XMINUS X- VALUE = 0. /BOUNdary T ID=XPLUS X- VALUE = 3. /BOUNdary T ID=YMINUS X- flux = 0. /BOUNdary T ID=YPLUS X- flux = 0. DIAGnostic element 1 print T RT every 10 steps OUTPUT OFF OUTPUT TABLES X Y T / MATRIX T redblack !face / SOLVE T STEADY 100 steps / SAVE T 'T103_SAV.TMP' / END / / ************************************************************************ TITLE T.104 Mixed Hex and Quad Based Pyramid - Linear X Conduction ************************************************************************ / GRID UNSTructured 7 elements THREED COORdinate 'ANST104.XYZ' CONNectivity HYBRid 'ANST104.CNC' / /LOCATE ID=XMINUS1 PAIR (5,3) /LOCATE ID=XMINUS2 PAIR (2,2) /LOCATE ID=XPLUS1 PAIR (3,3) /LOCATE ID=XPLUS2 PAIR (4,2) /LOCATE ID=YMINUS1 PAIR (2,3) /LOCATE ID=YMINUS2 PAIR (3,2) /LOCATE ID=YPLUS1 PAIR (4,3) /LOCATE ID=YPLUS2 PAIR (5,2) / LIMIT T MINIMUM -1.E30 MAX 1.E30 / **** ===> AKR: Note **** ===> bc's can be set explicitly or SET used or both; same results / /BOUNdary T LOCATE ID=XMINUS1 VALUE =-50. /BOUNdary T LOCATE ID=XMINUS2 VALUE = 0.0 /BOUNdary T LOCATE ID=XPLUS1 VALUE = 150. /BOUNdary T LOCATE ID=XPLUS2 VALUE = 100. /BOUNdary T LOCATE ID=YMINUS1 VALUE = 50. /BOUNdary T LOCATE ID=YMINUS2 VALUE = 150 /BOUNdary T LOCATE ID=YPLUS1 VALUE = 50. /BOUNdary T LOCATE ID=YPLUS2 VALUE =-50. / SET T LINEAR 0. +100 * X / **** ===> AKR: Note **** ===> If Diff Skew is not used error a little higher; still pretty good DIFFUSTION SKEW SECOND ORDER / /INLET ID=XMINUS X- /OUTLET ID=XPLUS X- /SYMMetry ID=YPLUS X- /WALL ID=YMINUS X- DENSity 1.0 VISCosity 0.1 /BOUNdary T ID=XMINUS X- VALUE = 0. /BOUNdary T ID=XPLUS X- VALUE = 3. /BOUNdary T ID=YMINUS X- flux = 0. /BOUNdary T ID=YPLUS X- flux = 0. DIAGnostic element 1 print T RT every 10 steps OUTPUT OFF OUTPUT TABLES X Y Z T / MATRIX T redblack !face / SOLVE T STEADY 100 steps / SAVE OFF T 'T104_SAV.TMP' / END / / ************************************************************************ TITLE T.105: Laminar Developing Flow in a Channel ************************************************************************ ** DEFINE TYPE_ELEMENTS = TRIANGLES ! Change TRIANGLES to QUAD for structured IF( TYPE_ELEMENTS = QUAD ) THEN ** ================================================================== ** ** This Section Active only for Quad Grid ** ================================================================== ** ** GRID NODEs 42 by 22 COOR X RANGE 10.0 1.05 COOR Y RANGE 0.5 LOCAte ID=INLET from ( 1,1) ( 1,999) LOCAte ID=OUTLET from (999,1) (999,999) LOCAte ID=SYMMETRY from (1,1) (999,1) LOCAte ID=OUTFLD from (41,2) to (41,21) DIAGnostic output U V P T RP at (41,2) EVERY 50 ** ELSE ** ================================================================== ** ** This Section Active only for Triangular Grid ** ================================================================== ** ** GRID UNSTructured elements = 1600 CONNectivity HYBRid in file 'ANST105.CNC' COORdinates X Y VERTices in file 'ANST105.XYZ' / LOCA ID=INLET PAIR 1 1 41 1 81 1 121 1 161 1 201 1 241 1 281 1 321 1 361 1 401 1 441 1 481 1 521 1 561 1 601 1 641 1 681 1 721 1 761 1 LoCAte ID=OUTLET PAIR 840 3 880 3 920 3 960 3 1000 3 1040 3 1080 3 1120 3 1160 3 1200 3 1240 3 1280 3 1320 3 1360 3 1400 3 1440 3 1480 3 1520 3 1560 3 1600 3 LOCAte ID=SYMMETRY PAIR 801 2 802 2 803 2 804 2 805 2 806 2 807 2 808 2 809 2 810 2 811 2 812 2 813 2 814 2 815 2 816 2 817 2 818 2 819 2 820 2 821 2 822 2 823 2 824 2 825 2 826 2 827 2 828 2 829 2 830 2 831 2 832 2 833 2 834 2 835 2 836 2 837 2 838 2 839 2 840 2 / LOCAte ID=OUTFLD 840 880 920 960 1000 1040 1080 1120 1160 1200 1240 1280 1320 1360 1400 1440 1480 1520 1560 1600 / DIFF SKEW U V P T DIAGnostic output U V P T RP at ELEMENT 40 EVERY 50 ** ENDIF ** ================================================================== ** ** Common Section Active for both Quad and Triangular Grid ** ================================================================== ** ** WALL at undefined outer boundaries INLEt at ID=INLET X- OUTLet at ID=OUTLET X+ SYMMetry at ID=SYMMETRY Y- SET U at ID=INLET 100.0 SET T at ID=INLET 1.0 / DENSity 1.0 VISC 1.0 SPEC 1.0 / DEBUG GEOMERTY OFF FLUX DEFAult output OFF / OUTPut OFF / PRAN EFFE 0.7 RELAX U=0.4 CONV GLOB 1.E-12 / SOLVE U V P T STEADY 5000 / SAVE TABLE U, V, P, T ID=OUTFLD 'T105_TBL.TMP' / END / / ************************************************************************ TITLE T.106 Square Cavity with Moving Grid (Triangular & Quad Grids) ************************************************************************ ** DEFINE TYPE_ELEMENTS = TRIANGLES ! Change TRIANGLES to QUAD for structured IF( TYPE_ELEMENTS = QUAD ) THEN ** ================================================================== ** ** This Section Active only for Quad Grid ** ================================================================== ** ** GRID NODEs BY 22 BY 22 COORDINATE X RANGE 1 COORDINATE Y RANGE 1 LOCAte ID=TOP from (1,999) to (999,999) DIAGNOSTIC NODE U V P RU RV AT ELEMENT (12,12) Every 20 STEPS ** ELSE ** ================================================================== ** ** This Section Active only for Triangular Grid ** ================================================================== ** ** GRID UNSTructured elements = 800 CONNectivity HYBRid in file 'ANST106.CNC' COORdinates X Y VERTices in file 'ANST106.XYZ' LOCAte PAIR LIST ID=TOP 381 3 382 3 383 3 384 3 385 3 386 3 387 3 388 3 389 3 390 3 391 3 392 3 393 3 394 3 395 3 396 3 397 3 398 3 399 3 800 3 LOCATE LIST ID = MID 10 30 50 70 90 110 130 150 170 190 210 230 250 270 290 310 330 350 370 390 DIAGNOSTIC NODE U V P RU RV AT ELEMENT 211 Every 20 STEPS ** ENDIF ** ================================================================== ** ** Common Section Active for both Quad and Triangular Grid ** ================================================================== ** ** WALL at all outer boundaries BOUNdary U = 1 at ID=TOP Y+ LAMINAR flow VISCOSITY 0.01 / DEBUG GEOMERTY OFF FLUX DEFAult output OFF SAVE OFF OUTPut OFF / CONVERGENCE REFERENCE for U in LOCAL mode 1.E-6 SOLVE FOR 1000 STEPS IN STEADY MODE LOCAte ID=H_MID COORDINATE (0., 0.47) TO (1., 0.53) LOCAte ID=V_MID COORDINATE (0.47, 0.) TO (0.53, 1.) SAVE TABLE U V ID=H_MID 'H_MID.TMP' NOW SAVE TABLE U V ID=V_MID 'V_MID.TMP' NOW / END / / ************************************************************************ TITLE T107: Backward Facing Step. Re = 2.7 x 10^4. !Driver and Seegmiller 1985, AIAA J., Vol. 23, No. 2, pp. 163 - 171 ************************************************************************ ** DEFINE TYPE_ELEMENTS = TRIANGLES ! Change TRIANGLES to QUAD for structured IF( TYPE_ELEMENTS = QUAD ) THEN ** ========================================================== ** ** This Section for Quad Elements ** ========================================================== ** ** GRID 74 by 52 COORdinate X CORNers