Nanosecond Pulse Plasma Simulation#

This example simulates a nanosecond-scale pulse discharge in a reactor. A Gaussian-shaped electric field pulse is applied over a short timescale. The plasma reaction mechanism used is based on Colin Pavan’s mechanism for methane-air plasmas which is described in his Ph.D. dissertation and the corresponding AIAA SciTech conference papers:

C. A. Pavan, “Nanosecond Pulsed Plasmas in Dynamic Combustion Environments,” Ph.D. Thesis, Massachusetts Institute of Technology, 2023. Chapter 5.

C. A. Pavan and C. Guerra-Garcia, “Modelling the Impact of a Repetitively Pulsed Nanosecond DBD Plasma on a Mesoscale Flame,” in AIAA SCITECH 2022 Forum, Reston, Virginia: American Institute of Aeronautics and Astronautics, Jan. 2022, pp. 1-15. DOI: 10.2514/6.2022-0975.

C. A. Pavan and C. Guerra-Garcia, “Modeling Flame Speed Modification by Nanosecond Pulsed Discharges to Inform Experimental Design,” in AIAA SCITECH 2023 Forum, Reston, Virginia: American Institute of Aeronautics and Astronautics, Jan. 2023, pp. 1-15. DOI: 10.2514/6.2023-2056.

Requires: cantera >= 3.2, matplotlib >= 2.0

Tags: Python plasma reactor network

