Mekanismien dynamiikan simuloinnissa sovellettavia numeerisia- ja mallinnusmenetelmiä

Koneet voidaan usein jakaa osajärjestelmiin, joita ovat ohjaus- ja säätöjärjestelmät, voimaa tuottavat toimilaitteet ja voiman välittävät mekanismit. Eri osajärjestelmiä on simuloitu tietokoneavusteisesti jo usean vuosikymmenen ajan. Osajärjestelmien yhdistäminen on kuitenkin uudempi ilmiö. Usein esimerkiksi mekanismien mallinnuksessa toimilaitteen tuottama voimaon kuvattu vakiona, tai ajan funktiona muuttuvana voimana. Vastaavasti toimilaitteiden analysoinnissa mekanismin toimilaitteeseen välittämä kuormitus on kuvattu vakiovoimana, tai ajan funktiona työkiertoa kuvaavana kuormituksena. Kun osajärjestelmät on erotettu toisistaan, on niiden välistenvuorovaikutuksien tarkastelu erittäin epätarkkaa. Samoin osajärjestelmän vaikutuksen huomioiminen koko järjestelmän käyttäytymissä on hankalaa. Mekanismien dynamiikan mallinnukseen on kehitetty erityisesti tietokoneille soveltuvia numeerisia mallinnusmenetelmiä. Useimmat menetelmistä perustuvat Lagrangen menetelmään, joka mahdollistaa vapaasti valittaviin koordinaattimuuttujiin perustuvan mallinnuksen. Numeerista ratkaisun mahdollistamiseksi menetelmän avulla muodostettua differentiaali-algebraaliyhtälöryhmää joudutaan muokkaamaan esim. derivoimalla rajoiteyhtälöitä kahteen kertaan. Menetelmän alkuperäisessä numeerisissa ratkaisuissa kaikki mekanismia kuvaavat yleistetyt koordinaatit integroidaan jokaisella aika-askeleella. Tästä perusmenetelmästä johdetuissa menetelmissä riippumattomat yleistetyt koordinaatit joko integroidaan ja riippuvat koordinaatit ratkaistaan rajoiteyhtälöiden perusteella tai yhtälöryhmän kokoa pienennetään esim. käyttämällä nopeus- ja kiihtyvyysanalyyseissä eri kiertymäkoordinaatteja kuin asema-analyysissä. Useimmat integrointimenetelmät on alun perin tarkoitettu differentiaaliyhtälöiden (ODE) ratkaisuunjolloin yhtälöryhmään liitetyt niveliä kuvaavat algebraaliset rajoiteyhtälöt saattavat aiheuttaa ongelmia. Nivelrajoitteiden virheiden korjaus, stabilointi, on erittäin tärkeää mekanismien dynamiikan simuloinnin onnistumisen ja tulosten oikeellisuuden kannalta. Mallinnusmenetelmien johtamisessa käytetyn virtuaalisen työn periaatteen oletuksena nimittäin on, etteivät rajoitevoimat tee työtä, eli rajoitteiden vastaista siirtymää ei tapahdu. Varsinkaan monimutkaisten järjestelmien pidemmissä analyyseissä nivelrajoitteet eivät toteudu tarkasti. Tällöin järjestelmän energiatasapainoei toteudu ja järjestelmään muodostuu virtuaalista energiaa, joka rikkoo virtuaalisen työn periaatetta, Tästä syystä tulokset eivät enää pidäpaikkaansa. Tässä raportissa tarkastellaan erityyppisiä mallinnus- ja ratkaisumenetelmiä, ja vertaillaan niiden toimivuutta yksinkertaisten mekanismien numeerisessa ratkaisussa. Menetelmien toimivuutta tarkastellaan ratkaisun tehokkuuden, nivelrajoitteiden toteutumisen ja energiatasapainon säilymisen kannalta.

Validation Studies of Thermal-hydraulic Code for Safety Analysis of Nuclear Power Plants

This thesis gives an overview of the validation process for thermal hydraulic system codes and it presents in more detail the assessment and validation of the French code CATHARE for VVER calculations. Three assessment cases are presented: loop seal clearing, core reflooding and flow in a horizontal steam generator. The experience gained during these assessment and validation calculations has been used to analyze the behavior of the horizontal steam generator and the natural circulation in the geometry of the Loviisa nuclear power plant. The cases presented are not exhaustive, but they give a good overview of the work performed by the personnel of Lappeenranta University of Technology (LUT). Large part of the work has been performed in co-operation with the CATHARE-team in Grenoble, France. The design of a Russian type pressurized water reactor, VVER, differs from that of a Western-type PWR. Most of thermal-hydraulic system codes are validated only for the Western-type PWRs. Thus, the codes should be assessed and validated also for VVER design in order to establish any weaknesses in the models. This information is needed before codes can be used for the safety analysis. Theresults of the assessment and validation calculations presented here show that the CATHARE code can be used also for the thermal-hydraulic safety studies for VVER type plants. However, some areas have been indicated which need to be reassessed after further experimental data become available. These areas are mostly connected to the horizontal stem generators, like condensation and phase separation in primary side tubes. The work presented in this thesis covers a large numberof the phenomena included in the CSNI code validation matrices for small and intermediate leaks and for transients. Also some of the phenomena included in the matrix for large break LOCAs are covered. The matrices for code validation for VVER applications should be used when future experimental programs are planned for code validation.

