Modelling a biochemical reaction with computational fluid dynamics

Earlier investigations (Cartland Glover et al., 2004) into the use of computational fluid dynamics (CFD) for the modelling of gas-liquid and gas-liquid-solid flow allowed a simple biochemical reaction model to be implemented. A single plane mesh was used to represent the transport and reaction of molasses, the mould Aspergillus niger and citric acid in a bubble column with a height to diameter aspect ratio of 20:1. Two specific growth rates were used to examine the impact that biomass growth had on the local solids concentration and the effect this had on the local hydrodynamics of the bubble column.

Kinetic study of the pyrolysis of miscanthus and its acid hydrolysis residue by thermogravimetric analysis

The kinetic parameters of the pyrolysis of miscanthus and its acid hydrolysis residue (AHR) were determined using thermogravimetric analysis (TGA). The AHR was produced at the University of Limerick by treating miscanthus with 5 wt.% sulphuric acid at 175 °C as representative of a lignocellulosic acid hydrolysis product. For the TGA experiments, 3 to 6 g of sample, milled and sieved to a particle size below 250 μm, were placed in the TGA ceramic crucible. The experiments were carried out under non-isothermal conditions heating the samples from 50 to 900 °C at heating rates of 2.5, 5, 10, 17 and 25 °C/min. The activation energy (EA) of the decomposition process was determined from the TGA data by differential analysis (Friedman) and three isoconversional methods of integral analysis (Kissinger–Akahira–Sunose, Ozawa–Flynn–Wall, Vyazovkin). The activation energy ranged from 129 to 156 kJ/mol for miscanthus and from 200 to 376 kJ/mol for AHR increasing with increasing conversion. The reaction model was selected using the non-linear least squares method and the pre-exponential factor was calculated from the Arrhenius approximation. The results showed that the best fitting reaction model was the third order reaction for both feedstocks. The pre-exponential factor was in the range of 5.6 × 1010 to 3.9 × 10+ 13 min− 1 for miscanthus and 2.1 × 1016 to 7.7 × 1025 min− 1 for AHR.

Thermal processing of miscanthus, sugarcane bagasse, sugarcane trash and their acid hydrolysis residues

The research presented in this thesis was developed as part of DIBANET, an EC funded project aiming to develop an energetically self-sustainable process for the production of diesel miscible biofuels (i.e. ethyl levulinate) via acid hydrolysis of selected biomass feedstocks. Three thermal conversion technologies, pyrolysis, gasification and combustion, were evaluated in the present work with the aim of recovering the energy stored in the acid hydrolysis solid residue (AHR). Mainly consisting of lignin and humins, the AHR can contain up to 80% of the energy in the original feedstock. Pyrolysis of AHR proved unsatisfactory, so attention focussed on gasification and combustion with the aim of producing heat and/or power to supply the energy demanded by the ethyl levulinate production process. A thermal processing rig consisting on a Laminar Entrained Flow Reactor (LEFR) equipped with solid and liquid collection and online gas analysis systems was designed and built to explore pyrolysis, gasification and air-blown combustion of AHR. Maximum liquid yield for pyrolysis of AHR was 30wt% with volatile conversion of 80%. Gas yield for AHR gasification was 78wt%, with 8wt% tar yields and conversion of volatiles close to 100%. 90wt% of the AHR was transformed into gas by combustion, with volatile conversions above 90%. 5volO2%-95vol%N2 gasification resulted in a nitrogen diluted, low heating value gas (2MJ/m3). Steam and oxygen-blown gasification of AHR were additionally investigated in a batch gasifier at KTH in Sweden. Steam promoted the formation of hydrogen (25vol%) and methane (14vol%) improving the gas heating value to 10MJ/m3, below the typical for steam gasification due to equipment limitations. Arrhenius kinetic parameters were calculated using data collected with the LEFR to provide reaction rate information for process design and optimisation. Activation energy (EA) and pre-exponential factor (ko in s-1) for pyrolysis (EA=80kJ/mol, lnko=14), gasification (EA=69kJ/mol, lnko=13) and combustion (EA=42kJ/mol, lnko=8) were calculated after linearly fitting the data using the random pore model. Kinetic parameters for pyrolysis and combustion were also determined by dynamic thermogravimetric analysis (TGA), including studies of the original biomass feedstocks for comparison. Results obtained by differential and integral isoconversional methods for activation energy determination were compared. Activation energy calculated by the Vyazovkin method was 103-204kJ/mol for pyrolysis of untreated feedstocks and 185-387kJ/mol for AHRs. Combustion activation energy was 138-163kJ/mol for biomass and 119-158 for AHRs. The non-linear least squares method was used to determine reaction model and pre-exponential factor. Pyrolysis and combustion of biomass were best modelled by a combination of third order reaction and 3 dimensional diffusion models, while AHR decomposed following the third order reaction for pyrolysis and the 3 dimensional diffusion for combustion.

