Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence modeling”

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,009 records · Page 56

RANS Modeling of Benchmark Shockwave / Boundary Layer Interaction Experiments

This presentation summarizes the computations of a set of shock wave / turbulent boundary layer interaction (SWTBLI) test cases using the Wind-US code, as part of the 2010 American Institute of Aeronautics and Astronautics (AIAA) shock / boundary layer interaction workshop. The experiments involve supersonic flows in wind tunnels with a shock generator that directs an oblique shock wave toward the boundary layer along one of the walls of the wind tunnel. The Wind-US calculations utilized structured grid computations performed in Reynolds-averaged Navier-Stokes mode. Three turbulence models were investigated: the Spalart-Allmaras one-equation model, the Menter Shear Stress Transport wavenumber-angular frequency two-equation model, and an explicit algebraic stress wavenumber-angular frequency formulation. Effects of grid resolution and upwinding scheme were also considered. The results from the CFD calculations are compared to particle image velocimetry (PIV) data from the experiments. As expected, turbulence model effects dominated the accuracy of the solutions with upwinding scheme selection indicating minimal effects.!

Georgiadis, Nick↗

Stability and Interaction of Coherent Structure in Supersonic Reactive Wakes

A theoretical formulation and analysis is presented for a study of the stability and interaction of coherent structure in reacting free shear layers. The physical problem under investigation is a premixed hydrogen-oxygen reacting shear layer in the wake of a thin flat plate. The coherent structure is modeled as a periodic disturbance and its stability is determined by the application of linearized hydrodynamic stability theory which results in a generalized eigenvalue problem for reactive flows. Detailed stability analysis of the reactive wake for neutral, symmetrical and antisymmetrical disturbance is presented. Reactive stability criteria is shown to be quite different from classical non-reactive stability. The interaction between the mean flow, coherent structure and fine-scale turbulence is theoretically formulated using the von-Kaman integral technique. Both time-averaging and conditional phase averaging are necessary to separate the three types of motion. The resulting integro-differential equations can then be solved subject to initial conditions with appropriate shape functions. In the laminar flow transition region of interest, the spatial interaction between the mean motion and coherent structure is calculated for both non-reactive and reactive conditions and compared with experimental data wherever available. The fine-scale turbulent motion determined by the application of integral analysis to the fluctuation equations. Since at present this turbulence model is still untested, turbulence is modeled in the interaction problem by a simple algebraic eddy viscosity model. The applicability of the integral turbulence model formulated here is studied parametrically by integrating these equations for the simple case of self-similar mean motion with assumed shape functions. The effect of the motion of the coherent structure is studied and very good agreement is obtained with previous experimental and theoretical works for non-reactive flow. For the reactive case, lack of experimental data made direct comparison difficult. It was determined that the growth rate of the disturbance amplitude is lower for reactive case. The results indicate that the reactive flow stability is in qualitative agreement with experimental observation.

Menon, Suresh↗

Molecular fluctuations inhibit intermittency in compressible turbulence

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

compressible turbulence↗

Turbulent Axial Odometer Model

A Galilean-invariant field equation is proposed and tested on standard turbulence model test cases. The field equation provides an additional non-dimensional outer scale which allows the turbulence models to reproduce the axial normal stress increase with Re seen in high Reynolds numbers experiments. The field equation provides a Reynolds number for every point based on the length of turbulent flow upstream of that point in the domain. This outer scale equation can be considered an odometer that gives a length scale conjectured to be related to the large stream wise structures that are seen in turbulent flow and that require run length to develop. A new RANS model using this additional scale is able to match the Reynolds number variation of the normal stresses seen at high Reynolds number. Furthermore, the good attached flow prediction capabilities of current RANS models appears to be attained. Using this scale equation, the entire Reynolds-stress state appears to be predicted correctly, over a large run length Reynolds number range such as experienced in aircraft design.

Olsen, Michael E.↗

Turbulent Axial Odometer Model

