Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulent flows”

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.

At least 919 records · Page 51

Contributions to the understanding of large-scale coherent structures in developing free turbulent shear flows

The physical problem of large-scale coherent structures in real, developing free turbulent shear flows are discussed from the point of view of a broader interpretation of the nonlinear aspects of hydrodynamic stability. Variations on the Amsden and Harlow problem are considered, and the role of linear theory in nonlinear problems is addressed. Spatially developing two-dimensional coherent structures and three-dimensional nonlinear effects in large-scale coherent mode interactions are considered.

Liu, J. T. C.↗

Numerical simulation of supersonic cone flow at high angle of attack

A conical symmetry assumption is applied to the full Navier-Stokes equations resulting in an equation set containing time and two coordinate directions as independent variables. The set is integrated by use of a finite difference technique for the particular case of sharp cones at incidence. Solutions are obtained and compared with experiment for auxiliary conditions corresponding to both laminar and turbulent flows. Closure for the turbulent flow case is provided by use of a scalor eddy viscosity model based on the mixing length hypothesis. Modifications to the eddy viscosity model were found which led to excellent surface pressure and surface flow direction agreement for turbulent flow at low supersonic Mach numbers. It is apparent that further work on the turbulence model is necessary for agreement at higher supersonic Mach numbers.

Mcrae, D. S.↗

Comparing Spray Characteristics from Reynolds Averaged Navier-Stokes (RANS) National Combustion Code (NCC) Calculations Against Experimental Data for a Turbulent Reacting Flow

Developing physics-based tools to aid in reducing harmful combustion emissions, like Nitrogen Oxides (NOx), Carbon Monoxide (CO), Unburnt Hydrocarbons (UHC s), and Sulfur Dioxides (SOx), is an important goal of aeronautics research at NASA. As part of that effort, NASA Glenn Research Center is performing a detailed assessment and validation of an in-house combustion CFD code known as the National Combustion Code (NCC) for turbulent reacting flows. To assess the current capabilities of NCC for simulating turbulent reacting flows with liquid jet fuel injection, a set of Single Swirler Lean Direct Injection (LDI) experiments performed at the University of Cincinnati was chosen as an initial validation data set. This Jet-A/air combustion experiment operates at a lean equivalence ratio of 0.75 at atmospheric pressure and has a 4 percent static pressure drop across the swirler. Detailed comparisons of NCC predictions for gas temperature and gaseous emissions (CO and NOx) against this experiment are considered in a previous work. The current paper is focused on detailed comparisons of the spray characteristics (radial profiles of drop size distribution and at several radial rakes) from NCC simulations against the experimental data. Comparisons against experimental data show that the use of the correlation for primary spray break-up implemented by Raju in the NCC produces most realistic results, but this result needs to be improved. Given the single or ten step chemical kinetics models, use of a spray size correlation gives similar, acceptable results

Iannetti, Anthony C.↗

Nonlinear Reynolds stress model for turbulent shear flows

A nonlinear algebraic Reynolds stress model, derived using the renormalization group, is applied to equilibrium homogeneous shear flow and fully developed flow in a square duct. The model, which is quadratically nonlinear in the velocity gradients, successfully captures the large-scale inhomogeneity and anisotropy of the flows studied. The ratios of normal stresses, as well as the actual magnitudes of the stresses are correctly predicted for equilibrium homogeneous shear flow. Reynolds normal stress anisotropy and attendant turbulence driven secondary flow are predicted for a square duct. Profiles of mean velocity and normal stresses are in good agreement with measurements. Very close to walls, agreement with measurements diminishes. The model has the benefit of containing no arbitrary constants; all values are determined directly from the theory. It seems that near wall behavior is influenced by more than the large scale anisotropy accommodated in the current model. More accurate near wall calculations may well require a model for anisotropic dissipation.

Barton, J. Michael↗

Low-Reynolds-number effects in a fully developed turbulent channel flow

