Second differential of the entropy as a criterion for the stability in low-dimensional climate models

The second differential of the entropy is used for analysing the stability of a thermodynamic climatic model. A delay time for the heat flux is introduced whereby it becomes an independent variable. Two different expressions for the second differential of the entropy are used: one follows classical irreversible thermodynamics theory; the second is related to the introduction of response time and is due to the extended irreversible thermodynamics theory. the second differential of the classical entropy leads to unstable solutions for high values of delay times. the extended expression always implies stable states for an ice-free earth. When the ice-albedo feedback is included, a discontinuous distribution of stable states is found for high response times. Following the thermodynamic analysis of the model, the maximum rates of entropy production at the steady state are obtained. A latitudinally isothermal earth produces the extremum in global entropy production. the material contribution to entropy production (by which we mean the production of entropy by material transport of heat) is a maximum when the latitudinal distribution of temperatures becomes less homogeneous than present values

New experimental apparatus for the study of the Bénard-Rayleigh problem

A new experimental system to measure the equivalent thermal conductivity of a liquid with regard to the Bénard-Rayleigh problem was constructed. The liquid is enclosed within walls of polymethylmethacrylate between two copper plates in which there are thermocouples to measure the difference in temperature between the lower and upper surfaces of the layer of liquid. Heat flux is measured by means of a linear heat fluxmeter consisting of 204 thermocouples in series. The fluxmeter was calibrated and the linear relationship that exists between the heat flux and the emf generated was verified. The thermal conductivity of the polymethylmethacrylate employed was measured and measurements of the equivalent conductivity in cylindrical boundaries of two silicone oils were made. The critical value of the temperature difference and the contribution of the convective process to the transmission of heat were determined.

Energy balance measurements over a banana orchard in the Semiarid region in the Northeast of Brazil

The objective of this work was to evaluate the reliability of eddy covariance measurements, analyzing the energy balance components, evapotranspiration and energy balance closure in dry and wet growing seasons, in a banana orchard. The experiment was carried out at a farm located within the irrigation district of Quixeré, in the Lower Jaguaribe basin, in Ceará state, Brazil. An eddy covariance system was used to measure the turbulent flux. An automatic weather station was installed in a grass field to obtain the reference evapotranspiration (ET0) from the combined FAO-Penman-Monteith method. Wind speed and vapor pressure deficit are the most important variables on the evaporative process in both growing seasons. In the dry season, the heat fluxes have a similar order of magnitude, and during the wet season the latent heat flux is the largest. The eddy covariance system had acceptable reliability in measuring heat flux, with actual evapotranspiration results comparing well with those obtained by using the water balance method. The energy balance closure had good results for the study area, with mean values of 0.93 and 0.86 for the dry and wet growing seasons respectively.

Analysis of thermal properties of phase change materials (PCM) using differential scanning calorimeter (DSC)

Commercially available PCM RT 20, RT 27, SP 22, A17 and SP 25 A8. Were analyzed using dynamic and step method of heat flux DSC. The results of the dinamic and step method were compared with commercial valures. It was found that RT 20 & RT 27 showed good conforming of results with commercial values while SP 22 A17 & SP 25 A8 did not show conformity.

Parametrización de kB-1 para estimar el flujo de calor sensible con imágenes AVHRR

La necesidad de evaluar la evapotranspiración a escala regional para la gestión de regadíos ha hecho que sean innumerables los intentos por aplicar imágenes AVHRR-NOAA en la determinación el flujo de calor sensible. La principal limitación de estos métodos es la estimación de la resistencia aerodinámica. El parámetro crítico en la expresión de la resistencia aerodinámica es kB-1. La parametrización de kB-1 ha sido infructuosa a escala regional por no disponer hasta ahora de medidas de flujo de calor sensible a escala del píxel AVHRR en superficies heterogéneas y durante toda una temporada de riegos. Para resolver esta medida de flujo se ha desarrollado el cintilómetro. En la primera parte de este trabajo se estudia la representatividad espacial de las medidas del cintilómetro. El núcleo de esta aportación consiste en la correlación entre el parámetro kB-1, el NDVI y la altura solar. Los buenos resultados obtenidos (r2=0.81) ofrecen una nueva metodología para determinar el flujo de calor sensible. La estimación de kB-1, las imágenes AVHRR y los datos meteorológicos permiten calcular el flujo de calor sensible durante toda la temporada de riegos con errores inferiores al 20%.

