Speaker: 

Professor Robert Krasny

Institution: 

University of Michigan

Time: 

Monday, June 1, 2009 - 4:00pm

Location: 

RH 306

We describe some recent developments involving particle methods in
electrostatics and fluid dynamics.

First we describe a Cartesian treecode for charged particle systems
undergoing screened Coulomb interactions in three dimensions. The
treecode reduces the operation count for computing potentials and
forces from $O(N^2)$ to $O(N\log N)$, where $N$ is the number of
particles in the system. The algorithm is especially well suited for
particles distributed on a surface. This is a step towards a treecode-
accelerated boundary integral solver for implicit solvent models in
biomolecular dynamics. (joint work with Peijun Li, Hans Johnston, and
Weihua Geng)

Second we describe a Lagrangian panel method for computing vortex
sheet motion in three dimensional fluid flow. Here too a treecode is
used, to evaluate the regularized Biot-Savart integral for the
sheet's self-induced velocity. The method is applied to compute the
azimuthal instability of a vortex ring. Details of the core dynamics
are clarified by tracking material lines on the sheet surface.
Results show the collapse of the vortex core in each wavelength and
the radial ejection of ringlets. These events are correlated with
local axial flow in the core of the ring. (joint work with Hualong
Feng and Leon Kaganovskiy)