Low-Reynolds-number effects are observed in the inner region of a fully developed turbulent channel flow, using data obtained either from experiments or by direct numerical simulations. The Reynolds-number influence is observed on the turbulence intensities and to a lesser degree on the average production and dissipation of the turbulent energy. In the near-wall region, the data confirm Wei and Willmarth's (1989) conclusion that the Reynolds stresses do not scale on wall variables. One of the reasons proposed to account for this behavior, namely, the 'geometry' effect or direct interaction between inner regions on opposite walls, was investigated in some detail by introducing temperature at one of the walls, both in experiment and simulation. Although the extent of penetration of thermal excursions into the opposite side of the channel can be significant at low Reynolds numbers, the contribution these excursions make to the Reynolds shear stress and the spanwise vorticity in the opposite wall region is negligible. In the inner region, spectra and cospectra of the velocity fluctuations u and v change rapidly with the Reynolds number, the variations being mainly confined to low wavenumbers in the u spectrum.

Antonia, R. A.↗

Analysis and computation of three-dimensional flow in strongly curved ducts

Numerical solutions are presented for three-dimensional laminar and turbulent flow in curved ducts of rectangular cross section and significant curvature. The analysis is based on a primary-secondary velocity decomposition in a given coordinate system, and leads to approximate governing equations which correct an a priori inviscid solution for viscous effects, secondary flows, total pressure distortion, heat transfer, and internal flow blockage and losses. Solution of the correction equations is accomplished as an initial-value problem in space using an implicit forward-marching technique. The overall solution procedure requires significantly less computational effort than Navier-Stokes algorithms. The present solution procedure is effective even with the extreme local mesh resolution which is necessary to resolve near-wall sublayer regions in turbulent flow calculations. Computed solutions for both laminar and turbulent flow compare very favorably with available analytical and experimental results.

Kreskovsky, J. P.↗

Large eddy simulation and direct numerical simulation of high speed turbulent reacting flows

The objective of this research is to make use of Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first phase of this research conducted within the past three years have been directed in several issues pertaining to intricate physics of turbulent reacting flows. In our previous 5 semi-annual reports submitted to NASA LaRC, as well as several technical papers in archival journals, the results of our investigations have been fully described. In this progress report which is different in format as compared to our previous documents, we focus only on the issue of LES. The reason for doing so is that LES is the primary issue of interest to our Technical Monitor and that our other findings were needed to support the activities conducted under this prime issue. The outcomes of our related investigations, nevertheless, are included in the appendices accompanying this report. The relevance of the materials in these appendices are, therefore, discussed only briefly within the body of the report. Here, results are presented of a priori and a posterior analyses for validity assessments of assumed Probability Density Function (PDF) methods as potential subgrid scale (SGS) closures for LES of turbulent reacting flows. Simple non-premixed reacting systems involving an isothermal reaction of the type A + B yields Products under both chemical equilibrium and non-equilibrium conditions are considered. A priori analyses are conducted of a homogeneous box flow, and a spatially developing planar mixing layer to investigate the performance of the Pearson Family of PDF's as SGS models. A posteriori analyses are conducted of the mixing layer using a hybrid one-equation Smagorinsky/PDF SGS closure. The Smagorinsky closure augmented by the solution of the subgrid turbulent kinetic energy (TKE) equation is employed to account for hydrodynamic fluctuations, and the PDF is employed for modeling the effects of scalar fluctuations. The implementation of the model requires the knowledge of the local values of the first two SGS moments. These are provided by additional modeled transport equations. In both a priori and a posteriori analyses, the predicted results are appraised by comparison with subgrid averaged results generated by DNS. Based on these results, the paths to be followed in future investigations are identified.

Adumitroaie, V.↗

Development of computational methods for heavy lift launch vehicles

