Showing 1–38 of 38 results for author: Lozinski, A

MsFEM for advection-dominated problems in heterogeneous media: Stabilization via nonconforming variants

Authors: Rutger A. Biezemans, Claude Le Bris, Frédéric Legoll, Alexei Lozinski

Abstract: We study the numerical approximation of advection-diffusion equations with highly oscillatory coefficients and possibly dominant advection terms by means of the Multiscale Finite Element Method. The latter method is a now classical, finite element type method that performs a Galerkin approximation on a problem-dependent basis set, itself pre-computed in an offline stage. The approach is implemente… ▽ More We study the numerical approximation of advection-diffusion equations with highly oscillatory coefficients and possibly dominant advection terms by means of the Multiscale Finite Element Method. The latter method is a now classical, finite element type method that performs a Galerkin approximation on a problem-dependent basis set, itself pre-computed in an offline stage. The approach is implemented here using basis functions that locally resolve both the diffusion and the advection terms. Variants with additional bubble functions and possibly weak inter-element continuity are proposed. Some theoretical arguments and a comprehensive set of numerical experiments allow to investigate and compare the stability and the accuracy of the approaches. The best approach constructed is shown to be adequate for both the diffusion- and advection-dominated regimes, and does not rely on an auxiliary stabilization parameter that would have to be properly adjusted. △ Less

Submitted 19 June, 2024; originally announced June 2024.

An a posteriori error estimate for a 0D/2D coupled model

Authors: Hussein Albazzal, Alexei Lozinski, Roberta Tittarelli

Abstract: This work is motivated by the need of efficient numerical simulations of gas flows in the serpentine channels used in proton-exchange membrane fuel cells. In particular, we consider the Poisson problem in a 2D domain composed of several long straight rectangular sections and of several bends corners. In order to speed up the resolution, we propose a 0D model in the rectangular parts of the channel… ▽ More This work is motivated by the need of efficient numerical simulations of gas flows in the serpentine channels used in proton-exchange membrane fuel cells. In particular, we consider the Poisson problem in a 2D domain composed of several long straight rectangular sections and of several bends corners. In order to speed up the resolution, we propose a 0D model in the rectangular parts of the channel and a Finite Element resolution in the bends. To find a good compromise between precision and time consuming, the challenge is double: how to choose a suitable position of the interface between the 0D and the 2D models and how to control the discretization error in the bends. We shall present an \textit{a posteriori} error estimator based on an equilibrated flux reconstruction in the subdomains where the Finite Element method is applied. The estimates give a global upper bound on the error measured in the energy norm of the difference between the exact and approximate solutions on the whole domain. They are guaranteed, meaning that they feature no undetermined constants. (global) Lower bounds for the error are also derived. An adaptive algorithm is proposed to use smartly the estimator for aforementioned double challenge. A numerical validation of the estimator and the algorithm completes the work. \end{abstract} △ Less

Submitted 13 December, 2023; originally announced December 2023.

Nonlinear reduced basis using mixture Wasserstein barycenters: application to an eigenvalue problem inspired from quantum chemistry

Authors: Maxime Dalery, Genevieve Dusson, Virginie Ehrlacher, Alexei Lozinski

Abstract: The aim of this article is to propose a new reduced-order modelling approach for parametric eigenvalue problems arising in electronic structure calculations. Namely, we develop nonlinear reduced basis techniques for the approximation of parametric eigenvalue problems inspired from quantum chemistry applications. More precisely, we consider here a one-dimensional model which is a toy model for the… ▽ More The aim of this article is to propose a new reduced-order modelling approach for parametric eigenvalue problems arising in electronic structure calculations. Namely, we develop nonlinear reduced basis techniques for the approximation of parametric eigenvalue problems inspired from quantum chemistry applications. More precisely, we consider here a one-dimensional model which is a toy model for the computation of the electronic ground state wavefunction of a system of electrons within a molecule, solution to the many-body electronic Schrödinger equation, where the varying parameters are the positions of the nuclei in the molecule. We estimate the decay rate of the Kolmogorov n-width of the set of solutions for this parametric problem in several settings, including the standard L2-norm as well as with distances based on optimal transport. The fact that the latter decays much faster than in the traditional L2-norm setting motivates us to propose a practical nonlinear reduced basis method, which is based on an offline greedy algorithm, and an efficient stochastic energy minimization in the online phase. We finally provide numerical results illustrating the capabilities of the method and good approximation properties, both in the offline and the online phase. △ Less

