Search NASA⌕ Search

SEARCH · Search NASA

Results for “COMPRESSIBLE TURBULENT FLOW”

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 127 records · Page 7

The analysis and simulation of compressible turbulence

Compressible turbulent flows at low turbulent Mach numbers are considered. Contrary to the general belief that such flows are almost incompressible (i.e., the divergence of the velocity field remains small for all times), it is shown that even if the divergence of the initial velocity field is negligibly small, it can grow rapidly on a nondimensional time scale which is the inverse of the fluctuating Mach number. An asymptotic theory which enables one to obtain a description of the flow in terms of its divergence-free and vorticity-free components has been developed to solve the initial-value problem. As a result, the various types of low Mach number turbulent regimes have been classified with respect to the initial conditions. Formulae are derived that accurately predict the level of compressibility after the initial transients have disappeared. These results are verified by extensive direct numerical simulations of isotropic turbulence.

Erlebacher, Gordon↗

Numerical solutions of Navier-Stokes equations for compressible turbulent two/three dimensional flows in terminal shock region of an inlet/diffuser

The multidimensional ensemble averaged compressible time dependent Navier Stokes equations in conjunction with mixing length turbulence model and shock capturing technique were used to study the terminal shock type of flows in various flight regimes occurring in a diffuser/inlet model. The numerical scheme for solving the governing equations is based on a linearized block implicit approach and the following high Reynolds number calculations were carried out: (1) 2 D, steady, subsonic; (2) 2 D, steady, transonic with normal shock; (3) 2 D, steady, supersonic with terminal shock; (4) 2 D, transient process of shock development and (5) 3 D, steady, transonic with normal shock. The numerical results obtained for the 2 D and 3 D transonic shocked flows were compared with corresponding experimental data; the calculated wall static pressure distributions agree well with the measured data.

Liu, N. S.↗

Adaptive unstructured triangular mesh generation and flow solvers for the Navier-Stokes equations at high Reynolds number

A method for generating high quality unstructured triangular grids for high Reynolds number Navier-Stokes calculations about complex geometries is described. Careful attention is paid in the mesh generation process to resolving efficiently the disparate length scales which arise in these flows. First the surface mesh is constructed in a way which ensures that the geometry is faithfully represented. The volume mesh generation then proceeds in two phases thus allowing the viscous and inviscid regions of the flow to be meshed optimally. A solution-adaptive remeshing procedure which allows the mesh to adapt itself to flow features is also described. The procedure for tracking wakes and refinement criteria appropriate for shock detection are described. Although at present it has only been implemented in two dimensions, the grid generation process has been designed with the extension to three dimensions in mind. An implicit, higher-order, upwind method is also presented for computing compressible turbulent flows on these meshes. Two recently developed one-equation turbulence models have been implemented to simulate the effects of the fluid turbulence. Results for flow about a RAE 2822 airfoil and a Douglas three-element airfoil are presented which clearly show the improved resolution obtainable.

Ashford, Gregory A.↗

Tetrahedral-Mesh Simulation of Turbulent Flows with the Space-Time Conservative Schemes

Direct numerical simulations of turbulent flows are predominantly carried out using structured, hexahedral meshes despite decades of development in unstructured mesh methods. Tetrahedral meshes offer ease of mesh generation around complex geometries and the potential of an orientation free grid that would provide un-biased small-scale dissipation and more accurate intermediate scale solutions. However, due to the lack of consistent multi-dimensional numerical formulations in conventional schemes for triangular and tetrahedral meshes at the cell interfaces, numerical issues exist when flow discontinuities or stagnation regions are present. The space-time conservative conservation element solution element (CESE) method - due to its Riemann-solver-free shock capturing capabilities, non-dissipative baseline schemes, and flux conservation in time as well as space - has the potential to more accurately simulate turbulent flows using unstructured tetrahedral meshes. To pave the way towards accurate simulation of shock/turbulent boundary-layer interaction, a series of wave and shock interaction benchmark problems that increase in complexity, are computed in this paper with triangular/tetrahedral meshes. Preliminary computations for the normal shock/turbulence interactions are carried out with a relatively coarse mesh, by direct numerical simulations standards, in order to assess other effects such as boundary conditions and the necessity of a buffer domain. The results indicate that qualitative agreement with previous studies can be obtained for flows where, strong shocks co-exist along with unsteady waves that display a broad range of scales, with a relatively compact computational domain and less stringent requirements for grid clustering near the shock. With the space-time conservation properties, stable solutions without any spurious wave reflections can be obtained without a need for buffer domains near the outflow/farfield boundaries. Computational results for the isotropic turbulent flow decay, at a relatively high turbulent Mach number, show a nicely behaved spectral decay rate for medium to high wave numbers. The high-order CESE schemes offer very robust solutions even with the presence of strong shocks or widespread shocklets. The explicit formulation in conjunction with a close to unity theoretical upper Courant number bound has the potential to offer an efficient numerical framework for general compressible turbulent flow simulations with unstructured meshes.

