Analysis and modelling of chemical looping combustion process with and without oxygen uncoupling

Effective control and limiting of carbon dioxide (CO₂) emissions in energy production are major challenges of science today. Current research activities include the development of new low-cost carbon capture technologies, and among the proposed concepts, chemical combustion (CLC) and chemical looping with oxygen uncoupling (CLOU) have attracted significant attention allowing intrinsic separation of pure CO₂ from a hydrocarbon fuel combustion process with a comparatively small energy penalty. Both CLC and CLOU utilize the well-established fluidized bed technology, but several technical challenges need to be overcome in order to commercialize the processes. Therefore, development of proper modelling and simulation tools is essential for the design, optimization, and scale-up of chemical looping-based combustion systems. The main objective of this work was to analyze the technological feasibility of CLC and CLOU processes at different scales using a computational modelling approach. A onedimensional fluidized bed model frame was constructed and applied for simulations of CLC and CLOU systems consisting of interconnected fluidized bed reactors. The model is based on the conservation of mass and energy, and semi-empirical correlations are used to describe the hydrodynamics, chemical reactions, and transfer of heat in the reactors. Another objective was to evaluate the viability of chemical looping-based energy production, and a flow sheet model representing a CLC-integrated steam power plant was developed. The 1D model frame was succesfully validated based on the operation of a 150 kWth laboratory-sized CLC unit fed by methane. By following certain scale-up criteria, a conceptual design for a CLC reactor system at a pre-commercial scale of 100 MWth was created, after which the validated model was used to predict the performance of the system. As a result, further understanding of the parameters affecting the operation of a large-scale CLC process was acquired, which will be useful for the practical design work in the future. The integration of the reactor system and steam turbine cycle for power production was studied resulting in a suggested plant layout including a CLC boiler system, a simple heat recovery setup, and an integrated steam cycle with a three pressure level steam turbine. Possible operational regions of a CLOU reactor system fed by bituminous coal were determined via mass, energy, and exergy balance analysis. Finally, the 1D fluidized bed model was modified suitable for CLOU, and the performance of a hypothetical 500 MWth CLOU fuel reactor was evaluated by extensive case simulations.

Att använda Choice Reaction Time (CRT) paradigmet för att differentiera mellan sexualförbrytare med barnoffer, sexualförbrytare med vuxna offer, förbrytare skyldiga till andra brott än sexualbrott, samt kontroller från normalpopulationen

Pro graduavhanlingens svenska sammanfattning

Quantitative refinement of reaction-based biomodels

In the field of molecular biology, scientists adopted for decades a reductionist perspective in their inquiries, being predominantly concerned with the intricate mechanistic details of subcellular regulatory systems. However, integrative thinking was still applied at a smaller scale in molecular biology to understand the underlying processes of cellular behaviour for at least half a century. It was not until the genomic revolution at the end of the previous century that we required model building to account for systemic properties of cellular activity. Our system-level understanding of cellular function is to this day hindered by drastic limitations in our capability of predicting cellular behaviour to reflect system dynamics and system structures. To this end, systems biology aims for a system-level understanding of functional intraand inter-cellular activity. Modern biology brings about a high volume of data, whose comprehension we cannot even aim for in the absence of computational support. Computational modelling, hence, bridges modern biology to computer science, enabling a number of assets, which prove to be invaluable in the analysis of complex biological systems, such as: a rigorous characterization of the system structure, simulation techniques, perturbations analysis, etc. Computational biomodels augmented in size considerably in the past years, major contributions being made towards the simulation and analysis of large-scale models, starting with signalling pathways and culminating with whole-cell models, tissue-level models, organ models and full-scale patient models. The simulation and analysis of models of such complexity very often requires, in fact, the integration of various sub-models, entwined at different levels of resolution and whose organization spans over several levels of hierarchy. This thesis revolves around the concept of quantitative model refinement in relation to the process of model building in computational systems biology. The thesis proposes a sound computational framework for the stepwise augmentation of a biomodel. One starts with an abstract, high-level representation of a biological phenomenon, which is materialised into an initial model that is validated against a set of existing data. Consequently, the model is refined to include more details regarding its species and/or reactions. The framework is employed in the development of two models, one for the heat shock response in eukaryotes and the second for the ErbB signalling pathway. The thesis spans over several formalisms used in computational systems biology, inherently quantitative: reaction-network models, rule-based models and Petri net models, as well as a recent formalism intrinsically qualitative: reaction systems. The choice of modelling formalism is, however, determined by the nature of the question the modeler aims to answer. Quantitative model refinement turns out to be not only essential in the model development cycle, but also beneficial for the compilation of large-scale models, whose development requires the integration of several sub-models across various levels of resolution and underlying formal representations.

