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

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↗

PV TOMCAT

PV TOMCAT predicts the outdoor operating temperature and output power of a PV module based on measured physical properties and outdoor conditions. It includes a ray tracing model of a PV module, which runs in the commercial SunSolve ray tracer. Output from this model is combined with meteorological time series data to produce a combined input file for a finite-element model of heat conduction and radiation. Finally, this finite-element model is run in the commercial COMSOL Multiphysics software, producing time series outputs of temperature and power.

Deceglie, Michael↗

Berkeley Lab Finite Element Framework (BELFEM) v0.1

The software program, referred to as BELFEM, is a specialized finite element code designed for the magnetodynamic modeling of high-temperature superconducting (HTS) tapes. It incorporates novel mixed finite element formulations, particularly the h-ϕ-formulation with thin-shell simplification, to efficiently simulate larger geometries. This methodology is extremely promising for predicting the electrodynamic performance of HTS tapes used in superconducting cables and magnets, offering the benefit of reduced computational cost. Compared to similar technologies like COMSOL Multiphysics and GetDP, BELFEM's performance benchmarking indicates superior efficiency in its thin-shell implementation. In the future, it will also support features like thermal coupling and inter-tape current sharing, enhancing its utility in research and development, particularly in nuclear fusion applications. The intent is to develop BELFEM as a robust and efficient tool for the HTS community, contributing to the analysis and design of superconducting cables and magnets.

Messe, Christian↗

Design and characterization of “flow-cell” integrated-flow active cooling for high-average-power ceramic optics

We used COMSOL Multiphysics to design a prototype actively cooled “flow-cell” substrate targeted at high-average-power applications, fabricated the prototype from cordierite ceramic, and tested the substrate under load in our thermal loading test stand. Sub-aperture testing revealed average-power handling up to 3.88-W/cm 2 absorbed power density, in excellent agreement with model predictions. Gratings fabricated on 2-in.-diam cordierite coupons were subjected to laser-damage testing and showed a damage threshold of 250 mJ/cm 2 .

47 OTHER INSTRUMENTATION↗

Conduction Heat Transfer through Solid in Porous Materials: A Comparative Study by Finite-Element Simulations and Effective Medium Approximations

The development of insulation materials with low effective thermal conductivity is essential for energy savings in various applications, including buildings, food services, pipe insulation, and refrigeration. Such materials can be developed by using micro- or nanoporous structures, as well as low-thermal conductivity gases and/or reduced pressure. A variety of effective medium approximation (EMA) models have been developed to study thermal transport through solid in porous structures. However, in many cases, the impacts of porosity on solid conduction are not well predicted by EMA models because of the assumptions made for simplification. Furthermore, the results can vary by a factor of up to 1.5, depending on the morphology of the pores. Hence, proper guidance is needed to choose the appropriate EMA model for a given morphology. This work presents a finite element method study using COMSOL Multiphysics software for various geometries, including hollow cubes, overlapping hollow spheres (normally stacked and tightly packed), and Voronoi structures, and compares the results to those obtained from various EMA models. Simulation results for the hollow cubes and Voronoi structures matched well with the Maxwell-Eucken and Russell model. The results for hollow spheres closely matched the Bauer model and the Glicksman model with fitting parameters. This work provides guidance on analyzing and designing insulation materials for energy savings in the future.

36 MATERIALS SCIENCE↗

Feasibility of Process Intensification of Water-Gas Shift Reaction Using a Microreactor with Integrated Cooling

