Turbulence modeling for flows around convex features giving rapid eddy distortion

Reynolds averaged Navier-Stokes model performances in the stagnation and wake regions for turbulent flows with relatively large Lagrangian length scales (generally larger than the scale of geometrical features) approaching small cylinders (both square and circular) is explored. The effective cylinder (or wire) diameter based Reynolds number, ReW ≤ 2.5 × 103. The following turbulence models are considered: a mixing-length; standard Spalart and Allmaras (SA) and streamline curvature (and rotation) corrected SA (SARC); Secundov's νt-92; Secundov et al.'s two equation νt-L; Wolfshtein's k-l model; the Explicit Algebraic Stress Model (EASM) of Abid et al.; the cubic model of Craft et al.; various linear k-ε models including those with wall distance based damping functions; Menter SST, k-ω and Spalding's LVEL model. The use of differential equation distance functions (Poisson and Hamilton-Jacobi equation based) for palliative turbulence modeling purposes is explored. The performance of SA with these distance functions is also considered in the sharp convex geometry region of an airfoil trailing edge. For the cylinder, with ReW ≈ 2.5 × 103 the mixing length and k-l models give strong turbulence production in the wake region. However, in agreement with eddy viscosity estimates, the LVEL and Secundov νt-92 models show relatively little cylinder influence on turbulence. On the other hand, two equation models (as does the one equation SA) suggest the cylinder gives a strong turbulence deficit in the wake region. Also, for SA, an order or magnitude cylinder diameter decrease from ReW = 2500 to 250 surprisingly strengthens the cylinder's disruptive influence. Importantly, results for ReW ≪ 250 are virtually identical to those for ReW = 250 i.e. no matter how small the cylinder/wire its influence does not, as it should, vanish. Similar tests for the Launder-Sharma k-ε, Menter SST and k-ω show, in accordance with physical reality, the cylinder's influence diminishing albeit slowly with size. Results suggest distance functions palliate the SA model's erroneous trait and improve its predictive performance in wire wake regions. Also, results suggest that, along the stagnation line, such functions improve the SA, mixing length, k-l and LVEL results. For the airfoil, with SA, the larger Poisson distance function increases the wake region turbulence levels by just under 5%. © 2007 Elsevier Inc. All rights reserved.

Thermal transport in high porosity cellular metal foams

The heat dissipation capability of highly porous cellular metal foams with open cells subject to forced air convection is studied using a combined experimental and analytical approach. The cellular morphologies of six FeCrAlY (an iron-based alloy) foams and six copper alloy foams with a range of pore sizes and porosities are quantified with the scanning electronic microscope and image analysis. Experimental measurements on pressure drop and heat transfer for copper foams are carried out. A numerical model for forced convection across open-celled metal foams is subsequently developed, and the predictions are compared with those measured. Reasonably good agreement with test data is obtained, given the complexity of the cellular foam morphology and the associated momentum/energy transport. The results show that cell size has a more significant effect on the overall heat transfer than porosity. An optimal porosity is obtained based on the balance between pressure drop and overall heat transfer, which decreases as the Reynolds number is increased.

Large eddy simulations for fan-OGV broadband noise prediction

The work in this paper forms part of a project on the use of large eddy simulation (LES) for broadband rotor-stator interaction noise prediction. Here we focus on LES of the flow field near a fan blade trailing edge. The first part of the paper aims to evaluate LES suitability for predicting the near-field velocity field for a blunt NACA-0012 airfoil at moderate Reynolds numbers (2× 10 5 and 4× 10 5). Preliminary computations of turbulent mean and root-mean-square velocities, as well as energy spectra at the trailing edge, are compared with those from a recent experiment.1 The second part of the paper describes preliminary progress on an LES calculation of the fan wakes on a fan rig. 2 The CFD code uses a mixed element unstructured mesh with a median dual control volume. A wall-adapting local eddy-viscosity sub-grid scale model is employed. A very small amount of numerical dissipation is added in the numerical scheme to keep the compressible solver stable. Further results for the fan turbulentmean and RMS velocity, and especially the aeroacoustics field will be presented at a later stage. Copyright © 2008 by Qinling LI, Nigel Peake & Mark Savill.

Forcing of self-excited round jet diffusion flames

In this experimental and numerical study, two types of round jet are examined under acoustic forcing. The first is a non-reacting low density jet (density ratio 0.14). The second is a buoyant jet diffusion flame at a Reynolds number of 1100 (density ratio of unburnt fluids 0.5). Both jets have regions of strong absolute instability at their base and this causes them to exhibit strong self-excited bulging oscillations at welldefined natural frequencies. This study particularly focuses on the heat release of the jet diffusion flame, which oscillates at the same natural frequency as the bulging mode, due to the absolutely unstable shear layer just outside the flame. The jets are forced at several amplitudes around their natural frequencies. In the non-reacting jet, the frequency of the bulging oscillation locks into the forcing frequency relatively easily. In the jet diffusion flame, however, very large forcing amplitudes are required to make the heat release lock into the forcing frequency. Even at these high forcing amplitudes, the natural mode takes over again from the forced mode in the downstream region of the flow, where the perturbation is beginning to saturate non-linearly and where the heat release is high. This raises the possibility that, in a flame with large regions of absolute instability, the strong natural mode could saturate before the forced mode, weakening the coupling between heat release and incident pressure perturbations, hence weakening the feedback loop that causes combustion instability. © 2009 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Convective heat dissipation with lattice-frame materials

