Numerical Study of Heat Transfer Mechanism in Turbulent Supercritical CO2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

On the thermal budget of Pahoehoe lava flows

In this thesis I investigate some aspects of the thermal budget of pahoehoe lava flows. This is done with a combination of general field observations, quantitative modeling, and specific field experiments. The results of this work apply to pahoehoe flows in general, even though the vast bulk of the work has been conducted on the lavas formed by the Pu'u 'O'o - Kupaianaha eruption of Kilauea Volcano on Hawai'i. The field observations rely heavily on discussions with the staff of the United States Geological Survey's Hawaiian Volcano Observatory (HVO), under whom I labored repeatedly in 1991-1993 for a period totaling about 10 months.

The quantitative models I have constructed are based on the physical processes observed by others and myself to be active on pahoehoe lava flows. By building up these models from the basic physical principles involved, this work avoids many of the pitfalls of earlier attempts to fit field observations with "intuitively appropriate" mathematical expressions. Unlike many earlier works, my model results can be analyzed in terms of the interactions between the different physical processes. I constructed models to: (1) describe the initial cooling of small pahoehoe flow lobes and (2) understand the thermal budget of lava tubes.

The field experiments were designed either to validate model results or to constrain key input parameters. In support of the cooling model for pahoehoe flow lobes, attempts were made to measure: (1) the cooling within the flow lobes, (2) the amount of heat transported away from the lava by wind, and (3) the growth of the crust on the lobes. Field data collected by Jones [1992], Hon et al. [1994b], and Denlinger [Keszthelyi and Denlinger, in prep.] were also particularly useful in constraining my cooling model for flow lobes. Most of the field observations I have used to constrain the thermal budget of lava tubes were collected by HVO (geological and geophysical monitoring) and the Jet Propulsion Laboratory (airborne infrared imagery [Realmuto et al., 1992]). I was able to assist HVO for part of their lava tube monitoring program and also to collect helicopterborne and ground-based IR video in collaboration with JPL [Keszthelyi et al., 1993].

The most significant results of this work are (1) the quantitative demonstration that the emplacement of pahoehoe and 'a'a flows are the fundamentally different, (2) confirmation that even the longest lava flows observed in our Solar System could have formed as low effusion rate, tube-fed pahoehoe flows, and (3) the recognition that the atmosphere plays a very important role throughout the cooling of history of pahoehoe lava flows. In addition to answering specific questions about the thermal budget of tube-fed pahoehoe lava flows, this thesis has led to some additional, more general, insights into the emplacement of these lava flows. This general understanding of the tube-fed pahoehoe lava flow as a system has suggested foci for future research in this part of physical volcanology.

Growth and energy budget of F-2 'all-fish' growth hormone gene transgenic common carp

The growth and energy budget for F-2 'all-fish' growth hormone gene transgenic common carp Cyprinus carpio of two body sizes were investigated at 29.2 degrees C for 21 days. Specific growth rate, feed intake, feed efficiency, digestibility coefficients of dry matter and protein, gross energy intake (I-E), and the proportion of I-E utilized for heat production (H-E) were significantly higher in the transgenics than in the controls. The proportion of I-E directed to waste products [faecal energy (F-E) and excretory energy loss (Z(E) + U-E) where Z(E) is through the gills and U-E through the kidney], and the proportion of metabolizable energy (M-E) for recovered energy (R-E) were significantly lower in the transgenics than in the controls. The average energy budget equation of transgenic fish was as follows: 100 I-E = 19.3 F-E + 6.0 (Z(E) + U-E) + 45.2 H-E + 29.5 R-E or 100 M-E = 60.5 H-E + 39.5 R-E. The average energy budget equation of the controls was: 100 I-E = 25.2 F-E + 7.4 (Z(E) + U-E) + 35.5 H-E + 31.9 R-E or 100 M-E = 52.7 H-E + 47.3 R-E. These findings indicate that the high growth rate of 'all-fish' transgenic common carp relative to their non-transgenic counterparts was due to their increased feed intake, reduced lose of waste productions and improved feed efficiency. The benefit of the increased energy intake by transgenic fish, however, was diminished by their increased metabolism.

