Steady state optimization of an ammonia reactor

A hybrid simulation technique for identification and steady state optimization of a tubular reactor used in ammonia synthesis is presented. The parameter identification program finds the catalyst activity factor and certain heat transfer coefficients that minimize the sum of squares of deviation from simulated and actual temperature measurements obtained from an operating plant. The optimization program finds the values of three flows to the reactor to maximize the ammonia yield using the estimated parameter values. Powell's direct method of optimization is used in both cases. The results obtained here are compared with the plant data.

Evaluation of residence time on nitrogen oxides removal in non-thermal plasma reactor

Non-thermal plasma (NTP) has been introduced over the last few years as a promising after- treatment system for nitrogen oxides and particulate matter removal from diesel exhaust. NTP technology has not been commercialised as yet, due to its high rate of energy consumption. Therefore, it is important to seek out new methods to improve NTP performance. Residence time is a crucial parameter in engine exhaust emissions treatment. In this paper, different electrode shapes are analysed and the corresponding residence time and NOx removal efficiency are studied. An axisymmetric laminar model is used for obtaining residence time distribution numerically using FLUENT software. If the mean residence time in a NTP plasma reactor increases, there will be a corresponding increase in the reaction time and consequently the pollutant removal efficiency increases. Three different screw thread electrodes and a rod electrode are examined. The results show the advantage of screw thread electrodes in comparison with the rod electrode. Furthermore, between the screw thread electrodes, the electrode with the thread width of 1 mm has the highest NOx removal due to higher residence time and a greater number of micro-discharges. The results show that the residence time of the screw thread electrode with a thread width of 1 mm is 21% more than for the rod electrode.

