MathSciDoc: An Archive for Mathematician ∫

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.

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

This paper deals with a degenerate diffusion Patlak–Keller–Segel system in $n\geq 3$ dimension. The main difference between the current work and many other recent studies on the same model is that we study the diffusion exponent $m=2n/(n+2)$, which is smaller than the usual exponent $m^{*}=2-2/n$ used in other studies. With the exponent $m=2n/(n+2)$, the associated free energy is conformal invariant, and there is a family of stationary solutions $U_{\lambda,x_0}(x)=C(\frac{\lambda} {\lambda^2+|x-x_0|^2})^{\frac{n+2}{2}}$ $\forall \lambda>0$, $x_0\in {\mathbb R}^n$. For radially symmetric solutions, we prove that if the initial data are strictly below $U_{\lambda,0}(x)$ for some $\lambda$, then the solution vanishes in $L^1_{loc}$ as $t\to\infty$; if the initial data are strictly above $U_{\lambda,0}(x)$ for some $\lambda$, then the solution either blows up at a finite time or has a mass concentration at $r=0$ as time goes to infinity. For general initial data, we prove that there is a global weak solution provided that the $L^m$ norm of initial density is less than a universal constant, and the weak solution vanishes as time goes to infinity. We also prove a finite time blow-up of the solution if the $L^m$ norm for initial data is larger than the $L^m$ norm of $U_{\lambda,x_0}(x)$, which is constant independent of $\lambda$ and $x_0$, and the free energy of initial data is smaller than that of $U_{\lambda,x_0}(x)$.

Motivated by a phenomenon of phase transition in a model of alignment of self-propelled particles, we obtain a kinetic mean-field equation which is nothing else than the Doi equation (also called Smoluchowski equation) with dipolar potential. In a self-contained article, using only basic tools, we analyze the dynamics of this equation in any dimension. We first prove global well-posedness of this equation, starting with an initial condition in any Sobolev space. We then compute all possible steady-states. There is a threshold for the noise parameter: over this threshold, the only equilibrium is the uniform distribution, and under this threshold, there is also a family of non-isotropic equilibria. We give a rigorous prove of convergence of the solution to a steady-state as time goes to infinity. In particular we show that in the supercritical case, the only initial conditions leading to the uniform distribution in large time are those with vanishing momentum. For any positive value of the noise parameter, and any initial condition, we give rates of convergence towards equilibrium, exponentially for both supercritical and subcritical cases and algebraically for the critical case.

The existence of solution for the 2D-Keller-Segel system in the subcritical case, i.e. when the initial mass is less than 8π, is reproved. Instead of using the entropy in the free energy and free energy dissipation, which was used in the proofs (Blanchet et al. in SIAM J. Numer. Anal. 46:691–721, 2008; Electron. J. Differ. Equ. Conf. 44:32, 2006 (electronic)), the potential energy term is fully utilized by adapting Delort’s theory on 2D incompressible Euler equation (Delort in J. Am. Math. Soc. 4:553–386, 1991).

We consider a kinetic model of self-propelled particles with alignment interaction and with precession about the alignment direction. We derive a hydrodynamic system for the local density and velocity orientation of the particles. The system consists of the conservative equation for the local density and a non-conservative equation for the orientation. First, we assume that the alignment interaction is purely local and derive a first order system. However, we show that this system may lose its hyperbolicity. Under the assumption of weakly non-local interaction, we derive diffusive corrections to the first order system which lead to the combination of a heat flow of the harmonic map and Landau-Lifschitz-Gilbert dynamics. In the particular case of zero self-propelling speed, the resulting model reduces to the phenomenological Landau-Lifschitz-Gilbert equations. Therefore the present theory provides a kinetic formulation of classical micromagnetization models and spin dynamics.

The purpose of this paper is to provide global existence and uniqueness results for a family of fluid transport equations by establishing convergence results for the particle method applied to these equations. The considered family of PDEs is a collection of strongly nonlinear equations which yield traveling wave solutions and can be used to model a variety of flows in fluid dynamics. We apply a particle method to the studied evolutionary equations and provide a new self-contained method for proving its convergence. The latter is accomplished by using the concept of space-time bounded variation and the associated compactness properties. From this result, we prove the existence of a unique global weak solution in some special cases and obtain stronger regularity properties of the solution than previously established.

The computation of compressible flows becomes more challenging when the Mach number has different orders of magnitude. When the Mach number is of order one, modern shock capturing methods are able to capture shocks and other complex structures with high numerical resolutions. However, if the Mach number is small, the acoustic waves lead to stiffness in time and exces-sively large numerical viscosity, thus demanding much smaller time step and mesh size than normally needed for incompressible flow simulation. In this paper, we develop an all-speed asymptotic preserving (AP) numerical scheme for the compressible isentropic Euler and Navier-Stokes equations that is uni-formly stable and accurate for all Mach numbers. Our idea is to split the system into two parts: one involves a slow, nonlinear and conservative hyper-bolic system adequate for the use of modern shock capturing methods, and the other a linear hyperbolic system which contains the stiff acoustic dynam-ics, to be solved implicitly. This implicit part is reformulated into a standard pressure Poisson projection system, and thus possesses sufficient structure for efficient fast Fourier transform solution techniques. In the zero Mach number limit, the scheme automatically becomes a projection method-like incompress-ible solver. We present numerical results in one and two dimensions in both compressible and incompressible regimes.

The purpose of this paper is to establish a new method for proving the convergence of the particle method applied to the Camassa-Holm (CH) equation. The CH equation is a strongly nonlinear, bi-Hamiltonian, completely integrable model in the context of shallow water waves. The equation admits solutions that are nonlinear superpositions of traveling waves that have a discontinuity in the first derivative at their peaks and therefore are called peakons. This behavior admits several diverse scientific applications, but introduce difficult numerical challenges. To accurately capture these solutions, one may apply the particle method to the CH equation. Using the concept of space-time bounded variation, we show that the particle solution converges to a global weak solution of the CH equation for positive Radon measure initial data. Read More: http://www.worldscientific.com/doi/abs/10.1142/9789814417099_0033