Modeling of Pressurizer Using APROS and TRACE Thermal Hydraulic Codes

Diplomityön tavoitteena on paineistimen yksityiskohtainen mallintaminen APROS- ja TRACE- termohydrauliikkaohjelmistoja käyttäen. Rakennetut paineistinmallit testattiin vertaamalla laskentatuloksia paineistimen täyttymistä, tyhjentymistä ja ruiskutusta käsittelevistä erilliskokeista saatuun mittausdataan. Tutkimuksen päätavoitteena on APROSin paineistinmallin validoiminen käyttäen vertailuaineistona PACTEL ATWS-koesarjan sopivia paineistinkokeita sekä MIT Pressurizer- ja Neptunus- erilliskokeita. Lisäksi rakennettiin malli Loviisan ydinvoimalaitoksen paineistimesta, jota käytettiin turbiinitrippitransientin simulointiin tarkoituksena selvittää mahdolliset voimalaitoksen ja koelaitteistojen mittakaavaerosta johtuvat vaikutukset APROSin paineistinlaskentaan. Kokeiden simuloinnissa testattiin erilaisia noodituksia ja mallinnusvaihtoehtoja, kuten entalpian ensimmäisen ja toisen kertaluvun diskretisointia, ja APROSin sekä TRACEn antamia tuloksia vertailtiin kattavasti toisiinsa. APROSin paineistinmallin lämmönsiirtokorrelaatioissa havaittiin merkittävä puute ja laskentatuloksiin saatiin huomattava parannus ottamalla käyttöön uusi seinämälauhtumismalli. Työssä tehdyt TRACE-simulaatiot ovat osa United States Nuclear Regulatory Commissionin kansainvälistä CAMP-koodinkehitys-ja validointiohjelmaa.

Regulation of bistable ICEclc transfer in pseudomonas knackmussii B13

