APPROXIMATION OF HYPERBOLIC MODELS FOR CHEMOSENSITIVE MOVEMENT.pdfVIP

SIAM J. SCI. COMPUT. c 2005 Society for Industrial and Applied Mathematics Vol. 27, No. 3, pp. 850–872 APPROXIMATION OF HYPERBOLIC MODELS FOR CHEMOSENSITIVE MOVEMENT∗ FRANCIS FILBET† AND CHI-WANG SHU‡ Abstract. Numerical methods with different orders of accuracy are proposed to approximate hyperbolic models for chemosensitive movements. On the one hand, first- and second-order well- balanced finite volume schemes are presented. This approach provides exact conservation of the steady state solutions. On the other hand, a high-order finite difference weighted essentially nonoscil- latory (WENO) scheme is constructed and the well-balanced reconstruction is adapted to this scheme to exactly preserve steady states and to retain high-order accuracy. Numerical simulations are per- formed to verify accuracy and the well-balanced property of the proposed schemes and to observe the formation of networks in the hyperbolic models similar to those observed in the experiments. Key words. chemotaxis, hyperbolic systems, finite volume methods, finite difference methods, weighted essentially nonoscillatory schemes, well-balanced schemes AMS subject classifications. 35L60, 92C15, 65M06 DOI. 10.1137/040604054 1. Introduction. Chemosensitive movement is a process by which cells change their direction when reacting to the presence of a chemical substance, approach- ing chemically favorable environments, and avoiding unfavorable environments. In the biological literature there is a distinction between chemotaxis and chemokine- sis. Chemokinesis describes nondirected bias in the movement behavior, which indi- rectly leads to an oriented movement