Developing discontinuous Galerkin methods for the resolution of the incompressible Navier-Stokes equations

Informe de investigación elaborado a partir de una estancia en el Laboratorio de Diseño Computacional en Aeroespacial en el Massachusetts Institute of Technology (MIT), Estados Unidos, entre noviembre de 2006 y agosto de 2007. La aerodinámica es una rama de la dinámica de fluidos referida al estudio de los movimientos de los líquidos o gases, cuya meta principal es predecir las fuerzas aerodinámicas en un avión o cualquier tipo de vehículo, incluyendo los automóviles. Las ecuaciones de Navier-Stokes representan un estado dinámico del equilibrio de las fuerzas que actúan en cualquier región dada del fluido. Son uno de los sistemas de ecuaciones más útiles porque describen la física de una gran cantidad de fenómenos como corrientes del océano, flujos alrededor de una superficie de sustentación, etc. En el contexto de una tesis doctoral, se está estudiando un flujo viscoso e incompresible, solucionando las ecuaciones de Navier- Stokes incompresibles de una manera eficiente. Durante la estancia en el MIT, se ha utilizado un método de Galerkin discontinuo para solucionar las ecuaciones de Navier-Stokes incompresibles usando, o bien un parámetro de penalti para asegurar la continuidad de los flujos entre elementos, o bien un método de Galerkin discontinuo compacto. Ambos métodos han dado buenos resultados y varios ejemplos numéricos se han simulado para validar el buen comportamiento de los métodos desarrollados. También se han estudiado elementos particulares, los elementos de Raviart y Thomas, que se podrían utilizar en una formulación mixta para obtener un algoritmo eficiente para solucionar problemas numéricos complejos.

Influència del context geològic i hidrogeològic en l'eficiència energètica dels sistemes geotèrmics de molt baixa entalpia

L’elevat consum energètic de les societats actuals, així com la impossibilitat de sostenir-lo a llarg termini implica la cerca de noves fonts d’energia. D’aquesta manera, en el camp de la climatització residencial, l’energia geotèrmica de molt baixa entalpia es posiciona com una alternativa als recursos energètics actuals. Així, els primers metres de subsòl presenten una temperatura adequada per al seu aprofitament calorífic, mitjançant els sistemes geotèrmics de bomba de calor. Si bé a nivell energètic, són sistemes, intrínsecament, molt eficients, el seu rendiment pot patir importants variacions davant dels canvis en les condicions del medi geològic i hidrogeològic. Especialment, els col·lectors de calor verticals, treballen, amb freqüència, en el si de les formacions hidrogeològiques. En aquest sentit, els canvis del nivell hidràulic i de la temperatura de l’aigua de l’aqüífer es manifesten amb variacions de la conductivitat tèrmica equivalent i del flux subterrani d’aigua, que alhora, aquestes, es tradueixen en alteracions del flux subterrani de calor. Davant d’aquest fet, l’avaluació quantitativa de l’efecte d’aquestes fluctuacions que es presenta en aquest treball mostra petites variacions del nivell hidràulic i de la temperatura de l’aigua comporten canvis molt notables en l’eficiència dels sistemes verticals de bomba de calor geotèrmica.

Simulium (Chirostilbia) bifenestratum (Diptera, Simuliidae), a new black-fly species from the Atlantic forest, State of São Paulo, Brazil

The larva, pupa, male, and female of Simulium bifenestratum n. sp. are described and illustrated. The pupae of the new species have 10 gill filaments, thick at their base and arranged in a three-dimensional way, surrounding the head and thorax. Its pupal cocoon is peculiar, not found in any of the known Brazilian black-fly species; it is very thick and hard with two openings in the anterior region. S. bifenestratum n. sp. was collected in one stream in the Bocaina mountain chain, Atlantic forest, in São José do Barreiro county, state of São Paulo, in a high (1500 m) natural grassland. Larvae and pupae were collected on the edges of small waterfalls and in places with-high speed laminar water flow, attached to the bedrock.

Numerical modeling of the effect of variation of boundary conditions on vadose zone hydraulic properties

