Modelling of combustion and sorbent reactions in three-dimensional flow environment of a circulating fluidized bed furnace

This thesis presents a three-dimensional, semi-empirical, steady state model for simulating the combustion, gasification, and formation of emissions in circulating fluidized bed (CFB) processes. In a large-scale CFB furnace, the local feeding of fuel, air, and other input materials, as well as the limited mixing rate of different reactants produce inhomogeneous process conditions. To simulate the real conditions, the furnace should be modelled three-dimensionally or the three-dimensional effects should be taken into account. The only available methods for simulating the large CFB furnaces three-dimensionally are semi-empirical models, which apply a relatively coarse calculation mesh and a combination of fundamental conservation equations, theoretical models and empirical correlations. The number of such models is extremely small. The main objective of this work was to achieve a model which can be applied to calculating industrial scale CFB boilers and which can simulate all the essential sub-phenomena: fluid dynamics, reactions, the attrition of particles, and heat transfer. The core of the work was to develop the model frame and the required sub-models for determining the combustion and sorbent reactions. The objective was reached, and the developed model was successfully used for studying various industrial scale CFB boilers combusting different types of fuel. The model for sorbent reactions, which includes the main reactions for calcitic limestones, was applied for studying the new possible phenomena occurring in the oxygen-fired combustion. The presented combustion and sorbent models and principles can be utilized in other model approaches as well, including other empirical and semi-empirical model approaches, and CFD based simulations. The main achievement is the overall model frame which can be utilized for the further development and testing of new sub-models and theories, and for concentrating the knowledge gathered from the experimental work carried out at bench scale, pilot scale and industrial scale apparatus, and from the computational work performed by other modelling methods.

Computational modeling for irrigated agriculture planning. Part I: general description and linear programming

Linear programming models are effective tools to support initial or periodic planning of agricultural enterprises, requiring, however, technical coefficients that can be determined using computer simulation models. This paper, presented in two parts, deals with the development, application and tests of a methodology and of a computational modeling tool to support planning of irrigated agriculture activities. Part I aimed at the development and application, including sensitivity analysis, of a multiyear linear programming model to optimize the financial return and water use, at farm level for Jaíba irrigation scheme, Minas Gerais State, Brazil, using data on crop irrigation requirement and yield, obtained from previous simulation with MCID model. The linear programming model outputted a crop pattern to which a maximum total net present value of R$ 372,723.00 for the four years period, was obtained. Constraints on monthly water availability, labor, land and production were critical in the optimal solution. In relation to the water use optimization, it was verified that an expressive reductions on the irrigation requirements may be achieved by small reductions on the maximum total net present value.

Computational modeling for irrigated agriculture planning. Part II: risk analysis

Techniques of evaluation of risks coming from inherent uncertainties to the agricultural activity should accompany planning studies. The risk analysis should be carried out by risk simulation using techniques as the Monte Carlo method. This study was carried out to develop a computer program so-called P-RISCO for the application of risky simulations on linear programming models, to apply to a case study, as well to test the results comparatively to the @RISK program. In the risk analysis it was observed that the average of the output variable total net present value, U, was considerably lower than the maximum U value obtained from the linear programming model. It was also verified that the enterprise will be front to expressive risk of shortage of water in the month of April, what doesn't happen for the cropping pattern obtained by the minimization of the irrigation requirement in the months of April in the four years. The scenario analysis indicated that the sale price of the passion fruit crop exercises expressive influence on the financial performance of the enterprise. In the comparative analysis it was verified the equivalence of P-RISCO and @RISK programs in the execution of the risk simulation for the considered scenario.

Computational methods for tactical simulations

