Orifice Contraction Coefficient for Inviscid Incompressible Flow

1985 ◽  
Vol 107 (1) ◽  
pp. 36-43 ◽  
Author(s):  
R. D. Grose

The theory for steady flow of an incompressible fluid through an orifice has been semi-empirically established for only certain flow conditions. In this paper, the development of a more rigorous theory for the prediction of the orifice flow contraction effect is presented. This theory is based on the conservation of momentum and mass principles applied to global control volumes for continuum flow. The control volumes are chosen to have a particular geometric construction which is based on certain characteristics of the Navier-Stokes equations for incompressible and, in the limit, inviscid flow. The treatment is restricted to steady incompressible, single phase, single component, inviscid Newtonian flow, but the principles that are developed hold for more general conditions. The resultant equations predict the orifice contraction coefficient as a function of the upstream geometry ratio for both axisymmetric and two-dimensional flow fields. The predicted contraction coefficient values agree with experimental orifice discharge coefficient data without the need for empirical adjustment.

Author(s):  
Markus Scholle ◽  
André Haas ◽  
Philip H. Gaskell

Although it is well known that Bernoulli's equation is obtained as the first integral of Euler's equations in the absence of vorticity and that in the case of non-vanishing vorticity a first integral of them can be found using the Clebsch transformation for inviscid flow, generalization of the procedure for viscous flow has remained elusive. Accordingly, in this paper, a first integral of the Navier–Stokes equations for steady flow is constructed. In the case of a two-dimensional flow, this is made possible by formulating the governing equations in terms of complex variables and introducing a new scalar potential. Associated boundary conditions are considered, and an extension of the theory to three dimensions is proposed. The capabilities of the new approach are demonstrated by calculating a Reynolds number correction to the laminar shear flow generated in the narrow gap between a flat moving and a stationary wavy wall, as is often encountered in lubrication problems. It highlights the first integral as a suitable tool for the development of new analytical and numerical methods in fluid dynamics.


Author(s):  
S. G. Rajeev

The initial value problem of the incompressible Navier–Stokes equations is explained. Leray’s classic study of it (using Picard iteration) is simplified and described in the language of physics. The ideas of Lebesgue and Sobolev norms are explained. The L2 norm being the energy, cannot increase. This gives sufficient control to establish existence, regularity and uniqueness in two-dimensional flow. The L3 norm is not guaranteed to decrease, so this strategy fails in three dimensions. Leray’s proof of regularity for a finite time is outlined. His attempts to construct a scale-invariant singular solution, and modern work showing this is impossible, are then explained. The physical consequences of a negative answer to the regularity of Navier–Stokes solutions are explained. This chapter is meant as an introduction, for physicists, to a difficult field of analysis.


1994 ◽  
Vol 61 (3) ◽  
pp. 629-633 ◽  
Author(s):  
S. H. Smith

When a stretching surface is moved quickly, for a short period of time, a pulse is transmitted to the surrounding fluid. Here we describe an exact solution in terms of a similarity variable for the Navier-Stokes equations which represents the effect of this pulse for two-dimensional flow. The unusual feature is that this solution is only valid for a limited range of the Reynolds number; outside this domain unbounded velocities result.


2010 ◽  
Vol 132 (10) ◽  
Author(s):  
Minoru Watari

Relationship between accuracy and number of velocity particles in velocity slip phenomena was investigated by numerical simulations and theoretical considerations. Two types of 2D models were used: the octagon family and the D2Q9 model. Models have to possess the following four prerequisites to accurately simulate the velocity slip phenomena: (a) equivalency to the Navier–Stokes equations in the N-S flow area, (b) conservation of momentum flow Pxy in the whole area, (c) appropriate relaxation process in the Knudsen layer, and (d) capability to properly express the mass and momentum flows on the wall. Both the octagon family and the D2Q9 model satisfy conditions (a) and (b). However, models with fewer velocity particles do not sufficiently satisfy conditions (c) and (d). The D2Q9 model fails to represent a relaxation process in the Knudsen layer and shows a considerable fluctuation in the velocity slip due to the model’s angle to the wall. To perform an accurate velocity slip simulation, models with sufficient velocity particles, such as the triple octagon model with moving particles of 24 directions, are desirable.


1991 ◽  
Vol 227 ◽  
pp. 1-33 ◽  
Author(s):  
Stephen M. Cox

We consider the flow of a viscous incompressible fluid in a parallel-walled channel, driven by steady uniform suction through the porous channel walls. A similarity transformation reduces the Navier-Stokes equations to a single partial differential equation (PDE) for the stream function, with two-point boundary conditions. We discuss the bifurcations of the steady solutions first, and show how a pitchfork bifurcation is unfolded when a symmetry of the problem is broken.Then we describe time-dependent solutions of the governing PDE, which we calculate numerically. We analyse these unsteady solutions when there is a high rate of suction through one wall, and the other wall is impermeable: there is a limit cycle composed of an explosive phase of inviscid growth, and a slow viscous decay. The inviscid phase ‘almost’ has a finite-time singularity. We discuss whether solutions of the governing PDE, which are exact solutions of the Navier-Stokes equations, may develop mathematical singularities in a finite time.When the rates of suction at the two walls are equal so that the problem is symmetrical, there is an abrupt transition to chaos, a ‘homoclinic explosion’, in the time-dependent solutions as the Reynolds number is increased. We unfold this transition by perturbing the symmetry, and compare direct numerical integrations of the governing PDE with a recent theory for ‘Lorenz-like’ dynamical systems. The chaos is found to be very sensitive to symmetry breaking.


