BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Adrian Muntean

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

ResearcherUnclaimed source record
math.APmath-phmath.MPcond-mat.stat-mechphysics.soc-phnlin.CGmath.NAmath.FAmath.PRNumerical Analysiscond-mat.mtrl-scicond-mat.softmath.DSnlin.AOphysics.chem-phphysics.comp-phSocial and Information Networks

Catalog footprint

What is connected

47works
17topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

47 published item(s)

preprint2022arXiv

A reduction scheme for coupled Brownian harmonic oscillators

We propose a reduction scheme for a system constituted by two coupled harmonically-bound Brownian oscillators. We reduce the description by constructing a lower dimensional model which inherits some of the basic features of the original dynamics and is written in terms of suitable transport coefficients. The proposed procedure is twofold: while the deterministic component of the dynamics is obtained by a direct application of the invariant manifold method, the diffusion terms are determined via the Fluctuation-Dissipation Theorem. We highlight the behavior of the coefficients up to a critical value of the coupling parameter, which marks the endpoint of the interval in which a contracted description is available. The study of the weak coupling regime is addressed and the commutativity of alternative reduction paths is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Coupled stochastic systems of Skorokhod type: well-posedness of a mathematical model and its applications

Population dynamics with complex biological interactions, accounting for uncertainty quantification, is critical for many application areas. However, due to the complexity of biological systems, the mathematical formulation of the corresponding problems faces the challenge that the corresponding stochastic processes should, in most cases, be considered in bounded domains. We propose a model based on a coupled system of reflecting Skorokhod-type stochastic differential equations with jump-like exit from a boundary. The setting describes the population dynamics of active and passive populations. As main working techniques, we use compactness methods and Skorokhod's representation of solutions to SDEs posed in bounded domains to prove the well-posedness of the system. This functional setting is a new point of view in the field of modelling and simulation of population dynamics. We provide the details of the model, as well as representative numerical examples, and discuss the applications of a Wilson-Cowan-type system, modelling the dynamics of two interacting populations of excitatory and inhibitory neurons. Furthermore, the presence of random input current, reflecting factors together with Poisson jumps, increases firing activity in neuronal systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Error estimates for semi-discrete finite element approximations for a moving boundary problem capturing the penetration of diffusants into rubber

We consider a moving boundary problem with kinetic condition that describes the diffusion of solvent into rubber and study semi-discrete finite element approximations of the corresponding weak solutions. We report on both a priori and a posteriori error estimates for the mass concentration of the diffusants, and respectively, for the a priori unknown position of the moving boundary. Our working techniques include integral and energy-based estimates for a nonlinear parabolic problem posed in a transformed fixed domain combined with a suitable use of the interpolation-trace inequality to handle the interface terms. Numerical illustrations of our FEM approximations are within the experimental range and show good agreement with our theoretical investigation. This work is a preliminary investigation necessary before extending the current moving boundary modeling to account explicitly for the mechanics of hyperelastic rods to capture a directional swelling of the underlying elastomer.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Identifying Processes Governing Damage Evolution in Quasi-Static Elasticity Part I -- Analysis

We present a quasi-static elasticity model that accounts for damage evolution based on the ideas of Kachanov 1958 and and Rabotnov 1968. We analyze the resulting strongly nonlinear system of differential equations in view of well-posedness. The specific feature is that the displacements are connected to the damage evolution equation via a Nemytskii- or superposition- operator. The novelty in this work is that we present an inverse problem in a parameter identification setting, in which we are able to identify the shape of this Nemytskii-operator. From the material modelling point of view, the relation of material damage and displacements is modelled by the shape of this operator. We establish the Fréchet-derivative of the forward operator as well as the adjoint of the derivative and characterize both via systems of linear differential equations. We prove ill-posedness of the inverse problem and provide a sufficient condition for the classical nonlinear Landweber method to converge.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Identifying Processes Governing Damage Evolution in Quasi-Static Elasticity. Part 2 -- Numerical Simulations