The research effort has been focused on the development of an advanced flow solver for complex viscous turbulent flows with shock waves. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. A new computer program named CENS3D has been developed for viscous turbulent flows with discontinuities. Details of the code are described in Appendix A and Appendix B. With the developments of the numerical algorithm and dissipation model, the simulation of three-dimensional viscous compressible flows has become more efficient and accurate. The results of the research are expected to yield a direct impact on the design process of future liquid fueled launch systems.

Yoon, Seokkwan↗

Some structural features of unsteady separating turbulent shear flows

Some physical features of unsteady separating turbulent boundary layers are presented for practical Reynolds numbers and reduced frequencies for helicopter and turbomachinery flows. Upstream of detachment in moderate amplitude flows, the flow is quasi steady, i.e., the phase averaged flow is described by the steady-free stream flow structure. Results show that oscillation waveform and amplitude strongly influence the detached flow behavior.

Simpson, R. L.↗

A method for calculating strut and splitter plate noise in exit ducts: Theory and verification

Portions of a four-year analytical and experimental investigation relative to noise radiation from engine internal components in turbulent flow are summarized. Spectra measured for such airfoils over a range of chord, thickness ratio, flow velocity, and turbulence level were compared with predictions made by an available rigorous thin-airfoil analytical method. This analysis included the effects of flow compressibility and source noncompactness. Generally good agreement was obtained. This noise calculation method for isolated airfoils in turbulent flow was combined with a method for calculating transmission of sound through a subsonic exit duct and with an empirical far-field directivity shape. These three elements were checked separately and were individually shown to give close agreement with data. This combination provides a method for predicting engine internally generated aft-radiated noise from radial struts and stators, and annular splitter rings. Calculated sound power spectra, directivity, and acoustic pressure spectra were compared with the best available data. These data were for noise caused by a fan exit duct annular splitter ring, larger-chord stator blades, and turbine exit struts.

Fink, M. R.↗

Implementation of wall boundary conditions for transpiration in F3D thin-layer Navier-Stokes code

Numerical boundary conditions for mass injection/suction at the wall are incorporated in the thin-layer Navier-Stokes code, F3D. The accuracy of the boundary conditions and the code is assessed by a detailed comparison of the predictions of velocity distributions and skin-friction coefficients with exact similarity solutions for laminar flow over a flat plate with variable blowing/suction, and measurements for turbulent flow past a flat plate with uniform blowing. In laminar flow, F3D predictions for friction coefficient compare well with exact similarity solution with and without suction, but produces large errors at moderate-to-large values of blowing. A slight Mach number dependence of skin-friction coefficient due to blowing in turbulent flow is computed by F3D code. Predicted surface pressures for turbulent flow past an airfoil with mass injection are in qualitative agreement with measurements for a flat plate.

Kandula, M.↗

Numerical simulation of shear flows over wavy boundaries

The full time-dependent two dimensional Navier-Stokes equations are solved using spectral methods to achieve high spatial accuracy and high order time-splitting methods together with conformal mapping methods to allow simulation of flow over steep waves, including general turbulent Reynolds stresses. The programs are used to study the properties of laminar and turbulent flow over a flat plate and in various wavy geometries. Results are compared with both available theory and experiment. For laminar flow, good agreement is achieved. For turbulent flow, the agreement is less satisfactory and may be explained by the simplified turbulent modelling use in the calculations.

Orszag, S. A.↗

A comprehensive model to determine the effects of temperature and species fluctuations on reactions in turbulent reacting flows

The joint probability distribution function (pdf), which is a modification of the bivariate Gaussian pdf, is discussed and results are presented for a global reaction model using the joint pdf. An alternative joint pdf is discussed. A criterion which permits the selection of temperature pdf's in different regions of turbulent, reacting flow fields is developed. Two principal approaches to the determination of reaction rates in computer programs containing detailed chemical kinetics are outlined. These models represent a practical solution to the modeling of species reaction rates in turbulent, reacting flows.

Antaki, P. J.↗

Turbulent-like flows in quasi two-dimensional dense suspensions of motile colloids

