Researcher profile

Rainer Helmig

Rainer Helmig contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

8 published item(s)

preprint2022arXiv

Flux-mortar mixed finite element methods with multipoint flux approximation

The flux-mortar mixed finite element method was recently developed for a general class of domain decomposition saddle point problems on non-matching grids. In this work we develop the method for Darcy flow using the multipoint flux approximation as the subdomain discretization. The subdomain problems involve solving positive definite cell-centered pressure systems. The normal flux on the subdomain interfaces is the mortar coupling variable, which plays the role of a Lagrange multiplier to impose weakly continuity of pressure. We present well-posedness and error analysis based on reformulating the method as a mixed finite element method with a quadrature rule. We develop a non-overlapping domain decomposition algorithm for the solution of the resulting algebraic system that reduces it to an interface problem for the flux-mortar, as well as an efficient interface preconditioner. A series of numerical experiments is presented illustrating the performance of the method on general grids, including applications to flow in complex porous media.

preprint2022arXiv

Modeling of two phase flow in a hydrophobic porous medium interacting with a hydrophilic structure

Fluid flow through layered materials with different wetting behavior is observed in a wide range of applications in biological, environmental and technical systems. Therefore, it is necessary to understand the occuring transport mechanisms of the fluids at the interface between the layered constituents. Of special interest is the water transport in polymer electrolyte membrane fuel cells (PEM FC). Here, it is necessary to understand the transport mechanisms of water throughout the cell constituents especially on the cathode side, where the excess water has to be removed. This is crucial to choose optimal operating conditions and improve the overall cell performance. Pore-scale modeling of gas diffusion layers (GDLs) and gas distributor has been established as a favorable technique to investigate the ongoing processes. Investigating the interface between the hydrophobic porous GDL and the hydrophilic gas distributor, a particular challenge is the combination and interaction of the different material structures and wetting properties at the interface and its influence on the flow. In this paper, a modeling approach is presented which captures the influence of a hydrophilic domain on the flow in a hydrophobic porous domain at the interface between the two domains. A pore-network model is used as the basis of the developed concept which is extended to allow the modeling of mixed-wet interactions at the interface. The functionality of the model is demonstrated using basic example configurations with one and several interface pores and it is applied to a realistic GDL representation in contact with a channel-land structured gas distributor.

preprint2022arXiv

Two-phase flow dynamics at the interface between GDL and gas distributor channel using a pore-network model

For improved operating conditions of a polymer electrolyte membrane (PEM) fuel cell, a sophisticated water management is crucial. Therefore, it is necessary to understand the transport mechanisms of water throughout the cell constituents especially on the cathode side, where the excess water has to be removed. Pore-scale modeling of diffusion layers and gas distributor has been established as a favorable technique to investigate the ongoing processes. Investigating the interface between the cathode layers, a particular challenge is the combination and interaction of the multi-phase flow in the porous material of the gas diffusion layer (GDL) with the free flow in the gas distributor channels. The formation, growth and detachment of water droplets on the hydrophobic, porous surface of the GDL have a major influence on the mass, momentum and energy exchange between the layers. A dynamic pore-network model is used to describe the flow through the porous GDL on the pore-scale. To capture the droplet occurrence and its influence on the flow, this dynamic two-phase pore-network model is extended to capture droplet formation and growth at the surface of the GDL as well as droplet detachment due to the gas flow in the gas distributor channels. In this article, the developed model is applied to single- and multi-tube systems to investigate the general drop behavior. These rather simple test-cases are compared to experimental and numerical data available in the literature. Finally, the model is applied to a GDL unit cell to analyse the interaction between two-phase flow through the GDL and drop formation at the interface between GDL and gas distributor channel.

preprint2021arXiv

Evaporation-driven density instabilities in saturated porous media

