# American Institute of Mathematical Sciences

December  2018, 11(6): 1475-1501. doi: 10.3934/krm.2018058

## Linear Boltzmann dynamics in a strip with large reflective obstacles: Stationary state and residence time

 Dipartimento di Scienze di Base e Applicate per l'Ingegneria, Sapienza Università di Roma, via A. Scarpa 16, Ⅰ - 00161, Roma, Italy

* Corresponding author: Alessandro Ciallella

Received  August 2017 Published  June 2018

The presence of obstacles modifies the way in which particles diffuse. In cells, for instance, it is observed that, due to the presence of macromolecules playing the role of obstacles, the mean-square displacement of biomolecules scales as a power law with exponent smaller than one. On the other hand, different situations in grain and pedestrian dynamics in which the presence of an obstacle accelerates the dynamics are known. We focus on the time, called the residence time, needed by particles to cross a strip assuming that the dynamics inside the strip follows the linear Boltzmann dynamics. We find that the residence time is not monotonic with respect to the size and the location of the obstacles, since the obstacle can force those particles that eventually cross the strip to spend a smaller time in the strip itself. We focus on the case of a rectangular strip with two open sides and two reflective sides and we consider reflective obstacles into the strip. We prove that the stationary state of the linear Boltzmann dynamics, in the diffusive regime, converges to the solution of the Laplace equation with Dirichlet boundary conditions on the open sides and homogeneous Neumann boundary conditions on the other sides and on the obstacle boundaries.

Citation: Alessandro Ciallella, Emilio N. M. Cirillo. Linear Boltzmann dynamics in a strip with large reflective obstacles: Stationary state and residence time. Kinetic & Related Models, 2018, 11 (6) : 1475-1501. doi: 10.3934/krm.2018058
##### References:

show all references

