
:html_theme.sidebar_secondary.remove:

.. py:currentmodule:: cantera


.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "examples/python/onedim/catalytic_combustion.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_examples_python_onedim_catalytic_combustion.py>`
        to download the full example code.

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_examples_python_onedim_catalytic_combustion.py:


Catalytic combustion of methane on platinum
===========================================

This script solves a catalytic combustion problem. A stagnation flow is set
up, with a gas inlet 10 cm from a platinum surface at 900 K. The lean,
premixed methane/air mixture enters at ~6 cm/s (0.06 kg/m2/s), and burns
catalytically on the platinum surface. Gas-phase chemistry is included too,
and has some effect very near the surface.

The catalytic combustion mechanism is from Deutschmann et al., 26th
Symp. (Intl.) on Combustion,1996 pp. 1747-1754

Requires: cantera >= 3.0, matplotlib >= 2.0

.. tags:: Python, catalysis, combustion, 1D flow, surface chemistry

.. GENERATED FROM PYTHON SOURCE LINES 18-23

.. code-block:: Python


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








.. GENERATED FROM PYTHON SOURCE LINES 24-28

Problem Definition
------------------

Parameter values are collected here to make it easier to modify them

.. GENERATED FROM PYTHON SOURCE LINES 28-34

.. code-block:: Python

    p = ct.one_atm  # pressure
    tinlet = 300.0  # inlet temperature
    tsurf = 900.0  # surface temperature
    mdot = 0.06  # kg/m^2/s
    transport = 'mixture-averaged'  # transport model








.. GENERATED FROM PYTHON SOURCE LINES 35-37

We will solve first for a hydrogen/air case to use as the initial estimate
for the methane/air case

.. GENERATED FROM PYTHON SOURCE LINES 37-49

.. code-block:: Python


    # composition of the inlet premixed gas for the hydrogen/air case
    comp1 = 'H2:0.05, O2:0.21, N2:0.78, AR:0.01'

    # composition of the inlet premixed gas for the methane/air case
    comp2 = 'CH4:0.095, O2:0.21, N2:0.78, AR:0.01'

    # The inlet/surface separation is 10 cm.
    width = 0.1  # m

    loglevel = 1  # amount of diagnostic output (0 to 5)








.. GENERATED FROM PYTHON SOURCE LINES 50-61

Create the phase objects
------------------------

The ``surf_phase`` object will be used to evaluate all surface chemical production
rates. It will be created from the interface definition ``Pt_surf`` in input file
``ptcombust.yaml``, which implements the reaction mechanism of Deutschmann et
al., 1995 for catalytic combustion on platinum.

This phase definition also references the phase ``gas`` in the same input file,
which will be created and used to evaluate all thermodynamic, kinetic, and
transport properties. It is a stripped-down version of GRI-Mech 3.0.

.. GENERATED FROM PYTHON SOURCE LINES 61-82

.. code-block:: Python

    surf_phase = ct.Interface("ptcombust.yaml", "Pt_surf")
    surf_phase.TP = tsurf, p
    gas = surf_phase.adjacent["gas"]
    gas.TPX = tinlet, p, comp1

    # integrate the coverage equations in time for 1 s, holding the gas
    # composition fixed to generate a good starting estimate for the coverages.
    surf_phase.coverages = {'PT(S)': 0.5, 'O(S)':0.5}
    surf_phase.advance_coverages(1.0)

    # create the object that simulates the stagnation flow, and specify an initial
    # grid
    sim = ct.ImpingingJet(gas=gas, width=width, surface=surf_phase)

    # Objects of class ImpingingJet have members that represent the gas inlet
    # ('inlet') and the surface ('surface'). Set some parameters of these objects.
    sim.inlet.mdot = mdot
    sim.inlet.T = tinlet
    sim.inlet.X = comp1
    sim.surface.T = tsurf








.. GENERATED FROM PYTHON SOURCE LINES 83-84

Show the initial solution estimate

.. GENERATED FROM PYTHON SOURCE LINES 84-86

.. code-block:: Python

    sim.show()





.. rst-class:: sphx-glr-script-out

 .. code-block:: none



    >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> inlet <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

        Mass Flux:         0.06 kg/m^2/s 
        Temperature:        300 K 
        Mass Fractions: 
                          H2    0.003467 
                          O2      0.2311 
                          AR     0.01374 
                          N2      0.7516 



    >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> flame <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

        Pressure:   1.013e+05 Pa

    -------------------------------------------------------------------------------
              z    velocity  spreadRate           T      Lambda      eField 
    -------------------------------------------------------------------------------
              0     0.05335           0         300           0           0 
           0.02     0.04268           0         420           0           0 
           0.04     0.03201           0         540           0           0 
           0.06     0.02134           0         660           0           0 
           0.08     0.01067           0         780           0           0 
            0.1           0           0         900           0           0 

    -------------------------------------------------------------------------------
              z          Uo          H2           H           O          O2 
    -------------------------------------------------------------------------------
              0           0    0.003467           0           0      0.2311 
           0.02           0    0.003467           0           0      0.2311 
           0.04           0    0.003467           0           0      0.2311 
           0.06           0    0.003467           0           0      0.2311 
           0.08           0    0.003467           0           0      0.2311 
            0.1           0    0.003467           0           0      0.2311 

    -------------------------------------------------------------------------------
              z          OH         H2O         HO2        H2O2           C 
    -------------------------------------------------------------------------------
              0           0           0           0           0           0 
           0.02           0           0           0           0           0 
           0.04           0           0           0           0           0 
           0.06           0           0           0           0           0 
           0.08           0           0           0           0           0 
            0.1           0           0           0           0           0 

    -------------------------------------------------------------------------------
              z          CH         CH2      CH2(S)         CH3         CH4 
    -------------------------------------------------------------------------------
              0           0           0           0           0           0 
           0.02           0           0           0           0           0 
           0.04           0           0           0           0           0 
           0.06           0           0           0           0           0 
           0.08           0           0           0           0           0 
            0.1           0           0           0           0           0 

    -------------------------------------------------------------------------------
              z          CO         CO2         HCO        CH2O       CH2OH 
    -------------------------------------------------------------------------------
              0           0           0           0           0           0 
           0.02           0           0           0           0           0 
           0.04           0           0           0           0           0 
           0.06           0           0           0           0           0 
           0.08           0           0           0           0           0 
            0.1           0           0           0           0           0 

    -------------------------------------------------------------------------------
              z        CH3O       CH3OH         C2H        C2H2        C2H3 
    -------------------------------------------------------------------------------
              0           0           0           0           0           0 
           0.02           0           0           0           0           0 
           0.04           0           0           0           0           0 
           0.06           0           0           0           0           0 
           0.08           0           0           0           0           0 
            0.1           0           0           0           0           0 

    -------------------------------------------------------------------------------
              z        C2H4        C2H5        C2H6        HCCO       CH2CO 
    -------------------------------------------------------------------------------
              0           0           0           0           0           0 
           0.02           0           0           0           0           0 
           0.04           0           0           0           0           0 
           0.06           0           0           0           0           0 
           0.08           0           0           0           0           0 
            0.1           0           0           0           0           0 

    -------------------------------------------------------------------------------
              z       HCCOH          AR          N2 
    -------------------------------------------------------------------------------
              0           0     0.01374      0.7516 
           0.02           0     0.01374      0.7516 
           0.04           0     0.01374      0.7516 
           0.06           0     0.01374      0.7516 
           0.08           0     0.01374      0.7516 
            0.1           0     0.01374      0.7516 


    >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> surface <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

        Temperature:        900 K 
        Coverages: 
                       PT(S)     0.4879 
                        H(S)  4.738e-06 
                      H2O(S)  2.755e-05 
                       OH(S)    0.01028 
                       CO(S)  4.549e-22 
                      CO2(S)   2.85e-27 
                      CH3(S)          0 
                     CH2(S)s          0 
                       CH(S)          0 
                        C(S)          0 
                        O(S)     0.5018 





