Wolfram Library Archive

All Collections Articles Books Conference Proceedings
Courseware Demos MathSource Technical Notes
Title Downloads

Three dimensional trajectories of hydrodynamically interacting spheres at low Reynolds numbers in shear flows of Newtonian fluids

Housam Binous
Organization: KFUPM
Department: Chemical Engineering
URL: http://sites.google.com/site/homepageofdrhousambinous/
Revision date


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.

*Engineering > Chemical Engineering
*Engineering > Mechanical and Structural Engineering
*Science > Physics > Fluid Mechanics

TwoSpheresShearFlow.zip (133.6 KB) - ZIP archive