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 109 records · Page 6

Application of viscous-inviscid interaction methods to transonic turbulent flows

Two different viscous-inviscid interaction schemes were developed for the analysis of steady, turbulent, transonic, separated flows over axisymmetric bodies. The viscous and inviscid solutions are coupled through the displacement concept using a transpiration velocity approach. In the semi-inverse interaction scheme, the viscous and inviscid equations are solved in an explicitly separate manner and the displacement thickness distribution is iteratively updated by a simple coupling algorithm. In the simultaneous interaction method, local solutions of viscous and inviscid equations are treated simultaneously, and the displacement thickness is treated as an unknown and is obtained as a part of the solution through a global iteration procedure. The inviscid flow region is described by a direct finite-difference solution of a velocity potential equation in conservative form. The potential equation is solved on a numerically generated mesh by an approximate factorization (AF2) scheme in the semi-inverse interaction method and by a successive line overrelaxation (SLOR) scheme in the simultaneous interaction method. The boundary-layer equations are used for the viscous flow region. The continuity and momentum equations are solved inversely in a coupled manner using a fully implicit finite-difference scheme.

Lee, D.↗

An eddy viscosity calculation method for a turbulent duct flow

The mean velocity profile across a fully developed turbulent duct flow is obtained from an eddy viscosity relation combined with an empirical outer region wake function. Results are in good agreement with experiments and with direct numerical simulations in the same flow at two Reynolds numbers. In particular, the near-wall trend of the Reynolds shear stress and its variation with Reynolds number are similar to those of the simulations. The eddy viscosity method is more accurate than previous mixing length or implicit function methods.

Antonia, R. A.↗

A dispersion relationship governing incompressible wall turbulence

The method of separation of variables is shown to make turbulent correlation equations of Karman-Howarth type tractable for shear turbulence as well under the condition of neglected triple correlation. The separated dependent variable obeys an Orr-Sommerfeld equation. A new analytical method is developed using a scaling law different from the classical one due to Heisenberg and Lin and more appropriate for wall turbulent profiles. A dispersion relationship between the wave number and the separation constant which has the dimension of a frequency is derived in support of experimental observations of wave or coherent structure of wall turbulence.

Tsuge, S.↗

Dry Deposition Methods Based on Turbulence Kinetic Energy: 1. Evaluation of Various Resistances and Sensitivity Studies Using a Single‐Point Model

Different functions are used to account for turbulence strength in the atmospheric boundary layer for different stability regimes. These functions are one of the sources for differences among different atmospheric models' predictions and associated biases. Also, turbulence strength is underrepresented in some of the resistance formulations. Here, to address these issues with dry deposition, firstly we take advantage of three‐dimensional (3‐D) turbulence information in estimating resistances by proposing and validating a 3‐D turbulence velocity scale that is relevant for different stability regimes of boundary layer. Secondly, we hypothesize and validate that friction velocity measured by 3‐D sonic anemometer can be effectively replaced by the new turbulence velocity scale multiplied by the von Karman constant. Finally, we (a) present a set of resistance formulations for ozone (O 3 ) based on the 3‐D turbulence velocity scale; (b) intercompare estimations of such resistances with those obtained using existing formulations; and, (c) evaluate simulated O 3 fluxes using a single‐point dry deposition model against long‐term observations of O 3 fluxes at the Harvard Forest (MA) site. Results indicate that the new resistance formulations work very well in simulating surface latent heat and O 3 fluxes when compared to respective existing formulations and measurements at a decadal time scale. Findings from this research may help to improve the capability of dry deposition schemes for better estimation of dry deposition fluxes and create opportunities for the development of a community dry deposition model for use in regional/global air quality models.

54 ENVIRONMENTAL SCIENCES↗

STC-SAB program users manual for the turbulent boundary layer and turbulent separation prediction methods employed in the NASA Langley streamtube curvature computer program

The streamtube curvature program (STC) has been developed to predict the inviscid flow field and the pressure distribution about nacelles at transonic speeds. The effects of boundary layer are to displace the inviscid flow and effectively change the body shape. Thus, the body shape must be corrected by the displacement thickness in order to calculate the correct pressure distribution. This report describes the coupling of the Stratford and Beavers boundary layer solution with the inviscid STC analysis so that all nacelle pressure forces, friction drag, and incipient separation may be predicted. The usage of the coupled STC-SAB computer program is outlined and the program input and output are defined. Included in this manual are descriptions of the principal boundary layer tables and other revisions to the STC program. The use of the viscous option is controlled by the engineer during program input definition.

Ferguson, D. R.↗

Essentially non-oscillatory shock capturing methods applied to turbulence amplification in shock wave calculations

ENO (essentially non-oscillatory) schemes can provide uniformly high order accuracy right up to discontinuities while keeping sharp, essentially non-oscillatory shock transitions. Recently, an efficient implementation of ENO schemes was obtained based on fluxes and TVD Runge-Kutta time discretizations. The resulting code is very simple to program for multi-dimensions. ENO schemes are especially suitable for computing problems with both discontinuities and fine structures in smooth regions, such as shock interaction with turbulence, for which results for 1-D and 2-D Euler equations are presented. Much better resolution is observed by using third order ENO schemes than by using second order TVD schemes for such problems.

Osher, Stanley↗

A defect stream function, law of the wall/wake method for compressible turbulent boundary layers