Bio and chemocatalysis for stereo and regioselective one-pot reaction applications

The development of cost efficient, selective and sustainable chemical processes for production of chiral building blocks is of great importance in synthetic and industrial organic chemistry. One way to reach these objectives is to carry out several reactions steps in one vessel at one time. Furthermore, when this kind of one-pot multi step reactions are catalyzed by heterogeneous chemo- and bio-catalysts, which can be separated from the reaction products by filtration, practical access to chiral small molecules for further utilization can be obtained. The initial reactions studied in this thesis are the two step dynamic kinetic resolution of rac-2-hydroxy-1-indanone and the regioselective hydrogenation of 1,2-indanedione. These reactions are then combined in a new heterogeneously catalyzed one-pot reaction sequence enabling simple recovery of the catalysts by filtration, facilitating simple reaction product isolation. Conclusively, the readily available 1,2-indanedione is by the presented one-pot sequence, utilizing heterogeneous enzyme and transition metal based catalysts, transferred with high regio- and stereoselectivity to a useful chiral vicinal hydroxyl ketone structure. Additional and complementary investigation of homogeneous half-sandwich ruthenium complexes for catalyzing the epimerization of chiral secondary alcohols of five natural products containing additional non-functionalized stereocenters was conducted. In principle, this kind of epimerization reactions of single stereocenters could be utilized for converting inexpensive starting materials, containing other stereogenic centers, into diastereomeric mixtures from which more valuable compounds can be isolated by traditional isolation techniques.

Implementation of Quantitative Reaction-based Model Refinement

Building a computational model for complex biological systems is an iterative process. It starts from an abstraction of the process and then incorporates more details regarding the specific biochemical reactions which results in the change of the model fit. Meanwhile, the model’s numerical properties such as its numerical fit and validation should be preserved. However, refitting the model after each refinement iteration is computationally expensive resource-wise. There is an alternative approach which ensures the model fit preservation without the need to refit the model after each refinement iteration. And this approach is known as quantitative model refinement. The aim of this thesis is to develop and implement a tool called ModelRef which does the quantitative model refinement automatically. It is both implemented as a stand-alone Java application and as one of Anduril framework components. ModelRef performs data refinement of a model and generates the results in two different well known formats (SBML and CPS formats). The development of this tool successfully reduces the time and resource needed and the errors generated as well by traditional reiteration of the whole model to perform the fitting procedure.

Computational models for and from biology : simple gene assembly and reaction systems

The advancement of science and technology makes it clear that no single perspective is any longer sufficient to describe the true nature of any phenomenon. That is why the interdisciplinary research is gaining more attention overtime. An excellent example of this type of research is natural computing which stands on the borderline between biology and computer science. The contribution of research done in natural computing is twofold: on one hand, it sheds light into how nature works and how it processes information and, on the other hand, it provides some guidelines on how to design bio-inspired technologies. The first direction in this thesis focuses on a nature-inspired process called gene assembly in ciliates. The second one studies reaction systems, as a modeling framework with its rationale built upon the biochemical interactions happening within a cell. The process of gene assembly in ciliates has attracted a lot of attention as a research topic in the past 15 years. Two main modelling frameworks have been initially proposed in the end of 1990s to capture ciliates’ gene assembly process, namely the intermolecular model and the intramolecular model. They were followed by other model proposals such as templatebased assembly and DNA rearrangement pathways recombination models. In this thesis we are interested in a variation of the intramolecular model called simple gene assembly model, which focuses on the simplest possible folds in the assembly process. We propose a new framework called directed overlap-inclusion (DOI) graphs to overcome the limitations that previously introduced models faced in capturing all the combinatorial details of the simple gene assembly process. We investigate a number of combinatorial properties of these graphs, including a necessary property in terms of forbidden induced subgraphs. We also introduce DOI graph-based rewriting rules that capture all the operations of the simple gene assembly model and prove that they are equivalent to the string-based formalization of the model. Reaction systems (RS) is another nature-inspired modeling framework that is studied in this thesis. Reaction systems’ rationale is based upon two main regulation mechanisms, facilitation and inhibition, which control the interactions between biochemical reactions. Reaction systems is a complementary modeling framework to traditional quantitative frameworks, focusing on explicit cause-effect relationships between reactions. The explicit formulation of facilitation and inhibition mechanisms behind reactions, as well as the focus on interactions between reactions (rather than dynamics of concentrations) makes their applicability potentially wide and useful beyond biological case studies. In this thesis, we construct a reaction system model corresponding to the heat shock response mechanism based on a novel concept of dominance graph that captures the competition on resources in the ODE model. We also introduce for RS various concepts inspired by biology, e.g., mass conservation, steady state, periodicity, etc., to do model checking of the reaction systems based models. We prove that the complexity of the decision problems related to these properties varies from P to NP- and coNP-complete to PSPACE-complete. We further focus on the mass conservation relation in an RS and introduce the conservation dependency graph to capture the relation between the species and also propose an algorithm to list the conserved sets of a given reaction system.

