Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence simulation”

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 1,297 records · Page 72

Evolution of a vortical structure associated with the bursting event in a channel flow

The temporal evolution of a horseshoe vortex structure in a channel is investigated by a numerical simulation. A spectral numerical method is employed to integrate the time-dependent, three-dimensional Navier-Stokes equations. The initial vortical structure is obtained by applying a conditional sampling technique to a data base generated from a direct simulation of a turbulent channel flow. The evolution of this vortical structure under the influence of the self-induced motion and the mean shear is presented. It is shown that the initial sheet-like vortical structure rolls up into a vortex tube as it is convected downstream. Turbulence structures associated with the vortex stretching are investigated. Production of vorticity stretching is high inside the vortex legs, although it is also substantial in the tip region and above the legs. High Reynolds shear stress is produced near the tip of the vortex.

Kim, J.↗

Plasma physics abstracts, 1 January - 31 December, 1986

Topics addressed include: ion-cyclotron waves; plasma waves; solar wind lithium releases; bow shock; Pi2 wave bursts; auroral kilometric radiation; ion energization; magnetic field corrections; electric fields; magnetospheric processes; electron acceleration; inner heliosphere; nightside auroral zone; computerized simulation; plasma wave turbulence; and magnetohydrodynamic waves in plasma sheets.

Gurnett, D. A.↗

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.↗

Geometry of premixed flames in three-dimensional turbulence

Constant density premixed flame propagation in three dimensional Navier-Stokes turbulence has been simulated. The zero-thickness flame model of Kerstein has been used. There are two aspects to this study: (1) adjustment of the large-scale strain rate in order to achieve a constant energy system; and (2) determination of flame curvature. The sampled distribution of curvature indicates that in most cases the flame has a cylindrical shape, with one curvature at least three times larger than the other. This implies that realistic chemical reactions in a flame-vortex interaction may be simulated in two dimensions.

Ashurst, William T.↗

Reynolds number influences in aeronautics

Reynolds number, a measure of the ratio of inertia to viscous forces, is a fundamental similarity parameter for fluid flows and therefore, would be expected to have a major influence in aerodynamics and aeronautics. Reynolds number influences are generally large, but monatomic, for attached laminar (continuum) flow; however, laminar flows are easily separated, inducing even stronger, non-monatomic, Reynolds number sensitivities. Probably the strongest Reynolds number influences occur in connection with transitional flow behavior. Transition can take place over a tremendous Reynolds number range, from the order of 20 x 10(exp 3) for 2-D free shear layers up to the order of 100 x 10(exp 6) for hypersonic boundary layers. This variability in transition behavior is especially important for complex configurations where various vehicle and flow field elements can undergo transition at various Reynolds numbers, causing often surprising changes in aerodynamics characteristics over wide ranges in Reynolds number. This is further compounded by the vast parameterization associated with transition, in that any parameter which influences mean viscous flow development (e.g., pressure gradient, flow curvature, wall temperature, Mach number, sweep, roughness, flow chemistry, shock interactions, etc.), and incident disturbance fields (acoustics, vorticity, particulates, temperature spottiness, even electro static discharges) can alter transition locations to first order. The usual method of dealing with the transition problem is to trip the flow in the generally lower Reynolds number wind tunnel to simulate the flight turbulent behavior. However, this is not wholly satisfactory as it results in incorrectly scaled viscous region thicknesses and cannot be utilized at all for applications such as turbine blades and helicopter rotors, nacelles, leading edge and nose regions, and High Altitude Long Endurance and hypersonic airbreathers where the transitional flow is an innately critical portion of the problem.

Bushnell, Dennis M.↗

Multiblock Navier-Stokes solutions about the F/A-18 wing-LEX-fuselage configuration

Three-dimensional thin-layer Navier-Stokes computations are presented for the F/A-18 configuration. The modeled configuration includes an accurate surface representation of the fuselage, leading-edge extension (LEX), and wing, both with and without leading-edge flap deflection. A multiblock structured grid strategy is employed to decompose the computational flowfield domain around the subject configuration. Steady-state solutions are obtained from an algorithm that solves the compressible Navier-Stokes equations with an upwind-biased, flux-difference splitting approach. The results presented are based on a fully turbulent flow assumption, simulating the high Reynolds number flow conditions that correspond to a recent FI A-18 flight experiment. Good agreements between the computations and the flight test results are obtained for both surface flow patterns as well as surface pressure distributions. Furthermore, a correlation between the computed LEX vortex-core and the flight test results, observed by way of smoke visualization, is also presented.

