Engineering PapersSearch

Engineering topics

Spalart, P. R.

Publications and source records attributed to Spalart, P. R..

At least 19 records

Numerical Study of the Effect of Mean Three-Dimensionality on Turbulence in Adverse-Pressure-Gradient Boundary Layers

Direct numerical simulation (DNS) is used to isolate the influence of sweep on a separating turbulent boundary layer. Attention here is limited to the behavior of the turbulence within the adverse-pressure-gradient (APG) region upstream of separation. Other regions and quantities are considered in Coleman, Rumsey & Spalart. The mean three-dimensionality and outer-layer inviscid skewing have only a slight effect upon the structure of the turbulence (measured by the relationship of the components of the Reynolds-stress tensor and the efficiency of the turbulence energy transfer) compared with that of the adverse pressure gradient, which dominates both the skewed and unskewed layers.

Coleman, G. N.

Numerical Study of Turbulent Separation Bubbles with Varying Pressure Gradient and Reynolds Number

A family of cases each containing a small separation bubble is treated by direct numerical simulation (DNS), varying two parameters: the severity of the pressure gradients, generated by suction and blowing across the opposite boundary, and the Reynolds number. Each flow contains a well-developed entry region with essentially zero pressure gradient, and all are adjusted to have the same value for the momentum thickness, extrapolated from the entry region to the centre of the separation bubble. Combined with fully defined boundary conditions this will make comparisons with other simulations and turbulence models rigorous; we present results for a set of eight Reynolds-averaged Navier–Stokes turbulence models. Even though the largest Reynolds number is approximately 5.5 times higher than in a similar DNS study we presented in 1997, the models have difficulties matching the DNS skin friction very closely even in the zero pressure gradient, which complicates their assessment. In the rest of the domain, the separation location per se is not particularly difficult to predict, and the most definite disagreement between DNS and models is near reattachment. Curiously, the better models tend to cluster together in their predictions of pressure and skin friction even when they deviate from the DNS, although their eddy-viscosity levels are widely different in the outer region near the bubble (or they do not rely on an eddy viscosity). Stratford’s square-root law is satisfied by the velocity profiles, both at separation and reattachment. The Reynolds-number range covers a factor of two, with the Reynolds number based on the extrapolated momentum thickness equal to approximately 1500 and 3000. This allows tentative estimates of the improvements that even higher values will bring to the model comparisons. The solutions are used to assess models through pressure, skin friction and other measures; the flow fields are also used to produce effective eddy-viscosity targets for the models, thus guiding turbulence-modelling work in each region of the flow.

Coleman, G. N.

Direct Numerical Simulation and Theories of Wall Turbulence with a Range of Pressure Gradients

A new Direct Numerical Simulation (DNS) of Couette-Poiseuille flow at a higher Reynolds number is presented and compared with DNS of other wall-bounded flows. It is analyzed in terms of testing semi-theoretical proposals for universal behavior of the velocity, mixing length, or eddy viscosity in pressure gradients, and in terms of assessing the accuracy of two turbulence models. These models are used in two modes, the traditional one with only a dependence on the wall-normal coordinate y, and a newer one in which a lateral dependence on z is added. For pure Couette flow and the Couette-Poiseuille case considered here, this z-dependence allows some models to generate steady streamwise vortices, which generally improves the agreement with DNS and experiment. On the other hand, it complicates the comparison between DNS and models.

Coleman, G. N.

Direct simulation of the stably stratified turbulent Ekman layer

The Navier-Stokes equations and the Boussinesq approximation were used to compute a 3D time-dependent turbulent flow in the stably stratified Ekman layer over a smooth surface. The simulation data are found to be in very good agreement with atmospheric measurements when nondimensionalized according to Nieuwstadt's local scaling scheme. Results suggest that, when Reynolds number effects are taken into account, the 'constant Froud number' stable layer model (Brost and Wyngaard, 1978) and the 'shearing length' stable layer model (Hunt, 1985) for the dissipitation rate of turbulent kinetic energy are both valid. It is concluded that there is good agreement between the direct numerical simulation results and large-eddy simulation results obtained by Mason and Derbyshire (1990).