An accurate estimation of hydraulic fluxes in the vadose zone is essential for the prediction of water, nutrient and contaminant transport in natural systems. The objective of this study was to simulate the effect of variation of boundary conditions on the estimation of hydraulic properties (i.e. water content, effective unsaturated hydraulic conductivity and hydraulic flux) in a one-dimensional unsaturated flow model domain. Unsaturated one-dimensional vertical water flow was simulated in a pure phase clay loam profile and in clay loam interlayered with silt loam distributed according to the third iteration of the Cantor Bar fractal object Simulations were performed using the numerical model Hydrus 1D. The upper and lower pressure heads were varied around average values of -55 cm for the near-saturation range. This resulted in combinations for the upper and lower constant head boundary conditions, respectively, of -50 and -60 cm, -40 and -70 cm, -30 and -80 cm, -20 and -90 cm, and -10 and -100 cm. For the drier range the average head between the upper and lower boundary conditions was set to -550 cm, resulting in the combinations -500 and -600 cm, -400 and -700 cm, -300 and -800 cm, -200 and -900 cm, and -100 and -1,000 cm, for upper and lower boundary conditions, respectively. There was an increase in water contents, fluxes and hydraulic conductivities with the increase in head difference between boundary conditions. Variation in boundary conditions in the pure phase and interlayered one-dimensional profiles caused significant deviations in fluxes, water contents and hydraulic conductivities compared to the simplest case (a head difference between the upper and lower constant head boundaries of 10 cm in the wetter range and 100 cm in the drier range).

Determination of hydraulic properties of a tropical soil of Hawaii using column experiments and inverse modeling

A method for determining soil hydraulic properties of a weathered tropical soil (Oxisol) using a medium-sized column with undisturbed soil is presented. The method was used to determine fitting parameters of the water retention curve and hydraulic conductivity functions of a soil column in support of a pesticide leaching study. The soil column was extracted from a continuously-used research plot in Central Oahu (Hawaii, USA) and its internal structure was examined by computed tomography. The experiment was based on tension infiltration into the soil column with free outflow at the lower end. Water flow through the soil core was mathematically modeled using a computer code that numerically solves the one-dimensional Richards equation. Measured soil hydraulic parameters were used for direct simulation, and the retention and soil hydraulic parameters were estimated by inverse modeling. The inverse modeling produced very good agreement between model outputs and measured flux and pressure head data for the relatively homogeneous column. The moisture content at a given pressure from the retention curve measured directly in small soil samples was lower than that obtained through parameter optimization based on experiments using a medium-sized undisturbed soil column.

Rill erosion on an oxisol influenced by a thin compacted layer

The presence of compacted layers in soils can induce subprocesses (e.g., discontinuity of water flow) and induces soil erosion and rill development. This study assesses how rill erosion in Oxisols is affected by a plow pan. The study shows that changes in hydraulic properties occur when the topsoil is eroded because the compacted layer lies close below the surface. The hydraulic properties that induce sediment transport and rill formation (i.e., hydraulic thresholds at which these processes occur) are not the same. Because of the resistance of the compacted layer, the hydraulic conditions leading to rill incision on the soil surface differed from the conditions inducing rill deepening. The Reynolds number was the best hydraulic predictor for both processes. The formed rills were shallow and could easily be removed by tillage between crops. However, during rill development, large amounts of soil and contaminants could also be transferred.

Contribution of microbial and invertebrate communities to leaf litter colonization in a Mediterranean stream

