Researcher profile

Sigal Gottlieb

Sigal Gottlieb contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 15 - UnverifiedVerification L1Unclaimed author
3works
0followers
2topics
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

3 published item(s)

preprint2021arXiv

An EIM-degradation free reduced basis method via over collocation and residual hyper reduction-based error estimation

The need for multiple interactive, real-time simulations using different parameter values has driven the design of fast numerical algorithms with certifiable accuracies. The reduced basis method (RBM) presents itself as such an option. RBM features a mathematically rigorous error estimator which drives the construction of a low-dimensional subspace. A surrogate solution is then sought in this low-dimensional space approximating the parameter-induced high fidelity solution manifold. However when the system is nonlinear or its parameter dependence nonaffine, this efficiency gain degrades tremendously, an inherent drawback of the application of the empirical interpolation method (EIM). In this paper, we augment and extend the EIM approach as a direct solver, as opposed to an assistant, for solving nonlinear partial differential equations on the reduced level. The resulting method, called Reduced Over-Collocation method (ROC), is stable and capable of avoiding the efficiency degradation. Two critical ingredients of the scheme are collocation at about twice as many locations as the number of basis elements for the reduced approximation space, and an efficient error indicator for the strategic building of the reduced solution space. The latter, the main contribution of this paper, results from an adaptive hyper reduction of the residuals for the reduced solution. Together, these two ingredients render the proposed R2-ROC scheme both offline- and online-efficient. A distinctive feature is that the efficiency degradation appearing in traditional RBM approaches that utilize EIM for nonlinear and nonaffine problems is circumvented, both in the offline and online stages. Numerical tests on different families of time-dependent and steady-state nonlinear problems demonstrate the high efficiency and accuracy of our R2-ROC and its superior stability performance.

preprint2020arXiv

Explicit and implicit error inhibiting schemes with post-processing

Efficient high order numerical methods for evolving the solution of an ordinary differential equation are widely used. The popular Runge--Kutta methods, linear multi-step methods, and more broadly general linear methods, all have a global error that is completely determined by analysis of the local truncation error. In prior work in we investigated the interplay between the local truncation error and the global error to construct {\em error inhibiting schemes} that control the accumulation of the local truncation error over time, resulting in a global error that is one order higher than expected from the local truncation error. In this work we extend our error inhibiting framework to include a broader class of time-discretization methods that allows an exact computation of the leading error term, which can then be post-processed to obtain a solution that is two orders higher than expected from truncation error analysis. We define sufficient conditions that result in a desired form of the error and describe the construction of the post-processor. A number of new explicit and implicit methods that have this property are given and tested on a variety of ordinary and partial differential equation. We show that these methods provide a solution that is two orders higher than expected from truncation error analysis alone.

preprint2020arXiv

IMEX error inhibiting schemes with post-processing

High order implicit-explicit (IMEX) methods are often desired when evolving the solution of an ordinary differential equation that has a stiff part that is linear and a non-stiff part that is nonlinear. This situation often arises in semi-discretization of partial differential equations and many such IMEX schemes have been considered in the literature. The methods considered usually have a a global error that is of the same order as the local truncation error. More recently, methods with global errors that are one order higher than predicted by the local truncation error have been devised (by Kulikov and Weiner, Ditkowski and Gottlieb). In prior work we investigated the interplay between the local truncation error and the global error to construct explicit and implicit {\em error inhibiting schemes} that control the accumulation of the local truncation error over time, resulting in a global error that is one order higher than expected from the local truncation error, and which can be post-processed to obtain a solution which is two orders higher than expected. In this work we extend our error inhibiting with post-processing framework introduced in our previous work to a class of additive general linear methods with multiple steps and stages. We provide sufficient conditions under which these methods with local truncation error of order p will produce solutions of order (p+1), which can be post-processed to order (p+2), and describe the construction of one such post-processor. We apply this approach to obtain implicit-explicit (IMEX) methods with multiple steps and stages. We present some of our new IMEX methods and show their linear stability properties, and investigate how these methods perform in practice on some numerical test cases.