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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer

The objective of this work is to compare a high-order solver with a low-order solver for performing Large-Eddy Simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the highorder method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.

large eddy simulations↗
NASA NTRS2015

Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer

The objective of this work is to compare a high-order solver with a low-order solver for performing large-eddy simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the high-order method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.

large eddy simulations↗
DOE OSTI2020

Turbulence model reduction by deep learning

A defining problem of turbulence theory is to produce a predictive model for turbulent fluxes. These have profound implications for virtually all aspects of the turbulence dynamics. In magnetic confinement devices, drift-wave turbulence produces anomalous fluxes via cross-correlations between fluctuations. In this work, we introduce an alternative, data-driven method for parametrizing these fluxes. The method uses deep supervised learning to infer a reduced mean-field model from a set of numerical simulations. We apply the method to a simple drift-wave turbulence system and find a significant new effect which couples the particle flux to the local gradient of vorticity. Notably, here, this effect is much stronger than the oft-invoked shear suppression effect. We also recover the result via a simple calculation. The vorticity gradient effect tends to modulate the density profile. In addition, our method recovers a model for spontaneous zonal flow generation by negative viscosity, stabilized by nonlinear and hyperviscous terms. We highlight the important role of symmetry to implementation of the new method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Transonic airfoil calculations including wind tunnel wall-interference effects

The results of Reynolds-averaged time-dependent inviscid and turbulent compressible Navier-Stokes computations using the implicit finite-difference approach of Steger (1978), modified by incorporating a pressure boundary condition, (PBC) to account for wall interference are compared with experimental data on a NACA 64A010 airfoil (Johnson and Bachalo, 1980) in graphs and briefly characterized. The computational approach is the same as that used by King and Johnson (1980), but a 137 x 50 mesh is used instead of a 97 x 35 mesh, and special care is taken in resolving the nose, shock, and trailing-edge regions. Imposition of PBC is shown to improve significantly the accuracy of the computations for the flowfield on the upper surface of the airfoil, shifting the shock forward to its experimentally measured position in the case of turbulent flow. The failure of the method, even with PBC, to match the experimental shock location in the case of a flow with a separation bubble is attributed to inadequacies in the algebraic turbulence model employed (Baldwin and Lomax, 1978).

King, L. S.↗
NASA NTRS1978

Critical study of higher order numerical methods for solving the boundary-layer equations

A fourth order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method, which is the natural extension of the second order box scheme to fourth order, was demonstrated with application to the incompressible, laminar and turbulent, boundary layer equations. The efficiency of the present method is compared with two point and three point higher order methods, namely, the Keller box scheme with Richardson extrapolation, the method of deferred corrections, a three point spline method, and a modified finite element method. For equivalent accuracy, numerical results show the present method to be more efficient than higher order methods for both laminar and turbulent flows.

Wornom, S. F.↗
NASA NTRS1978

Theoretical and experimental study of the drag of multielement airfoils

The viscous/potential flow past single-element and multielement airfoils is studied theoretically and experimentally. A computerized analysis, based on iteratively coupled potential-flow and boundary-layer analysis, is used to predict the flow field of the airfoil. The method yields detailed characteristics of conventional laminar and turbulent boundary layers, turbulent wakes, and confluent boundary layers. The viscous flows are analyzed with a method that uses finite-difference solutions of the boundary-layer equations. Reynolds stress in the boundary layers and wakes is simulated with eddy viscosity models for the various flow zones. The viscous calculations are carried into the wake of the airfoil where the drag is found from the defect in the wake momentum.

Olson, L. E.↗
NASA NTRS1985

Simultaneous measurements of temperature, density, and pressure in a supersonic turbulent flow using laser-induced fluorescence

A pulsed laser-induced fluorescence technique is described that provides simultaneous measurements of temperature, density, and pressure in low-temperature, turbulent flows. The measurements are made with spatial and temporal resolution comparable to that obtained with modern laser anemometer techniques used for turbulent boundary layer research. The capabilities of the method are briefly described and its demonstration in a simple two-dimensional turbulent boundary layer at Mach 2 is reported. The results are compared with conventional hot-wire anemometer data obtained in the same flow.

Mckenzie, R. L.↗
NASA NTRS1988

A multiple-time-scale turbulence model based on variable partitioning of turbulent kinetic energy spectrum

The paper presents a multiple-time-scale turbulence model of a single point closure and a simplified split-spectrum method. Consideration is given to a class of turbulent boundary layer flows and of separated and/or swirling elliptic turbulent flows. For the separated and/or swirling turbulent flows, the present turbulence model yielded significantly improved computational results over those obtained with the standard k-epsilon turbulence model.

Kim, S.-W.↗
NASA NTRS1987

Relations between two-point correlations and pressure strain terms

