Offdesign operation of kraft recovery boiler

A model to solve heat and mass balances during the offdesign load calculations was created. These equations are complex and nonlinear. The main new ideas used in the created offdesign model of a kraft recovery boiler are the use of heat flows as torn iteration variables instead of the current practice of using the mass flows, vectorizing equation solving, thus speeding up the process, using non dimensional variables for solving the multiple heat transfer surface problem and using a new procedure for calculating pressure losses. Recovery boiler heat and mass balances are reduced to vector form. It is shown that these vectorized equations can be solved virtually without iteration. The iteration speed is enhanced by the use of the derived method of calculating multiple heat transfer surfaces simultaneously. To achieve this quick convergence the heat flows were used as the torn iteration parameters. A new method to handle pressure loss calculations with linearization was presented. This method enabled less time to be spent calculating pressure losses. The derived vector representation of the steam generator was used to calculate offdesign operation parameters for a 3000 tds/d example recovery boiler. The model was used to study recovery boiler part load operation and the effect of the black liquor dry solids increase on recovery boiler dimensioning. Heat flows to surface elements for part load calculations can be closely approximated with a previously defined exponent function. The exponential method can be used for the prediction of fouling in kraft recovery boilers. For similar furnaces the firing of 80 % dry solids liquor produces lower hearth heat release rate than the 65 % dry solids liquor if we fire at constant steam flow. The furnace outlet temperatures show that capacity increase with firing rate increase produces higher loadings than capacity increase with dry solids increase. The economizers, boiler banks and furnaces can be dimensioned smaller if we increase the black liquor dry solids content. The main problem with increased black liquor dry solids content is the decrease in the heat available to superheat. Whenever possible the furnace exit temperature should be increased by decreasing the furnace height. The increase in the furnace exit temperature is usually opposed because of fear of increased corrosion.

Development of power flows around Baltic countries with new grid links

In the Thesis main focus is done on power flow development paths around the Baltic States as well as on market-based requirements for creation of the common Baltic electricity market. Current market regulations between the countries are presented; barriers for creating competitive common Baltic power market and for electricity trading with third countries are clarified; solutions are offered and corresponding road map is developed. Future power development paths around the Baltic States are analysed. For this purpose the 330 kV transmission grid of Estonia, Latvia and Lithuania is modelled in a power flow tool. Power flow calculations are carried out for winter and summer peak and off-peak load periods in 2020 with different combinations of interconnections. While carrying out power balance experiments several power flow patterns in the Baltic States are revealed. Conclusions are made about security of supply, grid congestion and transmission capacity availability for different scenarios.

Spottier Targets Are Less Attractive to Tabanid Flies: On the Tabanid-Repellency of Spotty Fur Patterns

During blood-sucking, female members of the family Tabanidae transmit pathogens of serious diseases and annoy their host animals so strongly that they cannot graze, thus the health of the hosts is drastically reduced. Consequently, a tabanid-resistant coat with appropriate brightness, colour and pattern is advantageous for the host. Spotty coats are widespread among mammals, especially in cattle (Bos primigenius). In field experiments we studied the influence of the size and number of spots on the attractiveness of test surfaces to tabanids that are attracted to linearly polarized light. We measured the reflection-polarization characteristics of living cattle, spotty cattle coats and the used test surfaces. We show here that the smaller and the more numerous the spots, the less attractive the target (host) is to tabanids. We demonstrate that the attractiveness of spotty patterns to tabanids is also reduced if the target exhibits spottiness only in the angle of polarization pattern, while being homogeneous grey with a constant high degree of polarization. Tabanid flies respond strongly to linearly polarized light, and we show that bright and dark parts of cattle coats reflect light with different degrees and angles of polarization that in combination with dark spots on a bright coat surface disrupt the attractiveness to tabanids. This could be one of the possible evolutionary benefits that explains why spotty coat patterns are so widespread in mammals, especially in ungulates, many species of which are tabanid hosts

Computational and Experimental Investigation of the Grooved Roll in Paper Machine Environment

