Parallelized Rigid Body Dynamics

Physics engines are collections of API-like software designed for video games, movies and scientific simulations. While physics engines often come in many shapes and designs, all engines can benefit from an increase in speed via parallelization. However, despite this need for increased speed, it is uncommon to encounter a parallelized physics engine today. Many engines are long-standing projects and changing them to support parallelization is too costly to consider as a practical matter. Parallelization needs to be considered from the design stages through completion to ensure adequate implementation. In this project we develop a realistic approach to simulate physics in a parallel environment. Utilizing many techniques we establish a practical approach to significantly reduce the run-time on a standard physics engine

Bounding stationary averages of polynomial diffusions via semidefinite programming

We introduce an algorithm based on semidefinite programming that yields increasing (resp. decreasing) sequences of lower (resp. upper) bounds on polynomial stationary averages of diffusions with polynomial drift vector and diffusion coefficients. The bounds are obtained by optimising an objective, determined by the stationary average of interest, over the set of real vectors defined by certain linear equalities and semidefinite inequalities which are satisfied by the moments of any stationary measure of the diffusion. We exemplify the use of the approach through several applications: a Bayesian inference problem; the computation of Lyapunov exponents of linear ordinary differential equations perturbed by multiplicative white noise; and a reliability problem from structural mechanics. Additionally, we prove that the bounds converge to the infimum and supremum of the set of stationary averages for certain SDEs associated with the computation of the Lyapunov exponents, and we provide numerical evidence of convergence in more general settings

Convergence of Entropic Schemes for Optimal Transport and Gradient Flows

Replacing positivity constraints by an entropy barrier is popular to approximate solutions of linear programs. In the special case of the optimal transport problem, this technique dates back to the early work of Schr\"odinger. This approach has recently been used successfully to solve optimal transport related problems in several applied fields such as imaging sciences, machine learning and social sciences. The main reason for this success is that, in contrast to linear programming solvers, the resulting algorithms are highly parallelizable and take advantage of the geometry of the computational grid (e.g. an image or a triangulated mesh). The first contribution of this article is the proof of the $\Gamma$-convergence of the entropic regularized optimal transport problem towards the Monge-Kantorovich problem for the squared Euclidean norm cost function. This implies in particular the convergence of the optimal entropic regularized transport plan towards an optimal transport plan as the entropy vanishes. Optimal transport distances are also useful to define gradient flows as a limit of implicit Euler steps according to the transportation distance. Our second contribution is a proof that implicit steps according to the entropic regularized distance converge towards the original gradient flow when both the step size and the entropic penalty vanish (in some controlled way)

SPH with the multiple boundary tangent method

In this article, we present an improved solid boundary treatment formulation for the smoothed particle hydrodynamics (SPH) method. Benchmark simulations using previously reported boundary treatments can suffer from particle penetration and may produce results that numerically blow up near solid boundaries. As well, current SPH boundary approaches do not properly treat curved boundaries in complicated flow domains. These drawbacks have been remedied in a new boundary treatment method presented in this article, called the multiple boundary tangent (MBT) approach. In this article we present two important benchmark problems to validate the developed algorithm and show that the multiple boundary tangent treatment produces results that agree with known numerical and experimental solutions. The two benchmark problems chosen are the lid-driven cavity problem, and flow over a cylinder. The SPH solutions using the MBT approach and the results from literature are in very good agreement. These solutions involved solid boundaries, but the approach presented herein should be extendable to time-evolving, free-surface boundaries

Verification of system properties of polynomial systems using discrete-time approximations and set-based analysis

Magdeburg, Univ., Fak. für Elektrotechnik und Informationstechnik, Diss., 2015von Philipp Rumschinsk

PT-symmetric operators and metastable states of the 1D relativistic oscillators

We consider the one-dimensional Dirac equation for the harmonic oscillator and the associated second order separated operators giving the resonances of the problem by complex dilation. The same operators have unique extensions as closed PT-symmetric operators defining infinite positive energy levels converging to the Schroedinger ones as c tends to infinity. Such energy levels and their eigenfunctions give directly a definite choice of metastable states of the problem. Precise numerical computations shows that these levels coincide with the positions of the resonances up to the order of the width. Similar results are found for the Klein-Gordon oscillators, and in this case there is an infinite number of dynamics and the eigenvalues and eigenvectors of the PT-symmetric operators give metastable states for each dynamics.Comment: 13 pages, 2 figure

Archipelagian Cosmology: Dynamics and Observables in a Universe with Discretized Matter Content

We consider a model of the Universe in which the matter content is in the form of discrete islands, rather than a continuous fluid. In the appropriate limits the resulting large-scale dynamics approach those of a Friedmann-Robertson-Walker (FRW) universe. The optical properties of such a space-time, however, do not. This illustrates the fact that the optical and `average' dynamical properties of a relativistic universe are not equivalent, and do not specify each other uniquely. We find the angular diameter distance, luminosity distance and redshifts that would be measured by observers in these space-times, using both analytic approximations and numerical simulations. While different from their counterparts in FRW, the effects found do not look like promising candidates to explain the observations usually attributed to the existence of Dark Energy. This incongruity with standard FRW cosmology is not due to the existence of any unexpectedly large structures or voids in the Universe, but only to the fact that the matter content of the Universe is not a continuous fluid.Comment: 49 pages, 15 figures. Corrections made to description of lattice constructio