Coleman, G. N.

Linear and nonlinear stability of the Blasius boundary layer

Two new techniques for the study of the linear and nonlinear instability in growing boundary layers are presented. The first technique employs partial differential equations of parabolic type exploiting the slow change of the mean flow, disturbance velocity profiles, wavelengths, and growth rates in the streamwise direction. The second technique solves the Navier-Stokes equation for spatially evolving disturbances using buffer zones adjacent to the inflow and outflow boundaries. Results of both techniques are in excellent agreement. The linear and nonlinear development of Tollmien-Schlichting (TS) waves in the Blasius boundary layer is investigated with both techniques and with a local procedure based on a system of ordinary differential equations. The results are compared with previous work and the effects of non-parallelism and nonlinearity are clarified. The effect of nonparallelism is confirmed to be weak and, consequently, not responsible for the discrepancies between measurements and theoretical results for parallel flow.

Bertolotti, F. P.

Quasi-Coherent Structures In Turbulent Boundary Layers

Two-part report reviews knowledge of coherent structures in turbulent boundary layers. Part I describes processes and status of cooperative project to summarize data from research on boundary-layer turbulence. Part II presents results of study of numerically simulated flat-plate canonical turbulent boundary layer.

Robinson, S. K.

Comment paper: Workshop on Engineering Turbulence Modeling

The speaker for this paper describes and evaluates a k-epsilon model for calculating Samuel-Joubert flow. He proceeds to present both Boeing's and his positions on the state-of-the-art in this area and future goals. Finally, presented is a one equation mathematical model for calculating Samuel-Joubert flow. All results are presented in viewgraph format.

Spalart, P. R.

Spectral Method For Simulation Of Vortex Rings

Method of computation relying on spectral basis functions developed especially for simulation of axisymmetric vortex rings in incompressible, viscous fluid with quiescent far field. Contributes to understanding of flows in and around vortex rings during long propagation times, including such theoretically and practically important phenomena as drift and expansion of ring, "leapfrogging" and coalescence of two rings, and shedding of vorticity into wake of propagating ring.

Stanaway, S. K.

A numerical study of the turbulent Ekman layer

The three-dimensional time-dependent turbulent flow in a neutrally stratified Ekman layer over a smooth flat surface was numerically simulated by directly solving the Navier-Stokes equations. Issues addressed using the direct numerical simulation (DNS) fields include the presence or absence of large-scale coherent structures ('longitudinal' or 'roll' vortices) in neutrally stratified Ekman-layer turbulence, the effects of the horizontal component of the angular velocity vector (i.e., latitude), and implications for models of the PBL. Experimental and DNS profiles are compared.

Coleman, G. N.

Computation of instability and transition

The state of transition calculations in 1989 and desirable future trends are analyzed within the framework of aeronautical applications with emphasis on wall-bounded flows. A shift to studies of the early stages of transition, particularly to receptivity, from studies of the late stages in excessively confined periodic domains is proposed. It is pointed out that early-stage studies require a larger computational domain and more difficult boundary conditions but not a very fine grid. Intermediate strategies, more involved than the classical Orr-Sommerfeld equation but much more economical than full Navier-Stokes equations, are emphasized, along with the need for compressible calculations.

Spalart, P. R.

Direct numerical study of crossflow instability

Disturbances in the swept Hiemenz flow are calculated by solving the Navier-Stokes equations. The spatially-evolving base flow is treated exactly, allowing a check of the 'local' stability theories. Different types of disturbances such as random noise, waves, and wave packets, are input near the attachment line, develop in space, and exit through an outflow boundary. They all generate streamwise vortices. The effect of the Reynolds number, of the time-dependence of the noise, and of nonlinearity, are investigated.

Spalart, P. R.

Understanding transition and turbulence through direct simulations

