BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ali Mani

Ali Mani appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynphysics.comp-phphysics.chem-phmath.NANumerical Analysiscond-mat.mes-hallphysics.ao-phphysics.app-ph

Catalog footprint

What is connected

15works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

15 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A multi-scale model for electrokinetic transport in networks of micro-scale and nano-scale pores

We present an efficient and robust numerical model for simulation of electrokinetic phenomena in porous networks over a wide range of applications including energy conversion, desalination, and lab-on-a-chip systems. Coupling between fluid flow and ion transport in these networks is governed by the Poisson-Nernst-Planck-Stokes equations. These equations describe a wide range of transport phenomena that can interact in complex and highly nonlinear ways in networks involving multiple pores with variable properties. Capturing these phenomena by direct simulation of the governing equations in multiple dimensions is prohibitively expensive. We present here a reduced order computational model that treats a network of many pores via solutions to 1D equations. Assuming that each pore in the network is long and thin, we derive a 1D model describing the transport in pore's longitudinal direction. We take into account the non-uniformity of potential and ion concentration profiles across the pore cross-section in the form of area-averaged coefficients in different flux terms representing fluid flow, electric current, and ion fluxes. Distinct advantages of the present framework include: a fully conservative discretization, fully bounded tabulated area-averaged coefficients without any singularity in the limit of infinitely thick electric double layers (EDLs), a flux discretization that exactly preserves equilibrium conditions, and extension to general network of pores with multiple intersections. By considering a hierarchy of canonical problems with increasing complexity, we demonstrate that the developed framework can capture a wide range of phenomena. Example demonstrations include, prediction of osmotic pressure built up in thin pores subject to concentration gradient, propagation of deionization shocks and induced recirculations for intersecting pores with varying properties.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Description of clustering of inertial particles in turbulent flows via finite-time Lyapunov exponents

An asymptotic solution is derived for the motion of inertial particles exposed to Stokes drag in an unsteady random flow. This solution provides the finite-time Lyapunov exponents as a function of Stokes number and Lagrangian strain- and rotation-rates autocovariances. The sum of these exponents, which corresponds to a concentration-weighted divergence of particle velocity field, is considered as a measure of clustering. For inertial particles dispersed in an isotropic turbulent flow our analysis predicts maximum clustering at an intermediate Stokes number and minimal clustering at small and large Stokes numbers. Direct numerical simulations are performed for quantitative validation of our analysis, showing a reasonable agreement between the two.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the coupling between buoyancy forces and electroconvective instability near ion-selective surfaces

Recent investigations have revealed that ion transport from aqueous electrolytes to ion-selective surfaces is subject to electroconvective instability that stems from coupling of hydrodynamics with electrostatic forces. Electroconvection is shown to enhance ion mixing and the net rate of transport. However, systems subject to electroconvection inherently involve fluid density variation set by salinity gradient in the bulk fluid. In this study we thoroughly examine the interplay of gravitational convection and chaotic electroconvection. Our results reveal that buoyant forces can significantly influence the transport rates, otherwise set by electroconvection, when the Rayleigh number $Ra$ of the system exceeds a value $Ra \sim 1000$. We show that buoyancy forces can significantly alter the flow patterns in these systems. When the buoyancy acts in the stabilizing direction, it limits the extent of penetration of electroconvection, but without eliminating it. When the buoyancy destabilizes the flow, it alters the electroconvective patterns by introducing upward and downward fingers of respectively light and heavy fluid.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Parallel variable-density particle-laden turbulence simulation

We have developed a fully parallel C++/MPI based simulation code for variable-density particle-laden turbulent flows. The fluid is represented through a uniform Eulerian staggered grid, while particles are modeled using a Lagrangian point-particle framework. Spatial discretization is second-order accurate, and time integration has a fourth-order accuracy. Two-way coupling of the particles with the background flow is considered in both momentum and energy equations. The code is fully modular and abstracted, and easily can be extended or modified. We have considered two different boundary conditions. We have also developed a novel parallel linear solver for the variable density Poisson equation that arises in the calculation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accurate calculation of Stokes drag for point-particle tracking in two-way coupled flows

In this work, we propose and test a method for calculating Stokes drag applicable to particle-laden fluid flows where two-way momentum coupling is important. In the point-particle formulation, particle dynamics are coupled to fluid dynamics via a source term that appears in the respective momentum equations. When the particle Reynolds number is small and the particle diameter is smaller than the fluid scales, it is common to approximate the momentum coupling source term as the Stokes drag. The Stokes drag force depends on the difference between the undisturbed fluid velocity evaluated at the particle location, and the particle velocity. However, owing to two-way coupling, the fluid velocity is modified in the neighborhood of a particle, relative to its undisturbed value. This causes the computed Stokes drag force to be underestimated in two-way coupled point-particle simulations. We develop estimates for the drag force error as function of the particle size relative to the grid size. We then develop a correction method that estimates the undisturbed fluid velocity from the computed disturbed velocity field. The correction scheme is tested for a particle settling in an otherwise quiescent fluid and is found to reduce the error in computed settling velocity by an order of magnitude compared with common interpolation schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electro-Osmotic Instability and Chaos near Ion-Selective Surfaces

We present a "fluid dynamics video" showing results from direct numerical simulations of electrokinetic transport in a binary electrolyte in two systems: (1) next to an ion-selective membrane and (2) around a metallic cylinder. The governing equations are the Poisson--Nernst--Planck and Navier--Stokes and our DNS resolves the time and length scales associated with the thin electric double layers. Our results show transition to chaotic dynamics when the applied voltage across the characteristic length scale is above a threshold of approximately 0.5 Volt. More details, in particular regarding system (1), are given in our paper [Druzgalski et al. POF 25, 110804 (2013)].

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Overlimiting Current and Shock Electrodialysis in Porous Media

Most electrochemical processes, such as electrodialysis, are limited by diffusion, but in porous media, surface conduction and electro-osmotic flow also contribute to ionic fluxes. In this paper, we report experimental evidence for surface-driven over-limiting current (faster than diffusion) and deionization shocks (propagating salt removal) in a porous medium. The apparatus consists of a silica glass frit (1 mm thick with 500 nm mean pore size) in an aqueous electrolyte (CuSO$_4$ or AgNO$_3$) passing ionic current from a reservoir to a cation-selective membrane (Nafion). The current-voltage relation of the whole system is consistent with a proposed theory based on the electro-osmotic flow mechanism over a broad range of reservoir salt concentrations (0.1 mM - 1.0 M), after accounting for (Cu) electrode polarization and pH-regulated silica charge. Above the limiting current, deionized water ($\approx 10 μ$ $M$) can be continuously extracted from the frit, which implies the existence of a stable shock propagating against the flow, bordering a depleted region that extends more than 0.5mm across the outlet. The results suggest the feasibility of "shock electrodialysis" as a new approach to water desalination and other electrochemical separations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Desalination shocks in microstructures

Salt transport in bulk electrolytes is limited by diffusion and convection, but in microstructures with charged surfaces (e.g. microfluidic devices, porous media, soils, or biological tissues) surface conduction and electro-osmotic flow also contribute to ionic fluxes. For small applied voltages, these effects lead to well known linear electrokinetic phenomena. In this paper, we predict some surprising nonlinear dynamics that can result from the competition between bulk and interfacial transport at higher voltages. When counter-ions are selectively removed by a membrane or electrode, a "desalination shock" can propagate through the microstructure, leaving in its wake an ultrapure solution, nearly devoid of co-ions and colloidal impurities. We elucidate the basic physics of desalination shocks and develop a mathematical theory of their existence, structure, and stability, allowing for slow variations in surface charge or channel geometry. Via asymptotic approximations and similarity solutions, we show that desalination shocks accelerate and sharpen in narrowing channels, while they decelerate and weaken, and sometimes disappear, in widening channels. These phenomena may find applications in separations (desalination, decontamination, biological assays) and energy storage (batteries, supercapacitors) involving electrolytes in microstructures.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Overlimiting current in a microchannel

We revisit the classical problem of diffusion-limited ion transport to a membrane (or electrode) by considering the effects of charged side walls. Using simple mathematical models and numerical simulations, we identify three basic mechanisms for over-limiting current in a microchannel: (i) {\it surface conduction} carried by excess counterions, which dominates for very thin channels, (ii) convection by {\it electro-osmotic flow} on the side walls, which dominates for thicker channels and transitions to (iii) {\it electro-osmotic instability} on the membrane end in very thick channels. These intriguing electrokinetic phenomena may find applications in biological separations, water desalination, and electrochemical energy storage.

0 citations0 reviewsNo code linkedOpen access