SEARCH · Search NASA

Results for “Biot model”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Robust block preconditioners for poroelasticity

In this paper we develop and analyze several preconditioners for the linear systems arising from discretized poroelasticity problems. The preconditioners include one block preconditioner for the two-field Biot model and several preconditioners for the classical three-field Biot model. We manage to analyze these different preconditioners under a same theoretical framework and show that all of them are uniformly optimal with respect to material and discretization parameters. Numerical tests demonstrate the robustness of these preconditioners.

42 ENGINEERING↗

DOE OSTI2021

Stepwise Dynamic Calibration of a Hydromechanical Simulation Using Time-Lapse Vertical Seismic Profile

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse Vertical Seismic Profiling (VSP) integration. The selected study site is the 13-10A injector well within the ongoing CO2-EOR operation of the Farnsworth Field Unit. The Time-lapse VSP dataset carries the combined effects of fluid substitution and mean effective stress changes, thereby providing a dataset amenable for the calibration of production and injection-induced stress changes. The concept is similar to calibrating a reservoir simulation model in that the process honor real field data to set up an inverse problem. The solution optimizes the independent and impactful geomechanical parameters that replicate the observed time-lapse seismic velocity changes. This stress calibration is enabled by 4D geomechanical modeling and the VSP Integration workflow. This calibration benefits from extensive geological, geophysical and geomechanical characterization through 3D seismic data, geophysical well logs, and core assessed as part of the 1D MEM conducted on the 13-10A subject well. These data are used to develop a site-specific rock physics model. The Biot Gassmann workflow combines rock physics and reservoir simulation outputs to determine the fluid substitution contribution to seismic velocity change. Additionally, modeled seismic velocity attributed to mean effective stress are determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from the ultrasonic seismic velocity measurements on the extracted Morrow B core. A penalty function is then formed between the modeled seismic velocities and the observed time-lapse VSP dataset. Four independent and impactful geomechanical parameters have been determined. These are the bulk modulus and shear modulus for zero porosity and the shear and compressional seismic velocity to mean effective stress derivatives. The dataset of numerous coupled hydromechanical- geomechanical simulation realizations is built by combining variations of the four stated geomechanical parameters. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer are used to converge on the optimal geomechanical parameters. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling.

02 PETROLEUM↗

DOE OSTI2022

Experimental and Computational Investigation of Integrated Internal and Film Cooling Designs Incorporating a Thermal Barrier Coating

Abstract Few studies in the open literature have studied the effect of thermal barrier coatings (TBC) when used in combination with shaped hole film cooling and enhanced internal cooling techniques. The current study presents Reynolds-averaged Navier–Stokes (RANS) conjugate heat transfer simulations that identify trends in cooling design performance as well as experimental measurements of overall effectiveness using a flat-plate matched-Biot number model with a simulated TBC layer of 0.42D thickness, where D is the film cooling hole diameter. Coolant is fed to the film cooling holes in a co-flow configuration, and the results of both the smooth and rib-turbulated channels are compared. At a constant coolant flow rate, enhanced internal cooling was found to provide a 44% increase in spatially-averaged overall effectiveness, $\bar{\bar{\phi}}$, without a TBC. The results show that the addition of a TBC can raise $\bar{\bar{\phi}}$ on a film-cooled component surface by 47%. The optimum velocity ratio was found to decrease with the addition of enhanced cooling techniques and a TBC as the film provided minimal benefit at the expense of reduced internal cooling. While the computational results closely identified trends in overall system performance without a TBC, the model over-predicted effectiveness on the metal–TBC interface. The results of this study will inform turbine component design as material science advances increase the reliability of the TBC.

Engineering↗

DOE OSTI2024

CO 2 rock physics modeling for reliable monitoring of geologic carbon storage

