Potentials of process intensification in the Finnish petrochemical industry

This research was motivated by the need to examine the potential application areas of process intensification technologies in Neste Oil Oyj. According to the company’s interest membrane reactor technology was chosen and applicability of this technology in refining industry was investigated. Moreover, Neste Oil suggested a project which is related to the CO2 capture from FCC unit flue gas stream. The flowrate of the flue gas is 180t/h and consist of approximately 14% by volume CO2. Membrane based absorption process (membrane contactor) was chosen as a potential technique to model CO2 capture from fluid catalytic cracking (FCC) unit effluent. In the design of membrane contactor, a mathematical model was developed to describe CO2 absorption from a gas mixture using monoethanole amine (MEA) aqueous solution. According to the results of literature survey, in the hollow fiber contactor for laminar flow conditions approximately 99 % percent of CO2 can be removed by using a 20 cm in length polyvinylidene fluoride (PDVF) membrane. Furthermore, the design of whole process was performed by using PRO/II simulation software and the CO2 removal efficiency of the whole process obtained as 97 %. The technical and economical comparisons among existing MEA absorption processes were performed to determine the advantages and disadvantages of membrane contactor technology.

Role of oxidation-reduction cycle of iron in direct leaching of zinc concentrate

Direct leaching is an alternative to conventional roast-leach-electrowin (RLE) zinc production method. The basic reaction of direct leach method is the oxidation of sphalerite concentrate in acidic liquid by ferric iron. The reaction mechanism and kinetics, mass transfer and current modifications of zinc concentrate direct leaching process are considered. Particular attention is paid to the oxidation-reduction cycle of iron and its role in direct leaching of zinc concentrate, since it can be one of the limiting factors of the leaching process under certain conditions. The oxidation-reduction cycle of iron was experimentally studied with goal of gaining new knowledge for developing the direct leaching of zinc concentrate. In order to obtain this aim, ferrous iron oxidation experiments were carried out. Affect of such parameters as temperature, pressure, sulfuric acid concentration, ferrous iron and copper concentrations was studied. Based on the experimental results, mathematical model of the ferrous iron oxidation rate was developed. According to results obtained during the study, the reaction rate orders for ferrous iron concentration, oxygen concentration and copper concentration are 0.777, 0.652 and 0.0951 respectively. Values predicted by model were in good concordance with the experimental results. The reliability of estimated parameters was evaluated by MCMC analysis which showed good parameters reliability.

Development of chemical looping combustion process

Chemical looping combustion (CLC) provides a promising technology to help cut carbon dioxide emissions. CLC is based on separated oxidation and reduction processes. Oxygen carrier, which is made from metal and supporting material, is in continuous recirculation between the air and fuel reactors. The CLC process does not require separation unit for carbon dioxide. The fuel reactor can produce an almost pure carbon dioxide feed which decrease costs of carbon capture and storage (CCS). The CLC method is one of the most promising ones for energy efficient carbon capture. A large amount of literature was examined for this study and from it the most promising methods and designs were chosen. These methods and designs were combined as reactor system design which was then sized during the making of this thesis. Sizing was done with a mathematical model that was further improved during the study.

Automated, adaptive methods for forest inventory

Forest inventories are used to estimate forest characteristics and the condition of forest for many different applications: operational tree logging for forest industry, forest health state estimation, carbon balance estimation, land-cover and land use analysis in order to avoid forest degradation etc. Recent inventory methods are strongly based on remote sensing data combined with field sample measurements, which are used to define estimates covering the whole area of interest. Remote sensing data from satellites, aerial photographs or aerial laser scannings are used, depending on the scale of inventory. To be applicable in operational use, forest inventory methods need to be easily adjusted to local conditions of the study area at hand. All the data handling and parameter tuning should be objective and automated as much as possible. The methods also need to be robust when applied to different forest types. Since there generally are no extensive direct physical models connecting the remote sensing data from different sources to the forest parameters that are estimated, mathematical estimation models are of "black-box" type, connecting the independent auxiliary data to dependent response data with linear or nonlinear arbitrary models. To avoid redundant complexity and over-fitting of the model, which is based on up to hundreds of possibly collinear variables extracted from the auxiliary data, variable selection is needed. To connect the auxiliary data to the inventory parameters that are estimated, field work must be performed. In larger study areas with dense forests, field work is expensive, and should therefore be minimized. To get cost-efficient inventories, field work could partly be replaced with information from formerly measured sites, databases. The work in this thesis is devoted to the development of automated, adaptive computation methods for aerial forest inventory. The mathematical model parameter definition steps are automated, and the cost-efficiency is improved by setting up a procedure that utilizes databases in the estimation of new area characteristics.

