A semi-discrete large-time behavior preserving scheme for the augmented Burgers equation

In this paper we analyze the large-time behavior of the augmented Burgers equation. We first study the well-posedness of the Cauchy problem and obtain $L^1$-$L^p$ decay rates. The asymptotic behavior of the solution is obtained by showing that the influence of the convolution term $K*u_{xx}$ is the same as $u_{xx}$ for large times. Then, we propose a semi-discrete numerical scheme that preserves this asymptotic behavior, by introducing two correcting factors in the discretization of the non-local term. Numerical experiments illustrating the accuracy of the results of the paper are also presented.Comment: Accepted for publication in ESAIM: Mathematical Modelling and Numerical Analysi

Large-time asymptotics, vanishing viscosity and numerics for 1-D scalar conservation laws

In this paper we analyze the large time asymptotic behavior of the discrete solutions of numerical approximation schemes for scalar hyperbolic conservation laws. We consider three monotone conservative schemes that are consistent with the one-sided Lipschitz condition (OSLC): Lax-Friedrichs, Engquist-Osher and Godunov. We mainly focus on the inviscid Burgers equation, for which we know that the large time behavior is of self-similar nature, described by a two-parameter family of N-waves. We prove that, at the numerical level, the large time dynamics depends on the amount of numerical viscosity introduced by the scheme: while Engquist-Osher and Godunov yield the same N-wave asymptotic behavior, the Lax-Friedrichs scheme leads to viscous self-similar profiles, corresponding to the asymptotic behavior of the solutions of the continuous viscous Burgers equation. The same problem is analyzed in the context of self-similar variables that lead to a better numerical performance but to the same dichotomy on the asymptotic behavior: N-waves versus viscous ones. We also give some hints to extend the results to more general fluxes. Some numerical experiments illustrating the accuracy of the results of the paper are also presented.Comment: Error corrected in main theorem in v3. Obtained results do not change in essence. Minor typos corrected in v

Lower and upper bounds for the first eigenvalue of nonlocal diffusion problems in the whole space

We find lower and upper bounds for the first eigenvalue of a nonlocal diffusion operator of the form T(u) = - \int_{\rr^d} K(x,y) (u(y)-u(x)) \, dy. Here we consider a kernel $K(x,y)=\psi (y-a(x))+\psi(x-a(y))$ where $\psi$ is a bounded, nonnegative function supported in the unit ball and $a$ means a diffeomorphism on \rr^d. A simple example being a linear function $a(x)= Ax$. The upper and lower bounds that we obtain are given in terms of the Jacobian of $a$ and the integral of $\psi$. Indeed, in the linear case $a(x) = Ax$ we obtain an explicit expression for the first eigenvalue in the whole \rr^d and it is positive when the the determinant of the matrix $A$ is different from one. As an application of our results, we observe that, when the first eigenvalue is positive, there is an exponential decay for the solutions to the associated evolution problem. As a tool to obtain the result, we also study the behaviour of the principal eigenvalue of the nonlocal Dirichlet problem in the ball $B_R$ and prove that it converges to the first eigenvalue in the whole space as $R\to \infty$