As a typical example of hyperbolic conservation laws, the system of Euler equations representing the conservation of mass, momentum and energy has three basic wave patterns. They are two nonlinear waves called shock and rarefaction wave, and one linearly degenerate wave called contact discontinuity. These basic wave patterns have their counterparts both in other physical models for gas motions in equilibrium involving viscosity and heat conductivity, and gas motion in non-equilibrium. The stability of these basic wave patterns in the system of Navier-Stokes equations and the Boltzmann equation has been an active research topic. Even though the stability of the two nonlinear wave patterns has been extensively studied, the stability of the linearly degenerate contact wave was not solved until recently. In this paper, we will briefly present our recent results in [32] on the stability of the contact wave patterns for both the Navier-Stokes equations and the Boltzmann equation.

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 prove the existence of a resonance-free region in scattering by a strictly convex obstacle $\mathcal{O}$ with the Robin boundary condition $\partial_{\nu}u+\gamma u|_{\partial\mathcal{O}}=0$ . More precisely, we show that the scattering resonances lie below a cubic curve ℑ$ζ$=−$S$|$ζ$|^{1/3}+$C$. The constant$S$is the same as in the case of the Neumann boundary condition$γ$=0. This generalizes earlier results on cubic pole-free regions obtained for the Dirichlet boundary condition.

Motivated by a model proposed by Peng et al. [Advances in Coll. and Interf. Sci. 206 (2014)] for break-up of tear films on human eyes, we study the dynamics of a generalized thin film model. The governing equations form a fourth-order coupled system of nonlinear parabolic PDE for the film thickness and salt concentration subject to nonconservative effects representing evaporation. We analytically prove the global existence of solutions to this model with mobility exponents in several different ranges and the results are then validated against PDE simulations. We also numerically capture other interesting dynamics of the model, including finite-time rupture-shock phenomenon due to the instabilities caused by locally elevated evaporation rates, convergence to equilibrium and infinite-time thinning.

No abstract uploaded!