Cyclic animation using Partial differential Equations

YesThis work presents an efficient and fast method for achieving cyclic animation using Partial Differential Equations (PDEs). The boundary-value nature associ- ated with elliptic PDEs offers a fast analytic solution technique for setting up a framework for this type of animation. The surface of a given character is thus cre- ated from a set of pre-determined curves, which are used as boundary conditions so that a number of PDEs can be solved. Two different approaches to cyclic ani- mation are presented here. The first consists of using attaching the set of curves to a skeletal system hold- ing the animation for cyclic motions linked to a set mathematical expressions, the second one exploits the spine associated with the analytic solution of the PDE as a driving mechanism to achieve cyclic animation, which is also manipulated mathematically. The first of these approaches is implemented within a framework related to cyclic motions inherent to human-like char- acters, whereas the spine-based approach is focused on modelling the undulatory movement observed in fish when swimming. The proposed method is fast and ac- curate. Additionally, the animation can be either used in the PDE-based surface representation of the model or transferred to the original mesh model by means of a point to point map. Thus, the user is offered with the choice of using either of these two animation repre- sentations of the same object, the selection depends on the computing resources such as storage and memory capacity associated with each particular application

An Annular Plate Model in Arbitrary-Lagrangian-Eulerian Description for the DLR FlexibleBodies Library

The bending deformation of rotating annular plates and the associated vibration behaviour is important in engineering applications which range from automotive or railway brake systems to discs that form essential components in turbomachinery. In order to extend the capabilities of the DLR FlexibleBodies library for such use cases, a new Modelica class has been implemented which is based on the analytical description of an annular Kirchhoff plate. In addition the so-called Arbitray Langrangian-Eulerian (ALE) representation has been adopted so that rotating and non-rotating external loads may be applied conventiently to rotating plates. Besides these particularities the new class AnnularPlate completely corresponds to the concept of FlexibleBodies library with the two already available model classes Beam and ModalBody. This paper gives an overview on the theoretical background of the new class AnnularPlate, explains the usage and presents application examples

Waltzing peakons and compacton pairs in a cross-coupled Camassa-Holm equation

We consider singular solutions of a system of two cross-coupled Camassa-Holm (CCCH) equations. This CCCH system admits peakon solutions, but it is not in the two-component CH integrable hierarchy. The system is a pair of coupled Hamiltonian partial differential equations for two types of solutions on the real line, each of which separately possesses exp(-|x|) peakon solutions with a discontinuity in the first derivative at the peak. However, there are no self-interactions, so each of the two types of peakon solutions moves only under the induced velocity of the other type. We analyse the `waltzing' solution behaviour of the cases with a single bound peakon pair (a peakon couple), as well as the over-taking collisions of peakon couples and the antisymmetric case of the head-on collision of a peakon couple and a peakon anti-couple. We then present numerical solutions of these collisions, which are inelastic because the waltzing peakon couples each possess an internal degree of freedom corresponding to their `tempo' -- that is, the period at which the two peakons of opposite type in the couple cycle around each other in phase space. Finally, we discuss compacton couple solutions of the cross-coupled Euler-Poincar\'e (CCEP) equations and illustrate the same types of collisions as for peakon couples, with triangular and parabolic compacton couples. We finish with a number of outstanding questions and challenges remaining for understanding couple dynamics of the CCCH and CCEP equations

Optimising Spatial and Tonal Data for PDE-based Inpainting

Some recent methods for lossy signal and image compression store only a few selected pixels and fill in the missing structures by inpainting with a partial differential equation (PDE). Suitable operators include the Laplacian, the biharmonic operator, and edge-enhancing anisotropic diffusion (EED). The quality of such approaches depends substantially on the selection of the data that is kept. Optimising this data in the domain and codomain gives rise to challenging mathematical problems that shall be addressed in our work. In the 1D case, we prove results that provide insights into the difficulty of this problem, and we give evidence that a splitting into spatial and tonal (i.e. function value) optimisation does hardly deteriorate the results. In the 2D setting, we present generic algorithms that achieve a high reconstruction quality even if the specified data is very sparse. To optimise the spatial data, we use a probabilistic sparsification, followed by a nonlocal pixel exchange that avoids getting trapped in bad local optima. After this spatial optimisation we perform a tonal optimisation that modifies the function values in order to reduce the global reconstruction error. For homogeneous diffusion inpainting, this comes down to a least squares problem for which we prove that it has a unique solution. We demonstrate that it can be found efficiently with a gradient descent approach that is accelerated with fast explicit diffusion (FED) cycles. Our framework allows to specify the desired density of the inpainting mask a priori. Moreover, is more generic than other data optimisation approaches for the sparse inpainting problem, since it can also be extended to nonlinear inpainting operators such as EED. This is exploited to achieve reconstructions with state-of-the-art quality. We also give an extensive literature survey on PDE-based image compression methods

Implementation of a partitioned algorithm for simulation of large CSI problems

The implementation of a partitioned numerical algorithm for determining the dynamic response of coupled structure/controller/estimator finite-dimensional systems is reviewed. The partitioned approach leads to a set of coupled first and second-order linear differential equations which are numerically integrated with extrapolation and implicit step methods. The present software implementation, ACSIS, utilizes parallel processing techniques at various levels to optimize performance on a shared-memory concurrent/vector processing system. A general procedure for the design of controller and filter gains is also implemented, which utilizes the vibration characteristics of the structure to be solved. Also presented are: example problems; a user's guide to the software; the procedures and algorithm scripts; a stability analysis for the algorithm; and the source code for the parallel implementation

Parametric Forcing of Confined and Stratified Flows

abstract: A continuously and stably stratified fluid contained in a square cavity subjected to harmonic body forcing is studied numerically by solving the Navier-Stokes equations under the Boussinesq approximation. Complex dynamics are observed near the onset of instability of the basic state, which is a flow configuration that is always an exact analytical solution of the governing equations. The instability of the basic state to perturbations is first studied with linear stability analysis (Floquet analysis), revealing a multitude of intersecting synchronous and subharmonic resonance tongues in parameter space. A modal reduction method for determining the locus of basic state instability is also shown, greatly simplifying the computational overhead normally required by a Floquet study. Then, a study of the nonlinear governing equations determines the criticality of the basic state's instability, and ultimately characterizes the dynamics of the lowest order spatial mode by the three discovered codimension-two bifurcation points within the resonance tongue. The rich dynamics include a homoclinic doubling cascade that resembles the logistic map and a multitude of gluing bifurcations. The numerical techniques and methodologies are first demonstrated on a homogeneous fluid contained within a three-dimensional lid-driven cavity. The edge state technique and linear stability analysis through Arnoldi iteration are used to resolve the complex dynamics of the canonical shear-driven benchmark problem. The techniques here lead to a dynamical description of an instability mechanism, and the work serves as a basis for the remainder of the dissertation.Dissertation/ThesisSupplemental Materials Description Filezip file containing 10 mp4 formatted video animations, as well as a text readme and the previously submitted Supplemental Materials Description FileDoctoral Dissertation Mathematics 201