Extent and economic effect of heat loads on dairy cattle production in Australia

Objective To investigate the extent of heat load problems, caused by the combination of excessive temperature and humidity, in Holstein-Friesian cows in Australia. Also, to outline how milk production losses and consequent costs from this can be estimated and minimised. Procedures Long-term meteorological data for Australia were analysed to determine the distribution of hot conditions over space and time. Fifteen dairy production regions were identified for higher-resolution data analysis. Both the raw meteorological data and their integration into a temperature-humidity thermal index were compiled onto a computer program. This mapping software displays the distribution of climatic patterns, both Australia-wide and within the selected dairying regions. Graphical displays of the variation in historical records for 200 locations in the 15 dairying regions are also available. As a separate study, production data from research stations, on-farm trials and milk factory records were statistically analysed and correlated with the climatic indices, to estimate production losses due to hot conditions. Results Both milk yields and milk constituents declined with increases in the temperature-humidity index. The onset and rate of this decline are dependent on a number of factors, including location, level of production, adaptation, and management regime. These results have been integrated into a farm-level economic analysis for managers of dairy properties. Conclusion By considering the historical patterns of hot conditions over time and space, along with expected production losses, managers of dairy farms can now conduct an economic evaluation of investment strategies to alleviate heat loads. These strategies include the provision of sprinklers, shade structures, or combinations of these.

Effects of environmental heat on conception rates in lactating dairy cows: Critical periods of exposure

Environmental heat can reduce conception rates (the proportion of services that result in pregnancy) in lactating dairy cows. The study objectives were to identify periods of exposure relative to the service date in which environmental heat is most closely associated with conception rates, and to assess whether the total time cows are exposed to high environmental heat within each 24-h period is more closely associated with conception rates than is the maximum environmental heat for each 24-h period. A retrospective observational study was conducted in 25 predominantly Holstein-Friesian commercial dairy herds located in Australia. Associations between weather and conception rates were assessed using 16,878 services performed over a 21-mo period. Services were classified as successful based on rectal palpation. Two measures of heat load were defined for each 24-h period: the maximum temperature-humidity index (THI) for the period, and the number of hours in the 24-h period when the THI was >72. Conception rates were reduced when cows were exposed to a high heat load from the day of service to 6 d after service, and in wk -1. Heat loads in wk -3 to -5 were also associated with reduced conception rates. Thus, management interventions to ameliorate the effects of heat load on conception rates should be implemented at least 5 wk before anticipated service and should continue until at least 1 wk after service. High autocorrelations existed between successive daily values in both measures, and associations between day of heat load relative to service day and conception rates differed substantially when ridge regression was used to account for this autocorrelation. This indicates that when assessing the effects of heat load on conception rates, the autocorrelation in heat load between days should be accounted for in analyses. The results suggest that either weekly averages or totals summarizing the daily heat load are adequate to describe heat load when assessing effects on conception rates in lactating dairy cows.

Theoretical and Experimental Studies on an Ammonia-Based Loop Heat Pipe With a Flat Evaporator

An ammonia loop heat pipe (LHP) with a flat plate evaporator is developed and tested. The device uses a nickel wick encased in an aluminum-stainless steel casing. The loop is tested for various heat loads and different sink temperatures, and it demonstrated reliable startup characteristics. Results with the analysis of the experimental observation indicate that the conductance between the compensation chamber and the heater plate can significantly influence the operating temperatures of the LHP. A mathematical model is also presented which is validated against the experimental observations.

Performance Analysis of Heat Sinks With Phase-Change Materials Subjected to Transient and Cyclic Heating

Phase-change cooling technique is a suitable method for thermal management of electronic equipment subjected to transient or cyclic heat loads. The thermal performance of a phase-change based heat sink under cyclic heat load depends on several design parameters, namely, applied heat flux, cooling heat transfer coefficient, thermophysical properties of phase-change materials (PCMs), and physical dimensions of phase-change storage system during melting and freezing processes. A one-dimensional conduction heat transfer model is formulated to evaluate the effectiveness of preliminary design of practical PCM-based energy storage units. In this model, the phase-change process of the PCM is divided into melting and solidification subprocesses, for which separate equations are written. The equations are solved sequentially and an explicit closed-form solution is obtained. The efficacy of analytical model is estimated by comparing with a finite-volume-based numerical solution for both transient and cyclic heat loads.

