Modeling study on the flow, heat transfer and energy conversion characteristics of low-power arc-heated hydrogen/nitrogen thrusters

A modeling study is conducted to investigate the effect of hydrogen content in propellants on the plasma flow, heat transfer and energy conversion characteristics of low-power (kW class) arc-heated hydrogen/nitrogen thrusters (arcjets). 1:0 (pure hydrogen), 3:1 (to simulate decomposed ammonia), 2:1 (to simulate decomposed hydrazine) and 0:1 (pure nitrogen) hydrogen/nitrogen mixtures are chosen as the propellants. Both the gas flow region inside the thruster nozzle and the anode-nozzle wall are included in the computational domain in order to better treat the conjugate heat transfer between the gas flow region and the solid wall region. The axial variations of the enthalpy flux, kinetic energy flux, directed kinetic-energy flux, and momentum flux, all normalized to the mass flow rate of the propellant, are used to investigate the energy conversion process inside the thruster nozzle. The modeling results show that the values of the arc voltage, the gas axial-velocity at the thruster exit, and the specific impulse of the arcjet thruster all increase with increasing hydrogen content in the propellant, but the gas temperature at the nitrogen thruster exit is significantly higher than that for other three propellants. The flow, heat transfer, and energy conversion processes taking place in the thruster nozzle have some common features for all the four propellants. The propellant is heated mainly in the near-cathode and constrictor region, accompanied with a rapid increase of the enthalpy flux, and after achieving its maximum value, the enthalpy flux decreases appreciably due to the conversion of gas internal energy into its kinetic energy in the divergent segment of the thruster nozzle. The kinetic energy flux, directed kinetic energy flux and momentum flux also increase at first due to the arc heating and the thermodynamic expansion, assume their maximum inside the nozzle and then decrease gradually as the propellant flows toward the thruster exit. It is found that a large energy loss (31-52%) occurs in the thruster nozzle due to the heat transfer to the nozzle wall and too long nozzle is not necessary. Modeling results for the NASA 1-kW class arcjet thruster with hydrogen or decomposed hydrazine as the propellant are found to compare favorably with available experimental data.

Numerical Simulation of a Low-power Hydrazine Arcjet Thruster

A modeling study is conducted to investigate the plasma flow and heat transfer characteristics of low-power (kW class) arc-heated thrusters (arcjets) with 2:1 hydrogen/nitrogen to simulate decomposed hydrazine as the propellant. The all-speed SIMPLE algorithm is employed to solve the governing equations, which take into account the effects of compressibility, the Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. Typical computed results about the temperature, velocity and Mach number distributions within arcjet thruster are presented for the case with arc current of 9 A and inlet stagnant pressure of 3.3×105 Pa to show the flow and heat transfer characteristics. It is found that the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appearing near the cathode tip, and the flow transition from the subsonic to supersonic regime occurs within the constrictor region. The effect of gas viscosity on the plasma flow within arcjet thruster is examined by an additional numerical test using artificially reduced values of gas viscosity. The test results show that the gas viscosity appreciably affects the plasma flow and the performance of the arcjet thruster for the cases with the hydrazine or hydrogen as the propellant. The integrated axial Lorentz force in the thruster nozzle is also calculated and compared with the thrust force of the arcjet thruster. It is found that the integrated axial Lorentz force is much smaller than the thrust force for the low-power arcjet thruster. Modeling results for the NASA 1-kW class arcjet thruster with simulated hydrazine as the propellant are found to be reasonably consistent with available experimental data.

川西亚高山森林生态系统辐射传输研究

