Source author record

David C. Del Rey Fernández

David C. Del Rey Fernández appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis

Catalog footprint

What is connected

6works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2021arXiv

Nonlinearly Stable Flux Reconstruction High-Order Methods in Split Form

The flux reconstruction (FR) method has gained popularity in the research community as it recovers promising high-order methods through modally filtered correction fields, such as the discontinuous Galerkin method, amongst others, on unstructured grids over complex geometries. Moreover, FR schemes, specifically energy stable FR (ESFR) schemes also known as Vincent-Castonguay-Jameson-Huynh schemes, have proven attractive as they allow for design flexibility as well as stability proofs for the linear advection problem on affine elements. Additionally, split forms have recently seen a resurgence in research activity due to their resultant nonlinear (entropy) stability proofs. This paper derives for the first time nonlinearly stable ESFR schemes in split form that enable nonlinear stability proofs for, uncollocated, modal, ESFR split forms with different volume and surface cubature nodes. The critical enabling technology is applying the splitting to the discrete stiffness operator. This naturally leads to appropriate surface and numerical fluxes, enabling both entropy stability and conservation proofs. When these schemes are recast in strong form, they differ from schemes found in the ESFR literature as the ESFR correction functions are incorporated on the volume integral. Furthermore, numerical experiments are conducted verifying that the new class of proposed ESFR split forms is nonlinearly stable in contrast to the standard split form ESFR approach. Lastly, the new ESFR split form is shown to obtain the correct orders of accuracy.

preprint2021arXiv

Provably Stable Flux Reconstruction High-Order Methods on Curvilinear Elements

Provably stable flux reconstruction (FR) schemes are derived for partial differential equations cast in curvilinear coordinates. Specifically, energy stable flux reconstruction (ESFR) schemes are considered as they allow for design flexibility as well as stability proofs for the linear advection problem on affine elements. Additionally, split forms are examined as they enable the development of energy stability proofs. The first critical step proves, that in curvilinear coordinates, the discontinuous Galerkin (DG) conservative and non-conservative forms are inherently different--even under exact integration and analytically exact metric terms. This analysis demonstrates that the split form is essential to developing provably stable DG schemes on curvilinear coordinates and motivates the construction of metric dependent ESFR correction functions in each element. Furthermore, the provably stable FR schemes differ from schemes in the literature that only apply the ESFR correction functions to surface terms or on the conservative form, and instead incorporate the ESFR correction functions on the full split form of the equations. It is demonstrated that the scheme is divergent when the correction functions are only used for surface reconstruction in curvilinear coordinates. We numerically verify the stability claims for our proposed FR split forms and compare them to ESFR schemes in the literature. Lastly, the newly proposed provably stable FR schemes are shown to obtain optimal orders of convergence. The scheme loses the orders of accuracy at the equivalent correction parameter value c as that of the one-dimensional ESFR scheme.

preprint2020arXiv

Mortar-based entropy-stable discontinuous Galerkin methods on non-conforming quadrilateral and hexahedral meshes

High-order entropy-stable discontinuous Galerkin (DG) methods for nonlinear conservation laws reproduce a discrete entropy inequality by combining entropy conservative finite volume fluxes with summation-by-parts (SBP) discretization matrices. In the DG context, on tensor product (quadrilateral and hexahedral) elements, SBP matrices are typically constructed by collocating at Lobatto quadrature points. Recent work has extended the construction of entropy-stable DG schemes to collocation at more accurate Gauss quadrature points. In this work, we extend entropy-stable Gauss collocation schemes to non-conforming meshes. Entropy-stable DG schemes require computing entropy conservative numerical fluxes between volume and surface quadrature nodes. On conforming tensor product meshes where volume and surface nodes are aligned, flux evaluations are required only between "lines" of nodes. However, on non-conforming meshes, volume and surface nodes are no longer aligned, resulting in a larger number of flux evaluations. We reduce this expense by introducing an entropy-stable mortar-based treatment of non-conforming interfaces via a face-local correction term, and provide necessary conditions for high-order accuracy. Numerical experiments in both two and three dimensions confirm the stability and accuracy of this approach.

