Development of a novel high-order flux correction method on strand grids is presented. The method uses a combination of flux correction in the unstructured plane and summation-by-parts operators in the strand direction to achieve high-fidelity solutions. Low-order truncation errors are cancelled with accurate flux and solution gradients in the flux correction method, thereby achieving a formal order of accuracy of 3, although higher orders are often obtained, especially for highly viscous flows. In this work, the scheme is extended to high-Reynolds number computations in both two and three dimensions. Turbulence closure is achieved with a robust version of the Spalart-Allmaras turbulence model that accommodates negative values of the turbulence working variable, and the Menter SST turbulence model, which blends the k-ε and k-ω turbulence models for better accuracy. A major advantage of this high-order formulation is the ability to implement traditional finite volume-like limiters to cleanly capture shocked and discontinuous flow. In this work, this approach is explored via a symmetric limited positive (SLIP) limiter. Extensive verification and validation is conducted in two and three dimensions to determine the accuracy and fidelity of the scheme for a number of different cases. Verification studies show that the scheme achieves better than third order accuracy for low and high-Reynolds number flow. Cost studies show that in three-dimensions, the third-order flux correction scheme requires only 30% more walltime than a traditional second-order scheme on strand grids to achieve the same level of convergence. In order to overcome meshing issues at sharp corners and other small-scale features, a unique approach to traditional geometry, coined "asymptotic geometry," is explored. Asymptotic geometry is achieved by filtering out small-scale features in a level set domain through min/max flow. This approach is combined with a curvature based strand shortening strategy in order to qualitatively improve strand grid mesh quality.

Standard eddy-viscosity models lack curvature and system rotation sensitized terms in their formulation. Hence they fail to capture the effects of curvature and system rotation on turbulence anisotropy. As part of this effort, an algebraic expression for a characteristic rotation term is developed and tuned with the help of rotating homogeneous shear flow. This formulation is primarily based upon the rotation and curvature sensitized eddy-viscosity coefficient developed by York et al. (2009). A new scalar transport equation loosely based on Durbin’s wall normal turbulent velocity scale (Durbin, 1991) is introduced to account for the modification in turbulence structure due to system rotation and curvature effects. The added transport equation also introduces history effects and stability in the solution with small increase in computational cost. The eddy-viscosity is redefined based on new turbulent velocity scale and hence the effects of rotation and streamline curvature are introduced into the mean momentum equation. A number of canonical test cases with significant curvature and rotation effects along with a cyclone flow, a representative of complex industrial flows, are considered for model validation. Hybrid modeling framework combines the strength of RANS in boundary layers and LES in separated shear layers to alleviate the weaknesses of RANS and limitations of LES model in some complex flows. A recently proposed hybrid RANS-LES modeling framework uses a weighing parameter that dynamically determines the RANS and LES regions based on solution statistics. The hybrid modeling methodology is implemented on a normal jet in crossflow, and a film cooling case for the purpose of model validation and evaluation. The final goal of the proposed effort is to combine advanced RANS modeling capability with LES using the new hybrid modeling framework. Specifically, the curvature and rotation sensitive RANS model developed here is coupled with commonly used LES models to produce a novel model for complex turbulent flows with the potential to improve accuracy of CFD predictions (versus existing RANS models) as well as significantly reduce the computational expense (versus existing LES models). Performance of the model form hence developed is evaluated on a cyclone flow case.