Hydrogen is an increasingly attractive low-carbon energy carrier for a variety of stationary and mobile applications. The water-gas shift (WGS) reaction is a key processing step used for large-scale hydrogen production via the steam methane reforming process. However, the thermodynamics and kinetics of the reaction are such that standard two-stage adiabatic reactors used in these systems are large, increasing catalyst volume and cost. To intensify the process and realize the economical distributed production of hydrogen, adiabatic WGS microreactors with integrated cooling directly regulate the reaction temperature via integrated coolant channels to promote higher conversion within a smaller reactor volume. This study investigates the conversion efficiency of a single WGS microchannel operating under such cooling conditions. A COMSOL Multiphysics model is developed and validated with isothermal experimental data from the literature. The model is then used to evaluate improvements in conversion efficiency when the reaction is cooled via a specified wall temperature profile. Lastly, the model is modified to include cooling channels with a secondary fluid that can practically achieve a similar conversion profile as the specified wall temperature profile previously applied. Initial results show that reactor conversion can be significantly increased by the inclusion of appropriate cooling and that there is a potential for the recovery of energy from the reaction stream that can be used for other applications within the overall process.

energy recovery↗

Evaluating the Effects of Nitrogen Doping and Oxygen Doping on SRF Cavity Performance

Superconducting radiofrequency (SRF) cavities are resonators with extremely low surface resistance that enable accelerating cavities to have extremely high quality factors (Q₀). High Q₀ decreases the capital required to keep the accelerators cold by reducing power loss. The performance of SRF cavities is largely governed by the surface composition of the first §I{100}{nm} of the cavity surface. Impurities such as oxygen and nitrogen have been observed to yield high Q₀, but their precise roles are still being studied. Here, we compare the performance of cavities doped with nitrogen and oxygen in terms of surface composition and heating behavior with field. A simulation of the diffusion of oxygen into the bulk of the cavity was built using COMSOL Multiphysics software. Simulated results were compared to the actual surface composition of the cavities as determined from secondary ion mass spectrometry analysis. Understanding how these impurities affects performance allows us to have further insight into the underlying mechanisms that enable these surface treatments to yield high Q₀.

43 PARTICLE ACCELERATORS↗

Simulation of the dynamics of gas mixtures during plasma processing in the C75 Cavity

Plasma processing using a mixture of noble gas and oxygen is a technique that is currently being used to reduce field emission and multipacting in accelerating cavities. Plasma is created inside the cavity when the gas mixture is exposed to an electromagnetic field that is generated by applying RF power through the fundamental power or higher-order mode couplers. Oxygen ions and atomic oxygen are created in the plasma which breaks down the hydrocarbons on the surface of the cavity and the residuals from this process are removed as part of the process gas flow. Removal of hydrocarbons from the surface increases the work function and reduces the secondary emission coefficient. This work describes the initial results of plasma simulation, which provides insight into the ignition process, distribution of different species, and interactions of free oxygen and oxygen ions with the cavity surfaces. The simulations have been done with an Ar/?2 plasma using COMSOL® multiphysics. These simulations help in understanding the dynamics and control of plasma inside the cavity and the exploration of different gas mixtures.

Dhakal, Pashupati↗

Modeling Continuous-Flow Reactor Improvements over Batch Reactions for Enzyme Catalyzed Microsphere Surface Reactions

Our research objective was to study enzymatic reactions occurring on microsphere surfaces. Specifically, our work involved enzymatic polymerization reactions (chain lengthening) taking place both in batch and in a continuous flow reactor. My research involved creating simulations using COMSOL Multiphysics 5.4a to further study variable interactions within the continuous flow reactor. Our goal was to use suitable software to create a simulation platform suite that could be used for advanced research studies of continuous flow enzymatic surfaces reactions. In addition, we studied various characteristics of our present enzymatic polymerization reactor assembly, modeled inlet species flow (Species M) and estimated outlet species production (Species P). Our simulations spanned three levels of detail: Single-Bead, Packed-Bead Cell and Full Reactor Column. Furthermore, the Packed-Bead Cells were studied using three possible regular arrangement of solid spheres: Face-Centered Cubic (FCC), Hexagonal Close Packed (HCP) and Cubic Close Packed (CCP).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulating Dopant Diffusion in a Detalied Porous Structure

