Chemical reaction in a turbulent flow field with uniform velocity gradient
Chemical reactions in turbulent flow fields with uniform velocity gradient
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Chemical reactions in turbulent flow fields with uniform velocity gradient
Measurements of turbulent flow in a duct with a cross sectional transition from square to round are presented. Laser Doppler velocimetry was used to measure the mean velocity components, turbulence levels and shear stresses. The boundary layers at the inlet and exit of the transition were approximately 13 percent and 20 percent of the hydraulic diameter respectively, becoming thicker near the corner fillets. The development of secondary flow, of magnitudes up to 7 percent of the bulk velocity is shown. This flow is directed away from the corner fillets and along the periphery of the duct and is associated with the longitudinal curvature of the wall and the related pressure gradients.
A detailed experimental investigation has been carried out on the effects of different types of screens on turbulent flow, in particular turbulent boundary layers. The effect of a screen on a turbulent boundary layer is to give it a 'new lease of life'. The boundary layer turbulence is reorganized and the thickness reduced, thus making it less susceptible to separation. The aerodynamic properties of plastic screens are found to differ significantly from those of the conventional metal screens, evidently because of differences in the weaving properties. The 'overshoot' in mean velocity profile near the boudnary layer edge is shown to be a result of the effect of screen inclination on pressure drop coefficient. A more accurate formulation for the deflection coefficient of a screen is also proposed.
Grid convergence studies are performed to establish reference solutions for benchmark three dimensional turbulent flows in support of the ongoing turbulence model verification and validation e ort at the Turbulence Modeling Resource website curated by NASA. The bench- mark cases are a subsonic flow around a hemisphere cylinder and a transonic flow around the ONERA M6 wing with a sharp trailing edge. The study applies widely-used computational fluid dynamics codes developed and supported at the NASA Langley Research Center: FUN3D, USM3D, and CFL3D. Reference steady-state solutions are computed for the Reynolds-Averaged Navier-Stokes equations with the Spalart-Allmaras turbulence model on families of consistently-refined grids composed of different types of cells. Coarse-to- ne and code-to-code solution variation is described in detail.
Predictions of swirl, confinement and nozzle effects on confined turbulent flow are exhibited and compared with five-hole pitot-probe time-mean velocity measurements. Two sets of computations are given, one using the standard k-epsilon turbulence model and the other using a C sub mu formulation model deduced from recent six-orientation single-wire hot-wire measurements. Results confirm that the accuracy of the latter model is superior. To highlight the effects of confinement and exit nozzle area on this flow, three expansion ratios and two contraction ratios are used. Predictions are given for a full range of swirl strengths using measured inlet conditions for axial, radial and swirl velocity profiles. The predicted velocity profiles illustrate the large-scale effects of inlet swirl on flowfields. It appears that a strong contraction nozzle has a pronounced effect, on swirl flow cases, with discouragement of central recirculation zones, and forward flow in highly swirled vortex core regions. The expansion ratio value has large-scale effects on the size and location of the recirculation zones.
Baldwin-Lomax and kappa-epsilon turbulence models were modified for use in Navier-Stokes numerical computations of Mach 2.9 supersonic turbulent boundary layer flows along compression ramps. The computational results of Reynolds shear stress profiles were compared with experimental data. The Baldwin-Lomax model was modified to account for the Reynolds shear stress amplification within the flow field. A hybrid kappa-epsilon model with viscous sublayer turbulence treatment was constructed to predict the Reynolds shear stress profiles within the entire flow field. These modified turbulence models were effective for the computations of the surface pressure and the skin friction factor variations along an 8 deg ramp surface. The hybrid kappa-epsilon model could improve the predictions of the Reynolds shear stress profile and the skin friction factor near the corner of a 16 deg ramp.
The physical aspects governing accurate numerical simulation of turbulent flows having large regions of crossflow separation are re-examined. Time-accurate, three-dimensional fine-grid Navier-Stokes solutions were obtained for turbulent subsonic and supersonic flows over a slender ogive-cylinder body of revolution at large angles of attack. These flowfields are complex and contain regions of crossflow separation and an organized leeward-side vortex structure. An algebraic eddy-viscosity turbulence model has been modified to correctly account for the effects of the vortices on the underlying viscous layers. The numerical results show that the vortical flow structure of both subsonic and supersonic high incidence flows are similar, and can be accurately treated by the turbulence model. As the angle of attack is increased, the effectiveness of the model increases since the influence of the vortices on the underlying boundary layer decreases.
A technique for improving the numerical predictions of turbulent flows with the effect of streamline curvature is developed. Separated flows, the flow in a curved duct, and swirling flows are examples of flow fields where streamline curvature plays a dominant role. A comprehensive literature review on the effect of streamline curvature was conducted. New algebraic formulations for the eddy viscosity incorporating the kappa-epsilon turbulence model are proposed to account for various effects of streamline curvature. The loci of flow reversal of the separated flows over various backward-facing steps are employed to test the capability of the proposed turbulence model in capturing the effect of local curvature. The inclusion of the effect of longitudinal curvature in the proposed turbulence model is validated by predicting the distributions of the static pressure coefficients in an S-bend duct and in 180 degree turn-around ducts. The proposed turbulence model embedded with transverse curvature modification is substantiated by predicting the decay of the axial velocities in the confined swirling flows. The numerical predictions of different curvature effects by the proposed turbulence models are also reported.
The effects of prescribed wall motion on turbulent channel flow were examined with the objective of understanding the drag mechanisms and obtaining a possible means for drag reduction. A computer program for large-eddy simulation (LES) of turbulent flow in a channel was adapted to treat the case of time-varying wall motion using linearized boundary conditions. The code was applied to flow in a channel with prescribed wall motion on one wall. It was found that the results can be explained by a model based on a Stokes-like layer near the moving wall. This provides a theory which can be used to guide the selection of parameters. Considering both viscous and pressure drag components, it was found that net drag reduction is possible under certain conditions of wall motion. However, the work required to drive the wall motion can be greater than the reduced through-flow energy requirement. The major effects of the wall motion are restricted to the viscous sublayer of the turbulent flow.
In this paper, the turbulent flow in a serpentine with oblique ribs is investigated experimentally and by numerical simulations. The measurements are carried out by using Magnetic Resonance Velocimetry (MRV) and the simulations using the Immersed Boundary (IB) technique. A brief description of these two approaches is reported in following sections. The results are reported in terms of velocity distributions in various planes in the serpentine; differences between measurements and simulations are presented qualitatively and quantitatively. The study of the discrepancy allows us to identify areas of needed improvements in the turbulence modeling.
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.
The physical properties of the commonly used second-order closure models are examined theoretically for rotating turbulent flows. Comparisons are made with results which are a rigorous consequence of the Navier-Stokes equations for the problem of fully-developed turbulent channel flow in a rapidly rotating framework. It is demonstrated that all existing second-order closures yield spurious physical results for this test problem of rotating channel flow. In fact, the results obtained are shown to be substantially more unphysical than those obtained from the simpler k-epsilon and k-l models. Modifications in the basic structure of these second-order closure models are proposed which can alleviate this problem.
The physical properties of the commonly used second-order closure models are examined theoretically for rotating turbulent flows. Comparisons are made with results which are a rigorous consequence of the Navier-Stokes equations for the problem of a fully-developed turbulent channel flow in a rapidly rotating framework. It is demonstrated that all existing second-order closures yield spurious physical results for this test problem of rotating channel flow. In fact, the results obtained are shown to be substantially more unphysical than those obtained from the simpler K-epsilon and K-l models. Modifications in the basic structure of these second-order closure models are proposed which can alleviate this problem.
Substantial advances made over the past decade in the prediction of turbulent flows are discussed. There was extensive work in the development of turbulence models, particularly for use in boundary layer calculations. Basic aspects of several important methods based on partial differential equations for the mean velocity field and turbulence quantities, including the relationship between the methods and suggestions for future development were reviewed. Work on three-dimensional time-dependent large eddy simulations is discussed. The emphasis is on the hydrodynamics of incompressible flows, but sources for consideration of heat transfer and compressibility are mentioned.
A grid convergence study is performed to establish benchmark solutions for turbulent flows in three dimensions (3D) in support of turbulence-model verification campaign at the Turbulence Modeling Resource (TMR) website. The three benchmark cases are subsonic flows around a 3D bump and a hemisphere-cylinder configuration and a supersonic internal flow through a square duct. Reference solutions are computed for Reynolds Averaged Navier Stokes equations with the Spalart-Allmaras turbulence model using a linear eddy-viscosity model for the external flows and a nonlinear eddy-viscosity model based on a quadratic constitutive relation for the internal flow. The study involves three widely-used practical computational fluid dynamics codes developed and supported at NASA Langley Research Center: FUN3D, USM3D, and CFL3D. Reference steady-state solutions computed with these three codes on families of consistently refined grids are presented. Grid-to-grid and code-to-code variations are described in detail.
Energy decay and spectra measurements of grid-injection turbulent flows
A system for measuring momentum flux in a turbulently flowing fluid includes: a sensing apparatus for dynamically sensing the mainstream and the cross velocity components of the fluid, a transducer operative to provide two electrical output signals representative of the velocity components in the mainstream and in the cross direction, and signal processors to derive the Reynolds stress wave and the Reynolds stress.