A Galilean-invariant field equation is proposed and tested on standard turbulence model test cases. The field equation provides an additional non-dimensional outer scale which allows the turbulence models to reproduce the axial normal stress increase with Re seen in high Reynolds numbers experiments. The field equation provides a Reynolds number for every point based on the length of turbulent flow upstream of that point in the domain. This outer scale equation can be considered an odometer that gives a length scale conjectured to be related to the large streamwise structures that are seen in turbulent flow and that require run length to develop. A new RANS model using this additional scale is able to match the Reynolds number variation of the normal stresses seen at high Reynolds number. Furthermore, the good attached flow prediction capabilities of current RANS models appears to be attained. Using this scale equation, the entire Reynolds-stress state appears to be predicted correctly, over a large run length Reynolds number range such as experienced in aircraft design.

Olsen, Michael E.↗

Compressible turbulent flows: Modeling and similarity considerations

With the recent revitalization of high speed flow research, compressibility presents a new set of challenging problems to turbulence researchers. Questions arise as to what extent compressibility affects turbulence dynamics, structures, the Reynolds stress-mean velocity (constitutive) relation, and the accompanying processes of heat transfer and mixing. In astrophysical applications, compressible turbulence is believed to play an important role in intergalactic gas cloud dynamics and in accretion disk convection. Understanding and modeling of the compressibility effects in free shear flows, boundary layers, and boundary layer/shock interactions is discussed.

Zeman, Otto↗

Features of a reattaching turbulent shear layer subject to an adverse pressure gradient

Experimental data have been obtained in an incompressible turbulent flow over a rearward-facing step with superimposed adverse pressure gradient. Mean velocities, Reynolds stresses and triple-products measured by a laser Doppler velocimeter are presented for two cases of adverse pressure gradient. Mixing lengths, eddy viscosities, production, convection, turbulent diffusion, and dissipation terms are extracted from the data. These data are compared with various mixing length and eddy-viscosity turbulence models. Numerical calculations incorporating the k-epsilon and the algebraic-stress turbulence models are compared with the data. When determining quantities of engineering interest, the modified algebraic-stress model (ASM) is a significant improvement over the unmodified ASM and the unmodified k-epsilon model

Driver, D. M.↗

Recent advances in Runge-Kutta schemes for solving 3-D Navier-Stokes equations

A thin-layer Navier-Stokes has been developed for solving high Reynolds number, turbulent flows past aircraft components under transonic flow conditions. The computer code has been validated through data comparisons for flow past isolated wings, wing-body configurations, prolate spheroids and wings mounted inside wind-tunnels. The basic code employs an explicit Runge-Kutta time-stepping scheme to obtain steady state solution to the unsteady governing equations. Significant gain in the efficiency of the code has been obtained by implementing a multigrid acceleration technique to achieve steady-state solutions. The improved efficiency of the code has made it feasible to conduct grid-refinement and turbulence model studies in a reasonable amount of computer time. The non-equilibrium turbulence model of Johnson and King has been extended to three-dimensional flows and excellent agreement with pressure data has been obtained for transonic separated flow over a transport type of wing.

Vatsa, Veer N.↗

High-Fidelity Analysis of Lift+Cruise VTOL Urban Air Mobility Concept Aircraft

This paper presents a high-fidelity multidisciplinary analysis for the NASA lift+cruise vertical takeoff and landing urban air mobility concept aircraft. The reported simulations couple comprehensive rotorcraft aeromechanics and high-fidelity computational fluid dynamics (CFD). Aerodynamic solutions are computed on dynamic, deforming, unstructured, overset grid systems by an unsteady CFD solver, FUN3D, developed at the NASA Langley Research Center. An integrated overset-grid assembler, Yoga, is used for communications between component grids. Two turbulence models are compared for loose-coupling simulations in hover and low-speed forward flight conditions. One model is the negative variant of the one-equation Spalart-Allmaras model (SA-neg), and the other model is the SA-neg-R model adding a simple rotation correction term to the SA-neg model. It demonstrated that the rotation correction in the SA-neg-R turbulence model significantly improves resolution of secondary vortices and wake interactions for rotorcraft simulations.

CFD↗

High-Fidelity Analysis of Lift+Cruise VTOL Urban Air Mobility Concept Aircraft

