Inverse convection problem of simultaneous estimation of two boundary heat fluxes in parallel plate channels

This paper deals with the use of the conjugate gradient method of function estimation for the simultaneous identification of two unknown boundary heat fluxes in parallel plate channels. The fluid flow is assumed to be laminar and hydrodynamically developed. Temperature measurements taken inside the channel are used in the inverse analysis. The accuracy of the present solution approach is examined by using simulated measurements containing random errors, for strict cases involving functional forms with discontinuities and sharp-corners for the unknown functions. Three different types of inverse problems are addressed in the paper, involving the estimation of: (i) Spatially dependent heat fluxes; (ii) Time-dependent heat fluxes; and (iii) Time and spatially dependent heat fluxes.

Determination of turbulent heat fluxes using a large aperture scintillometer over undulating mixed agricultural terrain

Scintillometry is an established technique for determining large areal average sensible heat fluxes. The scintillometer measurement is related to sensible heat flux via Monin–Obukhov similarity theory, which was developed for ideal homogeneous land surfaces. In this study it is shown that judicious application of scintillometry over heterogeneous mixed agriculture on undulating topography yields valid results when compared to eddy covariance (EC). A large aperture scintillometer (LAS) over a 2.4 km path was compared with four EC stations measuring sensible (H) and latent (LvE) heat fluxes over different vegetation (cereals and grass) which when aggregated were representative of the LAS source area. The partitioning of available energy into H and LvE varied strongly for different vegetation types, with H varying by a factor of three between senesced winter wheat and grass pasture. The LAS derived H agrees (one-to-one within the experimental uncertainty) with H aggregated from EC with a high coefficient of determination of 0.94. Chronological analysis shows individual fields may have a varying contribution to the areal average sensible heat flux on short (weekly) time scales due to phenological development and changing soil moisture conditions. Using spatially aggregated measurements of net radiation and soil heat flux with H from the LAS, the areal averaged latent heat flux (LvELAS) was calculated as the residual of the surface energy balance. The regression of LvELAS against aggregated LvE from the EC stations has a slope of 0.94, close to ideal, and demonstrates that this is an accurate method for the landscape-scale estimation of evaporation over heterogeneous complex topography.

Multi-scale sensible heat fluxes in the urban environment from large aperture scintillometry and eddy covariance

Sensible heat fluxes (QH) are determined using scintillometry and eddy covariance over a suburban area. Two large aperture scintillometers provide spatially integrated fluxes across path lengths of 2.8 km and 5.5 km over Swindon, UK. The shorter scintillometer path spans newly built residential areas and has an approximate source area of 2-4 km2, whilst the long path extends from the rural outskirts to the town centre and has a source area of around 5-10 km2. These large-scale heat fluxes are compared with local-scale eddy covariance measurements. Clear seasonal trends are revealed by the long duration of this dataset and variability in monthly QH is related to the meteorological conditions. At shorter time scales the response of QH to solar radiation often gives rise to close agreement between the measurements, but during times of rapidly changing cloud cover spatial differences in the net radiation (Q*) coincide with greater differences between heat fluxes. For clear days QH lags Q*, thus the ratio of QH to Q* increases throughout the day. In summer the observed energy partitioning is related to the vegetation fraction through use of a footprint model. The results demonstrate the value of scintillometry for integrating surface heterogeneity and offer improved understanding of the influence of anthropogenic materials on surface-atmosphere interactions.

Infrared and millimetre-wave scintillometry in the suburban environment – Part 2: Large-area sensible and latent heat fluxes

A millimetre-wave scintillometer was paired with an infrared scintillometer, enabling estimation of large-area evapotranspiration across northern Swindon, a suburban area in the UK. Both sensible and latent heat fluxes can be obtained using this "two-wavelength" technique, as it is able to provide both temperature and humidity structure parameters, offering a major advantage over conventional single-wavelength scintillometry. The first paper of this two-part series presented the measurement theory and structure parameters. In this second paper, heat fluxes are obtained and analysed. These fluxes, estimated using two-wavelength scintillometry over an urban area, are the first of their kind. Source area modelling suggests the scintillometric fluxes are representative of 5–10 km2. For comparison, local-scale (0.05–0.5 km2) fluxes were measured by an eddy covariance station. Similar responses to seasonal changes are evident at the different scales but the energy partitioning varies between source areas. The response to moisture availability is explored using data from 2 consecutive years with contrasting rainfall patterns (2011–2012). This extensive data set offers insight into urban surface-atmosphere interactions and demonstrates the potential for two-wavelength scintillometry to deliver fluxes over mixed land cover, typically representative of an area 1–2 orders of magnitude greater than for eddy covariance measurements. Fluxes at this scale are extremely valuable for hydro-meteorological model evaluation and assessment of satellite data products

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.

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.