Modeling of complex flows involving the combined effects of flow transition and streamline curvature using two advanced turbulence models, one in the Reynolds-averaged Navier-Stokes (RANS) category and the other in the hybrid RANS-Large eddy simulation (LES) category is considered in this research effort. In the first part of the research, a new scalar eddy-viscosity model (EVM) is proposed, designed to exhibit physically correct responses to flow transition, streamline curvature, and system rotation effects. The four equation model developed herein is a curvature-sensitized version of a commercially available three-equation transition-sensitive model. The physical effects of rotation and curvature (RC) enter the model through the added transport equation, analogous to a transverse turbulent velocity scale. The eddy-viscosity has been redefined such that the proposed model is constrained to reduce to the original transition-sensitive model definition in nonrotating flows or in regions with negligible RC effects. In the second part of the research, the developed four-equation model is combined with a LES technique using a new hybrid modeling framework, dynamic hybrid RANS-LES. The new framework is highly generalized, allowing coupling of any desired LES model with any given RANS model and addresses several deficiencies inherent in most current hybrid models. In the present research effort, the DHRL model comprises of the proposed four-equation model for RANS component and the MILES scheme for LES component. Both the models were implemented into a commercial computational fluid dynamics (CFD) solver and tested on a number of engineering and generic flow problems. Results from both the RANS and hybrid models show successful resolution of the combined effects of transition and curvature with reasonable engineering accuracy, and for only a small increase in computational cost. In addition, results from the hybrid model indicate significant levels of turbulent fluctuations in the flowfield, improved accuracy compared to RANS models predictions, and are obtained at a significant reduction of computational cost compared to full LES models. The results suggest that the advanced turbulence modeling techniques presented in this research effort have potential as practical tools for solving low/high Re flows over blunt/curved bodies for the prediction of transition and RC effects.

No description available.

In recent years, airports have experienced increasing airport congestion, partially due to the hub-and-spoke model on which airline operations are based. Current airline operations utilize large airports, focusing traffic to a small number of airports. One way to relieve such congestion is to transition to a more accessible and efficient point-to-point operation, which utilizes a large web of smaller airports. This expansion to regional airports propagates the need for next-generation low-noise aircraft with short take-off and landing capabilities. NASA has attacked this problem with a high-lift, low-noise concept dubbed the Cruise Efficient Short Take-Off and Landing (CESTOL) aircraft. The goal of the CESTOL project is to produce aircraft designs that can further expand the air travel industry to currently untapped regional airports. One method of obtaining a large lifting capability with low noise production is to utilize circulation control (CC) technology. CC is an active flow control approach that makes use of the Coanda effect. A high speed jet of air is blown over a wing flap and/or the leading edge of the wing, which entrains the freestream flow and effectively increases circulation around the wing. A promising tool for predicting CESTOL aircraft performance is computational fluid dynamics (CFD,) due to the relatively low cost and easy implementation in the design process. However, the unique flows that CC introduces are not well understood, and traditional turbulence modeling does not correctly resolve these complex flows (including high speed jet flow, complex shear flows and mixing phenomena, streamline curvature, and other challenging flow phenomena). The recent derivation of the v2-f turbulence model shows theoretical promise in increasing the accuracy of CFD predictions for CC flows, but this has not yet been assessed in great detail. This paper presents a methodical verification of several variations on the v2-f turbulence model. These models are verified using simple, well-understood flows. Results for CC flows are compared to those obtained with more traditional turbulence modeling techniques (including the Spalart-Allmaras, k-ε, and k-ω turbulence models). Wherever possible, computed results are compared to experimental data and more accurate numerical methods. Results indicate that the v2-f turbulence models predict some aspects of circulation control flow fields quite well, in particular the lift coefficient. The linear v2-f, nonlinear v2-f, and nonlinear v2-f-cc turbulence models have generated lift coefficients within 19%, 14%, and -26%, respectively of experimental values, whereas the Spalart-Allmaras, k-ε, and k-ω turbulence models produce errors as high as 85%, 36%, and 39%, respectively. The predicted stagnation points and pressure coefficient distributions match experimental data roughly as well as standard turbulence models do, though the modeling of these aspects of the flow do show some room for improvement. The nonlinear v2-f-cc turbulence model shows very non-physical skin friction coefficient profiles, pressure coefficient profiles, and stagnation points, indicating that the streamline curvature correction terms need attention. Regardless of the source of the discrepancies, the v2-f turbulence models show promise in the modeling of circulation control flow fields, but are not quite ready for application in the design of circulation control aircraft.

