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,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Viscous Fluxes - Part 1

Discretization of the viscous terms in current finite-volume unstructured-grid schemes are compared using node-centered and cell-centered approaches in two dimensions. Accuracy and efficiency are studied for six nominally second-order accurate schemes: a node-centered scheme, cell-centered node-averaging schemes with and without clipping, and cell-centered schemes with unweighted, weighted, and approximately mapped least-square face gradient reconstruction. The grids considered range from structured (regular) grids to irregular grids composed of arbitrary mixtures of triangles and quadrilaterals, including random perturbations of the grid points to bring out the worst possible behavior of the solution. Two classes of tests are considered. The first class of tests involves smooth manufactured solutions on both isotropic and highly anisotropic grids with discontinuous metrics, typical of those encountered in grid adaptation. The second class concerns solutions and grids varying strongly anisotropically over a curved body, typical of those encountered in high-Reynolds number turbulent flow simulations. Results from the first class indicate the face least-square methods, the node-averaging method without clipping, and the node-centered method demonstrate second-order convergence of discretization errors with very similar accuracies per degree of freedom. The second class of tests are more discriminating. The node-centered scheme is always second order with an accuracy and complexity in linearization comparable to the best of the cell-centered schemes. In comparison, the cell-centered node-averaging schemes are less accurate, have a higher complexity in linearization, and can fail to converge to the exact solution when clipping of the node-averaged values is used. The cell-centered schemes using least-square face gradient reconstruction have more compact stencils with a complexity similar to the complexity of the node-centered scheme. For simulations on highly anisotropic curved grids, the least-square methods have to be amended either by introducing a local mapping of the surface anisotropy or modifying the scheme stencil to reflect the direction of strong coupling.

Diskin, Boris↗
NASA NTRS2010

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Viscous Fluxes

Discretization of the viscous terms in current finite-volume unstructured-grid schemes are compared using node-centered and cell-centered approaches in two dimensions. Accuracy and complexity are studied for four nominally second-order accurate schemes: a node-centered scheme and three cell-centered schemes - a node-averaging scheme and two schemes with nearest-neighbor and adaptive compact stencils for least-square face gradient reconstruction. The grids considered range from structured (regular) grids to irregular grids composed of arbitrary mixtures of triangles and quadrilaterals, including random perturbations of the grid points to bring out the worst possible behavior of the solution. Two classes of tests are considered. The first class of tests involves smooth manufactured solutions on both isotropic and highly anisotropic grids with discontinuous metrics, typical of those encountered in grid adaptation. The second class concerns solutions and grids varying strongly anisotropically over a curved body, typical of those encountered in high-Reynolds number turbulent flow simulations. Tests from the first class indicate the face least-square methods, the node-averaging method without clipping, and the node-centered method demonstrate second-order convergence of discretization errors with very similar accuracies per degree of freedom. The tests of the second class are more discriminating. The node-centered scheme is always second order with an accuracy and complexity in linearization comparable to the best of the cell-centered schemes. In comparison, the cell-centered node-averaging schemes may degenerate on mixed grids, have a higher complexity in linearization, and can fail to converge to the exact solution when clipping of the node-averaged values is used. The cell-centered schemes using least-square face gradient reconstruction have more compact stencils with a complexity similar to that of the node-centered scheme. For simulations on highly anisotropic curved grids, the least-square methods have to be amended either by introducing a local mapping based on a distance function commonly available in practical schemes or modifying the scheme stencil to reflect the direction of strong coupling. The major conclusion is that accuracies of the node centered and the best cell-centered schemes are comparable at equivalent number of degrees of freedom.

Diskin, Boris↗
NASA NTRS2012

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗
NASA NTRS2010

Direct Measurement Sensor of the Boundary Shear Stress in Fluid Flow

The flow fields and boundary erosion that are associated with scour at bridge piers are very complex. Direct measurement of the boundary shear stress and boundary pressure fluctuations in experimental scour research has always been a challenge and high spatial resolution and fidelity have been almost impossible. Most researchers have applied an indirect process to determine shear stress using precise measured velocity profiles. Laser Doppler Anemometry and Particle Image Velocimetry are common techniques used to accurately measure velocity profiles. These methods are based on theoretical assumptions to estimate boundary shear stress. In addition, available turbulence models cannot very well account for the effect of bed roughness which is fundamentally important for any CFD simulation. The authors have taken on the challenge to advance the magnitude level to which direct measurements of the shear stress in water flow can be performed. This paper covered the challenges and the efforts to develop a higher accuracy and small spatial resolution sensor. Also, preliminary sensor designs and test results are presented.

Badescu, Mircea↗
NASA NTRSDate not supplied

Streamline-Traced, External-Compression Supersonic Inlets for Mach 2