辐射传输研究是贯穿森林生态系统的纽带，太阳辐射为植物的生长发育提供光合能量、适宜的环境温度以及发育信息。一方面，气候变化使到达地面辐射能的质和量发生变化，影响到植被的生长发育，改变森林的结构，而森林结构的变化又会影响林冠内辐射能的分配和质量，这些变化会进一步影响到林下土壤温度，改变森林根系活性以及土壤营养转化的效率；连锁反应的结果有可能会使森林生态系统的生产力发生变化，改变碳素和氮素源库的调节方向，从而反馈影响地球气候系统。另一方面，人类作为生态系统的成员，必然需要森林生态系统为其提供更多的原材料和更好的生态服务功能，如何实现这些目标，就需要人类适度调整干预方式和频度，达到预期的目的。本文在建立适合于川西亚高山森林的叶面积测量技术、光照辐射模型和土壤温度变化模型的基础上，对川西亚高山地带森林生态系统的辐射传输特征进行了分析，并从森林结构的角度探讨了林分内的辐射分布以及对土壤温度的影响。主要成果如下： 1. 提出了一种照相法测量叶面积的方法。通过对摆放在平面上的叶片照相，利用投影变化，把非正射图像转化为正射图像，然后经过计算机图像处理得到每一片叶片的面积、周长、长度、宽度等信息。这种方法可使用户以任意方向和距离拍摄处于平面上的叶片，能同时处理大量的叶片，适于野外离体或活体叶片测量。叶片面积分辨率可调，分辨率可以与常用的激光叶面积仪相近甚至更高，而且叶片图像可以存档查询。 2. 提出一种模拟林内光照变化的模型。利用林冠半球照片，记录视点以上半球内的林冠构件空间分布，作为林冠子模型；天空辐射子模型采用国际照明委员会(CIE)的标准晴天和阴天以及插值模型。该模型能够模拟林下某一位点处的实时光斑变化。 3. 提出一种土壤温度变化模型。把土壤视为具有容量和阻力性质的结构，利用电阻和电容器件构建土壤能量分布模型。外界太阳辐射能经过植被以及其它一些能量分配器后进入土壤，其中有一部分转化为土壤势能，即土壤温度。土壤温度的变化类似于电池的充放电过程。在已知模型参数的情况下，可以从太阳辐射计算土壤温度的变化。在模型参数未知的情况下，通过输入和输出值推算模型的参数，而模型参数中的时间常数与土壤组成和含水量有关，这样就可以知道土壤水分的变化情况。 4. 从王朗亚高山森林典型样地林分结构的测量获得林地三维结构图、树冠形态、叶面积密度等参数，这些参数输入到Brunner (1998)开发的tRAYci 模型中计算出一段时间内林分任意位置处的光照值。与林下辐射计测量值以及半球照片计算结果的比较，该模型基本上能够满足对林分光环境了解的要求。 5. 从川西亚高山森林生产力的角度，探讨了森林生产力研究的方法以及川西地区的研究历史和成果，发现了其中的一些规律和问题，特别是在叶面积测量上，还没有使用标准的叶面积指数定义。综合来看，川西地区针叶林叶面积指数(单位土地面积上植物冠层总叶面积的一半) 应在4-5 之间。降雨丰富的华西雨屏带是川西地区森林生产力最高的地区，而向西北森林生产力逐渐降低。川西地区云冷杉林森林生产力平均约为600 gDM m-2 a-1，但是根据辐射能计算的潜在生产力则达到1800 gDM m-2 a-1。实际与潜在森林生产力的巨大差异说明其它因子对生产力的影响。 6. 王朗亚高山3 个典型森林林分中，白桦林样地(BF) 林下草本以糙野青茅、牛至、紫菀等喜阳性物种为主，林下透光度较高；冷杉林样地(FF) 林下透光度最低，以喜阴性物种水金凤、蟹甲草、囊瓣芹等为主；而云杉林样地(SF)林分林龄最大，林下透光度介于冷杉林和白桦林之间，草本层仍然以喜阴性物种东方草莓、紫花碎米芥、酢浆草等为主。冷杉林和云杉林的灌木层也很丰富，卫矛属、五加属、茶藨子属、忍冬属植物很丰富，而在白桦林则以栒摘要子属、榛子属、鹅耳枥属等植物为主。藓类植物在云杉林中最丰富，并且形成毯状层，其它两个林分则很稀少。3 个样地林分结构与林下光环境有很强的相关性，从光环境特征可以在一定程度上推测林分的结构。各样地单纯从乔木层材积推算的NPP 排列顺序为BF>FF>SF，与林下辐射透射率和林分年龄的顺序相同，暗示辐射对群落演替过程的驱动作用。 7. 用半球照相法测得BF、FF 和SF 3 个样地的有效叶面积指数以SF 样地最高，BF 最低。如果考虑针叶树叶片在小枝上的丛聚分布，利用北方针叶林的数值进行校正，则SF 样地LAI 显著增加(达到89%)，其它样地的LAI 基本不变甚至有所下降。校正后的数值与文献中地面测量的结果较相近，说明在使用半球照相法测量川西亚高山针叶林LAI 时必须加以校正。 8. 在3 个样地中，白桦、岷江冷杉和方枝柏种群为丛聚分布，紫果云杉在FF和SF 样地中基本上为随机分布。3 个物种出现丛聚分布的最短距离约为2m，在最短距离以内则为随机分布。最短距离可能与树冠大小有关，种子传播特征以及对光照的需求状况可能是造成这种分布格局类型的原因。 Radiative transfer plays a key role in forest ecosystems. Solar radiation providesenergy for photosynthesis, appropriate ambient temperature and development informationfor plants. However, quality and quantity of radiation reaching land surface are affected byweather and subsequently influence the growth and development of plants, which in turnchanges the budget of radiation in forest. Soil temperature changes with the variation ofradiation under forest canopy and influences the activity of roots and rate of nutrientturnover. Thus, any changes of radiation will induce chain reactions in the entireecosystem and display in the value of net primary productivity which will possibly shiftthe relationship between carbon source and sink at local or regional scale and feed back tothe global climate system. On the other hand, as a component of ecosystems, humanbeings of course need to demand more materials and better service from ecosystems. Forthese purpose, man must adapt their pattern and frequency of interference to ecosystems.This paper aims to research on the canopy structure, the radiation distribution and theirinfluence on soil temperature from the process of radiative transfer in subalpine forestecosystem of western Sichuan. The main results are: 1 Present a new photogrammetric method for leaf area. The main idea is to convertnon-vertically taken images of planar leaves to orthoimages through projectivetransformation. The resultant images are used to get leaf morphological parametersthrough image processing. This method enables users to take photos at almost anyorientation and distance if only the leaves are placed on same plane, and to processlarge quantity of leaves in a short time, which is suitable for field measurement. Theresolution of leaf area is adjustable to fit for special requirement. 2 A model using hemispherical photos combining with solar tracks and radiation courseis provided to simulate light variation in forest. The hemispherical photos of canopyrecord the real spatial distribution of each element of plants viewed from a point. Skyradiance is simulated with CIE standard clear sky or cloudy sky model. This modelcan be used to simulate real time light variation under canopy. 3 Present a soil temperature model. Soil could be regarded as a body of resistor andcapacitor. Some of the budget of solar radiation in soil body is transformed into soilpotential energy, the soil temperature. Variation of soil temperature is driven by solarradiation, vegetation, soil properties, etc. This model has two parameters, one of whichis time constant and is related to soil water content. The inversed model can be used tosimulate the variation of soil water. 4 By using model tRAYci developed by Brunner (1998), the 3-D distribution of light inthree subalpine forest stands of Wanglang Nature Reserve has been simulated andvalidated with value of radiometers in these stands. This model can basically satisfythe need for understanding light regimes of these stands. 5 Present some principles and questions of NPP (net primary of productivity) researchesin western Sichuan. The standard leaf area index (LAI) defined by Chen and Black(1997) has not been used in this region. Total leaf area and projected leaf area indexare still used in NPP researches which may differ around 1-fold in magnitude. Thestandard LAI which is a half of total leaf area above unit land area should be between4 and 5 for typical subalpine coniferous forest of western Sichuan concluded fromliteratures. The maximum forest NPP occurs in West China rain belt and decreasesnorthwestwards. Average NPP of spruce-fir forest in western Sichuan is about600gDM m-2 a-1, which is below the potential NPP of 1800gDM m-2 a-1 based onmeasured radiation in this region. The significant difference between potential and realNPP suggests that other factors influence the growth of stands. 6 In the three subalpine forest stands of Wanglang Nature Reserve, herbage layer ofAbstractbirch stand (BF) with age of 40 is dominated by heliophytes of Deyeuxia scabrescens,Origanum vulgare, Aster tongoloa etc.. However, both of the other two stands aredominated by shade tolerent species, such as Impatiens noli-tangere, Impatiensdicentra, Cacalia deltophylla and Pternopetalum tanakae etc. in fir stand (FF) withage of 180 and Fragaria orientalis, Cardamine tangutorum and Oxalis corniculata etc.in spruce stand (SF) with age of 330. Shrub species in the latter two stands arerelatively rich, typical dominant genera being Euonymus, Acanthopanax, Ribes andLonicera. Birch stand has relatively sparse shrubs dominated by genera of Cotoneaster,Corylus and Carpinus. Mosses are significant only in spruce stand. The canopystructure controls the light regime of stand, which influence the composition of herblayers beneath the canopy. This light regime-community structure relationship can beused to infer the herb community from canopy structure. The NPP derived from timbervolume of arbor layer of the three stands decreases from BF to SF, which is in thesame order of transmitted total radiation under canopy and age of these stands,suggesting the driving effect of radiation in the succession of community. 7 The highest effective LAI of the three stands obtained by hemispherical photos is inplot SF and lowest in plot BF. After rectification of the clumping effect of leaves onshoot, the real LAI in plot SF increases significantly (89%) and approximate to theaverage LAI of coniferous forest in western Sichuan. Therefore, the LAI obtainedfrom hemispherical photos needs rectification for clumping effect. 8 Spatial distribution pattern for Betula platyphylla, Abies faxoniana and Sabinasaltuaria is clumpy, but Picea purpurea almost random in plot FF and SF. The shortestdistance for clumpy distribution for Betula platyphylla and Sabina saltuaria is 1.5m,and 2m for Abies faxoniana. And random pattern for these trees is exhibited within thisrange which almost coincides with the diameter of crown. Seed dispersalcharacteristics and light requirement may be the reason for different spatial pattern.

