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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Increasing Radon Daughter Plate-out onto Superconducting Qubits

Radon daughter plate-out may pose a fundamental limitation to realizing large-scale superconducting quantum computers. Radon daughters plated-out onto a superconducting qubit chip can decay, releasing energy into the substrate, producing phonons. Phonon induced quasiparticles generate bit flip errors and reduce coherence times. We seek to increase plate-out on chips to enable future study of the impact of Radon on coherence times. In this project, we utilize COMSOL Multiphysics simulations to finalize design for a device to increase radon plate-out.

Kay, Dax [Northwestern U.]↗
DOE OSTI2020

Iodine Effects on Medical Isotope Producer Delay & Guard Bed Temperature and Performance in Support of NA-23 Program

In support of NA-23’s objective to accelerate the establishment of a reliable commercial United States (U.S.) domestic supply of 99Mo produced without the use of highly enriched uranium (HEU), the Pacific Northwest National Laboratory (PNNL) was tasked with modeling the performance of activated carbon delay and guard beds that may be typical for non-HEU 99Mo production processes. The specific tasks accomplished by PNNL included calculations of long-term effects of iodine on delay bed performance and development of a finite element simulation to better understand delay bed performance under various scenarios specifically with radioxenon and radioiodine decay heating of the delay beds. The main issue of interest in this study is understanding the effect of iodine on the abatement performance of the beds if iodine is not sufficiently removed from gases prior to reaching the delay beds. Finite element simulations of delay bed performance were developed using COMSOL Multiphysics® Software. These simulations are based on partial differential equations that describe xenon and iodine adsorption, heat transfer and decay of radioxenon and radioiodine isotopes. Heat transfer in the activated carbon beds is simulated with a heat source tied to the decay of radioactive xenon and iodine species. These simulations were used to investigate delay bed temperature and hold up time using radioxenon and radioiodine activities and flowrates expected from a non-HEU production process. A bounding case scenario was simulated, where all the iodine was released into the delay beds due to insufficient guarding of the delay beds. In addition, preliminary investigations of the guard beds were conducted.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Aerodynamic Characterization of 3D Scanned Wind Turbine Blades Using Experimental and Computational Methods

This study presents an aerodynamic characterization of 3D scanned wind turbine blades using both experimental and computational methods. The research was conducted by Gulf Wind Technology and Sandia National Laboratories. The primary objective was to investigate the aerodynamic impacts of leading-edge manufacturing defects on wind turbine blades. The study utilized the Stratasys NEO 800 3D Printer for high-precision manufacturing and the GWT Accelerator Wind Tunnel for experimental testing. Computational simulations were performed using COMSOL Multiphysics to model the wind tunnel and analyze flow characteristics and OpenFOAM to study the aerodynamic impacts of leading-edge defects. OpenFAST was used to estimate how these defects can lead to revenue losses for wind farm operators as high as 6%. The results demonstrated significant aerodynamic performance variations due to defects, with detailed analysis provided through wind tunnel and CFD data. The findings contribute to the understanding of defect impacts on wind turbine blade performance and offer insights for future design improvements.

17 WIND ENERGY↗
DOE OSTI2025

Thermal Hydraulic Experimental Test Article: Second Year of Testing with Secondary Sodium System (Fiscal Year 2025 Final Report)

The Thermal Hydraulic Experimental Test Article (THETA) is currently installed in the Mechanisms Engineering Test Loop (METL) 28” test vessel #4. Both the primary and secondary sodium systems remain online to facilitate continued testing. This fiscal year, work was performed using a COMSOL Multiphysics magnetohydrodynamic model to characterize flow more accurately in the secondary electromagnetic flowmeters. Experimental campaigns were then performed to study the thermal hydraulic differences between sodium and water as a surrogate fluid in the THETA geometry as well as a study to better characterize and understand temperature oscillations that exist at the outlet of the core to the hot pool. A peer-reviewed article was published in the ASME Journal of Nuclear Engineering and Radiation Science detailing the THETA facility and providing an overview of a test that was performed with the primary and secondary system online [1]. Work continues to develop a database to house experimental THETA data to better facilitate collaboration with industry and laboratory partners for their use of the data for code benchmarking/validation. THETA remains fully operational and is positioned for continued testing in fiscal year 2026.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Enhancing Even Gas Distribution in Porous Media with Radial Flow

