Next Article in Journal
State-Space Compression for Efficient Policy Learning in Crude Oil Scheduling
Next Article in Special Issue
Predictability of Magnetic Field Reversals
Previous Article in Journal
On Global Solutions of Hyperbolic Equations with Positive Coefficients at Nonlocal Potentials
Previous Article in Special Issue
General Solutions for MHD Motions of Ordinary and Fractional Maxwell Fluids through Porous Medium When Differential Expressions of Shear Stress Are Prescribed on Boundary
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Velocity Field due to a Vertical Deformation of the Bottom of a Laminar Free-Surface Fluid Flow

by
Rodrigo González
1 and
Aldo Tamburrino
2,*
1
Doctoral Program in Engineering Sciences, Mention Fluid Dynamics, Faculty of Physical and Mathematical Sciences, University of Chile, Santiago 8370448, Chile
2
Department of Civil Engineering and Advanced Mining Technology Center, Faculty of Physical and Mathematical Sciences, University of Chile, Santiago 8370448, Chile
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(3), 394; https://doi.org/10.3390/math12030394
Submission received: 30 November 2023 / Revised: 18 January 2024 / Accepted: 22 January 2024 / Published: 25 January 2024
(This article belongs to the Special Issue Applications of Mathematics to Fluid Dynamics)

Abstract

:
This article investigates the velocity field of a free-surface flow subjected to harmonic deformation of the channel bottom, progressing asymptotically from a flat initial state to a maximum amplitude. Assuming a uniform main flow with the primary velocity component transverse to the bed undulation, analytical solutions are obtained for the three velocity components and free surface distortion using the method of perturbations. The perturbation components of the velocity field, streamlines, and surface deformation depend on a dimensionless parameter that reflects the fluid inertia induced by bed deformation relative to viscous resistance. When viscous effects dominate, a monotonic decay of the perturbations from the bed to the free surface is observed. In contrast, when inertia dominates, the perturbations can exhibit an oscillatory behavior and introduce circulation cells in the plane normal to the main flow. The interplay between inertia and viscosity reveals scenarios where surface and bed deformations are either in or out of phase, influencing vertical velocity components. Figures illustrate these phenomena, providing insights into the complex dynamics of free-surface flows with harmonic bed deformation in the direction normal to the main flow, and amplitude growing with time. The results are limited to small deformations of the channel bottom, as imposed by the linearization of the momentum equations. Even so, to the best of the authors’ knowledge, this problem has not been addressed before.

1. Introduction and Objective

The Navier–Stokes equations constitute the foundations of fluid mechanics. Published by Claude Navier in two articles in 1821 and 1823, they were definitively consolidated with the work of G.G. Stokes, who settled the no-slip condition at the fluid–solid boundary in 1845 [1]. Due to the nonlinearity of the equations, analytical solutions are limited to relatively few cases, typically those in which the nonlinear terms are absent, a condition that arises in specific flow scenarios, for instance, when parallel streamlines are present. This is the case of some flows solved in the seminal articles by Navier and Stokes. Thus, Navier derived the velocity distribution for a steady uniform free-surface flow in a rectangular channel [2]. Stokes, in his study on the motion of pendulums, solved “the extremely simple case of an oscillating plane” [3] (p. 20). In a note to the same study, Stokes also provided the solution of the velocity distribution due to the sudden motion of a plate initially at rest [3] (Note B, p. 100). Problems where the solid boundary is not a plane, but a deformed surface are much more challenging to solve analytically. In this case, analytical solutions are available only when the nonlinear terms can be addressed through linearization. Focusing on the velocity distribution in free-surface flows, perhaps the classic problem corresponds to a wavy bed with undulations transverse to the flow, as solved by Wang [4], who addressed the problem for a low Reynolds number film flow over a wavy bed with amplitude much higher than the wavelength. The propagation of water waves over rigid rippled beds was studied numerically by Huang and Dong [5].
In the case of unsteady fluid motion, aside from instances where the equations can be linearized, there is limited knowledge about the structure of solutions as the physical and geometric parameters of the problem undergo variations [6]. Wierschem et al. [7,8] determined the velocity distribution of film flows down an inclined wavy channel, considering the effect of surface tension. The flow over a moving wavy bottom was approached by Benjamin [9]. Another flow configuration for which an analytical result was obtained is that studied by Lyne [10], who, in addition to the velocity field, determined the steady stream generated by an oscillatory viscous flow over a wavy wall. Lyne’s solution, obtained by the method of perturbations was extended by Kaneko and Honji [11], and subsequently by Vittori [12].
The free surface behavior, including the stability, of the flow over a wavy bed has been studied by Pozrikidis [13], Bontozoglou and Papapolymerou [14], Scholle et al. [15], Trifonov [16,17,18], and Tseluiko et al. [19], among many others. Pozrikidis [13] and Scholle et al. [15] made the analysis considering the lubrication approximation. Trifonov [18] also computed the velocity field numerically. Dávalos-Orozco [20] studied the nonlinear instability of a thin film flowing down a smoothly deformed surface. Although his analysis is general as the surface can vary as much as in the direction of the main flow as in the direction normal to it, for simplicity, the numerical analysis is carried out considering a sinusoidal wall that deforms in the flow direction.
The oscillating flow over a wavy bottom problem has been approached by many authors due to its implication in sediment transport and ripple formation at the bottom of sea waves. In addition to the analytical approaches developed by Lyne [10], Kaneko and Honji [11], and Vittori [12], the works by Hara and Mei [21], Nishimura et al. [22], and Shen and Chan [23] can be mentioned. In all the above studies, the generatrix of the bed undulations is transverse to the direction of the main flow, as sketched in Figure 1.
The configuration in which the main flow is parallel to the bed deformation has also been studied, although all the literature found by the authors refers to the case in which the fluid is confined between two plates, and the flow is driven by a pressure gradient. Thus, Szumbarski [24] computed semi-analytically the velocity distribution and performed a linear stability analysis of the flow. Results of linear stability analysis have also been presented by Kowalewski et al. [25], Szumbarski and Błoński [26], and Moradi and Floryan [27]. Szumbarski et al. [28] studied numerically the effect of the wall corrugation on the hydraulic resistance of the flow. Among the results presented in their article, they show the circulation pattern induced by the velocity’s normal to the main flow direction, generating closed circulation cells. Gepner et al. [29] focused on the nonlinear solutions of the velocity field generated by the flow limited by two corrugated plates and the effect of the induced secondary currents enhancing diffusive transport. Yadav and Gepner [30] studied the effect on the flow instabilities for a geometry consisting in a fixed corrugated bottom plate and an upper moving plane plate. Later, Gepner and Floryan [31] approached numerically the flow driven by a pressure gradient of a fluid bounded by two sinusoidal walls which are in phase in flows at low and moderate sub-turbulent Reynolds numbers. As the fluid moves transverse to the plate deformations, it is prone to separation and generation of a recirculating region in the valleys of the walls. Previously, Smith [32] had studied the viscous Poiseuille 2-D flow, the walls of which had symmetrical and non-symmetrical cavities. His approach was theoretical, and the final equations were solved numerically. Bordner [33] had solved numerically the laminar boundary layer over a periodic stationary wavy surface, noticing the separation effect near the valleys of the wavy surface. A similar problem was solved by Sobey [6], who modeled and computed numerically the flow in an extracorporeal membrane oxygenator for both the steady and oscillatory flow conditions. Caponi et al. [34], also numerically, found that separation is obtained from the linearized Navier–Stokes equations for the laminar flow near a moving wavy surface. They also determined the range of validity of the linear approximation of the Navier–Stokes equations. Tsangaris and Leiter [35], using a perturbation method, and solving numerically some resulting equations, studied the separation of viscous laminar flows through wavy-walled channels. The separation characteristics of fluid flow inside two parallel wavy plates were also studied numerically by Mahmud et al. [36], who found that separation occurs for a critical Reynolds number which depends on dimensionless values of the amplitude and wavelength of the bed undulation.
Figure 1. Configuration analyzed in previous studies with free-surface or unconfined flows: (a) free-surface fluid flow over a fixed bottom [1]; (b) flow induced by the sudden motion of a plane plate [3]; (c) flow induced by the oscillatory motion of a plane plate [3]; (d) flow over a fixed wavy bed [4,5,9,13,14,15,16,17,18,19,20]; (e) flow resulting from the motion of a wavy bottom [9,34]; and (f) oscillatory flow over a wavy bottom [10,11,12,21,22,23].
Figure 1. Configuration analyzed in previous studies with free-surface or unconfined flows: (a) free-surface fluid flow over a fixed bottom [1]; (b) flow induced by the sudden motion of a plane plate [3]; (c) flow induced by the oscillatory motion of a plane plate [3]; (d) flow over a fixed wavy bed [4,5,9,13,14,15,16,17,18,19,20]; (e) flow resulting from the motion of a wavy bottom [9,34]; and (f) oscillatory flow over a wavy bottom [10,11,12,21,22,23].
Mathematics 12 00394 g001
So far, the cited references on flows over deformed walls correspond to pre-established geometries that do not change over time. However, there are studies in which the upper and lower walls confining the fluid fluctuate temporarily. Thus, Carlsson et al. [37,38] investigated a scenario involving the confinement of fluid within flexible walls, exhibiting transverse motion in the form of small amplitude standing waves. The velocity field was ascertained through a perturbation approach, revealing that streaming occurs at the second order, leading to the formation of cellular flow patterns in the channel. Lebbal et al. [39] analyzed the linear instabilities of the flow between flexible walls, considering that the base condition of the flow and the displacement of the walls are sinusoidal in both time and the direction of the main flow. The resulting equations were solved numerically. Vladimirov [40] considered that the flow domain is the half-space, with a boundary subject to tangential vibrations. The solution of the problem contains both oscillating and non-oscillating components and the streaming flows (either steady or unsteady) are obtained.
The goal of this article is to determine the velocity field of a free-surface flow when the bottom undergoes harmonic deformation in the direction transverse to the main flow (i.e., the undulations’ generatrix is parallel to the main flow), asymptotically approaching a maximum amplitude from a flat initial condition. The problem addressed here arose from the analysis of the generation of rills in their incipient stage due to overland flow over erodible soil [41]. The original problem is highly complex due to the interaction between the flow and the soil. As the soil erodes, it alters the geometry, subsequently changing the flow characteristics, and so on, until eventually reaching a final equilibrium condition. The primary objective of the article is not to solve the physical problem, as the equations accounting for soil erosion are not considered. The problem presented is merely an idealization of the bed deformation process, or a mimicry of it. Nevertheless, as far as the authors are aware, this has not been accomplished previously. Thus, the aim is to determine the flow characteristics of the rill formation process considering the depth of erosion to be small compared to the width of the rill. The main flow is assumed to be uniform (i.e., the flow depth does not change along the main flow direction), with the primary velocity component oriented in the direction transverse to the bed undulation (the main flow is parallel to the undulation generatrix), as illustrated in Figure 2. The specific goal of this paper is to present analytical solutions for the three velocity components and the distortion of the free surface caused by the deformation of the channel bottom. The bottom starts as initially flat and evolves to reach an equilibrium wavy shape.

