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 613 records · Page 34

Numerical simulation of turbomachinery flows with advanced turbulence models

A three dimensional full Navier-Stokes (FNS) code is used to simulate complex turbomachinery flows. The code incorporates an explicit multistep scheme and solves a conservative form of the density averaged continuity, momentum, and energy equations. A compressible low Reynolds number form of the k-epsilon turbulence model, and a q-omega model and an algebraic Reynolds stress model have been incorporated in a fully coupled manner to approximate Reynolds stresses. The code is used to predict the viscous flow field in a backswept transonic centrifugal compressor for which laser two focus data is available. The code is also used to simulate the tip clearance flow in a cascade. The code has been extended to include unsteady Euler solutions for predicting the unsteady flow through a cascade due to incoming wakes, simulating rotor-stator interactions.

Lakshminarayana, B.↗

Simulations using a 3D-moment turbulence closure in a two-dimensional cloud model - Comparisons with observations

A method for deep cumulus cloud simulations is presented, in which the plumes or clouds are simulated using a two-dimensional axisymmetric numerical cloud model with third-moment turbulence closure. This method is based on the assumption that the only uncertainty in the turbulence closure is in the turbulent length scale specification. Therefore, the length scale is varied, and, using different length scales, the simulated plumes are compared to the observed ones, making it possible to optimize the length scale prescription.

Krueger, Steven K.↗

Least Squares Shadowing Sensitivity Analysis of Chaotic Flow Around a Two-Dimensional Airfoil

Gradient-based sensitivity analysis has proven to be an enabling technology for many applications, including design of aerospace vehicles. However, conventional sensitivity analysis methods break down when applied to long-time averages of chaotic systems. This breakdown is a serious limitation because many aerospace applications involve physical phenomena that exhibit chaotic dynamics, most notably high-resolution large-eddy and direct numerical simulations of turbulent aerodynamic flows. A recently proposed methodology, Least Squares Shadowing (LSS), avoids this breakdown and advances the state of the art in sensitivity analysis for chaotic flows. The first application of LSS to a chaotic flow simulated with a large-scale computational fluid dynamics solver is presented. The LSS sensitivity computed for this chaotic flow is verified and shown to be accurate, but the computational cost of the current LSS implementation is high.

Blonigan, Patrick J.↗

Turbulence Model Effects on RANS Simulations of the HIFiRE Flight 2 Ground Test Configurations

The Wind-US Reynolds-averaged Navier-Stokes solver was applied to the Hypersonic International Flight Research Experimentation (HIFiRE) Flight 2 scramjet ground test configuration. Two test points corresponding to flight Mach numbers of 5.9 and 8.9 were examined. The emphasis was examining turbulence model effects on the prediction of flow path pressures. Three variants of the Menter k-omega turbulence model family were investigated. These include the baseline (BSL) and shear stress transport (SST) as well as a modified SST model where the shear stress limiter was altered. Variations in the turbulent Schmidt number were also considered. Choice of turbulence model had a substantial effect on prediction of the flow path pressures. The BSL model produced the highest pressures and the SST model produced the lowest pressures. As expected, the settings for the turbulent Schmidt number also had significant effects on predicted pressures. Small values for the turbulent Schmidt number enabled more rapid mass transfer, faster combustion, and in turn higher flowpath pressures. Optimal settings for turbulence model and turbulent Schmidt number were found to be rather case dependent, as has been concluded in other scramjet investigations.

Georgiadis, Nicholas J.↗

Atmospheric effects on the Space Shuttle approach and landing

The present work describes piloted simulations of Space Shuttle approach and landing where the pilot either flies the aircraft or monitors a computer-controlled flight with the option to take control. The simulated runs encompassed approach to the airport, acquisition of the glide slope and localizer, and accomplishing final flare and touchdown. The degree of difficulty associated with the flight task was imposed through modeled atmospheric winds and superimposed atmospheric turbulence. The simulator models the Orbiter vehicle with the aid of an all-digital computer complex and a fixed-base crew station. Comparison between piloted and computer-controlled responses to unexpected atmospheric conditions was obtained by recording data from instrument displays available to the pilot.

