Researcher profile

Ali Mani

Ali Mani contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2023arXiv

Direct numerical simulation of electrokinetic transport phenomena in fluids: variational multi-scale stabilization and octree-based mesh refinement

Computational modeling of charged species transport has enabled the analysis, design, and optimization of a diverse array of electrochemical and electrokinetic devices. These systems are represented by the Poisson-Nernst-Planck (PNP) equations coupled with the Navier-Stokes (NS) equation. Direct numerical simulation (DNS) to accurately capture the spatio-temporal variation of ion concentration and current flux remains challenging due to the (a) small critical dimension of the diffuse charge layer (DCL), (b) stiff coupling due to fast charge relaxation times, large advective effects, and steep gradients close to boundaries, and (c) complex geometries exhibited by electrochemical devices. In the current study, we address these challenges by presenting a direct numerical simulation framework that incorporates (a) a variational multiscale (VMS) treatment, (b) a block-iterative strategy in conjunction with semi-implicit (for NS) and implicit (for PNP) time integrators, and (c) octree based adaptive mesh refinement. The VMS formulation provides numerical stabilization critical for capturing the electro-convective flows often observed in engineered devices. The block-iterative strategy decouples the difficulty of non-linear coupling between the NS and PNP equations and allows the use of tailored numerical schemes separately for NS and PNP equations. The carefully designed second-order, hybrid implicit methods circumvent the harsh timestep requirements of explicit time steppers, thus enabling simulations over longer time horizons. Finally, the octree-based meshing allows efficient and targeted spatial resolution of the DCL. These features are incorporated into a massively parallel computational framework, enabling the simulation of realistic engineering electrochemical devices. The numerical framework is illustrated using several challenging canonical examples.

preprint2022arXiv

Direct 3D observation and unraveling of electroconvection phenomena during concentration polarization at ion-exchange membranes

A decade ago, two-dimensional microscopic flow visualization proved the theoretically predicted existence of electroconvection roles as well as their decisive role in destabilizing the concentration polarization layer at ion-selective fluid/membrane interfaces. Electroconvection induces chaotic flow vortices injecting volume having bulk concentration into the ion-depleted diffusion layer at the interface. Experimental quantification of these important flow patterns have so far only been carried out in 2D. Numerical direct simulations suggest 3D features, yet experimental proof is lacking. 3D simulations are also limited in covering extended spacial and temporal scales. This study presents a new comprehensive experimental method for the time-resolved recording of the 3D electroconvective velocity field near a cation-exchange membrane. For the first time, the spatio-temporal velocity field can be visualized in 3D at multiples of the overlimiting current density. In contrast to today's simulations, these experiments cover length and time scales typical for actual electrodialytic membrane processes. We visualize coherent vortex structures and reveal the changes in the velocity field and its statistics during the transition from vortex rolls to vortex rings with increasing current density. The transition is characterized by changes in the rotational direction, mean square velocity, and temporal energy spectrum with only little influence on the spatial spectrum. These findings indicate a more significant impact of EC's structural change on the mean square velocities and temporal spectra than on the spatial spectra. This knowledge is a prerequisite for engineering ion-selective surfaces that will enable the operation of electrically driven processes beyond the diffusion-limited Nernst regime.

preprint2022arXiv

Eddy diffusivity in homogeneous isotropic turbulence

We use the recently developed Macroscopic Forcing Method [Mani and Park, Physical Review Fluids, 6:054607, 2021] 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 wavenumber. Behavior at all scales is reasonably captured by a non-local eddy diffusivity operator modeled as $D/\sqrt{\mathcal{I}-l^2\nabla^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 the first direct measurement of eddy diffusivity in turbulence with implications in turbulence modeling.

preprint2020arXiv

A conservative diffuse-interface method for compressible two-phase flows

In this article, we propose a novel conservative diffuse-interface method for the simulation of immiscible compressible two-phase flows. The proposed method discretely conserves the mass of each phase, momentum and total energy of the system. We use the baseline five-equation model and propose interface-regularization (diffusion--sharpening) terms in such a way that the resulting model maintains the conservative property of the underlying baseline model; and lets us use a central-difference scheme for the discretization of all the operators in the model, which leads to a non-dissipative implementation that is crucial for the simulation of turbulent flows and acoustics. Furthermore, the provable strengths of the proposed model are: (a) the model maintains the boundedness property of the volume fraction field, which is a physical realizability requirement for the simulation of two-phase flows, (b) the proposed model is such that the transport of volume fraction field inherently satisfies the total-variation-diminishing property without having to add any flux limiters that destroy the non-dissipative nature of the scheme, (c) the proposed interface-regularization terms in the model do not spuriously contribute to the kinetic energy of the system and therefore do not affect the non-linear stability of the numerical simulation, and (d) the model is consistent with the second law of thermodynamics. Finally, we present numerical simulations using the model and assess (a) the accuracy of evolution of the interface shape, (b) implementation of surface tension effects, (c) propagation of acoustics and their interaction with material interfaces, (d) the accuracy and robustness of the numerical scheme for simulation of complex high-Reynolds-number flows, and (e) performance and scalability of the method.

preprint2020arXiv

Time-Accurate and highly-Stable Explicit operators for stiff differential equations

Unconditionally stable implicit time-marching methods are powerful in solving stiff differential equations efficiently. In this work, a novel framework to handle stiff physical terms implicitly is proposed. Both physical and numerical stiffness originating from convection, diffusion and source terms (typically related to reaction) can be handled by a set of predefined Time-Accurate and highly-Stable Explicit (TASE) operators in a unified framework. The proposed TASE operators act as preconditioners on the stiff terms and can be deployed to any existing explicit time-marching methods straightforwardly. The resulting time integration methods remain the original explicit time-marching schemes, yet with nearly unconditional stability. The TASE operators can be designed to be arbitrarily high-order accurate with Richardson extrapolation such that the accuracy order of original explicit time-marching method is preserved. Theoretical analyses and stability diagrams show that the $s$-stages $s$th-order explicit Runge-Kutta (RK) methods are unconditionally stable when preconditioned by the TASE operators with order $p \leq s$ and $p \leq 2$. On the other hand, the $s$th-order RK methods preconditioned by the TASE operators with order of $p \leq s$ and $p > 2$ are nearly unconditionally stable. The only free parameter in TASE operators can be determined a priori based on stability arguments. Unlike classical implicit methods, the TASE methodology allows for solving non-linear problems with arbitrary order without requiring solving a nonlinear system of equations. A set of benchmark problems with strong stiffness is simulated to assess the performance of the TASE method. Numerical results suggest that the proposed framework preserves the high-order accuracy of the explicit time-marching methods with very-large time steps for all the considered cases.

preprint2019arXiv

A modal analysis of the behavior of inertial particles in turbulence

The clustering of small heavy inertial particles subjected to Stokes drag in turbulence is known to be minimal at small and large Stokes number and substantial at $\rm St = \mathcal O(1)$. This non-monotonic 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 $τQ$ as $τ\to 0$ and $τ^{-1/2} |Q|^{1/4}$ as $τ\to \infty$, where $τ$ is the particle relaxation time and $Q(τ)$ is the difference between the norm of the rotation- and strain-rate tensors computed along the particle trajectory. Additionally, the trajectory crossing is predicted only in hyperbolic flows where $Q<0$ for sufficiently inertial particles with a $τ$ 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. These findings, which are stemmed from our analysis, are corroborated by the direct numerical simulations.