Multiphase flows with surface tension

Aphros on GitHubBubbles and drops are critical components of important industrial applications such as boiling and condensation, bubble column reactors, electrochemical cells and physical systems involving air entrainment in plunging jets and liquid jet atomization.

Effects of surface tension are important for these applications and present a challenge for numerical simulations.
We develop models, numerical methods and software for simulation of such flows.
People: Petr Karnakov, Sergey Litvinov, Petros Koumoutsakos

Code:  Aphros: Finite volume solver for incompressible multiphase flows with surface tension

Membrane-less electrochemical reactors

In traditional electrochemical reactors, the products forming on the electrodes are separated by a membrane. Flow based electrochemical instead rely on hydrodynamic mechanisms to keep the products separated. One of these mechanisms is the Segre-Silberberg effect of lateral migration of a particle in a channel flow: the particle finds a stable equilibrium position near the wall. With proper flow conditions, the gaseous products in the form of bubbles can be kept separated.

MFER membraneMFER flow

Removing the membrane can reduce the ionic resistance of the reactor leading to an increase in power conversion efficiency. This, however, requires strong understanding of the underlying processes. Combining numerical modelling with experiments done by our collaborators, we aim to fully investigate the operation of flow-based electrochemical reactors and provide tools for their optimal design.

Publications

2020

  • P. Karnakov, S. Litvinov, and P. Koumoutsakos, “A hybrid particle volume-of-fluid method for curvature estimation in multiphase flows," Int. J. Multiphas. Flow, vol. 125, p. 103209, 2020.

BibTeX

