BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Anna-Karin Tornberg

Anna-Karin Tornberg appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysisphysics.flu-dynphysics.comp-phcond-mat.softmath-phmath.APmath.MP

Catalog footprint

What is connected

14works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

An adaptive kernel-split quadrature method for parameter-dependent layer potentials

Panel-based, kernel-split quadrature is currently one of the most efficient methods available for accurate evaluation of singular and nearly singular layer potentials in two dimensions. However, it can fail completely for the layer potentials belonging to the modified Helmholtz, biharmonic and Stokes equations. These equations depend on a parameter, denoted $α$, and kernel-split quadrature loses its accuracy rapidly when this parameter grows beyond a certain threshold. This paper describes an algorithm that remedies this problem, using per-target adaptive sampling of the source geometry. The refinement is carried out through recursive bisection, with a carefully selected rule set. This maintains accuracy for a wide range of the parameter $α$, at an increased cost that scales as $\logα$. Using this algorithm allows kernel-split quadrature to be both accurate and efficient for a much wider range of problems than previously possible.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quadrature error estimates for layer potentials evaluated near curved surfaces in three dimensions

The quadrature error associated with a regular quadrature rule for evaluation of a layer potential increases rapidly when the evaluation point approaches the surface and the integral becomes nearly singular. Error estimates are needed to determine when the accuracy is insufficient and a more costly special quadrature method should be utilized. The final result of this paper are such quadrature error estimates for the composite Gauss-Legendre rule and the global trapezoidal rule, when applied to evaluate layer potentials defined over smooth curved surfaces in R^3. The estimates have no unknown coefficients and can be efficiently evaluated given the discretization of the surface, invoking a local one-dimensional root-finding procedure. They are derived starting with integrals over curves, using complex analysis involving contour integrals, residue calculus and branch cuts. By complexifying the parameter plane, the theory can be used to derive estimates also for curves in in R^3. These results are then used in the derivation of the estimates for integrals over surfaces. In this procedure, we also obtain error estimates for layer potentials evaluated over curves in R^2. Such estimates combined with a local root-finding procedure for their evaluation were earlier derived for the composite Gauss-Legendre rule for layer potentials written on complex form [4]. This is here extended to provide quadrature error estimates for both complex and real formulations of layer potentials, both for the Gauss-Legendre and the trapezoidal rule. Numerical examples are given to illustrate the performance of the quadrature error estimates. The estimates for integration over curves are in many cases remarkably precise, and the estimates for curved surfaces in R^3 are also sufficiently precise, with sufficiently low computational cost, to be practically useful.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fast Ewald summation for electrostatic potentials with arbitrary periodicity

A unified treatment for fast and spectrally accurate evaluation of electrostatic potentials subject to periodic boundary conditions in any or none of the three spatial 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 computational cost is furthermore reduced by employing an adaptive FFT for the doubly and singly periodic cases, allowing for different local upsampling factors. The SE method will always be most efficient for the triply periodic case as the cost of computing FFTs will then be the smallest, whereas the computational cost of the rest of the algorithm is essentially independent of periodicity. We show that the cost of removing periodic boundary conditions from one or two directions out of three will only moderately increase the total runtime. Our comparisons also show that the computational cost of the SE method in the free-space case is around four times that of the triply periodic case. The Gaussian window function previously used in the SE method, is here compared to a piecewise polynomial approximation of the Kaiser-Bessel window function. 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. Furthermore, we consider different methods for computing the force, and compare the runtime of the SE method with that of the Fast Multipole Method.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Parabolic velocity profile causes drift of inertial prolate spheroids -- but gravity is stronger

