Source author record

Jennifer K. Ryan

Jennifer K. Ryan 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

4works

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.

preprint2022arXiv

Enhanced Multi-Resolution Analysis for Multi-Dimensional Data Utilizing Line Filtering Techniques

In this article we introduce Line Smoothness-Increasing Accuracy-Conserving Multi-Resolution Analysis\linebreak (LSIAC-MRA). This is a procedure for exploiting convolution kernel post-processors for obtaining more accurate multi-dimensional multi-resolution analysis (MRA) in terms of error reduction. This filtering-projection tool allows for the transition of data between different resolutions while simultaneously decreasing errors in the fine grid approximation. It specifically allows for defining detail multi-wavelet coefficients when translating coarse data onto finer meshes. These coefficients are usually not defined in such cases. We show how to analytically evaluate the resulting convolutions and express the filtered approximation in a new basis. This is done by combining the filtering procedure with projection operators that allow for computational implementation of this scale transition procedure. Further, this procedure can be applied to piecewise constant approximations to functions, as it provides error reduction. We demonstrate the effectiveness of this technique in two and three dimensions.

preprint2016arXiv

Multi-dimensional filtering: Reducing the dimension through rotation

Over the past few decades there has been a strong effort towards the development of Smoothness-Increasing Accuracy-Conserving (SIAC) filters for Discontinuous Galerkin (DG) methods, designed to increase the smoothness and improve the convergence rate of the DG solution through this post-processor. These advantages can be exploited during flow visualization, for example by applying the SIAC filter to the DG data before streamline computations [Steffan {\it et al.}, IEEE-TVCG 14(3): 680-692]. However, introducing these filters in engineering applications can be challenging since a tensor product filter grows in support size as the field dimension increases, becoming computationally expensive. As an alternative, [Walfisch {\it et al.}, JOMP 38(2);164-184] proposed a univariate filter implemented along the streamline curves. Until now, this technique remained a numerical experiment. In this paper we introduce the SIAC line filter and explore how the orientation, structure and filter size affect the order of accuracy and global errors. We present theoretical error estimates showing how line filtering preserves the properties of traditional tensor product filtering, including smoothness and improvement in the convergence rate. Furthermore, numerical experiments are included, exhibiting how these filters achieve the same accuracy at significantly lower computational costs, becoming an attractive tool for the scientific visualization community.

preprint2016arXiv

Multiwavelet troubled-cell indicator for discontinuity detection of discontinuous Galerkin schemes

In this paper, we introduce a new global troubled-cell indicator for the discontinuous Galerkin (DG) method in one- and two-dimensions. This is done by taking advantage of the global expression of the DG method and re-expanding it in terms of a multiwavelet basis, which is a sum of the global average and finer details on different levels. Examining the higher level difference coefficients acts as a troubled-cell indicator, thus avoiding unnecessary increased computational cost of a new expansion. In two-dimensions the multiwavelet decomposition uses combinations of scaling functions and multiwavelets in the x- and y-directions for improved troubled-cell indication. By using such a troubled-cell indicator, we are able to reduce the computational cost by avoiding limiting in smooth regions. We present numerical examples in one- and two-dimensions and compare our troubled-cell indicator to the subcell resolution technique of Harten (1989) and the shock detector of Krivodonova, Xin, Remacle, Chevaugeon, and Flaherty (2004), which were previously investigated by Qiu and Shu (2005).

preprint2015arXiv

Automated parameters for troubled-cell indicators using outlier detection

In Vuik and Ryan (2014) we studied the use of troubled-cell indicators for discontinuity detection in nonlinear hyperbolic partial differential equations and introduced a new multiwavelet technique to detect troubled cells. We found that these methods perform well as long as a suitable, problem-dependent parameter is chosen. This parameter is used in a threshold which decides whether or not to detect an element as a troubled cell. Until now, these parameters could not be chosen automatically. The choice of the parameter has impact on the approximation: it determines the strictness of the troubled-cell indicator. An inappropriate choice of the parameter will result in detection (and limiting) of too few or too many elements. The optimal parameter is chosen such that the minimal number of troubled cells is detected and the resulting approximation is free of spurious oscillations. In this paper we will see that for each troubled-cell indicator the sudden increase or decrease of the indicator value with respect to the neighboring values is important for detection. Indication basically reduces to detecting the outliers of a vector (one dimension) or matrix (two dimensions). This is done using Tukey's boxplot approach to detect which coefficients in a vector are straying far beyond others (Tukey, 1977). We provide an algorithm that can be applied to various troubled-cell indication variables. Using this technique the problem-dependent parameter that the original indicator requires is no longer necessary as the parameter will be chosen automatically.