BZPEERReading, trust and collaboration for research
Menu

Discover

SearchFind papers, people, topics and opportunities from one place.GraphOpen maps around papers, topics and saved views.

Workspaces

HomeOpen the ranked research feed.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Source author record

Bryan Quaife

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

ResearcherUnclaimed source record
math.NAcond-mat.softphysics.flu-dyn

Catalog footprint

What is connected

9works
3topics
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 work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this map preview

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

Published work

9 published item(s)

preprint2022arXiv

Effects of Tunable Hydrophobicity on the Collective Hydrodynamics of Janus Particles under Flows

Active colloidal systems with non-equilibrium self-organization is a long-standing, challenging area in biology. To understand how hydrodynamic flow may be used to actively control self-assembly of Janus particles (JPs), we use a model recently developed for the many-body hydrodynamics of amphiphilic JPs suspended in a viscous background flow (JFM, 941, 2022). We investigate how various morphologies arise from tuning the hydrophobic distribution of the JP-solvent interface. We find JPs assembled into uni-lamella, multi-lamella and striated structures. To introduce dynamics, we include a linear shear flow and a steady Taylor-Green mixing flow, and measure the collective dynamics of JP particles in terms of their (a) free energy from the hydrophobic interactions between the JPs, (b) order parameter for the ordering of JPs in terms of alignment of their directors, and (c) strain parameter that captures the deformation in the assembly. We characterize the effective material properties of the JP structures and find that the uni-lamellar structures increases orientation order under shear flow, the multilamellar structure behaves as a shear thinning fluid, and the striated structure possesses a yield stress. These numerical results provide insights into dynamic control of non-equilibrium active biological systems with similar self-organization.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An efficient preconditioner for the fast simulation of a 2D Stokes flow in porous media

We consider an efficient preconditioner for boundary integral equation (BIE) formulations of the two-dimensional Stokes equations in porous media. While BIEs are well-suited for resolving the complex porous geometry, they lead to a dense linear system of equations that is computationally expensive to solve for large problems. This expense is further amplified when a significant number of iterations is required in an iterative Krylov solver such as GMRES. In this paper, we apply a fast inexact direct solver, the inverse fast multipole method (IFMM), as an efficient preconditioner for GMRES. This solver is based on the framework of $\mathcal{H}^{2}$-matrices and uses low-rank compressions to approximate certain matrix blocks. It has a tunable accuracy $\varepsilon$ and a computational cost that scales as $\mathcal{O} (N \log^2 1/\varepsilon)$. We discuss various numerical benchmarks that validate the accuracy and confirm the efficiency of the proposed method. We demonstrate with several types of boundary conditions that the preconditioner is capable of significantly accelerating the convergence of GMRES when compared to a simple block-diagonal preconditioner, especially for pipe flow problems involving many pores.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Adaptive Time Stepping for Vesicle Suspensions

We present an adaptive arbitrary-order accurate time-stepping numerical scheme for the flow of vesicles suspended in Stokesian fluids. Our scheme can be summarized as an approximate implicit spectral deferred correction (SDC) method. Applying a textbook fully implicit SDC scheme to vesicle flows is prohibitively expensive. For this reason we introduce several approximations. Our scheme is based on a semi-implicit linearized low-order time stepping method. (Our discretization is spectrally accurate in space.) We also use invariant properties of vesicle flows, constant area and boundary length in two dimensions, to reduce the computational cost of error estimation for adaptive time stepping. We present results in two dimensions for single-vesicle flows, constricted geometry flows, converging flows, and flows in a Couette apparatus. We experimentally demonstrate that the proposed scheme enables automatic selection of the step size and high-order accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-order adaptive time stepping for vesicle suspensions with viscosity contrast

We construct a high-order adaptive time stepping scheme for vesicle suspensions with viscosity contrast. The high-order accuracy is achieved using a spectral deferred correction (SDC) method, and adaptivity is achieved by estimating the local truncation error with the numerical error of physically constant values. Numerical examples demonstrate that our method can handle suspensions with vesicles that are tumbling, tank-treading, or both. Moreover, we demonstrate that a user-prescribed tolerance can be automatically achieved for simulations with long time horizons.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-volume fraction simulations of two-dimensional vesicle suspensions

