-
Convergence analysis of the Monte Carlo method for random Navier--Stokes--Fourier system
Authors:
Maria Lukacova -- Medvidova,
Bangwei She,
Yuhuan Yuan
Abstract:
In the present paper we consider the initial data, external force, viscosity coefficients, and heat conductivity coefficient as random data for the compressible Navier--Stokes--Fourier system. The Monte Carlo method, which is frequently used for the approximation of statistical moments, is combined with a suitable deterministic discretisation method in physical space and time. Under the assumption…
▽ More
In the present paper we consider the initial data, external force, viscosity coefficients, and heat conductivity coefficient as random data for the compressible Navier--Stokes--Fourier system. The Monte Carlo method, which is frequently used for the approximation of statistical moments, is combined with a suitable deterministic discretisation method in physical space and time. Under the assumption that numerical densities and temperatures are bounded in probability, we prove the convergence of random finite volume solutions to a statistical strong solution by applying genuine stochastic compactness arguments. Further, we show the convergence and error estimates for the Monte Carlo estimators of the expectation and deviation. We present several numerical results to illustrate the theoretical results.
△ Less
Submitted 2 April, 2023;
originally announced April 2023.
-
Convergence and error analysis of compressible fluid flows with random data: Monte Carlo method
Authors:
Eduard Feireisl,
Mária Lukáčová - Medviďová,
Bangwei She,
Yuhuan Yuan
Abstract:
The goal of this paper is to study convergence and error estimates of the Monte Carlo method for the Navier-Stokes equations with random data. To discretize in space and time, the Monte Carlo method is combined with a suitable deterministic discretization scheme, such as a finite volume method. We assume that the initial data, force and the viscosity coefficients are random variables and study bot…
▽ More
The goal of this paper is to study convergence and error estimates of the Monte Carlo method for the Navier-Stokes equations with random data. To discretize in space and time, the Monte Carlo method is combined with a suitable deterministic discretization scheme, such as a finite volume method. We assume that the initial data, force and the viscosity coefficients are random variables and study both, the statistical convergence rates as well as the approximation errors. Since the compressible Navier-Stokes equations are not known to be uniquely solvable in the class of global weak solutions, we cannot apply pathwise arguments to analyze the random Navier-Stokes equations. Instead we have to apply intrinsic stochastic compactness arguments via the Skorokhod representation theorem and the Gyöngy-Krylov method. Assuming that the numerical solutions are bounded in probability, we prove that the Monte Carlo finite volume method converges to a statistical strong solution. The convergence rates are discussed as well. Numerical experiments illustrate theoretical results.
△ Less
Submitted 9 May, 2022;
originally announced May 2022.
-
Finite Volume Evolution Galerkin Methods for the Shallow Water Equations with Dry Beds
Authors:
Andreas Bollermann,
Sebastian Noelle,
Maria Lukáčová - Medvidová
Abstract:
We present a new Finite Volume Evolution Galerkin (FVEG) scheme for the solution of the shallow water equations (SWE) with the bottom topography as a source term. Our new scheme will be based on the FVEG methods presented in (Lukáčová, Noelle and Kraft, J. Comp. Phys. 221, 2007), but adds the possibility to handle dry boundaries. The most important aspect is to preserve the positivity of the water…
▽ More
We present a new Finite Volume Evolution Galerkin (FVEG) scheme for the solution of the shallow water equations (SWE) with the bottom topography as a source term. Our new scheme will be based on the FVEG methods presented in (Lukáčová, Noelle and Kraft, J. Comp. Phys. 221, 2007), but adds the possibility to handle dry boundaries. The most important aspect is to preserve the positivity of the water height. We present a general approach to ensure this for arbitrary finite volume schemes. The main idea is to limit the outgoing fluxes of a cell whenever they would create negative water height. Physically, this corresponds to the absence of fluxes in the presence of vacuum. Well-balancing is then re-established by splitting gravitational and gravity driven parts of the flux. Moreover, a new entropy fix is introduced that improves the reproduction of sonic rarefaction waves.
△ Less
Submitted 15 January, 2015;
originally announced January 2015.
-
Well-balanced finite volume evolution Galerkin methods for the shallow water equations
Authors:
Maria Lukáčová - Medvidová,
Sebastian Noelle,
Marcus Kraft
Abstract:
We present a new well-balanced finite volume method within the framework of the finite volume evolution Galerkin (FVEG) schemes. The methodology will be illustrated for the shallow water equations with source terms modelling the bottom topography and Coriolis forces. Results can be generalized to more complex systems of balance laws. The FVEG methods couple a finite volume formulation with approxi…
▽ More
We present a new well-balanced finite volume method within the framework of the finite volume evolution Galerkin (FVEG) schemes. The methodology will be illustrated for the shallow water equations with source terms modelling the bottom topography and Coriolis forces. Results can be generalized to more complex systems of balance laws. The FVEG methods couple a finite volume formulation with approximate evolution operators. The latter are constructed using the bicharacteristics of multidimensional hyperbolic systems, such that all of the infinitely many directions of wave propagation are taken into account explicitly. We derive a well- balanced approximation of the integral equations and prove that the FVEG scheme is well-balanced for the stationary steady states as well as for the steady jets in the rotational frame. Several numerical experiments for stationary and quasi-stationary states as well as for steady jets confirm the reliability of the well-balanced FVEG scheme.
△ Less
Submitted 15 January, 2015;
originally announced January 2015.