Experimental investigation of regular fluids and nanofluids during flow boiling in a single microchannel at different heat fluxes and mass fluxes

The dissipation of high heat flux from integrated circuit chips and the maintenance of acceptable junction temperatures in high powered electronics require advanced cooling technologies. One such technology is two-phase cooling in microchannels under confined flow boiling conditions. In macroscale flow boiling bubbles will nucleate on the channel walls, grow, and depart from the surface. In microscale flow boiling bubbles can fill the channel diameter before the liquid drag force has a chance to sweep them off the channel wall. As a confined bubble elongates in a microchannel, it traps thin liquid films between the heated wall and the vapor core that are subject to large temperature gradients. The thin films evaporate rapidly, sometimes faster than the incoming mass flux can replenish bulk fluid in the microchannel. When the local vapor pressure spike exceeds the inlet pressure, it forces the upstream interface to travel back into the inlet plenum and create flow boiling instabilities. Flow boiling instabilities reduce the temperature at which critical heat flux occurs and create channel dryout. Dryout causes high surface temperatures that can destroy the electronic circuits that use two-phase micro heat exchangers for cooling. Flow boiling instability is characterized by periodic oscillation of flow regimes which induce oscillations in fluid temperature, wall temperatures, pressure drop, and mass flux. When nanofluids are used in flow boiling, the nanoparticles become deposited on the heated surface and change its thermal conductivity, roughness, capillarity, wettability, and nucleation site density. It also affects heat transfer by changing bubble departure diameter, bubble departure frequency, and the evaporation of the micro and macrolayer beneath the growing bubbles. Flow boiling was investigated in this study using degassed, deionized water, and 0.001 vol% aluminum oxide nanofluids in a single rectangular brass microchannel with a hydraulic diameter of 229 µm for one inlet fluid temperature of 63°C and two constant flow rates of 0.41 ml/min and 0.82 ml/min. The power input was adjusted for two average surface temperatures of 103°C and 119°C at each flow rate. High speed images were taken periodically for water and nanofluid flow boiling after durations of 25, 75, and 125 minutes from the start of flow. The change in regime timing revealed the effect of nanoparticle suspension and deposition on the Onset of Nucelate Boiling (ONB) and the Onset of Bubble Elongation (OBE). Cycle duration and bubble frequencies are reported for different nanofluid flow boiling durations. The addition of nanoparticles was found to stabilize bubble nucleation and growth and limit the recession rate of the upstream and downstream interfaces, mitigating the spreading of dry spots and elongating the thin film regions to increase thin film evaporation.

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

Sensible and latent heat fluxes in the atmospheric surface layer over the equatorial Atlantic (Appendix)

During a four weeks anchoring station of R.V. ,,Meteor" on the equator at 30° W longitude, vertical profiles of wind, temperature, and humidity were measured by means of a meteorological buoy carrying a mast of 10 m height. After eliminating periods of instrumental failure, 18 days are available for the investigation of the diurnal variations of the meteorological parameters and 9 days for the investigation of the vertical heat fluxes. The diurnal variations of the above mentioned quantities are caused essentially by two periodic processes: the 24-hourly changing solar energy supply and the 12-hourly oscillation of air pressure, which both originate in the daily rotation of the earth. While the temperature of the water and of the near water layers of the air show a 24 hours period in their diurnal course, the wind speed, as a consequence of the pressure wave, has a 12 hours period, which is also observable in evaporation and, consequently, in the water vapor content of the surface layer. Concerning the temperature, a weak dependence of the daily amplitude on height was determined. Further investigation of the profiles yields relations between the vertical gradients of wind, temperature, and water vapor and the wind speed, the difference between sea and air of temperature and water vapor, respectively, thus giving a contribution to the problem of parameterizing the vertical fluxes. Mean profile coefficients for the encountered stabilities, which were slightly unstable, are presented, and correction terms are given due to the fact that the conditions at the very surface are not sufficiently represented by measuring in a water depth of 20 cm and assuming water vapor saturation. This is especially true for the water vapor content, where the relation between the gradient and the air-sea difference suggests a reduction of relative humidity to appr. 96% at the very surface, if the gradients are high. This effect may result in an overestimation of the water vapor flux, if a ,,bulk"-formula is used. Finally sensible and latent heat fluxes are computed by means of a gradient-formula. The influence of stability on the transfer process is taken into account. As the air-sea temperature differences are small, sensible heat plays no important role in that region, but latent heat shows several interesting features. Within the measuring period of 18 days, a regular variation by a factor of ten is observed. Unperiodic short term variations are superposed by periodic diurnal variations. The mean diurnal course shows a 12-hours period caused by the vertical wind speed gradient superposed by a 24-hours period due to the changing stabilities. Mean values within the measuring period are 276 ly/day for latent heat and 9.41y/day for sensible heat.

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.