Calculation of breed direct and maternal genetic fractions and breed specific direct and maternal heterozygosity for crossbreeding data

Teaching, research, and herd breeding applications may require calculation of breed additive contributions for direct and maternal genetic effects and fractions of heterozygosity associated with breed specific direct and maternal heterosis effects. These coefficients can be obtained from the first NB rows of a pseudo numerator relationship matrix where the first NB rows represent fractional contributions by breed to each animal or group representing a specific breed cross. The table begins with an NB x NB identity matrix representing pure breeds. Initial animals or representative crosses must be purebreds or two-breed crosses. Parents of initial purebreds are represented by the corresponding column and initial two-breed cross progeny by the two corresponding columns of the identity matrix. After that, usual rules are used to calculate the NB column entries corresponding to breeds for each animal. The NB entries are fractions of genes expected to be contributed by each of the pure breeds and correspond to the breed additive direct fractions. Entries in the column corresponding to the dam represent breed additive maternal fractions. Breed specific direct heterozygosity coefficients are entries of an NB x NB matrix formed by the outer product of the two NB by 1 columns associated with sire and dam of the animal. One minus sum of the diagonals represents total direct heterozygosity. Similarly, the NB x NB matrix formed by the outer product of columns associated with sire of dam and dam of dam contains breed specific maternal heterozygosity coefficients. These steps can be programmed to create covariates to merge with data. If X represents these coefficients for all unique breed crosses, then the reduced row echelon form function of MATLAB or SAS can be used on X to determine estimable functions of additive breed direct and maternal effects and breed specific direct and maternal heterosis effects

Virtual Testing of Active Magnetic Bearing Systems based on Design Guidelines given by the Standards

Active Magnetic Bearings offer many advantages that have brought new applications to the industry. However, similarly to all new technology, active magnetic bearings also have downsides and one of those is the low standardization level. This thesis is studying mainly the ISO 14839 standard and more specifically the system verification methods. These verifying methods are conducted using a practical test with an existing active magnetic bearing system. The system is simulated with Matlab using rotor-bearing dynamics toolbox, but this study does not include the exact simulation code or a direct algebra calculation. However, this study provides the proof that standardized simulation methods can be applied in practical problems.

Yliopiston sähköverkon mallintaminen

Lappeenrannan teknillinen yliopisto tutkii älykkäiden sähköverkkojen kehittämistä. Yliopisto on hankkinut sähköverkkoonsa tuuliturbiinin ja aurinkopaneeleita, joilla pystytään tuottamaan sähköenergiaa sähköverkkoon. Näitä tuotantoja voidaan käyttää myös tutkimuksessa. Tässä työssä luodaan simulaatiomalli yliopiston sähköverkosta Matlab® Simulink® -ohjelmalla. Simulaatiomalliin mallinnetaan yliopiston sisäinen keskijänniteverkko ja osa pienjänniteverkosta. Simulaatiomalli toteutetaan ohjelman valmiilla komponenteilla, joihin lasketaan tarvittavat parametrit. Tuuliturbiinin ja aurinkopaneelien sähköntuotantotehot määritetään säätiladatojen avulla. Verkon komponenteille lasketaan arvot komponenttien tyyppitietojen perusteella ja asetetaan simulaatiomallin parametreiksi. Simulaatiomalli luodaan yliopiston sisäisen verkon tehonjaon tarkastelemiseksi. Työssä selvitetään myös mahdollisuuksia luodun simulaatiomallin käyttämiseen vikatilanteiden tarkastelussa.

Application and development of numerical methods for the modelling of innovative gas cooled fission reactors