Chang, Chau-Lyan↗

Turbulence modeling

The performance of existing two-equation eddy viscosity models was examined. An effort was made to develop better models for near-wall turbulence using direct numerical simulations of plane channel and boundary layer flows. The asymptotic near-wall behavior of turbulence was used to examine the problems of current second order closure models and develop new models with the correct near-wall behavior. Rapid Distortion Theory was used to analytically study the effects of mean deformation on turbulence, obtain analytical solutions for the spectrum tensor, Reynolds stress tensor, anisotropy tensor and its invariants, which can be used in the turbulence model development. The potential of the renormalization group theory in turbulence modeling was studied, as well as compressible turbulent flows, and modeling of bypass transition.

Shih, Tsan-Hsing↗

An experimental and analytical investigation of axisymmetric diffusers

A finite difference computer program for turbulent compressible flow was used to establish the performance of several diffuser shapes for experimental testing. The diffusers were designed to have a linear change in Mach number, a linear change in pressure, or a curvature fitted by a quadratic equation. Testing was performed with M = 0.1 to 0.9 with and without boundary layer bleed. Above M = 0.6, data were obtained with a normal shock upstream of the diffuser entrance. Peak static pressure recovery occurred with a diffuser inlet M0.75. The quadratic diffuser yielded the highest total pressure recovery.

Povinelli, L. A.↗

An experimental and analytical investigation of axisymmetric diffusers

A finite difference computer program for turbulent compressible flow was used to establish the performance of several diffuser shapes for experimental testing. The diffusers were designed to have a linear change in Mach number, a linear change in pressure or a curvature fitted by a quadratic equation. Testing was performed with M = 0.1 to 0.9 with and without boundary layer bleed. Above M = 0.6, data were obtained with a normal shock upstream of the diffuser entrance. Peak static pressure recovery occurred with a diffuser inlet M-0.75. The quadratic diffuser yielded the highest total pressure recovery.

Povinelli, L. A.↗

An analytical and experimental study of a short s-shaped subsonic diffuser of a supersonic inlet

A subscale HiMAT forebody and inlet was investigated over a range of Mach numbers to 1.4. The inlet exhibited a transitory separation within the diffuser but steady state data indicated reattachment at the diffuser exit. A finite difference procedure for turbulent compressible flow in axisymmetric ducts was used to successfully model the HiMAT duct flow. The analysis technique was further used to estimate the initiation of separation and delineate the steady and unsteady flow regimes in similar S-shaped ducts.

Neumann, H. E.↗

An analytical and experimental study of a short S-shaped subsonic diffuser of a supersonic inlet

An experimental investigation of a subscale HiMat forebody and inlet was conducted over a range of Mach numbers to 1.4. The inlet exhibited a transitory separation within the diffuser but steady state data indicated reattachment at the diffuser exit. A finite difference procedure for turbulent compressible flow in axisymmetric ducts was used to successfully model the HiMAT duct flow. The analysis technique was further used to estimate the initiation of separation and delineate the steady and unsteady flow regimes in similar S-shaped ducts.

Neumann, H. E.↗

Trends and pacing items in computational aerodynamics

A perspective is presented of trends in computational aerodynamics, and of important technology development items that pace future advanced applications. From a survey of AIAA Journal papers published during the past two decades, the growth trends and the progressively increasing emphasis on code development for viscous, compressible, turbulent flow are illustrated. These trends are reflected in the chronology of introduction by the aerospace industry of new computational methods in aircraft design. Key pacing items outlined are: automatic grid generation for nonlinear inviscid computations; advanced computers, improved efficiency of numerical methods, and improved turbulence models for Reynolds-averaged Navier-Stokes computations; advanced computers, time-dependent three-dimensional law-of-the-wall, code development, improved efficiency of numerical methods, and improved subgrid-scale turbulence modeling for large eddy simulations.

Chapman, D. R.↗

Three-dimensional unsteady transonic viscous-inviscid interaction using the Euler and boundary-layer equations

The objective of this study is the development of a numerical technique which can provide three-dimensional, time-accurate, compressible, turbulent flow solutions in a practical and relatively economical way. The approach taken is that of the method of viscous-inviscid interaction. The Euler equations are assumed to govern the outer inviscid portion of the flow, and the viscous layer close to the solid wall is described by a set of integral boundary-layer equations. The viscous solutions are obtained in a direct fashion with a weighted-average phase error scheme. The method of equivalent sources is used for viscous-inviscid coupling. Steady-state and unsteady computations for an AGARD airfoil and a wing show that satisfactory engineering solutions can be obtained for attached, high Reynolds number flows using this method. Quasi-unsteady interactions are shown to produce similar results to those provided by true-unsteady interactions. Considerable computer resources can be saved for some cases by using quasi-unsteady interactions.

Whitfield, David L.↗

Three-dimensional compressible turbulent computations for a nondiffusing S-duct

