# American Institute of Mathematical Sciences

September  2017, 10(3): 669-688. doi: 10.3934/krm.2017027

## Finite range method of approximation for balance laws in measure spaces

 Institute of Applied Mathematics and Mechanics, University of Warsaw, ul. Banacha 2, 02-097 Warsaw, Poland

* Corresponding author: Piotr Gwiazda

Received  April 2016 Revised  July 2016 Published  December 2016

In the following paper we reconsider a numerical scheme recently introduced in [10]. The method was designed for a wide class of size structured population models with a nonlocal term describing the birth process. Despite its numerous advantages it features the exponential growth in time of the number of particles constituting the numerical solution. We introduce a new algorithm free from this inconvenience. The improvement is based on the application the Finite Range Approximation to the nonlocal term. We prove the convergence of the derived method and provide the rate of its convergence. Moreover, the results are illustrated by numerical simulations applied to various test cases.

Citation: Piotr Gwiazda, Piotr Orlinski, Agnieszka Ulikowska. Finite range method of approximation for balance laws in measure spaces. Kinetic & Related Models, 2017, 10 (3) : 669-688. doi: 10.3934/krm.2017027
##### References:

show all references

##### References:
Function $f$ and its Finite Range Approximation
Function $f^{\varepsilon}$ and its approximation $\bar f^{\varepsilon}$
Order of convergence

Results for $\varepsilon \in \{0.1, 0.0125, 0.0015625, \Delta t\}$ presented in Tables 1-4

The error of the FRA method for $\varepsilon=0.1$
 $\Delta t$ ${\rm{Err}}(1,\Delta t,0.1)$ $\log({\frac{{\rm{Err}}(T,2\Delta t,\varepsilon)}{{\rm{Err}}(T,\Delta t,\varepsilon)}}) / \log 2$ $1.0000\cdot 10^{-1}$ $0.0008865973$ - $5.0000\cdot 10^{-2}$ $0.0005886514$ 0.59086543 $2.5000\cdot 10^{-2}$ $0.0003434678$ 0.77723882 $1.2500\cdot 10^{-2}$ $0.0002668756$ 0.36400752 $6.2500\cdot 10^{-3}$ $0.0002400045$ 0.15310595 $3.1250\cdot 10^{-3}$ $0.0002258416$ 0.08775007 $1.5625\cdot 10^{-3}$ $0.0002187803$ 0.04582869 $7.8125\cdot 10^{-4}$ $0.0002154824$ 0.02191237
 $\Delta t$ ${\rm{Err}}(1,\Delta t,0.1)$ $\log({\frac{{\rm{Err}}(T,2\Delta t,\varepsilon)}{{\rm{Err}}(T,\Delta t,\varepsilon)}}) / \log 2$ $1.0000\cdot 10^{-1}$ $0.0008865973$ - $5.0000\cdot 10^{-2}$ $0.0005886514$ 0.59086543 $2.5000\cdot 10^{-2}$ $0.0003434678$ 0.77723882 $1.2500\cdot 10^{-2}$ $0.0002668756$ 0.36400752 $6.2500\cdot 10^{-3}$ $0.0002400045$ 0.15310595 $3.1250\cdot 10^{-3}$ $0.0002258416$ 0.08775007 $1.5625\cdot 10^{-3}$ $0.0002187803$ 0.04582869 $7.8125\cdot 10^{-4}$ $0.0002154824$ 0.02191237
The error of the FRA method for $\varepsilon=0.0125$
 $\Delta t$ ${\rm{Err}}(1,\Delta t,0.0125)$ $\log({\frac{{\rm{Err}}(T,2\Delta t,\varepsilon)}{{\rm{Err}}(T,\Delta t,\varepsilon)}})/ \log 2$ $1.0000\cdot 10^{-1}$ $6.702793\cdot 10^{-4}$ - $5.0000\cdot 10^{-2}$ $4.255311\cdot 10^{-4}$ 0.6554979 $2.5000\cdot 10^{-2}$ $1.757573\cdot 10^{-4}$ 1.2756799 $1.2500\cdot 10^{-2}$ $9.330613\cdot 10^{-5}$ 0.9135408 $6.2500\cdot 10^{-3}$ $5.658711\cdot 10^{-5}$ 0.7214983 $3.1250\cdot 10^{-3}$ $4.055238\cdot 10^{-5}$ 0.4806869 $1.5625\cdot 10^{-3}$ $3.332588\cdot 10^{-5}$ 0.2831438 $7.8125\cdot 10^{-4}$ $2.973577\cdot 10^{-5}$ 0.1644434
 $\Delta t$ ${\rm{Err}}(1,\Delta t,0.0125)$ $\log({\frac{{\rm{Err}}(T,2\Delta t,\varepsilon)}{{\rm{Err}}(T,\Delta t,\varepsilon)}})/ \log 2$ $1.0000\cdot 10^{-1}$ $6.702793\cdot 10^{-4}$ - $5.0000\cdot 10^{-2}$ $4.255311\cdot 10^{-4}$ 0.6554979 $2.5000\cdot 10^{-2}$ $1.757573\cdot 10^{-4}$ 1.2756799 $1.2500\cdot 10^{-2}$ $9.330613\cdot 10^{-5}$ 0.9135408 $6.2500\cdot 10^{-3}$ $5.658711\cdot 10^{-5}$ 0.7214983 $3.1250\cdot 10^{-3}$ $4.055238\cdot 10^{-5}$ 0.4806869 $1.5625\cdot 10^{-3}$ $3.332588\cdot 10^{-5}$ 0.2831438 $7.8125\cdot 10^{-4}$ $2.973577\cdot 10^{-5}$ 0.1644434
