On the use of spatially discrete data to compute energy and mass balance

In many practical applications the state of field soils is monitored by recording the evolution of temperature and soil moisture at discrete depths. We theoretically investigate the systematic errors that arise when mass and energy balances are computed directly from these measurements. We show that, even with no measurement or model errors, large residuals might result when finite difference approximations are used to compute fluxes and storage term. To calculate the limits set by the use of spatially discrete measurements on the accuracy of balance closure, we derive an analytical solution to estimate the residual on the basis of the two key parameters: the penetration depth and the distance between the measurements. When the thickness of the control layer for which the balance is computed is comparable to the penetration depth of the forcing (which depends on the thermal diffusivity and on the forcing period) large residuals arise. The residual is also very sensitive to the distance between the measurements, which requires accurately controlling the position of the sensors in field experiments. We also demonstrate that, for the same experimental setup, mass residuals are sensitively larger than the energy residuals due to the nonlinearity of the moisture transport equation. Our analysis suggests that a careful assessment of the systematic mass error introduced by the use of spatially discrete data is required before using fluxes and residuals computed directly from field measurements.

Estimation of soil moisture in the root-zone from remote sensing data

Field-based soil moisture measurements are cumbersome. Thus, remote sensing techniques are needed because allows field and landscape-scale mapping of soil moisture depth-averaged through the root zone of existing vegetation. The objective of the study was to evaluate the accuracy of an empirical relationship to calculate soil moisture from remote sensing data of irrigated soils of the Apodi Plateau, in the Brazilian semiarid region. The empirical relationship had previously been tested for irrigated soils in Mexico, Egypt, and Pakistan, with promising results. In this study, the relationship was evaluated from experimental data collected from a cotton field. The experiment was carried out in an area of 5 ha with irrigated cotton. The energy balance and evaporative fraction (Λ) were measured by the Bowen ratio method. Soil moisture (θ) data were collected using a PR2 - Profile Probe (Delta-T Devices Ltd). The empirical relationship was tested using experimentally collected Λ and θ values and was applied using the Λ values obtained from the Surface Energy Balance Algorithm for Land (SEBAL) and three TM - Landsat 5 images. There was a close correlation between measured and estimated θ values (p<0.05, R² = 0.84) and there were no significant differences according to the Student t-test (p<0.01). The statistical analyses showed that the empirical relationship can be applied to estimate the root-zone soil moisture of irrigated soils, i.e. when the evaporative fraction is greater than 0.45.

Radiation and energy balance of lettuce culture inside a polyethylene greenhouse

The objective of this paper was to describe the radiation and energy balance, during the lettuce (Lactuca sativa, L. cv. Verônica) crop cycle inside a polyethylene greenhouse. The radiation and energy balance was made inside a tunnel greenhouse with polyethylene cover (100 mum) and in an external area, both areas with 35 m². Global, reflected and net radiation, soil heat flux and air temperature (dry and humid) were measured during the crop cycle. A Datalogger, which operated at 1 Hz frequency, storing 5 minutes averages was utilized. The global (K¯) and reflected (K­) radiations showed that the average transmission of global radiation (K¯in / K¯ex) was almost constant, near to 79.59%, while the average ratio of reflected radiation (K­in / K­ex) was 69.21% with 8.47% standard-deviation. The normalized curves of short-wave net radiation, in relation to the global radiation (K*/ K¯), found for both environments, were almost constant at the beginning of cycle; this relation decreased in the final stage of culture. The normalized relation (Rn/ K¯) was bigger in the external area, about 12%, when the green culture covered the soil surface. The long-wave radiation balance average (L*) was bigger outside, about 50%. The energy balance, estimated in terms of vertical fluxes, showed that, for the external area, in average, 83.07% of total net radiation was converted in latent heat evaporation (LE), and 18% in soil heat flux (G), and 9.96% in sensible heat (H), while inside of the greenhouse, 58.71% of total net radiation was converted in LE, 42.68% in H, and 28.79% in G.

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.

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.

Kerrostumien muodostumisen määritys sondin lämmönsiirron perusteella