Motion of elongated particles in shear is studied. In applications where particles are much heavier than the carrying fluid, e.g. aerosols, the influence of particle inertia dominates the particle dynamics. Assuming that the particle only experiences a local linear velocity profile, its rotational and translational motion are independent. However, we show that quadratic terms of the local velocity profile combined with particle inertia cause a lateral drift of prolate spheroidal particles. We find that this drift is maximal when particle inertial forces are of the same order of magnitude as viscous forces, and that both extremely light and extremely heavy particles have negligible drift. In the non-inertial case, the particle rotates according to the local linear velocity profile, with each instantaneous orientation corresponding to a velocity that gives zero force on the particle. This results in a translational motion in the flow direction with periodic velocity fluctuations. With added particle inertia, the particle is slow to react to the surrounding fluid motion and the particle will switch between slower and faster rotation compared to the zero-force solution. The final motion that gives zero integrated force over a rotational period, is a motion with a lateral drift. We show that this drift is purely an effect of the non-sphericity of the particle and its translational inertia, while rotational inertia is negligible. Finally, although this inertial drift will contribute to the lateral motion of heavy elongated particles in channel flow, sedimentation due to gravity will dominate in any practical application on earth.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An accurate integral equation method for Stokes flow with piecewise smooth boundaries

Two-dimensional Stokes flow through a periodic channel is considered. The channel walls need only be Lipschitz continuous, in other words they are allowed to have corners. Boundary integral methods are an attractive tool for numerically solving the Stokes equations, as the partial differential equation can be reformulated into an integral equation that must be solved only over the boundary of the domain. When the boundary is at least C 1 smooth, the boundary integral kernel is a compact operator, and traditional Nyström methods can be used to obtain highly accurate solutions. In the case of Lipschitz continuous boundaries, however, obtaining accurate solutions using the standard Nyström method can require high resolution. We adapt a technique known as recursively compressed inverse preconditioning to accurately solve the Stokes equations without requiring any more resolution than is needed to resolve the boundary. Combined with a periodic fast summation method we construct a method that is O(N log N ) where N is the number of quadrature points on the boundary. We demonstrate the robustness of this method by extending an existing boundary integral method for viscous drops to handle the movement of drops near corners.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Highly accurate special quadrature methods for Stokesian particle suspensions in confined geometries

Boundary integral methods are highly suited for problems with complicated geometries, but require special quadrature methods to accurately compute the singular and nearly singular layer potentials that appear in them. This paper presents a boundary integral method that can be used to study the motion of rigid particles in three-dimensional periodic Stokes flow with confining walls. A centrepiece of our method is the highly accurate special quadrature method, which is based on a combination of upsampled quadrature and quadrature by expansion (QBX), accelerated using a precomputation scheme. The method is demonstrated for rodlike and spheroidal particles, with the confining geometry given by a pipe or a pair of flat walls. A parameter selection strategy for the special quadrature method is presented and tested. Periodic interactions are computed using the Spectral Ewald (SE) fast summation method, which allows our method to run in O(n log n) time for n grid points, assuming the number of geometrical objects grows while the grid point concentration is kept fixed.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Adaptive quadrature by expansion for layer potential evaluation in two dimensions

When solving partial differential equations using boundary integral equation methods, accurate evaluation of singular and nearly singular integrals in layer potentials is crucial. A recent scheme for this is quadrature by expansion (QBX), which solves the problem by locally approximating the potential using a local expansion centered at some distance from the source boundary. In this paper we introduce an extension of the QBX scheme in 2D denoted AQBX - adaptive quadrature by expansion - which combines QBX with an algorithm for automated selection of parameters, based on a target error tolerance. A key component in this algorithm is the ability to accurately estimate the numerical errors in the coefficients of the expansion. Combining previous results for flat panels with a procedure for taking the panel shape into account, we derive such error estimates for arbitrarily shaped boundaries in 2D that are discretized using panel-based Gauss-Legendre quadrature. Applying our scheme to numerical solutions of Dirichlet problems for the Laplace and Helmholtz equations, and also for solving these equations, we find that the scheme is able to satisfy a given target tolerance to within an order of magnitude, making it useful for practical applications. This represents a significant simplification over the original QBX algorithm, in which choosing a good set of parameters can be hard.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A fast integral equation method for solid particles in viscous flow using quadrature by expansion

