Lagrangian ADER-WENO Finite Volume Schemes on Unstructured Triangular Meshes Based On Genuinely Multidimensional HLL Riemann Solvers

In this paper we use the genuinely multidimensional HLL Riemann solvers recently developed by Balsara et al. to construct a new class of computationally efficient high order Lagrangian ADER-WENO one-step ALE finite volume schemes on unstructured triangular meshes. A nonlinear WENO reconstruction operator allows the algorithm to achieve high order of accuracy in space, while high order of accuracy in time is obtained by the use of an ADER time-stepping technique based on a local space-time Galerkin predictor. The multidimensional HLL and HLLC Riemann solvers operate at each vertex of the grid, considering the entire Voronoi neighborhood of each node and allows for larger time steps than conventional one-dimensional Riemann solvers. The results produced by the multidimensional Riemann solver are then used twice in our one-step ALE algorithm: first, as a node solver that assigns a unique velocity vector to each vertex, in order to preserve the continuity of the computational mesh; second, as a building block for genuinely multidimensional numerical flux evaluation that allows the scheme to run with larger time steps compared to conventional finite volume schemes that use classical one-dimensional Riemann solvers in normal direction. A rezoning step may be necessary in order to overcome element overlapping or crossing-over. We apply the method presented in this article to two systems of hyperbolic conservation laws, namely the Euler equations of compressible gas dynamics and the equations of ideal classical magneto-hydrodynamics (MHD). Convergence studies up to fourth order of accuracy in space and time have been carried out. Several numerical test problems have been solved to validate the new approach

Simulating magnetized neutron stars with discontinuous Galerkin methods

Discontinuous Galerkin methods are popular because they can achieve high order where the solution is smooth, because they can capture shocks while needing only nearest-neighbor communication, and because they are relatively easy to formulate on complex meshes. We perform a detailed comparison of various limiting strategies presented in the literature applied to the equations of general relativistic magnetohydrodynamics. We compare the standard minmod/$\Lambda\Pi^N$ limiter, the hierarchical limiter of Krivodonova, the simple WENO limiter, the HWENO limiter, and a discontinuous Galerkin-finite-difference hybrid method. The ultimate goal is to understand what limiting strategies are able to robustly simulate magnetized TOV stars without any fine-tuning of parameters. Among the limiters explored here, the only limiting strategy we can endorse is a discontinuous Galerkin-finite-difference hybrid method

An explicit high-order single-stage single-step positivity-preserving finite difference WENO method for the compressible Euler equations

In this work we construct a high-order, single-stage, single-step positivity-preserving method for the compressible Euler equations. Space is discretized with the finite difference weighted essentially non-oscillatory (WENO) method. Time is discretized through a Lax-Wendroff procedure that is constructed from the Picard integral formulation (PIF) of the partial differential equation. The method can be viewed as a modified flux approach, where a linear combination of a low- and high-order flux defines the numerical flux used for a single-step update. The coefficients of the linear combination are constructed by solving a simple optimization problem at each time step. The high-order flux itself is constructed through the use of Taylor series and the Cauchy-Kowalewski procedure that incorporates higher-order terms. Numerical results in one- and two-dimensions are presented

Monotonicity-preserving finite element methods for hyperbolic problems