Developing an enhanced model for combined heat and air infiltration energy simulation

The need for efficient, sustainable, and planned utilization of resources is ever more critical. In the U.S. alone, buildings consume 34.8 Quadrillion (1015) BTU of energy annually at a cost of $1.4 Trillion. Of this energy 58% is utilized for heating and air conditioning. ^ Several building energy analysis tools have been developed to assess energy demands and lifecycle energy costs in buildings. Such analyses are also essential for an efficient HVAC design that overcomes the pitfalls of an under/over-designed system. DOE-2 is among the most widely known full building energy analysis models. It also constitutes the simulation engine of other prominent software such as eQUEST, EnergyPro, PowerDOE. Therefore, it is essential that DOE-2 energy simulations be characterized by high accuracy. ^ Infiltration is an uncontrolled process through which outside air leaks into a building. Studies have estimated infiltration to account for up to 50% of a building's energy demand. This, considered alongside the annual cost of buildings energy consumption, reveals the costs of air infiltration. It also stresses the need that prominent building energy simulation engines accurately account for its impact. ^ In this research the relative accuracy of current air infiltration calculation methods is evaluated against an intricate Multiphysics Hygrothermal CFD building envelope analysis. The full-scale CFD analysis is based on a meticulous representation of cracking in building envelopes and on real-life conditions. The research found that even the most advanced current infiltration methods, including in DOE-2, are at up to 96.13% relative error versus CFD analysis. ^ An Enhanced Model for Combined Heat and Air Infiltration Simulation was developed. The model resulted in 91.6% improvement in relative accuracy over current models. It reduces error versus CFD analysis to less than 4.5% while requiring less than 1% of the time required for such a complex hygrothermal analysis. The algorithm used in our model was demonstrated to be easy to integrate into DOE-2 and other engines as a standalone method for evaluating infiltration heat loads. This will vastly increase the accuracy of such simulation engines while maintaining their speed and ease of use characteristics that make them very widely used in building design.^

Developing an Enhanced Model for Combined Heat and Air Infiltration Energy Simulation

The need for efficient, sustainable, and planned utilization of resources is ever more critical. In the U.S. alone, buildings consume 34.8 Quadrillion (1015) BTU of energy annually at a cost of $1.4 Trillion. Of this energy 58% is utilized for heating and air conditioning. Several building energy analysis tools have been developed to assess energy demands and lifecycle energy costs in buildings. Such analyses are also essential for an efficient HVAC design that overcomes the pitfalls of an under/over-designed system. DOE-2 is among the most widely known full building energy analysis models. It also constitutes the simulation engine of other prominent software such as eQUEST, EnergyPro, PowerDOE. Therefore, it is essential that DOE-2 energy simulations be characterized by high accuracy. Infiltration is an uncontrolled process through which outside air leaks into a building. Studies have estimated infiltration to account for up to 50% of a building’s energy demand. This, considered alongside the annual cost of buildings energy consumption, reveals the costs of air infiltration. It also stresses the need that prominent building energy simulation engines accurately account for its impact. In this research the relative accuracy of current air infiltration calculation methods is evaluated against an intricate Multiphysics Hygrothermal CFD building envelope analysis. The full-scale CFD analysis is based on a meticulous representation of cracking in building envelopes and on real-life conditions. The research found that even the most advanced current infiltration methods, including in DOE-2, are at up to 96.13% relative error versus CFD analysis. An Enhanced Model for Combined Heat and Air Infiltration Simulation was developed. The model resulted in 91.6% improvement in relative accuracy over current models. It reduces error versus CFD analysis to less than 4.5% while requiring less than 1% of the time required for such a complex hygrothermal analysis. The algorithm used in our model was demonstrated to be easy to integrate into DOE-2 and other engines as a standalone method for evaluating infiltration heat loads. This will vastly increase the accuracy of such simulation engines while maintaining their speed and ease of use characteristics that make them very widely used in building design.

Solar thermal power system augmented with LHP