Boundary integral methods are advantageous when simulating viscous flow around rigid particles, due to the reduction in number of unknowns and straightforward handling of the geometry. In this work we present a fast and accurate framework for simulating spheroids in periodic Stokes flow, which is based on the completed double layer boundary integral formulation. The framework implements a new method known as quadrature by expansion (QBX), which uses surrogate local expansions of the layer potential to evaluate it to very high accuracy both on and off the particle surfaces. This quadrature method is accelerated through a newly developed precomputation scheme. The long range interactions are computed using the spectral Ewald (SE) fast summation method, which after integration with QBX allows the resulting system to be solved in M log M time, where M is the number of particles. This framework is suitable for simulations of large particle systems, and can be used for studying e.g. porous media models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An accurate integral equation method for simulating multi-phase Stokes flow

We introduce a numerical method based on an integral equation formulation for simulating drops in viscous fluids in the plane. It builds upon the method introduced by Kropinski in 2001, but improves on it by adding an interpolatory quadrature approach for handling near-singular integrals. Such integrals typically arise when drop boundaries come close to one another, and are difficult to compute accurately using standard quadrature rules. Adapting the interpolatory quadrature method introduced by Helsing and Ojala in 2008 to the current application, very general drop configurations can be handled while still maintaining stability and high accuracy. The performance of the new method is demonstrated by some challenging numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Ewald sums for singly, doubly and triply periodic electrostatic systems

When evaluating the electrostatic potential, periodic boundary conditions in one, two or three of the spatial dimensions are often needed for different applications. The triply periodic Ewald summation formula is classical, and Ewald summation formulas for the other two cases have also been derived. In this paper, derivations of the Ewald sums in the doubly and singly periodic cases are presented in a uniform framework based on Fourier analysis, which also yields a natural starting point for FFT-based fast summation methods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Interface tracking using patches

A new method for interface tracking is presented. The interface representation, based on domain decomposition, provides the interface location explicitly, yet is Eulerian. This allows for well established finite difference methods on uniform grids to be used for the numerics of advecting the interface and other computations. CFL-stable and second order accurate explicit and implicit time-stepping methods are derived. Numerical results are given to substantiate stated convergence properties, as well as convergence in interface curvature and mass conservation. Our method is applied to a boundary integral formulation for Stokes flow, and the resulting integrals are analyzed and treated numerically to second order accuracy. Finally, we embed our method in the familiar immersed boundary- and immersed interface methods for two-phase Navier-Stokes flow.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On diffuse interface modeling and simulation of surfactants in two-phase fluid flow

An existing phase-field model of two immiscible fluids with a single soluble surfactant present is discussed in detail. We analyze the well-posedness of the model and provide strong evidence that it is mathematically ill-posed for a large set of physically relevant parameters. As a consequence, critical modifications to the model are suggested that substantially increase the domain of validity. Carefully designed numerical simulations offer informative demonstrations as to the sharpness of our theoretical results and the qualities of the physical model. A fully coupled hydrodynamic test-case demonstrates the potential to capture also non-trivial effects on the overall flow.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fast and spectrally accurate Ewald summation for 2-periodic electrostatic systems

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 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fast and spectrally accurate summation of 2-periodic Stokes potentials

We derive a Ewald decomposition for the Stokeslet in planar periodicity and a novel PME-type O(N log N) method for the fast evaluation of the resulting sums. The decomposition is the natural 2P counterpart to the classical 3P decomposition by Hasimoto, and is given in an explicit form not found in the literature. Truncation error estimates are provided to aid in selecting parameters. The fast, PME-type, method appears to be the first fast method for computing Stokeslet Ewald sums in planar periodicity, and has three attractive properties: it is spectrally accurate; it uses the minimal amount of memory that a gridded Ewald method can use; and provides clarity regarding numerical errors and how to choose parameters. Analytical and numerical results are give to support this. We explore the practicalities of the proposed method, and survey the computational issues involved in applying it to 2-periodic boundary integral Stokes problems.

0 citations0 reviewsNo code linkedOpen access