The National Ignition Facility (NIF) at Lawrence Livermore National Laboratory (LLNL) uses highpowered lasers to explore cutting edge fusion technology. Specialized optics, known as fiber lasers, generate and amplify the NIF’s high-powered lasers that facilitate fusion experiments. Fiber lasers are glass fibers that are doped with heavy elements to stimulate light emission. The fiber laser doping processing uses a predeposited layer of porous silica nanoparticles on the surface of the glass to absorb a dopant rich salt solution. Our project simulated the diffusion of dopants into porous silica with prescribed colloidal stacking geometries using COMSOL Multiphysics. With these simulations we identified the porosity and tortuosity of the porous nanoparticle structures as the primary parameters that will dictate how the dopants will diffusive into the porous surface of the glass. Porosity and Tortuosity are geometric parameters that can easily be input into more detailed simulations to estimate an effective diffusivity within the porous layer. Future work should account for surface reactions to accurately simulate the full-scale diffusion of the solution doping process.

36 MATERIALS SCIENCE↗

Numerical and Analytical Modeling of the Effect of Cracks on the Self-Inductance of a COTS YJ-41003-TC Toroid

COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μ r is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μ core >> μ air ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.

36 MATERIALS SCIENCE↗

Computational Fluid Dynamics Modeling to Simulate a Combined Reforming Process for Syngas and Hydrogen Production

An Oxygen Transport Membrane (OTM) combined reforming technology for producing syngas and hydrogen integrates the advantages of multiple processes—steam methane reforming (SMR), autothermal reforming (ATR), an air separation unit (ASU)—into a single integrated technology. The OTM consists of a primary reforming tube, in which desulfurized natural gas is partially reformed by steam at high pressure in the presence of a metal catalyst. This process is followed in series by a ceramic OTM with a secondary reformer, in which residual methane reforms and O 2 - ions react with a portion of the CO and H 2 fuel to provide the heat to support both primary and secondary reforming. Although the OTM combined reformer technology for syngas and H 2 production has been substantially developed in the last decade, several challenges that affect the overall production efficiency and reliability are yet to be fully understood, addressed, and resolved. Therefore, developing Computational Fluid Dynamics (CFD) models that incorporate fluid dynamics, mass transport, kinetics, heat transport, and structural mechanics is critical to understanding and minimizing the probability of tube failures during the startup and operation. In this report, an exhaustive literature review was performed to survey the current state of technology for producing syngas and H 2 using either conventional or renewable energy sources. The feedstocks reviewed include natural gas and coal for the conventional technologies, whereas biomass, solar, wind, and nuclear energy for the renewable technologies. The existing industry-grade COMSOL multiphysics models of OTM were upgraded for the latest software release. In addition, they were improved to help achieve grid and solver independence and were successfully ported on the ORNL high-performance computing clusters to speed up their run times. A 42% reduction in the simulation run time was achieved. A new higher-fidelity CFD model of an OTM tube was developed in the StarCCM+ simulation platform. This new model was designed to simulate various physics using first principles, e.g., turbulent flow, heat transfer, and chemical reactions while avoiding unnecessary simplifications. The resulting predictions were qualitatively assessed and provided useful insights into the multiphysics complexity of an OTM tube.

08 HYDROGEN↗

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↗

