Search NASA⌕ Search

SEARCH · Search NASA

Results for “COMSOL multiphysics”

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.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Anomalous Thrust Production from an RF Test Device Measured on a Low-Thrust Torsion Pendulum

This paper describes the eight-day August 2013 test campaign designed to investigate and demonstrate viability of using classical magnetoplasmadynamics to obtain a propulsive momentum transfer via the quantum vacuum virtual plasma. This paper will not address the physics of the quantum vacuum plasma thruster, but instead will describe the test integration, test operations, and the results obtained from the test campaign. Approximately 30-50 micro-Newtons of thrust were recorded from an electric propulsion test article consisting primarily of a radio frequency (RF) resonant cavity excited at approximately 935 megahertz. Testing was performed on a low-thrust torsion pendulum that is capable of detecting force at a single-digit micronewton level, within a stainless steel vacuum chamber with the door closed but at ambient atmospheric pressure. Several different test configurations were used, including two different test articles as well as a reversal of the test article orientation. In addition, the test article was replaced by an RF load to verify that the force was not being generated by effects not associated with the test article. The two test articles were designed by Cannae LLC of Doylestown, Pennsylvania. The torsion pendulum was designed, built, and operated by Eagleworks Laboratories at the NASA Johnson Space Center of Houston, Texas. Approximately six days of test integration were required, followed by two days of test operations, during which, technical issues were discovered and resolved. Integration of the two test articles and their supporting equipment was performed in an iterative fashion between the test bench and the vacuum chamber. In other words, the test article was tested on the bench, then moved to the chamber, then moved back as needed to resolve issues. Manual frequency control was required throughout the test. Thrust was observed on both test articles, even though one of the test articles was designed with the expectation that it would not produce thrust. Specifically, one test article contained internal physical modifications that were designed to produce thrust, while the other did not (with the latter being referred to as the "null" test article). Test data gathered includes torsion pendulum displacement measurements which are used to calculate generated force, still imagery in the visible spectrum to document the physical configuration, still imagery in the infrared spectrum to characterize the thermal environment, and video imagery. Post-test data includes static and animated graphics produced during RF resonant cavity characterization using the COMSOL Multiphysics® software application. Excerpts from all of the above are included and discussed in this paper. Lessons learned from test integration and operations include identification of the need to replace manual control of the resonant cavity target frequency with an automated frequency control capability. Future test plans include the development of an automatic frequency control circuit. Test results indicate that the RF resonant cavity thruster design, which is unique as an electric propulsion device, is producing a force that is not attributable to any classical electromagnetic phenomenon and therefore is potentially demonstrating an interaction with the quantum vacuum virtual plasma. Future test plans include independent verification and validation at other test facilities.

Brady, David↗
NASA NTRS2015

Magnetostrictive Vibration Damper and Energy Harvester for Rotating Machinery

Vibrations generated by machine driveline components can cause excessive noise and structural damage. Magnetostrictive materials, including Galfenol (iron-gallium alloys) and Terfenol-D (terbium-iron-dysprosium alloys), are able to convert mechanical energy to magnetic energy. A magnetostrictive vibration ring is proposed, which generates electrical energy and dampens vibration, when installed in a machine driveline. A 2D axisymmetric finite element (FE) model incorporating magnetic, mechanical, and electrical dynamics is constructed in COMSOL Multiphysics. Based on the model, a parametric study considering magnetostrictive material geometry, pickup coil size, bias magnet strength, flux path design, and electrical load is conducted to maximize loss factor and average electrical output power. By connecting various resistive loads to the pickup coil, the maximum loss factors for Galfenol and Terfenol-D due to electrical energy loss are identified as 0.14 and 0.34, respectively. The maximum average electrical output power for Galfenol and Terfenol-D is 0.21 W and 0.58 W, respectively. The loss factors for Galfenol and Terfenol-D are increased to 0.59 and 1.83, respectively, by using an L-C resonant circuit.

Magnetostrictive↗
NASA NTRS2015

The Effect of Resistance on Rocket Injector Acoustics

