An Affine Semi-Lagrangian Advection Method

In computer graphics, the standard semi-Lagrangian advection as in the work of Stam (1999) is a widespread unconditionally stable transport scheme used in incompressible fluid solvers. Due to its stability, which disconnects the grid resolution from the time step required to prevent the numerical solution from blowing up, the method provides a good artistic control over the quality-performance trade-off. However, it is also notoriously known to include a great amount of artificial dissipation into the solution, hence destroying fine-scale details, and making simulated fluids appear overly viscous. Previous research efforts to counteract this unfortunate side effect have been spent notably on reinserting lost small-scale features, and on adapting different parts of the method to improve its accuracy. As part of the latter group, we present an affine semi-Lagrangian advection method, which we refer to as the ASLAM (pronounced "ay-slam"). This novel ASLAM adapts the locally affine descriptor of velocity from the affine particle-in-cell (APIC) method, a hybrid approach by Jiang et al. (2015), to the particle-free context of the Eulerian framework. We analyse the ASLAM's behaviour on a selection of testing scenarios, and evaluate it both qualitatively and quantitatively against a range of competing techniques, showing that it successfully reduces the artificial dissipation arising from standard semi-Lagrangian advection

Interactive freeform editing techniques for large-scale, multiresolution level set models

Level set methods provide a volumetric implicit surface representation with automatic smooth blending properties and no self-intersections. They can handle arbitrary topology changes easily, and the volumetric implicit representation does not require the surface to be re-adjusted after extreme deformations. Even though they have found some use in movie productions and some medical applications, level set models are not highly utilized in either special effects industry or medical science. Lack of interactive modeling tools makes working with level set models difficult for people in these application areas.This dissertation describes techniques and algorithms for interactive freeform editing of large-scale, multiresolution level set models. Algorithms are developed to map intuitive user interactions into level set speed functions producing specific, desired surface movements. Data structures for efficient representation of very high resolution volume datasets and associated algorithms for rapid access and processing of the information within the data structures are explained. A hierarchical, multiresolution representation of level set models that allows for rapid decomposition and reconstruction of the complete full-resolution model is created for an editing framework that allows level-of-detail editing. We have developed a framework that identifies surface details prior to editing and introduces them back afterwards. Combining these two features provides a detail-preserving level set editing capability that may be used for multi-resolution modeling and texture transfer. Given the complex data structures that are required to represent large-scale, multiresolution level set models and the compute-intensive numerical methods to evaluate them, optimization techniques and algorithms have been developed to evaluate and display the dynamic isosurface embedded in the volumetric data.Ph.D., Computer Science -- Drexel University, 201

Problems in cosmology and numerical relativity

Includes bibliographical references.A generic feature of most inflationary scenarios is the generation of primordial perturbations. Ordinarily, such perturbations can interact with a weak magnetic field in a plasma, resulting in a wide range of phenomena, such as the parametric excitation of plasma waves by gravitational waves. This mechanism has been studied in different contexts in the literature, such as the possibility of indirect detection of gravitational waves through electromagnetic signatures of the interaction. In this work, we consider this concept in the particular case of magnetic field amplification. Specifically, we use non-linear gauge-in variant perturbation theory to study the interaction of a primordial seed magnetic field with density and gravitational wave perturbations in an almost Friedmann-Lemaıtre-Robertson- Walker (FLRW) spacetime with zero spatial curvature. We compare the effects of this coupling under the assumptions of poor conductivity, perfect conductivity and the case where the electric field is sourced via the coupling of velocity perturbations to the seed field in the ideal magnetohydrodynamic (MHD) regime, thus generalizing, improving on and correcting previous results. We solve our equations for long wavelength limits and numerically integrate the resulting equations to generate power spectra for the electromagnetic field variables, showing where the modes cross the horizon. We find that the interaction can seed Electric fields with non-zero curl and that the curl of the electric field dominates the power spectrum on small scales, in agreement with previous arguments. The second focus area of the thesis is the development a stable high order mesh refinement scheme for the solution of hyperbolic partial differential equations. It has now become customary in the field of numerical relativity to couple high order finite difference schemes to mesh refinement algorithms. This approach combines the efficiency of local mesh refinement with the robustness and accuracy of higher order methods. To this end, different modifications of the standard Berger-Oliger adaptive mesh refinement a logarithm have been proposed. In this work we present a new fourth order convergent mesh refinement scheme with sub- cycling in time for numerical relativity applications. One of the distinctive features of our algorithm is that we do not use buffer zones to deal with refinement boundaries, as is currently done in the literature, but explicitly specify boundary data for refined grids instead. We argue that the incompatibility of the standard mesh refinement algorithm with higher order Runge Kutta methods is a manifestation of order reduction phenomena which is caused by inconsistent application of boundary data in the refined grids. Indeed, a peculiar feature of high order explicit Runge Kutta schemes is that they behave like low order schemes when applied to hyperbolic problems with time dependent Dirichlet boundary conditions. We present a new algorithm to deal with this phenomenon and through a series of examples demonstrate fourth order convergence. Our scheme also addresses the problem of spurious reflections that are generated when propagating waves cross mesh refinement boundaries. We introduce a transition zone on refined levels within which the phase velocity of propagating modes is allowed to decelerate in order to smoothly match the phase velocity of coarser grids. We apply the method to test problems involving propagating waves and show a significant reduction in spurious reflections

Particle based modeling and simulation of natural phenomena

Ankara : The Department of Computer Engineering and the Institute of Engineering and Science of Bilkent University, 2010.Thesis (Ph. D.) -- Bilkent University, 2010.Includes bibliographical references leaves 92-108.This thesis is about modeling and simulation of fluids and cloth-like deformable objects by the physically-based simulation paradigm. Simulated objects are modeled with particles and their interaction with each other and the environment is defined by particle-to-particle forces. We propose several improvements over the existing particle simulation techniques. Neighbor search algorithms are crucial for the performance efficiency and robustness of a particle system. We present a sorting-based neighbor search method which operates on a uniform grid, and can be parallelizable. We improve upon the existing fluid surface generation methods so that our method captures surface details better since we consider the relative position of fluid particles to the fluid surface. We investigate several alternatives of particle interaction schema (i.e. Smoothed Particle Hydrodynamics, the Discrete Element Method, and Lennard-Jones potential) for the purpose of defining fluid-fluid, fluid-cloth, fluid-boundary interaction forces. We also propose a practical way to simulate knitwear and its interaction with fluids. We employ capillary pressure–based forces to simulate the absorption of fluid particles by knitwear. We also propose a method to simulate the flow of miscible fluids. Our particle simulation system is implement to exploit parallel computing capabilities of the commodity computers. Specifically, we implemented the proposed methods on multicore CPUs and programmable graphics boards. The experiments show that our method is computationally efficient and produces realistic results.Bayraktar, SerkanPh.D

Task Allocation in Foraging Robot Swarms:The Role of Information Sharing

Autonomous task allocation is a desirable feature of robot swarms that collect and deliver items in scenarios where congestion, caused by accumulated items or robots, can temporarily interfere with swarm behaviour. In such settings, self-regulation of workforce can prevent unnecessary energy consumption. We explore two types of self-regulation: non-social, where robots become idle upon experiencing congestion, and social, where robots broadcast information about congestion to their team mates in order to socially inhibit foraging. We show that while both types of self-regulation can lead to improved energy efficiency and increase the amount of resource collected, the speed with which information about congestion flows through a swarm affects the scalability of these algorithms