Parametric pattern selection in a reaction-diffusion model

We compare spot patterns generated by Turing mechanisms with those generated by replication cascades, in a model one-dimensional reaction-diffusion system. We determine the stability region of spot solutions in parameter space as a function of a natural control parameter (feed-rate) where degenerate patterns with different numbers of spots coexist for a fixed feed-rate. While it is possible to generate identical patterns via both mechanisms, we show that replication cascades lead to a wider choice of pattern profiles that can be selected through a tuning of the feed-rate, exploiting hysteresis and directionality effects of the different pattern pathways.

Simultaneous chemical reaction and distillation of formaldehyde

The thesis is concerned with the development and testing of a mathematical model of a distillation process in which the components react chemically. The formaldehyde-methanol-water system was selected and only the reversible reactions between formaldehyde and water giving methylene glycol and between formaldehyde and methanol producing hemiformal were assumed to occur under the distillation conditions. Accordingly the system has been treated as a five component system. The vapour-liquid equilibrium calculations were performed by solving iteratively the thermodynamic relationships expressing the phase equilibria with the stoichiometric equations expressing the chemical equilibria. Using optimisation techniques, the Wilson single parameters and Henry's constants were calculated for binary systems containing formaldehyde which was assumed to be a supercritical component whilst Wilson binary parameters were calculated for the remaining binary systems. Thus the phase equilibria for the formaldehyde system could be calculated using these parameters and good accuracy was obtained when calculated values were compared with experimental values. The distillation process was modelled using the mass and energy balance equations together with the phase equilibria calculations. The plate efficiencies were obtained from a modified A.I.Ch.E. Bubble Tray method. The resulting equations were solved by an iterative plate to plate calculation based on the Newton Raphson method. Experiments were carried out in a 76mm I.D., eight sieve plate distillation column and the results were compared with the mathematical model calculations. Overall, good agreement was obtained but some discrepancies were observed in the concentration profiles and these may have been caused by the effect of limited physical property data and a limited understanding of the reactions mechanism. The model equations were solved in the form of modular computer programs. Although they were written to describe the steady state distillation with simultaneous chemical reaction of the formaldehyde system, the approach used may be of wider application.

Bio-inspired carbon electro-catalysis for the oxygen reduction reaction

We report the synthesis, characterisation and catalytic performance of two nature-inspired biomass-derived electro-catalysts for the oxygen reduction reaction in fuel cells. The catalysts were prepared via pyrolysis of a real food waste (lobster shells) or by mimicking the composition of lobster shells using chitin and CaCO3 particles followed by acid washing. The simplified model of artificial lobster was prepared for better reproducibility. The calcium carbonate in both samples acts as a pore agent, creating increased surface area and pore volume, though considerably higher in artificial lobster samples due to the better homogeneity of the components. Various characterisation techniques revealed the presence of a considerable amount of hydroxyapatite left in the real lobster samples after acid washing and a low content of carbon (23%), nitrogen and sulphur (<1%), limiting the surface area to 23 m2/g, and consequently resulting in rather poor catalytic activity. However, artificial lobster samples, with a surface area of ≈200 m2/g and a nitrogen doping of 2%, showed a promising onset potential, very similar to a commercially available platinum catalyst, with better methanol tolerance, though with lower stability in long time testing over 10,000 s.

