From 8501a2320603942371c54b07029aac429bc69f37 Mon Sep 17 00:00:00 2001 From: Harry Moffat Date: Tue, 24 Mar 2009 19:36:09 +0000 Subject: [PATCH] Adding example dir. --- .../examples/flames/stflame1/.cvsignore | 19 + .../examples/flames/stflame1/Makefile.in | 15 + .../python/examples/flames/stflame1/cleanup | 5 + .../flames/stflame1/output_blessed_0.txt | 999 ++++++++++++++++++ .../examples/flames/stflame1/runtest.in | 64 ++ .../examples/flames/stflame1/stflame1.py | 117 ++ .../flames/stflame1/stflame1_1_blessed_0.csv | 79 ++ 7 files changed, 1298 insertions(+) create mode 100644 Cantera/python/examples/flames/stflame1/.cvsignore create mode 100644 Cantera/python/examples/flames/stflame1/Makefile.in create mode 100755 Cantera/python/examples/flames/stflame1/cleanup create mode 100644 Cantera/python/examples/flames/stflame1/output_blessed_0.txt create mode 100755 Cantera/python/examples/flames/stflame1/runtest.in create mode 100644 Cantera/python/examples/flames/stflame1/stflame1.py create mode 100644 Cantera/python/examples/flames/stflame1/stflame1_1_blessed_0.csv diff --git a/Cantera/python/examples/flames/stflame1/.cvsignore b/Cantera/python/examples/flames/stflame1/.cvsignore new file mode 100644 index 000000000..c6fb8d066 --- /dev/null +++ b/Cantera/python/examples/flames/stflame1/.cvsignore @@ -0,0 +1,19 @@ +Makefile +runtest +Makefile +ct2ctml.log +diff_csv.txt +diff_out_0.txt +h2o2.xml +output_0.txt +runtest +stflame1.xml +stflame1_1.csv +stflame1_2.csv +stflame1_3.csv +stflame1_4.csv +stflame1_5.csv +stflame1_6.csv +stflame1_7.csv +transport_log.xml + diff --git a/Cantera/python/examples/flames/stflame1/Makefile.in b/Cantera/python/examples/flames/stflame1/Makefile.in new file mode 100644 index 000000000..4a4e84482 --- /dev/null +++ b/Cantera/python/examples/flames/stflame1/Makefile.in @@ -0,0 +1,15 @@ +#!/bin/sh + +PYTHON_CMD = @PYTHON_CMD@ + +run: + $(PYTHON_CMD) stflame1.py + +test: + ./runtest +clean: + rm -f *.log *.csv *.xml + ./cleanup + +# end of file + diff --git a/Cantera/python/examples/flames/stflame1/cleanup b/Cantera/python/examples/flames/stflame1/cleanup new file mode 100755 index 000000000..55b97eb53 --- /dev/null +++ b/Cantera/python/examples/flames/stflame1/cleanup @@ -0,0 +1,5 @@ +#!/bin/sh +# +/bin/rm -rf equilibrate_log*.html +/bin/rm -rf .cttmp* ct2ctml.log transport_log.xml vcs_equilibrate_res*.csv \ + stflame1_[1234567].csv output_0.txt diff* diff --git a/Cantera/python/examples/flames/stflame1/output_blessed_0.txt b/Cantera/python/examples/flames/stflame1/output_blessed_0.txt new file mode 100644 index 000000000..e0537ec89 --- /dev/null +++ b/Cantera/python/examples/flames/stflame1/output_blessed_0.txt @@ -0,0 +1,999 @@ + + +>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> inlet <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< + + Mass Flux: 0.06 kg/m^2/s + Temperature: 373 K + Mass Fractions: + H2 0.01151 + O2 0.1015 + AR 0.887 + + + +>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> flow <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< + + Pressure: 5066 Pa + +------------------------------------------------------------------------------- + z u V T lambda H2 +------------------------------------------------------------------------------- + 0 1.142 0 373 0 0.01151 + 0.02 1.028 0 966.2 0 0.007795 + 0.04 0.9134 0 1559 0 0.00408 + 0.06 0.7992 0 2153 0 0.000365 + 0.08 0.685 0 2153 0 0.000365 + 0.1 0.5709 0 2153 0 0.000365 + 0.15 0.2854 0 1894 0 0.000365 + 0.2 0 0 600 0 0.000365 + +------------------------------------------------------------------------------- + z H O O2 OH H2O +------------------------------------------------------------------------------- + 0 0 0 0.1015 0 0 + 0.02 1.461e-05 0.000152 0.07153 0.000956 0.03256 + 0.04 2.922e-05 0.000304 0.04156 0.001912 0.06512 + 0.06 4.383e-05 0.000456 0.01159 0.002868 0.09769 + 0.08 4.383e-05 0.000456 0.01159 0.002868 0.09769 + 0.1 4.383e-05 0.000456 0.01159 0.002868 0.09769 + 0.15 4.383e-05 0.000456 0.01159 0.002868 0.09769 + 0.2 4.383e-05 0.000456 0.01159 0.002868 0.09769 + +------------------------------------------------------------------------------- + z HO2 H2O2 AR +------------------------------------------------------------------------------- + 0 0 0 0.887 + 0.02 1.259e-07 4.647e-09 0.887 + 0.04 2.517e-07 9.293e-09 0.887 + 0.06 3.776e-07 1.394e-08 0.887 + 0.08 3.776e-07 1.394e-08 0.887 + 0.1 3.776e-07 1.394e-08 0.887 + 0.15 3.776e-07 1.394e-08 0.887 + 0.2 3.776e-07 1.394e-08 0.887 + + +>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> surface <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< + + Temperature: 600 K + +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 1 timesteps 1.5e-05 6.261 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 2 timesteps 1.688e-05 2.914 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 5 timesteps 3.204e-05 2.909 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.0001155 2.876 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.0004161 2.753 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [8] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 0 1 2 + to resolve H H2 H2O H2O2 HO2 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [11] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 0 + to resolve HO2 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [12] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 0 1 + to resolve H2O2 HO2 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [14] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 0 + to resolve HO2 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [15] point grid(s). + +no new points needed in flow +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 1 timesteps 5.859e-08 5.265 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 2 timesteps 1.006e-12 5.265 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 5 timesteps 3.819e-12 5.265 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 5.506e-11 5.265 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 7.937e-10 5.265 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 1.144e-08 5.265 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 1.65e-07 5.264 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 2.378e-06 5.254 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 3.428e-05 5.239 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.0002471 5.143 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.0008906 4.716 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.00428 4.059 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.00457 3.871 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.003856 2.499 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [15] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 0 1 2 3 4 5 6 7 8 10 11 12 13 + to resolve H H2 H2O H2O2 HO2 O O2 OH T V u +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [28] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 0 1 2 3 4 5 6 7 8 9 10 11 23 24 25 26 + to resolve H H2 H2O H2O2 HO2 O O2 OH T V u +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [44] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 0 1 2 3 4 5 6 8 20 38 39 40 41 42 + to resolve H2O H2O2 HO2 O O2 OH T V u +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [58] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 0 1 28 52 53 54 55 56 + to resolve H2O2 HO2 O OH T V +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [66] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 0 59 60 61 62 63 64 + to resolve HO2 OH T V +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [73] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 68 69 70 71 + to resolve T V +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [77] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 75 + to resolve V +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [78] point grid(s). + +no new points needed in flow +Solution saved to file stflame1.xml as solution mdot1_1. +solution saved to flame1.csv +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [78] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 20 21 22 23 24 33 34 35 74 + to resolve H H2 H2O H2O2 O O2 OH V +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [87] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 24 40 41 + to resolve H H2O2 O O2 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [90] point grid(s). + +no new points needed in flow +Solution saved to file stflame1.xml as solution mdot2_1. +solution saved to flame1.csv +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 1 timesteps 1.5e-05 5.337 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 2 timesteps 8.438e-06 4.411 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 5 timesteps 3.204e-05 4.402 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.0004618 4.315 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.003329 4.156 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.003999 4.088 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [90] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 31 32 33 39 40 41 42 45 46 47 48 49 50 + to resolve H H2 H2O H2O2 HO2 O O2 OH T u +refine: discarding point at 3.90625e-05 +refine: discarding point at 7.8125e-05 +refine: discarding point at 0.00015625 +refine: discarding point at 0.000234375 +refine: discarding point at 0.0003125 +refine: discarding point at 0.00046875 +refine: discarding point at 0.000625 +refine: discarding point at 0.00078125 +refine: discarding point at 0.0009375 +refine: discarding point at 0.00109375 +refine: discarding point at 0.00125 +refine: discarding point at 0.00140625 +refine: discarding point at 0.0015625 +refine: discarding point at 0.00171875 +refine: discarding point at 0.001875 +refine: discarding point at 0.00203125 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [87] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 0 21 22 23 24 34 35 36 40 41 42 43 44 45 46 + to resolve H H2 H2O H2O2 HO2 O O2 point 0 u +refine: discarding point at 0.0025 +refine: discarding point at 0.0028125 +refine: discarding point at 0.003125 +refine: discarding point at 0.0034375 +refine: discarding point at 0.00375 +refine: discarding point at 0.00453125 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [96] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 2 41 42 43 44 46 47 48 49 51 52 53 54 55 56 57 60 + to resolve H H2 H2O H2O2 HO2 O O2 OH point 2 +refine: discarding point at 0.004375 +refine: discarding point at 0.0046875 +refine: discarding point at 0.005 +refine: discarding point at 0.0053125 +refine: discarding point at 0.005625 +refine: discarding point at 0.0059375 +refine: discarding point at 0.00625 +refine: discarding point at 0.0065625 +refine: discarding point at 0.006875 +refine: discarding point at 0.0071875 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [103] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 4 40 41 47 48 49 53 54 55 56 57 60 61 62 63 64 65 + to resolve H H2 H2O H2O2 HO2 O O2 OH T point 4 u +refine: discarding point at 0.003125 +refine: discarding point at 0.0078125 +refine: discarding point at 0.008125 +refine: discarding point at 0.0084375 +refine: discarding point at 0.00875 +refine: discarding point at 0.0090625 +refine: discarding point at 0.009375 +refine: discarding point at 0.0096875 +refine: discarding point at 0.01 +refine: discarding point at 0.0103125 +refine: discarding point at 0.010625 +refine: discarding point at 0.01125 +refine: discarding point at 0.0115625 +refine: discarding point at 0.011875 +refine: discarding point at 0.0121875 +refine: discarding point at 0.0125 +refine: discarding point at 0.0128125 +refine: discarding point at 0.013125 +refine: discarding point at 0.0134375 +refine: discarding point at 0.01375 +refine: discarding point at 0.0140625 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [99] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 5 24 25 26 27 28 36 37 48 49 50 51 54 55 56 57 58 59 64 + to resolve H H2 H2O H2O2 HO2 O O2 OH point 5 +refine: discarding point at 0.0146875 +refine: discarding point at 0.015 +refine: discarding point at 0.0153125 +refine: discarding point at 0.015625 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [114] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 60 61 70 71 72 73 74 + to resolve H H2 H2O H2O2 O O2 +refine: discarding point at 0.016875 +refine: discarding point at 0.0171875 +refine: discarding point at 0.0175 +refine: discarding point at 0.0178125 +refine: discarding point at 0.018125 +refine: discarding point at 0.0184375 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [115] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 8 + to resolve point 8 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [116] point grid(s). + +no new points needed in flow +Solution saved to file stflame1.xml as solution mdot3_1. +solution saved to flame1.csv +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 1 timesteps 1.5e-05 5.331 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 2 timesteps 1.688e-05 4.34 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [116] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 59 60 61 62 68 69 70 74 75 76 77 78 79 80 82 + to resolve H H2 H2O H2O2 HO2 O O2 OH T u +refine: discarding point at 0.0190625 +refine: discarding point at 0.019375 +refine: discarding point at 0.0196875 +refine: discarding point at 0.02 +refine: discarding point at 0.0203125 +refine: discarding point at 0.020625 +refine: discarding point at 0.0209375 +refine: discarding point at 0.02125 +refine: discarding point at 0.0215625 +refine: discarding point at 0.021875 +refine: discarding point at 0.0221875 +refine: discarding point at 0.0225 +refine: discarding point at 0.0228125 +refine: discarding point at 0.023125 +refine: discarding point at 0.0234375 +refine: discarding point at 0.02375 +refine: discarding point at 0.0240625 +refine: discarding point at 0.024375 +refine: discarding point at 0.0246875 +refine: discarding point at 0.025 +refine: discarding point at 0.025625 +refine: discarding point at 0.02625 +refine: discarding point at 0.026875 +refine: discarding point at 0.0271875 +refine: discarding point at 0.0275 +refine: discarding point at 0.0278125 +refine: discarding point at 0.028125 +refine: discarding point at 0.0284375 +refine: discarding point at 0.02875 +refine: discarding point at 0.0290625 +refine: discarding point at 0.029375 +refine: discarding point at 0.0296875 +refine: discarding point at 0.03 +refine: discarding point at 0.0303125 +refine: discarding point at 0.030625 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [96] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 10 32 43 44 48 49 50 51 52 54 55 56 57 + to resolve H H2 H2O H2O2 HO2 O O2 OH T point 10 u +refine: discarding point at 0.03125 +refine: discarding point at 0.0315625 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [107] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 11 33 34 47 58 59 60 62 63 64 65 66 67 + to resolve H H2 H2O H2O2 HO2 O O2 point 11 +refine: discarding point at 0.031875 +refine: discarding point at 0.0321875 +refine: discarding point at 0.0325 +refine: discarding point at 0.0328125 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [116] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 13 34 35 62 63 64 77 78 + to resolve H2O2 HO2 O point 13 +refine: discarding point at 0.0334375 +refine: discarding point at 0.03375 +refine: discarding point at 0.0340625 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [121] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 66 67 68 82 85 + to resolve H2O2 O OH +refine: discarding point at 0.0346875 +refine: discarding point at 0.035 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [124] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 70 + to resolve H2O2 +refine: discarding point at 0.035625 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [124] point grid(s). + +no new points needed in flow +Solution saved to file stflame1.xml as solution mdot4_1. +solution saved to flame1.csv +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [124] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 73 74 75 76 79 80 81 82 83 86 87 88 89 90 91 92 + to resolve H H2 H2O H2O2 O O2 OH T u +refine: discarding point at 0.0353125 +refine: discarding point at 0.0359375 +refine: discarding point at 0.03625 +refine: discarding point at 0.0365625 +refine: discarding point at 0.036875 +refine: discarding point at 0.0371875 +refine: discarding point at 0.0375 +refine: discarding point at 0.0378125 +refine: discarding point at 0.038125 +refine: discarding point at 0.0384375 +refine: discarding point at 0.03875 +refine: discarding point at 0.0390625 +refine: discarding point at 0.039375 +refine: discarding point at 0.0396875 +refine: discarding point at 0.04 +refine: discarding point at 0.0403125 +refine: discarding point at 0.040625 +refine: discarding point at 0.0409375 +refine: discarding point at 0.04125 +refine: discarding point at 0.0415625 +refine: discarding point at 0.041875 +refine: discarding point at 0.0421875 +refine: discarding point at 0.0425 +refine: discarding point at 0.0428125 +refine: discarding point at 0.043125 +refine: discarding point at 0.0434375 +refine: discarding point at 0.04375 +refine: discarding point at 0.0440625 +refine: discarding point at 0.044375 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [111] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 16 52 53 64 65 66 67 68 69 70 71 72 73 + to resolve H H2 H2O H2O2 O O2 point 16 +refine: discarding point at 0.045 +refine: discarding point at 0.0453125 +refine: discarding point at 0.045625 +refine: discarding point at 0.0459375 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [120] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 39 40 65 66 67 68 69 79 80 81 + to resolve H H2O H2O2 HO2 O O2 +refine: discarding point at 0.0465625 +refine: discarding point at 0.046875 +refine: discarding point at 0.0471875 +refine: discarding point at 0.0475 +refine: discarding point at 0.0478125 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [125] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 19 72 + to resolve H2O2 point 19 +refine: discarding point at 0.0484375 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [126] point grid(s). + +no new points needed in flow +Solution saved to file stflame1.xml as solution mdot5_1. +solution saved to flame1.csv +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [126] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 75 76 77 78 82 83 84 85 86 88 89 90 91 92 95 + to resolve H H2 H2O H2O2 O O2 OH T u +refine: discarding point at 0.00109375 +refine: discarding point at 0.0320312 +refine: discarding point at 0.0471875 +refine: discarding point at 0.048125 +refine: discarding point at 0.04875 +refine: discarding point at 0.0490625 +refine: discarding point at 0.049375 +refine: discarding point at 0.0496875 +refine: discarding point at 0.05 +refine: discarding point at 0.0503125 +refine: discarding point at 0.050625 +refine: discarding point at 0.0509375 +refine: discarding point at 0.05125 +refine: discarding point at 0.0515625 +refine: discarding point at 0.051875 +refine: discarding point at 0.0521875 +refine: discarding point at 0.0525 +refine: discarding point at 0.0528125 +refine: discarding point at 0.053125 +refine: discarding point at 0.0534375 +refine: discarding point at 0.05375 +refine: discarding point at 0.0540625 +refine: discarding point at 0.054375 +refine: discarding point at 0.0546875 +refine: discarding point at 0.055 +refine: discarding point at 0.0553125 +refine: discarding point at 0.055625 +refine: discarding point at 0.0559375 +refine: discarding point at 0.05625 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [112] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 17 54 55 56 67 68 69 70 71 72 73 74 75 76 + to resolve H H2 H2O H2O2 O O2 OH point 17 +refine: discarding point at 0.056875 +refine: discarding point at 0.0571875 +refine: discarding point at 0.0575 +refine: discarding point at 0.0578125 +refine: discarding point at 0.058125 +refine: discarding point at 0.0584375 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [120] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 19 40 71 79 80 81 82 + to resolve H H2 H2O H2O2 HO2 O O2 point 19 +refine: discarding point at 0.0590625 +refine: discarding point at 0.059375 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [125] point grid(s). + +no new points needed in flow +Solution saved to file stflame1.xml as solution mdot6_1. +solution saved to flame1.csv +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 1 timesteps 7.5e-06 5.618 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 2 timesteps 8.438e-06 4.225 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 5 timesteps 1.602e-05 4.219 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.0001155 4.174 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.0008323 4.037 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.005999 3.803 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.005405 3.511 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.00974 3.018 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.0234 2.641 +.............................................................................. + +Attempt Newton solution of steady-state problem... failure. + +.............................................................................. +Take 10 timesteps 0.01111 2.98 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [125] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 78 79 80 81 84 85 86 87 88 90 91 92 93 94 95 96 + to resolve H H2 H2O H2O2 HO2 O O2 OH T u +refine: discarding point at 0.0576562 +refine: discarding point at 0.05875 +refine: discarding point at 0.0596875 +refine: discarding point at 0.06 +refine: discarding point at 0.0603125 +refine: discarding point at 0.060625 +refine: discarding point at 0.0609375 +refine: discarding point at 0.06125 +refine: discarding point at 0.0615625 +refine: discarding point at 0.061875 +refine: discarding point at 0.0621875 +refine: discarding point at 0.0625 +refine: discarding point at 0.0628125 +refine: discarding point at 0.063125 +refine: discarding point at 0.0634375 +refine: discarding point at 0.06375 +refine: discarding point at 0.0640625 +refine: discarding point at 0.064375 +refine: discarding point at 0.0646875 +refine: discarding point at 0.065 +refine: discarding point at 0.0653125 +refine: discarding point at 0.065625 +refine: discarding point at 0.0659375 +refine: discarding point at 0.06625 +refine: discarding point at 0.0665625 +refine: discarding point at 0.066875 +refine: discarding point at 0.0671875 +refine: discarding point at 0.0675 +refine: discarding point at 0.0678125 +refine: discarding point at 0.068125 +refine: discarding point at 0.0684375 +refine: discarding point at 0.06875 +refine: discarding point at 0.0690625 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [108] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 19 43 44 53 54 65 66 67 68 69 + to resolve H H2 H2O H2O2 HO2 O O2 OH point 19 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [118] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 20 76 77 78 + to resolve H H2O HO2 O O2 point 20 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [122] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 21 44 + to resolve H2O2 point 21 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [124] point grid(s). + +############################################################################## +Refining grid in flow. + New points inserted after grid points 22 + to resolve point 22 +.............................................................................. + +Attempt Newton solution of steady-state problem... success. + +Problem solved on [125] point grid(s). + +no new points needed in flow +Solution saved to file stflame1.xml as solution mdot7_1. +solution saved to flame1.csv + +Statistics: + + Grid Functions Time Jacobians Time + 10 401 0.3200 21 0.2800 + 13 16 0.0300 2 0.0400 + 14 14 0.0100 2 0.0500 + 16 14 0.0200 2 0.0600 + 17 10 0.0200 1 0.0300 + 17 4750 6.3400 129 3.1300 + 30 35 0.1400 2 0.1300 + 46 23 0.1200 2 0.2300 + 60 14 0.1100 2 0.3000 + 68 14 0.1400 2 0.3400 + 75 8 0.0800 1 0.2000 + 79 4 0.0500 1 0.2100 + 80 4 0.0400 1 0.1500 + 80 53 0.3800 5 0.7000 + 89 10 0.0700 1 0.1600 + 92 12 0.0900 1 0.1600 + 92 843 7.2400 44 7.4700 + 89 25 0.2000 2 0.3100 + 98 46 0.4000 4 0.7000 + 105 60 0.6000 6 1.1600 + 101 57 0.5600 5 0.9200 + 116 40 0.4300 4 0.8900 + 117 10 0.1100 1 0.2200 + 118 2 0.0200 1 0.2200 + 118 162 1.8200 21 5.2800 + 98 18 0.1800 2 0.3500 + 109 22 0.2200 2 0.4100 + 118 30 0.3100 3 0.6600 + 123 30 0.3500 3 0.7000 + 126 26 0.2900 3 0.7200 + 126 18 0.2200 2 0.4800 + 126 59 0.9200 8 2.2300 + 113 28 0.3900 3 0.8500 + 122 41 0.4700 4 0.9200 + 127 22 0.2400 2 0.4800 + 128 18 0.2200 2 0.4900 + 128 59 0.7000 7 1.7200 + 114 39 0.4100 4 0.8500 + 122 26 0.2900 3 0.6900 + 127 14 0.1500 2 0.5000 + 127 2580 29.7300 120 29.4600 + 110 53 0.5300 4 0.8200 + 120 18 0.1900 2 0.4500 + 124 14 0.1500 2 0.4700 + 126 10 0.1200 1 0.2500 + 127 2 0.0200 1 0.2500 diff --git a/Cantera/python/examples/flames/stflame1/runtest.in b/Cantera/python/examples/flames/stflame1/runtest.in new file mode 100755 index 000000000..457d2dff1 --- /dev/null +++ b/Cantera/python/examples/flames/stflame1/runtest.in @@ -0,0 +1,64 @@ +#!/bin/sh +# +# +temp_success="1" +/bin/rm -f output_0.txt stflame1_[1234567].csv diff_csv.txt diff_out_0.txt + +########################################################################## +PYTHON_CMD=@PYTHON_CMD@ +prog=stflame1.py +if test ! -f $prog ; then + echo $prog ' does not exist' + exit -1 +fi +################################################################# +# +CANTERA_DATA=${CANTERA_DATA:=../../../data/inputs}; export CANTERA_DATA +CANTERA_BIN=${CANTERA_BIN:=../../../bin} + +################################################################# + +$PYTHON_CMD $prog > output_0.txt <<+ +1.0 ++ +retnStat=$? +if [ $retnStat != "0" ] +then + temp_success="0" + echo "$prog returned with bad status, $retnStat, check output" +fi + +diff -w output_blessed_0.txt output_0.txt > diff_out_0.txt +retnStat_0=$? + +csvdiff -a 1.0E-50 stflame1_1_blessed_0.csv stflame1_1.csv > diff_csv.txt +retnStat_csv_0=$? + +retnTotal=1 +if test $retnStat_0 = "0" +then + retnTotal=0 +fi + +retnCSVTotal=1 +if test $retnStat_csv_0 = "1" +then + retnCSVTotal=0 +fi + +if test $retnCSVTotal = "0" +then + echo "Successful test comparison on "`pwd` + if test $retnTotal = "1" + then + echo " But text files show differences. See diff_out_0.txt" + fi +else + echo "Unsuccessful test comparison of csv files on "`pwd` " test" + echo " see diff_csv.txt " + if test $retnTotal != "0" + then + echo " ASCII files are different too - see diff_test*.txt" + fi +fi + diff --git a/Cantera/python/examples/flames/stflame1/stflame1.py b/Cantera/python/examples/flames/stflame1/stflame1.py new file mode 100644 index 000000000..88bd4fee2 --- /dev/null +++ b/Cantera/python/examples/flames/stflame1/stflame1.py @@ -0,0 +1,117 @@ +# +# STFLAME1 - A detached flat flame stabilized at a stagnation point +# + +# This script simulates a lean hydrogen-oxygen flame stabilized in +# a strained flowfield at an axisymmetric stagnation point on a +# non-reacting surface. The solution begins with a flame attached +# to the inlet (burner), and the mass flow rate is progressively +# increased, causing the flame to detach and move closer to the +# surface. This example illustrates use of the new 'prune' grid +# refinement parameter, which allows grid points to be removed if +# they are no longer required to resolve the solution. This is +# important here, since the flamefront moves as the mass flowrate +# is increased. Without using 'prune', a large number of grid +# points would be concentrated upsteam of the flame, where the +# flamefront had been previously. (To see this, try setting prune +# to zero.) + +from Cantera import * +from Cantera.OneD import * +from Cantera.OneD.StagnationFlow import StagnationFlow + +################################################################ +# +# parameter values +# +p = 0.05*OneAtm # pressure +tburner = 373.0 # burner temperature +tsurf = 600.0 + +# each mdot value will be solved to convergence, with grid refinement, +# and then that solution will be used for the next mdot +mdot = [0.06, 0.07, 0.08, 0.09, 0.1, 0.11, 0.12] # kg/m^2/s + +rxnmech = 'h2o2.cti' # reaction mechanism file +comp = 'H2:1.8, O2:1, AR:7' # premixed gas composition + +# The solution domain is chosen to be 50 cm, and a point very near the +# downstream boundary is added to help with the zero-gradient boundary +# condition at this boundary. +initial_grid = [0.0, 0.02, 0.04, 0.06, 0.08, 0.1, + 0.15, 0.2] # m + +tol_ss = [1.0e-5, 1.0e-13] # [rtol atol] for steady-state + # problem +tol_ts = [1.0e-4, 1.0e-9] # [rtol atol] for time stepping + +loglevel = 1 # amount of diagnostic output (0 + # to 5) + +refine_grid = 1 # 1 to enable refinement, 0 to + # disable +ratio = 5.0 +slope = 0.1 +curve = 0.2 +prune = 0.05 + + + +################ create the gas object ######################## +# +# This object will be used to evaluate all thermodynamic, kinetic, +# and transport properties +# +gas = IdealGasMix(rxnmech) + +# set its state to that of the unburned gas at the burner +gas.setState_TPX(tburner, p, comp) + +# Create the stagnation flow object with a non-reactive surface. (To +# make the surface reactive, supply a surface reaction mechanism. see +# example catcomb.py for how to do this.) +f = StagnationFlow(gas = gas, grid = initial_grid) + +# set the properties at the inlet +f.inlet.set(massflux = mdot[0], mole_fractions = comp, temperature = tburner) + +# set the surface state +f.surface.setTemperature(tsurf) + +f.set(tol = tol_ss, tol_time = tol_ts) +f.setMaxJacAge(5, 10) +f.set(energy = 'off') +f.init(products = 'equil') # assume adiabatic equilibrium products +f.showSolution() + +f.solve(loglevel, refine_grid) + +f.setRefineCriteria(ratio = ratio, slope = slope, + curve = curve, prune = prune) +f.set(energy = 'on') + +m = 0 +for md in mdot: + f.inlet.set(mdot = md) + f.solve(loglevel,refine_grid) + m = m + 1 + f.save('stflame1.xml','mdot'+`m`,'mdot = '+`md`+' kg/m2/s') + + + # write the velocity, temperature, and mole fractions to a CSV file + z = f.flow.grid() + T = f.T() + u = f.u() + V = f.V() + fcsv = open('stflame1_'+`m`+'.csv','w') + writeCSV(fcsv, ['z (m)', 'u (m/s)', 'V (1/s)', 'T (K)'] + + list(gas.speciesNames())) + for n in range(f.flow.nPoints()): + f.setGasState(n) + writeCSV(fcsv, [z[n], u[n], V[n], T[n]]+list(gas.moleFractions())) + fcsv.close() + + print 'solution saved to flame1.csv' + +f.showStats() + diff --git a/Cantera/python/examples/flames/stflame1/stflame1_1_blessed_0.csv b/Cantera/python/examples/flames/stflame1/stflame1_1_blessed_0.csv new file mode 100644 index 000000000..e29045f78 --- /dev/null +++ b/Cantera/python/examples/flames/stflame1/stflame1_1_blessed_0.csv @@ -0,0 +1,79 @@ +z (m), u (m/s), V (1/s), T (K), H2, H, O, O2, OH, H2O, HO2, H2O2, AR, +0.0, 1.1269138452534573, -1.1602530123599895e-18, 373.0, 0.16695382378294518, 0.0031594159146740831, 2.5749750841786712e-05, 0.10184082149446209, 0.00047317349057369264, 0.0035735664005197725, 1.697012788952471e-05, 9.7306818262049284e-06, 0.72394674835626749, +3.9062500000000001e-05, 1.1371219334494318, 0.0051978363943706724, 376.44446003546778, 0.16666059746237277, 0.0031933200289385876, 2.6797576976430301e-05, 0.10175171253386257, 0.00049282266137237784, 0.0037715912660259299, 1.8079695162514993e-05, 1.037184088742118e-05, 0.72407470693440135, +7.8125000000000002e-05, 1.1472568596118624, 0.010391665633537635, 379.8658962688084, 0.16636480638127599, 0.0032301164727841451, 2.7765942324678221e-05, 0.10166397263426888, 0.00050963971822362011, 0.0039721697412226538, 1.8418716079916611e-05, 1.0956317836311545e-05, 0.72420215407598387, +0.00015625, 1.1673010198488134, 0.020768200220209993, 386.63716433658993, 0.16576833577204059, 0.0033117656237137078, 2.9453637327520617e-05, 0.10149250221256256, 0.00053485392321793511, 0.0043784000693143045, 1.8096916918393308e-05, 1.1945581793614534e-05, 0.72445464626311151, +0.00023437499999999999, 1.1870556401592762, 0.031130846573563489, 393.31754748110751, 0.16516234581423489, 0.0034037169054290423, 3.085868075989289e-05, 0.10132396658386209, 0.00055107645284741935, 0.0047919652582060169, 1.7503268341397989e-05, 1.2707101360791019e-05, 0.72470585993495851, +0.00031250000000000001, 1.2065407315752454, 0.041480540894505202, 399.91360659446309, 0.16454722590756557, 0.0035051530228125734, 3.2027045942592047e-05, 0.10115747870235019, 0.00056022151250060225, 0.0052117376653196736, 1.6890616990866092e-05, 1.3251906181049915e-05, 0.7249560136203369, +0.00046874999999999998, 1.2447801361488315, 0.062146732436985701, 412.8759957637867, 0.1633007069323115, 0.0037323732654734264, 3.3781398600230934e-05, 0.10083132921230167, 0.00056216670269419982, 0.0060603514324214887, 1.5783891131805383e-05, 1.3744009753828619e-05, 0.72544976315531184, +0.00062500000000000001, 1.2822034733595999, 0.082773928548448103, 425.58729646738658, 0.16202181121934209, 0.0039897428458355892, 3.51114818680169e-05, 0.10050785781173724, 0.00055163300989821207, 0.006923794817287842, 1.4831023124096334e-05, 1.3656418091604356e-05, 0.7259415613728154, +0.00078125000000000004, 1.3189290864573651, 0.10336670825046008, 438.08660525944407, 0.16071183971991326, 0.0042741161281480361, 3.6256126722182145e-05, 0.10018519093843194, 0.00053392419679190457, 0.0077992574618779538, 1.3995449562465722e-05, 1.315106348625625e-05, 0.72643226891506607, +0.00093749999999999997, 1.3550519949336213, 0.12392878310180702, 450.40541185501814, 0.15937155246811427, 0.0045832025599421879, 3.7402367245515362e-05, 0.099862113649702114, 0.00051237409632229163, 0.0086851955787108654, 1.3250312512904096e-05, 1.2363377207843735e-05, 0.72692254559024183, +0.0010937500000000001, 1.3906491621445891, 0.14446319070039659, 462.56928594041131, 0.15800127843085149, 0.004915325820741878, 3.870541723425692e-05, 0.099537734536521053, 0.00048909529400325148, 0.0095810148201524916, 1.2578041898445439e-05, 1.1400153592753749e-05, 0.72741286748500433, +0.00125, 1.4257836469436014, 0.16497243614576237, 474.59921218758899, 0.15660100240395861, 0.0052692578161408756, 4.0305955873324369e-05, 0.099211285013011768, 0.00046542809094366388, 0.010486869081794787, 1.1966732763473041e-05, 1.0342907360713014e-05, 0.72790354199815277, +0.0014062499999999999, 1.4605077512516862, 0.18545859538804255, 486.5126039539013, 0.15517043473974171, 0.0056441017882332654, 4.2343893436038596e-05, 0.098881996452862136, 0.00044221473247492844, 0.011403527192190566, 1.1407838448741067e-05, 9.2524954597076296e-06, 0.72839472086715285, +0.0015625000000000001, 1.4948653750957082, 0.20592339098202836, 498.32406113397315, 0.15370906787826225, 0.0060392081003653479, 4.4969123438261168e-05, 0.098549022504770378, 0.00041997149898005242, 0.01233228142537926, 1.0894870152065874e-05, 8.1734286880518497e-06, 0.72888641116996422, +0.00171875, 1.5288937699844487, 0.22636824780704656, 510.04593255722614, 0.15221622245564551, 0.0064541123540089902, 4.8349934485048127e-05, 0.098211387512211742, 0.0003989996998082628, 0.013274883377285735, 1.0422664843561336e-05, 7.1374342572973523e-06, 0.7293784845674538, +0.0018749999999999999, 1.5626248438191674, 0.24679433411217191, 521.68873200347025, 0.15069108505805082, 0.0068884892026782642, 5.2679618539344149e-05, 0.097867949952274172, 0.00037945884522511684, 0.014233498886238396, 9.9869623094088169e-06, 6.1662696484876188e-06, 0.72987068520503595, +0.0020312500000000001, 1.5960861353083913, 0.267202591659643, 533.26144596188487, 0.14913273925655568, 0.0073421174822777373, 5.8181649144244276e-05, 0.0975173745259027, 0.00036141530955842801, 0.015210676976891241, 9.5841377410960237e-06, 5.2739262080305269e-06, 0.73036263673572088, +0.0021875000000000002, 1.6293015500889272, 0.28759375758954264, 544.77176314997155, 0.1475401912473219, 0.0078148536710352205, 6.5113685687838646e-05, 0.097158109365405476, 0.00034487419308606136, 0.016209329735635662, 9.2108359675111434e-06, 4.4683786325791908e-06, 0.73085384888722782, +0.0025000000000000001, 1.6950922130399506, 0.32833104578390071, 567.6283457517809, 0.1442839977695024, 0.0088066259006495725, 8.3919785909670629e-05, 0.096416255057566458, 0.00031646352376704065, 0.018258642112893603, 8.5415984272714952e-06, 3.1400744601070932e-06, 0.73182241417682392, +0.0028124999999999999, 1.7601622395233132, 0.36901366212278919, 590.31870313094339, 0.14088483344451352, 0.0098719750196863704, 0.00011163737997050631, 0.095623770953552775, 0.00029371184349899379, 0.020420866055861034, 7.9611546172885317e-06, 2.1320693115309563e-06, 0.73278311207898794, +0.0031250000000000002, 1.8246286300679793, 0.40964404501147378, 612.87705972675963, 0.13733474778187321, 0.011010645764518753, 0.00015114112299886392, 0.094763226121449565, 0.00027657370910097109, 0.022724582236932248, 7.4552291024611113e-06, 1.4002735722768633e-06, 0.7337302277604516, +0.0037499999999999999, 1.9521607688587745, 0.49077192568950184, 657.68309408614073, 0.12989962892314597, 0.013459705478942251, 0.00027251049278819247, 0.092815357633142387, 0.00025929489736134189, 0.02775188715726919, 6.615987752514378e-06, 6.0112368266695839e-07, 0.7355343983059156, +0.0043750000000000004, 2.0782697316183487, 0.57173796795836584, 702.2236384948236, 0.12184330358254464, 0.016177003587254928, 0.00046837284595409659, 0.090433025016398591, 0.00026685838871252378, 0.033568046781064134, 5.9527172183304049e-06, 2.5038304360501136e-07, 0.73723718669780902, +0.0050000000000000001, 2.2031577992780838, 0.6525386730449847, 746.49909985601039, 0.11317232017836959, 0.019138349511985477, 0.00075406185999996847, 0.087484913611943688, 0.00030408951716745367, 0.040339849873275067, 5.4131308141089158e-06, 1.1661327608584701e-07, 0.7388008857031686, +0.0056249999999999998, 2.3266045446486152, 0.73316170321133411, 790.35936761613914, 0.1039560084571382, 0.022296627215448003, 0.0011359165443020106, 0.083877529192241299, 0.0003760964167247532, 0.048153498250387908, 4.9567780860884836e-06, 7.4194975865213397e-08, 0.74019929295069586, +0.0062500000000000003, 2.4479214314061588, 0.81358991161259298, 833.49915753070206, 0.094342281079607085, 0.025576515598032364, 0.0016093379165598887, 0.079589020297861954, 0.00048650534096092017, 0.056968106191529824, 4.5523878389231024e-06, 6.7032298849189715e-08, 0.74142361415531022, +0.006875, 2.5659880801402859, 0.89380568407752115, 875.48144891531763, 0.084558504114132543, 0.02887564904789016, 0.0021587955911611337, 0.074692266304889551, 0.00063548478274258545, 0.066591533702783423, 4.1776136445993059e-06, 7.3383909261230379e-08, 0.74248351545884683, +0.0074999999999999997, 2.679397320457686, 0.97379492721512451, 915.7945483998285, 0.074891658980043241, 0.032074034490164896, 0.0027597011676134373, 0.069355990408358173, 0.00081848711048261211, 0.076693396035891495, 3.81958038131416e-06, 8.5179156847188386e-08, 0.74340282704790805, +0.0081250000000000003, 2.7866877324919019, 1.0535498699610373, 953.93280354121919, 0.065649346009905199, 0.035050331284822378, 0.0033820558165062698, 0.063819498245732242, 0.0010264911152755339, 0.086856929633320423, 3.4744641184932023e-06, 9.8937718133226469e-08, 0.74421177449260145, +0.0087500000000000008, 2.8865969655502388, 1.1330700792045947, 989.47817570227369, 0.057112015453486001, 0.037700270948181645, 0.0039952793491947487, 0.058347955566787769, 0.0012478902436022757, 0.096654188398265034, 3.1453332073885087e-06, 1.1269670724508324e-07, 0.74493914201056788, +0.0093750000000000014, 2.978255980578755, 1.2123615756592105, 1022.1581506196891, 0.04949150383519825, 0.039951274241403413, 0.004572962462803953, 0.053184021407520994, 0.0014712387392487179, 0.10571874839747665, 2.8387688304227487e-06, 1.2526018213844753e-07, 0.74560728688733546, +0.01, 3.0612837303484013, 1.2914342945805504, 1051.8698390915185, 0.042907963852313061, 0.041769137411653737, 0.0050962069769417999, 0.048511911927462799, 0.0016876111685787617, 0.11379359894307464, 2.5614991221790733e-06, 1.3606021498471971e-07, 0.74623087216063799, +0.010624999999999999, 3.1357463709104816, 1.3702984802990987, 1078.6604673289501, 0.037388614518914044, 0.043156275525769706, 0.0055548085787310315, 0.04444218060939626, 0.0018916553123911274, 0.12074576125966001, 2.3181795962567304e-06, 1.4503891580231648e-07, 0.74681824097662597, +0.01125, 3.2020415691409063, 1.448960567428119, 1102.6854217753801, 0.032883380517046756, 0.044144021520961363, 0.0059464574224023454, 0.041015676607536161, 0.0020813407191971642, 0.12655322433813762, 2.110570031301196e-06, 1.5246638696000219e-07, 0.74737363583830041, +0.012500000000000001, 3.3121878560158216, 1.6058500081811453, 1143.1575235307675, 0.026645013291340502, 0.045071139543787407, 0.0065352408653600912, 0.036086146168142229, 0.0024165220843660553, 0.13487019988675802, 1.8002115650996147e-06, 1.6398664135726015e-07, 0.74837377396203919, +0.014999999999999999, 3.4702707643220618, 1.918844109519684, 1203.094394898302, 0.021459053172201936, 0.044274863831088376, 0.007147633748197599, 0.031347684369649555, 0.0029540922952010031, 0.14283845879587859, 1.4614640368461466e-06, 1.8000888522746873e-07, 0.74997657231486081, +0.017500000000000002, 3.5867964469021296, 2.2311865115228948, 1249.1848940417526, 0.020006979500262648, 0.042317328181192077, 0.0073956784720450187, 0.02927945069501613, 0.0034003881190899214, 0.14632359267224729, 1.2812472513064534e-06, 1.8939004615031488e-07, 0.75127511172284955, +0.02, 3.6776605639689586, 2.5430579614213622, 1286.9080967562491, 0.019783353886530502, 0.040195001729350054, 0.0074791630594031227, 0.028246457542008514, 0.0037804917643500472, 0.14812446769993065, 1.1694833657180304e-06, 1.9264027655252762e-07, 0.75238970219478496, +0.025000000000000001, 3.8062218661088032, 3.1664641834273768, 1345.3375313534179, 0.020069469547008605, 0.03651060715338638, 0.0074468831539421894, 0.027138010755168095, 0.0043779084822032494, 0.15022672466350395, 1.0312082986825226e-06, 1.8631380371644787e-07, 0.75422917872268513, +0.029999999999999999, 3.8964523799099826, 3.7897282346058812, 1392.5425655885936, 0.020288096023165163, 0.033530769882025739, 0.0073296095403677644, 0.026317745985898339, 0.0048618842774350701, 0.15189880837836778, 9.3560224772902113e-07, 1.7143134876546083e-07, 0.75577197887914371, +0.035000000000000003, 3.9586752601947635, 4.4131032621018136, 1432.1310686638305, 0.020383361631021689, 0.031073965180691825, 0.0071783970425025514, 0.025637143310665959, 0.0052626043347806396, 0.15336027560540569, 8.6355205687248186e-07, 1.5425894260730429e-07, 0.75710323508393207, +0.040000000000000001, 3.9988577414527242, 5.0367609497864478, 1466.1996248048285, 0.020389347273240155, 0.028998808287582837, 0.0070131092893370092, 0.025050347570421145, 0.005599546789696843, 0.1546730881166766, 8.0648696809369336e-07, 1.3774350413201087e-07, 0.75827480844257322, +0.044999999999999998, 4.0207699805593204, 5.6608255423746447, 1496.0795406181303, 0.020334286089737273, 0.0272108712831849, 0.0068430370318006972, 0.024533761724802372, 0.0058861543287198007, 0.15586938150908675, 7.5970999949226152e-07, 1.2296533980778867e-07, 0.75932162535732894, +0.050000000000000003, 4.0268319709231264, 6.285375160384584, 1522.6273727602772, 0.020236190055544903, 0.025646423545325259, 0.0066728437840372929, 0.024072615160771714, 0.0061315752166984554, 0.15697079451993964, 7.2042625786976623e-07, 1.1016373624666346e-07, 0.76026872712768845, +0.059999999999999998, 3.9954715076305489, 7.5370961306618929, 1567.046270226188, 0.019954345924716552, 0.023065449653235648, 0.0063476810280637197, 0.023298938513315232, 0.0065164824393074494, 0.15890690323181508, 6.5900286262737154e-07, 9.030833660433104e-08, 0.76190944989834697, +0.070000000000000007, 3.9199086361114239, 8.7921414621501324, 1605.1278754471107, 0.019605495248628049, 0.020918183272324184, 0.0060358049694081524, 0.022631356937731825, 0.0068168665197777716, 0.16064815117619569, 6.1002837933327518e-07, 7.5363354185294861e-08, 0.76334345648420088, +0.080000000000000002, 3.8049375869106234, 10.051244722401661, 1638.6518722532867, 0.01921215737409027, 0.019078761538933977, 0.0057379738026155837, 0.022040084767416552, 0.0070538238712225072, 0.16224743328628899, 5.6938843443230727e-07, 6.3933003469375504e-08, 0.76462913203799421, +0.089999999999999997, 3.6536509178782928, 11.315225491227944, 1668.8125855596343, 0.01878458929576025, 0.017463381881796032, 0.0054524610901615653, 0.021505386096201135, 0.0072411499205829889, 0.16374349727944509, 5.346216791715588e-07, 5.4992511125819936e-08, 0.76580894482186268, +0.10000000000000001, 3.4680905100725519, 12.585101170668079, 1696.4504931194019, 0.018325156210624388, 0.016012300298996442, 0.0051763571439689508, 0.021013518129501477, 0.0073877137268993183, 0.1651673557628989, 5.0415621573349164e-07, 4.7836691263300553e-08, 0.76691704673420358, +0.1125, 3.1897572328325454, 14.183253752635501, 1728.2595546684843, 0.017702815552013099, 0.01436751073581297, 0.004840379923476584, 0.020446727014850764, 0.0075210885953569985, 0.16687912245690117, 4.7066897693984347e-07, 4.0705593822036701e-08, 0.76824184434701759, +0.125, 2.861891624737785, 15.795003950061801, 1758.1437154169494, 0.017008873446908601, 0.012834780716031167, 0.0045043196205411789, 0.019913007753461304, 0.0076043351617225618, 0.1685777645236769, 4.4065046699408111e-07, 3.4952746498433607e-08, 0.76955644317444472, +0.13750000000000001, 2.4847826693992552, 17.422178299315796, 1786.4307271301702, 0.016214106766358791, 0.011359025308892086, 0.0041573022595386448, 0.019397936207958346, 0.0076304747593489439, 0.17032860339458941, 4.1326202060478539e-07, 3.0217663503250393e-08, 0.77091210782362962, +0.14374999999999999, 2.2770558401831931, 18.238730311716395, 1799.4957684600345, 0.01576109307190152, 0.010624589461568681, 0.0039730304563995993, 0.019139495836953978, 0.0076104426812164004, 0.17126681503003738, 4.0054729635656477e-07, 2.8201819752163684e-08, 0.77162410471280629, +0.14999999999999999, 2.0552061620168458, 19.052921126670338, 1810.5919684430958, 0.015258249416852306, 0.0098889207849126963, 0.0037761882407849902, 0.018872610325050793, 0.0075502341824138326, 0.17229090540397235, 3.8886593343690788e-07, 2.6494945224781377e-08, 0.77236247628513444, +0.15625, 1.8174654217941273, 19.850433363629929, 1817.5024674967121, 0.014691821003212957, 0.0091514721965593789, 0.0035595308434884097, 0.018582851226749995, 0.0074213646753439437, 0.1734671658244227, 3.7883331302056433e-07, 2.5237768341404418e-08, 0.77312539015914128, +0.16249999999999998, 1.5614276192033629, 20.595797640087845, 1815.9634137766275, 0.014042603675041685, 0.0084160994397404851, 0.0033122846680085872, 0.018244182112174247, 0.0071743076809984716, 0.17490452992836494, 3.7169304991284128e-07, 2.472084022324292e-08, 0.77390559608178133, +0.16874999999999998, 1.2848439006398622, 21.209859407347761, 1798.5892811579283, 0.013283368976569751, 0.0076923251770785262, 0.0030196337193669595, 0.017812154998107566, 0.0067323576952801889, 0.17676852736487453, 3.7012790522676784e-07, 2.5558235689681428e-08, 0.77469123638258153, +0.171875, 1.1382164187730344, 21.407485658921317, 1779.7574029601565, 0.012850174665669191, 0.0073442347524166835, 0.0028499156000188529, 0.017537699148039199, 0.0063998132290477884, 0.17793893396867863, 3.7354604756653276e-07, 2.6968402038326325e-08, 0.77507882812168005, +0.17499999999999999, 0.98617863747333057, 21.472995814642747, 1750.8026865586155, 0.012368663794682943, 0.0070092161297527038, 0.0026586378356279498, 0.017200588887095837, 0.005964803929325484, 0.17933571405172963, 3.8192394076185316e-07, 2.9619106009921667e-08, 0.77546196382873878, +0.17812499999999998, 0.82987840405755042, 21.342693402582224, 1708.920457562846, 0.011827497030354329, 0.0066922628574613459, 0.0024421816092741406, 0.016777956314181237, 0.0054108532609216308, 0.18100812833678023, 3.9768060034609105e-07, 3.4403235968189336e-08, 0.77584068850719079, +0.18124999999999999, 0.67145766561153086, 20.930026669154678, 1650.9187320412996, 0.01121094714556985, 0.0063999884360225511, 0.0021968259265958976, 0.016238333774208673, 0.004726786471278833, 0.18301005080342692, 4.239167153159463e-07, 4.3171961213720126e-08, 0.7762166003542208, +0.18281249999999999, 0.59241465081008782, 20.576304998909894, 1614.6826570688215, 0.010869428332468421, 0.0062690462242506809, 0.0020631379373647345, 0.015914510501020544, 0.0043360226155424028, 0.18414452905987913, 4.4220424487953039e-07, 5.0231793717051413e-08, 0.77640283289343559, +0.18437500000000001, 0.51418210557136812, 20.096337097725328, 1572.7818220534091, 0.010499516219101527, 0.0061489010171829573, 0.0019203914552325223, 0.015543046151993561, 0.0039119586132670024, 0.18538576658981648, 4.6466342638672339e-07, 6.0268474295514317e-08, 0.77658989502150522, +0.18593750000000001, 0.43747526454577007, 19.467524920795679, 1524.6587790557048, 0.010097322355126625, 0.0060411761126855038, 0.0017680119657623023, 0.015116197179070714, 0.0034578496485199511, 0.18674039198962022, 4.9187277560889051e-07, 7.4777006286711021e-08, 0.77677848409943273, +0.1875, 0.36313147318575267, 18.663952131597618, 1469.7028961216049, 0.0096599499725037001, 0.0059471578843499923, 0.0016053677264801351, 0.014627317211017628, 0.0029787964720486, 0.18821165780621549, 5.2461874093902912e-07, 9.5899412933792073e-08, 0.77696913240923049, +0.18906249999999999, 0.29211558898522427, 17.655691683210282, 1407.2176055311322, 0.0091880878301557421, 0.0058668045314332245, 0.001431853013839542, 0.014074202077233233, 0.0024822813458974177, 0.18979426594627039, 5.6425602417920409e-07, 1.2591539600027082e-07, 0.77716181508375037, +0.19062499999999999, 0.22552511951869023, 16.407696691882055, 1336.3603094559869, 0.0086906989849737434, 0.0057969220050752036, 0.0012471953688583976, 0.013464940328317351, 0.0019790087550489222, 0.19146517539732993, 6.1356155222814421e-07, 1.6400687590482523e-07, 0.77735528159196821, +0.19218750000000001, 0.1645976230947685, 14.877911575022168, 1256.0367661476128, 0.0081909702178872418, 0.0057288379555029897, 0.0010519876682219241, 0.012825018359898287, 0.0014845636057194014, 0.19317153730876621, 6.786151771662939e-07, 1.9750624835605178e-07, 0.77754620876257841, +0.19375000000000001, 0.11072490825613812, 13.013855475512626, 1164.7223755227883, 0.0077290133943349916, 0.0056476349129977134, 0.00084754537479999519, 0.012201113865899027, 0.0010215005391175582, 0.19482341064247913, 7.7221127964360055e-07, 1.9880727301944028e-07, 0.77772881025181895, +0.19453124999999999, 0.086875428136627034, 11.93460634789748, 1114.2863211040678, 0.0075339788614498492, 0.0055978169470309027, 0.00074353045866084656, 0.011925550981397503, 0.00081492908920098573, 0.19557048118163842, 8.3746894298422668e-07, 1.8375684101845595e-07, 0.77781269125483743, +0.1953125, 0.065387567171209762, 10.742546989325882, 1060.1273214713146, 0.0073659767147021436, 0.0055393514115634184, 0.00063712102757782028, 0.011676583401477824, 0.00062998355453721723, 0.19625830348606668, 9.2114254811505824e-07, 1.6155744412982579e-07, 0.77789159770408267, +0.19609375000000001, 0.046491650612488229, 9.4240295569474029, 1001.7159138255283, 0.0072286581832338185, 0.0054717508677656398, 0.0005282439022776421, 0.011458866674890301, 0.00047344314576905464, 0.19687325153174595, 1.0309930113041476e-06, 1.3903662312876773e-07, 0.77796461566468311, +0.19687500000000002, 0.030433529836239048, 7.9616396821344377, 938.32479698019745, 0.0071232901356660409, 0.0053960875040308543, 0.00041770708847931298, 0.011274293816151337, 0.00035204766088763826, 0.19740448049877671, 1.1794647734858792e-06, 1.2565283648166381e-07, 0.7780307881783981, +0.19765625000000003, 0.017476039086012533, 6.3321208347116018, 868.8943893747105, 0.00704894629106303, 0.00531522695439662, 0.00030910135669143907, 0.011121153230818878, 0.00027249614925341632, 0.19784284912739869, 1.3879749381188476e-06, 1.3546455318674106e-07, 0.77808870345088665, +0.19843750000000002, 0.0079005409017527777, 4.5024856494690031, 791.76943205355201, 0.00700339874905344, 0.0052345400860735001, 0.00021113489053437893, 0.010995469037131909, 0.00024035735830448758, 0.1981773335928593, 1.6976823975704344e-06, 1.9258304053971736e-07, 0.77813587602060508, +0.19921875, 0.0020063660064339448, 2.421834060377221, 704.12243020995618, 0.0069839268363947743, 0.0051652738912929385, 0.0001390425432120124, 0.010899112629199329, 0.00025394096488660559, 0.19838825652559391, 2.2022302405485384e-06, 3.286281178233616e-07, 0.77816791575106214, +0.19960937500000001, 0.0004929236741182477, 1.2618846057426176, 654.65258598081732, 0.00698154145099243, 0.005146641929851938, 0.00012174475845906474, 0.010874085721525416, 0.00026462101081261327, 0.19843313490621356, 2.5852186402687997e-06, 3.9434196095075079e-07, 0.7781752506615438, +0.20000000000000001, 0.0, 0.0, 600.0, 0.0069815359087274671, 0.0051466374871214516, 0.00012174471378542655, 0.010874082917076219, 0.00026462086108068578, 0.19843304230836817, 2.5852177095946074e-06, 3.9434175287653015e-07, 0.77817535624437806,