Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence methods”

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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Statistical Methodologies to Integrate Experimental and Computational Research

Development of advanced algorithms for simulating engine flow paths requires the integration of fundamental experiments with the validation of enhanced mathematical models. In this paper, we provide an overview of statistical methods to strategically and efficiently conduct experiments and computational model refinement. Moreover, the integration of experimental and computational research efforts is emphasized. With a statistical engineering perspective, scientific and engineering expertise is combined with statistical sciences to gain deeper insights into experimental phenomenon and code development performance; supporting the overall research objectives. The particular statistical methods discussed are design of experiments, response surface methodology, and uncertainty analysis and planning. Their application is illustrated with a coaxial free jet experiment and a turbulence model refinement investigation. Our goal is to provide an overview, focusing on concepts rather than practice, to demonstrate the benefits of using statistical methods in research and development, thereby encouraging their broader and more systematic application.

Parker, P. A.↗
DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. II. Kinematics

This second part of our research explores the kinematic aspect of fine-scale bubble-induced turbulence (BIT) to (i) present the effect of bubble breakup and coalescence and (ii) compare it against the universal kinematic fine-scale turbulence characteristics reported in the literature. To this end, we simulate a dilute bubbly system of 0.5% void fraction using two distinct numerical simulations. In the volume-of-fluid (VoF) simulation, bubbles undergo breakup and coalescence. In the immersed boundary method (IBM) simulation, however, they act as rigid spheres. We also perform a simulation of classical homogeneous isotropic turbulence (HIT). The first important outcome of this study is that BIT is radically different from HIT in terms of its kinematic fine-scale characteristics. In the vorticity-dominating regions, BIT exhibits a weak vortex stretching. This weak vortex stretching is due to (a) the intermediate strain-rate eigenvalues skewed weakly to positive and (b) the extensive strain-rate eigenvector aligning perpendicular to the vorticity vector. The BIT has, on average, not only a weak enstrophy production but also a weak strain production in strain-dominating regions. The weak strain production is due to (a) the presence of vortex stretching in highly strained fluid elements and (b) the absolute magnitude of compressive strain-rate eigenvalue being as close to the extensive strain-rate eigenvalue. Thus, none of the kinematic fine-scale HIT characteristics is noted for BIT. The second important conclusion is that bubble breakup and coalescence play little to no influence on the kinematics of fine-scale BIT as VoF and IBM simulations produce similar results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1988

Turbulence production near walls: The role of flow structures with spanwise asymmetry

Space-time evolution of near wall flow structures is described by conditional sampling methods, in which conditional averages are formed at various stages of development of shear layer structures. The development of spanwise asymmetry of the structures was found to be important in the creation of the structures and for the process of turbulence production.

Alfredsson, P. Henrik↗
DOE OSTI2020

Experimental determination of the thermal, turbulent, and rotational ion motion and magnetic field profiles in imploding plasmas

A tutorial is presented on advances in spectroscopic diagnostic methods developed for measuring key plasma properties in pulsed-power systems such as Z-pinches, magnetized-plasma compression devices, ion and electron diodes, and plasma switches. The parameters measured include the true ion temperature in Z-pinch implosions, which led to a discovery that much of the ion kinetic energy at stagnation is stored in hydrodynamic rather than in thermal motion. This observation contributed a new important insight into the understanding of the ion thermalization at stagnation and stimulated further investigations of turbulence at stagnation, discussed here too. The second part of this tutorial is devoted to the development of measurements for magnetic-field distributions in Z-pinches and in other pulsed-power systems, as well as their use in studying the plasma dynamics, resistivity, and pressure and energy balance. The latter study raises intriguing questions on the implosion process. In particular, in Z-pinches, the current during stagnation was found to largely flow at relatively large radii, outside the stagnation region. The magnetic-field measurements also enable investigations into the compression of a pre-magnetized cylindrical plasma that uncover striking phenomena related to the current flow, where the current was found to redistribute toward the outer regions during the implosion. Observation of the rotation of the magnetized plasma is also discussed. Lastly, experimental and theoretical investigations of a non-diffusive fast penetration of magnetic field into a low-density plasma, including its effect on the plasma dynamics, are described

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Structure of nonevaporating sprays. II - Drop and turbulence properties

This is the second part of a study reporting structure measurements in the dilute portion of axisymmetric nonevaporating sprays. Measurements are compared with predictions of three typical methods for analyzing sprays: (1) locally homogeneous flow (LHF) analysis, where slip between the phases is neglected; (2) deterministic separated flow (DSF) analysis, where slip is considered but effects of drop interactions with turbulence are ignored; and (3) stochastic separated flow (SSF) analysis, where both slip and effects of drop interactions with turbulence are considered. This part of the study reports measurements of mean and fluctuating drop velocities, the variation of Sauter mean diameter, and gas-phase turbulence properties in the dilute portion of the sprays. Best agreement between predictions and measurements was obtained with the SSF analysis. For present measurements in the dilute region (void fraction greater than 99.1 percent), effects of drops on gas-phase turbulence properties (turbulence modulation) were small. However, as the dense spray regions near the injector were approached, the measurements indicated modification of turbulence properties by drop motion.

