The simulation of the opposing fluxes of latent heat and CO2 over various land-use types: coupling a gas exchange model to a mesoscale atmospheric model

A mesoscale meteorological model (FOOT3DK) is coupled with a gas exchange model to simulate surface fluxes of CO2 and H2O under field conditions. The gas exchange model consists of a C3 single leaf photosynthesis sub-model and an extended big leaf (sun/shade) sub-model that divides the canopy into sunlit and shaded fractions. Simulated CO2 fluxes of the stand-alone version of the gas exchange model correspond well to eddy-covariance measurements at a test site in a rural area in the west of Germany. The coupled FOOT3DK/gas exchange model is validated for the diurnal cycle at singular grid points, and delivers realistic fluxes with respect to their order of magnitude and to the general daily course. Compared to the Jarvis-based big leaf scheme, simulations of latent heat fluxes with a photosynthesis-based scheme for stomatal conductance are more realistic. As expected, flux averages are strongly influenced by the underlying land cover. While the simulated net ecosystem exchange is highly correlated with leaf area index, this correlation is much weaker for the latent heat flux. Photosynthetic CO2 uptake is associated with transpirational water loss via the stomata, and the resulting opposing surface fluxes of CO2 and H2O are reproduced with the model approach. Over vegetated surfaces it is shown that the coupling of a photosynthesis-based gas exchange model with the land-surface scheme of a mesoscale model results in more realistic simulated latent heat fluxes.

An assessment of air-sea heat fluxes from ocean and coupled reanalyses

Sixteen monthly air–sea heat flux products from global ocean/coupled reanalyses are compared over 1993–2009 as part of the Ocean Reanalysis Intercomparison Project (ORA-IP). Objectives include assessing the global heat closure, the consistency of temporal variability, comparison with other flux products, and documenting errors against in situ flux measurements at a number of OceanSITES moorings. The ensemble of 16 ORA-IP flux estimates has a global positive bias over 1993–2009 of 4.2 ± 1.1 W m−2. Residual heat gain (i.e., surface flux + assimilation increments) is reduced to a small positive imbalance (typically, +1–2 W m−2). This compensation between surface fluxes and assimilation increments is concentrated in the upper 100 m. Implied steady meridional heat transports also improve by including assimilation sources, except near the equator. The ensemble spread in surface heat fluxes is dominated by turbulent fluxes (>40 W m−2 over the western boundary currents). The mean seasonal cycle is highly consistent, with variability between products mostly <10 W m−2. The interannual variability has consistent signal-to-noise ratio (~2) throughout the equatorial Pacific, reflecting ENSO variability. Comparisons at tropical buoy sites (10°S–15°N) over 2007–2009 showed too little ocean heat gain (i.e., flux into the ocean) in ORA-IP (up to 1/3 smaller than buoy measurements) primarily due to latent heat flux errors in ORA-IP. Comparisons with the Stratus buoy (20°S, 85°W) over a longer period, 2001–2009, also show the ORA-IP ensemble has 16 W m−2 smaller net heat gain, nearly all of which is due to too much latent cooling caused by differences in surface winds imposed in ORA-IP.

Alternating heat transfer between objects of large geometrical difference

The heat transfer on the surface of an object in a gas fluidised bed is sequentially and alternately induced by particle-packet and gas bubble. This phenomenon is studied with computational simulation. The particle-packet and bubble are modelled by a double particle-layers and porous medium model and a hemispherical model, respectively. The heat transfer to and within the object is simulated concurrently. Different grid schemes are applied and different grid sizes are used in meshing the particle-packet and the object as there is a very large difference in their geometrical sizes. Based on theoretical analysis, an approximate method is developed to calculate the heat flux at the surface of the object. The simulation is implemented in a CFD package and the results are compared with experiments.

Study of heat flow through a rammed earth wall building