Lämmönsiirron numeerinen laskenta profiloidulle kanavalle.

Työn tavoitteena oli tutkia aaltomaisen profiloinnin vaikutusta suorien jäähdytyskanavien lämmönsiirtoon ja painehäviöön. Erilaisia profiileja oli kymmenen kappaletta ja ne olivat 5mm leveitä ja 30cm pitkiä kukin. Ne laskettiin kolmeulotteisina tapauksina FINFLO-virtausratkaisijalla kolmella eri Reynoldsin luvulla, jotka vastasivat laminaarista, osittain turbulenttista ja lähes kokonaan turbulenttista virtausta. Lämmönsiirtoaine oli kuiva +30°C ilma ja profiloinnin toteutustapa oli toisiaan sivuavat ympyräkaaret kolmella erilaisella säteen arvolla ja kolmella erilaisella aallonpituuden arvolla. Lisäksi laskettiin saman levyisen tasokanavan arvot jokaisella Reynoldsin luvulla kaksiulotteisina tapauksina. Näitä profiloimattomia kanavia pidettiin referenssitapauksina. Tuloksena havaittiin että profiloimalla saadaan yksiselitteisesti suurempi lämpöteho ulos samasta tilavuudesta. Lämmönsiirtokerroin kasvaa profiloinnin avulla parhaimmillaan n. 20% käytetystä turbulenssimallista tai lämmönsiirtokertoimen määritelmästä riippumatta. Painehäviö kasvaa myös aina, mutta kitkakerroin voi hieman pienentyä. Profiilin varsinaisena hyvyyskriteerinä pidettiin lämmönsiirtokertoimen ja kitkakertoimen suhdetta h/f. Se osoittautui riippuvaksi Reynoldsin luvusta ja turbulenssimallista; ASM ja Chien k-έ -mallit ennustavat transitioetäisyyden eri tavalla. Laminaarisilla virtauksilla h/f :n vaihtelu oli vähäistä; suhde vaihteli vain ±5% eri profiilien kesken. ASM-mallilla havaittiin sekundääripyörteilyä, ehkä siksi että se mallintaa anisotrooppisen turbulenssin. Chien k-έ malli ennusti suuremman ja aikaisemmin alkavan turbulenttisuuden kuin ASM. Lisäksi havaittiin mm. että tietyillä profiileilla muodostuu kanavan kapeimpaan kohtaan selvä virtausnopeuden paikallinen minimi seinämän läheisyyden takia.

Arinakattilan lämmönsiirron laskentamallin kehittäminen

Kattilalaitosten polttoaineen syötössä ilmenevät häiriöt ja biopolttoaineiden laatuvaihtelut aiheuttavat epävakaata palamista ja tekevät prosessin hallinnasta vaikeampaa. Polttoaineen laatuvaihtelut vaikuttavat koko prosessiin ja näkyvät lopulta myös höyryntuotannossa. Kompensoinnilla pyritään estämään häiriöiden suuret vaikutukset höyryn tuotantoon. Tarkoituksena on saada kattilan toiminta ja tehon tuotanto tasaisemmaksi ja helpommin hallittavaksi. Tässä diplomityössä tarkastellaan polttoaineen ominaisuuksien, säteilylämmönsiirron sekä säätöjen vaikutusta toisiinsa ja merkitystä kattilan toiminnan kannalta. Työssä muodostetaan säteilylämmönsiirron laskentamalli arinakattilan tulipesälle käyttäen hyväksi hyvin sekoittuneen tulipesän menetelmää. Menetelmällä voidaan määrittääsavukaasujen keskimääräinen lämpötila tulipesässä, lämpövirta tulipesän seiniin tai poltossa vapautuva lämpöteho. Mallin avulla voidaan paremmin ymmärtää prosessin käyttäytymistä polttoaineen laadun muuttuessa sekä helpottaa ja nopeuttaa kattilan käyttäytymisen ennustamista. Laskentamalli tehtiin Excel –laskentaohjelmaan, jossa se testataan. Verifioinnin jälkeen malli on tarkoitus siirtää toimimaan apros –simulointiympäristöön.

Modeling of water/steam circulation in Circulating Fluidized Bed boiler

