Cusped ripples at the plane surface of a viscous liquid

Summary We study the time-evolution of periodical ripples of a viscous liquid at the plane free surface under the action of a distant pure straining flow. We neglect inertial forces (Stokes flow) and include surface tension effects. The solutions for a contracting surface and constant strain rate show that the ripples may develop near-cusps during a stage of the evolution, though later the free surface inevitably asymptotically tends to a smooth plane with vanishing ripples due to the action of capillarity. We obtain the condition for cusp formation in this intermediate stage in terms of the initial capillary number and aspect ratio. If the capillary number is kept constant, the surface tends to shrink through a succession of self-similar trochoidal shapes, whose aspect ratio is given by the capillary number. Received 23 March 1998, accepted for publication 23 July 1998     

Numerical simulation of orbitally shaken viscous fluids with free surface

Orbitally shaken bioreactors are an emerging alternative to stirred‐tank bioreactors for large‐scale mammalian cell culture, but their fluid dynamics is still not well defined. Among the theoretical and practical issues that remain to be resolved, the characterization of the liquid free surface during orbital shaking remains a major challenge because it is an essential aspect of gas transfer and mixing in these reactors. To simulate the fluid behavior and the free surface shape, we developed a numerical method based on the finite element framework. We found that the large density ratio between the liquid and the gas phases induced unphysical results for the free surface shape. We therefore devised a new pressure correction scheme to deal with large density ratios. The simulations operated with this new scheme gave values of wave amplitude similar to the ones measured experimentally. These simulations were used to calculate the shear stress and to study the mixing principle in orbitally shaken bioreactors. Copyright © 2012 John Wiley & Sons, Ltd.     

Interaction of viscous wakes with a free surface

The interaction of laminar wakes with.free-surface waves generated by a moving body beneath the surface of an incompressible viscous fluid of infinite depth was investigated analytically. The analysis was based on the steady Oseen equations for disturbed flows.The kinematic and dynamic boundary conditions were linearized for the small-amplitude free-surface waves. The effect of the moving body was mathematically modeled as an Oseenlet.The disturbed flow was regarded as the sum of an unbounded singular Oseen flow which represents the effect of the viscous wake and a bounded regular Oseen flow which represents the influence of the free surface. The exact solution for the free-surface waves was obtained by the method of integral transforms. The asymptotic representation with additive corrections for the free-surface waves was derived by means of Lighthill‘s two-stage scheme. The symmetric solution obtained shows that the amplitudes of the free-surface waves are exponentially damped by the presences of viscosity and submergence depth.     

Similarity solutions for immiscible phase migration in porous media: an analysis of free boundaries

A fuel pollutant migrating in a water flow throughout a porous medium is distributed between the moving (continuous) and residual (discontinuous) phases. Usually, there is an equilibrium condition between these phases. In this study, the migration of a fuel slug confied within free boundaries moving in the porous medium is considered. This type of fuel migration pertains to circumstances in which convective fuel transport dominates fuel dispersion when fuel saturation approaches zero. A one-dimensional self-similar model is developed, describing the movement of fuel saturation fronts in a porous medium against and with the water flow direction. Several analytical solutions are found revealing the effects of the pore size, fuel viscosity, fuel mass, and the capillary number on the fuel migration in the porous medium.     

Interaction of an obliquely rising vortex ring with a free surface in a viscous fluid

A vortex ring, which approaches a free surface under various initial conditions and different values of flow parameters, is investigated by means of numerical solutions of the Navier-Stokes equations. The study focuses on the connection process of the vortex ring with the free surface at low Reynolds number. Discrepancies between a numerical solution and experimental observations found recently with regard to the shape of the reconnected vorticity tubes have been resolved. The crucial parameter, which determines the essential differences at the low Reynolds number selected, is the angle of inclination at which the vortex ring rises toward the free surface. At an angle of inclination of 20° circular vorticity tubes at the free surface evolve, while at an angle of 45° the reconnection takes place in the form of an almost circular sheet. At higher Reynolds number this form of reconnection may change to the other form through instability. The possibility of vortex-ring reflection at the free surface is discussed.

---

Sommario Usando soluzioni numeriche delle equazioni di Navier-Stokes, viene analizzato il comportamento di un vortice and anello che si avvicina ad una superficie libera, per diverse condizioni iniziali e valori dei parametri del flusso. Lo studio è incentrato sul processo di connessione del vortice adnanello con la superficie libera per bassi numeri di Reynolds. con riguardo alla forma dei tubi di vorticità riconnessi con la superficie, tutte le differenze tra osservazioni sperimentali e simulazioni numeriche sono state risolte. Per bassi numeri di Reynolds, il parametro cruciale dell' interazione è l'angolo di inclinazione tra vortice e superficie. Per un angolo di 20° si formano tubi circolari di vorticità alla superficie, mentre a 45° la riconnessione ha luogo in forma di strato sottile di vorticità. Per altri numeri di Reynolds la forma della riconnessione può cambiare in altre forme a causa dell'instabilità. Infine viene discussa la possibilità di riflessione del vortice ad anello alla superficie.

     