Ragsdale, W. A.↗

Low-dimensional chaos in turbulence

Direct numerical simulations are being performed on two different fluid flows in an attempt to discover the mechanism underlying the transition to turbulence in each. The first system is Taylor-Couette flow; the second, two-dimensional flow over an airfoil. Both flows exhibit a gradual transition to high-dimensional turbulence through low-dimensional chaos. The hope is that the instabilities leading to chaos will be easier to relate to physical processes in this case, and that the understanding of these mechanisms can then be applied to a wider array of turbulent systems.

Vastano, John A.↗

Turbulence flight director analysis and preliminary simulation

A control column and trottle flight director display system is synthesized for use during flight through severe turbulence. The column system is designed to minimize airspeed excursions without overdriving attitude. The throttle system is designed to augment the airspeed regulation and provide an indication of the trim thrust required for any desired flight path angle. Together they form an energy management system to provide harmonious display indications of current aircraft motions and required corrective action, minimize gust upset tendencies, minimize unsafe aircraft excursions, and maintain satisfactory ride qualities. A preliminary fixed-base piloted simulation verified the analysis and provided a shakedown for a more sophisticated moving-base simulation to be accomplished next. This preliminary simulation utilized a flight scenario concept combining piloting tasks, random turbulence, and discrete gusts to create a high but realistic pilot workload conducive to pilot error and potential upset. The turbulence director (energy management) system significantly reduced pilot workload and minimized unsafe aircraft excursions.

Johnson, D. E.↗

Large Eddy Simulation of Gravitational Effects on Transitional and Turbulent Gas-Jet Diffusion Flames

The basic objective of this work is to assess the influence of gravity on "the compositional and the spatial structures" of transitional and turbulent diffusion flames via large eddy simulation (LES), and direct numerical simulation (DNS). The DNS is conducted for appraisal of the various closures employed in LES, and to study the effect of buoyancy on the small scale flow features. The LES is based on our "filtered mass density function"' (FMDF) model. The novelty of the methodology is that it allows for reliable simulations with inclusion of "realistic physics." It also allows for detailed analysis of the unsteady large scale flow evolution and compositional flame structure which is not usually possible via Reynolds averaged simulations.

Givi, Peyman↗

Non-local transport in turbulent MHD convection

The nonlocal non-diffusive transport of passive scalars in turbulent magnetohydrodynamic (MHD) convection is investigated using transilient matrices. These matrices describe the probability that a tracer particle beginning at one position in a flow will be advected to another position after some time. A method for the calculation of these matrices from simulation data which involves following the trajectories of passive tracer particles and calculating their transport statistics, is presented. The method is applied to study the transport in several simulations of turbulent, rotating, three dimensional compressible, penetrative MDH convection. Transport coefficients and other diagnostics are used to quantify the transport, which is found to resemble advection more closely than diffusion. Some of the results are found to have direct relevance to other physical problems, such as the light element depletion in sun-type stars. The large kurtosis found for downward moving particles at the base of the convection zone implies several extreme events.

Miesch, Mark↗

Numerical simulation of the flow and fuel-air mixing in an axisymmetric piston-cylinder arrangement

The implicit factored method of Beam and Warming was employed to describe the flow and the fuel-air mixing in an axisymmetric piston-cylinder configuration during the intake and compression strokes. The governing equations were established on the basis of laminar flow. The increased mixing due to turbulence was simulated by appropriately chosen effective transport properties. Calculations were performed for single-component gases and for two-component gases and for two-component gas mixtures. The flow field was calculated as functions of time and position for different geometries, piston speeds, intake-charge-to-residual-gas-pressure ratios, and species mass fractions of the intake charge. Results are presented in graphical form which show the formation, growth, and break-up of those vortices which form during the intake stroke and the mixing of fuel and air throughout the intake and compression strokes. It is shown that at bore-to-stroke ratio of less than unity, the vortices may break-up during the intake stroke. It is also shown that vortices which do not break-up during the intake stroke coalesce during the compression stroke. The results generated were compared to existing numerical solutions and to available experimental data.

