A generalized 1D particle transport method for convection diffusion reaction model

This work is devoted to the development of numerical method to deal with convection diffusion dominated problem with reaction term, non - stiff chemical reaction and stiff chemical reaction. The technique is based on the unifying Eulerian - Lagrangian schemes (particle transport method) under the framework of operator splitting method. In the computational domain, the particle set is assigned to solve the convection reaction subproblem along the characteristic curves created by convective velocity. At each time step, convection, diffusion and reaction terms are solved separately by assuming that, each phenomenon occurs separately in a sequential fashion. Moreover, adaptivities and projection techniques are used to add particles in the regions of high gradients (steep fronts) and discontinuities and transfer a solution from particle set onto grid point respectively. The numerical results show that, the particle transport method has improved the solutions of CDR problems. Nevertheless, the method is time consumer when compared with other classical technique e.g., method of lines. Apart from this advantage, the particle transport method can be used to simulate problems that involve movingsteep/smooth fronts such as separation of two or more elements in the system.

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.

Synthesis of NiO-MgO-ZrO(2) catalysts and their performance in reforming of model biogas

Catalysts containing NiO/MgO/ZrO(2) mixtures were synthesized by the polymerization method in a single step. They were characterized by X-ray diffraction (XRD), temperature programmed reduction (TPR) and physisorption of N(2) (BET) and then tested in the reforming of a model biogas (1.5CH4:1CO(2)) in the presence of air (1.5CH(4) + 1CO(2) + 0.25O(2)) at 750 degrees C for 6h. It was observed that the catalyst Ni20MZ performed better in catalytic processes than the well known catalysts, Ni/ZrO(2) and Ni/MgO, synthesized under the same conditions. The formation of solid solutions, MgO-ZrO(2) and NiO-MgO, increased the rate of conversion of reactants (CH(4) and CO(2)) into synthesis gas (H(2) + CO). The formation of oxygen vacancies (in samples containing ZrO(2) and MgO) seems to promote removal of the coke deposited on the nickel surface. The values of the H(2)/CO ratio were generally found to be slightly lower than stoichiometric, owing to the reverse water gas shift reaction occurring in parallel. (C) 2011 Elsevier B.V. All rights reserved.

Local scale air quality model with several pollutant sources

[EN]This paper shows a finite element method for pollutant transport with several pollutant sources. An Eulerian convection–diffusion–reaction model to simulate the pollutant dispersion is used. The discretization of the different sources allows to impose the emissions as boundary conditions. The Eulerian description can deal with the coupling of several plumes. An adaptive stabilized finite element formulation, specifically Least-Squares, with a Crank-Nicolson temporal integration is proposed to solve the problem. An splitting scheme has been used to treat separately the transport and the reaction. A mass-consistent model has been used to compute the wind field of the problem…

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.

Nanohybrid materials as transducer surfaces for electrochemical sensing applications

A nanohybrid electrochemical transducer surface was developed using carbon and gold nanomaterials. The strategy relayed on casting multiwalled carbon nanotubes or carbon nanofibers onto a screen-printed carbon electrode surface, followed by in situ generation of gold nanoparticles by electrochemical deposition of ionic gold, in a reproducible manner. These transducers, so fabricated, were characterized using both electrochemical and microscopic techniques. Biofunctionality was evaluated using the streptavidin-biotin interaction system as the biological reaction model. These platforms allow to achieve low detection limits (in the order of pmoles), are reproducible and stable at least for a month after their preparation, being a perfect candidate to be used as transducer of different sensor devices.

Simulação de um processo de pirólise rápida de biomassa

Mestrado em Engenharia Química - Ramo Optimização Energética na Indústria Química

Comprehensive Study of Crystal Growth from Solution

