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,081 records · Page 60

CFD Validation Study of a Hypersonic Cone-Slice-Flap Configuration

Model validation is the process of determining the degree of accuracy between the real world and the model. The result of model validation can be used to either improve the model through calibration or quantify the model-form uncertainty. This work focuses on estimating the model-form uncertainty through an area metric for a hypersonic cone-slice-flap geometry configuration given the uncertainty in both the simulation and experimental data. A procedure that can give an accurate representation of the model-form uncertainty using a small number of high-fidelity runs is outlined. The work also assesses the impact of using different high-fidelity solvers and different turbulence models. The goal of performing this comparison is to provide a quantifiable measurement of the accuracy of each solver and turbulence model for this type of design. A low-fidelity analysis is also performed to get a model-form uncertainty for this type of analysis. The two high-fidelity CFD solvers used here are VULCAN-CFD and FUN3D, and the low-fidelity results comes from Cart3D.The experimental data comes from the 20-inch Mach 6 Tunnel located at NASA Langley Research Center. The work here shows that the models used tend to under-predict the drag and moment aerodynamic coefficients for this type of design leading to a model-form uncertainty estimate which can account for up to an 86% of the total uncertainty in the model predictions.

Laura M White↗

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The paper discusses the following: 1. No-cost Extension. The no-cost extension is required to complete the work on the unified model codes (both hydrodynamic and kinetic Monte Carlo) as described in the initial proposal and previous annual reports. 2. Scientific Accomplishments during the Report Period. We completed a comprehensive model of Alfvtn wave reflection that spans the full distance from the photosphere to the distant heliosphere. 3. Comparison of Accomplishments with Proposed Goals. The proposal contained two specific objectives for Year 3: (1) to complete the unified model code, and (2) to apply it to various kinds of coronal holes (and polar plumes within coronal holes). Although the anticipated route toward these two final goals has changed (see accomplishments 2a and 2b above), they remain the major milestones for the extended period of performance. Accomplishments la and IC were necessary prerequisites for the derivation of "physically relevant transport and mode-coupling terms" for the unified model codes (as stated in the proposal Year 3 goals). We have fulfilled the proposed "core work" to study 4 general types of physical processes; in previous years we studied turbulence, mode coupling (Le., non-WKB reflection), and kinetic wave damping, and accomplishment lb provides the fourth topic: nonlinear steepening.

Wagner, William↗

Mesh and model requirements for capturing deep-stall aerodynamics in low-Mach-number flows

Here, this work presents a comprehensive computational fluid dynamics investigation of the effects of grid resolution and turbulence-model choice for capturing the unsteady three-dimensional aerodynamic performance of NACA 0012 and 0021 airfoils, with specific focus on the deep-stall regime. At high angles of attack (α), wind turbine blades routinely experience vortex-induced vibrations, which can cause significant structural damages. Accurate predictions of post-stall aerodynamics can identify the frequencies at which such vibrations maybe triggered. In this context, the NACA 0012 airfoil simulations are conducted at a chord-based Reynolds number, Re c =2×10 6 , with the k-ω Shear-Stress Transport Reynolds-Averaged Navier-Stokes (RANS) and Improved Delayed Detached Eddy Simulation (IDDES) hybrid RANS-Large Eddy Simulation turbulence models. The effect of mesh resolution both in the wall-normal and spanwise directions is investigated. Only the IDDES model with a minimum spanwise resolution of 24 cells per chord length correctly predicts the aerodynamic forces. Spectral analysis shows the peak primary shedding frequency at α=30°, which signifies the end of the stall region. In the post-stall regime, both lift and drag frequencies drop asymptotically with increasing α. The Strouhal number, based on normalised chord length, remains nearly constant in this region. Based on this study, NACA 0021 airfoil runs are performed with IDDES for Re c =2.7×10 5 and 2.0×10 6 on the finest wall-normal mesh and three spanwise grids. Simulations conducted on the finer spanwise grids demonstrate grid independence and show good agreement with experiments. The effect of varying Rec on the airfoil frequency statistics is investigated. Additionally, comparison studies are presented to investigate the impact of airfoil thickness on the frequency content at Re c =2.0×10 6 . The results from the study provide guidance on the choice of mesh resolution with the IDDES model to accurately capture aerodynamic quantities for complex industrial applications.

