scholarly journals A class of exact Navier–Stokes solutions for homogeneous flat-plate boundary layers and their linear stability

2014 ◽  
Vol 752 ◽  
pp. 462-484 ◽  
Author(s):  
Michael O. John ◽  
Dominik Obrist ◽  
Leonhard Kleiser

AbstractWe introduce a new boundary layer formalism on the basis of which a class of exact solutions to the Navier–Stokes equations is derived. These solutions describe laminar boundary layer flows past a flat plate under the assumption of one homogeneous direction, such as the classical swept Hiemenz boundary layer (SHBL), the asymptotic suction boundary layer (ASBL) and the oblique impingement boundary layer. The linear stability of these new solutions is investigated, uncovering new results for the SHBL and the ASBL. Previously, each of these flows had been described with its own formalism and coordinate system, such that the solutions could not be transformed into each other. Using a new compound formalism, we are able to show that the ASBL is the physical limit of the SHBL with wall suction when the chordwise velocity component vanishes while the homogeneous sweep velocity is maintained. A corresponding non-dimensionalization is proposed, which allows conversion of the new Reynolds number definition to the classical ones. Linear stability analysis for the new class of solutions reveals a compound neutral surface which contains the classical neutral curves of the SHBL and the ASBL. It is shown that the linearly most unstable Görtler–Hämmerlin modes of the SHBL smoothly transform into Tollmien–Schlichting modes as the chordwise velocity vanishes. These results are useful for transition prediction of the attachment-line instability, especially concerning the use of suction to stabilize boundary layers of swept-wing aircraft.

1995 ◽  
Vol 298 ◽  
pp. 211-248 ◽  
Author(s):  
U. Rist ◽  
H. Fasel

The three-dimensional development of controlled transition in a flat-plate boundary layer is investigated by direct numerical simulation (DNS) using the complete Navier-Stokes equations. The numerical investigations are based on the so-called spatial model, thus allowing realistic simulations of spatially developing transition phenomena as observed in laboratory experiments. For solving the Navier-Stokes equations, an efficient and accurate numerical method was developed employing fourth-order finite differences in the downstream and wall-normal directions and treating the spanwise direction pseudo-spectrally. The present paper focuses on direct simulations of the wind-tunnel experiments by Kachanov et al. (1984, 1985) of fundamental breakdown in controlled transition. The numerical results agreed very well with the experimental measurements up to the second spike stage, in spite of relatively coarse spanwise resolution. Detailed analysis of the numerical data allowed identification of the essential breakdown mechanisms. In particular, from our numerical data, we could identify the dominant shear layers and vortical structures that are associated with this breakdown process.


2016 ◽  
Vol 792 ◽  
pp. 499-525 ◽  
Author(s):  
Hui Xu ◽  
Spencer J. Sherwin ◽  
Philip Hall ◽  
Xuesong Wu

This paper is concerned with the behaviour of Tollmien–Schlichting (TS) waves experiencing small localised distortions within an incompressible boundary layer developing over a flat plate. In particular, the distortion is produced by an isolated roughness element located at $\mathit{Re}_{x_{c}}=440\,000$. We considered the amplification of an incoming TS wave governed by the two-dimensional linearised Navier–Stokes equations, where the base flow is obtained from the two-dimensional nonlinear Navier–Stokes equations. We compare these solutions with asymptotic analyses which assume a linearised triple-deck theory for the base flow and determine the validity of this theory in terms of the height of the small-scale humps/indentations taken into account. The height of the humps/indentations is denoted by $h$, which is considered to be less than or equal to $x_{c}\mathit{Re}_{x_{c}}^{-5/8}$ (corresponding to $h/{\it\delta}_{99}<6\,\%$ for our choice of $\mathit{Re}_{x_{c}}$). The rescaled width $\hat{d}~(\equiv d/(x_{c}\mathit{Re}_{x_{c}}^{-3/8}))$ of the distortion is of order $\mathit{O}(1)$ and the width $d$ is shorter than the TS wavelength (${\it\lambda}_{TS}=11.3{\it\delta}_{99}$). We observe that, for distortions which are smaller than 0.1 of the inner deck height ($h/{\it\delta}_{99}<0.4\,\%$), the numerical simulations confirm the asymptotic theory in the vicinity of the distortion. For larger distortions which are still within the inner deck ($0.4\,\%<h/{\it\delta}_{99}<5.5\,\%$) and where the flow is still attached, the numerical solutions show that both humps and indentations are destabilising and deviate from the linear theory even in the vicinity of the distortion. We numerically determine the transmission coefficient which provides the relative amplification of the TS wave over the distortion as compared to the flat plate. We observe that for small distortions, $h/{\it\delta}_{99}<5.5\,\%$, where the width of the distortion is of the order of the boundary layer, a maximum amplification of only 2 % is achieved. This amplification can however be increased as the width of the distortion is increased or if multiple distortions are present. Increasing the height of the distortion so that the flow separates ($7.2\,\%<h/{\it\delta}_{99}<12.8\,\%$) leads to a substantial increase in the transmission coefficient of the hump up to 350 %.