Combustion instability, where unsteady heat release couples with acoustic modes, has long been an area of concern in liquid rocket engines. Accurate modeling of the acoustic normal modes of the combustion chamber is important to understanding and preventing combustion instability. This study evaluates the effect of injector resistance on the mode shapes and complex eigen-frequencies of an injector/combustion chamber system by defining a high Mach-flow form of the convective wave equation (see Eq. 1) in COMSOL Multiphysics' Coefficient Form PDE Mathematics Module.

Morgan, C. J.↗
NASA NTRS2016

My Internship at NASA

My name is Isaac Lopez and I am a junior at the University of Houston majoring in Mechanical Engineering Technology. I will be completing my first tour at the NASA-Johnson Space Center ("JSC") as a Mechanical Engineer within the Human Interfaces Branch. Throughout my tour, I was given the opportunity to work on multiple projects that have expanded my knowledge and interest in acoustics and engineering design. One of the projects I worked on at JSC consisted of doing acoustic simulation of the EVA comm. cap. While working on the comm. cap headset, my main duty consisted of simulating the acoustics of the headset to find a solution to the condensing water that can accumulate and block the acoustic tube, causing attenuation or complete loss of audio in one ear for an astronaut using the EVA. For this project, I had to create a Creo model of the comm. cap so that I would be able to import it into Comsol for acoustic simulation. I also had the opportunity to design a portable and lightweight beam degrader for the EEE Parts and Radiation team. With the help of Creo, I was able to make a CAD design and put together a small working prototype for the radiation team to demonstrate the capabilities that the beam degrader had. In addition to these projects, JSC allowed me to work closely on projects with other interns. I had the opportunity to help another intern with his acoustic diverter, intended to improve the sound quality in Node 1 of the ISS. During this project, I helped with some of the acoustic testing inside the anechoic chamber as well as helping record data during testing at the ISS mock up. During the course of my first tour, I was able to learn and continually improve on my CAD drafting skills. With each project I worked on, I acquired new ways to create and improve various designs with various constraints. Furthermore, I also had the opportunity to work with electrical engineers and learn about the electronic components that would provide control of the beam degrader functions. I was able to take the knowledge, skills and designs I had acquired during my tour at JSC (such as Finite Element Analysis) and apply them to my work on the Comsol Multiphysics software for acoustics simulation. I was also able to use it to do some stress analysis on a proposed design for a proposed Habitat structure. My time at JSC has been invaluable as it helped me grow and more fully-comprehend my major as well as spark an interest in acoustics and how much it relates to mechanical engineering. This internship has taught me that engineering is much more than a profession - it is embedded in our daily lives as a puzzle that requires critical thinking as an everyday challenge that makes life more interesting.

Lopez, Isaac↗
NASA NTRS2016

A Virtual Laboratory for the 4 Bed Molecular Sieve of the Carbon Dioxide Removal Assembly

Ongoing work to improve water and carbon dioxide separation systems to be used on crewed space vehicles combines sub-scale systems testing and multi-physics simulations. Thus, as part of NASA's Advanced Exploration Systems (AES) program and the Life Support Systems Project (LSSP), fully predictive COMSOL Multiphysics models of the Four Bed Molecular Sieve (4BMS) of the Carbon Dioxide Removal Assembly (CDRA) on the International Space Station (ISS) have been developed. This Virtual Laboratory is being used to help reduce mass, power, and volume requirements for exploration missions. In this paper we describe current and planned modeling developments in the area of carbon dioxide removal to support future missions as well as the resolution of anomalies observed in the ISS CDRA.

Coker, Robert↗
NASA NTRS2017

Investigation of Vapor Cooling Enhancements for Applications on Large Cryogenic Systems

