Fiber flows in global paper and pulp industry. A panel data approach.

The purpose of this thesis is to study factors that explain the bilateral fiber trade flows. This is done by analyzing bilateral trade flows during 1990-2006. It will be studied also, whether there are differences between fiber types. This thesis uses a gravity model approach to study the trade flows. Gravity model is mostly used to study the aggregate data between trading countries. In this thesis the gravity model is applied to single fibers. This model is then applied to panel data set. Results from the regression show clearly that there are benefits in studying different fibers in separate. The effects differ considerably from each other. Furthermore, this thesis speaks for the existence of Linder’s effect in certain fiber types.

Identification of new transformation products during enzymatic treatment of tetracycline and erythromycin antibiotics at laboratory scale by an on-line turbulent flow liquid-chromatography coupled to a high resolution mass spectrometer LTQ-Orbitrap

This work describes the formation of transformation products (TPs) by the enzymatic degradation at laboratory scale of two highly consumed antibiotics: tetracycline (Tc) and erythromycin (ERY). The analysis of the samples was carried out by a fast and simple method based on the novel configuration of the on-line turbulent flow system coupled to a hybrid linear ion trap – high resolution mass spectrometer. The method was optimized and validated for the complete analysis of ERY, Tc and their transformation products within 10 min without any other sample manipulation. Furthermore, the applicability of the on-line procedure was evaluated for 25 additional antibiotics, covering a wide range of chemical classes in different environmental waters with satisfactory quality parameters. Degradation rates obtained for Tc by laccase enzyme and ERY by EreB esterase enzyme without the presence of mediators were ∼78% and ∼50%, respectively. Concerning the identification of TPs, three suspected compounds for Tc and five of ERY have been proposed. In the case of Tc, the tentative molecular formulas with errors mass within 2 ppm have been based on the hypothesis of dehydroxylation, (bi)demethylation and oxidation of the rings A and C as major reactions. In contrast, the major TP detected for ERY has been identified as the “dehydration ERY-A”, with the same molecular formula of its parent compound. In addition, the evaluation of the antibiotic activity of the samples along the enzymatic treatments showed a decrease around 100% in both cases

Process intensification: From optimised flow patterns to microprocess technology

This thesis is focused on process intensification. Several significant problems and applications of this theme are covered. Process intensification is nowadays one of the most popular trends in chemical engineering and attempts have been made to develop a general, systematic methodology for intensification. This seems, however, to be very difficult, because intensified processes are often based on creativity and novel ideas. Monolith reactors and microreactors are successful examples of process intensification. They are usually multichannel devices in which a proper feed technique is important for creating even fluid distribution into the channels. Two different feed techniques were tested for monoliths. In the first technique a shower method was implemented by means of perforated plates. The second technique was a dispersion method using static mixers. Both techniques offered stable operation and uniform fluid distribution. The dispersion method enabled a wider operational range in terms of liquid superficial velocity. Using dispersion method, a volumetric gas-liquid mass transfer coefficient of 2 s-1 was reached. Flow patterns play a significant role in terms of the mixing performance of micromixers. Although the geometry of a T-mixer is simple, channel configurations and dimensions had a clear effect on mixing efficiency. The flow in the microchannel was laminar, but the formation of vortices promoted mixing in micro T-mixers. The generation of vortices was dependent on the channel dimensions, configurations and flow rate. Microreactors offer a high ratio of surface area to volume. Surface forces and interactions between fluids and surfaces are, therefore, often dominant factors. In certain cases, the interactions can be effectively utilised. Different wetting properties of solid materials (PTFE and stainless steel) were applied in the separation of immiscible liquid phases. A micro-scale plate coalescer with hydrophilic and hydrophobic surfaces was used for the continuous separation of organic and aqueous phases. Complete phase separation occurred in less than 20 seconds, whereas the separation time by settling exceeded 30 min. Fluid flows can be also intensified in suitable conditions. By adding certain additives into turbulent fluid flow, it was possible to reduce friction (drag) by 40 %. Drag reduction decreases frictional pressure drop in pipelines which leads to remarkable energy savings and decreases the size or number of pumping facilities required, e.g., in oil transport pipes. Process intensification enables operation often under more optimal conditions. The consequent cost savings from reduced use of raw materials and reduced waste lead to greater economic benefits in processing.

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.

LOW COST ANALYZER FOR THE DETERMINATION OF PHOSPHORUS BASED ON OPEN-SOURCE HARDWARE AND PULSED FLOWS

