In this article, we propose a macro-micro (two-scale) mathematical model for describing the macroscopic swelling of a rubber foam caused by the microscopic absorption of some liquid. In our modeling approach, we suppose that the material occupies a one-dimensional domain which swells as described by the standard beam equation including an additional term determined by the liquid pressure. As special feature of our model, the absorption takes place inside the rubber foam via a lower length scale, which is assumed to be inherently present in such a structured material. The liquid's absorption and transport inside the material is modeled by means of a nonlinear parabolic equation derived from Darcy's law posed in a non-cylindrical domain defined by the macroscopic deformation (which is a solution of the beam equation). Under suitable assumptions, we establish the existence and uniqueness of a suitable class of solutions to our evolution system coupling the nonlinear parabolic equation posed on the microscopic non-cylindrical domain with the beam equation posed on the macroscopic cylindrical domain. In order to guarantee the regularity of the non-cylindrical domain, we impose a singularity to the elastic response function appearing in the beam equation.
In many industrial applications, rubber-based materials are routinely used in conjunction with various penetrants or diluents in gaseous or liquid form. It is of interest to estimate theoretically the penetration depth as well as the amount of diffusants stored inside the material. In this framework, we prove the global solvability and explore the large time-behavior of solutions to a one-phase free boundary problem with nonlinear kinetic condition that is able to describe the migration of diffusants into rubber. The key idea in the proof of the large time behavior is to benefit from a contradiction argument, since it is difficult to obtain uniform estimates for the growth rate of the free boundary due to the use of a Robin boundary condition posed at the fixed boundary.
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 t-behavior of reaction penetration depths by including nonlinear effects due to deviations from the classical Henry's law and time-dependent Dirichlet data.
We prove the large time behavior of solutions to a coupled thermo-diffusion arising in the modeling 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.
In our previous works we studied a one-dimensional free-boundary model related to the aggressive penetration of gaseous carbon dioxide in unsaturated concrete. Essentially, global existence and uniqueness of weak solutions to the model were obtained when the initial functions are bounded on the domain. In this paper we investigate the well-posedness of the problem for the case when the initial functions belong to a class. Specifically, the uniqueness of weak solutions is proved by applying the dual equation method.
We study the weak solvability of a system of coupled Allen–Cahn-like equations resembling cross-diffusion which arises as a model for the consolidation of saturated porous media. Besides using energy-like estimates, we cast the special structure of the system in the framework of the Leray–Schauder fixed-point principle and ensure in this way the local existence of strong solutions to a regularized version of our system. Furthermore, weak convergence techniques ensure the existence of weak solutions to the original consolidation problem. The uniqueness of global-in-time solutions is guaranteed in a particular case. Moreover, we use a finite difference scheme to show the negativity of the vector of solutions.
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. © 2013 Elsevier Inc.
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.
We present modeling strategies that describe the motion and interaction of groups of pedestrians in obscured spaces. We start off with an approach based on balance equations in terms of measures and then we exploit the descriptive power of a probabilistic cellular automaton model. Based on a variation of the simple symmetric random walk on the square lattice, we test the interplay between population size and an interpersonal attraction parameter for the evacuation of confined and darkened spaces. We argue that information overload and coordination costs associated with information processing in small groups are two key processes that influence the evacuation rate. Our results show that substantial computational resources are necessary to compensate for incomplete information - the more individuals in (information processing) groups the higher the exit rate for low population size. For simple social systems, it is likely that the individual representations are not redundant and large group sizes ensure that this non-redundant information is actually available to a substantial number of individuals. For complex social systems, information redundancy makes information evaluation and transfer inefficient and, as such, group size becomes a drawback rather than a benefit. The effect of group sizes on outgoing fluxes, evacuation times and wall effects is carefully studied with a Monte Carlo framework accounting also for the presence of an internal obstacle.
We study the upscaling of a system of many interacting particles through a heterogenous thin elongated obstacle as modeled via a two-dimensional diffusion problem with a one-directional nonlinear convective drift. Assuming that the obstacle can be described well by a thin composite strip with periodically placed microstructures, we aim at deriving the upscaled model equations as well as the effective transport coefficients for suitable scalings in terms of both the inherent thickness at the strip and the typical length scales of the microscopic heterogeneities. Aiming at computable scenarios, we consider that the heterogeneity of the strip is made of an array of periodically arranged impenetrable solid rectangles and identify two scaling regimes what concerns the small asymptotics parameter for the upscaling procedure: the characteristic size of the microstructure is either significantly smaller than the thickness of the thin obstacle or it is of the same order of magnitude. We scale up the diffusion-polynomial drift model and list computable formulas for the effective diffusion and drift tensorial coefficients for both scaling regimes. Our upscaling procedure combines ideas of two-scale asymptotics homogenization with dimension reduction arguments. Consequences of these results for the construction of more general transmission boundary conditions are discussed. We illustrate numerically the concentration profile of the chemical species passing through the upscaled strip in the finite thickness regime and point out that trapping of concentration inside the strip is likely to occur in at least two conceptually different transport situations: (i) full diffusion/dispersion matrix and nonlinear horizontal drift, and (ii) diagonal diffusion matrix and oblique nonlinear drift.
We study the pedestrian escape from an obscure room using a lattice gas model with twospecies of particles. One species, called passive, performs a symmetric random walk on the lattice,whereas the second species, called active, is subject to a drift guiding the particles towards the exit.The drift mimics the awareness of some pedestrians of the geometry of the room and of the location ofthe exit. We provide numerical evidence that, in spite of the hard core interaction between particles –namely, there can be at most one particle of any species per site – adding a fraction of active particlesin the system enhances the evacuation rate of all particles from the room. A similar effect is alsoobserved when looking at the outgoing particle flux, when the system is in contact with an externalparticle reservoir that induces the onset of a steady state. We interpret this phenomenon as a discretespace counterpart of the drafting effect typically observed in a continuum set–up as the aerodynamicdrag experienced by pelotons of competing cyclists.
Particle diffusion is modified by the presence of barriers. In cells macro-molecules, behaving as obstacles, slow down the dynamics so that the mean-square displacement of molecules grows with time as a power law with exponent smaller than one. In different situations, such as grain and pedestrian dynamics, it can happen that an obstacle can accelerate the dynamics. In the framework of very basic models, we study the time needed by particles to cross a strip for different bulk dynamics and discuss the effect of obstacles. We find that in some regimes such a residence time is not monotonic with respect to the size and the position of the obstacles. We can then conclude that, even in very elementary systems where no interaction among particles is considered, obstacles can either slow down or accelerate the particle dynamics depending on their geometry and position.
We consider the setup 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 us to control the particle fraction in the condensed phase. This phenomenon results in a current versus blockage parameter curve characterized by two nonanalyticity points.
Human crowds base most of their behavioral decisions upon anticipated states oftheir walking environment. We explore a minimal version of a lattice model to study lanesformation in pedestrian counterow. Using the concept of horizon depth, our simulationresults suggest that the anticipation eect together with the presence of a small backgroundnoise play an important role in promoting collective behaviors in a counterow setup. Theseingredients facilitate the formation of seemingly stable lanes and ensure the ergodicity of the system.
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.
Stimulated by experimental evidence in the field of solution-born thin films, we study the morphology formation in a three state lattice system subjected to the evaporation of one component. The practical problem that we address is the understanding of the parameters that govern morphology formation from a ternary mixture upon evaporation, as is the case in the fabrication of thin films from solution for organic photovoltaics. We use, as a tool, a generalized version of the Potts and Blume-Capel models in 2D, with the Monte Carlo Kawasaki-Metropolis algorithm, to simulate the phase behaviour of a ternary mixture upon evaporation of one of its components. The components with spin 1, −1 and 0 in the Blume-Capel dynamics correspond to the electron-acceptor, electron-donor and solvent molecules, respectively, in a ternary mixture used in the preparation of the active layer films in an organic solar cell. Furthermore, we introduce parameters that account for the relative composition of the mixture, temperature, and interaction between the species in the system. We identify the parameter regions that are prone to facilitate the phase separation. Furthermore, we study qualitatively the types of formed configurations. We show that even a relatively simple model, as the present one, can generate key morphological features, similar to those observed in experiments, which proves the method valuable for the study of complex systems.
We study the motion of pedestrians through an obscure tunnel where the lack of visibility hides the exits. Using a lattice model, we explore the effects of communication on the effective transport properties of the crowd of pedestrians. More precisely, we study the effect of two thresholds on the structure of the effective nonlinear diffusion coefficient. One threshold models pedestrian communication efficiency in the dark, while the other one describes the tunnel capacity. Essentially, we note that if the evacuees show a maximum trust (leading to a fast communication), they tend to quickly find the exit and hence the collective action tends to prevent the occurrence of disasters.
Fundamental diagrams describing the relation between pedestrians' speed and density are key points in understanding pedestrian dynamics. Experimental data evidence the onset of complex behaviors in which the velocity decreases with the density, and different logistic regimes are identified. This paper addresses the issue of pedestrian transport and of fundamental diagrams for a scenario involving the motion of pedestrians escaping from an obscure tunnel. We capture the effects of communication efficiency and exit capacity by means of two thresholds controlling the rate at which particles (walkers, pedestrians) move on the lattice. Using a particle system model, we show that in the absence of limitation in communication among pedestrians, we reproduce with good accuracy the standard fundamental diagrams, whose basic behaviors can be interpreted in terms of exit capacity limitation. When the effect of limited communication ability is considered, then interesting nonintuitive phenomena occur. In particular, we shed light on the loss of monotonicity of the typical speed-density curves, revealing the existence of a pedestrian density optimizing the escape. We study both the discrete particle dynamics and the corresponding hydrodynamic limit (a porous medium equation and a transport (continuity) equation). We also point out the dependence of the effective transport coefficients on the two thresholds---the essence of the microstructure information.Read More: http://epubs.siam.org/doi/10.1137/15M1030960
Npoint 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 thresholdT. When the geometrical parameters of the billiard table are fixed, this mechanism gives rise to non-equilibrium phase transitions in the largeNlimit: lettingT/Ndecrease, 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.
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.
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. The research reported here opens many fundamental questions and should be seen therefore as a preliminary investigation of the very complex behavior of the people and their motion in dark regions.
We investigate the appearance of trapping states in pedestrian flows through bottlenecks as a result of the interplay between the geometry of the system and the microscopic stochastic dynamics. We model the flow through a bottleneck via a Zero Range Process on a one-dimensional periodic lattice. Particle are removed from the lattice sites with rates proportional to the local occupation numbers. The bottleneck is modeled by a particular site of the lattice whose updating rate saturates to a constant value as soon as the local occupation number exceeds a fixed threshold. We show that for any finite value of the threshold the stationary particle current saturates to the limiting bottleneck rate when the total particle density in the system exceeds a critical value corresponding to the bottleneck rate itself.
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.
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.
A procedure for model reduction of stochastic ordinary differential equations with additive noise was recently introduced in Colangeli et al (2022 J. Phys. A: Math. Theor.55 505002), based on the Invariant Manifold method and on the Fluctuation–Dissipation relation. A general question thus arises as to whether one can rigorously quantify the error entailed by the use of the reduced dynamics in place of the original one. In this work we provide explicit formulae and estimates of the error in terms of the Wasserstein distance, both in the presence or in the absence of a sharp time-scale separation between the variables to be retained or eliminated from the description, as well as in the long-time behavior.
Understanding and modeling the dynamics of pedestrian crowds can help with designing and increasing the safety of civil facilities. A key feature of a crowd is its intrinsic stochasticity, appearing even under very diluted conditions, due to the variability in individual behaviors. Individual stochasticity becomes even more important under densely crowded conditions, since it can be nonlinearly magnified and may lead to potentially dangerous collective behaviors. To understand quantitatively crowd stochasticity, we study the real-life dynamics of a large ensemble of pedestrians walking undisturbed, and we perform a statistical analysis of the fully resolved pedestrian trajectories obtained by a yearlong high-resolution measurement campaign. Our measurements have been carried out in a corridor of the Eindhoven University of Technology via a combination of Microsoft Kinect 3D range sensor and automatic head-tracking algorithms. The temporal homogeneity of our large database of trajectories allows us to robustly define and separate average walking behaviors from fluctuations parallel and orthogonal with respect to the average walking path. Fluctuations include rare events when individuals suddenly change their minds and invert their walking directions. Such tendency to invert direction has been poorly studied so far, even if it may have important implications on the functioning and safety of facilities. We propose a model for the dynamics of undisturbed pedestrians, based on stochastic differential equations, that provides a good agreement with our field observations, including the occurrence of rare events.
Focusing on a specific crowd dynamics situation, including real life experiments and measurements, our paper targets a twofold aim: (1) we present a Bayesian probabilistic method to estimate the value and the uncertainty (in the form of a probability density function) of parameters in crowd dynamic models from the experimental data; and (2) we introduce a fitness measure for the models to classify a couple of model structures (forces) according to their fitness to the experimental data, preparing the stage for a more general model-selection and validation strategy inspired by probabilistic data analysis. Finally, we review the essential aspects of our experimental setup and measurement technique.
We report the results of a simulation study in which we explore the joint effect of group absorptive capacity (as the average individual rationality of the group members) and cognitive distance (as the distance between the most rational group member and the rest of the group) on the emergence of collective rationality in groups. We start from empirical results reported in the literature on group rationality as collective group level competence and use data on real-life groups of four and five to validate a mathematical model. We then use this mathematical model to predict group level scores from a variety of possible group configurations (varying both in cognitive distance and average individual rationality). Our results show that both group competence and cognitive distance are necessary conditions for emergent group rationality. Group configurations, in which the groups become more rational than the most rational group member, are groups scoring low on cognitive distance and scoring high on absorptive capacity.
We discuss the existence of a class of weak solutions to a nonlinear parabolic system of reaction-diffusion type endowed with singular production terms by reaction. The singularity is due to a potential occurrence of quenching localized to the domain boundary. The kind of quenching we have in mind is due to a twofold contribution: (i) the choice of boundary conditions, modeling in our case the contact with an infinite reservoir filled with ready-to-react chemicals and (ii) the use of a particular nonlinear, non-Lipschitz structure of the reaction kinetics. Our working techniques use fine energy estimates for approximating non-singular problems and uniform control on the set where singularities are localizing.
We are interested in exploring interacting particle systemsthat can be seen as microscopic models for a particular structure ofcoupled transport flux arising when different populations are jointlyevolving. The scenarios we have in mind are inspired by the dynamicsof pedestrian flows in open spaces and are intimately connectedto cross-diffusion and thermo-diffusion problems holding a variationalstructure. The tools we use include a suitable structure of the relativeentropy controlling TV-norms, the construction of Lyapunov functionalsand particular closed-form solutions to nonlinear transport equations,a hydrodynamics limiting procedure due to Philipowski, as wellas the construction of numerical approximates to both the continuumlimit problem in 2D and to the original interacting particle systems.
We investigate corrector estimates for the solutions of a thermoelasticity problem posed in a highly heterogeneous two-phase medium and its corresponding two-scale thermoelasticity model which was derived in [11] by two-scale convergence arguments. The medium in question consists of a connected matrix with disconnected, initially periodically distributed inclusions separated by a sharp interface undergoing a priori known phase transformations. While such estimates seem not to be obtainable in the fully coupled setting, we show that for some simplified scenarios optimal convergence rates can be proven rigorously. The main technique for the proofs are energy estimates using special reconstructions of two-scale functions and particular operator estimates for periodic functions with zero average. Here, additional regularity results for the involved functions are necessary.
We investigate a linear, fully coupled thermoelasticity problem for a highly heterogeneous, two-phase medium. The medium in question consists of a connected matrix with disconnected, initially periodically distributed inclusions separated by a sharp interface undergoing an a prioriknown interface movement because of phase transformations. After transforming the moving geometry to an ϵ-periodic, fixed reference domain, we establish the well-posedness of the model and derive a number of ϵ-independent a priori estimates. Via a two-scale convergence argument, we then show that the ϵ-dependent solutions converge to solutions of a corresponding upscaled model with distributed time-dependent microstructures.
We study the weak solvability of a macroscopic, quasilinear reaction–diffusion system posed in a 2D porous medium which undergoes microstructural problems. The solid matrix of this porous medium is assumed to be made out of circles of not-necessarily uniform radius. The growth or shrinkage of these circles, which are governed by an ODE, has direct feedback to the macroscopic diffusivity via an additional elliptic cell problem. The reaction–diffusion system describes the macroscopic diffusion, aggregation, and deposition of populations of colloidal particles of various sizes inside a porous media made of prescribed arrangement of balls. The mathematical analysis of this two-scale problem relies on a suitable application of Schauder's fixed point theorem which also provides a convergent algorithm for an iteration method to compute finite difference approximations of smooth solutions to our multiscale model. Numerical simulations illustrate the behavior of the local concentration of the colloidal populations close to clogging situations.