Real-Time Particle-Based Simulation on GPUs (sap 0151) Takahiro Harada∗

Masayuki Tanaka†

Seiichi Koshizuka‡

Yoichiro Kawaguchi§

The University of Tokyo

Figure 1: Real-time simulation of glasses and liquid. Glass tower is filled with liquid and a glass is thrown into the scene. This simulation runs 17.1 frames per second on GeForce 8800GTX.

1

Introduction

The current trend in processor technology is to improve the efficiency of processors and not increase their frequency. Processors nowadays are equipped with parallel architecture. Cell Broadband Engine Architecture is a multi-core processor for general-purpose computation and Graphics Processing Units (GPUs) are specialized parallel processors for graphics tasks. Additionally, CPUs are also shifting to multi-core design. What we need to do now is adaptation to these platforms. Therefore, we need to develop data-parallel algorithms that exploit their computational powers.

particles (spheres) as the title of this skech implies. We call them rigid particles. The size of rigid particles is all the same and also the same to the size of fluid particles. An advantage of this shape representation is the computation speed is controllable by changing the accuracy, i.e., the resolution of particles. With this shape representation, not only collision between rigid bodies but also interaction between a rigid body and a fluid can be converted to the simple problem of computation of particle interactions. Thus, the computation is simple and it can be computed in parallel. However, the shape representation using particles increases the number of simulation entites because a rigid body is consists of a few rigid particles. A uniform grid is introduced to make the neighboring particle search efficient. The interaction between fluid particles and rigid particles is calculated by assuming rigid bodies as a fluid. The density is also computed on rigid particles and then the pressure and viscosity forces are calculated between fluid particles. The force on the rigid particle, which is the sum of the force from fluid and the force from collisions between rigid particles, is used to update the linear and anguler momenta of a rigid body.

In this sketch, we show that a particle-based simulation can be parallelized and implemented entirely on Graphics Processing Units (GPUs) as a parallel computation platform. As a result, we can obtain unprecedented performance with scalar processors. We also presents a particle-based method to interact fluids and rigid bodies. In this method, rigid bodies are represented by a set of particles. The benefits of this method are low computational cost and parallelism of its algorithm.

As described above, the force computation between rigid bodies and a fluid can be executed in parallel. When GPUs are used, a frame buffer is rendered with a fragment shader by assigning a pixel to a particle. The fragment shader compute the force with physical values of neighboring particles which are read from other textures. A data which stores an information about neighboring particles is generated in advance by a vertex shader. In this way, forces on rigid bodies and a fluid are computed in parallel.

As physical laws govern the motion of objects around us, a physically-based simulation plays an important role in computer graphics. For instance, the motion of a fluid, which is difficult to generate by hand, can be produced by solving the governing equations. Acceleration of a simulation is one of the most important research themes because the speed and stability of a simulation are essential for real-time applications.

2

Methods

Smoothed Particle Hydrodynamics (SPH) is employed to solve the governing equation of a fluid[M¨uller et al. 2003]. A characteristics of particle method including SPH is that there is no numerical dissipation caused by advection calculation and so mass loss does not occure even if the resolution of a simulation is low. Therefore, the particle methods are suited for a real-time application. As for the rigid body simulation, a rigid body is represented by a set of ∗ e-mail:

[email protected]

† e-mail:[email protected] ‡ e-mail:[email protected] § e-mail:[email protected]

3

Results

In Figure 1, 10 glasses are stacked and a fluid is poured from above them. Then, a glass is thrown onto them and the stacked glasses collapse. This simulation uses 49,153 particles and runs 17.1 frames per second on GeForce 8800GTX using a rendering in which point sprites are used to render particles. The simulator outputs simulation data and the surface of the fluid is constructed by Marching Cubes by assigning densities to fluid particles. The polygons are rendered after the simulation with an offline renderer. The accompanying video includes several examples which runs in real-time. A simulation which uses the largest particle number runs 3.85 frames per second with 245,760 particles. These examples show the capability of the present technique.

References ¨ M ULLER , M., C HARYPAR , D., AND G ROSS , M. 2003. Particle-based fluid simulation for interactive applications. In Proc. of SIGGRAPH Symposium on Computer Animation, 154–159.

Real-Time Particle-Based Simulation on GPUs - Semantic Scholar

tion to these platforms. Therefore, we need to ... (GPUs) as a parallel computation platform. As a result, we can ... ∗e-mail: [email protected].

359KB Sizes 0 Downloads 263 Views

Recommend Documents

Real-Time Particle-Based Simulation on GPUs - Semantic Scholar
†e-mail:[email protected]. ‡e-mail:[email protected]. §e-mail:[email protected] particles (spheres) as the title of this skech implies ...

On numerical simulation of high-speed CCD ... - Semantic Scholar
have a smaller photosensitive area that cause higher photon shot noise, higher ... that there are i interactions per pixel and PI is the number of interacting photons. .... to the random trapping and emission of mobile charge carriers resulting in ..

