Supercritical withdrawal from a two-layer fluid through a line sink

1995 ◽  
Vol 297 ◽  
pp. 37-47 ◽  
Author(s):  
G. C. Hocking

Accurate numerical solutions to the problem of finding the location of the interface between two unconfined regions of fluid of different density during the withdrawal process are presented. Supercritical flows are considered, in which the interface is drawn directly into the sink. As the flow rate is reduced, the interface enters the sink more steeply, until the solution method breaks down just before the interface enters the sink vertically from above, and becomes flow from the lower layer only. This lower bound on supercritical flow is compared with the upper bound on single-layer (free surface) flow with good agreement.

2005 ◽  
Vol 127 (3) ◽  
pp. 572-582 ◽  
Author(s):  
Shin Hyung Rhee

The present study is concerned with liquid tank sloshing at low filling level conditions. The volume of fluid method implemented in a Navier–Stokes computational fluid dynamics code is employed to handle the free-surface flow of liquid sloshing. The geometric reconstruction scheme for the interface representation is employed to ensure sharpness at the free surface. The governing equations are discretized by second order accurate schemes on unstructured grids. Several different computational approaches are verified and numerical uncertainties are assessed. The computational results are validated against existing experimental data, showing good agreement. The capability is demonstrated for a generic membrane-type liquefied natural gas carrier tank with a simplified pump tower inside. The validation results suggest that the present computational approach is both easy to apply and accurate enough for more realistic problems.


2012 ◽  
Vol 23 (4) ◽  
pp. 441-467 ◽  
Author(s):  
CHRISTOPHER J. LUSTRI ◽  
SCOTT W. MCCUE ◽  
BENJAMIN J. BINDER

The problem of steady subcritical free surface flow past a submerged inclined step is considered. The asymptotic limit of small Froude number is treated, with particular emphasis on the effect that changing the angle of the step face has on the surface waves. As demonstrated by Chapman & Vanden-Broeck, (2006) Exponential asymptotics and gravity waves. J. Fluid Mech.567, 299–326, the divergence of a power series expansion in powers of the square of the Froude number is caused by singularities in the analytic continuation of the free surface; for an inclined step, these singularities may correspond to either the corners or stagnation points of the step, or both, depending on the angle of inclination. Stokes lines emanate from these singularities, and exponentially small waves are switched on at the point the Stokes lines intersect with the free surface. Our results suggest that for a certain range of step angles, two wavetrains are switched on, but the exponentially subdominant one is switched on first, leading to an intermediate wavetrain not previously noted. We extend these ideas to the problem of flow over a submerged bump or trench, again with inclined sides. This time there may be two, three or four active Stokes lines, depending on the inclination angles. We demonstrate how to construct a base topography such that wave contributions from separate Stokes lines are of equal magnitude but opposite phase, thus cancelling out. Our asymptotic results are complemented by numerical solutions to the fully nonlinear equations.


2020 ◽  
Vol 61 ◽  
pp. C75-C88
Author(s):  
Shaymaa Mukhlif Shraida ◽  
Graeme Hocking