Kinetics of Extraction of Iron(II1) by Bis(2-ethylhexyl Phosphate in a Laminar Liquid-Liquid Jet Reactor

The kinetics of iron(II1) extraction by bis(Zethylhexy1) phosphate (HDEHP, HA) in kerosene from sulfuric acid solutions has been studied in a liquid-liquid laminar jet reactor. The contact time of the interface in this reacting device is of the same order of magnitude as the surface renewal time in dispersion mixing and much less than that obtained in the relatively quiescent condition of the Lewis cell. Yet the analysis of the data in this study suggested a rate-controlling step involving surface saturation quite in conformity with that obtained in the Lewis cell and not with that in dispersion mixing as reported in the literature. Further, the mechanism suggested a weaker dependence of the rate on hydrogen ion concentration which was reported by other workers.

Interatomic potentials for fusion reactor material simulations

Fusion energy is a clean and safe solution for the intricate question of how to produce non-polluting and sustainable energy for the constantly growing population. The fusion process does not result in any harmful waste or green-house gases, since small amounts of helium is the only bi-product that is produced when using the hydrogen isotopes deuterium and tritium as fuel. Moreover, deuterium is abundant in seawater and tritium can be bred from lithium, a common metal in the Earth's crust, rendering the fuel reservoirs practically bottomless. Due to its enormous mass, the Sun has been able to utilize fusion as its main energy source ever since it was born. But here on Earth, we must find other means to achieve the same. Inertial fusion involving powerful lasers and thermonuclear fusion employing extreme temperatures are examples of successful methods. However, these have yet to produce more energy than they consume. In thermonuclear fusion, the fuel is held inside a tokamak, which is a doughnut-shaped chamber with strong magnets wrapped around it. Once the fuel is heated up, it is controlled with the help of these magnets, since the required temperatures (over 100 million degrees C) will separate the electrons from the nuclei, forming a plasma. Once the fusion reactions occur, excess binding energy is released as energetic neutrons, which are absorbed in water in order to produce steam that runs turbines. Keeping the power losses from the plasma low, thus allowing for a high number of reactions, is a challenge. Another challenge is related to the reactor materials, since the confinement of the plasma particles is not perfect, resulting in particle bombardment of the reactor walls and structures. Material erosion and activation as well as plasma contamination are expected. Adding to this, the high energy neutrons will cause radiation damage in the materials, causing, for instance, swelling and embrittlement. In this thesis, the behaviour of a material situated in a fusion reactor was studied using molecular dynamics simulations. Simulations of processes in the next generation fusion reactor ITER include the reactor materials beryllium, carbon and tungsten as well as the plasma hydrogen isotopes. This means that interaction models, {\it i.e. interatomic potentials}, for this complicated quaternary system are needed. The task of finding such potentials is nonetheless nearly at its end, since models for the beryllium-carbon-hydrogen interactions were constructed in this thesis and as a continuation of that work, a beryllium-tungsten model is under development. These potentials are combinable with the earlier tungsten-carbon-hydrogen ones. The potentials were used to explain the chemical sputtering of beryllium due to deuterium plasma exposure. During experiments, a large fraction of the sputtered beryllium atoms were observed to be released as BeD molecules, and the simulations identified the swift chemical sputtering mechanism, previously not believed to be important in metals, as the underlying mechanism. Radiation damage in the reactor structural materials vanadium, iron and iron chromium, as well as in the wall material tungsten and the mixed alloy tungsten carbide, was also studied in this thesis. Interatomic potentials for vanadium, tungsten and iron were modified to be better suited for simulating collision cascades that are formed during particle irradiation, and the potential features affecting the resulting primary damage were identified. Including the often neglected electronic effects in the simulations was also shown to have an impact on the damage. With proper tuning of the electron-phonon interaction strength, experimentally measured quantities related to ion-beam mixing in iron could be reproduced. The damage in tungsten carbide alloys showed elemental asymmetry, as the major part of the damage consisted of carbon defects. On the other hand, modelling the damage in the iron chromium alloy, essentially representing steel, showed that small additions of chromium do not noticeably affect the primary damage in iron. Since a complete assessment of the response of a material in a future full-scale fusion reactor is not achievable using only experimental techniques, molecular dynamics simulations are of vital help. This thesis has not only provided insight into complicated reactor processes and improved current methods, but also offered tools for further simulations. It is therefore an important step towards making fusion energy more than a future goal.

Computer simulation of H and He effects in fusion reactor materials

Fusion power is an appealing source of clean and abundant energy. The radiation resistance of reactor materials is one of the greatest obstacles on the path towards commercial fusion power. These materials are subject to a harsh radiation environment, and cannot fail mechanically or contaminate the fusion plasma. Moreover, for a power plant to be economically viable, the reactor materials must withstand long operation times, with little maintenance. The fusion reactor materials will contain hydrogen and helium, due to deposition from the plasma and nuclear reactions because of energetic neutron irradiation. The first wall divertor materials, carbon and tungsten in existing and planned test reactors, will be subject to intense bombardment of low energy deuterium and helium, which erodes and modifies the surface. All reactor materials, including the structural steel, will suffer irradiation of high energy neutrons, causing displacement cascade damage. Molecular dynamics simulation is a valuable tool for studying irradiation phenomena, such as surface bombardment and the onset of primary damage due to displacement cascades. The governing mechanisms are on the atomic level, and hence not easily studied experimentally. In order to model materials, interatomic potentials are needed to describe the interaction between the atoms. In this thesis, new interatomic potentials were developed for the tungsten-carbon-hydrogen system and for iron-helium and chromium-helium. Thus, the study of previously inaccessible systems was made possible, in particular the effect of H and He on radiation damage. The potentials were based on experimental and ab initio data from the literature, as well as density-functional theory calculations performed in this work. As a model for ferritic steel, iron-chromium with 10% Cr was studied. The difference between Fe and FeCr was shown to be negligible for threshold displacement energies. The properties of small He and He-vacancy clusters in Fe and FeCr were also investigated. The clusters were found to be more mobile and dissociate more rapidly than previously assumed, and the effect of Cr was small. The primary damage formed by displacement cascades was found to be heavily influenced by the presence of He, both in FeCr and W. Many important issues with fusion reactor materials remain poorly understood, and will require a huge effort by the international community. The development of potential models for new materials and the simulations performed in this thesis reveal many interesting features, but also serve as a platform for further studies.

Oxidation of acetaldehyde to acetic acid in a sparger reactor

The kinetics of oxidation of acetaldehyde to acetic acid was studied in a sparger reactor using manganese acetate as the catalyst. Data obtained in a stirred tank reactor are used for analyzing the sparger reactor data. The rate of chemical reaction is extremely fast and can be neglected for the rate equation of the sparger reactor. A kinetic model applicable at any temperature and concentration within the range of the variables studied is developed which predicts the performance of the sparger reactor satisfactorily.

Acoustic emission technique for leak detection in an end shield of a pressurised heavy water reactor

This paper discusses a successful application of the Acoustic Emission Technique (AET) for the detection and location of leak paths present on an inaccessible side of an end shield of a Pressurised Heavy Water Reactor (PHWR). The methodology was based on the fact that air- and water-leak AE signals have different characteristic features. Baseline data was generated from a sound end shield of a PHWR for characterising the background noise. A mock-up end shield system with saw-cut leak paths was used to verify the validity of the methodology. It was found that air-leak signals under pressurisation (as low as 3 psi) could be detected by frequency domain analysis. Signals due to air leaks from various locations of defective end shield were acquired and analysed. It was possible to detect and locate leak paths. The presence of detected leak paths was further confirmed by an alternative test.

Optimization of an ammonia reactor using regression analysis

This paper presents an optimization algorithm for an ammonia reactor based on a regression model relating the yield to several parameters, control inputs and disturbances. This model is derived from the data generated by hybrid simulation of the steady-state equations describing the reactor behaviour. The simplicity of the optimization program along with its ability to take into account constraints on flow variables make it best suited in supervisory control applications.

Computer simulation of multi-elemental fusion reactor materials

Thermonuclear fusion is a sustainable energy solution, in which energy is produced using similar processes as in the sun. In this technology hydrogen isotopes are fused to gain energy and consequently to produce electricity. In a fusion reactor hydrogen isotopes are confined by magnetic fields as ionized gas, the plasma. Since the core plasma is millions of degrees hot, there are special needs for the plasma-facing materials. Moreover, in the plasma the fusion of hydrogen isotopes leads to the production of high energetic neutrons which sets demanding abilities for the structural materials of the reactor. This thesis investigates the irradiation response of materials to be used in future fusion reactors. Interactions of the plasma with the reactor wall leads to the removal of surface atoms, migration of them, and formation of co-deposited layers such as tungsten carbide. Sputtering of tungsten carbide and deuterium trapping in tungsten carbide was investigated in this thesis. As the second topic the primary interaction of the neutrons in the structural material steel was examined. As model materials for steel iron chromium and iron nickel were used. This study was performed theoretically by the means of computer simulations on the atomic level. In contrast to previous studies in the field, in which simulations were limited to pure elements, in this work more complex materials were used, i.e. they were multi-elemental including two or more atom species. The results of this thesis are in the microscale. One of the results is a catalogue of atom species, which were removed from tungsten carbide by the plasma. Another result is e.g. the atomic distributions of defects in iron chromium caused by the energetic neutrons. These microscopic results are used in data bases for multiscale modelling of fusion reactor materials, which has the aim to explain the macroscopic degradation in the materials. This thesis is therefore a relevant contribution to investigate the connection of microscopic and macroscopic radiation effects, which is one objective in fusion reactor materials research.

General model of a foam bed reactor

A general model of a foam bed reactor has been developed which rigorously accounts for the extent of gas absorption with chemical reaction occurring in both the storage and foam sections. Its applicability extends to a wide spectrum of reaction velocities. The possibilities of the predominance of the bulk-liquid reaction in the storage section or the absorption with reaction in the foam section can be handled as merely special cases of the general analysis. The importance of foam for carrying out a particular gas-liquid reaction is characterised by a criterion in terms of the fractional rate of reaction in the foam section. Trends of variations in the concentrations of dissolved free A, solute B, and gas-phase A with time of operation of the reactor are presented. The nature of the variation in the fractional rate of reaction in the foam section with time, at different reaction velocities, and the effect of the liquid flow rate (across the storage section) on the transience are also illustrated. Finally, the predictions of the general model have been validated using the available experimental data on the oxidation of sodium sulphide in a foam bed reactor. The agreement between the experimental and the present theoretical information is fairly good, apart from being more insightful than all the previous models of this reactor.

Modelling of a batch sonochemical reactor

Ultrasonication of aqueous KI solution is known to yield I2 due to reaction of iodide ions with hydroxyl radicals, which in turn are generated due to cavitation. Based on this conceptual framework, a model has been developed to predict the rate of iodine formation for KI solutions of various concentrations under different gas atmospheres. The model follows the growth and collapse of a gas—vapour cavity using the Rayleigh—Plesset bubble dynamics equation. The bubble is assumed to behave isothermally during its growth phase and a part of the collapse phase. Thereafter it is assumed to collapse adiabatically, yielding high temperatures and pressures. Thermodynamic equilibrium is assumed in the bubble at the end of collapse phase. The contents of the bubble are assumed to mix with the liquid, and the reactor contents are assumed to be well stirred. The model has been verified by conducting experiments with KI solutions of different concentrations and using different gas atmospheres. The model not only explains these results but also the existence of a maximum when Ar---O2 mixtures of different compositions are employed.

Modelling of a batch sonochemical reactor

Ultrasonication of aqueous KI solution is known to yield I2 due to reaction of iodide ions with hydroxyl radicals, which in turn are generated due to cavitation. Based on this conceptual framework, a model has been developed to predict the rate of iodine formation for KI solutions of various concentrations under different gas atmospheres. The model follows the growth and collapse of a gas-vapour cavity using the Rayleigh-Plesset bubble dynamics equation. The bubble is assumed to behave isothermally during its growth phase and a part of the collapse phase. Thereafter it is assumed to collapse adiabatically, yielding high temperatures and pressures. Thermodynamic equilibrium is assumed in the bubble at the end of collapse phase. The contents of the bubble are assumed to mix with the liquid, and the reactor contents are assumed to be well stirred. The model has been verified by conducting experiments with KI solutions of different concentrations and using different gas atmospheres. The model not only explains these results but also the existence of a maximum when Ar-O2 mixtures of different compositions are employed.

Temperature effect on hydrocarbon-enhanced nitric oxide conversion using a dielectric barrier discharge reactor

Experimental investigations into the effect of temperature on conversion of NO in the presence of hydrocarbons (ethylene, acetylene and n-hexane) are presented. An AC energized dielectric barrier discharge reactor was used as the plasma reactor. The experiments were carried out at different temperatures up to 200 degreesC. The discharge powers were measured at all the temperatures. The discharge power was found to increase with temperature. NO conversion in the presence of ethylene and n-hexane was better than that of acetylene at all temperatures. The addition of acetylene at room temperature showed no better conversion of NO compared to no additive case. While at higher temperatures, it could enhance the conversion of NO. A slight enhancement in NO and NOx removal was observed in the presence of water vapor. (C) 2003 Elsevier Science B.V. All rights reserved.

Catalyst-Packed Non-Thermal Plasma Reactor for Removal of Nitrogen Oxides

A single-stage plasma-catalytic reactor in which catalytic materials were packed was used to remove nitrogen oxides. The packing material was scoria being made of various metal oxides including Al2O3, MgO, TiO2, etc. Scoria was able to act not only as dielectric pellets but also as a catalyst in the presence of reducing agent such as ethylene and ammonia. Without plasma discharge, scoria did not work well as a catalyst in the temperature range of 100 °C to 200 °C, showing less than 10% of NOx removal efficiency. When plasma is produced inside the reactor, the NOx removal efficiency could be increased to 60% in this temperature range.