Source author record

Joar Bagge

Joar Bagge appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis cond-mat.soft physics.comp-ph physics.flu-dyn

Catalog footprint

What is connected

3works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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.

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.

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.