In the paper machine, it is not a desired feature for the boundary layer flows in the fabric and the roll surfaces to travel into the closing nips, creating overpressure. In this thesis, the aerodynamic behavior of the grooved roll and smooth rolls is compared in order to understand the nip flow phenomena, which is the main reason why vacuum and grooved roll constructions are designed. A common method to remove the boundary layer flow from the closing nip is to use the vacuum roll construction. The downside of the use of vacuum rolls is high operational costs due to pressure losses in the vacuum roll shell. The deep grooved roll has the same goal, to create a pressure difference over the paper web and keep the paper attached to the roll or fabric surface in the drying pocket of the paper machine. A literature review revealed that the aerodynamic functionality of the grooved roll is not very well known. In this thesis, the aerodynamic functionality of the grooved roll in interaction with a permeable or impermeable wall is studied by varying the groove properties. Computational fluid dynamics simulations are utilized as the research tool. The simulations have been performed with commercial fluid dynamics software, ANSYS Fluent. Simulation results made with 3- and 2-dimensional fluid dynamics models are compared to laboratory scale measurements. The measurements have been made with a grooved roll simulator designed for the research. The variables in the comparison are the paper or fabric wrap angle, surface velocities, groove geometry and wall permeability. Present-day computational and modeling resources limit grooved roll fluid dynamics simulations in the paper machine scale. Based on the analysis of the aerodynamic functionality of the grooved roll, a grooved roll simulation tool is proposed. The smooth roll simulations show that the closing nip pressure does not depend on the length of boundary layer development. The surface velocity increase affects the pressure distribution in the closing and opening nips. The 3D grooved roll model reveals the aerodynamic functionality of the grooved roll. With the optimal groove size it is possible to avoid closing nip overpressure and keep the web attached to the fabric surface in the area of the wrap angle. The groove flow friction and minor losses play a different role when the wrap angle is changed. The proposed 2D grooved roll simulation tool is able to replicate the grooved aerodynamic behavior with reasonable accuracy. A small wrap angle predicts the pressure distribution correctly with the chosen approach for calculating the groove friction losses. With a large wrap angle, the groove friction loss shows too large pressure gradients, and the way of calculating the air flow friction losses in the groove has to be reconsidered. The aerodynamic functionality of the grooved roll is based on minor and viscous losses in the closing and opening nips as well as in the grooves. The proposed 2D grooved roll model is a simplification in order to reduce computational and modeling efforts. The simulation tool makes it possible to simulate complex paper machine constructions in the paper machine scale. In order to use the grooved roll as a replacement for the vacuum roll, the grooved roll properties have to be considered on the basis of the web handling application.

Impacts of climate change on water resources in the Mediterranean basin : a case study in Catalonia, Spain

Most climate change projections show important decreases in water availability in the Mediterranean region by the end of this century. We assess those main climate change impacts on water resources in three medium-sized catchments with varying climatic conditions in north-eastern Spain. A combination of hydrological modelling and climate projections with B1 and A2 IPCC emission scenarios is performed to infer future stream flows. The largest reduction (22-48% for 2076-2100) of stream flows is expected in the headwaters of the two wettest catchments, while lower decreases (22-32% for 2076-2100) are expected in the drier one. In all three catchments, autumn and summer are the seasons with the most notable projected decreases in stream flow, 50% and 34%, respectively (2076-2100). Thus, ecological flows might be noticeably impacted by climate change in the catchments, especially in the headwaters of those wet catchments.

Immobilization of laccase on hybrid layered double hydroxide

Crystals of Mg/Al layered double hydroxide were synthesized by alkaline precipitation and treated in an aqueous solution of glutamic acid. The glutamate ions were not intercalated into the interlayer space, but were detected in the material by Fourier transform infrared spectroscopy, suggesting that only the external surfaces of crystals were modified with glutamate ions. The resulting hybrid material was tested as a support for immobilization of the enzyme laccase (Myceliophthora thermophila). The immobilized enzyme preparation was characterized by electronic paramagnetic resonance spectroscopy and by assays of catalytic activity. The activity of the immobilized laccase was 97% of the activity in the free enzyme. Layered double hydroxide is a suitable support for use in remediation of soil studies.

