Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Understanding the impact of the gas diffusion layer structure on catalyst utilization in the PEM water electrolyzer

A multiphysics half-cell model of a polymer electrolyte membrane water electrolyzer (PEMWE) was developed to probe impacts of the detailed 3-dimensional pore structure of the gas diffusion layer (GDL) on performance characteristics. We show that pores in the titanium GDL mesh led to significant underutilization of the catalyst layer (CL), with only 45% of the catalyst effectively utilized. This contradicts the assumption of uniform electron flow across the CL, as shown in graphical abstract (a), as near-zero current was observed near GDL pore regions and the current distribution in CL was influenced by GDL structure, as shown in graphical abstract (b). Instead, oxygen generation was primarily concentrated under the solid titanium regions, diffusing out around the pore walls. High current density peaks were also noted at the GDL-catalyst contact, correlating with degradation hotspots that were directly observed in companion experiments, as shown in graphical abstract (c). Collectively, these findings point to the critical importance of the heterogeneous GDL porous architecture not only for PEMWE efficiency but also for uneven degradation of the CL.

COMSOL↗
DOE OSTI2025

Spatially Patterned Architectures to Modulate CO 2 Reduction Cascade Catalysis Kinetics

Electrochemical CO 2 reduction using renewable sources of electrical energy holds promise for converting CO 2 into fuels and chemicals. The complex interactions among chemical/electrochemical reactions and mass transport make it difficult to analyze the effect of an individual process on electrode performance based only on experimental methods. Here, we developed a generalized steady-state simulation to describe an electrode surface in which sequential cascade catalysts are patterned in a periodic trench design. If appropriately constructed, this trench geometry is hypothesized to be able to yield a higher net current density for a CO 2 reduction (CO 2 R) cascade reaction. We have used realistic experimental reaction kinetics to investigate the role of trench geometry in mass transport, local microenvironments, and selectivity for a model CO 2 R cascade reaction. The model considers local concentration gradients of bicarbonate species at quasi-equilibrium and catalytic surface reactions based on concentration-dependent Butler–Volmer kinetics. Our results suggest that varying the spatial distribution of active sites plays a significant role in facilitating effective mass transport between active sites, modulating selectivity for the cascade reaction, and enhancing the yield of desirable cascade products. Moreover, we observe that this trench geometry significantly alters the cascade reaction rate by affecting the local pH, which can cause inadvertent depletion of available aqueous CO 2 to limit the CO 2 R cascade kinetics and modest suppression of the hydrogen evolution reaction (HER). The results highlight the trade-offs between mass transport, pH, and reaction kinetics that become apparent only when considering the coupled physics of all processes at the electrode surface. Here, this model can thus serve as a primary tool to build more selective and efficient patterned architectures for the CO 2 R cascade catalysis.

CO2 reduction↗
DOE OSTI2022

Spatiotemporally resolved measurements of electric field around a piezoelectric transformer using electric-field induced second harmonic (E-FISH) generation

When a piezoelectric transformer (PT) is actuated at its second harmonic frequency by a low input voltage, the generated electric field at the distal end can be sufficient to breakdown the surrounding gas, making them attractive power sources for non-equilibrium plasma generation. Understanding the potential and electric fields produced in the surrounding medium by the PT is important for effectively designing and using PT plasma devices. In this work, the spatiotemporally resolved characteristics of the electric field generated by a PT operating in open air have been investigated using the femtosecond electric field-induced second harmonic generation (E-FISH) method. Here, electric field components were determined by simultaneously conducting E-FISH measurements with the incident laser polarized in two orthogonal directions relative to the PT crystal. Results of this work demonstrate the spatial distribution of electric field around the PT’s output distal end and how it evolves as a function of time. Notably, the strongest electric field appears on the face of the PT’s distal surface, near the top and bottom edges and decreases by approximately 70% over 3 mm. The time delay between the PT’s input voltage and measured electric field indicates that there is an about 0.45 phase difference between the PT’s input voltage and output signal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Comprehensive model for evaluating voltage losses and performance improvements in thin-film photovoltaic devices