Packed beds, such as those used in thermal energy storage (TES) systems, typically use flow from one end to the other. This axial flow configuration leads to inefficiencies due to dispersion effects and high pressure drop. Instead of axial flow, this work proposes radial flow in packed beds, where a central tube provides flow that transports fluid from the center to the bed’s wall. Radial flow could be a promising solution to increase the efficiency of charging/discharging processes in TES systems. For instance, studies have been conducted on the thermal behavior during the charging process using radial flow and axial flow, and it was found that radial flow is better than axial flow in terms of thermal performance, where more energy can be stored in the storage tank during the charging process. However, the design of the radial pipe should be optimized to improve even flow distribution. In this work, numerical investigations on different designs were analyzed to enhance even flow gas distribution in the porous media when using this radial technique. This work shows the impact of different parameters on the even flow distribution into the packed bed for two different designs: 1) one radial tube at the center to provide the radial flow along with an annular tube at the wall to receive the flow, and 2) one radial tube at the center and four radial tubes at the wall. Air was used as fluid and 6 mm alumina beads as packing materials. Computational fluid dynamics (CFD) models in COMSOL Multiphysics were used to simulate the behavior of air flow through the piping and packed bed. It was found that the flow into the bed from the radial tube can be affected by different parameters: space between the holes, size/diameter of the holes, number of segments/zones, number of the holes in each zone and the length of each zone.

Beck, David M.↗
DOE OSTI2021

Improving the W Coating Uniformity by a COMSOL Model-Based CVD Parameter Study for Denser W f /W Composites

Tungsten (W) has the unique combination of excellent thermal properties, low sputter yield, low hydrogen retention, and acceptable activation. Therefore, W is presently the main candidate for the first wall and armor material for future fusion devices. However, its intrinsic brittleness and its embrittlement during operation bears the risk of a sudden and catastrophic component failure. As a countermeasure, tungsten fiber-reinforced tungsten (W f /W) composites exhibiting extrinsic toughening are being developed. A possible W f /W production route is chemical vapor deposition (CVD) by reducing WF 6 with H 2 on heated W fabrics. The challenge here is that the growing CVD-W can seal gaseous domains leading to strength reducing pores. In previous work, CVD models for W f /W synthesis were developed with COMSOL Multiphysics and validated experimentally. In the present article, these models were applied to conduct a parameter study to optimize the coating uniformity, the relative density, the WF 6 demand, and the process time. A low temperature and a low total pressure increase the process time, but in return lead to very uniform W layers at the micro and macro scales and thus to an optimized relative density of the W f /W composite. High H 2 and low WF 6 gas flow rates lead to a slightly shorter process time and an improved coating uniformity as long as WF 6 is not depleted, which can be avoided by applying the presented reactor model.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A Threshold Helium Leakage Detection Switch with Ultra Low Power Operation

Detecting helium leakage is important in many applications, such as in dry cask nuclear waste storage systems. This work develops a helium detection system based on the relative permittivity (dielectric constant) difference between air and helium. This difference changes the status of an electrostatic microelectromechanical system (MEMS) switch. The switch is a capacitive-based device and requires a very negligible amount of power. Exciting the switch’s electrical resonance enhances the MEMS switch sensitivity to detect low helium concentration. This work simulates two different MEMS switch configurations: a cantilever-based MEMS modeled as a single-degree-freedom model and a clamped-clamped beam MEMS molded using the COMSOL Multiphysics finite-element software. While both configurations demonstrate the switch’s simple operation concept, the clamped-clamped beam was selected for detailed parametric characterization due to its comprehensive modeling approach. The beam detects at least 5% helium concentration levels when excited at 3.8 MHz, near electrical resonance. The switch performance decreases at lower excitation frequencies or increases the circuit resistance. The MEMS sensor detection level was relatively immune to beam thickness and parasitic capacitance changes. However, higher parasitic capacitance increases the switch’s susceptibility to errors, fluctuations, and uncertainties.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Dielectric Resonator Design for Low Power and Low Temperature Microwave Plasma

