Search NASA⌕ Search

SEARCH · Search NASA

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

Design workflow of a symmetric traveling wave antenna for fast ion production on DD tokamaks

Initial computational plasma physics scoping and a finite element method antenna modeling design workflow for a symmetric center-fed high-field side high harmonic fast wave traveling wave array (TWA) antenna are reported here. The TWA is designed to generate a test population of fast deuterium ions in an existing D–D tokamak by heating neutral beam deuterium ions, accelerating them from 80 keV to several hundred keV. The resulting fast particles are tailored to mimic key reactor energetic particle parameters with regards to exciting Alfven eigenmode instabilities, allowing for a D–D tokamak like DIII-D or ASDEX-U to replicate reactor-relevant conditions experimentally. Initial scenario scoping for high single-pass absorption as well as good preferential fast ion damping relative to electron damping was completed using the ray-tracing/Fokker–Planck codes GENRAY and CQL3D. Python RF network analysis packages were used to create a custom TWA optimization tool to inform a COMSOL flat antenna design, and Petra-M was used to study cold plasma effects. The TWA produced by this workflow has several novel features when compared to previous TWA studies, including symmetric center feeding, and passive end straps for image current cancellation for reduced impurity production. We show here that the antenna design workflow can readily produce TWA antennas optimized for reflection coefficient, image current cancellation, and launched power spectrum shape; and that a population of fast ions can be generated in the correct region of parameter space, warranting future more detailed studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Integrated modeling of RF-induced tungsten erosion at ICRH antenna structures in the WEST tokamak *