A distinctive design feature of steam boiler with natural circulation is the presence of the steam drum which plays a role of the separator of vapor from the flow of water-and-steam mixture coming into steam drum from the furnace tubes. Steam drum with unheated downcomer tubes, deducing from it, and riser (screen/furnace tubes) inside the furnace is a closed circulation loop in which movement of water (downcomer tubes) and water-and-steam mixture (riser tubes) is organized. The movement of the working fluid is appears due to occurrence of the natural pressure, determined by the difference in hydrostatic pressure and the mass of water and water-and-steam mixtures in downcomer and riser tubes and called the driving pressure of the natural circulation: S drive = H steam (ρ down + ρ mix) g where: ρ down - density of water in downcomer tubes; ρ mix - density of water in riser tubes; H steam - height of steam content section; g - acceleration of gravity. In steam boilers with natural circulation the circulation rate is usually between 10 and 30. Thus, consumption of water in the circulation circuit “circulation rate times” more than steam output of the boiler. There are two aspects of the design of natural water circulation loops. One is to ensure a sufficient mass flux of circulating water to avoid burnout of evaporator tubes. The other is to avoid tube wall temperature fluctuation and tube vibration due to oscillation of circulation velocity. The design criteria are therefore reduced, in principle, to those of critical heat flux, critical flow rate for burnout, and flow instability. In practical design, however, the circulation velocity and the void fraction at the evaporator tube outlet are used as the design criteria (Seikan I., et. al., 1999). This study has been made with assumption that the heat flux in the furnace of the boiler even all the time. The target of the study was to define the circulation rate of the boiler, thus average heat flux do not change it. I would like to acknowledge professionals from “Foster Wheeler” company for good and comfortable cooperation.

Modeling of residual stresses and distortion due to welding in fillet welds

Welding has a growing role in modern world manufacturing. Welding joints are extensively used from pipes to aerospace industries. Prediction of welding residual stresses and distortions is necessary for accurate evaluation of fillet welds in relation to design and safety conditions. Residual stresses may be beneficial or detrimental, depending whether they are tensile or compressive and the loading. They directly affect the fatigue life of the weld by impacting crack growth rate. Beside theoretical background of residual stresses this study calculates residual stresses and deformations due to localized heating by welding process and subsequent rapid cooling in fillet welds. Validated methods are required for this purpose due to complexity of process, localized heating, temperature dependence of material properties and heat source. In this research both empirical and simulation methods were used for the analysis of welded joints. Finite element simulation has become a popular tool of prediction of welding residual stresses and distortion. Three different cases with and without preload have been modeled during this study. Thermal heat load set is used by calculating heat flux from the given heat input energy. First the linear and then nonlinear material behavior model is modeled for calculation of residual stresses. Experimental work is done to calculate the stresses empirically. The results from both the methods are compared to check their reliability. Residual stresses can have a significant effect on fatigue performance of the welded joints made of high strength steel. Both initial residual stress state and subsequent residual stress relaxation need to be considered for accurate description of fatigue behavior. Tensile residual stresses are detrimental and will reduce the fatigue life and compressive residual stresses will increase it. The residual stresses follow the yield strength of base or filler material and the components made of high strength steel are typically thin, where the role of distortion is emphasizing.

Two-phase flow measurements with spacer grid in a rod bundle geometry

Heat transfer effectiveness in nuclear rod bundles is of great importance to nuclear reactor safety and economics. An important design parameter is the Critical Heat Flux (CHF), which limits the transferred heat from the fuel to the coolant. The CHF is determined by flow behaviour, especially the turbulence created inside the fuel rod bundle. Adiabatic experiments can be used to characterize the flow behaviour separately from the heat transfer phenomena in diabatic flow. To enhance the turbulence, mixing vanes are attached to spacer grids, which hold the rods in place. The vanes either make the flow swirl around a single sub-channel or induce cross-mixing between adjacent sub-channels. In adiabatic two-phase conditions an important phenomenon that can be investigated is the effect of the spacer on canceling the lift force, which collects the small bubbles to the rod surfaces leading to decreased CHF in diabatic conditions and thus limits the reactor power. Computational Fluid Dynamics (CFD) can be used to simulate the flow numerically and to test how different spacer configurations affect the flow. Experimental data is needed to validate and verify the used CFD models. Especially the modeling of turbulence is challenging even for single-phase flow inside the complex sub-channel geometry. In two-phase flow other factors such as bubble dynamics further complicate the modeling. To investigate the spacer grid effect on two-phase flow, and to provide further experimental data for CFD validation, a series of experiments was run on an adiabatic sub-channel flow loop using a duct-type spacer grid with different configurations. Utilizing the wire-mesh sensor technology, the facility gives high resolution experimental data in both time and space. The experimental results indicate that the duct-type spacer grid is less effective in canceling the lift force effect than the egg-crate type spacer tested earlier.