Waveguide-based microwave plasmas generally operate at high temperatures (2000 - 6000K)[1], making it difficult to directly interface solid materials with the plasma without significant thermal damage. Dielectric microwave resonators (DMRs), long studied for wave-based manipulation of electromagnetic radiation for telecom and optics, can focus radiation to extremely small mode volumes, creating intense localized fields with low-power input.[2] This phenomenon can be used for applications ranging from efficient plasma electronics to near-ambient plasma-materials interactions. Such DMR-based plasmas have been demonstrated a handful of times in the literature, but the majority of research towards this utilize the lowest frequency resonance mode.[3], [4], [5] By carefully controlling the geometry of cylindrical resonators, a variety of electromagnetic modes can be excited. In this work, COMSOL Multiphysics simulations are used to study the electric field enhancement and absorption properties of CaTiO3 DMRs as a function of geometry and excitation frequency. Whereas previous studies have utilized the HEM111 resonance frequency to drive low power plasma excitation, we find that higher order resonance frequencies are more effective at field enhancement and result in less power loss within the dielectric material, hence less wasted heating. The effectiveness of these modes is also geometry dependent and can be computationally optimized for plasma generation. Complementing these computational efforts, we demonstrate a new closed-system reactor design built in a WR-650 waveguide and experimentally demonstrate the formation of atmospheric argon microwave plasma using < 30 W input power on DMR dimers. We observe a shifting resonance frequency as the DMRs heat in response to microwave excitation and develop a Python-based lock-in mechanism to effectively track the DMR resonance over time, leading to stable plasma operation. We use infrared thermal imaging to monitor the temperature of the DMR dimers and surrounding quartz chamber, demonstrating thermal temperatures < 60 degreesC. Finally, we utilize optical emission spectroscopy (OES) to probe the plasma properties as a function of the resonance mode.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗
DOE OSTI2022

Modeling Vapor Transport Deposition of Metal-halide Perovskite Thin Films for Photovoltaic and Optoelectronic Devices

Over the past decade, metal halide perovskites (MHPs) have emerged as a promising materials platform for high-efficiency solar cells and low-cost optoelectronics. However, there are challenges that frustrate the large-scale manufacturing of MHP devices, including difficulty in controlling film composition, interface formation and their device instability under ambient conditions. Vapor processing offers an attractive path to manufacturability, while also opening the door to new opportunities in device design that could favorably impact ultrahigh efficiency tandem solar cells or overall stability. Emphasis here is on the use of an alternate processing methodology, vapor transport deposition (VTD), to deposit the MHP layer. In experimental work, we found that our deposition system had significant run-to-run variations in film thickness and composition. To understand and resolve these issues, we use COMSOL Multiphysics software to model the precursor deposition rate numerically, trying to guide the broad range of parameters in the system. The result can be matched with experimental data and provides insight into system hydrodynamics and molar transport effects. Finally, we show that chamber pressure is a key factor to scale up the VTD technique.

Hsu, Wan-Ju↗
DOE OSTI2020

Development of the Intelligent, Preventive Infrared (IR) Inspection System Housed in Hybrid Robotic Platforms