Chemical changes in biomass during torrefaction

Torrefaction is moderate thermal treatment (~200-300 °C) of biomass in an inert atmosphere. The torrefied fuel offers advantages to traditional biomass, such as higher heating value, reduced hydrophilic nature, increased its resistance to biological decay, and improved grindability. These factors could, for instance, lead to better handling and storage of biomass and increased use of biomass in pulverized combustors. In this work, we look at several aspects of changes in the biomass during torrefaction. We investigate the fate of carboxylic groups during torrefaction and its dependency to equilibrium moisture content. The changes in the wood components including carbohydrates, lignin, extractable materials and ashforming matters are also studied. And at last, the effect of K on torrefaction is investigated and then modeled. In biomass, carboxylic sites are partially responsible for its hydrophilic characteristic. These sites are degraded to varying extents during torrefaction. In this work, methylene blue sorption and potentiometric titration were applied to measure the concentration of carboxylic groups in torrefied spruce wood. The results from both methods were applicable and the values agreed well. A decrease in the equilibrium moisture content at different humidity was also measured for the torrefied wood samples, which is in good agreement with the decrease in carboxylic group contents. Thus, both methods offer a means of directly measuring the decomposition of carboxylic groups in biomass during torrefaction as a valuable parameter in evaluating the extent of torrefaction. This provides new information to the chemical changes occurring during torrefaction. The effect of torrefaction temperature on the chemistry of birch wood was investigated. The samples were from a pilot plant at Energy research Center of the Netherlands (ECN). And in that way they were representative of industrially produced samples. Sugar analysis was applied to analyze the hemicellulose and cellulose content during torrefaction. The results show a significant degradation of hemicellulose already at 240 °C, while cellulose degradation becomes significant above 270 °C torrefaction. Several methods including Klason lignin method, solid state NMR and Py-GC-MS analyses were applied to measure the changes in lignin during torrefaction. The changes in the ratio of phenyl, guaiacyl and syringyl units show that lignin degrades already at 240 °C to a small extent. To investigate the changes in the extractives from acetone extraction during torrefaction, gravimetric method, HP-SEC and GC-FID followed by GC-MS analysis were performed. The content of acetone-extractable material increases already at 240 °C torrefaction through the degradation of carbohydrate and lignin. The molecular weight of the acetone-extractable material decreases with increasing the torrefaction temperature. The formation of some valuable materials like syringaresinol or vanillin is also observed which is important from biorefinery perspective. To investigate the change in the chemical association of ash-forming elements in birch wood during torrefaction, chemical fractionation was performed on the original and torrefied birch samples. These results give a first understanding of the changes in the association of ashforming elements during torrefaction. The most significant changes can be seen in the distribution of calcium, magnesium and manganese, with some change in water solubility seen in potassium. These changes may in part be due to the destruction of carboxylic groups. In addition to some changes in water and acid solubility of phosphorous, a clear decrease in the concentration of both chlorine and sulfur was observed. This would be a significant additional benefit for the combustion of torrefied biomass. Another objective of this work is studying the impact of organically bound K, Na, Ca and Mn on mass loss of biomass during torrefaction. These elements were of interest because they have been shown to be catalytically active in solid fuels during pyrolysis and/or gasification. The biomasses were first acid washed to remove the ash-forming matters and then organic sites were doped with K, Na, Ca or Mn. The results show that K and Na bound to organic sites can significantly increase the mass loss during torrefaction. It is also seen that Mn bound to organic sites increases the mass loss and Ca addition does not influence the mass loss rate on torrefaction. This increase in mass loss during torrefaction with alkali addition is unlike what has been found in the case of pyrolysis where alkali addition resulted in a reduced mass loss. These results are important for the future operation of torrefaction plants, which will likely be designed to handle various biomasses with significantly different contents of K. The results imply that shorter retention times are possible for high K-containing biomasses. The mass loss of spruce wood with different content of K was modeled using a two-step reaction model based on four kinetic rate constants. The results show that it is possible to model the mass loss of spruce wood doped with different levels of K using the same activation energies but different pre-exponential factors for the rate constants. Three of the pre-exponential factors increased linearly with increasing K content, while one of the preexponential factors decreased with increasing K content. Therefore, a new torrefaction model was formulated using the hemicellulose and cellulose content and K content. The new torrefaction model was validated against the mass loss during the torrefaction of aspen, miscanthus, straw and bark. There is good agreement between the model and the experimental data for the other biomasses, except bark. For bark, the mass loss of acetone extractable material is also needed to be taken into account. The new model can describe the kinetics of mass loss during torrefaction of different types of biomass. This is important for considering fuel flexibility in torrefaction plants.