Solomon, A. S. P.↗
NASA NTRS1992

Application of space-marching methods to hypersonic forebody flow fields

The heavy reliance on CFD in the design of hypersonic flight vehicles emphasizes the importance of choosing computationally efficient techniques. As a result, the development and application of space-marching flow solvers for the parabolized Navier-Stokes equations have received considerable attention in the past decade. This paper summarizes a space-marching approach used at NASA Ames Research Center for the analysis of hypersonic vehicle forebodies. Aspects of the simulation process are addressed that contain features unique to space-marching flow solvers, such as grid generation, initial condition specification, and finite-rate chemistry models. For completeness, sections are also included on turbulence modeling and equilibrium air chemistry modeling. The discussion includes an overview of the computational methods utilized as well as a number of relevant results.

Lawrence, Scott L.↗
NASA NTRS1975

Supersonic turbulent boundary-layer flows with mass injection through slots and/or porous walls

An implicit finite-difference method was used to solve the compressible boundary-layer equations, and to study the effects of mass transfer through porous plates, slots, and a combination of the two. The effects of the external pressure field were also included by using a global pressure interaction scheme. Two different eddy viscosity models were used for the slot and slot-porous combination cases: one was a two-layer model with inner and outer laws, and the other was a multi-layer model with as many as five separate layers. Results of the present method were compared with experimental data at a Mach number of 2.8. Comparisons of the skin friction reduction and Mach number profiles gave good to excellent agreement. Pressure interaction had little effect on the slot injection skin friction but increased the skin friction of the porous and slot-porous combination markedly.

Murray, A. L.↗
NASA NTRS1980

Drag reduction effects in turbulent boundary layers over wavy walls

Two dimensional incompressible flow over wavy surfaces are analyzed numerically by spectral methods. Algorithms for periodic flows (Fourier modes in the periodic flow direction and Chebycheff modes in the normal direction), and inflow-outflow boundary conditions (Chebycheff modes used in both directions) are described. Results obtained using both codes are reported for laminar flows. Comparisons with known theoretical and experimental results are made.

Balasubramanian, R.↗
NASA NTRS1993

Unstructured grid large eddy simulation of wall bounded turbulent flows

Historically, large eddy simulations (LES) have been restricted to simple geometries where spectral or finite difference methods have dominated due to their efficient use of structured grids. Structured grids, however, not only have difficulty representing complex domains and adapting to complicated flow features, but also are rather inefficient for simulating flows at high Reynolds numbers. The lack of efficiency stems from the need to resolve the viscous sub layer which requires very fine resolution in all three directions near the wall. Structured grids make use of a stretching to reduce the normal grid spacing but must carry the fine resolution in the streamwise and spanwise directions throughout the domain. The unnecessarily fine grid for much of the domain leads to disturbingly high grid estimates. Chapman (1979), and later Moin & Jimenez (1993), pointed out that, in order to advance the technology to airfoils at flight Reynolds numbers, structured grids must be abandoned in lieu of what are known as nested or unstructured grids. We illustrate the ability of an unstructured mesh to refine only the near wall region. Note the large number of points near the wall (where the fine vortical features need better resolution) and the coarseness in all directions away from the wall (where the scales are much larger). The important difference between this approach and the usual structured grid stretching is that the number of elements used to discretize the spanwise and streamwise features of the flow is reduced in each successive layer coming off the wall. This is due to the fact that the elements not only grow in the normal direction, but in the other directions as well. This greatly reduces the total number of points or elements required for a given Reynolds number flow.

Jansen, Kenneth↗
NASA NTRS1986

Advanced boundary layer transition measurement methods for flight applications

In modern laminar flow flight research, it is important to understand the specific cause(s) of laminar to turbulent boundary-layer transition. Such information is crucial to the exploration of the limits of practical application of laminar flow for drag reduction on aircraft. The transition modes of interest in current flight investigations include the viscous Tollmien-Schlichting instability, the inflectional instability at laminar separation, and the crossflow inflectional instability, as well as others. This paper presents the results to date of research on advanced devices and methods used for the study of laminar boundary-layer transition phenomena in the flight environment. Recent advancements in the development of arrayed hot-film devices and of a new flow visualization method are discussed. Arrayed hot-film devices have been designed to detect the presence of laminar separation, and of crossflow vorticity. The advanced flow visualization method utilizes color changes in liquid-crystal coatings to detect boundary-layer transition at high altitude flight conditions. Flight and wind tunnel data are presented to illustrate the design and operation of these advanced methods. These new research tools provide information on disturbance growth and transition mode which is essential to furthering our understanding of practical design limits for applications of laminar flow technology.

