Convection of water vapour in snowpacks

This paper studies numerically the convection of water vapour in snowpacks using an Eulerian-Eulerian two-phase approach. The convective water vapour transport in snow and its effects on snow density are often invoked to explain observed density profiles, e.g. of thin Arctic snow covers, but this process has never been numerically simulated and analysed in a systematic manner. Here, the impact of convection on the thermal and phase change regimes as a function of different snowpack depths, thermal boundary conditions and Rayleigh numbers is analysed. We find considerable impact of natural convection on the snow density distribution with a layer of significantly lower density at the bottom of the snowpack and a layer of higher density located higher in the snowpack or at the surface. We find that emergent heterogeneity in the snow porosity results in a feedback effect on the spatial organization of convection cells causing their horizontal displacement. Regions where the snowpack is most impacted by phase changes are found to be horizontally extended and vertically thin, `pancake'-like layers at the top and bottom of the snowpack. We further demonstrate that among the parameters important for natural convection, the snowpack depth has the strongest influence on the heat and mass transfer. Despite our simplifying assumptions, our study represents a significant improvement over the state of the art and a first step to accurately simulate snowpack dynamics in diverse regions of the cryosphere.

The Impact of Diffusive Water Vapor Transport on Snow Profiles in Deep and Shallow Snow Covers and on Sea Ice

Margaux Caroline Patricia Couttet, Isabelle Amandine Gouttevin, Mahdi Jafari, Michael Lehning, Adrien Michel, Varun Sharma, Nander Wever

Water vapor transport has been highlighted as a critical process in Arctic snowpacks, shaping the snow cover structure in terms of density, thermal conductivity, and temperature profile among others. Here, we present an attempt to describe the thermally-induced vertical diffusion of water vapor in the snow cover and its effects of the snowpack structure using the SNOWPACK model. Convection, that may also constitute a significant part of vapor transport, is not addressed. Assuming saturated conditions at the upper boundary of the snowpack and as initial condition, the vapor flux between snow layers is expressed by a 1-dimensional transient diffusion equation, which is solved with a finite difference routine. The implications on the snowpack of this vertical diffusive flux, are analyzed using metrics such as the cumulative density change due to diffusive vapor transport, the degree of over- or undersaturation, the instantaneous snow density change rate, and the percentage of snow density change. We present results for four different regions sampling the space of natural snow cover variability: Alpine, Subarctic, Arctic, and Antarctic sea ice. The largest impact of diffusive water vapor transport is observed in snow on sea ice in the Weddell Sea and the shallow Arctic snowpack. The simulations show significant density reductions upon inclusion of diffusive water vapor transport: cumulative density changes from diffusive vapor transport can reach −62 and −66 kg m−3 for the bottom layer in the shallow Arctic snowpack and snow on sea ice, respectively. For comparison, in deeper snow covers, they rarely exceed −40 kg m−3. This leads to changes in density for shallow snowpacks at the soil-snow interface in the range of −5 to −21%. Mirroring the density decease at depth is a thicker deposition layer above it with increase in density around 7.5%. Similarly, for the sea ice, the density decreased at the sea ice-snow interface by −20%. We acknowledge that vapor transport by diffusion may in some snow covers—such as in thin tundra snow—be small compared to convective transport, which will have to be addressed in future work.

2020

Thermo-Hydrodynamics in a Closed Loop Pulsating Heat Pipe

Giulia Spinato

Pulsating heat pipes (PHPs) represent a promising solution for passive on-chip, two-phase cooling of micro-electronics, providing advantages such as a simple construction and operation in any gravitational orientation. Unfortunately, the unique coupling of thermodynamics, hydrodynamics and heat transfer responsible for their operation has so far eluded comprehensive description or accurate prediction. The complexity of the self-sustained two-phase flow in PHPs presents many challenges to the understanding on the physical phenomena taking place. It is important to evaluate the heat and mass transfer mechanisms occurring during their operation in order to better describe their performance as a function of the operating conditions. In the present study, a new facility at the Laboratory of Heat and Mass Transfer (LTCM) was built to allow the synchronized thermal and visual investigation of a Closed Loop Pulsating Heat Pipe (CLPHP). A single-turn channel CLPHP was investigated using R245fa as the working fluid. The tests were carried out at filling ratios from 10 to 90 % and heat inputs from 2 to 60 W, for vertical and inclined orientation. Flow visualization was attained via the transparent front side of the test section which provided full optical access to the flow inside of the CLPHP channels. A novel time-strip image processing technique was applied to the high speed videos to extract qualitative details of the flow regimes and quantitative flow data concerning the liquid/vapor interface dynamics. Local temperature oscillations were also measured and their frequency spectra further helped in characterizing the self-sustained two-phase flow. Thermal resistance measurements were used to qualitatively and quantitatively assess the effect of the flow dynamics on the system thermal performance, which are presented as operational maps. The dynamics governing the PHP operation during oscillating and circulating flows and their effects on the system thermal performance were investigated and insight gained. Four distinct flow regimes and their thermal and flow-dynamics characteristics were identified, suggesting the strong coupling between the two-phase flow pattern and the system thermal behavior. Thin film evaporation was observed to be the most dominant thermal mechanism while heat transfer into the oscillating liquid slug was of second importance, together with localized nucleate boiling. Moreover, the net forces acting on the system could be identified through the novel time-strip technique, revealing new details on the mechanisms producing self-sustained two-phase flow oscillation and circulation, and the two-phase flow pattern transition. The role of gravity for the operation of this single-turn CLPHP was also assessed.

EPFL2014

Replicated microcellular aluminium structures for thermal management

Alexandros Athanasiou-Ioannou

Open-pore aluminium foams are of interest in thermal management of power electronics due to their large specific surface area, low weight and relatively high thermal conductivity. The replication process is a low-cost and simple foam manufacturing method that offers many opportunities for the investigation of the influence exerted by the microstructure (porosity, pore size) and geometry or their combination in the performance of metal foams in heat-exchange applications. In this study, we explored replicated aluminium foams as a part of an integrated heat sink manufactured in one step that seamlessly combines metal and ceramic phases. Replicated aluminium foams were manufactured and tested first under natural and then under forced convection. A dedicated test apparatus was built to test these samples under both natural and forced convection, using the guarding ring technique. In the case of natural convection, the aluminium foams had different pore sizes (125-180 µm, 400-450 µm and 5 mm), porosities (78 to 86%) and geometry (disk, cylinder and finned structures). In all cases, replicated aluminium foams were found to dissipate the same amount of heat as their bulk aluminium (alloy) equivalents. The thermal heat exchange performance of replicated microcellular aluminium was also measured under forced convection. The test rig was modified to measure the pressure drop and volumetric air flow in a cylindrically symmetric jet impingement configuration. Foams of different pore size (125-180, 400-450 and 900-1300 µm), porosity (between 75 and 85 %), height (1 and 2.8 cm) and infiltrated at various pressures (between 2.7 and 70 bars) were manufactured for thermal testing, all having a hollow central channel, through which air is injected at a range of flow rates. The results show that within the covered parameter range foams of higher porosities, larger pore sizes and lower infiltration pressures dissipate more heat at lower cost in terms of pumping power. Shorter heights also dissipate more heat for a given volumetric flow but do so at the cost of an increased pressure drop and hence greater pumping power expenditure. A numerical model, developed and written by Dr. D. Ingram, using the finite-volumes technique, was employed to elucidate how the various parameters influence the heat transfer. The model solves the Darcy-Forchheimer formulation of fluid flow in porous media, which is then coupled to a convective/conductive heat transfer model of the same structure. These simulations were used to perform parametric studies on the influence of height and pore size with regard to the thermal behaviour of replicated foams. Furthermore, the model was benchmarked against the experimental results; overall, despite some discrepancies with measured temperature profiles, the model captures data trends well, and provides insight on the physical phenomena that underlie the experimental data and observations. In particular, the presence and influence of recirculating air flow along the inner portion of the heat exchanger structures was identified. Finally, we successfully fabricated the proposed integrated multimaterial structure and tested it under forced convection, using the same conditions as the replicated aluminium foams. Results demonstrate that this type of structures can be effectively used to cool power electronics.

EPFL2017