非粮原料燃料乙醇发酵关键技术的研究

生物质燃料乙醇是一种高度清洁的交通液体燃料，是减少温室气体排放，缓解大气污染的最佳技术选择。以非粮原料生产燃料乙醇可以在进行能源生产的同时保证粮食安全，有利于产业的可持续发展。在众多的非粮原料中，甘薯是我国开发潜力最大的生物质能源作物之一。我国占世界甘薯种植总面积和产量的90％。同时，甘薯的单位面积燃料乙醇产量远大于玉米和小麦。其成本是目前酒精中最低廉的，因此利用甘薯生产乙醇是发展生物质燃料乙醇的首要选择。目前采用薯类全原料主要采用分批发酵生产乙醇，其技术水平低，发酵强度低，一般在0.7－2.5g/（L•h），乙醇浓度低，甘薯发酵乙醇为6－8％(v/v)，能耗高，环境负荷大，污染严重。针对上述问题，本文从菌株选育、原料预处理、中试放大、残糖成分分析等方面进行研究。 为了研究乙醇发酵生产规模扩大过程中，大型发酵罐底部高压条件下，CO2对酵母乙醇发酵的影响，我们通过CO2 加压的方法进行模拟试验，研究结果表明，发酵时间随压强的升高而逐渐延长，高压CO2 对乙醇发酵效率影响不大，在0.3 MPa 以下时，发酵效率均可达到90%以上。高压CO2 对发酵的抑制作用是高压和CO2 这两个因素联合作用的结果。高压CO2 条件下，酵母胞外酶和胞内重要酶类的酶活均表现出特征性。0.2 MPa 下，酶活性的变化趋势和0.1 MPa 条件下的较为一致。而0.3 MPa 下的酶活变化趋势与0.4 MPa 下的酶活更为接近。通过全基因表达分析发现在CO2 压力为0.3 MPa 下，乙醇发酵途径中多个基因表达量下调，同时海藻糖合成酶和热激蛋白基因表达量上调。 筛选耐高温的乙醇酵母菌株能够解决糖化温度和发酵温度不协调的矛盾，实现真正意义上的边糖化边发酵。高温发酵还能够降低发酵时的冷却成本，实现乙醇的周年生产。本研究筛选出一株高温发酵菌株Y-H1，进而我们对该菌株的胞外酶和胞内乙醇代谢重要酶类的酶活性进行了分析。结果表明Y-H1 能够在40 ℃条件下正常进行乙醇发酵，发酵33h，最终乙醇浓度达到10.7%（w/w），发酵效率达到90%以上。同时发酵液最终pH 在3.5 左右，显示菌株具有一定的耐酸性能力。同时观察到40 ℃下，菌株的胞外酶和胞内乙醇代谢重要酶类的酶活性发生了变化，乙醇发酵途径中关键酶基因表达下调，而海藻糖合成酶与热激蛋白基因表达量上调，这些结果为进一步研究酵母菌耐热调控机理提供了依据。 糖蜜是一种大规模工业生产乙醇的理想原料，本研究利用选育高浓度乙醇发酵菌株结合配套的发酵稳定剂，研究了糖蜜高浓度乙醇发酵情况。结果表明采用冷酸沉淀预处理糖蜜溶液，采用分批补料的发酵方式，乙醇浓度最高达到了10.26% (w/w)，发酵时间为42 h。同时观察到在糖蜜发酵中，乙醛含量与乙醇浓度存在一定的相关性。 快速乙醇发酵对于缩短乙醇生产周期、降低乙醇生产成本、减少原料腐烂损失具有重要意义。本研究诱变和筛选得到了一株快速乙醇发酵菌株10232B。在优化后的发酵条件下，采用10L 发酵罐进行分批乙醇发酵，经过18h,乙醇的最终浓度达到88.5g/L,发酵效率93.6%，平均乙醇生产速度达到4.92 g/L/h。此菌株在保持较高乙醇生产浓度的同时，拥有快速生产乙醇的能力，适合作为快速乙醇发酵生产菌种。 由于鲜甘薯具有粘度大的特点，传统液化糖化处理很难在短时间内充分糖化原料；高粘度的醪液也难以进行管道输送，容易堵塞管路；同时，也会降低后续的乙醇发酵效率。 本文采用了快速粘度分析法对鲜甘薯糊化粘度特性进行了分析，进而对预处理条件进行了研究，在最佳预处理条件下，糖化2h 后，醪液葡萄糖值最高可达99.3，粘度4.5×104 mPa.s，而采用传统糖化工艺，醪液DE 值仅为85.8，粘度大于1.0×105 mPa.s。 此预处理方法也可用于快速糖化不加水的醪液。后续的乙醇发酵试验表明，通过此预处理方法获得的糖化醪液对乙醇发酵无负面影响。 在前期已实现了实验室水平的鲜甘薯燃料乙醇快速乙醇发酵基础上，进一步将发酵规模扩大到500L，在中试水平上对甘薯乙醇发酵进行了研究。结果表明在500L 中试规模，采用边糖化边发酵（SSF）工艺，在料液比为3∶1，发酵醪液最高粘度为6×104mPa.s 条件下，发酵37h，乙醇浓度达到了12.7%（v/v），发酵效率91%，发酵强度为2.7 g/（L•h）。与目前国内的薯类乙醇发酵生产技术水平具有明显的优越性。 为研究甘薯、木薯乙醇发酵中残糖的组成，采用了高效液相色谱—蒸发光散射检测法，对乙醇发酵残糖进行了分析。结果表明，甘薯、木薯乙醇发酵残糖均为寡聚糖，主要由葡萄糖、木糖、半乳糖、阿拉伯糖和甘露糖构成。随着发酵时间延长，寡聚糖中的葡萄糖、半乳糖、甘露糖可被缓慢的水解释放。提高糖化酶量仅在一定程度上降低残糖，过量的糖化酶反而会导致残糖增加。同时发现3, 5-二硝基水杨酸法不能准确测定甘薯、木薯乙醇发酵中的残总糖含量。进一步筛选了两株残糖降解菌株，对甘薯乙醇发酵残糖的降解利用率均达到了40%以上，而且还能显著降低发酵醪液粘度。经形态学和rRNA ITS 序列分析，确定这两株菌分别属于为木霉属和曲霉属黑曲霉组。 通过对以甘薯原料为代表的非粮原料发酵技术研究开发，以期形成乙醇转化率高，能耗低，生产效率高、季节适应性好，原料适应性广，经济性强，符合清洁生产机制的燃料乙醇高效转化技术，为具有我国特色的燃料乙醇发展模式提供技术支持。 Sweet potato is one of the major feedstock for the fuel ethanol production in China. The planting area and the yield in China take 90% of the world. Sweet potato is an efficient kind of energy crops. The energy outcome per area is higher than corn or wheat. And the manufacture cost of ethanol is the lowest, compared with corn and wheat. So sweet potato is the favorable crop for the bioethanol production in China. However, the low-level fermentation technology restricts the development of ethanol production by sweet potato, including slow ethanol production rate, low ethanol concentration and high energy cost. To solve these problems, we conducted research on the strain breeding, pretreatment, pilot fermentation test and residual saccharides analysis. To study the impact of hyperbaric condition at bottom of the large fermentor on yeast fermentation, high pressure carbon dioxide (CO2) was adopted to simulate the situation. The results showed that the fermentation was prolonged with the increasing pressure. The pressure of CO2 had little impact on the ethanol yield which could reach 90% under the pressure below 0.3 MPa. The inhibition was combined by the high pressure and CO2. Under the high CO2 pressure, the extracellular and important intracellular enzyme activities were different from those under normal state. The changes under 0.1 MPa and 0.2 MPa were similar. The changes under 0.3 MPa were closer to those under 0.4 MPa. The application of thermotolerance yeast could solve the problem of the inconsistent temperature between fermentation and saccharificaton and fulfill the real simultaneous saccharification and fermentation. And it could reduce the cooling cost. A thermotolerance strain Y-H1 was isolated in our research. It gave high ethanol concentration of 10.7%（w/w）at 40 ℃ for 33 h. The ethanol yield efficiency was over 90%. At 40 ℃, the extracellular and important intracellular enzyme activities of Y-H1 showed the difference with normal state, which may indicate its physiological changes at the high temperature. Molasses is another feedstock for industrial ethanol production. By our ethanol-tolerance strain and the regulation reagents, the fermentation with high ethanol concentration was investigated. In fed-batch mode combined with cold acid deposition, the highest ethanol concentration was 10.26% (w/w) for 42h. The aldehyde concentration in fermentation was found to be related to ethanol concentration. The development of a rapid ethanol fermentation strain of Zymomonas mobilis is essential for reducing the cost of ethanol production and for the timely utilization of fresh material that is easily decayed in the Chinese bioethanol industry. A mutant Z. mobilis strain, 10232B, was generated by UV mutagenesis. Under these optimized conditions, fermentation of the mutant Z. mobilis 10232B strain was completed in just 18 h with a high ethanol production rate, at an average of 4.92 gL-1h-1 per batch. The final maximum ethanol concentration was 88.5 gL-1, with an ethanol yield efficiency of 93.6%. This result illustrated the potential use of the mutant Z. mobilis 10232B strain in rapid ethanol fermentation in order to help reduce the cost of industrial ethanol production. As fresh sweet potato syrup shows high viscosity, it is hard to be fully converted to glucose by enzymes in the traditional saccharification process. The high-viscosity syrup is difficult to be transmitted in pipes, which may be easily blocked. Meanwhile it could also reduce the later ethanol fermentation efficiency. To solve these problems, effects of the pretreatment conditions were investigated. The highest dextrose equivalent value of 99.3 and the lowest viscosity of 4.5×104 mPa.s were obtained by the most favorable pretreatment conditions, while those of 85.8 and over 1.0×105 mPa.s was produced by traditional treatment conditions. The pretreatment could also be applied on the material syrup without adding water. The later experiments showed that the pretreated syrup had no negative effect on the ethanol fermentation and exhibited lower viscosity. The fuel ethanol rapid production from fresh sweet potato was enlarged in the 500L pilot scale after its fulfillment on the laboratory level. The optimal ratio of material to water was 3 to 1 in 500L fermentor. With low-temperature-cooking (85 ℃) using SSF, the Saccharomyces cerevisiae was able to produce ethanol 97.44 g/kg for 37h, which reached 92% of theoretical yield. The average ethanol production rate was 4.06 g/kg/h. And the maximum viscosity of syrup reached 6×104mPa.s. The results showed its superiority over current industrial ethanol fermentation. The compositions of the residual saccharides in the ethanol fermentation by sweet potato and cassava were analyzed by high performance liquid chromatography coupled with evaporative light-scattering detector. The results showed that all the residual saccharides were oligosaccharides, mainly composed of glucose, xylose, galactose, arabinose and mannose. The glucose, galactose and mannose could be slowly hydrolyzed from oligosaccharides in syrup during a long period. To increase the glucoamylase dosage could lower the residual saccharides to a certain extent. However, excess glucoamylase dosage led to more residual saccharides. And the method of 3, 5-dinitrosalicylic acid could not accurately quantify the residual total saccharides content. Two residual saccharides degrading strains were isolated, which could utilize 40% of total residual saccharide and lower the syrup viscosity. With the analysis of morphology and internal transcribed spacer sequence, they were finally identified as species of Trichoderma and Aspergillus niger.

