A sharp interface immersed boundary method for simulating incompressible viscous flow past threedimensional immersed bodies is described. The method employs a multi-dimensional ghost-cell methodology to satisfy the boundary conditions on the immersed boundary and the method is designed to handle highly complex three-dimensional, stationary, moving and/or deforming bodies. The complex immersed surfaces are represented by grids consisting of unstructured triangular elements; while the flow is computed on non-uniform Cartesian grids. The paper describes the salient features of the methodology with special emphasis on the immersed boundary treatment for stationary and moving boundaries. Simulations of a number of canonical two-and three-dimensional flows are used to verify the accuracy and fidelity of the solver over a range of Reynolds numbers. Flow past suddenly accelerated bodies are used to validate the solver for moving boundary problems. Finally two cases inspired from biology with highly complex three-dimensional bodies are simulated in order to demonstrate the versatility of the method.

Numerical simulations are used to investigate the effect of aspect ratio on the wake topology and hydrodynamic performance of thin ellipsoidal flapping foils. The study is motivated by the quest to understand the hydrodynamics of fish pectoral fins. The simulations employ an immersed boundary method that allows us to simulate flows with complex moving boundaries on fixed Cartesian grids. A detailed analysis of the vortex topology shows that the wake of low-aspect-ratio flapping foils is dominated by two sets of interconnected vortex loops that evolve into distinct vortex rings as they convect downstream. The flow downstream of these flapping foils is characterized by two oblique jets and the implications of this characteristic on the hydrodynamic performance are examined. Simulations are also used to examine the thrust and propulsive efficiency of these foils over a range of Strouhal and Reynolds numbers as well as pitch-bias angles.

To understand and better model the hydrodynamic force acting on a finite-sized particle moving in a wall-bounded linear shear flow, here we consider the two limiting cases of (a) a rigid stationary spherical particle in a linear wall-bounded shear flow and (b) a rigid spherical particle in rectilinear motion parallel to a wall in a quiescent ambient flow. In the present computations, the particle Reynolds number ranges from 2 to 250 at separation distances to the wall from nearly sitting on the wall to far away from the wall. First we characterize the structure of the wake for a stationary particle in a linear shear flow and compare with those for a particle moving parallel to a wall in a quiescent ambient [see L. Zeng, S. Balachandar, and P. Fischer, J. Fluid Mech. 536, 1 (2005)]. For both these cases we present drag and lift results and obtain composite drag and lift correlations that are valid for a wide range of Re and distance from the wall. These correlations have been developed to be consistent with all available low Reynolds number theories and approach the appropriate uniform flow results at large distance from the wall. Particular attention is paid to the case of particle in contact with the wall and the computational results are compared with those from experiments.

Pulsatile flow in a planar channel with a one-sided semicircular constriction has been simulated using direct numerical simulation and large-eddy simulation. This configuration is intended as a simple model for studying blood flow in a constricted artery. Simulations have been carried out over a range of Reynolds numbers (based on channel height and peak bulk velocity) from 750 to 2000 and a fixed non-dimensional pulsation frequency of 0.024. The results indicate that despite the simplicity of the chosen geometry, the simulated flow exhibits a number of features that have been observed in previous experiments carried out in more realistic configurations. It is found that over the entire Reynolds number range studied here, the flow downstream of the constriction is dominated by the complex dynamics associated with two shear-layers, one of which separates from the lip of the constriction and other from the opposite wall. Computed statistics indicate that for Reynolds numbers higher than about 1000, the flow transitions to turbulence downstream of the region where the separated shear layers first reattach to the channel walls. Large fluctuations in wall pressure and shear stress have also been associated with this reattachment phenomenon. Frequency spectra corresponding to velocity and pressure fluctuations have been analysed in detail and these indicate the presence of a characteristic shearlayer frequency which increases monotonically with Reynolds number. For Reynolds numbers greater than 1000, this frequency is found to be associated with the periodic formation of vortex structures in the shear-layers and the impact of this characteristic shear-layer frequency on the dynamics of the flow is described in detail.

No abstract

No abstract

The separated flow past a zero-thickness flat plate held normal to a free stream at Re=250 has been investigated through numerical experiments. The long-time signatures of the drag and lift coefficients clearly capture a low-frequency unsteadiness with a period of approximately 10 times the primary shedding period. The amplitude and frequency of drag and lift variations during the shedding process are strongly modulated by the low frequency. A physical interpretation of the low-frequency behaviour is that the flow gradually varies between two different regimes: a regime H of high mean drag and a regime L of low mean drag. It is observed that in regime H the shear layer rolls up closer to the plate to form coherent spanwise vortices, while in regime L the shear layer extends farther downstream and the rolled-up Kármán vortices are less coherent. In the high-drag regime three-dimensionality is characterized by coherent Kármán vortices and reasonably well-organized streamwise vortices connecting the Kármán vortices. With a non-dimensional spanwise wavelength of about 1.2, the three-dimensionality in this regime is reminiscent of mode-B three-dimensionality. It is observed that the high degree of spanwise coherence that exists in regime H breaks down in regime L. Based on detailed numerical flow visualization we conjecture that the formation of streamwise and spanwise vortices is not in perfect synchronization and that the low-frequency unsteadiness is the result of this imbalance (or phase mismatch).

The properties of the time-and span-averaged mean wake recirculation region are investigated in separated flows over several different two-dimensional bluff bodies. Ten different cases are considered and they divide into two groups : cylindrical geometries of circular, elliptic and square cross-sections and the normal plate. A wide Reynolds number range from 250 to 140 000 is considered, but in all the cases the attached portion of the boundary layer remains laminar until separation. The lower Reynolds number data are from direct numerical simulations, while the data at the higher Reynolds number are obtained from large-eddy simulation and the experimental work of Cantwell & Coles (1983), Krothapalli (1996, personal communication), Leder (1991 and Lyn et al. (1995). Unlike supersonic and subsonic separations with a splitter plate in the wake, in all the cases considered here there is strong interaction between the shear layers resulting in Ka! rma! n vortex shedding. The impact of this fundamental difference on the distribution of Reynolds stress components and pressure in relation to the mean wake recirculation region (wake bubble) is considered. It is observed that in all cases the contribution from Reynolds normal stress to the force balance of the wake bubble is significant. In fact, in the cylinder geometries this contribution can outweigh the net force from the shear stress, so that the net pressure force tends to push the bubble away from the body. In contrast, in the case of normal plate, owing to the longer wake, the net contribution from shear stress outweighs that from the normal stress. At higher Reynolds numbers, separation of the Reynolds stress components into incoherent contributions provides more insight. The behaviour of the coherent contribution, arising from the dominant vortex shedding, is similar to that at lower Reynolds numbers. The incoherent contribution to Reynolds stress, arising from small-scale activity, is compared with that of a canonical free shear layer. Based on these observations a simple extension of the wake model (Sychev 1982 ; Roshko 1993 a, b) is proposed.