The structure of the two-point spatial correlations (velocity-velocity, velocity-scalar, and scalar-scalar) were studied with a view to improve turbulence closure models. The linear model for the two-point correlations proposed by Naot provides a method of including the information about the turbulence structure in the turbulence models. The assumptions and adequacy of this model were tested against the homogeneous shear flow simulation data base. The model performs poorly in some details and it is suggested how it may be improved. The models were also tested for rapid pressure-strain terms in a variety of flows including axisymmetric expansion and contraction flows, homogeneous shear flow, channel flow, and boundary layer.

Wolfshtein, Micha↗
NASA NTRS1995

Assessment of climate variability of the Greenland ice sheet: Integration of in situ and satellite data

The proposed research involves the application of multispectral satellite data in combination with ground truth measurements to monitor surface properties of the Greenland ice sheet which are essential for describing the energy and mass of the ice sheet. Several key components of the energy balance are parameterized using satellite data and in situ measurements. The analysis will be done for a ten year time period in order to get statistics on the seasonal and interannual variations of the surface processes and the climatology. Our goal is to investigate to what accuracy and over what geographic areas large scale snow properties and radiative fluxes can be derived based upon a combination of available remote sensing and meterological data sets. Data analysis showed the following results: (1)cloud classification based on longwave sky radiation revealed that overcast sky occurred for 25% of the time in winter, and for 15% in spring and summer respectively (winter and summer both show the same occurrence of clear sky of approximately 26%); (2) comparison of aerodynamic profile method with eddy correlation method to derive sensible and latent heat flux showed good agreement in the diurnal cycle and the turbulent fluxes were underestimated with the aerodynamic method by 10 - 30% as compared to the in situ eddy flux method; (3) the katabatic wind shows a distinct diurnal cycle with a maximum in the morning (7-9 h solar time) and a minimum in the later afternoon (18 h solar time); (4) snow grain size was modeled with a surface energy balance model (SNTHERM) and compared with in situ measurements. Sharp decreases in the modeled snow grain size, caused by accumulation events such as precipitation and deposition, could be verified with observational data; (4) radiative transfer modeling of firn supports our beliefs that the observed trends in 18 and 19 GHz passive microwave brightness temperatures are attributable to accumulation rate changes (modeling also indicates the above relationship is detectable because of the presence of depth hoar; (5) snow melt can be detected by a distinct signal in the passive microwave cross-polarized gradient ratio (19h-37v)/(19h+37v) and has been used for wet/dry snow classification; (6) top of the atmosphere (TOA) broadband albedos were derived from AVHRR visible and near infrared reflectances for the entire ice sheet from May 1990 - June 1991, and the highest albedo values are found along the southeast coast of the ice sheet which is consistent with the summer peak of precipitation due to onshore flow loaded with high water vapor content (TOA albedo values dropped to around 40% along the south-western coast during July and August due to bare ice surfaces); and (7) the net all-wave radiation balance at the top of the atmosphere is negative over the entire ice sheet except for the summer months - June-July-August (in June, the net radiation balance is slightly positive over the dry snow areas--15 W/m2).

Steffen, K.↗
NASA NTRS1978

A method for calculating strut and splitter plate noise in exit ducts: Theory and verification

Portions of a four-year analytical and experimental investigation relative to noise radiation from engine internal components in turbulent flow are summarized. Spectra measured for such airfoils over a range of chord, thickness ratio, flow velocity, and turbulence level were compared with predictions made by an available rigorous thin-airfoil analytical method. This analysis included the effects of flow compressibility and source noncompactness. Generally good agreement was obtained. This noise calculation method for isolated airfoils in turbulent flow was combined with a method for calculating transmission of sound through a subsonic exit duct and with an empirical far-field directivity shape. These three elements were checked separately and were individually shown to give close agreement with data. This combination provides a method for predicting engine internally generated aft-radiated noise from radial struts and stators, and annular splitter rings. Calculated sound power spectra, directivity, and acoustic pressure spectra were compared with the best available data. These data were for noise caused by a fan exit duct annular splitter ring, larger-chord stator blades, and turbine exit struts.

Fink, M. R.↗
NASA NTRS1991

Computational Fluid Dynamics Symposium on Aeropropulsion

Recognizing the considerable advances that have been made in computational fluid dynamics, the Internal Fluid Mechanics Division of NASA Lewis Research Center sponsored this symposium with the objective of providing a forum for exchanging information regarding recent developments in numerical methods, physical and chemical modeling, and applications. This conference publication is a compilation of 4 invited and 34 contributed papers presented in six sessions: algorithms one and two, turbomachinery, turbulence, components application, and combustors. Topics include numerical methods, grid generation, chemically reacting flows, turbulence modeling, inlets, nozzles, and unsteady flows.

Source record↗
DOE OSTI2023

Large Eddy Simulation of Gasoline Sprays in a Lagrangian–Eulerian Framework Using the High-Order Spectral Element Method

Predicting the spray evolution using simulations requires accurate modeling of the turbulent gas-phase flow field. Here, in this study, the high-order spectral-element method (SEM), implemented in the code Nek5000, was used to provide highly resolved solutions to the turbulent flow field. Spray modeling capabilities were implemented into the Nek5000 code. The spray is modeled in a Lagrangian–Eulerian (LE) framework, where the liquid is represented by discrete parcels of droplets. The method for coupling liquid and gas in the context of SEM is described, which allows for very fine meshes to be used without affecting the stability of the solution. Large-eddy simulations (LES) of the eight-hole ECN Spray G gasoline injector were conducted. Numerical results are compared against experimental data for liquid penetration, droplet size and gas velocity. The morphology of the multiplume spray is compared against experimental data. The effect of different spray injection inputs is analyzed. It was found that using a plume direction of 33 deg and an injection cone angle of 30 deg produced the best results overall. This work shows the applicability of SEM for spray modeling applications, where use of a high-order flow solver can help us understand the multiplume spray aerodynamics and how it leads to plume collapse under certain conditions. Results also highlight the need for tuning spray input parameters in the LE framework, even when high-fidelity gas flow solutions are possible.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1993

Surface hot-film method for the measurement of transition, separation and reattachment points

A real-time method of determining positions of laminar-to-turbulent transition region, separation and reattachment points and stagnation points using an array of simultaneously operated surface-mounted hot-film sensors has been developed and applied to a wind-tunnel test of a multielement airfoil model. Determination of various types of transitions and flow directions in various regimes of flows seen on multielement airfoils are possible without precise sensor calibration or laborious post-test data analysis. The results agree with established method and theoretical methods, but determination of turbulent reattachment points are not yet satisfactory.

Nakayama, Akihiko↗
NASA NTRS1990

Computation and turbulence modeling for three-dimensional boundary layers including turbomachinery rotor flows

A method is developed for predicting the behavior of three-dimensional, turbulent boundary layers occurring in internal flows, including those on turbomachinery rotor blades. These boundary layers are complex, turbulent, and subject to Coriolis and centrifugal forces. The major thrust of this paper is the development and use of an algebraic Reynolds stress model (ARSM) that captures the changes in turbulent flow structure arising from curvature, rotation, and three dimensionality. The prediction of pressure-driven secondary flow agrees well with the measured data, and all three turbulence models (k-epsilon, algebraic eddy viscosity, and ARSM) show the same level of agreement. The prediction of boundary-layer development on rotor blades shows much better agreement with measurements with the ARSM. It is essential to employ higher-order turbulence models to capture the effects of rotation, curvature, and three dimensionality on boundary layers in turbomachinery.

Zhang, J.↗
NASA NTRS1979

Large Eddy Simulation of turbulent shear flows

The conceptual foundation underlying Large Eddy Simulation (LES) is summarized, and the numerical methods developed for simulation of the time-developing turbulent mixing layer and turbulent plane Poiseuille flow are discussed. Computational results show that the average Reynolds stress profile nearly attains the equilibrium shape which balances the downstream mean pressure gradient in the regions away from the walls. In the vicinity of the walls, viscous stresses are shown to be significant; together with the Reynolds stresses, these stresses balance the mean pressure gradient. It is stressed that the subgrid scale contribution to the total Reynolds stress is significant only in the vicinity of the walls. The continued development of LES is urged.

Moin, P.↗
NASA NTRS2001

Observational Tests of Recent MHD Turbulence Perspectives

This grant seeks to analyze the Heliospheric Missions data to test current theories on the angular dependence (with respect to mean magnetic field direction) of magnetohydrodynamic (MHD) turbulence in the solar wind. Solar wind turbulence may be composed of two or more dynamically independent components. Such components include magnetic pressure-balanced structures, velocity shears, quasi-2D turbulence, and slab (Alfven) waves. We use a method, developed during the first two years of this grant, for extracting the individual reduced spectra of up to three separate turbulence components from a single spacecraft time series. The method has been used on ISEE-3 data, Pioneer Venus Orbiter, Ulysses, and Voyager data samples. The correlation of fluctuations as a function of angle between flow direction and magnetic-field direction is the focus of study during the third year.

Ghosh, Sanjoy↗
NASA NTRSDate not supplied

Validation of Two-Phase Flow Modeling Techniques in Loci-Stream Using Axial Jet Pressure Control and Autogenous Pressurization Experiments

Improvements to computational fluid dynamics modeling of cryogenic fluid management operations are presented with experimental validations. Modeling improvements includes utilizing improved delayed detaches eddy simulation to more effectively model turbulence compared to the Reynolds-averaged Navier Stokes modeling, modeling temperature dependent tank solid properties to capture large variation in thermodynamic properties at cryogenic conditions, modeling real gas equation of state for hydrogen instead of ideal gas behavior, and an on-going effort to improve interface turbulence transport for the volume of fluid method. Two flight-scale ground tests, five K-site tank axial jet pressure control experiments and an Engineering Design Unit tank autogenous pressurization experiment are validated and show the effectiveness of the improved modeling techniques. The successful validation establishes further trust in these models for flight project support and provides a pathway for maturation of in-development interface turbulence damping model for the volume of fluid method for more generalized applications.

CFD↗