粘红酵母对铀的吸附研究

Characteristic of uranium biosorption in water solution by Rhodotorula glutinis was investigated in the present study and the optimal pH for uranium adsorption was found to be 6-7.At the same time,maximum adsorption capacity of 149.4 mgU/(g dry cell)was identified,and Langmuir adsorption models can be used to simulate the isothermal biosorption process with high correlation coefficient of 0.99.According to Fourier transform infrared spectra,a new peak at wave number of 904 cm-1,which can be assigned to the ...中文摘要：研究了粘红酵母对水溶液中铀的吸附行为,发现其吸附铀的最佳pH值为6～7,最大吸附量为149.4mgU.g-1,其吸附等温线和Langmuir吸附等温方程符合较好,相关系数R2达到0.99;比较吸附铀前后粘红酵母的红外光谱图发现,吸附过铀的菌体的红外光谱在904cm-1处出现了一个新的峰,此峰为UO2的伸缩振动峰,说明粘红酵母确实对铀发生了吸附作用。此外,氨基或羟基的伸缩振动峰由3309移至3287cm-1,细胞壁中碳水化合物或醇中C—O键伸缩振动发生位移,由1068移至1080cm-1,说明这些基团可能参与了吸附过程;蛋白质的特征吸收峰(1653,1540,1237cm-1)在吸附前后基本无明显变化,表明粘红酵母的主要成分及结构仍保持完整。吸附后的菌体利用0.1mol.L-1的NaHCO3处理后可解吸出其中96%的铀,可见该菌在铀矿废水处理方面具有广阔的应用前景。

