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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Navier-Stokes computation of compressible turbulent flows with a second order closure, part 1

A second order closure turbulence model for compressible flows is developed and implemented in a 2D Reynolds-averaged Navier-Stokes solver. From the beginning where a kappa-epsilon turbulence model was implemented in the bidiagonal implicit method of MACCORMACK (referred to as the MAC3 code) to the final stage of implementing a full second order closure in the efficient line Gauss-Seidel algorithm, numerous work was done, individually and collectively. Besides the collaboration itself, the final product of this work is a second order closure derived from the Launder, Reece, and Rodi model to account for near wall effects, which has been called FRAME model, which stands for FRench-AMerican-Effort. During the reporting period, two different problems were worked out. The first was to provide Ames researchers with a reliable compressible boundary layer code including a wide collection of turbulence models for quick testing of new terms, both in two equations and in second order closure (LRR and FRAME). The second topic was to complete the implementation of the FRAME model in the MAC5 code. The work related to these two different contributions is reported. dilatation in presence of stron shocks. This work, which has been conducted during a work at the Center for Turbulence Research with Zeman aimed also to cros-check earlier assumptions by Rubesin and Vandromme.

Haminh, Hieu↗

NASA NTRS1984

Maximum likelihood estimation of turbulence spectrum parameters

Estimation of the integral scale and intensity of a generic turbulence record is treated as a statistical problem of parameter estimation. Properties of parameter estimators and the method of maximum likelihood are reviewed. Likelihood equations are derived for estimation of the integral scale and intensity applicable to a general class of turbulence spectra that includes the von Karman and Dryden transverse and longitudinal spectra as special cases. The method is extended to include the Bullen transverse and longitudinal spectra. Coefficients of variation are given for maximum likelihood estimates of the integral scale and intensity of the von Karman spectra. Application of the method is illustrated by estimating the integral scale and intensity of an atmospheric turbulence vertical velocity record assumed to be governed by the von Karman transverse spectrum.

Mark, W. D.↗

NASA NTRS1953

Notes on the Prediction of Shock-induced Boundary-layer Separation

The present status of available information relative to the prediction of shock-induced boundary-layer separation is discussed. Experimental results showing the effects of Reynolds number and Mach number on the separation of both laminar and turbulent boundary layer are given and compared with available methods for predicting separation. The flow phenomena associated with separation caused by forward-facing steps, wedges, and incident shock waves are discussed. Applications of the flat-plate data to problems of separation on spoilers, diffusers, and scoop inlets are indicated for turbulent boundary layers.

Lange, Roy H.↗

NASA NTRS1983

Psuedospectral calculation of shock turbulence interactions

A Chebyshev-Fourier discretization with shock fitting is used to solve the unsteady Euler equations. The method is applied to shock interactions with plane waves and with a simple model of homogeneous isotropic turbulence. The plane wave solutions are compared to linear theory.

Zang, T. A.↗

NASA NTRS1983

Pseudospectral calculation of shock turbulence interactions

Zang, T. A.↗

NASA NTRS1984

Formation and inflammation of a turbulent jet

The formation and inflammation of a planar, turbulent jet in an incompressible medium is modeled numerically by the use of the random vortex method amended by a flame propagation algorithm. The results demonstrate the dominant influence of turbulent eddies and their interactions upon the development of the jet. Its growth is shown to consist of three stages: formation of small eddies, pairing of eddies with the same sign of circulation, and pairing of eddies of opposite signs. On this basis a number of features of the jet mechanism are revealed, namely penetration, engulfment, entrainment, and intermittency. Two cases of inflammation are considered. In one, the jet is ignited at the center of the orifice, the solution tracing its own inflammation. In the other, combustion is initiated across its full cross section, the results modeling the action of a turbulent torch as it spreads the flame into the combustible surroundings. In both cases the flow field is still dominated by the turbulent eddies and their interactions. However, the coherence among them is encumbered as a consequence of expansion due to the exothermicity of the combustion process.

Ghoniem, A. F.↗

NASA NTRS2008