Analysis of Microwave Resonant Cavity Transducer Degradation Mechanisms

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with the flat wall of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the resonant frequency, which can be detected with a spectrum analyzer. Because the transducer is intended for immersion in a high temperature corrosive fluid, understanding of material degradation is crucial for estimation of transducer performance lifetime, and development of measurement interpretation algorithms. We conducted a preliminary computational investigation of relevant damage mechanisms of a stainless steel 316 cylindrical resonator in FLiBe salt. The two main damage mechanisms, creep and corrosion, were modeled using multiphysics COMSOL software. Degradation was modeled for a temperature range 500°C to 700°C. Coupling of the damage mechanisms was not considered. These models predict significant inelastic deformation at most temperatures due to creep, and qualitatively predict chromium depletion both along the liquid/solid interface and along the grain boundaries. An algorithmic approach for compensation of these degradation effects during fluid flow measurements will be developed in the future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Second Annual Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with the flat wall of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the resonant frequency, which can be detected with a spectrum analyzer. Because the transducer is intended for immersion in a high temperature corrosive fluid, understanding of material degradation is crucial for estimation of transducer performance lifetime, and development of measurement interpretation algorithms. We conducted a preliminary computational investigation of relevant damage mechanisms of a stainless steel 316 cylindrical resonator in FLiBe salt. The two main damage mechanisms, creep and corrosion, were modeled using multiphysics COMSOL software. Degradation was modeled for a temperature range 500°C to 700°C. Coupling of the damage mechanisms was not considered. These models predict significant inelastic deformation at most temperatures due to creep, and qualitatively predict chromium depletion both along the liquid/solid interface and along the grain boundaries. An algorithmic approach for compensation of these degradation effects during fluid flow measurements will be developed in the future work. To validate sensor physics, we have performed proof-of-principle test of flow sensing in water. For this test, we have developed a cylindrical resonator for K-band, which was machined from brass. The cavity was excited through WR-42 waveguide through a subwavelength hole on the side of the wall of the cylinder. To increase the spectral signal visibility, we developed a signal processing method for baseline subtraction. A flow loop for proof-of-principle test of transducer performance in water was assembled. A commercial flow meter was installed in the loop for reference measurements. Cylindrical cavity was excited in the TE 011 mode with resonant frequency f ≈ 17.8GHz. Frequency shift of cavity spectral response was obtained by gradually increasing water flow rate from 0 to 60gpm. Corresponding monotonic increase of resonant frequency shift by several MHz was observed. Approximate figure of merit of sensitivity to flow rate is 100KHz/GPM. In addition, we have identified an existing liquid sodium experimental setup for demonstration of flow sensing in environment similar to that of an advanced reactor. The setup consists of a cylindrical vessel and center feed line, where transducer inserted through the lid will measure velocity of the impinging liquid jet. As a calibration experiment, we have assembled a water vessel with center feed with the same dimensions as those of the liquid sodium setup. We have also developed and insertion probe consisting of a 50cm brass waveguide enclosed in protective SS316 tube. Using the water loop, we have demonstrated feasibility of sensing the impinging liquid jet in the vessel.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Report of METL Mutual Inductance Level Sensor Development for Use in Liquid Metals – FY2023

A robust mutual inductance level sensor (MILS) for use in high temperature liquid metals has been developed at Argonne National Laboratory (ANL). The current design utilizes mineral insulated cables with a 300-series stainless steel sheath, magnesium oxide insulation, and a single copper conductor. Two coils are wrapped on a common core made of a 300-series stainless steel tube. One coil is energized with an alternating current (AC) source, and this electromagnetically couples with the second coil to generate an induced voltage. The voltage is measured to determine the mutual induction between the coils, and the mutual inductance measurement can be used to determine the level of a nearby liquid metal volume. The MILS can be located in an isolating thimble that is fully sealed to a vessel containing liquid metal for ease of maintenance and replacement. When the primary coil is energized with the AC source, the coil not only electromagnetic (EM) couples with the second coil, but also the surrounding liquid metal. The EM coupling with the surrounding liquid metal reduces the EM coupling with the second coil, producing an inverse-linear relationship between liquid metal level and secondary coil voltage. This has all been demonstrated at the Mechanisms Engineering Test Loop (METL) liquid sodium facility at ANL. The most current iteration of the MILS system is the MILS-MK-II. This sensor system has been commissioned in a non-sodium test stand where calibrations were performed using a sodium analog. The calibrations proved to be highly linear and repeatable. The MILS-MK-II has been installed in the METL expansion tank where high temperature sodium tests have been performed. The MILS-MK-II has been calibrated against a known standard at temperature of 300°C, and the calibrations have proved to be highly linear and repeatable. The calibrated MILS-MK-II has been in operation in the METL facility for several thousands of hours at temperatures around 300°C. The experimental data has been used to develop and validate electromagnetic finite element models in COMSOL Multiphysics, and now these models can be used to advance the development of the sensor system. Next steps will include temperature compensation to allow for operation at various temperatures. Additionally, efforts to incorporate a self-calibration methodology are underway.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Design, Modeling, and Testing of a Novel Inductor for Electric Vehicles: Iron Nitride Soft Magnetic Composites