气相色谱过程的蒙特卡罗模拟

Based on Monte Carlo method the microscopic gas-phase adsorption chromatographic behaviors of volatile chemical compounds were simulated,and the migration process of those compounds followed by the carrier gases along the chromatographic surface and the process of adsorption-desorption were also described.The programs have been compiled to simulate the processes of thermochromatography and isothermal chromatography with different experimental conditions and properties.The calculated results agree well with ...中文摘要：用蒙特卡罗方法模拟了易挥发化合物的气相色谱微观动力学行为,描述了化合物随气流在固体色谱柱表面的迁移以及吸附-解吸等过程。按气相色谱微观动力学模型编译的计算机程序模拟了超重元素化合物在等温色谱及热色谱的动力学过程,并根据不同实验条件进行了大量计算。计算结果与实验数据较为吻合。讨论了超重元素的半衰期、吸附态的周期、载气流量以及化合物的质量密度等对实验结果的影响及该理论模型的优点和待改进的地方。

Measurement of neutron flow with the external target experimental device on CSR

With the construction of the neutron detection wall at the external target position on Heavy Ion Research Facility in Lanzhou-Cooling Storage Ring (HIRFL-CSR), it will be possible to detect high energy neutron. A BUU model is applied to simulate the flow in both symmetric (Ni+Ni, Pb+Pb) and asymmetric(Pb+Ni) systems. It is shown that at above several hundreds MeV/u, the flow signals are very obvious and depend clearly on the centrality of the collisions. Based on the products in the forward angle less than 20 degrees, the simulation also reveals that the determination of the reaction plane and the selection of the impact parameter, both of which are essential in the flow measurement, are well implemented. The double event and its influence on the determination of the neutron flow are also simulated.

