Thermal expansion, fluid flow, and thermal shock of cement and a cement/steel interface at elevated pressure and temperature .
Abstract not provided.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Abstract not provided.
Explore the source record for details and available documents.
A well cement composite and a method for making a well cement composite includes a mixture of calcium aluminate cement (CAC) and fly ash cenospheres (CS) in a weight ratio of from 30:70 to 80:20 CAC to CS; sodium metasilicate (SMS) in an amount of from 1 to 10% of the total weight of the mixture of CAC and CS; polymethylhydrosiloxane (PMHS) in an amount of from 0.5 to 6.0% of the total weight of the mixture of CAC and CS; and water in a weight ratio of from 0.5:1.0 to 1.2:1.0 of water to CAC and CS.
Purpose AlSi10Mg alloy is commonly used in laser powder bed fusion due to its printability, relatively high thermal conductivity, low density and good mechanical properties. However, the thermal conductivity of as-built materials as a function of processing (energy density, laser power, laser scanning speed, support structure) and build orientation, are not well explored in the literature. This study aims to elucidate the relationship between processing, microstructure, and thermal conductivity. Design/methodology/approach The thermal conductivity of laser powder bed fusion (L-PBF) AlSi10Mg samples are investigated by the flash diffusivity and frequency domain thermoreflectance (FDTR) techniques. Thermal conductivities are linked to the microstructure of L-PBF AlSi10Mg, which changes with processing conditions. The through-plane exceeded the in-plane thermal conductivity for all energy densities. A co-located thermal conductivity map by frequency domain thermoreflectance (FDTR) and crystallographic grain orientation map by electron backscattered diffraction (EBSD) was used to investigate the effect of microstructure on thermal conductivity. Findings The highest through-plane thermal conductivity (136 ± 2 W/m-K) was achieved at 59 J/mm 3 and exceeded the values reported previously. The in-plane thermal conductivity peaked at 117 ± 2 W/m-K at 50 J/mm 3 . The trend of thermal conductivity reducing with energy density at similar porosity was primarily due to the reduced grain size producing more Al-Si interfaces that pose thermal resistance. At these interfaces, thermal energy must convert from electrons in the aluminum to phonons in the silicon. The co-located thermal conductivity and crystallographic grain orientation maps confirmed that larger colonies of columnar grains have higher thermal conductivity compared to smaller columnar grains. Practical implications The thermal properties of AlSi10Mg are crucial to heat transfer applications including additively manufactured heatsinks, cold plates, vapor chambers, heat pipes, enclosures and heat exchangers. Additionally, thermal-based nondestructive testing methods require these properties for applications such as defect detection and simulation of L-PBF processes. Industrial standards for L-PBF processes and components can use the data for thermal applications. Originality/value To the best of the authors’ knowledge, this paper is the first to make coupled thermal conductivity maps that were matched to microstructure for L-PBF AlSi10Mg aluminum alloy. This was achieved by a unique in-house thermal conductivity mapping setup and relating the data to local SEM EBSD maps. This provides the first conclusive proof that larger grain sizes can achieve higher thermal conductivity for this processing method and material system. This study also shows that control of the solidification can result in higher thermal conductivity. It was also the first to find that the build substrate (with or without support) has a large effect on thermal conductivity.
Studying the thermal conductivity of 2D materials is important due to the applications of 2D materials in fields such as thermal management, thermoelectricity, renewable energy, and sensors. As such, measurements of the thermal conductivity of these 2D materials become important to measure. Thermal conductivity is often difficult to measure for 2D materials due to their atomically thin nature and many experimental methods for doing so requiring contact with the sample, which can alter the thermal properties. A non-contact method for calculating the thermal conductivity of 2D materials supported on substrates in order to model the thermal conductivity of 2D materials for devices, is proposed and experimentally performed in this dissertation. The optothermal Raman technique is a useful non-contact diagnostic technique useful in determining the thermal conductivity of 2D materials. The optothermal Raman typically does not account for heat losses due to convection or radiation or substrate resistance, which are shown to be important factors to consider when developing an optothermal Raman model. Additionally, the calculation of the interfacial thermal conductance between the bottom surface of the sample and the top surface of the substrate, plays an important role in determining the final value of the thermal conductivity of a supported sample, and will yield differing results based on whether or not the conductance is calculated using an approach such as the Diffuse Mismatch Model (DMM) or calculated directly by varying the laser heating profile (usually done by changing the laser objective). This is shown to be the case for both graphene on Ni, graphene on Cu, and SnSe 2 on Cu. In addition to experimentally calculating the thermal conductivity of a 2D material with the optothermal Raman technique, the thermal conductivity of 2D materials can also be calculated using computational methods. The three-phonon method is a method which can be used to simulate phonon scattering processes and determine the thermal conductivity of semiconductors, wherein phonon scattering is the dominant mechanism which determines the thermal conductivity. The three-phonon method uses relaxation times for phonon scattering with other phonons, electrons, and other material system elements, such as isotopes or material defects, in order to create a single-mode relaxation time approximation (SMRTA), which is used to calculate the final value of the thermal conductivity. An important consideration when determining the thermal conductivity of a 2D material using this method is the device geometry, which is reflected in this work as the phonon-boundary scattering relaxation time. This inclusion is important along with the inclusion of phonon-electron scattering in accurately determining the thermal conductivity of a 2D material. In both the optothermal Raman experiments and the three-phonon method computations, strain is shown to have a demonstrable effect on the thermal conductivity of 2D materials. When a 1.1% strain was applied to the mechanical properties of SnSe, the three-phonon processes yielded a lower thermal conductivity than the no-strain case. For the optothermal Raman experiments, the strain induced in the Cu substrate and transferred to a single-layer graphene (SLG) sample yields a trend where the thermal conductivity of the SLG decreases with respect to strain applied. In the case where the interfacial thermal conductance was calculated directly, the conductance increased with respect to strain applied. This presents strain as a reliable and viable method for tuning the thermal properties of 2D materials for device applications.
We investigated the impact of disc-shaped and linear thermal enhancements on performance of closed-loop geothermal systems using steady-state and transient COMSOL numerical simulations. Thermal enhancements refer to thermally conductive material introduced in the rocks surrounding a wellbore to compensate for the relatively low rock thermal conductivity and increase heat extraction. Materials proposed include composite cements with thermally conductive additives such as metals or graphite to obtain overall cement thermal conductivities of 10 to 100x the rock thermal conductivity. Our simulation approach was modeling a subsegment (e.g., 100 m long) of a wellbore with idealized (e.g., perfect disc-shaped) thermal enhancements, and we report results as relative increase in heat extraction with respect to systems without thermal enhancements. We did not evaluate technical feasibility or cost of installing thermal enhancements. Simulations indicate that for 5-mm thick disc-shaped thermal enhancements with 5-m radius, repeated every 1 m along the wellbore and with thermal conductivity of 100x the rock thermal conductivity, the thermal output increases roughly 20% with respect to a closed-loop system without thermal enhancements. For fishbone structure thermal enhancements, we estimate a roughly 7% increase in thermal output for 5-m long radially outward pipes with half the radius of the main wellbore, repeated every 2.5 m, with thermal conductivity of 100x the rock thermal conductivity.
The spreading rates of convective thermals are linked to their net entrainment, and previous literature has suggested differences in spreading rates between moist and dry thermals. Here, growth rates of idealized numerically simulated axisymmetric dry and moist convective thermals are directly compared. In an environment with dry-neutral stratification, the increase of thermal radius with height dR / dz is a larger by a factor of 1.7 for dry thermals relative to moist thermals. The fractional change in thermal volume ε is also greater for dry thermals within a distance of ~4 radii from the initial thermal height. Values of dR / dz are nearly constant with height for both moist and dry thermals consistent with classical theory based on dimensional analysis. The simulations are also consistent with theory relating impulse, circulation, and spreading rate for dry thermals proposed in previous papers and extended here to moist thermals, predicting they will spread less than dry thermals. Tests adding heating to dry thermals, either spread uniformly across the thermal volume or concentrated in the inner core, indicate that dR / dz and ε are smaller for moist thermals because latent heating is confined mostly to their cores. Additional axisymmetric moist simulations using modified lapse rates and large-eddy simulations support this analysis. Overall, these results indicate that slow spreading rates are a fundamental feature of moist thermals caused by latent heating that alters the spatial distribution of buoyancy within them relative to dry thermals.
In a recent report in Matter, Bergström and co-workers describe a novel thermal transport behavior that subverts this assumption of porous materials. In their work, the authors describe a nanocellulose-based foam that features ultrahigh porosity (> 99.6%) and aligned µm-scale pores (10–100 µm) whose radial thermal conductivity (i.e., perpendicular to the pore alignment) is close to that of free air in the dry state. Furthermore, the radial thermal conductivity of the cellulose-based foam can be reduced to ~14 mW/(m·K) when the relative humidity is ~35%. Super-thermal insulators, which demonstrate a thermal conductivity below that of stationary air (~ 25 mW/(m·K), 20 ºC, 1.0 atm), are needed to minimize heat loss in various applications (e.g., buildings, thermal energy storage tanks, cold chain packaging, etc.) to mitigate the energy crisis and reduce carbon emissions. Introducing pores into a material is a facile and effective way to achieve low thermal conductivity as pores can suppress thermal transport through solids by reducing the cross-sectional area and increasing the tortuosity of the heat transfer pathway. In porous structures there are two kinds of pores: open and closed. While increased porosity can reduce heat conduction through solids with an open porous structure, the improved gas conduction creates a competing effect that simultaneously elevates the heat transfer. Therefore, the thermal conductivity of a material with µm-scale open pores is usually larger than that of stationary air. Reducing the pore size to less than the mean free path of air (~ 70 nm, 20 ºC, 1.0 atm) can effectively reduce gas conduction, enabling the material to achieve a thermal conductivity below that of stationary air. However, high cost of the nanosized raw materials and time-consuming fabrication processes limit the large-scale applications of nanoporous thermal insulators. Meanwhile, closed pores can block heat transport only through the continuous gas phase. As a result, the thermal conductivity of a closed porous structure can theoretically be much smaller than that of stationary air if the solid conduction can also be suppressed by (1) lowering the solid content of the material, (2) reducing the thermal conductivity of the building blocks of the material, and/or (3) increasing interfacial thermal resistance between neighboring building blocks. Most processes used to generate porous structures (e.g., supercritical drying, freeze drying) involve a solvent that escapes the material. Therefore, it is difficult to create pores and isolate them simultaneously, preventing the fabrication of super-thermal insulators with closed pores.
Here this study examines the entrainment, detrainment, and dilution of dry and moist (cloud) atmospheric thermals in large-eddy simulations. In a neutrally stable environment (with respect to dry dynamics), moist thermals have an increase in radius R with thermal height z t (α ≡ dR/dz t ) about 4 times smaller compared to dry thermals when density stratification is considered and ~2.4 times smaller without density stratification (i.e., applying the Boussinesq approximation). An analytic expression relating α to several dimensionless parameters is derived from the thermal impulse–circulation relation to clarify the factors impacting α. This expression shows that the difference in buoyancy structure between moist and dry thermals, with buoyancy concentrated in the central cores of moist thermals owing to latent heating, explains their smaller spreading rates. Individual contributions of entrainment and detrainment are analyzed using a direct parcel-based approach in the simulations. Moist thermals have similar fractional detrainment but much smaller fractional entrainment rates compared to dry thermals, consistent with the differences in α. Despite having smaller α, moist thermals are similarly dilute (quantified by a passive tracer) as dry thermals because of their greater mixing efficiency with the environment. Thus, moist thermals are substantially dilute but expand much less in size/volume as they rise compared to dry thermals. The α values for moist thermals in (dry) neutral and statically stable environments are similar, but fractional entrainment and especially detrainment rates are greater in the stable environment. Large detrainment rates are associated with a breakdown of the broader thermal vortex ring structure, especially with low environmental relative humidity, attributed in part to evaporation and buoyancy reversal.
Carbon nanotubes (CNTs) are quasi-one dimensional nanostructures that display both high thermal conductivity for potential thermal management applications and intriguing low-dimensional phonon transport phenomena. In comparison to the advances made in the theoretical calculation of the lattice thermal conductivity of CNTs, thermal transport measurements of CNTs have been limited by either the poor temperature sensitivity of Raman thermometry technique or the presence of contact thermal resistance errors in sensitive two-probe resistance thermometry measurements. Here we report advances in a multi-probe measurement of the intrinsic thermal conductivity of individual multi-walled CNT samples that are transferred from the growth substrate onto the measurement device. The sample-thermometer thermal interface resistance is directly measured by this multi-probe method and used to model the temperature distribution along the contacted sample segment. The detailed temperature profile helps to eliminate the contact thermal resistance error in the obtained thermal conductivity of the suspended sample segment. A differential electro-thermal bridge measurement method is established to enhance the signal-to-noise ratio and reduce the measurement uncertainty by over 40%. The obtained thermal resistances of multiple suspended segments of the same MWCNT samples increase nearly linearly with increasing length, revealing diffusive phonon transport as a result of phonon-defect scattering in these MWCNT samples. The measured thermal conductivity increases with temperature and reaches up to 390 ± 20 W m –1 K –1 at room temperature for a 9-walled MWCNT. In conclusion, theoretical analysis of the measurement results suggests submicron phonon mean free paths due to extrinsic phonon scattering by extended defects such as grain boundaries. The obtained thermal conductivity is decreased by a factor of 3 upon electron beam damage and surface contamination of the CNT sample.
Vacuum-insulated glazing (VIG) with a low-emittance coating has a great market potential as an effective transparent insulator. The thermal insulating performance of VIG is determined by its design, including material selection and configuration of different components. Thermal conductance of the vacuum gap as a transport bottleneck is one of the primary factors controlling the thermal transport across VIG. In particular, because support pillars provide the main thermal transport channels across the vacuum gap, increasing the pillar thermal resistance is a key strategy for creating effective thermal insulation while maintaining the vacuum space. Here, the effects of various pillar design parameters, such as thermal conductivity, geometry, and arrangement, on the VIG thermal performance were comprehensively investigated via the finite element method. In addition, analytical models for thermal transport were examined and thermal conductance across the VIG unit was experimentally measured for validation. The pillar design parameters, especially the height, shape, spacing, and arrangement of the pillars, showed significant effects on the thermal performance of VIG. This research also shows that the smaller contact area for horizontal pillars can effectively decrease the heat loss by more than 30%. Because current VIG analytical equations of thermal performance are only applicable to cylindrical pillars, an analytical equation that can better describe the thermal performance of rectangular parallelepiped pillars is presented, along with a discussion about the mechanism of thermal transfer for different pillar shapes. Thermo-mechanical analyses based on 3D FEM simulations can provide valuable insights into the effect of various design parameters on the overall performance the VIG, allowing for the development of an optimal VIG design.
This report investigates numerical modeling methods for thermal ratcheting analysis of packed-bed thermal energy storage (TES) tank and discusses the validation results via comparison with experimental data. The experimental data obtained from various design characteristics of packed-bed thermocline tanks, including the Thermal Energy Distribution System (TEDS) TES tank at Idaho National Laboratory (INL), were used to validate thermal and mechanical models developed in this study to evaluate the thermal ratcheting potential. The thermal model was shown to predict the transient thermal propagation through the packed-bed thermocline tanks generally well. However, a larger discrepancy was observed during the comparison with the data from TEDS, presumably due to the uncertainty of boundary conditions given from the experiment. Based on the comparative study between the thermal model predictions and experimental data of various packed-bed thermocline tanks, potential improvements were suggested for the future TEDS experiments for more precise validation study. For mechanical (thermally induced stress) analysis, two different modeling approaches were tested to evaluate hoop stress applied to the packed-bed TES tank wall, which is a major cause of thermal ratcheting process: (i) infinite rigidity model and (ii) Drucker-Prager (DP) model. The ‘model (i)’ is a conservative method with infinite rigidity assumption of granular filler inside a TES tank, whereas the ‘model (ii)’ is a method that takes into account more realistic processes such as thermal expansion of filler and tank wall as well as inter-particle interactions during the cyclic operation of a packed-bed TES tank. The validity of each modeling method was examined by comparing the numerical simulation with the experimental data obtained from the packed-bed TES tank for Solar One pilot plant. Then, the effects of various model parameters were discussed to evaluate the thermal ratcheting potential of the TEDS TES tank. The preliminary thermal ratcheting analysis implies that the TEDS TES tank will hold its structural integrity during the normal operation cycles.
Abstract Liquid metals are being investigated as coolants in many advanced reactor designs because of their high thermal conductivity and effectiveness at high temperatures. However, they often pose challenges to reactor operation and safety because of the complex thermal mixing and stratification in the plenum of pool-type reactor designs. The advanced system analysis code System Analysis Module (SAM) currently under development at Argonne National Laboratory aims to develop and implement thermal mixing models to accurately capture these complex thermal fluid behaviors. In this study, the SAM thermal mixing model was compared against experimental data from the Gallium Thermal-Hydraulic Experiment facility, a scaled liquid metal test facility that uses gallium as a surrogate fluid to investigate the stratification and thermal mixing of low-Prandtl-number fluids in the upper plenum of a liquid metal-cooled reactor. Two cold shock transient cases were used: one with stable stratified flow (Ri = 32) and one with stronger thermal mixing (Ri = 0.5). The resultant temperatures were then compared with the experimental temperatures over the entire plenum to assess the ability of the mixing models to capture the thermal behavior and to better correspond mixing parameters to various flow scenarios. Generally, the zero-dimensional mixing model was more capable of capturing the bulk temperature of the component modeled assuming that an accurate mass flow rate was provided, but it was inherently unable to capture thermal gradients in space. The one-dimensional mixing model was capable of capturing that the thermal gradients provided accurate selection of the mixing coefficients. Further, the temperature at the outlet junction was compared over time for each of the mixing models with the recorded experimental temperature. The implemented mixing models demonstrated the ability to effectively capture the overall thermal behavior for stronger mixing scenarios but struggled with more stably stratified flows. It was found that a system analysis code's covering of the entire range of different operating conditions still remains a challenging task, and it is suggested that further model and closure improvements are necessary to accurately capture complex thermal mixing and stratification phenomena.
Abstract Understanding factors influencing carbon effluxes from soils to the atmosphere is important in a world experiencing climatic change. Two important uncertainties related to soil organic carbon (SOC) stock responses to a changing climate are (a) whether soil microbial communities acclimate or adapt to changes in soil temperature and (b) how to represent this process in SOC models. To further explore these issues, we included thermal adaptation of enzyme‐mediated processes in a mechanistic SOC model (ReSOM) using the macromolecular rate theory. Thermal adaptation is defined here to encompass all potential responses of soil microbes and microbial communities following a change in temperature. To assess the effects of thermal adaptation of enzyme‐mediated processes on simulated SOC losses, ReSOM was applied to data collected from a 13‐year soil warming experiment. Results show that a model omitting thermal adaptation of enzyme‐mediated processes substantially overestimates observed CO 2 effluxes during the initial years of soil warming. The bias against observed CO 2 effluxes was lower for models including thermal adaptation of enzyme‐mediated processes. In addition, for a simulated linear 3°C soil warming over 100 years, models including thermal adaptation of enzyme‐mediated processes simulated SOC losses of a factor of three smaller than models omitting this process. As thermal adaptation of microbial community characteristics is generally not included in models simulating feedback between the soil, biosphere and atmosphere, we encourage future studies to assess the potential impact that microbial adaptation has on soil carbon – climate feedback representations in models. Plain Language Summary A major uncertainty in projecting how much soil organic carbon (SOC) will be converted to CO 2 as a consequence of climate change is related to how soil microbes may adapt to increasing soil temperatures. While this “microbial thermal adaptation” has been shown to occur in short‐term lab incubation experiments, its effect on SOC cycling on a decadal timescale is not clear. To address this knowledge gap, a mechanistic SOC model was used to simulate data collected from a 13‐year soil warming experiment, to assess how microbial thermal adaptation affects predicted SOC losses upon soil warming. The model results show that incorporating microbial thermal adaptation into the model led to reduced CO 2 effluxes from the soil to the atmosphere compared to the common approach of omitting this mechanism. Our results imply that projected SOC losses for the decades to come may be reduced when this mechanism is incorporated in land models. We therefore advocate for more research on the mechanisms controlling microbial thermal adaptation, and how to implement this mechanism in SOC models. Key Points A crucial aspect of soil organic carbon (SOC) models is the representation of soil microbes Predicted soil CO 2 fluxes upon soil warming are reduced when accounting for microbial thermal adaptation On a centennial time scale, this thermal adaptation results in up to a factor of three lower predicted SOC loss
Abstract This study examines dynamic pressure drag on rising dry buoyant thermals. A theoretical expression for drag coefficient C d as a function of several other nondimensional parameters governing thermal dynamics is derived based on combining the thermal momentum budget with the similarity theory of Scorer. Using values for these nondimensional parameters from previous studies, the theory suggests drag on thermals is small relative to that on solid spheres in laminar or turbulent flow. Two sets of numerical simulations of thermals in an unstratified, neutrally stable environment using an LES configuration of the Cloud Model 1 (CM1) are analyzed. One set has a relatively low effective Reynolds number R e and the other has an order-of-magnitude-higher R e ; these produce laminar and turbulent resolved-scale flows, respectively. Consistent with the theoretical C d , the magnitude of drag is small in all simulations. However, whereas the laminar thermals have C d ≈ 0.01, the turbulent thermals have weakly negative drag ( C d ≈ −0.1). This difference is explained by the laminar thermals having near vertical symmetry but the turbulent thermals exhibiting considerable vertical asymmetry of their azimuthally averaged flows. In the laminar thermals, buoyancy rapidly becomes concentrated around the main centers of rotation located along the horizontal central axis, leading to expansion of thermals via baroclinic vorticity generation but doing little to break vertical symmetry of the flow. Vertical asymmetry of the azimuthally averaged flow of turbulent thermals is attributed mainly to small-scale resolved eddies that are concentrated in the upper part of the thermals.
Poly(vinyl chloride) (PVC) is widely used in various fields and requires the use of thermal stabilizers to enhance its thermal stability during processing because of its poor thermal stability. Layered double hydroxides (LDHs) are widely considered to be one kind of highly efficient and environmentally friendly PVC thermal stabilizer. To investigate the thermal stabilizing process of layered double hydroxides (LDHs) in PVC resin, PVC and MgAl-LDHs powders with different interlayer anions (CO 3 2- , Cl - , and NO 3 - ) were physically mixed and aged at 180 °C. The structure of LDHs at different aging times was studied using XRD, SEM, and FT-IR. The results show that the thermal stabilizing process of LDHs on PVC mainly has three stages. In the first stage, the layers of LDHs undergo a reaction with HCl, which is released during the thermal decomposition of PVC. Subsequently, the ion exchange process occurs between Cl - and interlayer CO 3 2- , resulting in the formation of MgAl-Cl-LDHs. Finally, the layers of MgAl-Cl-LDHs react with HCl slowly. During the thermal stabilizing process of MgAl-Cl-LDHs, the peak intensity of XRD reduces slightly, and no new XRD peak emerges. It indicates that only the first step happens for MgAl-Cl-LDHs. The TG-DTA analysis of LDHs indicates that the interaction of LDHs with different interlayer anions has the following order: NO 3 - < CO 3 2- < Cl - , according to the early coloring in the thermal aging test of PVC composites. The results of the thermal aging tests suggest that LDHs with a weak interaction between interlayer anions and layers can enhance the early stability of PVC significantly. Furthermore, the thermal aging test demonstrates that LDHs with high HCl absorption capacities exhibit superior long-term stabilizing effects on PVC resin. This finding provides a valuable hint for designing an LDHs/PVC resin with a novel structure and excellent thermal stability.
Global warming due to climate change is a threat to humankind. Nuclear energy is one of the promising solutions to reduce fossil fuel usage. Nuclear energy can handle the base load, compensating for the volatility of renewable energy. If nuclear energy could achieve load following capability, the combination with renewable energy would be more suitable. Thermal energy storage (TES) is one of the options for enabling load following of nuclear reactors. The TES makes it possible to store surplus nuclear thermal energy and release it later as needed. In Idaho National Laboratory (INL), a new concept of latent heat TES integrated with high-temperature heat pipe has been proposed and is under development, which is called Heat pipe-Integrated Thermal Battery (HITB). HITB exchanges thermal energy between the reactor system and TES via heat pipe. The heat transferred to TES medium, made of phase change material (PCM), stores energy as sensible heat and/or latent heat. As PCM typically has poor thermal conductivity, however, various heat transfer enhancement techniques are required to achieve a rapid charging cycle. There are many techniques to enhance the heat transfer ability of TES medium such as disk, fin, and metal foam. Among them, metal foam is an appropriate option to enhance the heat transfer because it maximizes the heat transfer area through metal wicks. Metal foam is a lightweight metal structure that has a high porosity of over 0.9. The typical materials for metal foam are Aluminum, Copper, Nickel, and Silicon Carbide (SiC). Metal foam not only enhances heat transfer via conduction but also increases contact surface area. In the HITB design , the metal foam is being considered as one of the options to enhance the heat transfer of TES medium (PCM) [1]. To predict the enhanced thermal performance of TES, one should properly estimate the effective thermal conductivity of metal foam combined with PCM material or calculate heat transfer in distributed model. There are many experimental works that provides effective thermal conductivity of metal foam with various PCM [2,3]. Also, many theoretical models were developed based on the unit cell model of metal foam [4,5]. With a distributed model, on the other hand, detail heat transfer characteristics between metal foam and PCM material can be analyzed considering the geometry or buoyancy effect. However, due to the complex geometry of metal foam pores, the computational cost for three-dimensional modeling highly increases. Therefore, if metal foam structure can be modeled in simple and repetitive design, the computational cost would decrease Among the various metal foam models [2], lattice model is one of the simple and extendable design. The porosity and pores per inch (PPI) can be characterized by the size and spatial distance of lattice structure. If the three-dimensional metal foam model consists of lattice structure could properly estimate the heat transfer, which is characterized by effective thermal conductivity, it would be a good option to assess the thermal performance of metal foam with PCM. In this study, a three-dimensional numerical model was developed to simulate conductive heat transfer between metal foam and PCM. The three-dimensional lattice structure of square pillars was selected as a basic structure of the metal foam. The calculation result was characterized by the effective thermal conductivity of the whole domain. A sensitivity study was conducted for mesh size, domain size, and PPI to check whether the calculation result gives a converged result or not. Lastly, the effective thermal conductivity from the lattice model was compared with existing experimental data to validate the model result
Thermal energy storage systems in buildings can store cooling/heating energy during non-peak load hours or when renewable energy sources are available for later use when demanded. Using the building envelope as a thermal battery can help to shave the peak load, reduce the burden in electric grids, and enhance the occupant's thermal comfort. Sensible energy storage on wall systems such as thermally activated building systems can provide an active thermal storage strategy. However, most of the stored energy is used through passive means directed by the thermal lag, which can impede the on-demand release of the stored energy. Integrating active insulation systems with building thermal storage systems can increase the flexibility of charging and discharging time and duration. In this study, a wall system equipped with an active insulation system and thermally activated storage system was designed, and its performance on active cooling energy contribution was studied. The performances of energy storage (charging), release (discharging) of the thermal energy storage energy, and the active insulation system were studied separately and together as an integrated system. Results showed that the thermal properties of the thermal energy storage core material and the pipe spacing of both embedded pipes in the thermal energy storage and hydronic pipes used in the active insulation system affected the wall performance the most out of all the tested wall system parameters. During discharging, the heat flux into the wall was as high as 81.92 W/m 2 . This active insulation system's dynamic R-value varied between less than 1ft 2 •°F•h/BTU (0.18 m 2 •K/W) to 98% of the equally thick foam insulation's R-value.