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.994548
1.660e-08 300.000 101325.000 -254430.994543
1.670e-08 300.000 101325.000 -254430.994538
1.680e-08 300.000 101325.000 -254430.994533
1.690e-08 300.000 101325.000 -254430.994528
1.700e-08 300.000 101325.000 -254430.994522
1.710e-08 300.000 101325.000 -254430.994517
1.720e-08 300.000 101325.000 -254430.994482
1.730e-08 300.000 101325.000 -254430.994448
1.740e-08 300.000 101325.000 -254430.994413
1.750e-08 300.000 101325.000 -254430.994379
1.760e-08 300.000 101325.000 -254430.994344
1.770e-08 300.000 101325.000 -254430.994310
1.780e-08 300.000 101325.000 -254430.994275
1.790e-08 300.000 101325.000 -254430.994240
1.800e-08 300.000 101325.000 -254430.994206
1.810e-08 300.000 101325.000 -254430.994171
1.820e-08 300.000 101325.000 -254430.993995
1.830e-08 300.000 101325.000 -254430.993819
1.840e-08 300.000 101325.000 -254430.993643
1.850e-08 300.000 101325.000 -254430.993467
1.860e-08 300.000 101325.000 -254430.993291
1.870e-08 300.000 101325.000 -254430.993115
1.880e-08 300.000 101325.000 -254430.992939
1.890e-08 300.000 101325.000 -254430.992764
1.900e-08 300.000 101325.000 -254430.992588
1.910e-08 300.000 101325.000 -254430.992412
1.920e-08 300.000 101325.000 -254430.991857
1.930e-08 300.000 101325.000 -254430.991303
1.940e-08 300.000 101325.000 -254430.990749
1.950e-08 300.000 101325.000 -254430.990195
1.960e-08 300.000 101325.000 -254430.989641
1.970e-08 300.000 101325.000 -254430.989087
1.980e-08 300.000 101325.000 -254430.988533
1.990e-08 300.000 101325.000 -254430.987979
2.000e-08 300.000 101325.000 -254430.987425
2.010e-08 300.000 101325.000 -254430.986871
2.020e-08 300.000 101325.000 -254430.985993
2.030e-08 300.000 101325.000 -254430.985114
2.040e-08 300.000 101325.000 -254430.984236
2.050e-08 300.000 101325.000 -254430.983358
2.060e-08 300.000 101325.000 -254430.982480
2.070e-08 300.000 101325.000 -254430.981602
2.080e-08 300.000 101325.000 -254430.980724
2.090e-08 300.000 101325.000 -254430.979846
2.100e-08 300.000 101325.000 -254430.978968
2.110e-08 300.000 101325.000 -254430.977733
2.120e-08 300.000 101325.000 -254430.976499
2.130e-08 300.000 101325.000 -254430.975265
2.140e-08 300.000 101325.000 -254430.974030
2.150e-08 300.000 101325.000 -254430.972796
2.160e-08 300.000 101325.000 -254430.971562
2.170e-08 300.000 101325.000 -254430.970327
2.180e-08 300.000 101325.000 -254430.969093
2.190e-08 300.000 101325.000 -254430.967859
2.200e-08 300.000 101325.000 -254430.966624
2.210e-08 300.000 101325.000 -254430.964756
2.220e-08 300.000 101325.000 -254430.962887
2.230e-08 300.000 101325.000 -254430.961019
2.240e-08 300.000 101325.000 -254430.959150
2.250e-08 300.000 101325.000 -254430.957282
2.260e-08 300.000 101325.000 -254430.955414
2.270e-08 300.000 101325.000 -254430.953545
2.280e-08 300.000 101325.000 -254430.951677
2.290e-08 300.000 101325.000 -254430.949808
2.300e-08 300.000 101325.000 -254430.947940
2.310e-08 300.000 101325.000 -254430.945239
2.320e-08 300.000 101325.000 -254430.942537
2.330e-08 300.000 101325.000 -254430.939836
2.340e-08 300.000 101325.000 -254430.937135
2.350e-08 300.000 101325.000 -254430.934434
2.360e-08 300.000 101325.000 -254430.931732
2.370e-08 300.000 101325.000 -254430.929031
2.380e-08 300.000 101325.000 -254430.926330
2.390e-08 300.000 101325.000 -254430.923629
2.400e-08 300.000 101325.000 -254430.920927
2.410e-08 300.000 101325.000 -254430.917483
2.420e-08 300.000 101325.000 -254430.914038
2.430e-08 300.000 101325.000 -254430.910594
2.440e-08 300.000 101325.000 -254430.907149
2.450e-08 300.000 101325.000 -254430.903704
2.460e-08 300.000 101325.000 -254430.900260
2.470e-08 300.000 101325.000 -254430.896815
2.480e-08 300.000 101325.000 -254430.893370
2.490e-08 300.000 101325.000 -254430.889926
2.500e-08 300.000 101325.000 -254430.886481
2.510e-08 300.000 101325.000 -254430.882735
2.520e-08 300.000 101325.000 -254430.878988
2.530e-08 300.000 101325.000 -254430.875242
2.540e-08 300.000 101325.000 -254430.871495
2.550e-08 300.000 101325.000 -254430.867749
2.560e-08 300.000 101325.000 -254430.864002
2.570e-08 300.000 101325.000 -254430.860256
2.580e-08 300.000 101325.000 -254430.856510
2.590e-08 300.000 101325.000 -254430.852763
2.600e-08 300.000 101325.000 -254430.849017
2.610e-08 300.000 101325.000 -254430.845572
2.620e-08 300.000 101325.000 -254430.842127
2.630e-08 300.000 101325.000 -254430.838683
2.640e-08 300.000 101325.000 -254430.835238
2.650e-08 300.000 101325.000 -254430.831794
2.660e-08 300.000 101325.000 -254430.828349
2.670e-08 300.000 101325.000 -254430.824904
2.680e-08 300.000 101325.000 -254430.821460
2.690e-08 300.000 101325.000 -254430.818015
2.700e-08 300.000 101325.000 -254430.814570
2.710e-08 300.000 101325.000 -254430.811869
2.720e-08 300.000 101325.000 -254430.809168
2.730e-08 300.000 101325.000 -254430.806467
2.740e-08 300.000 101325.000 -254430.803765
2.750e-08 300.000 101325.000 -254430.801064
2.760e-08 300.000 101325.000 -254430.798363
2.770e-08 300.000 101325.000 -254430.795662
2.780e-08 300.000 101325.000 -254430.792960
2.790e-08 300.000 101325.000 -254430.790259
2.800e-08 300.000 101325.000 -254430.787558
2.810e-08 300.000 101325.000 -254430.785690
2.820e-08 300.000 101325.000 -254430.783821
2.830e-08 300.000 101325.000 -254430.781953
2.840e-08 300.000 101325.000 -254430.780084
2.850e-08 300.000 101325.000 -254430.778216
2.860e-08 300.000 101325.000 -254430.776347
2.870e-08 300.000 101325.000 -254430.774479
2.880e-08 300.000 101325.000 -254430.772611
2.890e-08 300.000 101325.000 -254430.770742
2.900e-08 300.000 101325.000 -254430.768874
2.910e-08 300.000 101325.000 -254430.767639
2.920e-08 300.000 101325.000 -254430.766405
2.930e-08 300.000 101325.000 -254430.765171
2.940e-08 300.000 101325.000 -254430.763936
2.950e-08 300.000 101325.000 -254430.762702
2.960e-08 300.000 101325.000 -254430.761468
2.970e-08 300.000 101325.000 -254430.760233
2.980e-08 300.000 101325.000 -254430.758999
2.990e-08 300.000 101325.000 -254430.757765
3.000e-08 300.000 101325.000 -254430.756530
3.010e-08 300.000 101325.000 -254430.755652
3.020e-08 300.000 101325.000 -254430.754774
3.030e-08 300.000 101325.000 -254430.753896
3.040e-08 300.000 101325.000 -254430.753018
3.050e-08 300.000 101325.000 -254430.752140
3.060e-08 300.000 101325.000 -254430.751262
3.070e-08 300.000 101325.000 -254430.750384
3.080e-08 300.000 101325.000 -254430.749506
3.090e-08 300.000 101325.000 -254430.748628
3.100e-08 300.000 101325.000 -254430.747750
3.110e-08 300.000 101325.000 -254430.747196
3.120e-08 300.000 101325.000 -254430.746642
3.130e-08 300.000 101325.000 -254430.746088
3.140e-08 300.000 101325.000 -254430.745534
3.150e-08 300.000 101325.000 -254430.744981
3.160e-08 300.000 101325.000 -254430.744427
3.170e-08 300.000 101325.000 -254430.743873
3.180e-08 300.000 101325.000 -254430.743319
3.190e-08 300.000 101325.000 -254430.742765
3.200e-08 300.000 101325.000 -254430.742211
3.210e-08 300.000 101325.000 -254430.742035
3.220e-08 300.000 101325.000 -254430.741859
3.230e-08 300.000 101325.000 -254430.741684
3.240e-08 300.000 101325.000 -254430.741508
3.250e-08 300.000 101325.000 -254430.741332
3.260e-08 300.000 101325.000 -254430.741156
3.270e-08 300.000 101325.000 -254430.740980
3.280e-08 300.000 101325.000 -254430.740804
3.290e-08 300.000 101325.000 -254430.740628
3.300e-08 300.000 101325.000 -254430.740453
3.310e-08 300.000 101325.000 -254430.740418
3.320e-08 300.000 101325.000 -254430.740384
3.330e-08 300.000 101325.000 -254430.740349
3.340e-08 300.000 101325.000 -254430.740315
3.350e-08 300.000 101325.000 -254430.740280
3.360e-08 300.000 101325.000 -254430.740246
3.370e-08 300.000 101325.000 -254430.740211
3.380e-08 300.000 101325.000 -254430.740177
3.390e-08 300.000 101325.000 -254430.740142
3.400e-08 300.000 101325.000 -254430.740108
3.410e-08 300.000 101325.000 -254430.740103
3.420e-08 300.000 101325.000 -254430.740097
3.430e-08 300.000 101325.000 -254430.740092
3.440e-08 300.000 101325.000 -254430.740087
3.450e-08 300.000 101325.000 -254430.740082
3.460e-08 300.000 101325.000 -254430.740076
3.470e-08 300.000 101325.000 -254430.740071
3.480e-08 300.000 101325.000 -254430.740066
3.490e-08 300.000 101325.000 -254430.740061
3.500e-08 300.000 101325.000 -254430.740055
3.510e-08 300.000 101325.000 -254430.740055
3.520e-08 300.000 101325.000 -254430.740054
3.530e-08 300.000 101325.000 -254430.740054
3.540e-08 300.000 101325.000 -254430.740053
3.550e-08 300.000 101325.000 -254430.740052
3.560e-08 300.000 101325.000 -254430.740052
3.570e-08 300.000 101325.000 -254430.740051
3.580e-08 300.000 101325.000 -254430.740050
3.590e-08 300.000 101325.000 -254430.740050
3.600e-08 300.000 101325.000 -254430.740049
3.610e-08 300.000 101325.000 -254430.740049
3.620e-08 300.000 101325.000 -254430.740049
3.630e-08 300.000 101325.000 -254430.740049
3.640e-08 300.000 101325.000 -254430.740049
3.650e-08 300.000 101325.000 -254430.740049
3.660e-08 300.000 101325.000 -254430.740049
3.670e-08 300.000 101325.000 -254430.740049
3.680e-08 300.000 101325.000 -254430.740049
3.690e-08 300.000 101325.000 -254430.740049
3.700e-08 300.000 101325.000 -254430.740049
3.710e-08 300.000 101325.000 -254430.740049
3.720e-08 300.000 101325.000 -254430.740049
3.730e-08 300.000 101325.000 -254430.740049
3.740e-08 300.000 101325.000 -254430.740048
3.750e-08 300.000 101325.000 -254430.740048
3.760e-08 300.000 101325.000 -254430.740048
3.770e-08 300.000 101325.000 -254430.740048
3.780e-08 300.000 101325.000 -254430.740048
3.790e-08 300.000 101325.000 -254430.740048
3.800e-08 300.000 101325.000 -254430.740048
3.810e-08 300.000 101325.000 -254430.740048
3.820e-08 300.000 101325.000 -254430.740048
3.830e-08 300.000 101325.000 -254430.740048
3.840e-08 300.000 101325.000 -254430.740048
3.850e-08 300.000 101325.000 -254430.740048
3.860e-08 300.000 101325.000 -254430.740048
3.870e-08 300.000 101325.000 -254430.740048
3.880e-08 300.000 101325.000 -254430.740048
3.890e-08 300.000 101325.000 -254430.740048
3.900e-08 300.000 101325.000 -254430.740048
3.910e-08 300.000 101325.000 -254430.740048
3.920e-08 300.000 101325.000 -254430.740048
3.930e-08 300.000 101325.000 -254430.740048
3.940e-08 300.000 101325.000 -254430.740048
3.950e-08 300.000 101325.000 -254430.740048
3.960e-08 300.000 101325.000 -254430.740048
3.970e-08 300.000 101325.000 -254430.740048
3.980e-08 300.000 101325.000 -254430.740048
3.990e-08 300.000 101325.000 -254430.740048
4.000e-08 300.000 101325.000 -254430.740048
4.010e-08 300.000 101325.000 -254430.740048
4.020e-08 300.000 101325.000 -254430.740048
4.030e-08 300.000 101325.000 -254430.740048
4.040e-08 300.000 101325.000 -254430.740048
4.050e-08 300.000 101325.000 -254430.740048
4.060e-08 300.000 101325.000 -254430.740048
4.070e-08 300.000 101325.000 -254430.740048
4.080e-08 300.000 101325.000 -254430.740048
4.090e-08 300.000 101325.000 -254430.740048
4.100e-08 300.000 101325.000 -254430.740048
4.110e-08 300.000 101325.000 -254430.740048
4.120e-08 300.000 101325.000 -254430.740048
4.130e-08 300.000 101325.000 -254430.740048
4.140e-08 300.000 101325.000 -254430.740048
4.150e-08 300.000 101325.000 -254430.740048
4.160e-08 300.000 101325.000 -254430.740048
4.170e-08 300.000 101325.000 -254430.740048
4.180e-08 300.000 101325.000 -254430.740048
4.190e-08 300.000 101325.000 -254430.740048
4.200e-08 300.000 101325.000 -254430.740048
4.210e-08 300.000 101325.000 -254430.740048
4.220e-08 300.000 101325.000 -254430.740048
4.230e-08 300.000 101325.000 -254430.740048
4.240e-08 300.000 101325.000 -254430.740048
4.250e-08 300.000 101325.000 -254430.740048
4.260e-08 300.000 101325.000 -254430.740048
4.270e-08 300.000 101325.000 -254430.740048
4.280e-08 300.000 101325.000 -254430.740048
4.290e-08 300.000 101325.000 -254430.740048
4.300e-08 300.000 101325.000 -254430.740048
4.310e-08 300.000 101325.000 -254430.740048
4.320e-08 300.000 101325.000 -254430.740048
4.330e-08 300.000 101325.000 -254430.740048
4.340e-08 300.000 101325.000 -254430.740048
4.350e-08 300.000 101325.000 -254430.740048
4.360e-08 300.000 101325.000 -254430.740048
4.370e-08 300.000 101325.000 -254430.740048
4.380e-08 300.000 101325.000 -254430.740048
4.390e-08 300.000 101325.000 -254430.740048
4.400e-08 300.000 101325.000 -254430.740048
4.410e-08 300.000 101325.000 -254430.740048
4.420e-08 300.000 101325.000 -254430.740048
4.430e-08 300.000 101325.000 -254430.740048
4.440e-08 300.000 101325.000 -254430.740048
4.450e-08 300.000 101325.000 -254430.740048
4.460e-08 300.000 101325.000 -254430.740048
4.470e-08 300.000 101325.000 -254430.740048
4.480e-08 300.000 101325.000 -254430.740048
4.490e-08 300.000 101325.000 -254430.740048
4.500e-08 300.000 101325.000 -254430.740048
4.510e-08 300.000 101325.000 -254430.740048
4.520e-08 300.000 101325.000 -254430.740048
4.530e-08 300.000 101325.000 -254430.740048
4.540e-08 300.000 101325.000 -254430.740048
4.550e-08 300.000 101325.000 -254430.740048
4.560e-08 300.000 101325.000 -254430.740048
4.570e-08 300.000 101325.000 -254430.740048
4.580e-08 300.000 101325.000 -254430.740048
4.590e-08 300.000 101325.000 -254430.740048
4.600e-08 300.000 101325.000 -254430.740048
4.610e-08 300.000 101325.000 -254430.740048
4.620e-08 300.000 101325.000 -254430.740048
4.630e-08 300.000 101325.000 -254430.740048
4.640e-08 300.000 101325.000 -254430.740048
4.650e-08 300.000 101325.000 -254430.740048
4.660e-08 300.000 101325.000 -254430.740048
4.670e-08 300.000 101325.000 -254430.740048
4.680e-08 300.000 101325.000 -254430.740048
4.690e-08 300.000 101325.000 -254430.740048
4.700e-08 300.000 101325.000 -254430.740048
4.710e-08 300.000 101325.000 -254430.740048
4.720e-08 300.000 101325.000 -254430.740048
4.730e-08 300.000 101325.000 -254430.740048
4.740e-08 300.000 101325.000 -254430.740048
4.750e-08 300.000 101325.000 -254430.740048
4.760e-08 300.000 101325.000 -254430.740048
4.770e-08 300.000 101325.000 -254430.740048
4.780e-08 300.000 101325.000 -254430.740048
4.790e-08 300.000 101325.000 -254430.740048
4.800e-08 300.000 101325.000 -254430.740048
4.810e-08 300.000 101325.000 -254430.740048
4.820e-08 300.000 101325.000 -254430.740048
4.830e-08 300.000 101325.000 -254430.740048
4.840e-08 300.000 101325.000 -254430.740048
4.850e-08 300.000 101325.000 -254430.740048
4.860e-08 300.000 101325.000 -254430.740048
4.870e-08 300.000 101325.000 -254430.740048
4.880e-08 300.000 101325.000 -254430.740048
4.890e-08 300.000 101325.000 -254430.740048
4.900e-08 300.000 101325.000 -254430.740048
4.910e-08 300.000 101325.000 -254430.740048
4.920e-08 300.000 101325.000 -254430.740048
4.930e-08 300.000 101325.000 -254430.740048
4.940e-08 300.000 101325.000 -254430.740048
4.950e-08 300.000 101325.000 -254430.740048
4.960e-08 300.000 101325.000 -254430.740048
4.970e-08 300.000 101325.000 -254430.740048
4.980e-08 300.000 101325.000 -254430.740048
4.990e-08 300.000 101325.000 -254430.740048
5.000e-08 300.000 101325.000 -254430.740048
5.010e-08 300.000 101325.000 -254430.740048
5.020e-08 300.000 101325.000 -254430.740048
5.030e-08 300.000 101325.000 -254430.740048
5.040e-08 300.000 101325.000 -254430.740048
5.050e-08 300.000 101325.000 -254430.740048
5.060e-08 300.000 101325.000 -254430.740048
5.070e-08 300.000 101325.000 -254430.740048
5.080e-08 300.000 101325.000 -254430.740048
5.090e-08 300.000 101325.000 -254430.740048
5.100e-08 300.000 101325.000 -254430.740048
5.110e-08 300.000 101325.000 -254430.740048
5.120e-08 300.000 101325.000 -254430.740048
5.130e-08 300.000 101325.000 -254430.740048
5.140e-08 300.000 101325.000 -254430.740048
5.150e-08 300.000 101325.000 -254430.740048
5.160e-08 300.000 101325.000 -254430.740048
5.170e-08 300.000 101325.000 -254430.740048
5.180e-08 300.000 101325.000 -254430.740048
5.190e-08 300.000 101325.000 -254430.740048
5.200e-08 300.000 101325.000 -254430.740048
5.210e-08 300.000 101325.000 -254430.740048
5.220e-08 300.000 101325.000 -254430.740048
5.230e-08 300.000 101325.000 -254430.740048
5.240e-08 300.000 101325.000 -254430.740048
5.250e-08 300.000 101325.000 -254430.740048
5.260e-08 300.000 101325.000 -254430.740048
5.270e-08 300.000 101325.000 -254430.740048
5.280e-08 300.000 101325.000 -254430.740048
5.290e-08 300.000 101325.000 -254430.740048
5.300e-08 300.000 101325.000 -254430.740048
5.310e-08 300.000 101325.000 -254430.740048
5.320e-08 300.000 101325.000 -254430.740048
5.330e-08 300.000 101325.000 -254430.740048
5.340e-08 300.000 101325.000 -254430.740048
5.350e-08 300.000 101325.000 -254430.740048
5.360e-08 300.000 101325.000 -254430.740048
5.370e-08 300.000 101325.000 -254430.740048
5.380e-08 300.000 101325.000 -254430.740048
5.390e-08 300.000 101325.000 -254430.740048
5.400e-08 300.000 101325.000 -254430.740048
5.410e-08 300.000 101325.000 -254430.740048
5.420e-08 300.000 101325.000 -254430.740048
5.430e-08 300.000 101325.000 -254430.740048
5.440e-08 300.000 101325.000 -254430.740048
5.450e-08 300.000 101325.000 -254430.740048
5.460e-08 300.000 101325.000 -254430.740048
5.470e-08 300.000 101325.000 -254430.740048
5.480e-08 300.000 101325.000 -254430.740048
5.490e-08 300.000 101325.000 -254430.740048
5.500e-08 300.000 101325.000 -254430.740048
5.510e-08 300.000 101325.000 -254430.740048
5.520e-08 300.000 101325.000 -254430.740048
5.530e-08 300.000 101325.000 -254430.740048
5.540e-08 300.000 101325.000 -254430.740048
5.550e-08 300.000 101325.000 -254430.740048
5.560e-08 300.000 101325.000 -254430.740048
5.570e-08 300.000 101325.000 -254430.740048
5.580e-08 300.000 101325.000 -254430.740048
5.590e-08 300.000 101325.000 -254430.740048
5.600e-08 300.000 101325.000 -254430.740048
5.610e-08 300.000 101325.000 -254430.740048
5.620e-08 300.000 101325.000 -254430.740048
5.630e-08 300.000 101325.000 -254430.740048
5.640e-08 300.000 101325.000 -254430.740048
5.650e-08 300.000 101325.000 -254430.740048
5.660e-08 300.000 101325.000 -254430.740048
5.670e-08 300.000 101325.000 -254430.740048
5.680e-08 300.000 101325.000 -254430.740048
5.690e-08 300.000 101325.000 -254430.740048
5.700e-08 300.000 101325.000 -254430.740048
5.710e-08 300.000 101325.000 -254430.740048
5.720e-08 300.000 101325.000 -254430.740048
5.730e-08 300.000 101325.000 -254430.740048
5.740e-08 300.000 101325.000 -254430.740048
5.750e-08 300.000 101325.000 -254430.740048
5.760e-08 300.000 101325.000 -254430.740048
5.770e-08 300.000 101325.000 -254430.740048
5.780e-08 300.000 101325.000 -254430.740048
5.790e-08 300.000 101325.000 -254430.740048
5.800e-08 300.000 101325.000 -254430.740048
5.810e-08 300.000 101325.000 -254430.740048
5.820e-08 300.000 101325.000 -254430.740048
5.830e-08 300.000 101325.000 -254430.740048
5.840e-08 300.000 101325.000 -254430.740048
5.850e-08 300.000 101325.000 -254430.740048
5.860e-08 300.000 101325.000 -254430.740048
5.870e-08 300.000 101325.000 -254430.740048
5.880e-08 300.000 101325.000 -254430.740048
5.890e-08 300.000 101325.000 -254430.740048
5.900e-08 300.000 101325.000 -254430.740048
5.910e-08 300.000 101325.000 -254430.740048
5.920e-08 300.000 101325.000 -254430.740048
5.930e-08 300.000 101325.000 -254430.740048
5.940e-08 300.000 101325.000 -254430.740048
5.950e-08 300.000 101325.000 -254430.740048
5.960e-08 300.000 101325.000 -254430.740048
5.970e-08 300.000 101325.000 -254430.740048
5.980e-08 300.000 101325.000 -254430.740048
5.990e-08 300.000 101325.000 -254430.740048
6.000e-08 300.000 101325.000 -254430.740048
6.010e-08 300.000 101325.000 -254430.740048
6.020e-08 300.000 101325.000 -254430.740048
6.030e-08 300.000 101325.000 -254430.740048
6.040e-08 300.000 101325.000 -254430.740048
6.050e-08 300.000 101325.000 -254430.740048
6.060e-08 300.000 101325.000 -254430.740048
6.070e-08 300.000 101325.000 -254430.740048
6.080e-08 300.000 101325.000 -254430.740048
6.090e-08 300.000 101325.000 -254430.740048
6.100e-08 300.000 101325.000 -254430.740048
6.110e-08 300.000 101325.000 -254430.740048
6.120e-08 300.000 101325.000 -254430.740048
6.130e-08 300.000 101325.000 -254430.740048
6.140e-08 300.000 101325.000 -254430.740048
6.150e-08 300.000 101325.000 -254430.740048
6.160e-08 300.000 101325.000 -254430.740048
6.170e-08 300.000 101325.000 -254430.740048
6.180e-08 300.000 101325.000 -254430.740048
6.190e-08 300.000 101325.000 -254430.740048
6.200e-08 300.000 101325.000 -254430.740048
6.210e-08 300.000 101325.000 -254430.740048
6.220e-08 300.000 101325.000 -254430.740048
6.230e-08 300.000 101325.000 -254430.740048
6.240e-08 300.000 101325.000 -254430.740048
6.250e-08 300.000 101325.000 -254430.740048
6.260e-08 300.000 101325.000 -254430.740048
6.270e-08 300.000 101325.000 -254430.740048
6.280e-08 300.000 101325.000 -254430.740048
6.290e-08 300.000 101325.000 -254430.740048
6.300e-08 300.000 101325.000 -254430.740048
6.310e-08 300.000 101325.000 -254430.740048
6.320e-08 300.000 101325.000 -254430.740048
6.330e-08 300.000 101325.000 -254430.740048
6.340e-08 300.000 101325.000 -254430.740048
6.350e-08 300.000 101325.000 -254430.740048
6.360e-08 300.000 101325.000 -254430.740048
6.370e-08 300.000 101325.000 -254430.740048
6.380e-08 300.000 101325.000 -254430.740048
6.390e-08 300.000 101325.000 -254430.740048
6.400e-08 300.000 101325.000 -254430.740048
6.410e-08 300.000 101325.000 -254430.740048
6.420e-08 300.000 101325.000 -254430.740048
6.430e-08 300.000 101325.000 -254430.740048
6.440e-08 300.000 101325.000 -254430.740048
6.450e-08 300.000 101325.000 -254430.740048
6.460e-08 300.000 101325.000 -254430.740048
6.470e-08 300.000 101325.000 -254430.740048
6.480e-08 300.000 101325.000 -254430.740048
6.490e-08 300.000 101325.000 -254430.740048
6.500e-08 300.000 101325.000 -254430.740048
6.510e-08 300.000 101325.000 -254430.740048
6.520e-08 300.000 101325.000 -254430.740048
6.530e-08 300.000 101325.000 -254430.740048
6.540e-08 300.000 101325.000 -254430.740048
6.550e-08 300.000 101325.000 -254430.740048
6.560e-08 300.000 101325.000 -254430.740048
6.570e-08 300.000 101325.000 -254430.740048
6.580e-08 300.000 101325.000 -254430.740048
6.590e-08 300.000 101325.000 -254430.740048
6.600e-08 300.000 101325.000 -254430.740048
6.610e-08 300.000 101325.000 -254430.740048
6.620e-08 300.000 101325.000 -254430.740048
6.630e-08 300.000 101325.000 -254430.740048
6.640e-08 300.000 101325.000 -254430.740048
6.650e-08 300.000 101325.000 -254430.740048
6.660e-08 300.000 101325.000 -254430.740048
6.670e-08 300.000 101325.000 -254430.740048
6.680e-08 300.000 101325.000 -254430.740048
6.690e-08 300.000 101325.000 -254430.740048
6.700e-08 300.000 101325.000 -254430.740048
6.710e-08 300.000 101325.000 -254430.740048
6.720e-08 300.000 101325.000 -254430.740048
6.730e-08 300.000 101325.000 -254430.740048
6.740e-08 300.000 101325.000 -254430.740048
6.750e-08 300.000 101325.000 -254430.740048
6.760e-08 300.000 101325.000 -254430.740048
6.770e-08 300.000 101325.000 -254430.740048
6.780e-08 300.000 101325.000 -254430.740048
6.790e-08 300.000 101325.000 -254430.740048
6.800e-08 300.000 101325.000 -254430.740048
6.810e-08 300.000 101325.000 -254430.740048
6.820e-08 300.000 101325.000 -254430.740048
6.830e-08 300.000 101325.000 -254430.740048
6.840e-08 300.000 101325.000 -254430.740048
6.850e-08 300.000 101325.000 -254430.740048
6.860e-08 300.000 101325.000 -254430.740048
6.870e-08 300.000 101325.000 -254430.740048
6.880e-08 300.000 101325.000 -254430.740048
6.890e-08 300.000 101325.000 -254430.740048
6.900e-08 300.000 101325.000 -254430.740048
6.910e-08 300.000 101325.000 -254430.740048
6.920e-08 300.000 101325.000 -254430.740048
6.930e-08 300.000 101325.000 -254430.740048
6.940e-08 300.000 101325.000 -254430.740048
6.950e-08 300.000 101325.000 -254430.740048
6.960e-08 300.000 101325.000 -254430.740048
6.970e-08 300.000 101325.000 -254430.740048
6.980e-08 300.000 101325.000 -254430.740048
6.990e-08 300.000 101325.000 -254430.740048
7.000e-08 300.000 101325.000 -254430.740048
7.010e-08 300.000 101325.000 -254430.740048
7.020e-08 300.000 101325.000 -254430.740048
7.030e-08 300.000 101325.000 -254430.740048
7.040e-08 300.000 101325.000 -254430.740048
7.050e-08 300.000 101325.000 -254430.740048
7.060e-08 300.000 101325.000 -254430.740048
7.070e-08 300.000 101325.000 -254430.740048
7.080e-08 300.000 101325.000 -254430.740048
7.090e-08 300.000 101325.000 -254430.740048
7.100e-08 300.000 101325.000 -254430.740048
7.110e-08 300.000 101325.000 -254430.740048
7.120e-08 300.000 101325.000 -254430.740048
7.130e-08 300.000 101325.000 -254430.740048
7.140e-08 300.000 101325.000 -254430.740048
7.150e-08 300.000 101325.000 -254430.740048
7.160e-08 300.000 101325.000 -254430.740048
7.170e-08 300.000 101325.000 -254430.740048
7.180e-08 300.000 101325.000 -254430.740048
7.190e-08 300.000 101325.000 -254430.740048
7.200e-08 300.000 101325.000 -254430.740048
7.210e-08 300.000 101325.000 -254430.740048
7.220e-08 300.000 101325.000 -254430.740048
7.230e-08 300.000 101325.000 -254430.740048
7.240e-08 300.000 101325.000 -254430.740048
7.250e-08 300.000 101325.000 -254430.740048
7.260e-08 300.000 101325.000 -254430.740048
7.270e-08 300.000 101325.000 -254430.740048
7.280e-08 300.000 101325.000 -254430.740048
7.290e-08 300.000 101325.000 -254430.740048
7.300e-08 300.000 101325.000 -254430.740048
7.310e-08 300.000 101325.000 -254430.740048
7.320e-08 300.000 101325.000 -254430.740048
7.330e-08 300.000 101325.000 -254430.740048
7.340e-08 300.000 101325.000 -254430.740048
7.350e-08 300.000 101325.000 -254430.740048
7.360e-08 300.000 101325.000 -254430.740048
7.370e-08 300.000 101325.000 -254430.740048
7.380e-08 300.000 101325.000 -254430.740048
7.390e-08 300.000 101325.000 -254430.740048
7.400e-08 300.000 101325.000 -254430.740048
7.410e-08 300.000 101325.000 -254430.740048
7.420e-08 300.000 101325.000 -254430.740048
7.430e-08 300.000 101325.000 -254430.740048
7.440e-08 300.000 101325.000 -254430.740048
7.450e-08 300.000 101325.000 -254430.740048
7.460e-08 300.000 101325.000 -254430.740048
7.470e-08 300.000 101325.000 -254430.740048
7.480e-08 300.000 101325.000 -254430.740048
7.490e-08 300.000 101325.000 -254430.740048
7.500e-08 300.000 101325.000 -254430.740048
7.510e-08 300.000 101325.000 -254430.740048
7.520e-08 300.000 101325.000 -254430.740048
7.530e-08 300.000 101325.000 -254430.740048
7.540e-08 300.000 101325.000 -254430.740048
7.550e-08 300.000 101325.000 -254430.740048
7.560e-08 300.000 101325.000 -254430.740048
7.570e-08 300.000 101325.000 -254430.740048
7.580e-08 300.000 101325.000 -254430.740048
7.590e-08 300.000 101325.000 -254430.740048
7.600e-08 300.000 101325.000 -254430.740048
7.610e-08 300.000 101325.000 -254430.740048
7.620e-08 300.000 101325.000 -254430.740048
7.630e-08 300.000 101325.000 -254430.740048
7.640e-08 300.000 101325.000 -254430.740048
7.650e-08 300.000 101325.000 -254430.740048
7.660e-08 300.000 101325.000 -254430.740048
7.670e-08 300.000 101325.000 -254430.740048
7.680e-08 300.000 101325.000 -254430.740048
7.690e-08 300.000 101325.000 -254430.740048
7.700e-08 300.000 101325.000 -254430.740048
7.710e-08 300.000 101325.000 -254430.740048
7.720e-08 300.000 101325.000 -254430.740048
7.730e-08 300.000 101325.000 -254430.740048
7.740e-08 300.000 101325.000 -254430.740048
7.750e-08 300.000 101325.000 -254430.740048
7.760e-08 300.000 101325.000 -254430.740048
7.770e-08 300.000 101325.000 -254430.740048
7.780e-08 300.000 101325.000 -254430.740048
7.790e-08 300.000 101325.000 -254430.740048
7.800e-08 300.000 101325.000 -254430.740048
7.810e-08 300.000 101325.000 -254430.740048
7.820e-08 300.000 101325.000 -254430.740048
7.830e-08 300.000 101325.000 -254430.740048
7.840e-08 300.000 101325.000 -254430.740048
7.850e-08 300.000 101325.000 -254430.740048
7.860e-08 300.000 101325.000 -254430.740048
7.870e-08 300.000 101325.000 -254430.740048
7.880e-08 300.000 101325.000 -254430.740048
7.890e-08 300.000 101325.000 -254430.740048
7.900e-08 300.000 101325.000 -254430.740048
7.910e-08 300.000 101325.000 -254430.740048
7.920e-08 300.000 101325.000 -254430.740048
7.930e-08 300.000 101325.000 -254430.740048
7.940e-08 300.000 101325.000 -254430.740048
7.950e-08 300.000 101325.000 -254430.740048
7.960e-08 300.000 101325.000 -254430.740048
7.970e-08 300.000 101325.000 -254430.740048
7.980e-08 300.000 101325.000 -254430.740048
7.990e-08 300.000 101325.000 -254430.740048
8.000e-08 300.000 101325.000 -254430.740048
8.010e-08 300.000 101325.000 -254430.740048
8.020e-08 300.000 101325.000 -254430.740048
8.030e-08 300.000 101325.000 -254430.740048
8.040e-08 300.000 101325.000 -254430.740048
8.050e-08 300.000 101325.000 -254430.740048
8.060e-08 300.000 101325.000 -254430.740048
8.070e-08 300.000 101325.000 -254430.740048
8.080e-08 300.000 101325.000 -254430.740048
8.090e-08 300.000 101325.000 -254430.740048
8.100e-08 300.000 101325.000 -254430.740048
8.110e-08 300.000 101325.000 -254430.740048
8.120e-08 300.000 101325.000 -254430.740048
8.130e-08 300.000 101325.000 -254430.740048
8.140e-08 300.000 101325.000 -254430.740048
8.150e-08 300.000 101325.000 -254430.740048
8.160e-08 300.000 101325.000 -254430.740048
8.170e-08 300.000 101325.000 -254430.740048
8.180e-08 300.000 101325.000 -254430.740048
8.190e-08 300.000 101325.000 -254430.740048
8.200e-08 300.000 101325.000 -254430.740048
8.210e-08 300.000 101325.000 -254430.740048
8.220e-08 300.000 101325.000 -254430.740048
8.230e-08 300.000 101325.000 -254430.740048
8.240e-08 300.000 101325.000 -254430.740048
8.250e-08 300.000 101325.000 -254430.740048
8.260e-08 300.000 101325.000 -254430.740048
8.270e-08 300.000 101325.000 -254430.740048
8.280e-08 300.000 101325.000 -254430.740048
8.290e-08 300.000 101325.000 -254430.740048
8.300e-08 300.000 101325.000 -254430.740048
8.310e-08 300.000 101325.000 -254430.740048
8.320e-08 300.000 101325.000 -254430.740048
8.330e-08 300.000 101325.000 -254430.740048
8.340e-08 300.000 101325.000 -254430.740048
8.350e-08 300.000 101325.000 -254430.740048
8.360e-08 300.000 101325.000 -254430.740048
8.370e-08 300.000 101325.000 -254430.740048
8.380e-08 300.000 101325.000 -254430.740048
8.390e-08 300.000 101325.000 -254430.740048
8.400e-08 300.000 101325.000 -254430.740048
8.410e-08 300.000 101325.000 -254430.740048
8.420e-08 300.000 101325.000 -254430.740048
8.430e-08 300.000 101325.000 -254430.740048
8.440e-08 300.000 101325.000 -254430.740048
8.450e-08 300.000 101325.000 -254430.740048
8.460e-08 300.000 101325.000 -254430.740048
8.470e-08 300.000 101325.000 -254430.740048
8.480e-08 300.000 101325.000 -254430.740048
8.490e-08 300.000 101325.000 -254430.740048
8.500e-08 300.000 101325.000 -254430.740048
8.510e-08 300.000 101325.000 -254430.740048
8.520e-08 300.000 101325.000 -254430.740048
8.530e-08 300.000 101325.000 -254430.740048
8.540e-08 300.000 101325.000 -254430.740048
8.550e-08 300.000 101325.000 -254430.740048
8.560e-08 300.000 101325.000 -254430.740048
8.570e-08 300.000 101325.000 -254430.740048
8.580e-08 300.000 101325.000 -254430.740048
8.590e-08 300.000 101325.000 -254430.740048
8.600e-08 300.000 101325.000 -254430.740048
8.610e-08 300.000 101325.000 -254430.740048
8.620e-08 300.000 101325.000 -254430.740048
8.630e-08 300.000 101325.000 -254430.740048
8.640e-08 300.000 101325.000 -254430.740048
8.650e-08 300.000 101325.000 -254430.740048
8.660e-08 300.000 101325.000 -254430.740048
8.670e-08 300.000 101325.000 -254430.740048
8.680e-08 300.000 101325.000 -254430.740048
8.690e-08 300.000 101325.000 -254430.740048
8.700e-08 300.000 101325.000 -254430.740048
8.710e-08 300.000 101325.000 -254430.740048
8.720e-08 300.000 101325.000 -254430.740048
8.730e-08 300.000 101325.000 -254430.740048
8.740e-08 300.000 101325.000 -254430.740048
8.750e-08 300.000 101325.000 -254430.740048
8.760e-08 300.000 101325.000 -254430.740048
8.770e-08 300.000 101325.000 -254430.740048
8.780e-08 300.000 101325.000 -254430.740048
8.790e-08 300.000 101325.000 -254430.740048
8.800e-08 300.000 101325.000 -254430.740048
8.810e-08 300.000 101325.000 -254430.740048
8.820e-08 300.000 101325.000 -254430.740048
8.830e-08 300.000 101325.000 -254430.740048
8.840e-08 300.000 101325.000 -254430.740048
8.850e-08 300.000 101325.000 -254430.740048
8.860e-08 300.000 101325.000 -254430.740048
8.870e-08 300.000 101325.000 -254430.740048
8.880e-08 300.000 101325.000 -254430.740048
8.890e-08 300.000 101325.000 -254430.740048
8.900e-08 300.000 101325.000 -254430.740048
8.910e-08 300.000 101325.000 -254430.740048
8.920e-08 300.000 101325.000 -254430.740048
8.930e-08 300.000 101325.000 -254430.740048
8.940e-08 300.000 101325.000 -254430.740048
8.950e-08 300.000 101325.000 -254430.740048
8.960e-08 300.000 101325.000 -254430.740048
8.970e-08 300.000 101325.000 -254430.740048
8.980e-08 300.000 101325.000 -254430.740048
8.990e-08 300.000 101325.000 -254430.740048
9.000e-08 300.000 101325.000 -254430.740048
9.010e-08 300.000 101325.000 -254430.740048
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.431 seconds)