Toimitusvalvonnan tehostaminen projektiympäristössä

Työn tutkimusongelmana on selvittää mitä on toimitusvalvonta, ja miten toimitusvalvontaa voidaan tehostaa projektitoiminnassa valvottaessa ulkopuolelta tilattuja laitteita. Tutkimuksen tuloksena annetaan hyviä käytäntöjä ja neuvoja miten toimitusvalvonta tulisi järjestää taloudellisesti ja kuitenkin tehokkaasti. Tutkimus tehtiin konstruktiivisena case-tutkimuksena havainnoimalla ja haastatteluilla. Raporttia varten tehtiin kirjallisuustutkimus, jonka tuloksena on tutkimuksen teoriaosuus. Tutkimuksen teoriaosassa on käsitelty toimitusvalvonnan työympäristöä yleisemmin ja on käyty läpi yksityiskohtaisesti toimitusvalvojan tehtävät aikataulun, kustannusten ja toimituslaajuuden valvonnassa. Teoriaosassa selvitetään myös niitä perimmäisiä syitä, jotka aiheuttavat ongelmia asiakkaan ja toimittajan välillä koskien toimituksen taloutta, riskien jakoa ja teknisiä ratkaisuja. Tutkimuksen loppuosissa on yhdistetty alkuosan teoria ja tutkimuksesta saadut tehostamistoimenpiteet sekä hyvät käytännöt. Siinä annetaan malli, miten toimitusvalvonta on mahdollista organisoida projektissa niin, että se on taloudellisesti tehokasta ja resurssit on mahdollista suunnata niihin kohteisiin, joissa niistä on eniten hyötyä. Lisäksi käydään läpi yleisimpiä ongelmakohtia.

Prosessilähtöisen laadunhallintajärjestelmän sopeuttaminen perinteiseen linjaorganisaatioon.

Työn tavoitteena oli perehtyä prosessilähtöisen laadunhallintajärjestelmän sopeuttamiseen perinteiseen linjaorganisaatioon. Kirjallisuuslähteitä apuna käyttäen perehdyttiin laadunhallintajärjestelmästandardien vaatimuksiin ja organisaatioiden erilaisiin johtamismalleihin. Prosessilähtöisen toimintamallin sopeuttamiseen linjaorganisaatioon liittyy olennaisesti esitutkimus, jonka avulla saadaan tietoa organisaation senhetkisestä tilasta. Tämä esitutkimus suoritettiin kohdeyrityksen henkilöstön keskuudessa eri tutkimusmenetelmiä käyttäen. Saaduista esitutkimuksen tuloksista yrityksen johto voi päättää toimenpiteistä, joilla edesautetaan uuden toimintamallin ja samalla yrityksen vision, strategioiden sekä arvojen saattamista käytännön toimintaan. Uusi prosessilähtöinen laadunhallintajärjestelmä, ISO 9001 -standardi, asettaa niin kyseisen laadunhallintajärjestelmän omaaville kuin sitä tavoittelevillekin yrityksille uudet vaatimukset. Pukkila Oy Ab:n lähtiessä täyttämään näitä vaatimuksia se on kokonaan nykyisestä poikkeavan toimintamallin edessä, joka edellyttää muutosta koko organisaatiolta. Prosessilähtöinen toimintamalli vaatii kokonaan nykyisen toimintamallin kyseenalaistamista, jolloin muutos kohdistuu myös nykyiseen yrityksen johtamisjärjestelmään. Laadunhallintajärjestelmä painottaa tässä toimintamallissa asiakaslähtöistä ajattelutapaa sekä toiminnallisella että teknillisellä puolella. Toimenpidesuosituksena yrityksen johdolle esitetään perehtymistä organisaation eri tasojen yhteiseen sosiaaliseen todellisuuteen sekä muuttamaan osastopohjainen johtamisjärjestelmä prosessilähtöistä toimintaa tukevaksi.

Tuotannon ja ostotoiminnan tarkastelu SCOR-mallilla

