Hydrodynamics of Suspensions of Passive and Active Rigid Particles: A Rigid Multiblob Approach

We develop a rigid multiblob method for numerically solving the mobility problem for suspensions of passive and active rigid particles of complex shape in Stokes flow in unconfined, partially confined, and fully confined geometries. As in a number of existing methods, we discretize rigid bodies using a collection of minimally-resolved spherical blobs constrained to move as a rigid body, to arrive at a potentially large linear system of equations for the unknown Lagrange multipliers and rigid-body motions. Here we develop a block-diagonal preconditioner for this linear system and show that a standard Krylov solver converges in a modest number of iterations that is essentially independent of the number of particles. For unbounded suspensions and suspensions sedimented against a single no-slip boundary, we rely on existing analytical expressions for the Rotne-Prager tensor combined with a fast multipole method or a direct summation on a Graphical Processing Unit to obtain an simple yet efficient and scalable implementation. For fully confined domains, such as periodic suspensions or suspensions confined in slit and square channels, we extend a recently-developed rigid-body immersed boundary method to suspensions of freely-moving passive or active rigid particles at zero Reynolds number. We demonstrate that the iterative solver for the coupled fluid and rigid body equations converges in a bounded number of iterations regardless of the system size. We optimize a number of parameters in the iterative solvers and apply our method to a variety of benchmark problems to carefully assess the accuracy of the rigid multiblob approach as a function of the resolution. We also model the dynamics of colloidal particles studied in recent experiments, such as passive boomerangs in a slit channel, as well as a pair of non-Brownian active nanorods sedimented against a wall.Comment: Under revision in CAMCOS, Nov 201

Computation of Steady Incompressible Flows in Unbounded Domains

In this study we revisit the problem of computing steady Navier-Stokes flows in two-dimensional unbounded domains. Precise quantitative characterization of such flows in the high-Reynolds number limit remains an open problem of theoretical fluid dynamics. Following a review of key mathematical properties of such solutions related to the slow decay of the velocity field at large distances from the obstacle, we develop and carefully validate a spectrally-accurate computational approach which ensures the correct behavior of the solution at infinity. In the proposed method the numerical solution is defined on the entire unbounded domain without the need to truncate this domain to a finite box with some artificial boundary conditions prescribed at its boundaries. Since our approach relies on the streamfunction-vorticity formulation, the main complication is the presence of a discontinuity in the streamfunction field at infinity which is related to the slow decay of this field. We demonstrate how this difficulty can be overcome by reformulating the problem using a suitable background "skeleton" field expressed in terms of the corresponding Oseen flow combined with spectral filtering. The method is thoroughly validated for Reynolds numbers spanning two orders of magnitude with the results comparing favourably against known theoretical predictions and the data available in the literature.Comment: 39 pages, 12 figures, accepted for publication in "Computers and Fluids

Baraj-rezervuar Etkileşim Problemlerinde Sonsuz Rezervuar Modellemesi İçin Yari Analitik Bir Yaklaşim

Konferans Bildirisi -- Teorik ve Uygulamalı Mekanik Türk Milli Komitesi, 2008Conference Paper -- Theoretical and Applied Mechanical Turkish National Committee, 2008Bu çalışmada, baraj-rezervuar etkileşim problemlerinde sonsuz rezervuar modellemesi için gerçek bir sınır şartı türetilmiş ve bir sonlu eleman programına monte edilmiştir. Geliştirilen sınır şartı baraj sıvısının viskoz olmadığı ve sıkıştırılamaz olduğu, ayrıca rezervuarın da sonsuz olduğu durum için geliştirilmiştir. Bu sınır şartı baraj-rezervuar etkileşiminde rezervuarın kesim yüzeyinde uygulanmıştır. Yöntemin yarı analitik bir yaklaşım olması, analizi yapılan durumlarda gerçek sonuçlara çok yakın sonuçlar elde edilmesini sağlamıştır.In this study, an exact truncation boundary condition is derived and implemented in a finite element code for the analysis of dam-reservoir interaction for incompressible, inviscid and unbounded fluid domains. This boundary condition is used at the truncation surface for dam-reservoir interaction considered. Since the method is a semi-analytical approach, obtained results are close to the exact solutions for the cases analyzed

Multi-resolution lattice Green's function method for incompressible flows

We propose a multi-resolution strategy that is compatible with the lattice Green's function (LGF) technique for solving viscous, incompressible flows on unbounded domains. The LGF method exploits the regularity of a finite-volume scheme on a formally unbounded Cartesian mesh to yield robust and computationally efficient solutions. The original method is spatially adaptive, but challenging to integrate with embedded mesh refinement as the underlying LGF is only defined for a fixed resolution. We present an ansatz for adaptive mesh refinement, where the solutions to the pressure Poisson equation are approximated using the LGF technique on a composite mesh constructed from a series of infinite lattices of differing resolution. To solve the incompressible Navier-Stokes equations, this is further combined with an integrating factor for the viscous terms and an appropriate Runge-Kutta scheme for the resulting differential-algebraic equations. The parallelized algorithm is verified through with numerical simulations of vortex rings, and the collision of vortex rings at high Reynolds number is simulated to demonstrate the reduction in computational cells achievable with both spatial and refinement adaptivity