Microscopic simulation of pigment coating consolidation

The control of coating layer properties is becoming increasingly important as a result of an emerging demand for novel coated paper-based products and an increasing popularity of new coating application methods. The governing mechanisms of microstructure formation dynamics during consolidation and drying are nevertheless, still poorly understood. Some of the difficulties encountered by experimental methods can be overcome by the utilisation of numerical modelling and simulation-based studies of the consolidation process. The objective of this study was to improve the fundamental understanding of pigment coating consolidation and structure formation mechanisms taking place on the microscopic level. Furthermore, it is aimed to relate the impact of process and suspension properties to the microstructure of the coating layer. A mathematical model based on a modified Stokesian dynamics particle simulation technique was developed and applied in several studies of consolidation-related phenomena. The model includes particle-particle and particle-boundary hydrodynamics, colloidal interactions, Born repulsion, and a steric repulsion model. The Brownian motion and a free surface model were incorporated to enable the specific investigation of consolidation and drying. Filter cake stability was simulated in various particle systems, and subjected to a range of base substrate absorption rates and system temperatures. The stability of the filter cake was primarily affected by the absorption rate and size of particles. Temperature was also shown to have an influence. The consolidation of polydisperse systems, with varying wet coating thicknesses, was studied using imposed pilot trial and model-based drying conditions. The results show that drying methods have a clear influence on the microstructure development, on small particle distributions in the coating layer and also on the mobility of particles during consolidation. It is concluded that colloidal properties can significantly impact coating layer shrinkage as well as the internal solids concentration profile. Visualisations of particle system development in time and comparison of systems at different conditions are useful in illustrating coating layer structure formation mechanisms. The results aid in understanding the underlying mechanisms of pigment coating layer consolidation. Guidance is given regarding the relationship between coating process conditions and internal coating slurry properties and their effects on the microstructure of the coating.

Thermal Stress in a Roll-to-Roll Printed Electronics Substrate

The Roll-to-Roll process makes it possible to print electronic products continuously onto a uniform substrate. Printing components on flexible surfaces can bring down the costs of simple electronic devices such as RFID tags, antennas and transistors. The possibility of quickly printing flexible electronic components opens up a wide array of novel products previously too expensive to produce on a large scale. Several different printing methods can be used in Roll-to-Roll printing, such as gravure, spray, offset, flexographic and others. Most of the methods can also be mixed in one production line. Most of them still require years of research to reach a significant commercial level. The research for this thesis was carried out at the Konkuk University Flexible Display Research Center (KU-FDRC) in Seoul, Korea. A system using Roll-to-Roll printing requires that the motion of the web can be controlled in every direction in order to align different layers of ink properly. Between printers the ink is dried with hot air. The effects of thermal expansion on the tension of the web are studied in this work, and a mathematical model was constructed on Matlab and Simulink. Simulations and experiments lead to the conclusion that the thermal expansion of the web has a great influence on the tension of the web. Also, experimental evidence was gained that the particular printing machine used for these experiments at KU-FDRC may have a problem in controlling the speeds of the cylinders which pull the web.

Disruptive innovations at the bottom of the pyramid : can they impact the sustainability of today´s top companies?