There is increasing interest in design methods and performance prediction for turbine engines operating at low Reynolds numbers. In this regime, boundary layer separation may be more likely to occur in the turbine flow passages. For accurate CFD predictions of the flow, correct modeling of laminar-turbulent boundary layer transition is essential to capture the details of the flow. To investigate possible improvements in model fidelity, both two-dimensional and three-dimensional CFD models were created for the flow over several low pressure turbine blade designs. A new three-equation eddy-viscosity type turbulent transitional flow model originally developed by Walters and Leylek was employed for the current RANS CFD calculations. Flows over three low pressure turbine blade airfoils with different aerodynamic characteristics were simulated over a Reynolds number range of 15,000-100,000, and predictions were compared to experiments. The turbulent transitional flow model sensitivity to inlet turbulent flow parameters showed a dependence on free-stream turbulence intensity and turbulent length scale. Using the total pressure loss coefficient as a measurement of aerodynamic performance, the Walters and Leylek transitional flow model produced adequate prediction of the Reynolds number performance in the Lightly Loaded blade. Furthermore, the correct qualitative flow response to separated shear layers was observed for the Highly Loaded blade. The vortex shedding produced by the separated flow was largely two-dimensional with small spanwise variations in the separation region. The blade loading and separation location was sufficiently predicted for the Aft-Loaded L1A blade flowfield. Investigations of the unsteady flowfield of the Aft-Loaded L1A blade showed the shear layer produced a large separation region on the suction surface. This separation region was located more downstream and significantly reduced in size when impinged upon by the upstream wakes, thus improving the aerodynamic performance consistent with experiments. For all cases investigated, the Walters and Leylek transitional flow model was judged to be sufficient for understanding the separation and transition characteristics, and superior to other widely-used turbulence models in accuracy of describing the details of the transitional and separated flow. This research characterized and assessed a new model for low Reynolds number turbine aerodynamic flow prediction and design improvement. / Ph. D.

In this work, a three-dimensional turbulent boundary layer experiment was set up with alternating stream-wise and span-wise pressure gradients. The pressure gradients are generated as a result of the test section wavy side wall shape. Each side had six sine waves with a trough to peak magnitude to wavelength ratio of 0.25. Boundary layer control was used so that the flow over the side walls remains attached. The mean flow velocity components, static and total pressures were measured at six plane along the stream-wise direction. The alternating mean span-wise and stream-wise pressure gradients created alternating stream-wise and span-wise vorticity fluxes, respectively, along the test section. As the flow developed downstream the vorticity created at the tunnel floor and ceiling diffused away from the wall. The vorticity components in the stream-wise and span-wise directions are strengthened due to stretching and tilting terms in the vorticity transport equations. The positive-z half of the test section contains large areas that generate positive vorticity flux in the trough region and smaller areas generating negative vorticity around the wave peak. The opposite is true for the negative-z half of the test-section. This results in a large positive stream-wise vorticity in the positive-z half and negative stream-wise vorticity in the negative-z half of the test-section. The smaller regions of opposite sign vorticity in each half tend to mix the flow such that as they diffuse away from the wall, the turbulent stresses are more uniform. Turbulent fluctuating velocity components were measured using Laser Doppler Velocimetery. Mean velocities as well as Reynolds stresses and triple velocity component correlations were measured at thirty stations along the last wave in the test section. Profiles at the center of the test section showed three dimensionality, but exhibited high turbulence intensities in the outer layer. Profiles off the test section centerline are highly three dimensional with multiple peaks in the normal stress profiles. The flow also reaches a state where all the normal stresses have equal magnitudes while the shear stresses are non-zero. Flow angles, flow gradient angles and shear stress angles show very large differences between wall values and outer layer vlaues. The shear stress angle lagged the flow gradient angle indicating non-equilibrium. A turbulent kinetic energy transport budget is performed for all profiles and the turbulence kinetic energy dissipation rate is estimated. Spectral measurements were also made and an independent estimate of the kinetic energy dissipation rate is made. These estimates agree very well with those estimates made by balancing the turbulence kinetic energy transport equation. Multiple turbulent diffusion models are compared to measured quantities. The models varied in agreement with experimental data. However, fair agreement with turbulence kinetic energy turbulent diffusion is observed. A model for the dissipation rate tensor anisotropy is used to extract estimates of the pressure-strain tensor from the Reynolds stress transport equations. The pressure-strain estimates are compared with some of the models in the literature. The comparison showed poor agreement with estimated pressure-strain values extracted from experimental data. A tentative model for the turbulent Reynolds shear stress angle is developed that captures the shear stress angle near wall behavior to a very good extent. The model contains one constant that is related to mean flow variables. However, the developed expression needs modification so that the prediction is improved along the entire boundary layer thickness. / Ph. D.

No description available.

No description available.

No description available.