L'élément génétique intégratif et conjugatif auto-transférable de 103 kb qui se trouve dans le génome de Pseudomonas knackmussii B13 (ICEc/c) confère la capacité de dégrader le 3-chlorobenzoate et le 2-aminophénol. L'élément ICE c/c peut être transféré par conjugaison de la souche B13 à diverses bêta- et gamma- protéobactéries. Seule une sous-population de 3 à 5% des cellules transfère l'élément, les cellules dites "compétentes pour le transfert". L'acquisition de la compétence pour le transfert est vraisemblablement la conséquence d'une régulation bistable, conduisant une partie des cellules au transfert de l'élément ICE c/c tandis que, dans les autres, l'élément reste quiescent et ne se transfère pas. À ce jour, les mécanismes et les acteurs moléculaires qui régulent l'activation bistable de l'élément sont restés inconnus. Mon travail de doctorat visait à identifier les éléments bistables du régulon de la compétence pour le transfert et d'analyser les fondements moléculaires de la bistabilité de l'élément ICE c/c chez P. knackmussii. Le premier chapitre introduit le thème du transfert génétique horizontal avec un accent particulier sur les éléments intégratifs et conjugatifs (ICE) et ICEcIc. L'état actuel des connaissances sur l'organisation génétique, la régulation, l'intégration et le transfert de différents modèles de ICEs est exposé en détail. En outre, je m'étends sur les phénomènes d'hétérogénéité et de bistabilité phénotyplques, qu'on peut distinguer dans une population isogénique dans des conditions de culture homogènes, et qui sont susceptibles de jouer un rôle dans le transfert de l'élément ICE c/c, dans la mesure où il ne s'active et n'est transférable que dans une très petite sous-population de cellules. Dans le chapitre 2, je présente une analyse globale des régions promotrices minimales des gènes appartenant au régulon de la compétence pour le transfert de l'élément ICE c/c. Nous avons étudié les caractéristiques d'expression des promoteurs et, s'ils s'avéraient bistables, leur activation dans le temps par comparaison avec le mutant lntB13. Pour ce faire, nous avons utilisé des fusions de promoteurs avec des gènes rapporteurs et testé l'expression bistable chez P. knackmussii par microscopie à épifluorescence. Pour six promoteurs présentant une expression bistable, nous avons employé de la microscopie temporelle pour déterminer la chronologie de leur expression par rapport à Pint et PinR. Parmi eux, nous avons identifié deux gènes exprimés précocement et trois gènes exprimés tardivement dans le processus d'acquisition de la compétence de transfert. Dans le chapitre 3, j'expose une analyse d'expression génétique pour l'un des groupes de gènes dont la transcription est la plus élevée dans la région conservée de ICE c/c, les gènes orf81655-orf68241 contenus dans une région de 14 kb. Nous montrons d'abord que cet opéron fait partie du même régulon bistable que intB13 et inrR et analysons les caractéristiques génétiques qui conduisent à une transcription élevée. Nous étudions les fonctions biologiques de ce groupe de gènes par des délétlons ciblées et montrons que certaines d'entre elles empêchent le transfert de l'élément. Nous approfondissons la caractérlsatlon de I'orf8l655 en construisant une fusion transcrlptionnelle avec le gène codant pour la protéine fluorescente verte (egfp) (en utilisant le système minl-Tn5). L'expression de Vorf81655 dans des cellules individuelles est comparée au signal mesuré par hybridation in situ en fluorescence (FISH) sur le ARN messager du gène. En utilisant FISH, des délétlons du promoteur et de l'analyse directe de transcription, nous avons localisé la région promotrice du groupe de gènes. En outre, nous avons utilisé des mutations dirigées pour comprendre la bistabilité de cette région promotrice, caractérisée par une transcription très élevée et une traduction lente de l'ARN messager.  Dans le chapitre 4, nous nous efforçons de comprendre comment la bistabilité est générée au sein du régulon te de l'élément ICE c/c. Pour ce faire, nous avons tenté de reconstituer une expression bistable, dans un hôte qui ne présente pas de bistabilité naturellement, à partir d'éléments génétiques individuels. L'hôte choisi est Pseudomonas putida dans lequel nous avons introduit une copie unique de Pint, PinR ou PaipA fusionnés à la egfp, construits qui permettent d'observer l'apparition de bistabilité. Nous avons ensuite construit différents assemblages de composants génétiques de l'élément ICE c/c, en nous concentrant sur la région parA-inrR. En effet, nous avons pu démontrer qu'une expression bistable apparaît dans P. putida grâce à ces éléments en l'absence de l'élément ICE c/c complet. À noter que la plupart des construits génétiques activent PaipA ou P|,,R, mais qu'un seul recrée la bistabilité de Pint, ce qui suggère que la région parA-inrR permet à la fois d'engendrer la bistabilité et d'opérer la transition entre les promoteurs précoces et les promoteurs tardifs du régulon de la bistabilité. Dans le chapitre 5, nous concluons sur une discussion de la pertinence de nos résultats et sur de futures perspectives de recherche. -- The 103-kb self-transmissible integrative and conjugative element (ICE) of Pseudomonas knackmussii B13 (ICEc/c) confers the capacity to degrade 3- chlorobenzoate and 2-aminophenol. ICEc/c can be conjugated from strain B13 to a variety of Beta- and Gammaproteobacteria. Interestingly, ICE c/c transfer is observed in a subpopulatlon of cells (3-5%) only, the so-called 'transfer competent' cells. The formation of transfer competence (tc) is thought to be the consequence of a 'bistable' decision, which forces those cells to follow the developmental path which leads to ICEc/c transfer, whereas in others ICE c/c remains silent and does not transfer. So far, the mechanisms and molecular partners generating this bistable transfer activation in cells of P. knackmussii B13 remain mostly unidentified. This thesis aimed at understanding the extent of the tc bistability regulon and to dissect the molecular basis of bistabillty formation of ICEc/c in P. knackmussii. The first chapter is a general Introduction on horizontal gene transfer (HGT) with particular emphasis on ICEs and ICE c/c. The emphasis is made on the current knowledge about the HGT gene organization, regulation and specific integration and transfer aspects of the different ICEs models. Furthermore, I focus on the phenomena of phenotypic heterogeneity and bistability (the property of two distinguishable phenotypes existing within an isogenic population under homogeneous conditions), which may play a particular role in ICEc/c behaviour, since ICE activation and transfer only occurs in a very small subpopulation of cells. In Chapter Two, I focus on a global analysis of the different core promoters that might belong to the ICEc/c tc pathway regulon. We studied both expression patterns of ICEc/c promoters and, once being identified as "bistable", their temporal activation compared to that of intB13. In order to do this, we used promoter reporter fusions and tested blstability expression in P. knackmussii using epifluorescence microscopy. For the 6 promoters that showed bistable expression, we used time-lapse microscopy to study the timing of promoter expression in comparison to that of P,,,t or PlnR. We could establish two "early" and 3 "late" phase promoters in the process of transfer competence. In Chapter Three, I focused my attention on analysis of gene expression of one of the most highly transcribed gene clusters in the conserved core region of ICEc/c, a 14-kb gene cluster formed by the genes orf81655-orf68241. First we showed that this operon is part of the same bistability 'regulon' as intB13 and inrR, and analysed the genetic features that lead to high transcription. We studied the potential biological function of this cluster for ICE c/c by making specific gene deletions, showing that some interrupt ICEc/c transfer. We further analysed the orfdl655 promoter by constructing transcriptional egfp fusion reporter strains using the miniTn5 delivery system. Expression of the orf81655 promoter in single cells was compared to signals measured by Fluorescence In Situ Hybridization (FISH) on orfSl655 mRNA. We localized the promoter region of the gene cluster using FISH, promoter deletions, and by direct transcript analysis. We further used site-directed mutagenesis to understand the bistability character of the promoter region and the extremely high transcription but low translation from this mRNA. In Chapter Four, we set out to understand how bistability is generated in the tc pathway of ICEc/c. For this we tried rebuilding bistable expression from ICEc/c individual gene components in a host, which normally does not display bistability. As host we used P. putida without ICEc/c but with a single copy Pint-, PlnR- or PalpA- egfp fusion that enabled us to verify bistability formation. Subsequently, we built different assemblages of ICEc/c gene components, focusing on the parA-inrR region. Indeed, we found that bistable expression can be build from those components in P. putida without ICEc/c. Interestingly, most genetic constructs activated PaipA or PlnR, but only one resulted in bistable activation of PinT. This suggests that the parA-inrR region acts as a bistability "generator", but also as a bistability "relay" from early to late promoters in the tc pathway hierarchy. In the final fifth chapter, we conclude with a discussion of the relevance of the present thesis and the resulting perspectives for future studies.