Multiplicação in vitro dos porta-enxertos de Prunus sp. 'Barrier' e 'Cadaman'

Este trabalho foi realizado com o objetivo de estabelecer a melhor concentração de sais do meio MS e da citocinina BAP para a multiplicação dos porta-enxertos de Prunus sp. 'Barrier' e 'Cadaman'. Segmentos nodais foram introduzidos em tubos de ensaio contendo 10 mL de meio de cultura com variações na concentração de sais (MS; ½MS; e 2/3MS) combinadas com cinco concentrações de BAP (0; 1,5; 2,5; 3,5 e 4,5 miM). Utilizou-se um fatorial 2x3x5, distribuído em blocos casualizados, compostos por quatro repetições contendo cinco tubos de ensaio cada uma, sendo inoculado um segmento nodal por tubo. As avaliações foram realizadas após cinco semanas de cultivo em ambiente com intensidade luminosa de 20 miE m-2 s-1, fotoperíodo de 16 horas e temperatura de 24 ± 4ºC. Verificou-se maior número médio de gemas e de brotações para a cultivar Barrier. À medida que se reduziu a concentração de sais do meio de cultura, obteve-se maior número de brotações, porém com menor tamanho. As regressões polinomiais das variáveis número de gemas, brotações por explante e comprimento das brotações apresentaram um ajustamento quadrático para níveis de BAP, atingindo os pontos de máximo 31,2 gemas/explante; 4,6 brotações por explante, e 8,1 mm de comprimento nas concentrações 3,3; 3,1, e 3,1 miM de BAP, respectivamente.

Dynamic Model of a Small Once-Through Boiler

In this study, a model for the unsteady dynamic behaviour of a once-through counter flow boiler that uses an organic working fluid is presented. The boiler is a compact waste-heat boiler without a furnace and it has a preheater, a vaporiser and a superheater. The relative lengths of the boiler parts vary with the operating conditions since they are all parts of a single tube. The present research is a part of a study on the unsteady dynamics of an organic Rankine cycle power plant and it will be a part of a dynamic process model. The boiler model is presented using a selected example case that uses toluene as the process fluid and flue gas from natural gas combustion as the heat source. The dynamic behaviour of the boiler means transition from the steady initial state towards another steady state that corresponds to the changed process conditions. The solution method chosen was to find such a pressure of the process fluid that the mass of the process fluid in the boiler equals the mass calculated using the mass flows into and out of the boiler during a time step, using the finite difference method. A special method of fast calculation of the thermal properties has been used, because most of the calculation time is spent in calculating the fluid properties. The boiler was divided into elements. The values of the thermodynamic properties and mass flows were calculated in the nodes that connect the elements. Dynamic behaviour was limited to the process fluid and tube wall, and the heat source was regarded as to be steady. The elements that connect the preheater to thevaporiser and the vaporiser to the superheater were treated in a special way that takes into account a flexible change from one part to the other. The model consists of the calculation of the steady state initial distribution of the variables in the nodes, and the calculation of these nodal values in a dynamic state. The initial state of the boiler was received from a steady process model that isnot a part of the boiler model. The known boundary values that may vary during the dynamic calculation were the inlet temperature and mass flow rates of both the heat source and the process fluid. A brief examination of the oscillation around a steady state, the so-called Ledinegg instability, was done. This examination showed that the pressure drop in the boiler is a third degree polynomial of the mass flow rate, and the stability criterion is a second degree polynomial of the enthalpy change in the preheater. The numerical examination showed that oscillations did not exist in the example case. The dynamic boiler model was analysed for linear and step changes of the entering fluid temperatures and flow rates.The problem for verifying the correctness of the achieved results was that there was no possibility o compare them with measurements. This is why the only way was to determine whether the obtained results were intuitively reasonable and the results changed logically when the boundary conditions were changed. The numerical stability was checked in a test run in which there was no change in input values. The differences compared with the initial values were so small that the effects of numerical oscillations were negligible. The heat source side tests showed that the model gives results that are logical in the directions of the changes, and the order of magnitude of the timescale of changes is also as expected. The results of the tests on the process fluid side showed that the model gives reasonable results both on the temperature changes that cause small alterations in the process state and on mass flow rate changes causing very great alterations. The test runs showed that the dynamic model has no problems in calculating cases in which temperature of the entering heat source suddenly goes below that of the tube wall or the process fluid.

Estabelecimento in vitro de porta-enxertos de Prunus sp.