Holmes, B. J.↗
NASA NTRS1985

Simulation of turbulent mixing layers

A hybrid numerical scheme, combined with a scalar convection scheme, is used to compute turbulent mixing layers. The hybrid scheme is used to solve the vorticity equation by dividing the domain into two regions, applying a Lagrangian method in one region and an Eulerian method in the other. The scalar convection scheme is used to track an interface which separates a dyed fluid from a dye-free fluid, thereby distinguishing the fluid that starts on the high speed side from the fluid that starts on the low speed side. A two-dimensional mixing layer has been simulated using this new method; the roll up of a vortex sheet into structures that interact and pair is observed. The growth of the momentum thickness and the development of the mean turbulence statistics are found to be consistent with experimental observations. Forcing of the layer at moderate and low frequencies shows a dramatic effect on the growth of the momentum thickness and the entrainment ratio. Low frequency forcing enhances the growth of the layer; moderate frequency forcing inhibits the growth.

Mansour, N. N.↗
DOE OSTI2020

Elucidating plasma dynamics in Hasegawa–Wakatani turbulence by information geometry

In this work, the impact of adiabatic electrons on drift-wave turbulence, modeled by the Hasegawa–Wakatani equations, is studied using information length. Information length is a novel theoretical method for measuring distances between statistical states represented by different probability distribution functions (PDFs) along the path of a system and represents the total number of statistically different states that a system evolves through in time. Specifically, the time-dependent PDFs of turbulent fluctuations for a given adiabatic index $A$ are computed. The changes in fluctuation statistics are then quantified in time by using information length. The numerical results provide time traces exhibiting intermittent plasma dynamics, and such behavior is identified by a rapid change in the information length. The effects of $A$ are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2006

Increased Fidelity in Prediction Methods For Landing Gear Noise

An aeroacoustic prediction scheme has been developed for landing gear noise. The method is designed to handle the complex landing gear geometry of current and future aircraft. The gear is represented by a collection of subassemblies and simple components that are modeled using acoustic elements. These acoustic elements are generic, but generate noise representative of the physical components on a landing gear. The method sums the noise radiation from each component of the undercarriage in isolation accounting for interference with adjacent components through an estimate of the local upstream and downstream flows and turbulence intensities. The acoustic calculations are made in the code LGMAP, which computes the sound pressure levels at various observer locations. The method can calculate the noise from the undercarriage in isolation or installed on an aircraft for both main and nose landing gear. Comparisons with wind tunnel and flight data are used to initially calibrate the method, then it may be used to predict the noise of any landing gear. In this paper, noise predictions are compared with wind tunnel data for model landing gears of various scales and levels of fidelity, as well as with flight data on fullscale undercarriages. The present agreement between the calculations and measurements suggests the method has promise for future application in the prediction of airframe noise.

Lopes, Leonard V.↗
NASA NTRS1982

Comparison Between Experiment and Prediction for a Transonic Turbulent Separated Flow

Solutions of the time-dependent, mass-averaged Navier-Stokes equations are compared In detail with experimental results obtained on an axisymmetric "bump" model at a transonic Mach number that produced an extensive separated now region. In addition, an inverse boundary method is evaluated for this type of flow. The Cebeci-Smith algebraic and the Wilcox-Rubesin two-equation turbulence models used in the Navier-Stokes calculations both predict the maximum boundary-layer displacement thickness generated by the interaction reasonably well, with the details of the now best described with the two-equation formulation. However, both models predict a shock location substantially farther aft on the bump than observed experimentally. This error in shock location was slightly less with the two-equation model (0.12 chord compared with 0.16 chord). In the vicinity of the shock, the calculations predict a more rapid increase in turbulent shear stress than observed in the experimental results; this more rapid increase is believed to be the cause or the poor predictions in shock position.

Johnson, D. A.↗
NASA NTRS1992

Further wind tunnel investigation of the SM701 airfoil with aileron and turbulators

Wind tunnel tests were performed on a two-dimensional model of the SM701 airfoil designed for use on the World Class gliders. The test covered a range of Reynolds numbers from 500,000 to 1.7 million. Aerodynamic forces and moments were measured with an external balance. Momentum loss method measurements of the section drag coefficient were also made. Flow visualization techniques provided information on transition from laminar to turbulent flow. Lift, drag, and pitching moment were analyzed and comparisons were made with predicted and previously obtained experimental data. The effects of V-tape turbulators for use in turbulent drag reduction were studied. The performance of a 25 percent chord aileron deflected through plus or minus 20 degrees was researched. The model was designed, constructed, and tested by students at Texas A&M University.

