Ändra sökning
RefereraExporteraLänk till posten
Permanent länk

Direktlänk
Referera
Referensformat
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Annat format
Fler format
Språk
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Annat språk
Fler språk
Utmatningsformat
  • html
  • text
  • asciidoc
  • rtf
Spectral accuracy in fast Ewald-based methods for particle simulations
KTH, Skolan för datavetenskap och kommunikation (CSC), Numerisk analys, NA.
KTH, Skolan för datavetenskap och kommunikation (CSC), Numerisk analys, NA.
2011 (Engelska)Ingår i: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 230, nr 24, s. 8744-8761Artikel i tidskrift (Refereegranskat) Published
Abstract [en]

A spectrally accurate fast method for electrostatic calculations under periodic boundary conditions is presented. We follow the established framework of FFT-based Ewald summation, but obtain a method with an important decoupling of errors: it is shown, for the proposed method, that the error due to frequency domain truncation can be separated from the approximation error added by the fast method. This has the significance that the truncation of the underlying Ewald sum prescribes the size of the grid used in the FFT-based fast method, which clearly is the minimal grid. Both errors are of exponential-squared order, and the latter can be controlled independently of the grid size. We compare numerically to the established SPME method by Essmann et al. and see that the memory required can be reduced by orders of magnitude. We also benchmark efficiency (i.e. error as a function of computing time) against the SPME method, which indicates that our method is competitive. Analytical error estimates are proven and used to select parameters with a great degree of reliability and ease.

Ort, förlag, år, upplaga, sidor
2011. Vol. 230, nr 24, s. 8744-8761
Nyckelord [en]
Ewald summation, FFT, Molecular dynamics, PME, Spectral accuracy, SPME
Nationell ämneskategori
Beräkningsmatematik Annan fysik
Identifikatorer
URN: urn:nbn:se:kth:diva-48765DOI: 10.1016/j.jcp.2011.08.022ISI: 000297081700007Scopus ID: 2-s2.0-80053628839OAI: oai:DiVA.org:kth-48765DiVA, id: diva2:458547
Forskningsfinansiär
Knut och Alice Wallenbergs StiftelseSwedish e‐Science Research Center
Anmärkning
QC 20111124Tillgänglig från: 2011-11-23 Skapad: 2011-11-23 Senast uppdaterad: 2017-12-08Bibliografiskt granskad
Ingår i avhandling
1. Spectral Accuracy in Fast Ewald Methods and Topics in Fluid Interface Simulation
Öppna denna publikation i ny flik eller fönster >>Spectral Accuracy in Fast Ewald Methods and Topics in Fluid Interface Simulation
2011 (Engelska)Doktorsavhandling, sammanläggning (Övrigt vetenskapligt)
Abstract [en]

  This work contains two separate but related parts: one on spectrally  accurate and fast Ewald methods for electrostatics and viscous flow,  and one on micro- and complex fluid interface problems.  In Part I we are concerned with fast and spectrally accurate methods  to compute sums of slowly decaying potentials over periodic  lattices. We consider two PDEs: Laplace (electrostatics, the Coulomb  potential) and Stokes (viscous flow, the ``Stokeslet''  potential). Moreover, we consider both full and planar periodicity,  the latter meaning that periodicity applies in two dimensions and  the third is ``free''. These are major simulation tasks in current  molecular dynamics simulations and in many areas of computational  fluid mechanics involving e.g. particle suspensions.   For each of the four combinations of PDE and periodic structure, we  give spectrally accurate and O(N log N) fast methods based on  Ewald's or Ewald-like decompositions of the underlying potential  sums. In the plane-periodic cases we derive the decompositions in a  manner that lets us develop fast methods. Associated error estimates  are developed as needed throughout. All four methods can be placed  in the P3M/PME (Particle Mesh Ewald) family. We argue that they  have certain novel and attractive features: first, they are spectral  accurate; secondly, they use the minimal amount of memory possible  within the PME family; third, each has a clear and reliable view of  numerical errors, such that parameters can be chosen  wisely. Analytical and numerical results are given to support these  propositions. We benchmark accuracy and performance versus an  established (S)PME method.  Part II deals with free boundary problems, specifically numerical  methods for multiphase flow. We give an interface tracking method  based on a domain-decomposition idea that lets us split the  interface into overlapping patches. Each patch is discretized on a  uniform grid, and accurate and efficient numerical methods are given  for the equations that govern interface transport. We demonstrate  that the method is accurate and how it's used in immersed boundary,  and interface, Navier-Stokes methods, as well as in a boundary  integral Stokes setting.  Finally, we consider a problem in complex fluidics where there is a  concentration of surfactants \emph{on} the interface and the  interface itself is in contact with a solid boundary (the contact  line problem). We argue that the domain-decomposition framework is  attractive for formulating and treating complex models  (e.g. involving PDEs on a dynamic interface) and proceed with  developing various aspects of such a method.

Ort, förlag, år, upplaga, sidor
Stockholm: KTH Royal Institute of Technology, 2011. s. xv, 104
Serie
Trita-CSC-A, ISSN 1653-5723 ; 2011:19
Nationell ämneskategori
Beräkningsmatematik
Identifikatorer
urn:nbn:se:kth:diva-48805 (URN)978-91-7501-195-0 (ISBN)
Disputation
2011-12-16, Salongen, KTHB, Osquars backe 25, Stockholm, 10:00 (Engelska)
Opponent
Handledare
Forskningsfinansiär
Swedish e‐Science Research Center
Anmärkning
QC 20111125Tillgänglig från: 2011-11-25 Skapad: 2011-11-23 Senast uppdaterad: 2012-05-24Bibliografiskt granskad

Open Access i DiVA

Fulltext saknas i DiVA

Övriga länkar

Förlagets fulltextScopus

Sök vidare i DiVA

Av författaren/redaktören
Lindbo, DagTornberg, Anna-Karin
Av organisationen
Numerisk analys, NA
I samma tidskrift
Journal of Computational Physics
BeräkningsmatematikAnnan fysik

Sök vidare utanför DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetricpoäng

doi
urn-nbn
Totalt: 128 träffar
RefereraExporteraLänk till posten
Permanent länk

Direktlänk
Referera
Referensformat
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Annat format
Fler format
Språk
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Annat språk
Fler språk
Utmatningsformat
  • html
  • text
  • asciidoc
  • rtf