We investigate numerically a quasi-static elasticity system of Kachanov-type. To do so we propose an Euler time discretization combined with a suitable finite elements scheme (FEM) to handle the discretization is space. We use ODE-type arguments to prove the consistency of the scheme as well as its convergence rate. We rely on the computational platform FEniCS to perform the FEM discretizations in space needed to compute the model output. The simulation results show a good agreement with both the physics of the problem and with our previous qualitative mathematical analysis results obtained for precisely the same problem setting. Furthermore, our implementation recovers nicely the theoretically expected convergence rate. This is a preliminary study preparing the framework for the rigorous numerical identification of the damage process in Kachanov-type models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Local existence of a solution to a free boundary problem describing migration into rubber with a breaking effect

We consider a one-dimensional free boundary problem describing the migration of diffusants into rubber. In our setting, the free boundary represents the position of the front delimitating the diffusant region. The growth rate of this region is described by an ordinary differential equation that includes the effect of breaking the growth of the diffusant region. In this specific context, the breaking mechanism is assumed to be the hyperelastic response to a too fast diffusion penetration. In recent works, we considered a similar class of free boundary problems modeling diffusants penetration in rubbers, but without attempting to deal with the possibility of breaking or accelerating the occurring free boundaries. For simplified settings, we were able to show the global existence and uniqueness as well as the large-time behavior of the corresponding solutions to our formulations. Since here the breaking effect is contained in the free boundary condition, our previous results are not anymore applicable. The main mathematical obstacle in ensuring the existence of a solution is the non-monotonic structure of the free boundary. In this paper, we establish the existence and uniqueness of a weak solution to the free boundary problem with breaking effect and give explicitly the maximum value that the free boundary can reach.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reduced Markovian Descriptions of Brownian Dynamics: Toward an Exact Theory

We outline a reduction scheme for a class of Brownian dynamics which leads to meaningful corrections to the Smoluchowski equation in the overdamped regime. The mobility coefficient of the reduced dynamics is obtained by exploiting the Dynamic Invariance principle, whereas the diffusion coefficient fulfils the Fluctuation-Dissipation theorem. Explicit calculations are carried out in the case of a harmonically bound particle. A quantitative pointwise representation of the reduction error is also provided and connections to both the Maximum Entropy method and the linear response theory are highlighted. Our study paves the way to the development of reduction procedures applicable to a wider class of diffusion processes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Upscaling of a reaction-diffusion-convection problem with exploding non-linear drift

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Well-posedness of a coupled system of Skorohod-like stochastic differential equations

We study the well-posedness of a coupled system of Skorohod-like stochastic differential equations with reflecting boundary condition. The setting describes the evacuation dynamics of a mixed crowd composed of both active and passive pedestrians moving through a domain with obstacles, fire and smoke. As main working techniques, we use compactness methods and the Skorohod's representation of solutions to SDEs posed in bounded domains. This functional setting is a new point of view in the field of modeling and simulation pedestrian dynamics. The main challenge is to handle the coupling in the model equations together with the multiple-connectedness of the domain and the pedestrian-obstacle interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Anticipation decides on lane formation in pedestrian counterflow -- a simulation study

Human crowds base most of their behavioral decisions upon anticipated states of their walking environment. We explore a minimal version of a lattice model to study lanes formation in pedestrian counterflow. Using the concept of horizon depth, our simulation results suggest that the anticipation effect together with the presence of a small background noise play an important role in promoting collective behaviors in a counterflow setup. These ingredients facilitate the formation of seemingly stable lanes and ensure the ergodicity of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Deterministic reversible model of non-equilibrium phase transitions and stochastic counterpart

N point particles move within a billiard table made of two circular cavities connected by a straight channel. The usual billiard dynamics is modified so that it remains deterministic, phase space volumes preserving and time reversal invariant. Particles move in straight lines and are elastically reflected at the boundary of the table, as usual, but those in a channel that are moving away from a cavity invert their motion (rebound), if their number exceeds a given threshold T. When the geometrical parameters of the billiard table are fixed, this mechanism gives rise to non--equilibrium phase transitions in the large N limit: letting T/N decrease, the homogeneous particle distribution abruptly turns into a stationary inhomogeneous one. The equivalence with a modified Ehrenfest two urn model, motivated by the ergodicity of the billiard with no rebound, allows us to obtain analytical results that accurately describe the numerical billiard simulation results. Thus, a stochastic exactly solvable model that exhibits non-equilibrium phase transitions is also introduced.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Global weak solvability, Continuous dependence on data and large time growth of swelling moving interfaces

