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 73 records · Page 4

A thermodynamics-based damage model for the non-linear mechanical behavior of SiC/SiC ceramic matrix composites in irradiation and thermal environments

A damage model is developed and validated with experimental data for the non-linear mechanical behavior of SiC/SiC composite materials in nuclear applications. Cyclic thermal and mechanical loading associated with neutron irradiation effects of these composites leads to wide-spread and progressive micro-cracking that leads to loss of thermal conductivity and further enhancement of thermo-mechanical damage. A physics-based model of wide-spread micro-cracking is developed within the thermodynamic framework of continuum damage mechanics. Evolution equations for damage parameters that describe the growth of continuum damage are developed, where the material variables are obtained from experiments. The model novelty is in coupling mechanical, thermal, and irradiation damage through a consistent thermodynamic framework, including loss of thermal conductivity due to the evolution of mechanically induced micro-cracks. A number of thermo-mechanical experiments were conducted to confirm model assumptions. The model is shown to be validated with out-of-pile experiments, and then implemented using commercial finite element code COMSOL to the fuel cladding problem with normal and high radiation dose cases.

Materials Science↗

Multiphysics modeling of accelerators through code integration

This work aims to improve the ability of particle accelerator researchers to develop high-performance accelerator cavity designs by creating an overall multiphysics framework that integrates and couples existing application codes. This framework will allow accelerator researchers to build multiphysics models that will optimize cavity design, improve understanding of whole-device performance, and reduce the development and fabrication costs of accelerator research. We utilize the open-source VizSchema data standard as an intermediate data structure interface layer to standardize interfaces between individual application codes. VizScema is extensively documented online, and plugins for VizSchema are available for popular visualization packages, including VisIt and ParaView. Currently, the work focuses on coupling the EM field solver COMSOL and the electron gun code MICHELLE to allow COMSOL field-solve results to be seamlessly used by MICHELLE for particle-solve. Later work will extend this integration to include other fields, particles, and thermodynamics simulation codes.

43 PARTICLE ACCELERATORS↗

PyGDH: Python Grid Discretization Helper

Mathematical models expressed in the form of discretized equations play an important role inmany scientific disciplines. In our experience, few domain scientists have sufficient backgroundin numerical computing (or the time required to acquire such a background) to use manyflexible and powerful but complex open source packages, such as FEniCS (Alnæs et al., 2015)and OpenFOAM (The OpenFOAM Foundation Ltd, n.d.). Many user-friendly open sourcepackages, such as FiPy (J. E. Guyer & Warren, 2009), and many commercial packages, suchas COMSOL Multiphysics (COMSOL AB, n.d.) and Simcenter STAR-CCM+ (Siemens DigitalIndustries Software, n.d.), provide limited flexibility in the equations that users can express.Additionally, the use of commercial packages, which by nature do not perform calculationstransparently, can hinder reproducibility, which is vital to the scientific process. PyGDH(“pigged”) is a Python 2 / Python 3 (Python Software Foundation, 1991–2020) package thatis meant to be accessible to scientists who might not be specialists in scientific computing,while approaching the level of flexibility associated with writing dedicated programs tailoredto solving specific problems. The PyGDH User’s Guide provides detailed instructions forcreating numerical models, including a brief introduction to necessary command line andPython (Python Software Foundation, 1991–2020) skills, and discussions of discretization andvalidation. Note that PyGDH emphasizes flexibility and simplicity over performance, and wasnot designed for high-performance applications or models describing complex spatial domains.

97 MATHEMATICS AND COMPUTING↗

Fabrication and Preliminary Demonstration of Microwave Resonant Cavity Transducer Performance