.. GENERATED FROM PYTHON SOURCE LINES 87-90

Solving problems with stiff chemistry coupled to flow can require a
sequential approach where solutions are first obtained for simpler problems
and used as the initial guess for more difficult problems.

.. GENERATED FROM PYTHON SOURCE LINES 90-97

.. code-block:: Python


    # disable the surface coverage equations, and turn off all gas and surface
    # chemistry.
    sim.surface.coverage_enabled = False
    surf_phase.set_multiplier(0.0)
    gas.set_multiplier(0.0)








.. GENERATED FROM PYTHON SOURCE LINES 98-100

solve the problem, refining the grid if needed, to determine the non-
reacting velocity and temperature distributions

.. GENERATED FROM PYTHON SOURCE LINES 100-102

.. code-block:: Python

    sim.solve(loglevel, auto=True)





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    ************ Solving on 6 point grid with energy equation enabled ************

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [6] point grid(s).
    grid refinement disabled.

    ******************** Solving with grid refinement enabled ********************

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [6] point grid(s).
    ##############################################################################
    Refining grid in flame.
        New points inserted after grid points 3 4 
        to resolve T spreadRate 
    ##############################################################################

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [8] point grid(s).
    ##############################################################################
    Refining grid in flame.
        New points inserted after grid points 5 6 
        to resolve spreadRate 
    ##############################################################################

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [10] point grid(s).
    no new points needed in flame




.. GENERATED FROM PYTHON SOURCE LINES 103-104

now turn on the surface coverage equations, and turn the chemistry on slowly

.. GENERATED FROM PYTHON SOURCE LINES 104-111

.. code-block:: Python

    sim.surface.coverage_enabled = True
    for mult in np.logspace(-5, 0, 6):
        surf_phase.set_multiplier(mult)
        gas.set_multiplier(mult)
        print('Multiplier =', mult)
        sim.solve(loglevel)





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    Multiplier = 1e-05

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [10] point grid(s).
    no new points needed in flame
    Multiplier = 0.0001

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [10] point grid(s).
    no new points needed in flame
    Multiplier = 0.001

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [10] point grid(s).
    no new points needed in flame
    Multiplier = 0.01

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [10] point grid(s).
    no new points needed in flame
    Multiplier = 0.1

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [10] point grid(s).
    no new points needed in flame
    Multiplier = 1.0

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [10] point grid(s).
    no new points needed in flame




.. GENERATED FROM PYTHON SOURCE LINES 112-113

At this point, we should have the solution for the hydrogen/air problem.

.. GENERATED FROM PYTHON SOURCE LINES 113-115

.. code-block:: Python

    sim.show()