Innovative gas cooled reactors, such as the pebble bed reactor (PBR) and the gas cooled fast reactor (GFR) offer higher efficiency and new application areas for nuclear energy. Numerical methods were applied and developed to analyse the specific features of these reactor types with fully three dimensional calculation models. In the first part of this thesis, discrete element method (DEM) was used for a physically realistic modelling of the packing of fuel pebbles in PBR geometries and methods were developed for utilising the DEM results in subsequent reactor physics and thermal-hydraulics calculations. In the second part, the flow and heat transfer for a single gas cooled fuel rod of a GFR were investigated with computational fluid dynamics (CFD) methods. An in-house DEM implementation was validated and used for packing simulations, in which the effect of several parameters on the resulting average packing density was investigated. The restitution coefficient was found out to have the most significant effect. The results can be utilised in further work to obtain a pebble bed with a specific packing density. The packing structures of selected pebble beds were also analysed in detail and local variations in the packing density were observed, which should be taken into account especially in the reactor core thermal-hydraulic analyses. Two open source DEM codes were used to produce stochastic pebble bed configurations to add realism and improve the accuracy of criticality calculations performed with the Monte Carlo reactor physics code Serpent. Russian ASTRA criticality experiments were calculated. Pebble beds corresponding to the experimental specifications within measurement uncertainties were produced in DEM simulations and successfully exported into the subsequent reactor physics analysis. With the developed approach, two typical issues in Monte Carlo reactor physics calculations of pebble bed geometries were avoided. A novel method was developed and implemented as a MATLAB code to calculate porosities in the cells of a CFD calculation mesh constructed over a pebble bed obtained from DEM simulations. The code was further developed to distribute power and temperature data accurately between discrete based reactor physics and continuum based thermal-hydraulics models to enable coupled reactor core calculations. The developed method was also found useful for analysing sphere packings in general. CFD calculations were performed to investigate the pressure losses and heat transfer in three dimensional air cooled smooth and rib roughened rod geometries, housed inside a hexagonal flow channel representing a sub-channel of a single fuel rod of a GFR. The CFD geometry represented the test section of the L-STAR experimental facility at Karlsruhe Institute of Technology and the calculation results were compared to the corresponding experimental results. Knowledge was gained of the adequacy of various turbulence models and of the modelling requirements and issues related to the specific application. The obtained pressure loss results were in a relatively good agreement with the experimental data. Heat transfer in the smooth rod geometry was somewhat under predicted, which can partly be explained by unaccounted heat losses and uncertainties. In the rib roughened geometry heat transfer was severely under predicted by the used realisable k − epsilon turbulence model. An additional calculation with a v2 − f turbulence model showed significant improvement in the heat transfer results, which is most likely due to the better performance of the model in separated flow problems. Further investigations are suggested before using CFD to make conclusions of the heat transfer performance of rib roughened GFR fuel rod geometries. It is suggested that the viewpoints of numerical modelling are included in the planning of experiments to ease the challenging model construction and simulations and to avoid introducing additional sources of uncertainties. To facilitate the use of advanced calculation approaches, multi-physical aspects in experiments should also be considered and documented in a reasonable detail.

Tasaisemmin jakaantuneen CIELAB-tyyppisen väriavaruuden määrittäminen värieroellipsien pääakselien eroja minimoimalla

Tässä työssä johdetaan lineaarimuunnoksella CIE x y z-värinsovitusfunktioista uudet värinsovitusfunktiot. Tarvittava muunnosmatriisi etsitään optimoimalla CIE ja BFD-RIT värieroellipsejä Matlab-ympäristössä. Työn tuloksena saatiin muunnosmatriisi, ja sillä muunnetut uudet värinsovitusfunktiot ja CIELAB-tyyppinen väriavaruus. Euklidisella etäisyydellä mitattuna CIE ja BFD-RIT värieroellipsien muoto ja koko paranivat noin kolmanneksen, mikä oli myös tavoitteena.

Simulation of the Interaction Mean Free Path between Neutrons and TRISO Particles

The interaction mean free path between neutrons and TRISO particles is simulated using scripts written in MATLAB to solve the increasing error present with an increase in the packing factor in the reactor physics code Serpent. Their movement is tracked both in an unbounded and in a bounded space. Their track is calculated, depending on the program, linearly directly using the position vectors of the neutrons and the surface equations of all the fuel particles; by dividing the space in multiple subspaces, each of which contain a fraction of the total number of particles, and choosing the particles from those subspaces through which the neutron passes through; or by choosing the particles that lie within an infinite cylinder formed on the movement axis of the neutron. The estimate from the current analytical model, based on an exponential distribution, for the mean free path, utilized by Serpent, is used as a reference result. The results from the implicit model in Serpent imply a too long mean free path with high packing factors. The received results support this observation by producing, with a packing factor of 17 %, approximately 2.46 % shorter mean free path compared to the reference model. This is supported by the packing factor experienced by the neutron, the simulation of which resulted in a 17.29 % packing factor. It was also observed that the neutrons leaving from the surfaces of the fuel particles, in contrast to those starting inside the moderator, do not follow the exponential distribution. The current model, as it is, is thus not valid in the determination of the free path lengths of the neutrons.

Massataseiden laskentatyökalujen kehitys bench-kokoluokan reaktiivisuuskoelaitteelle

