/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      hTCReactionsDPLR; // Reactions used in the DPLR solver
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

species
(
    N2
    //O2
    //NO
    N2+
    //O2+
    //NO+
    N
    //O
    N+
    //O+
    e-
);

vibTempAssociativity (0);

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
/* REACTIONS NO 1 TO 15 AND 18 TO 32 -> DISSOCIATION
   REACTIONS NO 16 AND 17            -> EXCHANGE
   REACTION NO 33                    -> ELECTRON IMPACT DISSOCIATION OF N2
   REACTIONS NO 34 AND 35            -> ELECTRON IMPACT IONISATION
   REACTIONS NO 36 TO 38             -> ASSOCIATIVE IONISATION
   REACTIONS NO 39 TO 49             -> CHARGE EXCHANGE                      */
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// Units
// A: m^3 kmol^-1 s^-1
// Ta: K
  
reactions
{
	// Reaction no 9
	nitrogenReaction
	{
		type     reversibleArrheniusReaction;
		reaction "N2 + N2 = 2N + N2";
		controlT dissociation;
		A        7.0e18;
		beta     -1.6;
		Ta       113220;
		
		ni       (1e20 1e21 1e22 1e23 1e24 1e25);
		A0       (3.4907 2.0723 1.606 1.5351 1.4766 1.4766);
		A1       (0.83133 1.3897 1.5732 1.6061 1.6291 1.6291);
		A2       (4.0978 2.0617 1.3923 1.2993 1.2153 1.2153);
		A3       (-12.728 -11.828 -11.533 -11.494 -11.457 -11.457);
		A4       (0.07487 0.015105 -0.004543 -0.00698 -0.009444 -0.009444);
	}
	
	// Reaction no 10
	nitrogenAtomicNitrogenReaction
	{
		type     reversibleArrheniusReaction;
		reaction "N2 + N = 2N + N";
		controlT dissociation;
		A        3.0e19;
		beta     -1.6;
		Ta       113220;
		
		ni       (1e20 1e21 1e22 1e23 1e24 1e25);
		A0       (3.4907 2.0723 1.606 1.5351 1.4766 1.4766);
		A1       (0.83133 1.3897 1.5732 1.6061 1.6291 1.6291);
		A2       (4.0978 2.0617 1.3923 1.2993 1.2153 1.2153);
		A3       (-12.728 -11.828 -11.533 -11.494 -11.457 -11.457);
		A4       (0.07487 0.015105 -0.004543 -0.00698 -0.009444 -0.009444);
	}

// Reaction no 18
	nitrogenNitrogenCationReaction
	{
		type     reversibleArrheniusReaction;
		reaction "N2 + N2+ = 2N + N2+";
		controlT dissociation;
		A        7.0e18;
		beta     -1.6;
		Ta       113220;
		
		ni       (1e20 1e21 1e22 1e23 1e24 1e25);
		A0       (3.4907 2.0723 1.606 1.5351 1.4766 1.4766);
		A1       (0.83133 1.3897 1.5732 1.6061 1.6291 1.6291);
		A2       (4.0978 2.0617 1.3923 1.2993 1.2153 1.2153);
		A3       (-12.728 -11.828 -11.533 -11.494 -11.457 -11.457);
		A4       (0.07487 0.015105 -0.004543 -0.00698 -0.009444 -0.009444);
	}

// Reaction no 21
	nitrogenAtomicNitrogenCationReaction
	{
		type     reversibleArrheniusReaction;
		reaction "N2 + N+ = 2N + N+";
		controlT dissociation;
		A        3.0e19;
		beta     -1.6;
		Ta       113220;
		
		ni       (1e20 1e21 1e22 1e23 1e24 1e25);
		A0       (3.4907 2.0723 1.606 1.5351 1.4766 1.4766);
		A1       (0.83133 1.3897 1.5732 1.6061 1.6291 1.6291);
		A2       (4.0978 2.0617 1.3923 1.2993 1.2153 1.2153);
		A3       (-12.728 -11.828 -11.533 -11.494 -11.457 -11.457);
		A4       (0.07487 0.015105 -0.004543 -0.00698 -0.009444 -0.009444);
	}

// Reaction no 33
	nitrogenElectronImpactDissociationReaction
	{
		type     reversibleArrheniusReaction;
		reaction "N2 + e- = 2N + e-";
		controlT impactDissociation;
		A        3.0e21;
		beta     -1.6;
		Ta       113220;
		
		ni       (1e20 1e21 1e22 1e23 1e24 1e25);
		A0       (3.4907 2.0723 1.606 1.5351 1.4766 1.4766);
		A1       (0.83133 1.3897 1.5732 1.6061 1.6291 1.6291);
		A2       (4.0978 2.0617 1.3923 1.2993 1.2153 1.2153);
		A3       (-12.728 -11.828 -11.533 -11.494 -11.457 -11.457);
		A4       (0.07487 0.015105 -0.004543 -0.00698 -0.009444 -0.009444);
	}
	
	// Reaction no 34
	atomicNitrogenElectronImpactIonisationReaction
	{
		type     reversibleArrheniusReaction;
		reaction "N + e- = N+ + 2e-";
		controlT impactIonisation;
		A        2.5e31;
		beta     -3.82;
		Ta       168600;
		
		ni       (1e20 1e21 1e22 1e23 1e24 1e25);
		A0       (-1.9094 -1.2002 -0.96709 -0.93184 -0.9026 -0.9026 );
		A1       (-3.0267 -3.3059 -3.3976 -3.414 -3.4255 -3.4255 );
		A2       (-3.6935 -2.6755 -2.3408 -2.2946 -2.2526 -2.2526 );
		A3       (-16.044 -16.494 -16.642 -16.661 -16.679 -16.679 );
		A4       (-0.050183 -0.020301 -0.010477 -0.009269 -0.008037 -0.008037);
	}

// Reaction no 37
	atomicNitrogenAtomicNitrogenAssociativeIonisationReaction
	{
		type     reversibleArrheniusReaction;
		reaction "N + N = N2+ + e-";
		controlT recombinationIonisation;
		A        4.4e04;
		beta     1.5;
		Ta       67500;
		
		ni       (1e20 1e21 1e22 1e23 1e24 1e25);
		A0       (-4.3785 -2.9601 -2.4938 -2.4229 -2.3644 -2.3644);
		A1       (-4.2726 -4.831 -5.0145 -5.0474 -5.0704 -5.0704);
		A2       (-7.8709 -5.8348 -5.1654 -5.0724 -4.9885 -4.9885);
		A3       (-4.4628 -5.3621 -5.6577 -5.6961 -5.7332 -5.7332);
		A4       (-0.12402 -0.064252 -0.044602 -0.042167 -0.039703 -0.039703);
	}

// Reaction no 49
	atomicNitrogenCationNitrogenChargeExchangeReaction
	{
		type     reversibleArrheniusReaction;
		reaction "N+ + N2 = N2+ + N";
		controlT chargeExchange;
		A        1.0e09;
		beta     0.5;
		Ta       12200;
		
		ni       (1e20 1e21 1e22 1e23 1e24 1e25);
		A0       ("no data");
		A1       ();
		A2       ();
		A3       ();
		A4       ();
	}
}