An optimization-based construction procedure for function space based summation-by-parts operators on arbitrary grids

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Glaubitz, Jan, Nordström, Jan, Öffner, Philipp
Format: Preprint
Publié: 2024
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866914796177195008
author Glaubitz, Jan
Nordström, Jan
Öffner, Philipp
author_facet Glaubitz, Jan
Nordström, Jan
Öffner, Philipp
contents We introduce a novel construction procedure for one-dimensional summation-by-parts (SBP) operators. Existing construction procedures for FSBP operators of the form $D = P^{-1} Q$ proceed as follows: Given a boundary operator $B$, the norm matrix $P$ is first determined and then in a second step the complementary matrix $Q$ is calculated to finally get the FSBP operator $D$. In contrast, the approach proposed here determines the norm and complementary matrices, $P$ and $Q$, simultaneously by solving an optimization problem. The proposed construction procedure applies to classical SBP operators based on polynomial approximation and the broader class of function space SBP (FSBP) operators. According to our experiments, the presented approach yields a numerically stable construction procedure and FSBP operators with higher accuracy for diagonal norm difference operators at the boundaries than the traditional approach. Through numerical simulations, we highlight the advantages of our proposed technique.
format Preprint
id arxiv_https___arxiv_org_abs_2405_08770
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle An optimization-based construction procedure for function space based summation-by-parts operators on arbitrary grids
Glaubitz, Jan
Nordström, Jan
Öffner, Philipp
Numerical Analysis
Optimization and Control
65M12, 65M60, 65M70, 65D25
We introduce a novel construction procedure for one-dimensional summation-by-parts (SBP) operators. Existing construction procedures for FSBP operators of the form $D = P^{-1} Q$ proceed as follows: Given a boundary operator $B$, the norm matrix $P$ is first determined and then in a second step the complementary matrix $Q$ is calculated to finally get the FSBP operator $D$. In contrast, the approach proposed here determines the norm and complementary matrices, $P$ and $Q$, simultaneously by solving an optimization problem. The proposed construction procedure applies to classical SBP operators based on polynomial approximation and the broader class of function space SBP (FSBP) operators. According to our experiments, the presented approach yields a numerically stable construction procedure and FSBP operators with higher accuracy for diagonal norm difference operators at the boundaries than the traditional approach. Through numerical simulations, we highlight the advantages of our proposed technique.
title An optimization-based construction procedure for function space based summation-by-parts operators on arbitrary grids
topic Numerical Analysis
Optimization and Control
65M12, 65M60, 65M70, 65D25
url https://arxiv.org/abs/2405.08770