Note
Go to the end to download the full example code.
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
t [s] T [K] P [Pa] h [J/kg]
1.000e-10 300.000 101325.000 -254430.994576
2.000e-10 300.000 101325.000 -254430.994576
3.000e-10 300.000 101325.000 -254430.994576
4.000e-10 300.000 101325.000 -254430.994576
5.000e-10 300.000 101325.000 -254430.994576
6.000e-10 300.000 101325.000 -254430.994576
7.000e-10 300.000 101325.000 -254430.994576
8.000e-10 300.000 101325.000 -254430.994576
9.000e-10 300.000 101325.000 -254430.994576
1.000e-09 300.000 101325.000 -254430.994576
1.100e-09 300.000 101325.000 -254430.994576
1.200e-09 300.000 101325.000 -254430.994576
1.300e-09 300.000 101325.000 -254430.994576
1.400e-09 300.000 101325.000 -254430.994576
1.500e-09 300.000 101325.000 -254430.994576
1.600e-09 300.000 101325.000 -254430.994576
1.700e-09 300.000 101325.000 -254430.994576
1.800e-09 300.000 101325.000 -254430.994576
1.900e-09 300.000 101325.000 -254430.994576
2.000e-09 300.000 101325.000 -254430.994576
2.100e-09 300.000 101325.000 -254430.994576
2.200e-09 300.000 101325.000 -254430.994576
2.300e-09 300.000 101325.000 -254430.994576
2.400e-09 300.000 101325.000 -254430.994576
2.500e-09 300.000 101325.000 -254430.994576
2.600e-09 300.000 101325.000 -254430.994576
2.700e-09 300.000 101325.000 -254430.994576
2.800e-09 300.000 101325.000 -254430.994576
2.900e-09 300.000 101325.000 -254430.994576
3.000e-09 300.000 101325.000 -254430.994576
3.100e-09 300.000 101325.000 -254430.994576
3.200e-09 300.000 101325.000 -254430.994576
3.300e-09 300.000 101325.000 -254430.994576
3.400e-09 300.000 101325.000 -254430.994576
3.500e-09 300.000 101325.000 -254430.994576
3.600e-09 300.000 101325.000 -254430.994576
3.700e-09 300.000 101325.000 -254430.994576
3.800e-09 300.000 101325.000 -254430.994576
3.900e-09 300.000 101325.000 -254430.994576
4.000e-09 300.000 101325.000 -254430.994576
4.100e-09 300.000 101325.000 -254430.994576
4.200e-09 300.000 101325.000 -254430.994576
4.300e-09 300.000 101325.000 -254430.994576
4.400e-09 300.000 101325.000 -254430.994576
4.500e-09 300.000 101325.000 -254430.994576
4.600e-09 300.000 101325.000 -254430.994576
4.700e-09 300.000 101325.000 -254430.994576
4.800e-09 300.000 101325.000 -254430.994576
4.900e-09 300.000 101325.000 -254430.994576
5.000e-09 300.000 101325.000 -254430.994576
5.100e-09 300.000 101325.000 -254430.994576
5.200e-09 300.000 101325.000 -254430.994576
5.300e-09 300.000 101325.000 -254430.994576
5.400e-09 300.000 101325.000 -254430.994576
5.500e-09 300.000 101325.000 -254430.994576
5.600e-09 300.000 101325.000 -254430.994576
5.700e-09 300.000 101325.000 -254430.994576
5.800e-09 300.000 101325.000 -254430.994576
5.900e-09 300.000 101325.000 -254430.994576
6.000e-09 300.000 101325.000 -254430.994576
6.100e-09 300.000 101325.000 -254430.994576
6.200e-09 300.000 101325.000 -254430.994576
6.300e-09 300.000 101325.000 -254430.994576
6.400e-09 300.000 101325.000 -254430.994576
6.500e-09 300.000 101325.000 -254430.994576
6.600e-09 300.000 101325.000 -254430.994576
6.700e-09 300.000 101325.000 -254430.994576
6.800e-09 300.000 101325.000 -254430.994576
6.900e-09 300.000 101325.000 -254430.994576
7.000e-09 300.000 101325.000 -254430.994576
7.100e-09 300.000 101325.000 -254430.994576
7.200e-09 300.000 101325.000 -254430.994576
7.300e-09 300.000 101325.000 -254430.994576
7.400e-09 300.000 101325.000 -254430.994576
7.500e-09 300.000 101325.000 -254430.994576
7.600e-09 300.000 101325.000 -254430.994576
7.700e-09 300.000 101325.000 -254430.994576
7.800e-09 300.000 101325.000 -254430.994576
7.900e-09 300.000 101325.000 -254430.994576
8.000e-09 300.000 101325.000 -254430.994576
8.100e-09 300.000 101325.000 -254430.994576
8.200e-09 300.000 101325.000 -254430.994576
8.300e-09 300.000 101325.000 -254430.994576
8.400e-09 300.000 101325.000 -254430.994576
8.500e-09 300.000 101325.000 -254430.994576
8.600e-09 300.000 101325.000 -254430.994576
8.700e-09 300.000 101325.000 -254430.994576
8.800e-09 300.000 101325.000 -254430.994576
8.900e-09 300.000 101325.000 -254430.994576
9.000e-09 300.000 101325.000 -254430.994576
9.100e-09 300.000 101325.000 -254430.994576
9.200e-09 300.000 101325.000 -254430.994576
9.300e-09 300.000 101325.000 -254430.994576
9.400e-09 300.000 101325.000 -254430.994576
9.500e-09 300.000 101325.000 -254430.994576
9.600e-09 300.000 101325.000 -254430.994576
9.700e-09 300.000 101325.000 -254430.994576
9.800e-09 300.000 101325.000 -254430.994576
9.900e-09 300.000 101325.000 -254430.994576
1.000e-08 300.000 101325.000 -254430.994576
1.010e-08 300.000 101325.000 -254430.994576
1.020e-08 300.000 101325.000 -254430.994576
1.030e-08 300.000 101325.000 -254430.994576
1.040e-08 300.000 101325.000 -254430.994576
1.050e-08 300.000 101325.000 -254430.994576
1.060e-08 300.000 101325.000 -254430.994576
1.070e-08 300.000 101325.000 -254430.994576
1.080e-08 300.000 101325.000 -254430.994576
1.090e-08 300.000 101325.000 -254430.994576
1.100e-08 300.000 101325.000 -254430.994576
1.110e-08 300.000 101325.000 -254430.994576
1.120e-08 300.000 101325.000 -254430.994576
1.130e-08 300.000 101325.000 -254430.994576
1.140e-08 300.000 101325.000 -254430.994576
1.150e-08 300.000 101325.000 -254430.994576
1.160e-08 300.000 101325.000 -254430.994576
1.170e-08 300.000 101325.000 -254430.994576
1.180e-08 300.000 101325.000 -254430.994576
1.190e-08 300.000 101325.000 -254430.994576
1.200e-08 300.000 101325.000 -254430.994576
1.210e-08 300.000 101325.000 -254430.994576
1.220e-08 300.000 101325.000 -254430.994576
1.230e-08 300.000 101325.000 -254430.994576
1.240e-08 300.000 101325.000 -254430.994576
1.250e-08 300.000 101325.000 -254430.994576
1.260e-08 300.000 101325.000 -254430.994576
1.270e-08 300.000 101325.000 -254430.994576
1.280e-08 300.000 101325.000 -254430.994576
1.290e-08 300.000 101325.000 -254430.994576
1.300e-08 300.000 101325.000 -254430.994576
1.310e-08 300.000 101325.000 -254430.994576
1.320e-08 300.000 101325.000 -254430.994576
1.330e-08 300.000 101325.000 -254430.994576
1.340e-08 300.000 101325.000 -254430.994576
1.350e-08 300.000 101325.000 -254430.994576
1.360e-08 300.000 101325.000 -254430.994576
1.370e-08 300.000 101325.000 -254430.994576
1.380e-08 300.000 101325.000 -254430.994576
1.390e-08 300.000 101325.000 -254430.994576
1.400e-08 300.000 101325.000 -254430.994576
1.410e-08 300.000 101325.000 -254430.994576
1.420e-08 300.000 101325.000 -254430.994576
1.430e-08 300.000 101325.000 -254430.994576
1.440e-08 300.000 101325.000 -254430.994576
1.450e-08 300.000 101325.000 -254430.994576
1.460e-08 300.000 101325.000 -254430.994576
1.470e-08 300.000 101325.000 -254430.994576
1.480e-08 300.000 101325.000 -254430.994575
1.490e-08 300.000 101325.000 -254430.994575
1.500e-08 300.000 101325.000 -254430.994575
1.510e-08 300.000 101325.000 -254430.994575
1.520e-08 300.000 101325.000 -254430.994575
1.530e-08 300.000 101325.000 -254430.994574
1.540e-08 300.000 101325.000 -254430.994574
1.550e-08 300.000 101325.000 -254430.994573
1.560e-08 300.000 101325.000 -254430.994572
1.570e-08 300.000 101325.000 -254430.994572
1.580e-08 300.000 101325.000 -254430.994571
1.590e-08 300.000 101325.000 -254430.994571
1.600e-08 300.000 101325.000 -254430.994570
1.610e-08 300.000 101325.000 -254430.994569
1.620e-08 300.000 101325.000 -254430.994564
1.630e-08 300.000 101325.000 -254430.994559
1.640e-08 300.000 101325.000 -254430.994554
1.650e-08 300.000 101325.000 -254430.994549
1.660e-08 300.000 101325.000 -254430.994544
1.670e-08 300.000 101325.000 -254430.994539
1.680e-08 300.000 101325.000 -254430.994534
1.690e-08 300.000 101325.000 -254430.994529
1.700e-08 300.000 101325.000 -254430.994524
1.710e-08 300.000 101325.000 -254430.994519
1.720e-08 300.000 101325.000 -254430.994485
1.730e-08 300.000 101325.000 -254430.994451
1.740e-08 300.000 101325.000 -254430.994418
1.750e-08 300.000 101325.000 -254430.994384
1.760e-08 300.000 101325.000 -254430.994351
1.770e-08 300.000 101325.000 -254430.994317
1.780e-08 300.000 101325.000 -254430.994283
1.790e-08 300.000 101325.000 -254430.994250
1.800e-08 300.000 101325.000 -254430.994216
1.810e-08 300.000 101325.000 -254430.994183
1.820e-08 300.000 101325.000 -254430.994012
1.830e-08 300.000 101325.000 -254430.993840
1.840e-08 300.000 101325.000 -254430.993668
1.850e-08 300.000 101325.000 -254430.993496
1.860e-08 300.000 101325.000 -254430.993324
1.870e-08 300.000 101325.000 -254430.993152
1.880e-08 300.000 101325.000 -254430.992980
1.890e-08 300.000 101325.000 -254430.992809
1.900e-08 300.000 101325.000 -254430.992637
1.910e-08 300.000 101325.000 -254430.992465
1.920e-08 300.000 101325.000 -254430.991924
1.930e-08 300.000 101325.000 -254430.991382
1.940e-08 300.000 101325.000 -254430.990841
1.950e-08 300.000 101325.000 -254430.990300
1.960e-08 300.000 101325.000 -254430.989758
1.970e-08 300.000 101325.000 -254430.989217
1.980e-08 300.000 101325.000 -254430.988676
1.990e-08 300.000 101325.000 -254430.988135
2.000e-08 300.000 101325.000 -254430.987593
2.010e-08 300.000 101325.000 -254430.987052
2.020e-08 300.000 101325.000 -254430.986194
2.030e-08 300.000 101325.000 -254430.985336
2.040e-08 300.000 101325.000 -254430.984478
2.050e-08 300.000 101325.000 -254430.983620
2.060e-08 300.000 101325.000 -254430.982762
2.070e-08 300.000 101325.000 -254430.981904
2.080e-08 300.000 101325.000 -254430.981046
2.090e-08 300.000 101325.000 -254430.980188
2.100e-08 300.000 101325.000 -254430.979330
2.110e-08 300.000 101325.000 -254430.978124
2.120e-08 300.000 101325.000 -254430.976918
2.130e-08 300.000 101325.000 -254430.975711
2.140e-08 300.000 101325.000 -254430.974505
2.150e-08 300.000 101325.000 -254430.973299
2.160e-08 300.000 101325.000 -254430.972093
2.170e-08 300.000 101325.000 -254430.970887
2.180e-08 300.000 101325.000 -254430.969680
2.190e-08 300.000 101325.000 -254430.968474
2.200e-08 300.000 101325.000 -254430.967268
2.210e-08 300.000 101325.000 -254430.965442
2.220e-08 300.000 101325.000 -254430.963616
2.230e-08 300.000 101325.000 -254430.961790
2.240e-08 300.000 101325.000 -254430.959965
2.250e-08 300.000 101325.000 -254430.958139
2.260e-08 300.000 101325.000 -254430.956313
2.270e-08 300.000 101325.000 -254430.954487
2.280e-08 300.000 101325.000 -254430.952661
2.290e-08 300.000 101325.000 -254430.950835
2.300e-08 300.000 101325.000 -254430.949010
2.310e-08 300.000 101325.000 -254430.946370
2.320e-08 300.000 101325.000 -254430.943730
2.330e-08 300.000 101325.000 -254430.941091
2.340e-08 300.000 101325.000 -254430.938451
2.350e-08 300.000 101325.000 -254430.935811
2.360e-08 300.000 101325.000 -254430.933172
2.370e-08 300.000 101325.000 -254430.930532
2.380e-08 300.000 101325.000 -254430.927892
2.390e-08 300.000 101325.000 -254430.925253
2.400e-08 300.000 101325.000 -254430.922613
2.410e-08 300.000 101325.000 -254430.919247
2.420e-08 300.000 101325.000 -254430.915881
2.430e-08 300.000 101325.000 -254430.912515
2.440e-08 300.000 101325.000 -254430.909149
2.450e-08 300.000 101325.000 -254430.905782
2.460e-08 300.000 101325.000 -254430.902416
2.470e-08 300.000 101325.000 -254430.899050
2.480e-08 300.000 101325.000 -254430.895684
2.490e-08 300.000 101325.000 -254430.892318
2.500e-08 300.000 101325.000 -254430.888952
2.510e-08 300.000 101325.000 -254430.885291
2.520e-08 300.000 101325.000 -254430.881630
2.530e-08 300.000 101325.000 -254430.877969
2.540e-08 300.000 101325.000 -254430.874308
2.550e-08 300.000 101325.000 -254430.870647
2.560e-08 300.000 101325.000 -254430.866986
2.570e-08 300.000 101325.000 -254430.863325
2.580e-08 300.000 101325.000 -254430.859664
2.590e-08 300.000 101325.000 -254430.856003
2.600e-08 300.000 101325.000 -254430.852342
2.610e-08 300.000 101325.000 -254430.848976
2.620e-08 300.000 101325.000 -254430.845610
2.630e-08 300.000 101325.000 -254430.842243
2.640e-08 300.000 101325.000 -254430.838877
2.650e-08 300.000 101325.000 -254430.835511
2.660e-08 300.000 101325.000 -254430.832145
2.670e-08 300.000 101325.000 -254430.828779
2.680e-08 300.000 101325.000 -254430.825413
2.690e-08 300.000 101325.000 -254430.822047
2.700e-08 300.000 101325.000 -254430.818681
2.710e-08 300.000 101325.000 -254430.816041
2.720e-08 300.000 101325.000 -254430.813401
2.730e-08 300.000 101325.000 -254430.810762
2.740e-08 300.000 101325.000 -254430.808122
2.750e-08 300.000 101325.000 -254430.805482
2.760e-08 300.000 101325.000 -254430.802843
2.770e-08 300.000 101325.000 -254430.800203
2.780e-08 300.000 101325.000 -254430.797563
2.790e-08 300.000 101325.000 -254430.794924
2.800e-08 300.000 101325.000 -254430.792284
2.810e-08 300.000 101325.000 -254430.790458
2.820e-08 300.000 101325.000 -254430.788632
2.830e-08 300.000 101325.000 -254430.786807
2.840e-08 300.000 101325.000 -254430.784981
2.850e-08 300.000 101325.000 -254430.783155
2.860e-08 300.000 101325.000 -254430.781329
2.870e-08 300.000 101325.000 -254430.779503
2.880e-08 300.000 101325.000 -254430.777677
2.890e-08 300.000 101325.000 -254430.775852
2.900e-08 300.000 101325.000 -254430.774026
2.910e-08 300.000 101325.000 -254430.772820
2.920e-08 300.000 101325.000 -254430.771613
2.930e-08 300.000 101325.000 -254430.770407
2.940e-08 300.000 101325.000 -254430.769201
2.950e-08 300.000 101325.000 -254430.767995
2.960e-08 300.000 101325.000 -254430.766788
2.970e-08 300.000 101325.000 -254430.765582
2.980e-08 300.000 101325.000 -254430.764376
2.990e-08 300.000 101325.000 -254430.763170
3.000e-08 300.000 101325.000 -254430.761964
3.010e-08 300.000 101325.000 -254430.761106
3.020e-08 300.000 101325.000 -254430.760248
3.030e-08 300.000 101325.000 -254430.759390
3.040e-08 300.000 101325.000 -254430.758532
3.050e-08 300.000 101325.000 -254430.757674
3.060e-08 300.000 101325.000 -254430.756816
3.070e-08 300.000 101325.000 -254430.755958
3.080e-08 300.000 101325.000 -254430.755100
3.090e-08 300.000 101325.000 -254430.754242
3.100e-08 300.000 101325.000 -254430.753384
3.110e-08 300.000 101325.000 -254430.752842
3.120e-08 300.000 101325.000 -254430.752301
3.130e-08 300.000 101325.000 -254430.751760
3.140e-08 300.000 101325.000 -254430.751219
3.150e-08 300.000 101325.000 -254430.750677
3.160e-08 300.000 101325.000 -254430.750136
3.170e-08 300.000 101325.000 -254430.749595
3.180e-08 300.000 101325.000 -254430.749054
3.190e-08 300.000 101325.000 -254430.748512
3.200e-08 300.000 101325.000 -254430.747971
3.210e-08 300.000 101325.000 -254430.747799
3.220e-08 300.000 101325.000 -254430.747627
3.230e-08 300.000 101325.000 -254430.747456
3.240e-08 300.000 101325.000 -254430.747284
3.250e-08 300.000 101325.000 -254430.747112
3.260e-08 300.000 101325.000 -254430.746940
3.270e-08 300.000 101325.000 -254430.746768
3.280e-08 300.000 101325.000 -254430.746596
3.290e-08 300.000 101325.000 -254430.746424
3.300e-08 300.000 101325.000 -254430.746253
3.310e-08 300.000 101325.000 -254430.746219
3.320e-08 300.000 101325.000 -254430.746185
3.330e-08 300.000 101325.000 -254430.746151
3.340e-08 300.000 101325.000 -254430.746118
3.350e-08 300.000 101325.000 -254430.746084
3.360e-08 300.000 101325.000 -254430.746050
3.370e-08 300.000 101325.000 -254430.746017
3.380e-08 300.000 101325.000 -254430.745983
3.390e-08 300.000 101325.000 -254430.745949
3.400e-08 300.000 101325.000 -254430.745916
3.410e-08 300.000 101325.000 -254430.745910
3.420e-08 300.000 101325.000 -254430.745905
3.430e-08 300.000 101325.000 -254430.745900
3.440e-08 300.000 101325.000 -254430.745895
3.450e-08 300.000 101325.000 -254430.745890
3.460e-08 300.000 101325.000 -254430.745885
3.470e-08 300.000 101325.000 -254430.745880
3.480e-08 300.000 101325.000 -254430.745875
3.490e-08 300.000 101325.000 -254430.745870
3.500e-08 300.000 101325.000 -254430.745864
3.510e-08 300.000 101325.000 -254430.745864
3.520e-08 300.000 101325.000 -254430.745863
3.530e-08 300.000 101325.000 -254430.745863
3.540e-08 300.000 101325.000 -254430.745862
3.550e-08 300.000 101325.000 -254430.745861
3.560e-08 300.000 101325.000 -254430.745861
3.570e-08 300.000 101325.000 -254430.745860
3.580e-08 300.000 101325.000 -254430.745860
3.590e-08 300.000 101325.000 -254430.745859
3.600e-08 300.000 101325.000 -254430.745858
3.610e-08 300.000 101325.000 -254430.745858
3.620e-08 300.000 101325.000 -254430.745858
3.630e-08 300.000 101325.000 -254430.745858
3.640e-08 300.000 101325.000 -254430.745858
3.650e-08 300.000 101325.000 -254430.745858
3.660e-08 300.000 101325.000 -254430.745858
3.670e-08 300.000 101325.000 -254430.745858
3.680e-08 300.000 101325.000 -254430.745858
3.690e-08 300.000 101325.000 -254430.745858
3.700e-08 300.000 101325.000 -254430.745858
3.710e-08 300.000 101325.000 -254430.745858
3.720e-08 300.000 101325.000 -254430.745858
3.730e-08 300.000 101325.000 -254430.745858
3.740e-08 300.000 101325.000 -254430.745858
3.750e-08 300.000 101325.000 -254430.745858
3.760e-08 300.000 101325.000 -254430.745858
3.770e-08 300.000 101325.000 -254430.745858
3.780e-08 300.000 101325.000 -254430.745858
3.790e-08 300.000 101325.000 -254430.745858
3.800e-08 300.000 101325.000 -254430.745858
3.810e-08 300.000 101325.000 -254430.745858
3.820e-08 300.000 101325.000 -254430.745858
3.830e-08 300.000 101325.000 -254430.745858
3.840e-08 300.000 101325.000 -254430.745858
3.850e-08 300.000 101325.000 -254430.745858
3.860e-08 300.000 101325.000 -254430.745858
3.870e-08 300.000 101325.000 -254430.745858
3.880e-08 300.000 101325.000 -254430.745858
3.890e-08 300.000 101325.000 -254430.745858
3.900e-08 300.000 101325.000 -254430.745858
3.910e-08 300.000 101325.000 -254430.745858
3.920e-08 300.000 101325.000 -254430.745858
3.930e-08 300.000 101325.000 -254430.745858
3.940e-08 300.000 101325.000 -254430.745858
3.950e-08 300.000 101325.000 -254430.745858
3.960e-08 300.000 101325.000 -254430.745858
3.970e-08 300.000 101325.000 -254430.745858
3.980e-08 300.000 101325.000 -254430.745858
3.990e-08 300.000 101325.000 -254430.745858
4.000e-08 300.000 101325.000 -254430.745858
4.010e-08 300.000 101325.000 -254430.745858
4.020e-08 300.000 101325.000 -254430.745858
4.030e-08 300.000 101325.000 -254430.745858
4.040e-08 300.000 101325.000 -254430.745858
4.050e-08 300.000 101325.000 -254430.745858
4.060e-08 300.000 101325.000 -254430.745858
4.070e-08 300.000 101325.000 -254430.745858
4.080e-08 300.000 101325.000 -254430.745858
4.090e-08 300.000 101325.000 -254430.745858
4.100e-08 300.000 101325.000 -254430.745858
4.110e-08 300.000 101325.000 -254430.745858
4.120e-08 300.000 101325.000 -254430.745858
4.130e-08 300.000 101325.000 -254430.745858
4.140e-08 300.000 101325.000 -254430.745858
4.150e-08 300.000 101325.000 -254430.745858
4.160e-08 300.000 101325.000 -254430.745858
4.170e-08 300.000 101325.000 -254430.745858
4.180e-08 300.000 101325.000 -254430.745858
4.190e-08 300.000 101325.000 -254430.745858
4.200e-08 300.000 101325.000 -254430.745858
4.210e-08 300.000 101325.000 -254430.745858
4.220e-08 300.000 101325.000 -254430.745858
4.230e-08 300.000 101325.000 -254430.745858
4.240e-08 300.000 101325.000 -254430.745858
4.250e-08 300.000 101325.000 -254430.745858
4.260e-08 300.000 101325.000 -254430.745858
4.270e-08 300.000 101325.000 -254430.745858
4.280e-08 300.000 101325.000 -254430.745858
4.290e-08 300.000 101325.000 -254430.745858
4.300e-08 300.000 101325.000 -254430.745858
4.310e-08 300.000 101325.000 -254430.745858
4.320e-08 300.000 101325.000 -254430.745858
4.330e-08 300.000 101325.000 -254430.745858
4.340e-08 300.000 101325.000 -254430.745858
4.350e-08 300.000 101325.000 -254430.745858
4.360e-08 300.000 101325.000 -254430.745858
4.370e-08 300.000 101325.000 -254430.745858
4.380e-08 300.000 101325.000 -254430.745858
4.390e-08 300.000 101325.000 -254430.745858
4.400e-08 300.000 101325.000 -254430.745858
4.410e-08 300.000 101325.000 -254430.745858
4.420e-08 300.000 101325.000 -254430.745858
4.430e-08 300.000 101325.000 -254430.745858
4.440e-08 300.000 101325.000 -254430.745858
4.450e-08 300.000 101325.000 -254430.745858
4.460e-08 300.000 101325.000 -254430.745858
4.470e-08 300.000 101325.000 -254430.745858
4.480e-08 300.000 101325.000 -254430.745858
4.490e-08 300.000 101325.000 -254430.745858
4.500e-08 300.000 101325.000 -254430.745858
4.510e-08 300.000 101325.000 -254430.745858
4.520e-08 300.000 101325.000 -254430.745858
4.530e-08 300.000 101325.000 -254430.745858
4.540e-08 300.000 101325.000 -254430.745858
4.550e-08 300.000 101325.000 -254430.745858
4.560e-08 300.000 101325.000 -254430.745858
4.570e-08 300.000 101325.000 -254430.745858
4.580e-08 300.000 101325.000 -254430.745858
4.590e-08 300.000 101325.000 -254430.745858
4.600e-08 300.000 101325.000 -254430.745858
4.610e-08 300.000 101325.000 -254430.745858
4.620e-08 300.000 101325.000 -254430.745858
4.630e-08 300.000 101325.000 -254430.745858
4.640e-08 300.000 101325.000 -254430.745858
4.650e-08 300.000 101325.000 -254430.745858
4.660e-08 300.000 101325.000 -254430.745858
4.670e-08 300.000 101325.000 -254430.745858
4.680e-08 300.000 101325.000 -254430.745858
4.690e-08 300.000 101325.000 -254430.745858
4.700e-08 300.000 101325.000 -254430.745858
4.710e-08 300.000 101325.000 -254430.745858
4.720e-08 300.000 101325.000 -254430.745858
4.730e-08 300.000 101325.000 -254430.745858
4.740e-08 300.000 101325.000 -254430.745858
4.750e-08 300.000 101325.000 -254430.745858
4.760e-08 300.000 101325.000 -254430.745858
4.770e-08 300.000 101325.000 -254430.745858
4.780e-08 300.000 101325.000 -254430.745858
4.790e-08 300.000 101325.000 -254430.745858
4.800e-08 300.000 101325.000 -254430.745858
4.810e-08 300.000 101325.000 -254430.745858
4.820e-08 300.000 101325.000 -254430.745858
4.830e-08 300.000 101325.000 -254430.745858
4.840e-08 300.000 101325.000 -254430.745858
4.850e-08 300.000 101325.000 -254430.745858
4.860e-08 300.000 101325.000 -254430.745858
4.870e-08 300.000 101325.000 -254430.745858
4.880e-08 300.000 101325.000 -254430.745858
4.890e-08 300.000 101325.000 -254430.745858
4.900e-08 300.000 101325.000 -254430.745858
4.910e-08 300.000 101325.000 -254430.745858
4.920e-08 300.000 101325.000 -254430.745858
4.930e-08 300.000 101325.000 -254430.745858
4.940e-08 300.000 101325.000 -254430.745858
4.950e-08 300.000 101325.000 -254430.745858
4.960e-08 300.000 101325.000 -254430.745858
4.970e-08 300.000 101325.000 -254430.745858
4.980e-08 300.000 101325.000 -254430.745858
4.990e-08 300.000 101325.000 -254430.745858
5.000e-08 300.000 101325.000 -254430.745858
5.010e-08 300.000 101325.000 -254430.745858
5.020e-08 300.000 101325.000 -254430.745858
5.030e-08 300.000 101325.000 -254430.745858
5.040e-08 300.000 101325.000 -254430.745858
5.050e-08 300.000 101325.000 -254430.745858
5.060e-08 300.000 101325.000 -254430.745858
5.070e-08 300.000 101325.000 -254430.745858
5.080e-08 300.000 101325.000 -254430.745858
5.090e-08 300.000 101325.000 -254430.745858
5.100e-08 300.000 101325.000 -254430.745858
5.110e-08 300.000 101325.000 -254430.745858
5.120e-08 300.000 101325.000 -254430.745858
5.130e-08 300.000 101325.000 -254430.745858
5.140e-08 300.000 101325.000 -254430.745858
5.150e-08 300.000 101325.000 -254430.745858
5.160e-08 300.000 101325.000 -254430.745858
5.170e-08 300.000 101325.000 -254430.745858
5.180e-08 300.000 101325.000 -254430.745858
5.190e-08 300.000 101325.000 -254430.745858
5.200e-08 300.000 101325.000 -254430.745858
5.210e-08 300.000 101325.000 -254430.745858
5.220e-08 300.000 101325.000 -254430.745858
5.230e-08 300.000 101325.000 -254430.745858
5.240e-08 300.000 101325.000 -254430.745858
5.250e-08 300.000 101325.000 -254430.745858
5.260e-08 300.000 101325.000 -254430.745858
5.270e-08 300.000 101325.000 -254430.745858
5.280e-08 300.000 101325.000 -254430.745858
5.290e-08 300.000 101325.000 -254430.745858
5.300e-08 300.000 101325.000 -254430.745858
5.310e-08 300.000 101325.000 -254430.745858
5.320e-08 300.000 101325.000 -254430.745858
5.330e-08 300.000 101325.000 -254430.745858
5.340e-08 300.000 101325.000 -254430.745858
5.350e-08 300.000 101325.000 -254430.745858
5.360e-08 300.000 101325.000 -254430.745858
5.370e-08 300.000 101325.000 -254430.745858
5.380e-08 300.000 101325.000 -254430.745858
5.390e-08 300.000 101325.000 -254430.745858
5.400e-08 300.000 101325.000 -254430.745858
5.410e-08 300.000 101325.000 -254430.745858
5.420e-08 300.000 101325.000 -254430.745858
5.430e-08 300.000 101325.000 -254430.745858
5.440e-08 300.000 101325.000 -254430.745858
5.450e-08 300.000 101325.000 -254430.745858
5.460e-08 300.000 101325.000 -254430.745858
5.470e-08 300.000 101325.000 -254430.745858
5.480e-08 300.000 101325.000 -254430.745858
5.490e-08 300.000 101325.000 -254430.745858
5.500e-08 300.000 101325.000 -254430.745858
5.510e-08 300.000 101325.000 -254430.745858
5.520e-08 300.000 101325.000 -254430.745858
5.530e-08 300.000 101325.000 -254430.745858
5.540e-08 300.000 101325.000 -254430.745858
5.550e-08 300.000 101325.000 -254430.745858
5.560e-08 300.000 101325.000 -254430.745858
5.570e-08 300.000 101325.000 -254430.745858
5.580e-08 300.000 101325.000 -254430.745858
5.590e-08 300.000 101325.000 -254430.745858
5.600e-08 300.000 101325.000 -254430.745858
5.610e-08 300.000 101325.000 -254430.745858
5.620e-08 300.000 101325.000 -254430.745858
5.630e-08 300.000 101325.000 -254430.745858
5.640e-08 300.000 101325.000 -254430.745858
5.650e-08 300.000 101325.000 -254430.745858
5.660e-08 300.000 101325.000 -254430.745858
5.670e-08 300.000 101325.000 -254430.745858
5.680e-08 300.000 101325.000 -254430.745858
5.690e-08 300.000 101325.000 -254430.745858
5.700e-08 300.000 101325.000 -254430.745858
5.710e-08 300.000 101325.000 -254430.745858
5.720e-08 300.000 101325.000 -254430.745858
5.730e-08 300.000 101325.000 -254430.745858
5.740e-08 300.000 101325.000 -254430.745858
5.750e-08 300.000 101325.000 -254430.745858
5.760e-08 300.000 101325.000 -254430.745858
5.770e-08 300.000 101325.000 -254430.745858
5.780e-08 300.000 101325.000 -254430.745858
5.790e-08 300.000 101325.000 -254430.745858
5.800e-08 300.000 101325.000 -254430.745858
5.810e-08 300.000 101325.000 -254430.745858
5.820e-08 300.000 101325.000 -254430.745858
5.830e-08 300.000 101325.000 -254430.745858
5.840e-08 300.000 101325.000 -254430.745858
5.850e-08 300.000 101325.000 -254430.745858
5.860e-08 300.000 101325.000 -254430.745858
5.870e-08 300.000 101325.000 -254430.745858
5.880e-08 300.000 101325.000 -254430.745858
5.890e-08 300.000 101325.000 -254430.745858
5.900e-08 300.000 101325.000 -254430.745858
5.910e-08 300.000 101325.000 -254430.745858
5.920e-08 300.000 101325.000 -254430.745858
5.930e-08 300.000 101325.000 -254430.745858
5.940e-08 300.000 101325.000 -254430.745858
5.950e-08 300.000 101325.000 -254430.745858
5.960e-08 300.000 101325.000 -254430.745858
5.970e-08 300.000 101325.000 -254430.745858
5.980e-08 300.000 101325.000 -254430.745858
5.990e-08 300.000 101325.000 -254430.745858
6.000e-08 300.000 101325.000 -254430.745858
6.010e-08 300.000 101325.000 -254430.745858
6.020e-08 300.000 101325.000 -254430.745858
6.030e-08 300.000 101325.000 -254430.745858
6.040e-08 300.000 101325.000 -254430.745858
6.050e-08 300.000 101325.000 -254430.745858
6.060e-08 300.000 101325.000 -254430.745858
6.070e-08 300.000 101325.000 -254430.745858
6.080e-08 300.000 101325.000 -254430.745858
6.090e-08 300.000 101325.000 -254430.745858
6.100e-08 300.000 101325.000 -254430.745858
6.110e-08 300.000 101325.000 -254430.745858
6.120e-08 300.000 101325.000 -254430.745858
6.130e-08 300.000 101325.000 -254430.745858
6.140e-08 300.000 101325.000 -254430.745858
6.150e-08 300.000 101325.000 -254430.745858
6.160e-08 300.000 101325.000 -254430.745858
6.170e-08 300.000 101325.000 -254430.745858
6.180e-08 300.000 101325.000 -254430.745858
6.190e-08 300.000 101325.000 -254430.745858
6.200e-08 300.000 101325.000 -254430.745858
6.210e-08 300.000 101325.000 -254430.745858
6.220e-08 300.000 101325.000 -254430.745858
6.230e-08 300.000 101325.000 -254430.745858
6.240e-08 300.000 101325.000 -254430.745858
6.250e-08 300.000 101325.000 -254430.745858
6.260e-08 300.000 101325.000 -254430.745858
6.270e-08 300.000 101325.000 -254430.745858
6.280e-08 300.000 101325.000 -254430.745858
6.290e-08 300.000 101325.000 -254430.745858
6.300e-08 300.000 101325.000 -254430.745858
6.310e-08 300.000 101325.000 -254430.745858
6.320e-08 300.000 101325.000 -254430.745858
6.330e-08 300.000 101325.000 -254430.745858
6.340e-08 300.000 101325.000 -254430.745858
6.350e-08 300.000 101325.000 -254430.745858
6.360e-08 300.000 101325.000 -254430.745858
6.370e-08 300.000 101325.000 -254430.745858
6.380e-08 300.000 101325.000 -254430.745858
6.390e-08 300.000 101325.000 -254430.745858
6.400e-08 300.000 101325.000 -254430.745858
6.410e-08 300.000 101325.000 -254430.745858
6.420e-08 300.000 101325.000 -254430.745858
6.430e-08 300.000 101325.000 -254430.745858
6.440e-08 300.000 101325.000 -254430.745858
6.450e-08 300.000 101325.000 -254430.745858
6.460e-08 300.000 101325.000 -254430.745858
6.470e-08 300.000 101325.000 -254430.745858
6.480e-08 300.000 101325.000 -254430.745858
6.490e-08 300.000 101325.000 -254430.745858
6.500e-08 300.000 101325.000 -254430.745858
6.510e-08 300.000 101325.000 -254430.745858
6.520e-08 300.000 101325.000 -254430.745858
6.530e-08 300.000 101325.000 -254430.745858
6.540e-08 300.000 101325.000 -254430.745858
6.550e-08 300.000 101325.000 -254430.745858
6.560e-08 300.000 101325.000 -254430.745858
6.570e-08 300.000 101325.000 -254430.745858
6.580e-08 300.000 101325.000 -254430.745858
6.590e-08 300.000 101325.000 -254430.745858
6.600e-08 300.000 101325.000 -254430.745858
6.610e-08 300.000 101325.000 -254430.745858
6.620e-08 300.000 101325.000 -254430.745858
6.630e-08 300.000 101325.000 -254430.745858
6.640e-08 300.000 101325.000 -254430.745858
6.650e-08 300.000 101325.000 -254430.745858
6.660e-08 300.000 101325.000 -254430.745858
6.670e-08 300.000 101325.000 -254430.745858
6.680e-08 300.000 101325.000 -254430.745858
6.690e-08 300.000 101325.000 -254430.745858
6.700e-08 300.000 101325.000 -254430.745858
6.710e-08 300.000 101325.000 -254430.745858
6.720e-08 300.000 101325.000 -254430.745858
6.730e-08 300.000 101325.000 -254430.745858
6.740e-08 300.000 101325.000 -254430.745858
6.750e-08 300.000 101325.000 -254430.745858
6.760e-08 300.000 101325.000 -254430.745858
6.770e-08 300.000 101325.000 -254430.745858
6.780e-08 300.000 101325.000 -254430.745858
6.790e-08 300.000 101325.000 -254430.745858
6.800e-08 300.000 101325.000 -254430.745858
6.810e-08 300.000 101325.000 -254430.745858
6.820e-08 300.000 101325.000 -254430.745858
6.830e-08 300.000 101325.000 -254430.745858
6.840e-08 300.000 101325.000 -254430.745858
6.850e-08 300.000 101325.000 -254430.745858
6.860e-08 300.000 101325.000 -254430.745858
6.870e-08 300.000 101325.000 -254430.745858
6.880e-08 300.000 101325.000 -254430.745858
6.890e-08 300.000 101325.000 -254430.745858
6.900e-08 300.000 101325.000 -254430.745858
6.910e-08 300.000 101325.000 -254430.745858
6.920e-08 300.000 101325.000 -254430.745858
6.930e-08 300.000 101325.000 -254430.745858
6.940e-08 300.000 101325.000 -254430.745858
6.950e-08 300.000 101325.000 -254430.745858
6.960e-08 300.000 101325.000 -254430.745858
6.970e-08 300.000 101325.000 -254430.745858
6.980e-08 300.000 101325.000 -254430.745858
6.990e-08 300.000 101325.000 -254430.745858
7.000e-08 300.000 101325.000 -254430.745858
7.010e-08 300.000 101325.000 -254430.745858
7.020e-08 300.000 101325.000 -254430.745858
7.030e-08 300.000 101325.000 -254430.745858
7.040e-08 300.000 101325.000 -254430.745858
7.050e-08 300.000 101325.000 -254430.745858
7.060e-08 300.000 101325.000 -254430.745858
7.070e-08 300.000 101325.000 -254430.745858
7.080e-08 300.000 101325.000 -254430.745858
7.090e-08 300.000 101325.000 -254430.745858
7.100e-08 300.000 101325.000 -254430.745858
7.110e-08 300.000 101325.000 -254430.745858
7.120e-08 300.000 101325.000 -254430.745858
7.130e-08 300.000 101325.000 -254430.745858
7.140e-08 300.000 101325.000 -254430.745858
7.150e-08 300.000 101325.000 -254430.745858
7.160e-08 300.000 101325.000 -254430.745858
7.170e-08 300.000 101325.000 -254430.745858
7.180e-08 300.000 101325.000 -254430.745858
7.190e-08 300.000 101325.000 -254430.745858
7.200e-08 300.000 101325.000 -254430.745858
7.210e-08 300.000 101325.000 -254430.745858
7.220e-08 300.000 101325.000 -254430.745858
7.230e-08 300.000 101325.000 -254430.745858
7.240e-08 300.000 101325.000 -254430.745858
7.250e-08 300.000 101325.000 -254430.745858
7.260e-08 300.000 101325.000 -254430.745858
7.270e-08 300.000 101325.000 -254430.745858
7.280e-08 300.000 101325.000 -254430.745858
7.290e-08 300.000 101325.000 -254430.745858
7.300e-08 300.000 101325.000 -254430.745858
7.310e-08 300.000 101325.000 -254430.745858
7.320e-08 300.000 101325.000 -254430.745858
7.330e-08 300.000 101325.000 -254430.745858
7.340e-08 300.000 101325.000 -254430.745858
7.350e-08 300.000 101325.000 -254430.745858
7.360e-08 300.000 101325.000 -254430.745858
7.370e-08 300.000 101325.000 -254430.745858
7.380e-08 300.000 101325.000 -254430.745858
7.390e-08 300.000 101325.000 -254430.745858
7.400e-08 300.000 101325.000 -254430.745858
7.410e-08 300.000 101325.000 -254430.745858
7.420e-08 300.000 101325.000 -254430.745858
7.430e-08 300.000 101325.000 -254430.745858
7.440e-08 300.000 101325.000 -254430.745858
7.450e-08 300.000 101325.000 -254430.745858
7.460e-08 300.000 101325.000 -254430.745858
7.470e-08 300.000 101325.000 -254430.745858
7.480e-08 300.000 101325.000 -254430.745858
7.490e-08 300.000 101325.000 -254430.745858
7.500e-08 300.000 101325.000 -254430.745858
7.510e-08 300.000 101325.000 -254430.745858
7.520e-08 300.000 101325.000 -254430.745858
7.530e-08 300.000 101325.000 -254430.745858
7.540e-08 300.000 101325.000 -254430.745858
7.550e-08 300.000 101325.000 -254430.745858
7.560e-08 300.000 101325.000 -254430.745858
7.570e-08 300.000 101325.000 -254430.745858
7.580e-08 300.000 101325.000 -254430.745858
7.590e-08 300.000 101325.000 -254430.745858
7.600e-08 300.000 101325.000 -254430.745858
7.610e-08 300.000 101325.000 -254430.745858
7.620e-08 300.000 101325.000 -254430.745858
7.630e-08 300.000 101325.000 -254430.745858
7.640e-08 300.000 101325.000 -254430.745858
7.650e-08 300.000 101325.000 -254430.745858
7.660e-08 300.000 101325.000 -254430.745858
7.670e-08 300.000 101325.000 -254430.745858
7.680e-08 300.000 101325.000 -254430.745858
7.690e-08 300.000 101325.000 -254430.745858
7.700e-08 300.000 101325.000 -254430.745858
7.710e-08 300.000 101325.000 -254430.745858
7.720e-08 300.000 101325.000 -254430.745858
7.730e-08 300.000 101325.000 -254430.745858
7.740e-08 300.000 101325.000 -254430.745858
7.750e-08 300.000 101325.000 -254430.745858
7.760e-08 300.000 101325.000 -254430.745858
7.770e-08 300.000 101325.000 -254430.745858
7.780e-08 300.000 101325.000 -254430.745858
7.790e-08 300.000 101325.000 -254430.745858
7.800e-08 300.000 101325.000 -254430.745858
7.810e-08 300.000 101325.000 -254430.745858
7.820e-08 300.000 101325.000 -254430.745858
7.830e-08 300.000 101325.000 -254430.745858
7.840e-08 300.000 101325.000 -254430.745858
7.850e-08 300.000 101325.000 -254430.745858
7.860e-08 300.000 101325.000 -254430.745858
7.870e-08 300.000 101325.000 -254430.745858
7.880e-08 300.000 101325.000 -254430.745858
7.890e-08 300.000 101325.000 -254430.745858
7.900e-08 300.000 101325.000 -254430.745858
7.910e-08 300.000 101325.000 -254430.745858
7.920e-08 300.000 101325.000 -254430.745858
7.930e-08 300.000 101325.000 -254430.745858
7.940e-08 300.000 101325.000 -254430.745858
7.950e-08 300.000 101325.000 -254430.745858
7.960e-08 300.000 101325.000 -254430.745858
7.970e-08 300.000 101325.000 -254430.745858
7.980e-08 300.000 101325.000 -254430.745858
7.990e-08 300.000 101325.000 -254430.745858
8.000e-08 300.000 101325.000 -254430.745858
8.010e-08 300.000 101325.000 -254430.745858
8.020e-08 300.000 101325.000 -254430.745858
8.030e-08 300.000 101325.000 -254430.745858
8.040e-08 300.000 101325.000 -254430.745858
8.050e-08 300.000 101325.000 -254430.745858
8.060e-08 300.000 101325.000 -254430.745858
8.070e-08 300.000 101325.000 -254430.745858
8.080e-08 300.000 101325.000 -254430.745858
8.090e-08 300.000 101325.000 -254430.745858
8.100e-08 300.000 101325.000 -254430.745858
8.110e-08 300.000 101325.000 -254430.745858
8.120e-08 300.000 101325.000 -254430.745858
8.130e-08 300.000 101325.000 -254430.745858
8.140e-08 300.000 101325.000 -254430.745858
8.150e-08 300.000 101325.000 -254430.745858
8.160e-08 300.000 101325.000 -254430.745858
8.170e-08 300.000 101325.000 -254430.745858
8.180e-08 300.000 101325.000 -254430.745858
8.190e-08 300.000 101325.000 -254430.745858
8.200e-08 300.000 101325.000 -254430.745858
8.210e-08 300.000 101325.000 -254430.745858
8.220e-08 300.000 101325.000 -254430.745858
8.230e-08 300.000 101325.000 -254430.745858
8.240e-08 300.000 101325.000 -254430.745858
8.250e-08 300.000 101325.000 -254430.745858
8.260e-08 300.000 101325.000 -254430.745858
8.270e-08 300.000 101325.000 -254430.745858
8.280e-08 300.000 101325.000 -254430.745858
8.290e-08 300.000 101325.000 -254430.745858
8.300e-08 300.000 101325.000 -254430.745858
8.310e-08 300.000 101325.000 -254430.745858
8.320e-08 300.000 101325.000 -254430.745858
8.330e-08 300.000 101325.000 -254430.745858
8.340e-08 300.000 101325.000 -254430.745858
8.350e-08 300.000 101325.000 -254430.745858
8.360e-08 300.000 101325.000 -254430.745858
8.370e-08 300.000 101325.000 -254430.745858
8.380e-08 300.000 101325.000 -254430.745858
8.390e-08 300.000 101325.000 -254430.745858
8.400e-08 300.000 101325.000 -254430.745858
8.410e-08 300.000 101325.000 -254430.745858
8.420e-08 300.000 101325.000 -254430.745858
8.430e-08 300.000 101325.000 -254430.745858
8.440e-08 300.000 101325.000 -254430.745858
8.450e-08 300.000 101325.000 -254430.745858
8.460e-08 300.000 101325.000 -254430.745858
8.470e-08 300.000 101325.000 -254430.745858
8.480e-08 300.000 101325.000 -254430.745858
8.490e-08 300.000 101325.000 -254430.745858
8.500e-08 300.000 101325.000 -254430.745858
8.510e-08 300.000 101325.000 -254430.745858
8.520e-08 300.000 101325.000 -254430.745858
8.530e-08 300.000 101325.000 -254430.745858
8.540e-08 300.000 101325.000 -254430.745858
8.550e-08 300.000 101325.000 -254430.745858
8.560e-08 300.000 101325.000 -254430.745858
8.570e-08 300.000 101325.000 -254430.745858
8.580e-08 300.000 101325.000 -254430.745858
8.590e-08 300.000 101325.000 -254430.745858
8.600e-08 300.000 101325.000 -254430.745858
8.610e-08 300.000 101325.000 -254430.745858
8.620e-08 300.000 101325.000 -254430.745858
8.630e-08 300.000 101325.000 -254430.745858
8.640e-08 300.000 101325.000 -254430.745858
8.650e-08 300.000 101325.000 -254430.745858
8.660e-08 300.000 101325.000 -254430.745858
8.670e-08 300.000 101325.000 -254430.745858
8.680e-08 300.000 101325.000 -254430.745858
8.690e-08 300.000 101325.000 -254430.745858
8.700e-08 300.000 101325.000 -254430.745858
8.710e-08 300.000 101325.000 -254430.745858
8.720e-08 300.000 101325.000 -254430.745858
8.730e-08 300.000 101325.000 -254430.745858
8.740e-08 300.000 101325.000 -254430.745858
8.750e-08 300.000 101325.000 -254430.745858
8.760e-08 300.000 101325.000 -254430.745858
8.770e-08 300.000 101325.000 -254430.745858
8.780e-08 300.000 101325.000 -254430.745858
8.790e-08 300.000 101325.000 -254430.745858
8.800e-08 300.000 101325.000 -254430.745858
8.810e-08 300.000 101325.000 -254430.745858
8.820e-08 300.000 101325.000 -254430.745858
8.830e-08 300.000 101325.000 -254430.745858
8.840e-08 300.000 101325.000 -254430.745858
8.850e-08 300.000 101325.000 -254430.745858
8.860e-08 300.000 101325.000 -254430.745858
8.870e-08 300.000 101325.000 -254430.745858
8.880e-08 300.000 101325.000 -254430.745858
8.890e-08 300.000 101325.000 -254430.745858
8.900e-08 300.000 101325.000 -254430.745858
8.910e-08 300.000 101325.000 -254430.745858
8.920e-08 300.000 101325.000 -254430.745858
8.930e-08 300.000 101325.000 -254430.745858
8.940e-08 300.000 101325.000 -254430.745858
8.950e-08 300.000 101325.000 -254430.745858
8.960e-08 300.000 101325.000 -254430.745858
8.970e-08 300.000 101325.000 -254430.745858
8.980e-08 300.000 101325.000 -254430.745858
8.990e-08 300.000 101325.000 -254430.745858
9.000e-08 300.000 101325.000 -254430.745858
9.010e-08 300.000 101325.000 -254430.745858
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="on", clone=False)
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.phase.state, t=sim.time)
print(f"{sim.time:10.3e} {r.T:10.3f} {r.phase.P:10.3f} {r.phase.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 0.856 seconds)