Source author record

Gerard Awanou

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

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

16 published item(s)

preprint2021arXiv

On weak convergence of Monge-Ampere measures for discrete convex mesh functions

To a mesh function we associate the natural analogue of the Monge-Ampere measure. The latter is shown to be equivalent to the Monge-Ampere measure of the convex envelope. We prove that the uniform convergence to a bounded convex function of mesh functions implies the uniform convergence on compact subsets of their convex envelopes and hence the weak convergence of the associated Monge-Ampere measures. We also give conditions for mesh functions to have a subsequence which converges uniformly to a convex function. Our result can be used to give alternate proofs of the convergence of some discretizations for the second boundary value problem for the Monge-Ampere equation and was used for a recently proposed discretization of the latter. For mesh functions which are uniformly bounded and satisfy a convexity condition at the discrete level, we show that there is a subsequence which converges uniformly on compact subsets to a convex function. The convex envelopes of the mesh functions of the subsequence also converge uniformly on compact subsets. If in addition they agree with a continuous convex function on the boundary, the limit function is shown to satisfy the boundary condition strongly.

preprint2015arXiv

A variational method for computing numerical solutions of the Monge-Ampere equation

We present a numerical method for solving the Monge-Ampere equation based on the characterization of the solution of the Dirichlet problem as the minimizer of a convex functional of the gradient and under convexity and nonlinear constraints. When the equation is discretized with a certain monotone scheme, we prove that the unique minimizer of the discrete problem solves the finite difference equation. For the numerical results we use both the standard finite difference discretization and the monotone scheme. Results with standard tests confirm that the numerical approximations converge to the Aleksandrov solution.

preprint2015arXiv

On standard finite difference discretizations of the elliptic Monge-Ampere equation

Given an orthogonal lattice with mesh length h on a bounded convex domain, we propose to approximate the Aleksandrov solution of the Monge-Ampere equation by regularizing the data and discretizing the equation in a subdomain using the standard finite difference method. The Dirichlet data is used to approximate the solution in the remaining part of the domain. We prove the uniform convergence on compact subsets of the solution of the discrete problems to an approximate problem on the subdomain. The result explains the behavior of methods based on the standard finite difference method and designed to numerically converge to non-smooth solutions. We propose an algorithm which for smooth solutions appears faster than the popular Newton's method with a high accuracy for non smooth solutions. The convergence of the algorithm is independent of how close to the numerical solution the initial guess is, upon rescaling the equation and given a user's measure of the closeness of an initial guess.

preprint2015arXiv

Standard finite elements for the numerical resolution of the elliptic Monge-Ampere equation: Aleksandrov solutions

We prove a convergence result for a natural discretization of the Dirichlet problem of the elliptic Monge-Ampere equation using finite dimensional spaces of piecewise polynomial C0 or C1 functions. Standard discretizations of the type considered in this paper have been previous analyzed in the case the equation has a smooth solution and numerous numerical evidence of convergence were given in the case of non smooth solutions. Our convergence result is valid for non smooth solutions, is given in the setting of Aleksandrov solutions, and consists in discretizing the equation in a subdomain with the boundary data used as an approximation of the solution in the remaining part of the domain. Our result gives a theoretical validation for the use of a non monotone finite element method for the Monge-Ampère equation.

preprint2015arXiv

Standard finite elements for the numerical resolution of the elliptic Monge-Ampere equation: classical solutions

We propose a new variational formulation of the elliptic Monge-Ampere equation and show how classical Lagrange elements can be used for the numerical resolution of classical solutions of the equation. Error estimates are given for Lagrange elements of degree d >= 2 in dimensions 2 and 3. No jump term is used in the variational formulation. We propose to solve the discrete nonlinear system of equations by a time marching method and numerical evidence is given which indicates that one approximates weak solutions in two dimensions.

preprint2014arXiv

Mixed finite elements for elasticity on quadrilateral meshes

We present stable mixed finite elements for planar linear elasticity on general quadrilateral meshes. The symmetry of the stress tensor is imposed weakly and so there are three primary variables, the stress tensor, the displacement vector field, and the scalar rotation. We develop and analyze a stable family of methods, indexed by an integer $r \geq 2$ and with rate of convergence in the $L^2$ norm of order $r$ for all the variables. The methods use Raviart-Thomas elements for the stress, piecewise tensor product polynomials for the displacement, and piecewise polynomials for the rotation. We also present a simple first order element, not belonging to this family. It uses the lowest order BDM elements for the stress, and piecewise constants for the displacement and rotation, and achieves first order convergence for all three variables.