.. rst-class:: sphx-glr-script-out

 .. code-block:: none



    >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> inlet <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

        Mass Flux:         0.06 kg/m^2/s 
        Temperature:        300 K 
        Mass Fractions: 
                          H2    0.003467 
                          O2      0.2311 
                          AR     0.01374 
                          N2      0.7516 



    >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> flame <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

        Pressure:   1.013e+05 Pa

    -------------------------------------------------------------------------------
              z    velocity  spreadRate           T      Lambda      eField 
    -------------------------------------------------------------------------------
              0     0.05335   4.314e-21         300        -0.4           0 
           0.02     0.05068      0.1333         300        -0.4           0 
           0.04     0.04268      0.2666         300        -0.4           0 
           0.06     0.02939      0.4001       300.6        -0.4           0 
           0.07     0.02105      0.4687       305.4        -0.4           0 
           0.08       0.013      0.5514       347.5        -0.4           0 
          0.085    0.009565      0.5931       426.2        -0.4           0 
           0.09    0.005933       0.577         569        -0.4           0 
          0.095    0.001988      0.3977       740.3        -0.4           0 
            0.1           0           0         900        -0.4           0 

    -------------------------------------------------------------------------------
              z          Uo          H2           H           O          O2 
    -------------------------------------------------------------------------------
              0           0    0.003467  -1.108e-21  -2.371e-20      0.2311 
           0.02           0    0.003467  -1.418e-27  -5.312e-21      0.2311 
           0.04           0    0.003462  -1.645e-28  -8.805e-22      0.2311 
           0.06           0    0.003399   1.746e-27   1.081e-20      0.2311 
           0.07           0    0.003197   6.826e-26   4.787e-19      0.2308 
           0.08           0    0.002515   3.923e-24   1.008e-17      0.2273 
          0.085           0    0.001917   7.589e-22   4.991e-17      0.2199 
           0.09           0    0.001239   2.823e-18   1.565e-15      0.2059 
          0.095           0    0.000594   1.554e-14   4.715e-13      0.1886 
            0.1           0   1.728e-05   2.758e-17   8.599e-16      0.1721 

    -------------------------------------------------------------------------------
              z          OH         H2O         HO2        H2O2           C 
    -------------------------------------------------------------------------------
              0   1.435e-25   6.854e-10   3.751e-16   4.325e-16   1.587e-21 
           0.02  -7.765e-24    3.18e-08   1.903e-14   2.209e-14   2.316e-29 
           0.04  -1.284e-24   1.404e-06   9.186e-13   1.074e-12   2.303e-37 
           0.06    1.33e-23   5.255e-05   3.755e-11   4.423e-11   7.661e-46 
           0.07   3.677e-22   0.0005515   4.291e-10   5.069e-10   4.848e-42 
           0.08   7.809e-21    0.005051   4.417e-09   5.105e-09   1.142e-33 
          0.085    3.04e-19     0.01333   1.441e-08   1.379e-08  -1.261e-39 
           0.09   2.408e-16      0.0275   3.848e-08   2.821e-08  -3.585e-34 
          0.095   2.228e-12     0.04362   7.898e-08   4.094e-08  -1.141e-26 
            0.1   1.912e-08     0.05837   7.921e-08   4.107e-08  -1.768e-19 

    -------------------------------------------------------------------------------
              z          CH         CH2      CH2(S)         CH3         CH4 
    -------------------------------------------------------------------------------
              0  -3.658e-21   1.055e-22   2.638e-38   1.472e-20  -1.583e-20 
           0.02  -1.664e-29   2.924e-29   1.202e-36   1.427e-20   -1.56e-20 
           0.04  -4.686e-37   7.259e-36   3.477e-41   1.393e-20  -1.532e-20 
           0.06   2.836e-44  -1.607e-37   1.211e-39   1.169e-20  -1.506e-20 
           0.07  -1.196e-41  -2.031e-36   6.661e-39   6.424e-21  -1.611e-20 
           0.08  -6.058e-33  -2.646e-33  -1.033e-38   1.126e-21  -2.838e-20 
          0.085  -1.125e-39   5.607e-35  -7.214e-38   2.004e-22   -5.29e-20 
           0.09  -1.107e-34   6.812e-33  -2.039e-36   1.885e-23  -9.856e-20 
          0.095  -6.596e-27    -5.3e-27  -2.684e-30   5.479e-23  -1.544e-19 
            0.1  -2.043e-19  -5.641e-27  -2.664e-30   4.437e-23  -2.059e-19 

    -------------------------------------------------------------------------------
              z          CO         CO2         HCO        CH2O       CH2OH 
    -------------------------------------------------------------------------------
              0  -2.083e-21   2.071e-24   -3.69e-29  -6.506e-22   1.274e-41 
           0.02  -5.303e-20   1.376e-22  -2.452e-27  -4.359e-20   8.552e-42 
           0.04  -5.305e-20   1.313e-22  -2.785e-34  -4.349e-20   1.127e-40 
           0.06  -5.398e-20  -2.075e-22  -3.891e-37   -4.09e-20   4.949e-39 
           0.07  -6.389e-20  -4.956e-21    2.28e-36  -4.261e-20    6.82e-38 
           0.08  -1.634e-19  -6.727e-20  -5.516e-32  -1.781e-19   1.516e-36 
          0.085  -3.721e-19  -2.349e-19   2.953e-34  -5.393e-19  -5.695e-36 
           0.09  -7.806e-19  -6.261e-19  -6.428e-30  -1.321e-18  -7.202e-32 
          0.095  -1.297e-18  -1.169e-18  -8.291e-26   -2.35e-18  -1.314e-25 
            0.1  -7.775e-19  -1.676e-18  -2.024e-25  -2.384e-18  -1.016e-18 

    -------------------------------------------------------------------------------
              z        CH3O       CH3OH         C2H        C2H2        C2H3 
    -------------------------------------------------------------------------------
              0    4.64e-30   6.682e-21  -3.294e-21  -8.998e-21    -3.9e-20 
           0.02   3.132e-28   6.682e-21  -2.942e-29  -9.312e-21  -4.079e-27 
           0.04   1.478e-26   6.684e-21   1.814e-38  -9.312e-21  -3.604e-34 
           0.06   4.901e-25   6.751e-21   5.921e-45  -9.315e-21  -2.843e-41 
           0.07   2.385e-24   6.963e-21    1.27e-48  -9.353e-21   4.176e-40 
           0.08   1.468e-24   7.295e-21   7.304e-46   -9.81e-21   3.682e-37 
          0.085   2.626e-25   7.038e-21   4.507e-42  -1.093e-20   6.579e-35 
           0.09  -1.228e-26   6.154e-21   2.377e-34   -1.34e-20   1.051e-31 
          0.095  -1.578e-25   4.998e-21   8.403e-27  -1.673e-20   6.372e-28 
            0.1   -1.52e-25   3.462e-21   1.821e-19   -4.01e-20   6.394e-28 

    -------------------------------------------------------------------------------
              z        C2H4        C2H5        C2H6        HCCO       CH2CO 
    -------------------------------------------------------------------------------
              0   1.598e-21   5.863e-41  -2.393e-21  -6.741e-41   9.752e-21 
           0.02   1.599e-21   4.191e-39  -2.393e-21  -2.782e-39   9.752e-21 
           0.04   1.606e-21   8.558e-38  -2.393e-21   1.796e-38   9.752e-21 
           0.06    2.01e-21    3.58e-36  -2.391e-21   7.806e-37    9.75e-21 
           0.07   7.461e-21    4.65e-35  -2.373e-21   1.011e-35   9.729e-21 
           0.08   7.421e-20   1.541e-33  -2.123e-21   1.714e-34   9.426e-21 
          0.085   2.342e-19   3.443e-31  -1.441e-21   2.372e-33   8.555e-21 
           0.09   5.709e-19   3.663e-28   1.705e-22   1.173e-31   6.452e-21 
          0.095   1.013e-18   7.342e-24   2.424e-21   1.854e-29   3.484e-21 
            0.1    1.01e-18   4.889e-19   3.985e-21   2.726e-29   5.762e-22 

    -------------------------------------------------------------------------------
              z       HCCOH          AR          N2 
    -------------------------------------------------------------------------------
              0   6.498e-21     0.01374      0.7516 
           0.02   6.498e-21     0.01374      0.7516 
           0.04   6.498e-21     0.01374      0.7516 
           0.06   6.498e-21     0.01374      0.7517 
           0.07   6.501e-21     0.01374      0.7517 
           0.08   6.561e-21     0.01374      0.7514 
          0.085   6.742e-21     0.01373      0.7511 
           0.09   7.202e-21     0.01374      0.7516 
          0.095   7.876e-21     0.01378      0.7534 
            0.1   8.539e-21     0.01382      0.7557 


    >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> surface <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

        Temperature:        900 K 
        Coverages: 
                       PT(S)    0.08832 
                        H(S)  1.352e-07 
                      H2O(S)  1.677e-05 
                       OH(S)    0.00567 
                       CO(S) -8.207e-22 
                      CO2(S) -1.057e-26 
                      CH3(S)  2.168e-26 
                     CH2(S)s  2.168e-26 
                       CH(S)  2.168e-26 
                        C(S)  6.863e-25 
                        O(S)      0.906 





