MATLAB code for sedimenting problem using Lattice Boltzmann approach

I am looking for a programmer who is expert in CFD and MATLAB programming. Problem Setup: We have a rectangular fluid domain confined by two vertically placed rigid walls. An isolated particle(may be rigid or flexible) is located along the channel centerline. The assumptions are fluid is viscous, incompressible and newtonian in nature. we want to study the settling velocity and deformation experienced (in case of flexible particles) while sedimenting under the action of gravity. the fluid domain will be assumed to be doubly periodic. we want the MATLAB code for the same problem setup using Lattice Boltzmann approach. The data provide for our side will be: !) shape of the particle 2) Unstructured mesh of the particle with all the nodes and edge details (made with help of distmesh software) 3)Length of the fluid domain, position of the walls, viscosity of the fluid and density difference between the structure and the fluid will also be provided.