We consider the outflow of water from the peak of a triangular ridge into a channel of finite depth. Solutions are computed for different flow rates and bottom angles. A numerical method is used to compute the flow from the source for small values of flow rate and it is found that there is a maximum flow rate beyond which steady solutions do not seem to exist. Limiting flows are computed for each geometrical configuration. One application of this work is as a model of saline water being returned to the ocean after desalination. References Craya, A. ''Theoretical research on the flow of nonhomogeneous fluids''. La Houille Blanche, (1):22–55, 1949. doi:10.1051/lhb/1949017 Dun, C. R. and Hocking, G. C. ''Withdrawal of fluid through a line sink beneath a free surface above a sloping boundary''. J. Eng. Math. 29:1–10, 1995. doi:10.1007/bf00046379 Hocking, G. ''Cusp-like free-surface flows due to a submerged source or sink in the presence of a flat or sloping bottom''. ANZIAM J. 26:470–486, 1985. doi:10.1017/s0334270000004665 Hocking, G. C. and Forbes, L. K. ''Subcritical free-surface flow caused by a line source in a fluid of finite depth''. J. Eng. Math. 26:455-466, 1992. doi:10.1007/bf00042763 Hocking, G. C. ''Supercritical withdrawal from a two-layer fluid through a line sink", J. Fluid Mech. 297:37–47, 1995. doi:10.1017/s0022112095002990 Hocking, G. C., Nguyen, H. H. N., Forbes, L. K. and Stokes,T. E. ''The effect of surface tension on free surface flow induced by a point sink''. ANZIAM J., 57:417–428, 2016. doi:10.1017/S1446181116000018 Landrini, M. and Tyvand, P. A. ''Generation of water waves and bores by impulsive bottom flux'', J. Eng. Math. 39(1–4):131-170, 2001. doi:10.1023/A:1004857624937 Lustri, C. J., McCue, S. W. and Chapman, S. J. ''Exponential asymptotics of free surface flow due to a line source''. IMA J. Appl. Math., 78(4):697–713, 2013. doi:10.1093/imamat/hxt016 Stokes, T. E., Hocking, G. C. and Forbes, L.K. ''Unsteady free surface flow induced by a line sink in a fluid of finite depth'', Comp. Fluids, 37(3):236–249, 2008. doi:10.1016/j.compfluid.2007.06.002 Tuck, E. O. and Vanden-Broeck, J.-M. ''A cusp-like free-surface flow due to a submerged source or sink''. ANZIAM J. 25:443–450, 1984. doi:10.1017/s0334270000004197 Vanden-Broeck, J.-M., Schwartz, L. W. and Tuck, E. O. ''Divergent low-Froude-number series expansion of nonlinear free-surface flow problems". Proc. Roy. Soc. A., 361(1705):207–224, 1978. doi:10.1098/rspa.1978.0099 Vanden-Broeck, J.-M. and Keller, J. B. ''Free surface flow due to a sink'', J. Fluid Mech, 175:109–117, 1987. doi:10.1017/s0022112087000314 Yih, C.-S. Stratified flows. Academic Press, New York, 1980. doi:10.1016/B978-0-12-771050-1.X5001-3


1993 ◽  
Vol 254 ◽  
pp. 649-670 ◽  
Author(s):  
S. B. G. M. O'Brien

In the field of industrial drying, a recent innovation has exploited the occurrence of Marangoni effects in such a way that the resultant free-surface flow enhances the drying process. To this end, alcohol vapour, soluble in water, is introduced above a drying film and as a result of diffusion through the air and water phases a favourable concentration gradient gives rise to the required shear flow. We consider here a simple process driven by this mechanism, and by means of asymptotic simplification and the concepts of singular perturbation theory a leading-order approximation is obtained in which the alcohol concentration in the water is a specified function of space and time. The evolution of the free surface thus reduces to a single nonlinear partial differential equation of a similar form to the Korteweg–de Vries and Burgers equations, higher-derivative terms corresponding to surface tension and gravity effects. Numerical solutions of this equation are obtained and are compared to the application of first order nonlinear kinematic wave theory with corresponding shock solutions.


1999 ◽  
Vol 392 ◽  
pp. 73-100 ◽  
Author(s):  
M. WIELAND ◽  
J. M. N. T. GRAY ◽  
K. HUTTER

A series of laboratory experiments and numerical simulations have been performed to investigate the rapid fluid-like flow of a finite mass of granular material down a chute with partial lateral confinement. The chute consists of a section inclined at 40° to the horizontal, which is connected to a plane run-out zone by a smooth transition. The flow is confined on the inclined section by a shallow parabolic cross-slope profile. Photogrammetric techniques have been used to determine the position of the evolving boundary during the flow, and the free-surface height of the stationary granular deposit in the run-out zone. The results of three experiments with different granular materials are presented and shown to be in very good agreement with numerical simulations based on the Savage–Hutter theory for granular avalanches. The basal topography over which the avalanche flows has a strong channelizing effect on the inclined section of the chute. As the avalanche reaches the run-out zone, where the lateral confinement ceases, the head spreads out to give the avalanche a characteristic ‘tadpole’ shape. Sharp gradients in the avalanche thickness and velocity began to develop at the interface between the nose and tail of the avalanche as it came to rest, indicating that a shock wave develops close to the end of the experiments.


1995 ◽  
Vol 284 ◽  
pp. 63-96 ◽  
Author(s):  
Robert I. Bowles

Steady planar flow of a liquid layer over an obstacle is studied for favourable slopes. First, half-plane Poiseuille flow is found to be a non-unique solution on a uniformly sloping surface since eigensolutions exist which are initially exponentially small far upstream. These have their origin in a viscous–inviscid interaction between the retarding action of viscosity and the hydrostatic pressure from the free surface. The cross-stream pressure gradient caused by the curvature of the streamlines also comes into play as the slope increases. As the interaction becomes nonlinear, separation of the liquid layer can occur, of a breakaway type if the slope is sufficiently large. The breakaway represents a hydraulic jump in the sense of a localized relatively short-scaled increase in layer thickness, e.g. far upstream of a large obstacle. The solution properties give predictions for the shape and structure of hydraulic jumps on various slopes. Secondly, the possibility of standing waves downstream of the jump is addressed for various slope magnitudes. A limiting case of small gradient, governed by lubrication theory, allows the downstream boundary condition to be included explicitly. Numerical solutions showing the free-surface flow over an obstacle confirm the analytical conclusions. In addition the predictions are compared with the experimental and computational results of Pritchardet al.(1992), yielding good qualitative and quantitative agreement. The effects of surface tension on the jump are also discussed and in particular the free interaction on small slopes is examined for large Bond numbers.


