result6

关于有限容积法,由西安交通大学陶文铨院士课题组的基本程序(FORTRAN)

               COMPUTATION  IN  CARTESIAN  COORDINATES
              ****************************************


   ITER       SMAX           SSUM          V(4,7)         T(4,7)
     0      0.000E+00      0.000E+00      1.200E+01      3.000E+02
     1      2.925E+00      4.470E-07      1.325E+01      3.597E+02
     2      1.519E+00      1.550E-06      1.777E+01      3.639E+02
     3      2.161E-01     -7.600E-07      1.635E+01      3.659E+02
     4      9.251E-02     -8.866E-07      1.562E+01      3.671E+02
     5      5.600E-02      9.313E-08      1.521E+01      3.676E+02
     6      3.087E-02     -1.788E-07      1.496E+01      3.677E+02
     7      1.698E-02      5.681E-07      1.482E+01      3.677E+02
     8      9.649E-03     -6.016E-07      1.475E+01      3.677E+02
     9      5.631E-03     -3.311E-07      1.472E+01      3.676E+02
    10      3.318E-03     -1.956E-07      1.470E+01      3.676E+02
    11      1.949E-03      3.374E-07      1.469E+01      3.676E+02
    12      1.125E-03     -5.169E-08      1.469E+01      3.676E+02
    13      6.330E-04     -1.017E-07      1.469E+01      3.676E+02
    14      3.498E-04     -4.289E-07      1.469E+01      3.676E+02
    15      1.849E-04     -4.265E-07      1.469E+01      3.676E+02
    16      9.821E-05     -4.376E-08      1.469E+01      3.676E+02
    17      4.871E-05      3.356E-07      1.469E+01      3.676E+02
    18      3.304E-05     -4.528E-07      1.469E+01      3.676E+02
    19      1.214E-05     -7.810E-07      1.469E+01      3.676E+02
    20      6.672E-06      1.781E-07      1.469E+01      3.676E+02
    21      4.681E-06      2.758E-07      1.469E+01      3.676E+02
    22      2.861E-06     -3.890E-09      1.469E+01      3.676E+02
    23      3.099E-06      2.777E-07      1.469E+01      3.676E+02
    24      9.537E-07     -2.346E-08      1.469E+01      3.676E+02
    25      9.537E-07     -3.388E-07      1.469E+01      3.676E+02
    26      7.153E-07      9.528E-08      1.469E+01      3.676E+02
    27      1.156E-06      2.714E-07      1.469E+01      3.676E+02
    28      1.455E-06     -2.555E-08      1.469E+01      3.676E+02
    29      1.002E-06      6.653E-08      1.469E+01      3.676E+02
    30      2.742E-06     -7.924E-08      1.469E+01      3.676E+02
  
  
  I =     1        2        3        4        5        6        7
  X = 0.00E+00 5.00E-02 1.50E-01 2.50E-01 3.50E-01 4.50E-01 5.00E-01
  
  
  J =     1        2        3        4        5        6        7
  Y = 0.00E+00 1.00E-01 3.00E-01 5.00E-01 7.00E-01 9.00E-01 1.10E+00
  
  J =     8        9       10       11       12
  Y = 1.30E+00 1.50E+00 1.70E+00 1.90E+00 2.00E+00
  
 **************************       VEL U    **************************
  
  I =     2        3        4        5        6        7
  J
  12   0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
  11   0.00E+00 1.41E-02 3.39E-02 4.04E-02 2.71E-02 0.00E+00
  10   0.00E+00-6.73E-02-1.96E-01-2.78E-01-2.11E-01 0.00E+00
   9   0.00E+00-1.55E-01-4.33E-01-5.97E-01-4.48E-01 0.00E+00
   8   0.00E+00-3.26E-01-8.75E-01-1.19E+00-8.95E-01 0.00E+00
   7   0.00E+00-6.17E-01-1.61E+00-2.16E+00-1.65E+00 0.00E+00
   6   0.00E+00-1.03E+00-2.62E+00-3.53E+00-2.75E+00 0.00E+00
   5   0.00E+00-1.42E+00-3.67E+00-5.06E+00-4.10E+00 0.00E+00
   4   0.00E+00-1.35E+00-3.91E+00-6.02E+00-5.42E+00 0.00E+00
   3   0.00E+00 1.37E-01-1.24E+00-4.69E+00-6.33E+00 0.00E+00
   2   0.00E+00 2.64E+00 6.16E+00 6.03E+00-7.70E+00 0.00E+00
   1   0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
  
 **************************       VEL V    **************************
  
  I =     1        2        3        4        5        6        7
  J
  12   0.00E+00 3.73E+00 9.97E+00 1.50E+01 1.87E+01 2.07E+01 1.20E+01
  11   0.00E+00 3.76E+00 1.01E+01 1.52E+01 1.89E+01 2.09E+01 1.20E+01
  10   0.00E+00 3.65E+00 9.94E+00 1.53E+01 1.95E+01 2.19E+01 1.20E+01
   9   0.00E+00 3.37E+00 9.57E+00 1.54E+01 2.04E+01 2.35E+01 1.20E+01
   8   0.00E+00 2.76E+00 8.70E+00 1.52E+01 2.17E+01 2.61E+01 1.20E+01
   7   0.00E+00 1.59E+00 7.02E+00 1.47E+01 2.35E+01 3.03E+01 1.20E+01
   6   0.00E+00-3.65E-01 4.21E+00 1.36E+01 2.60E+01 3.70E+01 1.20E+01
   5   0.00E+00-3.06E+00 1.81E-01 1.15E+01 2.89E+01 4.66E+01 1.20E+01
   4   0.00E+00-5.60E+00-4.41E+00 8.01E+00 3.09E+01 5.93E+01 1.20E+01
   3   0.00E+00-5.24E+00-6.77E+00 1.43E+00 2.77E+01 7.51E+01 1.20E+01
   2   0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 1.00E+02 0.00E+00
  
 **************************      STR FN    **************************
  
  I =     2        3        4        5        6        7
  J
  12   0.00E+00-3.63E-01-1.29E+00-2.63E+00-4.24E+00-6.00E+00
  11   0.00E+00-3.66E-01-1.29E+00-2.63E+00-4.25E+00-6.00E+00
  10   0.00E+00-3.53E-01-1.26E+00-2.58E+00-4.21E+00-6.00E+00
   9   0.00E+00-3.24E-01-1.18E+00-2.48E+00-4.14E+00-6.00E+00
   8   0.00E+00-2.64E-01-1.03E+00-2.29E+00-4.00E+00-6.00E+00
   7   0.00E+00-1.51E-01-7.61E-01-1.95E+00-3.74E+00-6.00E+00
   6   0.00E+00 3.46E-02-3.26E-01-1.40E+00-3.34E+00-6.00E+00
   5   0.00E+00 2.89E-01 2.74E-01-6.28E-01-2.74E+00-6.00E+00
   4   0.00E+00 5.31E-01 9.10E-01 2.79E-01-1.97E+00-6.00E+00
   3   0.00E+00 5.06E-01 1.12E+00 9.96E-01-1.09E+00-6.00E+00
   2   0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00-6.00E+00
  
 **************************      TEMP.     **************************
  
  I =     1        2        3        4        5        6        7
  J
  12   3.00E+02 3.08E+02 3.23E+02 3.37E+02 3.48E+02 3.53E+02 3.00E+02
  11   3.00E+02 3.08E+02 3.23E+02 3.37E+02 3.48E+02 3.53E+02 3.00E+02
  10   3.00E+02 3.09E+02 3.27E+02 3.43E+02 3.55E+02 3.62E+02 3.00E+02
   9   3.00E+02 3.11E+02 3.31E+02 3.50E+02 3.65E+02 3.73E+02 3.00E+02
   8   3.00E+02 3.12E+02 3.37E+02 3.59E+02 3.75E+02 3.84E+02 3.00E+02
   7   3.00E+02 3.14E+02 3.43E+02 3.68E+02 3.87E+02 3.97E+02 3.00E+02
   6   3.00E+02 3.16E+02 3.48E+02 3.76E+02 3.98E+02 4.10E+02 3.00E+02
   5   3.00E+02 3.18E+02 3.53E+02 3.83E+02 4.07E+02 4.23E+02 3.00E+02
   4   3.00E+02 3.18E+02 3.53E+02 3.85E+02 4.12E+02 4.35E+02 3.00E+02
   3   3.00E+02 3.15E+02 3.45E+02 3.76E+02 4.10E+02 4.49E+02 3.00E+02
   2   3.00E+02 3.06E+02 3.21E+02 3.42E+02 3.88E+02 4.69E+02 3.00E+02
   1   3.00E+02 3.00E+02 3.00E+02 3.00E+02 3.00E+02 5.00E+02 3.00E+02
  
 **************************     PRESSURE   **************************
  
  I =     1        2        3        4        5        6        7
  J
  12   8.40E+02 8.40E+02 8.39E+02 8.38E+02 8.34E+02 8.31E+02 8.30E+02
  11   8.52E+02 8.52E+02 8.52E+02 8.50E+02 8.48E+02 8.45E+02 8.44E+02
  10   8.77E+02 8.77E+02 8.76E+02 8.76E+02 8.75E+02 8.74E+02 8.73E+02
   9   8.99E+02 8.98E+02 8.97E+02 8.95E+02 8.94E+02 8.92E+02 8.91E+02
   8   9.12E+02 9.10E+02 9.08E+02 9.06E+02 9.05E+02 9.02E+02 9.00E+02
   7   9.06E+02 9.04E+02 9.01E+02 8.99E+02 8.99E+02 8.96E+02 8.94E+02
   6   8.63E+02 8.61E+02 8.56E+02 8.56E+02 8.62E+02 8.59E+02 8.58E+02
   5   7.55E+02 7.52E+02 7.46E+02 7.50E+02 7.66E+02 7.69E+02 7.70E+02
   4   5.57E+02 5.53E+02 5.45E+02 5.50E+02 5.85E+02 6.02E+02 6.11E+02
   3   2.91E+02 2.84E+02 2.72E+02 2.55E+02 3.32E+02 3.56E+02 3.68E+02
   2   9.85E+01 8.74E+01 6.54E+01-3.27E+01-2.08E+02 9.08E+01 2.40E+02
   1   0.00E+00-1.10E+01-3.79E+01-1.77E+02-4.78E+02-4.18E+01 1.07E+02
  
 **************************     DENSITY    **************************
  
  I =     1        2        3        4        5        6        7
  J
  12   1.00E+00 9.74E-01 9.28E-01 8.90E-01 8.63E-01 8.49E-01 1.00E+00
  11   1.00E+00 9.74E-01 9.28E-01 8.90E-01 8.63E-01 8.49E-01 1.00E+00
  10   1.00E+00 9.71E-01 9.18E-01 8.75E-01 8.44E-01 8.29E-01 1.00E+00
   9   1.00E+00 9.66E-01 9.05E-01 8.56E-01 8.23E-01 8.05E-01 1.00E+00
   8   1.00E+00 9.60E-01 8.91E-01 8.37E-01 7.99E-01 7.80E-01 1.00E+00
   7   1.00E+00 9.54E-01 8.76E-01 8.16E-01 7.76E-01 7.55E-01 1.00E+00
   6   1.00E+00 9.48E-01 8.62E-01 7.97E-01 7.54E-01 7.32E-01 1.00E+00
   5   1.00E+00 9.44E-01 8.51E-01 7.83E-01 7.37E-01 7.10E-01 1.00E+00
   4   1.00E+00 9.43E-01 8.49E-01 7.79E-01 7.27E-01 6.90E-01 1.00E+00
   3   1.00E+00 9.52E-01 8.68E-01 7.97E-01 7.32E-01 6.68E-01 1.00E+00
   2   1.00E+00 9.79E-01 9.35E-01 8.77E-01 7.73E-01 6.40E-01 1.00E+00
   1   1.00E+00 1.00E+00 1.00E+00 1.00E+00 1.00E+00 6.00E-01 1.00E+00