Researcher profile

Manuel Quezada de Luna

Manuel Quezada de Luna contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
6topics
2close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 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.