Researcher profile

Jan Glaubitz

Jan Glaubitz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
3topics
4close 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

7 published item(s)

preprint2022arXiv

Construction and application of provable positive and exact cubature formulas

Many applications require multi-dimensional numerical integration, often in the form of a cubature formula. These cubature formulas are desired to be positive and exact for certain finite-dimensional function spaces (and weight functions). Although there are several efficient procedures to construct positive and exact cubature formulas for many standard cases, it remains a challenge to do so in a more general setting. Here, we show how the method of least squares can be used to derive provable positive and exact formulas in a general multi-dimensional setting. Thereby, the procedure only makes use of basic linear algebra operations, such as solving a least squares problem. In particular, it is proved that the resulting least squares cubature formulas are ensured to be positive and exact if a sufficiently large number of equidistributed data points is used. We also discuss the application of provable positive and exact least squares cubature formulas to construct nested stable high-order rules and positive interpolatory formulas. Finally, our findings shed new light on some existing methods for multivariate numerical integration and under which restrictions these are ensured to be successful.

preprint2022arXiv

Energy-stable global radial basis function methods on summation-by-parts form

Radial basis function methods are powerful tools in numerical analysis and have demonstrated good properties in many different simulations. However, for time-dependent partial differential equations, only a few stability results are known. In particular, if boundary conditions are included, stability issues frequently occur. The question we address in this paper is how provable stability for RBF methods can be obtained. We develop and construct energy-stable radial basis function methods using the general framework of summation-by-parts operators often used in the Finite Difference and Finite Element communities.

preprint2022arXiv

Sequential image recovery from noisy and under-sampled Fourier data

A new algorithm is developed to jointly recover a temporal sequence of images from noisy and under-sampled Fourier data. Specifically, we consider the case where each data set is missing vital information that prevents its (individual) accurate recovery. Our new method is designed to restore the missing information in each individual image by "borrowing" it from the other images in the sequence. As a result, {\em all} of the individual reconstructions yield improved accuracy. The use of high resolution Fourier edge detection methods is essential to our algorithm. In particular, edge information is obtained directly from the Fourier data which leads to an accurate coupling term between data sets. Moreover, data loss is largely avoided as coarse reconstructions are not required to process inter- and intra-image information. Numerical examples are provided to demonstrate the accuracy, efficiency and robustness of our new method.

preprint2021arXiv

Towards Stable Radial Basis Function Methods for Linear Advection Problems

In this work, we investigate (energy) stability of global radial basis function (RBF) methods for linear advection problems. Classically, boundary conditions (BC) are enforced strongly in RBF methods. By now it is well-known that this can lead to stability problems, however. Here, we follow a different path and propose two novel RBF approaches which are based on a weak enforcement of BCs. By using the concept of flux reconstruction and simultaneous approximation terms (SATs), respectively, we are able to prove that both new RBF schemes are strongly (energy) stable. Numerical results in one and two spatial dimensions for both scalar equations and systems are presented, supporting our theoretical analysis.

preprint2020arXiv

Stable discretisations of high-order discontinuous Galerkin methods on equidistant and scattered points

In this work, we propose and investigate stable high-order collocation-type discretisations of the discontinuous Galerkin method on equidistant and scattered collocation points. We do so by incorporating the concept of discrete least squares into the discontinuous Galerkin framework. Discrete least squares approximations allow us to construct stable and high-order accurate approximations on arbitrary collocation points, while discrete least squares quadrature rules allow us their stable and exact numerical integration. Both methods are computed efficiently by using bases of discrete orthogonal polynomials. Thus, the proposed discretisation generalises known classes of discretisations of the discontinuous Galerkin method, such as the discontinuous Galerkin collocation spectral element method. We are able to prove conservation and linear $L^2$-stability of the proposed discretisations. Finally, numerical tests investigate their accuracy and demonstrate their extension to nonlinear conservation laws, systems, longtime simulations, and a variable coefficient problem in two space dimensions.

preprint2020arXiv

Stable High Order Quadrature Rules for Scattered Data and General Weight Functions

Numerical integration is encountered in all fields of numerical analysis and the engineering sciences. By now, various efficient and accurate quadrature rules are known; for instance, Gauss-type quadrature rules. In many applications, however, it might be impractical---if not even impossible---to obtain data to fit known quadrature rules. Often, experimental measurements are performed at equidistant or even scattered points in space or time. In this work, we propose stable high order quadrature rules for experimental data, which can accurately handle general weight functions.

preprint2019arXiv

Analysis of Artificial Dissipation of Explicit and Implicit Time-Integration Methods

Stability is an important aspect of numerical methods for hyperbolic conservation laws and has received much interest. However, continuity in time is often assumed and only semidiscrete stability is studied. Thus, it is interesting to investigate the influence of explicit and implicit time integration methods on the stability of numerical schemes. If an explicit time integration method is applied, spacially stable numerical schemes for hyperbolic conservation laws can result in unstable fully discrete schemes. Focusing on the explicit Euler method (and convex combinations thereof), undesired terms in the energy balance trigger this phenomenon and introduce an erroneous growth of the energy over time. In this work, we study the influence of artificial dissipation and modal filtering in the context of discontinuous spectral element methods to remedy these issues. In particular, lower bounds on the strength of both artificial dissipation and modal filtering operators are given and an adaptive procedure to conserve the (discrete) $\mathbf{L}_2$ norm of the numerical solution in time is derived. This might be beneficial in regions where the solution is smooth and for long time simulations. Moreover, this approach is used to study the connections between explicit and implicit time integration methods and the associated energy production. By adjusting the adaptive procedure, we demonstrate that filtering in explicit time integration methods is able to mimic the dissipative behavior inherent in implicit time integration methods. This contribution leads to a better understanding of existing algorithms and numerical techniques, in particular the application of artificial dissipation as well as modal filtering in the context of numerical methods for hyperbolic conservation laws together with the selection of explicit or implicit time integration methods.