2. Problem Formulation

The problem to be studied consists of determining the velocity field of the laminar flow with a free surface as it develops over a bottom inclined at an angle θ with respect to the horizontal in the main flow direction. The bottom features a small-amplitude undulation in the direction transverse to the flow that increases over time. A schematic of the problem and the coordinate system is presented in Figure 3. The flow extends indefinitely in the x and y directions. The origin of the z axis is in the median plane that defines the bottom and is perpendicular to it. The x axis is oriented in the direction of the main flow, forming an angle θ with the horizontal. The condition of small amplitude undulation is given by a / h 1 , where a is the amplitude and h the flow depth. Denoting by λ the length of the undulation, the condition a / λ 1 is also satisfied.

2.1. Governing Equations

The flow is governed by the continuity and Navier–Stokes equations [42], given below, where u , v , w are the velocity components along x , y , z , t is the time, p is the fluid pressure, g is the acceleration due to gravity, and ρ and μ are the density and the viscosity of the fluid, respectively.
u x + v y + w z = 0
ρ u t + u u x + v u y + w u z = p x + μ 2 u x 2 + 2 u y 2 + 2 u z 2 + ρ g sin θ
ρ v t + u v x + v v y + w v z = p y + μ 2 v x 2 + 2 v y 2 + 2 v z 2
ρ w t + u w x + v w y + w w z = p z + μ 2 w x 2 + 2 w y 2 + 2 w z 2 ρ g cos θ

2.2. Boundary Conditions

On the bottom, the kinematic boundary condition, corresponding to no-slip, is applied to the components u and v , while the bottom deformation rate is applied to w :
z = 0 ,     x , y , t   :   u = v = 0
  η x , y , t   :   w = d η d t
where η is the bottom level at x , y for the time t , and it is a prescribed function.
In absence of capillarity ( a λ ) and wind effects, the dynamic boundary condition to be imposed at the free surface corresponds to a shear-free one:
z = h ,     x , y , t   :   τ x z = τ y z = 0
To determine the free-surface deformation, the kinematic boundary condition is used:
z = h ,   x , y , t   :   w = d h d t
Thus, the system of differential Equations (1)–(4) with boundary conditions (5)–(8) allows us, formally, to solve the problem.

3. Solution

3.1. Transverse Velocity Field