Combined free and forced convection in a porous medium between two vertical walls with viscous dissipation

Viscous dissipation effects in the problem of a fully-developed combined free and forced convection flow between two symmetrically and asymmetrically heated vertical parallel walls filled with a porous medium is analyzed. The equation of motion contains the modified Rayleigh number for a porous medium and the small-order viscous dissipation parameter. Particular attention is given to the solutions near the critical Rayleigh numbers at which infinite flow rates are predicted. Information concerning the multiplicity of solutions at critical Rayleigh numbers is also deduced from perturbation solutions of the governing equation.     

Simulation of three-dimensional unsteady viscous free surface flow around a ship model

We present here a numerical method for solving the free surface flow around a ship at forward speed in calm water. The fluid is assumed to be Newtonian and the Reynolds-averaged Navier-Stokes equations are solved by a finite difference method. Modelization of turbulence is achieved by the algebraic model proposed by Baldwin and Lomax. Fully non-linear free surface conditions are satisfied in the model and a method to avoid the incompatibility between free surface conditions and no-slip conditions at the waterline is proposed. Numerical results obtained for a Wigley hull are compared with experimental results.     

Effect of vertical grid variability on a free surface flow model

A previously developed numerical model that solves the incompressible, non‐hydrostatic, Navier–Stokes equations for free surface flow is analysed on a non‐uniform vertical grid. The equations are vertically transformed to the σ‐coordinate system and solved in a fractional step manner in which the pressure is computed implicitly by correcting the hydrostatic flow field to be divergence free. Numerical consistency, accuracy and efficiency are assessed with analytical methods and numerical experiments for a varying vertical grid discretization. Specific discretizations are proposed that attain greater accuracy and minimize computational effort when compared to a uniform vertical discretization. Published in 2007 by John Wiley & Sons, Ltd.     

Formation of the free surface of a viscous fluid volume inside a rotating horizontal cylinder

Two-dimensional viscous flow with a free surface in a horizontal cylinder rotating at a constant speed is investigated numerically using the boundary element method. It is shown that in the initial stage of rotation of the cylinder four different variants of the behavior of the free surface can be realized in the stage of transition from horizontal to steady-state form.     

The effect of the slip condition on Stokes and Couette flows due to an oscillating wall: exact solutions

Stokes and Couette flows produced by an oscillatory motion of a wall are analyzed under conditions where the no-slip assumption between the wall and the fluid is no longer valid. The motion of the wall is assumed to have a generic sinusoidal behavior. The exact solutions include both steady periodic and transient velocity profiles. It is found that slip conditions between the wall and the fluid produces lower amplitudes of oscillations in the flow near the oscillating wall than when no-slip assumption is utilized. Further, the relative velocity between the fluid layer at the wall and the speed of the wall is found to overshoot at a specific oscillating slip parameter or vibrational Reynolds number at certain times. In addition, it is found that wall slip reduces the transient velocity for Stokes flow while minimum transient effects for Couette flow is achieved only for large and small values of the wall slip coefficient and the gap thickness, respectively. The time needed to reach to steady periodic Stokes flow due to sine oscillations is greater than that for cosine oscillations with both wall slip and no-slip conditions.     

The stability of the flow of a layer of a viscous liquid with moment stresses along an inclined plane

The stability of the flow of a layer of liquid over an inclined plane, taking account of the spin of the molecules and the internal moment stresses, was discussed in [1], However, in [1], a number of errors were allowed to creep in, which led the authors to untrue qualitative and quantitative results. In the present work, the stability of the flow of a layer with respect to long-wave perturbations is investigated by the method of successive approximations [2, 3] under the assumption that the coefficient of rotational viscosity nr is considerably less than the coefficient of Newtonian viscosity . It is shown that, in a first approximation, internal moment stresses do not affect the stability of the motion, and that the rotation of the particles exerts a destabilizing effect on the flow of the layer with respect to three-dimensional periodic perturbations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, pp. 149–151, September–October, 1977.     

Evolution of the free surface of a plane channel during filling with a viscous fluid

The behavior of the free surface of a viscous incompressible fluid is studied in the process of filling of a plane channel inclined to the gravity direction. The problem is numerically solved within the framework of the Navier-Stokes equations under the assumption that the fluid is Newtonian and capillary effects can be neglected. The numerical technique is based on the joint use of the SIMPLE and invariant methods. Two filling regimes are found to exist, namely, that with a steady shape of the free surface without the formation of internal cavities and that with cavity formation in the initial filling stage characterized by spouting behavior of the free boundary.     

Interaction of an internal gravity current with a submerged circular cylinder

