We analyze two numerical schemes of Euler type in time and C 0finite-element type with $\mathbb{P}_1$ -approximation in space forsolving a phase-field model of a binary alloy with thermalproperties. This model is written as a highly non-linear parabolicsystem with three unknowns: phase-field, solute concentration andtemperature, where the diffusion for the temperature and soluteconcentration may degenerate. The first scheme is nonlinear, unconditionally stableand convergent. The other scheme is linear but conditionally stableand convergent. A maximum principle is avoided in both schemes,using a truncation operator on the L 2 projection onto the $\mathbb{P}_0$ finite element for the discrete concentration. Inaddition, for the model when the heat conductivity and solutediffusion coefficients are constants, optimal error estimates forboth schemes are shown based on stability estimates.

We developed a mimetic finite difference method for solving elliptic equations with tensor coefficients on polyhedral meshes. The first-order convergence estimates in a mesh-dependent H1 norm are derived.

We propose and analyze numerical schemes for viscosity solutions of time-dependent Hamilton-Jacobi equations on the Heisenberg group.The main idea is to construct a grid compatible with the noncommutative group geometry. Under suitable assumptions on the data, the Hamiltonian and the parameters for the discrete first order scheme,we prove that the error between the viscosity solution computed at the grid nodes and the solution of the discrete problem behaves like $\sqrt{h}$ where h is the mesh step. Such an estimate is similar to those available in the Euclidean geometrical setting.The theoretical results are tested numerically on some examples for which semi-analytical formulas for the computation of geodesics are known. Other simulations are presented, for both steady and unsteady problems.

In this paper, by use of affine biquadratic elements, we constructand analyze a finite volume element scheme for elliptic equations onquadrilateral meshes. The scheme is shown to be of second-order inH 1-norm, provided that each quadrilateral in partition is almosta parallelogram. Numerical experiments are presented to confirm theusefulness and efficiency of the method.

In order to handle the flow of a viscous incompressible fluid in a porous medium with cracks, the thickness of which cannot be neglected, we consider a model which couples the Darcy equations in the medium with the Stokes equations in the cracks by a new boundary condition at the interface, namely the continuity of the pressure. We prove that this model admits a unique solution and propose a mixed formulation of it. Relying on this formulation, we describe a finite element discretization and derive a priori and a posteriori error estimates. We present some numerical experiments that are in good agreement with the analysis.

This paper is concerned with the finite volume approximation of the p-Laplacian equation with homogeneous Dirichlet boundary conditions on rectangular meshes.A reconstruction of the norm of the gradient on the mesh's interfaces is needed in order to discretize the p-Laplacian operator. We give a detailed description of the possible nine points schemes ensuring that the solution of the resulting finite dimensional nonlinear system exists and is unique. These schemes, called admissible, are locally conservative and in addition derive from the minimization of a strictly convexe and coercive discrete functional.The convergence rate is analyzed when the solution lies in W2,p . Numerical results are given in order to compare different admissible and non-admissible schemes.

This paper is concerned with the computation of 3D vertex singularities of anisotropic elastic fields with Dirichlet boundary conditions, focusing on the derivation of error estimates for a finite element method on graded meshes. The singularities are described by eigenpairs of a corresponding operator pencil on spherical polygonal domains. The main idea is to introduce a modified quadratic variational boundary eigenvalue problem which consists of two self-adjoint, positive definite sesquilinear forms and a skew-Hermitean form. This eigenvalue problem is discretized by a finite element method on graded meshes. Based on regularity results for the eigensolutions estimates for the finite element error are derived both for the eigenvalues and the eigensolutions. Finally, some numerical results are presented.

The goal of this paper is to derive some error estimates for thenumerical discretization of some optimal control problems governedby semilinear elliptic equations with bound constraints on thecontrol and a finitely number of equality and inequality stateconstraints. We prove some error estimates for the optimalcontrols in the L ∞ norm and we also obtain error estimatesfor the Lagrange multipliers associated to the state constraintsas well as for the optimal states and optimal adjoint states.

In this paper, a phase field system of Penrose–Fife type withnon–conserved order parameter is considered. A class of time–discrete schemes for an initial–boundary value problem for this phase–field system is presented.In three space dimensions, convergence is proved and an error estimate linear with respect to the time–step sizeis derived.

We consider a finite element discretization bythe Taylor–Hood element for the stationaryStokes and Navier–Stokesequations with slip boundary condition. The slip boundary conditionis enforced pointwise for nodal values of the velocity inboundary nodes. We prove optimal error estimates in theH 1 and L 2 norms for the velocity and pressure respectively.