Dense bacterial suspensions exhibit turbulent-like flows at low Reynolds numbers, driven by the activity of the microswimmers. In this study, we develop a model system to examine these dynamics using motile colloids that mimic bacterial locomotion. The colloids are powered by the Quincke instability, which causes them to spontaneously roll in a random-walk pattern when exposed to a square-wave electric field. We experimentally investigate the flow dynamics in dense suspensions of these Quincke random walkers under quasi two-dimensional conditions, where the particle size is comparable to the gap between the electrodes. The results suggest a scaling regime in the energy spectrum ∼k −4 at high wavenumbers, observed consistently across activity levels and particle concentrations. We observe that velocity time correlations decay within a single period of the square-wave field, yet an anti-correlation appears between successive field applications, indicative of a dynamic structural memory of the ensemble.

Luo, Rui [Northwestern Univ., Evanston, IL (United↗

Statistical analysis of near-wall structures in turbulent channel flow

The two-point correlation tensor is used to analyze near-wall structures in turbulent channel flow. Simulation results have been utilized to obtain the correlation tensor with sufficient spatial resolution to permit accurate differentiation for computing the vorticity correlation, and the resolution of sharp features such as the shear layer and its corresponding region of spanwise vorticity. Proper orthogonal decomposition and stochastic estimation were observed to yield similar results when appropriate conditions were employed for the estimation. Resulting structures were such that they would be detected by standard conditional sampling methods. The results of decomposing the vorticity field and the velocity field were significantly different, indicating that if a coherent structure is found that dominates the velocity fields, the curl of that velocity structure will not dominate the vorticity field.

Moser, Robert D.↗

Scalar dissipation and mixing in turbulent reacting flows

The mixing term that appears in probability density function-based formulations of the statistical behavior of turbulent reacting flows was studied with the use of direct numerical simulations of an irreversible, second-order, isothermal chemical reaction between initially segregated reactants in decaying, homogeneous turbulence. Three-dimensional, time-dependent pseudospectral calculations with 64-cubed Fourier modes were made as well as were Monte Carlo simulations of the joint concentration pdf equations for coalescence-dispersion (C/D) models. It was found that the development of the scalar microscales depends largely on the initial concentration distributions (both isotropic and nonisotropic distributions were used) and is relatively independent of the reaction rate constant. Joint concentration probability density functions were compared to predictions of D/D models that make use of parameters evaluated from the nonreacting or inert system.

Leonard, Andy D.↗

Automated and efficient local adaptive regression for principal component-based reduced-order modeling of turbulent reacting flows

Principal Component Analysis can be used to reduce the cost of Computational Fluid Dynamics simulations of turbulent reacting flows by reducing the dimensionality of the transported variables through projection of the thermochemical state onto a lower-dimensional manifold. However, because of the nonlinearity of the principal component source terms, nonlinear regression techniques must be utilized for the source terms in terms of the principal components. Unfortunately, widely available and utilized nonlinear regression techniques can have prohibitive computational requirements and/or accuracy that is highly dependent on user experience in ad hoc tuning of model architecture and hyperparameters. Here, in this work, a new nonlinear regression approach is proposed that is both computationally efficient and automated so does not require any user input. The approach is evaluated through a priori prediction of principal component source terms using data from a Direct Numerical Simulation of a turbulent nonpremixed n-heptane/air jet flame. In particular, the proposed framework consists of local regressions whose complexity is adapted according to the local nonlinearity of the data: local linear regression when accurate enough and local Artificial Neural Networks when nonlinear regression is required. The number of local clusters for local regression is determined automatically using the Davies-Bouldin index. In addition, Bayesian optimization is utilized for model training (i.e., to select the best architectures and hyperparameters of the nonlinear regressions in an unsupervised fashion), eliminating ad hoc hand-tuning and/or expensive grid searches. Overall, compared to a single, global neural network, the new local adaptive regression approach is shown to have comparable accuracy but 69% less training time due to the utilization of local linear regression and faster training of local neural networks.

42 ENGINEERING↗