Page 307 - Introduction to Computational Fluid Dynamics
P. 307

P1: ICD
                                                                                                 15:43
                                                                                   May 11, 2005
                            CB908/Date
            0521853265appb
                                         0 521 85326 5
                     286
                            Table B.1 (continued )                APPENDIX B. 1D CONDUCTION CODE
                            Variable      Meaning
                            T             Array containing temperature ( CorK)
                                                                   ◦
                            T1            Temperature at node 1
                            T1O           Temperature at node 1 at old time
                            T1SPEC        Logical – refers to T -boundary condition at node 1
                            THOMAS        Logical – refers to TDMA
                            TIMEMX        Maximum allowable time
                            TINF          Temperature T ∞
                            TINFO         Temperature T ∞ at old time
                            TINF1         Temperature T ∞ near node 1
                            TINFN         Temperature T ∞ near node N
                            TINF1O        Temperature T ∞ near node 1 at old time
                            TINFNO        Temperature T ∞ near node N at old time
                            TN            Temperature at node N
                            TNO           Temperature at node N at old time
                            TNSPEC        Logical – refers to T -boundary condition at node N
                            TO            Array containing temperature at old time
                            TTIME         Total current time
                            UNSTEADY      Logical – refers to unsteady-state calculation
                                                                    3
                            VOL           Array containing cell volume (m )
                            X             Coordinate of node P (m)
                            XCELL         Logical – refers to cell-face coordinate specification
                            XCF           Coordinate of cell face at w
                            XNODE         Logical – refers to node coordinate specification



                            B.3 File USER1D.FOR
                            This is the main control file at the command of the user. The first routine
                            PROGRAM ONED is the command routinefromwhere subroutine MAIN is called.
                            The latter is the first subroutine of the LIB1D.FOR file. When all operations are
                            completed, PROGRAM ONED calls the RESULT subroutine, which is a part of
                            the USER1D.FOR file.
                               Following the listing of the COM1D.FOR file, listings of two USER1D.FOR
                            files are given. They correspond to the two solved problems in Chapter 2. The reader
                            is advised to refer to these files as well as to Table B.1 to understand the description
                            of each routine in USER1D.FOR file.

                            BLOCK DATA This routine at the end of the USER1D.FOR file specifies all the
                            problem-dependent data such as properties, boundary conditions, and other control
                            parameters. It is assumed that all data are given in consistent units. Here, SI units
                            are used except for the grid data XCF or X, which are dimensionless. The physical
                            coordinates in meters are then evaluated by multiplying by AL (the domain length)
                            in PROGRAM ONED. Dimensionless specification provides better appreciation of
   302   303   304   305   306   307   308   309   310   311   312