The PARC3D code was used to compute the compressible turbulent flow within a three dimensional, nondiffusing S-duct. A frame of reference is provided for future computational fluid dynamics studies of internal flows with strong secondary flows and provides an understanding of the performance characteristics of a typical S-duct with attached flow. The predicted results, obtained with both H- and O-grids, are compared with the experimental wall pressure, static and total pressure fields, and velocity vectors. Additionally, computed boundary layer thickness, velocity profiles in wall coordinates, and skin friction values are presented.

Harloff, G. J.↗

Pseudospectral simulation of compressible turbulence using logarithmic variables

The direct numerical simulation of dissipative, highly compressible turbulent flow is performed using a pseudospectral Fourier technique. The governing equations are cast in a form where the important physical variables are the fluid velocity and the natural logarithms of the fluid density and temperature. Bulk viscosity is utilized to model polyatomic gases more accurately and to ensure numerical stability in the presence of strong shocks. Numerical examples include three-dimensional supersonic homogeneous turbulence and two-dimensional shock-turbulence interactions.

Shebalin, John V.↗

Carrying the mass flux terms exactly in the first and second moment equations of compressible turbulence

In compressible turbulence models, it is assumed that the Favre-mean velocities are suitable approximations to the Reynolds-mean velocities in order to close unknown terms. This neglects, in the mean momentum and energy equations, the contribution to the stress and work terms by the mean of the fluctuating Favre velocity, a quantity proportional to the turbulent mass flux. As the stress and work terms do not introduce any new unknown correlations requiring closure in either k-epsilon or Reynolds stress closures and because the exact form of the terms can, with little additional work, be carried there is no need to make any modeling assumptions. In the Reynolds stress equations the viscous terms appear naturally in Reynolds variables while the problem is posed in Favre variables. In the process of splitting the viscous terms into the viscous transport terms, carried in Favre variables, and the dissipation terms, carried in Reynolds variables, important contributions from the mass flux appear. The accurate accounting of these terms is important for any consistent near wall modeling and the retention of the mass flux terms is important in complex compressible turbulent flows.

Ristorcelli, J. R.↗

Discontinuous Galerkin Methods for Turbulence Simulation

A discontinuous Galerkin (DG) method is formulated, implemented, and tested for simulation of compressible turbulent flows. The method is applied to turbulent channel flow at low Reynolds number, where it is found to successfully predict low-order statistics with fewer degrees of freedom than traditional numerical methods. This reduction is achieved by utilizing local hp-refinement such that the computational grid is refined simultaneously in all three spatial coordinates with decreasing distance from the wall. Another advantage of DG is that Dirichlet boundary conditions can be enforced weakly through integrals of the numerical fluxes. Both for a model advection-diffusion problem and for turbulent channel flow, weak enforcement of wall boundaries is found to improve results at low resolution. Such weak boundary conditions may play a pivotal role in wall modeling for large-eddy simulation.

Collis, S. Scott↗

Numerical solution of the Navier-Stokes equations for compressible turbulent two/three dimensional flows in the terminal shock region of an inlet/diffuser

The multidimensional, ensemble-averaged, compressible, time-dependent Navier-Stokes equations in conjunction with mixing length turbulence model and shock capturing technique have been used to study the terminal shock type of flows in various flight regimes occurring in a diffuser/inlet model. The numerical scheme for solving the governing equations is based on a linearized block implicit approach and the following high Reynolds number calculations have been carried out: (1) 2-D, steady, subsonic; (2) 2-D, steady, transonic with normal shock, (3) 2-D, steady, supersonic with terminal shock, (4) 2-D, transient process of shock development and (5) 3-D, steady, transonic with normal shock. The numerical results obtained for the 2-D and 3-D transonic shocked flows have been compared with corresponding experimental data; the calculated wall static pressure distributions agree well with the measured data. The predicted transient responses of the flow to externally applied impulsive disturbance are consistent with results obtained via asymptotic analysis.

Liu, N.-S.↗

Direct simulation of compressible turbulence in a shear flow

The purpose of this study is to investigate compressibility effects on the turbulence in homogeneous shear flow. It is found that the growth of the turbulent kinetic energy decreases with increasing Mach number, a phenomenon similar to the reduction of turbulent velocity intensities observed in experiments on supersonic free shear layers. An examination of the turbulent energy budget shows that both the compressible dissipation and the pressure-dilatation contribute to the decrease in the growth of kinetic energy. The pressure-dilatation is predominantly negative in homogeneous shear flow, in contrast to its predominantly positive behavior in isotropic turbulence. The different signs of the pressure-dilatation are explained by theoretical consideration of the equations for the pressure variance and density variance.

Sarkar, S.↗

A compressible Navier-Stokes code for turbulent flow modeling

An implicit, finite volume code for solving two dimensional, compressible turbulent flows is described. Second order upwind differencing of the inviscid terms of the equations is used to enhance stability and accuracy. A diagonal form of the implicit algorithm is used to improve efficiency. Several zero and two equation turbulence models are incorporated to study their impact on overall flow modeling accuracy. Applications to external and internal flows are discussed.

Coakley, T. J.↗