Fady M. Najjar

Wake topology and hydrodynamic performance of low-aspect-ratio flapping foils

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.

Modeling superheat removal during continuous casting of steel slabs

To investigate superheat dissipation in a continuous slab casting machine, mathematical models have been developed to compute fluid flow velocities, temperature distribution within the liquid pool, heat transfer to the inside of the solidifying shell, and its effect on growth of the shell. Three-dimensional (3-D) velocity and heat-transfer predictions compare reasonably with pre-vious experimental measurements and two-dimensional (2-D) calculations. The results indicate that the maximum heat input to the shell occurs near the impingement point on the narrow face and confirm that most of the superheat is dissipated in or just below the mold. Superheat tem-perature and casting speed have the most important and direct influence on heat flux. The effects of other variables, including mold width, nozzle jet angle, and submergence depth, are also investigated. Calculated heat flux profiles are then input to a one-dimensional (1-D) solidifi-cation model to calculate growth of the shell. Shell thickness profiles down the wide and narrow faces are compared with the predictions of conventional heat conduction models and available measurements.

Simulation of fluid flow inside a continuous slab-casting machine

A finite element model has been developed and applied to compute the fluid flow distribution inside the shell in the mold region of a continuous, steel slab-casting machine. The model was produced with the commercial program FIDAP, which allows this nonlinear, highly turbulent problem to be simulated using the K- ε turbulence model. It consists of separate two-dimensional (2-D) models of the nozzle and a section through the mold, facing the broad face. The predicted flow patterns and velocity fields show reasonable agreement with experimental observations and measurements conducted using a transparent plastic water model. The effects of nozzle angle, casting speed, mold width, and turbulence simulation parameters on the flow pattern have been investigated. The overall flow field is relatively insensitive to process parameters.

STUDY OF DISCRETE TEST FILTERS AND FINITE DIFFERENCE APPROXIMATIONS FOR THE DYNAMIC SUBGRID-SCALE STRESS MODEL

This paper investigates the effects of discrete test filters and finite‐difference approximations for large‐eddy simulations using the dynamic subgrid‐scale stress model. Discrete explicit test filters based on finite‐difference formulations have been constructed and the characteristics of their transfer function are studied. Several definitions of the scaling factor are investigated in the context of the discrete test filters. Two test filters, one based on a discrete representation of the top‐hat filter (A), and another based on a high‐order filtering operation (C) are evaluated in simulations of the turbulent channel flow at Reτ=180. It is found that filter A calculates a higher turbulent viscosity than filter C, which behaves more like a cutoff filter. For the same test filtering operation, the results are found to be sensitive to the ratio of the characteristic lengths of the test and grid filters. By testing two approximations to the convection terms based on second‐order central difference and a no...

Low-frequency unsteadiness in the wake of a normal flat plate

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 Karman vortices are less coherent. In the high-drag regime three-dimensionality is characterized by coherent Karman vortices and reasonably well-organized streamwise vortices connecting the Karman 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).

Heat transfer enhancement mechanisms in inline and staggered parallel-plate fin heat exchangers

Abstract Heat transfer enhancement mechanisms and performance of parallel-plate fin heat exchangers are studied. Geometry effects such as finite fin thickness and inline vs staggered arrangements have been investigated. The time-dependent flow behavior due to vortex shedding has been taken into consideration by solving unsteady equations in two-dimensions. Companion steady symmetrized flow calculations have also been performed to identify the effect of vortices on heat transfer and frictional loss. Additional comparisons have been made with theoretical results for fully developed flow between parallel plates and numerical :results for restarted boundary layers to quantify boundary layer restart mechanism.

Numerical study of pulsatile flow in a constricted channel

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 shear-layer 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.

Properties of the mean recirculation region in the wakes of two-dimensional bluff bodies

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 140000 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 Karman 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.

Numerical study of steady turbulent flow through bifurcated nozzles in continuous casting

Bifurcated nozzles are used in continuous casting of molten steel, where they influence the quality of the cast steel slabs. The present study performs two-dimensional (2-D) and three-dimensional (3-D) simulations of steady turbulent(K- ε) flow in bifurcated nozzles, using a finite-element (FIDAP) model, which has been verified previously with water model experiments. The effects of nozzle design and casting process operating variables on the jet characteristics exiting the nozzle are investigated. The nozzle design parameters studied include the shape, angle, height, width, and thickness of the ports and the bottom geometry. The process operating practices include inlet velocity profile and angle as well as port curvature caused by erosion or inclusion buildup. Results show that the jet angle is controlled mainly by the port angle but is steeper with larger port area and thinner walls. The degree of swirl is increased by larger or rounder ports. The effective port area, where there is no recirculation, is increased by smaller or curved ports. Flow asymmetry is more severe with skewed or angled inlet conditions or unequal port sizes. Turbulence levels in the jet are higher with higher casting speed and smaller ports.

Forces on a finite-sized particle located close to a wall in a linear shear flow

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 wit...