A non-boundary-conforming formulation for simulating transitional and turbulent flows with complex geometries and dynamically moving boundaries on fixed orthogonal grids is developed. The underlying finite-difference solver for the filtered incompressible Navier-Stokes equations in both Cartesian and cylindrical coordinates is based on a second-order fractional step method on staggered grid. To satisfy the boundary conditions on an arbitrary immersed interface, the velocity field at the grid points near the interface is reconstructed locally without smearing the sharp interface. The complications caused by the Eulerian grid points emerging from a moving solid body into the fluid phase are treated with a novel "field-extension" strategy. To treat the two-way interactions between the fluid and structure, a strong coupling scheme based on Hamming's fourth-order predictor-corrector method has • attack angle at Re = 10, 000. Then, the turbulent flow over a traveling wavy wall at Re = 10, 170 are simulated are compared with the detailed DNS using bodyfitted grid in the literature. Finally, the simulation of the transitional flow past a prosthetic mechanical heart valve with moving leaflets at Re = 4, 000 has been performed. All results are in good agreement with the available reference data.

No abstract

▪ Abstract  Because the cost of large-eddy simulations (LES) of wall-bounded flows that resolve all the important eddies depends strongly on the Reynolds number, methods to bypass the wall layer are required to perform high-Reynolds-number LES at a reasonable cost. In this paper the available methodologies are reviewed, and their ranges of applicability are highlighted. Various unresolved issues in wall-layer modeling are presented, mostly in the context of engineering applications.

No abstract

SUMMARYIn the present study, a computational investigation was carried out to understand the influence of flexibility on the aerodynamic performance of a hovering wing. A flexible, two-dimensional, two-link model moving within a viscous fluid was considered. The Navier-Stokes equations governing the fluid dynamics were solved together with the equations governing the structural dynamics by using a strongly coupled fluid-structure interaction scheme. Harmonic kinematics was used to prescribe the motions of one of the links, thus effectively reducing the wing to a single degree-of-freedom oscillator. The wing's flexibility was characterized by the ratio of the flapping frequency to the natural frequency of the structure. Apart from the rigid case, different values of this frequency ratio (only in the range of 1/2 to 1/6) were considered at the Reynolds numbers of 75, 250 and 1000. It was found that flexibility can enhance aerodynamic performance and that the best performance is realized when the wing is excited by a non-linear resonance at 1/3 of the natural frequency. Specifically, at Reynolds numbers of 75, 250 and 1000, the aerodynamic performance that is characterized by the ratio of lift coefficient to drag coefficient is respectively increased by 28%, 23% and 21% when compared with the corresponding ratios of a rigid wing driven with the same kinematics. For all Reynolds numbers, the lift generated per unit driving power is also enhanced in a similar manner. The wake capture mechanism is enhanced, due to a stronger flow around the wing at stroke reversal, resulting from a stronger end of stroke vortex at the trailing edge. The present study provides some clues about how flexibility affects the aerodynamic performance in low Reynolds number flapping flight. In addition, it points to the importance of considering non-linear resonances for enhancing aerodynamic performance.

Comparisons of inflow conditions for large-eddy simulations of turbulent, wall-bounded flows are carried out. Consistent with previous investigations, it is found that the spectral content of the inflow velocity is important. Inflow conditions based on random-noise, or small-scale eddies only, dissipate quickly. Temporal and spatial filtering of a time series obtained from a separate calculation indicates that it is important to capture eddies of dimensions equal to or larger than the integral length scale of the flow. Three methods for generating inflow velocity fields are tested in a simulation of spatially developing turbulent channel flow. Synthetic turbulence generation methods that introduce realistic length scales are more suitable than uncorrelated random noise, but still require fairly long development lengths before realistic turbulence is established. A recycling method based on the use of turbulent data obtained from a separate calculation, in different flow conditions, was found to result in more rapid transition. A forcing method that includes a control loop also appears to be effective by generating turbulence with the correct Reynolds stresses and correlations within less than ten channel half heights.

No abstract

This work focuses on the direct numerical simulation of the pulsatile now through a bileaflet mechanical heart valve under physiological conditions and in a realistic aortic root geometry. The motion of the valve leaflets has been computed from the forces exerted by the fluid on the structure both being considered as a single dynamical system. To this purpose the immersed boundary method, combined with a fluid structure interaction algorithm, has shown to be an inexpensive and accurate technique for such complex flows. Several complete flow cycles have been simulated in order to collect enough phase-averaged statistics, and the results are in good agreement with experimental data obtained for a similar configuration. The flow analysis, strongly relying on the data accessibility provided by the numerical simulation, shows how some features of the leaflets motion depend on the flow dynamics and that the criteria for the red cell damages caused by the valve need to be formulated using very detailed analysis. In particular, it is shown that the standard Eulerian Computation of the Reynolds stresses, usually employed to assess the risk of haemolysis, might not be adequate oil several counts: (i) Reynolds stresses are only one part of the solicitation, the other part being the viscous stresses, (ii) the characteristic scales or the two solicitations are very different and the Reynolds stresses act on lengths much larger than the red cells diameter and (iii) the Eulerian zonal assessment of the stresses completely misses the information of time exposure to the solicitation which is a fundamental ingredient for the phenomenon of haemolysis. Accordingly, the trajectories of several fluid particles have been tracked in a Lagrangian way and the pointwise instantaneous Viscous stress tensor has been computed along the paths. The tensor has been then reduced to an equivalent scalar using the von Mises criterion, and the blood damage index has been evaluated following Grigioni et al. (Biomech. Model Mechanobiol., vol. 4, 2005, p. 249)