|
BEAVIS
Boundary Element Analysis
of VIscous Suspensions |
|
|
BEM formulation
|
|
2-D simulations (movies)
|
|
3-D simulations (stills)
|
BEAVIS is based on a completed double layer boundary integral equation
formulation, with fast multipole acceleration.
This formulation is the most effective when one deals with
flows of suspended particles in a container, because the
system of linear equations that results from collocation and
discretization is nicely conditioned, due to the diagonal
dominance afforded by the second kind integral equations.
For a more complete description of the implementation, go
here.
2-D simulations
The simulations shown below were performed on the IBM SP2 system at
the Maui High Performance Computing Center. Smooth mpeg animations
may be downloaded for each case (right button).
|
Wide-gap Couette geometry
neutrally buoyant particles
30% area fraction
download mpeg |
|
|
Wide-gap Couette geometry
sedimenting particles
10% area fraction
download mpeg |
|
|
Detail of S portion of wide-gap Couette geometry
sedimenting particles
10% area fraction |
|
|
Driven cavity
neutrally buoyant particles
20% area fraction |
|
|
Detail of SE corner in driven cavity
neutrally buoyant particles
20% area fraction |
|
3-D simulations
The 3-D code is still under development. Currently, it is being
benchmarked with known analytical solutions. The figure shown
below is a mesh for two sedimenting particles, where the line
connecting their centers is at 45 degrees to the vertical.