2002 ◽  
Vol 124 (4) ◽  
pp. 868-875 ◽  
Author(s):  
H. Yan ◽  
D. Knight ◽  
A. A. Zheltovodov

A supersonic flat-plate boundary layer at a Reynolds number of 2×104 based on the inflow boundary layer thickness is investigated at different Mach numbers (M=2.88 and 4) using the monotonically integrated large-eddy simulation (MILES) technique. The inherent numerical dissipation is taken as an implicit subgrid scales (SGS) model to close the Favre-filtered compressible Navier-Stokes (NS) equations. A finite volume method with second-order accuracy in time and space is implemented for the solution of the Navier-Stokes equations on an unstructured grid of tetrahedra. The heat transfer coefficient is predicted by simulating both adiabatic and isothermal cases. The mean flowfield and turbulent stresses are in good agreement with experiment. The relationship between the predicted skin friction coefficient and heat transfer coefficient is in close agreement with the Reynolds analogy factor. The variation of turbulent Prandtl number cross the boundary layer falls within the experimental envelope. These are the first LES predictions of adiabatic and isothermal supersonic flat plate boundary layers using the MILES technique.


1992 ◽  
Vol 114 (3) ◽  
pp. 306-312 ◽  
Author(s):  
R. E. Spall ◽  
M. R. Malik

The linear stability of the low-speed three-dimensional flow over a flat plate with an attached cylinder is studied. The region of interest is upstream of the initial separation point and includes the effects of both adverse and favorable pressure gradients, as well as crossflow. The resulting boundary-layer is subject to both the Tollmien-Schlichting (TS) and crossflow instabilities. Linear stability calculations, using N-factor correlations, indicate that the transition process would be dominated by TS instabilities, although for low frequencies crossflow-type disturbances are important.


1991 ◽  
Vol 113 (4) ◽  
pp. 608-616 ◽  
Author(s):  
H. M. Jang ◽  
J. A. Ekaterinaris ◽  
M. F. Platzer ◽  
T. Cebeci

Two methods are described for calculating pressure distributions and boundary layers on blades subjected to low Reynolds numbers and ramp-type motion. The first is based on an interactive scheme in which the inviscid flow is computed by a panel method and the boundary layer flow by an inverse method that makes use of the Hilbert integral to couple the solutions of the inviscid and viscous flow equations. The second method is based on the solution of the compressible Navier–Stokes equations with an embedded grid technique that permits accurate calculation of boundary layer flows. Studies for the Eppler-387 and NACA-0012 airfoils indicate that both methods can be used to calculate the behavior of unsteady blade boundary layers at low Reynolds numbers provided that the location of transition is computed with the en method and the transitional region is modeled properly.


1995 ◽  
Vol 291 ◽  
pp. 369-392 ◽  
Author(s):  
Ronald D. Joslin

The spatial evolution of three-dimensional disturbances in an attachment-line boundary layer is computed by direct numerical simulation of the unsteady, incompressible Navier–Stokes equations. Disturbances are introduced into the boundary layer by harmonic sources that involve unsteady suction and blowing through the wall. Various harmonic-source generators are implemented on or near the attachment line, and the disturbance evolutions are compared. Previous two-dimensional simulation results and nonparallel theory are compared with the present results. The three-dimensional simulation results for disturbances with quasi-two-dimensional features indicate growth rates of only a few percent larger than pure two-dimensional results; however, the results are close enough to enable the use of the more computationally efficient, two-dimensional approach. However, true three-dimensional disturbances are more likely in practice and are more stable than two-dimensional disturbances. Disturbances generated off (but near) the attachment line spread both away from and toward the attachment line as they evolve. The evolution pattern is comparable to wave packets in flat-plate boundary-layer flows. Suction stabilizes the quasi-two-dimensional attachment-line instabilities, and blowing destabilizes these instabilities; these results qualitatively agree with the theory. Furthermore, suction stabilizes the disturbances that develop off the attachment line. Clearly, disturbances that are generated near the attachment line can supply energy to attachment-line instabilities, but suction can be used to stabilize these instabilities.


