Source author record

Peter Knabner

Peter Knabner appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

15 published item(s)

preprint2022arXiv

Local existence of strong solutions to micro-macro models for reactive transport in evolving porous media

Two-scale models pose a promising approach in simulating reactive flow and transport in evolving porous media. Classically, homogenized flow and transport equations are solved on the macroscopic scale, while effective parameters are obtained from auxiliary cell problems on possibly evolving reference geometries (micro-scale). Despite their perspective success in rendering lab/field-scale simulations computationally feasible, analytic results regarding the arising two-scale bilaterally coupled system often restrict to simplified models. In this paper, we first derive smooth-dependence results concerning the partial coupling from the underlying geometry to macroscopic quantities. Therefore, alterations of the representative fluid domain are described by smooth paths of diffeomorphisms. Exploiting the gained regularity of the effective space- and time-dependent macroscopic coefficients, we present local-in-time existence results for strong solutions to the partially coupled micro-macro system using fixed-point arguments. What is more, we extend our results to the bilaterally coupled diffusive transport model including a level-set description of the evolving geometry.

preprint2020arXiv

Benchmark for numerical solutions of flow in heterogeneous groundwater formations

This article presents numerical investigations on accuracy and convergence properties of several numerical approaches for simulating steady state flows in heterogeneous aquifers. Finite difference, finite element, discontinuous Galerkin, spectral, and random walk methods are tested on one- and two-dimensional benchmark flow problems. Realizations of log-normal hydraulic conductivity fields are generated by Kraichnan algorithms in closed form as finite sums of random periodic modes, which allow direct code verification by comparisons with manufactured reference solutions. The quality of the methods is assessed for increasing number of random modes and for increasing variance of the log-hydraulic conductivity fields with Gaussian and exponential correlation. Experimental orders of convergence are calculated from successive refinements of the grid. The numerical methods are further validated by comparisons between statistical inferences obtained from Monte Carlo ensembles of numerical solutions and theoretical first-order perturbation results. It is found that while for Gaussian correlation of the log-conductivity field all the methods perform well, in the exponential case their accuracy deteriorates and, for large variance and number of modes, the benchmark problems are practically not tractable with reasonably large computing resources, for all the methods considered in this study.

preprint2020arXiv

FESTUNG: A MATLAB /GNU Octave toolbox for the discontinuous Galerkin method. Part IV: Generic problem framework and model-coupling interface

This is the fourth installment in our series on implementing the discontinuous Galerkin (DG) method as an open source MATLAB /GNU Octave toolbox. Similarly to its predecessors, this part presents new features for application developers employing DG methods and follows our strategy of relying on fully vectorized constructs and supplying a comprehensive documentation. The specific focus of the current work is the newly added generic problem implementation framework and the highly customizable model-coupling interface for multi-domain and multi-physics simulation tools based on this framework. The functionality of the coupling interface in the FESTUNG toolbox is illustrated using a two-way coupled free-surface / groundwater flow system as an example application.

preprint2020arXiv

Homogenization of two-phase flow in porous media from Pore to Darcy Scale: A phase-field approach

We extend the two-scale expansion approach of periodic homogenization to include time scales and thus can tackle the full instationary Navier-Stokes-Cahn-Hilliard model at the pore scale as microscale. Time scale separation allows us to keep microscale dynamics, responsible e.g. for hysteresis, and arrive at a numerically tractable micro-macro model including coupled generalized Darcy's laws.

preprint2016arXiv

A thermodynamically consistent model for multicomponent electrolyte solutions

This paper presents a thermodynamically consistent model for multicomponent electrolyte solutions. The first part of this paper derives the general governing equations for nonequilibrium systems within the theory of nonequilibrium thermodynamics. Here, we consider electrolyte solutions as general mixtures of charged constituents. Furthermore, in this part of the paper we combine the general theory of nonequilibrium thermodynamics with the well-known splittings of the entropy and the energy into a pure substance part and a part due to mixing. Thereby, we successfully establish evolution equations for both parts. Furthermore, we derive for both parts explicit expressions of the respective entropy production rates. Hence, we provide an approach that allows to study the entropy of mixing independently of the pure substance entropy and vice versa. This is of great value, in particular for a better understanding of the complex phenomena due to mixing in multicomponent systems. In the second part of this paper, we close the system of general balance equations by applying constitutive laws. This is the crucial step in the modeling procedure. ...

preprint2016arXiv

Building a Bridge from Moments to PDF's: A New Approach to Finding PDF Mixing Models

Probability density function (PDF) methods are a promising alternative to predicting the transport of solutes in groundwater under uncertainty. They make it possible to derive the evolution equations of the mean concentration and the concentration variance, used in moment methods. A mixing model, also known as a dissipation model, is essential for both methods. Finding a satisfactory mixing model is still an open question and due to the rather elaborate PDF methods, a difficult undertaking. Both the PDF equation and the concentration variance equation depend on the same mixing model. This connection is used to find and test an improved mixing model for the much easier to handle concentration variance. Subsequently, this mixing model is transferred to the PDF equation and tested. The newly proposed mixing model yields significantly improved results for both variance modelling and PDF modelling.

preprint2016arXiv

