MathSciDoc: An Archive for Mathematician ∫

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.

AbstractWe introduce a sub-cell shock capturing method for scalar conservation laws built upon the Jin-Xin relaxation framework. Here, sub-cell shock capturing is achieved thanks to an original defect measure correction technique. The proposed correction exactly restores entropy shock solutions of the exact Riemann problem and moreover, it produces monotone and entropy satisfying approximate self-similar solutions. These solutions are then sampled using Glimm's random choice method to advance in time. The resulting scheme combines the simplicity of the Jin-Xin relaxation method with the resolution of the Glimm's scheme in the sharp capture of discontinuities. The strong benefit of using defect measure corrections over usual sub-cell shock capturing methods is that the scheme can be easily made entropy satisfying with respect to infinitely many entropy pairs. Consequently, under a classical CFL condition, the method is proved to converge to the unique entropy weak solution of the Cauchy problem for general non-linear flux functions.

This work considers the rigorous derivation of continuum models of step motion starting from a mesoscopic Burton-Cabrera-Frank (BCF) type model following the work [Xiang, SIAM J. Appl. Math. 2002]. We prove that as the lattice parameter goes to zero, for a finite time interval, a modified discrete model converges to the strong solution of the limiting PDE with first order convergence rate.

ABSTRACT We study coagulation-fragmentation equations inspired by a simple model proposed in fisheries science to explain data for the size distribution of schools of pelagic fish. Although the equations lack detailed balance and admit no $H$-theorem, we are able to develop a rather complete description of equilibrium profiles and large-time behavior, based on recent developments in complex function theory for Bernstein and Pick functions. In the large-population continuum limit, a scaling-invariant regime is reached in which all equilibria are determined by a single scaling profile. This universal profile exhibits power-law behavior crossing over from exponent $-\frac23$ for small size to $-\frac32$ for large size, with an exponential cut-off.

We establish an optimal error estimate for a random particle blob method for the Keller-Segel equation in Rd (d  2). With a blob size " = h

In this paper, we study existence of global entropy weak solutions to a critical-case unstable thin lm equation in one-dimensional case ht + @x(h n @xxxh) + @x(h n+2@xh) = 0; where n  1. There exists a critical mass Mc = 2 p 6 3 found by Witelski etal. (2004 Euro. J. of Appl. Math. 15, 223{256) for n = 1. We obtain global existence of a non-negative entropy weak solution if initial mass is less than Mc. For n  4, entropy weak solutions are positive and unique. For n = 1,a nite time blow-up occurs for solutions with initial mass larger than Mc.For the Cauchy problem with n = 1 and initial mass less than Mc, we show that at least one of the following long-time behavior holds: the second moment goes to innity as the time goes to innity or h(
; tk) * 0 in L1 (R) for some subsequence tk ! 1.

We investigate the long-time dynamics of an opinion formation model inspired by a work by Borghesi, Bouchaud and Jensen. First, we derive a Fokker–Planck-type equation under the assumption that interactions between individuals produce little consensus of opinion (grazing collision approximation). Second, we study conditions under which the Fokker–Planck equation has non-trivial equilibria and derive the macroscopic limit (corresponding to the long-time dynamics and spatially localized interactions) for the evolution of the mean opinion. Finally, we compare two different types of interaction rates: the original one given in the work of Borghesi, Bouchaud and Jensen (symmetric binary interactions) and one inspired from works by Motsch and Tadmor (non-symmetric binary interactions). We show that the first case leads to a conservative model for the density of the mean opinion whereas the second case leads to a non-conservative equation. We also show that the speed at which consensus is reached asymptotically for these two rates has fairly different density dependence.

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

This paper investigates the existence of a uniform in time $L^{\infty}$ bounded weak entropy solution for the quasilinear parabolic-parabolic Keller-Segel model with the supercritical diffusion exponent $0 < m < 2-\frac{2}{d}$ in the multi-dimensional space ${\mathbb{R}}^d$ under the condition that the $L^{\frac{d(2-m)}{2}}$ norm of initial data is smaller than a universal constant. Moreover, the weak entropy solution $u(x,t)$ satisfies mass conservation when $m > 1-\frac{2}{d}$. We also prove the local existence of weak entropy solutions and a blow-up criterion for general $L^1\cap L^{\infty}$ initial data.

Abstract In this paper, we provide an alternative proof for the classical Sz. Nagy inequality in one dimension by a variational method and generalize it to higher dimensions J(h):=(∫Rd|h| dx)a611∫Rd|63h|2 dx(∫Rd|h|m+1 dx)a+1m+1≥β0, where m > 0 for d = 1, 2, for , and . The Euler–Lagrange equation for critical points of in the non-negative radial decreasing function space is given by a free boundary problem for a generalized Lane–Emden equation, which has a unique solution (denoted by h c ) and the solution determines the best constant for the above generalized Sz. Nagy inequality. The connection between the critical mass for the thin-film equation and the best constant of the Sz. Nagy inequality in one dimension was first noted by Witelski et al (2004 Eur. J. Appl. Math. 15 223–56). For the following critical thin film equation in multi-dimension ht+6366(h 63Δh)+6366(h 63hm)=0,x∈Rd, where m = 1 + 2/d, the critical mass is also given by . A finite time blow-up occurs for solutions with the initial mass larger than M c . On the other hand, if the initial mass is less than M c and a global non-negative entropy weak solution exists, then the second moment goes to infinity as or in for some subsequence . This shows that a part of the mass spreads to infinity.