Stability and nonlinear wavy regimes in downward film flows on a corrugated surface

The linear and nonlinear stability of downward viscous film flows on a corrugated surface to freesurface perturbations is analyzed theoretically. The study is performed with the use of an integral approach in ranges of parameters where the calculated results and the corresponding solutions of Navier-Stokes equations (downward wavy flow on a smooth wall and waveless flow along a corrugated surface) are in good agreement. It is demonstrated that, for moderate Reynolds numbers, there is a range of corrugation parameters (amplitude and period) where all linear perturbations of the free surface decay. For high Reynolds numbers, the waveless downward flow is unstable. Various nonlinear wavy regimes induced by varying the corrugation amplitude are determined. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 110–120, January–February, 2007.     

Wavy flow of a liquid film in the presence of a cocurrent turbulent gas flow

Wavy downflow of viscous liquid films in the presence of a cocurrent turbulent gas flow is analyzed theoretically. The parameters of two-dimensional steady-state traveling waves are calculated for wide ranges of liquid Reynolds number and gas flow velocity. The hydrodynamic characteristics of the liquid flow are computed using the full Navier-Stokes equations. The wavy interface is regarded as a small perturbation, and the equations for the gas are linearized in the vicinity of the main turbulent flow. Various optimal film flow regimes are obtained for the calculated nonlinear waves branching from the plane-parallel flow. It is shown that for high velocities of the cocurrent gas flow, the calculated wave characteristics correspond to those of ripple waves observed in experiments.     

Waves on down-flowing fluid films: Calculation of resistance to arbitrary two-dimensional perturbations and optimal downflow conditions

Wavy downflow of viscous fluid films is studied. The full Navier-Stokes equations are used to calculate the hydrodynamic characteristics of the flow. The stability of calculated nonlinear waves to arbitrary two-dimensional perturbations is considered within the framework of the Floquet theory. It is shown that, for small values of the Kapitza number, the waves are stable over a wide range of wavelengths and values of the Reynolds number. It is found that, as the Kapitza number increases, the parameter range where nonlinear waves are calculated is divided into a series of alternating zones of stable and unstable solutions. A large number of narrow zones where the solutions are stable are revealed on the wavelength-Reynolds number parameter plane for large values of the Kapitza number. Optimal regimes of film downflow that correspond to the minimum value of average film thickness for nonlinear waves with different wavelengths are determined. The basic characteristics of these waves are calculated in a wide range of Reynolds and Kapitza numbers.     

Calculation of Linear Stability of a Stratified Gas–Liquid Flow in an Inclined Plane Channel

Linear stability of liquid and gas counterflows in an inclined channel is considered. The full Navier–Stokes equations for both phases are linearized, and the dynamics of periodic disturbances is determined by means of solving a spectral problem in wide ranges of Reynolds numbers for the liquid and vapor velocity. Two unstable modes are found in the examined ranges: surface mode (corresponding to the Kapitsa waves at small velocities of the gas) and shear mode in the gas phase. The wave length and the phase velocity of neutral disturbances of both modes are calculated as functions of the Reynolds number for the liquid. It is shown that these dependences for the surface mode are significantly affected by the gas velocity.     

Effect of a Thin Longitudinal Inviscid Vortex on a Two-Dimensional Pre-Separation Boundary Layer

For large Reynolds numbers, an asymptotic solution of the Navier-Stokes equations describing the effect of a thin longitudinal vortex with a constant circulation on the development of an incompressible steady two-dimensional laminar boundary layer on a flat plate is obtained. It is established that, in a narrow wall region extending along the vortex filament, the viscous flow is described by the 3-D boundary layer equations. A solution of these equations for small values of the vortex circulation is studied. It is found that the solution of the two-dimensional pre-separation boundary layer equations collapses. This is attributable to the singular behavior of the 3-D disturbances near the zero-longitudinal-friction points.     

Exact and asymptotic solutions of the Navier-Stokes equations in a fluid layer between plates approaching and moving apart from one another