APPLICATION OF PARAFAC AND OPLS-DA ANALYSES ON HPLC FINGERPRINTS FOR THE CHARACTERIZATION OF Hibiscus sabdariffaCALYXES

The calyxes of Hibiscus sabdariffa are used in traditional medicine around the world. However, quality assurance protocols and chemical variability have not been previously analyzed. In the present study, chemical characterization of a set of samples of H. sabdariffa calyxes commercialized in Colombia was accomplished with the aim to explore the chemical variability among them. Chemometrics-based analyses on the data obtained from the HPLC-UV-DAD-derived profiles were then performed. Thus, the pre-processed single-wavelength data were subjected to principal component analysis (PCA). The PCA-derived results evidenced different groups which were well-correlated to the corresponding total phenolic and total anthocyanin contents. Multi-wavelength chromatographic (HPLC-UV-DAD surfaces) data were additionally examined via parallel factor analysis (PARAFAC) as data reduction method and the obtained loadings were subsequently submitted to PCA and orthogonal partial least squares discriminant analysis (OPLS-DA). Results were thus consistent with those from single-wavelength data. PCA loadings were employed to determine those chemical components responsible for the data variance and OPLS-DA model, constructed from PARAFAC loadings, and indicated differentiation according total anthocyanin contents among samples. The present chemometric analysis therefore demonstrated to be an excellent tool for differentiation of H. sabdariffacalyxes according to their chemical composition.

CFD Simulation of Two-phase and Three-phase Flows in Internal-loop Airlift

Airlift reactors are pneumatically agitated reactors that have been widely used in chemical, petrochemical, and bioprocess industries, such as fermentation and wastewater treatment. Computational Fluid Dynamics (CFD) has become more popular approach for design, scale-up and performance evaluation of such reactors. In the present work numerical simulations for internal-loop airlift reactors were performed using the transient Eulerian model with CFD package, ANSYS Fluent 12.1. The turbulence in the liquid phase is described using κ- ε the model. Global hydrodynamic parameters like gas holdup, gas velocity and liquid velocity have been investigated for a range of superficial gas velocities, both with 2D and 3D simulations. Moreover, the study of geometry and scale influence on the reactor have been considered. The results suggest that both, geometry and scale have significant effects on the hydrodynamic parameters, which may have substantial effects on the reactor performance. Grid refinement and time-step size effect have been discussed. Numerical calculations with gas-liquid-solid three-phase flow system have been carried out to investigate the effect of solid loading, solid particle size and solid density on the hydrodynamic characteristics of internal loop airlift reactor with different superficial gas velocities. It was observed that averaged gas holdup is significantly decreased with increasing slurry concentration. Simulations show that the riser gas holdup decreases with increase in solid particle diameter. In addition, it was found that the averaged solid holdup increases in the riser section with the increase of solid density. These produced results reveal that CFD have excellent potential to simulate two-phase and three-phase flow system.

First principles modeling of metallic alloys and alloy surfaces

By alloying metals with other materials, one can modify the metal’s characteristics or compose an alloy which has certain desired characteristics that no pure metal has. The field is vast and complex, and phenomena that govern the behaviour of alloys are numerous. Theories cannot penetrate such complexity, and the scope of experiments is also limited. This is why the relatively new field of ab initio computational methods has much to give to this field. With these methods, one can extend the understanding given by theories, predict how some systems might behave, and be able to obtain information that is not there to see in physical experiments. This thesis pursues to contribute to the collective knowledge of this field in the light of two cases. The first part examines the oxidation of Ag/Cu, namely, the adsorption dynamics and oxygen induced segregation of the surface. Our results demonstrate that the presence of Ag on the Cu(100) surface layer strongly inhibits dissociative adsorption. Our results also confirmed that surface reconstruction does happen, as experiments predicted. Our studies indicate that 0.25 ML of oxygen is enough for Ag to diffuse towards the bulk, under the copper oxide layer. The other part elucidates the complex interplay of various energy and entropy contributions to the phase stability of paramagnetic duplex steel alloys. We were able to produce a phase stability map from first principles, and it agrees with experiments rather well. Our results also show that entropy contributions play a very important role on defining the phase stability. This is, to the author’s knowledge, the first ab initio study upon this subject.