Simulations to design an online motion compensation system for scanned particle beams

Respiration-induced target motion is a major problem in intensity-modulated radiation therapy. Beam segments are delivered serially to form the total dose distribution. In the presence of motion, the spatial relation between dose deposition from different segments will be lost. Usually, this results in over-and underdosage. Besides such interplay effects between target motion and dynamic beam delivery as known from photon therapy, changes in internal density have an impact on delivered dose for intensity-modulated charged particle therapy. In this study, we have analysed interplay effects between raster scanned carbon ion beams and target motion. Furthermore, the potential of an online motion strategy was assessed in several simulations. An extended version of the clinical treatment planning software was used to calculate dose distributions to moving targets with and without motion compensation. For motion compensation, each individual ion pencil beam tracked the planned target position in the lateral aswell as longitudinal direction. Target translations and rotations, including changes in internal density, were simulated. Target motion simulating breathing resulted in severe degradation of delivered dose distributions. For example, for motion amplitudes of +/- 15 mm, only 47% of the target volume received 80% of the planned dose. Unpredictability of resulting dose distributions was demonstrated by varying motion parameters. On the other hand, motion compensation allowed for dose distributions for moving targets comparable to those for static targets. Even limited compensation precision (standard deviation similar to 2 mm), introduced to simulate possible limitations of real-time target tracking, resulted in less than 3% loss in dose homogeneity.

