Source author record

Maxim A. Olshanskii

Maxim A. Olshanskii 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

15works
7topics
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

15 published item(s)

preprint2022arXiv

Interpolatory tensorial reduced order models for parametric dynamical systems

The paper introduces a reduced order model (ROM) for numerical integration of a dynamical system which depends on multiple parameters. The ROM is a projection of the dynamical system on a low dimensional space that is both problem-dependent and parameter-specific. The ROM exploits compressed tensor formats to find a low rank representation for a sample of high-fidelity snapshots of the system state. This tensorial representation provides ROM with an orthogonal basis in a universal space of all snapshots and encodes information about the state variation in parameter domain. During the online phase and for any incoming parameter, this information is used to find a reduced basis that spans a parameter-specific subspace in the universal space. The computational cost of the online phase then depends only on tensor compression ranks, but not on space or time resolution of high-fidelity computations. Moreover, certain compressed tensor formats enable to avoid the adverse effect of parameter space dimension on the online costs (known as the curse of dimension). The analysis of the approach includes an estimate for the representation power of the acquired ROM basis. We illustrate the performance and prediction properties of the ROM with several numerical experiments, where tensorial ROM's complexity and accuracy is compared to those of conventional POD-ROM.

preprint2020arXiv

A stable method for 4D CT-based CFD simulation in the right ventricle of a TGA patient

The paper discusses a stabilization of a finite element method for the equations of fluid motion in a time-dependent domain. After experimental convergence analysis, the method is applied to simulate a blood flow in the right ventricle of a post-surgery patient with the transposition of the great arteries disorder. The flow domain is reconstructed from a sequence of 4D CT images. The corresponding segmentation and triangulation algorithms are also addressed in brief.

preprint2020arXiv

Error analysis of higher order trace finite element methods for the surface Stokes equations

The paper studies a higher order unfitted finite element method for the Stokes system posed on a surface in three-dimensional space. The method employs generalized Taylor-Hood finite element pairs on tetrahedral bulk mesh to discretize the Stokes system on embedded surface. Stability and optimal order convergence results are proved. The proofs include a complete quantification of geometric errors stemming from approximate parametric representation of the surface. Numerical experiments include formal convergence studies and an example of the Kelvin-Helmholtz instability problem on the unit sphere.

preprint2020arXiv

Inf-sup stability of the trace P2-P1 Taylor-Hood elements for surface PDEs

The paper studies a geometrically unfitted finite element method (FEM), known as trace FEM or cut FEM, for the numerical solution of the Stokes system posed on a closed smooth surface. A trace FEM based on standard Taylor-Hood (continuous P2-P1) bulk elements is proposed. A so-called volume normal derivative stabilization, known from the literature on trace FEM, is an essential ingredient of this method. The key result proved in the paper is an inf-sup stability of the trace P2-P1 finite element pair, with the stability constant uniformly bounded with respect to the discretization parameter and the position of the surface in the bulk mesh. Optimal order convergence of a consistent variant of the finite element method follows from this new stability result and interpolation properties of the trace FEM. Properties of the method are illustrated with numerical examples.

preprint2016arXiv

Finite stopping times for freely oscillating drop of a yield stress fluid

The paper addresses the question if there exists a finite stopping time for an unforced motion of a yield stress fluid with free surface. A variation inequality formulation is deduced for the problem of yield stress fluid dynamics with a free surface. Free surface is assumed to evolve with a normal velocity the flow. We also consider capillary forces acting along the free surface. Based on the variational inequality formulation an energy equality is obtained, where kinetic and free energy rate of change is in a balance with the internal energy viscoplastic dissipation and the work of external forces. Further, the paper considers free small-amplitude oscillations of a droplet of Herschel-Bulkley fluid under the action of surface tension forces. Under certain assumptions it is shown that the finite stopping time $T_f$ of oscillations exists once the yield stress parameter is positive and the flow index $α$ satisfies ($α\ge1$). Results of several numerical experiments illustrate the analysis, reveal the dependence of $T_f$ on problem parameters and suggest an instantaneous transition of the whole drop from yielding state to the rigid one.