17 WIND ENERGY↗

Duct flow nonuniformities: Effect of struts in SSME HGM II(+)

A numerical study, using the INS3D flow solver, of laminar and turbulent flow around a two dimensional strut, and three dimensional flow around a strut in an annulus is presented. A multi-block procedure was used to calculate two dimensional laminar flow around two struts in parallel, with each strut represented by one computational block. Single block calculations were performed for turbulent flow around a two dimensional strut, using a Baldwin-Lomax turbulence model to parameterize the turbulent shear stresses. A modified Baldwin-Lomax model was applied to the case of a three dimensional strut in an annulus. The results displayed the essential features of wing-body flows, including the presence of a horseshoe vortex system at the junction of the strut and the lower annulus surface. A similar system was observed at the upper annulus surface. The test geometries discussed were useful in developing the capability to perform multiblock calculations, and to simulate turbulent flow around obstructions located between curved walls. Both of these skills will be necessary to model the three dimensional flow in the strut assembly of the SSME. Work is now in progress on performing a three dimensional two block turbulent calculation of the flow in the turnaround duct (TAD) and strut/fuel bowl juncture region.

Burke, Roger↗

Investigation of third-order closure model of turbulence for the computation of incompressible flows in a channel with a backward-facing step

To predict the diffusion process of the Reynolds stresses in reattaching shear flows, the transport model for the triple-velocity products has been developed and tested for the computation of the flow in a channel with a backward-facing step. Upon comparison of the results of uuv, uvv, and vvv with those obtained by using existing algebraic correlations, it was shown that the present model improved the prediction of the triple-velocity products.

Amano, R. S.↗

Validation Assessment of Turbulent Reacting Flow Model Using The Area-Validation Metric on Medium-Scale Methanol Pool Fire Results

Accident analysis and ensuring power plant safety are pivotal in the nuclear energy sector. Significant strides have been achieved over the past few decades regarding fire protection and safety, primarily centered on design and regulatory compliance. Yet, after the Fukushima accident a decade ago, the imperative to enhance measures against fire, internal flooding, and power loss has intensified. Hence, a comprehensive, multilayered protection strategy against severe accidents is needed. Consequently, gaining a deeper insight into pool fires and their behavior through extensive validated data can greatly aid in improving these measures using advanced validation techniques. A model validation study was performed at Sandia National Laboratories in which a 30-cm diameter methanol pool fire was modeled using the SIERRA/Fuego turbulent reacting flow code. This validation study used a standard validation experiment to compare model results against, and conclusions have been published. The fire was modeled with a Large Eddy Simulation (LES) turbulence model with subgrid turbulent kinetic energy closure. Combustion was modeled using a strained laminar flamelet library approach. Radiative heat transfer was accounted for with a model utilizing the gray-gas approximation. In the present study, additional validation analysis is performed using the area validation metric (AVM). These activities are done on multiple datasets involving different variables and temporal/spatial ranges and intervals. In conclusion, the results provide insight into the use of the area validation metric on such temporally varying datasets and the importance of physics-aware use of the metric for proper analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Predicted Turbine Heat Transfer for a Range of Test Conditions

Comparisons are shown between predictions and experimental data for blade and endwall heat transfer. The comparisons of computational domain parisons are given for both vane and rotor geometries over an extensive range of Reynolds and Mach numbers. Comparisons are made with experimental data from a variety of sources. A number of turbulence models are available for predicting blade surface heat transfer, as well as aerodynamic performance. The results of an investigation to determine the turbulence model which gives the best agreement with experimental data over a wide range of test conditions are presented.

Boyle, R. J.↗

An experimental investigation of a supercritical airfoil at transonic speeds

Benchmark experimental data obtained in the two-dimensional, transonic flow field surrounding a supercritical airfoil are presented. Airfoil surface and tunnel wall pressure and LDV measurements are used to describe the flow on the model, above the wing and in the wake. Comparisons are made with calculations using the Reynolds-averaged Navier-Stokes equations. The results illustrate the performance of two turbulence models in both separated and attached flows. The largest differences between theory and experiment occurred in separated flows with the Johnson and King turbulence model providing the best estimates.