Effect of ration on the growth and energy budget of Chinese longsnout catfish, Leiocassis longirostris Gunther

The effect of ration on growth and energy budget of Chinese longsnout catfish was investigated in a growth trial. Fish of initial body weight of 6.5 g were fed at six ration levels (RLs): starvation, 0.8%, 1.6%, 2.4%, 3.2% of body weight per day, and apparent satiation for 8 weeks. Fish were weighed biweekly to adjust feed amount. The results showed that specific growth rate in wet weight, protein and energy increased logarithmically with increased RLs. The relationship of specific growth rate in wet weight (SGRw, % day(-1)) and RL (%) was a decelerating curve: SGRw=-0.62+3.10 Ln(RL+1). The energy budget equation at satiation was: 100 IE=12.94 FE+5.50(ZE+UE)+40.07 HE+41.49 RE, where IE, FE, (ZE+UE), HE, RE are food energy, faecal energy, excretory energy, heat production and recovered energy respectively. Body composition was slightly but significantly affected by ration size except for protein content. The most efficient ration based on the relationship between RL and feed efficiency ratio in energy (FERe) was 1.8% of body weight per day.

Effect of body size on growth and energy budget of Nile tilapia, Oreochromis niloticus

A growth trial was conducted at 30 degrees C to investigate the effect of body size on growth and energy budget of Nile tilapia. The average initial body weights of the four size groups tested were 9.3, 34.1, 80.3 and 172.4 g, respectively. Fish were fed to satiation twice a day with a diet containing 35.6% crude protein. Food consumption (C-max: kJ/day) increased with body size (W: g) according to the relationship: Ln C-max = 1.45 + 0.42 LnW. The final body contents of dry matter, crude protein and ash per unit body weight increased with increasing body size while contents of fat and energy were independent of body size. Specific growth rates of wet weight, dry weight, protein and energy decreased as the fish increased in size. Feed efficiencies in wet weigh, dry weight and crude protein decreased with increasing body size, while that of energy remained unchanged. The proportions of energy intake allocated to the various components (faecal energy, excretory energy, heat production and recovered energy) of the energy budget were not significantly affected by body size, and the average budget was: 100IE-18.5(+/- 1.33)FE + 5.9 (+/- 3.09)(ZE + UE) + 49.3(+/- 3.77)HE + 26.3(+/- 6.23)RE, where IE, FE, (ZE + UE), HE and RE represent gross energy intake, faecal energy, excretory (non-faecal) energy loss, heat production and recovered energy (growth), respectively. It is suggested that the decrease in growth rate in larger fish is mainly due to the decrease in relative food intake. (C) 1997 Elsevier Science B.V.

Energy budget of Nile tilapia (Oreochromis niloticus) in relation to ration size

Nile tilapia weighing 8.29-11.02 g were fed a practical diet at seven ration levels (starvation, 0.5, 1, 2, 3, 4% body weight per day and satiation) twice a day at 30 degrees C. Feed consumption, apparent digestibility, nitrogenous excretion and growth were determined directly, and heat production was calculated by difference of energy budget. The relationship between specific growth rate in wet weight (SGR(w), percentage per day) and ration size (RL, percentage per day) was a decelerating curve described as SGR(w) = 2.98 (1 - e(-0.61(RL-0.43))). The apparent digestibility coefficients for dry matter and protein showed a decreasing pattern with increasing ration while the apparent digestibility coefficient of energy was not significantly affected by ration size. The proportion of gross energy intake lost in nitrogenous excretion tended to decrease with increasing ration. Feed efficiency was highest, and the proportion of gross energy intake channelled to heat production was lowest, at an intermediate ration level (2% per day). The energy budget at the satiation level was: 100IE = 16.9FE + 1.2(ZE + UE) + 62.3HE + 19.6RE, where IE, FE, (ZE + UE), HE and RE represent gross energy intake, faecal energy, excretory (non-faecal) energy loss, heat production and recovered energy (growth), respectively. (C) 1997 Elsevier Science B.V.