The present work reports the study of the bubble formation dynamics in the compensation chamber (CC) of the evaporator in Loop Heat Pipes. A series of experiments were conducted at different heat loads and bubbles in the CC were visualized. Bubbles diameter, frequency and velocity were measured and correlated against heat loads. Temperatures were measured at various locations and heat transfer coefficient was calculated. Performance of the LHP evaporator was evaluated at different heat loads. (C) 2013 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license

Two-Phase Flow With Phase Change in Porous Channels

Phase change heat transfer in porous media finds applications in various geological flows and modern heat pipes. We present a study to show the effect of phase change on heat transfer in a porous channel. We show that the ratio of Jakob numbers based on wall superheat and inlet fluid subcooling governs the liquid-vapor interface location in the porous channel and below a critical value of the ratio, the liquid penetrates all the way to the extent of the channel in the flow direction. In such cases, the Nusselt number is higher due to the proximity of the liquid-vapor interface to the heat loads. For higher heat loads or lower subcooling of the liquid, the liquid-vapor interface is pushed toward the inlet, and heat transfer occurs through a wider vapor region thus resulting in a lower Nusselt number. This study is relevant in the designing of efficient two-phase heat exchangers such as capillary suction based heat pipes where a prior estimation of the interface location for the maximum heat load is required to ensure that the liquid-vapor interface is always inside the porous block for its operation.

ELM-induced transient tungsten melting in the JET divertor

The original goals of the JET ITER-like wall included the study of the impact of an all W divertor on plasma operation (Coenen et al 2013 Nucl. Fusion 53 073043) and fuel retention (Brezinsek et al 2013 Nucl. Fusion 53 083023). ITER has recently decided to install a full-tungsten (W) divertor from the start of operations. One of the key inputs required in support of this decision was the study of the possibility of W melting and melt splashing during transients. Damage of this type can lead to modifications of surface topology which could lead to higher disruption frequency or compromise subsequent plasma operation. Although every effort will be made to avoid leading edges, ITER plasma stored energies are sufficient that transients can drive shallow melting on the top surfaces of components. JET is able to produce ELMs large enough to allow access to transient melting in a regime of relevance to ITER.

Transient W melt experiments were performed in JET using a dedicated divertor module and a sequence of I-P = 3.0 MA/B-T = 2.9 T H-mode pulses with an input power of P-IN = 23 MW, a stored energy of similar to 6 MJ and regular type I ELMs at Delta W-ELM = 0.3 MJ and f(ELM) similar to 30 Hz. By moving the outer strike point onto a dedicated leading edge in the W divertor the base temperature was raised within similar to 1 s to a level allowing transient, ELM-driven melting during the subsequent 0.5 s. Such ELMs (delta W similar to 300 kJ per ELM) are comparable to mitigated ELMs expected in ITER (Pitts et al 2011 J. Nucl. Mater. 415 (Suppl.) S957-64).

Although significant material losses in terms of ejections into the plasma were not observed, there is indirect evidence that some small droplets (similar to 80 mu m) were released. Almost 1 mm (similar to 6 mm(3)) of W was moved by similar to 150 ELMs within 7 subsequent discharges. The impact on the main plasma parameters was minor and no disruptions occurred. The W-melt gradually moved along the leading edge towards the high-field side, driven by j x B forces. The evaporation rate determined from spectroscopy is 100 times less than expected from steady state melting and is thus consistent only with transient melting during the individual ELMs. Analysis of IR data and spectroscopy together with modelling using the MEMOS code Bazylev et al 2009 J. Nucl. Mater. 390-391 810-13 point to transient melting as the main process. 3D MEMOS simulations on the consequences of multiple ELMs on damage of tungsten castellated armour have been performed.

These experiments provide the first experimental evidence for the absence of significant melt splashing at transient events resembling mitigated ELMs on ITER and establish a key experimental benchmark for the MEMOS code.

The spatial distribution of carbon dioxide in an environmental test chamber

