Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “COMSOL”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Thermally Anisotropic Composites for Heat Redirection and Thermal Management in Building Envelopes

In building envelope systems, thermal management is important from energy conservation and thermal comfort perspectives, and is typically done via insulation, thermal mass, solar control and shading, etc. Thermal anisotropy is commonly utilized in electronics for heat dissipation. This article investigates the feasibility of utilizing thermally anisotropic composites (TACs) to reduce unwanted heat gains and losses through the building envelope by redirecting the heat transfer to or from appropriate heat sinks and sources.Numerical simulations of building envelope components using COMSOL Multiphysics and whole building simulations using EnergyPlus (E+) were performed to estimate the energy savings potential of TACs used in external walls of residential buildings. Application of TACs was compared with baseline walls with cavity insulation as well as with and without exterior continuous insulation (CI). The COMSOL simulations indicated 75% or higher reductions in heat gains when using a TAC coupled with a heat sink compared to an exterior CI case. In terms of energy savings in a hot climate, E+ models estimated cooling energy savings of 19% with a "TAC + heat sink" compared to a baseline wall with only cavity insulation. Further, a small-scale test sample with a TAC-heat sink system was built and tested in a multi-transducer heat flow meter (HFM). The heat sink was a copper tube circulating cold water. The experiments verified the re-direction of heat to a sink and significant reduction in heat transfer through the test specimen. The experiments showed that the heat redirection only occurred in the presence of both thermal anisotropy and water circulation; if either element was missing, the system behaved like a traditional envelope component.

Biswas, Kaushik↗
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↗
DOE OSTI2021

Radiative Heat Transport and Optical Characterization of High Temperature Molten Salts

The goal of this project was to advance the field of molten fluoride radiative heat transfer (RHT), through optical property measurements, heat transfer loop construction and operation, and computational fluid dynamics simulations. The following outcomes have been accomplished. A comprehensive review has been conducted to compile all available RHT-relevant optical absorption data for molten fluorides. The various measurement methods have also been compiled and critiqued. Thermal emission measurement methods have been refined and documented. A unique thermal emissivity measurement facility has been constructed and used to measure the emissivity of flowing FLiBe-exposed SS316 samples. A slight increase has been observed in the emissivity compared to unexposed SS316. A dual reflectance-emission setup was partially demonstrated for the purpose of molten fluoride absorption measurement. Preliminary infrared edge data for FLiNaK was obtained. A visible near-infrared (NIR) measurement setup has been constructed and will be used to measure absorption of corrosion products in FLiBe. A new multi-band thermal radiation method using blackbody-weighted spectral bands for the absorption-coefficient was deployed to investigate the thermal performance of the UW Natural Circulation FLiBe Loop under various power levels and operating conditions. The results showed an increasing and decreasing trend in the overall heat transfer with optical thickness with a maximum at optical thickness τ D ≈Ο(1). A forced convection FLiBe loop was constructed for the purpose of studying corrosion and thermal hydraulic phenomena. A vertical cantilever pump is used to achieve up to 5 gpm and 5 m/s flow velocity in hot and cold test sections. A methodology has been established for modeling radiative heat transfer (RHT) in molten salt coolants for the purpose of quantifying the contribution of thermal radiation in conditions relevant to FHR or MSR systems. The comparison between the results obtained for forced and mixed convection showed that in the studied ranges of optical thickness, the maximum effect of RHT on overall heat transfer was 60% improvement observed at z/D = 250 for forced convection; whereas for mixed convection, that value was higher at 170% observed at z/D = 30. The highest effect of RHT on the overall energy distribution was observed at τ D ≈3.4. Negligible effect (<20%) was observed above τ D =18 and below τ D ≈0.002. A CFD code-to-code comparison was conducted using COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT, for a forced and mixed convection laminar flow for a grey medium in a vertical heated using both COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT. The relative error was the largest at τ D =40 for the radiative heat flux, bulk temperature, and Nusselt number. On the other hand, the relative error was the least at τ D =0.002 for the conductive heat flux, bulk temperature, and Nusselt number.

Anderson, Mark↗
DOE OSTI2022

A Machine Learning-Assisted Framework to Control Thermally Anisotropic Building Envelopes in Residential Buildings