Tässä työssä tarkastellaan kehitettyjä laskentatyökaluja VTT:n Jyväskylän toimipisteen bench-kokoluokan reaktiivisuuskoelaitteiston tulosten käsittelyn nopeuttamiseksi. Työkalut kehitettiin Matlab- ohjelmistossa. Toinen työkalu toimii Savcor Wedge prosessinhallintaohjelmassa reaaliaikaisena laskentatyökaluna. Tulosten varmentamiseen ja tarkistamiseen kehitetty työkalu toimii Matlabissa mittauksista tuodulla mittausdatalla. Työn aluksi tarkastellaan yleisesti palamista ja sen eri vaiheita. Tämän jälkeen esitellään koelaitteisto, jolle analyysityökalut kehitettiin. Koelaitteiston kanssa esitellään esimerkki polttoainepanoskokeista, joiden tutkimisen nopeuttamiseen laskentatyökalut soveltuu. Laskentatyökalujen toiminta esitetään pääpiirteittäin ja niiden välinen tarkkuusero selvennetään.

Analyzing the Effects of the New Renewable Energy Policy in Russia on Investments into Wind, Solar and Small Hydro Power

This thesis presents an analysis of recently enacted Russian renewable energy policy based on capacity mechanism. Considering its novelty and poor coverage by academic literature, the aim of the thesis is to analyze capacity mechanism influence on investors’ decision-making process. The current research introduces a number of approaches to investment analysis. Firstly, classical financial model was built with Microsoft Excel® and crisp efficiency indicators such as net present value were determined. Secondly, sensitivity analysis was performed to understand different factors influence on project profitability. Thirdly, Datar-Mathews method was applied that by means of Monte Carlo simulation realized with Matlab Simulink®, disclosed all possible outcomes of investment project and enabled real option thinking. Fourthly, previous analysis was duplicated by fuzzy pay-off method with Microsoft Excel®. Finally, decision-making process under capacity mechanism was illustrated with decision tree. Capacity remuneration paid within 15 years is calculated individually for each RE project as variable annuity that guarantees a particular return on investment adjusted on changes in national interest rates. Analysis results indicate that capacity mechanism creates a real option to invest in renewable energy project by ensuring project profitability regardless of market conditions if project-internal factors are managed properly. The latter includes keeping capital expenditures within set limits, production performance higher than 75% of target indicators, and fulfilling localization requirement, implying producing equipment and services within the country. Occurrence of real option shapes decision-making process in the following way. Initially, investor should define appropriate location for a planned power plant where high production performance can be achieved, and lock in this location in case of competition. After, investor should wait until capital cost limit and localization requirement can be met, after that decision to invest can be made without any risk to project profitability. With respect to technology kind, investment into solar PV power plant is more attractive than into wind or small hydro power, since it has higher weighted net present value and lower standard deviation. However, it does not change decision-making strategy that remains the same for each technology type. Fuzzy pay-method proved its ability to disclose the same patterns of information as Monte Carlo simulation. Being effective in investment analysis under uncertainty and easy in use, it can be recommended as sufficient analytical tool to investors and researchers. Apart from described results, this thesis contributes to the academic literature by detailed description of capacity price calculation for renewable energy that was not available in English before. With respect to methodology novelty, such advanced approaches as Datar-Mathews method and fuzzy pay-off method are applied on the top of investment profitability model that incorporates capacity remuneration calculation as well. Comparison of effects of two different RE supporting schemes, namely Russian capacity mechanism and feed-in premium, contributes to policy comparative studies and exhibits useful inferences for researchers and policymakers. Limitations of this research are simplification of assumptions to country-average level that restricts our ability to analyze renewable energy investment region wise and existing limitation of the studying policy to the wholesale power market that leaves retail markets and remote areas without our attention, taking away medium and small investment into renewable energy from the research focus. Elimination of these limitations would allow creating the full picture of Russian renewable energy investment profile.

A semi-automatic computerized method to measure baroreflex-mediated heart rate responses that reduces interobserver variability

Arterial baroreflex sensitivity estimated by pharmacological impulse stimuli depends on intrinsic signal variability and usually a subjective choice of blood pressure (BP) and heart rate (HR) values. We propose a semi-automatic method to estimate cardiovascular reflex sensitivity to bolus infusions of phenylephrine and nitroprusside. Beat-to-beat BP and HR time series for male Wistar rats (N = 13) were obtained from the digitized signal (sample frequency = 2 kHz) and analyzed by the proposed method (PRM) developed in Matlab language. In the PRM, time series were low-pass filtered with zero-phase distortion (3rd order Butterworth used in the forward and reverse direction) and presented graphically, and parameters were selected interactively. Differences between basal mean values and peak BP (deltaBP) and HR (deltaHR) values after drug infusions were used to calculate baroreflex sensitivity indexes, defined as the deltaHR/deltaBP ratio. The PRM was compared to the method traditionally (TDM) employed by seven independent observers using files for reflex bradycardia (N = 43) and tachycardia (N = 61). Agreement was assessed by Bland and Altman plots. Dispersion among users, measured as the standard deviation, was higher for TDM for reflex bradycardia (0.60 ± 0.46 vs 0.21 ± 0.26 bpm/mmHg for PRM, P < 0.001) and tachycardia (0.83 ± 0.62 vs 0.28 ± 0.28 bpm/mmHg for PRM, P < 0.001). The advantage of the present method is related to its objectivity, since the routine automatically calculates the desired parameters according to previous software instructions. This is an objective, robust and easy-to-use tool for cardiovascular reflex studies.

