# American Institute of Mathematical Sciences

• Previous Article
Low Mach number limit of strong solutions for 3-D full compressible MHD equations with Dirichlet boundary condition
• DCDS-B Home
• This Issue
• Next Article
Emergence of large densities and simultaneous blow-up in a two-species chemotaxis system with competitive kinetics
October  2019, 24(10): 5481-5502. doi: 10.3934/dcdsb.2019067

## Simulation of a simple particle system interacting through hitting times

 University of Oxford, Andrew Wiles Building, Woodstock Road, Oxford, UK, OX2 6GG

The frst author gratefully acknowledges support from the Economic and Social Research Council and Bank of America Merrill Lynch.

Received  June 2018 Published  April 2019

We develop an Euler-type particle method for the simulation of a McKean–Vlasov equation arising from a mean-field model with positive feedback from hitting a boundary. Under assumptions on the parameters which ensure differentiable solutions, we establish convergence of order $1/2$ in the time step. Moreover, we give a modification of the scheme using Brownian bridges and local mesh refinement, which improves the order to $1$. We confirm our theoretical results with numerical tests and empirically investigate cases with blow-up.

Citation: Vadim Kaushansky, Christoph Reisinger. Simulation of a simple particle system interacting through hitting times. Discrete & Continuous Dynamical Systems - B, 2019, 24 (10) : 5481-5502. doi: 10.3934/dcdsb.2019067
##### References:

show all references

