Saray Busto

This paper proposes a semi-implicit arbitrary Lagrangian-Eulerian (ALE) method for the solution of the unified Godunov-Peshkov-Romenski (GPR) model of continuum mechanics. To handle the curl free involutions arising in the solid limit of the model, the original system is augmented by adopting a thermodynamically compatible generalized Lagrangian multiplier (GLM) approach. Next, an operator splitting strategy decouples the computation of fast pressure waves from the bulk velocity of the medium yielding a transport subsystem, containing convective terms and non-conservative products, and a Poisson-type subsystem, for the pressure. A second splitting yields an ODE subsystem comprising only the potentially stiff source terms, responsible for the relaxation of the model between its fluid and solid limits. The mesh motion can be driven by two sources: the local fluid velocity and a prescribed boundary displacement. For the spatial discretization, we employ unstructured staggered grids, with the pressure defined on the primal mesh and all remaining variables on the dual grid. The transport subsystem is advanced via an explicit finite volume method, in which integration over closed space-time control volumes ensures verification of the geometric conservation law (GCL). On the other hand, implicit continuous finite elements are used for the discretization of the pressure subsystem and an implicit DIRK scheme is employed to solve the ODE subsystem. Consequently, the proposed approach is well suited to address all Mach number flows. A comprehensive set of benchmarks is employed to assess the accuracy and robustness of the proposed methodology in tackling both fluid and solid mechanics problems.

We present a novel high order semi-implicit hybrid finite volume/virtual element numerical scheme for the solution of compressible flows on Voronoi tessellations. The method relies on the flux splitting of the compressible Navier-Stokes equations into three sub-systems: a convective sub-system solved explicitly using a finite volume (FV) scheme, and the viscous and pressure sub-systems which are discretized implicitly at the aid of a virtual element method (VEM). Consequently, the time step restriction of the overall algorithm depends only on the mean flow velocity and not on the fast pressure waves nor on the viscous eigenvalues. As such, the proposed methodology is well suited for the solution of low Mach number flows at all Reynolds numbers. Moreover, the scheme is proven to be globally energy conserving so that shock capturing properties are retrieved in high Mach number flows. To reach high order of accuracy in time and space, an IMEX Runge-Kutta time stepping strategy is employed together with high order spatial reconstructions in terms of CWENO polynomials and virtual element space basis functions. The chosen discretization techniques allow the use of general polygonal grids, a useful tool when dealing with complex domain configurations. The new scheme is carefully validated in both the incompressible limit and the high Mach number regime through a large set of classical benchmarks for fluid dynamics, assessing robustness and accuracy.

In this paper, a new efficient, and at the same time, very simple and general class of thermodynamically compatible finite volume schemes is introduced for the discretization of nonlinear, overdetermined, and thermodynamically compatible first-order hyperbolic systems. By construction, the proposed semi-discrete method satisfies an entropy inequality and is nonlinearly stable in the energy norm. A very peculiar feature of our approach is that entropy is discretized directly, while total energy conservation is achieved as a mere consequence of the thermodynamically compatible discretization. The new schemes can be applied to a very general class of nonlinear systems of hyperbolic PDEs, including both, conservative and non-conservative products, as well as potentially stiff algebraic relaxation source terms, provided that the underlying system is overdetermined and therefore satisfies an additional extra conservation law, such as the conservation of total energy density. The proposed family of finite volume schemes is based on the seminal work of Abgrall [1], where for the first time a completely general methodology for the design of thermodynamically compatible numerical methods for overdetermined hyperbolic PDE was presented. We apply our new approach to three particular thermodynamically compatible systems: the equations of ideal magnetohydrodynamics (MHD) with thermodynamically compatible generalized Lagrangian multiplier (GLM) divergence cleaning, the unified first-order hyperbolic model of continuum mechanics proposed by Godunov, Peshkov, and Romenski (GPR model) and the first-order hyperbolic model for turbulent shallow water flows of Gavrilyuk et al. In addition to formal mathematical proofs of the properties of our new finite volume schemes, we also present a large set of numerical results in order to show their potential, efficiency, and practical applicability.

We present a novel fully implicit hybrid finite volume/finite element method for incompressible flows. Following previous works on semi-implicit hybrid FV/FE schemes, the incompressible Navier-Stokes equations are split into a pressure and a transport-diffusion subsystem. The first of them can be seen as a Poisson type problem and is thus solved efficiently using classical continuous Lagrange finite elements. On the other hand, finite volume methods are employed to solve the convective subsystem, in combination with Crouzeix-Raviart finite elements for the discretization of the viscous stress tensor. For some applications, the related CFL condition, even if depending only in the bulk velocity, may yield a severe time restriction in case explicit schemes are used. To overcome this issue an implicit approach is proposed. The system obtained from the implicit discretization of the transport-diffusion operator is solved using an inexact Newton-Krylov method, based either on the BiCStab or the GMRES algorithm. To improve the convergence properties of the linear solver a symmetric Gauss-Seidel (SGS) preconditioner is employed, together with a simple but efficient approach for the reordering of the grid elements that is compatible with MPI parallelization. Besides, considering the Ducros flux for the nonlinear convective terms we can prove that the discrete advection scheme is kinetic energy stable. The methodology is carefully assessed through a set of classical benchmarks for fluid mechanics. A last test shows the potential applicability of the method in the context of blood flow simulation in realistic vessel geometries.

We construct, analyse and assess various schemes of second order of accuracy in space and time for model advection-diffusion-reaction differential equations. The constructed schemes are meant to be of practical use in solving industrial problems and are derived following two related approaches, namely ADER and MUSCL-Hancock. Detailed analysis of linear stability and local truncation error are carried out. In addition, the schemes are implemented and assessed for various test problems. Empirical convergence rate studies confirm the theoretically expected accuracy in both space and time.