Toward including the effect of manufacturing processes in the pre-estimated losses of the switched reluctance motor

The estimation of losses plays a key role in the process of building any electrical machine. How to estimate those losses while designing any machine; by obtaining the characteristic of the electrical steel from the catalogue and calculate the losses. However, this way is inaccurate since the electrical steel performs several manufacturing processes during the process of building any machine, which affects directly the magnetic property of the electrical steel and accordingly the characteristic of the electrical steel will be affected. That means the B–H curve of the steel that was obtained from the catalogue will be changed. Moreover, during loading and rotating the machine, some important changes occur to the B–H characteristic of the electrical steel such as the stress on the laminated iron. Accordingly, the pre-estimated losses are completely far from the actual losses because they were estimated based on the data of the electrical steel obtained from the catalogue. So in order to estimate the losses precisely significant factors of the manufacturing processes must be included. The paper introduces the systematic estimation of the losses including the effect of one of the manufacturing factors. Similarly, any other manufacturing factor can be included in the pre-designed losses estimations.

Auringon asemaa seuraavan aurinkopaneelijärjestelmän kehittäminen

Tässä diplomityössä käsitellään Aurinkoa seuraavan aurinkopaneelijärjestelmän kehittämistä Suomen olosuhteisiin. Työ on tehty osana Lappeenrannan teknillisen yliopiston Green Campus-projektia, jossa tarkoituksena on tuottaa yliopistolle energiaa uusiutuvilla menetelmillä ja käyttää niitä apuna tutkimuksessa sekä opetuksessa. Tavoitteena työssä on ymmärtää Auringon seuraamisen hyödyt sekä mahdolliset haitat aurinkopaneeli sovellutuksissa. Aikaisemman tutkimustiedon ollessa vähäistä, on työssä pyritty löytämään laskentamalli tuottavuuden laskentaan riippumatta siitä, missä päin maapalloa aurinkopaneelijärjestelmä sijaitsee. Työ alkaa kirjallisuustutkimuksella, jossa käydään läpi aurinkopaneelien toimintaperiaate, aurinkoenergian ja auringonpaistetuntien suuruusluokat Suomessa, sekä Suomen sääoloista johtuvat vaatimukset aurinkopaneelijärjestelmille. Tämän jälkeen on vertailtu kaupallisia järjestelmiä. Lopuksi järjestelmällistä tuotesuunnittelua hyväksikäyttäen suunnitellaan oma versio Aurinkoa seuraavasta aurinkopaneelijärjestelmästä. Oman suunnitelman järkevyyttä simuloidaan pienoismallilla Matlab-Arduino ympäristössä ja pyritään löytämään mahdollisia heikkouksia. Suomessa aurinkoenergiasta 90 % saadaan maalis – syyskuun aikana. Nykyisillä akkujärjestelmillä aurinkoenergia Suomessa ei sovellu kuin täydentäväksi energianlähteeksi. Aurinkoa seuraamalla voidaan saavuttaa 25-30 % tuottavuuden lisäys kesäaikana verrattuna staattiseen järjestelmään. Talvella hyöty tippuu 0-10 % luokkaan. Pienoismallilla simuloidut ohjaustavat osoittivat, että Aurinkoa on mahdollista seurata ilman sensoreita laskemalla Auringon paikka tähtitieteen kaavoista.

Glukoosi-fruktoosisiirapin kromatografinen fraktiointi

