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.

No abstract uploaded!

In this paper, we introduce a nonlocal model for linear steady Stokes system with physical no-slip boundary condition. We use the idea of volume constraint to enforce the no-slip boundary condition. We prove that the nonlocal model is well-posed and the solution of the nonlocal system converges to the solution of the original Stokes system when the nonlocality vanishes.

We establish the global existence and uniqueness of classical solutions to the Cauchy problem for the isentropic compressible Navier-Stokes equations in three spatial dimensions with smooth initial data that are of small energy but possibly large oscillations with constant state as far field, which could be either vacuum or nonvacuum. The initial density is allowed to vanish, and the spatial measure of the set of vacuum can be arbitrarily large; in particular, the initial density can even have compact support. These results generalize previous results on classical solutions for initial densities being strictly away from vacuum and are the first for global classical solutions that may have large oscillations and can contain vacuum states. © 2012 Wiley Periodicals, Inc

In this paper we are concerned with the existence of a weak solution to the initial boundary value problem for the equation ∂ t u=Δ(Δu) −3 . This problem arises in the mathematical modeling of the evolution of a crystal surface. Existence of a weak solution u with Δu≥0 is obtained via a suitable substitution. Our investigations reveal the close connection between this problem and the equation ∂ t ρ+ρ 2 Δ 2 ρ 3 =0 , another crystal surface model first proposed by H. Al Hajj Shehadeh, R. V. Kohn and J. Weare in [1].