nanosecond pulse discharge
No existing EEDF. Using first guess method: maxwell
First guess EEDF maxwell
     t [s]      T [K]     P [Pa]       h [J/kg]
 1.000e-10    300.000 101325.000 -254430.994561
 2.000e-10    300.000 101325.000 -254430.994561
 3.000e-10    300.000 101325.000 -254430.994561
 4.000e-10    300.000 101325.000 -254430.994561
 5.000e-10    300.000 101325.000 -254430.994561
 6.000e-10    300.000 101325.000 -254430.994561
 7.000e-10    300.000 101325.000 -254430.994561
 8.000e-10    300.000 101325.000 -254430.994561
 9.000e-10    300.000 101325.000 -254430.994561
 1.000e-09    300.000 101325.000 -254430.994561
 1.100e-09    300.000 101325.000 -254430.994561
 1.200e-09    300.000 101325.000 -254430.994561
 1.300e-09    300.000 101325.000 -254430.994561
 1.400e-09    300.000 101325.000 -254430.994561
 1.500e-09    300.000 101325.000 -254430.994561
 1.600e-09    300.000 101325.000 -254430.994561
 1.700e-09    300.000 101325.000 -254430.994561
 1.800e-09    300.000 101325.000 -254430.994561
 1.900e-09    300.000 101325.000 -254430.994561
 2.000e-09    300.000 101325.000 -254430.994561
 2.100e-09    300.000 101325.000 -254430.994561
 2.200e-09    300.000 101325.000 -254430.994561
 2.300e-09    300.000 101325.000 -254430.994561
 2.400e-09    300.000 101325.000 -254430.994561
 2.500e-09    300.000 101325.000 -254430.994561
 2.600e-09    300.000 101325.000 -254430.994561
 2.700e-09    300.000 101325.000 -254430.994561
 2.800e-09    300.000 101325.000 -254430.994561
 2.900e-09    300.000 101325.000 -254430.994561
 3.000e-09    300.000 101325.000 -254430.994561
 3.100e-09    300.000 101325.000 -254430.994561
 3.200e-09    300.000 101325.000 -254430.994561
 3.300e-09    300.000 101325.000 -254430.994561
 3.400e-09    300.000 101325.000 -254430.994561
 3.500e-09    300.000 101325.000 -254430.994561
 3.600e-09    300.000 101325.000 -254430.994561
 3.700e-09    300.000 101325.000 -254430.994561
 3.800e-09    300.000 101325.000 -254430.994561
 3.900e-09    300.000 101325.000 -254430.994561
 4.000e-09    300.000 101325.000 -254430.994561
 4.100e-09    300.000 101325.000 -254430.994561
 4.200e-09    300.000 101325.000 -254430.994561
 4.300e-09    300.000 101325.000 -254430.994561
 4.400e-09    300.000 101325.000 -254430.994561
 4.500e-09    300.000 101325.000 -254430.994561
 4.600e-09    300.000 101325.000 -254430.994561
 4.700e-09    300.000 101325.000 -254430.994561
 4.800e-09    300.000 101325.000 -254430.994561
 4.900e-09    300.000 101325.000 -254430.994561
 5.000e-09    300.000 101325.000 -254430.994561
 5.100e-09    300.000 101325.000 -254430.994561
 5.200e-09    300.000 101325.000 -254430.994561
 5.300e-09    300.000 101325.000 -254430.994561
 5.400e-09    300.000 101325.000 -254430.994561
 5.500e-09    300.000 101325.000 -254430.994561
 5.600e-09    300.000 101325.000 -254430.994561
 5.700e-09    300.000 101325.000 -254430.994561
 5.800e-09    300.000 101325.000 -254430.994561
 5.900e-09    300.000 101325.000 -254430.994561
 6.000e-09    300.000 101325.000 -254430.994561
 6.100e-09    300.000 101325.000 -254430.994561
 6.200e-09    300.000 101325.000 -254430.994561
 6.300e-09    300.000 101325.000 -254430.994561
 6.400e-09    300.000 101325.000 -254430.994561
 6.500e-09    300.000 101325.000 -254430.994561
 6.600e-09    300.000 101325.000 -254430.994561
 6.700e-09    300.000 101325.000 -254430.994561
 6.800e-09    300.000 101325.000 -254430.994561
 6.900e-09    300.000 101325.000 -254430.994561
 7.000e-09    300.000 101325.000 -254430.994561
 7.100e-09    300.000 101325.000 -254430.994561
 7.200e-09    300.000 101325.000 -254430.994561
 7.300e-09    300.000 101325.000 -254430.994561
 7.400e-09    300.000 101325.000 -254430.994561
 7.500e-09    300.000 101325.000 -254430.994561
 7.600e-09    300.000 101325.000 -254430.994561
 7.700e-09    300.000 101325.000 -254430.994561
 7.800e-09    300.000 101325.000 -254430.994561
 7.900e-09    300.000 101325.000 -254430.994561
 8.000e-09    300.000 101325.000 -254430.994561
 8.100e-09    300.000 101325.000 -254430.994561
 8.200e-09    300.000 101325.000 -254430.994561
 8.300e-09    300.000 101325.000 -254430.994561
 8.400e-09    300.000 101325.000 -254430.994561
 8.500e-09    300.000 101325.000 -254430.994561
 8.600e-09    300.000 101325.000 -254430.994561
 8.700e-09    300.000 101325.000 -254430.994561
 8.800e-09    300.000 101325.000 -254430.994561
 8.900e-09    300.000 101325.000 -254430.994561
 9.000e-09    300.000 101325.000 -254430.994561
 9.100e-09    300.000 101325.000 -254430.994561
 9.200e-09    300.000 101325.000 -254430.994561
 9.300e-09    300.000 101325.000 -254430.994561
 9.400e-09    300.000 101325.000 -254430.994561
 9.500e-09    300.000 101325.000 -254430.994561
 9.600e-09    300.000 101325.000 -254430.994561
 9.700e-09    300.000 101325.000 -254430.994561
 9.800e-09    300.000 101325.000 -254430.994561
 9.900e-09    300.000 101325.000 -254430.994561
 1.000e-08    300.000 101325.000 -254430.994561
 1.010e-08    300.000 101325.000 -254430.994561
 1.020e-08    300.000 101325.000 -254430.994561
 1.030e-08    300.000 101325.000 -254430.994561
 1.040e-08    300.000 101325.000 -254430.994561
 1.050e-08    300.000 101325.000 -254430.994561
 1.060e-08    300.000 101325.000 -254430.994561
 1.070e-08    300.000 101325.000 -254430.994561
 1.080e-08    300.000 101325.000 -254430.994561
 1.090e-08    300.000 101325.000 -254430.994561
 1.100e-08    300.000 101325.000 -254430.994561
 1.110e-08    300.000 101325.000 -254430.994561
 1.120e-08    300.000 101325.000 -254430.994561
 1.130e-08    300.000 101325.000 -254430.994561
 1.140e-08    300.000 101325.000 -254430.994561
 1.150e-08    300.000 101325.000 -254430.994561
 1.160e-08    300.000 101325.000 -254430.994561
 1.170e-08    300.000 101325.000 -254430.994561
 1.180e-08    300.000 101325.000 -254430.994561
 1.190e-08    300.000 101325.000 -254430.994561
 1.200e-08    300.000 101325.000 -254430.994561
 1.210e-08    300.000 101325.000 -254430.994561
 1.220e-08    300.000 101325.000 -254430.994561
 1.230e-08    300.000 101325.000 -254430.994561
 1.240e-08    300.000 101325.000 -254430.994561
 1.250e-08    300.000 101325.000 -254430.994561
 1.260e-08    300.000 101325.000 -254430.994561
 1.270e-08    300.000 101325.000 -254430.994561
 1.280e-08    300.000 101325.000 -254430.994561
 1.290e-08    300.000 101325.000 -254430.994561
 1.300e-08    300.000 101325.000 -254430.994561
 1.310e-08    300.000 101325.000 -254430.994561
 1.320e-08    300.000 101325.000 -254430.994561
 1.330e-08    300.000 101325.000 -254430.994561
 1.340e-08    300.000 101325.000 -254430.994561
 1.350e-08    300.000 101325.000 -254430.994561
 1.360e-08    300.000 101325.000 -254430.994561
 1.370e-08    300.000 101325.000 -254430.994561
 1.380e-08    300.000 101325.000 -254430.994561
 1.390e-08    300.000 101325.000 -254430.994561
 1.400e-08    300.000 101325.000 -254430.994561
 1.410e-08    300.000 101325.000 -254430.994561
 1.420e-08    300.000 101325.000 -254430.994555
 1.430e-08    300.000 101325.000 -254430.994555
 1.440e-08    300.000 101325.000 -254430.994555
 1.450e-08    300.000 101325.000 -254430.994555
 1.460e-08    300.000 101325.000 -254430.994555
 1.470e-08    300.000 101325.000 -254430.994555
 1.480e-08    300.000 101325.000 -254430.994555
 1.490e-08    300.000 101325.000 -254430.994555
 1.500e-08    300.000 101325.000 -254430.994555
 1.510e-08    300.000 101325.000 -254430.994555
 1.520e-08    300.000 101325.000 -254430.994537
 1.530e-08    300.000 101325.000 -254430.994537
 1.540e-08    300.000 101325.000 -254430.994537
 1.550e-08    300.000 101325.000 -254430.994537
 1.560e-08    300.000 101325.000 -254430.994537
 1.570e-08    300.000 101325.000 -254430.994537
 1.580e-08    300.000 101325.000 -254430.994537
 1.590e-08    300.000 101325.000 -254430.994537
 1.600e-08    300.000 101325.000 -254430.994537
 1.610e-08    300.000 101325.000 -254430.994537
 1.620e-08    300.000 101325.000 -254430.994535
 1.630e-08    300.000 101325.000 -254430.994535
 1.640e-08    300.000 101325.000 -254430.994535
 1.650e-08    300.000 101325.000 -254430.994535
 1.660e-08    300.000 101325.000 -254430.994535
 1.670e-08    300.000 101325.000 -254430.994535
 1.680e-08    300.000 101325.000 -254430.994535
 1.690e-08    300.000 101325.000 -254430.994535
 1.700e-08    300.000 101325.000 -254430.994535
 1.710e-08    300.000 101325.000 -254430.994535
 1.720e-08    300.000 101325.000 -254430.994533
 1.730e-08    300.000 101325.000 -254430.994533
 1.740e-08    300.000 101325.000 -254430.994533
 1.750e-08    300.000 101325.000 -254430.994533
 1.760e-08    300.000 101325.000 -254430.994533
 1.770e-08    300.000 101325.000 -254430.994533
 1.780e-08    300.000 101325.000 -254430.994533
 1.790e-08    300.000 101325.000 -254430.994533
 1.800e-08    300.000 101325.000 -254430.994533
 1.810e-08    300.000 101325.000 -254430.994533
 1.820e-08    300.000 101325.000 -254430.994530
 1.830e-08    300.000 101325.000 -254430.994530
 1.840e-08    300.000 101325.000 -254430.994530
 1.850e-08    300.000 101325.000 -254430.994530
 1.860e-08    300.000 101325.000 -254430.994530
 1.870e-08    300.000 101325.000 -254430.994530
 1.880e-08    300.000 101325.000 -254430.994530
 1.890e-08    300.000 101325.000 -254430.994530
 1.900e-08    300.000 101325.000 -254430.994530
 1.910e-08    300.000 101325.000 -254430.994530
 1.920e-08    300.000 101325.000 -254430.994525
 1.930e-08    300.000 101325.000 -254430.994525
 1.940e-08    300.000 101325.000 -254430.994525
 1.950e-08    300.000 101325.000 -254430.994525
 1.960e-08    300.000 101325.000 -254430.994525
 1.970e-08    300.000 101325.000 -254430.994525
 1.980e-08    300.000 101325.000 -254430.994525
 1.990e-08    300.000 101325.000 -254430.994525
 2.000e-08    300.000 101325.000 -254430.994525
 2.010e-08    300.000 101325.000 -254430.994525
 2.020e-08    300.000 101325.000 -254430.994509
 2.030e-08    300.000 101325.000 -254430.994504
 2.040e-08    300.000 101325.000 -254430.994499
 2.050e-08    300.000 101325.000 -254430.994494
 2.060e-08    300.000 101325.000 -254430.994489
 2.070e-08    300.000 101325.000 -254430.994484
 2.080e-08    300.000 101325.000 -254430.994480
 2.090e-08    300.000 101325.000 -254430.994475
 2.100e-08    300.000 101325.000 -254430.994471
 2.110e-08    300.000 101325.000 -254430.994304
 2.120e-08    300.000 101325.000 -254430.994185
 2.130e-08    300.000 101325.000 -254430.994072
 2.140e-08    300.000 101325.000 -254430.993966
 2.150e-08    300.000 101325.000 -254430.993864
 2.160e-08    300.000 101325.000 -254430.993766
 2.170e-08    300.000 101325.000 -254430.993670
 2.180e-08    300.000 101325.000 -254430.993576
 2.190e-08    300.000 101325.000 -254430.993484
 2.200e-08    300.000 101325.000 -254430.993393
 2.210e-08    300.000 101325.000 -254430.992686
 2.220e-08    300.000 101325.000 -254430.992087
 2.230e-08    300.000 101325.000 -254430.991511
 2.240e-08    300.000 101325.000 -254430.990953
 2.250e-08    300.000 101325.000 -254430.990407
 2.260e-08    300.000 101325.000 -254430.989867
 2.270e-08    300.000 101325.000 -254430.989330
 2.280e-08    300.000 101325.000 -254430.988794
 2.290e-08    300.000 101325.000 -254430.988257
 2.300e-08    300.000 101325.000 -254430.987716
 2.310e-08    300.000 101325.000 -254430.985960
 2.320e-08    300.000 101325.000 -254430.984205
 2.330e-08    300.000 101325.000 -254430.982346
 2.340e-08    300.000 101325.000 -254430.980361
 2.350e-08    300.000 101325.000 -254430.978227
 2.360e-08    300.000 101325.000 -254430.975923
 2.370e-08    300.000 101325.000 -254430.973426
 2.380e-08    300.000 101325.000 -254430.970714
 2.390e-08    300.000 101325.000 -254430.967765
 2.400e-08    300.000 101325.000 -254430.964552
 2.410e-08    300.000 101325.000 -254430.957984
 2.420e-08    300.000 101325.000 -254430.950087
 2.430e-08    300.000 101325.000 -254430.940382
 2.440e-08    300.000 101325.000 -254430.928434
 2.450e-08    300.000 101325.000 -254430.913707
 2.460e-08    300.000 101325.000 -254430.895542
 2.470e-08    300.000 101325.000 -254430.873131
 2.480e-08    300.000 101325.000 -254430.845484
 2.490e-08    300.000 101325.000 -254430.811382
 2.500e-08    300.000 101325.000 -254430.769333
 2.510e-08    300.000 101325.000 -254430.705848
 2.520e-08    300.000 101325.000 -254430.622752
 2.530e-08    300.000 101325.000 -254430.513571
 2.540e-08    300.000 101325.000 -254430.370351
 2.550e-08    300.000 101325.000 -254430.182864
 2.560e-08    300.000 101325.000 -254429.938060
 2.570e-08    300.000 101325.000 -254429.619431
 2.580e-08    300.000 101325.000 -254429.206345
 2.590e-08    300.000 101325.000 -254428.673352
 2.600e-08    300.000 101325.000 -254427.989538
 2.610e-08    300.000 101325.000 -254427.280585
 2.620e-08    300.000 101325.000 -254426.434444
 2.630e-08    300.000 101325.000 -254425.435167
 2.640e-08    300.000 101325.000 -254424.259698
 2.650e-08    300.000 101325.000 -254422.882153
 2.660e-08    300.000 101325.000 -254421.273201
 2.670e-08    300.000 101325.000 -254419.399234
 2.680e-08    300.000 101325.000 -254417.221353
 2.690e-08    300.000 101325.000 -254414.694183
 2.700e-08    300.000 101325.000 -254411.764568
 2.710e-08    300.000 101325.000 -254410.091453
 2.720e-08    300.000 101325.000 -254408.252493
 2.730e-08    300.000 101325.000 -254406.289311
 2.740e-08    300.000 101325.000 -254404.186998
 2.750e-08    300.000 101325.000 -254401.936330
 2.760e-08    300.000 101325.000 -254399.530165
 2.770e-08    300.000 101325.000 -254396.961679
 2.780e-08    300.000 101325.000 -254394.223561
 2.790e-08    300.000 101325.000 -254391.307670
 2.800e-08    300.000 101325.000 -254388.204919
 2.810e-08    300.000 101325.000 -254387.812141
 2.820e-08    300.000 101325.000 -254387.072937
 2.830e-08    300.000 101325.000 -254386.204563
 2.840e-08    300.000 101325.000 -254385.232204
 2.850e-08    300.000 101325.000 -254384.176298
 2.860e-08    300.000 101325.000 -254383.052738
 2.870e-08    300.000 101325.000 -254381.873751
 2.880e-08    300.000 101325.000 -254380.648702
 2.890e-08    300.000 101325.000 -254379.384762
 2.900e-08    300.000 101325.000 -254378.087411
 2.910e-08    300.000 101325.000 -254378.690528
 2.920e-08    300.000 101325.000 -254378.922205
 2.930e-08    300.000 101325.000 -254379.043458
 2.940e-08    300.000 101325.000 -254379.079556
 2.950e-08    300.000 101325.000 -254379.049924
 2.960e-08    300.000 101325.000 -254378.969502
 2.970e-08    300.000 101325.000 -254378.849786
 2.980e-08    300.000 101325.000 -254378.699625
 2.990e-08    300.000 101325.000 -254378.525839
 3.000e-08    300.000 101325.000 -254378.333688
 3.010e-08    300.000 101325.000 -254378.644397
 3.020e-08    300.000 101325.000 -254378.792110
 3.030e-08    300.000 101325.000 -254378.907356
 3.040e-08    300.000 101325.000 -254378.997479
 3.050e-08    300.000 101325.000 -254379.068137
 3.060e-08    300.000 101325.000 -254379.123690
 3.070e-08    300.000 101325.000 -254379.167507
 3.080e-08    300.000 101325.000 -254379.202185
 3.090e-08    300.000 101325.000 -254379.229735
 3.100e-08    300.000 101325.000 -254379.251711
 3.110e-08    300.000 101325.000 -254379.319238
 3.120e-08    300.000 101325.000 -254379.350886
 3.130e-08    300.000 101325.000 -254379.378762
 3.140e-08    300.000 101325.000 -254379.403677
 3.150e-08    300.000 101325.000 -254379.426261
 3.160e-08    300.000 101325.000 -254379.447004
 3.170e-08    300.000 101325.000 -254379.466287
 3.180e-08    300.000 101325.000 -254379.484409
 3.190e-08    300.000 101325.000 -254379.501604
 3.200e-08    300.000 101325.000 -254379.518056
 3.210e-08    300.000 101325.000 -254379.548677
 3.220e-08    300.000 101325.000 -254379.564100
 3.230e-08    300.000 101325.000 -254379.579124
 3.240e-08    300.000 101325.000 -254379.593822
 3.250e-08    300.000 101325.000 -254379.608251
 3.260e-08    300.000 101325.000 -254379.622457
 3.270e-08    300.000 101325.000 -254379.636477
 3.280e-08    300.000 101325.000 -254379.650341
 3.290e-08    300.000 101325.000 -254379.664073
 3.300e-08    300.000 101325.000 -254379.677692
 3.310e-08    300.000 101325.000 -254379.702263
 3.320e-08    300.000 101325.000 -254379.715703
 3.330e-08    300.000 101325.000 -254379.729070
 3.340e-08    300.000 101325.000 -254379.742374
 3.350e-08    300.000 101325.000 -254379.755623
 3.360e-08    300.000 101325.000 -254379.768821
 3.370e-08    300.000 101325.000 -254379.781976
 3.380e-08    300.000 101325.000 -254379.795090
 3.390e-08    300.000 101325.000 -254379.808168
 3.400e-08    300.000 101325.000 -254379.821213
 3.410e-08    300.000 101325.000 -254379.840934
 3.420e-08    300.000 101325.000 -254379.853920
 3.430e-08    300.000 101325.000 -254379.866881
 3.440e-08    300.000 101325.000 -254379.879817
 3.450e-08    300.000 101325.000 -254379.892730
 3.460e-08    300.000 101325.000 -254379.905621
 3.470e-08    300.000 101325.000 -254379.918492
 3.480e-08    300.000 101325.000 -254379.931343
 3.490e-08    300.000 101325.000 -254379.944175
 3.500e-08    300.000 101325.000 -254379.956990
 3.510e-08    300.000 101325.000 -254379.975013
 3.520e-08    300.000 101325.000 -254379.987792
 3.530e-08    300.000 101325.000 -254380.000556
 3.540e-08    300.000 101325.000 -254380.013303
 3.550e-08    300.000 101325.000 -254380.026035
 3.560e-08    300.000 101325.000 -254380.038751
 3.570e-08    300.000 101325.000 -254380.051452
 3.580e-08    300.000 101325.000 -254380.064139
 3.590e-08    300.000 101325.000 -254380.076811
 3.600e-08    300.000 101325.000 -254380.089468
 3.610e-08    300.000 101325.000 -254380.159215
 3.620e-08    300.000 101325.000 -254380.171845
 3.630e-08    300.000 101325.000 -254380.184461
 3.640e-08    300.000 101325.000 -254380.197063
 3.650e-08    300.000 101325.000 -254380.209652
 3.660e-08    300.000 101325.000 -254380.222227
 3.670e-08    300.000 101325.000 -254380.234788
 3.680e-08    300.000 101325.000 -254380.247336
 3.690e-08    300.000 101325.000 -254380.259871
 3.700e-08    300.000 101325.000 -254380.272393
 3.710e-08    300.000 101325.000 -254380.305677
 3.720e-08    300.000 101325.000 -254380.318173
 3.730e-08    300.000 101325.000 -254380.330656
 3.740e-08    300.000 101325.000 -254380.343126
 3.750e-08    300.000 101325.000 -254380.355583
 3.760e-08    300.000 101325.000 -254380.368027
 3.770e-08    300.000 101325.000 -254380.380458
 3.780e-08    300.000 101325.000 -254380.392876
 3.790e-08    300.000 101325.000 -254380.405282
 3.800e-08    300.000 101325.000 -254380.417674
 3.810e-08    300.000 101325.000 -254380.430054
 3.820e-08    300.000 101325.000 -254380.442421
 3.830e-08    300.000 101325.000 -254380.454775
 3.840e-08    300.000 101325.000 -254380.467117
 3.850e-08    300.000 101325.000 -254380.479445
 3.860e-08    300.000 101325.000 -254380.491762
 3.870e-08    300.000 101325.000 -254380.504065
 3.880e-08    300.000 101325.000 -254380.516356
 3.890e-08    300.000 101325.000 -254380.528635
 3.900e-08    300.000 101325.000 -254380.540900
 3.910e-08    300.000 101325.000 -254380.553153
 3.920e-08    300.000 101325.000 -254380.565394
 3.930e-08    300.000 101325.000 -254380.577622
 3.940e-08    300.000 101325.000 -254380.589838
 3.950e-08    300.000 101325.000 -254380.602041
 3.960e-08    300.000 101325.000 -254380.614231
 3.970e-08    300.000 101325.000 -254380.626410
 3.980e-08    300.000 101325.000 -254380.638575
 3.990e-08    300.000 101325.000 -254380.650729
 4.000e-08    300.000 101325.000 -254380.662870
 4.010e-08    300.000 101325.000 -254380.674998
 4.020e-08    300.000 101325.000 -254380.687114
 4.030e-08    300.000 101325.000 -254380.699218
 4.040e-08    300.000 101325.000 -254380.711309
 4.050e-08    300.000 101325.000 -254380.723388
 4.060e-08    300.000 101325.000 -254380.735455
 4.070e-08    300.000 101325.000 -254380.747509
 4.080e-08    300.000 101325.000 -254380.759552
 4.090e-08    300.000 101325.000 -254380.771581
 4.100e-08    300.000 101325.000 -254380.783599
 4.110e-08    300.000 101325.000 -254380.795604
 4.120e-08    300.000 101325.000 -254380.807597
 4.130e-08    300.000 101325.000 -254380.819578
 4.140e-08    300.000 101325.000 -254380.831547
 4.150e-08    300.000 101325.000 -254380.843504
 4.160e-08    300.000 101325.000 -254380.855448
 4.170e-08    300.000 101325.000 -254380.867380
 4.180e-08    300.000 101325.000 -254380.879300
 4.190e-08    300.000 101325.000 -254380.891208
 4.200e-08    300.000 101325.000 -254380.903104
 4.210e-08    300.000 101325.000 -254380.914987
 4.220e-08    300.000 101325.000 -254380.926859
 4.230e-08    300.000 101325.000 -254380.938718
 4.240e-08    300.000 101325.000 -254380.950566
 4.250e-08    300.000 101325.000 -254380.962401
 4.260e-08    300.000 101325.000 -254380.974225
 4.270e-08    300.000 101325.000 -254380.986036
 4.280e-08    300.000 101325.000 -254380.997835
 4.290e-08    300.000 101325.000 -254381.009623
 4.300e-08    300.000 101325.000 -254381.021398
 4.310e-08    300.000 101325.000 -254381.033161
 4.320e-08    300.000 101325.000 -254381.044912
 4.330e-08    300.000 101325.000 -254381.056652
 4.340e-08    300.000 101325.000 -254381.068379
 4.350e-08    300.000 101325.000 -254381.080095
 4.360e-08    300.000 101325.000 -254381.091799
 4.370e-08    300.000 101325.000 -254381.103491
 4.380e-08    300.000 101325.000 -254381.115171
 4.390e-08    300.000 101325.000 -254381.126839
 4.400e-08    300.000 101325.000 -254381.138495
 4.410e-08    300.000 101325.000 -254381.150139
 4.420e-08    300.000 101325.000 -254381.161772
 4.430e-08    300.000 101325.000 -254381.173393
 4.440e-08    300.000 101325.000 -254381.185002
 4.450e-08    300.000 101325.000 -254381.196599
 4.460e-08    300.000 101325.000 -254381.208184
 4.470e-08    300.000 101325.000 -254381.219758
 4.480e-08    300.000 101325.000 -254381.231320
 4.490e-08    300.000 101325.000 -254381.242870
 4.500e-08    300.000 101325.000 -254381.254409
 4.510e-08    300.000 101325.000 -254381.265935
 4.520e-08    300.000 101325.000 -254381.277450
 4.530e-08    300.000 101325.000 -254381.288954
 4.540e-08    300.000 101325.000 -254381.300445
 4.550e-08    300.000 101325.000 -254381.311926
 4.560e-08    300.000 101325.000 -254381.323394
 4.570e-08    300.000 101325.000 -254381.334851
 4.580e-08    300.000 101325.000 -254381.346296
 4.590e-08    300.000 101325.000 -254381.357730
 4.600e-08    300.000 101325.000 -254381.369152
 4.610e-08    300.000 101325.000 -254381.380562
 4.620e-08    300.000 101325.000 -254381.391961
 4.630e-08    300.000 101325.000 -254381.403348
 4.640e-08    300.000 101325.000 -254381.414724
 4.650e-08    300.000 101325.000 -254381.426089
 4.660e-08    300.000 101325.000 -254381.437441
 4.670e-08    300.000 101325.000 -254381.448783
 4.680e-08    300.000 101325.000 -254381.460113
 4.690e-08    300.000 101325.000 -254381.471431
 4.700e-08    300.000 101325.000 -254381.482738
 4.710e-08    300.000 101325.000 -254381.494033
 4.720e-08    300.000 101325.000 -254381.505317
 4.730e-08    300.000 101325.000 -254381.516590
 4.740e-08    300.000 101325.000 -254381.527851
 4.750e-08    300.000 101325.000 -254381.539101
 4.760e-08    300.000 101325.000 -254381.550339
 4.770e-08    300.000 101325.000 -254381.561566
 4.780e-08    300.000 101325.000 -254381.572782
 4.790e-08    300.000 101325.000 -254381.583986
 4.800e-08    300.000 101325.000 -254381.595179
 4.810e-08    300.000 101325.000 -254381.606361
 4.820e-08    300.000 101325.000 -254381.617531
 4.830e-08    300.000 101325.000 -254381.628690
 4.840e-08    300.000 101325.000 -254381.639838
 4.850e-08    300.000 101325.000 -254381.650975
 4.860e-08    300.000 101325.000 -254381.662100
 4.870e-08    300.000 101325.000 -254381.673214
 4.880e-08    300.000 101325.000 -254381.684317
 4.890e-08    300.000 101325.000 -254381.695409
 4.900e-08    300.000 101325.000 -254381.706489
 4.910e-08    300.000 101325.000 -254381.717558
 4.920e-08    300.000 101325.000 -254381.728616
 4.930e-08    300.000 101325.000 -254381.739663
 4.940e-08    300.000 101325.000 -254381.750699
 4.950e-08    300.000 101325.000 -254381.761723
 4.960e-08    300.000 101325.000 -254381.772737
 4.970e-08    300.000 101325.000 -254381.783739
 4.980e-08    300.000 101325.000 -254381.794730
 4.990e-08    300.000 101325.000 -254381.805711
 5.000e-08    300.000 101325.000 -254381.816680
 5.010e-08    300.000 101325.000 -254381.827637
 5.020e-08    300.000 101325.000 -254381.838584
 5.030e-08    300.000 101325.000 -254381.849520
 5.040e-08    300.000 101325.000 -254381.860445
 5.050e-08    300.000 101325.000 -254381.871359
 5.060e-08    300.000 101325.000 -254381.882262
 5.070e-08    300.000 101325.000 -254381.893153
 5.080e-08    300.000 101325.000 -254381.904034
 5.090e-08    300.000 101325.000 -254381.914904
 5.100e-08    300.000 101325.000 -254381.925763
 5.110e-08    300.000 101325.000 -254381.936611
 5.120e-08    300.000 101325.000 -254381.947448
 5.130e-08    300.000 101325.000 -254381.958274
 5.140e-08    300.000 101325.000 -254381.969089
 5.150e-08    300.000 101325.000 -254381.979893
 5.160e-08    300.000 101325.000 -254381.990687
 5.170e-08    300.000 101325.000 -254382.001469
 5.180e-08    300.000 101325.000 -254382.012241
 5.190e-08    300.000 101325.000 -254382.023002
 5.200e-08    300.000 101325.000 -254382.033752
 5.210e-08    300.000 101325.000 -254382.044491
 5.220e-08    300.000 101325.000 -254382.055219
 5.230e-08    300.000 101325.000 -254382.065936
 5.240e-08    300.000 101325.000 -254382.076643
 5.250e-08    300.000 101325.000 -254382.087339
 5.260e-08    300.000 101325.000 -254382.098024
 5.270e-08    300.000 101325.000 -254382.108698
 5.280e-08    300.000 101325.000 -254382.119362
 5.290e-08    300.000 101325.000 -254382.130015
 5.300e-08    300.000 101325.000 -254382.140657
 5.310e-08    300.000 101325.000 -254382.151289
 5.320e-08    300.000 101325.000 -254382.161909
 5.330e-08    300.000 101325.000 -254382.172519
 5.340e-08    300.000 101325.000 -254382.183119
 5.350e-08    300.000 101325.000 -254382.193708
 5.360e-08    300.000 101325.000 -254382.204286
 5.370e-08    300.000 101325.000 -254382.214853
 5.380e-08    300.000 101325.000 -254382.225410
 5.390e-08    300.000 101325.000 -254382.235956
 5.400e-08    300.000 101325.000 -254382.246492
 5.410e-08    300.000 101325.000 -254382.257017
 5.420e-08    300.000 101325.000 -254382.267531
 5.430e-08    300.000 101325.000 -254382.278035
 5.440e-08    300.000 101325.000 -254382.288529
 5.450e-08    300.000 101325.000 -254382.299011
 5.460e-08    300.000 101325.000 -254382.309484
 5.470e-08    300.000 101325.000 -254382.319945
 5.480e-08    300.000 101325.000 -254382.330397
 5.490e-08    300.000 101325.000 -254382.340837
 5.500e-08    300.000 101325.000 -254382.351268
 5.510e-08    300.000 101325.000 -254382.361687
 5.520e-08    300.000 101325.000 -254382.372096
 5.530e-08    300.000 101325.000 -254382.382495
 5.540e-08    300.000 101325.000 -254382.392884
 5.550e-08    300.000 101325.000 -254382.403262
 5.560e-08    300.000 101325.000 -254382.413629
 5.570e-08    300.000 101325.000 -254382.423986
 5.580e-08    300.000 101325.000 -254382.434333
 5.590e-08    300.000 101325.000 -254382.444669
 5.600e-08    300.000 101325.000 -254382.454995
 5.610e-08    300.000 101325.000 -254382.465311
 5.620e-08    300.000 101325.000 -254382.475616
 5.630e-08    300.000 101325.000 -254382.485911
 5.640e-08    300.000 101325.000 -254382.496195
 5.650e-08    300.000 101325.000 -254382.506470
 5.660e-08    300.000 101325.000 -254382.516734
 5.670e-08    300.000 101325.000 -254382.526987
 5.680e-08    300.000 101325.000 -254382.537231
 5.690e-08    300.000 101325.000 -254382.547464
 5.700e-08    300.000 101325.000 -254382.557687
 5.710e-08    300.000 101325.000 -254382.567899
 5.720e-08    300.000 101325.000 -254382.578101
 5.730e-08    300.000 101325.000 -254382.588294
 5.740e-08    300.000 101325.000 -254382.598475
 5.750e-08    300.000 101325.000 -254382.608647
 5.760e-08    300.000 101325.000 -254382.618809
 5.770e-08    300.000 101325.000 -254382.628960
 5.780e-08    300.000 101325.000 -254382.639101
 5.790e-08    300.000 101325.000 -254382.649232
 5.800e-08    300.000 101325.000 -254382.659353
 5.810e-08    300.000 101325.000 -254382.669464
 5.820e-08    300.000 101325.000 -254382.679564
 5.830e-08    300.000 101325.000 -254382.689655
 5.840e-08    300.000 101325.000 -254382.699735
 5.850e-08    300.000 101325.000 -254382.709805
 5.860e-08    300.000 101325.000 -254382.719865
 5.870e-08    300.000 101325.000 -254382.729915
 5.880e-08    300.000 101325.000 -254382.739956
 5.890e-08    300.000 101325.000 -254382.749986
 5.900e-08    300.000 101325.000 -254382.760006
 5.910e-08    300.000 101325.000 -254382.770015
 5.920e-08    300.000 101325.000 -254382.780015
 5.930e-08    300.000 101325.000 -254382.790005
 5.940e-08    300.000 101325.000 -254382.799985
 5.950e-08    300.000 101325.000 -254382.809955
 5.960e-08    300.000 101325.000 -254382.819915
 5.970e-08    300.000 101325.000 -254382.829865
 5.980e-08    300.000 101325.000 -254382.839805
 5.990e-08    300.000 101325.000 -254382.849735
 6.000e-08    300.000 101325.000 -254382.859655
 6.010e-08    300.000 101325.000 -254382.869565
 6.020e-08    300.000 101325.000 -254382.879465
 6.030e-08    300.000 101325.000 -254382.889356
 6.040e-08    300.000 101325.000 -254382.899236
 6.050e-08    300.000 101325.000 -254382.909107
 6.060e-08    300.000 101325.000 -254382.918967
 6.070e-08    300.000 101325.000 -254382.928818
 6.080e-08    300.000 101325.000 -254382.938659
 6.090e-08    300.000 101325.000 -254382.948491
 6.100e-08    300.000 101325.000 -254382.958312
 6.110e-08    300.000 101325.000 -254382.968123
 6.120e-08    300.000 101325.000 -254382.977925
 6.130e-08    300.000 101325.000 -254382.987717
 6.140e-08    300.000 101325.000 -254382.997499
 6.150e-08    300.000 101325.000 -254383.007272
 6.160e-08    300.000 101325.000 -254383.017034
 6.170e-08    300.000 101325.000 -254383.026787
 6.180e-08    300.000 101325.000 -254383.036531
 6.190e-08    300.000 101325.000 -254383.046264
 6.200e-08    300.000 101325.000 -254383.055988
 6.210e-08    300.000 101325.000 -254383.065702
 6.220e-08    300.000 101325.000 -254383.075406
 6.230e-08    300.000 101325.000 -254383.085101
 6.240e-08    300.000 101325.000 -254383.094786
 6.250e-08    300.000 101325.000 -254383.104461
 6.260e-08    300.000 101325.000 -254383.114126
 6.270e-08    300.000 101325.000 -254383.123782
 6.280e-08    300.000 101325.000 -254383.133429
 6.290e-08    300.000 101325.000 -254383.143066
 6.300e-08    300.000 101325.000 -254383.152693
 6.310e-08    300.000 101325.000 -254383.162310
 6.320e-08    300.000 101325.000 -254383.171918
 6.330e-08    300.000 101325.000 -254383.181516
 6.340e-08    300.000 101325.000 -254383.191105
 6.350e-08    300.000 101325.000 -254383.200684
 6.360e-08    300.000 101325.000 -254383.210254
 6.370e-08    300.000 101325.000 -254383.219814
 6.380e-08    300.000 101325.000 -254383.229365
 6.390e-08    300.000 101325.000 -254383.238906
 6.400e-08    300.000 101325.000 -254383.248437
 6.410e-08    300.000 101325.000 -254383.257959
 6.420e-08    300.000 101325.000 -254383.267472
 6.430e-08    300.000 101325.000 -254383.276975
 6.440e-08    300.000 101325.000 -254383.286468
 6.450e-08    300.000 101325.000 -254383.295952
 6.460e-08    300.000 101325.000 -254383.305427
 6.470e-08    300.000 101325.000 -254383.314892
 6.480e-08    300.000 101325.000 -254383.324348
 6.490e-08    300.000 101325.000 -254383.333794
 6.500e-08    300.000 101325.000 -254383.343232
 6.510e-08    300.000 101325.000 -254383.352659
 6.520e-08    300.000 101325.000 -254383.362077
 6.530e-08    300.000 101325.000 -254383.371486
 6.540e-08    300.000 101325.000 -254383.380885
 6.550e-08    300.000 101325.000 -254383.390275
 6.560e-08    300.000 101325.000 -254383.399656
 6.570e-08    300.000 101325.000 -254383.409028
 6.580e-08    300.000 101325.000 -254383.418390
 6.590e-08    300.000 101325.000 -254383.427742
 6.600e-08    300.000 101325.000 -254383.437086
 6.610e-08    300.000 101325.000 -254383.446420
 6.620e-08    300.000 101325.000 -254383.455745
 6.630e-08    300.000 101325.000 -254383.465060
 6.640e-08    300.000 101325.000 -254383.474367
 6.650e-08    300.000 101325.000 -254383.483664
 6.660e-08    300.000 101325.000 -254383.492951
 6.670e-08    300.000 101325.000 -254383.502230
 6.680e-08    300.000 101325.000 -254383.511499
 6.690e-08    300.000 101325.000 -254383.520760
 6.700e-08    300.000 101325.000 -254383.530010
 6.710e-08    300.000 101325.000 -254383.539252
 6.720e-08    300.000 101325.000 -254383.548485
 6.730e-08    300.000 101325.000 -254383.557708
 6.740e-08    300.000 101325.000 -254383.566922
 6.750e-08    300.000 101325.000 -254383.576127
 6.760e-08    300.000 101325.000 -254383.585323
 6.770e-08    300.000 101325.000 -254383.594510
 6.780e-08    300.000 101325.000 -254383.603687
 6.790e-08    300.000 101325.000 -254383.612856
 6.800e-08    300.000 101325.000 -254383.622015
 6.810e-08    300.000 101325.000 -254383.631165
 6.820e-08    300.000 101325.000 -254383.640306
 6.830e-08    300.000 101325.000 -254383.649438
 6.840e-08    300.000 101325.000 -254383.658561
 6.850e-08    300.000 101325.000 -254383.667675
 6.860e-08    300.000 101325.000 -254383.676780
 6.870e-08    300.000 101325.000 -254383.685876
 6.880e-08    300.000 101325.000 -254383.694963
 6.890e-08    300.000 101325.000 -254383.704041
 6.900e-08    300.000 101325.000 -254383.713110
 6.910e-08    300.000 101325.000 -254383.722169
 6.920e-08    300.000 101325.000 -254383.731220
 6.930e-08    300.000 101325.000 -254383.740262
 6.940e-08    300.000 101325.000 -254383.749295
 6.950e-08    300.000 101325.000 -254383.758318
 6.960e-08    300.000 101325.000 -254383.767333
 6.970e-08    300.000 101325.000 -254383.776339
 6.980e-08    300.000 101325.000 -254383.785336
 6.990e-08    300.000 101325.000 -254383.794324
 7.000e-08    300.000 101325.000 -254383.803304
 7.010e-08    300.000 101325.000 -254383.812274
 7.020e-08    300.000 101325.000 -254383.821235
 7.030e-08    300.000 101325.000 -254383.830188
 7.040e-08    300.000 101325.000 -254383.839131
 7.050e-08    300.000 101325.000 -254383.848066
 7.060e-08    300.000 101325.000 -254383.856992
 7.070e-08    300.000 101325.000 -254383.865909
 7.080e-08    300.000 101325.000 -254383.874817
 7.090e-08    300.000 101325.000 -254383.883716
 7.100e-08    300.000 101325.000 -254383.892607
 7.110e-08    300.000 101325.000 -254383.901489
 7.120e-08    300.000 101325.000 -254383.910361
 7.130e-08    300.000 101325.000 -254383.919226
 7.140e-08    300.000 101325.000 -254383.928081
 7.150e-08    300.000 101325.000 -254383.936927
 7.160e-08    300.000 101325.000 -254383.945765
 7.170e-08    300.000 101325.000 -254383.954594
 7.180e-08    300.000 101325.000 -254383.963415
 7.190e-08    300.000 101325.000 -254383.972226
 7.200e-08    300.000 101325.000 -254383.981029
 7.210e-08    300.000 101325.000 -254383.989823
 7.220e-08    300.000 101325.000 -254383.998608
 7.230e-08    300.000 101325.000 -254384.007385
 7.240e-08    300.000 101325.000 -254384.016153
 7.250e-08    300.000 101325.000 -254384.024912
 7.260e-08    300.000 101325.000 -254384.033663
 7.270e-08    300.000 101325.000 -254384.042405
 7.280e-08    300.000 101325.000 -254384.051138
 7.290e-08    300.000 101325.000 -254384.059863
 7.300e-08    300.000 101325.000 -254384.068579
 7.310e-08    300.000 101325.000 -254384.077286
 7.320e-08    300.000 101325.000 -254384.085985
 7.330e-08    300.000 101325.000 -254384.094675
 7.340e-08    300.000 101325.000 -254384.103357
 7.350e-08    300.000 101325.000 -254384.112030
 7.360e-08    300.000 101325.000 -254384.120694
 7.370e-08    300.000 101325.000 -254384.129350
 7.380e-08    300.000 101325.000 -254384.137997
 7.390e-08    300.000 101325.000 -254384.146636
 7.400e-08    300.000 101325.000 -254384.155266
 7.410e-08    300.000 101325.000 -254384.163888
 7.420e-08    300.000 101325.000 -254384.172501
 7.430e-08    300.000 101325.000 -254384.181106
 7.440e-08    300.000 101325.000 -254384.189702
 7.450e-08    300.000 101325.000 -254384.198289
 7.460e-08    300.000 101325.000 -254384.206868
 7.470e-08    300.000 101325.000 -254384.215439
 7.480e-08    300.000 101325.000 -254384.224001
 7.490e-08    300.000 101325.000 -254384.232555
 7.500e-08    300.000 101325.000 -254384.241100
 7.510e-08    300.000 101325.000 -254384.249637
 7.520e-08    300.000 101325.000 -254384.258165
 7.530e-08    300.000 101325.000 -254384.266685
 7.540e-08    300.000 101325.000 -254384.275197
 7.550e-08    300.000 101325.000 -254384.283700
 7.560e-08    300.000 101325.000 -254384.292194
 7.570e-08    300.000 101325.000 -254384.300681
 7.580e-08    300.000 101325.000 -254384.309159
 7.590e-08    300.000 101325.000 -254384.317628
 7.600e-08    300.000 101325.000 -254384.326089
 7.610e-08    300.000 101325.000 -254384.334542
 7.620e-08    300.000 101325.000 -254384.342987
 7.630e-08    300.000 101325.000 -254384.351423
 7.640e-08    300.000 101325.000 -254384.359851
 7.650e-08    300.000 101325.000 -254384.368270
 7.660e-08    300.000 101325.000 -254384.376681
 7.670e-08    300.000 101325.000 -254384.385084
 7.680e-08    300.000 101325.000 -254384.393479
 7.690e-08    300.000 101325.000 -254384.401865
 7.700e-08    300.000 101325.000 -254384.410243
 7.710e-08    300.000 101325.000 -254384.418613
 7.720e-08    300.000 101325.000 -254384.426974
 7.730e-08    300.000 101325.000 -254384.435328
 7.740e-08    300.000 101325.000 -254384.443673
 7.750e-08    300.000 101325.000 -254384.452009
 7.760e-08    300.000 101325.000 -254384.460338
 7.770e-08    300.000 101325.000 -254384.468658
 7.780e-08    300.000 101325.000 -254384.476970
 7.790e-08    300.000 101325.000 -254384.485274
 7.800e-08    300.000 101325.000 -254384.493570
 7.810e-08    300.000 101325.000 -254384.501858
 7.820e-08    300.000 101325.000 -254384.510137
 7.830e-08    300.000 101325.000 -254384.518408
 7.840e-08    300.000 101325.000 -254384.526671
 7.850e-08    300.000 101325.000 -254384.534926
 7.860e-08    300.000 101325.000 -254384.543173
 7.870e-08    300.000 101325.000 -254384.551412
 7.880e-08    300.000 101325.000 -254384.559642
 7.890e-08    300.000 101325.000 -254384.567865
 7.900e-08    300.000 101325.000 -254384.576079
 7.910e-08    300.000 101325.000 -254384.584285
 7.920e-08    300.000 101325.000 -254384.592483
 7.930e-08    300.000 101325.000 -254384.600673
 7.940e-08    300.000 101325.000 -254384.608855
 7.950e-08    300.000 101325.000 -254384.617029
 7.960e-08    300.000 101325.000 -254384.625195
 7.970e-08    300.000 101325.000 -254384.633353
 7.980e-08    300.000 101325.000 -254384.641503
 7.990e-08    300.000 101325.000 -254384.649645
 8.000e-08    300.000 101325.000 -254384.657779
 8.010e-08    300.000 101325.000 -254384.665905
 8.020e-08    300.000 101325.000 -254384.674022
 8.030e-08    300.000 101325.000 -254384.682132
 8.040e-08    300.000 101325.000 -254384.690234
 8.050e-08    300.000 101325.000 -254384.698328
 8.060e-08    300.000 101325.000 -254384.706413
 8.070e-08    300.000 101325.000 -254384.714491
 8.080e-08    300.000 101325.000 -254384.722561
 8.090e-08    300.000 101325.000 -254384.730623
 8.100e-08    300.000 101325.000 -254384.738677
 8.110e-08    300.000 101325.000 -254384.746724
 8.120e-08    300.000 101325.000 -254384.754762
 8.130e-08    300.000 101325.000 -254384.762792
 8.140e-08    300.000 101325.000 -254384.770814
 8.150e-08    300.000 101325.000 -254384.778829
 8.160e-08    300.000 101325.000 -254384.786836
 8.170e-08    300.000 101325.000 -254384.794834
 8.180e-08    300.000 101325.000 -254384.802825
 8.190e-08    300.000 101325.000 -254384.810809
 8.200e-08    300.000 101325.000 -254384.818784
 8.210e-08    300.000 101325.000 -254384.826751
 8.220e-08    300.000 101325.000 -254384.834711
 8.230e-08    300.000 101325.000 -254384.842662
 8.240e-08    300.000 101325.000 -254384.850606
 8.250e-08    300.000 101325.000 -254384.858542
 8.260e-08    300.000 101325.000 -254384.866470
 8.270e-08    300.000 101325.000 -254384.874391
 8.280e-08    300.000 101325.000 -254384.882304
 8.290e-08    300.000 101325.000 -254384.890209
 8.300e-08    300.000 101325.000 -254384.898106
 8.310e-08    300.000 101325.000 -254384.905995
 8.320e-08    300.000 101325.000 -254384.913877
 8.330e-08    300.000 101325.000 -254384.921751
 8.340e-08    300.000 101325.000 -254384.929617
 8.350e-08    300.000 101325.000 -254384.937475
 8.360e-08    300.000 101325.000 -254384.945326
 8.370e-08    300.000 101325.000 -254384.953169
 8.380e-08    300.000 101325.000 -254384.961005
 8.390e-08    300.000 101325.000 -254384.968832
 8.400e-08    300.000 101325.000 -254384.976652
 8.410e-08    300.000 101325.000 -254384.984465
 8.420e-08    300.000 101325.000 -254384.992269
 8.430e-08    300.000 101325.000 -254385.000066
 8.440e-08    300.000 101325.000 -254385.007855
 8.450e-08    300.000 101325.000 -254385.015637
 8.460e-08    300.000 101325.000 -254385.023411
 8.470e-08    300.000 101325.000 -254385.031177
 8.480e-08    300.000 101325.000 -254385.038936
 8.490e-08    300.000 101325.000 -254385.046687
 8.500e-08    300.000 101325.000 -254385.054431
 8.510e-08    300.000 101325.000 -254385.062167
 8.520e-08    300.000 101325.000 -254385.069895
 8.530e-08    300.000 101325.000 -254385.077616
 8.540e-08    300.000 101325.000 -254385.085329
 8.550e-08    300.000 101325.000 -254385.093035
 8.560e-08    300.000 101325.000 -254385.100733
 8.570e-08    300.000 101325.000 -254385.108423
 8.580e-08    300.000 101325.000 -254385.116106
 8.590e-08    300.000 101325.000 -254385.123782
 8.600e-08    300.000 101325.000 -254385.131450
 8.610e-08    300.000 101325.000 -254385.139110
 8.620e-08    300.000 101325.000 -254385.146763
 8.630e-08    300.000 101325.000 -254385.154408
 8.640e-08    300.000 101325.000 -254385.162046
 8.650e-08    300.000 101325.000 -254385.169677
 8.660e-08    300.000 101325.000 -254385.177300
 8.670e-08    300.000 101325.000 -254385.184915
 8.680e-08    300.000 101325.000 -254385.192523
 8.690e-08    300.000 101325.000 -254385.200124
 8.700e-08    300.000 101325.000 -254385.207717
 8.710e-08    300.000 101325.000 -254385.215303
 8.720e-08    300.000 101325.000 -254385.222881
 8.730e-08    300.000 101325.000 -254385.230452
 8.740e-08    300.000 101325.000 -254385.238015
 8.750e-08    300.000 101325.000 -254385.245571
 8.760e-08    300.000 101325.000 -254385.253120
 8.770e-08    300.000 101325.000 -254385.260661
 8.780e-08    300.000 101325.000 -254385.268195
 8.790e-08    300.000 101325.000 -254385.275721
 8.800e-08    300.000 101325.000 -254385.283241
 8.810e-08    300.000 101325.000 -254385.290752
 8.820e-08    300.000 101325.000 -254385.298257
 8.830e-08    300.000 101325.000 -254385.305754
 8.840e-08    300.000 101325.000 -254385.313243
 8.850e-08    300.000 101325.000 -254385.320726
 8.860e-08    300.000 101325.000 -254385.328201
 8.870e-08    300.000 101325.000 -254385.335668
 8.880e-08    300.000 101325.000 -254385.343129
 8.890e-08    300.000 101325.000 -254385.350582
 8.900e-08    300.000 101325.000 -254385.358028
 8.910e-08    300.000 101325.000 -254385.365467
 8.920e-08    300.000 101325.000 -254385.372898
 8.930e-08    300.000 101325.000 -254385.380322
 8.940e-08    300.000 101325.000 -254385.387738
 8.950e-08    300.000 101325.000 -254385.395148
 8.960e-08    300.000 101325.000 -254385.402550
 8.970e-08    300.000 101325.000 -254385.409945
 8.980e-08    300.000 101325.000 -254385.417333
 8.990e-08    300.000 101325.000 -254385.424714
 9.000e-08    300.000 101325.000 -254385.432087
 9.010e-08    300.000 101325.000 -254385.439453