Kuluttajamarkkinoilla on suuri kysyntä vähäkalorisille ruuille ja juomille. Näitä tuotteita voidaan valmistaa mm. korvaamalla makeuttimena perinteisesti käytetty sakkaroosi fruktoosilla joka on 73% makeampi ja sisältää vähemmän kaloreita kuin sakkaroosi. Koska fruktoosi on makeampaa kuin sakkaroosi, sitä tarvitaan vähemmän elintarvikkeissa saman makeusasteen saavuttamiseksi. Glukoosi-fruktoosisiirappia valmistetaan maissista tai sakkaroosista hydrolyysin avulla. Teollisuudessa glukoosi-fruktoosisiirapin fraktiointi tehdään käyttämällä jatkuvatoimista simuloitu liikkuva peti-prosessia (Simulated Moving Bed, SMB-prosessi). Tässä työssä tutkitaan voidaanko kierrätyskromatografiaa (Steady State Recycling , SSR-prosessi) käyttäen fraktioida glukoosi-fruktoosisiirappia tehokkaasti. Fruktoosi erotetaan glukoosista käyttämällä erotusmateriaalina vahvoja kationinvaihtohartseja kalsium-muodossa. Työssä tehtiin kirjallisuusselvitys, jonka perusteella valittiin sopiva erotusmateriaali ja koeolosuhteet. Kokeellisessa osassa glukoosille ja fruktoosille määritettiin adsorptioisotermit Frontal analysis -menetelmällä. Kokeellisesti määritettyihin isotermeihin sovitetut mallit validoitiin ja SSR-prosessi suunniteltiin panoserotuskokeiden avulla. SSR-prosessia mallinnettiin käyttäen MATLAB-ohjelmaa.

Enzymatic Synthesis, Recovery, and Purification of Ethyl β-ᴅ-glucopyranoside

A method to synthesize ethyl β-ᴅ-glucopyranoside (BEG) was searched. Feasibility of different ion exchange resins was examined to purify the product from the synthetic binary solution of BEG and glucose. The target was to produce at least 50 grams of 99 % pure BEG with a scaled up process. Another target was to transfer the batch process into steady-state recycle chromatography process (SSR). BEG was synthesized enzymatically with reverse hydrolysis utilizing β-glucosidase as a catalyst. 65 % of glucose reacted with ethanol into BEG during the synthesis. Different ion exchanger based resins were examined to separate BEG from glucose. Based on batch chromatography experiments the best adsorbent was chosen between styrene based strong acid cation exchange resins (SAC) and acryl based weak acid cation exchange resins (WAC). CA10GC WAC resin in Na+ form was chosen for the further separation studies. To produce greater amounts of the product the batch process was scaled up. The adsorption isotherms for the components were linear. The target purity was possible to reach already in batch without recycle with flowrate and injection size small enough. 99 % pure product was produced with scaled-up batch process. Batch process was transferred to SSR process utilizing the data from design pulse chromatograms and Matlab simulations. The optimal operating conditions for the system were determined. Batch and SSR separation results were compared and by using SSR 98 % pure products were gained with 40 % higher productivity and 40 % lower eluent consumption compared to batch process producing as pure products.

Carbon footprint of recycled paper

Industrial production of pulp and paper is an intensive consumer of energy, natural resources, and chemicals that result in a big carbon footprint of the final product. At present companies and industries aspire to calculate their gas emissions into the atmosphere in order to afterwards reduce atmospheric contamination. One of the approaches allowing to increase carbon burden from the pulp and paper manufacture is paper recycling. The general purpose of the current paper is to establish methods of quantifying and minimizing the carbon footprint of paper. The first target of this research is to derive a mathematical relationship between virgin fibre requirements with respect to the amount of recycled paper used in the pulp. One more purpose is to establish a model to be used to clarify the contribution of recycling and transportation to decreasing carbon dioxide emissions. For this study sensitivity analysis is used to investigate the robustness of obtained results. The results of the present study show that an increasing of recycling rate does not always lead to minimizing the carbon footprint. Additionally, we derived that transportation of waste paper throughout distances longer than 5800 km has no sense because the use of that paper will only increase carbon dioxide emissions and it is better to reject recycling at all. Finally, we designed the model for organization of a new supply chain of paper product to a customer. The models were implemented as reusable MATLAB frameworks.

Akkuvarastot omakotitaloissa

Energian varastointi on noussut keskeiseksi energia-alan teemaksi viime vuosina. Erityi-sesti uusiutuvan tuotannon lisääntyminen ja energian käytön tehostaminen ovat edesautta-neet energiavarastoratkaisuiden mukaantuloa. Työssä tarkastellaan litiumrautafosfaattiak-kujen käytön kannattavuutta omakotitaloissa. Tavoitteena on selvittää, millä reunaehdoilla näiden akkujen käyttö energiavarastoina tulee kannattavaksi Lappeenrannan olosuhteissa. Kannattavuutta selvitetään litiumrautafosfaattiakkujen markkina-analyysin ja teknistalou-dellisen analyysin sekä Matlab-simulaation avulla.