We consider numerical algorithms for the simulation of the rheology of two-dimensional vesicles suspended in a viscous Stokesian fluid. The vesicle evolution dynamics is governed by hydrodynamic and elastic forces. The elastic forces are due to local inextensibility of the vesicle membrane and resistance to bending. Numerically resolving vesicle flows poses several challenges. For example, we need to resolve moving interfaces, address stiffness due to bending, enforce the inextensibility constraint, and efficiently compute the (non-negligible) long-range hydrodynamic interactions. Our method is based on the work of {\em Rahimian, Veerapaneni, and Biros, "Dynamic simulation of locally inextensible vesicles suspended in an arbitrary two-dimensional domain, a boundary integral method", Journal of Computational Physics, 229 (18), 2010}. It is a boundary integral formulation of the Stokes equations coupled to the interface mass continuity and force balance. We extend the algorithms presented in that paper to increase the robustness of the method and enable simulations with concentrated suspensions. In particular, we propose a scheme in which both intra-vesicle and inter-vesicle interactions are treated semi-implicitly. In addition we use special integration for near-singular integrals and we introduce a spectrally accurate collision detection scheme. We test the proposed methodologies on both unconfined and confined flows for vesicles whose internal fluid may have a viscosity contrast with the bulk medium. Our experiments demonstrate the importance of treating both intra-vesicle and inter-vesicle interactions accurately.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Integral equation methods for the Yukawa-Beltrami equation on the sphere

An integral equation method for solving the Yukawa-Beltrami equation on a multiply-connected sub-manifold of the unit sphere is presented. A fundamental solution for the Yukawa-Beltrami operator is constructed. This fundamental solution can be represented by conical functions. Using a suitable representation formula, a Fredholm equation of the second kind with a compact integral operator needs to be solved. The discretization of this integral equation leads to a linear system whose condition number is bounded independent of the size of the system. Several numerical examples exploring the properties of this integral equation are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On preconditioners for the Laplace double-layer in 2D

The discretization of the double-layer potential integral equation for the interior Dirichlet Laplace problem in a domain with smooth boundary results in a linear system that has a bounded condition number. Thus, the number of iterations required for the convergence of a Krylov method is, asymptotically, independent of the discretization size $N$. Using the Fast Multipole Method (FMM) to accelerate the matrix-vector products, we obtain an optimal $\mathcal{O}(N)$ solver. In practice, however, when the geometry is complicated, the number of Krylov iterations can be quite large---to the extend that necessitates the use of preconditioning. We summarize the different methodologies that have appeared in the literature (single-grid, multigrid, approximate sparse inverses) and we propose a new class of preconditioners based on an FMM-based spatial decomposition of the double-layer operator. We present an experimental study in which we compare the different approaches and we discuss the merits and shortcomings of our approach. Our method can be easily extended to other second-kind integral equations with non-oscillatory kernels in two and three dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fast integral equation methods for the modified Helmholtz equation

We present a collection of integral equation methods for the solution to the two-dimensional, modified Helmholtz equation, $u(\x) - α^2 Δu(\x) = 0$, in bounded or unbounded multiply-connected domains. We consider both Dirichlet and Neumann problems. We derive well-conditioned Fredholm integral equations of the second kind, which are discretized using high-order, hybrid Gauss-trapezoid rules. Our fast multipole-based iterative solution procedure requires only O(N) or $O(N\log N)$ operations, where N is the number of nodes in the discretization of the boundary. We demonstrate the performance of the methods on several numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fast integral equation methods for the heat equation and the modified Helmholtz equation in two dimensions

We present an efficient integral equation approach to solve the heat equation, $u_t (\x) - Δu(\x) = F(\x,t)$, in a two-dimensional, multiply connected domain, and with Dirichlet boundary conditions. Instead of using integral equations based on the heat kernel, we take the approach of discretizing in time, first. This leads to a non-homogeneous modified Helmholtz equation that is solved at each time step. The solution to this equation is formulated as a volume potential plus a double layer potential.The volume potential is evaluated using a fast multipole-accelerated solver. The boundary conditions are then satisfied by solving an integral equation for the homogeneous modified Helmholtz equation. The integral equation solver is also accelerated by the fast multipole method (FMM). For a total of $N$ points in the discretization of the boundary and the domain, the total computational cost per time step is $O(N)$ or $O(N\log N)$.

0 citations0 reviewsNo code linkedOpen access