Efficient Realization of the Mixed Finite Element Discretization for nonlinear Problems

We consider implementational aspects of the mixed finite element method for a special class of nonlinear problems. We establish the equivalence of the hybridized formulation of the mixed finite element method to a nonconforming finite element method with augmented Crouzeix-Raviart ansatz space. We discuss the reduction of unknowns by static condensation and propose Newton's method for the solution of local and global systems. Finally, we show, how such a nonlinear problem arises from the mixed formulation of Darcy-Forchheimer flow in porous media.

preprint2016arXiv

Global existence of weak solutions of a model for electrolyte solutions - Part 1: Two-component case

This paper analytically investigates the Darcy-Poisson-Nernst-Planck system. This system is a mathematical model for electrolyte solutions. In this paper, we consider electrolyte solutions, which consist of a neutral fluid and two suspended oppositely charged chemical species with arbitrary valencies z_1 > 0 > z_2. We prove global existence and uniqueness of weak solutions in two space dimensions and three space dimensions. So far, most of the existence results have been proven for symmetric electrolyte solutions. These solutions consist of a neutral fluid and two suspended charged chemical species with symmetric valencies +/-z. As many electrolyte solutions in biological applications and hydrodynamical applications are not symmetric, the presented extension of the previous existence results is an important step.

preprint2016arXiv

Global existence of weak solutions of a model for electrolyte solutions - Part 2: Multicomponent case

This paper analytically investigates the Darcy-Poisson-Nernst-Planck system. This system is a mathematical model for electrolyte solutions. In this paper, we consider electrolyte solutions, which consist of a neutral fluid and multiple suspended charged chemical species with arbitrary valencies. We prove global existence and uniqueness of weak solutions in two space dimensions.

preprint2016arXiv

Hybrid Discretization Methods for Transient Numerical Simulation of Combustion in Porous Media

We present an algorithm for the numerical solution of the equations governing combustion in porous inert media. The discretization of the flow problem is performed by the mixed finite element method, the transport problems are discretized by a cell-centered finite volume method. The resulting nonlinear equations are lineararized with Newton's method, the linearized systems are solved with a multigrid algorithm. Both subsystems are recoupled again in a Picard iteration. Numerical simulations based on a simplified model show how regions with different porosity stabilize the reaction zone inside the porous burner.

preprint2016arXiv

Including van der Waals Forces in Diffusion-Convection Equations - Modeling, Analysis, and Numerical Simulations

This paper presents a model of van der Waals forces in the framework of diffusion-convection equations. The model consists of a nonlinear and degenerated diffusion-convection equation, which furthermore can be considered as a model for slow perikinetic coagulation. For the analytical investigation, we transform the model to a porous medium equation, which provides us access to the comprehensive analytical results for porous medium equations. Additionally, this transformation reveals a new application for porous medium equations. Eventually, we present numerical simulations of the model by solving the porous medium equation. We note that we solve the porous medium equation without any further regularization, which is often applied in this context.

preprint2016arXiv

Modeling and simulation of coagulation according to DLVO-theory in a continuum model for electrolyte solutions

This paper presents a model of coagulation in electrolyte solutions. In this paper, the coagulation process is modeled according to DLVO-theory, which is an atomistic theory. On the other hand, we describe the dynamics in the electrolyte solutions by the Poisson-Nernst-Planck system, which is a continuum model. The contribution of this paper is to include the atomistic description of coagulation based on DLVO-theory in the continuum Poisson-Nernst-Planck system. Thereby, we involve information from different spatial scales. For this reason, the presented model accounts for the short-range interactions and the long-range interactions, which drive the coagulation process. Furthermore, many-body effects are naturally included as the resulting model is a continuum model.

preprint2014arXiv

Multiscale Modeling of Colloidal Dynamics in Porous Media: Capturing Aggregation and Deposition Effects

We investigate the influence of multiscale aggregation and deposition on the colloidal dynamics in a saturated porous medium. At the pore scale, the aggregation of colloids is modeled by the Smoluchowski equation. Essentially, the colloidal mass is distributed between different size clusters. We treat these clusters as different species involved in a diffusion-advection-reaction mechanism. This modeling procedure allows for different material properties to be varied between the different species, specifically the rates of diffusion, aggregation, deposition as well as the advection velocities. We apply the periodic homogenization procedure to give insight into the effective coefficients of the upscaled model equations. Benefiting from direct access to microstructural information, we capture by means of 2D numerical simulations the effect of aggregation on the deposition rates recovering this way both the blocking and ripening regimes reported in the literature.

preprint2011arXiv

Rigorous Homogenization of a Stokes-Nernst-Planck-Poisson Problem for various Boundary Conditions

We perform the periodic homogenization (i.e. $\eps\to 0$) of the non-stationary Stokes-Nernst-Planck-Poisson system using two-scale convergence, where $\eps$ is a suitable scale parameter. The objective is to investigate the influence of \textsl{different boundary conditions and variable choices of scalings in $\eps$} of the microscopic system of partial differential equations on the structure of the (upscaled) limit model equations. Due to the specific nonlinear coupling of the underlying equations, special attention has to be paid when passing to the limit in the electrostatic drift term. As a direct result of the homogenization procedure, various classes of upscaled model equations are obtained.