Leaf litter inputs and retention play an important role in ecosystem functioning in forested streams. We examined colonization of leaves by microbes (bacteria, fungi, and protozoa) and fauna in Fuirosos, an intermittent forested Mediterranean stream. Black poplar (Populus nigra) and plane (Platanus acerifolia) leaf packs were placed in the stream for 4 mo. We measured the biomasses and calculated the densities of bacteria, fungi, protozoa, meiofauna, and macroinvertebrates to determine their dynamics and potential interactions throughout the colonization process. Colonization was strongly correlated with hydrological variability (defined mainly by water temperature and discharge). The 1st week of colonization was characterized by hydrological stability and warm water temperatures, and allocation of C from microbial to invertebrate compartments on the leaf packs was rapid. Clumps of fine particulate organic matter (FPOM) were retained by the leaf packs, and enhanced rapid colonization by microfauna and meiofaunal collector-gatherers (ostracods and copepods). After 2 wk, an autumnal flood caused a 20-fold increase in water flow. Higher discharge and lower water temperature caused FPOM-related fauna to drift away from the packs and modified the subsequent colonization sequence. Fungi showed the highest biomass, with similar values to those recorded at the beginning of the experiment. After 70 d of postflood colonization, fungi decreased to nearly 40% of the total C in the leaf packs, whereas invertebrates became more abundant and accounted for 60% of the C. Natural flood occurrence in Mediterranean streams could be a key factor in the colonization and processing of organic matter.

Stream Stabilization in Western Iowa: Structure Evaluation and Design Manual, HR-385, 1998

Stream degradation is the action of deepening the stream bed and widening the banks due to the increasing velocity of water flow. Degradation is pervasive in channeled streams found within the deep to moderately deep loess regions of the central United States. Of all the streams, however, the most severe and widespread entrenchment occurs in western Iowa streams that are tributaries to the Missouri River. In September 1995 the Iowa Department of Transportation awarded a grant to Golden Hills Resource Conservation and Development, Inc. The purpose of the grant, HR-385 "Stream Stabilization in Western Iowa: Structure Evaluation and Design Manual", was to provide an assessment of the effectiveness and costs of various stabilization structures in controlling erosion on channeled streams. A review of literature, a survey of professionals, field observations and an analysis of the data recorded on fifty-two selected structures led to the conclusions presented in the project's publication, Design Manual, Streambed Degradation and Streambank Widening in Western Iowa. Technical standards and specifications for the design and construction of stream channel stabilization structures are included in the manual. Additional information on non-structural measures, monitoring and evaluation of structures, various permit requirements and further resources are also included. Findings of the research project and use and applications of the Design Manual were presented at two workshops in the Loess Hills region. Participants in these workshops included county engineers, private contractors, state and federal agency personnel, elected officials and others. The Design Manual continues to be available through Golden Hills Resource Conservation and Development.

Use of Longitudinal Subdrains in the 3R Program, HR-509, Progress Report #2, 1981

The main consideration for base construction under the pavement, in the design of Iowa's interstate, was structural capacity. The material was dense graded with the aim of supporting the pavement and distributing the load as it is transferred to the underlying grade. The drainage characteristics of the base was apparently not given adequate consideration. On jointed portland cement concrete pavement, the water that is trapped immediately beneath the pavement causes severe problems. The traffic causes rapid movement of the water resulting in the hydraulic pressures or "pumping" (movement and redeposit of base fine material), further resulting in faulting between individual slabs. The objective of this evaluation is to determine if longitudinal subdrains are effective in preventing or reducing pumping, faulting and related deterioration. Results suggest that, based upon the flow from the outlets observed during periodic checks and evidence of water flow at the outlets, it appears that to date the subdrains are effective in draining the subbase and subgrade. Because of the limited data available at this time, however, the pavement condition and faulting results are inconclusive.

Use of Longitudinal Subdrains in the 3R Program, Progress Report 1, HR-509, 1980

The main consideration for base construction under the pavement, in the design of Iowa's interstate, was structural capacity. The material was dense graded with the aim of supporting the pavement and distributing the load as it is transferred to the underlying grade. The drainage characteristics of the base was apparently not given adequate consideration. On jointed portland cement concrete pavement, the water that is trapped immediately beneath the pavement causes severe problems. The traffic causes rapid movement of the water resulting in the hydraulic pressures or "pumping" (movement and redeposit of base fine material), further resulting in faulting between individual slabs. The objective of this evaluation is to determine if longitudinal subdrains are effective in preventing or reducing pumping, faulting and related deterioration. Results suggest that, based upon the flow from the outlets observed during periodic checks and evidence of water flow at the outlets, it appears that to date the subdrains are effective in draining the subbase and subgrade. Because of the limited data available at this time, however, the pavement condition and faulting results are inconclusive.