Crystal growth is an essential phase in crystallization kinetics. The rate of crystal growth provides significant information for the design and control of crystallization processes; nevertheless, obtaining accurate growth rate data is still challenging due to a number of factors that prevail in crystal growth. In industrial crystallization, crystals are generally grown from multi-componentand multi-particle solutions under complicated hydrodynamic conditions; thus, it is crucial to increase the general understanding of the growth kinetics in these systems. The aim of this work is to develop a model of the crystal growth rate from solution. An extensive literature review of crystal growth focuses on themodelling of growth kinetics and thermodynamics, and new measuring techniques that have been introduced in the field of crystallization. The growth of a singlecrystal is investigated in binary and ternary systems. The binary system consists of potassium dihydrogen phosphate (KDP, crystallizing solute) and water (solvent), and the ternary system includes KDP, water and an organic admixture. The studied admixtures, urea, ethanol and 1-propanol, are employed at relatively highconcentrations (of up to 5.0 molal). The influence of the admixtures on the solution thermodynamics is studied using the Pitzer activity coefficient model. Theprediction method of the ternary solubility in the studied systems is introduced and verified. The growth rate of the KDP (101) face in the studied systems aremeasured in the growth cell as a function of supersaturation, the admixture concentration, the solution velocity over a crystal and temperature. In addition, the surface morphology of the KDP (101) face is studied using ex situ atomic force microscopy (AFM). The crystal growth rate in the ternary systems is modelled on the basis of the two-step growth model that contains the Maxwell-Stefan (MS) equations and a surface-reaction model. This model is used together with measuredcrystal growth rate data to develop a new method for the evaluation of the model parameters. The validation of the model is justified with experiments. The crystal growth rate in an imperfectly mixed suspension crystallizer is investigatedusing computational fluid dynamics (CFD). A solid-liquid suspension flow that includes multi-sized particles is described by the multi-fluid model as well as by a standard k-epsilon turbulence model and an interface momentum transfer model. The local crystal growth rate is determined from calculated flow information in a diffusion-controlled crystal growth regime. The calculated results are evaluated experimentally.

Bayesian methods for estimation, optimization and experimental design

Mathematical models often contain parameters that need to be calibrated from measured data. The emergence of efficient Markov Chain Monte Carlo (MCMC) methods has made the Bayesian approach a standard tool in quantifying the uncertainty in the parameters. With MCMC, the parameter estimation problem can be solved in a fully statistical manner, and the whole distribution of the parameters can be explored, instead of obtaining point estimates and using, e.g., Gaussian approximations. In this thesis, MCMC methods are applied to parameter estimation problems in chemical reaction engineering, population ecology, and climate modeling. Motivated by the climate model experiments, the methods are developed further to make them more suitable for problems where the model is computationally intensive. After the parameters are estimated, one can start to use the model for various tasks. Two such tasks are studied in this thesis: optimal design of experiments, where the task is to design the next measurements so that the parameter uncertainty is minimized, and model-based optimization, where a model-based quantity, such as the product yield in a chemical reaction model, is optimized. In this thesis, novel ways to perform these tasks are developed, based on the output of MCMC parameter estimation. A separate topic is dynamical state estimation, where the task is to estimate the dynamically changing model state, instead of static parameters. For example, in numerical weather prediction, an estimate of the state of the atmosphere must constantly be updated based on the recently obtained measurements. In this thesis, a novel hybrid state estimation method is developed, which combines elements from deterministic and random sampling methods.

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.

Influence of adsorption geometry in the heterogeneous enantioselective catalytic hydrogenation of a prototypical enone

Asymmetric catalysis is of paramount importance in organic synthesis and, in current practice, is achieved by means of homogeneous catalysts. The ability to catalyze such reactions heterogeneously would have a major impact both in the research laboratory and in the production of fine chemicals and pharmaceuticals, yet heterogeneous asymmetric hydrogenation of C═C bonds remains hardly explored. Very recently, we demonstrated how chiral ligands that anchor robustly to the surface of Pd nanoparticles promote asymmetric catalytic hydrogenation: ligand rigidity and stereochemistry emerged as key factors. Here, we address a complementary question: how does the enone reactant adsorb on the metal surface, and what implications does this have for the enantiodifferentiating interaction with the surface-tethered chiral modifiers? A reaction model is proposed, which correctly predicts the identity of the enantiomer experimentally observed in excess.

Electrochemical studies of the interface Fe/0.5 mol L(-1) H(2)SO(4) in the presence of benzotriazole and tolytriazole

The effect of benzotriazole (BTAH) and tolytriazole (TTAH) on the electrochemical behaviour of the Fe/0.5 mol L(-1) H(2)SO(4) interface at 25 degrees C was studied using cronopotentiometry, anodic and cathodic polarization curves and electrochemical impedance spectroscopy. BTAH and TTAH are inhibitors of anodic iron dissolution and the subsequent hydrogen evolution in 0.5 mol L(-1) H(2)SO(4) medium. Mass transport is an important step in the anodic process of inhibitive film formation. Electrochemical impedance spectroscopy was used to investigate the iron dissolution mechanism in the presence of the inhibitors and showed that BTAH and TTAH are adsorbed on the iron surface, thereby changing its dissolution mechanism in sulfate media. Starting from an iron dissolution model, it was possible to suggest two different mechanisms for iron dissolution in 0.5 mol L(-1) H(2)SO(4) containing BTAH or TTAH that involve a complex Fe(II)-inhibitor. (C) 2009 Elsevier B.V. All rights reserved