.. GENERATED FROM PYTHON SOURCE LINES 116-117

Now switch the inlet to the methane/air composition.

.. GENERATED FROM PYTHON SOURCE LINES 117-125

.. code-block:: Python

    sim.inlet.X = comp2

    # set more stringent grid refinement criteria
    sim.set_refine_criteria(100.0, 0.15, 0.2, 0.0)

    # solve the problem for the final time
    sim.solve(loglevel)





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    Attempt Newton solution of steady-state problem.
    Newton steady-state solve failed.

    Attempt 10 timesteps.
    Final timestep info: dt= 0.0003844  log(ss)= -0.8978   

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve failed.

    Attempt 10 timesteps.
    Final timestep info: dt= 0.01478    log(ss)= -0.9305   

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve failed.

    Attempt 10 timesteps.
    Final timestep info: dt= 0.5682     log(ss)= -0.3058   

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [10] point grid(s).
    ##############################################################################
    Refining grid in flame.
        New points inserted after grid points 0 1 2 3 4 5 6 7 8 
        to resolve AR C2H6 CH2O CH4 CO CO2 H2O N2 O2 OH T spreadRate velocity 
    ##############################################################################

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [19] point grid(s).
    ##############################################################################
    Refining grid in flame.
        New points inserted after grid points 9 10 11 12 13 14 15 16 17 
        to resolve AR CH2O CH4 CO CO2 H2O N2 O2 OH T spreadRate velocity 
    ##############################################################################

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [28] point grid(s).
    ##############################################################################
    Refining grid in flame.
        New points inserted after grid points 24 25 26 
        to resolve OH spreadRate velocity 
    ##############################################################################

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [31] point grid(s).
    ##############################################################################
    Refining grid in flame.
        New points inserted after grid points 29 
        to resolve OH 
    ##############################################################################

    Attempt Newton solution of steady-state problem.
    Newton steady-state solve succeeded.

    Problem solved on [32] point grid(s).
    no new points needed in flame




.. GENERATED FROM PYTHON SOURCE LINES 126-127

show the solution

.. GENERATED FROM PYTHON SOURCE LINES 127-129

.. code-block:: Python

    sim.show()