Exact solutions of the Navier-Stokes equations are investigated in the layer between parallel plates the distance between which changes proportionally to the square root of time. At the boundaries of the plates the no-slip condition is assigned. For approaching plates a countable family of exact solutions each of which continuously depends on the Reynolds number is obtained. At a sufficiently large Reynolds number, near the boundary a counterflow is formed: the velocity is directed oppositely to the average velocity. On the basis of the exact solution obtained, relative errors are calculated for the asymptotic theories of Reynolds lubricating layer and Prandtl boundary layer.     

Finite-Amplitude Equilibrium Solutions for Plane Poiseuille–Couette Flow

Two-dimensional nonlinear equilibrium solutions for the plane Poiseuille–Couette flow are computed by directly solving the full Navier–Stokes equations as a nonlinear eigenvalue problem. The equations are solved using the two-point fourth-order compact scheme and the Newton–Raphson iteration technique. The linear eigenvalue computations show that the combined Poiseuille–Couette flow is stable at all Reynolds numbers when the Couette velocity component σ₂ exceeds 0.34552. Starting with the neutral solution for the plane Poiseuille flow, the nonlinear neutral surfaces for the combined Poiseuille–Couette flow were mapped out by gradually increasing the velocity component σ₂. It is found that, for small σ₂, the neutral surfaces stay in the same family as that for the plane Poiseuille flow, and the nonlinear critical Reynolds number gradually increases with increasing σ₂. When the Couette velocity component is increased further, the neutral curve deviates from that for the Poiseuille flow with an appearance of a new loop at low wave numbers and at very low energy. By gradually increasing the σ₂ values at a constant Reynolds number, the nonlinear critical Reynolds numbers were determined as a function of σ₂. The results show that the nonlinear neutral curve is similar in shape to a linear case. The critical Reynolds number increases slowly up to σ₂∼ 0.2 and remains constant until σ₂∼ 0.58. Beyond σ₂ > 0.59, the critical Reynolds number increases sharply. From the computed results it is concluded that two-dimensional nonlinear equilibrium solutions do not exist beyond a critical σ₂ value of about 0.59. Received: 26 November 1996 and accepted 12 May 1997     

Stability of a viscous liquid film flowing down a periodic surface

The paper is devoted to a theoretical analysis of linear stability of the viscous liquid film flowing down a wavy surface. The study is based on the Navier–Stokes equations in their full statement. The developed numerical algorithm allows us to obtain pioneer results in the stability of the film flow down a corrugated surface without asymptotic approximations in a wide range over Reynolds and Kapitsa’s numbers. It is shown that in the case of moderate Reynolds numbers there is a region of the corrugation parameters (amplitude and period) where all disturbances decay in time and the wall corrugation demonstrates a stabilizing effect. At the same time, there exist corrugation parameters at which the steady-state solution is unstable with respect to perturbations of the same period as the period of corrugation. In this case the waveless solution cannot be observed in reality and the wall corrugation demonstrates a destabilizing effect.     

NUMERICAL SOLUTION OF THE SLOW TRANSLATION OF A SPHERE MOVING ALONG THE AXIS OF A ROTATING VISCOUS FLUID

The motion of a sphere along the axis of rotation of an incompressible viscous fluid that is rotating as a solid mass is investigated by means of numerical methods for small values of Reynolds numbers and moderate values of Taylor numbers. The Navier-Stokes equations governing the steady, axisymmetric, viscous flow can be written as three coupled, nonlinear, elliptic partial differential equations for the stream function, vorticity and rotational velocity component. Finite difference method is used for solving the governing equations. Second order derivatives are approximated by central differences and nonlinear terms are approximated by upwind differences. Results are presented mostly in the form of graphs of the streamlines and vorticity lines. When 1/ Ro > 2.2, separation occurs and reverse flow is obtained.     

The Onset and Nonlinear Development of Vortex Instabilities in a Horizontal Forced Convection Boundary Layer With Uniform Surface Suction