In 1999, a 2100 m2 (GFA) two-storey rammed earth building was built on the Thurgoona campus of Charles Sturt University. The climate at Thurgoona is considered Mediterranean – hot dry summers and cool winters. The internal and external walls of the building are constructed from 300-mm thick rammed earth (pise) and are load bearing. The thermal performance of the building has been investigated, both experimentally and theoretically over the summer and winter seasons of 2000/1. As part of these investigations heat flux sensors and thermistors were embedded in one of the external walls of a ground floor office, and data from the transducers has been used to determine the heat flow at the internal and external wall surfaces. The simulation software, TRNSYS, has been used to model the thermal performance of the same office. The programme allows the user to calculate the heat flow at the walls, which define any particular thermal zone. A comparison of measured and predicted values of heat flows and air temperatures has been used to validate the model. The model has then been used to simulate the effect of shading and added insulation on the thermal performance of the external walls in both summer and winter and these results are also presented in this paper.

Particle dynamics and heat transfer at workpiece surface in heat treatment fluidised beds

The particle behaviour is studied by the analysis of particle images taken with a high speed CCD digital video camera. The comparison of particle dynamics is performed for the fluidised beds without part, with single part and with multi-parts. The results show that there are significant differences in particle behaviours both in different beds and at different locations at part surfaces. The total and radiative heat transfer coefficients at different surfaces of a metallic component in a high temperature fluidised bed are measured by a heat transfer probe developed in the present work. The principle of the heat transfer probe is to measure the change in temperature of the heated metallic piece with time and, then, to extract the heat flux and heat transfer coefficients. The structure of the probe is optimized with numerical simulation of energy conservation for measuring the heat transfer coefficient of 150~600 W/m2 K. The relationship between the particle dynamics and the heat transfer is analysed to form the basis for future more rational designs of fluidised beds as well as for improved quality control.

Effects of cold fronts on MODIS-derived sensible and latent heat fluxes in Itumbiara reservoir (Central Brazil)

In this work we investigate the cold front passage effects on sensible and latent heat flux in a tropical hydroelectric reservoir. The study area, Itumbiara reservoir (Goiás State/Brazil) at the beginning of the austral winter, is characterized by the presence of a weak thermal stratification and the passage of several cold fronts from higher latitudes of South America. Sensible and latent heat fluxes were estimated considering the atmospheric boundary layer stability. In situ and MODIS water surface temperature data were used to adjust the coefficients for momentum and heat exchanges between water and atmosphere and spatialize the sensible and latent heat fluxes. The results showed that during a cold front event the sensible heat flux can be up to five times greater than the flux observed before. The latent heat flux tends to decrease during the cold front but increase again after the passage. The highest values of heat loss were observed at littoral zone and some Reservoir's embayment. The heat loss intensification can be separated in two moments: first, during the cold front passage, when the wind speed increases and the air temperature decreases; second, after the cold front passage, with air humidity decreasing. This can be considered a key process to understanding the heat loss in the Itumbiara reservoir. © 2013 COSPAR. Published by Elsevier Ltd. All rights reserved.

Influence of summertime mesoscale convective systems on the heat balance and surface mixed layer dynamics of a large Amazonian hydroelectric reservoir

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

(Table 1) Sea ice and snow characteristics and heat fluxes observed during R/V Aurora Australis cruise to East Antarctica in September/October 2007

The properties of snow on East Antarctic sea ice off Wilkes Land were examined during the Sea Ice Physics and Ecosystem Experiment (SIPEX) in late winter of 2007, focusing on the interaction with sea ice. This observation includes 11 transect lines for the measurement of ice thickness, freeboard, and snow depth, 50 snow pits on 13 ice floes, and diurnal variation of surface heat flux on three ice floes. The detailed profiling of topography along the transects and the d18O, salinity, and density datasets of snow made it possible to examine the snow-sea-ice interaction quantitatively for the first time in this area. In general, the snow displayed significant heterogeneity in types, thickness (mean: 0.14 +- 0.13 m), and density (325 +- 38 kg/m**3), as reported in other East Antarctic regions. High salinity was confined to the lowest 0.1 m. Salinity and d18O data within this layer revealed that saline water originated from the surface brine of sea ice in 20% of the total sites and from seawater in 80%. From the vertical profiles of snow density, bulk thermal conductivity of snow was estimated as 0.15 W/K/m on average, only half of the value used for numerical sea-ice models. Although the upward heat flux within snow estimated with this value was significantly lower than that within ice, it turned out that a higher value of thermal conductivity (0.3 to 0.4 W/K/m) is preferable for estimating ice growth amount in current numerical models. Diurnal measurements showed that upward conductive heat flux within the snow and net long-wave radiation at the surface seem to play important roles in the formation of snow ice from slush. The detailed surface topography allowed us to compare the air-ice drag coefficients of ice and snow surfaces under neutral conditions, and to examine the possibility of the retrieval of ice thickness distribution from satellite remote sensing. It was found that overall snow cover works to enhance the surface roughness of sea ice rather than moderate it, and increases the drag coefficient by about 10%. As for thickness retrieval, mean ice thickness had a higher correlation with ice surface roughness than mean freeboard or surface elevation, which indicates the potential usefulness of satellite L-band SAR in estimating the ice thickness distribution in the seasonal sea-ice zone.