Due to the different dynamics required for organizations to serve the emerging market which contains billions of people at the bottom of the pyramid (BOP) coupled with the increasing desire for organizations to grow and be more multinational, organizations need to continually innovate. However, the tendency for large and established companies to ignore the BOP market and rather focus on existing markets, gives an indication of the existence of a vulnerability that potentially disruptive innovations from the BOP will not be recognized in good time for a counter measure. This can be deduced from the fact that good management practice advocates that managers should learn and listen to their customers. Therefore majority of the large existing companies continually focus on their main customer/market with sustaining innovations which leaves aspiring new entrants with an underserved BOP market to experiment with. With the aid of research interviews and an agent-based model (ABM) simulation, this thesis examines the attributes of BOP innovations that can qualify them as disruptive and the possibilities of tangible disruptive innovations arising from the bottom of the pyramid and their underlying drivers. The thesis Furthermore, examines the associated impact of such innovations on the future sustainability of established large companies that are operating in the developed world, particularly those with a primary focus which is targeted towards the market at the top of the pyramid (TOP). Additionally, with the use of a scenario planning model, the research provides an evaluation of the possible evolution and potential sustainability impacts that could emerge, from the interplay of innovations at the two pyramidal market levels and the chosen market focus of organizations – TOP or BOP. Using four scenario quadrants, the thesis demonstrates the resulting possibilities from the interaction between the rate of innovations and the segment focused on by organizations with disruptive era characterizing the paradigm shift quadrant. Furthermore, a mathematical model and two theoretical propositions are developed for further research. As recommendations, the thesis also extends the ambidextrous organizational theory, business model innovation and portfolio diversification as plausible recommendations to limit a catastrophic impact, resulting from disruptive innovations.

Novel Methods for Error Modeling and Parameter Identification of a Redundant Serial-Parallel Hybrid Robot

To obtain the desirable accuracy of a robot, there are two techniques available. The first option would be to make the robot match the nominal mathematic model. In other words, the manufacturing and assembling tolerances of every part would be extremely tight so that all of the various parameters would match the “design” or “nominal” values as closely as possible. This method can satisfy most of the accuracy requirements， but the cost would increase dramatically as the accuracy requirement increases. Alternatively, a more cost-effective solution is to build a manipulator with relaxed manufacturing and assembling tolerances. By modifying the mathematical model in the controller, the actual errors of the robot can be compensated. This is the essence of robot calibration. Simply put, robot calibration is the process of defining an appropriate error model and then identifying the various parameter errors that make the error model match the robot as closely as possible. This work focuses on kinematic calibration of a 10 degree-of-freedom (DOF) redundant serial-parallel hybrid robot. The robot consists of a 4-DOF serial mechanism and a 6-DOF hexapod parallel manipulator. The redundant 4-DOF serial structure is used to enlarge workspace and the 6-DOF hexapod manipulator is used to provide high load capabilities and stiffness for the whole structure. The main objective of the study is to develop a suitable calibration method to improve the accuracy of the redundant serial-parallel hybrid robot. To this end, a Denavit–Hartenberg (DH) hybrid error model and a Product-of-Exponential (POE) error model are developed for error modeling of the proposed robot. Furthermore, two kinds of global optimization methods, i.e. the differential-evolution (DE) algorithm and the Markov Chain Monte Carlo (MCMC) algorithm, are employed to identify the parameter errors of the derived error model. A measurement method based on a 3-2-1 wire-based pose estimation system is proposed and implemented in a Solidworks environment to simulate the real experimental validations. Numerical simulations and Solidworks prototype-model validations are carried out on the hybrid robot to verify the effectiveness, accuracy and robustness of the calibration algorithms.

Modeling and Control of the Power Conversion Unit in a Solid Oxide Fuel Cell Environment

In this doctoral thesis, a power conversion unit for a 10 kWsolid oxide fuel cell is modeled, and a suitable control system is designed. The need for research was identified based on an observation that there was no information available about the characteristics of the solid oxide fuel cell from the perspective of power electronics and the control system, and suitable control methods had not previously been studied in the literature. In addition, because of the digital implementation of the control system, the inherent characteristics of the digital system had to be taken into account in the characteristics of the solid oxide fuel cell (SOFC). The characteristics of the solid oxide fuel cell as well the methods for the modeling and control of the DC/DC converter and the grid converter are studied by a literature survey. Based on the survey, the characteristics of the SOFC as an electrical power source are identified, and a solution to the interfacing of the SOFC in distributed generation is proposed. A mathematical model of the power conversion unit is provided, and the control design for the DC/DC converter and the grid converter is made based on the proposed interfacing solution. The limit cycling phenomenon is identified as a source of low-frequency current ripple, which is found to be insignificant when connected to a grid-tied converter. A method to mitigate a second harmonic originating from the grid interface is proposed, and practical considerations of the operation with the solid oxide fuel cell plant are presented. At the theoretical level, the thesis discusses and summarizes the methods to successfully derive a model for a DC/DC converter, a grid converter, and a power conversion unit. The results of this doctoral thesis can also be used in other applications, and the models and methods can be adopted to similar applications such as photovoltaic systems. When comparing the results with the objectives of the doctoral thesis, we may conclude that the objectives set for the work are met. In this doctoral thesis, theoretical and practical guidelines are presented for the successful control design to connect a SOFC-based distributed generation plant to the utility grid.

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.

