Kenneth W Childs

A model of heat transfer in sapwood and implications for sap flux density measurements using thermal dissipation probes

A variety of thermal approaches are used to estimate sap flux density in stems of woody plants. Models have proved valuable tools for interpreting the behavior of heat pulse, heat balance and heat field deformation techniques, but have seldom been used to describe heat transfer dynamics for the heat dissipation method. Therefore, to better understand the behavior of heat dissipation probes, a model was developed that takes into account the thermal properties of wood, the physical dimensions and thermal characteristics of the probes, and the conductive and convective heat transfer that occurs due to water flow in the sapwood. Probes were simulated as aluminum tubes 20 mm in length and 2 mm in diameter, whereas sapwood, heartwood and bark each had a density and water fraction that determined their thermal properties. Base simulations assumed a constant sap flux density with sapwood depth and no wounding or physical disruption of xylem beyond the 2 mm diameter hole drilled for probe installation. Simulations across a range of sap flux densities showed that the dimensionless quantity k [defined as (ΔT(m) -ΔT)/ΔT, where ΔT(m) is the temperature differential (ΔT) between the heated and unheated probe under zero-flow conditions] was dependent on the thermal conductivity of the sapwood. The relationship between sap flux density and k was also sensitive to radial gradients in sap flux density and to xylem disruption near the probe. Monte Carlo analysis in which 1000 simulations were conducted while simultaneously varying thermal conductivity and wound diameter revealed that sap flux density and k showed considerable departure from the original calibration equation used with this technique. The departure was greatest for variation in sap flux density typical of ring-porous species. Depending on the specific combination of thermal conductivity and wound diameter, use of the original calibration equation resulted in an 81% under- to 48% overestimation of sap flux density at modest flux rates. Future studies should verify these simulations and assess their utility in estimating sap flux density for this widely used technique.

Thermal behavior of mixtures of perlite and phase change material in a simulated climate

Carefully controlled and well documented experiments have been done for two candidate configurations to control the heat load on a conditioned space. The 2:1 PCM/perlite mixture and the 6:1 PCM/perlite mixture, both on a weight basis, accomplished thermal control. The 2:1 system seemed to have enough PCM to be effective and involve a much larger fraction of its PCM in diurnal freezing and melting than the 6:1 system. It is a good starting point for engineering design of an optimum thermal control system. The results from the 2:1 system were reproduced with the computer program HEATING to prove that we know the relevant mechanisms and thermophysical properties of the PCM used in the system. Even without a model for the supersaturation and hysteresis that this material exhibited, HEATING reproduced the heat fluxes to the conditioned space in the experiments accurately enough to mirror the good thermal control performance of the system. The modified sensible heat capacity that was used in HEATING is a handy way to account for phase change effects and could be used in a subroutine to compute hourly phase change effects for whole building models like DOE-2. The experiments were done with PCM/perlite mixtures sealed in small methylmethacrylate boxes and covered top and bottom by XPS. The boxes allowed precise placement of the instrumentation used to follow the phase change effects. The XPS gave high R-value per unit thickness. A more practical prototype configuration such as PCM/perlite hermetically sealed in plastic pouches between layers of batts or blown-in insulation should be tested over a larger cross section. A good candidate is the whole attic cavity of the manufactured home test section used in the present work. Use of a PCM that does not exhibit supersaturation and hysteresis would make interpretation of the results easier. If the results of the larger scale test areas are as encouraging as the test cell results, a whole house model with a phase change algorithm should be constructed to optimize the configuration for the climate in which it will perform.

In-Depth Analysis of Simulation Engine Codes for Comparison with DOE s Roof Savings Calculator and Measured Data

...................................................................................................................................................... vii