Monitoring, verification, and accounting (MVA) are crucial to ensure safe and long-term geologic carbon storage. Seismic monitoring is a key MVA technique that utilizes seismic data to infer elastic properties of CO 2 -saturated rocks. Reliable accounting of CO 2 in subsurface storage reservoirs and potential leakage zones requires an accurate rock physics model. However, the widely used CO 2 rock physics model based on the conventional Biot-Gassmann equation can substantially underestimate the influence of CO 2 saturation on seismic waves, leading to inaccurate accounting. We develop an accurate CO 2 rock physics model by accounting for both effects of the stress dependence of seismic velocities in porous rocks and CO 2 weakening on the rock framework. We validate our CO 2 rock physics model using the Kimberlina-1.2 model (a previously proposed geologic carbon storage site in California) and create time-lapse elastic property models with our new rock physics method. We compare the results with those obtained using the conventional Biot-Gassmann equation. Our innovative approach produces larger changes in elastic properties than the Biot-Gassmann results. Using our CO 2 rock physics model can replicate shear-wave speed reductions observed in the laboratory. Our rock physics model enhances the accuracy of time-lapse elastic-wave modeling and enables reliable CO 2 accounting using seismic monitoring.

58 GEOSCIENCES↗

DOE OSTI2024

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗

DOE OSTI2025

Thermal Modeling and Limitations for Power Electronics Embedded in Medium-Voltage Cables

As next-generation energy technologies gain traction and power demand increases, the existing electrical infrastructure faces significant stress, prompting innovative solutions to enhance the grid's capacity and lifespan. This work explores the possibility of embedding medium-voltage (MV) power electronics directly inline with the cable, and the resulting thermal challenges. Since the majority of power distribution cables installed in the U.S. are passively cooled, the work focuses primarily on passive cooling, with an emphasis on the limitations of axial heat spreading within the cable. To date, literature on axial spreading of high incident heat loads on cables and cable environments is limited, typically reporting cases with <10 W of incident heat load. This work will explore the considerations, limits, and tradeoffs of cable-embedded heat loads significantly larger than the cable losses. Both external and internal effects are modeled analytically in nondimensional terms via a Biot number analysis, allowing fundamental limits and tradeoffs to be derived. The work culminates in the design and experimental validation of a cable-embedded thermal system capable of passively dissipating 300 W of heat from a coaxial SiC mosfet switch module over a length of 20 cm, thus validating the possibility of MV cable-embedded power electronics from a thermal standpoint.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

A Multiscale Approach to Simulate Non‐Isothermal Multiphase Flow in Deformable Porous Materials

Coupled thermal, hydraulic, and mechanical processes in porous materials play important roles in several energy and environmental technologies. The Darcy-Brinkman-Biot (DBB) framework has proven effective in modeling multiphase fluid flow in deformable porous solids across both pore and Darcy scales, including in systems where fractures coexist with a porous matrix. In this study, we extend the DBB framework, originally designed for isothermal conditions, to address non-isothermal problems by incorporating an energy conservation equation. The resulting solver, hybridBiotThermalInterFoam, enables simulations of coupled multiphase fluid flow, heat transfer, and solid deformation in hybrid-scale systems containing both solid-free regions and ductile porous domains. The new solver is validated through comparisons with analytical solutions and, also, against established heat transfer solvers chtMultiRegionFoam and compressibleInterFoam. Further, a series of 2D and 3D case studies, including two-phase heat transfer in solid-free, static, or deformable porous media, highlights the solver's capacity to simulate complex flow dynamics and heat transport in systems involving high mobility ratios, viscous fingering, and fracture propagation. Our results establish the feasibility of incorporating thermal effects in simulations of a wide variety of energy geotechnics and environmental applications, including enhanced hydrocarbon recovery, soil remediation, and enhanced geothermal energy systems.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2026

Integration of the Biot–Gassmann Fluid Substitution Method and Machine Learning-Based Velocity–Stress Relationship for Estimating In Situ Stresses

