COPSS Features

General-geometry Ewald-like Method

The scheme provided by the General geometry Ewald-like Method (GgEM) is a platform for multi-scale problems in physics and chemistry. Applied to hydrodynamic equations, the method solves for the proper behavior (both Brownian and Hydrodynamic) of discrete entities such as polymers, cells and colloids, in continuous fluid media. The Navier-Stokes equations are solved using a Green's function approach for long-range interactions. At right is a demonstration of GGEM methods applied to the resolution of flow fields in suspensions of self-motile particles.

Fast Multipole Boundary-element Method

The Fast Multipole Boundary Element Method (FMBEM) is a code framework for various integral problems encountered in scientific and engineering communities. It runs in parallel and utilizes iterative solver GMRES accelerated by a kernel-independent Fast Multipole Method (FMM) with O(N) computational complexity. It applies to electrostatic polarization problem involving arbitrary-shaped polarizable dielectric objects embedded in a dielectric continuum, and also to the calculation of interactions between ferromagnetic/ferroelectric objects with arbitrary shapes. The code is being developed using libMesh (finite-element library) and ScalFMM (FMM library). At right is a demonstration of the FMBEM applied to the electrostatic polarization problem among interacting polarizable particles.

Parallel Finite-Element General-Geometry Ewald-like Method (pFE-GgEm)

Parallel finite element general geometry Ewald-like method (pFE-GgEm) is an efficient O(N) and fully scalable computation framework to model the dynamics of Brownian particles considering long-range hydrodynamic interactions (HI) in any geometries. The framework combines a FEM Stokes’ solver with GgEm, a mid-point integration scheme, and a Chebyshev polynomial approximation based on the fluctuation-dissipation theorem. It is built on open-source scientific packages, such as LIBMESH (framework of solving PDEs), PETSC (parallel linear, non-linear equation solver) and SLEPC (eigenvalue calculations). At right is a demonstration of the pFE-GgEm applied to study the behaviors of DNA chains under an elongational flow within a cross-channel geometry.