preprint2016arXiv

Simultaneous Approximation Terms for Multi-Dimensional Summation-by-Parts Operators

This paper is concerned with the accurate, conservative, and stable imposition of boundary conditions and inter-element coupling for multi-dimensional summation-by-parts (SBP) finite-difference operators. More precisely, the focus is on diagonal-norm SBP operators that are not based on tensor products and are applicable to unstructured grids composed of arbitrary elements. We show how penalty terms --- simultaneous approximation terms (SATs) --- can be adapted to discretizations based on multi-dimensional SBP operators to enforce boundary and interface conditions. A general SAT framework is presented that leads to conservative and stable discretizations of the variable-coefficient advection equation. This framework includes the case where there are no nodes on the boundary of the SBP element at which to apply penalties directly. This is an important generalization, because elements analogous to Legendre-Gauss collocation, \ie without boundary nodes, typically have higher accuracy for the same number of degrees of freedom. Symmetric and upwind examples of the general SAT framework are created using a decomposition of the symmetric part of an SBP operator; these particular SATs enable the pointwise imposition of boundary and inter-element conditions. We illustrate the proposed SATs using triangular-element SBP operators with and without nodes that lie on the boundary. The accuracy, conservation, and stability properties of the resulting SBP-SAT discretizations are verified using linear advection problems with spatially varying divergence-free velocity fields.

preprint2015arXiv

Multidimensional Summation-By-Parts Operators: General Theory and Application to Simplex Elements

Summation-by-parts (SBP) finite-difference discretizations share many attractive properties with Galerkin finite-element methods (FEMs), including time stability and superconvergent functionals; however, unlike FEMs, SBP operators are not completely determined by a basis, so the potential exists to tailor SBP operators to meet different objectives. To date, application of high-order SBP discretizations to multiple dimensions has been limited to tensor product domains. This paper presents a definition for multi-dimensional SBP finite-difference operators that is a natural extension of one-dimensional SBP operators. Theoretical implications of the definition are investigated for the special case of a diagonal norm (mass) matrix. In particular, a diagonal-norm SBP operator exists on a given domain if and only if there is a cubature rule with positive weights on that domain and the polynomial-basis matrix has full rank when evaluated at the cubature nodes. Appropriate simultaneous-approximation terms are developed to impose boundary conditions weakly, and the resulting discretizations are shown to be time stable. Concrete examples of multi-dimensional SBP operators are constructed for the triangle and tetrahedron; similarities and differences with spectral-element and spectral-difference methods are discussed. An assembly process is described that builds diagonal-norm SBP operators on a global domain from element-level operators. Numerical results of linear advection on a doubly periodic domain demonstrate the accuracy and time stability of the simplex operators.

preprint2014arXiv

Generalized Summation-by-Parts Operators for the Second Derivative with Variable Coefficients

The comprehensive generalization of summation-by-parts of Del Rey Fernández et al.\ (J. Comput. Phys., 266, 2014) is extended to approximations of second derivatives with variable coefficients. This enables the construction of second-derivative operators with one or more of the following characteristics: i) non-repeating interior stencil, ii) nonuniform nodal distributions, and iii) exclusion of one or both boundary nodes. Definitions are proposed that give rise to generalized SBP operators that result in consistent, conservative, and stable discretizations of PDEs with or without mixed derivatives. It is proven that such operators can be constructed using a correction to the application of the first-derivative operator twice that is the same as used for the constant-coefficient operator. Moreover, for operators with a repeating interior stencil, a decomposition is proposed that makes the application of such operators particularly simple. A number of novel operators are constructed, including operators on pseudo-spectral nodal distributions and operators that have a repeating interior stencil, but unequal nodal spacing near boundaries. The various operators are compared to the application of the first-derivative operator twice in the context of the linear convection-diffusion equation with constant and variable coefficients.