Radiative Heat Transfer in Boiler Furnaces

Better models are needed for radiative heat transfer in boiler furnaces. If the process is known better, combustion in the furnace can be optimized to produce low emissions. It makes the process to be environmental friendly. Furthermore, if there is a better model of the furnace it can more fully explain what is happening inside the furnace. Using of the model one can quickly and easily analyze how it operates with bio fuels, moist fuels or difficult fuels and improve the operation. Models helps with better estimation of furnace dimensions and result in more accurate understanding of operation. Key component lacking in these models is radiative heat transfer in particle laden gases. If there are no particles than radiative heat transfer can be calculated approximately. There are two problems with current models when used with flow modeling. The first one is a need to account for a particle laden gas and the second one is an absence of a fast algorithm. Fast calculation is needed if radiative heat transfer calculation is done for a large CDF model. Computations slow down if time is required for calculating radiative properties over and over again. This thesis presents a band model for radiative heat transfer in boiler furnaces. Advantage is a quickness of calculation and account of particles in the process.

Fischer-Tropsch Slurry Bubble Column Reactor Modeling

The literature part of the work reviews overall Fischer-Tropsch process, Fischer-Tropsch reactors and catalysts. Fundamentals of Fischer-Tropsch modeling are also presented. The emphasis is on the reactor unit. Comparison of the reactors and the catalysts is carried out to choose the suitable reactor setup for the modeling work. The effects of the operation conditions are also investigated. Slurry bubble column reactor model operating with cobalt catalyst is developed by taking into account the mass transfer of the reacting components (CO and H2) and the consumption of the reactants in the liquid phase. The effect of hydrostatic pressure and the change in total mole flow rate in gas phase are taken into account in calculation of the solubilities. The hydrodynamics, reaction kinetics and product composition are determined according to literature. The cooling system and furthermore the required heat transfer area and number of cooling tubes are also determined. The model is implemented in Matlab software. Commercial scale reactor setup is modeled and the behavior of the model is investigated. The possible inaccuraries are evaluated and the suggestions for the future work are presented. The model is also integrated to Aspen Plus process simulation software, which enables the usage of the model in more extensive Fischer-Tropsch process simulations. Commercial scale reactor of diameter of 7 m and height of 30 m was modeled. The capacity of the reactor was calculated to be about 9 800 barrels/day with CO conversion of 75 %. The behavior of the model was realistic and results were in the right range. The highest uncertainty to model was estimated to be caused by the determination of the kinetic rate.