Mateer, George C.↗

Ocean Turbulence: Mesoscale Modeling in Isopycnal Coordinates the role of the Spectrum of Vertical Shear - Part 4

We study the tracer subgrid term in isopycnal coordinates, S(sub I). We employ two ingredients: the experimental data on vertical spectra of ocean turbulence measured by Gargett et al.(1981) and the stochastic approach recently developed by Dukowicz and Smith (1997). Our result confirms that S(sub I) is made of two parts: an advection and a diffusion term. However, the tracer bolus velocity u** consists of two terms u** = u(sub 1) + u(sub 2) while in the GM model there is only a term related to u(sub 1) which is shown to be: u(sub 1) = k(bar-q)(sup -1)(delta)(sub rho) where bar-q is the thickness weighted average potential vorticity, a result in agreement with the recent suggestions by Treguier et al. (1997), Lee et al. (1997) and Greatbatch (1998). The second component u(sub 2) IS new. We compute it in the geostrophic approximation using the Gargett et al. data (1981) on ocean vertical turbulence. We find that u(sub 2) much greater than u(sub 1) and that u(sub 2) is orthogonal to u(sub 1).

Canuto, V. M.↗

Bayesian parameter estimation and evaluation of the K -ω shear stress transport model for plane impinging jets

Numerical simulations with semi-empirical turbulence models are commonly used to model impinging jets, often used for cooling solid surfaces. In this work, the constants in the k-ω shear stress transport model in ANSYS FLUENT are calibrated to experimental velocity and heat transfer data for a plane turbulent impinging air jet to determine if Kennedy-O'Hagan calibration (Kennedy and O'Hagan 2001 J. R. Stat. Soc. B 63 425–64) can improve predictions of near-surface velocities and surface Nusselt numbers for similar flows. Impinging jets have been proposed to cool the target plates of the divertor in future magnetic fusion energy reactors, where simulations are used to estimate divertor performance. The flat-plate divertor (Wang et al 2009 Fusion Sci. Technol .56 1023–7) uses a plane jet of helium issuing from a B = 0.5 mm slot to cool a surface with radius of curvature of 44 B at a distance 4 B from the slot. Predictions from the calibrated numerical model are compared with independent experimental data at different flow conditions, as well as surface temperature data for a flat plate divertor test section. The contribution of this work is evaluation of the accuracy of a calibrated turbulence model for modest extrapolations in flow geometry and flow conditions for a plane impinging jet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SST ..kappa..-..omega.. Simulations of the Atmospheric Boundary Layer Including the Coriolis Effect

For large-scale structures in the atmospheric boundary layer (ABL), the Coriolis effect and near-wall behavior can have a meaningful impact. For example, both the Coriolis effect and blade boundary layer impact how muchpower wind farms produce. RANS simulations of the ABL typically use the ..kappa..-epsilon turbulence model, which has been developed to accurately capture the Coriolis effect but typically does not perform well near walls. The SST ..kappa..-..omega.. turbulence model performs well near walls but does not accurately capture the Coriolis effect. This work modifies SST ..kappa..-..omega.. to accurately model the Coriolis effect. We discuss the similarities and differences in how to modify ..kappa..-epsilon and SST ..kappa..-..omega.. for the Coriolis effect. Finally, we compare ..kappa..-epsilon and SST ..kappa..-..omega.. simulations of the ABL, including the Coriolis effect, with large eddy simulations and measurements.

atmospheric boundary layer↗

Application of an invariant second-order closure model to compressible turbulent shear layers

A second-order closure model for two-dimensional, compressible shear flows is investigated using the invariant modeling technique developed by Donaldson. The invariant model parameters were originally selected by comparison of model predictions with critical experimental data on basic incompressible flows. Additional modeling for compressible flows has been introduced. A number of different shear flows such as the free shear layer, flat plate boundary layer and a simulation of the mixing region of a chemical laser are computed. The model is consistent with first-order closure turbulence models for equilibrium flows and is further capable of predicting nonequilibrium flows that cannot be correctly solved by eddy viscosity models. The results are in generally good agreement with experimental measurements but suggest the need for inclusion of a Mach number dependent model of the pressure diffusion terms in order to adequately represent high speed compressible flows.

