Invariance group properties and exact solutions of equations describing time-dependent free surface flows under gravity

Using the method of infinitesimal transformations, a 6-parameter family of exact solutions describing nonlinear sheared flows with a free surface are found. These solutions are a hybrid between the earlier self-propagating simple wave solutions of Freeman, and decaying solutions of Sachdev. Simple wave solutions are also derived via the method of infinitesimal transformations. Incomplete beta functions seem to characterize these (nonlinear) sheared flows in the absence of critical levels.

Exact, self-similar, time-dependent free surface flows under gravity

A class of exact, self-similar, time-dependent solutions describing free surface flows under gravity is found which extends the self-propagating class of solutions discovered earlier by Freeman (1972) to those which decay with time.

Exact free-surface flows for shallow-water equations .1. - the incompressible case

A comprehensive exact treatment of free surface flows governed by shallow water equations (in sigma variables) is given. Several new families of exact solutions of the governing PDEs are found and are shown to embed the well-known self-similar or traveling wave solutions which themselves are governed by reduced ODEs. The classes of solutions found here are explicit in contrast to those found earlier in an implicit form. The height of the free surface for each family of solutions is found explicitly. For the traveling or simple wave, the free surface is governed by a nonlinear wave equation, but is arbitrary otherwise. For other types of solutions, the height of the free surface is constant either on lines of constant acceleration or on lines of constant speed; in another case, the free surface is a horizontal plane while the flow underneath is a sine wave. The existence of simple waves on shear flows is analytically proved. The interaction of large amplitude progressive waves with shear flow is also studied.

Exact free-surface flows for shallow-water equations .2. the compressible case

Exact free surface flows with shear in a compressible barotropic medium are found, extending the authors' earlier work for the incompressible medium. The barotropic medium is of finite extent in the vertical direction, while it is infinite in the horizontal direction. The ''shallow water'' equations for a compressible barotropic medium, subject to boundary conditions at the free surface and at the bottom, are solved in terms of double psi-series, Simple wave and time-dependent solutions are found; for the former the free surface is of arbitrary shape while for the latter it is a damping traveling wave in the horizontal direction, For other types of solutions, the height of the free surface is constant either on lines of constant acceleration or on lines of constant speed. In the case of an isothermal medium, when gamma = 1, we again find simple wave and time-dependent solutions.

An improved free surface modeling for incompressible SPH

In Incompressible Smooth Particle Hydrodynamics (ISPH), a pressure Poisson equation (PPE) is solved to obtain a divergence free velocity field. When free surfaces are simulated using this method a Dirichlet boundary condition for pressure at the free surface has to be applied. In existing ISPH methods this is achieved by identifying free surface particles using heuristically chosen threshold of a parameter such as kernel sum, density or divergence of the position, and explicitly setting their pressure values. This often leads to clumping of particles near the free surface and spraying off of surface particles during splashes. Moreover, surface pressure gradients in flows where surface tension is important are not captured well using this approach. We propose a more accurate semi-analytical approach to impose Dirichlet boundary conditions on the free surface. We show the efficacy of the proposed algorithm by using test cases of elongation of a droplet and dam break. We perform two dimensional simulations of water entry and validate the proposed algorithm with experimental results. Further, a three dimensional simulation of droplet splash is shown to compare well with the Volume-of-Fluid simulations. (C) 2014 Elsevier Ltd. All rights reserved.

Energy identities in water wave theory for free-surface boundary condition with higher-order derivatives

A modified form of Green's integral theorem is employed to derive the energy identity in any water wave diffraction problem in a single-layer fluid for free-surface boundary condition with higher-order derivatives. For a two-layer fluid with free-surface boundary condition involving higher-order derivatives, two forms of energy identities involving transmission and reflection coefficients for any wave diffraction problem are also derived here by the same method. Based on this modified Green's theorem, hydrodynamic relations such as the energy-conservation principle and modified Haskind–Hanaoka relation are derived for radiation and diffraction problems in a single as well as two-layer fluid.

On the impact of a concave nosed axisymmetric body on a free surface

We report on an experimental study of the vertical impact of a concave nosed axisymmetric body on a free surface. Previous studies have shown that bodies with a convex nose, like a sphere, produce a well defined splash with a relatively large cavity behind the model. In contrast, we find that with a concave nose, there is hardly a splash and the cavity extent is greatly reduced. This may be explained by the fact that in the concave nosed case, the initial impact is between a confined air pocket and the free surface unlike in the convex nosed case. From measurements of the unsteady pressure in the concave nose portion, we show that in this case, the maximum pressures are significantly lower than the classically expected ``water hammer'' pressures and also lower than those generally measured on other geometries. Thus, the presence of an air pocket in the case of a concave nosed body adds an interesting dimension to the classical problem of impact of solid bodies on to a free surface. (C) 2015 AIP Publishing LLC.

Effects of mould filling on evolution of the solid-liquid interface during solidification

This work presents a numerical analysis of simultaneous mould filling and phase change for solidification in a two-dimensional rectangular cavity. The role of residual flow strength and temperature gradients within the solidifying domain, caused by the filling process, on the evolution of solidification interface are investigated. An implicit volume of fluid (VOF)-based algorithm has been employed for simulating the free surface flows during the filling process, while the model for solidification is based on a fixed-grid enthalpy-based control volume approach. Solidification modeling is coupled with VOF through User Defined Functions developed in the commercial computational fluid dynamics (CFD) code FLUENT 6.3.26. Comparison between results of the conventional analysis without filling effect and those of the present analysis shows that the residual flow resulting from the filling process significantly influences the progress of the solidification interface. A parametric study is also performed with variables such as cooling rate, filling velocity and filling configuration, in order to investigate the coupled effects of the buoyancy-driven flow and the residual flow on the solidification behavior.

Volume conservation issues in incompressible smoothed particle hydrodynamics

A divergence-free velocity field is usually sought in numerical simulations of incompressible fluids. We show that the particle methods that compute a divergence-free velocity field to achieve incompressibility suffer from a volume conservation issue when a finite time-step position update scheme is used. Further, we propose a deformation gradient based approach to arrive at a velocity field that reduces the volume conservation issues in free surface flows and maintains density uniformity in internal flows while retaining the simplicity of first order time updates. (C) 2015 Elsevier Inc. All rights reserved.

Finite Element Solution Of Surface-Tension Driven Flows In Laser Surface-Melting

Lasers are very efficient in heating localized regions and hence they find a wide application in surface treatment processes. The surface of a material can be selectively modified to give superior wear and corrosion resistance. In laser surface-melting and welding problems, the high temperature gradient prevailing in the free surface induces a surface-tension gradient which is the dominant driving force for convection (known as thermo-capillary or Marangoni convection). It has been reported that the surface-tension driven convection plays a dominant role in determining the melt pool shape. In most of the earlier works on laser-melting and related problems, the finite difference method (FDM) has been used to solve the Navier Stokes equations [1]. Since the Reynolds number is quite high in these cases, upwinding has been used. Though upwinding gives physically realistic solutions even on a coarse grid, the results are inaccurate. McLay and Carey have solved the thermo-capillary flow in welding problems by an implicit finite element method [2]. They used the conventional Galerkin finite element method (FEM) which requires that the pressure be interpolated by one order lower than velocity (mixed interpolation). This restricts the choice of elements to certain higher order elements which need numerical integration for evaluation of element matrices. The implicit algorithm yields a system of nonlinear, unsymmetric equations which are not positive definite. Computations would be possible only with large mainframe computers.Sluzalec [3] has modeled the pulsed laser-melting problem by an explicit method (FEM). He has used the six-node triangular element with mixed interpolation. Since he has considered the buoyancy induced flow only, the velocity values are small. In the present work, an equal order explicit FEM is used to compute the thermo-capillary flow in the laser surface-melting problem. As this method permits equal order interpolation, there is no restriction in the choice of elements. Even linear elements such as the three-node triangular elements can be used. As the governing equations are solved in a sequential manner, the computer memory requirement is less. The finite element formulation is discussed in this paper along with typical numerical results.

Nonlinear surface acoustic waves on an isotropic solid

A systematic derivation of the approximate coupled amplitude equations governing the propagation of a quasi-monochromatic Rayleigh surface wave on an isotropic solid is presented, starting from the non-linear governing differential equations and the non-linear free-surface boundary conditions, using the method of mulitple scales. An explicit solution of these equations for a signalling problem is obtained in terms of hyperbolic functions. In the case of monochromatic excitation, it is shown that the second harmonic amplitude grows initially at the expense of the fundamental and that the amplitudes of the fundamental and second harmonic remain bounded for all time.

Mass transfer with surface reaction in arrays of cylinders at low reynolds and high peclet numbers

Creeping flow hydrodynamics combined with diffusion boundary layer equation are solved in conjunction with free-surface cell model to obtain a solution of the problem of convective transfer with surface reaction for flow parallel to an array of cylindrical pellets at high Peclet numbers and under fast and intermediate kinetics regimes. Expressions are derived for surface concentration, boundary layer thickness, mass flux and Sherwood number in terms of Damkoehler number, Peclet number and void fraction of the array. The theoretical results are evaluated numerically.

Free vibrations of a reservoir-foundation system

The natural frequencies of a reservoir-foundation system are calculated by treating the foundation as a system of linear springs with inertia. The reservoir is treated as consisting of compressible liquid, and the influence of waves at the free surface is included. It is shown that the natural frequencies decrease monotonically as the depth of foundation participating in the motion increases. The influence of waves at the reservoir surface is negligible for the cases normally occurring in practice. It is also shown that the wavelength of motion along the reservoir has no influence on the frequencies when the foundation depth is large compared to the reservoir depth.

The stress in a slowly sheared granular column

We report measurements of the wall stress in a granular material sheared in a cylindrical Couette cell, as a function of the distance from the free surface. Our results shows that when the material is static, all components of the stress saturate to constant values within a short distance from the free surface, in conformity with earlier experiments and theoretical predictions. When the material is sheared by rotating the inner cylinder at a constant rate, the stresses are remarkably altered. The radial normal stress does not saturate, and increases even more rapidly with depth than the linear hydrostatic pressure profile. The axial shear stress changes sign on shearing, and its magnitude increases with depth. These results are discussed in the context of the predictions of the classical and Cosserat plasticity theories.

Transition to turbulence in the free convection boundary layers on an inclined heated plate

An experimental study has been made of transition to turbulence in the free convective flows on a heated plate. Observations have been made with the plate vertical and inclined at angles up to about 50° to the vertical, both above and below the plate. A fibre anemometer was used to survey the region of intermittent turbulence. Information has thus been obtained about the range of Grashof numbers over which transition takes place. Even when the plate is vertical the region of intermittent turbulence is long. When it is inclined, this region becomes still longer in the flow below the plate as a result of the stabilizing stratification, a Richardson number effect. It is possible to have a whole flow such that it should be described as transitional, not laminar or turbulent. It was noticed that in this flow and the vertical plate one, the velocity during the laminar periods could be either of two characteristic values, one of them close to zero. The behaviour above an inclined plate could be interpreted largely as a trend towards the behaviour described in a preceding paper.