Search NASA⌕ Search

Engineering topics

Mani, Ali

Publications and source records attributed to Mani, Ali.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Carbonate Management to Enable Energy- and Carbon-Efficient CO 2 Electrolysis (Final Technical Report)

The rapid growth and plummeting cost of solar energy have spurred growing interest in using CO 2 electrolysis to produce chemicals and fuels as an alternative to conventional petrochemical processes. High-temperature (>800 °C) solid oxide electrolyzers that convert CO 2 into CO and O 2 have recently become commercially available. Low-temperature electrolysis cells offer the prospect of more convenient and flexible operation, which is critical for utilizing intermittent solar energy, and provide access to more valuable C 2+ products such as ethylene, ethanol, and propanol. Over the past 10 years, research in this area has yielded substantial progress in both fundamental understanding of the requisite electrocatalytic reactions and design of prototype devices. Leveraging insights from fuel cells and membrane water electrolyzers, researchers have developed electrolysis cells with gas diffusion electrodes (GDE) that have demonstrated high CO 2 electrolysis current densities (>100 mA cm –2 ) as well as promising selectivity and stability. Despite these advances, the energy efficiency (electrical energy-to-product) and carbon efficiency (CO 2 -to-product) of low-temperature CO 2 electrolysis remain far too low for large-scale deployment. A preponderance of evidence indicates that the principal source of efficiency losses is the rapid and thermodynamically favorable reaction of CO 2 with hydroxide (OH – ) to form carbonate (CO 3 2– ). Carbonate formation imposes steady-state electrolysis conditions that result in large voltage and CO 2 losses for all known (photo)electrochemical CO 2 cells. While much current research remains focused on CO 2 reduction catalyst design, this largely overlooked CO 3 2– problem presents a fundamental scientific barrier to creating a viable electrochemical option for converting solar energy into chemicals and fuels. The project pursues an integrated, multi-PI research effort that establishes a fundamental science of CO 3 2– management. PI Kanan and Co-PI Mani’s contribution to the project is to evaluate strategies to mitigate the CO 3 2– problem by changing the properties of the electrolyte and the environment in which CO 2 reduction catalysis takes place. Experimental studies showed that electrolytes composed of a high concentration of both CO 3 2– and HCO 3 – , which serve as moderately alkaline buffers, improved the cell voltage by compared to all-HCO 3 – electrolytes, but these buffered systems still show substantial CO 2 uptake that reduces pH over time. Computational studies developed a homogenized model of a CO 2 reduction catalyst layer that permits a low-cost exploration of the high-dimensional parameter space associated with catalyst layers on gas diffusion electrodes. The model was validated by accurately reproducing experimental data for the related but simpler reaction of CO reduction and then used to probe the effects of catalyst layer architecture on CO 2 reduction. Minimizing the size of catalyst and hydrophobic domains in the catalyst layer is predicted to mitigate CO 3 2– formation and thereby enable prolonged operation at elevated pH. In support of the CO 2 electrolysis studies, a new method for rapidly prototyping electrochemical cells was developed and validated. The method uses a combination of 3D printing and electroless plating to generate conductive cell components for evaluating new cell designs. The carbonate problem encompasses mass transport processes and acid-base reactions that are relevant to many other electrochemical systems. Investigation of strategies to address the carbonate problem led to an additional line of inquiry into the physicochemical phenomena that determine the efficiency of electrochemical acid-base production, which has numerous applications in the broader field of carbon management. New strategies for using the supporting electrolyte to inhibit H + /OH – recombination in electrochemical acid-base production were evaluated, leading to the development of a novel acid-base producing system that eliminates the need for ion exchange membranes and exhibits promising efficiency and current densities for scalable applications.

25 ENERGY STORAGE↗
DOE OSTI2023

On the convergence of statistics in simulations of stationary incompressible turbulent flows

