Long time simulations of transport equations raise computational challenges since they require both a large domain of calculation and sufficient accuracy. It is therefore advantageous, in terms of computational costs, to use a time varying adaptive mesh, with small cells in the region of interest and coarser cells where the solution is smooth. Biological models involving cell dynamics fall for instance within this framework and are often non conservative to account for cell division. In that case the threshold controlling the spatial adaptivity may have to be time-dependent in order to keep up with the progression of the solution. In this article we tackle the difficulties arising when applying a Multiresolution method to a transport equation with discontinuous fluxes modeling localized mitosis. The analysis of the numerical method is performed on a simplified model and numerical scheme. An original threshold strategy is proposed and validated thanks to extensive numerical tests. It is then applied to a biological model in both cases of distributed and localized mitosis.
Mots-clés : adaptive finite volumes, non conservative PDE, discontinuous flux
@article{M2AN_2014__48_5_1381_0, author = {Aymard, Benjamin and Cl\'ement, Fr\'ed\'erique and Postel, Marie}, title = {Adaptive mesh refinement strategy for a non conservative transport problem}, journal = {ESAIM: Mathematical Modelling and Numerical Analysis }, pages = {1381--1412}, publisher = {EDP-Sciences}, volume = {48}, number = {5}, year = {2014}, doi = {10.1051/m2an/2014014}, mrnumber = {3264358}, language = {en}, url = {http://www.numdam.org/articles/10.1051/m2an/2014014/} }
TY - JOUR AU - Aymard, Benjamin AU - Clément, Frédérique AU - Postel, Marie TI - Adaptive mesh refinement strategy for a non conservative transport problem JO - ESAIM: Mathematical Modelling and Numerical Analysis PY - 2014 SP - 1381 EP - 1412 VL - 48 IS - 5 PB - EDP-Sciences UR - http://www.numdam.org/articles/10.1051/m2an/2014014/ DO - 10.1051/m2an/2014014 LA - en ID - M2AN_2014__48_5_1381_0 ER -
%0 Journal Article %A Aymard, Benjamin %A Clément, Frédérique %A Postel, Marie %T Adaptive mesh refinement strategy for a non conservative transport problem %J ESAIM: Mathematical Modelling and Numerical Analysis %D 2014 %P 1381-1412 %V 48 %N 5 %I EDP-Sciences %U http://www.numdam.org/articles/10.1051/m2an/2014014/ %R 10.1051/m2an/2014014 %G en %F M2AN_2014__48_5_1381_0
Aymard, Benjamin; Clément, Frédérique; Postel, Marie. Adaptive mesh refinement strategy for a non conservative transport problem. ESAIM: Mathematical Modelling and Numerical Analysis , Tome 48 (2014) no. 5, pp. 1381-1412. doi : 10.1051/m2an/2014014. http://www.numdam.org/articles/10.1051/m2an/2014014/
[1] Numerical simulation of the selection process of the ovarian follicles. ESAIM: Proc. 38 (2012) 99-117. | MR
, , and ,[2] A numerical method for transport equations with discontinuous flux functions: Application to mathematical modeling of cell dynamics. SIAM J. Sci. Comput. 35 (2013) 2442-2468. | MR | Zbl
, , and ,[3] Fully adaptive multiresolution schemes for strongly degenerate parabolic equations with discontinuous flux. J. Engrg. Math. 60 (2008) 365-385. | MR | Zbl
, , and ,[4] Multiscale modelling of ovarian follicular selection. Prog. Biophys. Mol. Biol. 113 (2013) 398-408.
and ,[5] Fully adaptive multiresolution finite volume schemes for conservation laws. Math. Comput. 72 (2003) 183-225. | MR | Zbl
, , and ,[6] Adaptive multiresolution for finite volume solutions of gas dynamics. Comput. Fluids 32 (2003) 31-38. | Zbl
, and ,[7] Multiresolution and Adaptive Methods for Convection-Dominated Problems, edited by F. Coquel, Y. Maday, S. Müller, M. Postel and Q. Tran, vol. 29. ESAIM: Proc. (2009) 1-108. | MR
[8] Local time stepping applied to implicit-explicit methods for hyperbolic systems. Multiscale Model. Simul. 8 (2010) 540-570. | MR | Zbl
, , and ,[9] Convergence of time-space adaptive algorithms for nonlinear conservation laws. IMA J. Numer. Anal. 32 (2012) 1440-1483. | MR | Zbl
, and ,[10] Multiscale modeling of follicular ovulation as a reachability problem. Multiscale Model. Simul. 6 (2007) 895-912. | MR | Zbl
, and ,[11] Multi-scale modeling of the follicle selection process in the ovary. Math. Biosci. 198 (2005) 57-79. | MR | Zbl
, , and ,[12] Multiresolution algorithms for the numerical solutions of hyperbolic conservation laws. Commun. Pure Appl. Math. 48 (1995) 1305-1342. | MR | Zbl
,[13] Summer school on multiresolution and adaptive mesh refinement methods, edited by V. Louvet and M. Massot, vol. 34. ESAIM: Proc. (2011) 1-290. | MR
[14] Multiscale modeling of follicular ovulation as a mass and maturity dynamical system. Multiscale Model. Simul. 9 (2011) 282-313. | MR | Zbl
,[15] Adaptive multiscale schemes for conservation laws, vol. 27 of Lect. Notes Comput. Sci. Eng. Springer-Verlag, Berlin (2003). | MR | Zbl
,[16] Transport Equations In Biology. Front. Math. Birkhauser Verlag, Basel (2007). | MR | Zbl
,[17] Visualizing spatiotemporal dynamics of multicellular cell-cycle progression. Cell 132 (2008) 487-498.
, , , , , , , , , , , , and ,[18] Cauchy problem for multiscale conservation laws: Applications to structured cell populations. J. Math. Anal. Appl. 401 (2013) 896-920. | MR
,[19] Contrasting quiescent G0 phase with mitotic cell cycling in the mouse immune system. PLoS ONE 8 (2013) e73801.
, , , , , , and ,Cité par Sources :