BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Alexei Lozinski

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

ResearcherUnclaimed source record
math.NAComputational Engineering, Finance, and Sciencemath-phmath.MPNumerical Analysisphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

13works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

13 published item(s)

preprint2022arXiv

Hierarchical a posteriori error estimation of Bank-Weiser type in the FEniCS Project

In the seminal paper of Bank and Weiser [Math. Comp., 44 (1985), pp.283-301] a new a posteriori estimator was introduced. This estimator requires the solution of a local Neumann problem on every cell of the finite element mesh. Despite the promise of Bank-Weiser type estimators, namely locality, computational efficiency, and asymptotic sharpness, they have seen little use in practical computational problems. The focus of this contribution is to describe a novel implementation of hierarchical estimators of the Bank-Weiser type in a modern high-level finite element software with automatic code generation capabilities. We show how to use the estimator to drive (goal-oriented) adaptive mesh refinement and to mixed approximations of the nearly-incompressible elasticity problems. We provide comparisons with various other used estimators. An open-source implementation based on the FEniCS Project finite element software is provided as supplementary material.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Quantifying discretization errors for soft-tissue simulation in computer assisted surgery: a preliminary study

Errors in biomechanics simulations arise from modeling and discretization. Modeling errors are due to the choice of the mathematical model whilst discretization errors measure the impact of the choice of the numerical method on the accuracy of the approximated solution to this specific mathematical model. A major source of discretization errors is mesh generation from medical images, that remains one of the major bottlenecks in the development of reliable, accurate, automatic and efficient personalized, clinically-relevant Finite Element (FE) models in biomechanics. The impact of mesh quality and density on the accuracy of the FE solution can be quantified with \emph{a posteriori} error estimates. Yet, to our knowledge, the relevance of such error estimates for practical biomechanics problems has seldom been addressed, see [25]. In this contribution, we propose an implementation of some a posteriori error estimates to quantify the discretization errors and to optimize the mesh. More precisely, we focus on error estimation for a user-defined quantity of interest with the Dual Weighted Residual (DWR) technique. We test its applicability and relevance in two situations, corresponding to computations for a tongue and an artery, using a simplified setting, i.e., plane linearized elasticity with contractility of the soft-tissue modeled as a pre-stress. Our results demonstrate the feasibility of such methodology to estimate the actual solution errors and to reduce them economically through mesh refinement.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical analysis of an asymptotic-preserving scheme for anisotropic elliptic equations

The main purpose of the present paper is to study from a numerical analysis point of view some robust methods designed to cope with stiff (highly anisotropic) elliptic problems. The so-called asymptotic-preserving schemes studied in this paper are very efficient in dealing with a wide range of $\varepsilon$-values, where $0 < \varepsilon \ll 1$ is the stiffness parameter, responsible for the high anisotropy of the problem. In particular, these schemes are even able to capture the macroscopic properties of the system, as $\varepsilon$ tends towards zero, while the discretization parameters remain fixed. The objective of this work shall be to prove some $\varepsilon$-independent convergence results for these numerical schemes and put hence some more rigor in the construction of such AP-methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A fictitious domain approach for Fluid-Structure Interactions based on the eXtended Finite Element Method

In this work we develop a fictitious domain method for the Stokes problem which allows computations in domains whose boundaries do not depend on the mesh. The method is based on the ideas of Xfem and has been first introduced for the Poisson problem. The fluid part is treated by a mixed finite element method, and a Dirichlet condition is imposed by a Lagrange multiplier on an immersed structure localized by a level-set function. A stabilization technique is carried out in order to get the convergence for this multiplier. The latter represents the forces that the fluid applies on the structure. The aim is to perform fluid-structure simulations for which these forces have a central role. We illustrate the capacities of the method by extending it to the incompressible Navier-Stokes equations coupled with a moving rigid solid.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A fictitious domain approach for the Stokes problem based on the extended finite element method

In the present work, we propose to extend to the Stokes problem a fictitious domain approach inspired by eXtended Finite Element Method and studied for Poisson problem in [Renard]. The method allows computations in domains whose boundaries do not match. A mixed finite element method is used for fluid flow. The interface between the fluid and the structure is localized by a level-set function. Dirichlet boundary conditions are taken into account using Lagrange multiplier. A stabilization term is introduced to improve the approximation of the normal trace of the Cauchy stress tensor at the interface and avoid the inf-sup condition between the spaces for velocity and the Lagrange multiplier. Convergence analysis is given and several numerical tests are performed to illustrate the capabilities of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-Conforming Multiscale Finite Element Method for Stokes Flows in Heterogeneous Media. Part I: Methodologies and Numerical Experiments

The Multiscale Finite Element Method (MsFEM) is developed in the vein of Crouzeix-Raviart element for solving viscous incompressible flows in genuine heterogeneous media. Such flows are relevant in many branches of engineering, often at multiple scales and at regions where analytical representations of the microscopic features of the flows are often unavailable. Full accounts to these problems heavily depend on the geometry of the system under consideration and are computationally expensive. Therefore, a method capable of solving multiscale features of the flow without confining itself to fine scale calculations is sought after. The approximation of boundary condition on coarse element edges when computing the multiscale basis functions critically influences the eventual accuracy of any MsFEM approaches. The weakly enforced continuity of Crouzeix - Raviart function space across element edges leads to a natural boundary condition for the multiscale basis functions which relaxes the sensitivity of our method to complex patterns of obstacles exempt from the needs of implementing any oversampling techniques. Additionally, the application of penalization method makes it possible to avoid complex unstructured domain and allows extensive use of simpler Cartesian meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An MsFEM type approach for perforated domains