The spatial distribution of CO2 level in a classroom carried out in previous field work research has demonstrated that there is some evidence of variations in CO2 concentration in a classroom space. Significant fluctuations in CO2 concentration were found at different sampling points depending on the ventilation strategies and environmental conditions prevailing in individual classrooms. However, how these variations are affected by the emitting sources and the room air movement remains unknown. Hence, it was concluded that detailed investigation of the CO2 distribution need to be performed on a smaller scale. As a result, it was decided to use an environmental chamber with various methods and rates of ventilation, for the same internal temperature and heat loads, to study the effect of ventilation strategy and air movement on the distribution of CO2 concentration in a room. The role of human exhalation and its interaction with the plume induced by the body's convective flow and room air movement due to different ventilation strategies were studied in a chamber at the University of Reading. These phenomena are considered to be important in understanding and predicting the flow patterns in a space and how these impact on the distribution of contaminants. This paper attempts to study the CO2 dispersion and distribution at the exhalation zone of two people sitting in a chamber as well as throughout the occupied zone of the chamber. The horizontal and vertical distributions of CO2 were sampled at locations with a probability that CO2 variation is considered high. Although the room size, source location, ventilation rate and location of air supply and extract devices all can have influence on the CO2 distribution, this article gives general guidelines on the optimum positioning of CO2 sensor in a room.

Evaluation of a high temperature solar thermal seasonal borehole storage

A solar thermal system with seasonal borehole storage for heating of a residential area in Anneberg, Sweden, approximately 10 km north of Stockholm, has been in operation since late 2002. Originally, the project was part of the EU THERMIE project “Large-scale Solar Heating Systems for Housing Developments” (REB/0061/97) and was the first solar heating plant in Europe with borehole storage in rock not utilizing a heat pump. Earlier evaluations of the system show lower performance than the preliminary simulation study, with residents complaining of a high use of electricity for domestic hot water (DHW) preparation and auxiliary heating. One explanation mentioned in the earlier evaluations is that the borehole storage had not yet reached “steady state” temperatures at the time of evaluation. Many years have passed since then and this paper presents results from a new evaluation. The main aim of this work is to evaluate the current performance of the system based on several key figures, as well as on system function based on available measurement data. The analysis show that though the borehole storage now has reached a quasi-steady state and operates as intended, the auxiliary electricity consumption is much higher than the original design values largely due to high losses in the distribution network, higher heat loads as well as lower solar gains.

Radiation properties of coil-coated steel in building envelope surfaces and the influence on building thermal performance

Recent studies have shown that the optical properties of building exterior surfaces are important in terms of energy use and thermal comfort. While the majority of the studies are related to exterior surfaces, the radiation properties of interior surfaces are less thoroughly investigated. Development in the coil-coating industries has now made it possible to allocate different optical properties for both exterior and interior surfaces of steel-clad buildings. The aim of this thesis is to investigate the influence of surface radiation properties with the focus on the thermal emittance of the interior surfaces, the modeling approaches and their consequences in the context of the building energy performance and indoor thermal environment. The study consists of both numerical and experimental investigations. The experimental investigations include parallel field measurements on three similar test cabins with different interior and exterior surface radiation properties in Borlänge, Sweden, and two ice rink arenas with normal and low emissive ceiling in Luleå, Sweden. The numerical methods include comparative simulations by the use of dynamic heat flux models, Building Energy Simulation (BES), Computational Fluid Dynamics (CFD) and a coupled model for BES and CFD. Several parametric studies and thermal performance analyses were carried out in combination with the different numerical methods. The parallel field measurements on the test cabins include the air, surface and radiation temperatures and energy use during passive and active (heating and cooling) measurements. Both measurement and comparative simulation results indicate an improvement in the indoor thermal environment when the interior surfaces have low emittance. In the ice rink arenas, surface and radiation temperature measurements indicate a considerable reduction in the ceiling-to-ice radiation by the use of low emittance surfaces, in agreement with a ceiling-toice radiation model using schematic dynamic heat flux calculations. The measurements in the test cabins indicate that the use of low emittance surfaces can increase the vertical indoor air temperature gradients depending on the time of day and outdoor conditions. This is in agreement with the transient CFD simulations having the boundary condition assigned on the exterior surfaces. The sensitivity analyses have been performed under different outdoor conditions and surface thermal radiation properties. The spatially resolved simulations indicate an increase in the air and surface temperature gradients by the use of low emittance coatings. This can allow for lower air temperature at the occupied zone during the summer. The combined effect of interior and exterior reflective coatings in terms of energy use has been investigated by the use of building energy simulation for different climates and internal heat loads. The results indicate possible energy savings by the smart choice of optical properties on interior and exterior surfaces of the building. Overall, it is concluded that the interior reflective coatings can contribute to building energy savings and improvement of the indoor thermal environment. This can be numerically investigated by the choice of appropriate models with respect to the level of detail and computational load. This thesis includes comparative simulations at different levels of detail.