When reporting statistics from simulations of statistically stationary chaotic phenomenon, it is important to verify that the simulations are time-converged. This condition is connected with the statistical error or number of digits with which statistics can be reliably reported. In this work we consider numerical experiments of low Reynolds number incompressible homogeneous and isotropic turbulence as a model problem to investigate statistical convergence over finite simulation times. Specifically, we investigate the time integration requirements that allow meaningful reporting of the statistical error associated with finiteness of the temporal domain. We address two key questions: (1) How long should a simulation be performed in terms of large eddy time, and (2) How should the simulation time be divided among temporal windows over which a quantity of interest is estimated so that its statistical error could be reliably reported? We find that reliable reporting of statistical errors requires simulations on the order of 10 4 large eddy times, which is orders of magnitude longer than typically performed. Additionally, data post-processing should employ windows of at least ten times the large eddy time scale, with the most robust computation of statistical error of the mean requiring window sizes of an additional factor of ten. For practical simulations, we demonstrate that it is possible to estimate the statistical error within a factor of two under a less stringent condition in which a minimum of four windows with size at least ten large eddy times are used. In conclusion, our observations for homogeneous isotropic turbulence are also shown to hold in turbulent channel flow.

42 ENGINEERING↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗
DOE OSTI2022

Eddy diffusivity operator in homogeneous isotropic turbulence

Here, we use the recently developed macroscopic forcing method to compute the scale-dependent eddy diffusivity characterizing ensemble-averaged scalar and momentum transport in incompressible homogeneous isotropic turbulence. For scales larger than the energy containing eddies, eddy diffusivity is found to be constant and consistent with the Boussinesq approximation. However, for small scales eddy diffusivity is found to vanish inversely proportional to the wave number. Behavior at all scales is reasonably captured by a nonlocal eddy diffusivity operator modeled as D/$\sqrt{I - l^2∇^2}$, where D is the eddy diffusivity in the Boussinesq limit, and l is a constant on the order of the large-eddy length. These results present a direct measurement of eddy diffusivity in turbulence with implications in turbulence modeling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A novel energy balance approach for a verifiable and accurate solution of radiation extinction in purely absorbing particle clouds

Here we consider the problem of radiation transport through purely absorbing particle clouds. The gold-standard solution of particle-resolved Monte Carlo ray-tracing method to this problem is computationally expensive and therefore solving the radiation transport equation (specifically, the Beer-Bouguer law (BB-law)) on a Eulerian mesh is often preferred. While the absorption coefficient in the real problem is infinite, the BB-law approximates it to be a finite number in the form of number density through a set of assumptions. For particle clouds that do not obey these assumptions, the BB-law predicts an incorrect exponential decay. Also, when the number density is computed using the nearest-neighbor approach, the BB-law solution diverges when the Eulerian mesh size becomes closer to or smaller than the particle size. This numerical divergence is due to the homogenization error. Although the filtering strategy for number density minimizes the homogenization error, it still converges to an incorrect exponential decay for particle clouds that break the BB-law constraints and the cost of the filtering is equivalent to that of the gold-standard solution. In this study, we develop a novel, highly accurate, verifiable and cost-effective solution to the radiation transport equation on a Eulerian domain using an energy balance approach where we derive an expression for the absorption coefficient as a function of particle and Eulerian mesh sizes. We apply our new method to Poisson and turbulent particle clouds that violate all the BB-law constraints and show that the solution is converging upon the mesh refinement, and eventually, we recover the same gold-standard solution for a much cheaper computational cost.

42 ENGINEERING↗
DOE OSTI2021

Settling of two-way momentum and energy coupled particles subject to Boussinesq and non-Boussinesq heating

This work establishes a procedure to accurately compute heat transfer between an Eulerian fluid and Lagrangian point-particles. Recent work has focused on accurately computing momentum transfer between fluid and particles. The coupling term for momentum involves the undisturbed fluid velocity at the particle location which is not directly accessible in the simulation. Analogously, in the context of thermal coupling, the undisturbed fluid temperature at the particle location is not directly accessible in simulations and must be estimated. In this paper, we develop a scheme to accurately estimate the undisturbed fluid temperature of a point-particle exchanging thermal energy with a surrounding fluid. Furthermore, the temperature disturbance is correlated with the enhanced temperature curvature in the vicinity of the particle and is formally valid in the low heating, low convection limit. We conduct extensive verification of the correction procedure for a settling particle subject to radiation. This setup allows the simultaneous testing of thermal and momentum corrections. By considering equations of drag and Nusselt number extended to finite Péclet and Boussinesq numbers, we establish a large range over which the correction procedure can be applied.