Effect of hydrodynamics on modelling, monitoring and control of crystallization

Crystallization is a purification method used to obtain crystalline product of a certain crystal size. It is one of the oldest industrial unit processes and commonly used in modern industry due to its good purification capability from rather impure solutions with reasonably low energy consumption. However, the process is extremely challenging to model and control because it involves inhomogeneous mixing and many simultaneous phenomena such as nucleation, crystal growth and agglomeration. All these phenomena are dependent on supersaturation, i.e. the difference between actual liquid phase concentration and solubility. Homogeneous mass and heat transfer in the crystallizer would greatly simplify modelling and control of crystallization processes, such conditions are, however, not the reality, especially in industrial scale processes. Consequently, the hydrodynamics of crystallizers, i.e. the combination of mixing, feed and product removal flows, and recycling of the suspension, needs to be thoroughly investigated. Understanding of hydrodynamics is important in crystallization, especially inlargerscale equipment where uniform flow conditions are difficult to attain. It is also important to understand different size scales of mixing; micro-, meso- and macromixing. Fast processes, like nucleation and chemical reactions, are typically highly dependent on micro- and mesomixing but macromixing, which equalizes the concentrations of all the species within the entire crystallizer, cannot be disregarded. This study investigates the influence of hydrodynamics on crystallization processes. Modelling of crystallizers with the mixed suspension mixed product removal (MSMPR) theory (ideal mixing), computational fluid dynamics (CFD), and a compartmental multiblock model is compared. The importance of proper verification of CFD and multiblock models is demonstrated. In addition, the influence of different hydrodynamic conditions on reactive crystallization process control is studied. Finally, the effect of extreme local supersaturation is studied using power ultrasound to initiate nucleation. The present work shows that mixing and chemical feeding conditions clearly affect induction time and cluster formation, nucleation, growth kinetics, and agglomeration. Consequently, the properties of crystalline end products, e.g. crystal size and crystal habit, can be influenced by management of mixing and feeding conditions. Impurities may have varying impacts on crystallization processes. As an example, manganese ions were shown to replace magnesium ions in the crystal lattice of magnesium sulphate heptahydrate, increasing the crystal growth rate significantly, whereas sodium ions showed no interaction at all. Modelling of continuous crystallization based on MSMPR theory showed that the model is feasible in a small laboratoryscale crystallizer, whereas in larger pilot- and industrial-scale crystallizers hydrodynamic effects should be taken into account. For that reason, CFD and multiblock modelling are shown to be effective tools for modelling crystallization with inhomogeneous mixing. The present work shows also that selection of the measurement point, or points in the case of multiprobe systems, is crucial when process analytical technology (PAT) is used to control larger scale crystallization. The thesis concludes by describing how control of local supersaturation by highly localized ultrasound was successfully applied to induce nucleation and to control polymorphism in reactive crystallization of L-glutamic acid.

Development of railway logistics’ information flows

Objective of this thesis was to develop the exchange of information and reduce the manual work done in the supply chain. In addition, the possibility to introduce electronic information exchange was studied between suppliers and Borealis. The aim was to create an accurate picture of Borealis’ current information flows and create from the basis of it short- and long-term improvement and development proposals. In this study the company's received and send information flows were mapped by interviewing persons who were responsible for the railroad imports and by examining documents that are used in the exchange of information. The data content of the information flows were prioritized and only the most important information contents were used for further development. Literature data was acquired concerning knowledge of electronic data interchange and information management to support the decisions and proposals. Long-term development proposals were compared with each other and the best one of them was recommended for further study. The final target of the proposal is to be able to receive electronic data and create an own database where to the information is stored and where from it is possible to follow up the rail tank cars and where from the needed reports can be retrieved.

MODELING AND SIMULATION OF RELIEF INFLUENCE ON EUCALYPTUS FORESTS: INTERACTION BETWEEN SOLAR IRRADIANCE AND PRODUCTIVITY