Application of FUN3D and CFL3D to the Third Workshop on CFD Uncertainty Analysis

Two Reynolds-averaged Navier-Stokes computer codes - one unstructured and one structured - are applied to two workshop cases (for the 3rd Workshop on CFD Uncertainty Analysis, held at Instituto Superior Tecnico, Lisbon, in October 2008) for the purpose of uncertainty analysis. The Spalart-Allmaras turbulence model is employed. The first case uses the method of manufactured solution and is intended as a verification case. In other words, the CFD solution is expected to approach the exact solution as the grid is refined. The second case is a validation case (comparison against experiment), for which modeling errors inherent in the turbulence model and errors/uncertainty in the experiment may prevent close agreement. The results from the two computer codes are also compared. This exercise verifies that the codes are consistent both with the exact manufactured solution and with each other. In terms of order property, both codes behave as expected for the manufactured solution. For the backward facing step, CFD uncertainty on the finest grid is computed and is generally very low for both codes (whose results are nearly identical). Agreement with experiment is good at some locations for particular variables, but there are also many areas where the CFD and experimental uncertainties do not overlap.

Rumsey, C. L.↗

NASA NTRS1996

Multigrid Acceleration of Time-Accurate DNS of Compressible Turbulent Flow

An efficient scheme for the direct numerical simulation of 3D transitional and developed turbulent flow is presented. Explicit and implicit time integration schemes for the compressible Navier-Stokes equations are compared. The nonlinear system resulting from the implicit time discretization is solved with an iterative method and accelerated by the application of a multigrid technique. Since we use central spatial discretizations and no artificial dissipation is added to the equations, the smoothing method is less effective than in the more traditional use of multigrid in steady-state calculations. Therefore, a special prolongation method is needed in order to obtain an effective multigrid method. This simulation scheme was studied in detail for compressible flow over a flat plate. In the laminar regime and in the first stages of turbulent flow the implicit method provides a speed-up of a factor 2 relative to the explicit method on a relatively coarse grid. At increased resolution this speed-up is enhanced correspondingly.

Broeze, Jan↗

NASA NTRS1982

Flow aerodynamics modeling of an MHD swirl combustor - Calculations and experimental verification

This paper describes a computer code for calculating the flow dynamics of constant density flow in the second stage trumpet shaped nozzle section of a two stage MHD swirl combustor for application to a disk generator. The primitive pressure-velocity variable, finite difference computer code has been developed to allow the computation of inert nonreacting turbulent swirling flows in an axisymmetric MHD model swirl combustor. The method and program involve a staggered grid system for axial and radial velocities, and a line relaxation technique for efficient solution of the equations. Turbulence simulation is by way of a two-equation Kappa-epsilon model. The code produces as output the flowfield map of the nondimensional stream function, axial, and swirl velocity. Good agreement was obtained between the theoretical predictions and the qualitative experimental results. The best seed injector location for uniform seed distribution at combustor exit is with injector located centrally on the combustor axis at entrance to the second stage combustor.

Gupta, A. K.↗

NASA NTRS1990

Numerical study of turbulent secondary flows in curved ducts

The pressure driven, fully-developed turbulent flow of an incompressible viscous fluid in curved ducts of square-section is studied numerically by making use of a finite volume method. A nonlinear Kappa - Iota model is used to represent the turbulence. The results for both straight and curved ducts are presented. For the case of fully-developed turbulent flow in straight and curved ducts, the secondary flow is characterized by an eight-vortex structure for which the computed flowfield is shown to be in good agreement with available experimental data. The introduction of moderate curvature is shown to cause a substantial increase in the strength of the secondary flow and to change the secondary flow pattern to either a double-vortex or a four-vortex configuration.

Hur, N.↗

NASA NTRS1985

Turbulent disruptions from the Strauss equations