Heat balances and thermally driven lagoon-ocean exchangeson a tropical coral reef system (Moorea, French Polynesia)

© 2015. American Geophysical Union. All Rights Reserved.The role of surface and advective heat fluxes on buoyancy-driven circulation was examined within a tropical coral reef system. Measurements of local meteorological conditions as well as water temperature and velocity were made at six lagoon locations for 2 months during the austral summer. We found that temperature rather than salinity dominated buoyancy in this system. The data were used to calculate diurnally phase-averaged thermal balances. A one-dimensional momentum balance developed for a portion of the lagoon indicates that the diurnal heating pattern and consistent spatial gradients in surface heat fluxes create a baroclinic pressure gradient that is dynamically important in driving the observed circulation. The baroclinic and barotropic pressure gradients make up 90% of the momentum budget in part of the system; thus, when the baroclinic pressure gradient decreases 20% during the day due to changes in temperature gradient, this substantially changes the circulation, with different flow patterns occurring during night and day. Thermal balances computed across the entire lagoon show that the spatial heating patterns and resulting buoyancy-driven circulation are important in maintaining a persistent advective export of heat from the lagoon and for enhancing ocean-lagoon exchange.

Climate entropy budget of the HadCM3 atmosphere-ocean general circulation model and of FAMOUS, its low resolution version

The entropy budget is calculated of the coupled atmosphere–ocean general circulation model HadCM3. Estimates of the different entropy sources and sinks of the climate system are obtained directly from the diabatic heating terms, and an approximate estimate of the planetary entropy production is also provided. The rate of material entropy production of the climate system is found to be ∼50 mW m−2 K−1, a value intermediate in the range 30–70 mW m−2 K−1 previously reported from different models. The largest part of this is due to sensible and latent heat transport (∼38 mW m−2 K−1). Another 13 mW m−2 K−1 is due to dissipation of kinetic energy in the atmosphere by friction and Reynolds stresses. Numerical entropy production in the atmosphere dynamical core is found to be about 0.7 mW m−2 K−1. The material entropy production within the ocean due to turbulent mixing is ∼1 mW m−2 K−1, a very small contribution to the material entropy production of the climate system. The rate of change of entropy of the model climate system is about 1 mW m−2 K−1 or less, which is comparable with the typical size of the fluctuations of the entropy sources due to interannual variability, and a more accurate closure of the budget than achieved by previous analyses. Results are similar for FAMOUS, which has a lower spatial resolution but similar formulation to HadCM3, while more substantial differences are found with respect to other models, suggesting that the formulation of the model has an important influence on the climate entropy budget. Since this is the first diagnosis of the entropy budget in a climate model of the type and complexity used for projection of twenty-first century climate change, it would be valuable if similar analyses were carried out for other such models.

Implementation of a new urban energy budget scheme into MetUM. Part II: Validation against observations and model intercomparison

In the first part of this article, we introduced a new urban surface scheme, the Met Office – Reading Urban Surface Exchange Scheme (MORUSES), into the Met Office Unified Model (MetUM) and compared its impact on the surface fluxes with respect to the current urban scheme. In this second part, we aim to analyze further the reasons behind the differences. This analysis is conducted by a comparison of the performance of the two schemes against observations and against a third model, the Single Column Reading Urban model (SCRUM). The key differences between the three models lie in how each model incorporates the heat stored in the urban fabric and how the surface-energy balance is coupled to the underlying substrate. The comparison of the models with observations from Mexico City reveals that the performance of MORUSES is improved if roof insulation is included by minimizing the roof thickness. A comparison of MORUSES and SCRUM reveals that, once insulation is included within MORUSES, these two models perform equally well against the observations overall, but that there are differences in the details of the simulations at the roof and canyon level. These differences are attributed to the different representations of the heat-storage term, specifically differences in the dominant frequencies captured by the urban canopy and substrate, between the models. These results strongly suggest a need for an urban model intercomparison exercise. Copyright © 2010 Royal Meteorological Society and Crown Copyright

The influence of Antarctica on the momentum budget of the southern extratropics