New technology for electric vehicles (EVs) must meet the requirements of higher energy usage, lower costs, and more sustainable source materials. One promising material for EV power system components is iron nitride (IN) soft magnetic composites (SMCs) because of their competitive magnetic properties and high abundance of the source materials. As part of an ongoing program at Sandia National Laboratories, this project focused on using computer modeling to optimize the prototyping process for an iron nitride SMC toroidal inductor to reach a target inductance of 600 μH. Four inductors with different combinations of wiring (26 AWG and 20 AWG) and vol% loading of iron nitride (65 vol% and 50 vol%) were fabricated at Cal Poly and characterized using an LCR meter. These inductors were also modeled using COMSOL Multiphysics™ with the Magnetic Fields module. The inductance data from the experiment and the model show that the 65 vol% IN prototypes and models agree with about 8% difference, while the 50 vol% IN samples show about a 9% difference between the prototype and the model. These results suggest that the model can predict inductance with both accuracy and precision with low confidence for the given sample size of four. An additional parameter of AC resistance is studied but the AC resistance results from the inductors and from the model generally do not agree closely, suggesting that the current model used in the project does not fully capture the mechanisms behind AC resistance of the inductor. With the focus of the project on inductance, the percent difference results of less than 9% across the four inductors that were tested increases confidence in the model’s predictive capabilities for inductance only. Using the inductance results from both the model and experiment, the final suggested inductor design is a 65 vol% core with 150 windings of 20 AWG wire that is 8 cm across and 1.5 cm tall to reach the inductance goal of 600 μH based on analysis using the optimized COMSOLTM model.

33 ADVANCED PROPULSION SYSTEMS↗

Novel Full-Ceramic Multi-Tubular Membrane Systems for Pre-Combustion CO2 Capture with Simultaneous H2 Production: Fabrication, Performance Testing, and 3D CFD Modeling

Inorganic membranes show promise for application in pre-combustion CO2 capture with simultaneous H2 production. State-of-the-art systems for use under high temperature and pressure conditions consist of multiple membrane tube bundles prepared in a "candle-filter" configuration, in which the membrane tubes are open at one end and sealed at the other. This configuration is used for practical reasons, specifically the need to minimize potential problems due to thermal expansion mismatch, at high temperatures, between the ceramic tube bundle and the steel housing. The primary technical problem with the candle-filter configuration for use in commercial-scale installations is the inability to purge the permeate side (typically the tube side), a feature that is crucial for high H2 recovery. In this study, we fabricated dual-end open, commercial-size ceramic multiple-tube bundles made of zeolite, palladium (Pd), and carbon molecular sieve (CMS) membranes that enable permeate-side (tube-side) purge for gas separation applications. Experimental gas separation data with these membrane bundles under harsh operating conditions (temperatures up to 350 and pressures up to 800 psig), to be presented at the meeting, manifest excellent performance. Parallel to the membrane bundle construction and testing efforts, we have also developed a detailed 3D CFD modeling package using COMSOL Multiphysics software to gain more insight into the effect on the H2 purity and recovery of the detailed geometry of the multi-tubular membrane system, including the number of tubes used, their dimensions, and placement in the bundle, as well as the number, type, and positioning of internal baffles and other flow-enhancement accessories. The CFD package is validated with experimental data from different systems (1-tube, 3-tube, and 19-tube bundles), and shows high accuracy in predicting the experimental results (<5 % error in all cases). The results of our study show that the detailed internal geometry of these multi-tubular membrane systems has a considerable impact on the performance of the system as the flow maldistribution within the shell-side can substantially decrease (>40%) the H2 recovery.

20 FOSSIL-FUELED POWER PLANTS↗