preprint2015arXiv

A narrow-band unfitted finite element method for elliptic PDEs posed on surfaces

The paper studies a method for solving elliptic partial differential equations posed on hypersurfaces in $\mathbb{R}^N$, $N=2,3$. The method allows a surface to be given implicitly as a zero level of a level set function. A surface equation is extended to a narrow-band neighborhood of the surface. The resulting extended equation is a non-degenerate PDE and it is solved on a bulk mesh that is unaligned to the surface. An unfitted finite element method is used to discretize extended equations. Error estimates are proved for finite element solutions in the bulk domain and restricted to the surface. The analysis admits finite elements of a higher order and gives sufficient conditions for archiving the optimal convergence order in the energy norm. Several numerical examples illustrate the properties of the method.

preprint2015arXiv

Unconditional long-time stability of a velocity-vorticity method for the 2D Navier-Stokes equations

We prove unconditional long-time stability for a particular velocity-vorticity discretization of the 2D Navier-Stokes equations. The scheme begins with a formulation that uses the Lamb vector to couple the usual velocity-pressure system to the vorticity dynamics equation, and then discretizes with the finite element method in space and implicit-explicit BDF2 in time, with the vorticity equation decoupling at each time step. We prove the method's vorticity and velocity are both long-time stable in the $L^2$ and $H^1$ norms, without any timestep restriction. Moreover, our analysis avoids the use of Gronwall-type estimates, which leads us to stability bounds with only polynomial (instead of exponential) dependence on the Reynolds number. Numerical experiments are given that demonstrate the effectiveness of the method.

preprint2014arXiv

A trace finite element method for a class of coupled bulk-interface transport problems

In this paper we study a system of advection-diffusion equations in a bulk domain coupled to an advection-diffusion equation on an embedded surface. Such systems of coupled partial differential equations arise in, for example, the modeling of transport and diffusion of surfactants in two-phase flows. The model considered here accounts for adsorption-desorption of the surfactants at a sharp interface between two fluids and their transport and diffusion in both fluid phases and along the interface. The paper gives a well-posedness analysis for the system of bulk-surface equations and introduces a finite element method for its numerical solution. The finite element method is unfitted, i.e., the mesh is not aligned to the interface. The method is based on taking traces of a standard finite element space both on the bulk domains and the embedded surface. The numerical approach allows an implicit definition of the surface as the zero level of a level-set function. Optimal order error estimates are proved for the finite element method both in the bulk-surface energy norm and the $L^2$-norm. The analysis is not restricted to linear finite elements and a piecewise planar reconstruction of the surface, but also covers the discretization with higher order elements and a higher order surface reconstruction.

preprint2014arXiv

An Eulerian space-time finite element method for diffusion problems on evolving surfaces

In this paper, we study numerical methods for the solution of partial differential equations on evolving surfaces. The evolving hypersurface in $\Bbb{R}^d$ defines a $d$-dimensional space-time manifold in the space-time continuum $\Bbb{R}^{d+1}$. We derive and analyze a variational formulation for a class of diffusion problems on the space-time manifold. For this variational formulation new well-posedness and stability results are derived. The analysis is based on an inf-sup condition and involves some natural, but non-standard, (anisotropic) function spaces. Based on this formulation a discrete in time variational formulation is introduced that is very suitable as a starting point for a discontinuous Galerkin (DG) space-time finite element discretization. This DG space-time method is explained and results of numerical experiments are presented that illustrate its properties.

preprint2014arXiv

Error analysis of a space-time finite element method for solving PDEs on evolving surfaces

