Design Principle for Mode-Consistent Galerkin Closure under a Physical Energy Metric for Hyperbolic Systems

Jul 16, 2026

Derives a design principle for Galerkin approximations of energy-conserving hyperbolic systems using Arakawa’s structure-preserving philosophy. The goal is to construct, within a resolved finite-mode space, the same modal energy-exchange structure as the continuous model so total energy conservation follows naturally. Introduces a state-dependent physical-energy metric H(U) and an energy-compatibility identity. A Galerkin projection with H-metric summation-by-parts and an energy-compatibility closure recover the desired modal exchange, including interface flux closure and an O(hp+1)-consistent defect analysis.

LanguageEnglish

CategoryResearch & Report

FormatPDF

Pages29

Reading Time1 h 13 min

Introduction

Structural preservation for resolved dynamics

Model hierarchy: continuous, finite-mode exact integration, and practical finite-quadrature

arXiv :2607 .05781v1 [math .NA] 7 Jul 2026 Design Principle for Mode-Consistent Galerkin Closure under a Physical Energy Metric for Hyperbolic Systems Hirofumi Tomitaa,∗ aRIKEN Center for Computational Science, Kobe, Japan Abstract This paper derives a design principle for Galerkin approximations of energy-conserving hyperbolic systems, based on Arakawa’s philosophy of structure preservation. The aim is to construct, within a resolved finite-mode space, the same modal-energy-exchange structure as in the continuous system, so that total energy conservation follows as its consequence. We introduce a state-dependent metric H( U) representing the physical energy density and derive the corresponding energy-compatibility identity. In the exact-integration infinite-mode reference model, H-orthogonalization makes the volume operator antisymmetric, so the volume contribution to the modal energy balance is expressed as pairwise exchange between modes. The boundary contribution is also represented as exchange with modes of adjacent elements, and internal interface exchanges cancel pairwise. To reproduce this structure in a semi-discrete finite-mode system, we combine two constructions. First, a Galerkin projection coupled with the physical energy metric guarantees the H-metric summation-by-parts identity. Second, an energy-compatibility closure cancels the compatibility residual by modifying the evolution of the H-metric mass matrix. As a result, the modal-energy-exchange structure is recovered in the finite-mode system. For discontinuous element-boundary traces, the interface contribution is closed by a shared numerical energy flux satisfying the same pairwise energy balance. We then compare the practical operator construction with the exact-integration finite-mode reference model. The defect in the antisymmetric modalenergy-exchange operator is decomposed into fixed-quadrature and projection-quadrature contributions, yielding an O(hp+1)-consistent estimate. Finally, transforming the equation back to the original Galerkin basis gives an equivalent fixed-basis coefficient equation that is directly implementable. Keywords: symmetric hyperbolic systems, Galerkin methods, structure-preserving discretization, energy-compatible closure, modal energy exchange 2020 MSC: 65M60, 65M12, 35L40 1. Introduction When a partial differential equation (PDE) is discretized, its dynamics are separated into resolved and unresolved components. This raises a fundamental question: what should be regarded as the correct resolved dynamics? The answer depends on which characteristics of the PDE are inherited by the resolved dynamics. In many cases, the discrete system is considered an approximation of the original PDE. This paper takes the position that, from an energy-structure perspective, the resolved dynamics should also reproduce the energy structure of the original PDE. To achieve this, it is necessary first to clarify which energy is to be represented. Then, by introducing a metric associated with the selected energy, we can describe how the energy is exchanged between modes. This paper treats physical energy as its target and focuses particularly on its modal dynamics. In this case, the resolved dynamics should retain the same modal-energy-exchange structure as that of the continuous system; that is, spurious energy transfer, artificial dissipation, or non-physical interaction pathways should not be introduced. This approach is based on Arakawa’s philosophy of structural preservation [1, 2] . Such structural preservation is also essential for the closure ∗ Corresponding author. Email address: [htomita@riken.jp](htomita@riken.jp) (Hirofumi Tomita) problem of unresolved dynamics, because unresolved effects should be modeled on the foundation that the resolved dynamics have the correct energy structure. Otherwise, it would be impossible to consistently discuss how unresolved effects should act on the discrete system. To represent

FAQs

What is the main objective of the paper for Galerkin discretization of hyperbolic systems?

To design a Galerkin closure so the resolved finite-mode dynamics reproduces the continuous system’s modal energy-exchange structure, preventing spurious energy transfer and non-physical interactions.

How does the paper define the energy used for the metric-based construction?

It treats physical energy as the target and introduces a state-dependent metric H(U) that represents physical energy density, enabling energy-compatibility relations for modal dynamics.

What mechanism restores the correct modal-energy-exchange structure in the semi-discrete finite-mode system?

A combination of an H-metric summation-by-parts-guaranteeing Galerkin projection and an energy-compatibility closure that modifies the evolution of the H-metric mass matrix to cancel the compatibility residual.

Jul 16, 2026

Upload by Theodore