The problem is solved by the perturbation method around a base flow [43,44]. The undisturbed condition (base flow) is denoted by the sub-index 0, and the perturbations are denoted by the sub-index 1. The base flow is defined by the velocity field u 0 ,   0 ,   0 , such that the velocity vector is given by
u , v , w = u 0 + ε u 1 , ε v 1 ,   ε w 1
where ε = a / h 0 is a small parameter and the velocity components u 1 , v 1 ,   w 1 , properly scaled, are of order 1. In the same way, the flow depth h is given by
h = h 0 + δ h 1
with δ 1 being a small parameter to be determined. The location η of the bottom is chosen to be
η = η R x   sin θ + a f η e i k y
where η R is an arbitrary reference level, a is the perturbation amplitude, f η is a function of time that describes the bottom deformation, k = 2 π / λ is the wavenumber of the lateral bed deformation, and i = 1 . To simplify the algebra, Euler’s formula is used, understanding that the real part of the results must be taken.
Substituting the velocities given by Equation (9) in the system of Equations (1)–(4), the solution for the base flow is obtained:
u 0 , v 0 ,   w 0 = ρ g sin θ h 0 2 μ z h 0 1 2 z 2 h 0 2   ,   0 ,   0
p = ρ g cos θ h 0 z
Furthermore, by substituting the velocities provided by Equation (9) into the continuity and momentum equations (Equations (1)–(4)) and subsequently subtracting the equations pertaining to the base flow, we derive the equations governing the disturbances:
v 1 y + w 1 z = 0
u 1 t + w 1 u 0 z = μ ρ 2 u 1 y 2 + 2 u 1 z 2
v 1 t = μ ρ 2 v 1 y 2 + 2 v 1 z 2
w 1 t = μ ρ 2 w 1 y 2 + 2 w 1 z 2
Note that, as the base flow is uniform and the bed is not perturbed along the x direction, the continuity and momentum equations are independent of x . The velocity perturbations v 1 and w 1 satisfy the continuity equation in the plane y , z , so it is feasible to look for the solution in terms of a stream function Ψ 1 :
v 1 = Ψ 1 z , w 1 = Ψ 1 y
Expressing Equations (16) or (17) in terms of the stream function, the solution for Ψ 1 is obtained. Using the momentum equation along y :
t ψ 1 z = μ ρ 2 y 2 ψ 1 z + 2 z 2 ψ 1 z
The following structure for the solution of ψ 1 is proposed:
ψ 1 = f   Φ e i k y
where f and Φ are functions of t and z, respectively, which must be determined. Substituting Equation (20) in Equation (19), the following ordinary differential equation for f is obtained:
ρ μ e i k y f t Φ z = f k 2 e i k y Φ z + f e i k y 3 Φ z 3
The above differential equation is of separable variables. Rearranging terms:
ρ μ 1 f f t = k 2 Φ z + 3 Φ z 3 Φ z = ζ
Thus, the temporal variation in the stream function is given by:
f = K 1 e μ ρ ζ t
The dependence of ψ on z is defined from the following differential equation:
3 Φ z 3 ζ + k 2 Φ z = 0
The ζ term of the above equations is not known and will be determined in the course of the solution. Although it can be a complex quantity, the final solution reveals that it is sufficient to perform the analysis assuming it is a real number. We will consider first that ζ is positive. Thus, the discriminant of ω = ζ + k 2 , is also positive and the solution of the above differential equation is
Φ = A sinh ( ω z ) + B cosh ( ω z ) + C
with A ,   B , and C as integration constants. Calling A ¯ = K 1 A ,   B ¯ = K 1 B , and C ¯ = K 1 C , after replacing Equations (23) and (25) in Equation (20), the stream function becomes:
ψ 1 = e μ ρ ζ t A ¯ sinh ( ω z ) + B ¯ cosh ( ω z ) + C ¯ e i k y
Thus, the velocity components v 1 and w 1 are:
v 1 = ω e μ ρ ζ t A ¯ cosh ( ω z ) + B ¯ sinh ( ω z ) e i k y
w 1 = i k e μ ρ ζ t A ¯ sinh ( ω z ) + B ¯ cosh ( ω z ) + C ¯ e i k y
The constants A ¯ ,   B ¯ , and C ¯ are determined from the boundary conditions, which correspond to two kinematic conditions at the bottom and one dynamic condition at the free surface. For v 1 , the condition of no slip on the bottom is imposed. Since the deformation of the bed has been assumed to be small, the condition v 1 = 0 can be applied at z = 0 , from which Equation (27) yields A ¯ = 0 . Thus:
v 1 = ω e μ ρ ζ t B ¯ sinh ( ω z ) e i k y
w 1 = i k e μ ρ ζ t B ¯ cosh ( ω z ) + C ¯ e i k y
Applying the same argument of small deformation as for the bottom, the dynamic condition at the free-surface is applied at z = h 0 , where the shear stress is zero. Replacing Equation (9) in the expression for the shear stress, it can be written as τ y z = τ y z 0 + ε τ y z 1 . With v 1 and w 1 from Equations (28) and (29), the shear stress associated with the perturbation is:
τ y z 1 = μ e μ ρ ζ t ω 2 + k 2 B ¯ cosh ( ω z ) + k 2 C ¯ e i k y
The boundary condition τ y z 1 = 0 at z = h 0 gives
B ¯ = k 2 C ¯ ω 2 + k 2 cosh ( ω h 0 )
In this way, the stream function and transverse velocity components associated to the perturbation are given by:
ψ 1 = e μ ρ ζ t C ¯ ω 2 + k 2 cosh ( ω h 0 ) k 2 cosh ( ω z ) ω 2 + k 2 cosh ( ω h 0 ) e i k y
v 1 = ω e μ ρ ζ t C ¯ k 2 sinh ( ω z ) ω 2 + k 2 cosh ( ω h 0 ) e i k y
w 1 = i k e μ ρ ζ t C ¯ ω 2 + k 2 cosh ( ω h 0 ) k 2 cosh ( ω z ) ω 2 + k 2 cosh ( ω h 0 ) e i k y
The second kinematic condition is enforced at the bottom ( z = 0 ). Its consideration, coupled with the temporal bottom deformation, facilitates the determination of the constants ζ and C ¯ . In effect,
w 1 z = 0 = d η 1 d t = η 1 t
Equation (11) can be written as η = η 0 + ε η 1 , with η 0 = η R x sin θ and η 1 given by:
η 1 = h 0 f η e i k y
where f η is a function of time. Given the structure imposed for the variation in time of the stream function ψ 1   ~   e μ ρ ζ t , the deformation of the bed cannot be arbitrary. In particular, some of the forms that the function f η in Equation (37) can take are an exponential dependence on time, or a constant plus an exponential dependence. For this article, the following function is chosen, which characterizes bed deformation that increases from zero until reaching asymptotically a maximum equilibrium amplitude equal to a .
f η = 1 e β t , β < 0
Thus,
η 1 t = h 0 β e β t e i k y
and the second kinematic condition is expressed as:
i k e μ ρ ζ t C ¯ ω 2 + k 2 cosh ( ω h 0 ) k 2 ω 2 + k 2 cosh ( ω h 0 ) = h 0 β e β t
To ensure the left side of Equation (40) is identical to the right side at any time, ζ and C ¯ must be:
ζ = β ρ μ
C ¯ = i h 0 β k ω 2 + k 2 cosh ω h 0 ω 2 + k 2 cosh ω h 0 k 2
Thus, the stream function associated to v 1 and w 1 is:
ψ 1 = i h 0 β k e β t ω 2 + k 2 cosh ( ω h 0 ) k 2 cosh ( ω z ) ω 2 + k 2 cosh ω h 0 k 2 e i k y
By construction, the stream function is composed of a base function ( ψ 0 ) and a disturbance, similar to velocities, shear stress, and bed deformation, i.e., ψ = ψ 0 + ε ψ 1 . As the base transverse velocities are zero, the stream function results solely from the perturbations, so ψ = a / h 0 ψ 1 . Thus, the stream function is given by:
ψ = i a β k e β t ω 2 + k 2 cosh ( ω h 0 ) k 2 cosh ( ω z ) ω 2 + k 2 cosh ω h 0 k 2 e i k y
From Equation (9), we have v = ε v 1 and w = ε w 1 . Introducing the expressions for the constants in Equations (34) and (35), the transverse velocities are obtained:
v = i a β e β t k ω sinh ( ω z ) ω 2 + k 2 cosh ω h 0 k 2 e i k y
w = a β e β t ω 2 + k 2 cosh ( ω h 0 ) k 2 cosh ( ω z ) ω 2 + k 2 cosh ω h 0 k 2 e i k y
Taking the real part of Equations (44)–(46), the stream function and transverse velocities result to be:
ψ = a β k e β t ω 2 + k 2 cosh ( ω h 0 ) k 2 cosh ( ω z ) ω 2 + k 2 cosh ω h 0 k 2 sin k y
v = a β e β t k ω sinh ( ω z ) ω 2 + k 2 cosh ( ω h 0 ) sin k y
w = a β e β t ω 2 + k 2 cosh ( ω h 0 ) k 2 cosh ( ω z ) ω 2 + k 2 cosh ( ω h 0 ) cos k y

3.2. Deformation of the Free Surface

The deformation of the free surface is obtained from the kinematic condition at z = h 0 :
w z = h 0 = d h d t = h t
Using Equation (46), a differential equation for h results:
h t = a β e β t ω 2 cosh ( ω h 0 ) ω 2 + k 2 cosh ω h 0 k 2 e i k y
Performing the integration with the initial condition t = 0 , h = h 0 , and taking the real part of the solution, the deformed free surface height is obtained.
h = h 0 + a 1 e β t ω 2 cosh ( ω h 0 ) ω 2 + k 2 cosh ω h 0 k 2 e i k y
Identifying terms in Equations (10) and (49), and taking the real part of η 1 (Equation (40)), the value of the amplitude of the free surface deformation δ is found:
δ = ε ω 2 cosh ( ω h 0 ) ω 2 + k 2 cosh ω h 0 k 2
Taking into account that ω h 0 is always a finite quantity, it is easy to see that
0 < ω 2 cosh ( ω h 0 ) ω 2 + k 2 cosh ω h 0 k 2 < 1
from which it results that the deformation of the free surface is always less than that of bottom. Moreover, based on the positive discriminant assumption of ω, it can be concluded from Equations (37) and (52) that both deformations are consistently in phase.

3.3. Longitudinal Velocity Distribution