In this paper we present an error analysis of an Eulerian finite element method for solving parabolic partial differential equations posed on evolving hypersurfaces in $\mathbb{R}^d$, $d=2,3$. The method employs discontinuous piecewise linear in time -- continuous piecewise linear in space finite elements and is based on a space-time weak formulation of a surface PDE problem. Trial and test surface finite element spaces consist of traces of standard volumetric elements on a space-time manifold resulting from the evolution of a surface. We prove first order convergence in space and time of the method in an energy norm and second order convergence in a weaker norm. Furthermore, we derive regularity results for solutions of parabolic PDEs on an evolving surface, which we need in a duality argument used in the proof of the second order convergence estimate.

preprint2013arXiv

A connection between filter stabilization and eddy viscosity models

Recently, a new approach for the stabilization of the incompressible Navier-Stokes equations for higher Reynolds numbers was introduced based on the nonlinear differential filtering of solutions on every time step of a discrete scheme. In this paper, the stabilization is shown to be equivalent to a certain eddy-viscosity model in LES. This allows a refined analysis and further understanding of desired filter properties. We also consider the application of the filtering in a projection (pressure correction) method, the standard splitting algorithm for time integration of the incompressible fluid equations. The paper proves an estimate on the convergence of the filtered numerical solution to the corresponding DNS solution.

preprint2013arXiv

A finite element solver and energy stable coupling for 3D and 1D fluid models

The paper develops a solver based on a conforming finite element method for a 3D--1D coupled incompressible flow problem. New coupling conditions are introduced to ensure a suitable bound for the cumulative energy of the model. We study the stability and accuracy of the discretization method, and the performance of some state-of-the-art linear algebraic solvers for such flow configurations. Motivated by the simulation of the flow over inferior vena cava (IVC) filter, we consider the coupling of a 1D fluid model and a 3D fluid model posed in a domain with anisotropic inclusions. The relevance of our approach to realistic cardiovascular simulations is demonstrated by computing a blood flow over a model IVC filter.

preprint2013arXiv

A stabilized finite element method for advection-diffusion equations on surfaces

A recently developed Eulerian finite element method is applied to solve advection-diffusion equations posed on hypersurfaces. When transport processes on a surface dominate over diffusion, finite element methods tend to be unstable unless the mesh is sufficiently fine. The paper introduces a stabilized finite element formulation based on the SUPG technique. An error analysis of the method is given. Results of numerical experiments are presented that illustrate the performance of the stabilized method.

preprint2013arXiv

Non-degenerate Eulerian finite element method for solving PDEs on surfaces

The paper studies a method for solving elliptic partial differential equations posed on hypersurfaces in $\mathbb{R}^N$, $N=2,3$. The method builds upon the formulation introduced in Bertalmio et al., J. Comput. Phys., 174 (2001), 759--780., where a surface equation is extended to a neighborhood of the surface. The resulting degenerate PDE is then solved in one dimension higher, but can be solved on a mesh that is unaligned to the surface. We introduce another extended formulation, which leads to uniformly elliptic (non-degenerate) equations in a bulk domain containing the surface. We apply a finite element method to solve this extended PDE and prove the convergence of finite element solutions restricted to the surface to the solution of the original surface problem. Several numerical examples illustrate the properties of the method.

preprint2013arXiv

On surface meshes induced by level set functions

The zero level set of a piecewise-affine function with respect to a consistent tetrahedral subdivision of a domain in $\mathbb{R}^3$ is a piecewise-planar hyper-surface. We prove that if a family of consistent tetrahedral subdivions satisfies the minimum angle condition, then after a simple postprocessing this zero level set becomes a consistent surface triangulation which satisfies the maximum angle condition. We treat an application of this result to the numerical solution of PDEs posed on surfaces, using a $P_1$ finite element space on such a surface triangulation. For this finite element space we derive optimal interpolation error bounds. We prove that the diagonally scaled mass matrix is well-conditioned, uniformly with respect to $h$. Furthermore, the issue of conditioning of the stiffness matrix is addressed.