Progress of state-of-the-art and next-generation thin-film photovoltaic devices is often stymied by open-circuit voltage (𝑉 oc ) that is significantly lower than theoretical and practical limits. Yet, effectively diagnosing the primary sources of voltage loss remains challenging. Herein, a sequence of device-level characterization techniques and simulations are employed to identify and rank loss mechanisms. For the research-based Cd⁡(Se,Te) device under study, most of the loss was at the front semiconductor heterointerface due to a clifflike conduction-band offset that lowered the recombination activation energy. Additional losses due to band tails were quantified by photoluminescence analysis. The latter provided the absorption coefficient and activation energy reduction associated with band tails as inputs to device models. Simulations showed that alleviating front-interface issues would improve 𝑉 oc , but it would then be limited by bulk recombination. Further improvement of the bulk would then lead to back-contact limitations. Reducing band tails is beneficial in any circumstance. In conclusion, this analysis provides guidance for reaching toward the radiative 𝑉 oc limit.

14 SOLAR ENERGY↗
DOE OSTI2026

Design and development of optical modules for the BUTTON-30 detector

BUTTON-30 is a neutrino detector demonstrator located in the STFC Boulby underground facility in the north-east of England. The main goal of the project is to deploy and test the performance of the gadolinium-loaded water-based liquid scintillator for neutrino detection in an underground environment. This will pave the way for a future large-volume neutrino observatory that can also perform remote monitoring of nuclear reactors for nonproliferation. This paper describes the design and construction of the watertight optical modules of the experiment.

COMSOL simulations↗
DOE OSTI2021

Scalable Technologies Achieving Risk-Informed Condition-Based Predictive Maintenance Enhancing the Economic Performance of Operating Nuclear Power Plants

The primary objective of the research presented in this report is to develop scalable technologies that are deployable across plant assets and across the nuclear fleet to achieve risk-informed predictive maintenance (PdM) strategies at commercial nuclear power plants (NPPs). Over the years, the nuclear fleet has relied on labor-intensive and time-consuming preventive maintenance (PM) programs, driving up operation and maintenance (O&M) costs to achieve high capacity factors. A well-constructed risk-informed PdM approach for an identified plant asset has been developed in this research, taking advantage of advancements in data analytics, machine learning (ML), artificial intelligence (AI), physics-informed modeling, and visualization. These technologies would allow commercial NPPs to reliably transition from current labor-intensive PM programs to a technology driven PdM program, eliminating unnecessary O&M costs. The work presented in the report is being developed as part of a collaborative research effort between Idaho National Laboratory and Public Service Enterprise Group Nuclear, LLC. This report (1) reflects the results of work by LWRS Program researchers with PSEG, Nuclear LLC-owned Salem and Hope Creek Nuclear Power Plants; (2) presents utilization of circulating water system (CWS) heterogeneous data and fault modes from both the Salem and Hope Creek nuclear power plant sites to develop salient fault signatures associated with each fault mode; (3) describes the integration of component-level predictive models into a robust system-level model enabled by the federated-transfer learning; (4) describes the development of physics-informed model of circulating water pump and motor; (5) develops a scalable risk and economic model; and (6) outlines the development of a user-centric visualization application. The outcomes presented in this report lays the foundation and provides a much-needed technical basis to focus on explainability and trustworthiness of ML and AI-based technologies, as part of future research.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Thermal Calculations to Determine Maximum Zircaloy Temperatures During Reactor Operations for Campaign 1 Fuel Assemblies