The application of the defect stream function to the solution of the two-dimensional, compressible boundary layer is examined. A law of the wall/law of the wake formulation for the inner part of the boundary layer is presented which greatly simplifies the computational task near the wall and eliminates the need for an eddy viscosity model in this region. The eddy viscosity model in the outer region is arbitrary. The modified Crocco temperature-velocity relationship is used as a simplification of the differential energy equation. Formulations for both equilibrium and nonequilibrium boundary layers are presented including a constrained zero-order form which significantly reduces the computational workload while retaining the significant physics of the flow. A formulation for primitive variables is also presented. Results are given for the constrained zero-order and second-order equilibrium formulations and are compared with experimental data. A compressible wake function valid near the wall has been developed from the present results.

Barnwell, Richard W.↗

More accurate predictions with transonic Navier-Stokes methods through improved turbulence modeling

Significant improvements in predictive accuracies for off-design conditions are achievable through better turbulence modeling; and, without necessarily adding any significant complication to the numerics. One well established fact about turbulence is it is slow to respond to changes in the mean strain field. With the 'equilibrium' algebraic turbulence models no attempt is made to model this characteristic and as a consequence these turbulence models exaggerate the turbulent boundary layer's ability to produce turbulent Reynolds shear stresses in regions of adverse pressure gradient. As a consequence, too little momentum loss within the boundary layer is predicted in the region of the shock wave and along the aft part of the airfoil where the surface pressure undergoes further increases. Recently, a 'nonequilibrium' algebraic turbulence model was formulated which attempts to capture this important characteristic of turbulence. This 'nonequilibrium' algebraic model employs an ordinary differential equation to model the slow response of the turbulence to changes in local flow conditions. In its original form, there was some question as to whether this 'nonequilibrium' model performed as well as the 'equilibrium' models for weak interaction cases. However, this turbulence model has since been further improved wherein it now appears that this turbulence model performs at least as well as the 'equilibrium' models for weak interaction cases and for strong interaction cases represents a very significant improvement. The performance of this turbulence model relative to popular 'equilibrium' models is illustrated for three airfoil test cases of the 1987 AIAA Viscous Transonic Airfoil Workshop, Reno, Nevada. A form of this 'nonequilibrium' turbulence model is currently being applied to wing flows for which similar improvements in predictive accuracy are being realized.

Johnson, Dennis A.↗

Turbulence modeling for high speed compressible flows

Turbulence models are necessary in numerical simulations because of the impracticality of computing all scales of turbulent motion. Approximate methods, or turbulence models, are introduced to simplify and make the computations practical. There are several approaches to turbulence modeling: subgrid scale modeling; modeling of all the scales of turbulent motion; eddy viscosity; and Reynolds stress turbulence closure models. In many instances, computations based on turbulence models are compared with available experimental data. The effort is directed at seeking models which have reasonable accuracy over a limited range of flow conditions. Future needs are briefly discussed. Work to date is outlined.

Chandra, Suresh↗

On the Yakhot-Orszag renormalization group method for deriving turbulence statistics and models

An independent, comprehensive, critical review of the 'renormalization group' (RNG) theory of turbulence developed by Yakhot and Orszag (1986) is provided. Their basic theory for the Navier-Stokes equations is confirmed, and approximations in the scale removal procedure are discussed. The YO derivations of the velocity-derivative skewness and the transport equation for the energy dissipation rate are examined. An algebraic error in the derivation of the skewness is corrected. The corrected RNG skewness value of -0.59 is in agreement with experiments at moderate Reynolds numbers. Several problems are identified in the derivation of the energy dissipation rate equations which suggest that the derivation should be reformulated.

Smith, L. M.↗

A Method for Computing Turbulent Heat Transfer in the Presence of a Streamwise Pressure Gradient for Bodies in High-Speed Flow

A modified Reynolds analogy between skin friction and heat transfer which depends upon local pressure gradient is derived. Exact and approximate solutions are derived from the differential equations; the exact solution is applicable for arbitrary initial (transition) conditions and the approximate solution requires fully developed turbulent flow from stagnation point or leading edge. The exact solution (restricted to stagnation initial conditions) and the approximate solutions are shown to agree with one another within 5 percent when applied to several blunt shapes. The present solutions generally predict the measured heating rates on these bodies within the accuracy of the measurements provided transition began upstream of the peak heating region. The present solutions appear to be sufficiently accurate for design purposes.

Cohen, Nathaniel B.↗

Exploring a Method for Improving Turbulent Separated-Flow Predictions with kappa-omega Models

A particular failing of Reynolds-averaged Navier-Stokes separated turbulent flow computations is addressed within the context of a kappa-omega two-equation turbulence model. The failing is the tendency for turbulence models to under-predict turbulent shear stress in the shear layers of some separation bubbles, yielding late boundary layer reattachment and recovery. Inspired by unpublished work of Volker, Langtry, and Menter, the author undertook an independent investigation in an attempt to improve the ability of the Menter shear stress transport (SST) model to predict flowfield characteristics in and downstream of separation bubbles. The fix is an ad hoc term that is a function of the local ratio of turbulent production to dissipation; it is used to multiply the omega-destruction term, increasing eddy viscosity in separated regions. With this fix, several flowfields are investigated. Results show that, although the "separation fix" can provide dramatic improvement in some cases, it is not consistently good for all flows. Thus, although it may prove helpful in many situations in its current form, this model may benefit from further refinements, including better sensitization to the energetics of turbulence in the separated region.

Rumsey, Christopher L.↗

A method for the analysis of nonlinearities in aircraft dynamic response to atmospheric turbulence

An analytical method is developed which combines the equivalent linearization technique for the analysis of the response of nonlinear dynamic systems with the amplitude modulated random process (Press model) for atmospheric turbulence. The method is initially applied to a bilinear spring system. The analysis of the response shows good agreement with exact results obtained by the Fokker-Planck equation. The method is then applied to an example of control-surface displacement limiting in an aircraft with a pitch-hold autopilot.

Sidwell, K.↗