Search NASA⌕ Search

SEARCH · Search NASA

Results for “COMSOL simulations”

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

Electromagnetic control of heat transport within a rectangular channel filled with flowing liquid metal

The behavior of free-surface, liquid-metal flows exposed to both magnetic fields and an injected electric current is investigated via experiment and numerical simulations. The purpose of this paper is to provide an experimental and theoretical proof-of-concept for enhanced thermal mixing within fast-flowing, free-surface, liquid-metal plasma facing components that could be used in next-generation fusion reactors. The enhanced hydrodynamic and thermal mixing induced by non-uniform current density near the electrodes appears to improve heat transfer through the thickness of the flowing metal. Also, the outflow heat flux profile is strongly affected by the impact of the J × B forces on flow velocity. The experimental results are compared to COMSOL simulations in order to lay the groundwork for future liquid-metal research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental and simulated heating in nanocrystalline diamond charge exchange injection films using an electron beam to mimic SNS conditions

Nanocrystalline diamond (NCD) foils have been in use at the Spallation Neutron Source as primary stripper foils for charge exchange injection since the machine’s inception. NCD’s low atomic number leads to reduced beam scattering while still being easy to handle, as opposed to even lighter elements. NCD foils also have the benefit of being free-standing, and rigid when grown under the right conditions to minimize residual stresses. This study overviews a method to mimic the Spallation Neutron Source (SNS) beam and characterize thermally driven failure of NCD foils. To do this, a foil test stand was developed with in situ diagnostics that tracks signs of foil sublimation and thinning. The foil test stand’s electron beam is equated to the SNS beamline by comparing peak deposited energy densities for both beams. Post-mortem testing is also conducted to help elucidate changes the NCD foils undergo during exposure to the beam. Further, a COMSOL simulation was also developed as a method to assess potential future changes to the NCD film and predict how changes to beam conditions will effect foil temperatures. This method allows for examination of the NCD foil’s thermal limits and can be used to assess future changes to SNS stripper foils.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Material Property Estimation in Thin Battery Components Using Guided Wave Measurement, Experimental Dispersion Curve Extraction and Finite Element Modeling

At NASA, we are investigating nondestructive evaluation (NDE) and structural health monitoring (SHM) techniques to detect precursors of thermal runaway failure in lithium metal based lithium ion batteries (LIB). The approach is centered on computational simulation models to guide inspection and aid in interpretation of results. Since lithium metal LIBs have combinations of solid and fluid-filled porous materials, obtaining accurate material properties is both challenging and critical for successful simulation of battery inspection. To this end, we investigated a multi-verification approach for material characterization of thin battery components. First, a laser Doppler vibrometer (LDV) was used to measure guided wave fields in thin (microns-thick) battery components subject to broadband excitation. The time-space wavefield data was converted to frequency-wave number data to extract guided wave dispersion curves. A data visualization and post processing graphics user interface (GUI) was developed at NASA to aid the data exploration and analysis. Due to the thinness of the samples, low frequency-thickness-product plate wave approximations allowed for the calculation of closed-form solutions for material elastic property estimation. These approximations were then verified by calculating the Lamb wave solutions using the previously obtained material properties. Finally, the estimated elastic material properties were implemented in a COMSOL simulation model, and dispersion curves were extracted from simulation results. The dispersion curves and derived material properties were then compared to the analysis results from the experimental data. These comparisons informed on the accuracy of the measured material properties and helped demonstrate the accuracy of the finite element analysis (FEA) computational models. This assessment will prove vital when we start simulating more complex multi-layer components and poroelastic media. This paper gives a brief background of the problem space, outlines the workflow for data analysis and verification, shows results from the workflow, and gives an overview of future plans for simulation of lithium metal LIB inspection.

Peter Juarez↗

Conception and realization of a parallel-plate free-air ionization chamber for the absolute dosimetry of an ultrasoft X-ray beam