The H-Canyon Facility plans on processing bundles from Non-Aluminum Spent Nuclear Fuel (NASNF) Risk Group 1 (UO 2 core fuels). Uranium-zirconium (UZr) alloys and zirconium clad UO 2 fuel containing UZr intermetallic compound have the potential to cause explosions during dissolution. UZr intermetallic typically only forms when uranium and zirconium (or zirconium alloy) cladding are in contact at high temperatures. Thermal modeling was performed to determine the maximum cladding temperature for Campaign 1 Zircoloy cladded fuel assembles during reactor operation to help determine if the UZr intermetallic compound could have formed during reactor operations. The temperature of 350°C was chosen as a conservative minimum value for the possible formation of intermetallic UZr compound. The modeling approach, inputs, and results are outlined in this document.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Hybrid Modeling of a Circulating Water Pump Motor

To achieve high capacity factors, the nuclear fleet has relied on labor-intensive and time-consuming operation and preventive maintenance programs for plant systems. Manually performed inspection, calibration, testing, and maintenance of plant assets at periodic frequencies, along with the time-based replacement of assets irrespective of condition, have resulted in a costly, labor-centric business model. Fortunately, there are technologies that can eliminate unnecessary preventive maintenance activities by deploying risk-informed predictive maintenance, enabling the transition to a technology-centric business model.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Robust Free-Space Optical Communication Utilizing Polarization for the Advancement of Quantum Communication

Free-space optical (FSO) communication can be subject to various types of distortion and loss as the signal propagates through non-uniform media. In experiment and simulation, we demonstrate that the state of polarization and degree of polarization of light passed though underwater bubbles, causing turbulence, is preserved. Our experimental setup serves as an efficient, low cost alternative approach to long distance atmospheric or underwater testing. We compare our experimental results with those of simulations, in which we model underwater bubbles, and separately, atmospheric turbulence. Our findings suggest potential improvements in polarization based FSO communication schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Reaction Diffusion Modelling of 3D Pillar Electrodes in Single-Catalyst CO 2 Reduction Cascades

Effective electrochemical CO 2 reduction to liquid fuels requires that the local catalytic environment facilitates the desired reactivity, yet a microscopic understanding of this environment is difficult to achieve from experiment alone. In this work, a 3D reaction-diffusion model was developed to explore the effects of electrode surface area and local geometry on the performance of a heterogeneous catalyst that performs a two-step CO 2 reduction cascade reaction to CO and then CH 3 OH under aqueous conditions. Kinetic parameters for the model were inspired by experimental results using a cobalt phthalocyanine (CoPc) catalyst. Three-dimensional architectures composed of arrays of square pillars with varying dimensions and either smooth or periodically modulated surfaces were tested, revealing the extent to which geometry modulates the performance of the cascade reactions. Although structural variations modulate local concentration gradients, we find that electrochemically active surface area predominantly governs the overall cascade reaction. Moreover, the results suggest that supersaturation of CO, with concentrations up to ten-fold higher than the equilibrium solubility limit, might be critical for more efficient conversion to CH 3 OH. For any given geometry, the spatially averaged ratio of [CO] to [CO 2 ] is dictated by the electrochemically active surface area and determines the yield of CH 3 OH. For a fixed surface area, geometries that spatially confine the electrolyte yield moderate local [CO] to [CO 2 ] ratios within small volumes. In contrast, less confining geometries result in a broader distribution of local ratios spread over larger volumes, with both configurations yielding the same spatially averaged [CO] to [CO 2 ] ratio. These insights provide valuable design principles—highlighting the critical importance of surface area and possibly CO supersaturation—for engineering advanced electrode architectures that leverage intermediate trapping and CO supersaturation to enhance overall performance in tandem CO 2 reduction systems.

COMSOL↗
DOE OSTI2022

Simulation and Modeling of Optical Properties of U, Th, Pb, and Co Nanoparticles of Interest to Nuclear Security Using Finite Element Analysis