The need to demonstrate and evaluate the effectiveness of heat interception methods for use on a relevant cryogenic propulsion stage at a system level has been identified. Evolvable Cryogenics (eCryo) Structural Heat Intercept, Insulation and Vibration Evaluation Rig (SHIIVER) will be designed with vehicle specific geometries (SLS Exploration Upper Stage (EUS) as guidance) and will be subjected to simulated space environments. One method of reducing structure-born heat leak being investigated utilizes vapor-based heat interception. Vapor-based heat interception could potentially reduce heat leak into liquid hydrogen propulsion tanks, increasing potential mission length or payload capability. Due to the high number of unknowns associated with the heat transfer mechanism and integration of vapor-based heat interception on a realistic large-scale skirt design, a sub-scale investigation was developed. The sub-project effort is known as the Small-scale Laboratory Investigation of Cooling Enhancements (SLICE). The SLICE aims to study, design, and test sub-scale multiple attachments and flow configuration concepts for vapor-based heat interception of structural skirts. SLICE will focus on understanding the efficiency of the heat transfer mechanism to the boil-off hydrogen vapor by varying the fluid network designs and configurations. Various analyses were completed in MATLAB, Excel VBA, and COMSOL Multiphysics to understand the optimum flow pattern for heat transfer and fluid dynamics. Results from these analyses were used to design and fabricate test article subsections of a large forward skirt with vapor cooling applied. The SLICE testing is currently being performed to collect thermal mechanical performance data on multiple skirt heat removal designs while varying inlet vapor conditions necessary to intercept a specified amount of heat for a given system. Initial results suggest that applying vapor-cooling provides a 50 heat reduction in conductive heat transmission along the skirt to the tank. The information obtained by SLICE will be used by the SHIIVER engineering team to design and implement vapor-based heat removal technology into the SHIIVER forward skirt hardware design.

small scale testing↗
NASA NTRS2017

Virtual Design of a 4-Bed Molecular Sieve for Exploration

Simulations of six new 4-Bed Molecular Sieve configurations have been performed using a COMSOL (COMSOL Multiphysics - commercial software) model. The preliminary results show that reductions in desiccant bed size and sorbent bed size when compared to the International Space Station configuration are feasible while still yielding a process that handles at least 4.0 kilograms a day CO2. The results also show that changes to the CO2 sorbent are likewise feasible. Decreasing the bed sizes was found to have very little negative effect on the adsorption process; breakthrough of CO2 in the sorbent bed was observed for two of the configurations, but a small degree of CO2 breakthrough is acceptable, and water breakthrough in the desiccant beds was not observed. Both configurations for which CO2 breakthrough was observed still yield relatively high CO2 efficiency, and future investigations will focus on bed size in order to find the optimum configuration.

Giesy, Timothy J.↗
NASA NTRS2019

Numerical Model of the Mars Electrostatic Precipitator

NASA's future human exploration missions will require chemical processing plants to convert local resources into consumables to support astronaut activities. The thin and mostly carbon dioxide atmosphere of Mars is estimated to have 1 - 10 particles/cu.cm with diameters of 1 - 10 m and up to 1000 particles/cu.cm during storms. The dust in the Martian atmosphere can foul chemical reactors and pose a risk to life support systems. Electrostatic precipitation (ESP) removes dust particles from the Martian atmosphere. The Electrostatics and Surface Physics Laboratory at NASA's Kennedy Space Center has developed a COMSOL Multiphysics(Registered Trademark) model of an ESP for dust filtration on Mars. The fundamental principles of an ESP can be simulated by four physics modules: plasma, AC/DC electromagnetics, computational fluid dynamics (CFD), and particle tracing. In the ESP model presented here, the plasma module was solved to estimate particle charge. The AC/DC and CFD module were solved for the electrostatic force and fluid force. The particle-tracing module was solved for particle collection efficiency.

Wang, Jerry J.↗
NASA NTRS2009

Understanding how molecular hydrogen impacts the total dose and dose rate response of linear bipolar circuits