Diplomityössä perehdyttiin toimitusketjun hallinnan ja johtamisen apuvälineenä käytettävään SCOR-malliin (Supply Chain Operations Refence-model). SCOR on standardimalli, josta jokaisen tulee poimia oman toiminnan kehittämisen ja tehostamisen kannalta tärkeät asiat. SCOR on hyvä instrumentti kasvun hallinnassa.Työn teoriaosa käsittelee logistiikkaa, toimitusketjun hallintaa, ostotoimintaa ja SCOR-mallia. Soveltavassa osassa käsitellään SCOR-mallin avulla sisäilma alalla toimivan yrityksen prosesseja. Työssä keskityttiin tarkastelemaan yrityksen ostotoiminta ja valmistusprosesseja SCOR-mallin avulla.SCOR-mallin käyttö oston ja tuotannon analysoinnissa toi esiin useita toiminnan kehitysajatuksia. Näitä olivat uusien menetelmien käyttöönotto varaston ja tuotannon hallittavuuden parantamiseksi, osto-organisaation selkiyttäminen ja toiminnanohjausjärjestelmän kehittäminen.

Metallien kompleksinmuodostus liuosfaasissa

Liuoksessa metallit muodostavat erilaisia koordinatioyhdisteitä epäorgaanisten ja orgaanisten anionien ja neutraalien molekyylien kanssa. Erityisesti siirtymämetalleilla on voimakas taipumus kompleksiyhdisteiden muodostamiseen elektroneja sisältävien 3-, 4-, ja 5d orbitaaliensa johdosta. Samassa liuoksessa voi samanaikaisesti esiintyä useita erilaisia, mutta samoista lähtöaineista muodostuneita, kompleksiyhdisteitä. Kompleksinmuodostusreaktiot ovat tasapainoreaktioita. Usein tasapainovakiot on esitetty termodynaamisina tasapainovakioina eli ne ovat päteviä standarditilassa. Standarditilan tasapainovakioista voidaan johtaa missä tahansa liuoksessa pätevät vakiot erilaisten Debye-Hückel-teoriasta johdettujen laskentamenetelmien avulla. Metalli-ligandiparin jakautuminen erilaisiksi kompleksiyhdisteiksi voidaan mallintaa kun tunnetaan muodostumisreaktioiden tasapainovakiot. Muodostumisreaktioiden tasapainovakioiden yhtälöistä voidaan johtaa epälineaarinen yhtälöryhmä, joka voidaan ratkaista jollakin numeerisella ratkaisimella. Esimerkiksi Matlab-ohjelmiston sisältämä fsolve-ratkaisin soveltuu tällaiseen tehtävään. Osana tätä työtä on kirjoitettu Matlab-sovellus, jolla voidaan mallintaa kationi-ligandiparin jakautumista koordinaatioyhdisteiksi tunnettujen tasapainovakioiden perusteella.

Service Productization and Product Servitization in The Contract Manufacturing Company

The role of contract manufacturing and subcontracting has been seen in black and white in product and service point of view. It used to be seen either as a product or a service. In the thesis product-service system, offering combining products and services, was discussed. Theory was created from two perspectives; Service productization via Business Model generation and product servitization via New Service Development process. Target for the case study was to point out new ways of service thinking and ways for changing business environment in contract manufacturing, especially in customer satisfaction and profitability point of view. The case study is following the New Service Development process phases. First ideas were collected from literature and via sales management interviews. Service offering and tool for service requirement evaluation was created. Last financial results of example service scenarios were calculated. It is recommended to take service offering into internal use and further develop it into modular service model. It is also recommended to take created customer service requirement evaluation tool into use for capturing customer service needs but also for communicating those internally.

Chemical changes in biomass during torrefaction