We consider the flow and heat transfer caused by a strong external flow passing over a hot surface with uniform surface suction. When the Péclet number based on the external velocity is sufficiently large, the resulting thermal boundary layer develops in a nonsimilar manner until it attains an asymptotic state which is independent of the streamwise coordinate, x, when it is dominated by the surface suction. For sufficiently large, but moderate, values of the Darcy–Rayleigh number this boundary layer becomes unstable to streamwise vortex disturbances. We employ a parabolic solver to determine how such disturbances, when placed very close to the leading edge, evolve with distance downstream. Neutral stability is then defined to be when a suitable energy functional ceases to decay/grow as x increases. Thus a neutral curve may be mapped out based upon the behaviour of this functional. Given that the uniform asymptotic state is well known to admit subcritical instabilities, our linearised analysis is extended into the nonlinear domain and the effect of different magnitudes of disturbance is ascertained. It is found that a surprisingly rich variety of vortex pattern emerges which is sometimes sensitively dependent on the values of the governing parameters. These patterns include wavy vortices and abrupt changes in perceived wavelength.     

Eine asymptotische Analyse des Wärmeüberganges im Einlaufbereich von turbulenten Kanal- und Rohrströmungen

The basic equations for turbulent entrance flow are deduced from an asymptotic expansion of the Navier-Stokes equations and the thermal energy equation forRe→∞. Together with a turbulence model they can be solved numerically. Solutions are independent of the Reynolds and Prandtl number. Based on theses solutions, the skin friction and heat transfer as well as velocity and temperature profiles can be determined for finite Reynolds numbers and Prandtl numbersO (1).     

Numerical study of Navier-Stokes equations

The problem of a viscous incompressible fluid flow around a body of revolution at incidence, which is described by Navier-Stokes equations, is considered. For low Reynolds numbers, the solutions of these equations are smooth functions. A numerical algorithm without saturation is constructed, which responds to solution smoothness. The calculations are performed on grids consisting of 900 (10 × 10 × 9) and 700 (10 × 10 × 7) nodes. On the grid consisting of 900 nodes, a system of 3600 nonlinear equations is solved by a standard code. The pressures on the shaded side of the body of revolution are compared. It is found that a numerical study (on this grid) is feasible for problems with Re ≈ 1. For high Reynolds numbers, the number of grid nodes has to be increased. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 5, pp. 43–52, September–October, 2007.     

Asymptotic analysis of a pulsating flow of viscous fluid in a symmetric deformed channel

The time-periodic flow of a viscous incompressible fluid in a two-dimensional symmetric channel with slightly deformed walls is considered. The solution of the Navier-Stokes equations is constructed by means of the method of matched asymptotic expansions [1] at large characteristic Reynolds numbers. It is shown that in an unsteady flow a region of nonlinear perturbations surrounds the line of zero velocity inside the fluid. The formation and development of such nonlinear zones with respect to time is considered. An alternation of the topological features of the streamline pattern in the nonlinear perturbation zone is discovered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 17–23, July–August, 1987.The author is deeply grateful to V. V. Sychev for his formulation of the problem and his attentive attitude to my work.     

Hyperbolic Approximation of the Navier-Stokes Equations for Viscous Mixed Flows

Simplified two-dimensional Navier-Stokes equations of the hyperbolic type are derived for viscous mixed (with transition through the sonic velocity) internal and external flows as a result of a special splitting of the pressure gradient in the predominant flow direction into hyperbolic and elliptic components. The application of these equations is illustrated with reference to the calculation of Laval nozzle flows and the problem of supersonic flow past blunt bodies. The hyperbolic approximation obtained adequately describes the interaction between the stream and surfaces for internal and external flows and can be used over a wide Mach number range at moderate and high Reynolds numbers. Examples of the calculation of viscous mixed flows in a Laval nozzle with large longitudinal throat curvature and in a shock layer in the neighborhood of a sphere and a large-aspect-ratio hemisphere-cylinder are given. The problem of determining the drag coefficient of cold and hot spheres is solved in a new formulation for supersonic air flow over a wide range of Reynolds numbers. In the case of low and moderate Reynolds numbers a drag reduction effect is detected when the surface of the sphere is cooled.     