The problem of hydrodynamic loads due to the interaction of a gravity current propagating over a bottom channel with a submerged circular cylinder is studied experimentally. It was shown that in the examined range of parameters, the hydrodynamic loads are simulated after Froude. The hydrodynamic loads are maximal if the cylinder lies on the bottom, and they decrease rapidly with increase in the distance from the cylinder to the channel bottom. The effects of mixing and entrainment on the nature of the hydrodynamic loads are considered.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 2, pp. 81–90, March–April, 2005.     

Doppler effects of an oscillating line source in shear flow with a free surface

The linearised water-wave radiation problem for the oscillating 2D submerged source in an inviscid shear flow with a free surface is investigated analytically. There is a nonzero surface velocity. The depth is infinite and the vorticity is uniform. The amplitudes radiated from the source are calculated analytically. Due to Doppler effects, there may be up to four different emitted waves, and there is resonance with zero group velocity and infinite amplitude.     

Influence of the capillarity on a creeping film flow down an inclined plane with an edge

Summary In creeping flows of thin films, the capillarity can play a dominant role. In this paper, the creeping film flow down an inclined plane with an edge is considered. The influence of the capillarity on the velocity and the film surface is studied analytically, numerically and experimentally. Received 12 April 1999; accepted for publication 9 May 1999     

Exact solutions for the incompressible viscous magnetohydrodynamic fluid of a rotating-disk flow with Hall current

The focus of the present study is to obtain exact solutions for the flow of a viscous hydromagnetic fluid due to the rotation of an infinite disk in the presence of an axial uniform steady magnetic field with the inclusion of Hall current effect. In place of the traditional von Karman's axisymmetric evolution of the flow, the rotational non-axisymmetric stationary conducting flow is taken into consideration here, whose governing equations allow an exact solution to develop bounded everywhere in the normal direction to the wall.The three-dimensional equations of motion are treated analytically yielding derivation of exact solutions, which differ from those of corresponding to the classical von Karman's conducting flow. Making use of this solution, analytical formulas for the angular velocity components, for the current density field as well as for the wall shear stresses are extracted. The critical peripheral locations at which extrema of the local skin friction occur are also determined. It is proved from the analytical results that for the specific flow the properly defined thicknesses decay as the magnetic field strength increases in magnitude, approaching their hydrodynamic value in the limit of large Hall numbers.Interaction of the resolved flow field with the surrounding temperature is further analyzed via the energy equation. The temperature field is shown to accord with the dissipation function. According to the Fourier's heat law, a constant heat transfer from the disk to the fluid occurs, though it increases by the presence of magnetic field, the increase is slowed down by the Hall effect eventually reaching its hydrodynamic limit.     

Effect of viscous dissipation on an MHD free convective flow past a semi-infinite vertical cone with a variable surface heat flux

An analysis is presented to investigate the influence of viscous dissipation on a free convection flow over a vertical cone with a variable surface heat flux under the action of a transverse magnetic field. The heat transfer characteristics of the free convection flow are investigated numerically. Numerical solutions for transformed governing equations with a variable surface heat flux are obtained. Velocity, temperature, local shear stress, and heat transfer coefficients are calculated for various values of the problem parameters and presented in the graphical form. The effects of the magnetic parameter, the dissipation number, the power-law index, the angle between the cone generatrix and the vertical line, and the Prandtl number on the flow are discussed. For validation of the present numerical results, they are compared with available experimental data and are found to agree well.     

Thermal criticality for a reactive gravity driven thin film flow of a third-grade fluid with adiabatic free surface down an inclined plane

This study is devoted to the investigation of thermal criticality for a reactive gravity driven thin film flow of a third-grade fluid with adiabatic free surface down an inclined isothermal plane. It is assumed that the reaction is exothermic under Arrhenius kinetics, neglecting the consumption of the material. The governing non-linear equations for conservation of momentum and energy are obtained and solved by using a new computational approach based on a special type of Hermite-Padé approximation technique implemented in MAPLE. This semi-numerical scheme offers some advantages over solutions obtained with traditional methods such as finite differences, spectral method, and shooting method. It reveals the analytical structure of the solution function. Important properties of overall flow structure including velocity field, temperature field, thermal criticality, and bifurcations are discussed.     

A computational method for simulating transient motions of an incompressible inviscid fluid with a free surface

A new numerical method has been developed for the analysis of unsteady free surface flow problems. The problem under consideration is formulated mathematically as a two-dimensional non-linear initial boundary value problem with unknown quantities of a velocity potential and a free surface profile. The basic equations are discretized spacewise with a boundary element method and timewise with a truncated forward-time Taylor series. The key feature of the present paper lies in the method used to compute the time derivatives of the unknown quantities in the Taylor series. The use of the Taylor series expansion has enabled us to employ a variable time-stepping method. The size of time increment is determined at each time step so that the remainders of the truncated Taylor series should be equal to a given small error limit. Such a variable time-stepping technique has made a great contribution to numerically stable computations. A wave-making problem in a two-dimensional rectangular water tank has been analysed. The computational accuracy has been verified by comparing the present numerical results with available experimental data. Good agreement is obtained.