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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Charge carrier motion and effect of fixed oxide charge in a microstructured silicon radiation detector

Signal formation in a microstructured semiconductor neutron detector is more complex than in planar diode geometry. Three-dimensional microstructures are laterally smaller than the ionization cloud length, and the electric fields may be weak enough to exhibit plasma time effects. This work is the first detailed treatment of charge carrier motion in these complex semiconductor devices to replicate the time profile and signal magnitude. Simulations were performed using COMSOL Multiphysics to investigate various parameters that affect the propagation of the charge cloud. It was observed that the size of the simulated three-dimensional structure had an impact on the induced current pulse, indicating the importance of simulation geometry optimization to accurately simulate charge cloud expansion. COMSOL Multiphysics was used to replicate accurate charge creation profiles using energy deposition information imported from radiation transport codes. A detailed simulation methodology is presented to benchmark preamplifier event pulses along with complexities in modeling the charge carrier motion along the etched microstructured trenches with Si–SiO 2 boundary conditions, including fixed oxide charge and interface trapping.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Lithium-ion battery physics and statistics-based state of health model

A pseudo-2d model using COMSOL Multiphysics® software is developed to simulate performance and performance degradation of Li-ion batteries consisting of layered and olivine cathodes with graphite anode when subjected to peak shaving grid service. Multiple degradation pathways are considered, including solid electrolyte interphase (SEI) formation and breakdown at the anode, cathode dissolution and its synergistic effect on SEI formation at the anode. The model is validated by simulating commercial cylindrical cell performance. A global model is developed to simulate performance across all chemistries, along with individual chemistry models using global model parameters as initial values. There is good agreement between these models for various optimization parameters such as SEI equilibrium potential, cathode dissolution exchange current density, solvent diffusivity in the SEI and SEI ionic conductivity. To circumvent time constraints related to the COMSOL model, a 0d global model is developed which fits data well and provides more clarity on differences in cathode dissolution exchange current density. Again, good agreement for various optimization parameters is obtained among the COMSOL global & individual chemistry models and the 0-d model. The lessons learned from the physics-based model is used to develop a top down statistics-based model using current, voltage and anode volumetric change per mole lithium intercalated, along with their interactions as degradation predictors. This model predicts out of sample degradation for multiple grid services and electric vehicle drive cycle with high accuracy and provides the pathway to develop an efficient battery management system combining machine learning and findings from physics-based computationally intensive algorithms.

Crawford, Aladsair J.↗
DOE OSTI2023

A COUPLED DETERMINISTIC TRANSPORT CALCULATION IN COMSOL USING PROPER ORTHOGONAL DECOMPOSITION

A coupled physics problem consisting of radiation transport and heat transfer was modeled using COMSOL Multiphysics to simulate a criticality accident. Reduced order models (ROM) were created to accelerate the radiation transport calculation through Proper Orthogonal Decomposition (POD). Additionally, on-the-fly neutron cross sections were generated from a non-linear function through a feed forward neural network. The Livelink for MATLAB module allowed for data transfer between an otherwise independent solver and in-house modules. A data driven approach to POD was developed to act as a surrogate model to couple with heat transfer, with the POD module written in MATLAB coupled to an isolated heat transfer solve in COMSOL. A coupled transient solution resulted in a good approximation of the flux and temperature, with the maximum error being 3E-2 and 2.5E-3 respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
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 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 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 OSTI2022

Fully Distributed Multi-Material Magnetic Sensing Structures for Multiparameter DAS Applications

