scholarly journals Natural convection in a shallow cavity with differentially heated end walls. Part 2. Numerical solutions

1974 ◽  
Vol 65 (2) ◽  
pp. 231-246 ◽  
Author(s):  
D. E. Cormack ◽  
L. G. Leal ◽  
J. H. Seinfeld

Numerical solutions of the full Navier-Stokes equations are obtained for the problem of natural convection in closed cavities of small aspect ratio with differentially heated end walls. These solutions cover the parameter range Pr = 6·983, 10 ≤ Gr 2 × 104 and 0·05 [les ] A [les ] 1. A comparison with the asymptotic theory of part 1 shows excellent agreement between the analytical and numerical solutions provided that A [lsim ] 0·1 and Gr2A3Pr2 [lsim ] 105. In addition, the numerical solutions demonstrate the transition between the shallow-cavity limit of part 1 and the boundary-layer limit; A fixed, Gr → ∞.

1990 ◽  
Vol 220 ◽  
pp. 397-411 ◽  
Author(s):  
Laura L. Pauley ◽  
Parviz Moin ◽  
William C. Reynolds

The separation of a two-dimensional laminar boundary layer under the influence of a suddenly imposed external adverse pressure gradient was studied by time-accurate numerical solutions of the Navier–Stokes equations. It was found that a strong adverse pressure gradient created periodic vortex shedding from the separation. The general features of the time-averaged results were similar to experimental results for laminar separation bubbles. Comparisons were made with the ‘steady’ separation experiments of Gaster (1966). It was found that his ‘bursting’ occurs under the same conditions as our periodic shedding, suggesting that bursting is actually periodic shedding which has been time-averaged. The Strouhal number based on the shedding frequency, local free-stream velocity, and boundary-layer momentum thickness at separation was independent of the Reynolds number and the pressure gradient. A criterion for onset of shedding was established. The shedding frequency was the same as that predicted for the most amplified linear inviscid instability of the separated shear layer.


According to Stewartson (1969, 1974) and to Messiter (1970), the flow near the trailing edge of a flat plate has a limit structure for Reynolds number Re →∞ consisting of three layers over a distance O (Re -3/8 ) from the trailing edge: the inner layer of thickness O ( Re -5/8 ) in which the usual boundary layer equations apply; an intermediate layer of thickness O ( Re -1/2 ) in which simplified inviscid equations hold, and the outer layer of thickness O ( Re -3/8 ) in which the full inviscid equations hold. These asymptotic equations have been solved numerically by means of a Cauchy-integral algorithm for the outer layer and a modified Crank-Nicholson boundary layer program for the displacement-thickness interaction between the layers. Results of the computation compare well with experimental data of Janour and with numerical solutions of the Navier-Stokes equations by Dennis & Chang (1969) and Dennis & Dunwoody (1966).


2002 ◽  
Vol 465 ◽  
pp. 99-130 ◽  
Author(s):  
A. V. OBABKO ◽  
K. W. CASSEL

Numerical solutions of the unsteady Navier–Stokes equations are considered for the flow induced by a thick-core vortex convecting along a surface in a two-dimensional incompressible flow. The presence of the vortex induces an adverse streamwise pressure gradient along the surface that leads to the formation of a secondary recirculation region followed by a narrow eruption of near-wall fluid in solutions of the unsteady boundary-layer equations. The locally thickening boundary layer in the vicinity of the eruption provokes an interaction between the viscous boundary layer and the outer inviscid flow. Numerical solutions of the Navier–Stokes equations show that the interaction occurs on two distinct streamwise length scales depending upon which of three Reynolds-number regimes is being considered. At high Reynolds numbers, the spike leads to a small-scale interaction; at moderate Reynolds numbers, the flow experiences a large-scale interaction followed by the small-scale interaction due to the spike; at low Reynolds numbers, large-scale interaction occurs, but there is no spike or subsequent small-scale interaction. The large-scale interaction is found to play an essential role in determining the overall evolution of unsteady separation in the moderate-Reynolds-number regime; it accelerates the spike formation process and leads to formation of secondary recirculation regions, splitting of the primary recirculation region into multiple corotating eddies and ejections of near-wall vorticity. These eddies later merge prior to being lifted away from the surface and causing detachment of the thick-core vortex.


1998 ◽  
Vol 120 (1) ◽  
pp. 70-75 ◽  
Author(s):  
K. O. Homan ◽  
S. L. Soo

This paper treats the steady flow of a wall jet into a large-width cavity for which the primary axis is normal to the direction of the jet inflow. Numerical solutions of the two-dimensional Navier-Stokes equations are computed for inlet Reynolds numbers of 10 to 50 and tank width to inlet height ratios of 16 to 128. The length and velocity scales of the wall jet boundary layer exhibit close agreement with the classic wall jet similarity solution and published experimental data but the width of the region for which the comparison proves to be favorable has a limited extent. This departure from a self-similar evolution of the wall jet is shown to result from the finite domain width and its influence on the large recirculation cell located immediately above the wall jet boundary layer.


1991 ◽  
Vol 113 (4) ◽  
pp. 574-578 ◽  
Author(s):  
V. C. Patel ◽  
J. Tyndall Chon ◽  
J. Y. Yoon