Self-potentials in partially saturated media: the importance of explicit modeling of electrode effects

Self-potential (SP) data are of interest to vadose zone hydrology because of their direct sensitivity to water flow and ionic transport. There is unfortunately little consensus in the literature about how to best model SP data under partially saturated conditions, and different approaches (often supported by one laboratory data set alone) have been proposed. We argue that this lack of agreement can largely be traced to electrode effects that have not been properly taken into account. A series of drainage and imbibition experiments were considered in which we found that previously proposed approaches to remove electrode effects were unlikely to provide adequate corrections. Instead, we explicitly modeled the electrode effects together with classical SP contributions using a flow and transport model. The simulated data agreed overall with the observed SP signals and allowed decomposing the different signal contributions to analyze them separately. After reviewing other published experimental data, we suggest that most of them include electrode effects that have not been properly taken into account. Our results suggest that previously presented SP theory works well when considering the modeling uncertainties presently associated with electrode effects. Additional work is warranted to not only develop suitable electrodes for laboratory experiments but also to assure that associated electrode effects that appear inevitable in longer term experiments are predictable, so that they can be incorporated into the modeling framework.

Étude des sédiments quaternaires, de la molasse et sa tectonique, dans le Grand Lac (Léman) à partir de données sismiques 2D et 3D