Modeling the transport phenomena within arterial wall: porous media approach

The transport of macromolecules, such as low-density lipoprotein (LDL), and their accumulation in the layers of the arterial wall play a critical role in the creation and development of atherosclerosis. Atherosclerosis is a disease of large arteries e.g., the aorta, coronary, carotid, and other proximal arteries that involves a distinctive accumulation of LDL and other lipid-bearing materials in the arterial wall. Over time, plaque hardens and narrows the arteries. The flow of oxygen-rich blood to organs and other parts of the body is reduced. This can lead to serious problems, including heart attack, stroke, or even death. It has been proven that the accumulation of macromolecules in the arterial wall depends not only on the ease with which materials enter the wall, but also on the hindrance to the passage of materials out of the wall posed by underlying layers. Therefore, attention was drawn to the fact that the wall structure of large arteries is different than other vessels which are disease-resistant. Atherosclerosis tends to be localized in regions of curvature and branching in arteries where fluid shear stress (shear rate) and other fluid mechanical characteristics deviate from their normal spatial and temporal distribution patterns in straight vessels. On the other hand, the smooth muscle cells (SMCs) residing in the media layer of the arterial wall respond to mechanical stimuli, such as shear stress. Shear stress may affect SMC proliferation and migration from the media layer to intima. This occurs in atherosclerosis and intimal hyperplasia. The study of blood flow and other body fluids and of heat transport through the arterial wall is one of the advanced applications of porous media in recent years. The arterial wall may be modeled in both macroscopic (as a continuous porous medium) and microscopic scales (as a heterogeneous porous medium). In the present study, the governing equations of mass, heat and momentum transport have been solved for different species and interstitial fluid within the arterial wall by means of computational fluid dynamics (CFD). Simulation models are based on the finite element (FE) and finite volume (FV) methods. The wall structure has been modeled by assuming the wall layers as porous media with different properties. In order to study the heat transport through human tissues, the simulations have been carried out for a non-homogeneous model of porous media. The tissue is composed of blood vessels, cells, and an interstitium. The interstitium consists of interstitial fluid and extracellular fibers. Numerical simulations are performed in a two-dimensional (2D) model to realize the effect of the shape and configuration of the discrete phase on the convective and conductive features of heat transfer, e.g. the interstitium of biological tissues. On the other hand, the governing equations of momentum and mass transport have been solved in the heterogeneous porous media model of the media layer, which has a major role in the transport and accumulation of solutes across the arterial wall. The transport of Adenosine 5´-triphosphate (ATP) is simulated across the media layer as a benchmark to observe how SMCs affect on the species mass transport. In addition, the transport of interstitial fluid has been simulated while the deformation of the media layer (due to high blood pressure) and its constituents such as SMCs are also involved in the model. In this context, the effect of pressure variation on shear stress is investigated over SMCs induced by the interstitial flow both in 2D and three-dimensional (3D) geometries for the media layer. The influence of hypertension (high pressure) on the transport of lowdensity lipoprotein (LDL) through deformable arterial wall layers is also studied. This is due to the pressure-driven convective flow across the arterial wall. The intima and media layers are assumed as homogeneous porous media. The results of the present study reveal that ATP concentration over the surface of SMCs and within the bulk of the media layer is significantly dependent on the distribution of cells. Moreover, the shear stress magnitude and distribution over the SMC surface are affected by transmural pressure and the deformation of the media layer of the aorta wall. This work reflects the fact that the second or even subsequent layers of SMCs may bear shear stresses of the same order of magnitude as the first layer does if cells are arranged in an arbitrary manner. This study has brought new insights into the simulation of the arterial wall, as the previous simplifications have been ignored. The configurations of SMCs used here with elliptic cross sections of SMCs closely resemble the physiological conditions of cells. Moreover, the deformation of SMCs with high transmural pressure which follows the media layer compaction has been studied for the first time. On the other hand, results demonstrate that LDL concentration through the intima and media layers changes significantly as wall layers compress with transmural pressure. It was also noticed that the fraction of leaky junctions across the endothelial cells and the area fraction of fenestral pores over the internal elastic lamina affect the LDL distribution dramatically through the thoracic aorta wall. The simulation techniques introduced in this work can also trigger new ideas for simulating porous media involved in any biomedical, biomechanical, chemical, and environmental engineering applications.