Methanol synthesis from CO2 using a silicone rubber/ceramic composite membrane reactor

A one-dimensional isothermal pseudo-homogeneous parallel flow model was developed for the methanol synthesis from CO2 in a silicone rubber/ceramic composite membrane reactor. The fourth-order Runge-Kutta method was adopted to simulate the process behaviors in the membrane reactor. How those parameters affect the reaction behaviors in the membrane reactor, such as Damkohler number Da, pressure ratio p(r), reaction temperature T, membrane separation factor alpha, membrane permeation parameter phi , as well as the non-uniform parameter of membrane permeation L-1, were discussed in detail. Parts of the theoretical results were tested and verified; the experimental results showed that the conversion of the main reaction in the membrane reactor increased by 22% against traditional fixed bed reactor, and the optimal non-uniform parameter of membrane permeation rate, L-1.opt ,does exist. (C) 2003 Elsevier B.V All rights reserved.

Geometries of the Halocarbene anions HCF- and CF2-: ab initio calculation and Franck-Condon analysis

A theoretical method to calculate multidimensional Franck-Condon factors including Duschinsky effects is described and used to simulate the photoelectron spectra of HCF- and CF2- radicals. Geometry optimization and harmonic vibrational frequency calculations have been performed on the (X) over tilde (1)A' state of HCF and (X) over tilde (2)A" state of HCF-, and (X) over tilde (1)A(1) state of CF2 and (X) over tilde B-2(1) state of CF2-. Franck-Condon analyses and spectral simulation were carried out on the first photoelectron band of HCF- and CF2- respectively. The theoretical spectra obtained by employing B3LYP/6-311 + G(2d,p) values are in excellent agreement with the observed ones. In addition, the equilibrium geometry parameters, R(CF) = 0.1475 +/- 0.0005 nm, of the (X) over tilde (2)A" state of HCF-, and r(FC) = 0.1425 +/- 0.0005 nm and angle(FCF) = 100.5 +/- 0.5degrees, of the (X) over tilde B-2(i) state of CF2-, are derived by employing an iterative Franck-Condon analysis procedure in the spectral simulation. (C) 2003 Elsevier B.V. All rights reserved.

不同森林管理预案下友好林业局森林景观的动态变化

Knowledge about cumulative effects of forest management alternatives on forest landscape is required to make forest management decision. In this paper,a spatially explicit landscape model,LANDIS,was applied to simulate forest landscape changes in 200 years under four management alternatives (no cutting,clearcutting,selective cutting I and II) in Youhao Forestry Bureau located in Small Khingan Mountains. APACK was used to calculate distribution area of the representative species and species age cohort for six species. The results showed:1) timber harvest decreased area percentage of representative conifer species,Pinus koraiensis,Picea koraiensis and Picea jezoensis,Larix gmelinii to some extent compared to no cutting. The most influencing cutting mode for the area percentage of Tilia amurensis and Quercus mongolica was selective cutting II,followed by selective cutting I and clearcutting. To the contrast,the change of area percentage of Betula phatyphylla was contrary to the management alternatives; 2) As to species age cohort composition,timber harvest significantly changed age structure,that is,it decreased over-matured age cohort of representative species,and increased seedling and middle-age cohort (B. phatyphylla was not included,because its area percentage of over-mature age cohort was the highest under clearcutting than other three scenarios).

An alternative approach to ascertain the rate-determining steps of TiO2 photoelectrocatalytic reaction by electrochemical impedance spectroscopy