##### References:
(a) $L_t$ for different $\alpha$ near the jump; (b) distribution of $Y_T$ for $Y_T > 0$ before and after the jump. Fitted by kernel density estimation with normal kernel for $N = 10^7$
Convergence rate for $d_1(L_t, \tilde{L}_t)$, $d_2(L_t, \tilde{L}_t)$, $d_3(L_t, \tilde{L}_t)$ for (a) Algorithm 1, (b) Algorithm 2
Error of the loss process at $t=T$ for $\frac{1}{Y_0} \sim \exp(1)$ : %depending on time, (a) for increasing number $n$ of timesteps; (b) for increasing number $N$ of samples, %with $\frac{1}{Y_0} \sim \exp(1)$ both for Algorithms 1 and 2.
(a) $L_t$ and (b) $L'_t$ for different values of $\alpha$
Error of the loss process at $t = T$ for ${{Y}_{0}}\tilde{\ }\text{Gammadistr}(1+\beta ,1/2)$: (a) for increasing number $n$ of timesteps; (b) for increasing number $N$ of samples, both for Algorithms 1 and 2
(a) $L_t$ and (b) $L'_t$ for different values of $\alpha$
(a) Loss process computed using Algorithm 1 for different $n$; (b) error as a function of $t$
Convergence rate: at (a) $T = 0.001$, (b) $T = 0.008$
 [1] Seung-Yeal Ha, Jinwook Jung, Jeongho Kim, Jinyeong Park, Xiongtao Zhang. A mean-field limit of the particle swarmalator model. Kinetic & Related Models, , () : -. doi: 10.3934/krm.2021011 [2] Marco Cirant, Diogo A. Gomes, Edgard A. Pimentel, Héctor Sánchez-Morgado. On some singular mean-field games. Journal of Dynamics & Games, 2021  doi: 10.3934/jdg.2021006 [3] Bin Pei, Yong Xu, Yuzhen Bai. Convergence of p-th mean in an averaging principle for stochastic partial differential equations driven by fractional Brownian motion. Discrete & Continuous Dynamical Systems - B, 2020, 25 (3) : 1141-1158. doi: 10.3934/dcdsb.2019213 [4] Tomáš Roubíček. An energy-conserving time-discretisation scheme for poroelastic media with phase-field fracture emitting waves and heat. Discrete & Continuous Dynamical Systems - S, 2017, 10 (4) : 867-893. doi: 10.3934/dcdss.2017044 [5] Qiang Guo, Dong Liang. An adaptive wavelet method and its analysis for parabolic equations. Numerical Algebra, Control & Optimization, 2013, 3 (2) : 327-345. doi: 10.3934/naco.2013.3.327 [6] Deren Han, Zehui Jia, Yongzhong Song, David Z. W. Wang. An efficient projection method for nonlinear inverse problems with sparsity constraints. Inverse Problems & Imaging, 2016, 10 (3) : 689-709. doi: 10.3934/ipi.2016017 [7] Petra Csomós, Hermann Mena. Fourier-splitting method for solving hyperbolic LQR problems. Numerical Algebra, Control & Optimization, 2018, 8 (1) : 17-46. doi: 10.3934/naco.2018002 [8] Christina Surulescu, Nicolae Surulescu. Modeling and simulation of some cell dispersion problems by a nonparametric method. Mathematical Biosciences & Engineering, 2011, 8 (2) : 263-277. doi: 10.3934/mbe.2011.8.263 [9] Boris Kramer, John R. Singler. A POD projection method for large-scale algebraic Riccati equations. Numerical Algebra, Control & Optimization, 2016, 6 (4) : 413-435. doi: 10.3934/naco.2016018 [10] Jiangxing Wang. Convergence analysis of an accurate and efficient method for nonlinear Maxwell's equations. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2429-2440. doi: 10.3934/dcdsb.2020185 [11] Xiaomao Deng, Xiao-Chuan Cai, Jun Zou. A parallel space-time domain decomposition method for unsteady source inversion problems. Inverse Problems & Imaging, 2015, 9 (4) : 1069-1091. doi: 10.3934/ipi.2015.9.1069 [12] Lunji Song, Wenya Qi, Kaifang Liu, Qingxian Gu. A new over-penalized weak galerkin finite element method. Part Ⅱ: Elliptic interface problems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2581-2598. doi: 10.3934/dcdsb.2020196 [13] Kaifang Liu, Lunji Song, Shan Zhao. A new over-penalized weak galerkin method. Part Ⅰ: Second-order elliptic problems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2411-2428. doi: 10.3934/dcdsb.2020184 [14] Marion Darbas, Jérémy Heleine, Stephanie Lohrengel. Numerical resolution by the quasi-reversibility method of a data completion problem for Maxwell's equations. Inverse Problems & Imaging, 2020, 14 (6) : 1107-1133. doi: 10.3934/ipi.2020056 [15] Shangzhi Li, Shangjiang Guo. Permanence and extinction of a stochastic SIS epidemic model with three independent Brownian motions. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2693-2719. doi: 10.3934/dcdsb.2020201 [16] Caifang Wang, Tie Zhou. The order of convergence for Landweber Scheme with $\alpha,\beta$-rule. Inverse Problems & Imaging, 2012, 6 (1) : 133-146. doi: 10.3934/ipi.2012.6.133 [17] Michael Grinfeld, Amy Novick-Cohen. Some remarks on stability for a phase field model with memory. Discrete & Continuous Dynamical Systems - A, 2006, 15 (4) : 1089-1117. doi: 10.3934/dcds.2006.15.1089 [18] Abdulrazzaq T. Abed, Azzam S. Y. Aladool. Applying particle swarm optimization based on Padé approximant to solve ordinary differential equation. Numerical Algebra, Control & Optimization, 2021  doi: 10.3934/naco.2021008 [19] Wolf-Jüergen Beyn, Janosch Rieger. The implicit Euler scheme for one-sided Lipschitz differential inclusions. Discrete & Continuous Dynamical Systems - B, 2010, 14 (2) : 409-428. doi: 10.3934/dcdsb.2010.14.409 [20] Alina Chertock, Alexander Kurganov, Mária Lukáčová-Medvi${\rm{\check{d}}}$ová, Șeyma Nur Özcan. An asymptotic preserving scheme for kinetic chemotaxis models in two space dimensions. Kinetic & Related Models, 2019, 12 (1) : 195-216. doi: 10.3934/krm.2019009

2019 Impact Factor: 1.27