An explicit structure‐preserving numerical scheme for EPDiff

We present a new structure‐preserving numerical scheme for solving the Euler‐Poincaré Differential (EPDiff) equation on arbitrary triangle meshes. Unlike existing techniques, our method solves the difficult non‐linear EPDiff equation by constructing energy preserving, yet fully explicit, update rule...

Full description

Saved in:
Bibliographic Details
Published in:Computer graphics forum Vol. 37; no. 5; pp. 107 - 119
Main Authors: Azencot, Omri, Vantzos, Orestis, Ben‐Chen, Mirela
Format: Journal Article
Language:English
Published: Oxford Blackwell Publishing Ltd 01.08.2018
Subjects:
Categories and Subject Descriptors (according to ACM CCS)
Computer simulation
Curvature
Differential equations
I.3.3 [Computer Graphics]: Three‐Dimensional Graphics and Realism—Animation
I.3.5 [Computer Graphics]: Computational Geometry and Object Modeling—Physically based modeling
Mathematical analysis
Matrix methods
Splitting
Test procedures
Wave fronts
ISSN:0167-7055, 1467-8659
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We present a new structure‐preserving numerical scheme for solving the Euler‐Poincaré Differential (EPDiff) equation on arbitrary triangle meshes. Unlike existing techniques, our method solves the difficult non‐linear EPDiff equation by constructing energy preserving, yet fully explicit, update rules. Our approach uses standard differential operators on triangle meshes, allowing for a simple and efficient implementation. Key to the structure‐preserving features that our method exhibits is a novel numerical splitting scheme. Namely, we break the integration into three steps which rely on linear solves with a fixed sparse matrix that is independent of the simulation and thus can be pre‐factored. We test our method in the context of simulating concentrated reconnecting wavefronts on flat and curved domains. In particular, EPDiff is known to generate geometrical fronts which exhibit wave‐like behavior when they interact with each other. In addition, we also show that at a small additional cost, we can produce globally‐supported periodic waves by using our simulated fronts with wavefronts tracking techniques. We provide quantitative graphs showing that our method exactly preserves the energy in practice. In addition, we demonstrate various interesting results including annihilation and recreation of a circular front, a wave splitting and merging when hitting an obstacle and two separate fronts propagating and bending due to the curvature of the domain.
Bibliography:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:0167-7055
1467-8659
DOI:10.1111/cgf.13495