Tämä taktiikan tutkimus keskittyy tietokoneavusteisen simuloinnin laskennallisiin menetelmiin, joita voidaan käyttää taktisen tason sotapeleissä. Työn tärkeimmät tuotokset ovat laskennalliset mallit todennäköisyyspohjaisen analyysin mahdollistaviin taktisen tason taistelusimulaattoreihin, joita voidaan käyttää vertailevaan analyysiin joukkue-prikaatitason tarkastelutilanteissa. Laskentamallit keskittyvät vaikuttamiseen. Mallit liittyvät vahingoittavan osuman todennäköisyyteen, jonka perusteella vaikutus joukossa on mallinnettu tilakoneina ja Markovin ketjuina. Edelleen näiden tulokset siirretään tapahtumapuuanalyysiin operaation onnistumisen todennäköisyyden osalta. Pienimmän laskentayksikön mallinnustaso on joukkue- tai ryhmätasolla, jotta laskenta-aika prikaatitason sotapelitarkasteluissa pysyisi riittävän lyhyenä samalla, kun tulokset ovat riittävän tarkkoja suomalaiseen maastoon. Joukkueiden mies- ja asejärjestelmävahvuudet ovat jakaumamuodossa, eivätkä yksittäisiä lukuja. Simuloinnin integroinnissa voidaan käyttää asejärjestelmäkohtaisia predictor corrector –parametreja, mikä mahdollistaa aika-askelta lyhytaikaisempien taistelukentän ilmiöiden mallintamisen. Asemallien pohjana ovat aiemmat tutkimukset ja kenttäkokeet, joista osa kuuluu tähän väitöstutkimukseen. Laskentamallien ohjelmoitavuus ja käytettävyys osana simulointityökalua on osoitettu tekijän johtaman tutkijaryhmän ohjelmoiman ”Sandis”- taistelusimulointiohjelmiston avulla, jota on kehitetty ja käytetty Puolustusvoimien Teknillisessä Tutkimuslaitoksessa. Sandikseen on ohjelmoitu karttakäyttöliittymä ja taistelun kulkua simuloivia laskennallisia malleja. Käyttäjä tai käyttäjäryhmä tekee taktiset päätökset ja syöttää nämä karttakäyttöliittymän avulla simulointiin, jonka tuloksena saadaan kunkin joukkuetason peliyksikön tappioiden jakauma, keskimääräisten tappioiden osalta kunkin asejärjestelmän aiheuttamat tappiot kuhunkin maaliin, ammuskulutus ja radioyhteydet ja niiden tila sekä haavoittuneiden evakuointi-tilanne joukkuetasolta evakuointisairaalaan asti. Tutkimuksen keskeisiä tuloksia (kontribuutio) ovat 1) uusi prikaatitason sotapelitilanteiden laskentamalli, jonka pienin yksikkö on joukkue tai ryhmä; 2) joukon murtumispisteen määritys tappioiden ja haavoittuneiden evakuointiin sitoutuvien taistelijoiden avulla; 3) todennäköisyyspohjaisen riskianalyysin käyttömahdollisuus vertailevassa tutkimuksessa sekä 4) kokeellisesti testatut tulen vaikutusmallit ja 5) toimivat integrointiratkaisut. Työ rajataan maavoimien taistelun joukkuetason todennäköisyysjakaumat luovaan laskentamalliin, kenttälääkinnän malliin ja epäsuoran tulen malliin integrointimenetelmineen sekä niiden antamien tulosten sovellettavuuteen. Ilmasta ja mereltä maahan -asevaikutusta voidaan tarkastella, mutta ei ilma- ja meritaistelua. Menetelmiä soveltavan Sandis -ohjelmiston malleja, käyttötapaa ja ohjelmistotekniikkaa kehitetään edelleen. Merkittäviä jatkotutkimuskohteita mallinnukseen osalta ovat muun muassa kaupunkitaistelu, vaunujen kaksintaistelu ja maaston vaikutus tykistön tuleen sekä materiaalikulutuksen arviointi.

Qualitative Characteristics and Quantitative Measures of Solution's Reliability in Discrete Optimization: Traditional Analytical Approaches, Innovative Computational Methods and Applicability