This thesis covers the development of monotonicity preserving finite element methods for hyperbolic problems. In particular, scalar convection-diffusion and Euler equations are used as model problems for the discussion in this dissertation. A novel artificial diffusion stabilization method has been proposed for scalar problems. This technique is proved to yield monotonic solutions, to be \ac{led}, Lipschitz continuous, and linearity preserving. These properties are satisfied in multiple dimensions and for general meshes. However, these results are limited to first order Lagrangian finite elements. A modification of this stabilization operator that is twice differentiable has been also proposed. With this regularized operator, nonlinear convergence is notably improved, while the stability properties remain unaltered (at least, in a weak sense). An extension of this stabilization method to high-order discretizations has also been proposed. In particular, arbitrary order space-time isogeometric analysis is used for this purpose. It has been proved that this scheme yields solutions that satisfy a global space-time discrete maximum principle unconditionally. A partitioned approach has also been proposed. This strategy reduces the computational cost of the scheme, while it preserves all stability properties. A regularization of this stabilization operator has also been developed. As for the first order finite element method, it improves the nonlinear convergence without harming the stability properties. An extension to Euler equations has also been pursued. In this case, instead of monotonicity-preserving, the developed scheme is local bounds preserving. Following the previous works, a regularized differentiable version has also been proposed. In addition, a continuation method using the parameters introduced for the regularization has been used. In this case, not only the nonlinear convergence is improved, but also the robustness of the method. However, the improvement in nonlinear convergence is limited to moderate tolerances and it is not as notable as for the scalar problem. Finally, the stabilized schemes proposed had been adapted to adaptive mesh refinement discretizations. In particular, nonconforming hierarchical octree-based meshes have been used. Using these settings, the efficiency of solving a monotonicity-preserving high-order stiff nonlinear problem has been assessed. Given a specific accuracy, the computational time required for solving the high-order problem is compared to the one required for solving a low-order problem (easy to converge) in a much finer adapted mesh. In addition, an error estimator based on the stabilization terms has been proposed and tested. The performance of all proposed schemes has been assessed using several numerical tests and solving various benchmark problems. The obtained results have been commented and included in the dissertation.La present tesi tracta sobre mètodes d'elements finits que preserven la monotonia per a problemes hiperbòlics. Concretament, els problemes que s'han utilitzat com a model en el desenvolupament d'aquesta tesi són l'equació escalar de convecció-difusió-reacció i les equacions d'Euler. Per a problemes escalars s'ha proposat un nou mètode d'estabilització mitjançant difusió artificial. S'ha provat que amb aquesta tècnica les solucions obtingudes són monòtones, l'esquema "disminueix els extrems locals", i preserva la linearitat. Aquestes propietats s'han pogut demostrar per múltiples dimensions i per malles generals. Per contra, aquests resultats només són vàlids per elements finits Lagrangians de primer ordre. També s'ha proposat una modificació de l'operador d'estabilització per tal de que aquest sigui diferenciable. Aquesta regularització ha permès millorar la convergència no-lineal notablement, mentre que les propietats d'estabilització no s'han vist alterades. L'anterior mètode d'estabilització s'ha adaptat a discretitzacions d'alt ordre. Concretament, s'ha utilitzat anàlisi isogeomètrica en espai i temps per a aquesta tasca. S'ha provat que les solucions obtingudes mitjançant aquest mètode satisfan el principi del màxim discret de forma global. També s'ha proposat un esquema particionat. Aquesta alternativa redueix el cost computacional, mentre preserva totes les propietats d'estabilitat. En aquest cas, també s'ha realitzat una regularització de l'operador d'estabilització per tal de que sigui diferenciable. Tal i com s'ha observat en els mètodes de primer ordre, aquesta regularització permet millorar la convergència no-lineal sense perdre les propietats d'estabilització. Posteriorment, s'ha estudiat l'adaptació dels mètodes anteriors a les equacions d'Euler. En aquest cas, en comptes de preservar la monotonia, l'esquema preserva "cotes locals". Seguint els desenvolupaments anteriors, s'ha proposat una versió diferenciable de l'estabilització. En aquest cas, també s'ha desenvolupat un mètode de continuació utilitzant els paràmetres introduïts per a la regularització. En aquest cas, no només ha millorat la convergència no-lineal sinó que l'esquema també esdevé més robust. Per contra, la millora en la convergència no-lineal només s'observa per a toleràncies moderades i no és tan notable com en el cas dels problemes escalars. Finalment, els esquemes d'estabilització proposat s'han adaptat a malles de refinament adaptatiu. Concretament, s'han utilitzat malles no-conformes basades en octrees. Utilitzant aquesta configuració, l'eficiència de resoldre un problema altament no-lineal ha estat avaluada de la següent forma. Donada una precisió determinada, el temps computacional requerit per resoldre el problema utilitzant un esquema d'alt ordre ha estat comparat amb el temps necessari per resoldre'l utilitzant un esquema de baix ordre en una malla adaptativa molt més refinada. Addicionalment, també s'ha proposat un estimador de l'error basat en l'operador d'estabilització. El comportament de tots els esquemes proposats anteriorment s'ha avaluat mitjançant varis tests numèrics. Els resultats s'han compilat i comentat en la present tesi.Postprint (published version

Hyperbolic Balance Laws: modeling, analysis, and numerics (hybrid meeting)

This workshop brought together leading experts, as well as the most promising young researchers, working on nonlinear hyperbolic balance laws. The meeting focused on addressing new cutting-edge research in modeling, analysis, and numerics. Particular topics included ill-/well-posedness, randomness and multiscale modeling, flows in a moving domain, free boundary problems, games and control