Adaptive control of Hydraulic Drive

Adaptive control systems are one of the most significant research directions of modern control theory. It is well known that every mechanical appliance’s behavior noticeably depends on environmental changes, functioning-mode parameter changes and changes in technical characteristics of internal functional devices. An adaptive controller involved in control process allows reducing an influence of such changes. In spite of this such type of control methods is applied seldom due to specifics of a controller designing. The work presented in this paper shows the design process of the adaptive controller built by Lyapunov’s function method for the Hydraulic Drive. The calculation needed and the modeling were conducting with MATLAB® software including Simulink® and Symbolic Math Toolbox™ etc. In the work there was applied the Jacobi matrix linearization of the object’s mathematical model and derivation of the suitable reference models based on Newton’s characteristic polynomial. The intelligent adaptive to nonlinearities algorithm for solving Lyapunov’s equation was developed. Developed algorithm works properly but considered plant is not met requirement of functioning with. The results showed confirmation that adaptive systems application significantly increases possibilities in use devices and might be used for correction a system’s behavior dynamics.

Effect of ultrastructural features on mechanical properties of softwood fiber

The main objective of this study was to develop mathematical model capable to describe the effect of ultrastructural features on the longitudinal modulus of elasticity of softwood fiber. Another objective was to identify, based on ultrastructural features, a potential explanatory factor for the mechanical difference between Norway spruce and Scots pine fibers and to demonstrate its influence utilizing developed modelling tools. According to the literature, the main difference between the pine and spruce fibers is the pit structure, which is clearly different in these fibers. The spruce fiber contains a lot of tiny pits, whereas the pits of the pine fiber are larger and the total number of them is smaller. The effect of the pits on the longitudinal modulus of elasticity of fiber is studied with both the analytical and the numerical model. The results show that, although the spruce fiber seems to contain clearly more pits, larger pits appearing in the pine fiber turn out to have a stronger influence on the longitudinal modulus of elasticity of the fiber. The effect of local variation of microfibril angle which occurs near the pits seems to be minor. Moreover, the results suggest that spruce fibers may have higher ultimate strength due to the more uniform straining behavior.

Surface studies of limestones and dolostones : characterisation using various techniques and batch dissolution experiments with hydrochloric acid solutions

