A Discrete Laplace–Beltrami Operator for Simplicial Surfaces

We define a discrete Laplace–Beltrami operator for simplicial surfaces (Definition 16). It depends only on the intrinsic geometry of the surface and its edge weights are positive. Our Laplace operator is similar to the well known finite-elements Laplacian (the so called “cotan formula”) except that it is based on the intrinsic Delaunay triangulation of the simplicial surface. This leads to new definitions of discrete harmonic functions, discrete mean curvature, and discrete minimal surfaces. The definition of the discrete Laplace–Beltrami operator depends on the existence and uniqueness of Delaunay tessellations in piecewise flat surfaces. While the existence is known, we prove the uniqueness. Using Rippa’s Theorem we show that, as claimed, Musin’s harmonic index provides an optimality criterion for Delaunay triangulations, and this can be used to prove that the edge flipping algorithm terminates also in the setting of piecewise flat surfaces. Research for this article was supported by the DFG Research Unit 565 “Polyhedral Surfaces” and the DFG Research Center Matheon “Mathematics for key technologies” in Berlin.     

The focal geometry of circular and conical meshes

Circular meshes are quadrilateral meshes all of whose faces possess a circumcircle, whereas conical meshes are planar quadrilateral meshes where the faces which meet in a vertex are tangent to a right circular cone. Both are amenable to geometric modeling – recently surface approximation and subdivision-like refinement processes have been studied. In this paper we extend the original defining property of conical meshes, namely the existence of face/face offset meshes at constant distance, to circular meshes. We study the close relation between circular and conical meshes, their vertex/vertex and face/face offsets, as well as their discrete normals and focal meshes. In particular we show how to construct a two-parameter family of circular (resp., conical) meshes from a given conical (resp., circular) mesh. We further discuss meshes which have both properties and their relation to discrete surfaces of negative Gaussian curvature. The offset properties of special quadrilateral meshes and the three-dimensional support structures derived from them are highly relevant for computational architectural design of freeform structures. Another aspect important for design is that both circular and conical meshes provide a discretization of the principal curvature lines of a smooth surface, so the mesh polylines represent principal features of the surface described by the mesh.      

Variational Laplacians for semidiscrete surfaces

We study a Laplace operator on semidiscrete surfaces that is defined by variation of the Dirichlet energy functional. We show existence and its relation to the mean curvature normal, which is itself defined via variation of area. We establish several core properties like linear precision (closely related to the mean curvature of flat surfaces), and pointwise convergence. It is interesting to observe how a certain freedom in choosing area measures yields different kinds of Laplacians: it turns out that using as a measure a simple numerical integration rule yields a Laplacian previously studied as the pointwise limit of geometrically meaningful Laplacians on polygonal meshes.     

基于Maxwell模型的线性粘弹性波传导问题的半离散/全离散杂交应力有限元法

1引言设ΩСR²为一个凸的有界开集,边界为ЭΩ;T为一个正常数.我们考虑如下基于Maxwell模型的二维粘弹性固体介质波传导问题。     

Application of a 14-point averaging operator in the grid method

The Dirichlet problem for Laplace’s equation in a rectangular parallelepiped is solved by applying the grid method. A 14-point averaging operator is used to specify the grid equations on the entire grid introduced in the parallelepiped. Given boundary values that are continuous on the parallelepiped edges and have first derivatives satisfying the Lipschitz condition on each parallelepiped face, the resulting discrete solution of the Dirichlet problem converges uniformly and quadratically with respect to the mesh size. Assuming that the boundary values on the faces have fourth derivatives satisfying the Hölder condition and the second derivatives on the edges obey an additional compatibility condition implied by Laplace’s equation, the discrete solution has uniform and quartic convergence with respect to the mesh size. The convergence of the method is also analyzed in certain cases when the boundary values are of intermediate smoothness.     

Numerical analysis for a locally damped wave equation

We consider a semi-discrete finite element formulation with artificial viscosity for the numerical approximation of a problem that models the damped vibrations of a string with fixed ends. The damping coefficient depends on the spatial variable and is effective only in a sub-interval of the domain. For this scheme, the energy of semi-discrete solutions decays exponentially and uniformly with respect to the mesh parameter to zero. We also introduce an implicit in time discretization. Error estimates for the semi-discrete and fully discrete schemes in the energy norm are provided and numerical experiments performed.     

Discrete Moving Frames on Lattice Varieties and Lattice-Based Multispaces