Varma, A. K.↗

Studying Turbulence Using Numerical Simulation Databases: Proceedings of the Summer Program - No. 7

The Seventh Summer Program of the Center for Turbulence Research took place in the four-week period, July 5 to July 31, 1998. This was the largest CTR Summer Program to date, involving thirty-six participants from the U. S. and nine other countries. Thirty-one Stanford and NASA-Ames staff members facilitated and contributed to most of the Summer projects. A new feature, and perhaps a preview of the future programs, was that many of the projects were executed on non-NASA computers. These included supercomputers located in Europe as well as those operated by the Departments of Defense and Energy in the United States. In addition, several simulation programs developed by the visiting participants at their home institutions were used. Another new feature was the prevalence of lap-top personal computers which were used by several participants to carry out some of the work that in the past were performed on desk-top workstations. We expect these trends to continue as computing power is enhanced and as more researchers (many of whom CTR alumni) use numerical simulations to study turbulent flows. CTR's main role continues to be in providing a forum for the study of turbulence for engineering analysis and in facilitating intellectual exchange among the leading researchers in the field. Once again the combustion group was the largest. Turbulent combustion has enjoyed remarkable progress in using simulations to address increasingly complex and practically more relevant questions. The combustion group's studies included such challenging topics as fuel evaporation, soot chemistry, and thermonuclear reactions. The latter study was one of three projects related to the Department of Energy's ASCI Program (www.llnl.gov/asci); the other two (rocket propulsion and fire safety) were carried out in the turbulence modeling group. The flow control and acoustics group demonstrated a successful application of the so-called evolution algorithms which actually led to a previously unknown forcing strategy for jets yielding increased spreading rate. A very efficient algorithm for flow in complex geometries with moving boundaries based on the immersed boundary forcing technique was tested with very encouraging results. Also a new strategy for the destruction of aircraft trailing vortices was introduced and tested. The Reynolds Averaged Modeling (RANS) group demonstrated that the elliptic relaxation concept for RANS calculations is also applicable to transonic flows with shocks; however, prediction of laminar/turbulent transition remains an important pacing item. A large fraction of the LES effort was devoted to the development and testing of a new algorithmic procedure (as opposed to phenomenological model) for subgrid scale modeling based on regularized de-filtering of the flow variables. This appears to be a very promising approach, and a significant effort is currently underway to assess its robustness in high Reynolds number flows and in conjunction with numerical methods for complex flows. As part of the Summer Program two review tutorials were given on Turbulent structures in hydrocarbon pool fires (Sheldon Tieszen), and Turbulent combustion modeling: from RANS to LES via DNS (Luc Vervisch); and two seminars entitled Assessment of turbulence models for engineering applications (Paul Durbin) and Subgrid-scale modeling for non-premixed, turbulent reacting flows (James Riley) were presented. A number of colleagues from universities, government agencies, and industry attended the final presentations of the participants on July 31 and participated in the discussions. There are twenty-six papers in this volume grouped in five areas. Each group is preceded with an overview by its coordinator.

Source record↗

Hypersonic Shock/Boundary-Layer Interaction Database

Turbulence modeling is generally recognized as the major problem obstructing further advances in computational fluid dynamics (CFD). A closed solution of the governing Navier-Stokes equations for turbulent flows of practical consequence is still far beyond grasp. At the same time, the simplified models of turbulence which are used to achieve closure of the Navier-Stokes equations are known to be rigorously incorrect. While these models serve a definite purpose, they are inadequate for the general prediction of hypersonic viscous/inviscid interactions, mixing problems, chemical nonequilibria, and a range of other phenomena which must be predicted in order to design a hypersonic vehicle computationally. Due to the complexity of turbulence, useful new turbulence models are synthesized only when great expertise is brought to bear and considerable intellectual energy is expended. Although this process is fundamentally theoretical, crucial guidance may be gained from carefully-executed basic experiments. Following the birth of a new model, its testing and validation once again demand comparisons with data of unimpeachable quality. This report concerns these issues which arise from the experimental aspects of hypersonic modeling and represents the results of the first phase of an effort to develop compressible turbulence models.