In this work, the optical characteristics of uranium (U), lead (Pb), cobalt (Co), and thorium (Th) nanoparticles are fashioned and simulated employing the finite element analysis (FEA) approach concerning multiple particle sizes. Applying finite element analysis, it was found that the simulated absorption peaks of electronic excitations of nuclear nanoparticles are red-shifted from 365 nm to 555 nm for U; from 355 nm to 550 nm for Pb; from 415 nm to 610 nm for Co; and from 350 nm to 540 nm for Th, comparing expanding particle sizes from 60 nm to 100 nm (except for Co, which varied from 70 nm to 100 nm). The FEA-simulated optical band gap energies and far-field radiation patterns were also obtained for nuclear materials. The simulation approach in this research enables the prediction of optical properties and design of nuclear materials before manufacture for nuclear security applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modeling the First Hydrogen Direct Reduction Pilot Reactor for Ironmaking in the USA Using Finite Element Analysis and Its Validation Using Pilot Plant Trial Data

Direct reduction of hematite pellets with hydrogen (H 2 ) was used to produce directly reduced iron (DRI) in a pilot scale reactor at a pellet feed rate of 21.4 kg/h. At a steady state, operational parameters of the pilot plant (gas recycling rate and inlet temperature) along with key reactor output parameters, the pellet metallization, and the internal temperature profile of the reactor were reported for two scenarios with high recycle and low recycle rate of H 2 . Scenario 1, with a high recycle rate of 400 L/min H 2 along with external heating of 870 °C, gave an average metallization of 91.8%, while Scenario 2, with low recycle rate of 100 L/min H 2 and external heating of 850 °C gave a metallization of 67.8% due to the higher moles of H 2 available for reduction and the external energy required for the endothermic reduction reaction in Scenario 1 as compared with Scenario 2. Finite element analysis was used to build a model of the shaft reactor, which was validated against the metallization and internal temperature profile data. The average metallization values predicted by the model were very close to the metallization values obtained from the pilot plant samples, with 90.9% average metallization for Scenario 1 and 65.6% average metallization for Scenario 2. The internal temperature profiles in the lower region of the reactor obtained from the model were very close to these pilot plant data, with a maximum difference of 52.7 °C and 67.6 °C for Scenarios 1 and 2, respectively. The pilot plant reactor model was used extensively in the commissioning of the pilot plant and to predict the startup outcomes for a given set of operating parameters.

08 HYDROGEN↗
NASA NTRSDate not supplied

Preliminary Considerations for Microwave Consolidation/Sintering of Lunar Regolith Simulant

As NASA prepares to establish permanent habitats on the Moon, a significant first step is to be able to land multiple times in the same area. As was consistently shown during the Apollo program, the very fine granular structure of the lunar regolith (the Moon’s “soil”)poses significant physical and health challenges [1]. One of the most concerning is the hyper-velocity lunar surface ejecta that results from the engine exhaust that exits the rocket as it lands and takes off [2]. It has been determined that to mitigate this, the regolith must be consolidated.

microwave↗
NASA NTRSDate not supplied

Comparing Acoustic Prediction Methods for Additively Manufactured Porous Strutures

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

impedance↗
NASA NTRSDate not supplied

Assessment of Acoustic Behavior for Perforate-Over-Large-Cell Liners

This study explores the effects of increasing the cell size for large-cell acoustic liners. Tests are conducted in the NASA Langley Grazing Flow Impedance Tube to evaluate liners with increasingly larger cell dimensions (up to 2” x 4” cross-section). Conventional impedance eduction confirms that liners with 2” x 3” cells (or larger) must be evaluated using nonlocally reacting assumptions. In addition, due to the sound propagation within cavities, the liners must be modeled using higher fidelity techniques. Thus, grazing flow duct and acoustic liner are modeled simultaneously using finite element methods. The facesheet is modeled using a transfer impedance, while the rest of the domain is modeled using the convected Helmholtz equation. The acoustic pressures predicted are shown to compare favorably with those measured in the Grazing Flow Impedance Tube.

large diameter↗