A reaction engineering approach to homogeneously catalyzed glycerol hydrochlorination

The production of biodiesel through transesterification has created a surplus of glycerol on the international market. In few years, glycerol has become an inexpensive and abundant raw material, subject to numerous plausible valorisation strategies. Glycerol hydrochlorination stands out as an economically attractive alternative to the production of biobased epichlorohydrin, an important raw material for the manufacturing of epoxy resins and plasticizers. Glycerol hydrochlorination using gaseous hydrogen chloride (HCl) was studied from a reaction engineering viewpoint. Firstly, a more general and rigorous kinetic model was derived based on a consistent reaction mechanism proposed in the literature. The model was validated with experimental data reported in the literature as well as with new data of our own. Semi-batch experiments were conducted in which the influence of the stirring speed, HCl partial pressure, catalyst concentration and temperature were thoroughly analysed and discussed. Acetic acid was used as a homogeneous catalyst for the experiments. For the first time, it was demonstrated that the liquid-phase volume undergoes a significant increase due to the accumulation of HCl in the liquid phase. Novel and relevant features concerning hydrochlorination kinetics, HCl solubility and mass transfer were investigated. An extended reaction mechanism was proposed and a new kinetic model was derived. The model was tested with the experimental data by means of regression analysis, in which kinetic and mass transfer parameters were successfully estimated. A dimensionless number, called Catalyst Modulus, was proposed as a tool for corroborating the kinetic model. Reactive flash distillation experiments were conducted to check the commonly accepted hypothesis that removal of water should enhance the glycerol hydrochlorination kinetics. The performance of the reactive flash distillation experiments were compared to the semi-batch data previously obtained. An unforeseen effect was observed once the water was let to be stripped out from the liquid phase, exposing a strong correlation between the HCl liquid uptake and the presence of water in the system. Water has revealed to play an important role also in the HCl dissociation: as water was removed, the dissociation of HCl was diminished, which had a retarding effect on the reaction kinetics. In order to obtain a further insight on the influence of water on the hydrochlorination reaction, extra semi-batch experiments were conducted in which initial amounts of water and the desired product were added. This study revealed the possibility to use the desired product as an ideal “solvent” for the glycerol hydrochlorination process. A co-current bubble column was used to investigate the glycerol hydrochlorination process under continuous operation. The influence of liquid flow rate, gas flow rate, temperature and catalyst concentration on the glycerol conversion and product distribution was studied. The fluid dynamics of the system showed a remarkable behaviour, which was carefully investigated and described. Highspeed camera images and residence time distribution experiments were conducted to collect relevant information about the flow conditions inside the tube. A model based on the axial dispersion concept was proposed and confronted with the experimental data. The kinetic and solubility parameters estimated from the semi-batch experiments were successfully used in the description of mass transfer and fluid dynamics of the bubble column reactor. In light of the results brought by the present work, the glycerol hydrochlorination reaction mechanism has been finally clarified. It has been demonstrated that the reactive distillation technology may cause drawbacks to the glycerol hydrochlorination reaction rate under certain conditions. Furthermore, continuous reactor technology showed a high selectivity towards monochlorohydrins, whilst semibatch technology was demonstrated to be more efficient towards the production of dichlorohydrins. Based on the novel and revealing discoveries brought by the present work, many insightful suggestions are made towards the improvement of the production of αγ-dichlorohydrin on an industrial scale.

Mass transfer and reaction kinetics in single droplet systems

As reactive extraction grown more and more popular in a variety of technological applications, optimizing its performance becomes more and more important. The process of complex formation is affected by a great number of both physical and chemical properties of all the components involved, and sometimes their interference with one another makes improving the effectiveness of such processes very difficult. In this Master’s Theses, the processes of complex formation between the aqueous phase - represented by copper sulfate water solution, and organic phase – represented by Acorga M5640 solvent extractor, were studied in order to establish the effect these components have on reactive extraction performance and to determine which step is bottlenecking the process the most.