https://hal.inria.fr/hal-03402701v2Kazolea, MariaMariaKazoleaCARDAMOM - Certified Adaptive discRete moDels for robust simulAtions of CoMplex flOws with Moving fronts - IMB - Institut de Mathématiques de Bordeaux - Université Bordeaux Segalen - Bordeaux 2 - UB - Université Sciences et Technologies - Bordeaux 1 - UB - Université de Bordeaux - Bordeaux INP - Institut Polytechnique de Bordeaux - CNRS - Centre National de la Recherche Scientifique - Inria Bordeaux - Sud-Ouest - Inria - Institut National de Recherche en Informatique et en AutomatiqueFilippini, AndreaAndreaFilippiniBRGM - Bureau de Recherches Géologiques et Minières (BRGM)Ricchiuto, MarioMarioRicchiutoCARDAMOM - Certified Adaptive discRete moDels for robust simulAtions of CoMplex flOws with Moving fronts - IMB - Institut de Mathématiques de Bordeaux - Université Bordeaux Segalen - Bordeaux 2 - UB - Université Sciences et Technologies - Bordeaux 1 - UB - Université de Bordeaux - Bordeaux INP - Institut Polytechnique de Bordeaux - CNRS - Centre National de la Recherche Scientifique - Inria Bordeaux - Sud-Ouest - Inria - Institut National de Recherche en Informatique et en AutomatiqueLow dispersion finite volume/element discretization of the enhanced Green-Naghdi equations for wave propagation, breaking and runup on unstructured meshesHAL CCSD2023Green-Naghdi equationsHybrid schemeFinite VolumesFinite ElementsHigh orderWave breaking[MATH.MATH-NA] Mathematics [math]/Numerical Analysis [math.NA][MATH.MATH-MP] Mathematics [math]/Mathematical Physics [math-ph]Kazolea, Maria2023-01-05 16:31:472023-02-01 18:44:102023-01-06 09:08:15enJournal articleshttps://hal.inria.fr/hal-03402701v2/document10.1016/j.ocemod.2022.102157https://hal.inria.fr/hal-03402701v1application/pdf2We study a hybrid approach combining a FV and FE method to solve a fully nonlinear and weakly-dispersive depth averaged wave propagation model. The FV method is used to solve the underlying hyperbolic shallow water system, while a standard P 1 finite element method is used to solve the elliptic system associated to the dispersive correction. We study the impact of several numerical aspects: the impact of the reconstruction used in the hyperbolic phase; the representation of the FV data in the FE method used in the elliptic phase and their impact on the theoretical accuracy of the method; the well-posedness of the overall method. For the first element we proposed a systematic implementation of an iterative reconstruction providing on arbitrary meshes up to third order solutions, full second order first derivatives, as well as a consistent approximation of the second derivatives. These properties are exploited to improve the assembly of the elliptic solver, showing dramatic improvement of the finale accuracy, if the FV representation is correctly accounted for. Concerning the elliptic step, the original problem is usually better suited for an approximation in H(div) spaces. However, it has been shown that perturbed problems involving similar operators with a small Laplace perturbation are well behaved in H 1. We show, based on both heuristic and strong numerical evidence, that numerical dissipation plays a major role in stabilizing the coupled method, and not only providing convergent results, but also providing the expected convergence rates. Finally, the full mode, coupling a wave breaking closure previously developed by the authors, is thoroughly tested on standard benchmarks using unstructured grids with sizes comparable or coarser than those usually proposed in literature.