Continuous chromatographic biochemical reaction-separation

Combined bioreaction separation studies have been carried out for the first time on a moving port semi-continuous counter-current chromatographic reactor-separator (SCCR-S1) consisting of twelve 5.4cm id x 75cm long columns packed with calcium charged cross-linked polystyrene resin (KORELA V07C). The inversion of sucrose to glucose and fructose in the presence of the enzyme invertase and the biochemIcal synthesis of dextran and fructose from sucrose in the presence of the enzyme dextransucrase were investigated. A dilute stream of the appropriate enzyme in deionised water was used as the eluent stream. The effect of switch time, feed concentration, enzyme activity, eluent rate and enzyme to feed concentration ratio on the combined bioreaction-separation were investigated. For the invertase reaction, at 20.77% w/v sucrose feed concentrations complete conversions were achieved. The enzyme usage was 34% of the theoretical enzyme amount needed to convert an equivalent amount of sucrose over the same time period when using a conventional fermenter. The fructose rich (FRP) and glucose rich (GRP) product purities obtained were over 90%. By operating at 35% w/v sucrose feed concentration and employing the product splitting and recycling techniques, the total concentration and purity of the GRP increased from 32% w/v to 4.6% and from 92.3% to 95% respectively. The FRP concentration also increased from 1.82% w/v to 2.88% w/v. A mathematical model was developed for the combined reaction-separation and used to simulate the continuous inversion of sucrose and product separation using the SCCR-S1. In the biosynthesis of dextran studies, 52% conversion of a 2% w/v sucrose concentration feed was achieved. An average dextran molecular weight of 4 millIon was obtained in the dextran rich (DRP) product stream. The enzyme dextransucrase was purifed successfully using centrifugation and ultrafiltration techniques.

Uncatalysed and potassium-catalysed pyrolysis of the cell-wall constituents of biomass and their model compounds

Cell-wall components (cellulose, hemicellulose (oat spelt xylan), lignin (Organosolv)), and model compounds (levoglucosan (an intermediate product of cellulose decomposition) and chlorogenic acid (structurally similar to lignin polymer units)) have been investigated to probe in detail the influence of potassium on their pyrolysis behaviours as well as their uncatalysed decomposition reaction. Cellulose and lignin were pretreated to remove salts and metals by hydrochloric acid, and this dematerialized sample was impregnated with 1% of potassium as potassium acetate. Levoglucosan, xylan and chlorogenic acid were mixed with CHCOOK to introduce 1% K. Characterisation was performed using thermogravimetric analysis (TGA) and differential thermal analysis (DTA). In addition to the TGA pyrolysis, pyrolysis-gas chromatography-mass spectrometry (PY-GC-MS) analysis was introduced to examine reaction products. Potassium-catalysed pyrolysis has a huge influence on the char formation stage and increases the char yields considerably (from 7.7% for raw cellulose to 27.7% for potassium impregnated cellulose; from 5.7% for raw levoglucosan to 20.8% for levoglucosan with CHCOOK added). Major changes in the pyrolytic decomposition pathways were observed for cellulose, levoglucosan and chlorogenic acid. The results for cellulose and levoglucosan are consistent with a base catalysed route in the presence of the potassium salt which promotes complete decomposition of glucosidic units by a heterolytic mechanism and favours its direct depolymerization and fragmentation to low molecular weight components (e.g. acetic acid, formic acid, glyoxal, hydroxyacetaldehyde and acetol). Base catalysed polymerization reactions increase the char yield. Potassium-catalysed lignin pyrolysis is very significant: the temperature of maximum conversion in pyrolysis shifts to lower temperature by 70 K and catalysed polymerization reactions increase the char yield from 37% to 51%. A similar trend is observed for the model compound, chlorogenic acid. The addition of potassium does not produce a dramatic change in the tar product distribution, although its addition to chlorogenic acid promoted the generation of cyclohexane and phenol derivatives. Postulated thermal decomposition schemes for chlorogenic acid are presented. © 2008 Elsevier B.V. All rights reserved.

