Moisture and energy dynamics in fields soils

Soil moisture is crucial to water-cycle as it affects infiltration, runoff and land-atmosphere interactions. In dry soils the role of water vapor fluxes becomes essential and causes a tight coupling between soil moisture dynamics and heat transfer, which makes measurements and modeling extremely difficult. Since advances in both theory and experimental techniques are required, we have focused on both aspects: we have performed theoretical-numerical investigations of soil moisture and energy dynamics at the land surface; tested novel measurements techniques, in laboratory and at field scale, to provide new insights into soil moisture dynamics and contribute to a comprehensive theory of coupled energy and moisture transport. As the dynamics of physical systems are dictated by conservation equations, it is common practice to estimate quantities that cannot be directly measured (e.g. vapor fluxes) from the residuals of balance equations. In laboratory as well as in field experiments, residuals are calculated from measurements that are sparse in space and discrete in time. We have theoretically and numerically shown how the use of discrete data may lead to large residuals even without measurements or model errors. We have demonstrated that residuals are very sensitive to the distance between the measurements and that they increase in presence of nonlinear processes, such as moisture transport. A careful error analysis is required to assess the reliability of residuals computed from discrete data and avoid misinterpretations of the underlying physical processes. Another crucial element dictating soil moisture dynamics and evaporation into the atmosphere is soil water-retention. Popular parametric models of the retention curve present a matric potential that becomes infinite at residual saturation, which is the water content below which liquid flow stops. Few extensions of the water-retention curves have been proposed, which remain rarely used in practice. We have studied the effects of different models on liquid and vapor transport. We have shown that parameterizations that allow vapor flux below residual saturation yield larger vapor fluxes that result in a drier soil surface, whereas deep soil remains relatively wet. When diurnal forcing is considered, potential evaporation is reduced in the energy-limited regime and enhanced in the moisture-limited regime. We have established that in this case extended water-retention models predict larger heat fluxes that might critically impact models at whole scales (possibly including global circulation models). Discrimination between different theoretical descriptions and parametric models requires experiments that monitor moisture and energy dynamics. This demands the use of accurate sensors and measurement techniques that are able to simultaneously monitor several quantities (e.g. soil moisture, temperature, soil thermal properties, solute concentration in soil water). As in general different probes are used, problems arise due to different footprints and sampled soil volumes. Multi functional heat pulse probes (MFHPPs) have been conceived to overcome these issues by allowing simultaneous measurements of all required variables. We have assessed MFHPPs ability to measure ground heat flux, which has to be accurately known to avoid significant imbalances in the surface energy budget. We have employed an experimental setup able to generate a non-uniform heat-flux field and we have demonstrated that MFHPPs can reliably capture the magnitudes and directions of the fluxes. This allows more complete information on ground heat flux than that obtained by the classic probes (e.g. ground heat flux plates). To monitor soil state in large-scale field applications we have investigated the possibility to infer distributed thermal conductivity and soil moisture profiles with an optical fiber. Information is obtained by monitoring soil response to an active heat pulse generated from the metal shield armoring the optic cable, and analyzing the temperature differences of the cooling phase, monitored each meter by a distributed temperature sensing system (DTS). We have shown that soil moisture measurements inferred from soil thermal response (recorded by the DTS) are in good agreement with independent measurements in wet and intermediately wet soils; however in dry soils moisture content is underestimated and further investigations are required to improve this technique for dry conditions. In this study we have significantly improved the knowledge of coupled energy and moisture transport on both theoretical-numerical and experimental sides. We have proven that reliable measurements can be obtained at different spatiotemporal scales, and provided both prognostic and diagnostic methods to optimize their use. Finally, we have suggested crucial improvements to numerical models which might have implications beyond field-soils dynamics.

Heterogeneous kinetics of some atmospherically relevant reactions

