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 163 records · Page 9

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↗

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↗

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↗

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↗

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↗

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↗

Joule-Heated Molten Regolith Electrolysis Reactor Concepts for Oxygen and Metals Production on the Moon and Mars

The maturation of Molten Regolith Electrolysis (MRE) as a viable technology for oxygen and metals production on explored planets relies on the realization of the self-heating mode for the reactor. Joule heat generated during regolith electrolysis creates thermal energy that should be able to maintain the molten phase (similar to electrolytic Hall-Heroult process for aluminum production). Self-heating via Joule heating offers many advantages: (1) The regolith itself is the crucible material, it protects the vessel walls (2) Simplifies the engineering of the reactor (3) Reduces power consumption (no external heating) (4) Extends the longevity of the reactor. Predictive modeling is a tool chosen to perform dimensional analysis of a self-heating reactor: (1) Multiphysics modeling (COMSOL) was selected for Joule heat generation and heat transfer (2) Objective is to identify critical dimensions for first reactor prototype.

Sibille, Laurent↗

DC Modeling of 4H-SiC nJFET Gate Length Reduction at 500°C

The development of robust, high-performance integrated circuits (ICs) will enable numerous potential NASA missions of current interest, including long-duration robotic missions exploring the 460°C surface of Venus. Currently, NASA is looking towards SiC-based devices to provide such a solution. However, the current NASA silicon carbide (SiC) JFET device with a channel length of 6 μm (recently fabricated Gen. 11 ICs) limits mission-relevant circuit capabilities. In this study, we combined experiments with simulations to explore two straightforward fabrication strategies (shallow n− and deep n+) to reduce the SiC JFET channel length while maintaining turn-off behavior needed to realize 500 °C circuit operation. First, the material properties for 4H-SiC were implemented in COMSOL and the Poisson equation was solved for twelve 2D device designs. Then, based on the insights gained from the distribution of the electron concentration, electrostatic potential, and electric field for twelve designs (with three fabrication strategies), simulation for a 1 μm gate length nJFET with the turn-off performance comparable to the state-of-the-art.

4H-SiC nJFET↗

DC Modeling of 4H-SiC nJFET Gate Length Reduction at 500°C

The development of robust, high-performance integrated circuits (ICs) will enable numerous potential NASA missions of current interest, including long-duration robotic missions exploring the 460°C surface of Venus. Currently, NASA is looking towards SiC-based devices to provide such a solution. However, the current NASA silicon carbide (SiC) JFET device with a channel length of 6 μm (recently fabricated Gen. 11 ICs) limits mission-relevant circuit capabilities. In this study, we combined experiments with simulations to explore two straightforward fabrication strategies (shallow n− and deep n+) to reduce the SiC JFET channel length while maintaining turn-off behavior needed to realize 500 °C circuit operation. First, the material properties for 4H-SiC were implemented in COMSOL and the Poisson equation was solved for twelve 2D device designs. Then, based on the insights gained from the distribution of the electron concentration, electrostatic potential, and electric field for twelve designs (with three fabrication strategies), simulation for a 1 μm gate length nJFET with the turn-off performance comparable to the state-of-the-art.

4H-SiC nJFET↗

Comparing Acoustic Prediction Methods for Additively Manufactured Porous Strutures

While macroscale methods for predicting the acoustic properties of porous structures have been popular in the past, they often require time-consuming manufacturing and testing workflows. Meanwhile, microscale approaches allow the prediction of transport parameters based exclusively on a periodic structure’s unit cell geometry. Here, we compare these methods to predict the characteristic impedance of additively manufactured porous structures. We use the microscale approach to estimate the geometry’s transport parameters, then predict the characteristic properties using the Johnson-Champoux-Allard (JCA) model. We measure the acoustic properties of the printed structures using a normal incidence impedance tube and estimate the transport parameters using an inverse characterization approach. We use the two-thickness method as a macroscale approach to predict the characteristic properties from the measured surface impedances of two sample thicknesses. Finally, we compare these characteristic prediction methods. Our results show that the inverse characterization and two-thickness methods offer the closest match to the measured values at low frequencies.

impedance↗

COMSOL Results for the Nominal Steady-State Operation of the Proposed 95-MW LEU Silicide Core for HFIR Conversion