The need for automated analyzers for industrial and environmental samples has triggered the research for new and cost-effective strategies of automation and control of analytical systems. The widespread availability of open-source hardware together with novel analytical methods based on pulsed flows have opened the possibility of implementing standalone automated analytical systems at low cost. Among the areas that can benefit from this approach are the analysis of industrial products and effluents and environmental analysis. In this work, a multi-pumping flow system is proposed for the determination of phosphorus in effluents and polluted water samples. The system employs photometric detection based on the formation of molybdovanadophosphoric acid, and the fluidic circuit is built using three solenoid micropumps. The detection is implemented with a low cost LED-photodiode photometric detection system and the whole system is controlled by an open-source Arduino Uno microcontroller board. The optimization of the timing to ensure the color development and the pumping cycle is discussed for the proposed implementation. Experimental results to evaluate the system behavior are presented verifying a linear relationship between the relative absorbance and the phosphorus concentrations for levels as high as 50 mg L-1.

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.

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.

CFD Modelling of Direct Contact Condensation in Suppression Pools by Applying Condensation Models of Separated Flow

The condensation rate has to be high in the safety pressure suppression pool systems of Boiling Water Reactors (BWR) in order to fulfill their safety function. The phenomena due to such a high direct contact condensation (DCC) rate turn out to be very challenging to be analysed either with experiments or numerical simulations. In this thesis, the suppression pool experiments carried out in the POOLEX facility of Lappeenranta University of Technology were simulated. Two different condensation modes were modelled by using the 2-phase CFD codes NEPTUNE CFD and TransAT. The DCC models applied were the typical ones to be used for separated flows in channels, and their applicability to the rapidly condensing flow in the condensation pool context had not been tested earlier. A low Reynolds number case was the first to be simulated. The POOLEX experiment STB-31 was operated near the conditions between the ’quasi-steady oscillatory interface condensation’ mode and the ’condensation within the blowdown pipe’ mode. The condensation models of Lakehal et al. and Coste & Lavi´eville predicted the condensation rate quite accurately, while the other tested ones overestimated it. It was possible to get the direct phase change solution to settle near to the measured values, but a very high resolution of calculation grid was needed. Secondly, a high Reynolds number case corresponding to the ’chugging’ mode was simulated. The POOLEX experiment STB-28 was chosen, because various standard and highspeed video samples of bubbles were recorded during it. In order to extract numerical information from the video material, a pattern recognition procedure was programmed. The bubble size distributions and the frequencies of chugging were calculated with this procedure. With the statistical data of the bubble sizes and temporal data of the bubble/jet appearance, it was possible to compare the condensation rates between the experiment and the CFD simulations. In the chugging simulations, a spherically curvilinear calculation grid at the blowdown pipe exit improved the convergence and decreased the required cell count. The compressible flow solver with complete steam-tables was beneficial for the numerical success of the simulations. The Hughes-Duffey model and, to some extent, the Coste & Lavi´eville model produced realistic chugging behavior. The initial level of the steam/water interface was an important factor to determine the initiation of the chugging. If the interface was initialized with a water level high enough inside the blowdown pipe, the vigorous penetration of a water plug into the pool created a turbulent wake which invoked the chugging that was self-sustaining. A 3D simulation with a suitable DCC model produced qualitatively very realistic shapes of the chugging bubbles and jets. The comparative FFT analysis of the bubble size data and the pool bottom pressure data gave useful information to distinguish the eigenmodes of chugging, bubbling, and pool structure oscillations.

Material embodiment and energy flows as efficiency indicators of soybean (Glycine max) production in Brazil

As the requirement for agriculture to be environmentally suitable there is a necessity to adopt indicators and methodologies approaching sustainability. In Brazil, biodiesel addition into diesel is mandatory and soybean oil is its main source. The material embodiment determines the convergence of inputs into the crop. Moreover, the material flows are necessary for any environmental analysis. This study evaluated distinct production scenarios, and also conventional versus GMO crops, through the material embodiment and energy analysis. GMO crops demanded less indirectly applied inputs. The energy balance showed linearity with yield, whereas for EROI, the increases in input and yield were not affected.

Analysis and validation of space averaged drag model for numerical simulations of gas-solid flows in fluidized beds

