Source author record

John R. Singler

John R. Singler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

3works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

3 published item(s)

preprint2022arXiv

A New Global Divergence Free and Pressure-Robust HDG Method for Tangential Boundary Control of Stokes Equations

In [ESAIM: M2AN, 54(2020), 2229-2264], we proposed an HDG method to approximate the solution of a tangential boundary control problem for the Stokes equations and obtained an optimal convergence rate for the optimal control {that reflects its global regularity}. However, the error estimates depend on the pressure, and the velocity is not divergence free. The importance of pressure-robust numerical methods for fluids was addressed by John et al. [SIAM Review, 59(2017), 492-544]. In this work, we devise a new HDG method to approximate the solution of the Stokes tangential boundary control problem; the HDG method is also of independent interest for solving the Stokes equations. This scheme yields a $\mathbf{H}(\mathrm{div})$ conforming, globally divergence free, and pressure-robust solution. To the best of our knowledge, this is the first time such a numerical scheme has been obtained for an optimal boundary control problem for the Stokes equations. We also provide numerical experiments to show the performance of the new HDG method and the advantage over the non pressure-robust scheme.

preprint2018arXiv

A note on incremental POD algorithms for continuous time data

In our earlier work [Fareed et al., Comput. Math. Appl. 75 (2018), no. 6, 1942-1960], we developed an incremental approach to compute the proper orthogonal decomposition (POD) of PDE simulation data. Specifically, we developed an incremental algorithm for the SVD with respect to a weighted inner product for the discrete time POD computations. For continuous time data, we used an approximate approach to arrive at a discrete time POD problem and then applied the incremental SVD algorithm. In this note, we analyze the continuous time case with simulation data that is piecewise constant in time such that each data snapshot is expanded in a finite collection of basis elements of a Hilbert space. We first show that the POD is determined by the SVD of two different data matrices with respect to weighted inner products. Next, we develop incremental algorithms for approximating the two matrix SVDs with respect to the different weighted inner products. Finally, we show neither approximate SVD is more accurate than the other; specifically, we show the incremental algorithms return equivalent results.

preprint2018arXiv

Error Analysis of an Incremental POD Algorithm for PDE Simulation Data

In our earlier work [Fareed et al., Comput. Math. Appl. 75 (2018), no. 6, 1942-1960], we proposed an incremental SVD algorithm with respect to a weighted inner product to compute the proper orthogonal decomposition (POD) of a set of simulation data for a partial differential equation (PDE) without storing the data. In this work, we perform an error analysis of the incremental SVD algorithm. We also modify the algorithm to incrementally update both the SVD and an error bound when a new column of data is added. We show the algorithm produces the exact SVD of an approximate data matrix, and the operator norm error between the approximate and exact data matrices is bounded above by the computed error bound. This error bound also allows us to bound the error in the incrementally computed singular values and singular vectors. We illustrate our analysis with numerical results for three simulation data sets from a 1D FitzHugh-Nagumo PDE system with various choices of the algorithm truncation tolerances.