To curb the energy consumption of buildings and their related CO2 emissions, Oak Ridge National Laboratory (ORNL) has developed the thermally anisotropic building envelope (TABE) —a multi-layer design comprising insulation materials and metal foils connected to thermal loops. In this study, a machine learning-assisted framework was developed to control the TABE in residential buildings to reduce the computation load for future optimal rule-based control and application. First, a 2D finite element model was established in COMSOL to calculate the hourly heat flux through exterior walls installed with the TABE. Then, TABE wall heat fluxes were simulated for various given indoor and outdoor boundary conditions, as well as thermal loops fluid temperatures and flow rates. Since the finite element simulations are computationally expensive, an artificial neural network (ANN) was then trained to use as a proxy of the finite element (COMSOL) modeling. Finally, the trained ANN model was coupled with the EnergyPlus model to predict the energy consumption of a US Department of Energy prototype single-family house installed with the TABE. An optimal simple rule-based control was determined from predefined rules for a case study. The results demonstrate that the developed machine learning–assisted framework can reduce 99.9% of the computation time while efficiently managing residential building energy for installed TABE walls.

Shen, Zhenglai↗
DOE OSTI2023

Simulation of gallium nitride vertical fin-shaped field effect transistor for use as thermal neutron detector

Through the use of a radiation detection system simulation framework, a gallium nitride vertical fin-shaped field effect transistor (FinFET) was studied for output response when utilized as a thermal neutron detector. The FinFET was assumed to had been backfilled with boron carbide, reactive to thermal neutrons. The GaN FinFET was modeled with radiation transport from MCNP, and the electronic transport from COMSOL Multiphysics. Fabricated FinFET devices (not neutron reactive) were tested to aid in the tuning of the COMSOL FinFET model. Through time-depenent studies, the drain current response pulse to simulated ionization due to single events lead to building of a database of device responses to radiation. By integrating the current pulses over time, the induced charge was calculated. Using the results of the radiation transport PTRAC file in combination with the induced charge database, an integrated charge spectrum was calculated.

Davidson, Bryce L.↗
DOE OSTI2026

Two-Way Coupled CFD/DEM Simulation of Electromagnetic Heating in a Fluidized Bed Reactor - Extended Abstract

A set of tools has been developed to couple a multiphase reacting computational fluid dynamics solver (MFIX) with two computational electromagnetic solvers (COMSOL and Elmer). In contrast to several recent investigations, which calculate a single static EM field prior to the CFD calculation, in the present work the electromagnetic field is periodically updated during the simulation. This enables the EM field and associated EM heating predictions to reflect the effects of evolving material properties (temperature, composition, size etc.) and the evolving distribution of material in the system. The present work uses effective medium theory (EMT) for the electromagnetic field calculations since the goal is to apply the tool for industrial scale reacting systems. Predictions using MFIX/COMSOL and MFIX/Elmer are compared to in bed thermocouple measurements from a fluidized bed experiment and recent simulation predictions and recent experiments performed at NETL. The NETL experiments also include packed bed operating conditions. Coupled MFIX/COMSOL simulations at these conditions are compared with similar simulations using only COMSOL to measure errors due to the solver coupling approach.

CFD↗
DOE OSTI2022

An SECM-Based Spot Analysis for Redoxmer-Electrode Kinetics: Identifying Redox Asymmetries on Model Graphitic Carbon Interfaces

