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