This dissertation demonstrates the first of its kind distributed magnetic field sensor based on a fiber optic distributed acoustic sensing (DAS) scheme. Ferromagnetic nickel and Metglas® were dispersed internally within a fiber optic preform and then drawn on an in-house fiber optic draw tower to lengths in the kilometers. Due to the close proximity of the ferromagnetic metals and fiber optic core, the magnetostrictive strain response of the ferromagnetic materials when exposed to a magnetic field would perturbate within the fiber cladding and transfer that strain, internally, to the fiber optic core. Strain resulting from the magnetostrictive effect allows the DAS based sensor to accurately translate strain into readable magnetic field data. Due to the high sensitivity seen in this sensor design, multiparameter sources, acoustic and magnetic fields, were tested and validated and a three dimensional magnetic-field vector sensor was proposed. Numerical analysis of the novel sensor design was first implemented using COMSOL Multiphysics, where inputs such as magnetostrictive element shape, size, distance, and number were first investigated. Upon optimizing system constraints, the sensor design was further modified such that single mode operation was consistent across multiple fiber draws while retaining high strain transfer from the ferromagnetic elements to the fiber optic core. Ferromagnetic material selection was evaluated as a function of the saturation magnetostriction constants and a total of 4 modules were used to fully characterize the complex physics involved in this sensor design. All fabrication and testing were performed in-house using a full scale 3-story fiber draw tower and custom environmental testing stations to imitate naturally occurring events such as magnetic or acoustic point sources. A unique stacking method was used to embed ferromagnetic nickel and Metglas® into a fiber optic preform which when combined with a custom fiber draw process resulted in consistent multi-material fibers drawn to lengths of 1-km. In-house testing facilities included different types of electromagnetic generators, in addition to a soil test bed, and an outdoor test bed which allowed 100 meters of fiber to be tested simultaneously. All tested sensors demonstrated high strain transfer capabilities on the order of 0.01-10 μϵ depending on the materials used, ferromagnetic rod number, and core to metal spacing. Due to the sensitivity of the system the difference between AC and DC was distinct, and directional magnetostriction was studied. Transverse and longitudinal magnetic wave propagation was controlled through a solenoid and rectangular Helmholtz coil, both built in-house. A three-dimensional magnetic field vector sensor was proposed due to the success of the magnetic field sensor, and a design was proposed and initially tested to validate direction as a function of field strength and distance. To summarize, this dissertation explores the first fully distributed magnetic field sensor using DAS based techniques and one of the first multi-material fiber draw processes which can produce consistent single mode fiber up to 1-km. Due to extensive FEA modeling, multiple iterations of the magnetic sensor were fully characterized and an equation describing the relationship between sensor design and strain transfer has been created and validated experimentally. Multi-parameter tests including acoustic and magnetic fields were implemented and an algorithm was developed to separate the mixed signals. Finally, a test was performed to demonstrate the feasibility of sensing magnetic fields directionally. Cumulative results demonstrate a high-quality sensor alternative to current designs which may surpass other magnetic sensors due to innate multi-parameter capabilities, in addition to the inexpensive production cost and extremely long operating lengths.

02 PETROLEUM↗
NASA NTRS2015

Simulation Apps Bring Us Closer to Mars

My work with mathematical modeling and computer simulations began in earnest in 1987 when I signed on with the National Aeronautics and Space Administration, better known as NASA. At the time I was developing computer simulations to support the design of carbon dioxide (CO2) removal systems for the International Space Station (ISS) life support system. CO2 is a human metabolic waste, produced at a rate of a kilogram per day per crewmember. CO2 must be removed from the crew cabin, as it will quickly become detrimental to crew health. The systems used to remove CO2 are typically based on sorption processes, which include complex interactions of heat transfer, mass transfer, and gas flow through porous media. At this time, there were no commercial options for solving a set of coupled mathematical models such as partial differential equations (PDEs); you either force-fit the physics into a thermal analysis package or wrote your own code, complete with discretization, meshing, and solution algorithms. Unfortunately, coding CO2 removal processes from scratch did not allow an appropriate focus on the challenging yet crucial task of understanding and capturing the underlying physics via appropriate mathematical models. Configuration control was often simply ignored due to tight schedules, resulting in a code that would be quicker for a second party to rewrite than modify. In the early 2000s, I decided to move CO2 removal simulations to a platform with built-in meshing, solvers, and postprocessing capabilities, and that could solve user-defined multiphysics PDEs. The program chosen to meet these needs was COMSOL Multiphysics®. Along with freeing the engineer to focus on the underlying physics, a degree of configuration control was automatically achieved via a consistent user interface, thus allowing COMSOL users to share computer models. My team has developed simulations that are already providing a valuable debugging capability for the ISS CO2 removal system, and will provide guidance in upgrades to that system.

Knox, Jim↗
NASA NTRS2008

A Thermal Precipitator for Fire Characterization Research

Characterization of the smoke from pyrolysis of common spacecraft materials provides insight for the design of future smoke detectors and post-fire clean-up equipment on the International Space Station. A thermal precipitator was designed to collect smoke aerosol particles for microscopic analysis in fire characterization research. Information on particle morphology, size and agglomerate structure obtained from these tests supplements additional aerosol data collected. Initial modeling for the thermal precipitator design was performed with the finite element software COMSOL Multiphysics, and includes the flow field and heat transfer in the device. The COMSOL Particle Tracing Module was used to determine particle deposition on SEM stubs which include TEM grids. Modeling provided optimized design parameters such as geometry, flow rate and temperatures. Microscopy results from fire characterization research using the thermal precipitator are presented.