import cantera as ct
import numpy as np
import matplotlib.pyplot as plt

# Gaussian pulse parameters
EN_peak = 190 * 1e-21  # 190 Td
pulse_center = 24e-9  # 24 ns
pulse_width = 3e-9  # standard deviation (3 ns)
pulse_fwhm = pulse_width * 2 * (2 * np.log(2))**.5
gaussian_EN = ct.Func1("Gaussian", [EN_peak, pulse_center, pulse_fwhm])

# setup
gas = ct.Solution('example_data/methane-plasma-pavan-2023.yaml')
gas.TPX = 300., 101325., 'CH4:0.095, O2:0.19, N2:0.715, e:1E-11'
gas.reduced_electric_field = gaussian_EN(0)
gas.update_electron_energy_distribution()

r = ct.ConstPressureReactor(gas, energy="off")

sim = ct.ReactorNet([r])
sim.verbose = False

# simulation parameters
t_total = 90e-9
dt_max = 1e-10
dt_chunk = 1e-9  # 1 ns chunk
states = ct.SolutionArray(gas, extra=['t'])

print(f"{'t [s]':>10} {'T [K]':>10} {'P [Pa]':>10} {'h [J/kg]':>14}")

# simulate in 1 ns chunks
t = 0.0
while t < t_total:

    # integrate over the next chunk
    t_end = min(t + dt_chunk, t_total)
    while sim.time < t_end:
        sim.advance(sim.time + dt_max) #use sim.step
        states.append(r.thermo.state, t=sim.time)
        print(f"{sim.time:10.3e} {r.T:10.3f} {r.thermo.P:10.3f} {r.thermo.h:14.6f}")

    EN_t = gaussian_EN(t)
    gas.reduced_electric_field = EN_t
    gas.update_electron_energy_distribution()

    # reinitialize integrator with new source terms
    sim.reinitialize()

    t = t_end

