Convergence of a stabilized discontinuous Galerkin method for incompressible nonlinear elasticity

In this paper we present a discontinuous Galerkin method applied to incompressible nonlinear elastostatics in a total Lagrangian deformation-pressure formulation, for which a suitable interior penalty stabilization is applied. We prove that the proposed discrete formulation for the linearized problem is well-posed, asymptotically consistent and that it converges to the corresponding weak solution. The derived convergence rates are optimal and further confirmed by a set of numerical examples in two and three spatial dimensions.

Object-oriented finite element programming

New technologies in computer science applied to numerical computations open the door to alternative approaches to mechanical problems using the finite element method. In classical approaches, theoretical developments often become cumbersome and the computer model which follows shows resemblance with the initial problem statement. The first step in the development consists usually in the analysis of the physics of the problem to simulate. The problem is generally described by a set of equations including partial differential equations. This first model is then replaced by successive equivalent or approximated models. The final result consists in a mathematical description of elemental matrices and algorithms describing the matrix form of the problem. The traditional approach consists then in constructing a computer model, generally complex and often quite different from the original mathematical description, thus making further corrections difficult. Therefore, the crucial problem of both the software architecture and the choice of the appropriate programming language is raised. Partially breaking with this approach, we propose a new approach to develop and program finite element formulations. The approach is based on a hybrid symbolic/numerical approach on the one hand, and on a high level software tool, object-oriented programming (supported here by the languages Smalltalk and C++) on the other hand. The aim of this work is to develop an appropriate environment for the algebraic manipulations needed for a finite element formulation applied to an initial boundary value problem, and also to perform efficient numerical computations. The new environment should make it possible to manage al1 the concepts necessary to solve a physical problem: manipulation of partial differential equations, variational formulations, integration by parts, weak forms, finite element approximations… The concepts manipulated therefore remain closely related to the original mathematical framework. The result of these symbolic manipulations is a set of elemental data (mass matrix, stiffness matrix, tangent stiffness matrix,…) to be introduced in a classical numerical code. The object-oriented paradigm is essential to the success of the implementation. In the context of the finite element codes, the object-oriented approach has already proved its capacity to represent and handle complex structures and phenomena. This is confirmed here with the symbolic environment for derivation of finite element formulations in which objects such as expression, integral and variational formulation appear. The link between both the numerical world and the symbolic world is based on an object-oriented concept for automatic programmation of matrix forms derived from the finite element method. As a result, a global environment in which the numerical is capable of evolving, using a language close to the natural mathematical one, is achieved. The potential of the approach is further demonstrated, on the one hand, by the wide range of problems solved in linear mechanics (electrodynamics in 1 and 2D, heat diffusion,…) as well as in nonlinear mechanics (advection dominated 1D problem, Navier Stokes problem), and, on the other hand by the diversity of the formulations manipulated (Galerkin formulations, space-time Galerkin formulations continuous in space and discontinuous in time, generalized Galerkin least-squares formulations).

EPFL1997

A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data

Fabio Nobile

We investigate the convergence rate of approximations by finite sums of rank-1 tensors of solutions of multiparametric elliptic PDEs. Such PDEs arise, for example, in the parametric, deterministic reformulation of elliptic PDEs with random field inputs, based, for example, on the M-term truncated Karhunen–Lo`eve expansion. Our approach could be regarded as either a class of compressed approximations of these solutions or as a new class of iterative elliptic problem solvers for high-dimensional, parametric, elliptic PDEs providing linear scaling complexity in the dimension M of the parameter space. It is based on rank-reduced, tensor-formatted separable approximations of the high-dimensional tensors and matrices involved in the iterative process, combined with the use of spectrally equivalent low-rank tensor-structured preconditioners to the parametric matrices resulting from a finite element discretization of the high-dimensional parametric, deterministic problems. Numerical illustrations for the M-dimensional parametric elliptic PDEs resulting from sPDEs on parameter spaces of dimensions M ≤ 100 indicate the advantages of employing low-rank tensor-structured matrix formats in the numerical solution of such problems.

Society for Industrial and Applied Mathematics2010

On the coupling of 1D and 3D Diffusion-Reaction Equations: Application to Tissue Pefusion Problems

Alfio Quarteroni

In this paper we consider the coupling between two diffusion-reaction problems, one taking place in a three-dimensional domain Ω, the other in a one-dimensional subdomain Λ. This coupled problem is the simplest model of fluid flow in a three-dimensional porous medium featuring fractures that can be described by one-dimensional manifolds. In particular this model can provide the basis for a multiscale analysis of blood flow through tissues, in which the capillary network is represented as a porous matrix, while the major blood vessels are described by thin tubular structures embedded into it: in this case, the model allows the computation of the 3D and 1D blood pressures respectively in the tissue and in the vessels. The mathematical analysis of the problem requires non-standard tools, since the mass conservation condition at the interface between the porous medium and the one-dimensional manifold has to be taken into account by means of a measure term in the 3D equation. In particular, the 3D solution is singular on Λ. In this work, suitable weighted Sobolev spaces are introduced to handle this singularity: the well-posedness of the coupled problem is established in the proposed functional setting. An advantage of such an approach is that it provides a Hilbertian framework which may be used for the convergence analysis of finite element approximation schemes. The investigation of the numerical approximation will be the subject of a forthcoming work.

2008

Object-oriented finite element programming

EPFL1997