RESUME L'Institut de Géophysique de l'Université de Lausanne a développé au cours de ces dernières années un système d'acquisition de sismique réflexion multitrace à haute résolution 2D et 3D. L'objectif de cette thèse était de poursuivre ce développement tout améliorant les connaissances de la géologie sous le lac Léman, en étudiant en particulier la configuration des grands accidents sous-lacustres dans la Molasse (Tertiaire) qui forme l'essentiel du substratum des formations quaternaires. En configuration 2D, notre système permet d'acquérir des profils sismiques avec une distance inter-CDP de 1,25 m. La couverture varie entre 6 et 18 selon le nombre de traces et la distance inter-tir. Le canon à air (15/15 eu. in.), offre une résolution verticale de 1,25 ni et une pénétration maximale de 300 m sous le fond de l'eau. Nous avons acquis au total plus de 400 km de sections 2D dans le Grand Lac et le Haut Lac entre octobre 2000 et juillet 2004. Une campagne de sismique 3D a fourni des données au large d'Evian sur une surface de 442,5 m sur 1450 m, soit 0,64 km2. La navigation ainsi que le positionnement des hydrophones et de la source ont été réalisés avec des GPS différentiels. Nous avons utilisé un traitement sismique conventionnel, sans appliquer d'AGC et en utilisant une migration post-stack. L'interprétation du substratum antéquaternaire est basée sur l'identification des sismofaciès, sur leurs relations avec les unités géologiques adjacentes au lac, ainsi que sur quelques données de forages. Nous obtenons ainsi une carte des unités géologiques dans le Grand Lac. Nous précisons la position du chevauchement subalpin entre la ville de Lausanne, sur la rive nord, et le bassin de Sciez, sur la rive sud. Dans la Molasse de Plateau, nous avons identifié les décrochements de Pontarlier et de St. Cergue ainsi que plusieurs failles non reconnues jusqu'ici. Nous avons cartographié les accidents qui affectent la Molasse subalpine ainsi que le plan de chevauchement du flysch sur la Molasse près de la rive sud du lac. Une nouvelle carte tectonique de la région lémanique a ainsi pu être dressée. L'analyse du substratum ne montre pas de failles suggérant une origine tectonique de la cuvette lémanique. Par contre, nous suggérons que la forme du creusement glaciaire, donc de la forme du lac Léman, a été influencée par la présence de failles dans le substratum antéquaternaire. L'analyse des sédiments quaternaires nous a permis de tracer des cartes des différentes interfaces ou unités qui les composent. La carte du toit du substratum antéquaternaire montre la présence de chenaux d'origine glaciaire dont la profondeur maximale atteint la cote -200 ni. Leur pente est dirigée vers le nord-est, à l'inverse du sens d'écoulement actuel des eaux. Nous expliquons cette observation par l'existence de circulations sous-glaciaires d'eau artésienne. Les sédiments glaciaires dont l'épaisseur maximale atteint 150 ni au centre du lac ont enregistré les différentes récurrences glaciaires. Dans la zone d'Evian, nous mettons en évidence la présence de lentilles de sédiments glaciolacustres perchées sur le flanc de la cuvette lémanique. Nous avons corrélé ces unités avec des données de forage et concluons qu'il s'agit du complexe inférieur de la pile sédimentaire d'Evian. Celui-ci, âgé de plus de 30 000 ans, serait un dépôt de Kame associé à un lac périglaciaire. La sismique réflexion 3D permet de préciser l'orientation de l'alimentation en matériel détritique de l'unité. La finesse des images obtenues nous permet également d'établir quels types d'érosion ont affecté certaines unités. Les sédiments lacustres, dont l'épaisseur maximale imagée atteint plus de 225 m et sans doute 400 ni sous le delta du Rhône, indiquent plusieurs mécanismes de dépôts. A la base, une mégaturbidite, épaisse d'une trentaine de mètres en moyenne, s'étend entre l'embouchure de la Dranse et le delta du Rhône. Au-dessus, la décantation des particules en suspension d'origine biologique et détritique fournit l'essentiel des sédiments. Dans la partie orientale du lac, les apports détritiques du Rhône forment un delta qui prograde vers l'ouest en s'imbriquant avec les sédiments déposés par décantation. La structure superficielle du delta a brutalement évolué, probablement à la suite de l'évènement catastrophique du Tauredunum (563 A.D.). Sa trace probable se marque par la présence d'une surface érosive que nous avons cartographiée. Le delta a ensuite changé de géométrie, avec notamment un déplacement des chenaux sous-lacustres. Sur l'ensemble de nos sections sismiques, nous n'observons aucune faille dans les sédiments quaternaires qui attesterait d'une tectonique postglaciaire du substratum. ABSTRACT During the last few years the institute of Geophysics of the University of Lausanne cleveloped a 2D and 3D high-resolution multichannel seismic reflection acquisition system. The objective of the present work was to carry on this development white improving our knowledge of the geology under Lake Geneva, in particular by studying the configuration of the large accidents affecting the Tertiary Molasse that makes up the basement of most Quaternary deposits. In its 2D configuration, our system makes it possible to acquire seismic profiles with a CDP interval of 1.25 m. The fold varies from 6 to 18 depending on the number of traces and the shooting interval. Our air gun (15/15 cu. in.) provides a vertical resolution of 1.25 m and a maximum penetration depth of approximately 300 m under water bottom. We acquired more than 400 km of 2D sections in the Grand Lac and the Haut Lac between October 2000 and July 2004. A 3D seismic survey off the city of Evian provided data on a surface of 442.5 m x 1450 m (0.64 km2). Ship's navigation as well as hydrophone- and source positioning were carried out with differential GPS. The seismic data were processed following a conventional sequence without .applying AGC and using post-stack migration. The interpretation of the pre-Quaternary substratum is based on sismofacies, on their relationships with terrestrial geological units and on some borehole data. We thus obtained a map of the geological units in the Grand Lac. We defined the location of the subalpine thrust from Lausanne, on the north shore, to the Sciez Basin, on the south shore. Within the Molasse de Plateau, we identified the already know Pontarlier and St Cergue transforms Fault as well as faults. We mapped faults that affect subalpine Molasse as well as the thrust fault plane between alpine flysch and Molasse near the lake's south shore. A new tectonic map of the Lake Geneva region could thus be drawn up. The substratum does not show faults indicating a tectonic origin for the Lake Geneva Basin. However, we suggest that the orientation of glacial erosion, and thus the shape of Lake Geneva, vas influenced by the presence of faults in the pre-Quaternary basement. The analysis of Quaternary sediments enabled us to draw up maps of various discontinuities or internal units. The top pre-Quaternary basement map shows channels of glacial origin, the deepest of them reaching an altitude of 200 m a.s.l. The channel's slopes are directed to the North-East, in opposite direction of the present water flow. We explain this observation by the presence of artesian subglacial water circulation. Glacial sediments, the maximum thickness of which reaches 150 m in the central part of the lake, record several glacial recurrences. In the Evian area, we found lenses of glacio-lacustrine sediments set high up on the flank of the Lake Geneva Bassin. We correlated these units with on-land borehole data and concluded that they represent the lower complex of the Evian sedimentary pile. The lower complex is aider than 30 000 years, and it could be a Kame deposit associated with a periglacial lake. Our 3D seismic reflexion survey enables us to specify the supply direction of detrital material in this unit. With detailed seismic images we established how some units were affected by different erosion types. The lacustrine sediments we imaged in Lake Geneva are thicker than 225 m and 400 m or more Linder the Rhone Delta. They indicate several depositional mechanisms. Their base is a major turbidite, thirty meters thick on average, that spreads between the Dranse mouth and the Rhone delta. Above this unit, settling of suspended biological and detrital particles provides most of the sediments. In the eastern part of the lake, detrital contribution from the Rhone builds a delta that progrades to the west and imbricates with the settling sediments. The shallow structure of the Rhone delta abruptly evolved, probably after the catastrophic Tauredunum event (563 A.D.). It probably coincides with an erosive surface that we mapped. As a result, the delta geometry changed, in particular associated with a displacement of water bottom channels. In all our seismic sections, we do not observe fault in the Quaternary sediments that would attest postglacial tectonic activity in the basement.

