Search NASA⌕ Search

Engineering topics

Ali Uzun

Publications and source records attributed to Ali Uzun.

Direct Numerical Simulation of Flow Past a Gaussian Bump at a High Reynolds Number

We present the results from a direct numerical simulation of a spanwise-periodic turbulent flow past a Gaussian bump. The problem setup is designed to investigate the interaction of an incoming turbulent boundary layer with the strong favorable and adverse pressure gradients generated by the Gaussian bump as the flow passes over it at a Reynolds number of 340000 based on the bump height, or 4 million based on the bump length. The statistical results from the present simulation are compared against our earlier results at a Reynolds number of 2 million. An internal layer, which forms beneath the strongly accelerated boundary layer over the windward side of the bump, is found to generate its near-wall turbulence stress peaks in closer proximity of the wall in the higher Reynolds-number case. Furthermore, the logarithmic layer of the higher Reynolds-number boundary layer appears more resistant to changes induced by strong acceleration and surface curvature effects over the same region. Despite a nearly identical flow separation point in the two flows, the detached shear layer grows at a faster rate and subsequently reattaches at an earlier point in the higher Reynolds-number flow. The surface pressure and skin-friction distributions over the attached flow region compare well against the corresponding experimental data for both flows. However, some differences appear in the separated flow region, which are attributed to the three-dimensionality of the experimental model setup that is not included in the simulation owing to the spanwise periodic assumption. Comparisons with the stereoscopic particle image velocimetry measurements on the central plane of the experimental model over the windward side of the bump show reasonable overall agreement in the mean velocity components, but the turbulence stress components do not agree well at some streamwise locations. Comparisons over the leeward side of the bump show that the mean separated shear layer in the simulation is tilted significantly more toward the wall than the experimental shear layer on the central plane. This mismatch in the mean shear layer orientation is due to the experimental model three-dimensionality and tunnel end-wall effects, which are not modeled in the present spanwise-periodic simulation.

Favorable and Adverse Pressure Gradients↗

Direct Numerical Simulation of Smooth-Body Flow Separation Around a Ramp

Spanwise-periodic computation of a turbulent flow past a two-dimensional smooth ramp geometry is performed in the form of a direct numerical simulation. The Reynolds number based on the ramp height is about 147,000. A straight section that precedes the smooth ramp allows the incoming turbulent boundary layer to grow under a weak favorable pressure gradient. The boundary layer introduced at the domain inlet has a momentum-thickness based Reynolds number of 2000. The turbulent boundary layer nearing the ramp first interacts with a relatively stronger favorable pressure gradient, followed by a strong adverse pressure gradient. Consequently, the boundary layer experiences a modest acceleration before decelerating and separating. Analysis of the data over this region hints at the formation of an internal layer beneath the accelerated boundary layer. The analysis also reveals that this internal layer forms the origin of the free shear layer that emerges in the deceleration region and separates. The streamwise extent of the separated region is comparable to the ramp length, while the viscous layer thickness near reattachment is about the same as the ramp height; hence, the boundary layer undergoing separation and subsequent reattachment in the present configuration experiences its thickness being amplified by about tenfold. The reattached flow continues to develop further under a diminishing pressure gradient in the recovery region in a similar fashion to a zero pressure gradient turbulent boundary layer.

Turbulence↗

A Dynamic Nonlinear Subgrid-Scale Model for Large-Eddy Simulation of Complex Turbulent Flows

We present a new dynamic nonlinear subgrid-scale (SGS) model for large-eddy simulations (LES) and apply it to compute a flow involving pressure gradients, surface curvature and separation, for which data from a direct numerical simulation are available for comparison. The model, inspired by the triple model idea of Bardina et al. (“Improved Turbulence Models Based on Large Eddy Simulation of Homogeneous, Incompressible, Turbulent Flows,” Report No. TF-19, Thermosciences Division, Department of Mechanical Engineering, Stanford University, 1983), includes a Galilean-invariant term called the modified Leonard stress tensor, and two nonlinear terms comprised of the products of the strain-rate and rotation-rate tensors for an improved representation of the subgrid-scale dissipation, backscatter and anisotropy effects. The model does not employ any ad hoc averaging or clipping procedures, and does not require the specification of a characteristic length scale; hence, it naturally avoids the ambiguities associated with defining a proper length scale for anisotropic grids. Results from the wall-resolved LES of flow past a Gaussian bump using the new model demonstrate improved prediction of skin-friction, flow separation, mean flow profiles and turbulent quantities when compared to implicit LES as well as explicit LES using the Vreman SGS model on the same grid.