We prove a global existence result for weak solutions to a one-dimensional free boundary problem with flux boundary conditions describing swelling along a halfline. Additionally, we show that solutions are not only unique but also depend continuously on data and parameters. The key observation is that the structure of our system of partial differential equations allows us to show that the moving a priori unknown interface never disappears. As main ingredients of the global existence proof, we rely on a local weak solvability result for our problem, uniform estimates of the solution, integral estimates on quantities defined at the free boundary, as well as a fine pointwise lower bound for the position of the moving boundary. Some of the estimates are time-independent. They allow us to explore the large time behavior of the position of the moving boundary. The approach is specific to one-dimensional settings.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

When diffusion faces drift: consequences of exclusion processes for bi--directional pedestrian flows

Stochastic particle--based models are useful tools for describing the collective movement of large crowds of pedestrians in crowded confined environments. Using descriptions based on the simple exclusion process, two populations of particles, mimicking pedestrians walking in a built environment, enter a room from two opposite sides. One population is passive -- being unaware of the local environment; particles belonging to this group perform a symmetric random walk. The other population has information on the local geometry in the sense that as soon as particles enter a visibility zone, a drift activates them. Their self-propulsion leads them towards the exit. This second type of species is referred here as active. The assumed crowdedness corresponds to a near--jammed scenario. The main question we ask in this paper is: Can we induce modifications of the dynamics of the active particles to improve the outgoing current of the passive particles? To address this question, we compute occupation number profiles and currents for both populations in selected parameter ranges. Besides observing the more classical faster--is--slower effect, new features appear as prominent like the non--monotonicity of currents, self--induced phase separation within the active population, as well as acceleration of passive particles for large--drift regimes of active particles.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A semidiscrete Galerkin scheme for a coupled two-scale elliptic-parabolic system: well-posedness and convergence approximation rates

In this paper, we study the numerical approximation of a coupled system of elliptic-parabolic equations posed on two separated spatial scales. The model equations describe the interplay between macroscopic and microscopic pressures in an unsaturated heterogeneous medium with distributed microstructures as they often arise in modeling reactive flow in cementitious-based materials. Besides ensuring the well-posedness of our two-scale model, we design two-scale convergent numerical approximations and prove a priori error estimates for the semidiscrete case. We complement our analysis with simulation results illustrating the expected behaviour of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A lattice model of reduced jamming by barrier

We study an asymmetric simple exclusion process in a strip in the presence of a solid impenetrable barrier. We focus on the effect of the barrier on the residence time of the particles, namely, the typical time needed by the particles to cross the whole strip. We explore the conditions for reduced jamming when varying the environment (different drifts, reservoir densities, horizontal diffusion walks, etc.). Particularly, we discover an interesting non--monotonic behavior of the residence time as a function of the barrier length. Besides recovering by means of both the lattice dynamics and mean-field model well-known aspects like faster-is-slower effect and the intermittence of the flow, we propose also a birth-and-death process and a reduced one-dimensional model with variable barrier permeability to capture qualitatively the behavior of the residence time with respect to the parameters. We report our first steps towards the understanding to which extent the presence of obstacles can fluidize pedestrian and biological transport in crowded heterogeneous environments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Note on Iterations-based Derivations of High-order Homogenization Correctors for Multiscale Semi-linear Elliptic Equations

This Note aims at presenting a simple and efficient procedure to derive the structure of high-order corrector estimates for the homogenization limit applied to a semi-linear elliptic equation posed in perforated domains. Our working technique relies on monotone iterations combined with formal two-scale homogenization asymptotics. It can be adapted to handle more complex scenarios including for instance nonlinearities posed at the boundary of perforations and the vectorial case, when the model equations are coupled only through the nonlinear production terms.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A priori feedback estimates for multiscale reaction-diffusion systems