In an attempt to ascertain the rate-determining steps (RDS) of TiO2 photoelectrocatalytic (PEC) reaction, the PEC oxidation of sulfosalicylic acid (SSA) solution in a TiO2-coated electrode reactor system was monitored by applying the electrochemical impedance spectroscopy (EIS) method. In the meantime, an EIS mathematical model was first established to theoretically simulate the PEC reaction. Based on the EIS model, the theoretical simulation indicates three typical reactions in a PEC oxidation process, which include the charge-transfer-dominated reaction, both the charge-transfer- and adsorption-dominated reaction, and the adsorption-dominated reaction. The experimental results of EIS measurement showed that there was only one arc/semicircle on the EIS plane display when the external bias applied was below 200 mV (vs SCE) in the SSA PEC degradation whereas there were two arcs/semicircles when the externally applied bias exceeded 200 mV (vs SCE). The experimental results have a good agreement with the model simulation. The EIS method in this study provides an easier way to determine the RDS in a PEC oxidation process, which would be helpful to better control the reaction in practice.

Transient state of chiral recognition in a binary mixture of cyclodextrins in capillary electrophoresis

The transient state (as the defined point where no enantioseparation is obtained in a dual chiral selector system) of chiral recognition of aminoglutethimide in a binary mixture of neutral cyclodextrins (CDs) was studied by capillary electrophoresis (CE). The following three dual selector systems were used: alpha-cyclodextrin (alpha-CD) and beta-cyclodextrin (beta-CD); alpha-CD and heptakis(di-O-methyl-beta-cyclodextrin) (DM-beta-CD); alpha-CD and heptakis(tri-O-methyl-beta-cyclodextrin) (TM-beta-CD). The S-(-) enantiomer of the analyte was more strongly retained in the presence of either alpha-CD or TM-beta-CD at pH 2.5, 100 mM phosphate buffer, while the R-(+) enantiomer was more strongly retained in the presence of either P-CD or DM-P-CD. In the more simple case, the elution order is invariably kept if the enantiomers have the same elution order in either one of the two hosts of the binary mixture. In contrast, the elution order may be switched by varying the concentration ratio of two hosts that produce opposite elution order for this particular analyte. In such a dual selector system, the enantioselectivity will disappear at the transient state at a certain ratio of host,:host, Moreover, the migration times of the two enantiomers with host, alone (diluted in buffer) is approximately equal to the migration times at the corresponding concentration of host, alone (diluted in buffer), where the ratio of concentrations of host,:host, is the same as in the binary mixture at the transient state. As found by nuclear magnetic resonance experiments, the analyte is forming a 1:1 complex with either one of the CDs applied. From this finding, a theoretical model based on the mobility difference of the two enantiomers was derived that was used to simulate the transient state. (C) 2000 Elsevier Science B.V. All rights reserved.

Stratigraphic Model Predictions of Geoacoustic Properties

Geoacoustic properties of the seabed have a controlling role in the propagation and reverberation of sound in shallow-water environments. Several techniques are available to quantify the important properties but are usually unable to adequately sample the region of interest. In this paper, we explore the potential for obtaining geotechnical properties from a process-based stratigraphic model. Grain-size predictions from the stratigraphic model are combined with two acoustic models to estimate sound speed with distance across the New Jersey continental shelf and with depth below the seabed. Model predictions are compared to two independent sets of data: 1) Surficial sound speeds obtained through direct measurement using in situ compressional wave probes, and 2) sound speed as a function of depth obtained through inversion of seabed reflection measurements. In water depths less than 100 m, the model predictions produce a trend of decreasing grain-size and sound speed with increasing water depth as similarly observed in the measured surficial data. In water depths between 100 and 130 m, the model predictions exhibit an increase in sound speed that was not observed in the measured surficial data. A closer comparison indicates that the grain-sizes predicted for the surficial sediments are generally too small producing sound speeds that are too slow. The predicted sound speeds also tend to be too slow for sediments 0.5-20 m below the seabed in water depths greater than 100 m. However, in water depths less than 100 m, the sound speeds between 0.5-20-m subbottom depth are generally too fast. There are several reasons for the discrepancies including the stratigraphic model was limited to two dimensions, the model was unable to simulate biologic processes responsible for the high sound-speed shell material common in the model area, and incomplete geological records necessary to accurately predict grain-size

Fast monostatic GPR modeling

We propose the exploding-reflector method to simulate a monostatic survey with a single simulation. The exploding reflector, used in seismic modeling, is adapted for ground-penetrating radar (GPR) modeling by using the analogy between acoustic and electromagnetic waves. The method can be used with ray tracing to obtain the location of the interfaces and estimate the properties of the medium on the basis of the traveltimes and reflection amplitudes. In particular, these can provide a better estimation of the conductivity and geometrical details. The modeling methodology is complemented with the use of the plane-wave method. The technique is illustrated with GPR data from an excavated tomb of the nineteenth century.