The fundamental process in non-aqueous redox flow battery (NRFB) operation revolves around electron transfer (ET) between a current collector electrode and redox-active organic molecules (redoxmers) in solution. Here, we present an approach utilizing scanning electrochemical microscopy (SECM) to evaluate interfacial ET kinetics between redoxmers and various electrode materials of interest at desired locations. This spot-analysis method relies on the measurement of heterogeneous electron transfer rate constants (k f or k b ) as a function of applied potential (E-E 0 '). As demonstrated by COMSOL simulations, this method enables the quantification of Butler-Volmer kinetic parameters, the standard heterogeneous rate constant, k 0 , and the transfer coefficient, α. Our method enabled the identification of inherent asymmetries in the ET kinetics arising during the reduction of ferrocene-based redoxmers, compared to their oxidation which displayed faster rate constants. Similar behavior was observed on a wide variety of carbon electrodes such as multi-layer graphene, highly ordered pyrolytic graphite, glassy carbon, and chemical vapor deposition-grown graphite films. However, aqueous systems and Pt do not exhibit such kinetic effects. Our analysis suggests that differential adsorption of the redoxmers is insufficient to account for our observations. Displaying a greater versatility than conventional electroanalytical methods, we demonstrate the operation of our spot analysis at concentrations up to 100 mM of redoxmer over graphite films. Looking forward, our method can be used to assess non-idealities in a variety of redoxmer/electrode/solvent systems with quantitative evaluation of kinetics for applications in redox-flow battery research.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Numerical modeling of electromagnetic field spatiotemporal evolution to evaluate the effects on calcium carbonate crystallization

Calcium carbonate (CaCO 3 ) scaling is a significant impediment to water systems. Electromagnetic field (EMF) treatment is a promising approach to control scaling owing to its simplicity and low or no energy requirements. However, the underlying mechanisms by which EMF impacts CaCO 3 crystallization remain unclear due to the challenges in measuring the EMFs in feed solutions and the lack of a fundamental understanding of the applied EMFs and the observed physicochemical phenomena. To fill this knowledge gap, a high-fidelity COMSOL model was first developed to simulate EMFs in bulk solutions for three alternating current-induced EMF devices with different configurations and properties. These were then integrated with experimental data to unveil the underlying mechanism by which applied EMFs alter the physicochemical processes. The study revealed that even low-strength EMFs (e.g., electric fields <0.15 V/m and magnetic fields <0.03 mT) promoted CaCO 3 precipitation in bulk solutions. The electric fields created by these EMF devices resulted in higher Lorentz force compared to their induced magnetic fields. In conclusion, the methodology of this study offers the capability to predict the effectiveness of different EMF devices in facilitating crystallization processes, and these mechanistic insights lay the foundation for the smart design of EMF devices for diverse water treatment applications.

COMSOL simulation↗
DOE OSTI2025

Enhancing EnergyPlus capabilities to model dynamic building envelopes using python plugin

Nearly half of the energy consumption in the United States is related to buildings, resulting in an urgent need to develop innovative technologies to improve building energy efficiency. Dynamic building envelopes, comprising switchable insulation and thermal energy storage materials, have been proposed recently as a promising solution to reduce buildings' heating and cooling loads by thermally coupling the indoor environment with the ambient environment when beneficial while decoupling them when outdoor conditions are not favorable. Although various related technologies are still underway, the whole-building energy modeling tools, like EnergyPlus, do not have the capability to simulate the transient and dynamic nature of dynamic envelope materials and components to accurately predict their impact on building energy use. The objective of this study is to formulate a method in EnergyPlus simulation engine to model multilayer envelopes, comprising dynamic building materials with variable thermophysical properties, and discuss the changes made to the program using a Python plugin. Furthermore, the thermal performance of the dynamic envelopes using the proposed method is compared and verified with the results from a well-established commercial code, COMSOL Multiphysics. A parametric assessment is also conducted to evaluate the energy efficiency benefits of dynamic envelopes in a single-family residential building, demonstrating total annual energy savings up to 11.6 %, when a dynamic envelope operates alone, and up to 18.2 % when it is combined with a thin layer of phase change material as a thermal storage medium. Finally, a United States wide energy efficiency assessment is presented to showcase the geographical spread of the energy savings. The method designed and implemented in this study provides the researchers with the ability to implement their dynamic insulation methods in EnergyPlus and evaluate the whole building energy impact.

25 ENERGY STORAGE↗
DOE OSTI2020

2-D moving mesh modeling of lithium dryout in open surface liquid metal flow applications

Liquid lithium displays increasing promise as a replacement for solid plasma facing components (PFCs) in fusion device applications. The Liquid Metal Infused Trench (LiMIT) system, developed at the University of Illinois (UIUC), has demonstrated how thermoelectric magnetohydrodynamics (TEMHD) can be harnessed to drive liquid lithium flow in an open surface PFC. However, in the highest heat flux applications, large local acceleration is created via TEMHD, and the sudden increase in velocity can cause the liquid level to expose the underlying solid, eliminating the protective benefits of the lithium. In order to study potential mitigation strategies, a 2-D COMSOL Multiphysics model was developed using the moving mesh module to capture free surface flow. The model depicts the development of the dryout phenomenon for 2 test cases – slow (1 cm/s) and medium (10 cm/s) flow in 5 mm deep trenches – including the liquid level reduction under the high heat flux and the pileup of slower flow both upstream and downstream of the heat stripe. The effectiveness of trench shaping dryout mitigation strategies is examined. For the slow flow case, it is shown that a 1.8 mm ledge placed under the heat stripe will stop dryout, and for the medium flow case, a 2.7 mm ledge is required to mitigate the effect. This model can be used to identify strategies for increasing the viable heat load for open surface liquid lithium PFCs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A hybrid numerical and machine learning framework for evaluating the performance of a 780 cm 2 aqueous organic redox flow battery

Aqueous organic redox flow battery (AORFB) is a promising cost-competitive technology for large-scale energy storage. Among existing work, the dihydroxyphenazine (DHP)-based AORFB has demonstrated high energy density and low capacity degradation in 10 cm2 cells during lab tests. However, its commercial-scale performance in more complex environments remains unknown, posing a barrier for commercialization. To address this gap, this work presents a comprehensive performance evaluation of a 780 cm 2 DHP-based AORFB by combining physics-based numerical model, machine learning (ML)-based surrogate models, and ML-derived sensitivity quantification. Specifically, we first select 12 key battery parameters that include 10 physicochemical quantities and 2 operation quantities, then select 6 performance metrics that include energy efficiency (EE), discharging capacity, charging energy, and power losses due to concentration, activation, and ohmic over-potentials. With such selection, 12800 combinations of the 12 parameters are subsequently generated using the Latin Hypercube Sampling method. These combinations, together with 38 pre-defined State of Charge, are then integrated to a validated AORFB model developed in COMSOL to compute the performance metrics. With both input parameters and performance metrics, 60 deep neural network (DNN) surrogate models are then trained to approximate the relationship between the 10 physicochemical quantities and 6 performance metrics at each flow rate and current density. Sensitivity scores are then calculated based on the DNN models. Two additional sensitivity analysis tools, i.e., MARS, and SHAP, are also used to cross-validate the sensitivity scores from the DNN. The results demonstrate that 1) the standard potential ranks the first in controlling EE and charging energy, 2) the membrane conductivity is most critical for power loss and EE, and 3) specific area and reaction rate control activation power loss.