This paper introduces STRIPE (Simulated Transport of RF Impurity Production and Emission), an advanced modeling framework developed to analyze material erosion and the global transport of eroded impurities originating from radio-frequency (RF) antenna structures in magnetic confinement fusion devices. STRIPE integrates multiple physics modules: SolEdge3x for scrape-off-layer plasma profiles, COMSOL for 3D RF rectified sheath potentials, RustBCA for erosion yields and surface interactions, and global impurity transport for 3D ion energy-angle distributions and impurity transport. The framework is applied to an ion cyclotron RF-heated L-mode discharge (#57877) in the WEST tokamak, where it predicts a thirty-fold increase in gross tungsten erosion at antenna limiters during the transition from ohmic to ICRH operation. Additionally, under ICRH conditions, a tenfold enhancement in erosion is observed when comparing RF sheath effects to purely thermal sheath conditions. High-charge-state oxygen ions ($\mathrm{O}$ 6+ and above) are identified as the dominant contributors to tungsten sputtering. To validate the model, a synthetic diagnostic tool based on inverse photon efficiency (S/XB coefficients) from the ColRadPy collisional-radiative model enables direct comparison with spectroscopic measurements. Model predictions using a plasma composition of 1% oxygen and 99% deuterium show good agreement with observed W − I (400.9 nm) emission for discharge #57877, supporting the accuracy of the STRIPE framework. This study focuses specifically on gross erosion calculations to demonstrate STRIPE’s capabilities. Future extensions of this work will incorporate net erosion, re-deposition, self-sputtering effects, and whole-device modeling of sputtered tungsten impurity transport. STRIPE is also being applied to other RF-heated linear and toroidal devices, offering valuable insights for antenna design, impurity control, and performance optimization in next-generation fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Optimizing the Internal Cell Structure of the Radon Monitor Based on Electrostatic Collection Method

Abstract As Rn-222 decays, an alpha particle is emitted and the residual polonium nucleus recoils in the opposite direction. At the end of the recoil path, 88% of the polonium atoms have a positive charge and 12% are neutral. The electric potential distribution in the 60 ml hemispherical internal cell of the radon monitor based on electrostatic collection is studied for reducing the combined probability of the positively charged Po-218 and the OH− produced by the ionization of water vapour in the air. The COMSOL software is used to simulate the electric potential distribution in the internal cell of the radon monitor based on the electrostatic collection method at 27°C, a pressure of 0.1 Mpa. For improving the collection efficiency of Po-218 ions, the average collection time along vertical and oblique lines is calculated when the upper surface of the internal cell is plastic, uncharged metal and charged metal, respectively. Assuming that the gas in the internal cell is uniformly distributed, the results show that if the upper surface of the hemispherical internal cell is plastic or uncharged metal, the electric potential formed in the internal cell is more uniform, and it is beneficial to reduce the total collection time of the positively charged Po-218 ions, thereby improving collection efficiency. The simulation results can be used as an effective reference for optimizing the design of the internal cell structure of the radon monitor based on electrostatic collection method.

Fan, Zhongkai (ORCID:0000000181247459)↗
DOE OSTI2021

Effect of Mechanical Deformation on the Dielectric Electric Field in Dynamic Umbilical Cables

Wave energy converter (WEC) systems have unique umbilical cable design requirements compared to conventional marine designs. Due to their dynamic nature, WEC umbilical cables are required to handle increased loads and larger motions over longer deployment scenarios. However, the life expectancy of dynamic umbilical cables is predominantly analyzed for mechanical fatigue, with little consideration given to the effect of mechanical stress on the response of the dielectric as it has been previously assumed that the non-polymeric components will fail prior to the polymeric components. In this work, we explore the effects of mechanical bending on the concentration of the electric field in the insulation of a conceptual three-core dynamic umbilical cable to determine the limiting requirements for life expectancy. A conceptual industrial schematic of a medium voltage dynamic umbilical cable was used as the basis for the WEC system under study. Simulations on the dynamic umbilical cable were performed using COMSOL Multiphysics®. Mechanical loads up to the calculated limit (2.64 m minimum bend radius or 388 kN transverse load) were applied to the dynamic umbilical in a step-wise fashion with one end of the umbilical cable fixed. One-way coupling in three-dimensions was conducted by first determining the engineering strains, then using the engineering strains as the initial condition for determination of the Green-Lagrange strains, and lastly calculating the electric field in the deformed dielectric. Based upon the input load, deformation of the dielectric at the minimum bend radius produced an increase in the electric field of more than 12% when compared to the undeformed umbilical cable, which may lead to dielectric breakdown of the insulation prior to mechanical failure of the umbilical cable.

Spencer, Mychal P.↗
DOE OSTI2022

3D Frequency Domain Reflectometry Digital Twin of an Electrical Cable: A First Glance

Electrical cables within nuclear power plants (NPPs) are critical components required for power, control, and instrumentation systems which may be exposed to stressors, such as elevated temperatures and gamma radiation. Such stressors can lead to a reduction in the remaining useful life of electrical cables, jeopardizing the safety of NPP systems. To evaluate the effect of stressors on the degradation of electrical cables, electrical reflectometry methods are commonly employed. Frequency domain reflectometry (FDR) is a non-destructive electrical reflectometry method that uses transmission line theory to detect degradation or impedance changes within electrical cables. However, in most cases FDR is only applied to de-energized cables, limiting the application in NPPs as the cable system must be taken offline. In this work, we explore the development of an FDR digital twin to predict the degradation of an electrical cable exposed to elevated temperature, which is expected to reduce the need for offline FDR. A 3-conductor low-voltage electrical cable was selected for evaluation of the digital twin. The fully three-dimensional digital twin was developed in COMSOL using the RF module. A cable length of 30-m and frequency bandwidth of 400 MHz was selected to mimic real-world application of FDR. Over a 1-m region, the permittivity of the insulation was varied by up to 20% to model thermal degradation. The results demonstrate accurate detection of the insulation damage region, supporting further investigation of the FDR digital twin using real-world data and machine learning for predictive damage estimation or remaining lifetime.

Spencer, Mychal P.↗
DOE OSTI2021

Investigating the Impact of Thickness, Calendering and Channel Structures of Printed Electrodes on the Energy Density of LIBs - 3D Simulation and Validation

Current lithium ion batteries (LIBs) are expensive and bulky, limited by relatively low charging rates. To increase the rate of charging and reduce weight, thin electrodes with high energy density are required. The increase in energy density can be achieved by several techniques including boosting electrolyte transport, high loading/utilization of active material, employing high conductive electrolytes and electrodes with advanced architectures, and increasing cell temperature. In this paper, a 3D physics-based electrochemical model of LIBs is developed in COMSOL simulation software for different thickness, calendering steps as well as channel structures (conical, cylindrical) to optimize the electrode design and in turn maximize volumetric energy density. The simulation results demonstrated that calendering the electrodes with high initial porosity increases the volumetric energy density of the cell. In addition, cylindrical channel structures with relatively lower edge-to-edge distance also results in increased volumetric energy density. The simulation results of the 3D model was validated by comparing it with experimental results.

improving volumetric energy density↗
DOE OSTI2024

Thermal Analysis of a Polypropylene Capacitor for Resonant Tuning Networks in WPT Applications

This study presents the thermal analysis of polypropylene capacitors operating at high frequencies and high currents used in the resonant tuning networks for electric vehicle (EV) wireless charging systems. The thermal analysis and performance results of the CELEM polypropylene capacitors are presented at different currents for a polyphase wireless power transfer (WPT) system while the study can be expanded to other capacitors by different manufacturers. Thermal equivalent circuit of the capacitor is derived and required cold plate size and design calculations are conducted by modeling the polypropene capacitor in MATLAB. The polypropene capacitor's finite element analysis (FEA) model is developed using COMSOL Multiphysics simulation software, and the thermal characteristics and analysis of the capacitor are presented in this paper.

Aktas, Ahmet↗
DOE OSTI2021

Simulation App for Time-Resolved Photoluminescence in Thin-Film Solar Cells

We develop a COMSOL Multiphysics-based model and application, or "app", for time-resolved photoluminescence (TRPL) measurements in thin-film solar cells. TRPL is commonly used to measure recombination in thin-film solar cells. However, simulations show that decay times extracted from TRPL data are not, in general, equal to the Shockley-Read-Hall lifetime parameter set in the model. Here, we present a model and app, using graded CdSe y Te 1-y solar cells as an example, that can be used to study impacts of a range of potentially important material and measurement parameters on TRPL results.

minority-carrier lifetime↗
DOE OSTI2022

Simulation for the fault current limiting operation of REBCO CORC superconducting cables with different core materials

With high power density and low loss, CORC superconducting cables composed of RE- Ba 2 Cu 3 O 7–δ (REBCO) coated conductors are of great interest for power transmission applications. They can also effectively be designed to serve as fault current limiters, thanks to their sharp superconducting-to-normal transition. Coupled electromagnetic-thermal finite element (FE) simulations implemented in COMSOL Multiphysics were developed and validated against the published results for the fault current limiting performance of two CORC cables of REBCO coasted conductors wound on either copper or stainless steel cores. For improved accuracy, temperature dependencies of electrical and thermal properties of all component materials were considered in the simulations. The simulations were performed in the cross-sections of the cables for every individual layer of each REBCO tape to deliver comprehensive understanding of evolution of the current distribution and temperature rise in those layers. Furthermore, the studies can suggest approaches to optimize cable design for fault current limiting (FCL) applications by assessing the role of individual components. In the paper, possible computational errors, challenges and improvements are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Cryogenic System Preliminary Design for a 0.5m-Long, Conduction-Cooled Nb3Sn Undulator Magnet Prototype

There are several NbTi superconducting undulator (SCU) magnets currently in operation at the Advanced Photon Source (APS) at Argonne National Laboratory (ANL). The development on Nb3Sn-based superconducting undulator magnets at APS is underway due to the potential to further enhance the performance of the SCUs. Superconducting undulator magnets need to keep temperature gradients minimized in order to retain thermal and operating current margin. We have designed the thermal links for efficient heat conduction using 3D finite element analysis (FEA) simulation in COMSOL Multiphysics software, which was later used for the calculation of the temperature distribution across a 0.5 m long, conduction-cooled, Nb3Sn undulator magnet prototype that includes both conductive heat transfer and radiative heating components. We have modelled the evolution of the thermal properties of the magnet winding as well as other cold parts during cool-down from ambient temperature with the operation of a SHI RDK-415D cryocooler to examine the estimated time that is needed for the cooldown as well as the baseline temperature we could achieve. A key result was that a maximum coil delta T = 0.11 K along the designed thermal links and a temperature range of 3.42 K-3.53 K for the winding were predicted at steady state. A pair of G10 support rods for carrying the magnet as well as a pair of current leads made of metals were also designed and optimized for the testing system, targeting for an operation temperature of 4 K and a coil current of 950 A DC.

Zhang, Danlu↗
DOE OSTI2023

Additive Manufacturing of Coreless Flyback Transformers Using Aerosol Jet Printing and Electrochemical Deposition

An additive manufacturing approach combining aerosol jet printing (AJP) and electrodeposition opens a new pathway to the production of lightweight coreless flyback transformer devices for power electronics. AJP of seed layers with resolution on the order of 30 μm is combined with electrodeposition of Cu and Ni for decreased resistance. This combined approach addresses known shortcomings of AJP and electrodeposition. Nanoparticle inks used in AJP of metals have low conductivity versus bulk materials due to their high grain boundary resistance. There is a lack of readily available high-resolution patterning techniques for electrodeposition outside of expensive clean-room-based lithography techniques. Combining these two techniques enables the patterning of high-resolution, high-conductivity components. Here in this manuscript, we report on the construction of coreless flyback transformers consisting of two-layer primary and two-layer secondary spiral inductors separated by layers of a printed UV-curable dielectric. An input voltage of 17 V at 400 kHz was amplified to an output of 1250 V corresponding to a gain of 73.5. COMSOL modeling at the individual inductor level and at the transformer level was used to compare expected inductance, equivalent series resistance (ESR), and coupling with experimentally measured values.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multiphysics Analysis of Li Cooled Divertor Substrate During Loss of Coolant Accident (LOCA)

In the ongoing study of potential designs for liquid metal (LM) plasma-facing components (PFCs), so-called “slow” and “fast” Li flow divertor concepts are under investigation. In the previous studies on design and analysis of the slow Li flow divertor and comparison with the fast Li flow divertor, the magnetohydrodynamics (MHD)/heat transfer effects of the Li flowing inside the substrate as a second coolant were comprehensively investigated under the normal steady-state operation conditions. Here, in the present study, the multiphysics analysis is extended to the unsteady abnormal divertor scenario where the Li layer on top of the substrate does not provide full coverage or even totally disappears for a certain period of time, so that the substrate becomes directly exposed to the incident high plasma heat flux. Such an unwanted event may happen regardless of the concept of the divertor and is worth detailed investigations, typically referred to as a loss of coolant accident (LOCA). To address this situation, a simplified scoping analysis is conducted first in 2-D, and then an integrated 3-D modeling is performed using a time-dependent multiphysics model in COMSOL Multiphysics that integrates LM MHD, heat transfer, and solid mechanics. The main goal is to evaluate conditions under which the major material limits, such as the maximum allowable temperature, stress, and displacement of the substrate, can still be met. It was shown that the maximum time over which the substrate of RAFM steel can retain structural integrity during the LOCA is around 0.2 ∼ 0.3 s. Any divertor concept that utilizes RAFM steel as a substrate material and liquid Li as a second coolant should take such a permitted time into consideration.

divertor↗
DOE OSTI2023

Numerical Study of Solar Receiver Tube With Modified Surface Roughness for Enhanced and Selective Absorptivity in Concentrated Solar Power Tower

Abstract Concentrated solar power (CSP) is a reliable renewable energy source that is progressively lowering its cost of energy. However, the heat loss due to reflected and emitted radiation hinders the maximum achievable thermal efficiency for solar receiver tubes on the solar tower. Current solar selective coatings cannot withstand the high temperatures that come with state-of-the-art CSP towers often needing to be recoated soon after initial operation. We intend to use Inconel 718 with different additive manufacturing (AM) practices to construct surfaces that allow for more light-trapping to occur. By adjusting printing parameters, we can tailor a surface to allow for more absorption while diminishing emitted radiation heat loss. By using COMSOL Multiphysics, we can generate these theoretical surfaces to emulate a printed surface, and using the coupled Multiphysics we can simulate how the surface dictates radiation properties. Our results show that by having a rougher surface we can enhance the absorptivity of Inconel 718 (IN718) by 38.8%. We expect this work to transform how solar absorber tubes are manufactured without using selective coatings and supplement the US Department of Energy (DOE) 2030 SunShot Initiative.

14 SOLAR ENERGY↗
DOE OSTI2024

Smart Adaptive Structures for an Ocean Wave Energy Converter

Ocean wave energy converters face significant challenges including cost-effectiveness, minimizing maintenance requirements, and withstanding extreme conditions. However, by utilizing smart materials, these converters could overcome these challenges. Such energy harvesters could use dielectric elastomer generators to convert ocean wave energy into electricity through their dynamic straining. Conversely, by applying electricity to these generators, they become actuators - dielectric elastomer actuators - thereby enabling them to alter their stiffness and adapt to the ever-changing ocean energy environments. Such active adaptation could enhance the converter's ability to: reach resonance with ocean waves and protect itself from dangerous waves. This study utilizes numerical analyses through the COMSOL software framework to evaluate the potential energy that could be harvested by a conceptual ocean wave energy converter based upon dielectric elastomer generators/actuators. The converter is composed of an external hull (that is a hollow cylinder), an inertial mass (that is a hollow cylinder and concentric with the hull), and 'spokes' - made of dielectric elastomer generators/actuators - that connect the hull to the inertial mass. Results of the numerical analyses include those outcomes arising from the conceptual converter being simulated via a sinusoidal motion analogous to ocean waves. That motion, therefore, causes relative motion between the converter's hull and inertial mass thereby dynamically stretching the corresponding connecting elastomers. The stretching of the elastomers enables them to 'gain elastic strain energy' and is, therefore, considered to be the theoretical limit of possible electrical energy conversion for the dielectric elastomer generator/actuator spokes. Additionally, the elastomer material properties of the spokes were altered to simulate the actuation of those same elastomers; with overall strain energy being subsequently investigated. Ultimately this is a preliminary study exploring the ability of such smart materials - electricity-generating and self-actuating elastomers - to actively adapt an ocean wave energy converter's structure to address and overcome the aforementioned challenges.

active materials↗
DOE OSTI2020

O · , H · , and · OH radical etching probability of polystyrene obtained for a radio frequency driven atmospheric pressure plasma jet

Atmospheric pressure plasma jets have great potential for the surface modification of polymers. In this work, the authors report on polystyrene etching by a radio frequency driven atmospheric pressure plasma jet with a focus on the role of H · , O · , and · OH radicals in this process. The absolute flux of H · , O · , and · OH radicals reaching the surface of the polymer was determined by a COMSOL MULTIPHYSICS reacting fluid dynamics model incorporating detailed transport phenomena in the boundary layer near the substrate. The simulated results of H · and · OH densities in the jet effluent were experimentally verified by two-photon absorption laser induced fluorescence and laser induced fluorescence, respectively. The carbon atom removal flux from the polystyrene surface was taken from previously reported measurements using the same plasma source. The authors show that the boundary layer effects in the interfacial region above the substrate can have a significant impact on the calculated etching probabilities. The reaction probability (β) has a significant uncertainty although a variation of 2 orders of magnitude in β leads to uncertainties of approximately 1 order of magnitude variation in the determined etching probability. The etching probability of polystyrene by · OH radicals was confirmed to be at least an order of magnitude larger than the polystyrene etching probability by O · radicals. The authors also confirmed the weak polystyrene etching probability by H · radicals. The model suggests that the presence of a 30 ppm O 2 impurity can lead to the production of · OH radicals in the far effluent of the Ar + 1% H 2 plasma jet close to the substrate at sufficient densities to enable effective etching.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Computational modeling of phononic pseudocrystal isolators

Methods for the efficient computational prediction of the performance of phononic pseudocrystals (structured materials capable of blocking extraordinary ranges of frequency) in COMSOL and other comparable finite element method codes are set forth. These methods include boundary conditions that make possible halving the size of the computational domain. Also included is an introduction of elastic energy density methods for assessing the extinction of elastic waves within the patterned region.

Swift, Stephen Hales↗
DOE OSTI2023

Direct Measurements of Ionic Transport Behavior of Dual-Layer Porous Electrodes

To improve power and cycling performance of lithium-ion batteries, dual-layer or porosity-gradient electrodes have been proposed. By engineering a higher porosity close to the separator, the intention is to improve ion transport where it is most impactful. In this research, MacMullin numbers of two dual-layer anode samples are tested using an impedance measurement technique developed previously. To characterize the microstructure of each layer independently, we developed an improved transmission-line model that accounts for each layer's properties.Virtual experiments using COMSOL Multiphysics to simulate impedance measurements are used to examine and improve the accuracy of the numerical inversion procedure. The results for the two dual-layer anodes studied show that MacMullin numbers follow expected trends, though the anodes are quite different from each other.

Dual-Layer Electrodes↗
DOE OSTI2024

In Situ Liquid Electron Microscope Cells Strongly Attenuate Electrochemical Behavior

In situ electrochemical Scanning Transmission Electron Microscopy (ec-S/TEM) probes the dynamics of electrochemical processes in real time and at high spatial resolution. This method utilizes microfabricated electrochemical cells, with enclosed dimensions of 800 μm wide, 800 μm long, and 1–2 μm tall, with electron-transparent silicon nitride windows. Working, counter, and reference electrodes in the configuration studied are deposited Pt with a planar surface area of about 1 μm 2 and a height of 50 nm. The cell confines the electrolyte within a 1 μm tall channel, and the microfabricated electrodes have a non-standard geometry, are coplanar, and are not front-facing. As such, standard assumptions of electrochemical experiments do not apply. COMSOL modeling was used to compare and contrast differences in electrochemical behavior between conventional experimental setups and the in situ ec-S/TEM cells. Cell height strongly affected voltammogram measurements: peak magnitudes and shifts correlate directly to cell height. Reactant concentration above the electrode quickly decreases during deposition, leading to a planar diffusion-dominated regime. The concentration gradient prioritizes particle growth along the thin edges of the electrode and not on the planar part of the electrode. Experiments done in the ec-S/TEM ex situ of a TEM verify the model’s accuracy.

Wittman, Reed M. [Oak Ridge National Laboratory (O↗