Author(s):  
Shin Hyung Rhee

The present study is concerned with the liquid tank sloshing at low filling level conditions. The volume of fluid method implemented in a Navier-Stokes computational fluid dynamics code is employed to handle the free-surface flow of liquid sloshing. A geometric reconstruction scheme for the interface representation is employed to ensure sharpness at the free-surface. The governing equations are discretized by second order accurate schemes on unstructured grids. Several different computational approaches are verified and numerical uncertainties are assessed. The computational results are validated against existing experimental data, showing good agreement. The capability is demonstrated for a generic membrane type LNG carrier tank with a simplified pump tower inside. The validation results suggest that the present computational approach is both easy to apply and accurate enough for more realistic problems.


2016 ◽  
Vol 808 ◽  
pp. 441-468 ◽  
Author(s):  
S. L. Gavrilyuk ◽  
V. Yu. Liapidevskii ◽  
A. A. Chesnokov

A two-layer long-wave approximation of the homogeneous Euler equations for a free-surface flow evolving over mild slopes is derived. The upper layer is turbulent and is described by depth-averaged equations for the layer thickness, average fluid velocity and fluid turbulent energy. The lower layer is almost potential and can be described by Serre–Su–Gardner–Green–Naghdi equations (a second-order shallow water approximation with respect to the parameter $H/L$, where $H$ is a characteristic water depth and $L$ is a characteristic wavelength). A simple model for vertical turbulent mixing is proposed governing the interaction between these layers. Stationary supercritical solutions to this model are first constructed, containing, in particular, a local turbulent subcritical zone at the forward slope of the wave. The non-stationary model was then numerically solved and compared with experimental data for the following two problems. The first one is the study of surface waves resulting from the interaction of a uniform free-surface flow with an immobile wall (the water hammer problem with a free surface). These waves are sometimes called ‘Favre waves’ in homage to Henry Favre and his contribution to the study of this phenomenon. When the Froude number is between 1 and approximately 1.3, an undular bore appears. The characteristics of the leading wave in an undular bore are in good agreement with experimental data by Favre (Ondes de Translation dans les Canaux Découverts, 1935, Dunod) and Treske (J. Hydraul Res., vol. 32 (3), 1994, pp. 355–370). When the Froude number is between 1.3 and 1.4, the transition from an undular bore to a breaking (monotone) bore occurs. The shoaling and breaking of a solitary wave propagating in a long channel (300 m) of mild slope (1/60) was then studied. Good agreement with experimental data by Hsiao et al. (Coast. Engng, vol. 55, 2008, pp. 975–988) for the wave profile evolution was found.


Author(s):  
Aggelos S. Dimakopoulos ◽  
Athanassios A. Dimas

A numerical model is presented for the simulation of the two-dimensional, inviscid, free-surface flow developing by the propagation and breaking of water waves over a flat bottom of steep slope. The simulation is based on the numerical solution of the unsteady, two-dimensional, Euler equations subject to the fully-nonlinear free-surface boundary conditions, the non-penetration condition at the bottom and appropriate inflow and outflow conditions. A boundary-fitted transformation, which includes both the time-dependent free surface and the arbitrary bottom shape, is applied. For the numerical solution of the Euler equations, a two-stage fractional time-step method is employed for the temporal discretization, while a hybrid scheme is used for the spatial discretization. Finite differences are used in the streamwise direction and a pseudo-spectral method in the vertical direction. An absorption zone is placed at the outflow region in order to minimize wave reflection by the outflow boundary. Wave breaking is modeled by a surface roller breaking model, which modifies the dynamic free-surface condition. The simulation results are in very good agreement with available experimental results for the wave propagation and breaking over bottom with slope 1:35. Results, from the simulations over bottom with steeper slopes of 1:15 and 1:10, which generate strong spilling and mild plunging breakers, respectively, are also in very good agreement with available predictions for the breaking depth and wave height. In all cases, a vortex is formed under the breaking wave front and convected in the surf zone.


Sign in / Sign up

Export Citation Format

Share Document