25 ENERGY STORAGE↗
DOE OSTI2023

Progress towards the FRIB-EDM 3 -Frontend: A tool to provide radioactive molecules from isotope harvesting for fundamental symmetry studies

The under-construction FRIB-EDM 3 -instrument was designed to study polar radioactive molecules (such as RaF) in transparent cryogenic solids by laser spectroscopy. The instrument is divided into a frontend- and a backend section. The frontend accepts an aqueous sample from isotope harvesting and provides a mass-separated molecular ion beam in an ultra-high vacuum environment. In the backend, the ions are guided into alkali-metal vapor and the resulting neutrals are co-deposited in a solid argon matrix to perform laser spectroscopy. Here, this work addresses the frontend of the instrument. The efficient ionization of harvested radioisotopes from aqueous samples is achieved with a spray-ionization method. Subsequently, the molecular ion beam is analyzed by mass-to-charge ratio by a quadrupole mass filter. To verify the feasibility of the approach, numerical simulations with the COMSOL and SIMION packages have been conducted. While the former was applied to study transport in ion funnels, the latter was used to investigate ion beam transmission through the lower pressure sections. Following promising simulation results, a first experimental setup is under construction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Multiphysics modeling of printed surface acoustic wave thermometers

Surface acoustic wave (SAW) transducers are a class of sensors and actuators that operate on the fundamental principle of piezoelectricity. Aerosol jet printing and other additive manufacturing techniques have resulted in the low-cost fabrication of low-power and small-footprint SAW devices that are suitable for sensing in high-temperature and radioactive environments. In this work, we developed a series of temperature-dependent finite element models for a SAW transducer consisting of printed silver interdigitated transducers (IDTs) deposited onto piezoelectric lithium niobate. Modeling accuracy was evaluated experimentally from room temperature to 200 °C using an aerosol-jet-printed SAW thermometer. A time-domain study enabled visualization of the wave propagation and successfully guided the denoising of the scattering parameter measurement. Additionally, frequency-domain models using traditional modal analysis or the unique port boundary condition feature in COMSOL Multiphysics accurately predicted the temperature-driven natural frequency drift in the SAW thermometer. The finite element models developed in this study serve to facilitate the computer-aided design of future SAW transducers for applications in harsh environments.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Predicting the effects of thermally-induced density gradients on the hydrodynamic behavior of PBX 9502