Heat fluxes and surface radiation measurements from Samoylov Island, Lena Delta, Siberia, April to September 2007

In this article, we present a study on the surface energy balance of a polygonal tundra landscape in northeast Siberia. The study was performed during half-year periods from April to September in each of 2007 and 2008. The surface energy balance is obtained from independent measurements of the net radiation, the turbulent heat fluxes, and the ground heat flux at several sites. Short-wave radiation is the dominant factor controlling the magnitude of all the other components of the surface energy balance during the entire observation period. About 50% of the available net radiation is consumed by the latent heat flux, while the sensible and the ground heat flux are each around 20 to 30%. The ground heat flux is mainly consumed by active layer thawing. About 60% of the energy storage in the ground is attributed to the phase change of soil water. The remainder is used for soil warming down to a depth of 15 m. In particular, the controlling factors for the surface energy partitioning are snow cover, cloud cover, and the temperature gradient in the soil. The thin snow cover melts within a few days, during which the equivalent of about 20% of the snow-water evaporates or sublimates. Surface temperature differences of the heterogeneous landscape indicate spatial variabilities of sensible and latent heat fluxes, which are verified by measurements. However, spatial differences in the partitioning between sensible and latent heat flux are only measured during conditions of high radiative forcing, which only occur occasionally.

Mixed layer heat and salinity budgets during the onset of the 2011 Atlantic cold tongue

The mixed layer (ML) temperature and salinity changes in the central tropical Atlantic have been studied by a dedicated experiment (Cold Tongue Experiment (CTE)) carried out from May to July 2011. The CTE was based on two successive research cruises, a glider swarm, and moored observations. The acquired in situ data sets together with satellite, reanalysis, and assimilation model data were used to evaluate box-averaged ML heat and salinity budgets for two subregions: (1) the western equatorial Atlantic cold tongue (ACT) (23°-10°W) and (2) the region north of the ACT. The strong ML heat loss in the ACT region during the CTE was found to be the result of the balance of warming due to net surface heat flux and cooling due to zonal advection and diapycnal mixing. The northern region was characterized by weak cooling and the dominant balance of net surface heat flux and zonal advection. A strong salinity increase occurred at the equator, 10°W, just before the CTE. During the CTE, ML salinity in the ACT region slightly increased. Largest contributions to the ML salinity budget were zonal advection and the net surface freshwater flux. While essential for the ML heat budget in the ACT region, diapycnal mixing played only a minor role for the ML salinity budget. In the region north of the ACT, the ML freshened at the beginning of the CTE due to precipitation, followed by a weak salinity increase. Zonal advection changed sign contributing to ML freshening at the beginning of the CTE and salinity increase afterward.

Difference schemes for numerical solutions of lagging models of heat conduction

Non-Fourier models of heat conduction are increasingly being considered in the modeling of microscale heat transfer in engineering and biomedical heat transfer problems. The dual-phase-lagging model, incorporating time lags in the heat flux and the temperature gradient, and some of its particular cases and approximations, result in heat conduction modeling equations in the form of delayed or hyperbolic partial differential equations. In this work, the application of difference schemes for the numerical solution of lagging models of heat conduction is considered. Numerical schemes for some DPL approximations are developed, characterizing their properties of convergence and stability. Examples of numerical computations are included.

