Influx, outflux, and no-flux boundary conditions (respectively, on the left, right, and remaining boundaries). $x\in(0,1)$, $y\in(0,1.5)$, $t\in(0,5)$. $\varepsilon = 0.2$. $\rho_{\textrm{In}} = 0.5$.
Second-order scheme. $M_x=512$, $M_y=256$, $\Delta x= \Delta y=2^{-9}$, $\Delta t=\Delta x$.
Pedro Aceves-Sánchez · Rafael Bailo · Pierre Degond · Zoé Mercier
Mathematical Models and Methods in Applied Sciences, 34 (6), 2024.