In this thesis, stepwise titration with hydrochloric acid was used to obtain chemical reactivities and dissolution rates of ground limestones and dolostones of varying geological backgrounds (sedimentary, metamorphic or magmatic). Two different ways of conducting the calculations were used: 1) a first order mathematical model was used to calculate extrapolated initial reactivities (and dissolution rates) at pH 4, and 2) a second order mathematical model was used to acquire integrated mean specific chemical reaction constants (and dissolution rates) at pH 5. The calculations of the reactivities and dissolution rates were based on rate of change of pH and particle size distributions of the sample powders obtained by laser diffraction. The initial dissolution rates at pH 4 were repeatedly higher than previously reported literature values, whereas the dissolution rates at pH 5 were consistent with former observations. Reactivities and dissolution rates varied substantially for dolostones, whereas for limestones and calcareous rocks, the variation can be primarily explained by relatively large sample standard deviations. A list of the dolostone samples in a decreasing order of initial reactivity at pH 4 is: 1) metamorphic dolostones with calcite/dolomite ratio higher than about 6% 2) sedimentary dolostones without calcite 3) metamorphic dolostones with calcite/dolomite ratio lower than about 6% The reactivities and dissolution rates were accompanied by a wide range of experimental techniques to characterise the samples, to reveal how different rocks changed during the dissolution process, and to find out which factors had an influence on their chemical reactivities. An emphasis was put on chemical and morphological changes taking place at the surfaces of the particles via X-ray Photoelectron Spectroscopy (XPS) and Scanning Electron Microscopy (SEM). Supporting chemical information was obtained with X-Ray Fluorescence (XRF) measurements of the samples, and Inductively Coupled Plasma-Mass Spectrometry (ICP-MS) and Inductively Coupled Plasma-Optical Emission Spectrometry (ICP-OES) measurements of the solutions used in the reactivity experiments. Information on mineral (modal) compositions and their occurrence was provided by X-Ray Diffraction (XRD), Energy Dispersive X-ray analysis (EDX) and studying thin sections with a petrographic microscope. BET (Brunauer, Emmet, Teller) surface areas were determined from nitrogen physisorption data. Factors increasing chemical reactivity of dolostones and calcareous rocks were found to be sedimentary origin, higher calcite concentration and smaller quartz concentration. Also, it is assumed that finer grain size and larger BET surface areas increase the reactivity although no certain correlation was found in this thesis. Atomic concentrations did not correlate with the reactivities. Sedimentary dolostones, unlike metamorphic ones, were found to have porous surface structures after dissolution. In addition, conventional (XPS) and synchrotron based (HRXPS) X-ray Photoelectron Spectroscopy were used to study bonding environments on calcite and dolomite surfaces. Both samples are insulators, which is why neutralisation measures such as electron flood gun and a conductive mask were used. Surface core level shifts of 0.7 ± 0.1 eV for Ca 2p spectrum of calcite and 0.75 ± 0.05 eV for Mg 2p and Ca 3s spectra of dolomite were obtained. Some satellite features of Ca 2p, C 1s and O 1s spectra have been suggested to be bulk plasmons. The origin of carbide bonds was suggested to be beam assisted interaction with hydrocarbons found on the surface. The results presented in this thesis are of particular importance for choosing raw materials for wet Flue Gas Desulphurisation (FGD) and construction industry. Wet FGD benefits from high reactivity, whereas construction industry can take advantage of slow reactivity of carbonate rocks often used in the facades of fine buildings. Information on chemical bonding environments may help to create more accurate models for water-rock interactions of carbonates.

Forecasting financial weather - can we foresee market sentiment? Spectrum of stock price behavior - NYSE case study

The desire to create a statistical or mathematical model, which would allow predicting the future changes in stock prices, was born many years ago. Economists and mathematicians are trying to solve this task by applying statistical analysis and physical laws, but there are still no satisfactory results. The main reason for this is that a stock exchange is a non-stationary, unstable and complex system, which is influenced by many factors. In this thesis the New York Stock Exchange was considered as the system to be explored. A topological analysis, basic statistical tools and singular value decomposition were conducted for understanding the behavior of the market. Two methods for normalization of initial daily closure prices by Dow Jones and S&P500 were introduced and applied for further analysis. As a result, some unexpected features were identified, such as a shape of distribution of correlation matrix, a bulk of which is shifted to the right hand side with respect to zero. Also non-ergodicity of NYSE was confirmed graphically. It was shown, that singular vectors differ from each other by a constant factor. There are for certain results no clear conclusions from this work, but it creates a good basis for the further analysis of market topology.

Modeling a Network of Heat Exchangers

This research work addresses the problem of building a mathematical model for the given system of heat exchangers and to determine the temperatures, pressures and velocities at the intermediate positions. Such model could be used in nding an optimal design for such a superstructure. To limit the size and computing time a reduced network model was used. The method can be generalized to larger network structures. A mathematical model which includes a system of non-linear equations has been built and solved according to the Newton-Raphson algorithm. The results obtained by the proposed mathematical model were compared with the results obtained by the Paterson approximation and Chen's Approximation. Results of this research work in collaboration with a current ongoing research at the department will optimize the valve positions and hence, minimize the pumping cost and maximize the heat transfer of the system of heat exchangers.