We propose a semi-discrete scheme for 2D Keller-Segel equations based on a symmetrization reformation, which is equivalent to the convex splitting method and is free of any nonlinear solver. We show that, this new scheme is unconditionally stable as long as the initial condition does not exceed certain threshold, and it asymptotically preserves the quasi-static limit in the transient regime. Furthermore, we prove that the fully discrete scheme is conservative and positivity preserving, which makes it ideal for simulations. The analogical schemes for the radial symmetric cases and the subcritical degenerate cases are also presented and analyzed. With extensive numerical tests, we verify the claimed properties of the methods and demonstrate their superiority in various challenging applications.

We consider an inverse time-dependent source problem governed by a distributed time-fractional diffusion equation using interior measurement data. Such a problem arises in some ultra-slowdiffusion phenomena in many applied areas. Based on the regularity result of the solution to the direct problem, we establish the solvability of this inverse problem as well as the conditional stability in suitable function space with a weak norm. By a variational identity connecting the unknown time-dependent source and the interior measurement data, the conjugate gradient method is also introduced to construct the inversion algorithm under the framework of regularizing scheme. We show the validity of the proposed scheme by several numerical examples.

No abstract uploaded!

We solve the Plateau problem for marginally outer trapped surfaces in general Cauchy data sets. We employ the Perron method and tools from geometric measure theory to force and control a blow-up of Jang’s equation. Substantial new geometric insights regarding the lower order properties of marginally outer trapped surfaces are gained in the process. The techniques developed in this paper are flexible and can be adapted to other non-variational existence problems.

We prove the Poisson geometric version of the Local Reeb Stability (from foliation theory) and of the Slice Theorem (from equivariant geometry), which is also a generalization of Conn’s linearization theorem.