We propose a simple finite difference scheme for Navier--Stokes equations in primitive formulation on curvilinear domains. With proper boundary treatment and interplay between covariant and contravariant components, the spatial discretization admits exact Hodge decomposition and energy identity. As a result, the pressure can be decoupled from the momentum equation with explicit time stepping. No artificial pressure boundary condition is needed. In addition, it can be shown that this spatially compatible discretization leads to uniform inf-sup condition, which plays a crucial role in the pressure approximation of both dynamic and steady state calculations. Numerical experiments demonstrate the robustness and efficiency of our scheme.

This paper is concerned with the derivation and analysis of hydrodynamic models for systems of self-propelled particles subject to alignment interaction and attraction-repulsion. The starting point is the kinetic model considered in earlier work of Degond & Motsch with the addition of an attraction-repulsion interaction potential. Introducing different scalings than in Degond & Motsch, the non-local effects of the alignment and attraction-repulsion interactions can be kept in the hydrodynamic limit and result in extra pressure, viscosity terms and capillary force. The systems are shown to be symmetrizable hyperbolic systems with viscosity terms. A local-in-time existence result is proved in the 2D case for the viscous model and in the 3D case for the inviscid model. The proof relies on the energy method.

This paper investigates infinite-time spreading and finite-time blow-up for the Keller-Segel system. For 0 < m ≤ 2 − 2 / d, the L p space for both dynamic and steady solutions are detected with ${p:=\frac{d(2-m)}{2} }$ . Firstly, the global existence of the weak solution is proved for small initial data in L p . Moreover, when m > 1 − 2 / d, the weak solution preserves mass and satisfies the hyper-contractive estimates in L q for any p < q < ∞. Furthermore, for slow diffusion 1 < m ≤ 2 − 2/d, this weak solution is also a weak entropy solution which blows up at finite time provided by the initial negative free energy. For m > 2 − 2/d, the hyper-contractive estimates are also obtained. Finally, we focus on the L p norm of the steady solutions, it is shown that the energy critical exponent m = 2d/(d + 2) is the critical exponent separating finite L p norm and infinite L p norm for the steady state solutions.

Abstract We prove the convergence of the vortex blob method for a family of nonlinear evolutionary partial differential equations (PDEs), the so-called b-equation. This kind of PDEs, including the Camassa-Holm equation and the Degasperis-Procesi equation, has many applications in diverse scientific fields. Our convergence analysis also provides a proof for the existence of the global weak solution to the b-equation when the initial data is a nonnegative Radon measure with compact support.

This paper establishes the hyper-contractivity in L∞ (ℝd) (it's known as ultra-contractivity) for the multi-dimensional Keller-Segel systems with the dif;fusion exponent m > 1 -2/d. The results show that for the supercritical and critical case 1 - 2/d < m ≤ 2 - 2/d, if || U0||d(2-m)/2) < Cd, m where Cd, m is a universal constant, then for any t > 0, || u (⋅, t)|| ... is bounded and decays as t goes to infinity, For the subcritical case m > 2 - 2/d, the solution u (⋅, t) ∈ L∞ (ℝd) with any initial data U0 ∈ L1+ (ℝd) for any positive time.