Search NASA⌕ Search

Engineering topics

Roy, Thomas

Publications and source records attributed to Roy, Thomas.

Topology optimization for the full-cell design of porous electrodes in electrochemical energy storage devices

In this paper, we introduce a density-based topology optimization framework to design porous electrodes for maximum energy storage. We simulate the full cell with a model that incorporates electronic potential, ionic potential, and electrolyte concentration. The system consists of three materials, namely pure liquid electrolyte and the porous solids of the anode and cathode, for which we determine the optimal placement. We use separate electronic potentials to model each electrode, which allows interdigitated designs. As a result, a penalization is required to ensure that the anode and cathode do not touch, i.e., causing a short circuit. We compare multiple 2D designs generated for different fixed conditions, e.g. material properties. A 3D design with complex channel and interlocked structure is also created. All optimized designs are far superior to the traditional monolithic electrode design with respect to energy storage metrics. We observe up to a 750% increase in energy storage for cases with slow effective ionic diffusion within the porous electrode.

25 ENERGY STORAGE↗

Shape Matters: Understanding the Effect of Electrode Geometry on Cell Resistance and Chemo-Mechanical Stress

Rechargeable batteries that incorporate shaped three-dimensional electrodes have been shown to have increased power and energy densities when compared to a conventional geometry, i.e. a planar cathode and anode that sandwich an electrolyte. Electrodes can be shaped to enable a higher active material loading, while keeping ion transport distances small. However, the relationship between electrical and mechanical performance of shaped electrodes remains poorly understood. Many electrode designs have been explored, where the electrodes are individually shaped or intertwined, and advances in manufacturing and shape/topology optimization have made such designs a reality. Here, we explore sinusoidal half cells and interdigitated full cells. First, we use a simple electrostatics model to understand the cell resistance as a function of shape. Here we focus on low-temperature conditions, where the electrolyte conductivity decreases relative to that of the electrode; here, LiPF 6 EC:DMC electrolyte and MnO 2 electrode are considered. Next, we use a chemo-mechanics model to examine the stress that arises due to intercalation-driven volume expansion. We show that shaped electrodes provide a significant reduction in resistance in low-temperature conditions, however, they exhibit unfavorable stress concentrations. Overall, we find that the fully interdigitated electrodes may provide the best balance with respect to this resistance-stress trade-off.

25 ENERGY STORAGE↗

EchemFEM: A Firedrake-based Python package for electrochemical transport

The transition from fossil fuels to renewable energy has brought about a rapid increase in the availability of clean electricity. However, electricity generated from sources such as wind and solar are limited to intermittent operation due to daily and seasonal variation. One solution is to utilize electrochemical devices in energy storage and electrochemical manufacturing applications, where they can harness surplus energy and decarbonize chemical industries traditionally reliant on petrochemical feedstocks. Managing the growing prevalence of renewable energy underscores the importance of developing and scaling up these technologies, which can in turn facilitate the achievement of carbon emission reduction commitments of companies and developed economies. Likewise, the electrification of transport creates an increasing need for energy-dense electrochemical energy storage devices such as batteries and supercapacitors. Naturally, simulation tools are required to assist in the design of efficient and industrial-scale electrochemical devices.

30 DIRECT ENERGY CONVERSION↗

Train small, model big: Scalable physics simulators via reduced order modeling and domain decomposition

Numerous cutting-edge scientific technologies originate at the laboratory scale, but transitioning them to practical industry applications is a formidable challenge. Traditional pilot projects at intermediate scales are costly and time-consuming. An alternative, the pilot-scale model, relies on high-fidelity numerical simulations, but even these simulations can be computationally prohibitive at larger scales. To overcome these limitations, we propose a scalable, physics-constrained reduced order model (ROM) method. The ROM identifies critical physics modes from small-scale unit components, projecting governing equations onto these modes to create a reduced model that retains essential physics details. We also employ Discontinuous Galerkin Domain Decomposition (DG-DD) to apply ROM to unit components and interfaces, enabling the construction of large-scale global systems without data at such large scales. Here this method is demonstrated on the Poisson and Stokes flow equations, showing that it can solve equations about 15–40 times faster with only ~1% relative error. Furthermore, ROM takes one order of magnitude less memory than the full order model, enabling larger scale predictions at a given memory limitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Topology Optimization for Full-Cell Electrochemical Energy Storage Device