Robots and robotic systems that are designed for inspection, environmental study, and health and safety aid are becoming an increasing necessity. However, there is a number of challenges that accompany robots that are designed for these specific applications. These challenges include: navigating compact, enclosed spaces, travelling over multiple terrains and large obstacles, using the proper sensing and detection methods to assess an environment, and the use of lightweight and durable materials. The robotic platforms currently in development look at all of these challenges and attempt to overcome them. These designs specific use of hybrid robotic platforms, or platform that utilizes soft and rigid materials, allows for a more flexible platform and makes environments more navigable. To further improve the navigation of the platforms and environmental assessment, a novel infrared detection system is housed in the platforms to create a robot that can be used for the applications listed above and more. Objectives: Further develop two types of robotic platforms that utilize additive manufacturing, soft materials, and rigid materials. Continue the development of an intelligent and inhibitory infrared detection system based on an artificial intelligence (AI) algorithm. Continued study and fabrication of active soft materials designed for both sensing and actuation in hybrid robotic systems. Improve additive manufacturing fabrication to design rigid and semi-rigid components for hybrid robotic platforms. Transformable Wheel Robotic Platform: The chassis, wheels, tires and inspection system housing use different additive manufacturing techniques for fabrication. Continued work with additive manufacturing has lead to studies in metal-based printing and modular design and manufacturing. The new platform design with integrated electrical component printed. This will allow integration of the sensor housing onto the platform. Electrical components are being tested for battery life and performance. To improve this performance, such as integration of Lithium Polymer (LiPo) batteries. Snake Robotic Platform: The main focus of the development has centered around liquid-based soft actuators. That act on the principles of electrostatic and hydraulic actuation. A liquid dielectric sits between two compliant electrodes, contained by a flexible polymer shell. The electrodes and film gradually collapse toward each other from one corner of the electrode to the other. When the electrodes and film close together, a majority of the fluid is pushed into the area not covered by an electrode. A thin layer of the liquid dielectric remains between the electrode. The actuators will be stacked to cause large displacement, and move the linkages. The chassis of this platform uses purely additively manufactured linkages. Intelligent, Preventive IR Inspection System: Development of the AI for the system has lead to using a Scikit-Learn which assists in creating predictive models based on Regression, clustering, classification etc. To improve the infrared thermometry for low emissivity sources, work on the fabrication of a tandem photoconductive infrared thermometer was a main focus. Distance-Voltage-Temperature response data has been collected in the range of 7 cm - 100 cm and 200-400 deg. C. Modifications were made to the existing test bench to have a better control over the data. Automated data collection was realized using a Python code and an Arduino controlled stepper motor to increase the sample rate. A protective enclosure has been built around the setup to minimize the effect of the environment on the measurements. To better predict temperature, different AI models are being optimized. The regression model, LARS showed a high accuracy but had convergence issues and only works for the current test set-up. Results: The Transformable Wheel Robot has developed into a more flexible and modular platform. With the improvements to the current work, effort on the tire or soft gripper design been a large focus. The soft grippers will be interchange able to allow for increased performance in identified terrain types. The development of the liquid-based actuators allows the snake robotic platform to achieve the goals of being flexible and able to navigate confined spaces. However, there is room for improvement. Optimization work is currently being done in COMSOL Multiphysics. With the current IR system set-up the LARS model perfectly predicts the data; however, considering the mobility aspect of the project other models will allow for an optimized system. Testing of other model types is currently being done.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Increasing Radon Daughter Plate-out onto Superconducting Qubits

Radon daughter plate-out may pose a fundamental limitation to realizing large-scale superconducting quantum computers. Radon daughters plated-out on a superconducting qubit can decay, releasing energy into the substrate, producing phonons. Phonon-induced quasiparticles then generate bit-flip errors and reduce coherence times. We seek to increase plate-out on chips to enable future study of the impact of radon on coherence times. In this project, we utilize COMSOL Multiphysics simulations to visualize electric field configuration and model ion transport to finalize design of a device to increase radon plate-out on qubit chips.

Kay, Dax [Fermilab; Northwestern U.]↗
DOE OSTI2024

Thermal performance of phase change materials with anisotropic carbon fiber inserts