Knowing the vertical perturbation of the component of the velocity along z (Equation (46)), and the solution at order 0 of u (Equation (12)), the perturbation of u can be obtained by integrating Equation (15), which becomes:
u 1 t ρ g h 0 2 sin θ μ 1 z h 0 β e β t e i k y ω 2 + k 2 cosh ω h 0 k 2 cosh ω z ω 2 + k 2 cosh ω h 0 k 2 = μ ρ 2 u 1 y 2 + 2 u 1 z 2
Consider that the solution u 1 has the following structure:
u 1 = e β t e i k y ρ g h 0 2 sin θ μ U
where U is a function of z . Replacing Equation (56) into Equation (55) and using the definitions of ω and ζ , the equation for u 1 is reduced to:
d 2 U d z 2 ω 2 U + ω 2 k 2 1 z h 0 ω 2 + k 2 cosh ω h 0 k 2 cosh ω z ω 2 + k 2 cosh ω h 0 k 2 = 0
The solution of Equation (57) is rather cumbersome. The following expression for U results after imposing the boundary conditions of zero x -velocity at z = 0 and no shear stress at z = h 0 :
U = A 1 ω 2 1 z h 0 cosh ω z + B 1 4 ω 2 z h 0 cosh ω z + + B 1 h 0 2 ω z h 0 1 2 z 2 h 0 2 sinh ω z + [ A 1 ω 2 cosh ω h 0 1 ω h 0 + sinh ω h 0 B 1 4 ω 2 1 ω h 0 + sinh ω h 0 cosh ω h 0 + ω h 0 ] sinh ω z
Taking the real part of Equation (58) and adding the perturbation to the zeroth-order solution, the velocity along the x-axis is obtained ( u = u 0 + ε u 1 ):
u = ρ g sin θ h 0 2 μ { z h 0 z 2 2 h 0 2 + a h 0 e β t cos ( k y ) ω 2 k 2 ω 2 { ω 2 + k 2 cosh ω h 0 ω 2 + k 2 cosh ω h 0 k 2 × 1 z h 0 cosh ( ω z ) + 1 4 k 2 ω 2 + k 2 cosh ω h 0 k 2 z h 0 cosh ( ω z ) + k 2 ω 2 + k 2 cosh ω h 0 k 2 h 0 ω 2 z h 0 z 2 2 h 0 2 sinh ( ω z ) + [ 1 cosh ( ω h 0 ) ω 2 + k 2 cosh ω h 0 ω 2 + k 2 cosh ω h 0 k 2 1 h 0 ω + sinh ω h 0 1 4 k 2 ω 2 + k 2 cosh ω h 0 k 2 × 1 h 0 ω + sinh ω h 0 cosh ( ω h 0 ) + h 0 ω ] sinh ( ω z ) } }
Thus, the stream function (Equation (47)), the velocity field (Equations (48), (49) and (59)), and the free surface deformation (real part of Equation (52)) were determined for the case ζ + k 2 > 0 .
When ζ + k 2 < 0 , and ω is an imaginary number, there is a slight change in the solutions. From the definition of ω , another variable ϖ is introduced such that ϖ 2 =   ζ + k 2 > 0 . Thus ω = i ϖ , and the expressions for the stream function, velocities, and free surface are:
ψ = a β k e β t ϖ 2 k 2 cos ( ϖ h 0 ) + k 2 cos ( ϖ z ) ϖ 2 k 2 cos ϖ h 0 + k 2 sin k y
v = a β e β t k ϖ sin ( ϖ z ) ϖ 2 k 2 cos ϖ h 0 + k 2 sin k y
w = a β e β t ϖ 2 k 2 cos ( ϖ h 0 ) + k 2 cos ( ϖ z ) ϖ 2 k 2 cos ϖ h 0 + k 2 cos k y
h = h 0 + a 1 e β t ϖ 2 cos ( ϖ h 0 ) ϖ 2 k 2 cos ϖ h 0 + k 2 cos k y
u = ρ g sin θ h 0 2 μ { z h 0 z 2 2 h 0 2 + a h 0 e β t cos ( k y ) ϖ 2 + k 2 ϖ 2 { ϖ 2 k 2 cos ϖ h 0 ϖ 2 k 2 cos ϖ h 0 + k 2 × 1 z h 0 cos ( ϖ z ) 1 4 k 2 ϖ 2 k 2 cos ϖ h 0 + k 2 z h 0 cos ( ϖ z ) + k 2 ϖ 2 k 2 cos ϖ h 0 + k 2 h 0 ϖ 2 z h 0 z 2 2 h 0 2 sin ( ϖ z ) + [ ϖ 2 k 2 ϖ 2 k 2 cos ϖ h 0 + k 2 1 h 0 ϖ sin ϖ h 0 1 4 k 2 ϖ 2 k 2 cos ϖ h 0 + k 2 × 1 h 0 ϖ + sin ϖ h 0 cos ( ϖ h 0 ) + h 0 ϖ ] sin ( ϖ z ) } }
The ratio between the amplitude of the free surface and the amplitude of the deformation of the bottom for ω 2 < 0 is given by
δ ε = ϖ 2 cos ( ϖ h 0 ) ϖ 2 k 2 cos ϖ h 0 + k 2

4. Analysis and Visualization of the Solution

4.1. Dimensionless Parameters

To generalize the results, h 0 , λ , and β are used as scales to make dimensionless z ,   y , and t , respectively:
z = z h 0 , y = y λ , t = β t
The minus sign in t appears because β < 0 (Equation (38)). The dimensionless bed perturbation (Equations (37) and (38)) is written as η 1 = 1 e t e i y . k ^ = 1 / λ is also defined, such that k ^ = λ k ^ = 1 . Thus, the following dimensionless parameters arise:
B = β ρ λ 2 μ , Ω = ω λ = 1 B ,   Λ = h 0 λ
B is always positive and it appears naturally from the scaling of the equations for the perturbed components of the transverse velocities v 1 or w 1 (Equations (16) or (17)). Let V be a scale for the velocity v 1 such that the dimensionless velocity is given by v ^ 1 = v 1 / V . In terms of the dimensionless variables, Equation (16) becomes:
β V v ^ 1 t = μ ρ V λ 2 2 v ^ 1 y 2 + V h 0 2 2 v ^ 1 z 2
β v ^ 1 t = μ ρ λ 2 2 v ^ 1 y 2 + λ h 0 2 2 v ^ 1 z 2
B v ^ 1 t = 2 v ^ 1 y 2 + 1 Λ 2 2 v ^ 1 z 2
Calling W the velocity scale to make dimensionless the perturbation of the velocity in the z direction, w ^ 1 = w 1 / W . Thus, the dimensionless version of Equation (17) is:
B w ^ 1 t = 2 w ^ 1 y 2 + 1 Λ 2 2 w ^ 1 z 2
As the momentum equations for disturbances in the y and z directions are linear and depend on only one velocity component, the final result is independent of the chosen velocity scale, depending on two parameters, B and Λ . From the above equations, clearly B is a kind of Reynolds number that represents the ratio between the fluid’s inertia due to the vertical deformation of the bed and the fluid’s viscous resistance. Note that from the continuity equation for the perturbed velocity components (Equation (14)), the following relationship between the velocity scales is obtained: V / W ~ λ / h 0 = Λ 1 .
For 0 B 1 , Ω is a real number, with 0 Ω 1 , and the solutions for the stream function, velocity field, flow depth, and free surface deformation involve hyperbolic functions (Equations (47)–(49), (52), (53), and (59)). When B > 1 , Ω is an imaginary number and the hyperbolic functions become trigonometric ones (Equations (60)–(65)), changing the behavior of the solutions, as will be seen below. Defining Ω ~ = B 1 > 0 , = Ω ~ i and Ω ~ = ω λ .

4.1.1. Behavior of the Solutions When 0 B 1 ( Ω Real)

Taking the real part of Equation (44) and rendering it dimensionless results in:
ψ = ψ a β = e t Ω 2 + 1 cosh ( Λ Ω ) cosh ( Λ Ω z ) Ω 2 + 1 cosh Λ Ω 1 sin y
Similarly, we can define v and w as new dimensionless velocity components, such that
v = v a β = e t Ω sinh ( Λ Ω z ) Ω 2 + 1 cosh Λ Ω 1 sin y
w = w a β = e t Ω 2 + 1 cosh ( Λ Ω ) cosh ( Λ Ω z ) Ω 2 + 1 cosh ( Λ Ω ) 1 cos y
The transverse velocities have been scaled with α β and not with the scales V and W , solely to facilitate the graphical presentation of the results. Indeed, by having the product α β appear as a coefficient in the transverse velocities, generating the new variables v and w reduces the number of parameters to be considered.
Denoting A = δ / ε as the ratio between the deformation of the free surface and the bed, Equation (53) in terms of the dimensionless parameters is
A = ε h 1 ε η 1 = Ω 2 cosh ( Λ Ω ) Ω 2 + 1 cosh ( Λ Ω ) 1  
For 0 B 1 (i.e., 0 Ω 1 ), A is in the range 0 A 1 for any value of Λ 0 . This means that the amplitude of the free surface deformation is always smaller than the bed deformation. Also, the term in brackets in Equations (73) and (74) is always less than 1, for any z 1 .
Finally, the dimensionless perturbation of the longitudinal component of the velocity is as follows:
u 1 = u 1 ρ g sin θ h 0 2 / μ = e t cos ( y ) Ω 2 1 Ω 2 { Ω 2 + 1 cosh Λ Ω Ω 2 + 1 cosh Λ Ω 1 × 1 z cosh ( Λ Ω z ) + 1 4 1 Ω 2 + 1 cosh Λ Ω 1 z cosh ( Λ Ω z ) + 1 Ω 2 + 1 cosh Λ Ω 1 1 2 Λ Ω z 1 2 z 2 sinh ( Λ Ω z ) + [ 1 cosh ( Λ Ω ) Ω 2 + 1 cosh Λ Ω Ω 2 + 1 cosh Λ Ω 1 1 Λ Ω + sinh Λ Ω 1 4 1 Ω 2 + 1 cosh Λ Ω 1 × 1 Λ Ω + sinh Λ Ω cosh ( Λ Ω ) + Λ Ω ] sinh ( Λ Ω z ) }
It is easier to realize the analytical results presented above if they are plotted as in Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8. It must be noted that the general behavior of the streamlines, velocity components, and deformation of the free surface essentially does not change with the values of Ω and Λ , and the general shapes are those presented in the figures. Arbitrary values of these parameters were used in the computations. The spanwise dimension ( y -axis) of the plots correspond to a wavelength λ of the bed deformation. The obtained results require small values of ε to maintain the validity of the linear approximation. However, using ε < 0.1 produces results with variations that are practically imperceptible in the graphs. Thus, to facilitate visualizing the variations introduced by the first-order term in the base solution, the amplitudes of the perturbations were greatly amplified in the figures depicting flow properties.

