Search

Thinking of flows crossing through regular porous media, we numerically explore the behavior of weak solutions to a two-scale elliptic-parabolic system that is strongly coupled by means of a suitable nonlinear dispersion term. The two-scale system of interest originates from the fast-drift periodic homogenization of a nonlinear convective-diffusion-reaction problem, where the structure of the non-linearity in the drift fits to the hydrodynamic limit of a totally asymmetric simple exclusion process for a population of particles. In this article, we focus exclusively on numerical simulations that employ two decoupled approximation schemes, viz. 'scheme 1' - a Picard-type iteration - and 'scheme 2' - a time discretization decoupling. Additionally, we describe a computational strategy which helps to drastically improve computation times. Finally, we provide several numerical experiments to illustrate what dispersion effects are introduced by a specific choice of microstructure and model ingredients., Comment: 27 pages, 10 figures, 3 tables

We investigate a two-scale system featuring an upscaled parabolic dispersion-reaction equation intimately linked to a family of elliptic cell problems. The system is strongly coupled through a dispersion tensor, which depends on the solutions to the cell problems, and via the cell problems themselves, where the solution of the parabolic problem interacts nonlinearly with the drift term. This particular mathematical structure is motivated by a rigorously derived upscaled reaction-diffusion-convection model that describes the evolution of a population of interacting particles pushed by a large drift through an array of periodically placed obstacles (i.e., through a regular porous medium). We prove the existence and uniqueness of weak solutions to our system by means of an iterative scheme, where particular care is needed to ensure the uniform positivity of the dispersion tensor. Additionally, we use finite element-based approximations for the same iteration scheme to perform multiple simulation studies. Finally, we highlight how the choice of micro-geometry (building the regular porous medium) and of the nonlinear drift coupling affects the macroscopic dispersion of particles.

We study the periodic homogenization of a reaction-diffusion problem with large nonlinear drift and Robin boundary condition posed in an unbounded perforated domain. The nonlinear problem is associated with the hydrodynamic limit of a totally asymmetric simple exclusion process (TASEP) governing a population of interacting particles crossing a domain with obstacle. We are interested in deriving rigorously the upscaled model equations and the corresponding effective coefficients for the case when the microscopic dynamics are linked to a particular choice of characteristic length and time scales that lead to an exploding nonlinear drift. The main mathematical difficulty lies in proving the two-scale compactness and strong convergence results needed for the passage to the homogenization limit. To cope with the situation, we use the concept of two-scale compactness with drift, which is similar to the more classical two-scale compactness result but it is defined now in moving coordinates. We provide as well a strong convergence result for the corrector function, starting this way the search for the order of the convergence rate of the homogenization process for our target nonlinear drift problem., Comment: 37 pages

We study a reaction-diffusion-convection problem with nonlinear drift posed in a domain with periodically arranged obstacles. The non-linearity in the drift is linked to the hydrodynamic limit of a totally asymmetric simple exclusion process (TASEP) governing a population of interacting particles crossing a domain with obstacle. Because of the imposed large drift scaling, this nonlinearity is expected to explode in the limit of a vanishing scaling parameter. As main working techniques, we employ two-scale formal homogenization asymptotics with drift to derive the corresponding upscaled model equations as well as the structure of the effective transport tensors. Finally, we use Schauder's fixed point theorem as well as monotonicity arguments to study the weak solvability of the upscaled model posed in an unbounded domain. This study wants to contribute with theoretical understanding needed when designing thin composite materials that are resistant to high velocity impacts., Comment: 27 pages, 3 figures

We study the question of periodic homogenization of a variably scaled reaction-diffusion problem with non-linear drift posed for a domain crossed by a flat composite thin layer. The structure of the non-linearity in the drift was obtained in earlier works as hydrodynamic limit of a totally asymmetric simple exclusion process (TASEP) process for a population of interacting particles crossing a domain with obstacle. Using energy-type estimates as well as concepts like thin-layer convergence and two-scale convergence, we derive the homogenized evolution equation and the corresponding effective model parameters for a regularized problem. Special attention is paid to the derivation of the effective transmission conditions across the separating limit interface in essentially two different situations: (i) finitely thin layer and (ii) infinitely thin layer. This study should be seen as a preliminary step needed for the investigation of averaging fast non-linear drifts across material interfaces -- a topic with direct applications in the design of thin composite materials meant to be impenetrable to high-velocity impacts., Comment: 41 pages, 5 figures