Modeling of Gold Cyanidation

The chemistry of gold dissolution in alkaline cyanide solution has continually received attention and new rate equations expressing the gold leaching are still developed. The effect of leaching parameters on gold gold cyanidation is studied in this work in order to optimize the leaching process. A gold leaching model, based on the well-known shrinking-core model, is presented in this work. It is proposed that the reaction takes place at the reacting particle surface which is continuously reduced as the reaction proceeds. The model parameters are estimated by comparing experimental data and simulations. The experimental data used in this work was obtained from Ling et al. (1996) and de Andrade Lima and Hodouin (2005). Two different rate equations, where the unreacted amount of gold is considered in one equation, are investigated. In this work, it is presented that the reaction at the surface is the rate controlling step since there is no internal diffusion limitation. The model considering the effect of non-reacting gold shows that the reaction orders are consistent with the experimental observations reported by Ling et al. (1996) and de Andrade Lima and Hodouin (2005). However, it should be noted that the model obtained in this work is based on assumptions of no side reactions, no solid-liquid mass transfer resistances and no effect from temperature.

A three dimensional study of heat transfer in a supercritical circulating fluidized bed combustor

From the boiler design point of view, it is imperative to know and understand the operation of the boiler. Since comprehensive measurement of a large furnace is impossible, the furnace can be modeled in order to study its behavior and phenomena. This requires the used model to be validated to correspond with the physical furnace behavior. In this thesis, a three dimensional furnace model is validated to match a bituminous coal utilizing, supercritical once-through circulating fluidized bed combustor based on measurement data. The validated model is used for analyzing the furnace heat transfer. Other heat transfer analysis methods are energy balance method based on tube surface temperature measurements and a method based on measured temperature difference between the tube crest and the fin. The latter method was developed in the thesis using Fluent-software. In the theory part, literature is reviewed and the fundamental aspects of circulating fluidized bed are discussed. These aspects are solid particle behavior in fluidization known as hydrodynamics, behavior of fuel and combustion and heat transfer. Fundamental aspects of modeling are also presented.

Leaf-Targeted Ferredoxin-NADP+-Oxidoreductase (FNR): Electron Transfer Properties and Thylakoid Association in Arabidopsis thaliana