.. rst-class:: sphx-glr-script-out

 .. code-block:: none



    >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> inlet <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

        Mass Flux:         0.06 kg/m^2/s 
        Temperature:        300 K 
        Mass Fractions: 
                          O2      0.2204 
                         CH4     0.04998 
                          AR      0.0131 
                          N2      0.7166 



    >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> flame <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

        Pressure:   1.013e+05 Pa

    -------------------------------------------------------------------------------
              z    velocity  spreadRate           T      Lambda      eField 
    -------------------------------------------------------------------------------
              0     0.05304  -1.269e-24         300     -0.3867           0 
           0.01     0.05239     0.06444         300     -0.3867           0 
           0.02     0.05046      0.1289         300     -0.3867           0 
           0.03     0.04724      0.1933         300     -0.3867           0 
           0.04     0.04273      0.2578         300     -0.3867           0 
           0.05     0.03693      0.3222         300     -0.3867           0 
           0.06     0.02984      0.3867         300     -0.3867           0 
          0.065     0.02582      0.4189       300.1     -0.3867           0 
           0.07     0.02149      0.4513       300.4     -0.3867           0 
          0.075     0.01693      0.4842       302.2     -0.3867           0 
         0.0775     0.01465      0.5017       305.6     -0.3867           0 
           0.08     0.01247       0.521       314.3     -0.3867           0 
        0.08125     0.01147      0.5318       322.7     -0.3867           0 
         0.0825     0.01056      0.5434       335.5     -0.3867           0 
        0.08375    0.009736      0.5556         354     -0.3867           0 
          0.085    0.008971      0.5673       378.9     -0.3867           0 
        0.08625    0.008228      0.5766       410.4     -0.3867           0 
         0.0875    0.007466      0.5813       447.9     -0.3867           0 
        0.08875    0.006649      0.5789       490.3     -0.3867           0 
           0.09    0.005765      0.5668         536     -0.3867           0 
        0.09125    0.004822      0.5433       583.7     -0.3867           0 
         0.0925    0.003848      0.5069       632.1     -0.3867           0 
        0.09375    0.002883      0.4567       680.2     -0.3867           0 
          0.095    0.001979      0.3925       727.3     -0.3867           0 
        0.09625    0.001186      0.3144         773     -0.3867           0 
        0.09688   0.0008485      0.2702       795.2     -0.3867           0 
         0.0975   0.0005586      0.2226         817     -0.3867           0 
        0.09813   0.0003229      0.1717       838.4     -0.3867           0 
        0.09875   0.0001473      0.1176       859.4     -0.3867           0 
        0.09938   3.774e-05     0.06034       879.9     -0.3867           0 
        0.09969   9.583e-06     0.03055         890     -0.3867           0 
            0.1   4.186e-68           0         900     -0.3867           0 

    -------------------------------------------------------------------------------
              z          Uo          H2           H           O          O2 
    -------------------------------------------------------------------------------
              0           0  -5.538e-18   6.016e-21   1.841e-20      0.2204 
           0.01           0  -4.315e-17     5.6e-27  -3.408e-22      0.2204 
           0.02           0  -3.332e-16   5.129e-33  -1.163e-22      0.2204 
           0.03           0  -2.497e-15  -1.464e-30  -3.808e-23      0.2204 
           0.04           0  -1.775e-14   2.673e-31  -1.169e-23      0.2204 
           0.05           0  -1.165e-13  -1.313e-31  -3.079e-24      0.2204 
           0.06           0  -6.821e-13   1.933e-31   7.151e-24      0.2204 
          0.065           0  -1.775e-12   1.029e-31   1.311e-22      0.2203 
           0.07           0  -4.674e-12   4.836e-29   1.776e-21      0.2203 
          0.075           0  -1.154e-11   1.214e-26   2.311e-20        0.22 
         0.0775           0   -1.77e-11   2.161e-25   1.243e-19      0.2194 
           0.08           0  -2.654e-11   3.667e-24   5.612e-19      0.2177 
        0.08125           0  -3.209e-11   1.466e-23   1.152e-18       0.216 
         0.0825           0  -3.841e-11   5.539e-23   2.097e-18      0.2132 
        0.08375           0  -4.539e-11   1.904e-22   3.184e-18       0.209 
          0.085           0  -5.287e-11   6.075e-22   3.869e-18      0.2032 
        0.08625           0  -6.068e-11   2.312e-21   3.859e-18      0.1957 
         0.0875           0  -6.861e-11   1.987e-20    3.72e-18      0.1864 
        0.08875           0  -7.652e-11   2.864e-19   6.284e-18      0.1757 
           0.09           0  -8.428e-11   3.575e-18   3.116e-17      0.1638 
        0.09125           0  -9.182e-11   3.271e-17   2.062e-16      0.1512 
         0.0925           0   -9.91e-11   2.154e-16   1.207e-15      0.1382 
        0.09375           0  -1.061e-10   1.033e-15   5.892e-15      0.1251 
          0.095           0  -1.129e-10   3.759e-15   2.403e-14      0.1122 
        0.09625           0  -1.198e-10   1.124e-14   8.456e-14     0.09946 
        0.09688           0  -1.238e-10   1.876e-14    1.54e-13     0.09326 
         0.0975           0  -1.283e-10   3.082e-14   2.753e-13     0.08715 
        0.09813           0   -1.34e-10   5.036e-14   4.869e-13     0.08115 
        0.09875           0  -1.417e-10   8.257e-14   8.591e-13     0.07525 
        0.09938           0   -1.53e-10   1.371e-13   1.524e-12     0.06945 
        0.09969           0  -1.611e-10   1.809e-13   2.006e-12     0.06659 
            0.1           0   -1.71e-10   3.422e-15     3.9e-14     0.06376 

    -------------------------------------------------------------------------------
              z          OH         H2O         HO2        H2O2           C 
    -------------------------------------------------------------------------------
              0  -4.101e-26   1.356e-14   2.236e-20  -1.543e-21  -1.406e-20 
           0.01  -7.208e-25   3.307e-13   5.916e-19   1.114e-21  -4.537e-28 
           0.02   -7.28e-26   7.976e-12   8.973e-19   1.454e-19  -1.413e-35 
           0.03  -2.377e-26   1.858e-10    8.46e-18   3.882e-18  -4.135e-43 
           0.04  -4.378e-27   4.068e-09   1.873e-16   9.302e-17  -1.102e-50 
           0.05   1.131e-24   8.112e-08   4.031e-15   2.023e-15  -2.563e-58 
           0.06   3.702e-22   1.413e-06   7.596e-14    3.84e-14  -5.061e-66 
          0.065   1.493e-20   8.652e-06    4.98e-13   2.534e-13  -3.157e-68 
           0.07   6.964e-19   5.707e-05   3.531e-12   1.808e-12  -2.876e-66 
          0.075   2.467e-17   0.0003335   2.212e-11    1.14e-11   1.278e-59 
         0.0775   1.588e-16   0.0008499   5.923e-11   3.061e-11   2.917e-58 
           0.08   9.223e-16    0.002149   1.584e-10     8.2e-11   2.131e-57 
        0.08125   2.076e-15    0.003384    2.58e-10   1.335e-10  -1.203e-57 
         0.0825   4.287e-15    0.005243   4.153e-10   2.143e-10  -1.203e-56 
        0.08375   7.865e-15    0.007865   6.487e-10   3.331e-10  -3.706e-56 
          0.085    1.28e-14     0.01132   9.735e-10   4.964e-10   6.377e-56 
        0.08625   1.888e-14     0.01558   1.398e-09    7.06e-10   1.199e-50 
         0.0875   2.616e-14     0.02052   1.921e-09   9.587e-10   5.362e-51 
        0.08875   3.535e-14     0.02596   2.534e-09   1.247e-09  -2.877e-44 
           0.09   4.856e-14     0.03173   3.223e-09   1.561e-09  -2.053e-38 
        0.09125   7.209e-14     0.03764   3.972e-09    1.89e-09   2.465e-38 
         0.0925   1.256e-13     0.04356   4.768e-09   2.224e-09   1.038e-36 
        0.09375   2.696e-13     0.04939   5.604e-09   2.555e-09   1.423e-30 
          0.095   6.824e-13     0.05507   6.475e-09    2.87e-09   5.921e-36 
        0.09625   1.882e-12     0.06056   7.385e-09   3.149e-09   3.808e-30 
        0.09688   3.208e-12     0.06322   7.852e-09   3.263e-09   7.173e-35 
         0.0975   5.507e-12     0.06584   8.322e-09   3.354e-09   5.037e-30 
        0.09813   9.543e-12      0.0684   8.782e-09   3.417e-09   4.925e-28 
        0.09875   1.677e-11     0.07091   9.206e-09   3.449e-09  -5.483e-28 
        0.09938    3.07e-11     0.07337   9.539e-09   3.455e-09   2.066e-29 
        0.09969   1.934e-10     0.07458    9.63e-09   3.453e-09   6.738e-25 
            0.1   1.749e-08     0.07579   9.628e-09   3.452e-09   6.902e-22 

    -------------------------------------------------------------------------------
              z          CH         CH2      CH2(S)         CH3         CH4 
    -------------------------------------------------------------------------------
              0   1.719e-34  -3.497e-21  -3.841e-36  -5.887e-21     0.04998 
           0.01   1.309e-42  -2.303e-27   -9.21e-35  -5.165e-21     0.04998 
           0.02  -1.402e-48  -1.463e-33   2.413e-40    2.87e-19     0.04998 
           0.03  -3.943e-47   7.392e-39   3.814e-40   7.376e-18     0.04998 
           0.04  -1.427e-47   4.273e-38   1.971e-39   1.632e-16     0.04998 
           0.05    9.23e-46   8.805e-37   4.182e-38   3.271e-15     0.04998 
           0.06  -6.268e-45   1.188e-33   6.165e-35   5.726e-14     0.04998 
          0.065   1.206e-43   2.919e-31   1.517e-32   3.521e-13     0.04998 
           0.07   9.606e-43   9.019e-29    4.69e-30   2.333e-12     0.04995 
          0.075  -5.014e-41   1.878e-26   9.808e-28   1.373e-11     0.04984 
         0.0775  -1.541e-40   3.117e-25   1.641e-26   3.544e-11     0.04962 
           0.08  -3.674e-40   4.671e-24   2.503e-25   9.187e-11     0.04909 
        0.08125  -5.149e-39   1.698e-23   9.245e-25   1.492e-10     0.04858 
         0.0825   6.622e-39    5.62e-23   3.127e-24   2.414e-10     0.04781 
        0.08375   1.143e-38   1.619e-22   9.249e-24   3.852e-10     0.04672 
          0.085   1.826e-37   4.034e-22   2.367e-23   6.033e-10     0.04524 
        0.08625   2.717e-37   8.907e-22   5.345e-23   9.276e-10     0.04339 
         0.0875  -7.111e-37   1.819e-21   1.104e-22   1.403e-09     0.04118 
        0.08875   1.514e-36   3.606e-21   2.181e-22   2.096e-09     0.03868 
           0.09   1.957e-35   7.335e-21   4.338e-22   3.104e-09     0.03596 
        0.09125   2.881e-34   1.652e-20   9.355e-22   4.575e-09     0.03309 
         0.0925   5.029e-33    4.51e-20   2.394e-21   6.748e-09     0.03016 
        0.09375   8.511e-32   1.548e-19   7.586e-21   1.002e-08     0.02722 
          0.095   1.362e-30   6.331e-19   2.839e-20   1.511e-08     0.02431 
        0.09625   2.046e-29   2.883e-18   1.175e-19   2.335e-08     0.02146 
        0.09688   8.137e-29   6.467e-18   2.506e-19   2.969e-08     0.02006 
         0.0975   3.266e-28    1.48e-17    5.44e-19   3.815e-08     0.01869 
        0.09813   1.347e-27   3.477e-17   1.211e-18   4.972e-08     0.01734 
        0.09875   5.817e-27   8.472e-17   2.787e-18   6.606e-08       0.016 
        0.09938   1.738e-26   2.242e-16   6.944e-18       9e-08     0.01469 
        0.09969  -1.275e-23    2.83e-15    8.77e-17   1.066e-07     0.01405 
            0.1  -1.275e-23   2.829e-15   8.767e-17   1.066e-07      0.0134 

    -------------------------------------------------------------------------------
              z          CO         CO2         HCO        CH2O       CH2OH 
    -------------------------------------------------------------------------------
              0   3.827e-19    2.25e-16   8.428e-28   1.044e-20  -1.952e-45 
           0.01   1.015e-17   7.764e-15   2.917e-26   3.644e-19  -6.861e-44 
           0.02   2.592e-16   2.648e-13    8.93e-33   2.879e-18  -5.569e-43 
           0.03   6.553e-15   8.717e-12  -1.128e-38   8.649e-17  -6.323e-42 
           0.04   1.559e-13   2.694e-10   8.354e-37   2.694e-15  -2.476e-41 
           0.05   3.373e-12   7.566e-09   5.544e-34   7.651e-14   2.179e-38 
           0.06   6.374e-11   1.851e-07   2.459e-31   1.892e-12   1.416e-35 
          0.065   4.188e-10   1.529e-06   1.421e-29   1.578e-11   2.055e-33 
           0.07   2.976e-09   1.382e-05   1.546e-27   1.442e-10   7.012e-31 
          0.075   1.869e-08   0.0001091   2.607e-25   1.149e-09   1.925e-28 
         0.0775   5.007e-08    0.000344   4.908e-24   3.647e-09   3.933e-27 
           0.08   1.339e-07    0.001094   8.891e-23   1.166e-08   7.286e-26 
        0.08125    2.18e-07    0.001963   3.611e-22   2.098e-08   2.944e-25 
         0.0825   3.502e-07    0.003486   1.351e-21   3.729e-08    1.08e-24 
        0.08375   5.455e-07    0.005963   4.445e-21   6.375e-08   3.415e-24 
          0.085   8.154e-07    0.009689   1.303e-20   1.033e-07    9.29e-24 
        0.08625   1.165e-06     0.01487   3.661e-20   1.578e-07   2.284e-23 
         0.0875   1.589e-06     0.02156   1.095e-19   2.274e-07   5.671e-23 
        0.08875   2.079e-06     0.02965   3.677e-19   3.105e-07    1.68e-22 
           0.09   2.619e-06     0.03891   1.322e-18   4.042e-07   6.525e-22 
        0.09125   3.191e-06     0.04904   4.723e-18    5.05e-07   2.982e-21 
         0.0925   3.781e-06     0.05971    1.61e-17   6.094e-07   1.361e-20 
        0.09375   4.377e-06     0.07065   5.212e-17   7.144e-07   5.671e-20 
          0.095   4.967e-06     0.08163   1.625e-16   8.169e-07   2.098e-19 
        0.09625   5.545e-06     0.09247   4.954e-16   9.137e-07   6.981e-19 
        0.09688   5.828e-06      0.0978    8.67e-16   9.582e-07   1.243e-18 
         0.0975   6.107e-06      0.1031   1.511e-15   9.991e-07   2.194e-18 
        0.09813   6.381e-06      0.1082   2.629e-15   1.035e-06   3.878e-18 
        0.09875   6.651e-06      0.1133   4.575e-15   1.064e-06   6.977e-18 
        0.09938   6.916e-06      0.1184   7.962e-15   1.084e-06   1.278e-17 
        0.09969   7.047e-06      0.1208    7.25e-15   1.088e-06  -3.031e-17 
            0.1   7.177e-06      0.1233   7.248e-15   1.088e-06   -3.03e-17 

    -------------------------------------------------------------------------------
              z        CH3O       CH3OH         C2H        C2H2        C2H3 
    -------------------------------------------------------------------------------
              0  -3.669e-33  -6.849e-22   1.438e-20  -1.056e-21    2.56e-19 
           0.01  -1.289e-31   -4.08e-22   2.814e-28   9.575e-22   5.857e-26 
           0.02  -2.372e-32    9.31e-21   5.309e-36   9.554e-22   1.292e-32 
           0.03   2.015e-29   3.381e-19   9.397e-44   8.889e-22   2.674e-39 
           0.04   9.897e-27   1.077e-17    1.51e-51  -1.054e-21   4.163e-41 
           0.05   4.263e-24   3.114e-16  -1.607e-55  -5.262e-20   2.316e-38 
           0.06   1.401e-21   7.836e-15  -2.928e-51  -1.242e-18   9.854e-36 
          0.065   5.655e-20   6.638e-14  -9.303e-49   -9.79e-18   5.204e-34 
           0.07   2.652e-18   6.163e-13  -3.784e-46  -8.423e-17   4.212e-32 
          0.075   9.678e-17    4.99e-12  -1.164e-43  -6.342e-16   5.119e-30 
         0.0775   6.582e-16   1.602e-11  -3.017e-42  -1.932e-15   9.259e-29 
           0.08   4.362e-15   5.189e-11  -1.096e-40  -5.922e-15   1.881e-27 
        0.08125   1.103e-14   9.402e-11  -8.346e-40   -1.04e-14   8.813e-27 
         0.0825   2.673e-14   1.685e-10  -7.721e-39  -1.807e-14   4.123e-26 
        0.08375   5.961e-14   2.903e-10    -8.3e-38  -3.029e-14    1.86e-25 
          0.085   1.196e-13   4.738e-10  -9.725e-37   -4.83e-14   8.187e-25 
        0.08625   2.143e-13   7.287e-10   -1.14e-35  -7.294e-14     3.7e-24 
         0.0875    3.46e-13   1.057e-09  -1.233e-34  -1.043e-13   1.798e-23 
        0.08875   5.124e-13   1.451e-09  -1.195e-33  -1.418e-13   9.244e-23 
           0.09   7.108e-13   1.898e-09  -1.132e-32  -1.842e-13   4.695e-22 
        0.09125   9.346e-13   2.381e-09  -1.244e-31  -2.302e-13   2.221e-21 
         0.0925   1.159e-12   2.882e-09  -1.528e-30  -2.784e-13   9.562e-21 
        0.09375    1.34e-12   3.382e-09  -1.719e-29  -3.275e-13   3.755e-20 
          0.095   1.455e-12   3.862e-09  -1.609e-28  -3.766e-13   1.369e-19 
        0.09625   1.555e-12   4.303e-09  -1.255e-27  -4.252e-13   4.752e-19 
        0.09688   1.645e-12   4.502e-09  -3.333e-27  -4.494e-13   8.852e-19 
         0.0975   1.771e-12   4.684e-09  -8.597e-27  -4.735e-13   1.645e-18 
        0.09813   1.953e-12   4.845e-09   -2.16e-26  -4.979e-13   3.069e-18 
        0.09875   2.218e-12   4.981e-09  -5.386e-26  -5.228e-13   5.781e-18 
        0.09938   2.606e-12   5.085e-09  -1.003e-26  -5.491e-13   9.995e-18 
        0.09969   2.904e-12    5.12e-09  -4.373e-23   -5.63e-13  -2.215e-16 
            0.1   2.904e-12   5.119e-09  -1.073e-21  -5.629e-13  -2.215e-16 

    -------------------------------------------------------------------------------
              z        C2H4        C2H5        C2H6        HCCO       CH2CO 
    -------------------------------------------------------------------------------
              0   7.906e-22   1.108e-33  -1.007e-19  -4.587e-42   2.653e-28 
           0.01   1.391e-21   4.039e-32  -4.656e-20  -1.005e-40   3.889e-28 
           0.02   2.144e-20   4.864e-32   1.597e-18   -3.43e-41   5.309e-27 
           0.03   6.692e-19   6.228e-32   5.878e-17  -1.112e-41   1.942e-25 
           0.04   2.031e-17   3.226e-31   1.937e-15  -2.562e-42   6.997e-24 
           0.05   5.607e-16   5.085e-30   5.791e-14  -4.363e-43   2.292e-22 
           0.06   1.349e-14   9.045e-29   1.508e-12   2.502e-41   6.533e-21 
          0.065   1.097e-13   1.128e-27   1.317e-11   5.345e-39   6.195e-20 
           0.07   9.761e-13   1.166e-25   1.262e-10   2.389e-36   6.471e-19 
          0.075   7.589e-12   2.653e-23   1.054e-09   7.639e-34   5.866e-18 
         0.0775   2.366e-11   5.323e-22   3.464e-09   1.757e-32    2.05e-17 
           0.08   7.432e-11   9.501e-21   1.151e-08   3.844e-31   7.272e-17 
        0.08125   1.324e-10   3.704e-20   2.119e-08   1.752e-30   1.391e-16 
         0.0825   2.332e-10   1.287e-19   3.865e-08   7.462e-30   2.638e-16 
        0.08375   3.958e-10   3.772e-19   6.779e-08   2.786e-29   4.797e-16 
          0.085   6.381e-10   9.368e-19   1.126e-07   8.936e-29   8.223e-16 
        0.08625   9.722e-10   2.089e-18   1.761e-07   2.485e-28    1.32e-15 
         0.0875     1.4e-09   4.652e-18   2.593e-07   6.156e-28   1.987e-15 
        0.08875   1.913e-09   1.166e-17   3.611e-07   1.402e-27   2.815e-15 
           0.09   2.496e-09   3.432e-17   4.783e-07   3.025e-27   3.782e-15 
        0.09125   3.129e-09   1.112e-16   6.066e-07   6.251e-27   4.854e-15 
         0.0925   3.792e-09    3.67e-16   7.415e-07   1.163e-26   5.993e-15 
        0.09375   4.466e-09    1.19e-15   8.785e-07   1.266e-26   7.166e-15 
          0.095   5.138e-09    3.78e-15   1.013e-06  -4.391e-26   8.344e-15 
        0.09625   5.791e-09   1.199e-14   1.141e-06  -4.384e-25   9.507e-15 
        0.09688   6.104e-09    2.18e-14   1.201e-06  -1.068e-24   1.007e-14 
         0.0975   6.404e-09   4.005e-14   1.256e-06  -2.446e-24   1.063e-14 
        0.09813   6.686e-09   7.498e-14   1.305e-06  -5.349e-24   1.117e-14 
        0.09875   6.939e-09   1.443e-13   1.347e-06  -1.148e-23   1.168e-14 
        0.09938   7.144e-09   2.973e-13   1.375e-06    -4.1e-23   1.215e-14 
        0.09969   7.216e-09   1.105e-12   1.381e-06  -3.455e-21   1.234e-14 
            0.1   7.214e-09   1.105e-12   1.381e-06  -3.454e-21   1.234e-14 

    -------------------------------------------------------------------------------
              z       HCCOH          AR          N2 
    -------------------------------------------------------------------------------
              0   2.376e-31      0.0131      0.7166 
           0.01   1.448e-31      0.0131      0.7166 
           0.02   -3.55e-30      0.0131      0.7166 
           0.03  -1.454e-28      0.0131      0.7166 
           0.04  -5.254e-27      0.0131      0.7166 
           0.05  -1.721e-25      0.0131      0.7166 
           0.06  -4.906e-24      0.0131      0.7166 
          0.065  -4.652e-23      0.0131      0.7166 
           0.07  -4.859e-22      0.0131      0.7166 
          0.075  -4.405e-21      0.0131      0.7166 
         0.0775   -1.54e-20      0.0131      0.7167 
           0.08   -5.46e-20     0.01311      0.7169 
        0.08125  -1.044e-19     0.01311       0.717 
         0.0825  -1.981e-19     0.01311      0.7172 
        0.08375  -3.602e-19     0.01311      0.7173 
          0.085  -6.174e-19     0.01312      0.7174 
        0.08625  -9.914e-19     0.01312      0.7174 
         0.0875  -1.492e-18     0.01311      0.7172 
        0.08875  -2.114e-18     0.01311      0.7169 
           0.09   -2.84e-18      0.0131      0.7165 
        0.09125  -3.645e-18     0.01309      0.7159 
         0.0925    -4.5e-18     0.01308      0.7153 
        0.09375  -5.381e-18     0.01306      0.7145 
          0.095  -6.266e-18     0.01305      0.7138 
        0.09625   -7.14e-18     0.01304       0.713 
        0.09688  -7.568e-18     0.01303      0.7126 
         0.0975   -7.99e-18     0.01302      0.7122 
        0.09813  -8.406e-18     0.01302      0.7119 
        0.09875  -8.813e-18     0.01301      0.7115 
        0.09938  -9.213e-18       0.013      0.7111 
        0.09969   -9.41e-18       0.013      0.7109 
            0.1  -9.409e-18       0.013      0.7108 


    >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> surface <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

        Temperature:        900 K 
        Coverages: 
                       PT(S)     0.1251 
                        H(S)  2.668e-07 
                      H2O(S)   3.12e-05 
                       OH(S)   0.007566 
                       CO(S)  3.863e-05 
                      CO2(S)  4.183e-10 
                      CH3(S)  3.124e-09 
                     CH2(S)s  3.124e-09 
                       CH(S)  3.124e-09 
                        C(S)  1.373e-07 
                        O(S)     0.8673 