Torrefaction is moderate thermal treatment (~200-300 °C) of biomass in an inert atmosphere. The torrefied fuel offers advantages to traditional biomass, such as higher heating value, reduced hydrophilic nature, increased its resistance to biological decay, and improved grindability. These factors could, for instance, lead to better handling and storage of biomass and increased use of biomass in pulverized combustors. In this work, we look at several aspects of changes in the biomass during torrefaction. We investigate the fate of carboxylic groups during torrefaction and its dependency to equilibrium moisture content. The changes in the wood components including carbohydrates, lignin, extractable materials and ashforming matters are also studied. And at last, the effect of K on torrefaction is investigated and then modeled. In biomass, carboxylic sites are partially responsible for its hydrophilic characteristic. These sites are degraded to varying extents during torrefaction. In this work, methylene blue sorption and potentiometric titration were applied to measure the concentration of carboxylic groups in torrefied spruce wood. The results from both methods were applicable and the values agreed well. A decrease in the equilibrium moisture content at different humidity was also measured for the torrefied wood samples, which is in good agreement with the decrease in carboxylic group contents. Thus, both methods offer a means of directly measuring the decomposition of carboxylic groups in biomass during torrefaction as a valuable parameter in evaluating the extent of torrefaction. This provides new information to the chemical changes occurring during torrefaction. The effect of torrefaction temperature on the chemistry of birch wood was investigated. The samples were from a pilot plant at Energy research Center of the Netherlands (ECN). And in that way they were representative of industrially produced samples. Sugar analysis was applied to analyze the hemicellulose and cellulose content during torrefaction. The results show a significant degradation of hemicellulose already at 240 °C, while cellulose degradation becomes significant above 270 °C torrefaction. Several methods including Klason lignin method, solid state NMR and Py-GC-MS analyses were applied to measure the changes in lignin during torrefaction. The changes in the ratio of phenyl, guaiacyl and syringyl units show that lignin degrades already at 240 °C to a small extent. To investigate the changes in the extractives from acetone extraction during torrefaction, gravimetric method, HP-SEC and GC-FID followed by GC-MS analysis were performed. The content of acetone-extractable material increases already at 240 °C torrefaction through the degradation of carbohydrate and lignin. The molecular weight of the acetone-extractable material decreases with increasing the torrefaction temperature. The formation of some valuable materials like syringaresinol or vanillin is also observed which is important from biorefinery perspective. To investigate the change in the chemical association of ash-forming elements in birch wood during torrefaction, chemical fractionation was performed on the original and torrefied birch samples. These results give a first understanding of the changes in the association of ashforming elements during torrefaction. The most significant changes can be seen in the distribution of calcium, magnesium and manganese, with some change in water solubility seen in potassium. These changes may in part be due to the destruction of carboxylic groups. In addition to some changes in water and acid solubility of phosphorous, a clear decrease in the concentration of both chlorine and sulfur was observed. This would be a significant additional benefit for the combustion of torrefied biomass. Another objective of this work is studying the impact of organically bound K, Na, Ca and Mn on mass loss of biomass during torrefaction. These elements were of interest because they have been shown to be catalytically active in solid fuels during pyrolysis and/or gasification. The biomasses were first acid washed to remove the ash-forming matters and then organic sites were doped with K, Na, Ca or Mn. The results show that K and Na bound to organic sites can significantly increase the mass loss during torrefaction. It is also seen that Mn bound to organic sites increases the mass loss and Ca addition does not influence the mass loss rate on torrefaction. This increase in mass loss during torrefaction with alkali addition is unlike what has been found in the case of pyrolysis where alkali addition resulted in a reduced mass loss. These results are important for the future operation of torrefaction plants, which will likely be designed to handle various biomasses with significantly different contents of K. The results imply that shorter retention times are possible for high K-containing biomasses. The mass loss of spruce wood with different content of K was modeled using a two-step reaction model based on four kinetic rate constants. The results show that it is possible to model the mass loss of spruce wood doped with different levels of K using the same activation energies but different pre-exponential factors for the rate constants. Three of the pre-exponential factors increased linearly with increasing K content, while one of the preexponential factors decreased with increasing K content. Therefore, a new torrefaction model was formulated using the hemicellulose and cellulose content and K content. The new torrefaction model was validated against the mass loss during the torrefaction of aspen, miscanthus, straw and bark. There is good agreement between the model and the experimental data for the other biomasses, except bark. For bark, the mass loss of acetone extractable material is also needed to be taken into account. The new model can describe the kinetics of mass loss during torrefaction of different types of biomass. This is important for considering fuel flexibility in torrefaction plants.

Computational and Experimental Analysis of Flow Field in the Diffusers of Centrifugal Compressors

