In this paper, we introduce a random particle blob method for the Keller-Segel equation (with dimension ) and establish a rigorous convergence analysis.

A key difficulty in the analysis and numerical approximation of the shallow water equations is the nonconservative product of measures due to the gravitational force acting on a sloped bottom. Solutions may be nonunique, and numerical schemes are not only consistent discretizations of the shallow water equations, but they also determine how to model the physics. Our derivation is based on a subcell reconstruction using infinitesimal singular layers at the cell boundaries, as inspired by S. Noelle, Y. Xing, and C.-W. Shu [J. Comput. Phys., 226 (2007), pp. 29-58]. One key step is to separate the singular measures. Another aspect is the reconstruction of the solution variables in the singular layers. We study three reconstructions. The first leads to the well-known scheme of Audusse et al., [SIAM J. Sci. Comput., 25 (2004), pp. 2050-2065], which introduces the hydrostatic reconstruction. The second is a modification proposed in [T. Morales de Luna, M. J. Castro D´ıaz, and C. Par´es, Appl. Math. Comput., 219 (2013), pp. 9012-9032], which analyzes whether a wave has enough energy to overcome a step. The third is our new scheme, which borrows its structure from the wet-dry front. For a number of cases discussed in recent years, where water runs down a hill, Audusse’s scheme converges slowly or fails. Morales’ scheme gives a visible improvement. Both schemes are clearly outperformed by our new scheme.

We studied coupled systems of the Fokker--Planck equation and the Navier--Stokes equation modeling the Hookean and the finitely extensible nonlinear elastic (FENE)-type polymeric flows. We proved the continuous embedding and compact embedding theorems in weighted spaces that naturally arise from related entropy estimates. These embedding estimates are shown to be sharp. For the Hookean polymeric system with a center-of-mass diffusion and a superlinear spring potential, we proved the existence of a global weak solution. Moreover, we were able to tackle the FENE model with $L^2$ initial data for the polymer density instead of the $L^\infty$ counterpart in the literature.

No abstract uploaded!

This paper deals with a degenerate diffusion Patlak–Keller–Segel system in $n\geq 3$ dimension. The main difference between the current work and many other recent studies on the same model is that we study the diffusion exponent $m=2n/(n+2)$, which is smaller than the usual exponent $m^{*}=2-2/n$ used in other studies. With the exponent $m=2n/(n+2)$, the associated free energy is conformal invariant, and there is a family of stationary solutions $U_{\lambda,x_0}(x)=C(\frac{\lambda} {\lambda^2+|x-x_0|^2})^{\frac{n+2}{2}}$ $\forall \lambda>0$, $x_0\in {\mathbb R}^n$. For radially symmetric solutions, we prove that if the initial data are strictly below $U_{\lambda,0}(x)$ for some $\lambda$, then the solution vanishes in $L^1_{loc}$ as $t\to\infty$; if the initial data are strictly above $U_{\lambda,0}(x)$ for some $\lambda$, then the solution either blows up at a finite time or has a mass concentration at $r=0$ as time goes to infinity. For general initial data, we prove that there is a global weak solution provided that the $L^m$ norm of initial density is less than a universal constant, and the weak solution vanishes as time goes to infinity. We also prove a finite time blow-up of the solution if the $L^m$ norm for initial data is larger than the $L^m$ norm of $U_{\lambda,x_0}(x)$, which is constant independent of $\lambda$ and $x_0$, and the free energy of initial data is smaller than that of $U_{\lambda,x_0}(x)$.