
In recent papers the author introduced a simple alternative to isoparametric
finite elements of the nsimplex type, to enhance the accuracy of
approximations of secondorder boundary value problems with Dirichlet
conditions, posed in smooth curved domains. This technique is based upon
trialfunctions consisting of piecewise polynomials defined on straightedged
triangular or tetrahedral meshes, interpolating the Dirichlet boundary
conditions at points of the true boundary. In contrast the testfunctions 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 fourthorder problems in twodimensional domains was carried out in
arxiv NA1701.00663 and in a submitted paper, this article is devoted to the
study of the threedimensional case, in which the method is nonconforming.
Wellposedness, uniform stability and optimal a priori error estimates in the
energy norm are demonstrated for a tetrahedronbased Lagrange family of finite
elements. Novel L2error 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 finiteelement 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 functionvalue 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 straightedged
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
convectiondiffusion 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 finiteelement or finitevolume 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 threedimensional 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 quasiuniform 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 convectiondiffusion problems based
on Hermite interpolation. We extend to the case of these equations a Hermite
finite element method providing flux continuity across interelement
boundaries, shown to be a welladapted tool for simulating pure diffusion
phenomena (cf. V. Ruas, J. Comput. Appl. Maths., 246 p. 234242, 2013). We
consider two methods that can be viewed as non trivial improved versions of the
lowest order RaviartThomas mixed method, corresponding to its extensions to
convectiondiffusion problems proposed by Douglas and Roberts (cf.
Computational and Applied Mathematics, 1, p. 91103, 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.