MAC schemes on triangular Delaunay meshes

We study two classical generalized MAC schemes on unstructured triangular Delaunay meshes for the incompressible Stokes and Navier-Stokes equations and prove their convergence for the first time. These generalizations use the duality between Voronoi and triangles of Delaunay meshes, in order to construct two staggered discretization schemes. Both schemes are especially interesting, since compatible finite volume discretizations for coupled convection-diffusion equations can be constructed which preserve discrete maximum principles. In the first scheme, called tangential velocity scheme, the pressures are defined at the vertices of the mesh, and the discrete velocities are tangential to the edges of the triangles. In the second scheme, called normal velocity scheme, the pressures are defined in the triangles, and the discrete velocities are normal to the edges of the triangles. For both schemes, we prove the convergence in $L^2$ for the velocities and the discrete rotations of the velocities for the Stokes and the Navier-Stokes problem. Further, for the normal velocity scheme, we also prove the strong convergence of the pressure in $L^2$. Linear and nonlinear numerical examples illustrate the theoretical predictions

A numerical method for mass conservative coupling between fluid flow and solute transport

We present a new coupled discretization approach for species transport in an incompressible fluid. The Navier-Stokes equations for the flow are discretized by the divergence-free Scott-Vogelius element on barycentrically refined meshes guaranteeing LBB stability. The convection-diffusion equation for species transport is discretized by the Voronoi finite volume method. In accordance to the continuous setting, due to the exact integration of the normal component of the flow through the Voronoi surfaces, the species concentration fulfills discrete global and local maximum principles. Besides of the the numerical scheme itself, we present important aspects of its implementation. Further, for the case of homogeneous Dirichlet boundary conditions, we give a convergence proof for the coupled scheme. We report results of the application of the scheme to the interpretation of limiting current measurements in an electrochemical flow cell with cylindrical shape

A model of an electrochemical flow cell with porous layer

In this paper we discuss three different mathematical models for fluid-porous interfaces in a simple channel geometry that appears e.g. in thin-layer channel flow cells. Here the difficulties arise from the possibly different orders of the corresponding differential operators in the different domains. A finite volume discretization of this model allows to calculate the limiting current of the H_2 oxidation in a porous electrode with platinum catalyst particles

A numerical method for mass conservative coupling between fluid flow and solute transport

We present a new coupled discretization approach for species transport in an incompressible fluid. The Navier-Stokes equations for the flow are discretized by the divergence-free Scott-Vogelius element on barycentrically refined meshes guaranteeing LBB stability. The convection-diffusion equation for species transport is discretized by the Voronoi finite volume method. In accordance to the continuous setting, due to the exact integration of the normal component of the flow through the Voronoi surfaces, the species concentration fulfills discrete global and local maximum principles. Besides of the the numerical scheme itself, we present important aspects of its implementation. Further, for the case of homogeneous Dirichlet boundary conditions, we give a convergence proof for the coupled scheme. We report results of the application of the scheme to the interpretation of limiting current measurements in an electrochemical flow cell with cylindrical shape

Complex THz and DC inverse spin Hall effect in YIG/Cu1−x_{1-x}Irx_{x} bilayers across a wide concentration range

We measure the inverse spin Hall effect of Cu$_{1-x}$Ir$_{x}$ thin films on yttrium iron garnet over a wide range of Ir concentrations ($0.05 \leqslant x \leqslant 0.7$). Spin currents are triggered through the spin Seebeck effect, either by a DC temperature gradient or by ultrafast optical heating of the metal layer. The spin Hall current is detected by, respectively, electrical contacts or measurement of the emitted THz radiation. With both approaches, we reveal the same Ir concentration dependence that follows a novel complex, non-monotonous behavior as compared to previous studies. For small Ir concentrations a signal minimum is observed, while a pronounced maximum appears near the equiatomic composition. We identify this behavior as originating from the interplay of different spin Hall mechanisms as well as a concentration-dependent variation of the integrated spin current density in Cu$_{1-x}$Ir$_{x}$. The coinciding results obtained for DC and ultrafast stimuli show that the studied material allows for efficient spin-to-charge conversion even on ultrafast timescales, thus enabling a transfer of established spintronic measurement schemes into the terahertz regime.Comment: 12 pages, 4 figure

Mathematical modeling of channel-porous layer interfaces in PEM fuel cells

In proton exchange membrane (PEM) fuel cells, the transport of the fuel to the active zones, and the removal of the reaction products are realized using a combination of channels and porous diffusion layers. In order to improve existing mathematical and numerical models of PEM fuel cells, a deeper understanding of the coupling of the flow processes in the channels and diffusion layers is necessary. After discussing different mathematical models for PEM fuel cells, the work will focus on the description of the coupling of the free flow in the channel region with the filtration velocity in the porous diffusion layer as well as interface conditions between them. The difficulty in finding effective coupling conditions at the interface between the channel flow and the membrane lies in the fact that often the orders of the corresponding differential operators are different, e.g., when using stationary (Navier-)Stokes and Darcy's equation. Alternatively, using the Brinkman model for the porous media this difficulty does not occur. We will review different interface conditions, including the well-known Beavers-Joseph-Saffman boundary condition and its recent improvement by Le Bars and Worster

Models and numerical methods for electrolyte flows

The most common mathematical models for electrolyte flows are based on the dilute solution assumption, leading to a coupled system of the Nernst--Planck--Poisson drift-diffusion equations for ion transport and the Stokes resp. Navier--Stokes equations for fluid flow. This contribution discusses historical and recent model developments beyond the dilute solution assumption and focuses on the effects of finite ion sizes and solvation. A novel numerical solution approach is presented and verified here which aims at preserving on the discrete level consistency with basic thermodynamic principles and structural properties like independence of flow velocities from gradient contributions to external forces