Submitted 28 July, 2023; originally announced July 2023.

A penalty-free Shifted Boundary Method of arbitrary order

Authors: J. Haydel Collins, Alexei Lozinski, Guglielmo Scovazzi

Abstract: We introduce and analyze a penalty-free formulation of the Shifted Boundary Method (SBM), inspired by the asymmetric version of the Nitsche method. We prove its stability and convergence for arbitrary order finite element interpolation spaces and we test its performance with a number of numerical experiments. Moreover, while the SBM was previously believed to be only asymptotically consistent (in… ▽ More We introduce and analyze a penalty-free formulation of the Shifted Boundary Method (SBM), inspired by the asymmetric version of the Nitsche method. We prove its stability and convergence for arbitrary order finite element interpolation spaces and we test its performance with a number of numerical experiments. Moreover, while the SBM was previously believed to be only asymptotically consistent (in the sense of Galerkin orthogonality), we prove here that it is indeed exactly consistent. △ Less

Submitted 22 June, 2023; originally announced June 2023.

phi-FEM for the heat equation: optimal convergence on unfitted meshes in space

Authors: Michel Duprez, Vanessa Lleras, Alexei Lozinski, Killian Vuillemot

Abstract: Thanks to a finite element method, we solve numerically parabolic partial differential equations on complex domains by avoiding the mesh generation, using a regular background mesh, not fitting the domain and its real boundary exactly. Our technique follows the phi-FEM paradigm, which supposes that the domain is given by a level-set function. In this paper, we prove a priori error estimates in l2(… ▽ More Thanks to a finite element method, we solve numerically parabolic partial differential equations on complex domains by avoiding the mesh generation, using a regular background mesh, not fitting the domain and its real boundary exactly. Our technique follows the phi-FEM paradigm, which supposes that the domain is given by a level-set function. In this paper, we prove a priori error estimates in l2(H1) and linf(L2) norms for an implicit Euler discretization in time. We give numerical illustrations to highlight the performances of phi-FEM, which combines optimal convergence accuracy, easy implementation process and fastness. △ Less

Submitted 21 March, 2023; originally announced March 2023.

MSC Class: 65M60

Non-intrusive implementation of a wide variety of Multiscale Finite Element Methods

Authors: Rutger A. Biezemans, Claude Le Bris, Frédéric Legoll, Alexei Lozinski

Abstract: Multiscale Finite Element Methods (MsFEMs) are now well-established finite element type approaches dedicated to multiscale problems. They first compute local, oscillatory, problem-dependent basis functions that generate a suitable discretization space, and next perform a Galerkin approximation of the problem on that space. We investigate here how these approaches can be implemented in a non-intrus… ▽ More Multiscale Finite Element Methods (MsFEMs) are now well-established finite element type approaches dedicated to multiscale problems. They first compute local, oscillatory, problem-dependent basis functions that generate a suitable discretization space, and next perform a Galerkin approximation of the problem on that space. We investigate here how these approaches can be implemented in a non-intrusive way, in order to facilitate their dissemination within industrial codes or non-academic environments. We develop an abstract framework that covers a wide variety of MsFEMs for linear second-order partial differential equations. Non-intrusive MsFEM approaches are developed within the full generality of this framework, which may moreover be beneficial to steering software development and improving the theoretical understanding and analysis of MsFEMs. △ Less

Submitted 2 August, 2023; v1 submitted 30 November, 2022; originally announced November 2022.

Comments: 50 pages, 4 figures; typos corrected (notably in the proof of Lemma 3), some clarifications added, numerical results added in Section 7

Phi-FEM: an optimally convergent and easily implementable immersed boundary method for particulate flows and Stokes equations

Authors: Michel Duprez, Vanessa Lleras, Alexei Lozinski

Abstract: We present an immersed boundary method to simulate the creeping motion of a rigid particle in a fluid described by the Stokes equations discretized thanks to a finite element strategy on unfitted meshes, called Phi-FEM, that uses the description of the solid with a level-set function. One of the advantages of our method is the use of standard finite element spaces and classical integration tools,… ▽ More We present an immersed boundary method to simulate the creeping motion of a rigid particle in a fluid described by the Stokes equations discretized thanks to a finite element strategy on unfitted meshes, called Phi-FEM, that uses the description of the solid with a level-set function. One of the advantages of our method is the use of standard finite element spaces and classical integration tools, while maintaining the optimal convergence (theoretically in the H1 norm for the velocity and L2 for pressure; numerically also in the L2 norm for the velocity). △ Less

Submitted 27 January, 2023; v1 submitted 13 November, 2022; originally announced November 2022.

MSC Class: 65N30; 65N85; 65N15

Non-intrusive implementation of Multiscale Finite Element Methods: an illustrative example

Authors: Rutger A. Biezemans, Claude Le Bris, Frederic Legoll, Alexei Lozinski

Abstract: Multiscale Finite Element Methods (MsFEM) are finite element type approaches dedicated to multiscale problems. They first compute local, oscillatory, problem-dependent basis functions which generate a specific discretization space, and next perform a Galerkin approximation of the problem on that space. We investigate here how these approaches can be implemented in a non-intrusive way, in order to… ▽ More Multiscale Finite Element Methods (MsFEM) are finite element type approaches dedicated to multiscale problems. They first compute local, oscillatory, problem-dependent basis functions which generate a specific discretization space, and next perform a Galerkin approximation of the problem on that space. We investigate here how these approaches can be implemented in a non-intrusive way, in order to facilitate their dissemination within industrial codes or non academic environments. △ Less

Submitted 14 April, 2022; originally announced April 2022.

$φ$-FEM: an efficient simulation tool using simple meshes for problems in structure mechanics and heat transfer

Authors: Stephane Cotin, Michel Duprez, Vanessa Lleras, Alexei Lozinski, Killian Vuillemot

Abstract: One of the major issues in the computational mechanics is to take into account the geometrical complexity. To overcome this difficulty and to avoid the expensive mesh generation, geometrically unfitted methods, i.e. the numerical methods using the simple computational meshes that do not fit the boundary of the domain, and/or the internal interfaces, have been widely developed. In the present work,… ▽ More One of the major issues in the computational mechanics is to take into account the geometrical complexity. To overcome this difficulty and to avoid the expensive mesh generation, geometrically unfitted methods, i.e. the numerical methods using the simple computational meshes that do not fit the boundary of the domain, and/or the internal interfaces, have been widely developed. In the present work, we investigate the performances of an unfitted method called $φ$-FEM that converges optimally and uses classical finite element spaces so that it can be easily implemented using general FEM libraries. The main idea is to take into account the geometry thanks to a level set function describing the boundary or the interface. Up to now, the $φ$-FEM approach has been proposed, tested and substantiated mathematically only in some simplest settings: Poisson equation with Dirichlet/Neumann/Robin boundary conditions. Our goal here is to demonstrate its applicability to some more sophisticated governing equations arising in the computational mechanics. We consider the linear elasticity equations accompanied by either pure Dirichlet boundary conditions or by the mixed ones (Dirichlet and Neumann boundary conditions co-existing on parts of the boundary), an interface problem (linear elasticity with material coefficients abruptly changing over an internal interface), a model of elastic structures with cracks, and finally the heat equation. In all these settings, we derive an appropriate variant of $φ$-FEM and then illustrate it by numerical tests on manufactured solutions. We also compare the accuracy and efficiency of $φ$-FEM with those of the standard fitted FEM on the meshes of similar size, revealing the substantial gains that can be achieved by $φ$-FEM in both the accuracy and the computational time. △ Less

Submitted 11 October, 2021; originally announced October 2021.

MSC Class: 65N30; 65N85; 65N15

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

Authors: Raphaël Bulle, Jack S. Hale, Alexei Lozinski, Stéphane P. A. Bordas, Franz Chouly

Abstract: 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 computationa… ▽ More 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. △ Less

Submitted 7 February, 2022; v1 submitted 8 February, 2021; originally announced February 2021.

A new $φ$-FEM approach for problems with natural boundary conditions

Authors: Michel Duprez, Vanessa Lleras, Alexei Lozinski

Abstract: We present a new finite element method, called $φ$-FEM, to solve numerically elliptic partial differential equations with natural (Neumann or Robin) boundary conditions using simple computational grids, not fitted to the boundary of the physical domain. The boundary data are taken into account using a level-set function, which is a popular tool to deal with complicated or evolving domains. Our app… ▽ More We present a new finite element method, called $φ$-FEM, to solve numerically elliptic partial differential equations with natural (Neumann or Robin) boundary conditions using simple computational grids, not fitted to the boundary of the physical domain. The boundary data are taken into account using a level-set function, which is a popular tool to deal with complicated or evolving domains. Our approach belongs to the family of fictitious domain methods (or immersed boundary methods) and is close to recent methods of cutFEM/XFEM type. Contrary to the latter, $φ$-FEM does not need any non-standard numerical integration on cut mesh elements or on the actual boundary, while assuring the optimal convergence orders with finite elements of any degree and providing reasonably well conditioned discrete problems. In the first version of $φ$-FEM, only essential (Dirichlet) boundary conditions was considered. Here, to deal with natural boundary conditions, we introduce the gradient of the primary solution as an auxiliary variable. This is done only on the mesh cells cut by the boundary, so that the size of the numerical system is only slightly increased. We prove theoretically the optimal convergence of our scheme and a bound on the discrete problem conditioning, independent of the mesh cuts. The numerical experiments confirm these results. △ Less

Submitted 5 December, 2020; v1 submitted 25 March, 2020; originally announced March 2020.

MSC Class: 65N30; 65N85

$φ$-FEM: a finite element method on domains defined by level-sets

Authors: Michel Duprez, Alexei Lozinski

Abstract: We propose a new fictitious domain finite element method, well suited for elliptic problems posed in a domain given by a level-set function without requiring a mesh fitting the boundary. To impose the Dirichlet boundary conditions, we search the approximation to the solution as a product of a finite element function with the given level-set function, which also approximated by finite elements. Unl… ▽ More We propose a new fictitious domain finite element method, well suited for elliptic problems posed in a domain given by a level-set function without requiring a mesh fitting the boundary. To impose the Dirichlet boundary conditions, we search the approximation to the solution as a product of a finite element function with the given level-set function, which also approximated by finite elements. Unlike other recent fictitious domain-type methods (XFEM, CutFEM), our approach does not need any non-standard numerical integration (on cut mesh elements or on the actual boundary). We consider the Poisson equation discretized with piecewise polynomial Lagrange finite elements of any order and prove the optimal convergence of our method in the $H^1$-norm. Moreover, the discrete problem is proven to be well conditioned, \textit{i.e.} the condition number of the associated finite element matrix is of the same order as that of a standard finite element method on a comparable conforming mesh. Numerical results confirm the optimal convergence in both $H^1$ and $L^2$ norms. △ Less

Submitted 5 July, 2019; v1 submitted 8 March, 2019; originally announced March 2019.

MSC Class: 65N30; 65N85; 65N15

CutFEM without cutting the mesh cells: a new way to impose Dirichlet and Neumann boundary conditions on unfitted meshes

Authors: Alexei Lozinski

Abstract: We present a method of CutFEM type for the Poisson problem with either Dirichlet or Neumann boundary conditions. The computational mesh is obtained from a background (typically uniform Cartesian) mesh by retaining only the elements intersecting the domain where the problem is posed. The resulting mesh does not thus fit the boundary of the problem domain. Several finite element methods (XFEM, CutFE… ▽ More We present a method of CutFEM type for the Poisson problem with either Dirichlet or Neumann boundary conditions. The computational mesh is obtained from a background (typically uniform Cartesian) mesh by retaining only the elements intersecting the domain where the problem is posed. The resulting mesh does not thus fit the boundary of the problem domain. Several finite element methods (XFEM, CutFEM) adapted to such meshes have been recently proposed. The originality of the present article consists in avoiding integration over the elements cut by the boundary of the problem domain, while preserving the optimal convergence rates, as confirmed by both the theoretical estimates and the numerical results. △ Less

Submitted 5 July, 2019; v1 submitted 13 January, 2019; originally announced January 2019.

MSC Class: 65N30; 65N85; 65N15

Finite element method with local damage on the mesh

Authors: Michel Duprez, Vanessa Lleras, Alexei Lozinski

Abstract: We consider the finite element method on locally damaged meshes allowing for some distorted cells which are isolated from one another. In the case of the Poisson equation and piecewise linear Lagrange finite elements, we show that the usual a priori error estimates remain valid on such meshes. We also propose an alternative finite element scheme which is optimally convergent and, moreover, well co… ▽ More We consider the finite element method on locally damaged meshes allowing for some distorted cells which are isolated from one another. In the case of the Poisson equation and piecewise linear Lagrange finite elements, we show that the usual a priori error estimates remain valid on such meshes. We also propose an alternative finite element scheme which is optimally convergent and, moreover, well conditioned, i.e. its conditioning is of the same order as that of a standard finite element method on a regular mesh of comparable size. △ Less

Submitted 20 August, 2018; originally announced August 2018.

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

Authors: Michel Duprez, Stéphane P. A. Bordas, Marek Bucki, Huu Phuoc Bui, Franz Chouly, Vanessa Lleras, Claudio Lobos, Alexei Lozinski, Pierre-Yves Rohan, Satyendra Tomar

Abstract: 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… ▽ More 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. △ Less

Submitted 18 June, 2018; originally announced June 2018.

Comments: 35 pages, 10 figures

Journal ref: Applied Mathematical Modelling Volume 77, Part 1, January 2020, Pages 709-723

A primal discontinuous Galerkin method with static condensation on very general meshes

Authors: Alexei Lozinski

Abstract: We propose an efficient variant of a primal Discontinuous Galerkin method with interior penalty for the second order elliptic equations on very general meshes (polytopes with eventually curved boundaries). Efficiency, especially when higher order polynomials are used, is achieved by static condensation, i.e. a local elimination of certain degrees of freedom element by element. This alters the orig… ▽ More We propose an efficient variant of a primal Discontinuous Galerkin method with interior penalty for the second order elliptic equations on very general meshes (polytopes with eventually curved boundaries). Efficiency, especially when higher order polynomials are used, is achieved by static condensation, i.e. a local elimination of certain degrees of freedom element by element. This alters the original method in a way that preserves the optimal error estimates. Numerical experiments confirm that the solutions produced by the new method are indeed very close to that produced by the classical one. △ Less

Submitted 15 March, 2019; v1 submitted 19 March, 2018; originally announced March 2018.

Non-Conforming Multiscale Finite Element Method for Stokes Flows in Heterogeneous Media. Part II: error estimates for periodic microstructure

Authors: Gaspard Jankowiak, Alexei Lozinski

Abstract: This paper is dedicated to the rigorous numerical analysis of a Multiscale Finite Element Method (MsFEM) for the Stokes system, when dealing with highly heterogeneous media, as proposed in [B.P.~Muljadi et al., arXiv:1404.2837]. The method is in the vein of the classical Crouzeix-Raviart approach. It is generalized here to arbitrary sets of weighting functions used to enforced continuity across th… ▽ More This paper is dedicated to the rigorous numerical analysis of a Multiscale Finite Element Method (MsFEM) for the Stokes system, when dealing with highly heterogeneous media, as proposed in [B.P.~Muljadi et al., arXiv:1404.2837]. The method is in the vein of the classical Crouzeix-Raviart approach. It is generalized here to arbitrary sets of weighting functions used to enforced continuity across the mesh edges. We provide error bounds for a particular set of weighting functions in a periodic setting, using an accurate estimate of the homogenization error. Numerical experiments demonstrate an improved accuracy of the present variant with respect to that of Part I, both in the periodic case and in a broader setting. △ Less

Submitted 12 February, 2018; originally announced February 2018.

MSC Class: 35J15; 65N12; 65N30

An easily computable error estimator in space and time for the wave equation

Authors: Olga Gorynina, Alexei Lozinski, Marco Picasso

Abstract: We propose a cheaper version of \textit{a posteriori} error estimator from arXiv:1707.00057 for the linear second-order wave equation discretized by the Newmark scheme in time and by the finite element method in space. The new estimator preserves all the properties of the previous one (reliability, optimality on smooth solutions and quasi-uniform meshes) but no longer requires an extra computation… ▽ More We propose a cheaper version of \textit{a posteriori} error estimator from arXiv:1707.00057 for the linear second-order wave equation discretized by the Newmark scheme in time and by the finite element method in space. The new estimator preserves all the properties of the previous one (reliability, optimality on smooth solutions and quasi-uniform meshes) but no longer requires an extra computation of the Laplacian of the discrete solution on each time step. △ Less

Submitted 22 October, 2017; originally announced October 2017.

Comments: arXiv admin note: substantial text overlap with arXiv:1707.00057

Stability and optimal convergence of unfitted extended finite element methods with Lagrange multipliers for the Stokes equations

Authors: Michel Fournié, Alexei Lozinski

Abstract: We study a fictitious domain approach with Lagrange multipliers to discretize Stokes equations on a mesh that does not fit the boundaries. A mixed finite element method is used for fluid flow. Several stabilization terms are added to improve the approximation of the normal trace of the stress tensor and to avoid the inf-sup conditions between the spaces of the velocity and the Lagrange multipliers… ▽ More We study a fictitious domain approach with Lagrange multipliers to discretize Stokes equations on a mesh that does not fit the boundaries. A mixed finite element method is used for fluid flow. Several stabilization terms are added to improve the approximation of the normal trace of the stress tensor and to avoid the inf-sup conditions between the spaces of the velocity and the Lagrange multipliers. We generalize first an approach based on eXtended Finite Element Method due to Haslinger-Renard involving a Barbosa-Hughes stabilization and a robust reconstruction on the badly cut elements. Secondly, we adapt the approach due to Burman-Hansbo involving a stabilization only on the Lagrange multiplier. Multiple choices for the finite elements for velocity, pressure and multiplier are considered. Additional stabilization on pressure (Brezzi-Pitkäranta, Interior Penalty) is added, if needed. We prove the stability and the optimal convergence of several variants of these methods under appropriate assumptions. Finally, we perform numerical tests to illustrate the capabilities of the methods. △ Less

Submitted 22 October, 2017; originally announced October 2017.

Time and space adaptivity of the wave equation discretized in time by a second order scheme

Authors: Olga Gorynina, Alexei Lozinski, Marco Picasso

Abstract: The aim of this paper is to obtain a posteriori error bounds of optimal order in time and space for the linear second-order wave equation discretized by the Newmark scheme in time and the finite element method in space. Error estimates are derived in the $L^{\infty}$-in-time/energy-in-space norm. Numerical experiments are reported for several test cases and confirm equivalence of the proposed esti… ▽ More The aim of this paper is to obtain a posteriori error bounds of optimal order in time and space for the linear second-order wave equation discretized by the Newmark scheme in time and the finite element method in space. Error estimates are derived in the $L^{\infty}$-in-time/energy-in-space norm. Numerical experiments are reported for several test cases and confirm equivalence of the proposed estimators and the true error. △ Less

Submitted 21 September, 2017; v1 submitted 30 June, 2017; originally announced July 2017.

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

Authors: Alexei Lozinski, Jacek Narski, Claudia Negulescu

Abstract: 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… ▽ More 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. △ Less

Submitted 3 July, 2015; originally announced July 2015.

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

Authors: Bagus Putra Muljadi, Jacek Narski, Alexei Lozinski, Pierre Degond

Abstract: 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 hea… ▽ More 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. △ Less

Submitted 10 April, 2014; originally announced April 2014.

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

Authors: Sébastien Court, Michel Fournié, Alexei Lozinski

Abstract: 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 lo… ▽ More 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. △ Less

Submitted 2 January, 2014; originally announced January 2014.

Comments: 10 pages, 4 figures

MSC Class: 65M60; 74F10; 65M85

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

Authors: Pierre Degond, Alexei Lozinski, Bagus Putra Muljadi, Jacek Narski

Abstract: 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.… ▽ More 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. △ Less

Submitted 29 September, 2013; originally announced October 2013.

An MsFEM type approach for perforated domains

Authors: Claude Le Bris, Frederic Legoll, Alexei Lozinski

Abstract: 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… ▽ More 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. △ Less

Submitted 2 July, 2013; originally announced July 2013.

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

Authors: Sébastien Court, Michel Fournié, Alexei Lozinski

Abstract: 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. Diri… ▽ More 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. △ Less

Submitted 2 January, 2014; v1 submitted 27 March, 2013; originally announced March 2013.

Comments: 28 pages, 12 figures

MSC Class: 65M60; 74F10; 65M85

Highly anisotropic temperature balance equation and its asymptotic-preserving resolution

Authors: Alexei Lozinski, Jacek Narski, Claudia Negulescu

Abstract: 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… ▽ More 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. △ Less

Submitted 30 March, 2012; originally announced March 2012.

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

Authors: Pierre Degond, Alexei Lozinski, Jacek Narski, Claudia Negulescu

Abstract: 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 si… ▽ More 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. △ Less

Submitted 4 February, 2011; originally announced February 2011.

Journal ref: Journal of Computational Physics,231 (2012), pp. 2724-2740

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

Authors: Pierre Degond, Fabrice Deluzet, Alexei Lozinski, Jacek Narski, Claudia Negulescu

Abstract: 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… ▽ More 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. △ Less

Submitted 19 August, 2010; originally announced August 2010.

Journal ref: Communications in Mathematical Sciences, 10 (2012), pp. 1-31

On discretization in time in simulations of particulate flows

Authors: Matthieu Hillairet, Alexei Lozinski, Marcela Szopos

Abstract: 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… ▽ More 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. △ Less

Submitted 24 March, 2010; originally announced March 2010.

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

Authors: Pierre Degond, Alexei Lozinski, Robert G. Owens

Abstract: 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… ▽ More 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. △ Less

Submitted 14 January, 2010; originally announced January 2010.

Journal ref: Journal of Non-Newtonian Fluid Mechanics, 165 (2010), pp. 509-518

Geometric phase induced by a cyclically evolving squeezed vacuum reservoir

Authors: Angelo Carollo, G. Massimo Palma, Artur Lozinski, Marcelo Franca Santos, Vlatko Vedral

Abstract: We propose a new way to generate an observable geometric phase by means of a completely incoherent phenomenon. We show how to imprint a geometric phase to a system by "adiabatically" manipulating the environment with which it interacts. As a specific scheme we analyse a multilevel atom interacting with a broad-band squeezed vacuum bosonic bath. As the squeezing parameters are smoothly changed in… ▽ More We propose a new way to generate an observable geometric phase by means of a completely incoherent phenomenon. We show how to imprint a geometric phase to a system by "adiabatically" manipulating the environment with which it interacts. As a specific scheme we analyse a multilevel atom interacting with a broad-band squeezed vacuum bosonic bath. As the squeezing parameters are smoothly changed in time along a closed loop, the ground state of the system acquires a geometric phase. We propose also a scheme to measure such geometric phase by means of a suitable polarization detection. △ Less

Submitted 21 April, 2006; v1 submitted 11 July, 2005; originally announced July 2005.

Comments: 4 pages, 1 figure

Journal ref: Phys. Rev. Lett. 96, 150403 (2006)

Berry phase in open quantum systems: a quantum Langevin equation approach

Authors: Gabriele De Chiara, Artur Lozinski, G. Massimo Palma

Abstract: The evolution of a two level system with a slowly varying Hamiltonian, modeled as s spin 1/2 in a slowly varying magnetic field, and interacting with a quantum environment, modeled as a bath of harmonic oscillators is analyzed using a quantum Langevin approach. This allows to easily obtain the dissipation time and the correction to the Berry phase in the case of an adiabatic cyclic evolution. The evolution of a two level system with a slowly varying Hamiltonian, modeled as s spin 1/2 in a slowly varying magnetic field, and interacting with a quantum environment, modeled as a bath of harmonic oscillators is analyzed using a quantum Langevin approach. This allows to easily obtain the dissipation time and the correction to the Berry phase in the case of an adiabatic cyclic evolution. △ Less

Submitted 8 January, 2007; v1 submitted 22 October, 2004; originally announced October 2004.

Comments: 6 pages, no figures. Published version

Journal ref: Eur. Phys. J. D 41, 179-183 (2007)

Quantum dynamical entropy and decoherence rate

Authors: Robert Alicki, Artur Lozinski, Prot Pakonski, Karol Zyczkowski

Abstract: We investigate quantum dynamical systems defined on a finite dimensional Hilbert space and subjected to an interaction with an environment. The rate of decoherence of initially pure states, measured by the increase of their von Neumann entropy, averaged over an ensemble of random pure states, is proved to be bounded from above by the partial entropy used to define the ALF dynamical entropy. The… ▽ More We investigate quantum dynamical systems defined on a finite dimensional Hilbert space and subjected to an interaction with an environment. The rate of decoherence of initially pure states, measured by the increase of their von Neumann entropy, averaged over an ensemble of random pure states, is proved to be bounded from above by the partial entropy used to define the ALF dynamical entropy. The rate of decoherence induced by the sequence of the von Neumann projectors measurements is shown to be maximal, if the measurements are performed in a randomly chosen basis. The numerically observed linear increase of entropies is attributed to free-independence of the measured observable and the unitary dynamical map. △ Less

Submitted 26 September, 2003; originally announced September 2003.

Comments: 13 pages, 3 figures

Journal ref: J. Phys. A 37, 5157 (2004)

Entanglement of 2xK quantum systems

Authors: Artur Lozinski, Andreas Buchleitner, Karol Zyczkowski, Thomas Wellens

Abstract: We derive an analytical expression for the lower bound of the concurrence of mixed quantum states of composite 2xK systems. In contrast to other, implicitly defined entanglement measures, the numerical evaluation of our bound is straightforward. We explicitly evaluate its tightness for general mixed states of 2x3 systems, and identify a large class of states where our expression gives the exact… ▽ More We derive an analytical expression for the lower bound of the concurrence of mixed quantum states of composite 2xK systems. In contrast to other, implicitly defined entanglement measures, the numerical evaluation of our bound is straightforward. We explicitly evaluate its tightness for general mixed states of 2x3 systems, and identify a large class of states where our expression gives the exact value of the concurrence. △ Less

Submitted 19 February, 2003; originally announced February 2003.

Comments: 7 pages, 1 figure, to be published in Europhysics Letter

Journal ref: Europhys. Lett. 62, 168 (2003)

Quantum Iterated Function Systems

Authors: Artur Lozinski, Karol Zyczkowski, Wojciech Slomczynski

Abstract: Iterated functions system (IFS) is defined by specifying a set of functions in a classical phase space, which act randomly on an initial point. In an analogous way, we define a quantum iterated functions system (QIFS), where functions act randomly with prescribed probabilities in the Hilbert space. In a more general setting a QIFS consists of completely positive maps acting in the space of densi… ▽ More Iterated functions system (IFS) is defined by specifying a set of functions in a classical phase space, which act randomly on an initial point. In an analogous way, we define a quantum iterated functions system (QIFS), where functions act randomly with prescribed probabilities in the Hilbert space. In a more general setting a QIFS consists of completely positive maps acting in the space of density operators. We present exemplary classical IFSs, the invariant measure of which exhibits fractal structure, and study properties of the corresponding QIFSs and their invariant states. △ Less

Submitted 26 November, 2003; v1 submitted 4 October, 2002; originally announced October 2002.

Comments: 12 pages, 1 figure included

Journal ref: Phys. Rev. E 68, 046110 (2003)

Euroattractor: a brief introduction to Iterated Function Systems

Authors: Karol Zyczkowski, Artur Lozinski

Abstract: In this work we propose a definition of an Euroattractor: an attracting invariant measure of a certain iterated functions system (IFS). An IFS is defined by specifying a set of functions, defined in subsets of R^N or in a classical phase space, which act randomly on the initial point, so it may be considered as a generalization of the notion of classical dynamical system. If the functions are su… ▽ More In this work we propose a definition of an Euroattractor: an attracting invariant measure of a certain iterated functions system (IFS). An IFS is defined by specifying a set of functions, defined in subsets of R^N or in a classical phase space, which act randomly on the initial point, so it may be considered as a generalization of the notion of classical dynamical system. If the functions are sufficiently contracting, there exists an invariant measure of the system, often concentrated on a fractal set. We investigate invariant measures of a certain class of generalized or weakly contracting IFS's. △ Less

Submitted 30 October, 2002; originally announced October 2002.

Comments: 8 pages, 4 figures, to be published in Proceedings of the EUROATTRACTOR conference, Warsaw, June 2002

Journal ref: pages 231-241 in "Simplicity behind Complexity", (ed. W. Klonowski), Pabst Publishers, Langerich, 2004.

Irreversible Quantum Baker Map

Authors: Artur Lozinski, Prot Pakonski, Karol Zyczkowski

Abstract: We propose a generalization of the model of classical baker map on the torus, in which the images of two parts of the phase space do overlap. This transformation is irreversible and cannot be quantized by means of a unitary Floquet operator. A corresponding quantum system is constructed as a completely positive map acting in the space of density matrices. We investigate spectral properties of th… ▽ More We propose a generalization of the model of classical baker map on the torus, in which the images of two parts of the phase space do overlap. This transformation is irreversible and cannot be quantized by means of a unitary Floquet operator. A corresponding quantum system is constructed as a completely positive map acting in the space of density matrices. We investigate spectral properties of this super-operator and their link with the increase of the entropy of initially pure states. △ Less

Submitted 26 February, 2002; originally announced February 2002.

Comments: 4 pages, 3 figures included

Journal ref: Phys. Rev. E 66, 065201(2002)