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 trajectories of hydrodynamically interacting spheres in simple and oscillatory shear flows of a Newtonian fluid. Movies of particle 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 main programs (Simple_Shear_Flow.m and Oscillatory_Shear_Flow.m) call 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.dat and data2.dat .
Housam Binous (2020). Hydrodynamically Interacting Spheres at low Reynolds numbers in shear flows of a Newtonian fluid (https://www.mathworks.com/matlabcentral/fileexchange/8704-hydrodynamically-interacting-spheres-at-low-reynolds-numbers-in-shear-flows-of-a-newtonian-fluid), MATLAB Central File Exchange. Retrieved .