The purpose of this thesis is twofold. The first and major part is devoted to sensitivity analysis of various discrete optimization problems while the second part addresses methods applied for calculating measures of solution stability and solving multicriteria discrete optimization problems. Despite numerous approaches to stability analysis of discrete optimization problems two major directions can be single out: quantitative and qualitative. Qualitative sensitivity analysis is conducted for multicriteria discrete optimization problems with minisum, minimax and minimin partial criteria. The main results obtained here are necessary and sufficient conditions for different stability types of optimal solutions (or a set of optimal solutions) of the considered problems. Within the framework of quantitative direction various measures of solution stability are investigated. A formula for a quantitative characteristic called stability radius is obtained for the generalized equilibrium situation invariant to changes of game parameters in the case of the H¨older metric. Quality of the problem solution can also be described in terms of robustness analysis. In this work the concepts of accuracy and robustness tolerances are presented for a strategic game with a finite number of players where initial coefficients (costs) of linear payoff functions are subject to perturbations. Investigation of stability radius also aims to devise methods for its calculation. A new metaheuristic approach is derived for calculation of stability radius of an optimal solution to the shortest path problem. The main advantage of the developed method is that it can be potentially applicable for calculating stability radii of NP-hard problems. The last chapter of the thesis focuses on deriving innovative methods based on interactive optimization approach for solving multicriteria combinatorial optimization problems. The key idea of the proposed approach is to utilize a parameterized achievement scalarizing function for solution calculation and to direct interactive procedure by changing weighting coefficients of this function. In order to illustrate the introduced ideas a decision making process is simulated for three objective median location problem. The concepts, models, and ideas collected and analyzed in this thesis create a good and relevant grounds for developing more complicated and integrated models of postoptimal analysis and solving the most computationally challenging problems related to it.

Predictors of fluid intravasation during operative hysteroscopy: a preplanned prospective observational study with 200 cases

PURPOSE: To verify the predictors of intravasation rate during hysteroscopy.METHODS: Prospective observational study (Canadian Task Force classification II-1). All cases (n=200 women; 22 to 86 years old) were treated in an operating room setting. Considering respective bag overfill to calculate water balance, we tested two multiple linear regression models: one for total intravasation (mL) and the other for absorption rate (mL.min-1). The predictors tested (independent variables) were energy (mono/bipolar), tube patency (with/without tubal ligation), hysterometry (cm), age≤50 years, body surface area (m2), surgical complexity (with/without myomectomy) and duration (min).RESULTS: Mean intravasation was significantly higher when myomectomy was performed (442±616 versus 223±332 mL; p<0.01). In the proposed multiple linear regression models for total intravasation (adjusted R2=0.44; p<0.01), the only significant predictors were myomectomy and duration (p<0.01).In the proposed model for intravasation rate (R2=0.39; p<0.01), only myomectomy and hysterometry were significant predictors (p=0.02 and p<0.01, respectively).CONCLUSIONS: Not only myomectomy but also hysterometry were significant predictors of intravasation rate during operative hysteroscopy.

Peritoneal fluid changes in horses subjected to small colon distension

Intestinal devitalization in cases of small colon obstruction may be difficult to detect based only in clinical signs. The purpose was to serially evaluate blood and peritoneal fluid of horses subjected to small colon distension. Seventeen adult horses were allotted in three groups. In the small colon-distended group (DG, n=7) a surgically-implanted latex balloon was inflated to promote intraluminal small colon distension. In the shamoperated group (SG, n=5), the balloon was implanted but not inflated, and no surgery was done in the control group (CG, n=5). Blood and peritoneal fluid were sampled before and after (6 samples with a 30-minute interval) intestinal obstruction for cytological and biochemical analyses. No significant changes in clinical signs occurred within groups or across time during the experimental period. There were no statistical differences among SG and SG groups in hematologic and blood chemistry variables. Although total protein concentration and lactate dehydrogenase (LDH) activity in peritoneal fluid remained most of the time within reference values during the experimental period in all groups, increases from baseline values were detected in SG and DG groups. Such increases occurred earlier, progressively and with greater magnitude in the DG when compared with the SG (P<0.05). Increases from baselines values were also observed in total nucleated cells and neutrophils counts in the DG (P<0.05). In conclusion, distension of the equine small colon induced progressive subtle increases in total protein and LDH concentrations in the peritoneal fluid during the first hours. Serial evaluation of these variables in peritoneal fluid may be useful for early detection of intestinal devitalization in clinical cases of equine small colon obstruction.

Computational chemistry studies of wood-derived lignans