A computational study was performed to explore the aerodynamic design and performance of streamline-traced, external-compression (STEX) inlets for Mach 2.0. The performance metrics included inlet flow rates, total pressure recovery, and total pressure distortion. The study explored the use of round and flattened capture cross-sections and a bleed slot. The design of the inlet and generation of the inlet geometry was performed using the Supersonic Inlet Design and Analysis(SUPIN) Tool. Computational grids were generated, and methods of computational fluid dynamics (CFD) were applied to solve the three-dimensional, turbulent flow through the inlets using the Wind-US CFD flow solver. It was found that the use of a bleed slot with about 5% bleed resulted in an inlet total pressure recovery of 95% at Mach 2 with acceptable radial and circumferential total pressure distortion. It was also found that there was only a slight decrease in performance between a round and flattened capture cross-section. This suggests that the use of streamline tracing offers flexibility in shaping the inlet capture cross-section for a more favorable integration with an aircraft wing or fuselage. This flexibility and the good performance of the inlet provides a novel inlet design for future commercial supersonic aircraft.

Computational Fluid Dynamics↗
NASA NTRS1996

Biases In A Magnitude Limited Versus A Distant Limited Planet Search

All methods utilized to conduct a search or survey inevitably have some built-in biases. These biases are often traced to some limitation of the instrument used or some inherent character of the signal being recorded. We address these limitations for various methods used or proposed for planet detection: spectroscopy, astrometry, interferometry, and photometry. For spectroscopy, the turbulence in the photosphere limits the minimum measurable dossier velocity to 3 m/s and hence the minimum planet to star mass ratio, thereby favoring massive close-in planets. Limited available observing time will necessarily introduce additional selection biases in the targets observed and telescope aperture will limit the faintest magnitude stars to be measured; For astrometry, the angular resolution of the instrument along with motions in the photometric center of the star limit the furthest distance for which giant planets can be detected to about 10 pc and favor massive outer orbit planets around low mass non-solar like stars; For imaging interferometry, the minimum angular size of the central null limits both the distance to the star and closeness of the planet to the host star to about 1 AU at 10 pc and thus to the very few solar-like stars within 10 pc. Solar and extra-solar zodiacal emission will limit the minimum size of the detectable planet. For photometry, the inherent variability of the star does not limit the minimum planet size until earth-sized or smaller planets are considered around solar-like stars. The telescope aperture limits the faintest stars that can be monitored. As with spectroscopy, there is no inherent distant limit to the method. After addressing the limiting factors of each method, an estimate is made of the number of planets of various sizes that could be found for each stellar type based on the detection probability and the number of stars that can be searched for planets.

Koch, David↗
NASA NTRS2001

Turbulence Modeling Workshop

Advances in turbulence modeling are needed in order to calculate high Reynolds number flows near the onset of separation and beyond. To this end, the participants in this workshop made the following recommendations. (1) A national/international database and standards for turbulence modeling assessment should be established. Existing experimental data sets should be reviewed and categorized. Advantage should be taken of other efforts already under-way, such as that of the European Research Community on Flow, Turbulence, and Combustion (ERCOFTAC) consortium. Carefully selected "unit" experiments will be needed, as well as advances in instrumentation, to fill the gaps in existing datasets. A high priority should be given to document existing turbulence model capabilities in a standard form, including numerical implementation issues such as grid quality and resolution. (2) NASA should support long-term research on Algebraic Stress Models and Reynolds Stress Models. The emphasis should be placed on improving the length-scale equation, since it is the least understood and is a key component of two-equation and higher models. Second priority should be given to the development of improved near-wall models. Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) would provide valuable guidance in developing and validating new Reynolds-averaged Navier-Stokes (RANS) models. Although not the focus of this workshop, DNS, LES, and hybrid methods currently represent viable approaches for analysis on a limited basis. Therefore, although computer limitations require the use of RANS methods for realistic configurations at high Reynolds number in the foreseeable future, a balanced effort in turbulence modeling development, validation, and implementation should include these approaches as well.

Rubinstein, R.↗
NASA NTRS1983

Helicopter blade tips

Methods of improving helicopter performance and vibration level by proper shaping of helicopter blade tips are considered. The principle involved consists of reducing the extent of the supersonic zone above the advancing tip and of the turbulent interaction. For stationary and advancing flight, the influence of the rotor and the problems posed by blade tips are reviewed. The theoretical methods of dealing with the two types of flight are briefly stated, and the experimental apparatus is described, including model triple and quadruple rotors. Different blade tip shapes are shown and briefly discussed. The theoretical results include an advancing speed of 309 km/H and a blade tip rotational speed of 215 m/s. The experimental values are advancing speed of 302 km/h and blade tip Mach number 0.86 for both types of rotors.

Lyothier, R.↗
NASA NTRS1977

A numerical study of the unsteady leading edge separation bubble on an oscillating airfoil

