• In recent papers the author introduced a simple alternative to isoparametric finite elements of the n-simplex type, to enhance the accuracy of approximations of second-order boundary value problems with Dirichlet conditions, posed in smooth curved domains. This technique is based upon trial-functions consisting of piecewise polynomials defined on straight-edged triangular or tetrahedral meshes, interpolating the Dirichlet boundary conditions at points of the true boundary. In contrast the test-functions are defined upon the standard degrees of freedom associated with the underlying method for polytopic domains. While method's mathematical analysis for both second- and fourth-order problems in two-dimensional domains was carried out in arxiv NA-1701.00663 and in a submitted paper, this article is devoted to the study of the three-dimensional case, in which the method is nonconforming. Well-posedness, uniform stability and optimal a priori error estimates in the energy norm are demonstrated for a tetrahedron-based Lagrange family of finite elements. Novel L2-error estimates for the class of problems considered in this work are also proved. A series of numerical examples illustrates the potential of the new technique. In particular its better accuracy at equivalent cost as compared to the isoparametric technique is highlighted. Moreover the great generality of the new approach is exemplified through a method with degrees of freedom other than nodal values.
  • Since the 1960's the finite element method emerged as a powerful tool for the numerical simulation of countless physical phenomena or processes in applied sciences. One of the reasons for this undeniable success is the great versatility of the finite-element approach to deal with different types of geometries. This is particularly true of problems posed in curved domains of arbitrary shape. In the case of function-value Dirichlet conditions prescribed on curvilinear boundaries method's isoparametric version for meshes consisting of curved triangles or tetrahedra has been mostly employed to recover the optimal approximation properties known to hold for standard straight-edged elements in the case of polygonal or polyhedral domains. However, besides obvious algebraic and geometric inconveniences, the isoparametric technique is helplessly limited in scope and simplicity, since its extension to degrees of freedom other than function values is not straightforward if not unknown. The purpose of this paper is to propose, study and test a simple alternative that bypasses all the above drawbacks, without eroding qualitative approximation properties. More specifically this technique can do without curved elements and is based only on polynomial algebra. REMARKS : First submission (dated Jan. 3, 2017) updated on Jan. 11, 2017 with the addition of a footnote on author's research grant. A third version with several improvements was posted on Nov. 2, 2017. The fourth version incorporated some findings during the revision of a related paper. In the fifth version, besides minor changes, the convection-diffusion equation became the model problem; the text was reviewed in order to highlight the advantages of the new approach over classical techniques, more particularly by means of additional numerical examples.
  • This paper deals with a simple and straightforward procedure for automatic generation of finite-element or finite-volume meshes of spheroidal domains, consisting of tetrahedra. Besides the equation of the boundary, the generated meshes depend only on an integer parameter, whose value is associated with the degree of refinement. More specifically the procedure applies to the case where the boundary of a curved three-dimensional domain not so irregular can be expressed in spherical coordinates, with origin placed at a suitable location in its interior. An optimal numbering of mesh elements and nodes can be accomplished very easily. Several examples indicate that the generated meshes form a quasi-uniform family of partitions, as the corresponding value of the integer parameter increases, as long as the domain is not too distorted.
  • This work addresses techniques to solve convection-diffusion problems based on Hermite interpolation. We extend to the case of these equations a Hermite finite element method providing flux continuity across inter-element boundaries, shown to be a well-adapted tool for simulating pure diffusion phenomena (cf. V. Ruas, J. Comput. Appl. Maths., 246 p. 234-242, 2013). We consider two methods that can be viewed as non trivial improved versions of the lowest order Raviart-Thomas mixed method, corresponding to its extensions to convection-diffusion problems proposed by Douglas and Roberts (cf. Computational and Applied Mathematics, 1, p. 91-103, 1982) . A detailed convergence study is carried out for one of the methods, and numerical results illustrate the performance of both of them, as compared to each other and to the corresponding mixed methods.