Lämpövuoanturin mittauselektroniikan suunnittelu ja toteutus

Energian kulutuksen vähentäminen ja sen tutkiminen on kasvavan kiinnostuksen kohteena. Syntyneen lämmön mittaaminen on yksi tapa mitata energian siirtymistä. Lämpötilan mittaaminen on yleistä, vaikka usein on merkittävämpää selvittää missä ja miten lämpöenergia on siirtynyt. Tästä syystä tarvitaan lämpövuoantureita, jotka reagoivat suoraan lämpövuohon eli lämpöenergian siirtymiseen. Tässä tutkimuksessa suunnitellaan ja toteutetaan lämpövuoanturin mittauselektroniikka vaativaan käyttöympäristöön. Työssä käytettävän gradienttilämpövuoanturin tuottama jännitesignaali on mikrovolttiluokkaa ja ympäristön aiheuttama kohina voi olla huomattavasti suurempi. Tämän takia anturin tuottamaa signaalia on vahvistettava, jotta sitä voidaan mitata luotettavasti. Tutkimuksessa keskitytään vahvistimen suunnitteluun, mutta suunnittelussa on otettava huomioon koko järjestelmä. Anturin sähköiset ominaisuudet ja ympäristö asettavat rajoitteita vahvistimelle. Tavoitteena on selvittää miten voidaan mitata mikrovolttien jännitesignaalia mahdollisimman suurella taajuuskaistalla vaativassa käyttöympäristössä. Työn tuloksena syntyi mittalaite, jota voidaan käyttää vaativassa ympäristössä lämpövuon mittaamiseen. Suunnitteluparametrien mukainen vahvistus ja päästökaista sekä offset-jännitteen ryömintä saavutettiin suunnitellulla mittalaitteella, mutta offsetjännite ja kohina olivat hieman suunniteltua suuremmat. Mittalaitteella ja lämpövuoanturilla havaittiin selvästi lämpövuon muutoksia keinotekoisilla herätteillä.

Computational studies of deformation in stainless steel rings due to torch heating

A support ring of AISI 304L stainless steel that holds vertical, parallel wires arranged in a circle forming a cylinder is studied. The wires are attached to the ring with heat-induced shrinkage. When the ring is heated with a torch the heat affected zone tries to expand while the adjacent cool structure obstructs the expansion causing upsetting. During cooling, the ring shrinks smaller than its original size clamping the wires. The most important requirement for the ring is that it should be as round as possible and the deformations should occur as overall shrinkage in the ring diameter. A three-dimensional nonlinear transient sequential thermo-structural Abaqus model is used together with a Fortran code that enters the heat flux to each affected element. The local and overall deformations in one ring inflicted by the heating are studied with a small amount of inspection on residual stresses. A variety of different cases are chosen to be studied with the model constructed to provide directional knowledge; torch flux with the means of speed, location of the wires, heating location and structural factors. The decrease of heating speed increases heat flux that rises the temperature increasing shrinkage. In a single progressive heating uneven distribution of shrinkage appears to the start/end region that can be partially fixed with using speeded heating’s to strengthen the heating of that region. Location of the wires affect greatly to the caused shrinkage unlike heating location. The ring structure affects also greatly to the shrinkage; smaller diameter, bigger ring height, thinner thickness and greater number of wires increase shrinkage.

Electromagnetic and thermal design of a multilevel converter with high power density and reliability