Shih, T. I. P.↗

Simulation and modeling of homogeneous, compressed turbulence

Low Reynolds number homogeneous turbulence undergoing low Mach number isotropic and one-dimensional compression was simulated by numerically solving the Navier-Stokes equations. The numerical simulations were performed on a CYBER 205 computer using a 64 x 64 x 64 mesh. A spectral method was used for spatial differencing and the second-order Runge-Kutta method for time advancement. A variety of statistical information was extracted from the computed flow fields. These include three-dimensional energy and dissipation spectra, two-point velocity correlations, one-dimensional energy spectra, turbulent kinetic energy and its dissipation rate, integral length scales, Taylor microscales, and Kolmogorov length scale. Results from the simulated flow fields were used to test one-point closure, two-equation models. A new one-point-closure, three-equation turbulence model which accounts for the effect of compression is proposed. The new model accurately calculates four types of flows (isotropic decay, isotropic compression, one-dimensional compression, and axisymmetric expansion flows) for a wide range of strain rates.

Wu, C. T.↗

An explicit Runge-Kutta method for 3D turbulent incompressible flows

A computer code has been developed to solve for the steady-state solution of the 3D incompressible Reynolds-averaged Navier-Stokes equations. The approach is based on the cell-center, central-difference, finite-volume formulation and an explicit one-step, multistage Runge-Kutta time-stepping scheme. The Baldwin-Lomax turbulence model is used. Techniques to accelerate the rate of convergence to a steady-state solution include the preconditioned method, the local time stepping, and the implicit residual smoothing. Improvements in computational efficiency have been demonstrated in several areas. This numerical procedure has been used to simulate the turbulent horseshoe vortex flow around an airfoil/flat-plate juncture.

Sung, Chao-Ho↗

Laboratory and theoretical models of planetary-scale instabilities and waves

Research work is proceeding in theoretical, numerical, and experimental geophysical fluid dynamics leading up to a reflight of the GFFC (Geophysical Fluid Flow Cell Experiment) on USML-2. The work is intended not only to generate ideas for future space experiments, but to provide fundamental results concerned with nonlinear and chaotic properties of thermal convection and baroclinic waves in terrestrial and planetary atmospheres. The major efforts are focussed on thermal convection in a rapidly rotating annulus relevant to Jovian atmospheric dynamics, and on the chaotic behavior of baroclinic waves relevant to the Earth's atmosphere. The approach, in preparation for USML-2, is primarily theoretical and numerical. Mechanistic process models are solved numerically in order to identify physical mechanisms that may be observed in the GFFC, and which are important in real geophysical applications. The results from numerical simulations of geophysical fluid flow (subject to rotation and stratification) are compared with previous GFFC experiments on Spacelab-3 and with existing and proposed terrestrial laboratory experiments of various types. Pattern recognition algorithms have been employed to generate low-dimensional descriptions of the highly nonlinear and turbulent numerical simulations. Such empirically truncated descriptions provide for simplified but robust physical interpretations of the dynamics, as well as yielding highly efficient computations of these chaotic flows.

Hart, John E.↗

High-Order Simulation of Induced Disturbance in a Mach 6 Boundary Layer

