No Arabic abstract
The main purpose of this article is to develop a novel refinement strategy for four-dimensional hybrid meshes based on cubic pyramids. This optimal refinement strategy subdivides a given cubic pyramid into a conforming set of congruent cubic pyramids and invariant bipentatopes. The theoretical properties of the refinement strategy are rigorously analyzed and evaluated. In addition, a new class of fully symmetric quadrature rules with positive weights are generated for the cubic pyramid. These rules are capable of exactly integrating polynomials with degrees up to 12. Their effectiveness is successfully demonstrated on polynomial and transcendental functions. Broadly speaking, the refinement strategy and quadrature rules in this paper open new avenues for four-dimensional hybrid meshing, and space-time finite element methods.
The unsigned p-Willmore functional introduced in cite{mondino2011} generalizes important geometric functionals which measure the area and Willmore energy of immersed surfaces. Presently, techniques from cite{dziuk2008} are adapted to compute the first variation of this functional as a weak-form system of equations, which are subsequently used to develop a model for the p-Willmore flow of closed surfaces in $mathbb{R}^3$. This model is amenable to constraints on surface area and enclosed volume, and is shown to decrease the p-Willmore energy monotonically over time. In addition, a penalty-based regularization procedure is formulated to prevent artificial mesh degeneration along the flow; inspired by a conformality condition derived in cite{kamberov1996}, this procedure encourages angle-preservation in a closed and oriented surface immersion as it evolves. Following this, a finite-element discretization of both systems is discussed, and an application to mesh editing is presented.
We investigate moduli of planar circular quadrilaterals symmetric with respect to both the coordinate axes. First we develop an analytic approach which reduces this problem to ODEs and devise a numeric method to find out the accessory parameters. This method uses the Schwarz equation to determine conformal mapping of the unit disk onto a given circular quadrilateral. We also give an example of a circular quadrilateral for which the value of the conformal modulus can be found in the analytic form; this example is used to validate the numeric calculations. We also use another method, so called hpFEM, for the numeric calculation of the moduli. These two different approaches provide results agreeing with high accuracy.
We present a class of discretisation spaces and H(div)-conformal elements that can be built on any polytope. Bridging the flexibility of the Virtual Element spaces towards the elements shape with the divergence properties of the Raviart-Thomas elements on the boundaries, the designed frameworks offer a wide range of H(div)-conformal discretisations. As those elements are set up through degrees of freedom, their definitions are easily amenable to the properties the approximated quantities are wished to fulfil. Furthermore, we show that one straightforward restriction of this general setting share its properties with the classical Raviart-Thomas elements at each interface, for any order and any polytopial shape. Then, we investigate the shape of the basis functions corresponding to particular elements in the two dimensional case.
This paper presents an enhanced version of our previous work, hybrid non-uniform subdivision surfaces [19], to achieve optimal convergence rates in isogeometric analysis. We introduce a parameter $lambda$ ($frac{1}{4}<lambda<1$) to control the rate of shrinkage of irregular regions, so the method is called tuned hybrid non-uniform subdivision (tHNUS). Our previous work corresponds to the case when $lambda=frac{1}{2}$. While introducing $lambda$ in hybrid subdivision significantly complicates the theoretical proof of $G^1$ continuity around extraordinary vertices, reducing $lambda$ can recover the optimal convergence rates when tuned hybrid subdivision functions are used as a basis in isogeometric analysis. From the geometric point of view, the tHNUS retains comparable shape quality as [19] under non-uniform parameterization. Its basis functions are refinable and the geometric mapping stays invariant during refinement. Moreover, we prove that a tuned hybrid subdivision surface is globally $G^1$-continuous. From the analysis point of view, tHNUS basis functions form a non-negative partition of unity, are globally linearly independent, and their spline spaces are nested. We numerically demonstrate that tHNUS basis functions can achieve optimal convergence rates for the Poissons problem with non-uniform parameterization around extraordinary vertices.
We design a Hybrid High-Order (HHO) scheme for the Poisson problem that is fully robust on polytopal meshes in the presence of small edges/faces. We state general assumptions on the stabilisation terms involved in the scheme, under which optimal error estimates (in discrete and continuous energy norms, as well as $L^2$-norm) are established with multiplicative constants that do not depend on the maximum number of faces in each element, or the relative size between an element and its faces. We illustrate the error estimates through numerical simulations in 2D and 3D on meshes designed by agglomeration techniques (such meshes naturally have elements with a very large numbers of faces, and very small faces).