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

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Kim, Sungu, Saurabh, Kumar, Khanwale, Makrand A., Mani, Ali, Anand, Robbyn K., Ganapathysubramanian, Baskar
Natura: Preprint
Pubblicazione: 2023
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866913181299900416
author Kim, Sungu
Saurabh, Kumar
Khanwale, Makrand A.
Mani, Ali
Anand, Robbyn K.
Ganapathysubramanian, Baskar
author_facet Kim, Sungu
Saurabh, Kumar
Khanwale, Makrand A.
Mani, Ali
Anand, Robbyn K.
Ganapathysubramanian, Baskar
contents 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.
format Preprint
id arxiv_https___arxiv_org_abs_2301_05985
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle Direct numerical simulation of electrokinetic transport phenomena in fluids: variational multi-scale stabilization and octree-based mesh refinement
Kim, Sungu
Saurabh, Kumar
Khanwale, Makrand A.
Mani, Ali
Anand, Robbyn K.
Ganapathysubramanian, Baskar
Numerical Analysis
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.
title Direct numerical simulation of electrokinetic transport phenomena in fluids: variational multi-scale stabilization and octree-based mesh refinement
topic Numerical Analysis
url https://arxiv.org/abs/2301.05985