4.1.2. Behavior of the Solutions When B 1 ( Ω Imaginary)

By employing relationships that transform hyperbolic functions of a complex variable into trigonometric ones and taking their real parts, we obtain the dimensionless versions of the equations for the streamline, velocity components, perturbed deformation of the free surface, and the ratio between the deformation of the free surface and the bed. They are presented below:
ψ * = ψ a β = e t * 1 Ω ~ 2 cos ( Λ Ω ~ ) cos ( Λ Ω ~ z ) 1 Ω ~ 2 cos Λ Ω ~ 1 sin y
v = v a β = e t Ω ~ sin ( Λ Ω ~ z ) 1 Ω ~ 2 cos Λ Ω ~ 1 sin y
w = w a β = e t 1 Ω ~ 2 cos ( Λ Ω ~ ) cos ( Λ Ω ~ z ) 1 Ω ~ 2 cos Λ Ω ~ 1 cos y
ε h 1 = ε 1 e t cos ( y ) Ω ~ 2 cos ( Λ Ω ~ ) 1 Ω ~ 2 cos ( Λ Ω ~ ) 1  
A = ε h 1 ε η 1 = Ω ~ 2 cos ( Λ Ω ~ ) 1 Ω ~ 2 cos ( Λ Ω ~ ) 1  
It is worth discussing the behavior of A before analyzing the stream function, velocities, and free surface deformation. A becomes indeterminate when 1 Ω ~ 2 cos ( Λ Ω ~ ) 1 = 0 , fluctuating between positive and negative values, as it is observed in Figure 9. Although the figure corresponds to the particular case Λ = 2 , the general behavior is the same for any Λ . The change in sign of A indicates that the free surface deformation is not always in phase with the bed deformation. Furthermore, due to viscous dissipation, it is not possible to have values of the amplitude of the free surface greater than the amplitude of the bed, limiting the range of A to the interval 1 < A < 1 . This last condition is important because it imposes a relation between Λ and Ω ~ . For the upper limit A = 1 , Λ Ω ~ = 2 π n with n = 0,1 , 2 ,   must be satisfied. For the lower limit A = 1 , the relationship between Λ and Ω ~ is given by 2 Ω ~ 2 1 cos Ω ~ = 1 . For a given interval of Ω ~ , the larger the value of Λ , the greater the number of branches of A .
The behavior of ψ when Ω is imaginary is completely different from the case when it is real. The presence of an oscillatory function like cosine could lead to a change in the sign of ψ along z for a given transversal location y . In contrast to the case of real Ω , many configurations of streamlines arise depending on the values of Λ and Ω ~ , as shown in the examples in Figure 10. This sometimes generates cells to make compatible the velocities induced simultaneously by both bottom and free surface deformations.
The dimensionless transverse velocities v and w are presented in Figure 11 and Figure 12, respectively, for the three values of Ω ~ corresponding to Λ = 10 shown in Figure 10. The change in sign of the velocities is a result of the presence of circulation cells. The deformations of the free surface and bottom over time for Λ = 10 are shown in Figure 13. As it was indicated, both deformations can be in phase or out of phase by π , depending on the value of Ω ~ ( A can be positive or negative). Finally, the dimensionless main flow velocity u is presented in Figure 14, where the effect of the perturbation on the base velocity has been exaggerated for easier visualization.

4.2. Flow Separation

The presence of a flow over protrusions on the channel bed may lead to flow separation. Currently, there is a lack of studies determining the separation conditions for the specific configuration explored in this article—convergent flows over a sinusoidally undulated bed with a free surface. Nevertheless, as an initial approximation, extrapolating results from studies on pressure-driven flows in two-dimensional pipes with sinusoidally varying walls along the main flow direction is feasible. To estimate the separation condition for transverse flow with a free surface, we rely on a previous result obtained for flows constrained by wavy walls. This estimation is made despite differences in flow configurations, such as the constant net flow existing in the x -direction in pressure-driven cases (where separation occurs in the x z plane), which contrasts with the variable cross-sectional flow ( y z plane) of the problem addressed in this article.
In this regard, the analytical results of Tsangaris and Leiter [35] based on the separation of flow from a laminar boundary layer over a sinusoidally varying surface is applied. These authors determined that the separation is conditioned by a critical Reynolds number ( R e c ) that is related to the dimensionless amplitude of bed deformation ( ε ). In general, this relationship is not simple and must be determined by numerically integrating the vorticity equation in terms of the stream function, obtaining a relationship R e c = φ ( ε ) , which they graphically present in their article. The critical Reynolds number at which separation occurs is defined as R e c = U b / ν , where b is half of the mean separation between walls, and U is the average velocity. ε = a / b , with a being the amplitude of the wall deformation. For large Reynolds numbers, the authors provide the approximation φ     7.563 / ε ^ 2 .
To apply the separation point condition of Tsangaris and Leiter to the cases generated in our study, the velocity U is taken as the maximum transverse velocity v m a x , and b is set equal to the flow height coordinate, z m a x . Thus, the critical Reynolds number is given by R e c = ρ v m a x z m a x / μ   , which, expressed in terms of the dimensionless variables given by Equations (66) and (67), is written as R e c = ε B Λ 2 v m a x z V m a x . For 0 < B     1 and considering the most unfavorable scenario ( sin y 1   ,   z m a x = 1   ,   e t 1 ) , an upper bound for R e c is obtained:
R e c < ε 2 1 Ω 2 Λ 2 Ω sinh ( Λ Ω ) Ω 2 + 1 cosh Λ Ω 1 < φ ε
Similarly, for the case B > 1 , the flow separation condition is bounded by
R e c < ε 2 1 + Ω ~ 2 Λ 2 Ω ~ 1 Ω ~ 2 cos Λ Ω ~ 1 < φ ε
As an example, for the case 0 < B 1 , if ε = 0.1 ,   Λ = 10 , and Ω = 0.6 , the result is R e c < 0.2 . Similarly, when B > 1 , considering ε = 0.1 ,   Λ = 100 , and Ω ~ 0.1 , the result is R e c < 5.5. In both cases, the values of R e c are much lower than φ ( ε ) , obtained from the graphical relationship presented in the article by Tsangaris and Leiter [35].
Choosing ε = 0.1 represents merely an arbitrary value for the dimensionless amplitude of the bed deformation ( ε = a / h 0 ). While the preceding results confirm the absence of flow separation, they do not guarantee the extent to which the linear solution deviates from the one obtained by accounting for the non-linear terms in the momentum equations. In strict terms, linearizing the problem requires ε 1 . Therefore, it is reasonable to expect that the results obtained from the linear approach are valid for ε < 0.01 . However, this a priori assumption needs validation through a nonlinear analytical solution of the problem or via numerical modeling. Both approaches are beyond the scope of this work, and the authors are not aware of this problem being addressed previously.

4.3. Conditions for the Validity of the Results Obtained from the Linearized Equations

In order to determine the limitations of the applicability of the first-order results obtained in this article, the conditions that the parameters must satisfy when non-dimensionalizing the linearized equations are established below. The analysis can be carried out interchangeably for the momentum equations of perturbations in the y or z directions, taking into account that the velocity scales for both components are related as V / W ~ λ / h 0 .
The momentum equation for the z-component of the disturbance for a uniform flow is given by:
w 1 t + v 1 w 1 y + w 1 w 1 z = μ ρ 2 w 1 y 2 + 2 w 1 z 2
Considering the scales mentioned before, v 1 ~ V ,   w 1 ~ W ,   y ~ λ ,   z ~ h 0 , and t ~ 1 / β , the scalings for the different terms of Equation (79) are obtained:
w 1 t ~ β W , v 1 w 1 y ~ V W λ ~ W 2 h 0 , w 1 w 1 z ~ W 2 h 0 2 w 1 y 2 ~ W λ 2 , 2 w 1 z 2 ~ W h 0 2
To neglect the nonlinear terms means:
v 1 w 1 y + w 1 w 1 z w 1 t
That is to say:
v 1 w 1 y + w 1 w 1 z w 1 t 1
In terms of the scales:
W 2 / h 0 β W 1
The velocity in the z -direction is induced by the bed deformation (Equation (36)), so it is permissible to consider W ~ a β . Then, inequality (88) is reduced to:
a h 0 1
The previous relationship corresponds to ε 1 and is nothing more than the condition of small amplitude of bed deformation that has been imposed from the beginning and is the basis for the linearization of the equations. The nonlinear terms of the momentum flow are also non-negligible compared to the forces of viscous origin, that is:
v 1 w 1 y + w 1 w 1 z μ ρ 2 w 1 y 2 + 2 w 1 z 2
The scaling of relation (90) results in:
W 2 / h 0 μ W ρ λ 2 1 + λ / h 0 2 = ρ β λ 2 μ a h 0 1 1 + λ / h 0 2 1
B a h 0 1 1 + λ / h 0 2 1
To simplify the writing, we define G = 1 / 1 + λ / h 0 2 . Given condition (89), inequality (92) can be satisfied in various ways:
i.  
B ~ 1 and G ~ 1 . The second condition implies λ / h 0 1 or even λ / h 0 1 , a condition that corresponds to deep waters.
ii. 
B > 1 and G 1 , such that B ε ~ 1 . The condition G 1 means that λ / h 0 1 , an approximation of long waves. The condition B > 1 includes B 1 , which is interpreted as a transverse flow with a large Reynolds number but within the laminar regime.
iii. 
B < 1 and G 1 , such that G ε ~ 1 . This case cannot occur because, to have G 1 , its denominator must be less than one.
iv. 
B < 1 and G ~ 1 . This case implies λ / h 0 1 and that viscous effects dominate in the transverse flow.
Note that the conditions B < 1 and B > 1 correspond to the cases analyzed in Section 4.1.1 and Section 4.1.2, respectively. The condition for deep waters λ / h 0 1 , when applied to shallow depths means very short wavelengths and the effect of surface tension should be included in the dynamic boundary condition at the free surface.

