We establish an optimal error estimate for a random particle blob method for the Keller-Segel equation in Rd (d  2). With a blob size " = h

Linear elliptic equations in composite media with anisotropic fibres are concerned. The media consist of a periodic set of anisotropic fibres with low conductivity, included in a connected matrix with high conductivity. Inside the anisotropic fibres, the conductivity in the longitudinal direction is relatively high compared with that in the transverse directions. The coefficients of the elliptic equations depend on the conductivity. This work is to derive the H\"older and the gradient $L^p$ estimates (uniformly in the period size of the set of anisotropic fibres as well as in the conductivity ratio of the fibres in the transverse directions to the connected matrix) for the solutions of the elliptic equations. Furthermore, it is shown that, inside the fibres, the solutions have higher regularity along the fibres than in the transverse directions.

We study the asymptotic spreading of KolmogorovPetrovskyPiskunov (KPP) fronts in spacetime random incompressible flows in dimension d> 1. We prove that if the flow field is stationary, ergodic, and obeys a suitable moment condition, the large time front speeds (spreading rates) are deterministic in all directions for compactly supported initial data. The flow field can become unbounded at large times. The front speeds are characterized by the convex rate function governing large deviations of the associated diffusion in the random flow. Our proofs are based on the Harnack inequality, an application of the sub-additive ergodic theorem, and the construction of comparison functions. Using the variational principles for the front speed, we obtain general lower and upper bounds of front speeds in terms of flow statistics. The bounds show that front speed enhancement in incompressible flows can grow at most linearly in the root mean square amplitude of the flows, and may have much slower growth due to rapid temporal decorrelation of the flows.

In this work, we consider the numerical solution of the nonlinear Schrödinger equation with a highly oscillatory potential (NLSE-OP). The NLSE-OP is a model problem which frequently occurs in recent studies of some multiscale dynamical systems, where the potential introduces wide temporal oscillations to the solution and causes numerical difficulties. We aim to analyze rigorously the error bounds of the splitting schemes for solving the NLSE-OP to a fixed time. Our theoretical results show that the Lie–Trotter splitting scheme is uniformly and optimally accurate at the first order provided that the oscillatory potential is integrated exactly, while the Strang splitting scheme is not. Our results apply to general dispersive or wave equations with an oscillatory potential. The error estimates are confirmed by numerical results.

In this paper, we investigate the initial-boundary value problem to the heat conductive compressible Navier-Stokes equations. Local existence and uniqueness of strong solutions is established with any such initial data that the initial density $\rho_0$, velocity $u_0$, and temperature $\theta_0$ satisfy $\rho_0\in W^{1,q}$, with $q\in(3,6)$, $u_0\in H^1$, and $\sqrt{\rho_0}\theta_0\in L^2$. The initial density is assumed to be only nonnegative and thus the initial vacuum is allowed. In addition to the necessary regularity assumptions, we do not require any initial compatibility conditions such as those proposed in (Y. Cho and H. Kim, \emph{Existence results for viscous polytropic fluids with vacuum}, J. Differential Equations {\bf 228} (2006), no.~2, 377--411.), which although are widely used in many previous works but put some inconvenient constraints on the initial data. Due to the weaker regularities of the initial data and the absence of the initial compatibility conditions, leading to weaker regularities of the solutions compared with those in the previous works, the uniqueness of solutions obtained in the current paper does not follow from the arguments used in the existing literatures. Our proof of the uniqueness of solutions is based on the following new idea of two-stages argument: (i) showing that the difference of two solutions (or part of their components) with the same initial data is controlled by some power function of the time variable; (ii) carrying out some singular-in-time weighted energy differential inequalities fulfilling the structure of the Gr\"onwall inequality. The existence is established in the Euler coordinates, while the uniqueness is proved in the Lagrangian coordinates first and then transformed back to the Euler coordinates.