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.

We studied coupled systems of the Fokker--Planck equation and the Navier--Stokes equation modeling the Hookean and the finitely extensible nonlinear elastic (FENE)-type polymeric flows. We proved the continuous embedding and compact embedding theorems in weighted spaces that naturally arise from related entropy estimates. These embedding estimates are shown to be sharp. For the Hookean polymeric system with a center-of-mass diffusion and a superlinear spring potential, we proved the existence of a global weak solution. Moreover, we were able to tackle the FENE model with $L^2$ initial data for the polymer density instead of the $L^\infty$ counterpart in the literature.