It is shown that the boundary layer approximation to the flow of a viscous fluid past a flat plate of length l , generally valid near the plate when the Reynolds number Re is large, fails within a distance O( lRe -3/4 ) of the trailing edge. The appropriate governing equations in this neighbourhood are the full Navier- Stokes equations. On the basis of Imai (1966) these equations are linearized with respect to a uniform shear and are then completely solved by means of a Wiener-Hopf integral equation. The solution so obtained joins smoothly on to that of the boundary layer for a flat plate upstream of the trailing edge and for a wake downstream of the trailing edge. The contribution to the drag coefficient is found to be O ( Re -3/4 ) and the multiplicative constant is explicitly worked out for the linearized equations.


2015 ◽  
Vol 771 ◽  
pp. 520-546 ◽  
Author(s):  
Nicola De Tullio ◽  
Anatoly I. Ruban

The capabilities of the triple-deck theory of receptivity for subsonic compressible boundary layers have been thoroughly investigated through comparisons with numerical simulations of the compressible Navier–Stokes equations. The analysis focused on the two Tollmien–Schlichting wave linear receptivity problems arising due to the interaction between a low-amplitude acoustic wave and a small isolated roughness element, and the low-amplitude time-periodic vibrations of a ribbon placed on the wall of a flat plate. A parametric study was carried out to look at the effects of roughness element and vibrating ribbon longitudinal dimensions, Reynolds number, Mach number and Tollmien–Schlichting wave frequency. The flat plate is considered isothermal, with a temperature equal to the laminar adiabatic-wall temperature. Numerical simulations of the full and the linearised compressible Navier–Stokes equations have been carried out using high-order finite differences to obtain, respectively, the steady basic flows and the unsteady disturbance fields for the different flow configurations analysed. The results show that the asymptotic theory and the Navier–Stokes simulations are in good agreement. The initial Tollmien–Schlichting wave amplitudes and, in particular, the trends indicated by the theory across the whole parameter space are in excellent agreement with the numerical results. An important finding of the present study is that the behaviour of the theoretical solutions obtained for $\mathit{Re}\rightarrow \infty$ holds at finite Reynolds numbers and the only conditions needed for the theoretical predictions to be accurate are that the receptivity process be linear and the free-stream Mach number be subsonic.


1990 ◽  
Vol 221 ◽  
pp. 311-347 ◽  
Author(s):  
H. Fasel ◽  
U. Konzelmann

Non-parallel effects which are due to the growing boundary layer are investigated by direct numerical integration of the complete Navier-Stokes equations for incompressible flows. The problem formulation is spatial, i.e. disturbances may grow or decay in the downstream direction as in the physical experiments. In the past various non-parallel theories were published that differ considerably from each other in both approach and interpretation of the results. In this paper a detailed comparison of the Navier-Stokes calculation with the various non-parallel theories is provided. It is shown, that the good agreement of some of the theories with experiments is fortuitous and that the difference between experiments and theories concerning the branch I neutral location cannot be explained by non-parallel effects.


Author(s):  
H. M. Jang ◽  
M. F. Platzer ◽  
J. A. Ekaterinaris ◽  
T. Cebeci

Two methods are described for calculating pressure distributions and boundary layers on blades subjected to low Reynolds numbers and ramp–type motion. The first is based on an interactive scheme in which the inviscid flow is computed by a panel method and the boundary layer flow by an inverse method that makes use of the Hilbert integral to couple the solutions of the inviscid and viscous flow equations. The second method is based on the solution of the compressible Navier–Stokes equations with an embedded grid technique that permits accurate calculation of boundary layer flows. Studies for the Eppler and NACA–0012 airfoils indicate that both methods can be used to calculate the behavior of unsteady blade boundary layers at low Reynolds numbers provided that the location of transition is computed with the en–method and the transitional region is modelled properly.


Sign in / Sign up

Export Citation Format

Share Document