On Knowledge - Semantic Scholar
Rhizomatic Education: Community as Curriculum by Dave Cormier. The truths .... Couros's graduate-level course in educational technology offered at the University of Regina provides an .... Techknowledge: Literate practice and digital worlds.

On Knowledge - Semantic Scholar
Rhizomatic Education: Community as Curriculum .... articles (Nichol 2007). ... Couros's graduate-level course in educational technology offered at the University ...

Adiabatic Quantum Simulation of Quantum ... - Semantic Scholar
Oct 13, 2014 - quantum adiabatic algorithm to combinatorial optimization problems. ... applied to structured and unstructured search20,21, search engine ...... License. The images or other third party material in this article are included in the.

Parallel generation of samples for simulation ... - Semantic Scholar
Analytical modeling of complex systems is crucial to de- tect error conditions or ... The current SAN solver, PEPS software tool [4], works with less than 65 million ...

Parallel generation of samples for simulation ... - Semantic Scholar
This advantage justifies its usage in several contexts where .... The main advantage of. SAN is due to ..... ular Analytical Performance Models for Ad Hoc Wireless.

Scalable Quantum Simulation of Molecular Energies - Semantic Scholar
Jul 18, 2016 - Errors in our simulation as a function of R are shown in ..... P. J. J. O. compile quantum software and analyze data. R. Babbush, P. J. J. O., and ...

There are several different simulation paradigms ... - Semantic Scholar
A Comparison of System Dynamics and Agent-Based Simulation. Applied to the Study of ... based; 2) traditional stochastic, discrete time, discrete entity Monte Carlo ..... unbound molecules over time may be saved to a file for subsequent ...

Exponentially more precise quantum simulation of ... - Semantic Scholar
Mar 24, 2016 - Keywords: quantum algorithms, quantum simulation, electronic structure ..... Using atomic units in which the electron mass, electron charge, ...

Agent Based Modelling and Simulation of the ... - Semantic Scholar
guage. Independently of the programming language, ImmSim has a very rigid ... C-ImmSim and the correspondent parallel variant, ParImm, are versions of Imm-.

High-accuracy simulation of density driven flow in ... - Semantic Scholar
software tools and computing resources. In this paper a recently .... analytical form; in these cases highly accurate reference solutions have to be employed for ...

Simulation of 3D neuro-musculo-skeletal systems ... - Semantic Scholar
between them, is taken into account to determine the action of the forces generated in ... A graphics-based software system to develop and analyze models.

Exponentially more precise quantum simulation of ... - Semantic Scholar
Mar 24, 2016 - significantly more practical Trotter decompositions, the best known gate complexity ... The ancilla register is then put in a superposition state with .... integral in equation (4) usingμ grid points where the domain of the integral, 

High-accuracy simulation of density driven flow in ... - Semantic Scholar
software tools and computing resources. In this paper a recently .... analytical form; in these cases highly accurate reference solutions have to be employed for ...

building blocks of a simulation environment of the ... - Semantic Scholar
INTRODUCTION. The global Internet, wireless communication sys- .... the delay time of the packet (time of delivery at destination minus .... form a distributed routing table update of the whole network. ... We are not content with one sequence of ...

Agent Based Modelling and Simulation of the ... - Semantic Scholar
Generally, analytic treatment does not yield a complex system's ... based solutions is performed by comparing its perfomance with other approaches to the same ...

Online Video Recommendation Based on ... - Semantic Scholar
Department of Computer Science and Technology, Tsinghua University, Beijing 100084, P. R. ... precedented level, video recommendation has become a very.

On Deconstructing Ensemble Models - Semantic Scholar
Oct 1, 2015 - metrics. Our response measures a shift in user behavior observable only after a longer ..... Center G : G0 ≡ G − 1 ¯gT , where ¯g(m) = ave G(-, m). .... importance, we call the calculations based on ∆R2(j) an analysis of trees.

Article On Being Systematically Connectionist - Semantic Scholar
Address for correspondence: Lars F. Niklasson, Department of Computer Science, University of. Skövde, S-54128, SWEDEN. Email: [email protected]. Tim van Gelder, Philosophy Program,. Research School of Social Sciences, Australian National University, C

On Decomposability of Multilinear Sets - Semantic Scholar
May 28, 2016 - Examples include quadratic programs, polynomial programs, and multiplicative ... †Graduate Program in Operations Research and Industrial ...

Language Recognition Based on Score ... - Semantic Scholar
1School of Electrical and Computer Engineering. Georgia Institute of ... over all competing classes, and have been demonstrated to be effective in isolated word ...

Language Recognition Based on Score ... - Semantic Scholar
1School of Electrical and Computer Engineering. Georgia Institute ... NIST (National Institute of Standards and Technology) has ..... the best procedure to follow.

International Evidence on Sticky Consumption ... - Semantic Scholar
Oct 21, 2009 - We compare the random walk model of. 1Facts that have been ..... one-time events like the hosting of the Olympics). Standard theoretical ...