Recent Enhanced Low Dose Rate Sensitivity (ELDRS) investigations carried out by RLP Research, Arizona State University (ASU) and Jet Propulsion Laboratory (JPL) have shown significant differences in the degradation of bipolar micro-circuits with total dose in the presence of molecular hydrogen (H2) in packages. This has a significant impact on radiation hardness assurance and opens up opportunities to improve device performance. The general objectives of this program are to: 1. Determine the extent to which hydrogen contamination affects the total dose and dose rate response of linear bipolar circuits; 2. Develop a model that will enable the prediction of high dose rate (HDR) and low dose rate (LDR) response asymptotes and transition dose rates as a function of total dose, temperature, pressurized hydrogen, defect precursors, and other process dependent variables; 3. Explore the possibility of an accelerated hardness assurance method and possible hardening approaches; and 4. Extend the work to other technologies that have total dose response affected by hydrogen contamination.In this document, we xperimentally demonstrate with test transistors and circuits that hydrogen is correlated with ELDRS in bipolar linear circuits. We show that the amount of hydrogen determines: 1) the total dose response versus dose rate and 2) the transition dose rates between the high and low dose rate responses. The experimental results are supported with a steady state drift/diffusion analytical model as well as modeling calculations using COMSOL Multiphysics.

Pease, Ron↗
NASA NTRSDate not supplied

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

Additive Manufacturing↗
NASA NTRSDate not supplied

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

pragmatic method↗
NASA NTRS2021

Analysis of Arc Welding Process in Space

This work is motivated by NASA plans to conduct welding experiments on the ISS. It is expected that deployment of welding and additive manufacturing technologies in the space environment has the potential to revolutionise how orbiting platforms are designed, manufactured, and assembled. However, the structure, composition and quality of a weld is extremely dependent on the environment and can be difficult to control in space. Shielding gases would also be tough to manage as gases behave differently in zero gravity and airless environments. Additional points of concern are related to the spatter and sparks dynamics in space. Therefore, there is a need for a more basic understanding of welding processes by computational modelling. To provide such an insight we developed state of art models of the ARC torch, droplet detachment and transfer, and the melt pool build up using magnetohydrodynamics approximation and level set method for two-phase liquid metal/-gas flow modelling. Two finite element models were built using 2D axisymmetric geometry in COMSOL Multiphysics®: (i) stationary model of the ARC torch and (ii) dynamical model of droplet detachment and transfer. Both models demonstrate reasonable agreement with earlier experimental observations and high sensitivity to the temperature dependence of the thermophysical parameters on the system materials. The models were used to provide physical insight into ARC welding in various environments and geometries. It was shown, in particular, that there is a significant probability of gas bubbles trapping in the meltpool and a possibility of unbounded wondering of sparks in the welding chamber in zero gravity. In addition, we estimated metal evaporation rate that may be hazardous in the confined environment of the ISS. These issues raise concern of quality of the weld and safety for ISS applications.

ARC welding↗
NASA NTRSDate not supplied

Numerical Modelling of Sublimation Losses During Sample Handling Within the Light WAVE Lunar Water Prospecting Instrument

NASA’s focus of exploring the Moon’s south pole under the Artemis program has placed new emphasis on characterizing the quantity and distribution of water ice and other volatiles present within the Moon’s permanently shadowed regions (PSRs), first measured during the 2009 LCROSS mission. NASA’s VIPER rover is the first of a series of prospecting missions to characterize this water ice resource. NASA’s Johnson Space Center (JSC) has developed the Light Water Analysis and Volatiles Extraction (Light WAVE) instrument shown in Figure 1 for a follow-on prospecting mission, designed to be a rover payload for water ice prospecting at various depths and locations within PSRs. A key source of error for lunar water prospecting instruments, including Light WAVE, is the loss of water in an icy regolith sample to the lunar vacuum via sublimation due to heating above its ice-stable temperature. This abstract and conference presentation discuss a sublimation model for Light WAVE’s sample crucible using COMSOL Multiphysics to assess the sublimation losses that would occur between sample collection from an external drill and sealing the sample within Light WAVE’s boiler.

Lunar Polar Water↗
NASA NTRSDate not supplied

Design and Analysis of Battery Thermal Management Systems

Thermal management of battery cell packs is a critically needed technology. The purpose of this work is to design new and improved Battery Thermal Management Systems (BTMS) for use in electric airplanes. The BTMS should be 3D printable and hold twelve to sixteen 18650 batteries. To minimize aircraft’s weight, an actively air-cooled battery pack was chosen to avoid the excess weight of water-cooled and phase change material (PCM) cooled packs. As a result, significant changes were made to remove all metal in the pack and replace it with lower density polymer matrix composites (PMCs) which can be additively manufactured. Two air cooled battery pack configurations (traditional propeller fans & bladeless fans) were designed, modeled, and compared. These packs were first modeled in SolidWorks 2021 3D CAD, then imported into COMSOL MultiPhysics to be studied using the “Heat Transfer in Solids and Fluids” module. The design with the bladeless fans eliminated the need to use high conductivity heavy metal to remove unwanted thermal energy. These bladeless fans were designed entirely out of PMCs. This thermal pack design weighs 0.04 kg less than that of the traditional propeller fans design and has increased the battery pack energy density by 8.25 Wh/kg.

Thermal Management System for Battery Packs↗
NASA NTRS2024

University of South Alabama: Evaluation of Using Heat Pipes to Cool Adsorbent Beds in AC-TSAC Systems

This report details the results of the research and work of undergraduate students at the University of South Alabama across the 2023-24 academic year. These results examine the use of heat pipes for air revitalization in space exploration and habitation units. The current mechanism for removing carbon dioxide in the International Space Station is the Carbon Dioxide Removal Assembly (CDRA), which involves 2 cooperating adsorbent beds that continuously adsorb and desorb carbon dioxide. Recently, NASA has investigated the use of air-cooled temperature swing adsorption compressors (AC-TSAC) to provide both adsorption and compression of CO 2 . To effectively adsorb and desorb CO 2 , the beds must operate across a range of temperatures and require efficient heating and cooling to reach the necessary operating temperatures. Heat pipes, which are fluid-filled rods that use a wicking mechanism to transfer heat, may provide one alternative approach to moving heat into and out of adsorption beds. To test the cooling efficiency of packed beds equipped with heat pipes, a packed bed was constructed to examine the amount of heat transferred with this technology. In addition to a physical apparatus, a model was built using COMSOL Multiphysics® to obtain an initial assessment of the effectiveness of using heat pipes for temperature control in porous adsorption beds.

T. Grant Glover↗
NASA NTRSDate not supplied

System Modeling of a Lunar Molten Regolith Electrolysis Plant

Introduction: In-Situ Resource Utilization (ISRU) is the process of extracting local resources to produce commodities for propulsion, life support systems, and off-planet construction rather than transporting consumables from Earth. Molten Regolith Electrolysis (MRE) is a novel ISRU method of extracting oxygen gas and metal alloy from lunar regolith. The MRE process involves placing lunar regolith between two electrodes, through which current is passed, to melt the regolith and reduce the metal oxide constituents by direct electrolysis (e.g. FeO, SiO2, MgO, Al2O3) into oxygen gas and metal alloys. The oxygen is liquefied and used as propellant for landers, while the metals (e.g. Ferro-alloys) are further processed and used in structural building materials and parts manufacturing. A system model was developed that accounted for the major processes of an MRE plant (from excavation of raw materials to storage of products) to assess the feasibility of a lunar MRE plant. The System Engineering and Integration (SE&I) ISRU Modeling and Analysis (SIMA) team utilized its previously documented system sizing model, the Mission Analysis and Integration Tool (MAIT) [1] as framework of the system model. MAIT uses MATLAB/Simulink to integrate subsystem models into a complete system model of the MRE plant. Total mass, volume, and power requirements were computed for numerous iterations of a MRE plant. System Model: Figure 1: MRE Plant Block Diagram The regolith excavation model determines the mass and power needed to excavate sufficient regolith. The preheating auger initiates the regolith heating process before regolith enters the MRE re-actor to reduce the energy required to turn the solid into a molten liquid. The MRE reactor is modeled in COMSOL Multiphysics and based on the research by Dominguez, Sibille, and Schreiner [2, 3, 4]. This preliminary reactor model provides an accurate calculation of thermal equilibrium during electrochemical operation of the reactor system to assess the optimal mass and power required to process the inlet flow of regolith. The model also computes the outlet flowrates of oxygen and molten products. For this analysis, the primary components of the metal alloy considered were iron and silicon. The oxygen is then purified using an Yttrium Stabilized Zirconia (YSZ) electrode, followed by liquefaction using a 90K cryocooler to be stored as liquid oxygen in insulated cylindrical tanks. In future iterations of the system model, the molten metal tapped from the MRE reactor will undergo additional processing or refinement. However, downstream handling of metals is currently a technology gap that is missing a high TRL subsystem model. Therefore, for this analysis, the accumulated metal alloy stream terminates after leaving the MRE reactor. Study Goals: This analysis investigates multiple input variables to the system to determine the sensitivity of a (near) complete plant at full-scale. This preliminary investigation ran parametric sweeps on the MRE reactor geometry, electrical current supply, layers of multi-layer insulation (MLI) on the reactor, size of the electrodes in the oxygen purification model, and regolith composition (based on landing site location). Three production targets of oxygen (1,000, 10,000, and 50,000 kg/yr) were investigated for this analysis. The parametric sweeps conducted in this analysis provide valuable insight into the expected impact of the various model inputs on plant size. This information can be used to identify the most critical components of the plant and guide future decisions on allocating funding for research and development, providing subsystem developers with appropriate interfaces with downstream and upstream processes, and assessing the overall feasibility of MRE when compared to other ISRU plants. References: [1] Carlson, A. et al. (2024) ICES, ICES-2024-53. [2] Dominguez, D.A., and Sibille, L. (2011) AIAA, AIAA-2011-700. [3] Schreiner, S.S. (2015) MIT, Dissertation. [4] Schreiner, S.S. et al. (2016) ASR, 57(7), pp.1585-1603.

ISRU↗
DOE OSTI2021

An HPC-Based Hydrothermal Finite Element Simulator for Modeling Underground Response to Community-Scale Geothermal Energy Production

Geothermal heat, as renewable energy, shows great advantage with respect to its environmental impact due to its significantly lower CO2 emissions than conventional fossil fuel. Open and closed-loop geothermal heat pumps, which utilize shallow geothermal systems, are an efficient technology for cooling and heating buildings, especially in urban areas. Integrated use of geothermal energy technologies for district heating, cooling, and thermal energy storage can be applied to optimize the subsurface for communities to provide them with multiple sustainable energy and community resilience benefits. The utilization of the subsurface resources may lead to a variation in the underground environment, which might further impact local environmental conditions. However, very few simulators can handle such a highly complex set of coupled computations on a regional or city scale. We have developed high-performance computing (HPC) based hydrothermal finite element (FE) simulator that can simulate the subsurface and its hydrothermal conditions at a scale of tens of km. The HPC simulator enables us to investigate the subsurface thermal and hydrologic response to the built underground environment (such as basements and subways) at the community scale. In this study, a coupled hydrothermal simulator is developed based on the open-source finite element library deal.II. The HPC simulator was validated by comparing the results of a benchmark case study against COMSOL Multiphysics, in which Aquifer Thermal Energy Storage (ATES) is modeled and a process of heat injection into ATES is simulated. The use of an energy pile system at the Treasure Island redevelopment site (San Francisco, CA, USA) was selected as a case study to demonstrate the HPC capability of the developed simulator. The simulator is capable of modeling multiple city-scale geothermal scenarios in a reasonable amount of time.

Sun, Xiang↗
DOE OSTI2022

Experimental and numerical study on the impacts of various long carbon fiber inserts on the thermal performance of phase change materials

In thermal energy storage systems, phase change materials (PCMs) are widely used to provide efficient storing and extraction of thermal energy. However, most PCMs have low thermal conductivities, which limit the heat transfer of the latent energy inward/outward. In this study, the impacts of various fiber inserts on the solidification and melting processes with PCM were investigated by experiments and simulations. It was found that long anisotropic carbon fibers with high thermal conductivities can not only increase the heat transfer within PCM by enlarging hot-cold interface areas, but also act as crystallization nuclei and increase the solidification of PCM during the charging process. Recalescence of PCM during solidification was observed. COMSOL models were used to understand the thermal transfer mechanisms. COMSOL Multiphysics models of a single carbon fiber filament in PCM showed that the anisotropy in thermal conductivities controlled the radial and axial heat transfer ratios along its fiber axis and can enhance both the solidification and melting processes.

Zhou, Xinzhang↗