This code is used to conduct topology optimization of a general full cell energy storage device. The anode and cathode designs are optimized simultaneously, allowing possibility for interdigitation. The charging mechanism can consist of both electric double layer charging (EDLC) typically observed in capacitors and redox reaction typically observed in batteries.

Roy, Thomas↗

Simplified Models of the Bicarbonate Buffer for Scaled Simulations of CO 2 Electrolyzers

Bicarbonate electrolytes are used in a range of chemical processes; however, resolved simulation of these electrolytes is difficult, as disparate reaction time scales lead to numerical stiffness and the formation of fine boundary layers. Based on several physically motivated approximations, we reduce the full set of chemical reactions within a bicarbonate electrolyte to a simpler subset, eliminating the numerical stiffness. We supported this simplification via a two-variable singular perturbation expansion and demonstrated that under neutral conditions (6 < pH < 10) typically found within a fluid-fed CO 2 electrolyzer, the simplifications lead to negligible error. Here, we also discuss two alternative simplifications, one valid at high pH and another valid at arbitrary pH. These simplifications reduce the condition number of the matrices resulting from spatiotemporal discretization by up to 10 orders of magnitude and enable three-dimensional (3D) simulation of CO 2 electrolyzers containing carbonate solutions.

10 SYNTHETIC FUELS↗

Scalable Reduced Order Model with Discontinuous Galerkin Domain Decomposition

scaleupROM is a scalable, physics-constrained reduced order model (ROM). It aims to provide robust, accelerated physics predictions at extrapolated scales, based on the small, component-level data. This is implemented by combining projection-based ROM with discontinuous Galerkin domain decomposition, in the framework of MFEM and libROM. It currently supports the Poisson equation and Stokes flow equation, and more work is in progress toward general, nonlinear physics systems.

Chung, Seung Whan↗

FireCat

FireCat performs multiscale simulations by coupling ab initio microkinetic simulations from CatMAP with two-dimensional continuum transport simulations from EchemFEM (Firedrake). This coupling is performed in an iterative fashion. The examples provided here are specific to electrochemical CO2 reduction to CO on Au electrodes.

Roy, Thomas↗

Design and additive manufacturing of optimized electrodes for energy storage applications

Supercapacitors exhibit fast charging/discharging ability and have attracted considerable attention within the automotive, aerospace, and telecommunication industries. Porous carbons, prized for their high electrical conductivity and high surface area, have been attractive candidates for supercapacitor electrodes. Moving to thick electrodes is one strategy to further increase energy density due to a higher volume fraction of active material. However, thick electrodes suffer from sluggish charged species transport, which is why thin electrodes are currently favored. In this work, we investigate the use of computational optimization and additive manufacturing to design and fabricate thick porous electrodes with improved performance. Electrode performance was maximized by designing their morphologies via topology optimization and printing by projection micro stereolithography (PμSL) using commercial resin (PR48). The PR48 resin was then pyrolyzed (PR48-P) to create the final conductive electrode. The optimized PR48-P electrodes exhibited 99% improvement in capacitance compared to control electrodes printed with cubic lattice morphologies. To further improve performance, we formulated a resin combining graphene oxide (GO) and trimethylolpropane triacrylate (TMPTA). Electrodes printed with 3 wt% GO in TMPTA exhibited improved capacitance retention after pyrolysis compared to the PR48-P electrodes. Finally, this work demonstrates the benefits of using topology optimization to design electrodes and material development to improve functional properties of 3D printable electrodes.

25 ENERGY STORAGE↗

A scalable DG solver for the electroneutral Nernst-Planck equations

The robust, scalable simulation of flowing electrochemical systems is increasingly important due to the synergy between intermittent renewable energy and electrochemical technologies such as energy storage and chemical manufacturing. The high Péclet regime of many such applications prevents the use of off-the-shelf discretization methods. In this work, we present a high-order Discontinuous Galerkin scheme for the electroneutral Nernst-Planck equations. The chosen charge conservation formulation allows for the specific treatment of the different physics: upwinding for advection and migration, and interior penalty for diffusion of ionic species as well the electric potential. Similarly, the formulation enables different treatments in the preconditioner: AMG for the potential blocks and ILU-based methods for the advection-dominated concentration blocks. Here we evaluate the convergence rate of the discretization scheme through numerical tests. Strong scaling results for two preconditioning approaches are shown for a large 3D flow-plate reactor example.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Topology optimization for the design of porous electrodes

