Analysis of pipe flow with free surface. Part I. Numerical computation

1994 ◽  
Vol 13 (5) ◽  
pp. 229-247 ◽  
Author(s):  
Amane Tanaka ◽  
Ryuji Takaki
2016 ◽  
Vol 138 (11) ◽  
Author(s):  
Samuel Irvine ◽  
Luke Fullard

In this work, we examine the effect of wall slip for a gravity-driven flow of a Newtonian fluid in a partially filled circular pipe. An analytical solution is available for the no-slip case, while a numerical method is used for the case of flow with wall slip. We note that the partially filled circular pipe flow contains a free surface. The solution to the Navier–Stokes equations in such a case is a symmetry of a pipe flow (with no free surface) with the free surface as the symmetry plane. Therefore, we note that the analytical solution for the partially filled case is also the exact solution for fully filled lens and figure 8 shaped pipes, depending on the fill level. We find that the presence of wall slip increases the optimal fill height for maximum volumetric flow rate, brings the “velocity dip” closer to the free surface, and increases the overall flow rate for any fill. The applications of the work are twofold; the analytical solution may be used to verify numerical schemes for flows with a free surface in partially filled circular pipes, or for pipe flows in lens and figure 8 shaped pipes. Second, the work suggests that flows in pipes, particularly shallow filled pipes, can be greatly enhanced in the presence of wall slip, and optimal fill levels must account for the slip phenomenon when present.


1971 ◽  
Vol 50 (4) ◽  
pp. 645-655 ◽  
Author(s):  
M. D. Hughes ◽  
J. H. Gerrard

Flow visualization has been used quantitatively to determine the flow relative to a piston and a free surface started from rest. The discharge of water from a cylindrical reservoir was investigated. Flow with a free surface started from rest was found to have a critical Reynolds number (based on tube diameter and surface speed) of about 450 above which a ring vortex was produced just below the surface.Measurements at Reynolds numbers of 525 and 1200 were compared with computations made by the methods described in Part 1. The computed drift of tracer particles agreed well with observed values. The largest discrepancies occurred in the radial component of the drift in the early stages of the motion and amounted to 2½% of the tube diameter.


2018 ◽  
Vol 848 ◽  
pp. 467-507 ◽  
Author(s):  
Henry C.-H. Ng ◽  
Hope L. F. Cregan ◽  
Jonathan M. Dodds ◽  
Robert J. Poole ◽  
David J. C. Dennis

Pressure-driven laminar and turbulent flow in a horizontal partially filled pipe was investigated using stereoscopic particle imaging velocimetry (S-PIV) in the cross-stream plane. Laminar flow velocity measurements are in excellent agreement with a recent theoretical solution in the literature. For turbulent flow, the flow depth was varied independently of a nominally constant Reynolds number (based on hydraulic diameter, $D_{H}$; bulk velocity, $U_{b}$ and kinematic viscosity $\unicode[STIX]{x1D708}$) of $Re_{H}=U_{b}D_{H}/\unicode[STIX]{x1D708}\approx 30\,000\pm 5\,\%$. When running partially full, the inferred friction factor is no longer a simple function of Reynolds number, but also depends on the Froude number $Fr=U_{b}/\sqrt{gD_{m}}$ where $g$ is gravitational acceleration and $D_{m}$ is hydraulic mean depth. S-PIV measurements in turbulent flow reveal the presence of secondary currents which causes the maximum streamwise velocity to occur below the free surface consistent with results reported in the literature for rectangular cross-section open channel flows. Unlike square duct and rectangular open channel flow the mean secondary motion observed here manifests only as a single pair of vortices mirrored about the vertical bisector and these rollers, which fill the half-width of the pipe, remain at a constant distance from the free surface even with decreasing flow depth for the range of depths tested. Spatial distributions of streamwise Reynolds normal stress and turbulent kinetic energy exhibit preferential arrangement rather than having the same profile around the azimuth of the pipe as in a full pipe flow. Instantaneous fields reveal the signatures of elements of canonical wall-bounded turbulent flows near the pipe wall such as large-scale and very-large-scale motions and associated hairpin packets whilst near the free surface, the signatures of free surface turbulence in the absence of imposed mean shear such as ‘upwellings’, ‘downdrafts’ and ‘whirlpools’ are present. Two-point spatio-temporal correlations of streamwise velocity fluctuation suggest that the large-scale coherent motions present in full pipe flow persist in partially filled pipes but are compressed and distorted by the presence of the free surface and mean secondary motion.


2016 ◽  
Author(s):  
Hao Chen ◽  
Erik Damgaard Christensen

This paper presents numerical computation of added mass and damping coefficients of a slender horizontal cylinder in the free surface zone, which typically serves as a fish cage floater. A fully viscous two phase flow solver in OpenFOAM was employed in the numerical computation. The purpose was to validate the capability of this solver and dynamic mesh functionality. A two dimensional numerical wave tank was set up, and two wave relaxation zones were used to reduce the size of the computational domain. Harmonic forced oscillations of the cylinder were performed at different frequencies and amplitudes. The mesh at free surface zone was refined based on the radiated wave heights at different oscillation frequencies in order to properly resolve the radiated waves. The result shows that in most frequency ranges, the numerical computation agreed well with the experimental data and analytical solution. However at low frequency range for added mass coefficient in heave motion, deviations were observed, and it was due to the effect of finite water depth. In addition for sway motion at high frequency range, the damping coefficient was underestimated comparing with analytical solution. This was believed to be as a result of high steepness of the radiated waves.


Sign in / Sign up

Export Citation Format

Share Document