skip to main content


Title: Quadrilateral Mesh Generation III : Optimizing Singularity Configuration Based on Abel-Jacobi Theory
This work proposes a rigorous and practical algorithm for quad-mesh generation based the Abel-Jacobi theory of algebraic \textcolor{red}{curves}. We prove sufficient and necessary conditions for a flat metric with cone singularities to be compatible with a quad-mesh, in terms of the deck-transformation, then develop an algorithm based on the theorem. The algorithm has two stages: first, a meromorphic quartic differential is generated to induce a T-mesh; second, the edge lengths of the T-mesh are adjusted by solving a linear system to satisfy the deck transformation condition, which produces a quad-mesh. In the first stage, the algorithm pipeline can be summarized as follows: calculate the homology group; compute the holomorphic differential group; construct the period matrix of the surface and Jacobi variety; calculate the Abel-Jacobi map for a given divisor; optimize the divisor to satisfy the Abel-Jacobi condition by integer programming; compute \textcolor{red}{a} flat Riemannian metric with cone singularities at the divisor by Ricci flow; \textcolor{red}{isometrically} immerse the surface punctured at the divisor onto the complex plane and pull back the canonical holomorphic differential to the surface to obtain the meromorphic quartic differential; construct a motorcycle graph to generate a T-Mesh. In the second stage, the deck transformation constraints are formulated as a linear equation system of the edge lengths of the T-mesh. The solution provides a flat metric with integral deck transformations, which leads to the final quad-mesh. The proposed method is rigorous and practical. The T-mesh and quad-mesh results can be applied for constructing Splines directly. The efficiency and efficacy of the proposed algorithm are demonstrated by experimental results on surfaces with complicated topologies and geometries.  more » « less
Award ID(s):
1762287
NSF-PAR ID:
10291675
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Computer methods in applied mechanics and engineering
ISSN:
1879-2138
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    This work discovers the equivalence relation between quadrilateral meshes and meromorphic quartic differentials. Each quad-mesh induces a conformal structure of the surface, and a meromorphic quartic differential, where the configuration of singular vertices corresponds to the configurations of the poles and zeros (divisor) of the meromorphic differential. Due to Riemann surface theory, the configuration of singularities of a quad-mesh satisfies the Abel–Jacobi condition. Inversely, if a divisor satisfies the Abel–Jacobi condition, then there exists a meromorphic quartic differential whose divisor equals the given one. Furthermore, if the meromorphic quartic differential is with finite trajectories, then it also induces a quad-mesh, the poles and zeros of the meromorphic differential correspond to the singular vertices of the quad-mesh. Besides the theoretic proofs, the computational algorithm for verification of Abel–Jacobi condition is also explained in detail. Furthermore, constructive algorithm of meromorphic quartic differential on genus zero surfaces is proposed, which is based on the global algebraic representation of meromorphic differentials. Our experimental results demonstrate the efficiency and efficacy of the algorithm. This opens up a novel direction for quad-mesh generation using algebraic geometric approach. 
    more » « less
  2. A spherical conical metric g g on a surface Σ \Sigma is a metric of constant curvature 1 1 with finitely many isolated conical singularities. The uniformization problem for such metrics remains largely open when at least one of the cone angles exceeds 2 π 2\pi . The eigenfunctions of the Friedrichs Laplacian Δ g \Delta _g with eigenvalue λ = 2 \lambda =2 play a special role in this problem, as they represent local obstructions to deformations of the metric g g in the class of spherical conical metrics. In the present paper we apply the theory of multivalued harmonic maps to spheres to the question of existence of such eigenfunctions. In the first part we establish a new criterion for the existence of 2 2 -eigenfunctions, given in terms of a certain meromorphic data on Σ \Sigma . As an application we give a description of all 2 2 -eigenfunctions for metrics on the sphere with at most three conical singularities. The second part is an algebraic construction of metrics with large number of 2 2 -eigenfunctions via the deformation of multivalued harmonic maps. We provide new explicit examples of metrics with many 2 2 -eigenfunctions via both approaches, and describe the general algorithm to find metrics with arbitrarily large number of 2 2 -eigenfunctions. 
    more » « less
  3. A bstract We study Euclidean D3-branes wrapping divisors D in Calabi-Yau orientifold compactifications of type IIB string theory. Witten’s counting of fermion zero modes in terms of the cohomology of the structure sheaf $$ {\mathcal{O}}_D $$ O D applies when D is smooth, but we argue that effective divisors of Calabi-Yau threefolds typically have singularities along rational curves. We generalize the counting of fermion zero modes to such singular divisors, in terms of the cohomology of the structure sheaf $$ {\mathcal{O}}_{\overline{D}} $$ O D ¯ of the normalization $$ \overline{D} $$ D ¯ of D . We establish this by detailing compactifications in which the singularities can be unwound by passing through flop transitions, giving a physical incarnation of the normalization process. Analytically continuing the superpotential through the flops, we find that singular divisors whose normalizations are rigid can contribute to the superpotential: specifically, $$ {h}_{+}^{\bullet}\left({\mathcal{O}}_{\overline{D}}\right)=\left(1,0,0\right) $$ h + • O D ¯ = 1 0 0 and $$ {h}_{-}^{\bullet}\left({\mathcal{O}}_{\overline{D}}\right)=\left(0,0,0\right) $$ h − • O D ¯ = 0 0 0 give a sufficient condition for a contribution. The examples that we present feature infinitely many isomorphic geometric phases, with corresponding infinite-order monodromy groups Γ. We use the action of Γ on effective divisors to determine the exact effective cones, which have infinitely many generators. The resulting nonperturbative superpotentials are Jacobi theta functions, whose modular symmetries suggest the existence of strong-weak coupling dualities involving inversion of divisor volumes. 
    more » « less
  4. Abstract

    This paper presents a fast sweeping method (FSM) to calculate the first‐arrival traveltimes of the qP, qSV, and qSH waves in two‐dimensional (2D) transversely isotropic media, whose symmetry axis may have an arbitrary orientation (tilted transverse isotropy [TTI]). The method discretizes the anisotropic eikonal equation with finite difference approximations on a rectangular mesh and solves the discretized system iteratively with the Gauss‐Seidel iterations along alternating sweeping orderings. At each mesh point, a highly nonlinear equation is solved to update the numerical solution until its convergence. For solving the nonlinear equation, an interval that contains the solutions is first determined and partitioned into few subintervals such that each subinterval contains one solution; then, the false position method is applied on these subintervals to compute the solutions; after that, among all possible solutions for the discretized equation, a causality condition is imposed, and the minimum solution satisfying the causality condition is chosen to update the solution. For problems with a point‐source condition, the FSM is extended for solving the anisotropic eikonal equation after a factorization technique is applied to resolve the source singularities, which yields clean first‐order accuracy. When dealing with the triplication of the qSV wave, solutions corresponding to the minimal group velocity are chosen such that continuous solutions are computed. The accuracy, efficiency, and capability of the proposed method are demonstrated with numerical experiments.

     
    more » « less
  5. Recent advancements in observational techniques have led to new tests of the general relativistic predictions for black-hole spacetimes in the strong-field regime. One of the key ingredients for several tests is a metric that allows for deviations from the Kerr solution but remains free of pathologies outside its event horizon. Existing metrics that have been used in the literature often do not satisfy the null convergence condition that is necessary to apply the strong rigidity theorem and would have allowed us to calculate the location of the event horizon by identifying it with an appropriate Killing horizon. This has led earlier calculations of event horizons of parametrically deformed metrics to either follow numerical techniques or simply search heuristically for coordinate singularities. We show that several of these metrics, almost by construction, are circular. We can, therefore, use the weak rigidity and Carter’s rotosurface theorem and calculate algebraically the locations of their event horizons, without relying on expansions or numerical techniques. We apply this approach to a number of parametrically deformed metrics, calculate the locations of their event horizons, and place constraints on the deviation parameters such that the metrics remain regular outside their horizons. We show that calculating the angular velocity of the horizon and the effective gravity there offers new insights into the observational signatures of deformed metrics, such as the sizes and shapes of the predicted black-hole shadows. 
    more » « less