High explosives are often exposed to thermal and mechanical stimuli within normal and abnormal environments. These stimuli can introduce density gradients in the explosive charge but the effect of such gradients is not well-understood. The goal of this work is to develop a comprehensive methodology to model the evolution of density gradients under various thermal conditions and then to predict shock initiation and detonation behavior of the gradated sample. PBX 9502 was chosen for these scoping studies because it has a large coefficient of thermal expansion and because the available shock initiation data indicates that initial temperature strongly affects the shock-to-detonation transition (SDT). The multiphysics code COMSOL was first used to model induced density gradients in a slab of PBX 9502 with asymmetric heating. Then, the gradated sample was constructed in the Pagosa hydrocode to simulate steel flyer impact scenarios, with SDT and detonation propagation following a calibrated Scaled Uniform Reactive Front (SURF) burn model. Furthermore, the gradated material exhibited very different initiation and corner turning behavior compared to homogeneous materials, and the orientation of the gradation was also found to have some effect. The methodology used here can estimate thermomechanical response of the material for non-uniform thermal environments and compute the resulting detonation properties of the PBX.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Investigation of Anchor Link Magnetohydrodynamic Effects in the Toroidally Symmetric Lead-Lithium (TSLL) Blanket Concept

This paper presents a comprehensive study on the magnetohydrodynamic (MHD) flow in a slotted channel with a cylindrical anchor link, which is an imperative component of the toroidally symmetric lead-lithium (TSLL) blanket concept. Following the validation of the MHD solver implemented in COMSOL and a proposed subtraction approach to compute the anchor link pressure drop, the effects of computational domain size on the pressure drop and velocity distribution are examined. The results show that the pressure drop associated with the anchor link and the maximum flow velocity follow an asymptotic trend as the domain width increases, with wall-induced pressure drop being more dominant than that of the anchor link. The velocity distribution analysis revealed the formation of an internal boundary layer extending along the magnetic field direction, which is a unique feature of the investigated MHD flow. An estimation of the total MHD pressure drop associated with the array of anchor links under the TSLL blanket conditions suggests ~0.3 MPa, which is significantly lower than the recommended maximum allowable blanket pressure drop of 2 MPa. In conclusion, this work offers valuable insights into the anchor link–associated MHD phenomena and serves as a foundation for further development of the TSLL blanket concept for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Structured large-pore foams improve thermal performance of LiMIT-style liquid lithium PFC

As magnetically confined fusion devices improve, the conditions at the walls become increasingly intense. Plasma facing components (PFCs) must withstand these extreme heat and particle loads without damage or degradation. Liquid lithium PFCs are known to be quite resilient, and the presence of lithium also serves to improve plasma properties. The liquid metal infused trench (LiMIT) concept is an open surface liquid lithium PFC design that has been tested extensively at the University of Illinois and in fusion devices around the world. LiMIT utilizes thermoelectric magnetohydrodynamics (TEMHD) to passively drive liquid lithium flow. This work demonstrates an extension of the LiMIT trench geometry to three dimensions. Additively manufactured large pore metallic foams maintain TEMHD drive while drastically improving heat flux handling and resistance to lithium dryout, a phenomenon where locally high TEMHD forces depresses the lithium level and exposes underlying solid structure. COMSOL multiphysics modeling of the system yields insight into the forces at play in dryout development, and shows the 3D structures can eliminate dryout. Here, low heat proof-of-concept experimental testing of the system matches computational results, and high heat flux electron beam tests more than double the proven operational range of a LiMIT-style PFC, to 6.8 MW m –2 , with no indications of dryout or impending damage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

On the Inductance of a Single-Channel Spark-Gap Switch for Pulsed Power Applications

In this paper, we derive an alternative expression for calculating the inductance of a single-channel spark-gap switch. We believe this new inductance expression provides greater clarity compared to other approaches. Interestingly, both the Rosa formula and the coaxial inductance formula, which are commonly used in such calculations, emerge as special cases of our new expression. To validate the proposed formula, we conducted COMSOL simulations that show excellent agreement with our theoretical predictions.

COMSOL↗
DOE OSTI2020

A Cold Moderator For Sub-Thermal Neutron Flux Enhancement at the RPI-LINAC

The importance of high-quality thermal neutron scattering kernels coupled with a lack of accurate sub-thermal total cross section data prompted researchers at Rensselaer Polytechnic Institute (RPI) to design a cold moderation system to enhance sub-thermal neutron flux (defined as flux below 5 meV), for nuclear data measurements at the RPI Gaerttner LINAC facility. Using a series of neutronic and heat transfer simulations obtained with MCNP and COMSOL, an optimized final design was produced. The final design exceeded current target sub-thermal neutron producing capabilities by a factor of 8 in the sub-thermal region and is easily configurable with existing neutron producing targets.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