We report the design of a millimeter-sized parallel plate free-air ionization chamber (IC) aimed at determining the absolute air kerma rate of an ultra-soft X-ray beam (E = 1.5 keV). The size of the IC was determined so that the measurement volume satisfies the condition of charged-particle equilibrium. The correction factors necessary to properly measure the absolute kerma using the IC have been established. Particular attention was given to the determination of the effective mean energy for the 1.5 keV photons using the PENELOPE code. Other correction factors were determined by means of computer simulation (COMSOL™and FLUKA). Measurements of air kerma rates under specific operating parameters of the lab-bench X-ray source have been performed at various distances from that source and compared to Monte Carlo calculations. We show that the developed ionization chamber makes it possible to determine accurate photon fluence rates in routine work and will constitute substantial time-savings for future radiobiological experiments based on the use of ultra-soft X-rays.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

COMSOL-based Multiphysics Simulations to Support HFIR s Conversion to LEU Fuel

In this paper, development of at least one form of the COMSOL-based modeling framework for the HFIR is presented, key simulation steps are identified and several milestones achieved towards a coupled multi-physics capability are highlighted. COMSOL-based multi-physics simulation capability is able to answer the need for predictive 3D simulations of HFIR s involute plate and channels. Step-by-step development and analyses of the COMSOL models for the single and multi-channels will lead towards the desired full-core simulation capability for the HFIR. With very few experiments planned to support the conversion process, these 3D simulations will become the basis for the nuclear safety analysis of the HFIR s LEU fuel core.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An SECM-Based Spot Analysis for Redoxmer-Electrode Kinetics: Identifying Redox Asymmetries on Model Graphitic Carbon Interfaces