We are investigating a microwave cavity-based transducer for in-core high-temperature fluid flow sensing in molten salt cooled reactors (MSCR) and sodium fast reactors (SFR). This sensor is a hollow metallic cylindrical cavity, which can be fabricated from stainless steel, and as such is expected to be resilient to radiation, high temperature and corrosive environment of MSCR and SFR. The principle of sensing consists of making one wall of the cylindrical cavity flexible enough so that dynamic pressure, which is proportional to fluid velocity, will cause membrane deflection. Membrane deflection causes cavity volume change, which leads to a shift in the resonant frequency. We have developed an initial design for proof-of-principle testing of the flow sensor performance in microwave K-band. A cylindrical resonator prototype was fabricated from brass for initial tests in water. The cavity size is matched to the flange of a standard WR-42 waveguide. Microwave field is coupled into the resonant cavity through a subwavelength-size aperture. A test article was developed, consisting of a piping Tee with bulkhead WR-42 microwave waveguide installed in leak-proof design. In the test article, the cylindrical cavity is positioned in the center of the pipe. A microwave waveguide circulator was installed in the setup to suppress microwave reflections at the cavity entrance. Preliminary spectral characterization of cavity spectral response was performed with microwave VNA. Applying mechanical pressure to cavity membrane showed a measurable shift in the microwave resonant frequency. We also investigate mechanical integrity of the flowmeter’s membrane through computer simulations. By calculating the stress on the plate due to deflection and, comparing the stress to the material ultimate tensile strength and yield strength, it can be estimated if the plate will fail. The stress on the plate was calculated with an analytic closed form solution model, and with COMSOL Structural Mechanics Module which does not involve any approximations. Both the analytic model and COMSOL model showed that maximum stresses on the plate, which are at the radial boundary of the plate, are three orders of magnitude smaller than the yield strength and ultimate tensile strength. This indicates that the sensor is at a low risk of mechanical failure.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

PHOENIX Electrostatic Design

PHOENIX (Portable, High-efficiency, Optimal ENergy Imaging X-rays) is a quasi-DC, electrostatic, vacuum-diode designed as a portable x-ray source with national defense and commercial applications. The patent-pending PHOENIX concept combines a megavoltage, Cockroft-Walton voltage multiplier with a Van de-Graaff electrostatic charge-storage dome to create a vacuum-diode suitable for x-ray production. Naturally this structure must minimize internal electric fields to reduce electrical breakdown while simultaneously reducing size and weight to enhance portability. In this paper we describe the optimization process and model results obtained using the COMSOL multi-physics code. We describe three models: a prototype model built to demonstrate the PHOENIX concept as part of Laboratory Directed Research and Development (LDRD) Mission Foundation Research (MFR) Phase-I , a “back-of-the-envelope” design used as a starting point for further COMSOL optimization, and finally, the optimized geometry implemented in the MFR Phase-II. In all cases compromises resulting from cost, schedule, and manufacturing constraints were taken into account as the design progressed.

42 ENGINEERING↗

Finite Element Simulation of the Acoustic Pressure Inside a Beverage Container for Non-Thermal, Ultrasound-based Pasteurization

The purpose of this effort is to investigate whether large acoustic pressure waves can be transmitted inside beverage containers to enable pasteurization. Acoustic waves are known to induce large nonlinear compressive forces and shock waves in fluids, suggesting that compression waves may be capable of damaging bacteria inside beverage containers without appreciably increasingly the temperature or altering the freshness and flavor of the beverage contents. Although a combined process such as thermosonication (e.g., sonication with heating) is likely more efficient, it is instructive to compute the acoustic pressure field distribution inside the beverage container. The COMSOL simulations used two and three-dimensional models of beverage containers placed in a water bath to compute the acoustic pressure field. A limitation of these COMSOL models is that they cannot determine the bacterial lysis efficiency, rather the models provide an indirect metric of bacterial lysis based on the magnitude of the pressure field and its distribution.

42 ENGINEERING↗

Microwave-Assisted Heating for Gasification

Compared to traditional heating, microwave-assisted heating enhances catalyst activity, selectivity and stability. These features are essential in gasification to improve the process outcome while mitigating the generation undesirable by-products such as char. To study the effect of microwaves on fluidized bed reactors, the electromagnetics module of COMSOL and MFiX are coupled using file input/output wherein external calls are made to COMSOL from MFiX. The particle bed heating rates and temperature distributions are investigateed in the case of a monodisperse bed of magnetite particles subjected to single mode transverse electric and magnetic fields in reacting and non-reacting conditions. The effect of inlet gas velocity on the heating rate of particle is also explored in this work.

Koneru, Rahul Babu↗

Thermal Analysis of the Bus Bars for the 800MeV Injection Era

❑ Introduction ❑COMSOL thermal simulations for present conditions - Results - Validation ❑COMSOL thermal simulations for PIP-II conditions ❑Proposed changes - New busbars section - Results ❑Peak temperature tables ❑Conclusions

Picchi, Silvia↗

Two-Way Coupled CFD/DEM Simulation of Electromagnetic Heating in a Fluidized Bed Reactor