preprint2014arXiv

Pseudo transient continuation and time marching methods for Monge-Ampere type equations

We present two numerical methods for the fully nonlinear elliptic Monge-Ampere equation. The first is a pseudo transient continuation method and the second is a pure pseudo time marching method. The methods are proven to converge to a strictly convex solution of a natural discrete variational formulation with $C^1$ conforming approximations. The assumption of existence of a strictly convex solution to the discrete problem is proven for smooth solutions of the continuous problem and supported by numerical evidence for non smooth solutions.

preprint2013arXiv

Nonconforming tetrahedral mixed finite elements for elasticity

This paper presents a nonconforming finite element approximation of the space of symmetric tensors with square integrable divergence, on tetrahedral meshes. Used for stress approximation together with the full space of piecewise linear vector fields for displacement, this gives a stable mixed finite element method which is shown to be linearly convergent for both the stress and displacement, and which is significantly simpler than any stable conforming mixed finite element method. The method may be viewed as the three-dimensional analogue of a previously developed element in two dimensions. As in that case, a variant of the method is proposed as well, in which the displacement approximation is reduced to piecewise rigid motions and the stress space is reduced accordingly, but the linear convergence is retained.

preprint2013arXiv

Spline element method for the Monge-Ampere equation

We analyze the convergence of an iterative method for solving the nonlinear system resulting from a natural discretization of the Monge-Ampère equation with $C^1$ conforming approximations. We make the assumption, supported by numerical experiments for the two dimensional problem, that the discrete problem has a convex solution. The method we analyze is the discrete version of Newton's method in the vanishing moment methodology. Numerical experiments are given in the framework of the spline element method.

preprint2011arXiv

Rectangular Mixed Elements for Elasticity with Weakly Imposed symmetry Condition

We present new rectangular mixed finite elements for linear elasticity. The approach is based on a modification of the Hellinger-Reissner functional in which the symmetry of the stress field is enforced weakly through the introduction of a Lagrange multiplier. The elements are analogues of the lowest order elements described in Arnold, Falk and Winther [ Mixed finite element methods for linear elasticity with weakly imposed symmetry. Mathematics of Computation 76 (2007), pp. 1699--1723]. Piecewise constants are used to approximate the displacement and the rotation. The first order BDM elements are used to approximate each row of the stress field.

preprint2011arXiv

The serendipity family of finite elements

We give a new, simple, dimension-independent definition of the serendipity finite element family. The shape functions are the span of all monomials which are linear in at least s-r of the variables where s is the degree of the monomial or, equivalently, whose superlinear degree (total degree with respect to variables entering at least quadratically) is at most r. The degrees of freedom are given by moments of degree at most r-2d on each face of dimension d. We establish unisolvence and a geometric decomposition of the space.

preprint2011arXiv

Two remarks on rectangular mixed finite elements for elasticity

The lowest order nonconforming rectangular element in three dimensions involves 54 degrees of freedom for the stress and 12 degrees of freedom for the displacement. With a modest increase in the number of degrees of freedom (24 for the stress), we obtain a conforming rectangular element for linear elasticity in three dimensions. Moreover, unlike the conforming plane rectangular or simplicial elements, this element does not involve any vertex degrees of freedom. Second, we remark that further low order elements can be constructed by approximating the displacement with rigid body motions. This results in a pair of conforming elements with 72 degrees of freedom for the stress and 6 degrees of freedom for the displacement.

preprint2007arXiv

Finite elements for symmetric tensors in three dimension

We construct finite element subspaces of the space of symmetric tensors with square-integrable divergence on a three-dimensional domain. These spaces can be used to approximate the stress field in the classical Hellinger--Reissner mixed formulation of the elasticty equations, when standard discontinous finite element spaces are used to approximate the displacement field. These finite element spaces are defined with respect to an arbitrary simplicial triangulation of the domain, and there is one for each positive value of the polynomial degree used for the displacements. For each degree, these provide a stable finite element discretization. The construction of the spaces is closely tied to discretizations of the elasticity complex, and can be viewed as the three-dimensional analogue of the triangular element family for plane elasticity previously proposed by Arnold and Winther.