Modeling competitive reaction effects

In this study I critically review models that specify competitive reaction effects. I discuss different model structures and summarize my findings on competitive reaction effects and factors that explain competitive reactions. I discuss the many models of competitive market response that have been developed and classify them into twelve sets of models that are related to each other in a logical manner through the evolutionary model-building concept.

Multiscale modelling of heterogeneously catalysed transesterification reaction process:an overview

Biodiesel is fast becoming one of the key transport fuels as the world endeavours to reduce its carbon footprint and find viable alternatives to oil derived fuels. Research in the field is currently focusing on more efficient ways to produce biodiesel, with the most promising avenue of research looking into the use of heterogeneous catalysis. This article presents a framework for kinetic reaction and diffusive transport modelling of the heterogeneously catalysed transesterification of triglycerides into fatty acid methyl esters (FAMEs), unveiled by a model system of tributyrin transesterification in the presence of MgO catalysts. In particular, the paper makes recommendations on multicomponent diffusion calculations such as the diffusion coefficients and molar fluxes from infinite dilution diffusion coefficients using the Wilke and Chang correlation, intrinsic reaction kinetic studies using the Eley-Rideal kinetic mechanism with methanol adsorption as the rate determining steps and multiscale reaction-diffusion process simulation between catalytic porous and bulk reactor scales. © 2013 The Royal Society of Chemistry.

Hydrogen production by steam reforming of DME in a large scale CFB reactor. Part I:computational model and predictions

This study presents a computational fluid dynamic (CFD) study of Dimethyl Ether steam reforming (DME-SR) in a large scale Circulating Fluidized Bed (CFB) reactor. The CFD model is based on Eulerian-Eulerian dispersed flow and solved using commercial software (ANSYS FLUENT). The DME-SR reactions scheme and kinetics in the presence of a bifunctional catalyst of CuO/ZnO/Al2O3+ZSM-5 were incorporated in the model using in-house developed user-defined function. The model was validated by comparing the predictions with experimental data from the literature. The results revealed for the first time detailed CFB reactor hydrodynamics, gas residence time, temperature distribution and product gas composition at a selected operating condition of 300 °C and steam to DME mass ratio of 3 (molar ratio of 7.62). The spatial variation in the gas species concentrations suggests the existence of three distinct reaction zones but limited temperature variations. The DME conversion and hydrogen yield were found to be 87% and 59% respectively, resulting in a product gas consisting of 72 mol% hydrogen. In part II of this study, the model presented here will be used to optimize the reactor design and study the effect of operating conditions on the reactor performance and products.

Influence of reaction atmosphere (H2O, N2, H2, CO2, CO) on fluidized-bed fast pyrolysis of biomass using detailed tar vapor chemistry in computational fluid dynamics

Secondary pyrolysis in fluidized bed fast pyrolysis of biomass is the focus of this work. A novel computational fluid dynamics (CFD) model coupled with a comprehensive chemistry scheme (134 species and 4169 reactions, in CHEMKIN format) has been developed to investigate this complex phenomenon. Previous results from a transient three-dimensional model of primary pyrolysis were used for the source terms of primary products in this model. A parametric study of reaction atmospheres (H2O, N2, H2, CO2, CO) has been performed. For the N2 and H2O atmosphere, results of the model compared favorably to experimentally obtained yields after the temperature was adjusted to a value higher than that used in experiments. One notable deviation versus experiments is pyrolytic water yield and yield of higher hydrocarbons. The model suggests a not overly strong impact of the reaction atmosphere. However, both chemical and physical effects were observed. Most notably, effects could be seen on the yield of various compounds, temperature profile throughout the reactor system, residence time, radical concentration, and turbulent intensity. At the investigated temperature (873 K), turbulent intensity appeared to have the strongest influence on liquid yield. With the aid of acceleration techniques, most importantly dimension reduction, chemistry agglomeration, and in-situ tabulation, a converged solution could be obtained within a reasonable time (∼30 h). As such, a new potentially useful method has been suggested for numerical analysis of fast pyrolysis.

