**Are you an EPFL student looking for a semester project?**

Work with us on data science and visualisation projects, and deploy your project as an app on top of GraphSearch.

Publication# An accuracy study of resistive methods for elliptic PDEs

Abstract

We numerically study the resistive method for the numerical approximation of elliptic PDEs. In particular we focus on the resistive method for weakly setting solution values in specific subdomains or interfaces in the domain.

Official source

This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.

Related concepts

Loading

Related publications

Loading

Related concepts (8)

Related publications (3)

Numerical analysis

Numerical analysis is the study of algorithms that use numerical approximation (as opposed to symbolic manipulations) for the problems of mathematical analysis (as distinguished from discrete mathema

Accuracy and precision

Accuracy and precision are two measures of observational error.
Accuracy is how close a given set of measurements (observations or readings) are to their true value, while precision is how close the

Buffer solution

A buffer solution (more precisely, pH buffer or hydrogen ion buffer) is an acid or a base aqueous solution consisting of a mixture of a weak acid and its conjugate base, or vice versa. Its pH change

Loading

Loading

Loading

Simone Deparis, Davide Forti, Antonello Gerbi, Alfio Quarteroni

We consider the numerical simulation of integrated heart models for the study of the cardiac functioning with particular emphasis on the coupling of the muscle contraction, the fluid dynamics of the left ventricle, and the interaction with the valves. We address the numerical approximation of the blood flows inside the left ventricle by using the Finite Element method for the spatial approximation of the fluid–structure interaction problem together with the aortic valve and the cardiac muscle, specifically when the contraction of the latter is driven by electromechanical models. We present and discuss numerical results for the fluid dynamics of the left ventricle by analyzing the role of different coupling strategies.

The goal of this project is to numerically solve the Navier-Stokes equations by using different numerical methods with particular emphasis on solving the problem of the flow past a square cylinder. In particular, we use the finite element method based on P2/P1 elements for the velocity and pressure fields for the spatial approximation, while the backward Euler method (with semi-implicit treatment of the nonlinear term) for the time discretization. We firstly test the numerical schemes on a benchmark problem with known exact solution. Then, we discuss in detail the advantages, in terms of computational costs, in using the algebraic Chorin-Temam method with additional implementation improvements. We finally investigate the problem of the two-dimensional flow past a square cylinder, focusing our attention on the range 0.1-300 for the Reynolds number (Re). We describe the two different regimes associated to the steady and the unsteady flows and we remark as the latter is in fact due to a Hopf bifurcation of the system. We also discuss the relation between the Strouhal and Reynolds numbers concluding that the Strouhal number attains its maximum value in the range 169-170 for the Reynolds number. In particular, a cubic model is proposed, showing very good matching with observed data and a better fitting than other models available in literature.

2013Claudia Maria Colciago, Simone Deparis

This paper deals with fast simulations of the hemodynamics in large arteries by considering a reduced model of the associated fluid-structure interaction problem, which in turn allows an additional reduction in terms of the numerical discretisation. The resulting method is both accurate and computationally cheap. This goal is achieved by means of two levels of reduction: first, we describe the model equations with a reduced mathematical formulation which allows to write the fluid-structure interaction problem as a Navier-Stokes system with non-standard boundary conditions; second, we employ numerical reduction techniques to further and drastically lower the computational costs. The non standard boundary condition is of a generalized Robin type, with a boundary mass and boundary stiffness terms accounting for the arterial wall compliance. The numerical reduction is obtained coupling two well-known techniques: the proper orthogonal decomposition and the reduced basis method, in particular the greedy algorithm. We start by reducing the numerical dimension of the problem at hand with a proper orthogonal decomposition and we measure the system energy with specific norms; this allows to take into account the different orders of magnitude of the state variables, the velocity and the pressure. Then, we introduce a strategy based on a greedy procedure which aims at enriching the reduced discretization space with low offline computational costs. As application, we consider a realistic hemodynamics problem with a perturbation in the boundary conditions and we show the good performances of the reduction techniques presented in the paper. The results obtained with the numerical reduction algorithm are compared with the one obtained by a standard finite element method. The gains obtained in term of CPU time are of three orders of magnitude.

2018