Root niche separation can explain avoidance of seasonal drought stress and vulnerability of overstory trees to extended drought in a mature Amazonian forest

Large areas of Amazonian evergreen forest experience seasonal droughts extending for three or more months, yet show maximum rates of photosynthesis and evapotranspiration during dry intervals. This apparent resilience is belied by disproportionate mortality of the large trees in manipulations that reduce wet season rainfall, occurring after 2-3 years of treatment. The goal of this study is to characterize the mechanisms that produce these contrasting ecosystem responses. A mechanistic model is developed based on the ecohydrological framework of TIN (Triangulated Irregular Network)-based Real Time Integrated Basin Simulator + Vegetation Generator for Interactive Evolution (tRIBS+VEGGIE). The model is used to test the roles of deep roots and soil capillary flux to provide water to the forest during the dry season. Also examined is the importance of "root niche separation," in which roots of overstory trees extend to depth, where during the dry season they use water stored from wet season precipitation, while roots of understory trees are concentrated in shallow layers that access dry season precipitation directly. Observational data from the Tapajo's National Forest, Brazil, were used as meteorological forcing and provided comprehensive observational constraints on the model. Results strongly suggest that deep roots with root niche separation adaptations explain both the observed resilience during seasonal drought and the vulnerability of canopy-dominant trees to extended deficits of wet season rainfall. These mechanisms appear to provide an adaptive strategy that enhances productivity of the largest trees in the face of their disproportionate heat loads and water demand in the dry season. A sensitivity analysis exploring how wet season rainfall affects the stability of the rainforest system is presented. Citation: Ivanov, V. Y., L. R. Hutyra, S. C. Wofsy, J. W. Munger, S. R. Saleska, R. C. de Oliveira Jr., and P. B. de Camargo (2012), Root niche separation can explain avoidance of seasonal drought stress and vulnerability of overstory trees to extended drought in a mature Amazonian forest, Water Resour. Res., 48, W12507, doi:10.1029/2012WR011972.

Laminar-turbulent transition induced by a discrete roughness element in a supersonic boundary layer

The linear instability and breakdown to turbulence induced by an isolated roughness element in a boundary layer at Mach 2:5, over an isothermal flat plate with laminar adiabatic wall temperature, have been analysed by means of direct numerical simulations, aided by spatial BiGlobal and three-dimensional parabolized (PSE-3D) stability analyses. It is important to understand transition in this flow regime since the process can be slower than in incompressible flow and is crucial to prediction of local heat loads on next-generation flight vehicles. The results show that the roughness element, with a height of the order of the boundary layer displacement thickness, generates a highly unstable wake, which is composed of a low-velocity streak surrounded by a three-dimensional high-shear layer and is able to sustain the rapid growth of a number of instability modes. The most unstable of these modes are associated with varicose or sinuous deformations of the low-velocity streak; they are a consequence of the instability developing in the three-dimensional shear layer as a whole (the varicose mode) or in the lateral shear layers (the sinuous mode). The most unstable wake mode is of the varicose type and grows on average 17% faster tan the most unstable sinuous mode and 30 times faster than the most unstable boundary layer mode occurring in the absence of a roughness element. Due to the high growthrates registered in the presence of the roughness element, an amplification factor of N D 9 is reached within 50 roughness heights from the roughness trailing edge. The independently performed Navier–Stokes, spatial BiGlobal and PSE-3D stability results are in excellent agreement with each other, validating the use of simplified theories for roughness-induced transition involving wake instabilities. Following the linear stages of the laminar–turbulent transition process, the roll-up of the three-dimensional shear layer leads to the formation of a wedge of turbulence, which spreads laterally at a rate similar to that observed in the case of compressible turbulent spots for the same Mach number.