Steen, Gregory↗
NASA NTRS1997

Direct Numerical Simulation of a Temporally Evolving Incompressible Plane Wake: Effect of Initial Conditions on Evolution and Topology

Direct numerical simulations have been used to examine the effect of the initial disturbance field on the development of three-dimensionality and the transition to turbulence in the incompressible plane wake. The simulations were performed using a new numerical method for solving the time-dependent, three-dimensional, incompressible Navier-Stokes equations in flows with one infinite and two periodic directions. The method uses standard Fast Fourier Transforms and is applicable to cases where the vorticity field is compact in the infinite direction. Initial disturbances fields examined were combinations of two-dimensional waves and symmetric pairs of 60 deg oblique waves at the fundamental, subharmonic, and sub-subharmonic wavelengths. The results of these simulations indicate that the presence of 60 deg disturbances at the subharmonic streamwise wavelength results in the development of strong coherent three-dimensional structures. The resulting strong three-dimensional rate-of-strain triggers the growth of intense fine scale motions. Wakes initiated with 60 deg disturbances at the fundamental streamwise wavelength develop weak coherent streamwise structures, and do not develop significant fine scale motions, even at high Reynolds numbers. The wakes which develop strong three-dimensional structures exhibit growth rates on par with experimentally observed turbulent plane wakes. Wakes which develop only weak three-dimensional structures exhibit significantly lower late time growth rates. Preliminary studies of wakes initiated with an oblique fundamental and a two-dimensional subharmonic, which develop asymmetric coherent oblique structures at the subharmonic wavelength, indicate that significant fine scale motions only develop if the resulting oblique structures are above an angle of approximately 45 deg.

Sondergaard, R.↗
NASA NTRS2016

Direct Numerical Simulation of an Airfoil with Sand Grain Roughness on the Leading Edge

As part of a computational study of acoustic radiation due to the passage of turbulent boundary layer eddies over the trailing edge of an airfoil, the Lattice-Boltzmann method is used to perform direct numerical simulations of compressible, low Mach number flow past an NACA 0012 airfoil at zero degrees angle of attack. The chord Reynolds number of approximately 0.657 million models one of the test conditions from a previous experiment by Brooks, Pope, and Marcolini at NASA Langley Research Center. A unique feature of these simulations involves direct modeling of the sand grain roughness on the leading edge, which was used in the abovementioned experiment to trip the boundary layer to fully turbulent flow. This report documents the findings of preliminary, proof-of-concept simulations based on a narrow spanwise domain and a limited time interval. The inclusion of fully-resolved leading edge roughness in this simulation leads to significantly earlier transition than that in the absence of any roughness. The simulation data is used in conjunction with both the Ffowcs Williams-Hawkings acoustic analogy and a semi-analytical model by Roger and Moreau to predict the farfield noise. The encouraging agreement between the computed noise spectrum and that measured in the experiment indicates the potential payoff from a full-fledged numerical investigation based on the current approach. Analysis of the computed data is used to identify the required improvements to the preliminary simulations described herein.

Ribeiro, Andre F. P.↗
DOE OSTI2023

Temporal Subsampling Diminishes Small Spatial Scales in Recurrent Neural Network Emulators of Geophysical Turbulence

The immense computational cost of traditional numerical weather and climate models has sparked the development of machine learning (ML) based emulators. Because ML methods benefit from long records of training data, it is common to use data sets that are temporally subsampled relative to the time steps required for the numerical integration of differential equations. Here, we investigate how this often overlooked processing step affects the quality of an emulator's predictions. We implement two ML architectures from a class of methods called reservoir computing: (a) a form of Nonlinear Vector Autoregression (NVAR), and (b) an Echo State Network (ESN). Despite their simplicity, it is well documented that these architectures excel at predicting low dimensional chaotic dynamics. We are therefore motivated to test these architectures in an idealized setting of predicting high dimensional geophysical turbulence as represented by Surface Quasi-Geostrophic dynamics. In all cases, subsampling the training data consistently leads to an increased bias at small spatial scales that resembles numerical diffusion. Interestingly, the NVAR architecture becomes unstable when the temporal resolution is increased, indicating that the polynomial based interactions are insufficient at capturing the detailed nonlinearities of the turbulent flow. The ESN architecture is found to be more robust, suggesting a benefit to the more expensive but more general structure. Spectral errors are reduced by including a penalty on the kinetic energy density spectrum during training, although the subsampling related errors persist. Future work is warranted to understand how the temporal resolution of training data affects other ML architectures.

58 GEOSCIENCES↗