Automation of DEM Cutting for Hydrologic/Hydraulic Modeling, TR-631, 2015

Hydrologic analysis is a critical part of transportation design because it helps ensure that hydraulic structures are able to accommodate the flow regimes they are likely to see. This analysis is currently conducted using computer simulations of water flow patterns, and continuing developments in elevation survey techniques result in higher and higher resolution surveys. Current survey techniques now resolve many natural and anthropogenic features that were not practical to map and, thus, require new methods for dealing with depressions and flow discontinuities. A method for depressional analysis is proposed that uses the fact that most anthropogenically constructed embankments are roughly more symmetrical with greater slopes than natural depressions. An enforcement method for draining depressions is then analyzed on those depressions that should be drained. This procedure has been evaluated on a small watershed in central Iowa, Walnut Creek of the South Skunk River, HUC12 # 070801050901, and was found to accurately identify 88 of 92 drained depressions and place enforcements within two pixels, although the method often tries to drain prairie pothole depressions that are bisected by anthropogenic features.

Heat Transfer Analysis of the EPR Core Catcher Test Facility Volley

Uusi EPR-reaktorikonsepti on suunniteltu selviytymään tapauksista, joissa reaktorinsydän sulaa ja sula puhkaisee paineastian. Suojarakennuksen sisälle on suunniteltu alue, jolle sula passiivisesti kerätään, pidätetään ja jäähdytetään. Alueelle laaditaan valurautaelementeistä ns.sydänsieppari, joka tulvitetaan vedellä. Sydänsulan tuottama jälkilämpö siirtyyveteen, mistä se poistetaan suojarakennuksen jälkilämmönpoistojärjestelmän kautta. Suuri osa lämmöstä poistuu sydänsulasta sen yläpuolella olevaan veteen, mutta lämmönsiirron tehostamiseksi myös sydänsiepparin alapuolelle on sijoitettu vedellä täytettävät jäähdytyskanavat. Jotta sydänsiepparin toiminta voitaisiin todentaa, on Lappeenrannan Teknillisellä Yliopistolla rakennettu Volley-koelaitteisto tätä tarkoitusta varten. Koelaitteisto koostuu kahdesta täysimittaisesta valuraudasta tehdystä jäähdytyskanavasta. Sydänsulan tuottamaa jälkilämpöä simuloidaan koelaitteistossa sähkövastuksilla. Tässä työssä kuvataan simulaatioiden suorittaminen ja vertaillaan saatuja arvoja mittaustuloksiin. Työ keskittyy sydänsiepparista jäähdytyskanaviin tapahtuvan lämmönsiirron teoriaan jamekanismeihin. Työssä esitetään kolme erilaista korrelaatiota lämmönsiirtokertoimille allaskiehumisen tapauksessa. Nämä korrelaatiot soveltuvat erityisesti tapauksiin, joissa vain muutamia mittausparametreja on tiedossa. Työn toinen osa onVolley 04 -kokeiden simulointi. Ensin käytettyä simulointitapaa on kelpoistettuvertaamalla tuloksia Volley 04 ja 05 -kokeisiin, joissa koetta voitiin jatkaa tasapainotilaan ja joissa jäähdytteen käyttäytyminen jäähdytyskanavassa on tallennettu myös videokameralla. Näiden simulaatioiden tulokset ovat hyvin samanlaisiakuin mittaustulokset. Korkeammilla lämmitystehoilla kokeissa esiintyi vesi-iskuja, jotka rikkoivat videoinnin mahdollistavia ikkunoita. Tämän johdosta osassa Volley 04 -kokeita ikkunat peitettiin metallilevyillä. Joitakin kokeita jouduttiin keskeyttämään laitteiston suurten lämpöjännitysten johdosta. Tällaisten testien simulaatiot eivät ole yksinkertaisia suorittaa. Veden pinnan korkeudesta ei ole visuaalista havaintoa. Myöskään jäähdytteen tasapainotilanlämpötiloista ei ole tarkkaa tietoa, mutta joitakin oletuksia voidaan tehdä samoilla parametreilla tehtyjen Volley 05 -kokeiden perusteella. Mittaustulokset Volley 04 ja 05 -kokeista, jotka on videoitu ja voitu ajaa tasapainotilaan saakka, antoivat simulaatioiden kanssa hyvin samankaltaisia lämpötilojen arvoja. Keskeytettyjen kokeiden ekstrapolointi tasapainotilaan ei onnistunut kovin hyvin. Kokeet jouduttiin keskeyttämään niin paljon ennen termohydraulista tasapainoa, ettei tasapainotilan reunaehtoja voitu ennustaa. Videonauhoituksen puuttuessa ei veden pinnan korkeudesta saatu lisätietoa. Tuloksista voidaan lähinnä esittää arvioita siitä, mitä suuruusluokkaa mittapisteiden lämpötilat tulevat olemaan. Nämä lämpötilat ovat kuitenkin selvästi alle sydänsiepparissa käytettävän valuraudan sulamislämpötilan. Joten simulaatioiden perusteella voidaan sanoa, etteivät jäähdytyskanavien rakenteet sula, mikäli niissä on pienikin jäähdytevirtaus, eikä useampia kuin muutama vierekkäinen kanava ole täysin kuivana.

Morphological and functional variability in the root system of Quercus ilex L. subject to confinement: consequences for afforestation

We examined root morphological and functional differences caused by restrictions imposed to vertical growth in the root system of holm oak (Quercus ilex L.) seedlings to assess the consequences of using nursery containers in the development of a confined root system for this species. Thus, root morphological, topological and functional parameters, including hydraulic conductance per leaf unit surface area (K $_{\rm RL})$, were investigated in one-year seedlings cultivated in three PVC tubes differing in length (20, 60 and 100 cm). Longer tubes showed greater projected root area, root volume, total and fine root lengths, specific root length (SRL) and K$_{\rm RL}$ values than did shorter tubes. On the other hand, the length of coarse roots (diameter > 4.5 mm) and the average root diameter were greater in shorter tubes. The strong positive correlation found between K$_{\rm RL}$ and SRL (r=+0.69; P<0.001) indicated that root thickness was inversely related to water flow through the root system. We concluded that root systems developed in longer tubes are more efficient for plant water uptake and, therefore, changes in root pattern produced in standard forest containers (i.e. about 20 cm length) may in fact prevent a proper establishment of the holm oak in the field, particularly in xeric environments.