Centrifugal compressors are widely used for example in process industry, oil and gas industry, in small gas turbines and turbochargers. In order to achieve lower consumption of energy and operation costs the efficiency of the compressor needs to be improve. In the present work different pinches and low solidity vaned diffusers were utilized in order to improve the efficiency of a medium size centrifugal compressor. In this study, pinch means the decrement of the diffuser flow passage height. First different geometries were analyzed using computational fluid dynamics. The flow solver Finflo was used to solve the flow field. Finflo is a Navier-Stokes solver. The solver is capable to solve compressible, incompressible, steady and unsteady flow fields. Chien's k-e turbulence model was used. One of the numerically investigated pinched diffuser and one low solidity vaned diffuser were studied experimentally. The overall performance of the compressor and the static pressure distribution before and after the diffuser were measured. The flow entering and leaving the diffuser was measured using a three-hole Cobra-probe and Kiel-probes. The pinch and the low solidity vaned diffuser increased the efficiency of the compressor. Highest isentropic efficiency increment obtained was 3\% of the design isentropic efficiency of the original geometry. It was noticed in the numerical results that the pinch made to the hub and the shroud wall was most beneficial to the operation of the compressor. Also the pinch made to the hub was better than the pinchmade to the shroud. The pinch did not affect the operation range of the compressor, but the low solidity vaned diffuser slightly decreased the operation range.The unsteady phenomena in the vaneless diffuser were studied experimentally andnumerically. The unsteady static pressure was measured at the diffuser inlet and outlet, and time-accurate numerical simulation was conducted. The unsteady static pressure showed that most of the pressure variations lay at the passing frequency of every second blade. The pressure variations did not vanish in the diffuser and were visible at the diffuser outlet. However, the amplitude of the pressure variations decreased in the diffuser. The time-accurate calculations showed quite a good agreement with the measured data. Agreement was very good at the design operation point, even though the computational grid was not dense enough inthe volute and in the exit cone. The time-accurate calculation over-predicted the amplitude of the pressure variations at high flow.

Selective catalytic reduction of nitrogen oxides with ammonia in forced unsteady state reactors - Case based reasoning and mathematical model simulation reasoning

