42.9NAJun 3
Parametric Reduced Order Models for the Generalized Kuramoto--Sivashinsky EquationsMd Rezwan Bin Mizan, Maxim Olshanskii, Ilya Timofeyev
The paper studies parametric Reduced Order Models (ROMs) for the Kuramoto--Sivashinsky (KS) and generalized Kuramoto--Sivashinsky (gKS) equations. We consider several POD and POD-DEIM projection ROMs with various strategies for parameter sampling and snapshot collection. The aim is to identify an approach for constructing a ROM that is efficient across a range of parameters, encompassing several regimes exhibited by the KS and gKS solutions: weakly chaotic, transitional, and quasi-periodic dynamics. We describe such an approach and demonstrate that it is essential to develop ROMs that adequately represent the short-time transient behavior of the gKS model.
NAApr 20, 2017
Inf-sup stability of geometrically unfitted Stokes finite elementsJohnny Guzmán, Maxim Olshanskii
The paper shows an inf-sup stability property for several well-known 2D and 3D Stokes elements on triangulations which are not fitted to a given smooth or polygonal domain. The property implies stability and optimal error estimates for a class of unfitted finite element methods for the Stokes and Stokes interface problems, such as Nitsche-XFEM or cutFEM. The error analysis is presented for the Stokes problem. All assumptions made in the paper are satisfied once the background mesh is shape-regular and fine enough.
NAJan 22, 2016
Numerical integration over implicitly defined domains for higher order unfitted finite element methodsMaxim Olshanskii, Danil Safin
The paper studies several approaches to numerical integration over a domain defined implicitly by an indicator function such as the level set function. The integration methods are based on subdivision, moment--fitting, local quasi-parametrization and Monte-Carlo techniques. As an application of these techniques, the paper addresses numerical solution of elliptic PDEs posed on domains and manifolds defined implicitly. A higher order unfitted finite element method (FEM) is assumed for the discretization. In such a method the underlying mesh is not fitted to the geometry, and hence the errors of numerical integration over curvilinear elements affect the accuracy of the finite element solution together with approximation errors. The paper studies the numerical complexity of the integration procedures and the performance of unfitted FEMs which employ these tools.
NAJan 6, 2019
A computational study of lateral phase separation in biological membranesVladimir Yushutin, Annalisa Quaini, Sheereen Majd et al.
Conservative and non-conservative phase-field models are considered for the numerical simulation of lateral phase separation and coarsening in biological membranes. An unfitted finite element method is devised for these models to allow for a flexible treatment of complex shapes in the absence of an explicit surface parametrization. For a set of biologically relevant shapes and parameter values, the paper compares the dynamic coarsening produced by conservative and non-conservative numerical models, its dependence on certain geometric characteristics and convergence to the final equilibrium
4.5FLU-DYNMay 7
Reduced-Order Modeling of Parameterized Visco-Plastic Shallow FlowsMd Rezwan Bin Mizan, Ilya Timofeyev, Maxim Olshanskii
We propose a non-intrusive reduced-order modeling framework for parametrized visco-plastic free-surface flows governed by a shallow-water formulation of Herschel--Bulkley fluids. These flows exhibit strong nonlinearities, non-smooth rheology, moving fronts, and yield surfaces, making efficient surrogate modeling particularly challenging. To address this challenge, we employ a tensor-based approach in which the solution manifold is approximated using a low-rank representation obtained via higher-order singular value decomposition of snapshot data over a structured parameter space. The resulting tensorial reduced-order model (TROM) enables rapid online evaluation by directly reconstructing solution trajectories from the compressed representation, thereby avoiding the need to perform time integration of a reduced dynamical system. The proposed non-intrusive framework can be interpreted as an encoder--decoder architecture with a compressed latent representation and efficient multilinear decoding. Numerical experiments demonstrate that the proposed approach accurately captures key flow features, including front propagation, plug and shear regions, and near-stopping dynamics, while achieving substantial computational speedups relative to full-order simulations.
29.5NAMar 14
Tensorial Reduced-Order Models for Parametric Coupled Reaction-Diffusion Systems: Application to Brain Tumor Growth ModelingAsikul Islam, Md Rezwan Bin Mizan, Maxim Olshanskii et al.
We construct efficient surrogate models for parametric forward operators arising in brain tumor growth simulations, governed by coupled semilinear parabolic reaction-diffusion systems on heterogeneous two- and three-dimensional domains. We consider two models of increasing complexity: a scalar single-species formulation and a six-state, nine-parameter multi-species go-or-grow model. The governing equations are discretized using a finite volume method and integrated in time via an operator-splitting strategy. We develop tensorial reduced-order model (TROM) surrogates based on the Higher-Order Singular Value Decomposition in Tucker format and the Tensor Train decomposition, each in intrusive and non-intrusive variants. The models are compared against a classical proper orthogonal decomposition (POD) ROM baseline. Numerical experiments with up to $m=9$ model parameters demonstrate speedups of $85\times$-$120\times$ relative to the full-order solver while maintaining excellent accuracy, establishing tensorial surrogates as a rigorous and efficient computational foundation for many-query workflows.