The finite-difference PHYSALIS method can be used to perform direct numerical simulations of fluid-partilces interactions in complex geometry. Noteworthy features of the method are:
- The no-slip condition at the particle surface is satisfied exactly.
- The lower-order underdetermined coefficients give the force, torque, and stresslet on the sphere directly, thus avoiding the need to integrate the calculated fluid stress over the particle surface.
- If the Stokes approximation is used beyond its approximate domain of validity, the iterative process will not converge. In this sense, the method contains a built-in accuracy test.
The current solver has been implemented using a locally refined mesh for the spheres in the PHYSALIS method.
Download the tar file
and uncompress it.
First run the command make main
. Then, run the command ./main
The initial conditions for three velocity components can be set in the file initCond.cpp
. The boundary conditions for three velocity components can be set in the file boundaryCond.cpp
. They may be Dirichlet, Neumann, periodic, or other boundary conditions.
- Qianlong Liu <Qianlong.Liu@jhu.edu>
Advisor: Prof. Prosperetti email@example.com
- April 2008 - March 2009