Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEnergy consistent DG methods for the Navier-Stokes-Korteweg system
187
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We design consistent discontinuous Galerkin finite element schemes for the approximation of the Euler-Korteweg and the Navier-Stokes-Korteweg systems. We show that the scheme for the Euler-Korteweg system is energy and mass conservative and that the scheme for the Navier-Stokes-Korteweg system is mass conservative and monotonically energy dissipative. In this case the dissipation is isolated to viscous effects, that is, there is no numerical dissipation. In this sense the methods is consistent with the energy dissipation of the continuous PDE systems.
rate research
Read More
We propose an efficient, accurate and robust implicit solver for the incompressible Navier-Stokes equations, based on a DG spatial discretization and on the TR-BDF2 method for time discretization. The effectiveness of the method is demonstrated in a number of classical benchmarks, which highlight its superior efficiency with respect to other widely used implicit approaches. The parallel implementation of the proposed method in the framework of the deal.II software package allows for accurate and efficient adaptive simulations in complex geometries, which makes the proposed solver attractive for large scale industrial applications.
This paper presents an arbitrary h.o. accurate ADER DG method on space-time adaptive meshes (AMR) for the solution of two important families of non-linear time dependent PDE for compr. dissipative flows: the compr. Navier-Stokes equations and the equations of visc. and res. MHD in 2 and 3 space-dimensions. The work continues a recent series of papers concerning the development and application of a proper a posteriori subcell FV limiting procedure suitable for DG methods. It is a well known fact that a major weakness of h.o. DG methods lies in the difficulty of limiting discontinuous solutions, which generate spurious oscillations, namely the so-called Gibbs phenomenon. In the present work the main benefits of the MOOD paradigm, i.e. the computational robustness even in the presence of strong shocks, are preserved and the numerical diffusion is considerably reduced also for the limited cells by resorting to a proper sub-grid. An important feature of our new scheme is its ability to cure even floating point errors that may occur during a simulation, for example when taking real roots of negative numbers or after divisions by zero. We apply the whole approach for the first time to the equations of compr. gas dynamics and MHD in the presence of viscosity, thermal conductivity and magnetic resistivity, therefore extending our family of adaptive ADER-DG schemes to cases for which the numerical fluxes also depend on the gradient of the state vector. The distinguished high-resolution properties of the presented numerical scheme stands out against a wide number of non-trivial test cases both for the compr. Navier-Stokes and the viscous and resistive MHD equations. The present results show clearly that the shock-capturing capability of the news schemes are significantly enhanced within a cell-by-cell Adaptive Mesh Refinement implementation together with time accurate local time stepping (LTS).
In this paper, the main objective is to generalize to the Navier-Stokes-Korteweg (with density dependent viscosities satisfying the BD relation) and Euler-Korteweg systems a recent relative entropy [proposed by D. Bresch, P. Noble and J.--P. Vila, (2016)] introduced for the compressible Navier-Stokes equations with a linear density dependent shear viscosity and a zero bulk viscosity. As a concrete application, this helps to justify mathematically the convergence between global weak solutions of the quantum Navier-Stokes system [recently obtained simultaneously by I. Lacroix-Violet and A. Vasseur (2017)] and dissipative solutions of the quantum Euler system when the viscosity coefficient tends to zero: This selects a dissipative solution as the limit of a viscous system. We also get weak-strong uniqueness for the Quantum-Euler and for the Quantum-Navier-Stokes equations. Our results are based on the fact that Euler-Korteweg systems and corresponding Navier--Stokes-Korteweg systems can be reformulated through an augmented system such as the compressible Navier-Stokes system with density dependent viscosities satisfying the BD algebraic relation. This was also observed recently [by D. Bresch, F. Couderc, P. Noble and J.--P. Vila, (2016)] for the Euler-Korteweg system for numerical purposes. As a by-product of our analysis, we show that this augmented formulation helps to define relative entropy estimates for the Euler Korteweg systems in a simplest way compared to recent works [See D. Donatelli, E. Feireisl, P. Marcati (2015) and J. Giesselmann, C. Lattanzio, A.-E. Tzavaras (2017)] with less hypothesis required on the capillary coefficient.
Discontinuous Galerkin (DG) methods have a long history in computational physics and engineering to approximate solutions of partial differential equations due to their high-order accuracy and geometric flexibility. However, DG is not perfect and there remain some issues. Concerning robustness, DG has undergone an extensive transformation over the past seven years into its modern form that provides statements on solution boundedness for linear and nonlinear problems. This chapter takes a constructive approach to introduce a modern incarnation of the DG spectral element method for the compressible Navier-Stokes equations in a three-dimensional curvilinear context. The groundwork of the numerical scheme comes from classic principles of spectral methods including polynomial approximations and Gauss-type quadratures. We identify aliasing as one underlying cause of the robustness issues for classical DG spectral methods. Removing said aliasing errors requires a particular differentiation matrix and careful discretization of the advective flux terms in the governing equations.
Hydrodynamics coupled phase field models have intricate difficulties to solve numerically as they feature high nonlinearity and great complexity in coupling. In this paper, we propose two second order, linear, unconditionally stable decoupling methods based on the Crank--Nicolson leap-frog time discretization for solving the Allen--Cahn--Navier--Stokes (ACNS) phase field model of two-phase incompressible flows. The ACNS system is decoupled via the artificial compression method and a splitting approach by introducing an exponential scalar auxiliary variable. We prove both algorithms are unconditionally long time stable. Numerical examples are provided to verify the convergence rate and unconditional stability.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