Electric energy demand has been growing constantly as the global population increases. To avoid electric energy shortage, renewable energy sources and energy conservation are emphasized all over the world. The role of power electronics in energy saving and development of renewable energy systems is significant. Power electronics is applied in wind, solar, fuel cell, and micro turbine energy systems for the energy conversion and control. The use of power electronics introduces an energy saving potential in such applications as motors, lighting, home appliances, and consumer electronics. Despite the advantages of power converters, their penetration into the market requires that they have a set of characteristics such as high reliability and power density, cost effectiveness, and low weight, which are dictated by the emerging applications. In association with the increasing requirements, the design of the power converter is becoming more complicated, and thus, a multidisciplinary approach to the modelling of the converter is required. In this doctoral dissertation, methods and models are developed for the design of a multilevel power converter and the analysis of the related electromagnetic, thermal, and reliability issues. The focus is on the design of the main circuit. The electromagnetic model of the laminated busbar system and the IGBT modules is established with the aim of minimizing the stray inductance of the commutation loops that degrade the converter power capability. The circular busbar system is proposed to achieve equal current sharing among parallel-connected devices and implemented in the non-destructive test set-up. In addition to the electromagnetic model, a thermal model of the laminated busbar system is developed based on a lumped parameter thermal model. The temperature and temperature-dependent power losses of the busbars are estimated by the proposed algorithm. The Joule losses produced by non-sinusoidal currents flowing through the busbars in the converter are estimated taking into account the skin and proximity effects, which have a strong influence on the AC resistance of the busbars. The lifetime estimation algorithm was implemented to investigate the influence of the cooling solution on the reliability of the IGBT modules. As efficient cooling solutions have a low thermal inertia, they cause excessive temperature cycling of the IGBTs. Thus, a reliability analysis is required when selecting the cooling solutions for a particular application. The control of the cooling solution based on the use of a heat flux sensor is proposed to reduce the amplitude of the temperature cycles. The developed methods and models are verified experimentally by a laboratory prototype.

Coolability of porous core debris beds : Effects of bed geometry and multi-dimensional flooding

This thesis addresses the coolability of porous debris beds in the context of severe accident management of nuclear power reactors. In a hypothetical severe accident at a Nordic-type boiling water reactor, the lower drywell of the containment is flooded, for the purpose of cooling the core melt discharged from the reactor pressure vessel in a water pool. The melt is fragmented and solidified in the pool, ultimately forming a porous debris bed that generates decay heat. The properties of the bed determine the limiting value for the heat flux that can be removed from the debris to the surrounding water without the risk of re-melting. The coolability of porous debris beds has been investigated experimentally by measuring the dryout power in electrically heated test beds that have different geometries. The geometries represent the debris bed shapes that may form in an accident scenario. The focus is especially on heap-like, realistic geometries which facilitate the multi-dimensional infiltration (flooding) of coolant into the bed. Spherical and irregular particles have been used to simulate the debris. The experiments have been modeled using 2D and 3D simulation codes applicable to fluid flow and heat transfer in porous media. Based on the experimental and simulation results, an interpretation of the dryout behavior in complex debris bed geometries is presented, and the validity of the codes and models for dryout predictions is evaluated. According to the experimental and simulation results, the coolability of the debris bed depends on both the flooding mode and the height of the bed. In the experiments, it was found that multi-dimensional flooding increases the dryout heat flux and coolability in a heap-shaped debris bed by 47–58% compared to the dryout heat flux of a classical, top-flooded bed of the same height. However, heap-like beds are higher than flat, top-flooded beds, which results in the formation of larger steam flux at the top of the bed. This counteracts the effect of the multi-dimensional flooding. Based on the measured dryout heat fluxes, the maximum height of a heap-like bed can only be about 1.5 times the height of a top-flooded, cylindrical bed in order to preserve the direct benefit from the multi-dimensional flooding. In addition, studies were conducted to evaluate the hydrodynamically representative effective particle diameter, which is applied in simulation models to describe debris beds that consist of irregular particles with considerable size variation. The results suggest that the effective diameter is small, closest to the mean diameter based on the number or length of particles.

Evaporative respiratory cooling augments pit organ thermal detection in rattlesnakes

Rattlesnakes use their facial pit organs to sense external thermal fluctuations. A temperature decrease in the heat-sensing membrane of the pit organ has the potential to enhance heat flux between their endothermic prey and the thermal sensors, affect the optimal functioning of thermal sensors in the pit membrane and reduce the formation of thermal ‘‘afterimages’’, improving thermal detection. We examined the potential for respiratory cooling to improve strike behaviour, capture, and consumption of endothermic prey in the South American rattlesnake, as behavioural indicators of thermal detection. Snakes with a higher degree of rostral cooling were more accurate during the strike, attacking warmer regions of their prey, and relocated and consumed their prey faster. These findings reveal that by cooling their pit organs, rattlesnakes increase their ability to detect endothermic prey; disabling the pit organs caused these differences to disappear. Rattlesnakes also modify the degree of rostral cooling by altering their breathing pattern in response to biologically relevant stimuli, such as a mouse odour. Our findings reveal that low humidity increases their ability to detect endothermic prey, suggesting that habitat and ambush sites election in the wild may be influenced by external humidity levels as well as temperature.