Direct simulations consist in solving the full Navier-Stokes equations, without any turbulence model, and describing all the detailed features of the flow. Usually the flows are three-dimensional and time-dependent and contain both coarse and fine structures, which makes the numerical task very challenging in terms of both the algorithm and the computational effort. Most of the work until now has involved spectral methods, which are highly accurate but not very flexible in terms of geometry or complex equations. For that reason, future work will also rely on high-order finite-difference or other methods. Direct simulations complement experimental work, and both contribute to the theory and the empirical knowledge of turbulence. Once such a simulation has been shown to be accurate, the flow field is completely known in three dimensions and time, including the pressure, the vorticity and any other quantity. On the other hand, most simulations to date solved the incompressible equations in rather simple geometries, and direct simulations will always be limited to moderate Reynolds numbers. Extensive simulations have been conducted in homogeneous turbulence, channel flows, boundary layers, and mixing layers. Much effort is devoted to addressing flows with compressibility and chemical reactions, and to new geometries such as a backward-facing step.

Spalart, P. R.

Numerical and experimental evaluations of the flow past nested chevrons

An effort is made to contribute to the development of CFD by relating the successful use of vortex dynamics in the computation of the pressure drop past a planar array of chevron-shaped obstructions. An ensemble of results was used to compute the loss coefficient k, stimulating an experimental program for the assessment of the measured loss coefficient for the same geometry. The most provocative result of this study has been the representation of kinetic energy production in terms of vorticity source terms.

Foss, J. F.

A review of quasi-coherent structures in a numerically simulated turbulent boundary layer

Preliminary results of a comprehensive study of the structural aspects of a numerically simulated number turbulent boundary layer are presented. A direct Navier-Stokes simulation of a flat-plate, zero pressure gradient boundary layer at Re0 = 670 was used. Most of the known nonrandom, coherent features of turbulent boundary layers are confirmed in the simulation, and several new aspects of their spatial character are reported. The spatial relationships between many of the various structures are described, forming the basis for a more complete kinematical picture of boundary layer physics than has been previously known. In particular, the importance of vortex structures of various forms to the generation of Reynolds shear stress is investigated.

Robinson, S. K.

Direct numerical study of leading-edge contamination

Instability, turbulence, and relaminarization in the attachment-line region of swept and unswept cylindrical bodies are studied by numerical solution of the full Navier-Stokes equations. The flow is simulated over a strip containing the attachment-line and treated as homogeneous in the spanwise direction; the disturbances decay exponentially upstream. Transpiration through the wall may be prescribed. The new method, which admits completely general disturbance, agrees with published linear-stability results, which were limited to an apparently restrictive form of disturbance. Fully developed turbulent solutions with sweep are generated and compare well with the experiment. The turbulence is subcritical (except for blowing), resulting in large hysteresis loops. By lowering the sweep Reynolds number, or increasing the suction, the turbulent flow is made to relaminarize. The relaminarization Reynolds number is much less sensitive to suction than the linear-stability Reynolds number. Extensive attempts to detect the postulated nonlinear instability of the unswept flow failed, suggesting that this flow is linearly and nonlinearly stable.

Spalart, P. R.

A numerical study of bifurcations in a barotropic shear flow

In the last few years, more and more evidence has emerged suggesting that transition to turbulence may be viewed as a succession of bifurcations to deterministic chaos. Most experimental and numerical observations have been restricted to Rayleigh-Benard convection and Taylor-Couette flow between concentric cylinders. An attempt is made to accurately describe the bifurcation sequence leading to chaos in a 2-D temporal free shear layer on the beta-plane. The beta-plane is a locally Cartesian reduction of the equations describing the dynamicss of a shallow layer of fluid on a rotating spherical planet. It is a valid model for large scale flows of interest in meteorology and oceanography.

Huerre, P.

Navier-Stokes simulations of axisymmetric vortex rings

The solution of the incompressible Navier-Stokes equations by means of a spectral method is used in the present numerical investigation of viscous axisymmetric vortex rings, in the cases of a single ring over the 0.001-1000 Re number range, and a pair of interacting rings. The propagation speeds of vortex rings of varying Re are computed and compared with the results of Saffman's (1970) theory; the present error estimate is found to be smaller than Saffman's, and to decrease with increasing Re at fixed core/ring radius ratio. The 'leapfrogging' of two vortex rings of equal sign and Re=1000 each is observed, indicating severe strain of the vortex core through the first pass and merging during the second pass.

Stanaway, S. K.