Source author record

Maxim Olshanskii

Maxim Olshanskii 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 physics.comp-ph Biological Physics cond-mat.soft

Catalog footprint

What is connected

8works

5topics

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.

preprint2022arXiv

A CutFEM divergence-free discretization for the Stokes problem

We construct and analyze a CutFEM discretization for the Stokes problem based on the Scott-Vogelius pair. The discrete piecewise polynomial spaces are defined on macro-element triangulations which are not fitted to the smooth physical domain. Boundary conditions are imposed via penalization through the help of a Nitsche-type discretization, whereas stability with respect to small and anisotropic cuts of the bulk elements is ensured by adding local ghost penalty stabilization terms. We show stability of the scheme as well as a divergence--free property of the discrete velocity outside an $O(h)$ neighborhood of the boundary. To mitigate the error caused by the violation of the divergence-free condition, we introduce local grad-div stabilization. The error analysis shows that the grad-div parameter can scale like $O(h^{-1})$, allowing a rather heavy penalty for the violation of mass conservation, while still ensuring optimal order error estimates.

preprint2022arXiv

Recycling augmented Lagrangian preconditioner in an incompressible fluid solver

The paper discusses a reuse of matrix factorization as a building block in the Augmented Lagrangian (AL) and modified AL preconditioners for non-symmetric saddle point linear algebraic systems. The strategy is applied to solve two-dimensional incompressible fluid problems with efficiency rates independent of the Reynolds number. The solver is then tested to simulate motion of a surface fluid, an example of a 2D flow motivated by an interest in lateral fluidity of inextensible viscous membranes. Numerical examples include the Kelvin--Helmholtz instability problem posed on the sphere and on the torus. Some new eigenvalue estimates for the AL preconditioner are derived.

preprint2021arXiv

An unfitted finite element method for two-phase Stokes problems with slip between phases

We present an isoparametric unfitted finite element approach of the CutFEM or Nitsche-XFEM family for the simulation of two-phase Stokes problems with slip between phases. For the unfitted generalized Taylor--Hood finite element pair $\mathbf{P}_{k+1}-P_k$, $k\ge1$, we show an inf-sup stability property with a stability constant that is independent of the viscosity ratio, slip coefficient, position of the interface with respect to the background mesh and, of course, mesh size. In addition, we prove stability and optimal error estimates that follow from this inf-sup property. We provide numerical results in two and three dimensions to corroborate the theoretical findings and demonstrate the robustness of our approach with respect to the contrast in viscosity, slip coefficient value, and position of the interface relative to the fixed computational mesh.

preprint2020arXiv

Experimental validation of a phase-field model to predict coarsening dynamics of lipid domains in multicomponent membranes

Membrane phase-separation is a mechanism that biological membranes often use to locally concentrate specific lipid species in order to organize diverse membrane processes. Phase separation has also been explored as a tool for the design of liposomes with heterogeneous and spatially organized surfaces. These "patchy" liposomes are promising platforms for delivery purposes, however their design and optimization through experimentation can be expensive and time-consuming. We developed a computationally efficient method based on the surface Cahn-Hilliard phase-field model to complement experimental investigations in the design of patchy liposomes. The method relies on thermodynamic considerations to set the initial state for numerical simulations. We show that our computational approach delivers not only qualitative pictures, but also accurate quantitative information about the dynamics of the membrane organization. In particular, the computational and experimental results are in excellent agreement in terms of raft area fraction, total raft perimeter over time and total number of rafts over time for two different membrane compositions (DOPC:DPPC with a 2:1 molar ratio with 20% Chol and DOPC:DPPC with a 3:1 molar ratio with 20% Chol). Thus, the computational phase-field model informed by experiments has a considerable potential to assist in the design of liposomes with spatially organized surfaces, thereby containing the cost and time required by the design process.

preprint2020arXiv

New stability estimates for an unfitted finite element method for two-phase Stokes problem

The paper addresses stability and finite element analysis of the stationary two-phase Stokes problem with a piecewise constant viscosity coefficient experiencing a jump across the interface between two fluid phases. We first prove a priori estimates for the individual terms of the Cauchy stress tensor with stability constants independent of the viscosity coefficient. Next, this stability result is extended to the approximation of the two-phase Stokes problem by a finite element method. In the method considered, the interface between the phases does not respect the underlying triangulation, which put the finite element method into the class of unfitted discretizations. The finite element error estimates are proved with constants independent of viscosity. Numerical experiments supporting the theoretical results are provided.

preprint2019arXiv

Numerical modelling of phase separation on dynamic surfaces

The paper presents a model of lateral phase separation in a two component material surface. The resulting fourth order nonlinear PDE can be seen as a Cahn-Hilliard equation posed on a time-dependent surface. Only elementary tangential calculus and the embedding of the surface in $\mathbb{R}^3$ are used to formulate the model, thereby facilitating the development of a fully Eulerian discretization method to solve the problem numerically. A hybrid method, finite difference in time and trace finite element in space, is introduced and stability of its semi-discrete version is proved. The method avoids any triangulation of the surface and uses a surface-independent background mesh to discretize the equation. Thus, the method is capable of solving the Cahn-Hilliard equation numerically on implicitly defined surfaces and surfaces undergoing strong deformations and topological transitions. We assess the approach on a set of test problems and apply it to model spinodal decomposition and pattern formation on colliding surfaces. Finally, we consider the phase separation on a sphere splitting into two droplets.

preprint2016arXiv

Numerical integration over implicitly defined domains for higher order unfitted finite element methods

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.

preprint2014arXiv

A space-time FEM for PDEs on evolving surfaces

The paper studies a finite element method for computing transport and diffusion along evolving surfaces. The method does not require a parametrization of a surface or an extension of a PDE from a surface into a bulk outer domain. The surface and its evolution may be given implicitly, e.g., as the solution of a level set equation. This approach naturally allows a surface to undergo topological changes and experience local geometric singularities. The numerical method uses space-time finite elements and is provably second order accurate. The paper reviews the method, error estimates and shows results for computing the diffusion of a surfactant on surfaces of two colliding droplets.