A WKB-based fixed-grid method for capturing trait concentration in a dispersal evolution model
It provides an efficient numerical tool for evolutionary ecologists studying trait evolution under rare mutations, where traditional methods fail due to extreme concentration.
The paper develops a WKB-based numerical method for capturing trait concentration in a dispersal evolution model, achieving accurate resolution of the selected trait and concentration structure on fixed grids, with significant computational advantages over direct discretizations in the small-mutation regime.
The evolution of dispersal traits is a fundamental topic in evolutionary ecology, where natural selection may drive the trait distribution toward concentration in the rare-mutation regime. This singular behavior poses a serious numerical difficulty, since direct discretizations of the population density require very fine trait grids to identify the fittest trait and to resolve the concentrated profile accurately. In this paper, we develop a WKB-based numerical framework for a dispersal evolution model. By separating the exponentially concentrated trait dependence from a smoother amplitude variable and combining this WKB representation with dual trait-grid implementation and other specially designed techniques, the method recovers the selected trait and the associated concentration structure accurately and efficiently on fixed trait grids. We establish a semi-discrete stationary fixed-grid asymptotic-preserving structure for the rare-mutation limit of the steady-state problem. Numerical experiments compare the proposed method with direct density discretizations and confirm its advantage in the small-mutation regime.