PHYSALIS
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.
- Author:
- Qianlong Liu <Qianlong.Liu@jhu.edu>
Advisor: Prof. Prosperetti prosperetti@jhu.edu
- Date:
- April 2008 - March 2009
- Version:
- 1.0