Perturbation solutions for asymmetric laminar flow in porous channel with expanding and contracting walls

The cases of large Reynolds number and small expansion ratio for the asym- metric laminar flow through a two-dimensional porous expanding channel are considered. The Navier-Stokes equations are reduced to a nonlinear fourth-order ordinary differential equation by introducing a time and space similar transformation. A singular perturbation method is used for the large suction Reynolds case to obtain an asymptotic solution by matching outer and inner solutions. For the case of small expansion ratios, we are able to obtain asymptotic solutions by double parameter expansion in either a small Reynolds number or a small asymmetric parameter. The asymptotic solutions indicate that the Reynolds number and expansion ratio play an important role in the flow behavior. Nu- merical methods are also designed to confirm the correctness of the present asymptotic solutions.     

Effect of Inertial Terms on Fluid–Porous Medium Flow Coupling

The study considers an effect of the nonlinear inertial terms in the Brinkman filtration equation on the characteristics of coupled flows in a pure fluid and porous medium in the frameworks of two independent problems. The first problem is the forced boundary-layer flow overlying the Darcy–Brinkman porous medium. The Prandtl theory is used, and the self-similar equations are built to describe it. It is shown that the inertial terms have a valuable effect on the boundary-layer structure because of the large velocity gradient in the transition zone. The boundary-layer thickness in a porous medium rapidly grows at large Reynolds numbers. The velocity magnitude and gradient at the interface also change. The second independent problem is an analysis of the inertial terms effect on the flow stability. The neutral curves of the full and linearized flow models are built using the shooting method. They have different short-wave asymptotic, but there are no significant changes in the critical Reynolds numbers and corresponding wave numbers.     

On the flow of a viscous fluid driven along a channel by suction at porous walls

This paper is a theoretical treatment of the flow of a viscous incompressible fluid driven along a channel by steady uniform suction through porous parallel rigid walls. Many authors have found such flows when they are symmetric, steady and two-dimensional, by assuming a similarity form of solution due to Berman in order to reduce the Navier-Stokes equations to a nonlinear ordinary differential equation. We generalise their work by considering asymmetric flows, unsteady flows and three-dimensional perturbations. By use of numerical calculations, matched asymptotic expansions for large values of the Reynolds number, and the theory of dynamical systems, we find many more exact solutions of the Navier-Stokes equations, examine their stability, and interpret them. In particular, we show that most previously found steady solutions are unstable to antisymmetric two-dimensional disturbances. This leads to a pitchfork bifurcation, stable asymmetric steady solutions, a Hopf bifurcation, stable time-periodic solutions, stable quasi-periodic solutions, phase locking and chaos in succession as the Reynolds number increases.     

Wave formation on vertical falling liquid films

The method of integral relations is used to derive a nonlinear “two-wave” structure equation for long waves on the surface of vertical falling liquid films. This equation is valid in a wide range of Reynolds numbers and reduces to the known equations for high and low Re. Theoretical data for the fastest growing waves are compared with the experimental results on velocities, wave numbers and growth rates of the waves in the inception region. The validity of theoretical assumptions is also confirmed by the direct measurements of the instantaneous velocity profiles in a wave liquid film.     

Asymptotic solutions of the two-dimensional equations for a thin viscous shock layer in a rarefied gas

Two-dimensional hypersonic rarefied gas flow around blunt bodies is investigated for the continuum to free-molecular transition regime. In [1], as a result of an asymptotic analysis, three rarefied gas flow regimes, depending on the relationship between the problem parameters, were detected and one of these regimes was investigated. In the present study, asymptotic solutions of the thin viscous shock layer equations at small Reynolds numbers are obtained for the other two flow regimes. Analytical expressions for the heat transfer, friction and pressure coefficients are obtained as functions of the incident flow parameters and the body geometry and temperature. As the Reynolds number tends to zero, the values of these coefficients approach their values in free-molecular flow. The scaling parameters of hypersonic rarefied gas flow around bodies are determined for different regimes. The asymptotic solutions are compared with the results of direct Monte Carlo simulation.