Photosynthetic reactions are divided in two parts: light-driven electron transfer reactions and carbon fixation reactions. Electron transfer reactions capture solar energy and split water molecules to form reducing energy (NADPH) and energy-carrying molecules (ATP). These end-products are used for fixation of inorganic carbon dioxide into organic sugar molecules. Ferredoxin-NADP⁺ oxidoreductase (FNR) is an enzyme that acts at the branch point between the electron transfer reactions and reductive metabolism by catalyzing reduction of NADP+ at the last step of the electron transfer chain. In this thesis, two isoforms of FNR from A rabidopsis thaliana, FNR1 and FNR2, were characterized using the reverse genetics approach. The fnr1 and fnr2 mutant plants resembled each other in many respects. Downregulation of photosynthesis protected the single fnr mutant plants from excess formation of reactive oxygen species (ROS), even without significant upregulation of antioxidative mechanisms. Adverse growth conditions, however, resulted in phenotypic differences between fnr1 and fnr2. While fnr2 plants showed downregulation of photosynthetic complexes and upregulation of antioxidative mechanisms under low-temperature growth conditions, fnr1 plants had the wild-type phenotype, indicating that FNR2 may have a specific role in redistribution of electrons under unfavorable conditions. The heterozygotic double mutant (fnr1xfnr2) was severely devoid of chloroplastic FNR, which clearly restricted photosynthesis. The fnr1xfnr2 plants used several photoprotective mechanisms to avoid oxidative stress. In wild-type chloroplasts, both FNR isoforms were found from the stroma, the thylakoid membrane, and the inner envelope membrane. In the absence of the FNR1 isoform, FNR2 was found only in the stroma, suggesting that FNR1 and FNR2 form a dimer, by which FNR1 anchors FNR2 to the thylakoid membrane. Structural modeling predicted formation of an FNR dimer in complex with ferredoxin. In this thesis work, Tic62 was found to be the main protein that binds FNR to the thylakoid membrane, where Tic62 and FNR formed high molecular weight complexes. The formation of such complexes was shown to be regulated by the redox state of the chloroplast. The accumulation of Tic62-FNR complexes in darkness and dissociation of complexes from the membranes in light provide evidence that the complexes may have roles unrelated to photosynthesis. This and the high viability of fnr1 mutant plants lacking thylakoid-bound FNR indicate that the stromal pool of FNR is photosynthetically active.

Kinetic modeling of mechanisms of industrially important organic reactions in gas and liquid phase

This dissertation is based on 5 articles which deal with reaction mechanisms of the following selected industrially important organic reactions: 1. dehydrocyclization of n-butylbenzene to produce naphthalene 2. dehydrocyclization of 1-(p-tolyl)-2-methylbutane (MB) to produce 2,6-dimethylnaphthalene 3. esterification of neopentyl glycol (NPG) with different carboxylic acids to produce monoesters 4. skeletal isomerization of 1-pentene to produce 2-methyl-1-butene and 2-methyl-2-butene The results of initial- and integral-rate experiments of n-butylbenzene dehydrocyclization over selfmade chromia/alumina catalyst were applied when investigating reaction 2. Reaction 2 was performed using commercial chromia/alumina of different acidity, platina on silica and vanadium/calcium/alumina as catalysts. On all catalysts used for the dehydrocyclization, major reactions were fragmentation of MB and 1-(p-tolyl)-2-methylbutenes (MBes), dehydrogenation of MB, double bond transfer, hydrogenation and 1,6-cyclization of MBes. Minor reactions were 1,5-cyclization of MBes and methyl group fragmentation of 1,6- cyclization products. Esterification reactions of NPG were performed using three different carboxylic acids: propionic, isobutyric and 2-ethylhexanoic acid. Commercial heterogeneous gellular (Dowex 50WX2), macroreticular (Amberlyst 15) type resins and homogeneous para-toluene sulfonic acid were used as catalysts. At first NPG reacted with carboxylic acids to form corresponding monoester and water. Then monoester esterified with carboxylic acid to form corresponding diester. In disproportionation reaction two monoester molecules formed NPG and corresponding diester. All these three reactions can attain equilibrium. Concerning esterification, water was removed from the reactor in order to prevent backward reaction. Skeletal isomerization experiments of 1-pentene were performed over HZSM-22 catalyst. Isomerization reactions of three different kind were detected: double bond, cis-trans and skeletal isomerization. Minor side reaction were dimerization and fragmentation. Monomolecular and bimolecular reaction mechanisms for skeletal isomerization explained experimental results almost equally well. Pseudohomogeneous kinetic parameters of reactions 1 and 2 were estimated by usual least squares fitting. Concerning reactions 3 and 4 kinetic parameters were estimated by the leastsquares method, but also the possible cross-correlation and identifiability of parameters were determined using Markov chain Monte Carlo (MCMC) method. Finally using MCMC method, the estimation of model parameters and predictions were performed according to the Bayesian paradigm. According to the fitting results suggested reaction mechanisms explained experimental results rather well. When the possible cross-correlation and identifiability of parameters (Reactions 3 and 4) were determined using MCMC method, the parameters identified well, and no pathological cross-correlation could be seen between any parameter pair.