The boundary layer over a wavy wall and fully-developed flow in a duct with a wavy wall are considered. Numerical solutions of the Navier-Stokes equations have been obtained to provide insights into the various steady flow regimes that are possible, and to illustrate the nuances of predicting flows containing multiple separation and reattachment points.


1975 ◽  
Vol 69 (4) ◽  
pp. 631-656 ◽  
Author(s):  
W. Roger Briley ◽  
Henry Mcdonald

A method is presented for performing detailed computations of thin incompressible separation bubbles on smooth surfaces. The analysis consists of finite-difference solutions to the time-dependent boundary-layer or Navier-Stokes equations for the flow in the immediate vicinity of the bubble. The method employs the McDonald-Fish turbulence model, to predict the development of the time-mean flow field, as influenced by the free-stream turbulence level. It also employs a viscous-inviscid interaction model, which accounts for the elliptic interaction between the shear layer and inviscid free stream. The numerical method is based on an alternating-direction implicit scheme for the vorticity equation. It employs transformations, to allow the free-stream boundary to change in time with the shape of the computed shear layer, and to ensure an adequate resolution of the sublayer region. Numerical solutions are presented for transitional bubbles on an NACA 663-018 airfoil at zero angle of incidence with chordal Reynolds numbers of 2·0 × 106 and 1·7 × 106. These have a qualitative behaviour similar to that observed in numerous experiments; they are also in reasonable quantitative agreement with available experimental data. Little difference is found between steady solutions of the boundary-layer and Navier-Stokes equations for these flow conditions. Numerical studies based on mesh refinement suggest that the well-known singularity at separation, which is present in conventional solutions of the steady boundary-layer equations when the free-stream velocity is specified, is effectively removed when viscous-inviscid interaction is allowed to influence the imposed velocity distribution.


Numerical solutions of the Navier-Stokes equations are given for the steady, two-dimensional, laminar flow of an incompressible fluid through a channel with a symmetric constriction in the form of a semi-infinite step change in width. The flow proceeds from a steady Poiseuille velocity distribution far enough upstream of the step in the wider part of the channel to a corresponding distribution downstream in the narrower part and is assumed to remain symmetrical about the centre line of the channel. The numerical scheme involves an accurate and efficient centred difference treatment developed by Dennis & Hudson (1978) and results are obtained for Reynolds numbers, based on half the volumetric flow rate, up to 2000. For a step that halves the width of the channel it is found that very fine uniform grids, with 80 intervals spaced across half of the wider channel upstream, are necessary for resolution of the solution for the flow downstream of the onset of the step. Slightly less refined grids are adequate to resolve the flow upstream. The calculated flow ahead of the step exhibits very good agreement with the asymptotic theory of Smith (1979 b)for Reynolds numbers greater than about 100; indeed, comparisons of the upstream separation position and of the wall vorticity nearby are believed to yield the best agreement between numerical and asymptotic solutions yet found. Downstream there is also qualitative agreement regarding separation and reattachment as the grid size is refined in the numerical results.


Author(s):  
Marcel Escudier

This chapter starts by introducing the concept of a boundary layer and the associated boundary-layer approximations. The laminar boundary-layer equations are then derived from the Navier-Stokes equations. The assumption of velocity-profile similarity is shown to reduce the partial differential boundary-layer equations to ordinary differential equations. The results of numerical solutions to these equations are discussed: Blasius’ equation, for zero-pressure gradient, and the Falkner-Skan equation for wedge flows. Von Kármán’s momentum-integral equation is derived and used to obtain useful results for the zero-pressure-gradient boundary layer. Pohlhausen’s quartic-profile method is then discussed, followed by the approximate method of Thwaites. The chapter concludes with a qualitative account of the way in which aerodynamic lift is generated.


1971 ◽  
Vol 47 (4) ◽  
pp. 713-736 ◽  
Author(s):  
W. Roger Briley

The flow in a two-dimensional laminar separation bubble is analyzed by means of finite-difference solutions to the Navier-Stokes equations for incompressible flow. The study was motivated by the need to analyze high-Reynolds-number flow fields having viscous regions in which the boundary-layer assumptions are questionable. The approach adopted in the present study is to analyze the flow in the immediate vicinity of the separation bubble using the Navier-Stokes equations. It is assumed that the resulting solutions can then be patched to the remainder of the flow field, which is analyzed using boundary-layer theory and inviscid-flow analysis. Some of the difficulties associated with patching the numerical solutions to the remainder of the flow field are discussed, and a suggestion for treating boundary conditions is made which would permit a separation bubble to be computed from the Navier-Stokes equations using boundary conditions from inviscid and boundary-layer solutions without accounting for interaction between individual flow regions. Numerical solutions are presented for separation bubbles having Reynolds numbers (based on momentum thickness) of the order of 50. In these numerical solutions, separation was found to occur without any evidence of the singular behaviour at separation found in solutions to the boundary-layer equations. The numerical solutions indicate that predictions of separation by boundary-layer theory are not reliable for this range of Reynolds number. The accuracy and validity of the numerical solutions are briefly examined. Included in this examination are comparisons between the Howarth solution of the boundary-layer equations for a linearly retarded freestream velocity and the corresponding numerical solutions of the Navier-Stokes equations for various Reynolds numbers.


Sign in / Sign up

Export Citation Format

Share Document