This paper presents experimental results on heat transfer and pressure drop for a compact heat sink made of fully triangulated, lightweight (porosity∼0.938), aluminum lattice-frame materials (LFMs). Due to the inherent structural anisotropy of the LFMs, two mutually perpendicular orientations were selected for the measurements. Constant heat flux was applied to the heat sink under steady state conditions, and dissipated by forced air convection. The experimental data were compared with those predicted from an analytical model based on fin analogy. The experimental results revealed that pressure drop is strongly dependent upon the orientation of the structure, due mainly to the flow blockage effect. For heat transfer measurements, typical local temperature distributions on the substrate under constant heat flux conditions were captured with infrared camera. The thermal behavior of LFMs was found to follow closely that of cylinder banks, with early transition Reynolds number (based on strut diameter) equal to about 300. The Nusselt number prediction from the fin-analogy correlates well with experimental measurements, except at low Reynolds numbers where a slightly underestimation is observed. Comparisons with empty channels and commonly used heat exchanger media show that the present LFM heat sink can remove heat approximately seven times more efficient than an empty channel and as efficient as a bank of cylinders at the same porosity level. The aluminum LFMs are extremely stiff and strong, making them ideal candidates for multifunctional structures requiring both heat dissipation and mechanical load carrying capabilities. © 2003 Elsevier Ltd. All rights reserved.

Curvature and wrinkling of premixed flame Kernels comparison between planar laser induced fluorescence (PLIF) and 3D direct numerical simulations (DNS)

A direct comparison between time resolved PLIF measurements of OH and two dimensional slices from a full three dimensional DNS data set of turbulent premixed flame kernels in lean methane/air mixture was presented. The local flame structure and the degree of flame wrinkling were examined in response to differing turbulence intensities and turbulent Reynolds numbers. Simulations were performed using the SEGA DNS code, which is based on the solution of the compressible Navier Stokes, species, and energy equations for a lean hydrocarbon mixture. For the OH PLIF measurements, a cluster of four Nd:YAG laser was fired sequentially at high repetition rates and used to pump a dye laser. The frequency doubled laser beam was formed into a sheet of 40 mm height using a cylindrical telescope. The combination of PLIF and DNS has been demonstrated as a powerful tool for flame analysis. This research will form the basis for the development of sub-grid-scale (SGS) models for LES of lean-premixed combustion systems such as gas turbines. This is an abstract of a paper presented at the 30th International Symposium on Combustion (Chicago, IL 7/25-30/2004).

Normal shock boundary layer control with various vortex generator geometries

Various vortex generators which include ramp, split-ramp and a new hybrid concept "ramped-vane" are investigated under normal shock conditions with a diffuser at Mach number of 1.3. The dimensions of the computational domain were designed using Reynolds Average Navier-Stokes studies to be representative of the flow in an external-compression supersonic inlet. Using this flow geometry, various vortex generator concepts were studied with Implicit Large Eddy Simulation. In general, the ramped-vane provided increased vorticity compared to the other devices and reduced the separation length downstream of the device centerline. In addition, the size, edge gap and streamwise position respect to the shock were studied for the ramped-vane and it was found that a height of about half the boundary thickness and a large trailing edge gap yielded a fully attached flow downstream of the device. This ramped-vane also provided the largest reduction in the turbulent kinetic energy and pressure fluctuations. Additional benefits include negligible drag while the reductions in boundary layer displacement thickness and shape factor were seen compared to other devices. © 2010 by Sang Lee.

Normal shock boundary layer control with various vortex generator geometries

Various vortex generators which include ramp, split-ramp and a new hybrid concept "ramped-vane" are investigated under normal shock conditions with a diffuser at Mach number of 1.3. The dimensions of the computational domain were designed using Reynolds Average Navier-Stokes studies to be representative of the flow in an external-compression supersonic inlet. Using this flow geometry, various vortex generator concepts were studied with Implicit Large Eddy Simulation. In general, the ramped-vane provided increased vorticity compared to the other devices and reduced the separation length downstream of the device centerline. In addition, the size, edge gap and streamwise position respect to the shock were studied for the ramped-vane and it was found that a height of about half the boundary thickness and a large trailing edge gap yielded a fully attached flow downstream of the device. This ramped-vane also provided the largest reduction in the turbulent kinetic energy and pressure fluctuations. Additional benefits include negligible drag while the reductions in boundary layer displacement thickness and shape factor were seen compared to other devices. © 2011 Elsevier Ltd.