Bernaise: A flexible framework for simulating two-phase electrohydrodynamic flows in complex domains
arXiv:1805.11642 · doi:10.3389/fphy.2019.00021
Abstract
Bernaise (Binary ElectRohydrodyNAmIc SolvEr) is a flexible high-level finite element solver of two-phase electrohydrodynamic flow in complex geometries. Two-phase flow with electrolytes is relevant across a broad range of systems and scales, from 'lab-on-a-chip' devices for medical diagnostics to enhanced oil recovery at the reservoir scale. For the strongly coupled multi-physics problem, we employ a recently developed thermodynamically consistent model which combines a generalized Nernst-Planck equation for ion transport, the Poisson equation for electrostatics, the Cahn-Hilliard equation for the phase field (describing the interface separating the phases), and the Navier-Stokes equations for fluid flow. As an efficient alternative to solving the coupled system of partial differential equations in a monolithic manner, we present a linear, decoupled numerical scheme which sequentially solves the three sets of equations. The scheme is validated by comparison to limiting cases where analytical solutions are available, benchmark cases, and by the method of manufactured solution. The solver operates on unstructured meshes and is therefore well suited to handle arbitrarily shaped domains and problem set-ups where, e.g., very different resolutions are required in different parts of the domain. Bernaise is implemented in Python via the FEniCS framework, which effectively utilizes MPI and domain decomposition, and should therefore be suitable for large-scale/high-performance computing. Further, new solvers and problem set-ups can be specified and added with ease to the Bernaise framework by experienced Python users.
arXiv comment62 pages, 13 figures, 10 tables, 4 listings
References in corpus (12)
- DOLFIN: Automated Finite Element Computing
- Molecular scale contact line hydrodynamics of immiscible flows
- A variational approach to the moving contact line hydrodynamics
- Non-coalescence of oppositely charged drops
- Universality in dynamic wetting dominated by contact line friction
- Oasis: a high-level/high-performance open source Navier-Stokes solver
- Concentration polarization, surface currents, and bulk advection in a microchannel
- Wetting Boundary Conditions in Phase-Field-Based Simulation of Binary Fluids: Some Comparative Studies and New Development
- Adaptive and Iterative Methods for Simulations of Nanopores with the PNP-Stokes Equations
- A multiscale method for simulating fluid interfaces covered with large molecules such as asphaltenes
- Electrohydrodynamic channeling effects in narrow fractures and pores
- Controlling wetting with electrolytic solutions: phase-field simulations of a droplet-conductor system