Change search
Refine search result
1 - 12 of 12
CiteExportLink to result list
Permanent 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
Rows per page
  • 5
  • 10
  • 20
  • 50
  • 100
  • 250
Sort
  • Standard (Relevance)
  • Author A-Ö
  • Author Ö-A
  • Title A-Ö
  • Title Ö-A
  • Publication type A-Ö
  • Publication type Ö-A
  • Issued (Oldest first)
  • Issued (Newest first)
  • Created (Oldest first)
  • Created (Newest first)
  • Last updated (Oldest first)
  • Last updated (Newest first)
  • Disputation date (earliest first)
  • Disputation date (latest first)
  • Standard (Relevance)
  • Author A-Ö
  • Author Ö-A
  • Title A-Ö
  • Title Ö-A
  • Publication type A-Ö
  • Publication type Ö-A
  • Issued (Oldest first)
  • Issued (Newest first)
  • Created (Oldest first)
  • Created (Newest first)
  • Last updated (Oldest first)
  • Last updated (Newest first)
  • Disputation date (earliest first)
  • Disputation date (latest first)
Select
The maximal number of hits you can export is 250. When you want to export more records please use the Create feeds function.
  • 1.
    af Klinteberg, Ludvig
    et al.
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Saffar Shamshirgar, Davoud
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Tornberg, Anna-Karin
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Fast Ewald summation for free-space Stokes potentials2017In: Research in the Mathematical Sciences, ISSN 2197-9847, Vol. 4, no 1Article in journal (Refereed)
    Abstract [en]

    We present a spectrally accurate method for the rapid evaluation of free-space Stokes potentials, i.e., sums involving a large number of free space Green’s functions. We consider sums involving stokeslets, stresslets and rotlets that appear in boundary integral methods and potential methods for solving Stokes equations. The method combines the framework of the Spectral Ewald method for periodic problems (Lindbo and Tornberg in J Comput Phys 229(23):8994–9010, 2010. doi: 10.1016/j.jcp.2010.08.026 ), with a very recent approach to solving the free-space harmonic and biharmonic equations using fast Fourier transforms (FFTs) on a uniform grid (Vico et al. in J Comput Phys 323:191–203, 2016. doi: 10.1016/j.jcp.2016.07.028 ). Convolution with a truncated Gaussian function is used to place point sources on a grid. With precomputation of a scalar grid quantity that does not depend on these sources, the amount of oversampling of the grids with Gaussians can be kept at a factor of two, the minimum for aperiodic convolutions by FFTs. The resulting algorithm has a computational complexity of $$O(N \log N)$$ O ( N log N ) for problems with N sources and targets. Comparison is made with a fast multipole method to show that the performance of the new method is competitive.

  • 2.
    af Klinteberg, Ludvig
    et al.
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Tornberg, Anna-Karin
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Error estimation for quadrature by expansion in layer potential evaluation2017In: Advances in Computational Mathematics, ISSN 1019-7168, E-ISSN 1572-9044, Vol. 43, no 1, p. 195-234Article in journal (Refereed)
    Abstract [en]

    In boundary integral methods it is often necessary to evaluate layer potentials on or close to the boundary, where the underlying integral is difficult to evaluate numerically. Quadrature by expansion (QBX) is a new method for dealing with such integrals, and it is based on forming a local expansion of the layer potential close to the boundary. In doing so, one introduces a new quadrature error due to nearly singular integration in the evaluation of expansion coefficients. Using a method based on contour integration and calculus of residues, the quadrature error of nearly singular integrals can be accurately estimated. This makes it possible to derive accurate estimates for the quadrature errors related to QBX, when applied to layer potentials in two and three dimensions. As examples we derive estimates for the Laplace and Helmholtz single layer potentials. These results can be used for parameter selection in practical applications.

  • 3.
    Fryklund, Fredrik
    et al.
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.).
    Lehto, Erik
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Tornberg, Anna-Karin
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Partition of unity extension of functions on complex domains2018In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 375, p. 57-79Article in journal (Refereed)
    Abstract [en]

    We introduce an efficient algorithm, called partition of unity extension or PUX, to construct an extension of desired regularity of a function given on a complex multiply connected domain in 2D. Function extension plays a fundamental role in extending the applicability of boundary integral methods to inhomogeneous partial differential equations with embedded domain techniques. Overlapping partitions are placed along the boundaries, and a local extension of the function is computed on each patch using smooth radial basis functions; a trivially parallel process. A partition of unity method blends the local extrapolations into a global one, where weight functions impose compact support. The regularity of the extended function can be controlled by the construction of the partition of unity function. We evaluate the performance of the PUX method in the context of solving the Poisson equation on multiply connected domains using a boundary integral method and a spectral solver. With a suitable choice of parameters the error converges as a tenth order method down to 10−14.

  • 4.
    Pålsson, Sara
    et al.
    KTH, School of Engineering Sciences (SCI), Centres, Linné Flow Center, FLOW.
    Siegel, Michael
    New Jersey Inst Technol, Dept Math Sci, Newark, NJ 07102 USA..
    Tornberg, Anna-Karin
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Simulation and validation of surfactant-laden drops in two-dimensional Stokes flow2019In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 386, p. 218-247Article in journal (Refereed)
    Abstract [en]

    Performing highly accurate simulations of droplet systems is a challenging problem. This is primarily due to the interface dynamics which is complicated further by the addition of surfactants. This paper presents a boundary integral method for computing the evolution of surfactant-covered droplets in 2D Stokes flow. The method has spectral accuracy in space and the adaptive time-stepping scheme allows for control of the temporal errors. Previously available semi-analytical solutions (based on conformal-mapping techniques) are extended to include surfactants, and a set of algorithms is introduced to detail their evaluation. These semi-analytical solutions are used to validate and assess the accuracy of the boundary integral method, and it is demonstrated that the presented method maintains its high accuracy even when droplets are in close proximity. 

  • 5.
    Saffar Shamshirgar, Davood
    et al.
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Hess, Berk
    KTH, Centres, Science for Life Laboratory, SciLifeLab.
    Tornberg, Anna-Karin
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    A comparison of the Spectral Ewald and Smooth Particle Mesh Ewald methods in GROMACSManuscript (preprint) (Other academic)
    Abstract [en]

    The smooth particle mesh Ewald (SPME) method is an FFT based methodfor the fast evaluation of electrostatic interactions under periodic boundaryconditions. A highly optimized implementation of this method is availablein GROMACS, a widely used software for molecular dynamics simulations.In this article, we compare a more recent method from the same family ofmethods, the spectral Ewald (SE) method, to the SPME method in termsof performance and efficiency. We consider serial and parallel implementa-tions of both methods for single and multiple core computations on a desktopmachine as well as the Beskow supercomputer at KTH Royal Institute ofTechnology. The implementation of the SE method has been well optimized,however not yet comparable to the level of the SPME implementation thathas been improved upon for many years. We show that the SE method isvery efficient whenever used to achieve high accuracy and that it already atthis level of optimization can be competitive for low accuracy demands.

  • 6.
    Saffar Shamshirgar, Davood
    et al.
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Hess, Berk
    KTH, Centres, Science for Life Laboratory, SciLifeLab.
    Yokota, Rio
    Tornberg, Anna-Karin
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Regularized FMM for MD simulationsManuscript (preprint) (Other academic)
    Abstract [en]

    A regularized fast multipole method (FMM) which approximately conserves the total energy in Molecular dynamics (MD) simulations is presented. The new algorithm introduces a regularization which eliminates the discontinuity inherent in the FMM. This allows us to use FMM in simulations as a substitute for widely used FFT based methods. For a system of N particles, the computational complexity of the resulting method is still of order N though with a larger constant compared to the plain FMM. Numerical examples are provided to confirm that the new algorithm improves the accuracy and approximately conserves the long term total energy.

  • 7.
    Saffar Shamshirgar, Davood
    et al.
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Tornberg, Anna-Karin
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    A fast multipole method for evaluating exponential integral type kernelsManuscript (preprint) (Other academic)
    Abstract [en]

    We present a fast multipole method for evaluation of sums with exponential  integral type kernels. These sums appear while solving free space Poisson problems in two dimensions and in the derivation of 1d-periodic Ewald sums. The presented method uses recurrence relations to derive multipole expansions for computing interactions between particles and far clusters.

  • 8.
    Saffar Shamshirgar, Davood
    et al.
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Tornberg, Anna-Karin
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Fast Ewald summation for electrostatic potentials with arbitrary periodicityManuscript (preprint) (Other academic)
    Abstract [en]

    A unified treatment for fast and spectrally accurate evaluation of electrostatic potentials subject to periodic boundary conditions in any or none of the three space dimensions is presented. Ewald decomposition is used to split the problem into a real space and a Fourier space part, and the FFT based Spectral Ewald (SE) method is used to accelerate the computation of the latter. A key component in the unified treatment is an FFT based solution technique for the free-space Poisson problem in three, two or one dimensions, depending on the number of non-periodic directions. The cost of calculations is furthermore reduced by employing an adaptive FFT for the doubly and singly periodic cases, allowing for different local upsampling rates. The SE method will always be most efficient for the triply periodic case as the cost for computing FFTs will be the smallest, whereas the computational cost for the rest of the algorithm is essentially independent of the periodicity. We show that the cost of removing periodic boundary conditions from one or two directions out of three will only marginally increase the total run time. Our comparisons also show that the computational cost of the SE method for the free-space case is typically about four times more expensive as compared to the triply periodic case.

    The Gaussian window function previously used in the SE method, is here compared to an approximation of the Kaiser-Bessel window function, recently introduced. With a carefully tuned shape parameter that is selected based on an error estimate for this new window function, runtimes for the SE method can be further reduced.

  • 9. Siegel, M.
    et al.
    Tornberg, Anna-Karin
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.). KTH, School of Engineering Sciences (SCI), Centres, Linné Flow Center, FLOW. KTH, Centres, SeRC - Swedish e-Science Research Centre.
    A local target specific quadrature by expansion method for evaluation of layer potentials in 3D2018In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 364, p. 365-392Article in journal (Refereed)
    Abstract [en]

    Accurate evaluation of layer potentials is crucial when boundary integral equation methods are used to solve partial differential equations. Quadrature by expansion (QBX) is a recently introduced method that can offer high accuracy for singular and nearly singular integrals, using truncated expansions to locally represent the potential. The QBX method is typically based on a spherical harmonics expansion which when truncated at order p has O(p2) terms. This expansion can equivalently be written with p terms, however paying the price that the expansion coefficients will depend on the evaluation/target point. Based on this observation, we develop a target specific QBX method, and apply it to Laplace's equation on multiply-connected domains. The method is local in that the QBX expansions only involve information from a neighborhood of the target point. An analysis of the truncation error in the QBX expansions is presented, practical parameter choices are discussed and the method is validated and tested on various problems.

  • 10.
    Sorgentone, Chiara
    et al.
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.). KTH, School of Engineering Sciences (SCI), Centres, Linné Flow Center, FLOW. KTH, Centres, SeRC - Swedish e-Science Research Centre.
    Tornberg, Anna-Karin
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.). KTH, School of Engineering Sciences (SCI), Centres, Linné Flow Center, FLOW. KTH, Centres, SeRC - Swedish e-Science Research Centre.
    A highly accurate boundary integral equation method for surfactant-laden drops in 3D2018In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 360, p. 167-191Article in journal (Refereed)
    Abstract [en]

    The presence of surfactants alters the dynamics of viscous drops immersed in an ambient viscous fluid. This is specifically true at small scales, such as in applications of droplet based microfluidics, where the interface dynamics become of increased importance. At such small scales, viscous forces dominate and inertial effects are often negligible. Considering Stokes flow, a numerical method based on a boundary integral formulation is presented for simulating 3D drops covered by an insoluble surfactant. The method is able to simulate drops with different viscosities and close interactions, automatically controlling the time step size and maintaining high accuracy also when substantial drop deformation appears. To achieve this, the drop surfaces as well as the surfactant concentration on each surface are represented by spherical harmonics expansions. A novel reparameterization method is introduced to ensure a high-quality representation of the drops also under deformation, specialized quadrature methods for singular and nearly singular integrals that appear in the formulation are evoked and the adaptive time stepping scheme for the coupled drop and surfactant evolution is designed with a preconditioned implicit treatment of the surfactant diffusion.

  • 11.
    Sorgentone, Chiara
    et al.
    KTH, School of Engineering Sciences (SCI), Centres, Linné Flow Center, FLOW. KTH, School of Engineering Sciences (SCI), Mathematics (Dept.).
    Tornberg, Anna-Karin
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.). KTH, School of Engineering Sciences (SCI), Centres, Linné Flow Center, FLOW.
    Vlahovska, Petia M.
    Northwestern Univ, Engn Sci & Appl Math, Evanston, IL 60208 USA..
    A 3D boundary integral method for the electrohydrodynamics of surfactant-covered drops2019In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 389, p. 111-127Article in journal (Refereed)
    Abstract [en]

    We present a highly accurate numerical method based on a boundary integral formulation and the leaky dielectric model to study the dynamics of surfactant-covered drops in the presence of an applied electric field. The method can simulate interacting 3D drops (no axisymmetric simplification) in close proximity, can consider different viscosities, is adaptive in time and able to handle substantial drop deformation. For each drop global representations of the variables based on spherical harmonics expansions are used and the spectral accuracy is achieved by designing specific numerical tools: a specialized quadrature method for the singular and nearly singular integrals that appear in the formulation, a general preconditioner for the implicit treatment of the surfactant diffusion and a reparametrization procedure able to ensure a high-quality representation of the drops also under deformation. Our numerical method is validated against theoretical, numerical and experimental results available in the literature, as well as a new second-order theory developed for a surfactant-laden drop placed in a quadrupole electric field.

  • 12.
    Srinivasan, Shriram
    et al.
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Tornberg, Anna-Karin
    KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.
    Fast Ewald summation for Green's functions of Stokes flow in a half-space2018In: RESEARCH IN THE MATHEMATICAL SCIENCES, ISSN 2197-9847, Vol. 5, article id 35Article in journal (Refereed)
    Abstract [en]

    Recently, Gimbutas et al. (J Fluid Mech, 2015. https://doi.org/10.1017/jfm.2015.302) derived an elegant representation for the Green's functions of Stokes flow in a half-space. We present a fast summation method for sums involving these half-space Green's functions (stokeslets, stresslets and rotlets) that consolidates and builds on the work by Klinteberg et al. (Res Math Sci 4(1): 1, 2017. https://doi.org/10.1186/s40687-016-0092-7) for the corresponding free-space Green's functions. The fast method is based on two main ingredients: The Ewald decomposition and subsequent use of FFTs. The Ewald decomposition recasts the sum into a sum of two exponentially decaying series: one in real space (short-range interactions) and one in Fourier space (long-range interactions) with the convergence of each series controlled by a common parameter. The evaluation of short-range interactions is accelerated by restricting computations to neighbours within a specified distance, while the use of FFTs accelerates the computations in Fourier space thus accelerating the overall sum. We demonstrate that while the method incurs extra costs for the half-space in comparison with the free-space evaluation, greater computational savings is also achieved when compared to their respective direct sums.

1 - 12 of 12
CiteExportLink to result list
Permanent 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