The application of forced unsteady-state reactors in case of selective catalytic reduction of nitrogen oxides (NOx) with ammonia (NH3) is sustained by the fact that favorable temperature and composition distributions which cannot be achieved in any steady-state regime can be obtained by means of unsteady-state operations. In a normal way of operation the low exothermicity of the selective catalytic reduction (SCR) reaction (usually carried out in the range of 280-350°C) is not enough to maintain by itself the chemical reaction. A normal mode of operation usually requires supply of supplementary heat increasing in this way the overall process operation cost. Through forced unsteady-state operation, the main advantage that can be obtained when exothermic reactions take place is the possibility of trapping, beside the ammonia, the moving heat wave inside the catalytic bed. The unsteady state-operation enables the exploitation of the thermal storage capacity of the catalyticbed. The catalytic bed acts as a regenerative heat exchanger allowing auto-thermal behaviour when the adiabatic temperature rise is low. Finding the optimum reactor configuration, employing the most suitable operation model and identifying the reactor behavior are highly important steps in order to configure a proper device for industrial applications. The Reverse Flow Reactor (RFR) - a forced unsteady state reactor - corresponds to the above mentioned characteristics and may be employed as an efficient device for the treatment of dilute pollutant mixtures. As a main disadvantage, beside its advantages, the RFR presents the 'wash out' phenomena. This phenomenon represents emissions of unconverted reactants at every switch of the flow direction. As a consequence our attention was focused on finding an alternative reactor configuration for RFR which is not affected by the incontrollable emissions of unconverted reactants. In this respect the Reactor Network (RN) was investigated. Its configuration consists of several reactors connected in a closed sequence, simulating a moving bed by changing the reactants feeding position. In the RN the flow direction is maintained in the same way ensuring uniformcatalyst exploitation and in the same time the 'wash out' phenomena is annulated. The simulated moving bed (SMB) can operate in transient mode giving practically constant exit concentration and high conversion levels. The main advantage of the reactor network operation is emphasizedby the possibility to obtain auto-thermal behavior with nearly uniformcatalyst utilization. However, the reactor network presents only a small range of switching times which allow to reach and to maintain an ignited state. Even so a proper study of the complex behavior of the RN may give the necessary information to overcome all the difficulties that can appear in the RN operation. The unsteady-state reactors complexity arises from the fact that these reactor types are characterized by short contact times and complex interaction between heat and mass transportphenomena. Such complex interactions can give rise to a remarkable complex dynamic behavior characterized by a set of spatial-temporal patterns, chaotic changes in concentration and traveling waves of heat or chemical reactivity. The main efforts of the current research studies concern the improvement of contact modalities between reactants, the possibility of thermal wave storage inside the reactor and the improvement of the kinetic activity of the catalyst used. Paying attention to the above mentioned aspects is important when higher activity even at low feeding temperatures and low emissions of unconverted reactants are the main operation concerns. Also, the prediction of the reactor pseudo or steady-state performance (regarding the conversion, selectivity and thermal behavior) and the dynamicreactor response during exploitation are important aspects in finding the optimal control strategy for the forced unsteady state catalytic tubular reactors. The design of an adapted reactor requires knowledge about the influence of its operating conditions on the overall process performance and a precise evaluation of the operating parameters rage for which a sustained dynamic behavior is obtained. An apriori estimation of the system parameters result in diminution of the computational efforts. Usually the convergence of unsteady state reactor systems requires integration over hundreds of cycles depending on the initial guess of the parameter values. The investigation of various operation models and thermal transfer strategies give reliable means to obtain recuperative and regenerative devices which are capable to maintain an auto-thermal behavior in case of low exothermic reactions. In the present research work a gradual analysis of the SCR of NOx with ammonia process in forced unsteady-state reactors was realized. The investigation covers the presentationof the general problematic related to the effect of noxious emissions in the environment, the analysis of the suitable catalysts types for the process, the mathematical analysis approach for modeling and finding the system solutions and the experimental investigation of the device found to be more suitable for the present process. In order to gain information about the forced unsteady state reactor design, operation, important system parameters and their values, mathematical description, mathematicalmethod for solving systems of partial differential equations and other specific aspects, in a fast and easy way, and a case based reasoning (CBR) approach has been used. This approach, using the experience of past similarproblems and their adapted solutions, may provide a method for gaining informations and solutions for new problems related to the forced unsteady state reactors technology. As a consequence a CBR system was implemented and a corresponding tool was developed. Further on, grooving up the hypothesis of isothermal operation, the investigation by means of numerical simulation of the feasibility of the SCR of NOx with ammonia in the RFRand in the RN with variable feeding position was realized. The hypothesis of non-isothermal operation was taken into account because in our opinion ifa commercial catalyst is considered, is not possible to modify the chemical activity and its adsorptive capacity to improve the operation butis possible to change the operation regime. In order to identify the most suitable device for the unsteady state reduction of NOx with ammonia, considering the perspective of recuperative and regenerative devices, a comparative analysis of the above mentioned two devices performance was realized. The assumption of isothermal conditions in the beginningof the forced unsteadystate investigation allowed the simplification of the analysis enabling to focus on the impact of the conditions and mode of operation on the dynamic features caused by the trapping of one reactant in the reactor, without considering the impact of thermal effect on overall reactor performance. The non-isothermal system approach has been investigated in order to point out the important influence of the thermal effect on overall reactor performance, studying the possibility of RFR and RN utilization as recuperative and regenerative devices and the possibility of achieving a sustained auto-thermal behavior in case of lowexothermic reaction of SCR of NOx with ammonia and low temperature gasfeeding. Beside the influence of the thermal effect, the influence of the principal operating parameters, as switching time, inlet flow rate and initial catalyst temperature have been stressed. This analysis is important not only because it allows a comparison between the two devices and optimisation of the operation, but also the switching time is the main operating parameter. An appropriate choice of this parameter enables the fulfilment of the process constraints. The level of the conversions achieved, the more uniform temperature profiles, the uniformity ofcatalyst exploitation and the much simpler mode of operation imposed the RN as a much more suitable device for SCR of NOx with ammonia, in usual operation and also in the perspective of control strategy implementation. Theoretical simplified models have also been proposed in order to describe the forced unsteady state reactors performance and to estimate their internal temperature and concentration profiles. The general idea was to extend the study of catalytic reactor dynamics taking into account the perspectives that haven't been analyzed yet. The experimental investigation ofRN revealed a good agreement between the data obtained by model simulation and the ones obtained experimentally.

Cost allocation of corporate IT-services

Cost allocation is an inescapable problem in nearly every organization and in nearly every facet of accounting. Within large corporations there are several different types of units, like profit-making business units and non-profit service units. In order to evaluate the performance of the business units and to fund the operations of service units, the expenses of service production need to be allocated to the business units benefiting from the services.The objective of this thesis was to find good and fair allocating factors for the costs of corporate wide IT services. In order to reach this objective, the cost allocation process was studied in general and an overview of cost structure was established. All possible cost driver candidates were mapped and their good and bad properties were weighed. The cost allocation problem was handled separately according to organizational division of corporate IT department: infrastructure, administrative systems, sales system and e-business. The emphasis was on two largest cost groups: infrastructure costs and sales system costs. As a result of the study an allocation model is presented. It contains categorization of the costs, selected cost drivers and cost distributions for the current year.

Developing A Knowledge Management Process Model for Internal IT Organization, Case: Stora Enso Oyj, Corporate IT

The objectives of this Master’s Thesis were to find out what kind of knowledge management strategy would fit best an IT organization that uses ITIL (Information Technology Infrastructure Library) framework for IT Service Management and to create a knowledge management process model to support chosen strategy. The empirical material for this research was collected through qualitative semi-structured interviews of a case organization Stora Enso Corporate IT. The results of the qualitative interviews indicate that codification knowledge management strategy would fit best for the case organization. The knowledge management process model was created based on earlier studies and a literature of knowledge management. The model was evaluated in the interview research and the results showed that the created process model is realistic, useful, and it responds to a real life phenomenon.

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.

Methods for construction and analysis of computational models in systems biology : applications to the modelling of the heat shock response and the self-assembly of intermediate filaments

Systems biology is a new, emerging and rapidly developing, multidisciplinary research field that aims to study biochemical and biological systems from a holistic perspective, with the goal of providing a comprehensive, system- level understanding of cellular behaviour. In this way, it addresses one of the greatest challenges faced by contemporary biology, which is to compre- hend the function of complex biological systems. Systems biology combines various methods that originate from scientific disciplines such as molecu- lar biology, chemistry, engineering sciences, mathematics, computer science and systems theory. Systems biology, unlike “traditional” biology, focuses on high-level concepts such as: network, component, robustness, efficiency, control, regulation, hierarchical design, synchronization, concurrency, and many others. The very terminology of systems biology is “foreign” to “tra- ditional” biology, marks its drastic shift in the research paradigm and it indicates close linkage of systems biology to computer science. One of the basic tools utilized in systems biology is the mathematical modelling of life processes tightly linked to experimental practice. The stud- ies contained in this thesis revolve around a number of challenges commonly encountered in the computational modelling in systems biology. The re- search comprises of the development and application of a broad range of methods originating in the fields of computer science and mathematics for construction and analysis of computational models in systems biology. In particular, the performed research is setup in the context of two biolog- ical phenomena chosen as modelling case studies: 1) the eukaryotic heat shock response and 2) the in vitro self-assembly of intermediate filaments, one of the main constituents of the cytoskeleton. The range of presented approaches spans from heuristic, through numerical and statistical to ana- lytical methods applied in the effort to formally describe and analyse the two biological processes. We notice however, that although applied to cer- tain case studies, the presented methods are not limited to them and can be utilized in the analysis of other biological mechanisms as well as com- plex systems in general. The full range of developed and applied modelling techniques as well as model analysis methodologies constitutes a rich mod- elling framework. Moreover, the presentation of the developed methods, their application to the two case studies and the discussions concerning their potentials and limitations point to the difficulties and challenges one encounters in computational modelling of biological systems. The problems of model identifiability, model comparison, model refinement, model inte- gration and extension, choice of the proper modelling framework and level of abstraction, or the choice of the proper scope of the model run through this thesis.

Application of Computational Fluid Dynamics in Modeling Blood Flow in Human Thoracic Aorta

