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
Fast and spectrally accurate Ewald summation for 2-periodic electrostatic systems
KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).
KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).
2012 (English)In: Journal of Chemical Physics, ISSN 0021-9606, E-ISSN 1089-7690, Vol. 136, no 16, 164111-1-164111-16 p.Article in journal (Refereed) Published
Abstract [en]

A new method for Ewald summation in planar/slablike geometry, i.e., systems where periodicity applies in two dimensions and the last dimension is "free" (2P), is presented. We employ a spectral representation in terms of both Fourier series and integrals. This allows us to concisely derive both the 2P Ewald sum and a fast particle mesh Ewald (PME)-type method suitable for large-scale computations. The primary results are: (i) close and illuminating connections between the 2P problem and the standard Ewald sum and associated fast methods for full periodicity; (ii) a fast, O(N log N), and spectrally accurate PME-type method for the 2P k-space Ewald sum that uses vastly less memory than traditional PME methods; (iii) errors that decouple, such that parameter selection is simplified. We give analytical and numerical results to support this.

Place, publisher, year, edition, pages
2012. Vol. 136, no 16, 164111-1-164111-16 p.
Keyword [en]
Ewald summations, Fast methods, Fast particle, K-space, Numerical results, Parameter selection, Spectral representations, Two-dimension, Type methods
National Category
Physical Sciences Computational Mathematics
Identifiers
URN: urn:nbn:se:kth:diva-48766DOI: 10.1063/1.4704177ISI: 000303602200013Scopus ID: 2-s2.0-84860487992OAI: oai:DiVA.org:kth-48766DiVA: diva2:458562
Funder
Knut and Alice Wallenberg FoundationSwedish e‐Science Research Center
Note

Updated from manuscript to article in journal. QC 20120605

Available from: 2011-11-23 Created: 2011-11-23 Last updated: 2017-12-08Bibliographically approved
In thesis
1. Spectral Accuracy in Fast Ewald Methods and Topics in Fluid Interface Simulation
Open this publication in new window or tab >>Spectral Accuracy in Fast Ewald Methods and Topics in Fluid Interface Simulation
2011 (English)Doctoral thesis, comprehensive summary (Other academic)
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.

Place, publisher, year, edition, pages
Stockholm: KTH Royal Institute of Technology, 2011. xv, 104 p.
Series
Trita-CSC-A, ISSN 1653-5723 ; 2011:19
National Category
Computational Mathematics
Identifiers
urn:nbn:se:kth:diva-48805 (URN)978-91-7501-195-0 (ISBN)
Public defence
2011-12-16, Salongen, KTHB, Osquars backe 25, Stockholm, 10:00 (English)
Opponent
Supervisors
Funder
Swedish e‐Science Research Center
Note
QC 20111125Available from: 2011-11-25 Created: 2011-11-23 Last updated: 2012-05-24Bibliographically approved

Open Access in DiVA

No full text

Other links

Publisher's full textScopusarXiv:1109.1667

Search in DiVA

By author/editor
Lindbo, DagTornberg, Anna-Karin
By organisation
Numerical Analysis, NA (closed 2012-06-30)
In the same journal
Journal of Chemical Physics
Physical SciencesComputational Mathematics

Search outside of DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 57 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