@article{karnakov2020a,
author = {Petr Karnakov and Sergey Litvinov and Petros Koumoutsakos},
doi = {10.1016/j.ijmultiphaseflow.2020.103209},
journal = {{Int. J. Multiphas. Flow}},
month = {apr},
pages = {103209},
publisher = {Elsevier {BV}},
title = {A hybrid particle volume-of-fluid method for curvature estimation in multiphase flows},
url = {https://cse-lab.seas.harvard.edu/files/cse-lab/files/karnakov2020a.pdf},
volume = {125},
year = {2020}
}

Abstract

We present a particle method for estimating the curvature of interfaces in volume-of-fluid simulations of multiphase flows. The method is well suited for under-resolved interfaces, and it is shown to be more accurate than the parabolic fitting that is employed in such cases. The curvature is computed from the equilibrium positions of particles constrained to circular arcs and attracted to the interface. The pro- posed particle method is combined with the method of height functions at higher resolutions, and it is shown to outperform the current combinations of height functions and parabolic fitting. The algorithm is conceptually simple and straightforward to implement on new and existing software frameworks for multiphase flow simulations thus enhancing their capabilities in challenging flow problems. We evaluate the proposed hybrid method on a number of two- and three-dimensional benchmark flow problems and illustrate its capabilities on simulations of flows involving bubble coalescence and turbulent multiphase flows.
  • Z. Y. Wan, P. Karnakov, P. Koumoutsakos, and T. P. Sapsis, “Bubbles in turbulent flows: data-driven, kinematic models with history terms," Int. J. Multiphas. Flow, vol. 129, p. 103286, 2020.

BibTeX

@article{wan2020a,
author = {Zhong Yi Wan and Petr Karnakov and Petros Koumoutsakos and Themistoklis P. Sapsis},
doi = {10.1016/j.ijmultiphaseflow.2020.103286},
journal = {{Int. J. Multiphas. Flow}},
month = {aug},
pages = {103286},
publisher = {Elsevier {BV}},
title = {Bubbles in turbulent flows: Data-driven, kinematic models with history terms},
url = {https://cse-lab.seas.harvard.edu/files/cse-lab/files/wan2020a.pdf},
volume = {129},
year = {2020}
}

Abstract

We present data driven kinematic models for the motion of bubbles in high-Re turbulent fluid flows based on recurrent neural networks with long-short term memory enhancements. The models extend empirical relations, such as Maxey-Riley (MR) and its variants, whose applicability is limited when either the bubble size is large or the flow is very complex. The recurrent neural networks are trained on the trajectories of bubbles obtained by Direct Numerical Simulations (DNS) of the Navier Stokes equations for a two-component incompressible flow model. Long short term memory components exploit the time history of the flow field that the bubbles have encountered along their trajectories and the networks are further augmented by imposing rotational invariance to their structure. We first train and validate the formulated model using DNS data for a turbulent Taylor-Green vortex. Then we examine the model pre- dictive capabilities and its generalization to Reynolds numbers that are different from those of the train- ing data on benchmark problems, including a steady (Hill’s spherical vortex) and an unsteady (Gaussian vortex ring) flow field. We find that the predictions of the developed model are significantly improved compared with those obtained by the MR equation. Our results indicate that data-driven models with history terms are well suited in capturing the trajectories of bubbles in turbulent flows.

2019

  • S. M. H. Hashemi, P. Karnakov, P. Hadikhani, E. Chinello, S. Litvinov, C. Moser, P. Koumoutsakos, and D. Psaltis, “A versatile and membrane-less electrochemical reactor for the electrolysis of water and brine," Energ. Environ. Sci., 2019.

BibTeX

@article{hashemi2019a,
author = {Seyyed Mohammad Hosseini Hashemi and Petr Karnakov and Pooria Hadikhani and Enrico Chinello and Sergey Litvinov and Christophe Moser and Petros Koumoutsakos and Demetri Psaltis},
doi = {10.1039/c9ee00219g},
journal = {{Energ. Environ. Sci.}},
publisher = {Royal Society of Chemistry ({RSC})},
title = {A versatile and membrane-less electrochemical reactor for the electrolysis of water and brine},
url = {https://cse-lab.seas.harvard.edu/files/cse-lab/files/hashemi2019a.pdf},
year = {2019}
}

Abstract

Renewables challenge the management of energy supply and demand due to their intermittency. A promising solution is the direct conversion of the excess electrical energy into valuable chemicals in electrochemical reactors that are inexpensive, scalable, and compatible with irregular availability of electrical power. Membrane-less electrolyzers, deployed on a microfluidic platform, were recently shown to hold great promise for efficient electrolysis and cost-effective operation. The elimination of the membrane increases the reactor lifetime, reduces fabrication costs, and enables the deployment of liquid electrolytes with ionic conductivities that surpass those allowed by solid membranes. Here, we demonstrate a membrane-less architecture that enables unprecedented throughput by 3D printing a device that combines components such as the flow plates and the fluidic ports in a monolithic part while at the same time providing tight tolerances and smooth surfaces for precise flow conditioning. We show that inertial fluidic forces are effective even in milifluidic regimes and, therefore, are utilized to control the two-phase flows inside the device and prevent cross-contamination of the products. Simulations provide insight on governing fluid dynamics of coalescing bubbles and their rapid jumps away from the electrodes and help identify three key mechanisms for their fast and intriguing return towards the electrodes. Experiments and simulations are used to demonstrate the efficiency of the inertial separation mechanism in milichannels and at higher flow rates than in microchannels. We analyze the performance of the present device for two reactions: water splitting and the Chlor-Alkali process and find product purities of more than 99% and Faradaic efficiencies of more than 90%. The present membrane-less reactor - containing more efficient catalysts - provides close to 40 times higher throughput than its microfluidic counterpart and paves the way for realization of cost-effective and scalable electrochemical stacks that meet the performance and price targets of the renewable energy sector.
  • P. Karnakov, F. Wermelinger, M. Chatzimanolakis, S. Litvinov, and P. Koumoutsakos, “A high performance computing framework for multiphase, turbulent flows on structured grids," in Proceedings of the platform for advanced scientific computing conference – PASC ’19, 2019.

BibTeX

@inproceedings{karnakov2019a,
author = {Petr Karnakov and Fabian Wermelinger and Michail Chatzimanolakis and Sergey Litvinov and Petros Koumoutsakos},
booktitle = {Proceedings of the Platform for Advanced Scientific Computing Conference - {PASC} {\textquotesingle}19},
doi = {10.1145/3324989.3325727},
publisher = {{ACM} Press},
title = {A High Performance Computing Framework for Multiphase, Turbulent Flows on Structured Grids},
url = {https://cse-lab.seas.harvard.edu/files/cse-lab/files/karnakov2019a.pdf},
year = {2019}
}

Abstract

We present a high performance computing framework for mul- tiphase, turbulent flows on structured grids. The computational methods are validated on a number of benchmark problems such as the Taylor-Green vortex that are extended by the inclusion of bubbles in the flow field. We examine the effect of bubbles on the turbulent kinetic energy dissipation rate and provide extensive data for bubble trajectories and velocities that may assist the develop- ment of engineering models. The implementation of the present solver on massively parallel, GPU enhanced architectures allows for large scale and high throughput simulations of multiphase flows.

Videos

Select thumbnails to play videos