Soil salinization is a major cause of soil degradation and hampers plant growth. For soils saturated with saline water, the evaporation of water induces accumulation of salt near the top of the soil. The remaining liquid gets an increasingly larger density due to the accumulation of salt, giving a gravitationally unstable situation, where instabilities in the form of fingers can form. These fingers can hence lead to a net downwards transport of salt. We here investigate the appearance of these fingers through a linear stability analysis and through numerical simulations. The linear stability analysis gives criteria for onset of instabilities for a large range of parameters. Simulations using a set of parameters give information also about the development of the fingers after onset. With this knowledge we can predict whether and when the instabilities occur, and their effect on the salt concentration development near the top boundary.

preprint2020arXiv

Flux-mortar mixed finite element methods on non-matching grids

We investigate a mortar technique for mixed finite element approximations of Darcy flow on non-matching grids in which the normal flux is chosen as the coupling variable. It plays the role of a Lagrange multiplier to impose weakly continuity of pressure. In the mixed formulation of the problem, the normal flux is an essential boundary condition and it is incorporated with the use of suitable extension operators. Two such extension operators are considered and we analyze the resulting formulations with respect to stability and convergence. We further generalize the theoretical results, showing that the same domain decomposition technique is applicable to a class of saddle point problems satisfying mild assumptions. An example of coupled Stokes-Darcy flows is presented.

preprint2019arXiv

A new and consistent well model for one-phase flow in anisotropic porous media using a distributed source model

A new well model for one-phase flow in anisotropic porous media is introduced, where the mass exchange between well and a porous medium is modeled by spatially distributed source terms over a small neighborhood region. To this end, we first present a compact derivation of the exact analytical solution for an arbitrarily oriented, infinite well cylinder in an infinite porous medium with anisotropic permeability tensor in R3 , for constant well pressure and a given injection rate, using a conformal map. The analytical solution motivates the choice of a kernel function to distribute the sources. The presented model is independent from the discretization method and the choice of computational grids. In numerical experiments, the new well model is shown to be consistent and robust with respect to rotation of the well axis, rotation of the permeability tensor, and different anisotropy ratios. Finally, a comparison with a Peaceman-type well model suggests that the new scheme leads to an increased accuracy for injection (and production) rates for arbitrarily-oriented pressure-controlled wells.

preprint2019arXiv

Coupling Staggered-Grid and MPFA Finite Volume Methods for Free Flow/Porous-Medium Flow Problems

A discretization is proposed for models coupling free flow with anisotropic porous medium flow. Our approach employs a staggered grid finite volume method for the Navier-Stokes equations in the free flow subdomain and a MPFA finite volume method to solve Darcy flow in the porous medium. After appropriate spatial refinement in the free flow domain, the degrees of freedom are conveniently located to allow for a natural coupling of the two discretization schemes. In turn, we automatically obtain a more accurate description of the flow field surrounding the porous medium. Numerical experiments highlight the stability and applicability of the scheme in the presence of anisotropy and show good agreement with existing methods, verifying our approach.

preprint2019arXiv

Modeling tissue perfusion in terms of 1d-3d embedded mixed-dimension coupled problems with distributed sources

We present a new method for modeling tissue perfusion on the capillary scale. The microvasculature is represented by a network of one-dimensional vessel segments embedded in the extra-vascular space. Vascular and extra-vascular space exchange fluid over the vessel walls. This exchange is modeled by distributed sources using smooth kernel functions for the extra-vascular domain. It is shown that the proposed method may significantly improve the approximation of the exchange flux, in comparison with existing methods for mixed-dimension embedded problems. Furthermore, the method exhibits better convergence rates of the relevant quantities due to the increased regularity of the extra-vascular pressure solution. Numerical experiments with a vascular network from the rat cortex show that the error in the approximation of the exchange flux for coarse grid resolution may be decreased by a factor of $3$. This may open the way for computing on larger network domains, where a fine grid resolution cannot be achieved in practical simulations due to constraints in computational resources, for example in the context of uncertainty quantification.