Turbulence↗

High-Fidelity Simulation of Turbulent Flow Past a Gaussian Bump

A spanwise-periodic computation of a turbulent flow past a Gaussian bump is performed in the form of a hybrid direct numerical simulation and wall-resolved large-eddy simulation. A fourth-order spatially-accurate flow solver is employed to perform the simulation, using 10.2 billion grid points for a Reynolds number of 170000 based on the bump height. The key findings from the simulation are reported in the acceleration and deceleration flow regions associated with the bump shape. Significant anisotropy in the normal Reynolds stresses, along both the wall-normal and streamwise directions, is observed within the acceleration region. The ratio between the Reynolds shear stress and turbulent kinetic energy in that region also experiences significant deviations from the norms of a zero pressure gradient turbulent boundary layer. The chosen Reynolds number generates strong flow separation in the adverse pressure gradient region, which is in contrast with a previous simulation at half the Reynolds number that only indicated incipient separation. An internal layer generated in the acceleration region evolves into a free shear layer that develops in the deceleration region and separates. Proper modeling of this inner layer appears crucial to predict the flow separation. Surface curvature effects on the attached flow development are also discussed.

Turbulent Boundary Layer↗

Effect of Spatial Filtering in Implicit Large-Eddy Simulations of Separated Flows

The relatively high Reynolds number of turbulent flows encountered in various applications puts these problems well beyond the reach of direct numerical simulation (DNS) at present. Meanwhile, lower-fidelity Reynolds-averaged Navier-Stokes (RANS) calculations are known to be not accurate enough in complex problems, such as smooth-body flow separation and other flows involving highly-unsteady phenomena. Hence, given the current infeasibility of DNS and the unsatisfactory performance of RANS, intermediate techniques such as large-eddy simulation (LES) and hybrid RANS-LES, whose fidelity lie between RANS and DNS, have received much attention for application to various problems of practical importance. Modeling of the effect of missing scales on resolved scales, also known as subgrid-scale (SGS) modeling, is an important subject for LES. SGS models can be broadly categorized as explicit or implicit approaches. The explicit approach is based on an SGS model that explicitly appears in the governing equations expressed in the form of so-called “filtered Navier-Stokes equations”, which describe the evolution of the turbulence scales resolved by the LES grid. The effect of the scales unresolved by the grid is represented by the SGS model. The implicit modeling approach, on the other hand, does not employ an explicit model but instead treats the intrinsic dissipation of the numerical discretization scheme as an implicit SGS model. An LES without an explicit SGS model is commonly termed as an implicit LES (ILES). The relative merits of one SGS modeling approach over another is a subject of ongoing debate. We have opted to employ an ILES methodology, based on high-order compact finite-difference and spatial filtering schemes, in our recent investigations of separated flow problems [1, 2]. Further discussion of our preference of ILES over explicit LES is provided in Uzun and Malik [2]. The spatial filtering operation, described in the next section, is treated as an implicit SGS model for the ILES. Some observations made during the course of our recent investigations, which pointed out to excessive numerical dissipation in certain parts of the flowfield, prompted us to take a closer look at the potential effect of the spatial filter on ILES predictions. This technical note is therefore devoted to spatial filter effects in the context of a high Reynolds number, transonic shock-induced separated flow.

Ali Uzun↗

Simulations of a Turbulent Flow Subjected to Favorable and Adverse Pressure Gradients

This paper reports the results from a direct numerical simulation of an initially turbulent boundary layer passing over a wall-mounted “speed bump” geometry. The speed bump, represented in the form of a Gaussian distribution profile, generates a favorable pressure gradient region over the upstream half of the geometry, followed by an adverse pressure gradient over the downstream half. The boundary layer approaching the bump undergoes strong acceleration in the favorable pressure gradient region before experiencing incipient or very weak separation within the adverse pressure gradient region. These types of flows have proven to be particularly challenging to predict using lower-fidelity simulation tools based on various turbulence modeling approaches and warrant the use of the highest-fidelity simulation techniques. Simulation results are utilized to examine the key phenomena present in the flowfield, such as relaminarization/stabilization in the strong acceleration region succeeded by retransition to turbulence near the onset of adverse pressure gradient, incipient/weak separation, and development of internal layers where the sense of streamwise pressure gradient changes at the foot, apex and tail of the bump. The present direct numerical simulation is performed using a flow solver developed exclusively for graphics processing units, which is found to provide a significant speedup compared to an earlier solver optimized for central processing unit architectures.

Ali Uzun↗