The fundamental process in non-aqueous redox flow battery (NRFB) operation revolves around electron transfer (ET) between a current collector electrode and redox-active organic molecules (redoxmers) in solution. Here, we present an approach utilizing scanning electrochemical microscopy (SECM) to evaluate interfacial ET kinetics between redoxmers and various electrode materials of interest at desired locations. This spot-analysis method relies on the measurement of heterogeneous electron transfer rate constants (k f or k b ) as a function of applied potential (E-E 0 '). As demonstrated by COMSOL simulations, this method enables the quantification of Butler-Volmer kinetic parameters, the standard heterogeneous rate constant, k 0 , and the transfer coefficient, α. Our method enabled the identification of inherent asymmetries in the ET kinetics arising during the reduction of ferrocene-based redoxmers, compared to their oxidation which displayed faster rate constants. Similar behavior was observed on a wide variety of carbon electrodes such as multi-layer graphene, highly ordered pyrolytic graphite, glassy carbon, and chemical vapor deposition-grown graphite films. However, aqueous systems and Pt do not exhibit such kinetic effects. Our analysis suggests that differential adsorption of the redoxmers is insufficient to account for our observations. Displaying a greater versatility than conventional electroanalytical methods, we demonstrate the operation of our spot analysis at concentrations up to 100 mM of redoxmer over graphite films. Looking forward, our method can be used to assess non-idealities in a variety of redoxmer/electrode/solvent systems with quantitative evaluation of kinetics for applications in redox-flow battery research.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A compact light shield for the Psyche gamma-ray spectrometer

A compact, low-mass light shield has been designed for the Psyche space mission high-purity germanium gamma-ray spectrometer. The spectrometer will be vented to space, which can allow solar photons to enter the cryostat and cause charge-induced noise. The light shield was designed to block more than 99.9% of solar photons while maintaining high molecular conductance to space vacuum. COMSOL simulations and experiments were performed to characterize light shield photon attenuation and molecular conductance. Finally, the light shield was integrated onto the gamma-ray spectrometer and photon suppression and vacuum performance were demonstrated.

(16) Psyche↗

Finite element modeling of electropolishing of radio frequency dipole Nb crab cavity in hydrofluoric-sulfuric acid electrolyte

The superior performance of niobium superconducting radio frequency (SRF) cavities is achieved through meticulous surface treatment, notably via chemical electropolishing, ensuring exceptionally smooth surfaces. While this technique has been extensively employed for cylindrically symmetric structures, addressing more intricate geometries poses significant challenges in achieving uniform polishing and controlled material removal, especially when moving away from retractable cathodes. This presents a multifaceted electrochemical, thermal, and fluid dynamics conundrum. A prime example is the 197 MHz radio frequency dipole (RFD) crabbing cavity proposed for the Electron Ion Collider (EIC) project, exemplifying such complex structures. Our groundbreaking work integrates the localized oxide thickness variation, considering its impact on current distribution and Joule heating, within the framework of multi-physics modeling using the COMSOL® simulation suite. This approach was applied to a comprehensive model of the RFD cavity, allowing us to investigate optimal external cooling water flow conditions necessary for achieving desirable outcomes. In conclusion, this illustrates the increasing utility of such multi-physics codes to greatly aid in designing solutions to challenging processing requirements for increasingly complex accelerator cavities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

High-Temperature 2D Optical Relaxation Visualizes Enhanced Oxygen Exchange Kinetics at Metal-Mixed Conducting Oxide Interfaces

Solid-state heterointerfaces are of interest for emergent local behavior that is distinct from either bulk parent compound. One technologically relevant example is the case of mixed ionic/electronic conductor (MIEC)–metal interfaces, which play an important role in electrochemistry. Metal–MIEC composite electrodes can demonstrate improved catalytic activity vs single-phase MIECs, improving fuel cell efficiency. Similarly, MIEC surface reaction kinetics are often evaluated using techniques that place metal current collectors in contact with the surface under evaluation, potentially altering the response vs the native surface. Here, techniques enabling direct and local in situ observation of the behavior at and around such heterointerfaces are needed. Here, we develop a spatially resolved optical transmission relaxation (2D-OTR) method providing continuous evaluation of local, high-temperature, controlled atmosphere defect kinetics across a ~1 cm2 sample area simultaneously in a contact-free manner. We apply it to observe the spatial variance of oxygen incorporation and evolution rates at ~525–620 °C, in response to step changes in oxygen partial pressure, on MIEC SrTi 0.65 Fe 0.35 O 3–x films as a function of distance from porous Pt and Au layers. Using this model geometry, we find significant enhancements in kinetics adjacent to the metals that decay over a few millimeter distance. To extract kinetic parameters, we fit the short-term optical data (initial portion of relaxations) with an exponential decay function appropriate for surface-exchange-limited kinetics, yielding apparent surface exchange coefficients (k chem ) with spatial resolution, decreasing with distance from the metal. To understand the kinetic processes governing the complete (long-term) optical relaxations, we performed COMSOL simulations, which demonstrated that a combination of laterally varying k chem and in-plane diffusion controls the observed kinetics over the full time range. Further support for spatially varying kchem comes from demonstrations of changing surface and bulk chemistry vs distance from the metal–MIEC interface, by X-ray photoelectron and optical absorption spectroscopies, respectively. Although microporous Pt and Au are not excellent electrodes in isolation, both metals exert a synergistic effect on the oxygen surface exchange rate in the presence of the mixed conducting film.

36 MATERIALS SCIENCE↗

Integration of CeO 2 -Based Memristor with Vertically Aligned Nanocomposite Thin Film: Enabling Selective Conductive Filament Formation for High-Performance Electronic Synapses

The CeO 2 -based memristor has attracted significant attention due to its intrinsic resistive switching (RS) properties, large on/off ratio, and great plasticity, making it a promising candidate for artificial synapses. However, significant challenges such as high power consumption and poor device reliability hinder its broad application in neuromorphic microchips. To tackle these issues, in this work, we design a novel bilayer (BL) memristor by integrating a CeO 2 -based memristor with a Co-CeO 2 vertically aligned nanocomposite (VAN) layer and compare it with the single layer (SL) memristor. Preliminary electrical testing reveals that the BL memristor offers a reduced set/reset voltage (~67% lower), a higher on/off ratio (~5 × 10 2 ), enhanced device reliability, and improved device-to-device variation compared to the SL memristor. Insight from COMSOL simulation, coupled with microstructural analysis, provides a comprehensive elucidation on how the VAN layer facilitates the selective conductive filament (CF) formation. Subsequently, the plasticity of the BL memristor is evaluated through long-term potentiation/depression (LTP/LTD), paired-pulse facilitation (PPF), and spike-time-dependent plasticity (STDP). The spiking neural network (SNN) built upon the BL memristor achieves remarkable accuracy (~94%) after only 12 iterations, underscoring its potential for high-performance neural networks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bias-field controlled phasing and power combination of gyromagnetic nonlinear transmission lines

Gyromagnetic Nonlinear Transmission Lines (NLTLs) generate microwaves through the damped gyromagnetic precession of the magnetic moments in ferrimagnetic material, and are thus utilized as compact, solid-state, frequency agile, high power microwave (HPM) sources. The output frequency of a NLTL can be adjusted by control of the externally applied bias field and incident voltage pulse without physical alteration to the structure of the device. This property provides a frequency tuning capability not seen in many conventional e-beam based HPM sources. The NLTLs developed and tested are mesoband sources capable of generating MW power levels in the L, S, and C bands of the microwave spectrum. For an individual NLTL the output power at a given frequency is determined by several factors including the intrinsic properties of the ferrimagnetic material and the transmission line structure. Hence, if higher power levels are to be achieved, it is necessary to combine the outputs of multiple NLTLs. This can be accomplished in free space using antennas or in a transmission line via a power combiner. Using a bias-field controlled delay, a transient, high voltage, coaxial, three port, power combiner was designed and tested. Experimental results are compared with the results of a transient COMSOL simulation to evaluate combiner performance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Transport and photo-conduction in carbon nanotube fibers

We have characterized the conductivity of carbon nanotube (CNT) fibers enriched in semiconducting species as a function of temperature and pulsed laser irradiation of 266 nm wavelength. At high temperatures, the response approaches Arrhenius law behavior, but from room temperature down to 4.2 K, the response can be quantitatively framed within the predictions of the fluctuation induced tunneling model, with tunneling occurring between the inner fibrils (bundles) of the samples and/or the elementary CNTs constituting the fibers. Fiber conductivity is enhanced under laser irradiation, and analysis of the resulting data confirms the (exponential) dependence of the potential barrier upon temperature as expected from the fluctuation induced tunneling model. A thermal map of the laser-irradiated fibers is also obtained via COMSOL simulations in order to rule out bare heating phenomena as the background of our experiments.

36 MATERIALS SCIENCE↗

Spatially programmed alignment and actuation in printed liquid crystal elastomers

Liquid crystal elastomers (LCEs) exhibit reversible shape morphing behavior when cycled above their nematic-to-isotropic transition temperature. During extrusion-based 3D printing, LCE inks are subjected to coupled shear and extensional flows that can be harnessed to spatially control the alignment of their nematic director along prescribed print paths. Here, we combine experiment and modeling to elucidate the effects of ink composition, nozzle geometry, and printing parameters on director alignment. From rheological measurements, we quantify the dimensionless Weissenberg number ( Wi ) for the flow field each ink experiences as a function of printing conditions and demonstrate that Wi is a strong predictor of LCE alignment. We find that director alignment in LCE filaments printed through a tapered nozzle varies radially when Wi < 1, while it is uniform when Wi ≫ 1. Based on COMSOL simulations and in operando X-ray measurements, we show that LCE inks printed through nozzles with an internal hyperbolic geometry exhibit a more uniform director alignment for a given Wi compared to those through tapered nozzles. Concomitantly, the stiffness along the print direction and actuation strain of printed LCEs increases substantially under such conditions. By varying Wi during printing through adjusting the flow rate “on the fly”, LCE architectures with uniform composition, yet locally encoded shape morphing transitions can be realized.

36 MATERIALS SCIENCE↗

Experimental evidence of lower hybrid wave scattering in Alcator C-Mod due to scrape off layer density fluctuations

Novel experimental measurements of the lower hybrid (LH) wave electric field vector, E LH , obtained in Alcator C-Mod using dynamic Stark spectroscopy are presented. Three key results were obtained by comparison with 3D full-wave COMSOL simulations using the cold plasma dielectric tensor and reflectometry-measured density profiles. (1) As a function of space and LH power, the magnitude of E LH was found to agree with the full-wave simulation result. (2) The direction of E LH was found to have a substantial poloidal component and is in strong disagreement with the nearly radial full-wave simulation result. (3) The addition of scrape off layer density fluctuations to the density profile implemented in the full-wave simulations can be used to explain the E LH direction discrepancy. It is concluded that diffraction and scattering by turbulence-driven scrape off layer density fluctuations is the responsible mechanism behind the modification in the direction of E LH . Furthermore, these results further support the previously reported proposition that edge turbulence is an important mechanism behind the strong anomalous drop in LH current drive efficiency for diverted plasmas in the high-density regime.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

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↗

Thickness Dependency of Thin Film Samaria Doped Ceria for Oxygen Sensing

High temperature oxygen sensors are widely used for exhaust gas monitoring in automobiles. This particular study explores the use of thin film single crystalline samaria doped ceria as the oxygen sensing material. Desired signal to noise ratio can be achieved in a material system with high conductivity. From previous studies it is established that 6 atomic percent samarium doping is the optimum concentration for thin film samaria doped ceria to achieve high ionic conductivity. In this study, the conductivity of the 6 atomic percent samaria doped ceria thin film is measured as a function of the sensing film thickness. Hysteresis and dynamic response of this sensing platform is tested for a range of oxygen pressures from 0.001 Torr to 100 Torr for temperatures above 673 K. An attempt has been made to understand the physics behind the thickness dependent conductivity behavior of this sensing platform by developing a hypothetical operating model and through COMSOL simulations. This study can be used to identify the parameters required to construct a fast, reliable and compact high temperature oxygen sensor.

33 ADVANCED PROPULSION SYSTEMS↗

Investigating the dielectric constant of barium titanate in a polymer-matrix nanocomposite

Barium titanate (BTO) is a ferroelectric material used in capacitors because of its high bulk dielectric constant. However, the impact of the size of BTO on its dielectric constant is not yet fully understood and is highly contested. Here, we present an investigation into the dielectric constant of BTO nanoparticles with diameters ranging between 50 and 500 nm. BTO nanoparticles were incorporated into acrylonitrile butadiene styrene and injection molded into parallel plate capacitors, which were used to determine nanocomposite dielectric constants. The dielectric constants of BTO nanoparticles were obtained by combining experimental measurements with computational results from COMSOL simulations of ABS-matrix nanocomposites containing BTO. The dielectric constant of BTO was observed to be relatively constant at nanoparticle diameters as small as 200 nm but sharply declined at smaller nanoparticle sizes. Overall, these results will be useful in the development of improved energy storage and power conditioning systems utilizing BTO nanoparticles.

36 MATERIALS SCIENCE↗

Composite Structured Surfaces for Durable Dropwise Condensation

Condensation heat transfer is a vital process for a plethora of industrial application. Dropwise condensation of steam on hydrophobic polymers, typified by the formation of discrete water droplets which shed and clear the surface for re-nucleation, has a 10X higher heat transfer coefficient when compared to filmwise condensation on bare metallic surfaces. However, 8 decades of research have conclusively demonstrated a lack of durability of the hydrophobic coatings during dropwise condensation to quantify industrial implementation. They key bottleneck is the coating conundrum, which states that thick hydrophobic coatings (> 10 µm) are required for higher durability, while thin hydrophobic coatings (< 1 µm) are required for low parasitic thermal resistance due to the low intrinsic thermal conductivity of polymeric materials (k ~ 0.1 W/m·K). Here, we attempt to solve the conundrum by developing rationally-designed metal-polymer composite coatings to increase the effective thermal conductivity of the hydrophobic coating, thereby enabling thicker coatings for higher durability. We study nanowire, inverse opal, and sintered sphere metallic structures infilled with hydrophobic polymers. To quantify the heat transfer performance, we analyze the vapor-to-coolant and vapor-to-surface heat transfer coefficients of composite coated flat surfaces having varying wire/sphere diameters (D = 100 nm, 1 μm, 10 μm) and coating thickness (10 nm < t < 100 μm) using three-dimensional (3D), steady state, Comsol simulations. We show that analytical models such as Maxwell or effective medium theory fail to predict the composite thermal conductivity. By conducting simulation with two independent sets of boundary conditions; isothermal surfaces and convection, as well as by varying the interfacial metal fractions, we demonstrate the importance of 3D heat spreading near the coating-vapor and coating-substrate interface. To predict the wetting characteristics to ensure stable dropwise condensation, we coupled our 3D heat transfer simulations to a unified model for contact angle hysteresis on heterogeneous surfaces. Our wetting predictions identify the key tradeoff between low spreading resistance (high k, filmwise condensation) and low contact angle hysteresis (low k, dropwise condensation) on high and low surface metal fractions, respectively, and identify the optimum metallic structure, volume fraction, interfacial solid fraction, and thickness for enhanced durability. Our work provides a novel approach to utilize rational metallic nanostructuring to create composite hydrophobic coatings which promise enhanced durability.

Chang, Ho Chan↗