We follow up on our previous work [C. Le Bris, F. Legoll and A. Lozinski, Chinese Annals of Mathematics 2013] where we have studied a multiscale finite element (MsFEM) type method in the vein of the classical Crouzeix-Raviart finite element method that is specifically adapted for highly oscillatory elliptic problems. We adapt the approach to address here a multiscale problem on a perforated domain. An additional ingredient of our approach is the enrichment of the multiscale finite element space using bubble functions. We first establish a theoretical error estimate. We next show that, on the problem we consider, the approach we propose outperforms all dedicated existing variants of MsFEM we are aware of.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Crouzeix-Raviart MsFEM with Bubble Functions for Diffusion and Advection-Diffusion in Perforated Media

The adaptation of Crouzeix - Raviart finite element in the context of multiscale finite element method (MsFEM) is studied and implemented on diffusion and advection-diffusion problems in perforated media. It is known that the approximation of boundary condition on coarse element edges when computing the multiscale basis functions critically influences the eventual accuracy of any MsFEM approaches. The weakly enforced continuity of Crouzeix - Raviart function space across element edges leads to a natural boundary condition for the multiscale basis functions which relaxes the sensitivity of our method to complex patterns of perforations. Another ingredient to our method is the application of bubble functions which is shown to be instrumental in maintaining high accuracy amid dense perforations. Additionally, the application of penalization method makes it possible to avoid complex unstructured domain and allows extensive use of simpler Cartesian meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Highly anisotropic temperature balance equation and its asymptotic-preserving resolution

This paper deals with the numerical study of a nonlinear, strongly anisotropic heat equation. The use of standard schemes in this situation leads to poor results, due to the high anisotropy. An Asymptotic-Preserving method is introduced in this paper, which is second-order accurate in both, temporal and spacial variables. The discretization in time is done using an L-stable Runge-Kutta scheme. The convergence of the method is shown to be independent of the anisotropy parameter $0 < \eps <1$, and this for fixed coarse Cartesian grids and for variable anisotropy directions. The context of this work are magnetically confined fusion plasmas.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An Asymptotic-Preserving method for highly anisotropic elliptic equations based on a micro-macro decomposition

The concern of the present work is the introduction of a very efficient Asymptotic Preserving scheme for the resolution of highly anisotropic diffusion equations. The characteristic features of this scheme are the uniform convergence with respect to the anisotropy parameter $0<\eps <<1$, the applicability (on cartesian grids) to cases of non-uniform and non-aligned anisotropy fields $b$ and the simple extension to the case of a non-constant anisotropy intensity $1/\eps$. The mathematical approach and the numerical scheme are different from those presented in the previous work [Degond et al. (2010), arXiv:1008.3405v1] and its considerable advantages are pointed out.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Duality-based Asymptotic-Preserving method for highly anisotropic diffusion equations

The present paper introduces an efficient and accurate numerical scheme for the solution of a highly anisotropic elliptic equation, the anisotropy direction being given by a variable vector field. This scheme is based on an asymptotic preserving reformulation of the original system, permitting an accurate resolution independently of the anisotropy strength and without the need of a mesh adapted to this anisotropy. The counterpart of this original procedure is the larger system size, enlarged by adding auxiliary variables and Lagrange multipliers. This Asymptotic-Preserving method generalizes the method investigated in a previous paper [arXiv:0903.4984v2] to the case of an arbitrary anisotropy direction field.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Kinetic models for dilute solutions of dumbbells in non-homogeneous flows revisited

We propose a two fluid theory to model a dilute polymer solution assuming that it consists of two phases, polymer and solvent, with two distinct macroscopic velocities. The solvent phase velocity is governed by the macroscopic Navier-Stokes equations with the addition of a force term describing the interaction between the two phases. The polymer phase is described on the mesoscopic level using a dumbbell model and its macroscopic velocity is obtained through averaging. We start by writing down the full phase-space distribution function for the dumbbells and then obtain the inertialess limits for the Fokker-Planck equation and for the averaged friction force acting between the phases from a rigorous asymptotic analysis. The resulting equations are relevant to the modelling of strongly non-homogeneous flows, while the standard kinetic model is recovered in the locally homogeneous case.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On discretization in time in simulations of particulate flows

We propose a time discretization scheme for a class of ordinary differential equations arising in simulations of fluid/particle flows. The scheme is intended to work robustly in the lubrication regime when the distance between two particles immersed in the fluid or between a particle and the wall tends to zero. The idea consists in introducing a small threshold for the particle-wall distance below which the real trajectory of the particle is replaced by an approximated one where the distance is kept equal to the threshold value. The error of this approximation is estimated both theoretically and by numerical experiments. Our time marching scheme can be easily incorporated into a full simulation method where the velocity of the fluid is obtained by a numerical solution to Stokes or Navier-Stokes equations. We also provide a derivation of the asymptotic expansion for the lubrication force (used in our numerical experiments) acting on a disk immersed in a Newtonian fluid and approaching the wall. The method of this derivation is new and can be easily adapted to other cases.

0 citations0 reviewsNo code linkedOpen access