Understanding animal locomotion requires modeling the interaction of the organism with its environment. Locomotion within granular media like sand, soil, and debris that display both solid and fluid-like behavior in response to stress is less studied than locomotion within fluids or on solid ground. To begin to reveal the secrets of movement in sand, I developed models to explain the subsurface locomotion of the sand-swimming sandfish lizard. I developed a resistive force theory (RFT) with empirical force laws to explain the swimming speed observed in animal experiments. By varying the amplitude of the undulation in the RFT, I found that the range of amplitude used by the animal predicted the optimal swimming speed. I developed a numerical model of the sandfish coupled to a discrete element method simulation of the granular medium to test assumptions in the RFT and to study more detailed mechanics of sand-swimming. Inspired by the shovel-shaped head of the sandfish lizard, I used the simulation to study lift forces in granular media: I found that when a submerged intruder moved at a constant speed within a granular medium it experienced a lift force whose sign and magnitude depended on the intruder shape. The principles learned from the models guided the development of a biologically inspired robot that swam within granular media with similar performance to the lizard.

At low Reynolds number flows, the effect of inertia becomes negligible and the fluid motion is dominated by the effect of viscous forces. Understanding of the behaviour of low Reynolds number flows underpins the prediction of the motion of microorganisms and particle sedimentation as well as the development of micro-robots that could potentially swim inside the human body to perform targeted drug/cell delivery and non-invasive microsurgery. The work in this thesis focuses on developing an understanding in the mathematical analysis of objects moving at low Reynolds numbers. A boundary element implementation of the Method of regularized Stokeslets (MRS) is applied to analyse the low Reynolds number flow field around an object of simple shape (sphere and cube). It also showed that the results obtained by a boundary element implementation for an unbounded cube, where singularities are presented in the corners of the cube, agrees with more complex solutions methods such as a GBEM and FEM.A methodology for analysing the effect of walls by locating collocation points on the surface of the walls and the object is presented. First at all, this methodology is validated with a boundary element implementation of the method of images for a sphere at different locations. Then, the method is extended when more than one wall is presented. This methodology is applied to predict the velocity filed of a cube moving in a tow tank at low Reynolds numbers for two different cases with a supporting rod similar to an experimental set-up, and without the supporting rod as in the CFD simulations based on the FVM. The results indicate a good match between CFD and the MRS, and an excellent approximation between the MRS and experimental data from PIV measurements. The drag, thrust and torque generated by helices moving at low Reynolds numbers in an unbounded medium is analysed by the resistive force theory, a slender body theory, and a boundary element method of the MRS. The results show that the resistive force theory predict accurately the drag, thrust and torque of moving helices when the resistive force coefficients are calculated from a slender body theory approximation by calculating independently the resistive force coefficients for translation and rotation, because it is observed that the resistive force coefficients depend also of the nature of motion. Moreover, the thrust generated by helices of different pitch angles is analysed calculated by a CFD numerical simulation based on the FVM and a boundary element implementation, an compared with experimental data. The results also show an excellent prediction between the boundary element implementation, the CFD results and the experimental data. Finally, a boundary element implementation of the MRS is applied to predict swimming of a biomimetic swimmer that mimics the motion of E.coli bacteria in an unbounded medium. The results are compared with the propulsive velocity and induced angular velocity measurement by recording the motion of the biomimetic swimmer in a square tank. It is observed that special care needs to be taken when the biomimetic swimmer is modelled inside the tank, as there is an apparent increment in the calculate thrust propulsion which does not represent a real situation of the biometic swimmer which propels by a power supply. However, this increment does not represent the condition of the biomimetic swimmer and a suggested methodology based on the solution from an unbounded case and when the swimmer is moving inside the tank is presented. In addition, the prediction of the free-swimming velocity for the biomimetic swimmer agrees with the results obtained by the MRS when the resistive force coefficients are calculated from a SBT implementation. The results obtained in this work have showed that a boundary element implementation of the MRS produces results comparable with more complex numerical implementations such as GBEM, FEM, FVM, and also an excellent agreement with results obtained from experimentation. Therefore, it is a suitable and easy to apply methodology to analyse the motion of swimmers at low Reynolds numbers, such as the biomimetic swimmer modelled in this work.