We study the approximation of a multiscale reaction-diffusion system posed on both macroscopic and microscopic space scales. The coupling between the scales is done via micro-macro flux conditions. Our target system has a typical structure for reaction-diffusion-flow problems in media with distributed microstructures (also called, double porosity materials). Besides ensuring basic estimates for the convergence of two-scale semi-discrete Galerkin approximations, we provide a set of {\em a priori} feedback estimates and a local feedback error estimator that help in designing a distributed-high-errors strategy to allow for a computationally efficient zooming in and out from microscopic structures. The error control on the feedback estimates relies on two-scale-energy, regularity, and interpolation estimates as well as on a fine bookeeping of the sources responsible with the propagation of the (multiscale) approximation errors. The working technique based on {\em a priori } feedback estimates is in principle applicable to a large class of systems of PDEs with dual structure admitting strong solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Blockage induced condensation controlled by a local reaction

We consider the set-up of stationary Zero Range models and discuss the onset of condensation induced by a local blockage on the lattice. We show that the introduction of a local feedback on the hopping rates allows to control the particle fraction in the condensed phase. This phenomenon results in a current vs. blockage parameter curve characterized by two non--analyticity points.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Large-time behavior of solutions to a thermo-diffusion system with Smoluchowski interactions

We prove the large time behavior of solutions to a coupled thermo-diffusion arising in the modelling of the motion of hot colloidal particles in porous media. Additionally, we also ensure the uniqueness of solutions of the target problem. The main mathematical difficulty is due to the presence in the right-hand side of the equations of products between temperature and concentration gradients. Such terms mimic the so-called thermodynamic Soret and Dufour effects. These are cross-coupling terms emphasizing in this context a strong interplay between heat conduction and molecular diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Measure-valued mass evolution problems with flux boundary conditions and solution-dependent velocities

In this paper we prove well-posedness for a measure-valued continuity equation with solution-dependent velocity and flux boundary conditions, posed on a bounded one-dimensional domain. We generalize the results of [Evers, Hille and Muntean. Journal of Differential Equations, 259:1068-1097, 2015] to settings where the dynamics are driven by interactions. In a forward-Euler-like approach, we construct a time-discretized version of the original problem and employ the results of [Evers, Hille and Muntean. Journal of Differential Equations, 259:1068-1097, 2015] as a building block within each subinterval. A limit solution is obtained as the mesh size of the time discretization goes to zero. Moreover, the limit is independent of the specific way of partitioning the time interval $[0,T]$. This paper is partially based on results presented in Chapter 5 of [Evers. PhD thesis, Eindhoven University of Technology, 2015], while a number of issues that were still open there, are now resolved.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Asymmetric pedestrian dynamics on a staircase landing from continuous measurements

We investigate via extensive experimental data the dynamics of pedestrians walking in a corridor-shaped landing in a building at Eindhoven University of Technology. With year-long automatic measurements employing a Microsoft KinectTM 3D-range sensor and ad hoc tracking techniques, we acquired few hundreds of thousands pedestrian trajectories in real-life conditions. Here we discuss the asymmetric features of the dynamics in the two walking directions with respect to the flights of stairs (i.e. ascending or descending). We provide a detailed analysis of position and speed fields for the cases of pedestrians walking alone undisturbed and for couple of pedestrians in counter-flow. Then, we show average walking velocities exploring all the observed combinations in terms of numbers of pedestrians and walking directions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Asymptotic analysis of a semi-linear elliptic system in perforated domains: well-posedness and correctors for the homogenization limit

In this study, we prove results on the weak solvability and homogenization of a microscopic semi-linear elliptic system posed in perforated media. The model presented here explores the interplay between stationary diffusion and both surface and volume chemical reactions in porous media. Our interest lies in deriving homogenization limits (upscaling) for alike systems and particularly in justifying rigorously the obtained averaged descriptions. Essentially, we prove the well-posedness of the microscopic problem ensuring also the positivity and boundedness of the involved concentrations and then use the structure of the two scale expansions to derive corrector estimates delimitating this way the convergence rate of the asymptotic approximates to the macroscopic limit concentrations. Our techniques include Moser-like iteration techniques, a variational formulation, two-scale asymptotic expansions as well as energy-like estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Residence time estimates for asymmetric simple exclusion dynamics on stripes

