In this paper, we apply a simple finite element numerical scheme, proposed in an earlier work (Liu in Math Comput 70(234):579-593, 2000), to perform a high resolution numerical simulation of incompressible flow over an irregular domain and analyze its boundary layer separation. Compared with many classical finite element fluid solvers, this numerical method avoids a Stokes solver, and only two Poisson-like equations need to be solved at each time step/stage. In addition, its combination with the fully explicit fourth order Runge-Kutta (RK4) time discretization enables us to compute high Reynolds number flow in a very efficient way. As an application of this robust numerical solver, the dynamical mechanism of the boundary layer separation for a triangular cavity flow with Reynolds numbers \(Re=10^4\) and \(Re=10^5\) , including the precise values of bifurcation location and critical time, are reported in this paper. In addition, we provide a super-convergence analysis for the simple finite element numerical scheme, using linear elements over a uniform triangulation with right triangles.

We consider a system coupling the incompressible Navier–Stokes equations to the Vlasov–Fokker–Planck equation. Such a problem arises in the description of particulate flows. We design a numerical scheme to simulate the behavior of the system. This scheme is asymptotic-preserving, thus efficient in both the kinetic and hydrodynamic regimes. It has a numerical stability condition controlled by the non-stiff convection operator, with an implicit treatment of the stiff drag term and the Fokker–Planck operator. Yet, consistent to a standard asymptotic-preserving Fokker–Planck solver or an incompressible Navier–Stokes solver, only the conjugate–gradient method and fast Poisson and Helmholtz solvers are needed. Numerical experiments are presented to demonstrate the accuracy and asymptotic behavior of the scheme, with several interesting applications.

We consider a multidimensional version of an inequality due to Leray as a substitute for Hardy’s inequality in the case $p = n ≥ 2$. In this paper we provide an optimal Sobolev-type improvement of this substitute, analogous to the corresponding improvements obtained for $p = 2 < n$ in S. Filippas, A. Tertikas, Optimizing improved Hardy inequalities, J. Funct. Anal. 192 (1) (2002) 186–233, and for $p > n ≥ 1$ in G. Psaradakis, An optimal Hardy-Morrey inequality, Calc. Var. Partial Differential Equations 45 (3-4) (2012) 421–441.

We study the structures of Fourier coefficients of automorphic forms on symplectic groups based on their local and global structures related to Arthur parameters. This is a first step towards the general conjecture on the relation between the structure of Fourier coefficients and Arthur parameters for automorphic forms occurring in the discrete spectrum, given by the first named author.

The authors proved in [5] that in a complete, connected, and simply connected Riemannian manifold, the volume of the intersection of two small geodesic balls depends only on the distance between the centers and the radii if and only if the space is harmonic. In this paper, we show that this proposition remains true, if the condition is restricted to balls of equal radii.