Preliminary results are reported from application of a three-dimensional spectral method model to the solution of the Strauss (1976) reduced MHD equations. The investigation was focused on describing MHD turbulence in a current-carrying bounded magnetofluid. A cylindrical geometry with a square cross-section was considered, with the walls being rigid perfect conductors with free-slip boundary conditions. A uniform magnetic field and the electric current density both point in the z-direction. Initial conditions are specified which feature small amounts of random noise expressed as Fourier modes. Linearized equations are defined for tracing the movement to equlibrium conditions or other temporal development. The model is further refined with nonlinear equations to examine the effects of the appearance of disruptions. Comparisons are drawn between solutions obtained with linear and nonlinear equations, with an eye to the associated physical realities.

Dahlburg, J. P.↗

NASA NTRSDate not supplied

Turbulent Simulations of Cooling Jets in Crossflow

Turbulent flow simulations were completed of the THX-III single hole film cooling experiment performed at NASA Glenn Research Center (GRC). Two different Large-Eddy Simulation (LES) approaches were completed along with a RANS based study. The first LES approach utilizes the Flux Reconstruction (FR) based finite element code GFR, while the second approach makes use of the compact difference solver FDL3DI. Both LES codes were run on grids larger than 1.5 Billion degrees of freedom (DoF). The RANS simulations used FUN3D and included a grid convergence as a turbulence closure study. A summary of the computational approaches and key findings are presented along with the experimental measurements. Significant progress has been made towards the completion of this work, as all but one simulation have been completed. The preliminary findings have shown the LES approaches have excellent agreement with each other, and closely follow the trends of the experiment. The Reynolds stress RANS closure, while the best of the RANS results, had difficulty matching the results of the LES and experimental measurements.

High-Order Methods↗

NASA NTRS2019

Modeling of an Atmospheric-Boundary-Layer Profile in Support of Experiments in the NASA Langley Transonic Dynamics Tunnel

This paper presents the results of the Computational Fluid Dynamics (CFD) analyses to model two atmospheric-boundary-layer (ABL) profiles inside the NASA Langley Research Center Transonic Dynamics Tunnel (TDT). The CFD models include tunnel walls and all additional hardware needed to simulate the ABL profiles. This hardware includes Irwin spires and floor-mounted roughness elements. The numerical simulations show that the application of higher-fidelity numerical methods is necessary to compute boundary-layer and turbulenceintensity profiles that match experimental data. The ABL profiles are computed both inside the numerical model of the TDT and in a classical free-air model. Both Unsteady Reynolds Averaged Navier-Stokes (URANS) with Spalart-Allmaras (SA) turbulence model and Modified Delayed Detached Eddy (MDDES) simulation methods are used. The results show that the MDDES-simulation results match the experimental data very well while URANS-SA does not. The results also show that the wind-tunnel walls have a significant effect on ABL prediction.

Chwalowski, Pawel↗

NASA NTRS2011

An Upwind Solver for the National Combustion Code

An upwind solver is presented for the unstructured grid National Combustion Code (NCC). The compressible Navier-Stokes equations with time-derivative preconditioning and preconditioned flux-difference splitting of the inviscid terms are used. First order derivatives are computed on cell faces and used to evaluate the shear stresses and heat fluxes. A new flux limiter uses these same first order derivatives in the evaluation of left and right states used in the flux-difference splitting. The k-epsilon turbulence equations are solved with the same second-order method. The new solver has been installed in a recent version of NCC and the resulting code has been tested successfully in 2D on two laminar cases with known solutions and one turbulent case with experimental data.

Sockol, Peter M.↗

NASA NTRS1989

Numerical study of turbulent secondary flows in curved ducts

The pressure driven, fully-developed turbulent flow of an incompressible viscous fluid in curved ducts of square cross-section is studied numerically by making use of a finite volume method. A nonlinear Kappa - Iota model is used to represent the turbulence. The results for both straight and curved ducts are presented. For the case of fully-developed turbulent flow in straight ducts, the secondary flow is characterized by an eight-vortex structure for which the computed flowfield is shown to be in good agreement with available experimental data. The introduction of moderate curvature is shown to cause a substantial increase in the strength of the secondary flow and to change the secondary flow pattern to either a double-vortex or a four-vortex configuration.

Hur, N.↗

NASA NTRS2001