5. Discussion

The steady and unsteady flow in ducts with wavy walls driven by pressure gradients as well as the free-surface flow with a bed deformed in the downstream direction have been widely studied. In this study, we have examined free-surface flow, determining the characteristics of the velocity field and the surface when the channel bed deforms over time and in the direction transverse to the flow. To the best of the authors’ knowledge, this has not been studied before.
The pattern of the perturbation component of the velocity field, streamlines, and deformation of the free surface depends on the discriminant of Ω , given by ( 1 B ) , where B = ( β ρ λ 2 ) / μ > 0 . The parameter B is a kind of Reynolds number that represents the fluid inertia induced by bed deformation relative to viscous resistance to motion. Thus, when the viscous effect dominates, B < 1 , and the root of Ω is real. This results in a monotonic decay from the source of motion (bed deformation) to the free surface, generating simple and generalizable patterns of streamlines, velocity components, and surface deformation, as observed in Figure 4, Figure 5, Figure 6 and Figure 7. In this case, the deformation of the free surface is always in phase with the bed deformation.
More interesting and cumbersome to analyze is the case when fluid inertia dominates over viscous resistance, as indicated by B > 1 and the root of Ω being imaginary. This incorporates trigonometric functions into the solutions, leading to an oscillatory behavior of the perturbation velocity components. Consequently, circulation cells of flow transverse to the main flow direction can be observed, as illustrated in Figure 10, Figure 11, Figure 12 and Figure 13. With an imaginary Ω , the solution for the ratio of the amplitudes of free surface and bed deformations ( A ) is unbounded. However, due to viscous dissipation, the free surface deformation cannot exceed that of the bed, necessitating that 1 < A < 1 . This limitation constrains the range of Λ and Ω to ensure physically meaningful solutions. Depending on the relative importance of inertia to viscosity, the deformation of the free surface may be in phase with the bed deformation or out of phase by π . In the latter case, circulation cells are generated, serving as a way to reconcile the different signs imposed on the vertical velocity component ( w ) by the deformation of both the free surface and the bed.
In this work, the bed deformation is described by the perturbation η 1 = h 0 f η e i k y with f η = 1 e β t (Equations (37) and (38)). However, the perturbation η 1 can be generalized for any function whose variation in the y -direction can be expressed in terms of Fourier series, with f η being a bounded function, while maintaining the restriction that ε η 1 is small.
Finally, it is emphasized that the velocity fields and streamlines obtained here are first-order approximations of the nonlinear equations of motion. Therefore, the results are valid to the extent that the disturbances are small. The removal of this restriction can now be done without major difficulties if the problem is approached through proper numerical modeling, which is free from the constraints imposed by the linearization of equations, enabling the calculation of the velocity field and deformation of the free surface (including the effect of surface tension) and the identification of conditions leading to flow separation in real cases. Furthermore, numerical modeling is not limited to laminar flow, as is the case with most analytical approaches, and can also address turbulent flows. However, in the authors’ judgment, the approximate and linearized analytical solution can offer a clearer and more conceptual initial understanding of the system’s behavior. This is evident in this work through the role played by the parameter B , interpreted as a Reynolds number associated with the flow inertia induced by bed deformation. This parameter arises easily and naturally during the nondimensionalization of the (linear) equations governing first-order perturbations. Depending on whether B is greater or less than 1, the behavior of the solutions can be very different.

Author Contributions

Conceptualization, R.G. and A.T.; methodology, R.G. and A.T.; formal analysis, R.G. and A.T.; investigation, R.G. and A.T.; writing—original draft preparation and editing, A.T.; writing—review, R.G.; figures, R.G.; supervision, A.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

