Numerical Integration Procedure for the Steady State Navier-Stokes Equations

1969 ◽  
Vol 11 (5) ◽  
pp. 445-453 ◽  
Author(s):  
A. K. Runchal ◽  
M. Wolfshtein

A procedure is proposed for the integration of the full Navier-Stokes equations for constant-property two-dimensional flows. In contrast with earlier procedures, the present one is capable of dealing with cases of very high Reynolds number. The power of the new procedure is demonstrated in two cases: (1) the square recirculating eddy, for which no solution was previously available for Reynolds numbers larger than about 400, and (2) an impinging jet, for which no solution was available previously. The procedure has also been applied to the square cavity at Reynolds numbers below 400; it gives results of an accuracy comparable with that of previous solutions, but with a smaller computing time.

Author(s):  
Carl E. Rathmann

For well over 150 years now, theoreticians and practitioners have been developing and teaching students easily visualized models of fluid behavior that distinguish between the laminar and turbulent fluid regimes. Because of an emphasis on applications, perhaps insufficient attention has been paid to actually understanding the mechanisms by which fluids transition between these regimes. Summarized in this paper is the product of four decades of research into the sources of these mechanisms, at least one of which is a direct consequence of the non-linear terms of the Navier-Stokes equation. A scheme utilizing chaotic dynamic effects that become dominant only for sufficiently high Reynolds numbers is explored. This paper is designed to be of interest to faculty in the engineering, chemistry, physics, biology and mathematics disciplines as well as to practitioners in these and related applications.


2013 ◽  
Vol 734 ◽  
pp. 275-316 ◽  
Author(s):  
Rashad Moarref ◽  
Ati S. Sharma ◽  
Joel A. Tropp ◽  
Beverley J. McKeon

AbstractWe study the Reynolds-number scaling and the geometric self-similarity of a gain-based, low-rank approximation to turbulent channel flows, determined by the resolvent formulation of McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382), in order to obtain a description of the streamwise turbulence intensity from direct consideration of the Navier–Stokes equations. Under this formulation, the velocity field is decomposed into propagating waves (with single streamwise and spanwise wavelengths and wave speed) whose wall-normal shapes are determined from the principal singular function of the corresponding resolvent operator. Using the accepted scalings of the mean velocity in wall-bounded turbulent flows, we establish that the resolvent operator admits three classes of wave parameters that induce universal behaviour with Reynolds number in the low-rank model, and which are consistent with scalings proposed throughout the wall turbulence literature. In addition, it is shown that a necessary condition for geometrically self-similar resolvent modes is the presence of a logarithmic turbulent mean velocity. Under the practical assumption that the mean velocity consists of a logarithmic region, we identify the scalings that constitute hierarchies of self-similar modes that are parameterized by the critical wall-normal location where the speed of the mode equals the local turbulent mean velocity. For the rank-1 model subject to broadband forcing, the integrated streamwise energy density takes a universal form which is consistent with the dominant near-wall turbulent motions. When the shape of the forcing is optimized to enforce matching with results from direct numerical simulations at low turbulent Reynolds numbers, further similarity appears. Representation of these weight functions using similarity laws enables prediction of the Reynolds number and wall-normal variations of the streamwise energy intensity at high Reynolds numbers (${Re}_{\tau } \approx 1{0}^{3} {\unicode{x2013}} 1{0}^{10} $). Results from this low-rank model of the Navier–Stokes equations compare favourably with experimental results in the literature.


1997 ◽  
Vol 119 (3) ◽  
pp. 556-567 ◽  
Author(s):  
Mihai Arghir ◽  
Jean Freˆne

The paper presents a method to calculate the rotordynamic coefficients of circumferentially-grooved liquid seals operating in centered position and turbulent flow regimes. The method is based on the integration of the averaged Navier-Stokes equations and uses a coordinate transformation proposed by Dietzen and Nordmann (1987). The effect of the coordinate transformation on the components of the stress tensor is included in the first order transport equations. To ensure grid independent solutions, numerical boundary conditions for the first-order velocities were formulated using the logarithmic law. The perturbation of the turbulent viscosity was also considered. A pressure recovery effect at the exit section was included in the first order mathematical model. The method is validated by calculations for straight and circumferentially-grooved seals. Comparisons with experimental and theoretical results show a good agreement for straight seals and for seals with few grooves, and a reasonable agreement for severe industrial cases (high Reynolds numbers and large number of grooves).


Sign in / Sign up

Export Citation Format

Share Document