The target of our study is to approximate numerically and, in some particular physically relevant cases, also analytically, the residence time of particles undergoing an asymmetric simple exclusion dynamics on a stripe. The source of asymmetry is twofold: (i) the choice of boundary conditions (different reservoir levels) and (ii) the strong anisotropy from a nonlinear drift with prescribed directionality. We focus on the effect of the choice of anisotropy in the flux on the asymptotic behavior of the residence time with respect to the length of the stripe. The topic is relevant for situations occurring in pedestrian flows or biological transport in crowded environments, where lateral displacements of the particles occur predominantly affecting therefore in an essentially way the efficiency of the overall transport mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A thermo-diffusion system with Smoluchowski interactions: well-posedness and homogenization

We study the solvability and homogenization of a thermal-diffusion reaction problem posed in a periodically perforated domain. The system describes the motion of populations of hot colloidal particles interacting together via Smoluchowski production terms. The upscaled system, obtained via two-scale convergence techniques, allows the investigation of deposition effects in porous materials in the presence of thermal gradients.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Error control for the FEM approximation of an upscaled thermo-diffusion system with Smoluchowski interactions

We analyze a coupled system of evolution equations that describes the effect of thermal gradients on the motion and deposition of $N$ populations of colloidal species diffusing and interacting together through Smoluchowski production terms. This class of systems is particularly useful in studying drug delivery, contaminant transportin complex media, as well as heat shocks thorough permeable media. The particularity lies in the modeling of the nonlinear and nonlocal coupling between diffusion and thermal conduction. We investigate the semidiscrete as well as the fully discrete em a priori error analysis of the finite elements approximation of the weak solution to a thermo-diffusion reaction system posed in a macroscopic domain. The mathematical techniques include energy-like estimates and compactness arguments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mild solutions to a measure-valued mass evolution problem with flux boundary conditions

We investigate the well-posedness and approximation of mild solutions to a class of linear transport equations on the unit interval $[0,1]$ endowed with a linear discontinuous production term, formulated in the space $\mathcal{M}([0,1])$ of finite Borel measures. Our working technique includes a detailed boundary layer analysis in terms of a semigroup representation of solutions in spaces of measures able to cope with the passage to the singular limit where thickness of the layer vanishes. We obtain not only a suitable concept of solutions to the chosen measure-valued evolution problem, but also derive convergence rates for the approximation procedure and get insight in the structure of flux boundary conditions for the limit problem.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modelling with measures: Approximation of a mass-emitting object by a point source

We consider a linear diffusion equation on $Ω:=\mathbb{R}^2\setminus\bar{Ω_\mathcal{O}}$, where $Ω_\mathcal{O}$ is a bounded domain. The time-dependent flux on the boundary $Γ:=\partialΩ_\mathcal{O}$ is prescribed. The aim of the paper is to approximate the dynamics by the solution of the diffusion equation on the whole of $\mathbb{R}^2$ with a measure-valued point source in the origin and provide estimates for the quality of approximation. For all time $t$, we derive an $L^2([0,t];L^2(Γ))$-bound on the difference in flux on the boundary. Moreover, we derive for all $t>0$ an $L^2(Ω)$-bound and an $L^2([0,t];H^1(Ω))$-bound for the difference of the solutions to the two models.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Upscaling of the dynamics of dislocation walls

We perform the discrete-to-continuum limit passage for a microscopic model describing the time evolution of dislocations in a one dimensional setting. This answers the related open question raised by Geers et al. in [GPPS13]. The proof of the upscaling procedure (i.e. the discrete-to-continuum passage) relies on the gradient flow structure of both the discrete and continuous energies of dislocations set in a suitable evolutionary variational inequality framework. Moreover, the convexity and $Γ$-convergence of the respective energies are properties of paramount importance for our arguments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Crowds reaching targets by maximizing entropy: a Clausius-Duhem inequality approach