O objetivo deste trabalho foi avaliar o estabelecimento in vitro de cinco cultivares de porta-enxerto de Prunus sp. em variações de dois meios de cultura. O material vegetal foi obtido de plantas-matrizes cultivadas em telado, sendo a desinfestação realizada em soluções à base de álcool e hipoclorito de sódio. Utilizou-se um esquema fatorial 5 x 2 (cultivares Aldrighi, Barrier, Capdeboscq, Flordaguard e GF677 x formulações salinas ½MS e QL), distribuído em blocos casualizados, compostos por quatro repetições, com cinco tubos de ensaio cada uma, sendo inoculado um segmento nodal por tubo. O cultivo dos explantes foi realizado por quatro semanas, em ambiente com 20 µEm-2s-1, 24 ± 4ºC e fotoperíodo de 16 horas. Independentemente do meio de cultivo, os porta-enxertos 'Barrier', 'Capdebosq' e 'Flordaguard' apresentaram maiores porcentagem de explantes estabelecidos e número de brotações, sendo os mais responsivos in vitro. A formulação salina ½MS foi a mais indicada para o estabelecimento da cultivar GF677, enquanto a QL para a 'Aldrighi'. Não foram observadas plântulas com sintomas de clorose, vitrificação ou encarquilhamento, indicando que as condições de cultura utilizadas foram satisfatórias. As maiores porcentagens de contaminação foram verificadas nas cultivares Aldrighi (72%) e Capdeboscq (41%), as quais são as mais utilizadas na região Sul do Brasil.

Replicabilidad de los factores de segundo orden del 16PF-5 en muestras americanas y españolas

El objetivo de este estudio fue replicar las estructuras de 5 y 6 factores de segundo orden del 16PF-5. Para la estructura de 5 factores se toma como referencia teórica la estructura obtenida por Russell y Karol (1995), y para la estructura de 6 factores (incluyendo un factor adicional de Razonamiento) la obtenida en muestras americanas por Cattell y Cattell (1995). Se utilizan tres procedimientos para el estudio de la replicabilidad, a) análisis factorial exploratorio, b) análisis de la estructura ortogonal Procrustes, y c) análisis de los índices de congruencia entre las tres matrices factoriales. Las matrices factoriales obtenidas en el presente estudio son similares a las informadas en los estudios de referencia, aunque la solución Procrustes se revela ligeramente más parecida. Los índices de congruencia en ge- neral son aceptables por lo que se concluye que el 16PF-5 demuestra buena replicabilidad.

Clonal growth and demography of a hemicryptophyte alpine plant: Leontopodium alpinum Cassini

Clonally reproducing hemicryptophytic rosette plants are common in the alpine belt. However, their demography, and indirectly their growth and reproductive strategy in these harsh conditions, was rarely studied. We analysed the morphology, clonal reproduction and demography of one such species, Leontopodium alpinum, in two populations of the Swiss Alps. The species forms small colonies of 1-5 (maximum 30) sterile rosettes with a few flowering stalks. After flowering, the apical meristem dies and one or two new axillary buds grow below the previous rosette in the following year, developing into short rhizomes (<2 cm), which decay after four years. The new stalk produces sterile rosettes before flowering after two to four years, depending on climatic conditions. The apical meristem often dies during the sterile stage, and is replaced by a new axillary bud. Levkovitch matrices on two stages (sterile and flowering rosettes) showed that rosette survival and clonal reproduction maintain long-lived populations (λ = 0.96). Elasticities indicated that a change in the survival of sterile rosettes had the strongest effect on population dynamics, and this stage lasts, on average, 6.8 years at 2480 m. Altogether, L. alpinum is following Tomlinson's architectural model. This growth form appears perfectly adapted to harsh alpine conditions: the clonal ramification ensures longevity to genets and the semelparous behaviour of the rosettes allows an efficient flowering, whatever the climatic conditions. L. alpinum appears to follow a common growth model among rosette possessing hemicryptophytes in the alpine belt.

Spatially varying selection shapes life history clines among populations of Drosophila melanogaster from sub-Saharan Africa.

Clines in life history traits, presumably driven by spatially varying selection, are widespread. Major latitudinal clines have been observed, for example, in Drosophila melanogaster, an ancestrally tropical insect from Africa that has colonized temperate habitats on multiple continents. Yet, how geographic factors other than latitude, such as altitude or longitude, affect life history in this species remains poorly understood. Moreover, most previous work has been performed on derived European, American and Australian populations, but whether life history also varies predictably with geography in the ancestral Afro-tropical range has not been investigated systematically. Here, we have examined life history variation among populations of D. melanogaster from sub-Saharan Africa. Viability and reproductive diapause did not vary with geography, but body size increased with altitude, latitude and longitude. Early fecundity covaried positively with altitude and latitude, whereas lifespan showed the opposite trend. Examination of genetic variance-covariance matrices revealed geographic differentiation also in trade-off structure, and QST -FST analysis showed that life history differentiation among populations is likely shaped by selection. Together, our results suggest that geographic and/or climatic factors drive adaptive phenotypic differentiation among ancestral African populations and confirm the widely held notion that latitude and altitude represent parallel gradients.