Meyer, Marit↗
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 OSTI2025

Design and Analysis of a Mutual Inductance Level Sensor for Liquid Metals

Here, this article describes the design and analysis of an electromagnetic level sensor for use in high-temperature liquid metal systems. The mutual inductance level sensor (MILS) described in this work was fabricated using two single-conductor mineral insulated cables wrapped in a bifilar fashion around a stainless steel tube core and was housed in an isolating thimble that preserved the pressure boundary of the test vessel. Two sensor variations were fabricated that differ only in active length, 1016 and 1778 mm. Experimental data were collected using the 1016-mm sensor (MILS-MKII-040) that demonstrated a sensitivity of 9.2 μ V/mm in a room temperature testing stand that used solid aluminum as a surrogate for liquid metal. Experimental data were collected using the 1778-mm sensor (MILS-MKII-070) that demonstrated a sensitivity of 6.9 μ V/mm in the high-temperature (300 ° C) liquid sodium environment at the mechanisms engineering test loop (METL) of Argonne National Laboratory. The sensor performance was found to be repeatable over the course of several months, with roughly ±1% deviation from nominal output. Finite element models were developed in COMSOL Multiphysics that fully describe each test setup, and the models were validated using experimental data. The validated COMSOL models were used to perform an array of analyses that examined the performance of the sensor in differing environments. Maximizing the coil diameter inside the isolating thimble was found to maximize the signal and sensitivity of the sensor. An optimal operating frequency was found near 1000 Hz using both experimental data and COMSOL. The influence of a metallic thimble surrounding the sensor and a metallic sensor core was quantified and found to be negligible at the optimal operating frequency. The sensitivity of the sensor was quantified when monitoring the level of additional liquid metals. These include lead, lead-bismuth eutectic (LBE), sodium-potassium alloy (NaK), and lithium (in addition to sodium). The sensitivities were quantified using liquid metal properties at 350 ° C and 650 ° C. The geometry of the test stand model, all material properties used in the model, and the results are presented in a manner that allows the reader can replicate the model and perform additional analyses.

COMSOL↗
NASA NTRS2015

Design of a Thermal Precipitator for the Characterization of Smoke Particles from Common Spacecraft Materials

A thermal precipitator (TP) was designed to collect smoke aerosol particles for microscopic analysis in fire characterization research. Information on particle morphology, size and agglomerate structure obtained from these tests supplements additional aerosol data collected. Modeling of the thermal precipitator throughout the design process was performed with the COMSOL Multiphysics finite element software package, including the Eulerian flow field and thermal gradients in the fluid. The COMSOL Particle Tracing Module was subsequently used to determine particle deposition. Modeling provided optimized design parameters such as geometry, flow rate and temperatures. The thermal precipitator was built and testing verified the performance of the first iteration of the device. The thermal precipitator was successfully operated and provided quality particle samples for microscopic analysis, which furthered the body of knowledge on smoke particulates. This information is a key element of smoke characterization and will be useful for future spacecraft fire detection research.

Smoke↗
NASA NTRS2023

Multiphysics Modeling of Electrode Dissolution under Hydrodynamic Conditions in Low Conductivity Water

Efficient plug-flow electrochemical reactors are characterized by their ability to effectively convert species with significant ionic strength and/or with the assistance of a well-behaved electrolyte. Because of the classic electroconductivity requirement, electrochemical reactors are not particularly attractive options for applications with low-conductivity conditions, especially for low-pressure feeds. However, these reactors can still find application in dedicated processes such as ultra-clean water conditioning and biological-coupled separation systems. Consequently, a multi-physics model was developed using COMSOL Multiphysics to characterize different reactor-design conditions based on the Tertiary Current Distribution which accounts for the effects of expected variations in electrolyte composition and ionic strength on the electrochemical process, as well as solution resistance and electrode kinetics. The modelling framework presented here employs single-phase laminar fluid flow, the Nernst-Planck equation, the water-based electroneutrality condition, and concentration-dependent overpotentials. This investigation particularly explores the electrolysis process of “anodic dissolution,” in which elemental species dissociate from the solid matrix of the electrode and enter the liquid phase of the electrolyte as soluble ions. The reactor configuration of interest has a rectangular parallel-electrode design and is operated galvanostatically. A series of parametric studies are carried out to inspect the response of the system when hydrodynamic, kinetic, and geometric conditions are changed. The parametric sweep resulted in variable effluent concentrations and system voltage, which reveal underlaying optimization patterns for electrode dissolution in low conductivity water.

