Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
GPU Simulation of Rigid Fibers
KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
2015 (English)Independent thesis Advanced level (degree of Master (Two Years)), 20 credits / 30 HE creditsStudent thesisAlternative title
GPU simulering av stela fibrer (Swedish)
Abstract [en]

The major objective of this Master's thesis is to accelerate a serial implementation of a numerical algorithm for the simulation of slender fiber dynamics by using Graphical Processing Units (GPU). We focus on rigid fibers sedimenting due to gravity in a free-space Stokes flow. The ability to simulate a large number of fibers in a reasonable computational time on a high-performance parallel platform opens exciting new research opportunities.

The previous serial implementation is rewritten for parallel execution. The algorithm is implemented in single precision using the Compute Unified Device Architecture (CUDA) on NVIDIA GPUs. In addition, we develop an OpenMP version of the parallel implementation to run on multi-core CPUs. Using both implementations, we perform a number of benchmarks to determine the fastest variant of the algorithm. We observe a speedup of 20x to 40x on the NVIDIA GTX 970 compared to an Intel Core i7 4770. The GPU implementation can simulate up to 2000 fibers on a desktop computer and it takes only in the order of 8 seconds to advance one time step.

Furthermore, we have performed a number of simulations of known experiments for sedimenting fibers to validate the algorithm and to explore the numerical precision of the results. The results show an excellent agreement with results from prior experiments in the field.

Abstract [sv]

Huvudsyftet med detta arbete är att, med hjälp av grafikprocessorer (GPU:er), accelerera en seriell implementation av en numerisk algoritm för simulering av stela fibrer i en vätska. Fibrerna är modellerade på partikelnivå och ett stort antal fibrer krävs i beräkningarna för att kunna göra realistiska simuleringar. Trots förenklingar i modellen blir ekvationerna så komplexa att de måste lösas numeriskt och det krävs kraftfulla metoder för beräkningarna vilket innebär att snabba och effektiva algoritmer måste utvecklas.

Utgående från den seriella implementationen har en algoritm för parallell exekvering på GPU:er utvecklats. Algoritmen är implementerad i enkel precision på NVIDIA-GPU:er (NVIDIA GTX 970) med hjälp av Compute Unified Device Architecture, (CUDA). För att kunna göra en rättvis prestandajämförelse mellan GPU-implementationen och en CPU-implementering har dessutom en OpenMP-version av algoritmen implementerats för en CPU (Intel Core i7 4770) med flera kärnor. Båda versionerna har optimerats med avseende på respektive arkitektur. Ett antal prestandatester visar att den nya GPU-versionen av algoritmen ger en uppsnabbning av exekveringstiden på 20 till 40 gånger jämfört med CPU-OpenMP-versionen. Detta gör att vi nu kan simulera upp till 2000 fibrer på en bordsdator med en exekveringstid på ca 8 sekunder per tidssteg.

Vi har även utfört ett antal simuleringar av redan kända experiment med fibersuspensioner för att validera algoritmen samt för att utvärdera den numeriska precisionen i resultaten. Överensstämmelsen med tidigare rapporterade resultat är väldigt god.

Place, publisher, year, edition, pages
2015.
Series
TRITA-MAT-E, 2015:07
National Category
Computational Mathematics
Identifiers
URN: urn:nbn:se:kth:diva-162118OAI: oai:DiVA.org:kth-162118DiVA: diva2:797045
Subject / course
Scientific Computing
Educational program
Master of Science - Mathematics
Supervisors
Examiners
Available from: 2015-03-22 Created: 2015-03-22 Last updated: 2015-03-22Bibliographically approved

Open Access in DiVA

fulltext(3006 kB)117 downloads
File information
File name FULLTEXT01.pdfFile size 3006 kBChecksum SHA-512
1e141039513d4cbfcf35b4a8a88ff4b3c395c8d56ff65470b72d369e4c4ea09bed24ef76715d98d2149500823f045d970cf6dd0fe7c0e1a3ec25dfd3229e1845
Type fulltextMimetype application/pdf

By organisation
Numerical Analysis, NA
Computational Mathematics

Search outside of DiVA

GoogleGoogle Scholar
Total: 117 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

urn-nbn

Altmetric score

urn-nbn
Total: 106 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf