scholarly journals A Potential Field Description for Gravity-Driven Film Flow over Piece-Wise Planar Topography

Fluids ◽  
2019 ◽  
Vol 4 (2) ◽  
pp. 82 ◽  
Author(s):  
Markus Scholle ◽  
Philip H. Gaskell ◽  
Florian Marner
Keyword(s):  
Free Surface ◽  
Potential Field ◽  
Numerical Solutions ◽  
Film Flow ◽  
First Integral ◽  
Characteristic Parameter ◽  
Numerical Predictions ◽  
Planar Substrate ◽  
Surface Disturbance ◽  
Gravity Driven

Models based on a potential field description and corresponding first integral formulation, embodying a reduction of the associated dynamic boundary condition at a free surface to one of a standard Dirichlet-Neumann type, are used to explore the problem of continuous gravity-driven film flow down an inclined piece-wise planar substrate in the absence of inertia. Numerical solutions of the first integral equations are compared with analytical ones from a linearised form of a reduced equation set resulting from application of the long-wave approximation. The results obtained are shown to: (i) be in very close agreement with existing, comparable experimental data and complementary numerical predictions for isolated step-like topography available in the open literature; (ii) exhibit the same qualitative behaviour for a range of Capillary numbers and step heights/depths, becoming quantitively similar when both are small. A novel outcome of the formulation adopted is identification of an analytic criteria enabling a simple classification procedure for specifying the characteristic nature of the free surface disturbance formed; leading subsequently to the generation of a related, practically relevant, characteristic parameter map in terms of the substrate inclination angle and the Capillary number of the associated flow.

Gravity-driven creeping flow of two adjacent layers through a channel and down a plane wall

1998 ◽  
Vol 371 ◽  
pp. 345-376 ◽  
Author(s):  
C. POZRIKIDIS
Keyword(s):  
Free Surface ◽  
Numerical Simulations ◽  
Channel Flow ◽  
Film Flow ◽  
Finite Amplitude ◽  
Creeping Flow ◽  
Plane Wall ◽  
Surface Velocity ◽  
Inclined Plane ◽  
Gravity Driven

We study the stability of the interface between (a) two adjacent viscous layers flowing due to gravity through an inclined or vertical channel that is confined between two parallel plane walls, and (b) two superimposed liquid films flowing down an inclined or vertical plane wall, in the limit of Stokes flow. In the case of channel flow, linear stability analysis predicts that, when the fluids are stably stratified, the flow is neutrally stable when the surface tension vanishes and the channel is vertical, and stable otherwise. This behaviour contrasts with that of the gravity-driven flow of two superimposed films flowing down an inclined plane, where an instability has been identified when the viscosity of the fluid next to the plane is less than that of the top fluid, even in the absence of fluid inertia. We investigate the nonlinear stages of the motion subject to finite-amplitude two-dimensional perturbations by numerical simulations based on boundary-integral methods. In both cases of channel and film flow, the mathematical formulation results in integral equations for the unknown interface and free-surface velocity. The properties of the integral equation for multi-film flow are investigated with reference to the feasibility of computing a solution by the method of successive substitutions, and a deflation strategy that allows an iterative procedure is developed. In the case of channel flow, the numerical simulations show that disturbances of sufficiently large amplitude may cause permanent deformation in which the interface folds or develops elongated fingers. The ratio of the viscosities and densities of the two fluids plays an important role in determining the morphology of the emerging interfacial patterns. Comparing the numerical results with the predictions of a model based on the lubrication approximation shows that the simplified approach can only describe a limited range of motions. In the case of film flow down an inclined plane, we develop a method for extracting the properties of the normal modes, including the ratio of the amplitudes of the free-surface and interfacial waves and their relative phase lag, from the results of a numerical simulation for small deformations. The numerical procedure employs an adaptation of Prony's method for fitting a signal described by a time series to a sum of complex exponentials; in the present case, the signal is identified with the cosine or sine Fourier coefficients of the interface and free-surface waves. Numerical simulations of the nonlinear motion confirm that the deformability of the free surface is necessary for the growth of small-amplitude perturbations, and show that the morphology of the interfacial patterns developing subject to finite-amplitude perturbations is qualitatively similar to that for channel flow.

scholarly journals Virtual Replica of a Towing Tank Experiment to Determine the Kelvin Half-Angle of a Ship in Restricted Water

2020 ◽  
Vol 8 (4) ◽  
pp. 258 ◽  
Author(s):  
Momchil Terziev ◽  
Guangwei Zhao ◽  
Tahsin Tezdogan ◽  
Zhiming Yuan ◽  
Atilla Incecik
Keyword(s):  
Free Surface ◽  
Theoretical Method ◽  
Free Surface Flows ◽  
Surrounding Water ◽  
Towing Tank ◽  
Numerical Predictions ◽  
Galilean Relativity ◽  
Elevation Data ◽  
Half Angle ◽  
Surface Disturbance

The numerical simulation of ship flows has evolved into a highly practical approach in naval architecture. In typical virtual towing tanks, the principle of Galilean relativity is invoked to maintain the ship as fixed, while the surrounding water is prescribed to flow past it. This assumption may be identified, at least partly, as being responsible for the wide-scale adoption of computational solutions within practitioners’ toolkits. However, it carries several assumptions, such as the levels of inlet turbulence and their effect on flow properties. This study presents an alternative virtual towing tank, where the ship is simulated to advance over a stationary fluid. To supplement the present work, the free surface disturbance is processed into Fourier space to determine the Kelvin half-angle for an example case. The results suggest that it is possible to construct a fully unsteady virtual towing tank using the overset method, without relying on Galilean relativity. Differences between theoretical and numerical predictions for the Kelvin half-angle are predominantly attributed to the assumptions used by the theoretical method. The methods presented in this work can potentially be used to validate free-surface flows, even when one does not have access to experimental wave elevation data.

scholarly journals Weak, strong and detached oblique shocks in gravity-driven granular free-surface flows

2007 ◽  
Vol 579 ◽  
pp. 113-136 ◽  
Author(s):  
J. M. N. T. GRAY ◽  
X. CUI
Keyword(s):  
Free Surface ◽  
Numerical Solutions ◽  
Rapid Change ◽  
Free Surface Flows ◽  
Gravitational Acceleration ◽  
Source Terms ◽  
Surface Flows ◽  
Gravity Driven ◽  
Downstream Flow ◽  
Weak Shocks

Hazardous natural flows such as snow-slab avalanches, debris flows, pyroclastic flows and lahars are part of a much wider class of dense gravity-driven granular free-surface flows that frequently occur in industrial processes as well as in foodstuffs in our kitchens! This paper investigates the formation of oblique granular shocks, when the oncoming flow is deflected by a wall or obstacle in such a way as to cause a rapid change in the flow height and velocity. The theory for non-accelerative slopes is qualitatively similar to that of gasdynamics. For a given deflection angle there are three possibilities: a weak shock may form close to the wall; a strong shock may extend across the chute; or the shock may detach from the tip. Weak shocks have been observed in both dense granular free-surface flows and granular gases. This paper shows how strong shocks can be triggered in chute experiments by careful control of the downstream boundary conditions. The resulting downstream flow height is much thicker than that of weak shocks and there is a marked decrease in the downstream velocity. Strong shocks therefore dissipate much more energy than weak shocks. An exact solution for the angle at which the flow detaches from the wedge is derived and this is shown to be in excellent agreement with experiment. It therefore provides a very useful criterion for determining whether the flow will detach. In experimental, industrial and geophysical flows the avalanche is usually accelerated, or decelerated, by the net effect of the gravitational acceleration and basal sliding friction as the slope inclination angle changes. The presence of these source terms necessarily leads to gradual changes in the flow height and velocity away from the shocks, and this in turn modifies the local Froude number of the flow. A shock-capturing non-oscillating central method is used to compute numerical solutions to the full problem. This shows that the experiments can be matched very closely when the source terms are included and explains the deviations away from the classical oblique-shock theory. We show that weak shocks bend towards the wedge on accelerative slopes and away from it on decelerative slopes. In both cases the presence of the source terms leads to a gradual increase in the downstream flow thickness along the wedge, which suggests that defensive dams should increase in height further down the slope, contrary to current design criteria but in accordance with field observations of snow-avalanche deposits from a defensive dam in Northwestern Iceland. Movies are available with the online version of the paper.

