Source author record

Maciej Balajewicz

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

ResearcherUnclaimed source record

physics.flu-dyn math.NA physics.comp-ph Numerical Analysis Computational Engineering, Finance, and Science Machine Learning math.DS q-fin.CP

Catalog footprint

What is connected

7works

8topics

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.

preprint2020arXiv

Physics-aware registration based auto-encoder for convection dominated PDEs

We design a physics-aware auto-encoder to specifically reduce the dimensionality of solutions arising from convection-dominated nonlinear physical systems. Although existing nonlinear manifold learning methods seem to be compelling tools to reduce the dimensionality of data characterized by a large Kolmogorov n-width, they typically lack a straightforward mapping from the latent space to the high-dimensional physical space. Moreover, the realized latent variables are often hard to interpret. Therefore, many of these methods are often dismissed in the reduced order modeling of dynamical systems governed by the partial differential equations (PDEs). Accordingly, we propose an auto-encoder type nonlinear dimensionality reduction algorithm. The unsupervised learning problem trains a diffeomorphic spatio-temporal grid, that registers the output sequence of the PDEs on a non-uniform parameter/time-varying grid, such that the Kolmogorov n-width of the mapped data on the learned grid is minimized. We demonstrate the efficacy and interpretability of our approach to separate convection/advection from diffusion/scaling on various manufactured and physical systems.

preprint2016arXiv

Reduced Order Models for Pricing European and American Options under Stochastic Volatility and Jump-Diffusion Models

European options can be priced by solving parabolic partial(-integro) differential equations under stochastic volatility and jump-diffusion models like Heston, Merton, and Bates models. American option prices can be obtained by solving linear complementary problems (LCPs) with the same operators. A finite difference discretization leads to a so-called full order model (FOM). Reduced order models (ROMs) are derived employing proper orthogonal decomposition (POD). The early exercise constraint of American options is enforced by a penalty on subset of grid points. The presented numerical experiments demonstrate that pricing with ROMs can be orders of magnitude faster within a given model parameter variation range.

preprint2015arXiv

Minimal subspace rotation on the Stiefel manifold for stabilization and enhancement of projection-based reduced order models for the compressible Navier-Stokes equations

For a projection-based reduced order model (ROM) of a fluid flow to be stable and accurate, the dynamics of the truncated subspace must be taken into account. This paper proposes an approach for stabilizing and enhancing projection-based fluid ROMs in which truncated modes are accounted for a priori via a minimal rotation of the projection subspace. Attention is focused on the full non-linear compressible Navier-Stokes equations in specific volume form as a step toward a more general formulation for problems with generic non-linearities. Unlike traditional approaches, no empirical turbulence modeling terms are required, and consistency between the ROM and the full order model from which the ROM is derived is maintained. Mathematically, the approach is formulated as a trace minimization problem on the Stiefel manifold. The reproductive as well as predictive capabilities of the method are evaluated on several compressible flow problems, including a problem involving laminar flow over an airfoil with a high angle of attack, and a channel-driven cavity flow problem.

preprint2015arXiv

Projection-based model reduction for contact problems

To be feasible for computationally intensive applications such as parametric studies, optimization and control design, large-scale finite element analysis requires model order reduction. This is particularly true in nonlinear settings that tend to dramatically increase computational complexity. Although significant progress has been achieved in the development of computational approaches for the reduction of nonlinear computational mechanics models, addressing the issue of contact remains a major hurdle. To this effect, this paper introduces a projection-based model reduction approach for both static and dynamic contact problems. It features the application of a non-negative matrix factorization scheme to the construction of a positive reduced-order basis for the contact forces, and a greedy sampling algorithm coupled with an error indicator for achieving robustness with respect to model parameter variations. The proposed approach is successfully demonstrated for the reduction of several two-dimensional, simple, but representative contact and self contact computational models.

preprint2014arXiv

Reduction of Nonlinear Embedded Boundary Models for Problems with Evolving Interfaces

Embedded boundary methods alleviate many computational challenges, including those associated with meshing complex geometries and solving problems with evolving domains and interfaces. Developing model reduction methods for computational frameworks based on such methods seems however to be challenging. Indeed, most popular model reduction techniques are projection-based, and rely on basis functions obtained from the compression of simulation snapshots. In a traditional interface-fitted computational framework, the computation of such basis functions is straightforward, primarily because the computational domain does not contain in this case a fictitious region. This is not the case however for an embedded computational framework because the computational domain typically contains in this case both real and ghost regions whose definitions complicate the collection and compression of simulation snapshots. The problem is exacerbated when the interface separating both regions evolves in time. This paper addresses this issue by formulating the snapshot compression problem as a weighted low-rank approximation problem where the binary weighting identifies the evolving component of the individual simulation snapshots. The proposed approach is application independent and therefore comprehensive. It is successfully demonstrated for the model reduction of several two-dimensional, vortex-dominated, fluid-structure interaction problems.

preprint2013arXiv

Lyapunov stable Galerkin models of post-transient incompressible flows

A method for deriving provably stable low-dimensional Galerkin models of post-transient incompressible flows is introduced. The proposed approach involves an iterative procedure for expansion modes that satisfy Lyapunov stability in the neighborhood of a fixed point. The approach is demonstrated using two prototypical flow configurations: a two-dimensional mixing layer, and two-dimensional flow inside a lid-driven cavity. This new methodology can be a building block in an effort to develop more accurate and more robust low-dimensional models of incompressible flows.

preprint2012arXiv

A Novel Model Order Reduction Approach for Navier-Stokes Equations at High Reynolds Number

A new approach to model order reduction of the Navier-Stokes equations at high Reynolds number is proposed. Unlike traditional approaches, this method does not rely on empirical turbulence modeling or modification of the Navier-Stokes equations. It provides spatial basis functions different from the usual proper orthogonal decomposition basis function in that, in addition to optimally representing the training data set, the new basis functions also provide stable and accurate reduced-order models. The proposed approach is illustrated with two test cases: two-dimensional flow inside a square lid-driven cavity and a two-dimensional mixing layer.