Modeling study of limestone reactions in fluidized bed conditions

Traditionally limestone has been used for the flue gas desulfurization in fluidized bed combustion. Recently, several studies have been carried out to examine the use of limestone in applications which enable the removal of carbon dioxide from the combustion gases, such as calcium looping technology and oxy-fuel combustion. In these processes interlinked limestone reactions occur but the reaction mechanisms and kinetics are not yet fully understood. To examine these phenomena, analytical and numerical models have been created. In this work, the limestone reactions were studied with aid of one-dimensional numerical particle model. The model describes a single limestone particle in the process as a function of time, the progress of the reactions and the mass and energy transfer in the particle. The model-based results were compared with experimental laboratory scale BFB results. It was observed that by increasing the temperature from 850 °C to 950 °C the calcination was enhanced but the sulfate conversion was no more improved. A higher sulfur dioxide concentration accelerated the sulfation reaction and based on the modeling, the sulfation is first order with respect to SO2. The reaction order of O2 seems to become zero at high oxygen concentrations.

Modeling and simulation of forced-air cooling of strawberries using variable convective coefficient

In the forced-air cooling process of fruits occurs, besides the convective heat transfer, the mass transfer by evaporation. The energy need in the evaporation is taken from fruit that has its temperature lowered. In this study it has been proposed the use of empirical correlations for calculating the convective heat transfer coefficient as a function of surface temperature of the strawberry during the cooling process. The aim of this variation of the convective coefficient is to compensate the effect of evaporation in the heat transfer process. Linear and exponential correlations are tested, both with two adjustable parameters. The simulations are performed using experimental conditions reported in the literature for the cooling of strawberries. The results confirm the suitability of the proposed methodology.

Fluid dynamic modelling and simulation of circulating fluidized bed reactors: importance of the interface turbulence transfer

The main objective of this work is to analyze the importance of the gas-solid interface transfer of the kinetic energy of the turbulent motion on the accuracy of prediction of the fluid dynamic of Circulating Fluidized Bed (CFB) reactors. CFB reactors are used in a variety of industrial applications related to combustion, incineration and catalytic cracking. In this work a two-dimensional fluid dynamic model for gas-particle flow has been used to compute the porosity, the pressure, and the velocity fields of both phases in 2-D axisymmetrical cylindrical co-ordinates. The fluid dynamic model is based on the two fluid model approach in which both phases are considered to be continuous and fully interpenetrating. CFB processes are essentially turbulent. The model of effective stress on each phase is that of a Newtonian fluid, where the effective gas viscosity was calculated from the standard k-epsilon turbulence model and the transport coefficients of the particulate phase were calculated from the kinetic theory of granular flow (KTGF). This work shows that the turbulence transfer between the phases is very important for a better representation of the fluid dynamics of CFB reactors, especially for systems with internal recirculation and high gradients of particle concentration. Two systems with different characteristics were analyzed. The results were compared with experimental data available in the literature. The results were obtained by using a computer code developed by the authors. The finite volume method with collocated grid, the hybrid interpolation scheme, the false time step strategy and SIMPLEC (Semi-Implicit Method for Pressure Linked Equations - Consistent) algorithm were used to obtain the numerical solution.

Substitution-Newton-Raphson method applied to the modeling of a vapour compression refrigeration system using different representations of the thermodynamic properties of R-134a