This thesis is based on computational chemistry studies on lignans, focusing on the naturally occurring lignan hydroxymatairesinol (HMR) (Papers I II) and on TADDOL-like conidendrin-based chiral 1,4-diol ligands (LIGNOLs) (Papers III V). A complete quantum chemical conformational analysis on HMR was previously conducted by Dr. Antti Taskinen. In the works reported in this thesis, HMR was further studied by classical molecular dynamics (MD) simulations in aqueous solution including torsional angle analysis, quantum chemical solvation e ect study by the COnductorlike Screening MOdel (COSMO), and hydrogen bond analysis (Paper I), as well as from a catalytic point of view including protonation and deprotonation studies at di erent levels of theory (Paper II). The computational LIGNOL studies in this thesis constitute a multi-level deterministic structural optimization of the following molecules: 1,1-diphenyl (2Ph), two diastereomers of 1,1,4-triphenyl (3PhR, 3PhS), 1,1,4,4-tetraphenyl (4Ph) and 1,1,4,4-tetramethyl (4Met) 1,4-diol (Paper IV) and a conformational solvation study applying MD and COSMO (Paper V). Furthermore, a computational study on hemiketals in connection with problems in the experimental work by Docent Patrik Eklund's group synthesizing the LIGNOLs based on natural products starting from HMR, is shortly described (Paper III).

Enteral fluid therapy through nasogastric tube in rumen cannulated goats

This paper reports the effects of fluid therapy in goats through nasogastric route with an electrolyte solution composed by concentrations of sodium, potassium and chloride similar to goat plasma (140mmol/L of Na+, 4.5mmol/L of K+, 110mmol/L of Cl-). Four Alpine Chamoisee goats, two of them with evident leakage of the rumen cannulas, were used in a crossover experimental design of two periods and two groups. In one group the two goats were submitted to a treatment protocol to induce dehydration before the fluid therapy, whereas the other group was not. Fluid therapy consisted supplying 10mL/kg/h of the electrolyte solution during 8 hours. No signs of discomfort or stress were observed. The dehydration model employed caused a mild dehydration indicated by decrease in feces humidity, body weight and abdominal circumference, and increase in plasma total solids concentration. During fluid therapy globular volume and plasma total solids decreased, whereas % body weight and abdominal circumference increased. No signs of hyperhydration were observed and serum electrolytes (Na+, Cl-, K+) presented no significant alterations in both groups. Fluid therapy proposed in this study was efficient to treat dehydration, even for rumen cannulated animals with evident leakage, and can be administrated safely with no electrolyte imbalance.

Porcine follicular fluid concentration of free insulin-like growth factor-I collected from different diameter ovarian follicles

The study aimed to quantify the concentrations of free IGF-I in serum and fluid of ovarian follicles in pre-pubertal gilts and describe the ovarian morphology by measuring the size of the ovaries and counting the number of surface follicles. Ovaries (n=1,000) from pre-pubertal gilts were obtained immediately after slaughter. A total of 10 samplings were performed, with ovaries obtained from 50 females for each collection. The follicles situated on the surface of each ovary were classified as small (SFs, 2 to 5mm in diameter) or large (LFs 6 to 10mm in diameter) and the follicular fluid was obtained by follicle aspiration. The collection of serum samples was performed after the gilts exsanguination using sterile tubes. From the pool of serum and follicular fluid obtained from 50 females, the concentration of free IGF-I was determined in each sample using an enzyme immunoassay kit (ELISA). The description of ovarian morphometry was performed in 100 ovaries from randomly selected gilts. The larger and smaller lengths of ovaries were measured, and the total number of SFs and LFs present on the surface of each ovary were also counted. The IGF-I concentration was greater (P<0.05) in LFs (170.92±88.29 ng/mL) compared with SFs (67.39±49.90ng/mL) and serum (73.48±34.63ng/mL). The largest and smallest length of the ovaries was 26.0±3.0 and 19.0mm ±2.0mm, respectively. The number of SFs (70.86±25.76) was greater (P<0.01) than LFs (6.54±5.26). The study concluded that LFs present greater levels of IGF-I when compared with SFs and blood, which is related to increased activity of the LFs and its differentiation to ovulation. In addition, ovaries of pre-pubertal gilts have a higher number of SFs compared to LFs. Therefore, our study demonstrated unique data regarding the physiological concentration of free IGF-I in ovarian follicles, that can be used in future research to evaluate the addition of this hormone in the in vitro production media of porcine embryos with the goal to improve the technique efficiency.