This paper presents a high-fidelity multidisciplinary analysis for the NASA Lift+Cruise vertical takeoff and landing urban air mobility concept aircraft. The reported simulations couple comprehensive rotorcraft aeromechanics and high-fidelity computational fluid dynamics (CFD). Aerodynamic solutions are computed on dynamic, deforming, unstructured, overset grid systems by an unsteady CFD solver, FUN3D, developed at the NASA Langley Research Center. An integrated overset-grid assembler, Yoga, is used for communications between component grids. Two turbulence models are compared for loose-coupling simulations in hover and low-speed forward flight conditions. One model is the negative variant of the one-equation Spalart-Allmaras model (SA-neg), and the other model is the SA-neg-R model adding a simple rotation correction term to the SA-neg model. It demonstrated that the rotation correction in the SA-neg-R turbulence model significantly improves resolution of secondary vortices and wake interactions for rotorcraft simulations.

CFD↗

Advanced in turbulence physics and modeling by direct numerical simulations

The advent of direct numerical simulations of turbulence has opened avenues for research on turbulence physics and turbulence modeling. Direct numerical simulation provides values for anything that the scientist or modeler would like to know about the flow. An overview of some recent advances in the physical understanding of turbulence and in turbulence modeling obtained through such simulations is presented.

Reynolds, W. C.↗

An algebraic model for dissipation in supersonic boundary layers

With increasing mainstream Mach number, viscous dissipation becomes a progressively important influence in high-speed compressible turbulent boundary layers. An asymptotic analysis is carried out for high Reynolds numbers and Mach numbers of order 1, and it is shown that viscous dissipation gives rise to important terms in the solution of the total enthalpy equation. For simplicity, the case of supersonic flow over an adiabatic wall is considered. An expression for the adiabatic wall temperature is derived. It is shown that the asymptotic analysis constrains the types of turbulence models that can be used to represent the effects of viscous dissipation. A simple algebraic turbulence model is proposed and comparisons with measured total enthalpy profile data show good agreement.

He, J.↗

Observations, theoretical ideas and modeling of turbulent flows: Past, present and future

Turbulence was analyzed in a historical context featuring the interactions between observations, theoretical ideas, and modeling within three successive movements. These are identified as predominantly statistical, structural and deterministic. The statistical movement is criticized for its failure to deal with the structural elements observed in turbulent flows. The structural movement is criticized for its failure to embody observed structural elements within a formal theory. The deterministic movement is described as having the potential of overcoming these deficiencies by allowing structural elements to exhibit chaotic behavior that is nevertheless embodied within a theory. Four major ideas of this movement are described: bifurcation theory, strange attractors, fractals, and the renormalization group. A framework for the future study of turbulent flows is proposed, based on the premises of the deterministic movement.

Chapman, G. T.↗

High-Fidelity Simulation of the Light-to-Dense Stratification Transient in the HiRJET Facility

Density stratification in a large enclosure is a crucial phenomenon to heat transfer and sustainable passive heat removal of a sodium fast reactor during reactivity transients. However, engineering turbulence models were identified to have unsatisfactory performance in predicting propagation of a stratified front. Yet, the scarcity of high-resolution data for stratification hampers the development of models. To explor e applications of leveraging direct numerical simulation (DNS) data to support turbulence model development, this work conducted DNS using NekRS to study a long stratification transient in the High-Resolution Jet (HiRJET) experimental facility. Here, this work considers an experiment run where light fluid is injected into a tank containing a denser fluid with a relative density difference of 1.5%. Formation of the stratified layer is identified as impingement of the buoyant jet promoting mixing of the two fluids. Based on the transient statistics, transport of the concentration can be characterized by regions with dominating effects of turbulent mixing, buoyant dissipation, and molecular diffusion, respectively, as moving away from the elevation of jet impingement. Concentration near the stratified front also exhibits oscillation at Brunt-Väisälä frequency. Preliminary validation of the simulation showed encouraging agreement of the concentration distribution with the reference experiment.

buoyant mixing↗

Numerical Simulations of Vortex Generator Vanes and Jets on a Flat Plate

