Embedded function approach for turbulent flow prediction
Explore the source record for details and available documents.
Engineering topics
Publications and source records attributed to Werle, M. J..
Explore the source record for details and available documents.
In conventional prediction methods for turbulent flow influence on surface properties, very small mesh sizes and considerable computational effort is required to adequately resolve the intense velocity and temperature profile variations that occur in the wall-layer region. In this study, an approach is described wherein an outer region numerical solution is smoothly matched to a set of embedded analytic profile functions in the near-wall region; these wall-layer profile functions having been derived through consideration of the coherent structure of the time-dependent near-wall flow provide analytic expressions for the mean velocity and enthalpy profiles in the wall layer. The concept is demonstrated through example applications to turbulent boundary-layer flows in two dimensions. The technique is shown to be very efficient and it is demonstrated that a reduction of approximately half the mesh points across the 2-D layer may be realized (as compared to conventional methods) with no degradation in accuracy.
A general program was conducted to develop and assess a computational method for predicting the flow properties in a turbofan forced mixed duct. The detail assessment of the resulting computer code is presented. It was found that the code provided excellent predictions of the kinematics of the mixing process throughout the entire length of the mixer nozzle. The thermal mixing process between the hot core and cold fan flows was found to be well represented in the low speed portion of the flowfield.
The purpose of this paper is to review progress made in the solution of the interacting boundary-layer equations for subsonic flow. The interrelationship of triple deck theory and the interacting boundary-layer approach is discussed with emphasis placed on the development of efficient and reliable algorithms for the solution of the interacting boundary-layer equations. Example studies are presented for laminar and turbulent finite flat plate flow, laminar flow past a flat plate with a separation causing depression, and laminar and turbulent flow past a blunt based trailing edge.
The transformation permits a uniform mesh to be used in the computational coordinate which extends across the layer. This coordinate transformation uses the local value of the skin friction coefficient to scale the thickness of the wall layer region, and the local maximum value of turbulent viscosity to scale the boundary-layer thickness. Results are presented for two dimensional boundary layers in both positive and negative pressure gradients and comparisons are made with experimental data and conventional variable-grid results for low speed turbulent boundary-layers. The cases chosen illustrate the capability of this new transformation to capture the boundary layer growth over the full extent of laminar, transitional, and turbulent flow with no grid adjustment as well as its ability to consistently enlarge the wall layer region for accurate shear stress representation. Results of mesh refinement studies using the new coordinate transformation are presented.
The paper describes a new numerical scheme based on exponential difference operator concepts combined with Keller's (1968) box scheme approach to produce a stable second-order accurate finite-difference scheme for convection-diffusion problems arising in boundary layer flows in the presence of massive injection through a porous surface. The technique is demonstrated by application to the self-similar boundary layer equations with massive blowing at the surface.
A prediction method is developed for calculating distributions of surface heating rates, pressure and skin friction over a wavy wall in a two-dimensional supersonic flow. Of particular interest is the flow of thick turbulent boundary layers. The surface geometry and the flow conditions considered are such that there exists a strong interaction between the viscous and inviscid flow. First, using the interacting turbulent boundary layer equations, the problem is formulated in physical coordinates and then a reformulation of the governing equations in terms of Levy-Lees variables is given. Next, a numerical scheme for solving interacting boundary layer equations is adapted. A number of modifications which led to the improvement of the numerical algorithm are discussed. Finally, results are presented for flow over a train of up to six waves at various flow conditions.
This paper is concerned with the two-dimensional supersonic flow of a thick turbulent boundary layer over a train of relatively small wave-like protuberances. The flow conditions and the geometry are such that there exists a strong interaction between the viscous and inviscid flow. The problem cannot be solved without inclusion of interaction effects due to the occurrence of the separation singularity in classical boundary layer methods. Here the interacting boundary layer equations are solved numerically using a time-like relaxation method with turbulence effects represented by the inclusion of the eddy viscosity model of Cebeci and Smith. Results are presented for flow over a train of up to six waves for Mach numbers of 2.5 and 3.5, Reynolds numbers of 10,000,000/m and 32,000,000/m, and wall temperature ratios of 0.4 and 0.8. Limited comparisons with independent experimental and analytical results are also given.
A study of numerical schemes for solving viscous fluid flow problems with sizable regions of predominantly inviscid flow is presented. Difficulties associated with the familiar central difference approach for such problems were analyzed and alternative finite difference approaches employing windward concepts are presented. In addition, difference relations based on exponential operators were developed. All such schemes were demonstrated and evaluated through application to the case of Falkner Skan flow with blowing - a problem in which a sizable region of predominantly inviscid flow develops near the injection surface that traditionally causes numerical difficulty.
The two dimensional supersonic flow of a thick turbulent boundary layer over a train of relatively small wave-like protuberances is considered. The flow conditions and the geometry are such that there exists a strong interaction between the viscous and inviscid flow. The problem cannot be solved without inclusion of interaction effects due to the occurrence of the separation singularity in classical boundary layer methods. The interacting boundary layer equations are solved numerically using a time-like relaxation method with turbulence effects represented by the inclusion of the eddy viscosity model. Results are presented for flow over a train of up to six waves for Mach numbers of 10 and 32 million/meter, and wall temperature rations (T sub w/T sub 0) of 0.4 and 0.8. Limited comparisons with independent experimental and analytical results are also given. Detailed results on the influence of small protuberances on surface heating by boundary layers are presented.
A study and application of the fourth order spline collocation procedure, numerical solution of boundary layer like differential equations, is presented. A simple inversion algorithm for the simultaneous solution of the resulting difference equations is given. Particular attention is focused on the boundary condition representation for the spline second derivative approximations. Solutions using the spline procedure, as well as the three point finite difference method, are presented for several model problems in order to assess and improve the spline numerical scheme. Application of the resulting algorithm to the incompressible laminar self similar boundary layer equations is presented.
A numerical algorithm is presented for solving laminar, steady, supersonic interacting boundary-layer flows for quasi-three-dimensional configurations. The interaction problem is treated as a boundary-value problem and a salient feature of the scheme is the direct implementation of the downstream boundary condition. Solutions are presented for axisymmetric and swept (yawed) compression ramps for both adiabatic and heat transfer conditions over a Mach number range of 2-6. The results are in good agreement with experimental data and existing theories for axisymmetric cases. For the swept (yawed) configurations, lack of experimental data makes a direct comparison impossible, but the present solutions are found to be in qualitative agreement with earlier studies. In addition it is shown that the trends obtained here for the sweep effects are well predicted by a simple extension of the two-dimensional asymptotic theory.
Supersonic turbulent boundary layers over two-dimensional protuberances are investigated, using the numerical finite difference alternating direction implicit (ADI) method. The turbulence is modeled mathematically. The turbulence is represented here by the eddy viscosity approach. The turbulent boundary layer structure as well as an interest in thick boundary layers and much larger protuberance heights than in the laminar case lead to new difficulties. The problems encountered and the means to remove them are discussed.
The title problem is considered for the case of flow past a circular cylinder placed normal to a uniform mainstream with Reynolds numbers from 40 to 200. Implicit finite difference numerical solutions are obtained for a set of boundary-layer equations that account for the second order effects associated with surface curvature and displacement speed. It was found that both of these contributors have a significant influence on the internal structure of the viscous region and that an accurate estimate of the surface pressure distribution is essential for estimating the surface shear stress.
The title problem was studied employing implicit finite-difference methods to obtain numerical solutions to a composite set of boundary-layer equations valid to second order. Results are given for flow up a two-dimensional cubic compression ramp for free-stream Mach numbers of 6, 8, and 12.25 and for free-stream Re/inch equal to 85,800 and 25,800 at a wall-to-stagnation temperature ratio of 0.223. Comparisons with independent theories and experimental results are given. Nonsingular separation was produced at a free-stream Mach number of 12.25. For all cases considered, displacement and curvature effects canceled one another when a consistent treatment of inviscid and viscous curvature corrections was employed - the second-order theory virtually reproducing the first-order results.