Spectral Numerical Exterior Calculus Methods for Differential Equations on Radial Manifolds

We develop exterior calculus approaches for partial differential equations on radial manifolds. We introduce numerical methods that approximate with spectral accuracy the exterior derivative $\mathbf{d}$, Hodge star $\star$, and their compositions. To achieve discretizations with high precision and symmetry, we develop hyperinterpolation methods based on spherical harmonics and Lebedev quadrature. We perform convergence studies of our numerical exterior derivative operator $\overline{\mathbf{d}}$ and Hodge star operator $\overline{\star}$ showing each converge spectrally to $\mathbf{d}$ and $\star$. We show how the numerical operators can be naturally composed to formulate general numerical approximations for solving differential equations on manifolds. We present results for the Laplace-Beltrami equations demonstrating our approach.Comment: 22 pages, 13 figure

Parametric Deformation of Discrete Geometry for Aerodynamic Shape Design

We present a versatile discrete geometry manipulation platform for aerospace vehicle shape optimization. The platform is based on the geometry kernel of an open-source modeling tool called Blender and offers access to four parametric deformation techniques: lattice, cage-based, skeletal, and direct manipulation. Custom deformation methods are implemented as plugins, and the kernel is controlled through a scripting interface. Surface sensitivities are provided to support gradient-based optimization. The platform architecture allows the use of geometry pipelines, where multiple modelers are used in sequence, enabling manipulation difficult or impossible to achieve with a constructive modeler or deformer alone. We implement an intuitive custom deformation method in which a set of surface points serve as the design variables and user-specified constraints are intrinsically satisfied. We test our geometry platform on several design examples using an aerodynamic design framework based on Cartesian grids. We examine inverse airfoil design and shape matching and perform lift-constrained drag minimization on an airfoil with thickness constraints. A transport wing-fuselage integration problem demonstrates the approach in 3D. In a final example, our platform is pipelined with a constructive modeler to parabolically sweep a wingtip while applying a 1-G loading deformation across the wingspan. This work is an important first step towards the larger goal of leveraging the investment of the graphics industry to improve the state-of-the-art in aerospace geometry tools

Stable pair compactification of moduli of K3 surfaces of degree 2

We prove that the universal family of polarized K3 surfaces of degree 2 can be extended to a flat family of stable slc pairs $(X,\epsilon R)$ over the toroidal compactification associated to the Coxeter fan. One-parameter degenerations of K3 surfaces in this family are described by integral-affine structures on a sphere with 24 singularities.Comment: v2: Updated reference

Shape theory and mathematical design of a general geometric kernel through regular stratified objects

This thesis was submitted for the degree of Doctor of Philosophy and awarded by Brunel University.This dissertation focuses on the mathematical design of a unified shape kernel for geometric computing, with possible applications to computer aided design (CAM) and manufacturing (CAM), solid geometric modelling, free-form modelling of curves and surfaces, feature-based modelling, finite element meshing, computer animation, etc. The generality of such a unified shape kernel grounds on a shape theory for objects in some Euclidean space. Shape does not mean herein only geometry as usual in geometric modelling, but has been extended to other contexts, e. g. topology, homotopy, convexity theory, etc. This shape theory has enabled to make a shape analysis of the current geometric kernels. Significant deficiencies have been then identified in how these geometric kernels represent shapes from different applications. This thesis concludes that it is possible to construct a general shape kernel capable of representing and manipulating general specifications of shape for objects even in higher-dimensional Euclidean spaces, regardless whether such objects are implicitly or parametrically defined, they have ‘incomplete boundaries’ or not, they are structured with more or less detail or subcomplexes, which design sequence has been followed in a modelling session, etc. For this end, the basic constituents of such a general geometric kernel, say a combinatorial data structure and respective Euler operators for n-dimensional regular stratified objects, have been introduced and discussed

Surface parameterization over regular domains

Surface parameterization has been widely studied and it has been playing a critical role in many geometric processing tasks in graphics, computer-aided design, visualization, vision, physical simulation and etc. Regular domains, such as polycubes, are favored due to their structural regularity and geometric simplicity. This thesis focuses on studying the surface parameterization over regular domains, i.e. polycubes, and develops effective computation algorithms. Firstly, the motivation for surface parameterization and polycube mapping is introduced. Secondly, we briefly review existing surface parameterization techniques, especially for extensively studied parameterization algorithms for topological disk surfaces and parameterizations over regular domains for closed surfaces. Then we propose a polycube parameterization algorithm for closed surfaces with general topology. We develop an efficient optimization framework to minimize the angle and area distortion of the mapping. Its applications on surface meshing, inter-shape morphing and volumetric polycube mapping are also discussed

Formal study of plane Delaunay triangulation

This article presents the formal proof of correctness for a plane Delaunay triangulation algorithm. It consists in repeating a sequence of edge flippings from an initial triangulation until the Delaunay property is achieved. To describe triangulations, we rely on a combinatorial hypermap specification framework we have been developing for years. We embed hypermaps in the plane by attaching coordinates to elements in a consistent way. We then describe what are legal and illegal Delaunay edges and a flipping operation which we show preserves hypermap, triangulation, and embedding invariants. To prove the termination of the algorithm, we use a generic approach expressing that any non-cyclic relation is well-founded when working on a finite set