.. GENERATED FROM PYTHON SOURCE LINES 130-132

Save the full solution to HDF or YAML container files. The ``restore`` method can be
used to restore or restart a simulation from a solution stored in this form.

.. GENERATED FROM PYTHON SOURCE LINES 132-144

.. code-block:: Python

    if "native" in ct.hdf_support():
        filename = "catalytic_combustion.h5"
    else:
        filename = "catalytic_combustion.yaml"
    sim.save(filename, "soln1", description="catalytic combustion example",
             overwrite=True)

    # save selected solution components in a CSV file for plotting in Excel or MATLAB.
    sim.save('catalytic_combustion.csv', basis="mole", overwrite=True)

    sim.show_stats()





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    Statistics:

     Grid   Timesteps  Functions      Time  Jacobians      Time
        8           0         87    0.0268          2    0.0418
       10           0          6    0.0025          1    0.0291
       12           0          6    0.0032          1    0.0369
       12           0         10    0.0055          1    0.0374
       12           0          6    0.0033          1    0.0383
       12           0         12    0.0065          1    0.0382
       12           0         20    0.0106          1    0.0381
       12           0         14    0.0077          1    0.0383
       12           0          6    0.0032          1    0.0381
       12          30        394    0.2139         12    0.4571
       21           0         32    0.0327          2    0.1500
       30           0         32    0.0479          2    0.2269
       33           0         23    0.0380          1    0.1229
       34           0         19    0.0323          1    0.1281