42 ENGINEERING↗
DOE OSTI2020

Modal analysis of the behavior of inertial particles in turbulence subjected to Stokes drag

The clustering of small heavy inertial particles subjected to Stokes drag in turbulence, which is the concern of present study, is known to be minimal at small and large Stokes number and substantial at St = $\textit{O}$(1). This nonmonotonic trend, which has been shown computationally and experimentally, is yet to be explained analytically. In this study, we obtain an analytical expression for the Lyapunov exponents that quantitatively predicts this trend. The sum of the exponents, which is the normalized rate of change of the signed volume of a small cloud of particles, is correctly predicted to be negative and positive at small and large Stokes numbers, respectively, asymptoting to $\textit{τQ}$ as $\textit{τ}$ → 0 and $τ^{–1/2}|Q|^{1/4}$ as $\textit{τ} → ∞$, where $\textit{τ}$ is the particle relaxation time and $\textit{Q(τ)}$ is the difference between the norm of the rotation- and strain-rate tensors computed along the particle trajectory, which remains negative at all $\textit{τ}$ in turbulent flows. Additionally, the trajectory crossing is predicted only in hyperbolic flows where Q < 0 for sufficiently inertial particles with a $\textit{τ}$ that scales with $|Q|^{–1/2}$. Following the onset of crossovers, a transition from clustering to dispersion is predicted correctly. We show these behaviors are not unique to three-dimensional isotropic turbulence and can be reproduced closely by a one-dimensional mono-harmonic flow, which appears as a fundamental canonical problem in the study of particle clustering. Analysis of this one-dimensional canonical flow shows that the rate of clustering, quantified as the product of the Lyapunov exponent and particle relaxation time, is bounded by –1/2, behaving with extreme nonlinearity in the hyperbolic flows and always remaining positive in the elliptic flows. Finally, this model problem also confirms that dispersion in a hyperbolic flow will occur only if particles' trajectory cross. These findings, which stem from our analysis, are corroborated by the direct numerical simulations.

42 ENGINEERING↗
DOE OSTI2020

A filtering strategy for the numerical convergence of radiation transport through purely absorbing particle clouds

Radiation transport through particle clouds plays a major role in many engineering applications. Here in this paper, we study this problem by means of solving the radiative transport equation (RTE) on a Eulerian continuous domain. The particle clouds are generated through direct numerical simulation of Navier-Stokes equations coupled with Lagrangian particle transport at different Stokes numbers. We confirm the earlier observation noted in the literature that the solution to the RTE on Eulerian mesh diverges when the Eulerian mesh size is of the order of the particle diameter. This observation is often called the homogenization error stemming from relegation of number density onto the Eulerian domain. In order to circumvent this divergence problem, we propose a filtering strategy that spreads the information of the particles to the neighboring cells in a way that the representation of the particles remains the same even when the mesh size is smaller than that of the particle size. We show through our simulations that our filtering strategy solves the issue of homogenization error for the cases where the Beer-Bouger law is valid.

42 ENGINEERING↗
DOE OSTI2018

Particle-to-fluid heat transfer in particle-laden turbulence

Preferential concentration of inertial particles by turbulence is a well-recognized phenomenon. This work investigates how this phenomenon impacts the mean heat transfer between the fluid phase and the particle phase. Using direct numerical simulations of homogeneous and isotropic turbulent flows coupled with Lagrangian point particle tracking, we explore this phenomenon over a wide range of input parameters. Among the nine independent dimensionless numbers defining this problem, we show that the particle Stokes number, defined based on a large-eddy time, and an identified number called the heat-mixing parameter have the most significant effect on particle-to-gas heat transfer, while variation in other nondimensional numbers can be ignored. An investigation of regimes with significant particle mass loading suggests that the mean heat transfer from particles to gas is hardly affected by momentum two-way coupling. Using our numerical results, we propose an algebraic reduced-order model for heat transfer in particle-laden turbulence.

42 ENGINEERING↗