This thesis presents an approach for formulating and validating a space averaged drag model for coarse mesh simulations of gas-solid flows in fluidized beds using the two-fluid model. Proper modeling for fluid dynamics is central in understanding any industrial multiphase flow. The gas-solid flows in fluidized beds are heterogeneous and usually simulated with the Eulerian description of phases. Such a description requires the usage of fine meshes and small time steps for the proper prediction of its hydrodynamics. Such constraint on the mesh and time step size results in a large number of control volumes and long computational times which are unaffordable for simulations of large scale fluidized beds. If proper closure models are not included, coarse mesh simulations for fluidized beds do not give reasonable results. The coarse mesh simulation fails to resolve the mesoscale structures and results in uniform solids concentration profiles. For a circulating fluidized bed riser, such predicted profiles result in a higher drag force between the gas and solid phase and also overestimated solids mass flux at the outlet. Thus, there is a need to formulate the closure correlations which can accurately predict the hydrodynamics using coarse meshes. This thesis uses the space averaging modeling approach in the formulation of closure models for coarse mesh simulations of the gas-solid flow in fluidized beds using Geldart group B particles. In the analysis of formulating the closure correlation for space averaged drag model, the main parameters for the modeling were found to be the averaging size, solid volume fraction, and distance from the wall. The closure model for the gas-solid drag force was formulated and validated for coarse mesh simulations of the riser, which showed the verification of this modeling approach. Coarse mesh simulations using the corrected drag model resulted in lowered values of solids mass flux. Such an approach is a promising tool in the formulation of appropriate closure models which can be used in coarse mesh simulations of large scale fluidized beds.

Model-based estimation of liquid flows in the blast furnace hearth and taphole

This thesis presents a set of methods and models for estimation of iron and slag flows in the blast furnace hearth and taphole. The main focus was put on predicting taphole flow patterns and estimating the effects of various taphole conditions on the drainage behavior of the blast furnace hearth. All models were based on a general understanding of the typical tap cycle of an industrial blast furnace. Some of the models were evaluated on short-term process data from the reference furnace. A computational fluid dynamics (CFD) model was built and applied to simulate the complicated hearth flows and thus to predict the regions of the hearth exerted to erosion under various operating conditions. Key boundary variables of the CFD model were provided by a simplified drainage model based on the first principles. By examining the evolutions of liquid outflow rates measured from the furnace studied, the drainage model was improved to include the effects of taphole diameter and length. The estimated slag delays showed good agreement with the observed ones. The liquid flows in the taphole were further studied using two different models and the results of both models indicated that it is more likely that separated flow of iron and slag occurs in the taphole when the liquid outflow rates are comparable during tapping. The drainage process was simulated with an integrated model based on an overall balance analysis: The high in-furnace overpressure can compensate for the resistances induced by the liquid flows in the hearth and through the taphole. Finally, a recently developed multiphase CFD model including interfacial forces between immiscible liquids was developed and both the actual iron-slag system and a water-oil system in laboratory scale were simulated. The model was demonstrated to be a useful tool for simulating hearth flows for gaining understanding of the complex phenomena in the drainage of the blast furnace.

Effect of different surface roughnesses on a turbulent boundary layer

The classical treatment of rough wall turbulent boundary layers consists in determining the effect the roughness has on the mean velocity profile. This effect is usually described in terms of the roughness function delta U+. The general implication is that different roughness geometries with the same delta U+ will have similar turbulence characteristics, at least at a sufficient distance from the roughness elements. Measurements over two different surface geometries (a mesh roughness and spanwise circular rods regularly spaced in the streamwise direction) with nominally the same delta U+ indicate significant differences in the Reynolds stresses, especially those involving the wall-normal velocity fluctuation, over the outer region. The differences are such that the Reynolds stress anisotropy is smaller over the mesh roughness than the rod roughness. The Reynolds stress anisotropy is largest for a smooth wall. The small-scale anisotropy and interniittency exhibit much smaller differences when the Taylor microscale Reynolds number and the Kolmogorov-normalized mean shear are nominally the same. There is nonetheless evidence that the small-scale structure over the three-dimensional mesh roughness conforms more closely with isotropy than that over the rod-roughened and smooth walls.

Modeling of turbulent flow through radial diffuser

The work considers the modeling of turbulent flow in radial diffuser with axial feeding. Due to its claimed capability to predict flow including features such as separation, curvature and adverse pressure gradient, the RNG k-epsilon model of Orzag et al. (1993) is applied in the present analysis. The governing equations are numerically solved using the finite volume methodology. Experiments were conducted to assess the turbulence model. Numerical results of pressure distribution on the front disk surface for different flow conditions when compared to the experimental data indicated that the RNG k-epsilon model is adequate to predict this class of flow.