In this paper we propose the use of concepts from thermodynamics in the study of crowd dynamics. Our continuous model consists of the continuity equation for the density of the crowd and a kinetic equation for the velocity field. The latter includes a nonlocal term that models interactions between individuals. To support our modelling assumptions, we introduce an inequality that resembles the Second Law of Thermodynamics, containing an entropy-like functional. We show that its time derivative equals a positive dissipation term minus a corrector term. The latter term should be small for the time derivative of the entropy to be positive. In case of isotropic interactions the corrector term is absent. For the anisotropic case, we support the claim that the corrector term is small by simulations for the corresponding particle system. They reveal that this term is sufficiently small for the entropy still to increase. Moreover, we show that the entropy converges in time towards a limit value.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Is adding charcoal to soil a good method for CO2 sequestration? - Modeling a spatially homogeneous soil

Carbon sequestration is the process of capture and long-term storage of atmospheric carbon dioxide (CO2) with the aim to avoid dangerous climate change. In this paper, we propose a simple mathematical model (a coupled system of nonlinear ODEs) to capture some of the dynamical effects produced by adding charcoal to fertile soils. The main goal is to understand to which extent charcoal is able to lock up carbon in soils. Our results are preliminary in the sense that we do not solve the CO2 sequestration problem. Instead, we do set up a flexible modeling framework in which the interaction between charcoal and soil can be tackled by means of mathematical tools. We show that our model is well-posed and has interesting large-time behaviour. Depending on the reference parameter range (e.g. type of soil) and chosen time scale, numerical simulations suggest that adding charcoal typically postpones the release of CO2.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Large-time asymptotics of moving-reaction interfaces involving nonlinear Henry's law and time-dependent Dirichlet data

We study the large-time behavior of the free boundary position capturing the one-dimensional motion of the carbonation reaction front in concrete-based materials. We extend here our rigorous justification of the $\sqrt{t}$-behavior of reaction penetration depths by including non-linear effects due to deviations from the classical Henry's law and time-dependent Dirichlet data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Large-time behavior of a two-scale semilinear reaction-diffusion system for concrete sulfatation

We study the large-time behavior of (weak) solutions to a two-scale reaction-diffusion system coupled with a nonlinear ordinary differential equations modeling the partly dissipative corrosion of concrete (/cement)-based materials with sulfates. We prove that as $t\to\infty$ the solution to the original two-scale system converges to the corresponding two-scale stationary system. To obtain the main result we make use essentially of the theory of evolution equations governed by subdifferential operators of time-dependent convex functions developed combined with a series of two-scale energy-like time-independent estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pedestrians moving in dark: Balancing measures and playing games on lattices

We present two conceptually new modeling approaches aimed at describing the motion of pedestrians in obscured corridors: * a Becker-Döring-type dynamics * a probabilistic cellular automaton model. In both models the group formation is affected by a threshold. The pedestrians are supposed to have very limited knowledge about their current position and their neighborhood; they can form groups up to a certain size and they can leave them. Their main goal is to find the exit of the corridor. Although being of mathematically different character, the discussion of both models shows that it seems to be a disadvantage for the individual to adhere to larger groups. We illustrate this effect numerically by solving both model systems. Finally we list some of our main open questions and conjectures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The effect of perception anisotropy on particle systems describing pedestrian flows in corridors

We consider a microscopic model (a system of self-propelled particles) to study the behaviour of a large group of pedestrians walking in a corridor. Our point of interest is the effect of anisotropic interactions on the global behaviour of the crowd. The anisotropy we have in mind reflects the fact that people do not perceive (i.e. see, hear, feel or smell) their environment equally well in all directions. The dynamics of the individuals in our model follow from a system of Newton-like equations in the overdamped limit. The instantaneous velocity is modelled in such a way that it accounts for the angle under which an individual perceives another individual. We investigate the effects of this perception anisotropy by means of simulations, very much in the spirit of molecular dynamics. We define a number of characteristic quantifiers (including the polarization index and Morisita index) that serve as measures for e.g. organization and clustering, and we use these indices to investigate the influence of anisotropy on the global behaviour of the crowd. The goal of the paper is to investigate the potentiality of this model; extensive statistical analysis of simulation data, or reproducing any specific real-life situation are beyond its scope.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Can cooperation slow down emergency evacuations?