Recent advancements have shown that in situ stresses can be reliably estimated through an integrated machine/deep learning (ML/DL)-based framework, which relies on models trained and validated using true triaxial ultrasonic velocity (TUV) experimental data that involve measurements of ultrasonic velocity in saturated rocks under varying stress configurations. However, when the goal is to interpret lower frequency measurements, it may be more appropriate to run experiments on dry rocks and then obtain Biot–Gassmann-derived equivalent saturated velocities (low-frequency approximation) and employ these quantities for training ML/DL models to predict in situ stress. Whether the dispersion effect of frequency on the velocity–stress relationship substantially impacts in situ stress prediction is an important and unresolved question. This work presents an enhancement of ML/DL-based workflow by training and implementing ML/DL models using equivalent saturated acoustic velocities (low-frequency) obtained by applying Biot–Gassmann fluid substitution on the ultrasonic velocities of dry cores. The models were trained on TUV data sets derived from three subsurface cores extracted from the geothermal well 16B(78)-32 at the Utah FORGE site. Each core was subjected to 75 unique stress configurations for velocity measurement in the dry state. The ML/DL trained on the TUV data set with equivalent saturated velocities demonstrated promising performance to predict in situ stress in subsurface geological rocks using velocity–stress relationships with R 2 of 0.86, 0.971, and 0.975 and root mean squared error (RMSE) of 2.59, 1.92, and 1.80 for validation/testing phases of vertical, minimum horizontal, and maximum horizontal stress models, respectively. Additionally, interpretation and explanation by Shapley additive explanations (SHAP) analysis further improved scientific validation and model reliability for estimating in situ stresses.

colloids↗

DOE OSTI2021

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗

DOE OSTI2022

Current distribution monitoring enables quench and damage detection in superconducting fusion magnets

Abstract Fusion magnets made from high temperature superconducting ReBCO CORC ® cables are typically protected with quench detection systems that use voltage or temperature measurements to trigger current extraction processes. Although small coils with low inductances have been demonstrated, magnet protection remains a challenge and magnets are typically operated with little knowledge of the intrinsic performance parameters. We propose a protection framework based on current distribution monitoring in fusion cables with limited inter-cable current sharing. By employing inverse Biot-Savart techniques to distributed Hall probe arrays around CORC ® Cable-In-Conduit-Conductor (CICC) terminations, individual cable currents are recreated and used to extract the parameters of a predictive model. These parameters are shown to be of value for detecting conductor damage and defining safe magnet operating limits. The trained model is then used to predict cable current distributions in real-time, and departures between predictions and inverse Biot-Savart recreated current distributions are used to generate quench triggers. The methodology shows promise for quality control, operational planning and real-time quench detection in bundled CORC ® cables for compact fusion reactors.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗

DOE OSTI2023

Time-lapse VSP integration and calibration of subsurface stress field utilizing machine learning approaches: A case study of the morrow B formation, FWU

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse vertical seismic profiling (VSP) integration. The selected study site is a region around the injector well located within Farnsworth field unit (FWU), where there is an ongoing CO 2 -enhanced oil recovery (EOR) operation. In our study, a site-specific rock physics model was created from extensive geological, geophysical, and geomechanical characterization through 3D seismic data, well logs, and core assessed as part of the 1D MEM conducted on the characterization well within the study area. Here, the Biot-Gassmann workflow was utilized to combine the rock physics and reservoir simulation outputs to determine the seismic velocity change due to fluid substitution. Modeled seismic velocities attributed to mean effective stress were determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from ultrasonic seismic velocity measurements. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer (PSO) was utilized to minimize a penalty function created between the modeled seismic velocities and the observed time-lapse VSP dataset. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.

58 GEOSCIENCES↗

DOE OSTI2026

Flash Heating of a Hollow Rod to Determine Thermal Properties

In a nuclear reactor, prolonged exposure to high neutron flux levels can change material properties of reactor components. Due to the danger posed to personnel by highly radioactive components such as spent fuel, noncontact methods must be used to determine these properties. The objective of the current research is to test the viability of using a flash-heating method to determine the thermal diffusivity of a silicon carbide (SiC) tube, which could be used as cladding for a nuclear fuel rod. The flash-heating experiment was conducted at Oak Ridge National Laboratory and temperature results were analyzed as a two-Dimensional cross section which assumed axial symmetry. The model was fitted iteratively to the temperature measurements using nonlinear regression which required values for thermal diffusivity, Biot number and a term to determine the magnitude of the heat absorbed by the flash. The method could be expanded in the future to estimate the thermal diffusivity of spent nuclear fuel rods in an effort to increase the certainty of new reactor designs using precise thermal properties throughout the core life. A statistical analysis of the results of this work is provided as part of the analysis.

McMasters, Robert L [Virginia Military Institute, ↗

DOE OSTI2021

Impacts of Anisotropic Porosity on Heat Transfer and Off-Gassing during Biomass Pyrolysis

The pore structure of biogenic materials imbues the ability to deliver water and nutrients through a plant from root to leaf. This anisotropic pore granularity can also play a significant role in processes such as biomass pyrolysis that are used to convert these materials into useful products like heat, fuel, and chemicals. Evolutions in modeling of biomass pyrolysis as well as imaging of pore structures allow for further insights into the concerted physics of phase change-induced off-gassing, heat transfer, and chemical reactions. In this work, we report a biomass single particle model which incorporates these physics to explore the impact of implementing anisotropic permeability and diffusivity on the conversion time and yields predicted for pyrolysis of oak and pine particles. Simulation results showed that anisotropic permeability impacts predicted conversion time more than 2 times when the Biot number is above 0.1 and pyrolysis numbers (Py 1 , Py 2 ) are less than 20. Pore structure significantly impacts predicted pyrolytic conversion time (>8 times) when the Biot number is above 1 and the pyrolysis number is below 1, i.e., the “conduction controlled” regime. Therefore, these nondimensional numbers reflect that when internal heat conduction limits pyrolysis performance, internal pyrolysis off-gassing further retards effective heat transfer rates as a closely coupled phenomenon. Overall, this study highlights physically meaningful opportunities to improve particle-scale pyrolysis modeling and experimental validation relevant to a variety of feedstock identities and preparations, guiding the future design of pyrolyzers for efficient biomass conversion.

09 BIOMASS FUELS↗

DOE OSTI2021

Biomass pyrolysis devolatilization kinetics of herbaceous and woody feedstocks

Devolatilization kinetics were determined using a modified micropyrolyzer reactor for several biomass feedstocks: switchgrass, corn stover, red oak, and pine. Here, the micropyrolyzer was directly coupled to a flame ionization detector (FID) to track the release of volatiles from the biomass. Time series data from these experiments was analyzed to determine apparent devolatilization rates. Care was taken to assure the experiments were isothermal and kinetically limited calculating a Biot number less than 0.1 and pyrolysis numbers greater than 10, which simplifies the derivation of devolatilization rates. A single, first order reaction was able to model devolatilization rates at temperatures up to 500 °C. No correlation was found between the inorganic content of the biomass and its rate of devolatilization. Apparent activation energies were in the range of 54.9–88.4 kJ mol –1 . The rate coefficient at 500 °C was calculated as 1.90–5.14 s –1 for the four feedstocks.

09 BIOMASS FUELS↗

DOE OSTI2020

Thermal conductivity of two-dimensional disordered fibrous materials defined by interfiber thermal contact conductance and intrinsic conductivity of fibers

A general framework for the theoretical analysis and numerical calculations of the effective thermal conductivity of two-dimensional homogeneous and isotropic disordered fibrous materials is developed in this work based on the model of random contacts between soft-core spherocylinders. The analysis accounts for the interfiber contact conductance and intrinsic conductivity of fibers and is performed in a wide space of governing parameters that includes the fiber aspect ratio, Biot number calculated for a single thermal contact between fibers, and material density ranging from values corresponding to the percolation threshold up to those characteristic of dense fiber networks. For dense networks, exact theoretical equations for the thermal conductivity of materials composed of spherocylinders with an arbitrary aspect ratio and Biot number are derived. The effect of the intrinsic conductivity of fibers on the thermal transport in fibrous materials is found to depend on the density of contacts and can be significant in sufficiently dense fiber networks even if the Biot number for a single thermal contact is small. Semiempirical corrections to the theoretical equations are derived for small and moderate fiber densities. The power law exponent describing the approximate dependence of the conductivity on the density parameter is found to vary from values close to 1 up to values exceeding 2 when evaluated within different finite ranges of the density parameter. This finding explains the variability of scaling laws for thermal conductivity of fibrous materials suggested in the literature based on numerical simulations performed in different regions of the space of material parameters.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Simulations of biomass pyrolysis using glued-sphere CFD-DEM with 3-D intra-particle models

The shape of biomass particles influences intra-particle heat transfer, particle-particle collisions, interphase momentum, heat transfer, and pyrolysis yield. In this research, the effects of these multi-scale phenomena are captured using a glued-sphere computational fluid dynamic discrete element method (CFD-DEM) with 3-D intra-particle models. The glued-spheres resolved the shapes of particles and intra-particle temperature and species distributions. The implementation was validated using the pyrolysis data of a spherical particle and a cylindrical particle. Then, the influences of biomass shapes and sizes on pyrolysis were investigated. The results revealed the significant influence of surface areas and mass distributions on interphase heat transfer and chemical reactions. The intra-particle heat transfer is neglectable for a Biot number smaller than 0.41. Finally, the pyrolysis of irregularly shaped biomass particles in a fixed bed, fluidized bed, and spouted bed was simulated. Analysis of conversion distributions showed the largest deviations in the fixed bed and the slowest conversion rate in the spouted bed. A good balance between uniform and fast conversion is achieved in the fluidized bed. The multi-scale method proposed in this research provides a tool for the evaluation of large irregular-shaped biomass particle pyrolysis in different types of reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Interpreting T-History Data on Supercooling and Phase Change Enthalpy of a Salt-Hydrogel Complex for Building Thermal Management

Phase change materials (PCMs) are promising materials for storing and discharging latent heat in buildings to regulate the thermal environment. Measurement of salt hydrate PCM thermal performance is often challenging because of supercooling and phase segregation. Though differential scanning calorimetry (DSC) has been widely used to measure thermal properties of various phase change materials, DSC uses milligrams of sample and is unrealistic for predicting large-scale materials performance. The T-history method [1], in contrast, handles higher volumes but the analysis typically requires a lumped capacitance assumption. Here, we report and analyze T-History data on a novel salt-hydrogel complex with high viscosity, that renders low Biot number sample preparation extremely challenging. We re-examine the standard T-history method and compared against DSC data in this complex material. We report details of an experimental setup, built with necessary additions to a controlled environmental chamber. We carefully examine the validity of the lumped capacitance assumption and suggest a modified analysis using a computational approach for reliable data. Experimental parameters such as ramp rate and temperature range of the environmental chamber are considered via computational modeling to understand the effect of ramp rate on supercooling [2]. Finally, we report thermal cycling experiments on the salt-hydrogel complex. Our modified approach to measuring supercooling and enthalpy of fusion at large scales is important in better understanding the performance of phase-change materials at scale for building thermal storage. References: [1] Marin, Jose M, Belen Zalba, Luisa F Cabeza, and Harald Mehling. “Determination of Enthalpy Temperature Curves of Phase Change Materials with the Temperature-History Method: Improvement to Temperature Dependent Properties.” Measurement Science and Technology 14, no. 2 (February 1, 2003): 184–89. https://doi.org/10.1088/0957-0233/14/2/305. [2] Safari, A., R. Saidur, F.A. Sulaiman, Yan Xu, and Joe Dong. “A Review on Supercooling of Phase Change Materials in Thermal Energy Storage Systems.” Renewable and Sustainable Energy Reviews 70 (April 2017): 905–19. https://doi.org/10.1016/j.rser.2016.11.272

thermodynamics↗