The antarctic plateau acts as a strong heat sink for the global climate, cooling the atmosphere and radiating energy to space. A cold dense atmospheric boundary layer is formed. Strong surface winds are formed as the boundary layer drains off the plateau. These drainage winds and the eddy fluxes necessary to maintain them are analysed in a general circulation model (GCM). The drainage flow is well represented in the GCM. The associated mean meridional circulation is analysed in isentropic coordinates. The momentum budget over Antarctica reveals a balance between the Eliassen-Palm flux convergence and the Coriolis torque exerted by the mean meridional mass flux. Both vertical and horizontal components of the Eliassen-Palm flux contribute, the vertical component being the greater.

Observational constraints on atmospheric and oceanic cross-equatorial heat transports: revisiting the precipitation asymmetry problem in climate models

Satellite based top-of-atmosphere (TOA) and surface radiation budget observations are combined with mass corrected vertically integrated atmospheric energy divergence and tendency from reanalysis to infer the regional distribution of the TOA, atmospheric and surface energy budget terms over the globe. Hemispheric contrasts in the energy budget terms are used to determine the radiative and combined sensible and latent heat contributions to the cross-equatorial heat transports in the atmosphere (AHT_EQ) and ocean (OHT_EQ). The contrast in net atmospheric radiation implies an AHT_EQ from the northern hemisphere (NH) to the southern hemisphere (SH) (0.75 PW), while the hemispheric difference in sensible and latent heat implies an AHT_EQ in the opposite direction (0.51 PW), resulting in a net NH to SH AHT_EQ (0.24 PW). At the surface, the hemispheric contrast in the radiative component (0.95 PW) dominates, implying a 0.44 PW SH to NH OHT_EQ. Coupled model intercomparison project phase 5 (CMIP5) models with excessive net downward surface radiation and surface-to-atmosphere sensible and latent heat transport in the SH relative to the NH exhibit anomalous northward AHT_EQ and overestimate SH tropical precipitation. The hemispheric bias in net surface radiative flux is due to too much longwave surface radiative cooling in the NH tropics in both clear and all-sky conditions and excessive shortwave surface radiation in the SH subtropics and extratropics due to an underestimation in reflection by clouds.

Ocean heat content variability and change in an ensemble of ocean reanalyses

Accurate knowledge of the location and magnitude of ocean heat content (OHC) variability and change is essential for understanding the processes that govern decadal variations in surface temperature, quantifying changes in the planetary energy budget, and developing constraints on the transient climate response to external forcings. We present an overview of the temporal and spatial characteristics of OHC variability and change as represented by an ensemble of dynamical and statistical ocean reanalyses (ORAs). Spatial maps of the 0–300 m layer show large regions of the Pacific and Indian Oceans where the interannual variability of the ensemble mean exceeds ensemble spread, indicating that OHC variations are well-constrained by the available observations over the period 1993–2009. At deeper levels, the ORAs are less well-constrained by observations with the largest differences across the ensemble mostly associated with areas of high eddy kinetic energy, such as the Southern Ocean and boundary current regions. Spatial patterns of OHC change for the period 1997–2009 show good agreement in the upper 300 m and are characterized by a strong dipole pattern in the Pacific Ocean. There is less agreement in the patterns of change at deeper levels, potentially linked to differences in the representation of ocean dynamics, such as water mass formation processes. However, the Atlantic and Southern Oceans are regions in which many ORAs show widespread warming below 700 m over the period 1997–2009. Annual time series of global and hemispheric OHC change for 0–700 m show the largest spread for the data sparse Southern Hemisphere and a number of ORAs seem to be subject to large initialization ‘shock’ over the first few years. In agreement with previous studies, a number of ORAs exhibit enhanced ocean heat uptake below 300 and 700 m during the mid-1990s or early 2000s. The ORA ensemble mean (±1 standard deviation) of rolling 5-year trends in full-depth OHC shows a relatively steady heat uptake of approximately 0.9 ± 0.8 W m−2 (expressed relative to Earth’s surface area) between 1995 and 2002, which reduces to about 0.2 ± 0.6 W m−2 between 2004 and 2006, in qualitative agreement with recent analysis of Earth’s energy imbalance. There is a marked reduction in the ensemble spread of OHC trends below 300 m as the Argo profiling float observations become available in the early 2000s. In general, we suggest that ORAs should be treated with caution when employed to understand past ocean warming trends—especially when considering the deeper ocean where there is little in the way of observational constraints. The current work emphasizes the need to better observe the deep ocean, both for providing observational constraints for future ocean state estimation efforts and also to develop improved models and data assimilation methods.