Simulations of the early evolution of induced unstable disturbances in a Mach 6 hypersonic boundary layer are presented for the purposes of validating a high-order discontinuous Galerkin Navier-Stokes simulation code. The simulations are modeled after experiments performed in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. In these experiments, a forcing mechanism was employed to induce reproducible disturbances in a hypersonic boundary layer providing for the controlled study of the growth and breakdown of these disturbances into turbulent spots. Simulations revealed that the form and strength of the excitation can greatly influence the growth of the disturbance. In particular, at large forcing amplitude, the simulated forcing produces large advecting transients that appear to enhance the growth of the wave packet, relative to that of low amplitude forcing. Simulation results with large amplitude forcing agree well with experimental results while results from low amplitude forcing agree with linear stability theory.

Atkins, Harold L.↗

Airplane wing vibrations due to atmospheric turbulence

The magnitude of error introduced due to wing vibration when measuring atmospheric turbulence with a wind probe mounted at the wing tip was studied. It was also determined whether accelerometers mounted on the wing tip are needed to correct this error. A spectrum analysis approach is used to determine the error. Estimates of the B-57 wing characteristics are used to simulate the airplane wing, and von Karman's cross spectrum function is used to simulate atmospheric turbulence. It was found that wing vibration introduces large error in measured spectra of turbulence in the frequency's range close to the natural frequencies of the wing.

Pastel, R. L.↗

Implicit Large-Eddy Simulation of Axisymmetric and Asymmetric Turbulent Shock Wave Boundary Layer Interaction at Mach 2.5

High-speed inlet flows feature turbulent shock wave boundary layer interactions at supersonic Mach numbers. A canonical Mach 2.5 shock wave boundary layer interaction experiment at the NASA Glenn Research Center was designed for axisymmetric and asymmetric interactions depending on the placement of the shock generator. Four different wall-resolved implicit large-eddy simulations of the experiment were carried out at one-tenth of the Reynolds number of the experiment: (1) Two axisymmetric interactions with shock generator radii of 9.2% and 14.7% of the test section diameter and (2) two asymmetric interactions with 14.7% shock generator radius and shock generator centerline offsets of one-third and one-fourth of the test section radius. For the two axisymmetric interactions, various sub-domains upstream of the interaction, within the interaction, and downstream of the interaction were analyzed with the proper orthogonal decomposition. For all three regions, the dominant modes are streamwise coherent structures with Strouhal numbers based on separation length of 0.15 and below. Fourier spectra of the wall-pressure coefficient for the interaction region have peaks at Strouhal numbers based on separation length of 0.15 and 0.04. For the two asymmetric interactions, the streamwise extent of the separation varies considerably in the azimuthal direction. For both shock generator offsets, the interaction is strongest for the azimuthal location that is closest to the shock generator. Compared to the axisymmetric interaction, the streamwise extent of the separated flow region is noticeably reduced and a substantial crossflow is observed. Fourier analyses of the unsteady flow indicate low-frequency content f or the upstream part of the interaction which is more two-dimensional. Away from the upstream region, with increasing sweep angle and crossflow, the low-frequency content is diminished.

Computational Fluid Dynamics↗

Turbulent Mixing of Multiphase Flow

Thus we conduct numerical simulations of multiphase fluids stirred by two-dimensional turbulence to assess the possibility of self-similar drop size distribution in turbulence. In our turbulence simulations, we also explore the non-diffusive limit, where molecular mobility for the interface is vanishing. Special care is needed to transport the non-diffusive interface. Numerically, we use the particle level set method to evolve the interface. Instead of using the usual methods to calculate the surface tension force from the level set function, we reconstruct the interface based on phase- field modeling, and calculate the continuum surface tension forcing from the reconstructed interface.

Young, Y.-N.↗

Holocinematographic Velocimeter

Holographic movies used to study turbulent flows. High-speed movie composed of single-exposure holograms of tracer particles in each of two simultaneous orthogonal views of flow. Time resolution set by movie frame rate, and total time limited by length of movie. Data from system used to verify results of theoretical models and numerical simulations of turbulence, visualize and measure coherent structures in "non-simple" turbulent flows, to examine mechanisms operative in various turbulence-control/drag-reduction concepts.

Weinstein, Leonard M.↗