Engineering topics
Publications and source records attributed to Weinberg, B. C..
Work performed in the numerical simulation of turbine passage flows via a Navier-Stokes approach is discussed. Both laminar and turbulent simulations in both two and three dimensions are discussed. An outline of the approach, background, and an overview of the results are given.
The existence of large scale coherent structures in turbulent shear flows has been well documented. Discrepancies between experimental and computational data suggest a necessity to understand the roles they play in mass and momentum transport. Using conditional sampling and averaging on coincident two-component velocity and concentration velocity experimental data for swirling and nonswirling coaxial jets, triggers for identifying the structures were examined. Concentration fluctuation was found to be an adequate trigger or indicator for the concentration-velocity data, but no suitable detector was located for the two-component velocity data. The large scale structures are found in the region where the largest discrepancies exist between model and experiment. The traditional gradient transport model does not fit in this region as a result of these structures. The large scale motion was found to be responsible for a large percentage of the axial mass transport. The large scale structures were found to convect downstream at approximately the mean velocity of the overall flow in the axial direction. The radial mean velocity of the structures was found to be substantially greater than that of the overall flow.
An investigation of the influence of large scale structures on the flow development for coaxial jets with sudden expansion (with and without swirl) is presented. Both an experimental study and numerical predictions were performed for a configuration corresponding to that considered by Johnson and Bennett and Roback and Johnson. The effects of large scale structures on the swirling and nonswirling flows are documented, in particular their influence on turbulence modeling and the numerical simulation. The ensemble-averaged, time-dependent Navier-Stokes equations are solved by an LBI procedure to predict the turbulent flow field. Effects of artificial dissipation and placement of the upstream boundary in the numerical computation are also discussed.
The application of the time-dependent ensemble-averaged Navier-Stokes equations to transonic turbine cascade flow fields was examined. In particular, efforts focused on an assessment of the procedure in conjunction with a suitable turbulence model to calculate steady turbine flow fields using an O-type coordinate system. Three cascade configurations were considered. Comparisons were made between the predicted and measured surface pressures and heat transfer distributions wherever available. In general, the pressure predictions were in good agreement with the data. Heat transfer calculations also showed good agreement when an empirical transition model was used. However, further work in the development of laminar-turbulent transitional models is indicated. The calculations showed most of the known features associated with turbine cascade flow fields. These results indicate the ability of the Navier-Stokes analysis to predict, in reasonable amounts of computation time, the surface pressure distribution, heat transfer rates, and viscous flow development for turbine cascades operating at realistic conditions.
The multidimensional, ensemble-averaged, compressible, time-dependent Navier-Stokes equations have been used to study the turbulent flow field in two and three-dimensional turbine cascades. The viscous regions of the flow were resolved and non-slip boundary conditions were utilized on solid surfaces. The calculations were performed in a constructive 'O'-type grid which allows representation of the blade rounded trailing edge. Converged solutions were obtained in relatively few time steps (about 80-150) and comparisons for both surface pressure and heat transfer showed good agreement with data. The three-dimensional turbine cascade calculation showed many of the expected flow-field features.
An efficient computational procedure for solving three-dimensional unsteady turbulent flows is described. The consistently split Linearized Block Implicit (LBI) scheme is used in conjunction with the QR operator scheme to solve an approximate form of the Navier-Stokes equations in generalized nonorthogonal coordinates employing physical velocity components. As a demonstration calculation the turbulent oscillating flow over a flat plate corresponding to the experiment of Karlsson is considered in both two and three dimensions. New inflow boundary conditions are proposed which yield physically plausible solutions near the upstream boundary. The results obtained agree both qualitatively and quantitatively with Karlsson's data and shed new light on the controversy concerning the interpretation of the skin friction phase angle as a function of reduced frequency.
A Navier-Stokes analysis employing the time-dependent Linearized Block Implicit scheme (LBI) was applied to two-dimensional and three-dimensional transonic turbulent cascade flows. In general, the geometrical configuration of the turbine blade impacts both the grid construction procedure and the implementation of the numerical algorithm. Since modern turbine blades of interest are characterized by very blunt leading edges, rounded trailing edges and high stacking angles, a robust grid construction procedure is required that can accommodate the severe body shape while resolving regions of large flow gradients. A constructive O-type grid generation technique, suitable for cascades with rounded trailing edges, was developed and used to construct the C3X turbine cascade coordinate grid. Two-dimensional calculations were performed employing the Navier-Stokes procedure for the C3X turbine cascade, and the predicted pressure coefficients and heat transfer rates were compared with the experimental data. Three-dimensional Navier-Stokes calculations were also performed.
A three-dimensional solution procedure for the approximate form of the Navier-Stokes equation was exercised in the two- and three-dimensional modes to compute the unsteady turbulent boundary layer on a flat plate corresponding to the data of Karlsson. The procedure is based on the use of a consistently split Linearized Block Implicit technique in conjunction with a QR operator scheme. New time-dependent upstream boundary conditions were developed that yielded realistic solutions for the interior in the vicinity of the upstream boundary. Comparisons of the computation employing these boundary conditions with the data indicate that both qualitative and quantitative agreement was obtained for the mean velocity and the in phase and out of phase components of the first harmonic of the velocity. In addition, the calculation gave results for the skin friction phase angle that had expected physical behavior for large distances downstream of the inflow boundary. For the three-dimensional case, the two-dimensional data of Karlsson was considered, but in a coordinate system skewed at 45 deg to the free stream direction. The results of the calculations were in excellent agreement with the data and the two-dimensional computations.
A numerical scheme is developed for solving the time dependent, three dimensional compressible viscous flow equations to be used as an aid in the design of helicopter rotors. In order to further investigate the numerical procedure, the computer code developed to solve an approximate form of the three dimensional unsteady Navier-Stokes equations employing a linearized block implicit technique in conjunction with a QR operator scheme is tested. Results of calculations are presented for several two dimensional boundary layer flows including steady turbulent and unsteady laminar cases. A comparison of fourth order and second order solutions indicate that increased accuracy can be obtained without any significant increases in cost (run time). The results of the computations also indicate that the computer code can be applied to more complex flows such as those encountered on rotating airfoils. The geometry of a symmetric NACA four digit airfoil is considered and the appropriate geometrical properties are computed.
A procedure for solving three-dimensional, time-dependent turbulent flows is presented. The consistently split Linearized Block Implicit (LBI) scheme is used in conjunction with the QR Operator scheme to solve an approximate form of the Navier-Stokes equations in generalized nonorthogonal coordinates employing physical velocity components. Results of computations for both second order finite differences and the fourth order generalized Operator Compact Implicit (OCI) schemes are presented.
There is considerable interest in developing a numerical scheme for solving the time dependent viscous compressible three dimensional flow equations to aid in the design of helicopter rotors. The development of a computer code to solve a three dimensional unsteady approximate form of the Navier-Stokes equations employing a linearized block emplicit technique in conjunction with a QR operator scheme is described. Results of calculations of several Cartesian test cases are presented. The computer code can be applied to more complex flow fields such as these encountered on rotating airfoils.
There is considerable interest in developing a numerical scheme for solving the time dependent viscous compressible three dimensional flow to aid in the design of helicopter rotors. Numerical algorithms are examined to determine their overall suitability for the efficient and routine solution of an appropriate system of partial differential equations. It is concluded that a consistently split time linearized block implicit scheme using either quintic B spline collocation or the generalized operator compact implicit approach to generate a fourth order accurate algorithm is particularly well suited for use on the present problem. High cell Reynolds number behavior leads to favoring the generalized operator compact implicit approach over the quintic B spline collocation method.