ABSTRACT This study aimed to verify the differences in radiation intensity as a function of distinct relief exposure surfaces and to quantify these effects on the leaf area index (LAI) and other variables expressing eucalyptus forest productivity for simulations in a process-based growth model. The study was carried out at two contrasting edaphoclimatic locations in the Rio Doce basin in Minas Gerais, Brazil. Two stands with 32-year-old plantations were used, allocating fixed plots in locations with northern and southern exposure surfaces. The meteorological data were obtained from two automated weather stations located near the study sites. Solar radiation was corrected for terrain inclination and exposure surfaces, as it is measured based on the plane, perpendicularly to the vertical location. The LAI values collected in the field were used. For the comparative simulations in productivity variation, the mechanistic 3PG model was used, considering the relief exposure surfaces. It was verified that during most of the year, the southern surfaces showed lower availability of incident solar radiation, resulting in up to 66% losses, compared to the same surface considered plane, probably related to its geographical location and higher declivity. Higher values were obtained for the plantings located on the northern surface for the variables LAI, volume and mean annual wood increase, with this tendency being repeated in the 3PG model simulations.

Annual evolution of global, direct and diffuse radiation and fractions in tilted surfaces

It was evaluated the annual evolution of global, direct and diffuse components of incident solar radiation on tilted surfaces to 12.85, 22.85 and 32.85º, facing north, in Botucatu, state of São Paulo, Brazil. The radiometric fractions were obtained for each component of the radiation in the aforementioned surfaces, through the ratio with the global and top of the atmosphere radiations. Seasonality was evaluated based on monthly averages of daily values. The measures occurred between 04/1998 and 07/2001 at 22.85º; 08/2001 and 02/2003 at 12.85º; and from 03/2003 to 12/2007 at 32.85º, with concomitant measures in the horizontal surface (reference). The levels of global and direct radiation on tilted surfaces were lower in summer and higher in the equinoxes when compared with the horizontal. The diffuse radiation on tilted surfaces was lower in most months, with losses of up to 65%. A trend of increasing differences occurred between horizontal and tilted surfaces with the increase of the angle in all the components and fractions of incident radiation. The annual evolution of rainfall and cloud cover ratio directly affected the atmospheric transmissivity of direct and diffuse components in the region.

Influence of urbanization on land surface temperature in Recife city

Urbanization has caused significant environmental impacts, replacing natural surfaces by buildings, decreasing green vegetated areas, soil sealing and atmospheric pollution which contribute to increase the land surface temperature in such areas. Thus, this study aimed to analyze the influence of urbanization on land surface temperature (Ts) in Recife city - Pernambuco (PE), in Brazil, using the Thematic Mapper (TM) sensor images from Landsat 5 satellite. To perform the study, images of August 4, 1998 and September 6, 2010 were obtained and processed to generate Ts thematic maps of Recife-PE and of two districts of this city (Curado and Casa Amarela), in order to analyze the transformation dynamics that has occurred in the area. Through the profile produced for the study area, a spatial and temporal increase of the Ts surface was noticeable in the suburb-downtown direction: 6°C of difference between these areas. The Casa Amarela district, with high urban concentration, presented the highest Ts values observed (>27°C).

Outlining precision boundaries among areas with different variability standards using magnetic susceptibility and geomorphic surfaces

There is an increasing demand for detailed maps that represent in a simplified way the knowledge of the variability of a particular area or region maps. The objective was to outline precision boundaries among areas with different accuracy variability standards using magnetic susceptibility and geomorphic surfaces. The study was conducted in an area of 110 ha, which identified three compartment landscapes based on the geomorphic surfaces model. To determinate pH, organic matter, phosphorus, potassium and magnesium, the total sand and clay, 514 soil samples were collected at depths of 0-0.20 m and 0.60-0.80 m. The sum of base, cationic exchange capacity and base saturation were calculated and the magnetic susceptibility was evaluated in the laboratory using a system based on a balance of analytical precision method. Geomorphic surfaces identification allowed setting specific management areas (locations with maximum homogeneity of soil attributes). The map of spatial variability of magnetic susceptibility can be used to validate the precise boundaries among geomorphic surfaces identified in the field and infer the variability of clay content and soil base saturation.