Multiphysics↗
NASA NTRS2014

Multiphysics Modeling of an Annular Linear Induction Pump With Applications to Space Nuclear Power Systems

An annular linear induction pump (ALIP) that could be used for circulating liquid-metal coolant in a fission surface power reactor system is modeled in the present work using the computational COMSOL Multiphysics package. The pump is modeled using a two-dimensional, axisymmetric geometry and solved under conditions similar to those used during experimental pump testing. Real, nonlinear, temperature-dependent material properties can be incorporated into the model for both the electrically-conducting working fluid in the pump (NaK-78) and structural components of the pump. The intricate three-phase coil configuration of the pump is implemented in the model to produce an axially-traveling magnetic wave that is qualitatively similar to the measured magnetic wave. The model qualitatively captures the expected feature of a peak in efficiency as a function of flow rate.

Kilbane, J.↗
DOE OSTI2025

Precise magnetic field mapping of the EMPHATIC Phase 1 magnet with COMSOL

A compact Halbach array magnet is used to measure the momentum of the secondary particles in EMPHATIC(Experiment to Measure the Production ofHadrons At a Test beam InChicagoland). Hall probe data was taken for the central cylindrical bore of the magnet and a field map was constructed. COMSOL Multiphysics® Software is used for modeling the magnet and constructing the corresponding magnetic field map. Here, we present a fitting approach where the hall probe data is used to determine a 1mm-spacing map of the entire volume of the magnet using COMSOL. The new map will allow for linear interpolation within the volume, and expand the map to outside the measurement volume, thus increasing the acceptance and precision of EMPHATIC's tracking system.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Multiphysics simulations of a steady-state lower hybrid current drive antenna for the FSNF

The Fusion Nuclear Science Facility (FNSF) is a proposed tokamak reactor with the mission to investigate operation of a fusion reactor in a nuclear environment. The high neutron fluence component of the FNSF mission requires steady-state operation for extremely long pulses (t_{pulse} ∼ months) at full power. Plasma sustainment and current drive will be critical components of a successful FNSF. COMSOL Multiphysics software is used for combined radiofrequency (RF) and thermal simulations of the lower hybrid current drive (LHCD) antenna system. These simu- lations consider the resistive RF losses in the antenna including realistic surface roughness and a range of potential materials. The thermal analysis adds volumetric nuclear heating, plasma heat flux on leading edges, and electromagnetic radiation from the plasma to the RF heating calculated by COMSOL. Additional neutronics calculations have been performed to determine the impact of these antenna designs on activated waste disposal for the materials considered. The simulations show that it is technically feasible to implement a fully-active multi-junction (FAM) rather than a passive-active multi-junction (PAM) style of antenna if the septum between adjacent waveguides is sufficiently wide and the thermal conductivity of the structural material is sufficiently high. The FAM has the benefit of higher achievable power density with respect to the PAM, which results in a more compact antenna with potentially lower impact on neutron shielding and tritium breeding. These considerations point to tungsten rather than steel as the preferred structural material in constructing the antenna.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2026

Calculating the space-charge-limited current density for nonplanar geometries by simulating the charge-free electric field

Calculating the space-charge-limited-current density (SCLCD) for a complicated diode geometry often requires computationally expensive particle-in-cell (PIC) codes. Here, this paper addresses this issue by using the charge-free electric field $E_0$ calculated using COMSOL Multiphysics to determine local and global SCLCD. The SCLCD obtained by using the surface average of $|E_0|^2$ on the cathode recovers theoretical results for one-dimensional (1D) planar, cylindrical, and tip-to-tip geometries in appropriate limits. We next compared tip-to-tip calculations with the SCLCD obtained using the PIC code Empire. The SCLCD calculated using COMSOL agreed well with Empire for flatter 1D tip-to-tip geometries and diverged with increasing sharpness. Physically, Empire predicts lower SCLCD than COMSOL because the electrons spread due to concentrated space-charge at the tip, whereas theory assumes that the electrons follow the charge-free electric field lines. We further assess the behavior of the SCLCD for tips protruding from the centers of flat, circular plates of various areas. Larger plate areas with constant tip size recover the 1D planar SCLCD globally and 1D tip-to-tip SCLCD locally, while reducing the difference between Empire and COMSOL calculations since larger plates capture more of the emitted electrons, reducing SCLCD suppression due to beam spreading. These results show that charge-free electric field simulations can be used to determine the SCLCD without needing to simulate particle dynamics in PIC.

Wright, Jack K. [Purdue Univ., West Lafayette, IN ↗