In thermal energy storage systems, phase change materials (PCMs) are widely used for thermal energy management. Most PCMs have low thermal conductivities, which limits the heat transfer rate within PCM and thus makes the phase-changing process very slow. However, thermal conductivities of PCMs can be altered by inserting targeted additives. We hypothesize that the size and shape of these additive inserts play a key role in thermal management efficiency. To this end, the impacts of carbon fiber (CF) inserts on the phase change behavior and consequent heat transfer efficiencies of inorganic and organic PCMs were investigated using experiments and simulations. Long, anisotropic CFs with high thermal conductivities formed continuous fast heat flux tunnels inside PCMs to enhance the heat transfer. Such CFs could extend the phase change fronts from the limited container-shaped interface to the larger surface of numerous CF inserts inside the PCM. These special CF inserts work with a new heat transfer mechanism, different from conventional small additives or long isotropic CF inserts. The thermal energy release rate increased by 2.5 times with 1 wt.% anisotropic CF inserts in inorganic PCM. However, CF inserts in liquid organic PCM hindered the natural convection and compromised the improved heat conduction. In conclusion, the lab-scale multiphysics simulations support these experimental observations and indicate that CF inserts have potential to enhance heat transfer in inorganic PCMs, but they are less effective in organic PCMs.

25 ENERGY STORAGE↗
DOE OSTI2025

Modeling and Characterizing the electron backscatter in a cylindrical anode-based distributed X-ray source

Upcoming advancements in computed tomography architectures warrants the investigation of new X-ray source designs and the impacts that electron backscatter can have on these designs. One such design being investigated is a distributed, cylindrical anode-based X-ray source. For such a distributed X-ray source, we developed a modeling pipeline for simulating electron optics and transport to characterize the quality of the primary X-ray beam and the electron backscatter behavior. We report our results on the energy distributions of the bremsstrahlung spectra; electron backscatter ratio; and spatial, temporal, and energy distributions of backscattered electrons that return to the anode.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
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 OSTI2022

Spatiotemporally resolved measurements of electric field around a piezoelectric transformer using electric-field induced second harmonic (E-FISH) generation

When a piezoelectric transformer (PT) is actuated at its second harmonic frequency by a low input voltage, the generated electric field at the distal end can be sufficient to breakdown the surrounding gas, making them attractive power sources for non-equilibrium plasma generation. Understanding the potential and electric fields produced in the surrounding medium by the PT is important for effectively designing and using PT plasma devices. In this work, the spatiotemporally resolved characteristics of the electric field generated by a PT operating in open air have been investigated using the femtosecond electric field-induced second harmonic generation (E-FISH) method. Here, electric field components were determined by simultaneously conducting E-FISH measurements with the incident laser polarized in two orthogonal directions relative to the PT crystal. Results of this work demonstrate the spatial distribution of electric field around the PT’s output distal end and how it evolves as a function of time. Notably, the strongest electric field appears on the face of the PT’s distal surface, near the top and bottom edges and decreases by approximately 70% over 3 mm. The time delay between the PT’s input voltage and measured electric field indicates that there is an about 0.45 phase difference between the PT’s input voltage and output signal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Modelling the brittle failure of graphite induced by the controlled impact of runaway electrons in DIII-D

The thermo-mechanical response of an ATJ graphite sample to controlled runaway electron (RE) dissipation, realized in DIII-D, is modelled with a novel work-flow that features the RE orbit code KORC, the Monte Carlo particle transport code Geant4 and the finite element multiphysics software COMSOL. KORC provides the RE striking positions and momenta, Geant4 calculates the volumetric energy deposition and COMSOL simulates the thermoelastic response. Brittle failure is predicted according to the maximum normal stress criterion, which is suitable for ATJ graphite owing to its linear elastic behavior up to fracture and its isotropic mechanical properties. Measurements of the conducted energy, damage topology, explosion timing and blown-off material volume, impose a number of empirical constraints that suffice to distinguish between different RE impact scenarios and to identify RE parameters which provide the best match to the observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