A preliminary study of the unsteady viscous flow in the region of an airfoil leading edge was performed, in which the interaction between the viscous and inviscid flow fields is neglected. The solution method uses the finite difference form of the governing equations throughout the separated flow field and incorporates a transition model based on the integral turbulence kinetic energy equation. The validity of the numerical procedure is verified by making comparisons with analytical solutions to several test problems, including unsteady flow over a plate oscillating in its own plane. The method was then applied to the problem of unsteady viscous flow over a NACA 0012 airfoil oscillating sinusoidally in pitch. The flow field characteristics were in qualitative agreement with experimental results. The bubble moved forward on the airfoil and decreased in size as incidence was increased. Viscous flow in the leading edge region was found to be quasi-steady, while bubble height varied inversely with Reynolds number.

Kreskovsky, J. P.↗
NASA NTRS1992

Reactant conversion in homogeneous turbulence: Mathematical modeling, computational validations and practical applications

Closed form analytical expressions are obtained for predicting the limited rate of reactant conversion in a binary reaction of the type F + rO yields (1 + r) Product in unpremixed homogeneous turbulence. These relations are obtained by means of a single point Probability Density Function (PDF) method based on the Amplitude Mapping Closure. It is demonstrated that with this model, the maximum rate of the reactants' decay can be conveniently expressed in terms of definite integrals of the Parabolic Cylinder Functions. For the cases with complete initial segregation, it is shown that the results agree very closely with those predicted by employing a Beta density of the first kind for an appropriately defined Shvab-Zeldovich scalar variable. With this assumption, the final results can also be expressed in terms of closed form analytical expressions which are based on the Incomplete Beta Functions. With both models, the dependence of the results on the stoichiometric coefficient and the equivalence ratio can be expressed in an explicit manner. For a stoichiometric mixture, the analytical results simplify significantly. In the mapping closure, these results are expressed in terms of simple trigonometric functions. For the Beta density model, they are in the form of Gamma Functions. In all the cases considered, the results are shown to agree well with data generated by Direct Numerical Simulations (DNS). Due to the simplicity of these expressions and because of nice mathematical features of the Parabolic Cylinder and the Incomplete Beta Functions, these models are recommended for estimating the limiting rate of reactant conversion in homogeneous reacting flows. These results also provide useful insights in assessing the extent of validity of turbulence closures in the modeling of unpremixed reacting flows. Some discussions are provided on the extension of the model for treating more complicated reacting systems including realistic kinetics schemes and multi-scalar mixing with finite rate chemical reactions in more complex configurations.

Madnia, C. K.↗
NASA NTRS1996

Multiple Scale Simulation for Transitional and Turbulent Flow

A new concept, Multiple Scale Simulation (MSS), is presented in this paper. The basic idea is that the flow is decomposed into several component groups according to spatial and temporal length scales. Each group has its own subdomain, governing system, mesh size, and discretization method. The simulation is then performed groupwise. This approach has been successfully applied in combination with the intergrid dissipation technique for simulation of transitional and turbulent flow in 3-D boundary layers, and it is feasible for 3-D airfoils and other more complex configurations. MSS should prove to ameliorate the scale problems associated with conventional direct numerical simulation.

Liu, Chaoqun↗
NASA NTRS1971

Gust Alleviation

A discussion is presented of the response of an airplane to turbulence, including techniques for the alleviation of lateral and longitudinal motions. A simplified theoretical treatment of the longitudinal response to vertical gusts, which allows calculation of the effects of various types of gust-alleviation systems, is presented. Flight experiments on various gust-alleviation systems are reviewed. A brief discussion is given of the application of these systems to different types of airplanes. Finally, methods of refinement of the simplified theory to include effects of horizontal gusts, unsteady lift effects, and isotropic turbulence are described.

William H Phillips↗
NASA NTRS1946

On laminar and turbulent friction

Report deals, first with the theory of the laminar friction flow, where the basic concepts of Prandtl's boundary layer theory are represented from mathematical and physical points of view, and a method is indicated by means of which even more complicated cases can be treated with simple mathematical means, at least approximately. An attempt is also made to secure a basis for the computation of the turbulent friction by means of formulas through which the empirical laws of the turbulent pipe resistance can be applied to other problems on friction drag. (author)

Von Karman, TH↗
NASA NTRS1977

Response of periodic beam to supersonic boundary-layer pressure fluctuations

The response of a periodic beam (modeling a periodic fuselage) to supersonic boundary-layer pressure fluctuations is analyzed on the basis of a scheme in which a decaying turbulence is treated as a superposition of frozen-pattern components, thus allowing the structural response to be similarly superposed and the advantage of frozen-pattern analysis to be maximally utilized. The fundamental solution required for the construction of the total response is one corresponding to the excitation of a frozen-pattern sinusoid. To obtain this fundamental solution, the formulation follows Mead's wave-propagation method (1971), but also takes into account the effect of freestream velocity on the same side of the turbulence excitation and the effect of a cavity on the opposite side of the excitation. As a numerical example, the spectral density of the structural response is computed and the results are compared with experimental data.

