r /N-3 f c// t_-C_ NASA Technical Memo_ra__ndum 1063.4.,!................ Proteus Three-Dimensiona!Navier-Stokes Computer Code-Version 1.0 ............................ Volume 3-Programmer's Reference - -- Charles E. Towne, John R. Schwab, and Trong T. Bui Lewis Research Center Cleveland, Ohio N94-15865 _ _---( (NASA-TM- I06361) PROTEUS _'-n THREE-DIMENSIONAL NAVIER-STOKES )_ COMPUTER COPE, VERSION I.O. VOLUME Uncl as _" 3: PROGRAMMER'S REFERENCE (NASA) m_>_ October 1993 :_ 282 p -L=__ G3134 0191156 N L L I L -- I l L ---- N/ A i ......................... _j w r k _ CONTENTS SUMMARY ................................................................. 3 !.0 INTRODUCTION .......................................................... 5 2.0 PROGRAM STRUCTURE ................................................... 7 2.1 FLOWCHART ........................................................ 7 2.2 SUBPROGRAM CALLING TREE ........................................ 10 2.3 PROGRAMMING CONWENTIONS AND NOTES ........................... 14 2.3.1 Computer & Language .............................................. 14 2.3.2 Fortran Variables ................... • .............................. 15 3.0 COMMON BLOCKS ...................................................... 19 3.1 COMMON BLOCK SUMMARY ......................................... 19 3.2 COMMON VARIABLES LISTED ALPHABETICALLY ........................ 19 3.3 COMMON VARIABLES LISTED SYMBOLICALLY ................... , ...... 39 4.0 PROTEUS SUBPROGRAMS ................................................ 51 51 4.1 SUBPROGRAM SUMMARY ............................................ 54 4.2 SUBPROGRAM DETAILS .............................................. 55 Subroutine ADI ....................................................... Subroutine AVISC1 .................................................... 56 59 Subroutine AVISC2 ......................................... ........... Subroutine BCDENS ................................................... 62 Subroutine BCELIM ................................................... 65 Subroutine BCF ....................................................... 66 Subroutine BCFLIN .................................................... 71 Subroutine BCGEN .................................................... 73 Subrt ,tine BCGRAD .................................................. 75 76 Subroutine BCIMET .... , .............................................. Subroutine BCMET .................................................... 77 Subroutine BCNVEL ................................................... 79 Subroutine BCPRES .................................................... 85 Subroutine BCQ . ...................................................... 92 Subroutine BCSET ......................................... "............. 96 Subroutine BCTEMP .................................................... 98 105 Subroutine BCLrVEL ....................... ........................... Subroutine BCVN .................................................... 109 Subroutine BCWEL .................................................. 111 Subroutine BCVI ................. •.................................... 115 Subroutine BCV2 .. ................................................. .. 116 Subroutine BCV3 ..................................................... 117 118 Subroutine BCWVEL ............................ ,...................... Subroutine BCIVEL ................................................... 122 Subroutine BC2VEL ................................................... 127 Subroutine BC3VEL ................................................... 132 r Subroutine BLIN ..................................................... 137 Subroutine BLKOUT ...................................... ... ......... 139 Subroutine BLK4 ..................................................... 140 142 Subroutine BLK4P ............................................... ...... Subroutine BLK5 ..................................................... 144 Subroutine BLK5P ..................................................... 145 Proteus 3-D Programmer's Reference Contents ! BLOCK DATA ...................................................... 146 Subroutine BLOUT ................................................... 148 Subroutine BVUP ...................................................... 151 Subroutine COEFC ............................................. .. .... 158 Subroutine COEFE1 ................................................... 161 Subroutine COEFE2 ................................................... 167 Subroutine COEFX ................................................... 170 Subroutine COEFY ' 175 Subroutine COEFZ .................................................... 180 Subroutine CONV 185 Subroutine CUBIC .................................................... 187 Subroutine EQSTAT .................................................. 189 Subroutine EXEC ..................................................... 191 Subroutine EXECT ................................................... 196 Subroutine FILTER ................................................... 197 Subroutine FTEMP ................................................... 199 Function GATHER ................................................... 202 Subroutine GEOM .................................................... 203 Subroutine INIT ...................................................... 206 Subroutine INITC .................................................... 207 Subroutine INPUT .................................................... 212 Function ISAMAX .................................................... 214 Function ISAMIN .................................................... 216 Function ISRCHEQ ................................................... 217 Function ISRCHFGT .................................................. 219 Function ISRCHFLT .................................................. 220 Subroutine KEINIT ................................................... 221 MAIN Program ...................................................... 223 Subroutine METS ..................................................... 226 Subroutine OUTPUT .................................................. 228 Subroutine OUTW 230 Subroutine PAK ...................................................... 234 Subroutine PERIOD ................................................... 236 Subroutine PLOT ..................................................... 238 Subroutine PRODCT .................................................. 240 Subroutine, PRTHST .................................................. 241 Subroutine PRTOUT .................................................. 242 Subroutine RESID .................................................... 244 Subroutine REST ..................................................... 247 Subroutine ROBTS .................................................... 250 Function SASUM ..................................................... 252 Subroutine SGEFA .................................................... 253 Subroutine SGESL .................................................... 254 Function SNRM2 ..................................................... 255 Subroutine SWDOWN .................................................. 257 Subroutine SWUP .-................................................... 259 Subroutine TBC ...................................................... 262 Subroutine TIMSTP ................................................... 264 Subroutine TREMAIN ................................................. 268 Subroutine TURBBL .................................................. 269 Subroutine TURBCH .................................................. 271 Subroutine UPDATE .................................................. 273 Subroutine UPDTKE .................................................. 275 Subroutine VORTEX .................................................. 277 Subroutine WHENFLT ................................................ 278 Subroutine YPLUSN " 279 REFERENCES ............................................................. 281 2 Contents Proteus 3-D Programmer's Reference PROTEUSTHREE-DLMENSIONAL NAVIER-STOKES COMPUTER CODE - VERSION 1.0 Volume 3- Programmer's Reference Charles E. Towne, John R. Schwab, Trong T. Bui National Aeronautics and Space Administration Lewis Research Center Cleveland, Ohio SUMMARY A computer code called Proteus has been developed to solve the three-dimensional, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this ef- fort has been to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation have been emphasized. The governing equations are written in Cartesian coordinates and transformed into generalized nonorthogonal body-fitted coordinates. They are solved by marching in time using a fully-'coupled alternating-direction-implicit solution procedure with generalized first- or second-order time differencing. The boundary conditions are also treated implicitly, and may be steady or unsteady. Spatially periodic boundary conditions are also available. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The program contains many operating options. The thin-layer or Euler equations may be solved as subsets of the Navier-Stokes equations. The energy equation may be elLrninated by the assumption of constant