nthispaperwewillpresentthelocalstabilityanalysisandlocalerrorestimatefor the local discontinuous Galerkin (LDG) method, when solving the time-dependent singularly perturbed problems in one dimensional space with a stationary outflow boundary layer. Based on a general framework on the local stability, we obtain the optimal error estimate out of the local subdomain, which is nearby the outflow boundary point and has the width of O(h log(1/h)), for the semi-discrete LDG scheme and the fully-discrete LDG scheme with the second order explicit Runge–Kutta time-marching. Here h is the maximum mesh length. The numerical experiments are given also.

We present a domain decomposition method on a semilinear hyperbolic system with multiple relaxation times. In the region where the relaxation time is small, an asymptotic equilibrium equation can be used for computational efficiency. An interface condition based on the sign of the characteristic speed at the interface is provided to couple the two systems in a domain decomposition setting. A rigorous analysis, based on the Laplace Transform, on the $ L^2$ error estimate is presented for the linear case, which shows how the error of the domain decomposition method depends on the smaller relaxation time, and the boundary and interface layer effects. The given convergence rate is optimal. We present a numerical implementation of this domain decomposition method, and give some numerical results in order to study the performance of this method.

We prove the existence of the global weak entropy solution to the Doi–Saintillan–Shelley model for active and passive rod-like particle suspensions, which couples a Fokker–Planck equation with the incompressible Navier–Stokes or Stokes equation, under the no-flux boundary conditions, L2(Ω;L1(Sd611))L2(Ω;L1(Sd611)) initial data, and finite initial entropy for the particle distribution function in two and three dimensions. Furthermore, for the model with the Stokes equation, we obtain the global L2(Ω×Sd611)L2(Ω×Sd611) weak solution in two and three dimensions and the uniqueness in two dimension.

We investigate systems of self-propelled particles with alignment interaction. Compared to previous work (Degond and Motsch, Math. Models Methods Appl. Sci. 18:1193–1215, 2008a; Frouvelle, Math. Models Methods Appl. Sci., 2012), the force acting on the particles is not normalized, and this modification gives rise to phase transitions from disordered states at low density to aligned states at high densities. This model is the space-inhomogeneous extension of (Frouvelle and Liu, Dynamics in a kinetic model of oriented particles with phase transition, 2012), in which the existence and stability of the equilibrium states were investigated. When the density is lower than a threshold value, the dynamics is described by a nonlinear diffusion equation. By contrast, when the density is larger than this threshold value, the dynamics is described by a similar hydrodynamic model for self-alignment interactions as derived in (Degond and Motsch, Math. Models Methods Appl. Sci. 18:1193–1215, 2008a; Frouvelle, Math. Models Methods Appl. Sci., 2012). However, the modified normalization of the force gives rise to different convection speeds, and the resulting model may lose its hyperbolicity in some regions of the state space.

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.