We establish two nonlinear compactness theorems in Lp(0,T;B) with hypothesis on time translations, which are nonlinear counterparts of two results by Simon (1987) [1]. The first theorem sharpens a result by Maitre (2003) [10] and is important in the study of doubly nonlinear elliptic–parabolic equations. Based on this theorem, we then obtain a time translation counterpart of a result by Dubinskii˘ (1965) [5], which is supposed to be useful in the study of some nonlinear kinetic equations (e.g. the FENE-type bead–spring chains model).

In this paper we study the Euler-Poincar\'{e} equations in $\Bbb R^N$. We prove local existence of weak solutions in $W^{2,p}(\Bbb R^N),$ $p>N$, and local existence of unique classical solutions in $H^k (\Bbb R^N)$, $k>N/2+3$, as well as a blow-up criterion. For the zero dispersion equation($\alpha=0$) we prove a finite time blow-up of the classical solution. We also prove that as the dispersion parameter vanishes, the weak solution converges to a solution of the zero dispersion equation with sharp rate as $\alpha\to0$, provided that the limiting solution belongs to $C([0, T);H^k(\Bbb R^N))$ with $k>N/2 +3$. For the {\em stationary weak solutions} of the Euler-Poincar\'{e} equations we prove a Liouville type theorem. Namely, for $\alpha>0$ any weak solution $\mathbf{u}\in H^1(\Bbb R^N)$ is $\mathbf{u}=0$; for $\alpha=0$ any weak solution $\mathbf{u}\in L^2(\Bbb R^N)$ is $\mathbf{u}=0$.

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 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.