.. GENERATED FROM PYTHON SOURCE LINES 145-147

Temperature Profile
-------------------

.. GENERATED FROM PYTHON SOURCE LINES 147-153

.. code-block:: Python

    fig, ax = plt.subplots()
    ax.plot(sim.grid, sim.T, color='C3')
    ax.set_ylabel('heat release rate [MW/m³]')
    ax.set(xlabel='distance from inlet [m]')
    plt.show()




.. image-sg:: /examples/python/onedim/images/sphx_glr_catalytic_combustion_001.png
   :alt: catalytic combustion
   :srcset: /examples/python/onedim/images/sphx_glr_catalytic_combustion_001.png, /examples/python/onedim/images/sphx_glr_catalytic_combustion_001_2_00x.png 2.00x
   :class: sphx-glr-single-img





.. GENERATED FROM PYTHON SOURCE LINES 154-156

Major Species Profiles
----------------------

.. GENERATED FROM PYTHON SOURCE LINES 156-163

.. code-block:: Python

    fig, ax = plt.subplots()
    major = ('O2', 'CH4', 'H2O', 'CO2')
    states = sim.to_array()
    ax.plot(states.grid, states(*major).X, label=major)
    ax.set(xlabel='distance from inlet [m]', ylabel='mole fractions')
    ax.legend()
    plt.show()



.. image-sg:: /examples/python/onedim/images/sphx_glr_catalytic_combustion_002.png
   :alt: catalytic combustion
   :srcset: /examples/python/onedim/images/sphx_glr_catalytic_combustion_002.png, /examples/python/onedim/images/sphx_glr_catalytic_combustion_002_2_00x.png 2.00x
   :class: sphx-glr-single-img






.. rst-class:: sphx-glr-timing

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


.. _sphx_glr_download_examples_python_onedim_catalytic_combustion.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: catalytic_combustion.ipynb <catalytic_combustion.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: catalytic_combustion.py <catalytic_combustion.py>`

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: catalytic_combustion.zip <catalytic_combustion.zip>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_