We study the motion of pedestrians through obscure corridors where the lack of visibility hides the precise position of the exits. Using a lattice model, we explore the effects of cooperation on the overall exit flux (evacuation rate). More precisely, we study the effect of the buddying threshold (of no--exclusion per site) on the dynamics of the crowd. In some cases, we note that if the evacuees tend to cooperate and act altruistically, then their collective action tends to favor the occurrence of disasters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Distributed space scales in a semilinear reaction-diffusion system including a parabolic variational inequality: A well-posedness study

This paper treats the solvability of a semilinear reaction-diffusion system, which incorporates transport (diffusion) and reaction effects emerging from two separated spatial scales: $x$ - macro and $y$ - micro. The system's origin connects to the modeling of concrete corrosion in sewer concrete pipes. It consists of three partial differential equations which are mass-balances of concentrations, as well as, one ordinary differential equation tracking the damage-by-corrosion. The system is semilinear, partially dissipative, and coupled via the solid-water interface at the microstructure (pore) level. The structure of the model equations is obtained in \cite{tasnim1} by upscaling of the physical and chemical processes taking place within the microstructure of the concrete. Herein we ensure the positivity and $L^\infty-$bounds on concentrations, and then prove the global-in-time existence and uniqueness of a suitable class of positive and bounded solutions that are stable with respect to the two-scale data and model parameters. The main ingredient to prove existence include fixed-point arguments and convergent two-scale Galerkin approximations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamics of pedestrians in regions with no visibility - a lattice model without exclusion

We investigate the motion of pedestrians through obscure corridors where the lack of visibility (due to smoke, fog, darkness, etc.) hides the precise position of the exits. We focus our attention on a set of basic mechanisms, which we assume to be governing the dynamics at the individual level. Using a lattice model, we explore the effects of non-exclusion on the overall exit flux (evacuation rate). More precisely, we study the effect of the buddying threshold (of no-exclusion per site) on the dynamics of the crowd and investigate to which extent our model confirms the following pattern revealed by investigations on real emergencies: If the evacuees tend to cooperate and act altruistically, then their collective action tends to favor the occurrence of disasters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Macroscopic corrosion front computations of sulfate attack in sewer pipes based on a micro-macro reaction-diffusion model

We consider a two-scale reaction diffusion system able to capture the corrosion of concrete with sulfates. Our aim here is to define and compute two macroscopic corrosion indicators: typical pH drop and gypsum profiles. Mathematically, the system is coupled, endowed with micro-macro transmission conditions, and posed on two different spatially-separated scales: one microscopic (pore scale) and one macroscopic (sewer pipe scale). We use a logarithmic expression to compute values of pH from the volume averaged concentration of sulfuric acid which is obtained by resolving numerically the two-scale system (microscopic equations with direct feedback with the macroscopic diffusion of one of the reactants). Furthermore, we also evaluate the content of the main sulfatation reaction (corrosion) product---the gypsum---and point out numerically a persistent kink in gypsum's concentration profile. Finally, we illustrate numerically the position of the free boundary separating corroded from not-yet-corroded regions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Unfolding-based corrector estimates for a reaction-diffusion system predicting concrete corrosion

We use the periodic unfolding technique to derive corrector estimates for a reaction-diffusion system describing concrete corrosion penetration in the sewer pipes. The system, defined in a periodically-perforated domain, is semi-linear, partially dissipative, and coupled via a non-linear ordinary differential equation posed on the solid-water interface at the pore level. After discussing the solvability of the pore scale model, we apply the periodic unfolding techniques (adapted to treat the presence of perforations) not only to get upscaled model equations, but also to prepare a proper framework for getting a convergence rate (corrector estimates) of the averaging procedure.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A free-boundary problem for concrete carbonation: Rigorous justification of the $\sqrt{t}$-law of propagation