The aim of this study was to simulate blood flow in thoracic human aorta and understand the role of flow dynamics in the initialization and localization of atherosclerotic plaque in human thoracic aorta. The blood flow dynamics in idealized and realistic models of human thoracic aorta were numerically simulated in three idealized and two realistic thoracic aorta models. The idealized models of thoracic aorta were reconstructed with measurements available from literature, and the realistic models of thoracic aorta were constructed by image processing Computed Tomographic (CT) images. The CT images were made available by South Karelia Central Hospital in Lappeenranta. The reconstruction of thoracic aorta consisted of operations, such as contrast adjustment, image segmentations, and 3D surface rendering. Additional design operations were performed to make the aorta model compatible for the numerical method based computer code. The image processing and design operations were performed with specialized medical image processing software. Pulsatile pressure and velocity boundary conditions were deployed as inlet boundary conditions. The blood flow was assumed homogeneous and incompressible. The blood was assumed to be a Newtonian fluid. The simulations with idealized models of thoracic aorta were carried out with Finite Element Method based computer code, while the simulations with realistic models of thoracic aorta were carried out with Finite Volume Method based computer code. Simulations were carried out for four cardiac cycles. The distribution of flow, pressure and Wall Shear Stress (WSS) observed during the fourth cardiac cycle were extensively analyzed. The aim of carrying out the simulations with idealized model was to get an estimate of flow dynamics in a realistic aorta model. The motive behind the choice of three aorta models with distinct features was to understand the dependence of flow dynamics on aorta anatomy. Highly disturbed and nonuniform distribution of velocity and WSS was observed in aortic arch, near brachiocephalic, left common artery, and left subclavian artery. On the other hand, the WSS profiles at the roots of branches show significant differences with geometry variation of aorta and branches. The comparison of instantaneous WSS profiles revealed that the model with straight branching arteries had relatively lower WSS compared to that in the aorta model with curved branches. In addition to this, significant differences were observed in the spatial and temporal profiles of WSS, flow, and pressure. The study with idealized model was extended to study blood flow in thoracic aorta under the effects of hypertension and hypotension. One of the idealized aorta models was modified along with the boundary conditions to mimic the thoracic aorta under the effects of hypertension and hypotension. The results of simulations with realistic models extracted from CT scans demonstrated more realistic flow dynamics than that in the idealized models. During systole, the velocity in ascending aorta was skewed towards the outer wall of aortic arch. The flow develops secondary flow patterns as it moves downstream towards aortic arch. Unlike idealized models, the distribution of flow was nonplanar and heavily guided by the artery anatomy. Flow cavitation was observed in the aorta model which was imaged giving longer branches. This could not be properly observed in the model with imaging containing a shorter length for aortic branches. The flow circulation was also observed in the inner wall of the aortic arch. However, during the diastole, the flow profiles were almost flat and regular due the acceleration of flow at the inlet. The flow profiles were weakly turbulent during the flow reversal. The complex flow patterns caused a non-uniform distribution of WSS. High WSS was distributed at the junction of branches and aortic arch. Low WSS was distributed at the proximal part of the junction, while intermedium WSS was distributed in the distal part of the junction. The pulsatile nature of the inflow caused oscillating WSS at the branch entry region and inner curvature of aortic arch. Based on the WSS distribution in the realistic model, one of the aorta models was altered to induce artificial atherosclerotic plaque at the branch entry region and inner curvature of aortic arch. Atherosclerotic plaque causing 50% blockage of lumen was introduced in brachiocephalic artery, common carotid artery, left subclavian artery, and aortic arch. The aim of this part of the study was first to study the effect of stenosis on flow and WSS distribution, understand the effect of shape of atherosclerotic plaque on flow and WSS distribution, and finally to investigate the effect of lumen blockage severity on flow and WSS distributions. The results revealed that the distribution of WSS is significantly affected by plaque with mere 50% stenosis. The asymmetric shape of stenosis causes higher WSS in branching arteries than in the cases with symmetric plaque. The flow dynamics within thoracic aorta models has been extensively studied and reported here. The effects of pressure and arterial anatomy on the flow dynamic were investigated. The distribution of complex flow and WSS is correlated with the localization of atherosclerosis. With the available results we can conclude that the thoracic aorta, with complex anatomy is the most vulnerable artery for the localization and development of atherosclerosis. The flow dynamics and arterial anatomy play a role in the localization of atherosclerosis. The patient specific image based models can be used to diagnose the locations in the aorta vulnerable to the development of arterial diseases such as atherosclerosis.

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).

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.