Search

We present a finite difference-(multi-level) Monte Carlo algorithm to efficiently compute statistical solutions of the two-dimensional incompressible Navier-Stokes equations (NSE), with periodic boundary conditions and for arbitrary high Reynolds number. We propose a reformulation of statistical solutions in the sense of Foiaş and Prodi in the vorticity-stream function form. The vorticity-stream function formulation of the NSE in two-space dimensions is discretized with a finite difference scheme. We obtain a convergence rate error estimate for this approximation which is explicit in the viscosity parameter , under realistic assumptions on the solution regularity. We also prove convergence and complexity estimates, for the (multi-level) Monte Carlo finite difference algorithm to compute statistical solutions. Numerical experiments illustrating the validity of our estimates are presented. They show that the multi-level Monte Carlo algorithm can significantly accelerate the computation of statistical solutions in the sense of Foiaş and Prodi, even for very high Reynolds numbers. [ABSTRACT FROM AUTHOR]

We consider a sub-model of the Hall-MHD equations: the so-called magnetic induction equations with Hall effect. These equations are non-linear and include third-order spatial and spatio-temporal mixed derivatives. We show that the energy of the solutions is bounded and design finite difference schemes that preserve the energy bounds for the continuous problem. We design both divergence preserving schemes and schemes with bounded divergence. We present a set of numerical experiments that demonstrate the robustness of the proposed schemes. [ABSTRACT FROM AUTHOR]

We consider the scalar conservation law with flux function discontinuous in the space variable, i.e., \begin{eqnarray} \label{eq1} u_t+(H(x)f(u)+(1-H(x))g(u))_{x} &=& 0 \quad \mbox{in } \R \times \R_{+}, \nonumber \\ u(0, x) &=& u_{0}(x) \quad \mbox{in } \R, \label{0.1} \end{eqnarray} where $H$ is the Heaviside function and $f$ and $g$ are smooth with the assumptions that either $f$ is convex and $g$ is concave or $f$ is concave and $g$ is convex. The existence of a weak solution of (\ref{eq1}) is proved by showing that upwind finite difference schemes of Godunov and Enquist--Osher type converge to a weak solution. Uniqueness follows from a Kruzkhov-type entropy condition. We also provide explicit solutions to the Riemann problem for (\ref{eq1}). At the level of numerics, we give easy-to-implement numerical schemes of Godunov and Enquist--Osher type. The central feature of this paper is the modification of the singular mapping technique (the main analytical tool for these types of equations) which allows us to show that the numerical schemes converge. Equations of type (\ref{eq1}) with the above hypothesis on the flux may occur when considering the following scalar conservation law with discontinuous flux: \begin{equation} \label{eq2} \begin{array}{r@{\;}l} u_t + (k (x) f (u))_x &= 0, \\ u (0, x) &= u_0 (x), \end{array} \end{equation} with $f$ convex and $k$ of indefinite sign. [ABSTRACT FROM AUTHOR]