Fluctuation-driven traffic congestion in a scale-free model of the Internet

In studies of complex heterogeneous networks, particularly of the Internet, significant attention was paid to analysing network failures caused by hardware faults or overload. There network reaction was modelled as rerouting of traffic away from failed or congested elements. Here we model network reaction to congestion on much shorter time scales when the input traffic rate through congested routes is reduced. As an example we consider the Internet where local mismatch between demand and capacity results in traffic losses. We describe the onset of congestion as a phase transition characterised by strong, albeit relatively short-lived, fluctuations of losses caused by noise in input traffic and exacerbated by the heterogeneous nature of the network manifested in a power-law load distribution. The fluctuations may result in the network strongly overreacting to the first signs of congestion by significantly reducing input traffic along the communication paths where congestion is utterly negligible. © 2013 IEEE.

Effect of Cu and Sn promotion on the catalytic deoxygenation of model and algal lipids to fuel-like hydrocarbons over supported Ni catalysts

The ability of Cu and Sn to promote the performance of a 20% Ni/Al2O3 catalyst in the deoxygenation of lipids to fuel-like hydrocarbons was investigated using model triglyceride and fatty acid feeds, as well as algal lipids. In the semi-batch deoxygenation of tristearin at 260 °C a pronounced promotional effect was observed, a 20% Ni-5% Cu/Al2O3 catalyst affording both higher conversion (97%) and selectivity to C10-C17 alkanes (99%) in comparison with unpromoted 20% Ni/Al2O3 (27% conversion and 87% selectivity to C10-C17). In the same reaction at 350 °C, a 20% Ni-1% Sn/Al2O3 catalyst afforded the best results, giving yields of C10-C17 and C17 of 97% and 55%, respectively, which contrasts with the corresponding values of 87 and 21% obtained over 20% Ni/Al2O3. Equally encouraging results were obtained in the semi-batch deoxygenation of stearic acid at 300 °C, in which the 20% Ni-5% Cu/Al2O3 catalyst afforded the highest yields of C10-C17 and C17. Experiments were also conducted at 260 °C in a fixed bed reactor using triolein − a model unsaturated triglyceride − as the feed. While both 20% Ni/Al2O3 and 20% Ni-5% Cu/Al2O3 achieved quantitative yields of diesel-like hydrocarbons at all reaction times sampled, the Cu-promoted catalyst exhibited higher selectivity to longer chain hydrocarbons, a phenomenon which was also observed in experiments involving algal lipids as the feed. Characterization of fresh and spent catalysts indicates that Cu enhances the reducibility of Ni and suppresses both cracking reactions and coke-induced deactivation.

Experimental validation of a predictive pyrolysis model in AspenPlus®

A novel simulation model for pyrolysis processes oflignocellulosicbiomassin AspenPlus (R) was presented at the BC&E 2013. Based on kinetic reaction mechanisms, the simulation calculates product compositions and yields depending on reactor conditions (temperature, residence time, flue gas flow rate) and feedstock composition (biochemical composition, atomic composition, ash and alkali metal content). The simulation model was found to show good correlation with existing publications. In order to further verify the model, own pyrolysis experiments in a 1 kg/h continuously fed fluidized bed fast pyrolysis reactor are performed. Two types of biomass with different characteristics are processed in order to evaluate the influence of the feedstock composition on the yields of the pyrolysis products and their composition. One wood and one straw-like feedstock are used due to their different characteristics. Furthermore, the temperature response of yields and product compositions is evaluated by varying the reactor temperature between 450 and 550 degrees C for one of the feedstocks. The yields of the pyrolysis products (gas, oil, char) are determined and their detailed composition is analysed. The experimental runs are reproduced with the corresponding reactor conditions in the AspenPlus model and the results compared with the experimental findings.