Recalibration of the Shear Stress Transport Model to Improve Calculation of Shock Separated Flows
No abstract available
SEARCH · Engineering Papers
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.
No abstract available
The performance of the new LagRST models is compared with curvature corrected eddyviscosity RANS models. The Reynolds-stress models provide very good predictions of thisow without any additional curvature corrections. Along with comparison of predictionsof the vortex itself, some of the other measurements included in the oweld are comparedwith the predictions in the wake and over the wing. Dierences in predictions between thevarious RANS models in the wake region are small - not substantial, and all predictionsare in line with experimental measurements from 0.005c to 0.67c downstream of the trailingedge, at odds with what self-similar wake analysis would predict.
Explore the source record for details and available documents.
A full Reynolds stress turbulence model is applied to calculate the flow in various channels with complex cross-section. The model is shown to be more robust and to produce more consistent results than algebraic stress models. The relationship between various versions of the latter is provided. It is shown that differential or algebraic stress models derived from a direct application of a popular linear approximation for the pressure-strain terms in the Reynolds stress equations will always underpredict the secondary motion, and hence the accompanying distortion of other flow properties. A simple approximation which leads to better predictions is presented and related to previous practice.
A variational analysis of the stress state in microcracked cross-ply laminates has been used to investigate the phenomenon of curved microcracking in /(S)/90n/s laminates. Previous investigators proposed that the initiation and orientation of curved microcracks are controlled by local maxima and stress trajectories of the principal stresses. We have implemented a principal stress model using a variational mechanics stress analysis and we were able to make predictions about curved microcracks. The predictions agree well with experimental observations and therefore support the assertion that the variational analysis gives an accurate stress state that is useful for modeling the microcracking properties of cross-ply laminates. An important prediction about curved microcracks is that they are a late stage of microcracking damage. They occur only when the crack density of straight microcracks exceeds the critical crack density for curved microcracking. The predicted critical crack density for curved microcracking agrees well with experimental observations.
The predictive accuracy of typical hybrid RANS/LES models is limited by several shortcomings, such as modeled-stress depletion and a dependence on scalar grid measures. The "active model-split" (AMS) hybrid RANS/LES model was developed by Haering, Oliver, and Moser to address these shortcomings. In this model, the mean and fluctuating portions of the stress are modeled separately; this allows for a consistent treatment of the mean flow even in the presence of resolved fluctuations. Part of this hybrid model is an active forcing, which generates turbulent fluctuations in regions that support LES-type fluctuations. In this presentation, some recent results are shown using the active model-split for shock-separated flows. The effects of the active forcing are examined, including the effects on anisotropy and acoustics. Superior predictive accuracy over RANS is shown for some quantities of interest, such as the shock location.
The Brief Report demonstrates the relationship of eddy-viscosity models to subgrid stresses. A formula that determines the relative importance of the cross stress and the Reynolds stress for the net eddy-damping and backscatter contributions is derived. The cross-stress term with sharp-cut filtering is identified as an important quantity to model. These concepts could prove useful as a basis for constructing specific models for the Reynolds and cross stresses.
A three dimensional full Navier-Stokes (FNS) code is used to simulate complex turbomachinery flows. The code incorporates an explicit multistep scheme and solves a conservative form of the density averaged continuity, momentum, and energy equations. A compressible low Reynolds number form of the k-epsilon turbulence model, and a q-omega model and an algebraic Reynolds stress model have been incorporated in a fully coupled manner to approximate Reynolds stresses. The code is used to predict the viscous flow field in a backswept transonic centrifugal compressor for which laser two focus data is available. The code is also used to simulate the tip clearance flow in a cascade. The code has been extended to include unsteady Euler solutions for predicting the unsteady flow through a cascade due to incoming wakes, simulating rotor-stator interactions.
Reynolds stress closure models based on the recursion renormalization group theory are developed for the prediction of turbulent separated flows. The proposed model uses a finite wavenumber truncation scheme to account for the spectral distribution of energy. In particular, the model incorporates effects of both local and nonlocal interactions. The nonlocal interactions are shown to yield a contribution identical to that from the epsilon-renormalization group (RNG), while the local interactions introduce higher order dispersive effects. A formal analysis of the model is presented and its ability to accurately predict separated flows is analyzed from a combined theoretical and computational stand point. Turbulent flow past a backward facing step is chosen as a test case and the results obtained based on detailed computations demonstrate that the proposed recursion -RNG model with finite cut-off wavenumber can yield very good predictions for the backstep problem.
The use of computational fluid dynamics in scramjet engine component development is widespread in the existing literature. Unfortunately, the quantification of model-form uncertainties is rarely addressed with anything other than sensitivity studies, requiring that the computational results be intimately tied to and calibrated against existing test data. This practice must be replaced with a formal uncertainty quantification process for computational fluid dynamics to play an expanded role in the system design, development, and flight certification process. Due to ground test facility limitations, this expanded role is believed to be a requirement by some in the test and evaluation community if scramjet engines are to be given serious consideration as a viable propulsion device. An effort has been initiated at the NASA Langley Research Center to validate several turbulence closure models used for Reynolds-averaged simulations of scramjet isolator flows. The turbulence models considered were the Menter BSL, Menter SST, Wilcox 1998, Wilcox 2006, and the Gatski-Speziale explicit algebraic Reynolds stress models. The simulations were carried out using the VULCAN computational fluid dynamics package developed at the NASA Langley Research Center. A procedure to quantify the numerical errors was developed to account for discretization errors in the validation process. This procedure utilized the grid convergence index defined by Roache as a bounding estimate for the numerical error. The validation data was collected from a mechanically back-pressured constant area (1 2 inch) isolator model with an isolator entrance Mach number of 2.5. As expected, the model-form uncertainty was substantial for the shock-dominated, massively separated flowfield within the isolator as evidenced by a 6 duct height variation in shock train length depending on the turbulence model employed. Generally speaking, the turbulence models that did not include an explicit stress limiter more closely matched the measured surface pressures. This observation is somewhat surprising, given that stress-limiting models have generally been developed to better predict shock-separated flows. All of the models considered also failed to properly predict the shape and extent of the separated flow region caused by the shock boundary layer interactions. However, the best performing models were able to predict the isolator shock train length (an important metric for isolator operability margin) to within 1 isolator duct height.
The determination of the stress distribution in the inner lap of double-lap, double-bolt joints using photoelastic models of the joint is discussed. The principal idea is to fabricate the inner lap of a photoelastic material and to use a photoelastically sensitive material for the two outer laps. With this setup, polarized light transmitted through the stressed model responds principally to the stressed inner lap. The model geometry, the procedures for making and testing the model, and test results are described.
The objective here is to test the Dynamic Localization (DL) model in a wall-bounded channel flow for numerical stability and accuracy of results. Algebraic stress models suggest that the model for the residual subgrid-scale (SGS) Reynolds stress and scalar flux should generally have terms comprising most of the unique products of the resolved strain (S) and rotation (R) tensors with S and the resolved scalar gradient. The standard dynamic SGS model uses a simple (Smagorinsky) base model for the residual Reynolds stress, which is made proportional to S, and down-gradient base models for residual scalar fluxes; these correspond to the lowest, 'first-order' terms in algebraic stress models. Temporal scaling terms in these base models are formed from the magnitude of the resolved strain rate. While this is appropriate for simple shear flows, it may not be appropriate for more complicated flows (relevant to geophysical and astrophysical problems) that include any combination of shear, rotation, buoyancy, etc. On the other hand, the coefficient in the dynamic SGS model readily adjusts itself to different flow conditions and may adequately take account of these effects without the need for more complicated base models. Cabot (1993) has begun to test the dynamic SGS model in buoyant flows (Rayleigh-Benard and internally heated convection) with and without buoyancy terms explicitly included in the scaling terms of the base model; no great differences were found in large eddy simulation (LES) results for the different base model scalings. The second objective in this work is to test base models with additional, 'second-order' terms (e.g., S(sup 2) and RS for the residual Reynolds stress). These terms have been found to improve large-scale flow predictions by kappa-epsilon models in the presence of rotation and shear. Second-order base models will be tested here in the LES of channel flow with and without solid-body rotation and compared with results from the standard first-order base models to determine if there are significant differences or improvements in results that would warrant the added complexity of the second-order base models.
Computations, employing several turbulence models, are compared with a series of attached supersonic turbulent boundary-layer experiments over an extensive range of Reynolds numbers (11.7 x 10 to the 6th to 314 x 10 to the 6th). These experiments included measurements of surface pressure and skin friction for adverse pressure gradients ranging in magnitude from those of previous investigations to an order of magnitude greater. The turbulence models evaluated include algebraic and two-equation eddy-viscosity models and two full Reynolds stress models. In general all the models tested performed well independently of the magnitude of the pressure gradient or Reynolds number and could predict the measured skin friction for most cases with sufficient accuracy for engineering purposes.
The wakes of rotor blades are turbulent, three-dimensional, and are subjected to curvature and rotation effects. The objective of this study is to predict the development of such wakes and compare the predictions with the existing data. A finite difference procedure is employed in the numerical analysis of the wake utilizing the continuity, momentum, and turbulence closure equations in the rotating curvilinear and non-orthogonal coordinate system. The turbulence closure is affected by the two equations model and the simplified Reynolds stress model. The predictions from both the models agree well with the mean velocity and the turbulence data.
The space marching code was modified in order to be able to predict the flow field inside a rotor passage, including the blade and hub wall boundary layers. The basic changes incorporated are modifications of the equations so that the code can handle three dimensional configurations with changes in the radial direction (for example changes in stagger angle, blade camber and thickness), extensions and modifications in order to implement a physically realistic turbulence model such as a k sigma model and an algebraic Reynolds stress model.
The compressive flow strength-strain rate behavior of the oxide-dispersion-strengthened alloy MA 6000E has been studied in the temperature range 1144-1366 K, with strain rates ranging from 2.1 x 10 to the -5th to 2.1 x 10 to the -7th per s. It is found that the inherent strength of the alloy is essentially the same in all test directions and that the low strength observed in tensile tests results from the inability of grain boundaries to support high tensile stresses. The failure of MA 6000E under high-temperature, slow plastic flow conditions is shown to be the result of concentrated slip. Slow plastic deformation in MA 6000E can be described by a threshold stress model of creep where threshold stresses are calculated from relatively fast testing procedures and the effective stress exponent for creep is assumed to be 3.5.
The focus was on the use of databases obtained from direct numerical simulations of turbulent flows, for study of turbulence physics and modeling. Topics addressed included: stochastic decomposition/chaos/bifurcation; two-point closure (or k-space) modeling; scalar transport/reacting flows; Reynolds stress modeling; and structure of turbulent boundary layers.
The effect of free-stream turbulence level on stall conditions in plane diffusers is investigated with a finite-volume numerical method which utilizes a k-epsilon turbulence model or a differential Reynolds stress model for closure. Computed results show qualitative agreement with measured data. The ability to eliminate stall and thereby increase pressure recovery in diffusers at large included angles by utilizing high free-stream turbulence is correctly predicted.