A low pressure reactor (Knudsen cell) coupled to molecular beam modulated mass spectrometry was used to study the heterogeneous kinetics of the following three systems: (1) NO2 on amorphous carbon, (2) H2O vapor on ice and (3) gaseous HCl on ice. The first system was investigated at ambient temperature while the two last ones were investigated at temperatures typical of the stratosphere, that is roughly between 160 and 220K. The solid phase samples correspond to model surfaces whose physicochemical properties come very close to the ones of atmospheric particulates. Firstly, we studied the interaction of gaseous NO2 with three different commercial samples of amorphous carbon having widely differing physicochemical properties. Using in situ laser-induced fluorescence, the only major product was unambiguously determined to be NO with the oxidation product of carbon apparently remaining on the surface. Probing the gas phase with mass spectrometry (MS), both pulsed valve dosing and steady state experiments were performed and revealed a complex reaction mechanism for both the uptake as well as product formation. The initial uptake coefficient γo for NO2 was (6.4±2.0)·10-2 and proved to be identical within experimental uncertainty for all three types of amorphous carbon. The initial NO2 uptake rate was independent of the mass of the carbon sample and scaled with its external geometrical surface. Applying a simple chemical kinetic model to both pulsed dosing and steady state experiments resulted in an effective surface for uptake on amorphous carbon ranging between a factor of 2.8 to 8.4 larger than the geometrical surface area, but inversely proportional to the measured BET surface area of the three types of amorphous carbon. This means that the amorphous carbon with the highest BET surface (FW2) had the lowest external surface for NO2 adsorption. The chemical kinetic model included competing processes between Langmuir-type adsorption and inhibition controlling the adsorption kinetics of NO2 and revealed that the NO generation rate differed greatly between the three carbon samples examined. All samples showed saturation effects of differing degree that were partially reversible through prolonged pumping at 10-4 Torr and/or heating. Virgin amorphous carbon samples did not take up H2O vapor at 20 mTorr, and no HONO and/or HNO3 was detected in simultaneous NO2/H2O exposure experiments. CO and CO2 were detected when a sample previously exposed to NO2 was heated by an incandescent lamp. Moreover, upon heating such a sample, a MS signal m/e 62 originating from NO3 and/or N2O5 was detected. In addition to the experiments conducted on the three commercial carbon samples, we carried out uptake experiments of NO2 on acetylene soot. Originating from the flame of an acetylene burner, this soot was freshly deposited on glass dishes before each experiment. The initial uptake coefficient γo for NO2 was measured to be (3.0±1.1)·10-2 is smaller by a factor of two compared to γo observed in the uptake on commercial samples whose mass was a factor of 10 higher. In addition to NO, a net formation of HONO (m/e 47) was observed. This HONO formation is strongly dependent on the instantaneous water desorption from the sample as well as on the NO2 partial pressure in the reactor. On the other hand, it is not enhanced by an external water flow which did not show any interaction with the carbon sample. Upon heating the sample using an incandescent lamp, significant amounts of water desorbed, suggesting that the water vapor had condensed into the micropores of the soot during combustion. No quantitative results are presented yet. However, as no correlation between HONO formation and NO partial pressure has been observed so far, we suggest a reaction mechanism that only involves NO2 and H2O but no NO: 2 NO2 (g) + H2O HONO + HNO3 In a second study, the kinetics of condensation and evaporation of water on ice was determined in the temperature range of 170 to 220K. The rate of evaporation Fev corresponds to the evaporation of 70±10 monolayers of H2O per second at 200K. The condensation rate constant kc has a negative activation energy of -3.l±1.5kcal/mol which is significantly larger than the one recently measured by Haynes et al. We report the first real time kinetic measurements of water condensation on ice at stratospherically relevant temperatures. At temperatures above 160K, pulsed dosing experiments show a dose dependence of the condensation rate coefficient. For example, the condensation rate coefficient increases by a factor of five with a 1000 fold increase of the H2O dose from 1·1015 to 1·1018 molecules per pulse at 180K. This pressure dependence of the condensation rate constant may be explained by an autocatalytic mechanism involving the competition of two condensation channels. These two channels feed two different surface species, one of which leads to autocatalytic condensation. This mechanism gives insight into the manner in which a gas phase molecule is stepwise incorporated into tetrahedrally coordinated bulk ice. In a third and ongoing study, the uptake kinetics of HCl on ice has been investigated by pulsed valve experiments in the temperature range of 180 to 210K and are discussed as a function of the state of the surface. In fact, some of the pulsed valve experiments present a relaxation to a steady state level. It is shown that this level originates from the evaporation of a liquid solution corresponding to a temperature specific HCl concentration. The threshold of the injected dose leading to this solution layer is estimated to be 8·1014 molecules. This threshold corresponding to a fraction of a monolayer may be explained by the formation of small domains along grain boundaries or other lattice imperfections. Within detection limit, no difference in keff has been observed for pulses on the solution or on a previously exposed sample without formation of the solution. In the case of pulses forming a solution layer, a fraction of 20±20% of the injected molecules remains permanently adsorbed on the surface. In the temperature range from 210 to 180K, the pulsed valve experiments yield values of keff roughly between 5 and 25s-1 and an activation energy of -1.6±0.7kcal/mol which are consistent with the values observed by Flueckiger in steady state experiments using the 15mm escape aperture (kuni varying from 12 to 23s-1 and an activation energy of -1.9±0.25 kcal/mol).