We study a one-dimensional free-boundary problem describing the penetration of carbonation fronts (free reaction-triggered interfaces) in concrete. A couple of decades ago, it was observed experimentally that the penetration depth versus time curve (say $s(t)$ vs. $t$) behaves like $s(t)=C\sqrt{t}$ for sufficiently large times $t > 0$ (with $C$ a positive constant). Consequently, many fitting arguments solely based on this experimental law were used to predict the large-time behavior of carbonation fronts in real structures, a theoretical justification of the $\sqrt{t}$-law being lacking until now. %This is the place where our paper contributes: The aim of this paper is to fill this gap by justifying rigorously the experimentally guessed asymptotic behavior. We have previously proven the upper bound $s(t)\leq C'\sqrt{t}$ for some constant $C'$; now we show the optimality of the rate by proving the right nontrivial lower estimate, i.e. there exists $C">0$ such that $s(t)\geq C"\sqrt{t}$. Additionally, we obtain weak solutions to the free-boundary problem for the case when the measure of the initial domain vanishes. In this way, our mathematical model is now allowing for the appearance of a moving carbonation front -- a scenario that until was hard to handle from the analysis point of view.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Corrector estimates for the homogenization of a locally-periodic medium with areas of low and high diffusivity

We prove an upper bound for the convergence rate of the homogenization limit $ε\to 0$ for a linear transmission problem for a advection-diffusion(-reaction) system posed in areas with low and high diffusivity, where $ε$ is a suitable scale parameter. On this way, we justify the formal homogenization asymptotics obtained by us earlier by proving an upper bound for the convergence rate (a corrector estimate). The main ingredients of the proof of the corrector estimate include integral estimates for rapidly oscillating functions with prescribed average, properties of the macroscopic reconstruction operators, energy bounds and extra two-scale regularity estimates. The whole procedure essentially relies on a good understanding of the analysis of the limit two-scale problem.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Semi-discrete finite difference multiscale scheme for a concrete corrosion model: approximation estimates and convergence

We propose a semi-discrete finite difference multiscale scheme for a concrete corrosion model consisting of a system of two-scale reaction-diffusion equations coupled with an ode. We prove energy and regularity estimates and use them to get the necessary compactness of the approximation estimates. Finally, we illustrate numerically the behavior of the two-scale finite difference approximation of the weak solution.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Single and two-scale sharp-interface models for concrete carbonation -- Asymptotics and numerical approximation

We investigate the fast-reaction asymptotics for a one-dimensional reaction-diffusion (RD) system describing the penetration of the carbonation reaction in concrete. The technique of matched-asymptotics is used to show that the RD system leads to two distinct classes of sharp-interface models, that correspond to different scalings in a small parameter εrepresenting the fast-reaction. Here εis the ratio between the characteristic scale of the diffusion of the fastest species and the one of the carbonation reaction. We explore three conceptually different scaling regimes (in terms of ε) of the effective diffusivities of the driving chemical species. The limiting models include one-phase and two-phase generalised Stefan moving-boundary problems as well as a nonstandard two-scale (micro-macro) moving-boundary problem -- the main result of the paper. Numerical results, supporting the asymptotics, illustrate the behavior of the concentration profiles for relevant parameter regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling micro-macro pedestrian counterflow in heterogeneous domains

We present a micro-macro strategy able to describe the dynamics of crowds in heterogeneous media. Herein we focus on the example of pedestrian counterflow. The main working tools include the use of mass and porosity measures together with their transport as well as suitable application of a version of Radon-Nikodym Theorem formulated for finite measures. Finally, we illustrate numerically our microscopic model and emphasize the effects produced by an implicitly defined social velocity. Keywords: Crowd dynamics; mass measures; porosity measure; social networks

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Sulfate attack in sewer pipes: Derivation of a concrete corrosion model via two-scale convergence

We explore the homogenization limit and rigorously derive upscaled equations for a microscopic reaction-diffusion system modeling sulfate corrosion in sewer pipes made of concrete. The system, defined in a periodically-perforated domain, is semi-linear, partially dissipative and weakly coupled via a non-linear ordinary differential equation posed on the solid-water interface at the pore level. Firstly, we show the well-posedness of the microscopic model. We then apply homogenization techniques based on two-scale convergence for an uniformly periodic domain and derive upscaled equations together with explicit formulae for the effective diffusion coefficients and reaction constants. We use a boundary unfolding method to pass to the homogenization limit in the non-linear ordinary differential equation. Finally, besides giving its strong formulation, we also prove that the upscaled two-scale model admits a unique solution.

0 citations0 reviewsNo code linkedOpen access