Filmes finos dielétricos para dispositivos microeletrônicos avançados

Apresentamos mecanismos de formação e de degradação térmica de filmes fi- nos (espessura da ordem de 10 nm) de diferentes dielétricos sobre substrato de silício monocristalino. Tendo em vista a aplicação dessas estruturas em MOSFETs (transistores de efeito de campo metal-óxido-semicondutor), estudamos o consagrado óxido de silício (SiO2), os atuais substitutos oxinitretos de silício (SiOxNy) e o possível substituto futuro óxido de alumínio (Al2O3). Nossos resultados experimentais baseiam-se em técnicas preparativas de substituição isotópica e de caracterização física com feixes de íons (análise com reações nucleares) ou raios- X (espectroscopia de fotoelétrons). Observamos que: (a) átomos de silício não apresentam difusão de longo alcance (além de ~ 2 nm) durante o crescimento de SiO2 por oxidação térmica do silício em O2; (b) nitretação hipertérmica é capaz de produzir filmes finos de oxinitreto de silício com até dez vezes mais nitrogênio que o resultante do processamento térmico usual, sendo que esse nitrogênio tende a se acumular na interface SiOxNy/Si; e (c) átomos de oxigênio, alumínio e silício migram e promovem reações químicas durante o recozimento térmico de estruturas Al2O3/SiO2/Si em presença de O2. Desenvolvemos um modelo de difusão-reação que poderá vir a permitir o estabelecimento de condições ótimas de processamento térmico para filmes finos de Al2O3 sobre silício a serem empregados na fabricação de MOSFETs.

Síntese, caracterização e avaliação catalítica do silicoaluminofosfato SAPO-11 sulfatado

Due to environmental restrictions around the world, clean catalytic technology are of fundamental importance in the petrochemical industry and refineries. Creating the face of this a great interest in replacing the liquid acids for solid acids, so as molecular sieves have been extensively studied in reactions involving the acid catalysis to produce chemical substances with a high potential of quality. Being the activity of the catalysts involved in the reaction attributed to the acid character of them involved for the Lewis and Brönsted acid sites. Based on this context, this study aimed to prepare catalysts acids using a molecular sieve silicoalumino-phosphate (SAPO-11) synthesized in hidrotermical conditions and sulphated with sulphuric acid at different concentrations, using to it the method of controlled impregnating. The samples resulting from this process were characterized by x-ray difratometry (DRX), scanning electron microscopy (SEM), Fourier transform infrared spectroscopy (FT-IR), thermal analysis (TG-DTG) and determination of total acidity (by n-butilamin adsorption). The results show that the synthesis method used was efficient in the formation of AEL structure of SAPO-11 and when being incorporated the sulfate groups in this structure the acidity of the material was increased, pointing out that to very high concentrations of acid there is a trend of decrease the main peaks that form the structure. Finally they were tested catalytictly by the reaction model of conversion of m-xylene which showed favorable results of conversion for this catalyst, showing to be more selective of cracking products than isomerization, as expected, in order that for the o-xylene selectivity there was no positive change when to sulfate a sample of SAPO-11, while for light gases of C1-C4 this selectivity was remarkably observed

Clavulanic acid production processes in a tower bioreactor with immobilised cells

The feasibility of using Streptomyces clavuligerus ATCC 27064 bioparticles supported on alginate gel containing alumina to produce clavulanic acid (CA) was investigated. To this end, effectiveness factors for spherical bioparticles, relating respiration rates of immobilised and free cells, were experimentally determined for various dissolved oxygen (DO) levels and bioparticle radii. Monod kinetics was assumed as representative of the oxygen consuming reaction, while internal oxygen diffusion was considered the limiting step. A comparison was made of the results from a tower bioreactor operating under batch, repeated-batch and continuous conditions with immobilised bioparticles. The theoretical curve of the effectiveness factor for the zero-order reaction model, considering an inert nucleus - the dead core model - was very well fitted to the experimental data. The results of the bioprocess indicated that the batch operation was the most efficient and productive, requiring a do concentration in the reactor above 60% of the saturation value. (C) 2007 Elsevier B.V. All rights reserved.