result8

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

               COMPUTATION FOR AXISYMMETRIC SITUATION
              ****************************************


   ITER       SMAX           SSUM          U(4,4)         V(4,4)
     0     0.0000E+00     0.0000E+00     0.0000E+00     0.0000E+00
     1     3.7273E-01     0.0000E+00     4.4573E+00     7.6537E+00
     2     5.7905E-02    -2.7940E-08     2.4668E+00     9.4167E+00
     3     5.6399E-02     2.6077E-08     1.9613E+00     7.6476E+00
     4     3.4764E-02    -1.5367E-08     1.9209E+00     5.6635E+00
     5     1.5171E-02    -4.2492E-08     1.9350E+00     4.5638E+00
     6     5.0682E-03    -3.7617E-08     1.9302E+00     4.1275E+00
     7     1.3521E-03    -1.3410E-08     1.9164E+00     3.9963E+00
     8     5.1792E-04    -6.1473E-08     1.9051E+00     3.9703E+00
     9     2.0999E-04     6.4450E-08     1.8986E+00     3.9708E+00
    10     7.0803E-05     2.6450E-08     1.8956E+00     3.9745E+00
    11     2.7210E-05     2.1607E-08     1.8945E+00     3.9768E+00
    12     1.0341E-05     6.7428E-08     1.8941E+00     3.9777E+00
    13     3.2187E-06    -3.0920E-08     1.8940E+00     3.9781E+00
    14     8.0466E-07     4.8056E-08     1.8939E+00     3.9782E+00
    15     3.8743E-07     1.7881E-08     1.8939E+00     3.9782E+00
    16     2.0862E-07    -1.0431E-08     1.8939E+00     3.9782E+00
    17     9.6858E-08     1.8626E-09     1.8939E+00     3.9782E+00
    18     8.9407E-08     6.3702E-08     1.8939E+00     3.9782E+00
    19     8.9407E-08    -4.1723E-08     1.8939E+00     3.9782E+00
    20     9.8325E-08    -4.5449E-08     1.8939E+00     3.9782E+00
    21     6.3330E-08    -2.4959E-08     1.8939E+00     3.9782E+00
    22     5.9605E-08    -3.9861E-08     1.8939E+00     3.9782E+00
    23     5.9605E-08    -6.0350E-08     1.8939E+00     3.9782E+00
    24     4.4703E-08    -1.6019E-08     1.8939E+00     3.9782E+00
    25     5.9605E-08    -3.4273E-08     1.8939E+00     3.9782E+00
  
  
  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 5.00E-02 1.50E-01 2.50E-01 3.50E-01 4.50E-01 5.00E-01
  
 **************************      VEL U     **************************
  
  I =     2        3        4        5        6        7
  J
   7   0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
   6   0.00E+00-1.33E+00-2.67E+00-2.12E+00-8.36E-01 0.00E+00
   5   0.00E+00-1.86E+00-2.70E+00-1.86E+00-6.39E-01 0.00E+00
   4   0.00E+00-2.18E-01 1.89E+00 2.90E+00 1.65E+00 0.00E+00
   3   0.00E+00 1.33E+01 1.97E+01 1.92E+01 1.04E+01 0.00E+00
   2   1.00E+02 8.63E+01 7.43E+01 5.99E+01 3.27E+01 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
   7   0.00E+00 0.00E+00 0.00E+00 0.00E+00 3.99E+00 6.01E+00 0.00E+00
   6   0.00E+00-1.50E+00-1.50E+00 6.18E-01 6.43E+00 8.45E+00 0.00E+00
   5   0.00E+00-4.17E+00-2.98E+00 1.81E+00 1.00E+01 1.20E+01 0.00E+00
   4   0.00E+00-6.53E+00-1.84E+00 3.98E+00 1.34E+01 1.60E+01 0.00E+00
   3   0.00E+00 6.87E+00 5.96E+00 7.22E+00 1.36E+01 1.63E+01 0.00E+00
   2   0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
  
 **************************      STR FN    **************************
  
  I =     2        3        4        5        6        7
  J
   7   5.00E-01 5.00E-01 5.00E-01 5.00E-01 3.00E-01 0.00E+00
   6   5.00E-01 5.60E-01 6.20E-01 5.95E-01 3.38E-01 0.00E+00
   5   5.00E-01 6.25E-01 7.15E-01 6.60E-01 3.60E-01 0.00E+00
   4   5.00E-01 6.31E-01 6.67E-01 5.88E-01 3.19E-01 0.00E+00
   3   5.00E-01 4.31E-01 3.72E-01 2.99E-01 1.63E-01 0.00E+00
   2   0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
  
 **************************      R.VTH     **************************
  
  I =     1        2        3        4        5        6        7
  J
   7   0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 2.50E+01
   6   0.00E+00 1.24E-01 5.24E-01 1.64E+00 5.76E+00 1.26E+01 2.02E+01
   5   0.00E+00 2.02E-01 7.28E-01 1.69E+00 3.66E+00 7.75E+00 1.23E+01
   4   0.00E+00 1.40E-01 4.46E-01 8.49E-01 1.53E+00 3.54E+00 6.25E+00
   3   0.00E+00 5.15E-02 1.49E-01 2.47E-01 3.84E-01 1.09E+00 2.25E+00
   2   0.00E+00 4.66E-03 1.84E-02 3.72E-02 5.53E-02 1.55E-01 2.50E-01
   1   0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
  
 **************************     PRESSURE   **************************
  
  I =     1        2        3        4        5        6        7
  J
   7  -4.93E+02-4.81E+02-4.56E+02-3.68E+02-3.47E+02-3.61E+02-3.61E+02
   6  -5.08E+02-4.96E+02-4.72E+02-3.94E+02-3.61E+02-3.61E+02-3.61E+02
   5  -5.38E+02-5.26E+02-5.02E+02-4.46E+02-3.89E+02-3.61E+02-3.47E+02
   4  -6.85E+02-6.47E+02-5.72E+02-4.92E+02-3.60E+02-2.41E+02-1.81E+02
   3  -1.15E+03-9.63E+02-5.97E+02-4.57E+02-1.85E+02 1.02E+02 2.46E+02
   2  -3.01E+02-3.62E+02-4.84E+02-3.04E+02 1.83E+02 6.20E+02 8.39E+02
   1   0.00E+00-6.11E+01-4.27E+02-2.28E+02 3.67E+02 8.79E+02 1.10E+03