Browsing by Author "Rybak, Iryna"

Now showing 1 - 6 of 6

Open Access
Global sensitivity analysis using multi-resolution polynomial chaos expansion for coupled Stokes-Darcy flow problems
(2023) Kröker, Ilja; Oladyshkin, Sergey; Rybak, Iryna
Determination of relevant model parameters is crucial for accurate mathematical modelling and efficient numerical simulation of a wide spectrum of applications in geosciences. The conventional method of choice is the global sensitivity analysis (GSA). Unfortunately, at least the classical Monte-Carlo based GSA requires a high number of model runs. Response surfaces based techniques, e.g. arbitrary Polynomial Chaos (aPC) expansion, can reduce computational effort, however, they suffer from the Gibbs phenomena and high hardware requirements for higher accuracy. We introduce GSA for arbitrary Multi-Resolution Polynomial Chaos (aMR-PC) which is a localized aPC based data-driven polynomial discretization. The aMR-PC allows to reduce the Gibbs phenomena by construction and to achieve higher accuracy by means of localization also for lower polynomial degrees. We apply these techniques to perform the sensitivity analysis for the Stokes-Darcy problem which describes fluid flow in coupled free-flow and porous-medium systems. We consider the Stokes equations in the free-flow region, Darcy’s law in the porous-medium domain and the classical interface conditions across the fluid–porous interface including the conservation of mass, the balance of normal forces and the Beavers–Joseph condition for the tangential velocity. This coupled problem formulation contains four uncertain parameters: the exact location of the interface, the permeability, the Beavers-Joseph slip coefficient and the uncertainty in the boundary conditions. We carry out the sensitivity analysis of the coupled model with respect to these parameters using the Sobol indices on the aMR-PC expansion and conduct the corresponding numerical simulations.
Open Access
Mathematical modeling of coupled free flow and porous medium systems
(2016) Rybak, Iryna; Rohde, Christian (Prof. Dr.)
Different classes of physical systems with common interfaces arise in a variety of environmental and industrial problems. Striking examples originate from terrestrial-atmospheric contact zones, surface water-groundwater interaction, filters and fuel cells, where a free fluid system is in contact with a porous medium. Flow and transport processes in these systems evolve on multiple length and time scales, contributing to the complexity of these systems both from the modeling and the numerical side. An additional contributing factor to this complexity is the existence of multiple classes of entities, which include phases, interfaces between phases, and common curves that form at the boundary between three phases. Classical coupling approaches and traditional porous medium model formulations lead to reliable results in limited cases, however, applications require more realistic settings. The focus of the thesis is on derivation of mathematical models for multiphase multi-component porous medium systems that take into account lower dimensional entities, formulation of coupling conditions at the sharp interface and the transition region between porous medium and free flow systems, computation of effective parameters for the macroscale models, development and analysis of efficient numerical algorithms for coupled problems, and numerical simulation of applications.
Open Access
A modification of the Beavers-Joseph condition for arbitrary flows to the fluid-porous interface
(2023) Strohbeck, Paula; Eggenweiler, Elissa; Rybak, Iryna
Physically consistent coupling conditions at the fluid-porous interface with correctly determined effective parameters are necessary for accurate modeling and simulation of various applications. To describe single-fluid-phase flows in coupled free-flow and porous-medium systems, the Stokes/Darcy equations are typically used together with the conservation of mass across the interface, the balance of normal forces and the Beavers-Joseph condition on the tangential velocity. The latter condition is suitable for flows parallel to the interface but not applicable for arbitrary flow directions. Moreover, the value of the Beavers-Joseph slip coefficient is uncertain. In the literature, it is routinely set equal to one that is not correct for many applications, even if the flow is parallel to the porous layer. In this paper, we reformulate the generalized interface condition on the tangential velocity component, recently developed for arbitrary flows in Stokes/Darcy systems, such that it has the same analytical form as the Beavers-Joseph condition. We compute the effective coefficients appearing in this modified condition using theory of homogenization with boundary layers. We demonstrate that the modified Beavers-Joseph condition is applicable for arbitrary flow directions to the fluid-porous interface. In addition, we propose an efficient two-level numerical algorithm based on simulated annealing to compute the optimal Beavers-Joseph parameter. Article Highlights A modification of the Beavers-Joseph condition is proposed based on recently developed generalized coupling conditions. The Beavers-Joseph parameter can be found only for unidirectional flows. An efficient numerical algorithm to determine the optimal Beavers-Joseph parameter is developed.
Open Access
Permeability estimation of regular porous structures : a benchmark for comparison of methods
(2021) Wagner, Arndt; Eggenweiler, Elissa; Weinhardt, Felix; Trivedi, Zubin; Krach, David; Lohrmann, Christoph; Jain, Kartik; Karadimitriou, Nikolaos; Bringedal, Carina; Voland, Paul; Holm, Christian; Class, Holger; Steeb, Holger; Rybak, Iryna
The intrinsic permeability is a crucial parameter to characterise and quantify fluid flow through porous media. However, this parameter is typically uncertain, even if the geometry of the pore structure is available. In this paper, we perform a comparative study of experimental, semi-analytical and numerical methods to calculate the permeability of a regular porous structure. In particular, we use the Kozeny-Carman relation, different homogenisation approaches (3D, 2D, very thin porous media and pseudo 2D/3D), pore-scale simulations (lattice Boltzmann method, Smoothed Particle Hydrodynamics and finite-element method) and pore-scale experiments (microfluidics). A conceptual design of a periodic porous structure with regularly positioned solid cylinders is set up as a benchmark problem and treated with all considered methods. The results are discussed with regard to the individual strengths and limitations of the used methods. The applicable homogenisation approaches as well as all considered pore-scale models prove their ability to predict the permeability of the benchmark problem. The underestimation obtained by the microfluidic experiments is analysed in detail using the lattice Boltzmann method, which makes it possible to quantify the influence of experimental setup restrictions.
Open Access
A surrogate-assisted uncertainty-aware Bayesian validation framework and its application to coupling free flow and porous-medium flow
(2023) Mohammadi, Farid; Eggenweiler, Elissa; Flemisch, Bernd; Oladyshkin, Sergey; Rybak, Iryna; Schneider, Martin; Weishaupt, Kilian
Existing model validation studies in geoscience often disregard or partly account for uncertainties in observations, model choices, and input parameters. In this work, we develop a statistical framework that incorporates a probabilistic modeling technique using a fully Bayesian approach to perform a quantitative uncertainty-aware validation. A Bayesian perspective on a validation task yields an optimal bias-variance trade-off against the reference data. It provides an integrative metric for model validation that incorporates parameter and conceptual uncertainty. Additionally, a surrogate modeling technique, namely Bayesian Sparse Polynomial Chaos Expansion, is employed to accelerate the computationally demanding Bayesian calibration and validation. We apply this validation framework to perform a comparative evaluation of models for coupling a free flow with a porous-medium flow. The correct choice of interface conditions and proper model parameters for such coupled flow systems is crucial for physically consistent modeling and accurate numerical simulations of applications. We develop a benchmark scenario that uses the Stokes equations to describe the free flow and considers different models for the porous-medium compartment and the coupling at the fluid-porous interface. These models include a porous-medium model using Darcy’s law at the representative elementary volume scale with classical or generalized interface conditions and a pore-network model with its related coupling approach. We study the coupled flow problems’ behaviors considering a benchmark case, where a pore-scale resolved model provides the reference solution. With the suggested framework, we perform sensitivity analysis, quantify the parametric uncertainties, demonstrate each model’s predictive capabilities, and make a probabilistic model comparison.
Open Access
Validation and calibration of coupled porous-medium and free-flow problems using pore-scale resolved models
(2020) Rybak, Iryna; Schwarzmeier, Christoph; Eggenweiler, Elissa; Rüde, Ulrich
The correct choice of interface conditions and effective parameters for coupled macroscale free-flow and porous-medium models is crucial for a complete mathematical description of the problem under consideration and for accurate numerical simulation of applications. We consider single-fluid-phase systems described by the Stokes–Darcy model. Different sets of coupling conditions for this model are available. However, the choice of these conditions and effective model parameters is often arbitrary. We use large-scale lattice Boltzmann simulations to validate coupling conditions by comparison of the macroscale simulations against pore-scale resolved models. We analyse three settings (lid-driven cavity over a porous bed, infiltration problem and general filtration problem) with different geometrical configurations (channelised and staggered distributions of solid grains) and different sets of interface conditions. Effective parameters for the macroscale models (permeability tensor, boundary layer constants) are computed numerically for each geometrical configuration. Numerical simulation results demonstrate the sensitivity of the coupled Stokes–Darcy problem to the location of the sharp fluid–porous interface, the effective model parameters and the interface conditions.