Numerical simulations of a single low-profile vortex generator vane, which is only a small fraction of the boundary-layer thickness, and a vortex generating jet have been performed for flows over a flat plate. The numerical simulations were computed by solving the steady-state solution to the Reynolds-averaged Navier-Stokes equations. The vortex generating vane results were evaluated by comparing the strength and trajectory of the streamwise vortex to experimental particle image velocimetry measurements. From the numerical simulations of the vane case, it was observed that the Shear-Stress Transport (SST) turbulence model resulted in a better prediction of the streamwise peak vorticity and trajectory when compared to the Spalart-Allmaras (SA) turbulence model. It is shown in this investigation that the estimation of the turbulent eddy viscosity near the vortex core, for both the vane and jet simulations, was higher for the SA model when compared to the SST model. Even though the numerical simulations of the vortex generating vane were able to predict the trajectory of the stream-wise vortex, the initial magnitude and decay of the peak streamwise vorticity were significantly under predicted. A comparison of the positive circulation associated with the streamwise vortex showed that while the numerical simulations produced a more diffused vortex, the vortex strength compared very well to the experimental observations. A grid resolution study for the vortex generating vane was also performed showing that the diffusion of the vortex was not a result of insufficient grid resolution. Comparisons were also made between a fully modeled trapezoidal vane with finite thickness to a simply modeled rectangular thin vane. The comparisons showed that the simply modeled rectangular vane produced a streamwise vortex which had a strength and trajectory very similar to the fully modeled trapezoidal vane.

Allan, Brian G.↗

Ares I Vehicle Computed Turbulent Ascent Aerodynamic Data Development and Analysis

An overview of the computational ascent aerodynamic data development from a Reynolds-averaged Navier-Stokes flow solver for the Ares I vehicle design is presented. The computed results are assessed for grid and turbulence model effects and verified against the results obtained from other flow solvers. The numerical predictions are analyzed for the surface pressure, sectional line loads, longitudinal aerodynamic and rolling moment coefficients, and their trends with respect to angle of attack, Mach number, and the vehicle's roll angle. The results revealed that the solution development with the Spalart-Allmaras turbulence model was the most robust, stable, and efficient. These predictions were generally found to compare well with experimental data. Relative to wind-tunnel flow conditions, the results at flight Reynolds number showed the largest reduction of about 7% in the computed longitudinal aerodynamic coefficients. Protuberance size and relative position were found to have a significant effect on the vehicle's force and moment coefficients. The strake concept was proposed and shown to reduce the vehicle's maximum rolling moment coefficient Results and analyses have demonstrated the feasibility of the present numerical method, as an engineering tool, for predicting the external flow aerodynamic characteristics of this class of vehicle designs.

Aerodynamic characteristics↗

LAVA Simulations for the 3rd AIAA CFD High Lift Prediction Workshop with Body Fitted Grids

In response to the 3rd AIAA CFD High Lift Prediction Workshop, the workshop cases were analyzed using Reynolds-averaged Navier-Stokes flow solvers within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework. For the workshop cases the advantages and limitations of both overset-structured an unstructured polyhedral meshes were assessed. The workshop included 3 cases: a 2D airfoil validation case, a mesh convergence study using the High Lift Common Research Model, and a nacelle/pylon integration study using the JAXA (Japan Aerospace Exploration Agency) Standard Model. The 2D airfoil case from the workshop is used to verify the implementation of the Spalart-Allmaras turbulence model along with some of its variants within the solver. The High Lift Common Research Model case is used to assess solver performance and accuracy at varying mesh resolutions, as well as identify the minimum mesh fidelity required for LAVA on this class of problem. The JAXA Standard Model case is used to assess the solver's sensitivity to the turbulence model and to compare the structured and unstructured mesh paradigms. These workshop cases have helped establish best practices for high lift flow configurations for the LAVA solver.

LAVA↗

Low Reynolds number modeling of turbulent flows with and without wall transpiration

A full Reynolds-stress closure that is capable of describing the flow all the way to the wall is formulated. The closure is based on the conventional high Reynolds number form of the redistribution model, the inclusion of molecular diffusion, and a modified dissipation model to account for viscous effects near a wall. Two dissipation models are investigated along with two gradient diffusion and two redistribution models. Their respective effects on the calculated flow properties are assessed by comparing them with the data of fully developed turbulent flows and a developing pipe flow with wall transpiration. The near-wall behavior is very well predicted; however, the wall correction to the redistribution modeling is found to have little effect on the calculated results. The overall behavior of the fully developed turbulent flows is best described by a nonisotropic gradient diffusion model, a return-to-isotropy redistribution model, and a dissipation model that accounts for viscous behavior near a wall. This same closure also gives the best prediction of the axial pressure drop behavior along a pipe with a uniform wall suction. Furthermore, the near-wall behavior of such a flow is very well predicted by this closure.

So, Ronald M. C.↗