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 combine the spectral (viscosity) method and ensemble averaging to propose an algorithm that computes admissible measure-valued solutions of the incompressible Euler equations. The resulting approximate young measures are proved to converge (with increasing numerical resolution) to a measure-valued solution. We present numerical experiments demonstrating the robustness and efficiency of the proposed algorithm, as well as the appropriateness of measure-valued solutions as a solution framework for the Euler equations. Furthermore, we report an extensive computational study of the two-dimensional vortex sheet, which indicates that the computed measure-valued solution is non-atomic and implies possible non-uniqueness of weak solutions constructed by Delort. [ABSTRACT FROM AUTHOR]

Solutions of initial-boundary value problems for systems of conservation laws depend on the underlying viscous mechanism, namely different viscosity operators lead to different limit solutions. Standard numerical schemes for approximating conservation laws do not take into account this fact and converge to solutions that are not necessarily physically relevant. We design numerical schemes that incorporate explicit information about the underlying viscosity mechanism and approximate the physical-viscosity solution. Numerical experiments illustrating the robust performance of these schemes are presented. [ABSTRACT FROM AUTHOR]

We deal with a single conservation law in one space dimension whose flux function is discontinuous in the space variable and we introduce a proper framework of entropy solutions. We consider a large class of fluxes, namely, fluxes of the convex-convex type and of the concave-convex (mixed) type. The alternative entropy framework that is proposed here is based on a two step approach. In the first step, infinitely many classes of entropy solutions are defined, each associated with an interface connection. We show that each of these class of entropy solutions form a contractive semigroup in L1 and is hence unique. Godunov type schemes based on solutions of the Riemann problem are designed and shown to converge to each class of these entropy solutions. The second step is to choose one of these classes of solutions. This choice depends on the Physics of the problem being considered and we concentrate on the model of two-phase flows in a heterogeneous porous medium. We define an optimization problem on the set of admissible interface connections and show the existence of an unique optimal connection and its corresponding optimal entropy solution. The optimal entropy solution is consistent with the expected solutions for two-phase flows in heterogeneous porous media. [ABSTRACT FROM AUTHOR]