The error of the FRA method for $\varepsilon=0.0015625$
 $\Delta t$ ${\rm{Err}}(1,\Delta t,0.0015625)$ $\log({\frac{{\rm{Err}}(T,2\Delta t,\varepsilon)}{{\rm{Err}}(T,\Delta t,\varepsilon)}}) / \log 2$ $1.0000\cdot 10^{-1}$ $6.616854\cdot 10^{-4}$ - $5.0000\cdot 10^{-2}$ $4.139056\cdot 10^{-4}$ 0.6768439 $2.5000\cdot 10^{-2}$ $1.570582\cdot 10^{-4}$ 1.3980025 $1.2500\cdot 10^{-2}$ $7.418803\cdot 10^{-5}$ 1.0820407 $6.2500\cdot 10^{-3}$ $3.649949\cdot 10^{-5}$ 1.0820407 $3.1250\cdot 10^{-3}$ $1.883518\cdot 10^{-5}$ 0.9544464 $1.5625\cdot 10^{-3}$ $1.032763\cdot 10^{-5}$ 0.8669200 $7.8125\cdot 10^{-4}$ $6.277910\cdot 10^{-6}$ 0.7181537
 $\Delta t$ ${\rm{Err}}(1,\Delta t,0.0015625)$ $\log({\frac{{\rm{Err}}(T,2\Delta t,\varepsilon)}{{\rm{Err}}(T,\Delta t,\varepsilon)}}) / \log 2$ $1.0000\cdot 10^{-1}$ $6.616854\cdot 10^{-4}$ - $5.0000\cdot 10^{-2}$ $4.139056\cdot 10^{-4}$ 0.6768439 $2.5000\cdot 10^{-2}$ $1.570582\cdot 10^{-4}$ 1.3980025 $1.2500\cdot 10^{-2}$ $7.418803\cdot 10^{-5}$ 1.0820407 $6.2500\cdot 10^{-3}$ $3.649949\cdot 10^{-5}$ 1.0820407 $3.1250\cdot 10^{-3}$ $1.883518\cdot 10^{-5}$ 0.9544464 $1.5625\cdot 10^{-3}$ $1.032763\cdot 10^{-5}$ 0.8669200 $7.8125\cdot 10^{-4}$ $6.277910\cdot 10^{-6}$ 0.7181537