# Plotting
fig, ax = plt.subplots(2, layout="constrained")

ax[0].plot(states.t, states.X[:, gas.species_index('e')], label='e')
ax[0].plot(states.t, states.X[:, gas.species_index('O2+')], label='O2+')
ax[0].plot(states.t, states.X[:, gas.species_index('N2+')], label='N2+')
ax[0].plot(states.t, states.X[:, gas.species_index('H2O+')], label='H2O+')
ax[0].plot(states.t, states.X[:, gas.species_index('CH4+')], label='CH4+')
ax[0].plot(states.t, states.X[:, gas.species_index('O')], label='O')
ax[0].plot(states.t, states.X[:, gas.species_index('N2(A)')], label='N2(A)')
ax[0].plot(states.t, states.X[:, gas.species_index('N2(B)')], label='N2(B)')
ax[0].plot(states.t, states.X[:, gas.species_index('N2(C)')], label='N2(C)')
ax[0].plot(states.t, states.X[:, gas.species_index("N2(a')")], label="N2(a')")
ax[0].plot(states.t, states.X[:, gas.species_index('CH3')], label='CH3', linestyle='--')
ax[0].plot(states.t, states.X[:, gas.species_index('CO2')], label='CO2', linestyle='--')
ax[0].plot(states.t, states.X[:, gas.species_index('CO')], label='CO', linestyle='--')
ax[0].plot(states.t, states.X[:, gas.species_index('H2O')], label='H2O', linestyle='--')
ax[0].plot(states.t, states.X[:, gas.species_index('H')], label='H', linestyle='--')
ax[0].plot(states.t, states.X[:, gas.species_index('OH')], label='OH', linestyle='--')
# N2 vibrational states
""" ax[0].plot(states.t, states.X[:, gas.species_index('N2(v1)')], label='N2(v1)')
ax[0].plot(states.t, states.X[:, gas.species_index('N2(v2)')], label='N2(v2)')
ax[0].plot(states.t, states.X[:, gas.species_index('N2(v3)')], label='N2(v3)')
ax[0].plot(states.t, states.X[:, gas.species_index('N2(v4)')], label='N2(v4)')
ax[0].plot(states.t, states.X[:, gas.species_index('N2(v5)')], label='N2(v5)')
ax[0].plot(states.t, states.X[:, gas.species_index('N2(v6)')], label='N2(v6)')
ax[0].plot(states.t, states.X[:, gas.species_index('N2(v7)')], label='N2(v7)')
ax[0].plot(states.t, states.X[:, gas.species_index('N2(v8)')], label='N2(v8)') """

ax[0].set_yscale('log')
ax[0].set_ylim([1e-14, 1e-3])

ax[1].plot(states.t, states.T, label='T')
ax2 = ax[1].twinx()
EN_values = [gaussian_EN(t) for t in states.t]
ax2.plot(states.t, EN_values, label='E/N', color='tab:red', linestyle='--')
ax2.set_ylabel('E/N', color='tab:red')
ax2.tick_params(axis='y', labelcolor='tab:red')

for axx in ax:
    axx.legend(loc='lower right', ncol=2)
    axx.set_xlabel('Time [s]')

ax[0].set_ylabel('Mole fraction [-]')
ax[1].set_ylabel('Temperature [K]')

plt.show()

Total running time of the script: (0 minutes 3.004 seconds)

Gallery generated by Sphinx-Gallery