2005 ◽  
Vol 47 (1) ◽  
pp. 39-50 ◽  
Author(s):  
G. D. McBain

AbstractWe continue our study of the adaptation from spherical to doubly periodic slot domains of the poloidal-toroidal representation of vector fields. Building on the successful construction of an orthogonal quinquepartite decomposition of doubly periodic vector fields of arbitrary divergence with integral representations for the projections of known vector fields and equivalent scalar representations for unknown vector fields (Part 1), we now present a decomposition of vector field equations into an equivalent set of scalar field equations. The Stokes equations for slow viscous incompressible fluid flow in an arbitrary force field are treated as an example, and for them the application of the decomposition uncouples the conservation of momentum equation from the conservation of mass constraint. The resulting scalar equations are then solved by elementary methods. The extension to generalised Stokes equations resulting from the application of various time discretisation schemes to the Navier-Stokes equations is also solved.


1985 ◽  
Vol 154 ◽  
pp. 357-375 ◽  
Author(s):  
J. A. C. Humphrey ◽  
H. Iacovides ◽  
B. E. Launder

The paper reports numerical solutions to a semi-elliptic truncation of the Navier–Stokes equations for the case of developing laminar flow in circular-sectioned bends over a range of Dean numbers. The ratios of bend radius to pipe radius are 7:1 and 20:1, corresponding with the configurations examined experimentally by Talbot and his co-workers in recent years. The semi-elliptic treatment facilitates a much finer grid than has been possible in earlier studies. Numerical accuracy has been further improved by assuming radial equilibrium over a thin sublayer immediately adjacent to the wall and by re-formulating the boundary conditions at the pipe centre.Streamwise velocity profiles at Dean numbers of 183 and 565 are in excellent agreement with laser-Doppler measurements by Agrawal, Talbot & Gong (1978). Good, albeit less complete, accord is found with the secondary velocities, though the differences that exist may be mainly due to the difficulty of making these measurements. The paper provides new information on the behaviour of the streamwise shear stress around the inner line of symmetry. Upstream of the point of minimum shear stress, our numerical predictions display a progressive shift towards the result of Stewartson, Cebici & Chang (1980) as the Dean number is successively raised. Downstream of the minimum, however, in contrast with the monotonic approach to an asymptotic level reported by Stewartson, the numerical solutions display a damped oscillatory behaviour reminiscent of those from Hawthorne's (1951) inviscid-flow calculations. The amplitude of the oscillation grows as the Dean number is raised.


Author(s):  
Angelo Frisani ◽  
Yassin A. Hassan

A projection approach is presented for the coupled system of time-dependent incompressible Navier-Stokes equations in conjunction with the Immersed Boundary Method (IBM) for solving fluid flow problems in the presence of rigid objects not represented by the underlying mesh. The IBM allows solving the flow for geometries with complex objects without the need of generating a body fitted mesh. The no-slip boundary constraint is satisfied applying a boundary force at the immersed body surface. Using projection and interpolation operators from the fluid volume mesh to the solid surface mesh (i.e., the “immersed” boundary) and vice versa, it is possible to impose the extra constraint to the incompressible Navier-Stokes equations as a Lagrange multiplier in a fashion very similar to the effect pressure has on the momentum equations to satisfy the divergence-free constraint. The projection operation removes the immersed boundary surface slip and non-divergence-free components of the velocity field. The boundary force is determined implicitly at the inner iterations of the fractional step method implemented. No constitutive relations for the immersed boundary objects fluid interaction are required, allowing the formulation introduced to use larger CFL numbers compared to previous methodologies. An overview of the immersed boundary approach is presented showing third order accuracy in space and second order accuracy in time when the simulation results for the Taylor-Green decaying vortex are compared to the analytical solution using the Immersed Finite Element Method (IFEM). For the Immersed Finite Volume Method (IFVM) a ghost-cell approach is used. Second order accuracy in space and first order accuracy in time are obtained when the Taylor-Green decaying vortex test case is compared to the analytical solution. The numerical results are compared with the analytical solution also for adaptive mesh refinement (for the IFEM) showing an excellent error reduction. Computations were performed using IFEM and IFVM approaches for the time-dependent incompressible Navier-Stokes equations in a two-dimensional flow past a stationary circular cylinder at Re = 20, and 40, where shedding effects are not present. The drag coefficient and the recirculation length error compared to the experimental data is less than 3–4%. Simulations for the two-dimensional flow past a stationary circular cylinder at Re = 100 were also performed. For Re numbers above 46, unsteadiness generates vortex shedding, and an unsteady flow regime is present. The results shown are in excellent quantitative and qualitative agreement with the flow pattern expected. The numerical results obtained with the discussed IFEM and IFVM were also compared against other immersed boundary methodologies available in literature and simulation performed with the commercial computational fluid dynamics code STAR-CCM+/V5.02.009 for which a body fitted finite volume numerical discretization was used. The benchmark showed that the numerical results obtained with the implemented immersed boundary methods are very close to those obtained from STAR-CCM+ with a very fine mesh and in a good agreement with the other IBM techniques. The IBM based of finite element approach is numerically more accurate than the IBM based on finite volume discretization. In contrast, the latter is computationally more efficient than the former.


Sign in / Sign up

Export Citation Format

Share Document