Engineering design studies are being performed to determine the feasibility of converting the High Flux Isotope Reactor (HFIR) from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel at Oak Ridge National Laboratory. This activity is sponsored by the Office of Reactor Conversion and Uranium Supply (ORCUS) under the auspices of the US Department of Energy National Nuclear Security Administration’s Office of Material Management and Minimization. HFIR is a very high flux, pressurized, light water–cooled and moderated, flux trap–type research reactor with a core made of involute shaped U 3 O 8 /Al cermet fuel plates and coolant channels. HFIR currently operates at a thermal power of 85 MW and supports key national and international missions in neutron scattering, isotope production, materials/fuels irradiation, neutron activation analysis, gamma irradiation, and neutrino research. Advanced multiphysics computational fluid dynamics models have been developed in the COMSOL Multiphysics software to simulate the steady-state operating conditions for the proposed low-and high-density LEU U 3 Si 2 -Al (uranium silicide dispersion) fuel designs. The COMSOL models for HFIR inner and outer fuel element models incorporate various essential inputs and physics such as spatially dependent nuclear heat deposition, multilayer heat conduction, conjugate heat transfer, turbulent flows (using Reynolds-averaged Navier Stokes turbulence models), structural mechanics (thermal–structural interactions and fuel swelling), and oxide layer build-up. This report presents the best-estimate thermal hydraulics results for the low- and high-density optimized silicide LEU core designs at 95 MW steady-state nominal operation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Design and Analysis of the Open-Surface Slow Li Flow Divertor and Comparison to the Fast Li Flow Divertor

In the ongoing U.S. project, “Liquid Metal Plasma Facing Components,” sponsored by the U.S. Department of Energy, efforts have been taken to develop two open-surface divertor designs for the Fusion Nuclear Science Facility using liquid lithium (Li) as a heat and particle flux removal media. The main focus of this study is the design and analysis of a slow (~1 mm/s) and thin (<1 mm) open-surface Li flow divertor with a Li-cooled substrate, which is then compared with an earlier design of a fast (up to 10 m/s) and thick (~0.5 cm) Li flow divertor with the substrate cooled with helium. The slow Li flow divertor design is based on the original LiWall concept developed at the Princeton Plasma Physics Laboratory. Such a thin and slow Li layer can remove the particle flux by reducing the recycling flux, while the heat flux is removed mainly through the heat sink located beneath. In the present study, the heat sink is provided through a Li cooling flow inside the substrate of reduced activation ferritic/martensitic steel. By performing a multiphysics analysis with COMSOL that included liquid-metal magnetohydrodynamics (MHD), heat transfer, and structural mechanics, the impact of various factors on the divertor heat removal capability, such as Li flow velocity, MHD effects, and inlet velocity boundary condition, were examined. Based on comparisons of the two divertor designs, it was shown that the fast-flow divertor significantly outperformed the slow-flow design, whose heat removal capability was limited to ~1 to 2 MW/m2.

Jiang, Yuchen [ORNL] (ORCID:0000000283280284)↗

3D multiphysics simulation and analysis of a low temperature liquid metal magnetohydrodynamic power generator prototype

Magnetohydrodynamic (MHD) power generation is an interdisciplinary process, which was studied heavily in the late 1960’s but interest soon subsided. However, the current socioeconomic climate demands more diversity in renewable energy sources. MHD generators using liquid metal as the working fluid represent a viable option for this augmentation. Additionally, the steady increase in accessible computational capabilities provides a powerful new tool for interdisciplinary modeling and multiphysics simulations. These factors suggest that this is an opportune time to reevaluate MHD power generation. This paper provides mathematical background as well as a detailed multiphysics model of a low temperature liquid metal MHD power generator built in COMSOL. Simulation details are provided in addition to a comparison to experimental results. Finally, the prototype studied in this paper is the first proof of concept of a permanent magnet excited liquid metal MHD generator using Gallium as its working fluid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

ExoToComsol v.0.1

SAND2024-08641O ExoToComsol models capability-sharing. For example, the software can develop a multiphysics modeling chain using physics modeling capabilities from different research institutions using either SIERRA or COMSOL. The converter program allows exporting user-defined, region-of-interest from a large FE model, which results in a finer mesh of the submodel that can be used for subsequent analysis. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Rahman, Farhan↗

Simulating Thermoelectric Devices Using the MOOSE Framework

Thermoelectric generators (TEG) are devices that generate energy by converting heat into electricity or provide cooling via the Peltier effect. This feature of thermoelectric devices originates from the Seebeck, Peltier, Thomson, and Joule heating effects. TEGs can be applied in energy and thermal management systems such as waste heat recovery and refrigeration, respectively. Thermoelectric device design is influenced by the material selection and the device's geometry operating conditions. Therefore, predicting, verifying, and validating thermoelectric device performance using simulations tools is essential to deploying thermoelectric devices in industry. The Multiphysics Object-Oriented Simulation Environment (MOOSE) Framework is an open-source simulation tool capable of modeling simple to complex systems. In this work, we demonstrate MOOSE's thermoelectric device modeling capabilities by simulating a unicouple, module, and exhaust gas recovery system. The Seebeck, Peltier, Thomson, and Joule heating physics are implemented into MOOSE. The MOOSE thermoelectric physics were thoroughly verified and validated using published COMSOL® results and experimental data. In addition, thermoelectric modules were integrated into an exhaust gas recovery system using the MOOSE MultiApp function as a demonstration of the model's ability. The verification and validation results and exhaust gas heat recovery system showcases MOOSE's capability to model thermoelectric devices and integrate these devices into practical energy systems.

42 - ENGINEERING↗