Source author record

Manuel Quezada de Luna

Manuel Quezada de Luna appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
7topics
2close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2021arXiv

Numerical simulation and entropy dissipative cure of the carbuncle instability for the shallow water circular hydraulic jump

We investigate the numerical artifact known as a carbuncle, in the solution of the shallow water equations. We propose a new Riemann solver that is based on a local measure of the entropy residual and aims to avoid carbuncles while maintaining high accuracy. We propose a new challenging test problem for shallow water codes, consisting of a steady circular hydraulic jump that can be physically unstable. We show that numerical methods are prone to either suppress the instability completely or form carbuncles. We test existing cures for the carbuncle. In our experiments, only the proposed method is able to avoid unphysical carbuncles without suppressing the physical instability.

preprint2020arXiv

Algebraic entropy fixes and convex limiting for continuous finite element discretizations of scalar hyperbolic conservation laws

In this work, we modify a continuous Galerkin discretization of a scalar hyperbolic conservation law using new algebraic correction procedures. Discrete entropy conditions are used to determine the minimal amount of entropy stabilization and constrain antidiffusive corrections of a property-preserving low-order scheme. The addition of a second-order entropy dissipative component to the antidiffusive part of a nearly entropy conservative numerical flux is generally insufficient to prevent violations of local bounds in shock regions. Our monolithic convex limiting technique adjusts a given target flux in a manner which guarantees preservation of invariant domains, validity of local maximum principles, and entropy stability. The new methodology combines the advantages of modern entropy stable/entropy conservative schemes and their local extremum diminishing counterparts. The process of algebraic flux correction is based on inequality constraints which provably provide the desired properties. No free parameters are involved. The proposed algebraic fixes are readily applicable to unstructured meshes, finite element methods, general time discretizations, and steady-state residuals. Numerical studies of explicit entropy-constrained schemes are performed for linear and nonlinear test problems.

preprint2020arXiv

Effective Rankine-Hugoniot conditions for shock waves in periodic media

Solutions of first-order nonlinear hyperbolic conservation laws typically develop shocks in finite time even with smooth initial conditions. However, in heterogeneous media with rapid spatial variation, shock formation may be delayed or avoided. When shocks do form in such media, their speed of propagation depends on the material structure. We investigate conditions for shock formation and propagation in heterogeneous media. We focus on the propagation of plane waves in two-dimensional periodic media with material variation in only one direction. We propose an estimate for the speed of the shocks that is based on the Rankine-Hugoniot conditions applied to a leading-order homogenized (constant coefficient) system. We verify this estimate via numerical simulations using different nonlinear constitutive relations and layered and smoothly varying periodic media. In addition, we discuss conditions and regimes under which shocks form in this type of media.

preprint2020arXiv

Entropy conservation property and entropy stabilization of high-order continuous Galerkin approximations to scalar conservation laws

This paper addresses the design of linear and nonlinear stabilization procedures for high-order continuous Galerkin (CG) finite element discretizations of scalar conservation laws. We prove that the standard CG method is entropy conservative for the square entropy. In general, the rate of entropy production/dissipation depends on the residual of the governing equation and on the accuracy of the finite element approximation to the entropy variable. The inclusion of linear high-order stabilization generates an additional source/sink in the entropy budget equation. To balance the amount of entropy production in each cell, we construct entropy-dissipative element contributions using a coercive bilinear form and a parameter-free entropy viscosity coefficient. The entropy stabilization term is high-order consistent, and optimal convergence behavior is achieved in practice. To enforce preservation of local bounds in addition to entropy stability, we use the Bernstein basis representation of the finite element solution and a new subcell flux limiting procedure. The underlying inequality constraints ensure the validity of localized entropy conditions and local maximum principles. The benefits of the proposed modifications are illustrated by numerical results for linear and nonlinear test problems.

preprint2020arXiv

Subcell flux limiting for high-order Bernstein finite element discretizations of scalar hyperbolic conservation laws

This work extends the concepts of algebraic flux correction and convex limiting to continuous high-order Bernstein finite element discretizations of scalar hyperbolic problems. Using an array of adjustable diffusive fluxes, the standard Galerkin approximation is transformed into a nonlinear high-resolution scheme which has the compact sparsity pattern of the piecewise-linear or multilinear subcell discretization. The representation of this scheme in terms of invariant domain preserving states makes it possible to prove the validity of local discrete maximum principles under CFL-like conditions. In contrast to predictor-corrector approaches based on the flux-corrected transport methodology, the proposed flux limiting strategy is monolithic; i.e., limited antidiffusive terms are incorporated into the well-defined residual of a nonlinear (semi-)discrete problem. A stabilized high-order Galerkin discretization is recovered if no limiting is performed. In the limited version, the compact stencil property prevents direct mass exchange between nodes that are not nearest neighbors. A formal proof of sparsity is provided for simplicial and box elements. The involved element contributions can be calculated efficiently making use of matrix-free algorithms and precomputed element matrices of the reference element. Numerical studies for $\mathbb{Q}_2$ discretizations of linear and nonlinear two-dimensional test problems illustrate the virtues of monolithic convex limiting based on subcell flux decompositions.

preprint2014arXiv

Two-dimensional wave propagation in layered periodic media

We study two-dimensional wave propagation in materials whose properties vary periodically in one direction only. High order homogenization is carried out to derive a dispersive effective medium approximation. One-dimensional materials with constant impedance exhibit no effective dispersion. We show that a new kind of effective dispersion may arise in two dimensions, even in materials with constant impedance. This dispersion is a macroscopic effect of microscopic diffraction caused by spatial variation in the sound speed. We analyze this dispersive effect by using high-order homogenization to derive an anisotropic, dispersive effective medium. We generalize to two dimensions a homogenization approach that has been used previously for one-dimensional problems. Pseudospectral solutions of the effective medium equations agree to high accuracy with finite volume direct numerical simulations of the variable-coefficient equations.

preprint2013arXiv

Diffractons: solitary waves created by diffraction in periodic media

A new class of solitary waves arises in the solution of nonlinear wave equations with constant impedance and no dispersive terms. They depend on a balance between nonlinearity and a dispersion-like effect due to spatial variation in the sound speed of the medium. A high-order homogenized model confirms this effective dispersive behavior and its solutions agree well with those obtained by direct simulation of the variable-coefficient system. These waves are observed to be long-time stable, globally attracting solutions that arise in general as solutions to nonlinear wave problems with periodically-varying sound speed. They share some properties with known classes of solitary waves, but possess important differences as well.