This paper gives a detailed presentation of the Substitution-Newton-Raphson method, suitable for large sparse non-linear systems. It combines the Successive Substitution method and the Newton-Raphson method in such way as to take the best advantages of both, keeping the convergence features of the Newton-Raphson with the low requirements of memory and time of the Successive Substitution schemes. The large system is solved employing few effective variables, using the greatest possible part of the model equations in substitution fashion to fix the remaining variables, but maintaining the convergence characteristics of the Newton-Raphson. The methodology is exemplified through a simple algebraic system, and applied to a simple thermodynamic, mechanical and heat transfer modeling of a single-stage vapor compression refrigeration system. Three distinct approaches for reproducing the thermodynamic properties of the refrigerant R-134a are compared: the linear interpolation from tabulated data, the use of polynomial fitted curves and the use of functions derived from the Helmholtz free energy.

Modeling of high pressure pretreatment process for gold leaching

Investigation of high pressure pretreatment process for gold leaching is the objective of the present master's thesis. The gold ores and concentrates which cannot be easily treated by leaching process are called "refractory". These types of ores or concentrates often have high content of sulfur and arsenic that renders the precious metal inaccessible to the leaching agents. Since the refractory ores in gold manufacturing industry take a considerable share, the pressure oxidation method (autoclave method) is considered as one of the possible ways to overcome the related problems. Mathematical modeling is the main approach in this thesis which was used for investigation of high pressure oxidation process. For this task, available information from literature concerning this phenomenon, including chemistry, mass transfer and kinetics, reaction conditions, applied apparatus and application, was collected and studied. The modeling part includes investigation of pyrite oxidation kinetics in order to create a descriptive mathematical model. The following major steps are completed: creation of process model by using the available knowledge; estimation of unknown parameters and determination of goodness of the fit; study of the reliability of the model and its parameters.

Modeling of fireside deposit formation in two industrial furnaces

Fireside deposits can be found in many types of utility and industrial furnaces. The deposits in furnaces are problematic because they can reduce heat transfer, block gas paths and cause corrosion. To tackle these problems, it is vital to estimate the influence of deposits on heat transfer, to minimize deposit formation and to optimize deposit removal. It is beneficial to have a good understanding of the mechanisms of fireside deposit formation. Numerical modeling is a powerful tool for investigating the heat transfer in furnaces, and it can provide valuable information for understanding the mechanisms of deposit formation. In addition, a sub-model of deposit formation is generally an essential part of a comprehensive furnace model. This work investigates two specific processes of fireside deposit formation in two industrial furnaces. The first process is the slagging wall found in furnaces with molten deposits running on the wall. A slagging wall model is developed to take into account the two-layer structure of the deposits. With the slagging wall model, the thickness and the surface temperature of the molten deposit layer can be calculated. The slagging wall model is used to predict the surface temperature and the heat transfer to a specific section of a super-heater tube panel with the boundary condition obtained from a Kraft recovery furnace model. The slagging wall model is also incorporated into the computational fluid dynamics (CFD)-based Kraft recovery furnace model and applied on the lower furnace walls. The implementation of the slagging wall model includes a grid simplification scheme. The wall surface temperature calculated with the slagging wall model is used as the heat transfer boundary condition. Simulation of a Kraft recovery furnace is performed, and it is compared with two other cases and measurements. In the two other cases, a uniform wall surface temperature and a wall surface temperature calculated with a char bed burning model are used as the heat transfer boundary conditions. In this particular furnace, the wall surface temperatures from the three cases are similar and are in the correct range of the measurements. Nevertheless, the wall surface temperature profiles with the slagging wall model and the char bed burning model are different because the deposits are represented differently in the two models. In addition, the slagging wall model is proven to be computationally efficient. The second process is deposit formation due to thermophoresis of fine particles to the heat transfer surface. This process is considered in the simulation of a heat recovery boiler of the flash smelting process. In order to determine if the small dust particles stay on the wall, a criterion based on the analysis of forces acting on the particle is applied. Time-dependent simulation of deposit formation in the heat recovery boiler is carried out and the influence of deposits on heat transfer is investigated. The locations prone to deposit formation are also identified in the heat recovery boiler. Modeling of the two processes in the two industrial furnaces enhances the overall understanding of the processes. The sub-models developed in this work can be applied in other similar deposit formation processes with carefully-defined boundary conditions.