In this paper, we develop the theory of the discrete moving frame in two different ways. In the first half of the paper, we consider a discrete moving frame defined on a lattice variety and the equivalence classes of global syzygies that result from the first fundamental group of the variety. In the second half, we consider the continuum limit of discrete moving frames as a local lattice coalesces to a point. To achieve a well-defined limit of discrete frames, we construct multispace, a generalisation of the jet bundle that also generalises Olver’s one-dimensional construction. Using interpolation to provide coordinates, we prove that it is a manifold containing the usual jet bundle as a submanifold. We show that continuity of a multispace moving frame ensures that the discrete moving frame converges to a continuous one as lattices coalesce. The smooth frame is, at the same time, the restriction of the multispace frame to the embedded jet bundle. We prove further that the discrete invariants and syzygies approximate their smooth counterparts. In effect, a frame on multispace allows smooth frames and their discretisations to be studied simultaneously. In our last chapter we discuss two important applications, one to the discrete variational calculus, and the second to discrete integrable systems. Finally, in an appendix, we discuss a more general result concerning equicontinuous families of discretisations of moving frames, which are consistent with a smooth frame.     

An operator splitting method for nonlinear convection-diffusion equations

Summary. We present a semi-discrete method for constructing approximate solutions to the initial value problem for the -dimensional convection-diffusion equation . The method is based on the use of operator splitting to isolate the convection part and the diffusion part of the equation. In the case , dimensional splitting is used to reduce the -dimensional convection problem to a series of one-dimensional problems. We show that the method produces a compact sequence of approximate solutions which converges to the exact solution. Finally, a fully discrete method is analyzed, and demonstrated in the case of one and two space dimensions. ReceivedFebruary 1, 1996 / Revised version received June 24, 1996     

Biharmonic Green Functions on Homogeneous Trees

The study of biharmonic functions under the ordinary (Euclidean) Laplace operator on the open unit disk \mathbbD{\mathbb{D}} in \mathbbC{\mathbb{C}} arises in connection with plate theory, and in particular, with the biharmonic Green functions which measure, subject to various boundary conditions, the deflection at one point due to a load placed at another point. A homogeneous tree T is widely considered as a discrete analogue of the unit disk endowed with the Poincaré metric. The usual Laplace operator on T corresponds to the hyperbolic Laplacian. In this work, we consider a bounded metric on T for which T is relatively compact and use it to define a flat Laplacian which plays the same role as the ordinary Laplace operator on \mathbbD{\mathbb{D}}. We then study the simply-supported and the clamped biharmonic Green functions with respect to both Laplacians.     

Extremal eigenvalues of the Laplacian on Euclidean domains and closed surfaces

We investigate properties of the sequences of extremal values that could be achieved by the eigenvalues of the Laplacian on Euclidean domains of unit volume, under Dirichlet and Neumann boundary conditions, respectively. In a second part, we study sequences of extremal eigenvalues of the Laplace–Beltrami operator on closed surfaces of unit area.     

Carleman estimates for semi-discrete parabolic operators and application to the controllability of semi-linear semi-discrete parabolic equations

In arbitrary dimension, in the discrete setting of finite-differences we prove a Carleman estimate for a semi-discrete parabolic operator, in which the large parameter is connected to the mesh size. This estimate is applied for the derivation of a (relaxed) observability estimate, that yield some controlability results for semi-linear semi-discrete parabolic equations. Sub-linear and super-linear cases are considered.     

Uniform boundary controllability of a semi-discrete 1-D wave equation

Summary. A numerical scheme for the controlled semi-discrete 1-D wave equation is considered. We analyze the convergence of the boundary controls of the semi-discrete equations to a control of the continuous wave equation when the mesh size tends to zero. We prove that, if the high modes of the discrete initial data have been filtered out, there exists a sequence of uniformly bounded controls and any weak limit of this sequence is a control for the continuous problem. The number of the eliminated frequencies depends on the mesh size and the regularity of the continuous initial data. The case of the HUM controls is also discussed. Received March 3, 2001 / Published online October 17, 2001     

Semi-Discrete Isothermic Surfaces

We study mappings of the form ${x : \mathbb{Z}\times\mathbb{R}\to\mathbb{R}^3}$ which can be seen as a limit case of purely discrete surfaces, or as a semi-discretization of smooth surfaces. In particular we discuss circular surfaces, isothermic surfaces, conformal mappings, and dualizability in the sense of Christoffel. We arrive at a semi-discrete version of Koenigs nets and show that in the setting of circular surfaces, isothermicity is the same as dualizability. We show that minimal surfaces constructed as a dual of a sphere have vanishing mean curvature in a certain well-defined sense, and we also give an incidence-geometric characterization of isothermic surfaces.     

A Thomson’s Principle and a Rayleigh’s Monotonicity Law for Nonlinear Networks