A 7.5-Year Dataset of SSM/I-Derived Surface Turbulent Fluxes Over Global Oceans

The global air-sea turbulent fluxes are needed for driving ocean models and validating coupled ocean-atmosphere global models. A method was developed to retrieve surface air humidity from the radiances measured by the Special Sensor Microwave/Imager (SSM/I) Using both SSM/I-retrieved surface wind and air humidity, they computed daily turbulent fluxes over global oceans with a stability-dependent bulk scheme. Based on this method, we have produced Version 1 of Goddard Satellite-Based Surface Turbulent Fluxes (GSSTF) dataset from the SSM/I data and other data. It provides daily- and monthly-mean surface turbulent fluxes and some relevant parameters over global oceans for individual F8, F10, and F11 satellites covering the period July 1987-December 1994. It also provides 1988-94 annual- and monthly-mean climatologies of the same variables, using only F8 and F1 1 satellite data. It has a spatial resolution of 2.0 degrees x 2.5 degrees lat-long and is archived at the NASA/GSFC DAAC. The purpose of this paper is to present an updated assessment of the GSSTF 1.0 dataset.

Chou, Shu-Hsien↗

NASA NTRS1977

Theory of viscous transonic flow over airfoils at high Reynolds number

This paper considers viscous flows with unseparated turbulent boundary layers over two-dimensional airfoils at transonic speeds. Conventional theoretical methods are based on boundary layer formulations which do not account for the effect of the curved wake and static pressure variations across the boundary layer in the trailing edge region. In this investigation an extended viscous theory is developed that accounts for both effects. The theory is based on a rational analysis of the strong turbulent interaction at airfoil trailing edges. The method of matched asymptotic expansions is employed to develop formal series solutions of the full Reynolds equations in the limit of Reynolds numbers tending to infinity. Procedures are developed for combining the local trailing edge solution with numerical methods for solving the full potential flow and boundary layer equations. Theoretical results indicate that conventional boundary layer methods account for only about 50% of the viscous effect on lift, the remaining contribution arising from wake curvature and normal pressure gradient effects.

Melnik, R. E.↗

NASA NTRS2009

Motion Cueing Algorithm Modification for Improved Turbulence Simulation

Atmospheric turbulence cueing produced by flight simulator motion systems has been less than satisfactory because the turbulence profiles have been attenuated by the motion cueing algorithms. Cardullo and Ellor initially addressed this problem by directly porting the turbulence model output to the motion system. Reid and Robinson addressed the problem by employing a parallel aircraft model, which is only stimulated by the turbulence inputs and adding a filter specially designed to pass the higher turbulence frequencies. There have been advances in motion cueing algorithm development at the Man-Machine Systems Laboratory, at SUNY Binghamton. In particular, the system used to generate turbulence cues has been studied. The Reid approach, implemented by Telban and Cardullo, was employed to augment the optimal motion cueing algorithm installed at the NASA LaRC Simulation Laboratory, driving the Visual Motion Simulator. In this implementation, the output of the primary flight channel was added to the output of the turbulence channel and then sent through a non-linear cueing filter. The cueing filter is an adaptive filter; therefore, it is not desirable for the output of the turbulence channel to be augmented by this type of filter. The likelihood of the signal becoming divergent was also an issue in this design. After testing on-site it became apparent that the architecture of the turbulence algorithm was generating unacceptable cues. As mentioned above, this cueing algorithm comprised a filter that was designed to operate at low bandwidth. Therefore, the turbulence was also filtered, augmenting the cues generated by the model. If any filtering is to be done to the turbulence, it will utilize a filter with a much higher bandwidth, above the frequencies produced by the aircraft response to turbulence. The authors have developed an implementation wherein only the signal from the primary flight channel passes through the nonlinear cueing filter. This paper discusses three new algorithms. Testing shows that the new methods provide the pilot with a more realistic sensation of turbulence; the cues are not attenuated by algorithm. Results of offline testing show the credibility of the models. Offline test verification was based primarily on the evaluation of the power spectral density of the outputs and the time response.

Ercole, Anthony V.↗