R.G. acknowledges the fellowship awarded by the Chilean Agency for Research and Development (ANID) to pursue doctoral studies through the grant ANID-Subdirección de Capital Humano/Doctorado Nacional/2022-21221030. A.T. acknowledges the project ANID AFB230001 and the Department of Civil Engineering, University of Chile.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Tamburrino, A. From Navier to Stokes: Commemorating the Bicentenary of Navier’s Equation on the Lay of Fluid Motion. Fluids 2024, 9, 15. [Google Scholar] [CrossRef]
  2. Navier, C.J.M.H. Mémoire sur les lois du movement des fluides. Mem. Acd. R. Sci. Paris 1823, 6, 389–416. [Google Scholar]
  3. Stokes, G.G. On the Effect of the Internal Friction of Fluids on the Motion of Pendulums. Trans. Cambridge Phil. Soc. 1851, 9 Pt II, 8–106. [Google Scholar]
  4. Wang, C.Y. Liquid Film Flowing Slowly Down a Wavy Incline. AlChE J. 1981, 27, 207–212. [Google Scholar] [CrossRef]
  5. Huang, C.-J.; Dong, C.M. Propagation of Water Waves over Rigid Rippled Beds. J. Waterw. Port Coast. Ocean Eng. 2002, 128, 190–201. [Google Scholar] [CrossRef]
  6. Sobey, I.A. On flow through furrowed channels: Calculated flow patterns. J. Fluid Mech. 1980, 96 Pt 1, 1–26. [Google Scholar] [CrossRef]
  7. Wierschem, A.; Scholle, M.; Aksel, N. Comparison of different theoretical approaches to experiments on film flow down an inclined wavy channel. Exp. Fluids 2002, 33, 429–442. [Google Scholar] [CrossRef]
  8. Wierschem, A.; Lepski, C.; Aksel, N. Effect of long undulated bottoms on thin gravity-driven films. Acta Mech. 2005, 179, 41–66. [Google Scholar] [CrossRef]
  9. Benjamin, T.B. Shearing flow over a wavy boundary. J. Fluid Mech. 1959, 6, 161–205. [Google Scholar] [CrossRef]
  10. Lyne, W.H. Unsteady viscous flow over a wavy wall. J. Fluid Mech. 1971, 50 Pt 1, 33–48. [Google Scholar] [CrossRef]
  11. Kaneko, A.; Honji, H. Bouble structures of steady streaming in the oscillatory viscous flow over a wavy wall. J. Fluid Mech. 1979, 93 Pt 4, 727–736. [Google Scholar] [CrossRef]
  12. Vittori, G. Non-linear viscous oscillatory flow over a small amplitude wavy wall. J. Hy. Res. 1989, 27, 267–280. [Google Scholar] [CrossRef]
  13. Pozrikidis, C. The flow of a liquid film along a periodic wall. J. Fluid Mech. 1988, 188, 275–300. [Google Scholar] [CrossRef]
  14. Bontozoglou, V.; Papapolymerou, G. Laminar film flow down a wavy incline. Int. J. Multiph. Flow. 1996, 23, 69–79. [Google Scholar] [CrossRef]
  15. Scholle, M.; Wierschem, A.; Aksel, N. Creeping films with vortices over strongly undulated bottoms. Acta Mech. 2004, 168, 167–193. [Google Scholar] [CrossRef]
  16. Trifonov, Y.Y. Stability of a viscous liquid film flowing down a periodic surface. Int. J. Multiph. Flow. 2007, 33, 1186–1204. [Google Scholar] [CrossRef]
  17. Trifonov, Y.Y. Stability of a viscous liquid film flowing down an inclined corrugated plate: The direct Navier-Stokes computations and Floquet theory. Phys. Fluids 2014, 26, 114101. [Google Scholar] [CrossRef]
  18. Trifonov, Y. Nonlinear waves on a liquid film falling down an inclined corrugated surface. Phys. Fluids 2017, 29, 054104. [Google Scholar] [CrossRef]
  19. Tseluiko, D.; Blyth, M.G.; Papageorgiou, D.T. Stability film flow over inclined topography based on a long-wave nonlinear model. J. Fluid Mech. 2013, 729, 638–671. [Google Scholar] [CrossRef]
  20. Dávalos-Orozco, L.A. Nonlinear instability of a thin film flowing down a smoothy deformed surface. Phys. Fluids. 2007, 19, 074103. [Google Scholar] [CrossRef]
  21. Hara, T.; Mei, C.C. Oscillating flows over periodic ripples. J. Fluid Mech. 1990, 211, 183–209. [Google Scholar] [CrossRef]
  22. Nishimura, T.; Miyashita, H.; Murakami, S.; Kawamura, S. Oscillatory flow in a symmetric sinusoidal wavy-walled channel at intermediate Strouhal numbers. Chem. Eng. Sci. 1991, 46, 757–771. [Google Scholar] [CrossRef]
  23. Shen, L.; Chan, E.-S. Numerical simulation of oscillatory flows over a rippled bed by immersed boundary method. Appl. Ocean Res. 2013, 43, 27–36. [Google Scholar] [CrossRef]
  24. Szumbarski, J. Instability of viscous incompressible flow in a channel with transversely corrugated walls. J. Theor. Appl. Mech. 2007, 45, 659–683. [Google Scholar]
  25. Kowaleski, T.A.; Szumbarski, J.; Błoński, S. Low-Reynolds-number instability of the laminar flow between wavy walls. In Proceedings of the Sixth International ASME Conference on Nanochannels, Microchannels and Minichannels, ICNMM 2008, Darmstadt, Germany, 23–25 June 2008. [Google Scholar]
  26. Szumbarski, J.; Błoński, S. Destabilization of a laminar flow in a rectangular channel by transversely oriented wall corrugation. Arch. Mech. 2011, 63, 393–428. [Google Scholar]
  27. Moradi, H.V.; Floryan, J.M. Stability of flow in a channel with longitudinal grooves. J. Fluid Mech. 2014, 757, 613–648. [Google Scholar] [CrossRef]
  28. Szumbarski, J.; Błoński, S.; Kowaleski, T.A. Impact of transversely oriented wall corrugation on hydraulic resistance of a channel flow. Arch. Mech. Eng. 2011, 58, 441–466. [Google Scholar] [CrossRef]
  29. Gepner, S.W.; Yadav, N.; Szumbarski, J. Secondary flows in a longitudinal grooved channel and enhancement of diffusive transport. Int. J. Heat Mass Transf. 2020, 153, 119523. [Google Scholar] [CrossRef]
  30. Yadav, N.; Gepner, S.W. Slowing down convective instabilities in corrugated Couette-Poiseuille flow. J. Fluid Mech. 2022, 959, A5. [Google Scholar] [CrossRef]
  31. Gepner, S.W.; Floryan, J.M. Flow dynamics in sinusoidal channels at moderated Reynolds numbers. J. Fluid Mech. 2023, 972, A22. [Google Scholar] [CrossRef]
  32. Smith, F.T. flow through constricted or dilated pipes and channels: Part 1. Q. J. Mech. Appl. Math. 1976, 29 Pt 3, 343–364. [Google Scholar] [CrossRef]
  33. Bordner, G.L. Nonlinear analysis of laminar boundary layer flow over a periodic wavy surface. Phys. Fluids 1978, 21, 1471–1474. [Google Scholar] [CrossRef]
  34. Caponi, E.A.; Fornberg, B.; Knight, D.D.; McLean, J.W.; Saffman, P.G.; Yuen, H.C. Calculations of laminar viscous flow over a moving wavy surface. J. Fluid Mech. 1982, 124, 347–362. [Google Scholar] [CrossRef]
  35. Tsangaris, S.; Leiter, E. On laminar steady flow in sinusoidal channels. J. Eng. Math. 1984, 18, 89–103. [Google Scholar] [CrossRef]
  36. Mahmud, S.; Sadrul Islam, A.K.M.; Mamun, M.A.H. Separation characteristics of fluid flow inside two parallel plates with wavy surface. Int. J. Eng. Sci. 2002, 40, 1495–1509. [Google Scholar] [CrossRef]
  37. Carlsson, F.; Sen, M.; Löfdahl, L. Steady streaming due to vibrating walls. Phys. Fluids. 2004, 16, 1822–1825. [Google Scholar] [CrossRef]
  38. Carlsson, F.; Sen, M.; Löfdahl, L. Fluid mixing induced by vibrating walls. Eur. J. Mech. B Fluids 2005, 24, 366–378. [Google Scholar] [CrossRef]
  39. Lebbal, S.; Alizard, F.; Pier, B. Revisiting the linear instabilities of plane channel flow between compliant walls. Phys. Rev. Fluids 2022, 7, 023903. [Google Scholar] [CrossRef]
  40. Vladimirov, V.A. Viscous Flows in a Half Space Caused by Tangential Vibrations on its Boundary. Stud. Appl. Math. 2008, 121, 337–367. [Google Scholar] [CrossRef]
  41. Izumi, N.; Fuji, K. Channelization on plateaus composed of weakly cohesive fine sediment. J. Geophys. Res. 2006, 111, F01012. [Google Scholar] [CrossRef]
  42. Granger, R.A. Fluid Mechanics; Dover Publications Inc.: New York, NY, USA, 1995. [Google Scholar]
  43. Nayfeh, A.H. Perturbation Methods; WILEY-VCH Verlag GmbH & CO. KGaA: Weinheim, Germany, 2004; Chapter 2. [Google Scholar]
  44. Shivamoggi, B.K. Perturbation Methods for Differential Equations; Birkhauser: Boston, MA, USA, 2003. [Google Scholar]