Nonlinear Sloshing in Fixed and Vertically Excited Containers

2002 ◽  
Author(s):  
Jannette B. Frandsen ◽  
Alistair G. L. Borthwick
Keyword(s):  
Perturbation Theory ◽  
Free Surface ◽  
Small Perturbation ◽  
Numerical Solutions ◽  
Vertical Direction ◽  
Nonlinear Effects ◽  
Breaking Waves ◽  
Surface Flow ◽  
Nonlinear Behaviour ◽  
Computational Domain

Nonlinear effects of standing wave motions in fixed and vertically excited tanks are numerically investigated. The present fully nonlinear model analyses two-dimensional waves in stable and unstable regions of the free-surface flow. Numerical solutions of the governing nonlinear potential flow equations are obtained using a finite-difference time-stepping scheme on adaptively mapped grids. A σ-transformation in the vertical direction that stretches directly between the free-surface and bed boundary is applied to map the moving free surface physical domain onto a fixed computational domain. A horizontal linear mapping is also applied, so that the resulting computational domain is rectangular, and consists of unit square cells. The small-amplitude free-surface predictions in the fixed and vertically excited tanks compare well with 2nd order small perturbation theory. For stable steep waves in the vertically excited tank, the free-surface exhibits nonlinear behaviour. Parametric resonance is evident in the instability zones, as the amplitudes grow exponentially, even for small forcing amplitudes. For steep initial amplitudes the predictions differ considerably from the small perturbation theory solution, demonstrating the importance of nonlinear effects. The present numerical model provides a simple way of simulating steep non-breaking waves. It is computationally quick and accurate, and there is no need for free surface smoothing because of the σ-transformation.

NUMERICAL SOLUTIONS OF 2D AND 3D SLAMMING PROBLEMS

2021 ◽  
Vol 153 (A2) ◽  
Author(s):  
Q Yang ◽  
W Qiu
Keyword(s):  
Free Surface ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Type Equation ◽  
The Body ◽  
Navier Stokes ◽  
Time Step ◽  
Navier Stokes Equations ◽  
Highly Nonlinear ◽  
2D And 3D

Slamming forces on 2D and 3D bodies have been computed based on a CIP method. The highly nonlinear water entry problem governed by the Navier-Stokes equations was solved by a CIP based finite difference method on a fixed Cartesian grid. In the computation, a compact upwind scheme was employed for the advection calculations and a pressure-based algorithm was applied to treat the multiple phases. The free surface and the body boundaries were captured using density functions. For the pressure calculation, a Poisson-type equation was solved at each time step by the conjugate gradient iterative method. Validation studies were carried out for 2D wedges with various deadrise angles ranging from 0 to 60 degrees at constant vertical velocity. In the cases of wedges with small deadrise angles, the compressibility of air between the bottom of the wedge and the free surface was modelled. Studies were also extended to 3D bodies, such as a sphere, a cylinder and a catamaran, entering calm water. Computed pressures, free surface elevations and hydrodynamic forces were compared with experimental data and the numerical solutions by other methods.

scholarly journals Explicit Solutions of a Gravity-Induced Film Flow along a Convectively Heated Vertical Wall

2013 ◽  
Vol 2013 ◽  
pp. 1-7 ◽  
Author(s):  
Ammarah Raees ◽  
Hang Xu
Keyword(s):  
Film Flow ◽  
Vertical Wall ◽  
Boussinesq Approximation ◽  
Convective Boundary Condition ◽  
Convective Boundary ◽  
Similarity Transformations ◽  
Reduced Equations ◽  
Homotopy Analysis ◽  
Gravity Driven ◽  
The Mathematical Model

The gravity-driven film flow has been analyzed along a vertical wall subjected to a convective boundary condition. The Boussinesq approximation is applied to simplify the buoyancy term, and similarity transformations are used on the mathematical model of the problem under consideration, to obtain a set of coupled ordinary differential equations. Then the reduced equations are solved explicitly by using homotopy analysis method (HAM). The resulting solutions are investigated for heat transfer effects on velocity and temperature profiles.

Sign in / Sign up

Export Citation Format

Share Document