Related data repo.
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

166 lines
7.5KB

  1. units(length='cm', time='s', quantity='mol', act_energy='K')
  2. #ideal_gas(name='gas',
  3. # elements="N E Ar",
  4. # species="""e- N N+ N2 N2+ AR""",
  5. # reactions='all',
  6. # initial_state=state(temperature=300.0, pressure=OneAtm))
  7. #
  8. #ideal_gas(name='transport',
  9. # elements="N E Ar",
  10. # species="""N N+ N2 N2+ AR""",
  11. # transport='Mix',
  12. # initial_state=state(temperature=300.0, pressure=OneAtm))
  13. ideal_gas(name='gas',
  14. elements="N E",
  15. species="""e- N N+ N2 N2+ """,
  16. reactions='all',
  17. initial_state=state(temperature=300.0, pressure=OneAtm))
  18. ideal_gas(name='transport',
  19. elements="N E ",
  20. species="""N N+ N2 N2+ """,
  21. transport='Mix',
  22. initial_state=state(temperature=300.0, pressure=OneAtm))
  23. #-------------------------------------------------------------------------------
  24. # Species data
  25. #-------------------------------------------------------------------------------
  26. species(name='e-',
  27. atoms='E:1',
  28. thermo=(NASA9([298.15, 1000.00],
  29. [ 0.000000000E+00, 0.000000000E+00, 2.500000000E+00,
  30. 0.000000000E+00, 0.000000000E+00, 0.000000000E+00,
  31. 0.000000000E+00, -7.453750000E+02, -1.172081224E+01]),
  32. NASA9([1000.00, 6000.00],
  33. [ 0.000000000E+00, 0.000000000E+00, 2.500000000E+00,
  34. 0.000000000E+00, 0.000000000E+00, 0.000000000E+00,
  35. 0.000000000E+00, -7.453750000E+02, -1.172081224E+01]),
  36. NASA9([6000.00, 20000.00],
  37. [ 0.000000000E+00, 0.000000000E+00, 2.500000000E+00,
  38. 0.000000000E+00, 0.000000000E+00, 0.000000000E+00,
  39. 0.000000000E+00, -7.453750000E+02, -1.172081224E+01])),
  40. note='Ref-Species. Chase,1998 3/82. [g12/98]')
  41. #species(name = "AR",
  42. # atoms = " Ar:1 ",
  43. # thermo = (
  44. # NASA( [ 300.00, 1000.00], [ 2.500000000E+00, 0.000000000E+00,
  45. # 0.000000000E+00, 0.000000000E+00, 0.000000000E+00,
  46. # -7.453750000E+02, 4.366000000E+00] ),
  47. # NASA( [ 1000.00, 5000.00], [ 2.500000000E+00, 0.000000000E+00,
  48. # 0.000000000E+00, 0.000000000E+00, 0.000000000E+00,
  49. # -7.453750000E+02, 4.366000000E+00] )
  50. # ),
  51. # transport = gas_transport(
  52. # geom = "atom",
  53. # diam = 3.33,
  54. # well_depth = 136.50),
  55. # note = "120186"
  56. # )
  57. species(name='N',
  58. atoms='N:1',
  59. thermo=(NASA9([200.00, 1000.00],
  60. [ 0.000000000E+00, 0.000000000E+00, 2.500000000E+00,
  61. 0.000000000E+00, 0.000000000E+00, 0.000000000E+00,
  62. 0.000000000E+00, 5.610463780E+04, 4.193905036E+00]),
  63. NASA9([1000.00, 6000.00],
  64. [ 8.876501380E+04, -1.071231500E+02, 2.362188287E+00,
  65. 2.916720081E-04, -1.729515100E-07, 4.012657880E-11,
  66. -2.677227571E-15, 5.697351330E+04, 4.865231506E+00]),
  67. NASA9([6000.00, 20000.00],
  68. [ 5.475181050E+08, -3.107574980E+05, 6.916782740E+01,
  69. -6.847988130E-03, 3.827572400E-07, -1.098367709E-11,
  70. 1.277986024E-16, 2.550585618E+06, -5.848769753E+02])),
  71. transport = gas_transport(
  72. geom = "atom",
  73. diam = 3.30,
  74. well_depth = 71.40),
  75. note='Hf:Cox,1989. Moore,1975. Gordon,1999. [g 5/97]')
  76. species(name='N+',
  77. atoms='E:-1 N:1',
  78. thermo=(NASA9([298.15, 1000.00],
  79. [ 5.237079210E+03, 2.299958315E+00, 2.487488821E+00,
  80. 2.737490756E-05, -3.134447576E-08, 1.850111332E-11,
  81. -4.447350984E-15, 2.256284738E+05, 5.076830786E+00]),
  82. NASA9([1000.00, 6000.00],
  83. [ 2.904970374E+05, -8.557908610E+02, 3.477389290E+00,
  84. -5.288267190E-04, 1.352350307E-07, -1.389834122E-11,
  85. 5.046166279E-16, 2.310809984E+05, -1.994146545E+00]),
  86. NASA9([6000.00, 20000.00],
  87. [ 1.646092148E+07, -1.113165218E+04, 4.976986640E+00,
  88. -2.005393583E-04, 1.022481356E-08, -2.691430863E-13,
  89. 3.539931593E-18, 3.136284696E+05, -1.706646380E+01])),
  90. transport = gas_transport(
  91. geom = "atom",
  92. diam = 3.30,
  93. well_depth = 71.40),
  94. note='Moore,1975. Gordon,1999. [g 6/97]')
  95. species(name='N2',
  96. atoms='N:2',
  97. thermo=(NASA9([200.00, 1000.00],
  98. [ 2.210371497E+04, -3.818461820E+02, 6.082738360E+00,
  99. -8.530914410E-03, 1.384646189E-05, -9.625793620E-09,
  100. 2.519705809E-12, 7.108460860E+02, -1.076003744E+01]),
  101. NASA9([1000.00, 6000.00],
  102. [ 5.877124060E+05, -2.239249073E+03, 6.066949220E+00,
  103. -6.139685500E-04, 1.491806679E-07, -1.923105485E-11,
  104. 1.061954386E-15, 1.283210415E+04, -1.586640027E+01]),
  105. NASA9([6000.00, 20000.00],
  106. [ 8.310139160E+08, -6.420733540E+05, 2.020264635E+02,
  107. -3.065092046E-02, 2.486903333E-06, -9.705954110E-11,
  108. 1.437538881E-15, 4.938707040E+06, -1.672099740E+03])),
  109. transport = gas_transport(
  110. geom = "linear",
  111. diam = 3.62,
  112. well_depth = 97.53,
  113. polar = 1.76,
  114. rot_relax = 4.00),
  115. note='Ref-Elm. Gurvich,1978 pt1 p280 pt2 p207. [tpis78]')
  116. species(name='N2+',
  117. atoms='E:-1 N:2',
  118. thermo=(NASA9([298.15, 1000.00],
  119. [-3.474047470E+04, 2.696222703E+02, 3.164916370E+00,
  120. -2.132239781E-03, 6.730476400E-06, -5.637304970E-09,
  121. 1.621756000E-12, 1.790004424E+05, 6.832974166E+00]),
  122. NASA9([1000.00, 6000.00],
  123. [-2.845599002E+06, 7.058893030E+03, -2.884886385E+00,
  124. 3.068677059E-03, -4.361652310E-07, 2.102514545E-11,
  125. 5.411996470E-16, 1.340388483E+05, 5.090897022E+01]),
  126. NASA9([6000.00, 20000.00],
  127. [-3.712829770E+08, 3.139287234E+05, -9.603518050E+01,
  128. 1.571193286E-02, -1.175065525E-06, 4.144441230E-11,
  129. -5.621893090E-16, -2.217361867E+06, 8.436270947E+02])),
  130. transport = gas_transport(
  131. geom = "linear",
  132. diam = 3.62,
  133. well_depth = 97.53,
  134. polar = 1.76,
  135. rot_relax = 4.00),
  136. note='Gurvich,1989 pt1 p323 pt2 p200. [tpis89]')
  137. #-------------------------------------------------------------------------------
  138. # Reaction data
  139. #-------------------------------------------------------------------------------
  140. # Reaction 1
  141. reaction('N2 + N2 <=> N2 + N + N', [7.000000e+21, -1.6, 113200.0])
  142. # Reaction 2
  143. reaction('N2 + N <=> N + N + N', [3.000000e+22, -1.6, 113200.0])
  144. # Reaction 3
  145. reaction('N2 + e- <=> N + N + e-', [3.000000e+24, -1.6, 113200.0])
  146. # Reaction 4
  147. reaction('N + N <=> N2+ + e-', [4.400000e+07, 1.5, 67500.0])
  148. # Reaction 5
  149. reaction('N + e- <=> N+ + e- + e-', [2.500000e+34, -3.82, 168600.0])