##### References:
Domain $\Omega$: strip with large fixed obstacles, where $\partial\Omega_L$ and $\partial\Omega_R$ are the vertical open boundaries and $\partial\Omega_E$ are reflective boundaries.
Elastic collision with a scatterers: impact parameter $\delta$ and angle of incidence $\alpha$.
Plot of the simulated solutions $h_{t_m}$ in a $3D$ plot and in a $2D$ plot constructed by averaging on the $x_2$ variable: in dark gray $t_m = 2\cdot 10^{-1}$, in gray $t_m = 10^{-1}$, in light gray $t_m = 2\cdot10^{-2}$. In black (grid and dashed line) the analytic solution $\rho$ of the associated Laplace problem.
Simulation parameter $t_m = 10^{-2}$: relative error ${|h_{t_m}-\rho|}/{\rho}$.
Simulation parameter $t_m = 10^{-2}$: on the left in gray the numerical solution $h_{t_m}$ and in black the solution $\rho$ of the associated Laplace problem; on the right the relative error ${|h_{t_m}-\rho|}/{\rho}$. Into the strip there is a square obstacle with side $8\cdot10^{-1}$.
Simulation parameter $t_m = 10^{-2}$: on the left in gray the numerical solution $h_{t_m}$ and in black the solution $\rho$ of the associated Laplace problem; on the right the relative error ${|h_{t_m}-\rho|}/{\rho}$. In the strip is placed a very thin obstacle with height of $0.8$.
Simulation parameter $t_m = 10^{-2}$: on the left in gray the numerical solution $h_{t_m}$ and in black the solution $\rho$ of the associated Laplace problem; on the right the relative error ${|h_{t_m}-\rho|}/{\rho}$. In the first line we show the case of two squared obstacles with side $6\cdot 10^{-1}$, in the second one a couple of rectangular obstacles, taller and thinner than the squares.
Residence time vs. height of a centered rectangular obstacle with fixed width $4\cdot 10^{-2}$ (on the left) and $4 \cdot 10^{-1}$ (on the right). Simulation parameters: $L_1 = 4$, $L_2 = 1$, $t_m = 2\cdot 10^{-2}$, total number of inserted particles $10^8$, the total number of particles exiting through the right boundary varies from $5.3\cdot 10^{5}$ to $3.6\cdot 10^{5}$ (on the left) and from $5.3\cdot 10^{5}$ to $2.1\cdot 10^{5}$ (on the right) depending on the obstacle height. The solid lines represent the value of the residence time measured for the empty strip (no obstacle).
Residence time vs. height of a centered rectangular obstacle with fixed width $8\cdot 10^{-1}$ (on the left) and $12 \cdot 10^{-1}$ (on the right). Simulation parameters: $L_1 = 4$, $L_2 = 1$, $t_m = 2\cdot 10^{-2}$, total number of inserted particles $10^8$, the total number of particles exiting through the right boundary varies from $5.2\cdot 10^{5}$ to $1.4\cdot 10^{5}$ (on the left) and from $5.2 \cdot 10^{5}$ to $1.1 \cdot 10^{5}$ (on the right) depending on the obstacle height. The solid lines represent the value of the residence time measured for the empty strip (no obstacle).
Residence time vs. width of a centered rectangular obstacle with fixed height $0.8$ (on the left) and vs. the side length of a centered squared obstacle (on the right). Simulation parameters: $L_1 = 4$, $L_2 = 1$, $t_m = 2\cdot 10^{-2}$, total number of inserted particles $10^{8}$, the total number of particles exiting through the right boundary varies from $4.2 \cdot 10^{5}$ to $1.1\cdot 10^{5}$ (on the left) and from $5.3\cdot 10^5$ to $1.3\cdot 10^5$ (on the right) depending on the obstacle width. The solid lines represent the value of the residence time measured for the empty strip (no obstacle).
Residence time vs. position of the center of the obstacle. The obstacle is a square of side length $0.8$ on the left and a rectangle of side lengths $0.04$ and $0.8$ on the right. Simulation parameters: $L_1 = 4$, $L_2 = 1$, $t_m = 2\cdot 10^{-2}$, total number of inserted particles $10^8$, the total number of particles exiting through the right boundary is stable at the order of $2.6\cdot 10^5$ (on the left) and of $4\cdot 10^5$ (on the right) not depending on the obstacle position. The solid lines represent the value of the residence time measured for the empty strip (no obstacle).
As in the right panel in Figure 8. In the left panel the height of the obstacle is equal to $0.8$. Left panel: the mean time spent by particles crossing the strip in each point of the strip ($0.02\times0.02$ cells have been considered) for the empty strip case (black) and in presence of the obstacle (gray). Right panel: residence time in regions L (circles), C (squares), and R (triangles) in presence of the obstacle (gray) and for the empty strip case (black).
As in Figure 12 for the geometry in the left panel in Figure 10. In the left panel the width of the obstacle is $2.28$.
As in Figure 12 for the geometry in the left panel in Figure 11. In the left panel the position of the center of the obstacle is $0.8$.
Domain $\Lambda$: infinite strip with big fixed obstacles: the whole boundaries of $\Lambda$ is a specular reflective boundary.
 [1] Giacomo Dimarco. The moment guided Monte Carlo method for the Boltzmann equation. Kinetic & Related Models, 2013, 6 (2) : 291-315. doi: 10.3934/krm.2013.6.291 [2] Guillaume Bal, Ian Langmore, Youssef Marzouk. Bayesian inverse problems with Monte Carlo forward models. Inverse Problems & Imaging, 2013, 7 (1) : 81-105. doi: 10.3934/ipi.2013.7.81 [3] Daewa Kim, Annalisa Quaini. A kinetic theory approach to model pedestrian dynamics in bounded domains with obstacles. Kinetic & Related Models, 2019, 12 (6) : 1273-1296. doi: 10.3934/krm.2019049 [4] Claude Bardos, François Golse, Ivan Moyano. Linear Boltzmann equation and fractional diffusion. Kinetic & Related Models, 2018, 11 (4) : 1011-1036. doi: 10.3934/krm.2018039 [5] Jiakou Wang, Margaret J. Slattery, Meghan Henty Hoskins, Shile Liang, Cheng Dong, Qiang Du. Monte carlo simulation of heterotypic cell aggregation in nonlinear shear flow. Mathematical Biosciences & Engineering, 2006, 3 (4) : 683-696. doi: 10.3934/mbe.2006.3.683 [6] Michael B. Giles, Kristian Debrabant, Andreas Rössler. Analysis of multilevel Monte Carlo path simulation using the Milstein discretisation. Discrete & Continuous Dynamical Systems - B, 2019, 24 (8) : 3881-3903. doi: 10.3934/dcdsb.2018335 [7] Joseph Nebus. The Dirichlet quotient of point vortex interactions on the surface of the sphere examined by Monte Carlo experiments. Discrete & Continuous Dynamical Systems - B, 2005, 5 (1) : 125-136. doi: 10.3934/dcdsb.2005.5.125 [8] Chjan C. Lim, Joseph Nebus, Syed M. Assad. Monte-Carlo and polyhedron-based simulations I: extremal states of the logarithmic N-body problem on a sphere. Discrete & Continuous Dynamical Systems - B, 2003, 3 (3) : 313-342. doi: 10.3934/dcdsb.2003.3.313 [9] Olli-Pekka Tossavainen, Daniel B. Work. Markov Chain Monte Carlo based inverse modeling of traffic flows using GPS data. Networks & Heterogeneous Media, 2013, 8 (3) : 803-824. doi: 10.3934/nhm.2013.8.803 [10] Mazyar Zahedi-Seresht, Gholam-Reza Jahanshahloo, Josef Jablonsky, Sedighe Asghariniya. A new Monte Carlo based procedure for complete ranking efficient units in DEA models. Numerical Algebra, Control & Optimization, 2017, 7 (4) : 403-416. doi: 10.3934/naco.2017025 [11] Pedro Aceves-Sánchez, Christian Schmeiser. Fractional diffusion limit of a linear kinetic equation in a bounded domain. Kinetic & Related Models, 2017, 10 (3) : 541-551. doi: 10.3934/krm.2017021 [12] Tomasz Komorowski. Long time asymptotics of a degenerate linear kinetic transport equation. Kinetic & Related Models, 2014, 7 (1) : 79-108. doi: 10.3934/krm.2014.7.79 [13] Michael Herty, Giuseppe Visconti. Kinetic methods for inverse problems. Kinetic & Related Models, 2019, 12 (5) : 1109-1130. doi: 10.3934/krm.2019042 [14] Hélène Hivert. Numerical schemes for kinetic equation with diffusion limit and anomalous time scale. Kinetic & Related Models, 2018, 11 (2) : 409-439. doi: 10.3934/krm.2018019 [15] Robert M. Strain. Coordinates in the relativistic Boltzmann theory. Kinetic & Related Models, 2011, 4 (1) : 345-359. doi: 10.3934/krm.2011.4.345 [16] Luisa Arlotti, Bertrand Lods, Mustapha Mokhtar-Kharroubi. Non-autonomous Honesty theory in abstract state spaces with applications to linear kinetic equations. Communications on Pure & Applied Analysis, 2014, 13 (2) : 729-771. doi: 10.3934/cpaa.2014.13.729 [17] Qingguang Guan, Max Gunzburger. Stability and convergence of time-stepping methods for a nonlocal model for diffusion. Discrete & Continuous Dynamical Systems - B, 2015, 20 (5) : 1315-1335. doi: 10.3934/dcdsb.2015.20.1315 [18] Jean Dolbeault, Giuseppe Toscani. Fast diffusion equations: Matching large time asymptotics by relative entropy methods. Kinetic & Related Models, 2011, 4 (3) : 701-716. doi: 10.3934/krm.2011.4.701 [19] Darryl D. Holm, Vakhtang Putkaradze, Cesare Tronci. Collisionless kinetic theory of rolling molecules. Kinetic & Related Models, 2013, 6 (2) : 429-458. doi: 10.3934/krm.2013.6.429 [20] Emmanuel Frénod, Mathieu Lutz. On the Geometrical Gyro-Kinetic theory. Kinetic & Related Models, 2014, 7 (4) : 621-659. doi: 10.3934/krm.2014.7.621

2019 Impact Factor: 1.311