Modeling study of the aspect ratio influence on urban canopy energy fluxes with a modified wall-canyon energy budget scheme

The influence of the aspect ratio (building height/street canyon width) and the mean building height of cities on local energy fluxes and temperatures is studied by means of an Urban Canopy Model (UCM) coupled with a one-dimensional second-order turbulence closure model. The UCM presented is similar to the Town Energy Balance (TEB) model in most of its features but differs in a few important aspects. In particular, the street canyon walls are treated separately which leads to a different budget of radiation within the street canyon walls. The UCM has been calibrated using observations of incoming global and diffuse solar radiation, incoming long-wave radiation and air temperature at a site in So Paulo, Brazil. Sensitivity studies with various aspect ratios have been performed to assess their impact on urban temperatures and energy fluxes at the top of the canopy layer. In these simulations, it is assumed that the anthropogenic heat flux and latent heat fluxes are negligible. Results show that the simulated net radiation and sensible heat fluxes at the top of the canopy decrease and the stored heat increases as the aspect ratio increases. The simulated air temperature follows the behavior of the sensible heat flux. (C) 2010 Elsevier Ltd. All rights reserved.

Formulation of a tropical town energy budget (t-TEB) scheme

This work describes the tropical town energy budget (t-TEB) scheme addressed to simulate the diurnal occurrence of the urban heat island (UHI) as observed in the Metropolitan Area of Rio de Janeiro (MARJ; -22A degrees S; -44A degrees W) in Brazil. Reasoning about the tropical urban climate have guided the scheme implementation, starting from the original equations from Masson (Bound-Lay Meteorol 94:357-397, 2000). The modifications include (a) local scaling approaches for obtaining flux-gradient relationships in the roughness sub-layer, (b) the Monin-Obukhov similarity framework in the inertial sub-layer, (c) increasing aerodynamic conductance toward more unstable conditions, and (d) a modified urban subsurface drainage system to transfer the intercepted rainwater by roofs to the roads. Simulations along 2007 for the MARJ are obtained and compared with the climatology. The t-TEB simulation is consistent with the observations, suggesting that the timing and dynamics of the UHI in tropical cities could vary significantly from the familiar patterns observed in mid-latitude cities-with the peak heat island intensity occurring in the morning than at night. The simulations are suggesting that the thermal phase shift of this tropical diurnal UHI is a response of the surface energy budget to the large amount of solar radiation, intense evapotranspiration, and thermal response of the vegetated surfaces over a very humid soil layer.

Upwelling and associated heat flux in the equatorial Atlantic inferred from helium isotope disequilibrium

Sea-to-air and diapycnal fluxes of nitrous oxide (N2O) into the mixed layer were determined during three cruises to the upwelling region off Mauritania. Sea-to-air fluxes as well as diapycnal fluxes were elevated close to the shelf break, but elevated sea-to-air fluxes reached further offshore as a result of the offshore transport of upwelled water masses. To calculate a mixed layer budget for N2O we compared the regionally averaged sea-to-air and diapycnal fluxes and estimated the potential contribution of other processes, such as vertical advection and biological N2O production in the mixed layer. Using common parameterizations for the gas transfer velocity, the comparison of the average sea-toair and diapycnal N2O fluxes indicated that the mean sea-toair flux is about three to four times larger than the diapycnal flux. Neither vertical and horizontal advection nor biological production were found sufficient to close the mixed layer budget. Instead, the sea-to-air flux, calculated using a parameterization that takes into account the attenuating effect of surfactants on gas exchange, is in the same range as the diapycnal flux. From our observations we conclude that common parameterizations for the gas transfer velocity likely overestimate the air-sea gas exchange within highly productive upwelling zones.