The sharp-interface CURVIB approach of Ge and Sotiropoulos [L. Ge, F. Sotiropoulos, A Numerical Method for Solving the 3D Unsteady Incompressible Navier-Stokes Equations in Curvilinear Domains with Complex Immersed Boundaries, Journal of Computational Physics 225 (2007) 1782-1809] is extended to simulate fluid structure interaction (FSI) problems involving complex 3D rigid bodies undergoing large structural displacements. The FSI solver adopts the partitioned FSI solution approach and both loose and strong coupling strategies are implemented. The interfaces between immersed bodies and the fluid are discretized with a Lagrangian grid and tracked with an explicit front-tracking approach. An efficient ray-tracing algorithm is developed to quickly identify the relationship between the background grid and the moving bodies. Numerical experiments are carried out for two FSI problems: vortex induced vibration of elastically mounted cylinders and flow through a bileaflet mechanical heart valve at physiologic conditions. For both cases the computed results are in excellent agreement with benchmark simulations and experimental measurements. The numerical experiments suggest that both the properties of the structure (mass, geometry) and the local flow conditions can play an important role in determining the stability of the FSI algorithm. Under certain conditions unconditionally unstable iteration schemes result even when strong coupling FSI is employed. For such cases, however, combining the strong-coupling iteration with under-relaxation in conjunction with the Aitken's acceleration technique is shown to effectively resolve the stability problems. A theoretical analysis is presented to explain the findings of the numerical experiments. It is shown that the ratio of the added mass to the mass of the structure as well as the sign of the local time rate of change of the force or moment imparted on the structure by the fluid determine the stability and convergence of the FSI algorithm. The stabilizing role of under-relaxation is also clarified and an upper bound of the required for stability under-relaxation coefficient is derived.

SUMMARYFor all Re, however, the swimming power is shown to be significantly greater than that required to tow the rigid body at the same speed. We also show that the variation of the total drag and its viscous and form components with St depend on the Re. For Re=300, body undulations increase the drag over the rigid body level, while significant drag reduction is observed for Re=4000. This difference is shown to be due to the fact that at sufficiently high Re the drag force variation with St is dominated by its form component variation, which is reduced by undulatory swimming for St>0.2. Finally, our simulations clarify the 3D structure of various wake patterns observed in experiments -single and double row vortices -and suggest that the wake structure depends primarily on the St. Our numerical findings help elucidate the results of previous experiments with live fish, underscore the importance of scale (Re) effects on the hydrodynamic performance of carangiform swimming, and help explain why in nature this mode of swimming is typically preferred by fast swimmers.

SUMMARYWe carry out fluid-structure interaction simulations of self-propelled virtual swimmers to investigate the effects of body shape (form) and kinematics on the hydrodynamics of undulatory swimming. To separate the effects of form and kinematics, we employ four different virtual swimmers: a carangiform swimmer (i.e. a mackerel swimming like mackerel do in nature); an anguilliform swimmer (i.e. a lamprey swimming like lampreys do in nature); a hybrid swimmer with anguilliform kinematics but carangiform body shape (a mackerel swimming like a lamprey); and another hybrid swimmer with carangiform kinematics but anguilliform body shape (a lamprey swimming like a mackerel). By comparing the performance of swimmers with different kinematics but similar body shapes we study the effects of kinematics whereas by comparing swimmers with similar kinematics but different body shapes we study the effects of form. We show that the anguilliform kinematics not only reaches higher velocities but is also more efficient in the viscous (Re~10 2 ) and transitional (Re~10 3 ) regimes. However, in the inertial regime (Reϱ) carangiform kinematics achieves higher velocities and is also more efficient than the anguilliform kinematics. The mackerel body achieves higher swimming speeds in all cases but is more efficient in the inertial regime only whereas the lamprey body is more efficient in the transitional regime. We also show that form and kinematics have little overall effect on the 3-D structure of the wake (i.e. single vs double row vortex streets), which mainly depends on the Strouhal number. Nevertheless, body shape is found to somewhat affect the small-scale features and complexity of the vortex rings shed by the various swimmers.