Porous electrodes are an integral part of many electrochemical devices since they have high porosity to maximize electrochemical transport and high surface area to maximize activity. Traditional porous electrode materials are typically homogeneous, stochastic collections of small-scale particles and offer few opportunities to engineer higher performance. Fortunately, recent breakthroughs in advanced and additive manufacturing are yielding new methods to structure and pattern porous electrodes across length scales. These architected electrodes are emerging as a promising new technology to continue to drive improvement; however, it is still unclear which structures to employ and few tools are available to guide their design. Here we address this gap by applying topology optimization to the design of porous electrodes. We demonstrate our framework on two applications: a porous electrode driving a steady Faradaic reaction and a transiently operated electrode in a supercapacitor. We present computationally designed electrodes that minimize energy losses in a half-cell. For low-conductivity materials, the optimization algorithm creates electrode designs with a hierarchy of length scales. Further, the designed electrodes are found to outperform undesigned, homogeneous electrodes. Finally, we present three-dimensional porous electrode designs. We thus establish a topology optimization framework for designing porous electrodes.

25 ENERGY STORAGE↗

Neural Networks for Supercapacitor Modeling

This software implements a one-dimensional model of a supercapacitor and calibrates it to experimental data. A system of coupled partial differential equation models the salt concentration and the ionic and electronic potentials. The term in the equation modeling the capacitance contains a neural network, whose parameters are optimized to fit the provided experimental data. The software is written in python using the open-source library JAX to calculate the gradients needed for the optimization.

Salazar de Troya, Miquel↗

Topology Optimization for Porous Electrodes

This code is for the topology optimization of porous electrodes. It was used to produce the results in Roy, T., Salazar de Troya, M.A., Worsley, M.A. and Beck, V.A., 2021. Topology optimization for the design of porous electrodes. arXiv preprint arXiv:2111.12059.

Roy, Thomas↗

Electrospray device

An electrospray device includes an electrospray emitter adapted to receive electrospray fluid; an extractor plate spaced from the electrospray emitter and having at least one aperture; and a power supply for applying a first voltage between the extractor plate and emitter for generating at least one Taylor cone emission through the aperture to create an electrospray plume from the electrospray fluid, the extractor plate as well as accelerator and shaping plates may include a porous, conductive medium for transporting and storing excess, accumulated electrospray fluid away from the aperture.

Demmons, Nathaniel↗

Delivery of Colloid Micro-Newton Thrusters for the Space Technology 7 Mission

Two flight-qualified clusters of four Colloid Micro-Newton Thruster (CMNT) systems have been delivered to the Jet Propulsion Laboratory (JPL). The clusters will provide precise spacecraft control for the drag-free technology demonstration mission, Space Technology 7 (ST7). The ST7 mission is sponsored by the NASA New Millennium Program and will demonstrate precision formation flying technologies for future missions such as the Laser Interferometer Space Antenna (LISA) mission. The ST7 disturbance reduction system (DRS) will be on the ESA LISA Pathfinder spacecraft using the European gravitational reference sensor (GRS) as part of the ESA LISA Technology Package (LTP). Developed by Busek Co. Inc., with support from JPL in design and testing, the CMNT has been developed over the last six years into a flight-ready and flight-qualified microthruster system, the first of its kind. Recent flight-unit qualification tests have included vibration and thermal vacuum environmental testing, as well as performance verification and acceptance tests. All tests have been completed successfully prior to delivery to JPL. Delivery of the first flight unit occurred in February of 2008 with the second unit following in May of 2008. Since arrival at JPL, the units have successfully passed through mass distribution, magnetic, and EMI/EMC measurements and tests as part of the integration and test (I&T) activities including the integrated avionics unit (IAU). Flight software sequences have been tested and validated with the full flight DRS instrument successfully to the extent possible in ground testing, including full functional and 72 hour autonomous operations tests. Delivery of the cluster assemblies along with the IAU to ESA for integration into the LISA Pathfinder spacecraft is planned for the summer of 2008 with a planned launch and flight demonstration in late 2010.

Precision Formation Flying↗