Diplomityön tavoitteena oli tutkia kerrostumien muodostumista voimalaitoskattiloiden lämpöpinnoille ja niiden vaikutusta lämmönsiirtoon. Kerrostumien vaikutusta lämpöpintojen lämpötiloihin tutkittiin kerrostumasondin avulla kierto- ja kerrosleijukattiloissa. Saadun mittausaineiston muodostettiin mittaustilannetta vastaava laskentamalli, jonka avulla pystytään selvittämään voimalaitoksen lämpöpintojen likaantumista käyttöperiodin aikana. Työn alussa tarkasteltiin kerrostumien muodostumismekanismeja lämmönsiirtopinnoille, kerrostumien ominaisuuksia ja virtausympäristön vaikutusta lämmönsiirtoon. Poltossa syntyvä tuhka kiinnittyy eri mekanismeilla kattilan seinille ja lämpöpinnoille riippuen polttoaineesta ja ympäröivistä olosuhteista. Likakerrosten muodostuminen lämmönsiirtimen pinnalle pienentää putkeen siirtyvää lämpövirtaa ja alentaa kattilan hyötysuhdetta. Kerrostumien lämmönjohtumisesta ja mikrorakenteesta on hyvin vähän kokeellisen tietoa, mikä vaikeuttaa reaalisten mallien muodostamista kattilan käyttäytymisestä sen likaantuessa.Työn kokeellisessa osassa tarkasteltiin kahdella eri kattilalla tehtyjä mittauksia. Mittaukset tehtiin kerrostumasondilla, jonka lämmönsiirrosta luotiin laskentamalli SIMULINK-simulointiohjelmalla. Mittaustuloksina saatiin kolmen eri pinnan lämpötilat, jotka muuttuivat kerrostuman ja jäähdytyksen vaikutuksesta. Laskentamallista muokattiin mittaustilanteita vastaava, jolloin lämpötilamuutoksista nähdään likakerroksen ominaisuuksien vaikutus lämmönsiirtoon. Sondista muodostettiin myös FLUENT-malli, jolla tarkasteltiin yksittäisen putken virtauskenttää sekä kahden lähekkäin olevan putken vaikutusta virtaukseen.

Muutosilmiöt soodakattilakeossa

Soodakattilan sulakeon epästationaarinen käyttäytyminen sekä keon pitkä jäähtymisaika alasajon jälkeen ovat aiheuttaneet ongelmia kattilan taloudellisessa käytettävyydessä. Keon käyttäytymisestä on luotu CFD-malleja, joiden tavoitteena on havainnollistaa keon lämpötilajakaumaa ja rakennetta. Mallien ongelmana on se, että niissä huomioidaan vain keon aktiivinen pintakerros. Keon sisäosan rakennetta ja siinä tapahtuvia prosesseja ei toistaiseksi tunneta kunnolla luotettavan, koko keon kattavan mallin luomiseen. Tässä työssä tutkittiin sulakeon käytön aikana havaittuja muutosilmiöitä, jotka vaikuttavat keon rakenteeseen ja ominaisuuksiin sekä tutkittiin ilmiöiden taustalla olevia tekijöitä. Näitä tekijöitä ovat keon sisässä tapahtuvat kemialliset ja fyysiset prosessit, jotka aiheuttavat muutoksia niin lämpöteknisesti kuin fyysisesti sekä ulkoapäin tulevat tekijät, jotka aiheutuvat ajotilanteiden seurauksena tapahtuvista muutoksista. Työn kokeellisena osana luotiin sulakeon jäähtymismalli käyttäen 1-dimensionaalista ADL-mallia. Mallin pohjana käytettiin StoraEnso Oy:ltä Oulun soodakattilan sulakeosta saatua mittausraporttia. ADL-mallin avulla luotiin keon jäähtymiskäyrät lämpötilan ja syvyyden funktiona. Saadut käyrät täsmäsivät hyvin mittausraportin tuloksiin. Mallin avulla keolle saatiin muodostettua energiatase, jonka tuloksena keosta 12 tunnin aikana poistuva lämpövirta pinnalla oli noin 9.8kW/m2 ja pinnan lämmönsiirtokerroin 58.3W/m²°C. Pohjan poistuvaksi lämpövirraksi saatiin 14.1kW/m2 ja lämmönsiirtokertoimeksi 75.4W/m²°C. Termiseksi diffuusiokertoimeksi saatiin 3.9•10-7m²/s.

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

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.