Aircraft configurations↗

Automated pattern eduction from turbulent flow diagnostics

The development of an automated technique for the eduction of 3-D spatial patterns in vector or scalar diagnostics was completed. The method is based on an iterative convolution between a trial pattern and the data field. It was applied to the analysis of low Reynolds number turbulent channel flow and homogeneous shear flow. The results yielded new information on the dominant flow structures in these flows, particularly with respect to the spatial relationships between various forms of organized motion. A particular application of the pattern eduction method, which is tentatively referred to as an 'adaptive wavelet transformation', is proposed with the objective of investigating the way turbulence structure changes with scale. Preliminary results using data from homogeneous turbulent shear flow simulations are presented. At the low Reynolds numbers of the simulations, there is no evidence of scale similarity. The small scales appear to be associated with the edges of the larger scale vortical structures.

Stretch, Derek D.↗

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.↗

On the upstream boundary of electron foreshocks in the solar wind

The upstream boundary of electron foreshocks is defined as the path of the fastest electrons reflected by collisionless shocks and moving along the magnetic field in the solar wind. Considerable levels of magnetic fluctuations are found in these regions of the solar wind, and their effect is to create both a broadening and a fine structure of the electron foreshock boundary. The magnetic structure is studied by means of a 3-D numerical simulation of a turbulent magnetic field. Enhanced, anomalous diffusion is found, (Delta x(exp 2)) varies as s(sup alpha), where alpha is greater than 1 for typical values of the parameters (here, Delta x(exp 2) is the mean square width of the tangent magnetic surface and s is the field line length). This corresponds to a Levy flight regime for the magnetic field line random walk, and allows very efficient electron propagation perpendicular to the magnetic field. Implications on the observations of planetary foreshocks and of the termination shock foreshock are considered.

Zimbardo, G.↗

Dynamic model with scale-dependent coefficients in the viscous range