We give a representation of the solution of the Neumann problem for the Laplace operator on the n-dimensional unit ball in terms of the solution of an associated Dirichlet problem. The representation is extended to other operators besides the Laplacian, to smooth simply connected planar domains, and to the infinite-dimensional Laplacian on the unit ball of an abstract Wiener space, providing in particular an explicit solution for the Neumann problem in this case. As an application, we derive an explicit formula for the Dirichlet-to-Neumann operator, which may be of independent interest.     

First Order Approach and Index Theorems for Discrete and Metric Graphs

The aim of the present paper is to introduce a unified notion of Laplacians on discrete and metric graphs. In order to cover all self-adjoint vertex conditions for the associated metric graph Laplacian, we develop systematically a new type of discrete graph operators acting on a decorated graph. The decoration at each vertex of degree d is given by a subspace of , generalising the fact that a function on the standard vertex space has only a scalar value. We illustrate the abstract concept by giving classical examples throughout the article. Our approach includes infinite graphs as well. We develop the notion of exterior derivative, differential forms, Dirac and Laplace operators in the discrete and metric case, using a supersymmetric framework. We calculate the (supersymmetric) index of the discrete Dirac operator generalising the standard index formula involving the Euler characteristic of a graph. Finally, we show that for finite graphs, the corresponding index for the metric Dirac operator agrees with the discrete one.     

CONSTRAINED QUADRILATERAL NONCONFORMING ROTATED Q1 ELEMENT

In this paper, we define a new nonconforming quadrilateral finite element based on the nonconforming rotated Q1 element by enforcing a constraint on each element, which has only three degrees of freedom. We investigate the consistency, approximation, superclose property, discrete Green＇s function and superconvergence of this element. Moreover, we propose a new postprocessing technique and apply it to this element. It is proved that the postprocessed discrete solution is superconvergent under a mild assumption on the mesh.     

Conformal cochains

In this paper we define holomorphic cochains and an associated period matrix for triangulated closed topological surfaces. We use the combinatorial Hodge star operator introduced in the author's paper of 2007, which depends on the choice of an inner product on the simplicial 1-cochains.

We prove that for a triangulated Riemannian 2-manifold (or a Riemann surface), and a particularly nice choice of inner product, the combinatorial period matrix converges to the (conformal) Riemann period matrix as the mesh of the triangulation tends to zero.

     

CONVERGENCE OF LAPLACIAN SPECTRA FROM RANDOM SAMPLES

Eigenvectors and eigenvalues of discrete Laplacians are often used for manifold learning and nonlinear dimensionality reduction. Graph Laplacian is one widely used discrete laplacian on point cloud. It was previously proved by Belkin and Niyogithat the eigenvectors and eigenvalues of the graph Laplacian converge to the eigenfunctions and eigenvalues of the Laplace-Beltrami operator of the manifold in the limit of infinitely many data points sampled independently from the uniform distribution over the manifold. Recently, we introduced Point Integral method (PIM) to solve elliptic equations and corresponding eigenvalue problem on point clouds. In this paper, we prove that the eigenvectors and eigenvalues obtained by PIM converge in the limit of infinitely many random samples. Moreover, estimation of the convergence rate is also given.     

Uniformly Distributed Circular Porous Pattern Generation on Surface for 3D Printing

We present an algorithm for uniformly distributed circular porous pattern generation on surface for three-dimensional (3D) printing using a phase-field model. The algorithm is based on the narrow band domain method for the nonlocal Cahn–Hilliard (CH) equation on surfaces. Surfaces are embedded in 3D grid and the narrow band domain is defined as the neighborhood of surface. It allows one can perform numerical computation using the standard discrete Laplacian in 3D instead of the discrete surface Laplacian. For complex surfaces, we reconstruct them from point cloud data and represent them as the zero-level set of their discrete signed distance functions. Using the proposed algorithm, we can generate uniformly distributed circular porous patterns on surfaces in 3D and print the resulting 3D models. Furthermore, we provide the test of accuracy and energy stability of the proposed method.     

Numerical analysis and applications of Fokker-Planck equations for stochastic dynamical systems with multiplicative α-stable noises

In this paper, we study the nonlocal Fokker-Planck equations (FPEs) associated with Lévy-driven scalar stochastic dynamical systems. We first derive the Fokker-Planck equation for the case of multiplicative symmetric α-stable noises, by the adjoint operator method. Then we construct a finite difference scheme to simulate the nonlocal FPE on either bounded or infinite domain. It is shown that the semi-discrete scheme satisfies the discrete maximum principle and converges. Some experiments are conducted to validate the numerical method. Finally, we extend the results to the asymmetric case and present an application to the nonlinear filtering problem.