EPFL1996

Analyse de la migration d’huile dans les pompes à chaleur mono- et biétagées

Daniel Favrat, Michele Zehnder

The Laboratory of Industrial Energy systems has evaluated a series of high performance air-water heat pump concepts with the aim to retrofit existing fuel fired boilers for residential hydronic heating systems. The results are published within the framework of the Swiss Retrofit Heat Pump project, which is conducted by the Swiss Federal Office of Energy. The two-stage heat pump cycle has shown the best efficiency prospects. Oil migration between the upper stage and lower stage compressors is the main obstacle for commercialization in the low heat rate range. The oil mass concentration is measured by the Fourier Transform Infrared Spectrospopy (FTIR) method with an attenuated total reflection cell (ATR). This method was calibrated on a separated test loop and compared to measurements of mass-sampling and of density by a Coriolis effect densitymeter. Infrared analysis gives a precision of < 0.1% of oil concentration and the detection limit was identified at about 0.3%. Density measurements cannot be used in the ordinary temperature range of the liquid line in the heat pump. Mass sampling and the FT-IR measurements give corresponding results. An air-water heat pump (10 kWth, R-407C) was build in the laboratory to perform comparison tests in the single-stage and two-stage heating modes. The main elements of the heat pump are: two scroll type compressors, a brazed plate condenser, a finned tube evaporator, an economizer heat exchanger and two electronic expansion valves. Defrosting is done by cycle inversion with a 4-way reversing valve. Measured improvements on the two-stage heating mode of 30% for the COP and 24% for the heat output comparing to the corresponding test point (A-7/W50) in the single stage configuration have been shown. Important reductions of the maximum system temperature at the outlet of the compressors enable an extended range of application (including the point A-12/W65). In the two-stage heating mode an oil transfer rate of 5 - 10 g/min from the high pressure stage to the low pressure stage compressor was detected in all tested conditions, which is strongly limiting the continuous service of the heat pump. For long term running tests an additional balancing system composed of an oil pump and an oil vessel has been installed. The circulating oil concentration is between 0.2% and 0.4% for the two-stage configuration as measured with the FT-IR method. In single-stage mode the concentration is lower than 0.2% and only the mass sampling technique provided reliable results. Sensibility analysis in changing the regulation parameters did not give concluding results in order to design a stable configuration of the two-stage heat pump. In absence of integrated two-stage scroll compressors intermediate solutions have to be considered. These consist in the implementation of a booster compressor or the implementation of an oil separator at the outlet of the high pressure stage compressor coupled with a oil return management or a periodic shut down of the heat pump with external balancing of the oil level in the crankcase. An adapted heat transfer model for refrigerant-oil mixtures was developed, in order to investigate the theoretical impact of oil in the evaporator.

2003

Production d'hydrogène dans un réacteur microstructuré

Pierre Reuse

