In this paper, we apply the local discontinuous Galerkin (LDG) method to 2D Keller--Segel (KS) chemotaxis model. We improve the results upon (Y. Epshteyn and A. Kurganov, SIAM Journal on Numerical Analysis, 47 (2008), 368-408) and give optimal rate of convergence under special finite element spaces before the blow-up occurs (the exact solutions are smooth). Moreover, to construct physically relevant numerical approximations, we consider $P^1$ LDG scheme and develop a positivity-preserving limiter to the scheme, extending the idea in (Y. Zhang, X. Zhang and C.-W. Shu, Journal of Computational Physics, 229 (2010), 8918-8934). With this limiter, we can prove the $L^1$-stability of the numerical scheme. Numerical experiments are performed to demonstrate the good performance of the positivity-preserving LDG scheme. Moreover, it is known that the chemotaxis model will yield blow-up solutions under certain initial conditions. We numerically demonstrate how to find the approximate blow-up time by using the $L^2$-norm of the $L^1$-stable numerical solution.

In this paper, we develop a second-order asymptotic-preserving and positivity-preserving discontinuous Galerkin (DG) scheme for the Kerr-Debye model. By using the approach first introduced by Zhang and Shu in Ref.~\refcite{ZhSh:04} with an energy estimate and Taylor expansion, the asymptotic-preserving property of the semi-discrete DG methods is proved rigorously. In addition, we propose a class of unconditional positivity-preserving implicit-explicit (IMEX) Runge-Kutta methods for the system of ordinary differential equations arising from the semi-discretization of the Kerr-Debye model. The new IMEX Runge-Kutta methods are based on the modification of the strong-stability-preserving (SSP) implicit Runge-Kutta method and have second-order accuracy. The numerical results validate our analysis.

In this paper, we analyze the stability of the fourth order Runge-Kutta method for integrating semi-discrete approximations of time-dependent partial differential equations. Our study focuses on linear problems and covers general semi-bounded spatial discretizations. A counter example is given to show that the classical four-stage fourth order Runge-Kutta method can not preserve the one-step strong stability, even though the ordinary differential equation system is energy-decaying. But with an energy argument, we show that the strong stability property holds in two steps under an appropriate time step constraint. Based on this fact, the stability extends to general well-posed linear systems. As an application, we utilize the results to examine the stability of the fourth order Runge-Kutta approximations of several specific method of lines schemes for hyperbolic problems, including the spectral Galerkin method and the discontinuous Galerkin method.

To a torus action on a complex vector space, Gelfand, Kapranov and Zelevinsky introduce a system of differential equations, called the GKZ hypergeometric system. Its solutions are GKZ hypergeometric functions. We study the $p$-adic counterpart of the GKZ hypergeometric system. In the language of dagger spaces introduced by Grosse-Kl\"onne, the $p$-adic GKZ hypergeometric complex is a twisted relative de Rham complex of meromorphic differential forms with logarithmic poles for an affinoid toric dagger space over the dagger unit polydisc. It is a complex of ${\mathcal O}^\dagger$-modules with integrable connections and with Frobenius structures defined on the dagger unit polydisc such that traces of Frobenius on fibers at Techm\"uller points define the hypergeometric function over the finite field introduced by Gelfand and Graev.

We develop a general theory of canonical bases for quantum symmetric pairs (U, U^i) with parameters of arbitrary finite type. We construct new canonical bases for the simple integrable U-modules and their tensor products regarded as U^i-modules. We also construct a canonical basis for the modified form of the i-quantum group U^i. To that end, we establish several new structural results on quantum symmetric pairs, such as bilinear forms, braid group actions, integral forms, Levi subalgebras (of real rank one), and integrality of the intertwiners.