求问，CompositionReader是个python的第三方库吗？

13691255250

这个是cantera的一个事例，其中用到了import CompositionReader，求问CompositionReader怎么安装下载，我用pip失败了

1. # Clint Dunn, cdunn6754@gmail.com, 07-5-17
   
2. # written in Python 3.6.1 on cantera 2.3.0
   
3. 
   
4. ## Description
   
5. # A PSR combustor. Two separate streams - one coal volatiles and the other air, both at
   
6. # 300 K and 1 atm flow into an adiabatic and constant pressure
   
7. # combustor where they mix and burn.
   
8. 
   
9. # Writen so that we can write out a csv file for a particular combination of
   
10. # chemical mechanism and fuel composition. It writes these to the csv directory.
    
11. # The ppPsr.py script then reads these csv's and plots/ post processes them.
    
12. 
    
13. # From the original Cantera file starter thing that I kept:
    
14. # We are interested in the steady-state burning solution. Since at 300 K no
    
15. # reaction will occur between fuel and air, we need to use an 'igniter' to
    
16. # initiate the chemistry. A simple igniter is a pulsed flow of atomic hydrogen.
    
17. # After the igniter is turned off, the system approaches the steady burning
    
18. # solution.
    
19. 
    
20. 
    
21. # let me get at my utility functions
    
22. import sys
    
23. import os
    
24. sys.path.append(os.path.abspath("/home/cdunn3/Python_3.6.1/Cantera/utilities"))
    
25. 
    
26. import math
    
27. import csv
    
28. import matplotlib.pyplot as plt
    
29. import cantera as ct
    
30. import numpy as np
    
31. import CompositionReader as cr
    
32. 
    
33. #..........................................................................#
    
34. ## setup stuff
    
35. 
    
36. # user inputs
    
37. mechanism = 'drm'
    
38. fuel_name = 'heavy'
    
39. 
    
40. # Get the specified fuel compositions
    
41. fuel_dict = {'light':'light_fuel.txt',
    
42.              'heavy':'heavy_fuel.txt'
    
43.          }
    
44. composition = cr.read_compositions(fuel_dict[fuel_name])
    
45. 
    
46. # define the mechanism and gas
    
47. mech_dir = 'mechanisms/'
    
48. mechanism_dict = {'reduced':'r_creck_52.cti',
    
49.                   'full':'POLIMI_TOT_1412.cti',
    
50.                   'drm':'DRM_22_benzene.cti'
    
51.               }
    
52. gas = ct.Solution(mech_dir +  mechanism_dict[mechanism])
    
53. 
    
54. # ....................................................................#
    
55. ## run the simulation
    
56. 
    
57. # create a reservoir for the fuel inlet, and set to pure methane.
    
58. gas.TPX = 300.0, ct.one_atm, composition
    
59. fuel_in = ct.Reservoir(gas)
    
60. fuel_mw = gas.mean_molecular_weight
    
61. 
    
62. # use predefined function Air() for the air inlet
    
63. air = ct.Solution('air.xml')
    
64. air_in = ct.Reservoir(air)
    
65. air_mw = air.mean_molecular_weight
    
66. 
    
67. # to ignite the fuel/air mixture, we'll introduce a pulse of radicals. The
    
68. # steady-state behavior is independent of how we do this, so we'll just use a
    
69. # stream of pure atomic hydrogen.
    
70. gas.TPX = 300.0, ct.one_atm, 'H:1.0'
    
71. igniter = ct.Reservoir(gas)
    
72. 
    
73. # create the combustor, and fill it in initially with N2
    
74. gas.TPX = 300.0, ct.one_atm, 'N2:1.0'
    
75. combustor = ct.IdealGasReactor(gas)
    
76. combustor.volume = 1.0
    
77. 
    
78. # create a reservoir for the exhaust
    
79. exhaust = ct.Reservoir(gas)
    
80. 
    
81. # adjust the ratio coming in
    
82. Z  = 0.15
    
83. 
    
84. # compute fuel and air mass flow rates (light crashes at 2.5)
    
85. total_mdot = lambda t:  1 + (2.*t)**2.0
    
86. air_mdot = lambda t: total_mdot(t) * (1.0 -Z)
    
87. fuel_mdot = lambda t: total_mdot(t) * Z
    
88. 
    
89. # create and install the mass flow controllers. Controllers m1 and m2 provide
    
90. # constant mass flow rates, and m3 provides a short Gaussian pulse only to
    
91. # ignite the mixture
    
92. m1 = ct.MassFlowController(fuel_in, combustor, mdot=fuel_mdot)
    
93. 
    
94. # note that this connects two reactors with different reaction mechanisms and
    
95. # different numbers of species. Downstream and upstream species are matched by
    
96. # name.
    
97. m2 = ct.MassFlowController(air_in, combustor, mdot=air_mdot)
    
98. 
    
99. # The igniter will use a Gaussian time-dependent mass flow rate.
    
100. fwhm = 0.2
     
101. amplitude = 0.01
     
102. t0 = 0.5
     
103. igniter_mdot = lambda t: amplitude * math.exp(-(t-t0)**2 * 4 * math.log(2) / fwhm**2)
     
104. m3 = ct.MassFlowController(igniter, combustor, mdot=igniter_mdot)
     
105. 
     
106. # put a valve on the exhaust line to regulate the pressure
     
107. v = ct.Valve(combustor, exhaust, K=1.0)
     
108. 
     
109. # the simulation only contains one reactor
     
110. sim = ct.ReactorNet([combustor])
     
111. 
     
112. # take single steps to 6 s, writing the results to a CSV file for later
     
113. # plotting.
     
114. tfinal = 40.0
     
115. dt = 0.0001
     
116. times = np.arange(dt,tfinal,dt)
     
117. 
     
118. # data output arrays
     
119. states = ct.SolutionArray(gas, extra=['t', 'tres', 'fuel_mdot', 'Z'])
     
120. for i, time in enumerate(times):
     
121.     sim.advance(time)
     
122.     tres = combustor.mass/(m1.mdot(time) + m2.mdot(time))    #v.mdot(time)
     
123.     states.append(gas.state, t=time, tres=tres, fuel_mdot=m1.mdot(time), Z=Z)
     
124. 
     
125. csv_dir = 'csv_files/'
     
126. states.write_csv(csv_dir + fuel_name + '_' + mechanism + '_states.csv')
     

复制代码

Twilight6

应该不是，是自己写的代码吧