The standard dynamic procedure is based on the scale-invariance assumption that the model coefficient C is the same at the grid and test-filter levels. In many applications this condition is not met, e.g. when the filter-length, delta, approaches the Kolmogorov scale, and C(delta approaches eta) approaches O. Using a priori tests, we show that the standard dynamic model yields the coefficient corresponding to the test-filter scale (alpha delta) instead of the grid-scale (delta). Several approaches to account for scale dependence are examined and/or tested in large eddy simulation of isotropic turbulence: (a) take the limit alpha approaches 1; (b) solve for two unknown coefficients C(Delta) and C(alpha delta) in the least-square-error formulation; (c) the 'bi-dynamic model', in which two test-filters (e.g. at scales 2(delta) and 4(delta) are employed to gain additional information on possible scale-dependence of the coefficient, and an improved estimate for the grid-level coefficient is obtained by extrapolation, (d) use theoretical predictions for the ratio C(alpha delta)/C(delta) and dynamically solve for C(delta). None of these options is found to be entirely satisfactory, although the last approach appears applicable to the viscous range.

Meneveau, C.↗

Friction Factor for Flow in Rectangular Ducts with One Side Rib-Roughened

Numerical simulations of incompressible turbulent flow through rectangular ducts with one side rib-roughened were performed to determine pressure drop. The "PHOENICS " software package was used for the computations, which required provision of a wall function for transverse rib-roughened surfaces. The present study was conducted in the range of 10(exp 5) less than or equal to Reynolds number less than or equal to 10(exp 7), 0.01 less than or equal to rib height to hydraulic diameter ratio less than or equal to 0.04, 10 less than or equal to pitch to rib height ratio less than or equal to 40. Using the numerical results, friction factor charts for various aspect ratios were generated. The numerical results agreed well with experimental data that was obtained for 10(exp 5) less than Reynolds less than 2 x 10(exp 5). In addition, a scheme for predicting friction factor using existing correlations for smooth and rough walls was developed.

Youn, B.↗

Convection in Slab and Spheroidal Geometries

Three-dimensional numerical simulations of compressible turbulent thermally driven convection, in both slab and spheroidal geometries, are reviewed and analyzed in terms of velocity spectra and mixing-length theory. The same ideal gas model is used in both geometries, and resulting flows are compared. The piecewise-parabolic method (PPM), with either thermal conductivity or photospheric boundary conditions, is used to solve the fluid equations of motion. Fluid motions in both geometries exhibit a Kolmogorov-like k(sup -5/3) range in their velocity spectra. The longest wavelength modes are energetically dominant in both geometries, typically leading to one convection cell dominating the flow. In spheroidal geometry, a dipolar flow dominates the largest scale convective motions. Downflows are intensely turbulent and up drafts are relatively laminar in both geometries. In slab geometry, correlations between temperature and velocity fluctuations, which lead to the enthalpy flux, are fairly independent of depth. In spheroidal geometry this same correlation increases linearly with radius over the inner 70 percent by radius, in which the local pressure scale heights are a sizable fraction of the radius. The effects from the impenetrable boundary conditions in the slab geometry models are confused with the effects from non-local convection. In spheroidal geometry nonlocal effects, due to coherent plumes, are seen as far as several pressure scale heights from the lower boundary and are clearly distinguishable from boundary effects.

Porter, David H.↗

An Engine Research Program Focused on Low Pressure Turbine Aerodynamic Performance

A comprehensive test program was performed in the Propulsion Systems Laboratory at the NASA Glenn Research Center, Cleveland Ohio using a highly instrumented Pratt and Whitney Canada PW 545 turbofan engine. A key objective of this program was the development of a high-altitude database on small, high-bypass ratio engine performance and operability. In particular, the program documents the impact of altitude (Reynolds Number) on the aero-performance of the low-pressure turbine (fan turbine). A second objective was to assess the ability of a state-of-the-art CFD code to predict the effect of Reynolds number on the efficiency of the low-pressure turbine. CFD simulation performed prior and after the engine tests will be presented and discussed. Key findings are the ability of a state-of-the art CFD code to accurately predict the impact of Reynolds Number on the efficiency and flow capacity of the low-pressure turbine. In addition the CFD simulations showed the turbulent intensity exiting the low-pressure turbine to be high (9%). The level is consistent with measurements taken within an engine.

Castner, Raymond↗

Numerical Study of Outlet Boundary Conditions for Unsteady Turbulent Internal Flows Using the NCC

This paper presents the results of studies on the outlet boundary conditions for turbulent internal flow simulations. Several outlet boundary conditions have been investigated by applying the National Combustion Code (NCC) to the configuration of a LM6000 single injector flame tube. First of all, very large eddy simulations (VLES) have been performed using the partially resolved numerical simulation (PRNS) approach, in which both the nonlinear and linear dynamic subscale models were employed. Secondly, unsteady Reynolds averaged Navier- Stokes (URANS) simulations have also been performed for the same configuration to investigate the effects of different outlet boundary conditions in the context of URANS. Thirdly, the possible role of the initial condition is inspected by using three different initial flow fields for both the PRNS/VLES simulation and the URANS simulation. The same grid is used for all the simulations and the number of mesh element is about 0.5 million. The main purpose of this study is to examine the long-time behavior of the solution as determined by the imposed outlet boundary conditions. For a particular simulation to be considered as successful under the given initial and boundary conditions, the solution must be sustainable in a physically meaningful manner over a sufficiently long period of time. The commonly used outlet boundary condition for steady Reynolds averaged Navier-Stokes (RANS) simulation is a fixed pressure at the outlet with all the other dependent variables being extrapolated from the interior. The results of the present study suggest that this is also workable for the URANS simulation of the LM6000 injector flame tube. However, it does not work for the PRNS/VLES simulation due to the unphysical reflections of the pressure disturbances at the outlet boundary. This undesirable situation can be practically alleviated by applying a simple unsteady convection equation for the pressure disturbances at the outlet boundary. The numerical results presented in this paper suggest that this unsteady convection of pressure disturbances at the outlet works very well for all the unsteady simulations (both PRNS/VLES and URANS) of the LM6000 single injector flame tube.

Liu, Nan-Suey↗

Analysis of Instabilities in Non-Axisymmetric Hypersonic Boundary Layers Over Cones

Hypersonic flows over circular cones constitute one of the most important generic configurations for fundamental aerodynamic and aerothermodynamic studies. In this paper, numerical computations are carried out for Mach 6 flows over a 7-degree half-angle cone with two different flow incidence angles and a compression cone with a large concave curvature. Instability wave and transition-related flow physics are investigated using a series of advanced stability methods ranging from conventional linear stability theory (LST) and a higher-fidelity linear and nonlinear parabolized stability equations (PSE), to the 2D eigenvalue analysis based on partial differential equations. Computed N factor distribution pertinent to various instability mechanisms over the cone surface provides initial assessments of possible transition fronts and a guide to corresponding disturbance characteristics such as frequency and azimuthal wave numbers. It is also shown that strong secondary instability that eventually leads to transition to turbulence can be simulated very efficiently using a combination of advanced stability methods described above.

Li, Fei↗

An Assessment of CFD/CSD Prediction State-of-the-Art by Using the HART II International Workshop Data

Over the past decade, there have been significant advancements in the accuracy of rotor aeroelastic simulations with the application of computational fluid dynamics methods coupled with computational structural dynamics codes (CFD/CSD). The HART II International Workshop database, which includes descent operating conditions with strong blade-vortex interactions (BVI), provides a unique opportunity to assess the ability of CFD/CSD to capture these physics. In addition to a baseline case with BVI, two additional cases with 3/rev higher harmonic blade root pitch control (HHC) are available for comparison. The collaboration during the workshop permits assessment of structured, unstructured, and hybrid overset CFD/CSD methods from across the globe on the dynamics, aerodynamics, and wake structure. Evaluation of the plethora of CFD/CSD methods indicate that the most important numerical variables associated with most accurately capturing BVI are a two-equation or detached eddy simulation (DES)-based turbulence model and a sufficiently small time step. An appropriate trade-off between grid fidelity and spatial accuracy schemes also appears to be pertinent for capturing BVI on the advancing rotor disk. Overall, the CFD/CSD methods generally fall within the same accuracy; cost-effective hybrid Navier-Stokes/Lagrangian wake methods provide accuracies within 50% the full CFD/CSD methods for most parameters of interest, except for those highly influenced by torsion. The importance of modeling the fuselage is observed, and other computational requirements are discussed.

Smith, Marilyn J.↗

Flow Phenomena in the Very Near Wake of a Flat Plate with a Circular Trailing Edge

The very near wake of a flat plate with a circular trailing edge, exhibiting pronounced shedding of wake vortices, is investigated with data from a direct numerical simulation. The separating boundary layers are turbulent and statistically identical thus resulting in a wake that is symmetric in the mean. The focus here is on the instability of the detached shear layers, the evolution of rib-vortex induced localized regions of reverse flow that detach from the main body of reverse flow in the trailing edge region and convect downstream, and phaseaveraged velocity statistics in the very near wake. The detached shear layers are found to exhibit unstable behavior intermittently, including the development of shear layer vortices as in earlier cylinder flow investigations with laminar separating boundary layers. Only a small fraction of the separated turbulent boundary layers undergo this instability, and form the initial shed vortices. Pressure spectra within the shear layers show a broadband peak at a multiple of shedding frequency. Phase-averaged intensity and shear stress distributions of the randomly fluctuating component of velocity are compared with those obtained in the near wake. The distributions of the production terms in the transport equations for the turbulent stresses are also provided.

Turbulence↗

Tensor-Product Preconditioners for Higher-Order Space-Time Discontinuous Galerkin Methods

space-time discontinuous-Galerkin spectral-element discretization is presented for direct numerical simulation of the compressible Navier-Stokes equat ions. An efficient solution technique based on a matrix-free Newton-Krylov method is developed in order to overcome the stiffness associated with high solution order. The use of tensor-product basis functions is key to maintaining efficiency at high order. Efficient preconditioning methods are presented which can take advantage of the tensor-product formulation. A diagonalized Alternating-Direction-Implicit (ADI) scheme is extended to the space-time discontinuous Galerkin discretization. A new preconditioner for the compressible Euler/Navier-Stokes equations based on the fast-diagonalization method is also presented. Numerical results demonstrate the effectiveness of these preconditioners for the direct numerical simulation of subsonic turbulent flows.

Diosady, Laslo T.↗