Lactated Ringer's solution or 0.9% sodium chloride as fluid therapy in pigeons (Columba livia) submitted to humerus osteosynthesis

The study aimed to compare the effects of intraosseous infusion of lactated Ringer's and 0.9% sodium chloride solutions on the electrolytes and acid-base balance in pigeons submitted to humerus osteosynthesis. Eighteen pigeons were undergoing to isoflurane anesthesia by an avalvular circuit system. They were randomly assigned into two groups (n=9) receiving lactated Ringer's solution (LR) or 0.9% sodium chloride (SC), in a continuous infusion rate of 20mL/kg/h, by using an intraosseous catheter into the tibiotarsus during 60-minute anesthetic procedure. Heart rate (HR), and respiratory rate (RR) were measured every 10 min. Venous blood samples were collected at 0, 30 and 60 minutes to analyze blood pH, PvCO2, HCO3 -, Na+ and K+. Blood gases and electrolytes showed respiratory acidosis in both groups during induction, under physical restraint. This acidosis was evidenced by a decrease of pH since 0 min, associated with a compensatory response, observed by increasing of HCO3 - concentration, at 30 and 60 min. It was not observed any changes on Na+ and K+ serum concentrations. According to the results, there is no reason for choosing one of the two solutions, and it could be concluded that both fluid therapy solutions do not promote any impact on acid-base balance and electrolyte concentrations in pigeons submitted to humerus osteosynthesis.

An alternative finite element formulation for determination of streamlines in two-dimensional problems

It is well known that the numerical solutions of incompressible viscous flows are of great importance in Fluid Dynamics. The graphics output capabilities of their computational codes have revolutionized the communication of ideas to the non-specialist public. In general those codes include, in their hydrodynamic features, the visualization of flow streamlines - essentially a form of contour plot showing the line patterns of the flow - and the magnitudes and orientations of their velocity vectors. However, the standard finite element formulation to compute streamlines suffers from the disadvantage of requiring the determination of boundary integrals, leading to cumbersome implementations at the construction of the finite element code. In this article, we introduce an efficient way - via an alternative variational formulation - to determine the streamlines for fluid flows, which does not need the computation of contour integrals. In order to illustrate the good performance of the alternative formulation proposed, we capture the streamlines of three viscous models: Stokes, Navier-Stokes and Viscoelastic flows.

The Shadow Moiré method using the phase shifting technique and digital image processing: computational implementation and application to the 3D-reconstruction of a buckled plate

The Shadow Moiré fringe patterns are level lines of equal depth generated by interference between a master grid and its shadow projected on the surface. In simplistic approach, the minimum error is about the order of the master grid pitch, that is, always larger than 0,1 mm, resulting in an experimental technique of low precision. The use of a phase shift increases the accuracy of the Shadow Moiré technique. The current work uses the phase shifting method to determine the surfaces three-dimensional shape using isothamic fringe patterns and digital image processing. The current study presents the method and applies it to images obtained by simulation for error evaluation, as well as to a buckled plate, obtaining excellent results. The method hands itself particularly useful to decrease the errors in the interpretation of the Moiré fringes that can adversely affect the calculations of displacements in pieces containing many concave and convex regions in relatively small areas.

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.

Modeling of the Fluid-Structure Interaction in Inelastic Piping Systems

In many engineering applications, compliant piping systems conveying liquids are subjected to inelastic deformations due to severe pressure surges such as plastic tubes in modern water supply transmission lines and metallic pipings in nuclear power plants. In these cases the design of such systems may require an adequate modeling of the interactions between the fluid dynamics and the inelastic structural pipe motions. The reliability of the prediction of fluid-pipe behavior depends mainly on the adequacy of the constitutive equations employed in the analysis. In this paper it is proposed a systematic and general approach to consistently incorporate different kinds of inelastic behaviors of the pipe material in a fluid-structure interaction analysis. The main feature of the constitutive equations considered in this work is that a very simple numerical technique can be used for solving the coupled equations describing the dynamics of the fluid and pipe wall. Numerical examples concerning the analysis of polyethylene and stainless steel pipe networks are presented to illustrate the versatility of the proposed approach.