Difference schemes for time-dependent heat conduction models with delay

Different non-Fourier models of heat conduction, that incorporate time lags in the heat flux and/or the temperature gradient, have been increasingly considered in the last years to model microscale heat transfer problems in engineering. Numerical schemes to obtain approximate solutions of constant coefficients lagging models of heat conduction have already been proposed. In this work, an explicit finite difference scheme for a model with coefficients variable in time is developed, and their properties of convergence and stability are studied. Numerical computations showing examples of applications of the scheme are presented.

A compact difference scheme for numerical solutions of second order dual-phase-lagging models of microscale heat transfer

Dual-phase-lagging (DPL) models constitute a family of non-Fourier models of heat conduction that allow for the presence of time lags in the heat flux and the temperature gradient. These lags may need to be considered when modeling microscale heat transfer, and thus DPL models have found application in the last years in a wide range of theoretical and technical heat transfer problems. Consequently, analytical solutions and methods for computing numerical approximations have been proposed for particular DPL models in different settings. In this work, a compact difference scheme for second order DPL models is developed, providing higher order precision than a previously proposed method. The scheme is shown to be unconditionally stable and convergent, and its accuracy is illustrated with numerical examples.

Heat fluxes and surface radiation measurements from Samoylov Island, Lena Delta, Siberia, October 2007 to March 2008

In this study, we present the winter time surface energy balance at a polygonal tundra site in northern Siberia based on independent measurements of the net radiation, the sensible heat flux and the ground heat flux from two winter seasons. The latent heat flux is inferred from measurements of the atmospheric turbulence characteristics and a model approach. The long-wave radiation is found to be the dominant factor in the surface energy balance. The radiative losses are balanced to about 60 % by the ground heat flux and almost 40 % by the sensible heat fluxes, whereas the contribution of the latent heat flux is small. The main controlling factors of the surface energy budget are the snow cover, the cloudiness and the soil temperature gradient. Large spatial differences in the surface energy balance are observed between tundra soils and a small pond. The ground heat flux released at a freezing pond is by a factor of two higher compared to the freezing soil, whereas large differences in net radiation between the pond and soil are only observed at the end of the winter period. Differences in the surface energy balance between the two winter seasons are found to be related to differences in snow depth and cloud cover which strongly affect the temperature evolution and the freeze-up at the investigated pond.

TWO-PHASE HEAT TRANSFER MECHANISMS WITHIN PLATE HEAT EXCHANGERS: EXPERIMENTS AND MODELING

Two-phase flow heat exchangers have been shown to have very high efficiencies, but the lack of a dependable model and data precludes them from use in many cases. Herein a new method for the measurement of local convective heat transfer coefficients from the outside of a heat transferring wall has been developed, which results in accurate local measurements of heat flux during two-phase flow. This novel technique uses a chevron-pattern corrugated plate heat exchanger consisting of a specially machined Calcium Fluoride plate and the refrigerant HFE7100, with heat flux values up to 1 W cm-2 and flow rates up to 300 kg m-2s-1. As Calcium Fluoride is largely transparent to infra-red radiation, the measurement of the surface temperature of PHE that is in direct contact with the liquid is accomplished through use of a mid-range (3.0-5.1 µm) infra-red camera. The objective of this study is to develop, validate, and use a unique infrared thermometry method to quantify the heat transfer characteristics of flow boiling within different Plate Heat Exchanger geometries. This new method allows high spatial and temporal resolution measurements. Furthermore quasi-local pressure measurements enable us to characterize the performance of each geometry. Validation of this technique will be demonstrated by comparison to accepted single and two-phase data. The results can be used to come up with new heat transfer correlations and optimization tools for heat exchanger designers. The scientific contribution of this thesis is, to give PHE developers further tools to allow them to identify the heat transfer and pressure drop performance of any corrugated plate pattern directly without the need to account for typical error sources due to inlet and outlet distribution systems. Furthermore, the designers will now gain information on the local heat transfer distribution within one plate heat exchanger cell which will help to choose the correct corrugation geometry for a given task.