The error of the FRA method for $\varepsilon=\Delta t$
 $\Delta t$ ${\rm{Err}}(1,\Delta t,\Delta t)$ $\log({\frac{{\rm{Err}}(T,2\Delta t,2\varepsilon)}{{\rm{Err}}(T,\Delta t,\varepsilon)}}) / \log 2$ $1.0000\cdot 10^{-1}$ $8.865973\cdot 10^{-4}$ - $5.0000\cdot 10^{-2}$ $4.797352\cdot 10^{-4}$ 0.8860407 $2.5000\cdot 10^{-2}$ $1.991736\cdot 10^{-4}$ 1.2682122 $1.2500\cdot 10^{-2}$ $9.330613\cdot 10^{-5}$ 1.0939823 $6.2500\cdot 10^{-3}$ $4.477800\cdot 10^{-5}$ 1.0591816 $3.1250\cdot 10^{-3}$ $2.162411\cdot 10^{-5}$ 1.0501496 $1.5625\cdot 10^{-3}$ $1.032763\cdot 10^{-5}$ 1.0661308 $7.8125\cdot 10^{-4}$ $4.763338\cdot 10^{-6}$ 1.1164650
 $\Delta t$ ${\rm{Err}}(1,\Delta t,\Delta t)$ $\log({\frac{{\rm{Err}}(T,2\Delta t,2\varepsilon)}{{\rm{Err}}(T,\Delta t,\varepsilon)}}) / \log 2$ $1.0000\cdot 10^{-1}$ $8.865973\cdot 10^{-4}$ - $5.0000\cdot 10^{-2}$ $4.797352\cdot 10^{-4}$ 0.8860407 $2.5000\cdot 10^{-2}$ $1.991736\cdot 10^{-4}$ 1.2682122 $1.2500\cdot 10^{-2}$ $9.330613\cdot 10^{-5}$ 1.0939823 $6.2500\cdot 10^{-3}$ $4.477800\cdot 10^{-5}$ 1.0591816 $3.1250\cdot 10^{-3}$ $2.162411\cdot 10^{-5}$ 1.0501496 $1.5625\cdot 10^{-3}$ $1.032763\cdot 10^{-5}$ 1.0661308 $7.8125\cdot 10^{-4}$ $4.763338\cdot 10^{-6}$ 1.1164650
 [1] Giulia Luise, Giuseppe Savaré. Contraction and regularizing properties of heat flows in metric measure spaces. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 273-297. doi: 10.3934/dcdss.2020327 [2] Anna Abbatiello, Eduard Feireisl, Antoní Novotný. Generalized solutions to models of compressible viscous fluids. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 1-28. doi: 10.3934/dcds.2020345 [3] Laurent Di Menza, Virginie Joanne-Fabre. An age group model for the study of a population of trees. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020464 [4] Ebraheem O. Alzahrani, Muhammad Altaf Khan. Androgen driven evolutionary population dynamics in prostate cancer growth. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020426 [5] Fabio Camilli, Giulia Cavagnari, Raul De Maio, Benedetto Piccoli. Superposition principle and schemes for measure differential equations. Kinetic & Related Models, , () : -. doi: 10.3934/krm.2020050 [6] Shasha Hu, Yihong Xu, Yuhan Zhang. Second-Order characterizations for set-valued equilibrium problems with variable ordering structures. Journal of Industrial & Management Optimization, 2020  doi: 10.3934/jimo.2020164 [7] Giuseppina Guatteri, Federica Masiero. Stochastic maximum principle for problems with delay with dependence on the past through general measures. Mathematical Control & Related Fields, 2020  doi: 10.3934/mcrf.2020048 [8] Fanni M. Sélley. A self-consistent dynamical system with multiple absolutely continuous invariant measures. Journal of Computational Dynamics, 2021, 8 (1) : 9-32. doi: 10.3934/jcd.2021002 [9] Haodong Yu, Jie Sun. Robust stochastic optimization with convex risk measures: A discretized subgradient scheme. Journal of Industrial & Management Optimization, 2021, 17 (1) : 81-99. doi: 10.3934/jimo.2019100 [10] Xin Guo, Lexin Li, Qiang Wu. Modeling interactive components by coordinate kernel polynomial models. Mathematical Foundations of Computing, 2020, 3 (4) : 263-277. doi: 10.3934/mfc.2020010 [11] Wei Ouyang, Li Li. Hölder strong metric subregularity and its applications to convergence analysis of inexact Newton methods. Journal of Industrial & Management Optimization, 2021, 17 (1) : 169-184. doi: 10.3934/jimo.2019105 [12] Harrison Bray. Ergodicity of Bowen–Margulis measure for the Benoist 3-manifolds. Journal of Modern Dynamics, 2020, 16: 305-329. doi: 10.3934/jmd.2020011 [13] Mark F. Demers. Uniqueness and exponential mixing for the measure of maximal entropy for piecewise hyperbolic maps. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 217-256. doi: 10.3934/dcds.2020217 [14] Wei Feng, Michael Freeze, Xin Lu. On competition models under allee effect: Asymptotic behavior and traveling waves. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5609-5626. doi: 10.3934/cpaa.2020256 [15] Mehdi Bastani, Davod Khojasteh Salkuyeh. On the GSOR iteration method for image restoration. Numerical Algebra, Control & Optimization, 2021, 11 (1) : 27-43. doi: 10.3934/naco.2020013 [16] Russell Ricks. The unique measure of maximal entropy for a compact rank one locally CAT(0) space. Discrete & Continuous Dynamical Systems - A, 2021, 41 (2) : 507-523. doi: 10.3934/dcds.2020266 [17] Yongge Tian, Pengyang Xie. Simultaneous optimal predictions under two seemingly unrelated linear random-effects models. Journal of Industrial & Management Optimization, 2020  doi: 10.3934/jimo.2020168 [18] Annegret Glitzky, Matthias Liero, Grigor Nika. Dimension reduction of thermistor models for large-area organic light-emitting diodes. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020460 [19] Yen-Luan Chen, Chin-Chih Chang, Zhe George Zhang, Xiaofeng Chen. Optimal preventive "maintenance-first or -last" policies with generalized imperfect maintenance models. Journal of Industrial & Management Optimization, 2021, 17 (1) : 501-516. doi: 10.3934/jimo.2020149 [20] Hong Niu, Zhijiang Feng, Qijin Xiao, Yajun Zhang. A PID control method based on optimal control strategy. Numerical Algebra, Control & Optimization, 2021, 11 (1) : 117-126. doi: 10.3934/naco.2020019

2019 Impact Factor: 1.311