Figure 2. Sketch of the problem addressed in this article: determination of the velocity field and free-surface distortion due to the temporal variation of the bottom. The main flow is perpendicular to the sheet.
Figure 2. Sketch of the problem addressed in this article: determination of the velocity field and free-surface distortion due to the temporal variation of the bottom. The main flow is perpendicular to the sheet.
Mathematics 12 00394 g002
Figure 3. Sketch of the problem and definition of the coordinate system.
Figure 3. Sketch of the problem and definition of the coordinate system.
Mathematics 12 00394 g003
Figure 4. Streamlines for an arbitrary time t > 0 . Although the plot corresponds to the parameters Λ = 10 and Ω = 0.1 (representing the ratio of free surface to bed deformation A = 0.0276 ), the general shape of the streamlines is the same for any value of the triplet ( Λ , Ω , t ) .
Figure 4. Streamlines for an arbitrary time t > 0 . Although the plot corresponds to the parameters Λ = 10 and Ω = 0.1 (representing the ratio of free surface to bed deformation A = 0.0276 ), the general shape of the streamlines is the same for any value of the triplet ( Λ , Ω , t ) .
Mathematics 12 00394 g004
Figure 5. Distribution of the dimensionless transverse velocity v along z , at different y locations for an arbitrary time t , for Λ = 10 and three values of Ω . The general behavior of the transverse velocity is the same for other values of the triplet ( Λ , Ω , t ) . For y = 0 ,   π , and 2 π the component v is zero along the z -axis (from the bottom to the free surface). The locations y = 0 and y = 2 π correspond to the maximum bed deformation. y = π corresponds to the lowest point of the bottom. For t , the velocity v is zero in all the flow domain. Dashed color lines correspond to zero velocity at the location where the line is drawn, while continuous color lines represent the velocity at time t .
Figure 5. Distribution of the dimensionless transverse velocity v along z , at different y locations for an arbitrary time t , for Λ = 10 and three values of Ω . The general behavior of the transverse velocity is the same for other values of the triplet ( Λ , Ω , t ) . For y = 0 ,   π , and 2 π the component v is zero along the z -axis (from the bottom to the free surface). The locations y = 0 and y = 2 π correspond to the maximum bed deformation. y = π corresponds to the lowest point of the bottom. For t , the velocity v is zero in all the flow domain. Dashed color lines correspond to zero velocity at the location where the line is drawn, while continuous color lines represent the velocity at time t .
Mathematics 12 00394 g005
Figure 6. Distribution of the dimensionless transverse velocity w along y , at different height locations z for an arbitrary time t , for Λ = 10 and three values of Ω . This velocity component is in phase with the bottom deformation. The locations y = 0 and y = 2 π correspond to the maximum bed deformation and the maximum value of velocity w , whereas y = π corresponds to the lowest point of the bottom and the minimum value of w . The amplitude of the velocity w decreases as the location gets closer to the free surface. For t , when the bottom has reached its equilibrium amplitude, w is zero in all the flow domain. Dashed color lines correspond to zero velocity at the location where the line is drawn, while continuous color lines represent the velocity at time t .
Figure 6. Distribution of the dimensionless transverse velocity w along y , at different height locations z for an arbitrary time t , for Λ = 10 and three values of Ω . This velocity component is in phase with the bottom deformation. The locations y = 0 and y = 2 π correspond to the maximum bed deformation and the maximum value of velocity w , whereas y = π corresponds to the lowest point of the bottom and the minimum value of w . The amplitude of the velocity w decreases as the location gets closer to the free surface. For t , when the bottom has reached its equilibrium amplitude, w is zero in all the flow domain. Dashed color lines correspond to zero velocity at the location where the line is drawn, while continuous color lines represent the velocity at time t .
Mathematics 12 00394 g006
Figure 7. Deformation of the bottom and the free surface as functions of time. For t = 0 , both the bottom and the free surface are plane surfaces ( z = 0 and z = 1 , respectively). As time increases, so do the deformation amplitudes, tending asymptotically to a = ε h 0 for the bottom, and to δ h 0 < a for the free surface. The different curves show the deformations for the dimensionless times t = 0 ,   0.2 ,   0.6 ,   2 , and 10. Since the deformations of the bottom and the free surface are on the same scale, the deformation of the free surface is not noticeable in the plot for Ω = 0.1 .
Figure 7. Deformation of the bottom and the free surface as functions of time. For t = 0 , both the bottom and the free surface are plane surfaces ( z = 0 and z = 1 , respectively). As time increases, so do the deformation amplitudes, tending asymptotically to a = ε h 0 for the bottom, and to δ h 0 < a for the free surface. The different curves show the deformations for the dimensionless times t = 0 ,   0.2 ,   0.6 ,   2 , and 10. Since the deformations of the bottom and the free surface are on the same scale, the deformation of the free surface is not noticeable in the plot for Ω = 0.1 .
Mathematics 12 00394 g007
Figure 8. Dimensionless velocity component in the x direction ( u ). The effect of the perturbation has been amplified to better visualize the distortion of the velocity distribution in the y and z directions. The plot corresponds to Λ = 10 and Ω = 0.6 , but the general shape of u is similar for other values of Λ and Ω . The distortion due to the perturbation decreases as the dimensionless time t increases.
Figure 8. Dimensionless velocity component in the x direction ( u ). The effect of the perturbation has been amplified to better visualize the distortion of the velocity distribution in the y and z directions. The plot corresponds to Λ = 10 and Ω = 0.6 , but the general shape of u is similar for other values of Λ and Ω . The distortion due to the perturbation decreases as the dimensionless time t increases.
Mathematics 12 00394 g008
Figure 9. Ratio between the amplitude of the free surface and the bed according to Equation (77). Values of A with physical meaning lie within the range 1 < A < 1 , indicated by the horizontal red lines. This signifies that not every value of Ω ~ > 0 is valid.
Figure 9. Ratio between the amplitude of the free surface and the bed according to Equation (77). Values of A with physical meaning lie within the range 1 < A < 1 , indicated by the horizontal red lines. This signifies that not every value of Ω ~ > 0 is valid.
Mathematics 12 00394 g009
Figure 10. Streamlines for different values of Λ ,   Ω ~ , and A . Notice that A is not an independent parameter, but it is related to Λ and Ω ~ according to Equation (77).
Figure 10. Streamlines for different values of Λ ,   Ω ~ , and A . Notice that A is not an independent parameter, but it is related to Λ and Ω ~ according to Equation (77).
Mathematics 12 00394 g010
Figure 11. Distribution of the dimensionless transverse velocity v along z at different y locations for an arbitrary time t for Λ = 10 . The pattern reflects the presence of circulation cells arising from different values of Ω ~ , as shown in Figure 10. Dashed color lines correspond to zero velocity at the location where the line is drawn, while continuous color lines represent the velocity at time t .
Figure 11. Distribution of the dimensionless transverse velocity v along z at different y locations for an arbitrary time t for Λ = 10 . The pattern reflects the presence of circulation cells arising from different values of Ω ~ , as shown in Figure 10. Dashed color lines correspond to zero velocity at the location where the line is drawn, while continuous color lines represent the velocity at time t .
Mathematics 12 00394 g011
Figure 12. Distribution of the dimensionless transverse velocity w along y at different height locations for an arbitrary time t for Λ = 10 . Unlike the case of real Ω , depending on the value of Ω ~ , this velocity component can be in phase with the bottom deformation ( Ω ~ = 0.1 ), or out of phase ( Ω ~ = 0.6 ,   1.0 ). The amplitude of the velocity w decreases as the location gets closer to the free surface. For t , when the bottom has reached its equilibrium amplitude, w is zero throughout the entire flow domain. Dashed color lines correspond to zero velocity at the location where the line is drawn, while continuous color lines represent the velocity at time t .
Figure 12. Distribution of the dimensionless transverse velocity w along y at different height locations for an arbitrary time t for Λ = 10 . Unlike the case of real Ω , depending on the value of Ω ~ , this velocity component can be in phase with the bottom deformation ( Ω ~ = 0.1 ), or out of phase ( Ω ~ = 0.6 ,   1.0 ). The amplitude of the velocity w decreases as the location gets closer to the free surface. For t , when the bottom has reached its equilibrium amplitude, w is zero throughout the entire flow domain. Dashed color lines correspond to zero velocity at the location where the line is drawn, while continuous color lines represent the velocity at time t .
Mathematics 12 00394 g012
Figure 13. Deformation of the bottom and the free surface in function of time for Λ = 10 and the three values of Ω ~ presented in Figure 10. For t = 0 , both the bottom and the free surface are plane surfaces ( z = 0 and z = 1 , respectively). The different curves show the deformations for the dimensionless times t = 0 ,   0.2 ,   0.6 ,   2 , and 10. As time increases, so do the deformations. Distortion of the free surface is not noticeable in the plot for Ω = 0.1 . Depending on the value of Ω ~ , both deformations can be in phase or out of phase by π .
Figure 13. Deformation of the bottom and the free surface in function of time for Λ = 10 and the three values of Ω ~ presented in Figure 10. For t = 0 , both the bottom and the free surface are plane surfaces ( z = 0 and z = 1 , respectively). The different curves show the deformations for the dimensionless times t = 0 ,   0.2 ,   0.6 ,   2 , and 10. As time increases, so do the deformations. Distortion of the free surface is not noticeable in the plot for Ω = 0.1 . Depending on the value of Ω ~ , both deformations can be in phase or out of phase by π .
Mathematics 12 00394 g013
Figure 14. Dimensionless velocity component in the x direction ( u ) for Λ = 10 and Ω ~ = 1 . The effect of the perturbation has been amplified to better visualize the distortion of the velocity distribution in the y and z directions. The wavy shape of the velocity distribution depends on the number of circulation cells, which is determined by Ω ~ . The distortion due to the perturbation decreases as the dimensionless time t increases.
Figure 14. Dimensionless velocity component in the x direction ( u ) for Λ = 10 and Ω ~ = 1 . The effect of the perturbation has been amplified to better visualize the distortion of the velocity distribution in the y and z directions. The wavy shape of the velocity distribution depends on the number of circulation cells, which is determined by Ω ~ . The distortion due to the perturbation decreases as the dimensionless time t increases.
Mathematics 12 00394 g014
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

González, R.; Tamburrino, A. Velocity Field due to a Vertical Deformation of the Bottom of a Laminar Free-Surface Fluid Flow. Mathematics 2024, 12, 394. https://doi.org/10.3390/math12030394

AMA Style

González R, Tamburrino A. Velocity Field due to a Vertical Deformation of the Bottom of a Laminar Free-Surface Fluid Flow. Mathematics. 2024; 12(3):394. https://doi.org/10.3390/math12030394

Chicago/Turabian Style

González, Rodrigo, and Aldo Tamburrino. 2024. "Velocity Field due to a Vertical Deformation of the Bottom of a Laminar Free-Surface Fluid Flow" Mathematics 12, no. 3: 394. https://doi.org/10.3390/math12030394

APA Style

González, R., & Tamburrino, A. (2024). Velocity Field due to a Vertical Deformation of the Bottom of a Laminar Free-Surface Fluid Flow. Mathematics, 12(3), 394. https://doi.org/10.3390/math12030394

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop