|
Stokesian Dynamics, a method developed by Brady and Bossis in the 80s, simulate the 3D motion of hydrodynamically interacting spheres at low Reynolds numbers. The present contribution computes the 3D trajectories of hydrodynamically interacting spheres in simple and oscillatory shear flows of a Newtonian fluid. Animations of particles' trajectories are obtained as a result of the simulations. The FTS program was originally written by Professor Ron Phillips during his PhD thesis at MIT in 1989. The code was adapted from Fortran 77 to Mathematica by Housam Binous. Using Euler's integration scheme was one simplification, of the original code, that was performed by Housam Binous. This code is posted with permission from Professor Ron Phillips. The notebook Main_Stokesian_Dynamics_FTS_TWO_SPHERES.nb calls a .m file called VEL5.m, which needs four files TEMPM.m, FTS.m, FTSMOB.m and SSI.m as well as two data file called data1.m and data2.m. All files should be placed in the C: drive in the folder TwoSpheresShearFlow. If not, some changes to the Get instructions (<<) should be done by users.
|
|