Source author record

Neeraj Sarna

Neeraj Sarna 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

3works

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.

preprint2021arXiv

Hyper-reduction for parametrized transport dominated problems via online-adaptive reduced meshes

We propose an efficient residual minimization technique for the nonlinear model-order reduction of parameterized hyperbolic partial differential equations. Our nonlinear approximation space is a span of snapshots evaluated on a shifted spatial domain, and we compute our reduced approximation via residual minimization. To speed-up the residual minimization, we compute and minimize the residual on a (preferably small) subset of the mesh, the so-called reduced mesh. Due to the nonlinearity of our approximation space we show that, similar to the solution, the residual also exhibits transport-type behaviour. To account for this behaviour, we introduce online-adaptivity in the reduced mesh by "moving" it along the spatial domain with parameter dependent shifts. We also present an extension of our method to spatial transforms different from shifting. Numerical experiments showcase the effectiveness of our method and the inaccuracies resulting from a non-adaptive reduced mesh.

preprint2021arXiv

Stable Boundary Conditions and Discretization for PN Equations

A solution to the linear Boltzmann equation satisfies an energy bound, which reflects a natural fact: The energy of particles in a finite volume is bounded in time by the energy of particles initially occupying the volume augmented by the energy transported into the volume by particles entering the volume over time. In this paper, we present boundary conditions (BCs) for the spherical harmonic (PN) approximation, which ensure that this fundamental energy bound is satisfied by the PN approximation. Our BCs are compatible with the characteristic waves of PN equations and determine the incoming waves uniquely. Both, energy bound and compatibility, are shown based on abstract formulations of PN equations and BCs to isolate the necessary structures and properties. The BCs are derived from a Marshak type formulation of BC and base on a non-classical even/odd-classification of spherical harmonic functions and a stabilization step, which is similar to the truncation of the series expansion in the PN method. We show that summation by parts (SBP) finite differences on staggered grids in space and the method of simultaneous approximation terms (SAT) allows to maintain the energy bound also on the semi-discrete level.

preprint2020arXiv

Data-Driven Snapshot Calibration via Monotonic Feature Matching

Snapshot matrices of hyperbolic equations have a slow singular value decay, resulting in inefficient reduced-order models. We develop on the idea of inducing a faster singular value decay by computing snapshots on a transformed spatial domain, or the so-called snapshot calibration/transformation. We are particularly interested in problems involving shock collision, shock rarefaction-fan collision, shock formation, etc. For such problems, we propose a realizable algorithm to compute the spatial transform using monotonic feature matching. We consider discontinuities and kinks as features, and by carefully partitioning the parameter domain, we ensure that the spatial transform has properties that are desirable both from a theoretical and an implementation standpoint. We use these properties to prove that our method results in a fast m-width decay of a so-called calibrated manifold. A crucial observation we make is that due to calibration, the m-width does not only depend on m but also on the accuracy of the full order model, which is in contrast to elliptic and parabolic problems that do not need calibration. The method we propose only requires the solution snapshots and not the underlying partial differential equation (PDE) and is therefore, data-driven. We perform several numerical experiments to demonstrate the effectiveness of our method.