A set of tools has been developed to couple a multiphase reacting computational fluid dynamics solver (MFIX) with two computational electromagnetic solvers (COMSOL and Elmer). In contrast to several recent investigations, which calculate a single static EM field prior to the CFD calculation, in the present work the electromagnetic field is periodically updated during the simulation. This enables the EM field and associated EM heating predictions to reflect the effects of evolving material properties (temperature, composition, size etc.) and the evolving distribution of material in the system. The present work uses effective medium theory (EMT) for the electromagnetic field calculations since the goal is to apply the tool for industrial scale reacting systems. Predictions using MFIX/COMSOL and MFIX/Elmer are compared to in bed thermocouple measurements from a fluidized bed experiment and recent simulation predictions and recent experiments performed at NETL. The NETL experiments also include packed bed operating conditions. Coupled MFIX/COMSOL simulations at these conditions are compared with similar simulations using only COMSOL to measure errors due to the solver coupling approach.

CFD↗

Mapping the performance of a versatile water-based condensation particle counter (vWCPC) with numerical simulation and experimental study

Accurate airborne aerosol instrumentation is required to determine the spatial distribution of ambient aerosol particles, particularly when dealing with the complex vertical profiles and horizontal variations of atmospheric aerosols. A versatile water-based condensation particle counter (vWCPC) has been developed to provide aerosol concentration measurements under various environments with the advantage of reducing the health and safety concerns associated with using butanol or other chemicals as the working fluid. However, the airborne deployment of vWCPCs is relatively limited due to the lack of characterization of vWCPC performance at reduced pressures. Given the complex combinations of operating parameters in vWCPCs, modeling studies have advantages in mapping vWCPC performance. In this work, we thoroughly investigated the performance of a laminar-flow vWCPC using COMSOL Multiphysics® simulation coupled with MATLAB™. We compared it against a modified vWCPC (vWCPC model 3789, TSI, Shoreview, MN, USA). Our simulation determined the performance of particle activation and droplet growth in the vWCPC growth tube, including the supersaturation, $D$ p,kel,0 (smallest size of particle that can be activated), $D$ p,kel,50 (particle size activated with 50 % efficiency) profile, and final growth particle size D d under wide operating temperatures, inlet pressures $P$ (30–101 kPa), and growth tube geometry (diameter $D$ and initiator length $L$ ini ). The effect of inlet pressure and conditioner temperature on vWCPC 3789 performance was also examined and compared with laboratory experiments. The COMSOL simulation result showed that increasing the temperature difference (Δ$T$) between conditioner temperature $T$ con and initiator $T$ ini will reduce $D$ p,kel,0 and the cut-off size $D$ p,kel,50 of the vWCPC. In addition, lowering the temperature midpoint ($T$ mid = $\frac{T_{con}+T_{ini}}{2}$ increases the supersaturation and slightly decreases the $D$ p,kel . The droplet size at the end of the growth tube is not significantly dependent on raising or lowering the temperature midpoint but significantly decreases at reduced inlet pressure, which indirectly alters the vWCPC empirical cut-off size. Our study shows that the current simulated growth tube geometry ($D$=6.3 mm and $L$ ini =30 mm) is an optimized choice for current vWCPC flow and temperature settings. The current simulation can more realistically represent the $D$ p,kel for 7 nm vWCPC and also achieved good agreement with the 2 nm setting. Using the new simulation approach, we provide an optimized operation setting for the 7 nm setting. This study will guide further vWCPC performance optimization for applications requiring precise particle detection and atmospheric aerosol monitoring.

47 OTHER INSTRUMENTATION↗

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↗

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↗

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↗

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↗

My Internship at NASA

My name is Isaac Lopez and I am a junior at the University of Houston majoring in Mechanical Engineering Technology. I will be completing my first tour at the NASA-Johnson Space Center ("JSC") as a Mechanical Engineer within the Human Interfaces Branch. Throughout my tour, I was given the opportunity to work on multiple projects that have expanded my knowledge and interest in acoustics and engineering design. One of the projects I worked on at JSC consisted of doing acoustic simulation of the EVA comm. cap. While working on the comm. cap headset, my main duty consisted of simulating the acoustics of the headset to find a solution to the condensing water that can accumulate and block the acoustic tube, causing attenuation or complete loss of audio in one ear for an astronaut using the EVA. For this project, I had to create a Creo model of the comm. cap so that I would be able to import it into Comsol for acoustic simulation. I also had the opportunity to design a portable and lightweight beam degrader for the EEE Parts and Radiation team. With the help of Creo, I was able to make a CAD design and put together a small working prototype for the radiation team to demonstrate the capabilities that the beam degrader had. In addition to these projects, JSC allowed me to work closely on projects with other interns. I had the opportunity to help another intern with his acoustic diverter, intended to improve the sound quality in Node 1 of the ISS. During this project, I helped with some of the acoustic testing inside the anechoic chamber as well as helping record data during testing at the ISS mock up. During the course of my first tour, I was able to learn and continually improve on my CAD drafting skills. With each project I worked on, I acquired new ways to create and improve various designs with various constraints. Furthermore, I also had the opportunity to work with electrical engineers and learn about the electronic components that would provide control of the beam degrader functions. I was able to take the knowledge, skills and designs I had acquired during my tour at JSC (such as Finite Element Analysis) and apply them to my work on the Comsol Multiphysics software for acoustics simulation. I was also able to use it to do some stress analysis on a proposed design for a proposed Habitat structure. My time at JSC has been invaluable as it helped me grow and more fully-comprehend my major as well as spark an interest in acoustics and how much it relates to mechanical engineering. This internship has taught me that engineering is much more than a profession - it is embedded in our daily lives as a puzzle that requires critical thinking as an everyday challenge that makes life more interesting.

Lopez, Isaac↗

Extracting Damping Ratio from Dynamic Data and Numerical Solutions

There are many ways to extract damping parameters from data or models. This Technical Memorandum provides a quick reference for some of the more common approaches used in dynamics analysis. Described are six methods of extracting damping from data: the half-power method, logarithmic decrement (decay rate) method, an autocorrelation/power spectral density fitting method, a frequency response fitting method, a random decrement fitting method, and a newly developed half-quadratic gain method. Additionally, state-space models and finite element method modeling tools, such as COMSOL Multiphysics (COMSOL), provide a theoretical damping via complex frequency. Each method has its advantages which are briefly noted. There are also likely many other advanced techniques in extracting damping within the operational modal analysis discipline, where an input excitation is unknown; however, these approaches discussed here are objective, direct, and can be implemented in a consistent manner.

Casiano, M. J.↗

Virtual Design of a 4-Bed Molecular Sieve for Exploration

Simulations of six new 4-Bed Molecular Sieve configurations have been performed using a COMSOL (COMSOL Multiphysics - commercial software) model. The preliminary results show that reductions in desiccant bed size and sorbent bed size when compared to the International Space Station configuration are feasible while still yielding a process that handles at least 4.0 kilograms a day CO2. The results also show that changes to the CO2 sorbent are likewise feasible. Decreasing the bed sizes was found to have very little negative effect on the adsorption process; breakthrough of CO2 in the sorbent bed was observed for two of the configurations, but a small degree of CO2 breakthrough is acceptable, and water breakthrough in the desiccant beds was not observed. Both configurations for which CO2 breakthrough was observed still yield relatively high CO2 efficiency, and future investigations will focus on bed size in order to find the optimum configuration.

Giesy, Timothy J.↗

Oxidation of Silicon Particles Suspended in Mullite Matrix With Dual-Mode Vacancy Diffusion of Oxygen

Mullite bond coats filled with silicon particles have recently been proposed for the next generation of environmental barrier coatings (EBCs). The sacrificial oxidation of the silicon particles eliminates some of the environmental durability issues associated with traditional EBC designs. To investigate this novel design approach, a mathematical formulation and resulting numerical method have been developed for simulating the oxidation of silicon particles embedded in a mullite matrix. The focus is on the oxidation of the silicon in a dry oxygen environment where the oxidant is presumably oxygen. Oxygen transport from the exterior environment through the mullite is assumed to occur via two lattice defect diffusion mechanisms. As a result, the rate of oxygen transport is dictated by the self-diffusivity of the lattice defects. The numerical method is implemented using the COMSOL Multiphysics® program (COMSOL, Inc.). The method is applied to simulate the oxidation of silicon particles in a mullite pellet in recent oxidation experiments. Numerical results compare favorably with the oxidation behavior observed in the oxidation experiments. The numerical method is also used to investigate the effects of various features of the mullite-silicon particle system on the oxidation behavior.

silicon oxidation↗