Settles, G. S.↗

Computation of crossing shock/turbulent boundary layer interaction at Mach 8.3

A three-dimensional (3D) hypersonic crossing shock wave/turbulent boundary-layer interaction is examined numerically at Mach 8.3. The test geometry consists of a pair of opposing sharp fins of angle alpha = 15 deg, mounted on a flat plate. Two theoretical models are evaluated. The full (3D) Reynolds-averaged Navier-Stokes equations are solved using the Baldwin-Lomax and the Rodi (modified k-epsilon) turbulence models. Computed results for both cases show good agreement with experiment for flat plate surface pressure and for flowfield profiles of pitot pressure and yaw angle, indicating that the flowfield is primarily rotational and inviscid. Fair to poor agreement is obtained for surface heat transfer, indicating a need for more accurate turbulence models. The overall flowfield structure is similar to that observed in previous crossing shock interaction studies.

Narayanswami, N.↗

High speed nozzles task

Supersonic cruise exhaust nozzles for advanced applications are optimized for a high nozzle pressure ratio (NPR) at design supersonic cruise Mach number and altitude. The performance of these nozzles with large expansion ratios are severely degraded for operations at subsonic speeds near sea level for NPR significantly less than the design values. The prediction of over-expanded 2DCD nozzles performance is critical to evaluating the internal losses and to the optimization of the integrated vehicle and propulsion system performance. The reported research work was aimed at validating and assessing existing computational methods and turbulence models for predicting the flow characteristics and nozzle performance at over-expanded conditions. Flow simulations in 2DCD nozzles were performed using five different turbulence models. The results are compared with the experimental data for the wall pressure distribution and thrust and flow coefficients at over-expanded static conditions.

Hamed, Awatef↗

Reynolds-Averaged Navier-Stokes Simulation of a 2D Circulation Control Wind Tunnel Experiment

Numerical simulations are performed using a Reynolds-averaged Navier-Stokes (RANS) flow solver for a circulation control airfoil. 2D and 3D simulation results are compared to a circulation control wind tunnel test conducted at the NASA Langley Basic Aerodynamics Research Tunnel (BART). The RANS simulations are compared to a low blowing case with a jet momentum coefficient, C(sub u), of 0:047 and a higher blowing case of 0.115. Three dimensional simulations of the model and tunnel walls show wall effects on the lift and airfoil surface pressures. These wall effects include a 4% decrease of the midspan sectional lift for the C(sub u) 0.115 blowing condition. Simulations comparing the performance of the Spalart Allmaras (SA) and Shear Stress Transport (SST) turbulence models are also made, showing the SST model compares best to the experimental data. A Rotational/Curvature Correction (RCC) to the turbulence model is also evaluated demonstrating an improvement in the CFD predictions.

Allan, Brian G.↗

Numerical Simulations For the F-16XL Aircraft Configuration

Numerical simulations of flow around the F-16XL are presented as a contribution to the Cranked Arrow Wing Aerodynamic Project International II (CAWAPI-II). The NASA Tetrahedral Unstructured Software System (TetrUSS) is used to perform numerical simulations. This CFD suite, developed and maintained by NASA Langley Research Center, includes an unstructured grid generation program called VGRID, a postprocessor named POSTGRID, and the flow solver USM3D. The CRISP CFD package is utilized to provide error estimates and grid adaption for verification of USM3D results. A subsonic high angle-of-attack case flight condition (FC) 25 is computed and analyzed. Three turbulence models are used in the calculations: the one-equation Spalart-Allmaras (SA), the two-equation shear stress transport (SST) and the kε turbulence models. Computational results, and surface static pressure profiles are presented and compared with flight data. Solution verification is performed using formal grid refinement studies, the solution of Error Transport Equations, and adaptive mesh refinement. The current study shows that the USM3D solver coupled with CRISP CFD can be used in an engineering environment in predicting vortex-flow physics on a complex configuration at flight Reynolds numbers.

Elmiligui, Alaa A.↗