The aim of this work was to operate the steam reforming of methanol and the total oxidation of methanol in a microstructured two-passage reactor. The steam reforming enables the production of hydrogen at relatively low temperatures (250 °C) with only a small amount of carbon monoxide (~ 1%). After a cleaning step (CO preferential oxidation, for example), this hydrogen would be a good feed for a fuel cell to produce electricity for mobile applications. Steam reforming is an endothermic reaction. The energy needed to run this reaction is produced by the total oxidation, a strongly exothermic reaction. To allow an efficient heat transfer between these two reactions, a microstructured reactor was used. These structures, with dimensions in the micron range, are widely known to have far better heat transfer capacity than traditional heat exchangers. Furthermore, these structures have some unique features with respect to classical reactors. A narrow residence time distribution, close to that of a plug flow reactor, can be obtained, keeping selectivity high for complex reactions. Microstructured reactors also have small residence time, allowing very quick variations in the feed and experimental conditions. In chapter 4, the concepts relating to the channel coatings will be presented. A Cu/Zn/Al catalyst provided by Süd-Chemie was transformed into a suspension of nanoparticles in isopropanol by using an attrition grinding apparatus. A thin layer (~5 µm) of catalyst was formed on the microchannel walls. Several tests were performed to characterize the adhesion of this layer. The influence of the layer on the residence time distribution was also investigated; it was possible to maintain Bodenstein numbers as high as 50. Chapter 5 is focused on the steam reforming of methanol, the reaction that generates hydrogen. A kinetic analysis in differential and integral modes established that the rate of the reaction is inhibited by hydrogen, the main product of the reaction; the partial order is - 0.2. The influence of water on the rate of reaction is negligible (partial order of 0.1), and methanol has a positive partial order of 0.7. A kinetic constant of 5.57 · 10-6 mol0.4 · (m3)0.6 · gcat-1 · s-1 at 200 °C and an activation energy of 73.9 kJ · mol-1 complete the description of the kinetics of the steam reforming. Throughout the comparisons, we showed that the catalyst coated in the microchannels has an initial rate for the steam reforming 30% higher than the original catalyst in fixed bed. We can therefore conclude that the catalyst modification procedure does not diminish the catalyst activity. This was also shown by comparing our results with results available in the scientific literature for similar catalysts. CO2 selectivity remains high and constant then decreases as 90% conversion is approached. Working with an excess of water with respect to methanol contributes to keeping the selectivity high. A molar ratio of 1.2 is a good compromise between selectivity and increased energy costs due to the vaporization of the additional water. It's also beneficial to work at conversions below 90% to be in the high selectivity domain and to diminish the inhibiting influence of hydrogen. Finally, we showed that this catalyst deactivates over 100 hours at a temperature of 300 °C. A kinetic model for the deactivation at temperatures between 300 °C and 330 °C was established. In chapter 6 the total oxidation of methanol was studied. A cobalt based catalyst was developed to allow the synchronization of the two reactions temperatures. The minimum temperature for complete conversion for the total oxidation reaction must clearly not exceed an acceptable working temperature for the steam reforming catalyst (i.e. a temperature that results in an acceptable catalyst lifetime). A kinetic analysis determined partial orders for methanol and oxygen as 0.67, the kinetic constant at 230 °C as 5.45 · 10-3 mmol-0.33 · m4 · (gcat · s)-1 and the activation energy as 130 kJ · mol-1. The copper based catalyst used for the steam reforming of methanol was also active for the total oxidation of methanol. Unfortunately the minimal temperature to ensure a total conversion was not compatible with the working temperature of the steam reforming. In chapter 7, the results concerning the coupling of the two previously studied reactions are presented. The numerous fittings and cables necessary for the operation of the reactor didn't allow the reactor to be run in an autothermal mode. Therefore constant, not regulated, heating was used to compensate the thermal losses. By coupling the two reactions and by studying the dynamic response time of the reactor (temperature and conversion of the steam reforming reaction) at the time of stopping and starting the oxidation reaction, we've seen that the temperature increase was around +10 °C within 60 seconds and that the conversion for the steam reforming of the methanol fitted well with the temperature increase. This reactor also has good dynamics and responds quickly to feed variations. When the coupling mode was changed from co-current to counter-current we noticed a significant increase in conversion. Furthermore, selectivity in counter-current mode remained high, even for conversion approaching 100%.

EPFL2003