Lin, Y. K.↗
NASA NTRS1991

Simulating boundary layer transition with low-Reynolds-number k-epsilon turbulence models. I - An evaluation of prediction characteristics. II - An approach to improving the predictions

The capability of two k-epsilon low-Reynolds number (LRN) turbulence models, those of Jones and Launder (1972) and Lam and Bremhorst (1981), to predict transition in external boundary-layer flows subject to free-stream turbulence is analyzed. Both models correctly predict the basic qualitative aspects of boundary-layer transition with free stream turbulence, but for calculations started at low values of certain defined Reynolds numbers, the transition is generally predicted at unrealistically early locations. Also, the methods predict transition lengths significantly shorter than those found experimentally. An approach to overcoming these deficiencies without abandoning the basic LRN k-epsilon framework is developed. This approach limits the production term in the turbulent kinetic energy equation and is based on a simple stability criterion. It is correlated to the free-stream turbulence value. The modification is shown to improve the qualitative and quantitative characteristics of the transition predictions.

Schmidt, R. C.↗
NASA NTRS1980

Skin friction measurements by a new nonintrusive double-laser-beam oil viscosity balance technique

A portable dual-laser-beam interferometer that nonintrusively measures skin friction by monitoring the thickness change of an oil film subject to shear stress is described. The method is an advance over past versions in that the troublesome and error-introducing need to measure the distance to the oil leading edge and the starting time for the oil flow has been eliminated. The validity of the method was verified by measuring oil viscosity in the laboratory, and then using those results to measure skin friction beneath the turbulent boundary layer in a low-speed wind tunnel. The dual-laser-beam skin friction measurements are compared with Preston tube measurements, with mean velocity profile data in a 'law-of-the-wall' coordinate system, and with computations based on turbulent boundary-layer theory. Excellent agreement is found in all cases. This validation and the aforementioned improvements appear to make the present form of the instrument usable to measure skin friction reliably and nonintrusively in a wide range of flow situations in which previous methods are not practical.

Monson, D. J.↗
NASA NTRS1996

Effects of Gravity on Sheared Turbulence Laden with Bubbles or Droplets

This is a new project which started in May 1996. The main objective of the experimental/numerical study is to improve the understanding of the physics of two-way coupling between the dispersed phase and turbulence in a prototypical turbulent shear flow - homogeneous shear, laden with small liquid droplets (in gas) or gaseous bubbles (in liquid). The method of direct numerical simulation (DNS) is used to solve the full three-dimensional, time-dependent Navier-Stokes equations including the terms describing the two-way coupling between the dispersed phase and the carrier flow. The results include the temporal evolution of the three-dimensional energy and dissipation spectra and the rate of energy transfer across the energy spectrum to understand the fundamental physics of turbulence modulation, especially the effects of varying the magnitude of gravitational acceleration. The mean-square displacement and diffusivity of the droplets (or bubbles) of a given size and the preferential accumulation of droplets in low vorticity regions and bubbles in high vorticity regions will be examined in detail for different magnitudes of gravitational acceleration. These numerical results which will be compared with their corresponding measured data will provide a data base from which a subgrid-scale (SGS) model can be developed and validated for use in large-eddy simulation (LES) of particle-laden shear flows. Two parallel sets of experiments will be conducted: bubbles in an immiscible liquid and droplets in air. In both experiments homogeneous shear will be imposed on the turbulent carrier flow. The instantaneous velocities of the fluid and polydispersed-size particles (droplets or bubbles) will be measured simultaneously using a two-component Phase-Doppler Particle Analyzer (PDPA). Also, the velocity statistics and energy spectra for the carrier flow will be measured.

Elghobashi, Said↗
NASA NTRS1981

Skin Friction Measurements by a Dual-Laser-Beam Interferometer Technique

A portable dual-laser-beam interferometer that nonintrusively measures skin friction by monitoring the thickness change of an oil film subject to shear stress is described. The method is an advance over past versions in that the troublesome and error-introducing need to measure the distance to the oil leading edge and the starting time for the oil flow has been eliminated. The validity of the method was verified by measuring oil viscosity in the laboratory, and then using those results to measure skin friction beneath the turbulent boundary layer in a low speed wind tunnel. The dual-laser-beam skin friction measurements are compared with Preston tube measurements, with mean velocity profile data in a "law-of-the-well" coordinate system, and with computations based on turbulent boundary-layer theory. Excellent agreement is found in all cases. (This validation and the aforementioned improvements appear to make the present form of the instrument usable to measure skin friction reliably and nonintrusively in a wide range of flow situations in which previous methods are not practical.)

Monson, D. J.↗