**Evolution of Small Perturbations of the Free Surface of Viscous Fluid in the Stokes Approximation**

[17] Morrison, P.J.; (2009). Thoughts on Brackets and Dissipation: Old and New, *Journal of*

[18] Morrison, P. J. & Greene, J. M.; (1980). Noncanonical Hamiltonian Density Formula‐ tion of Hydrodynamics and Ideal Magnetohydrodynamics, *Physical Review Letters*,

[19] Padhye, N. S. (May 1998). Topics in Lagrangian and Hamiltonian fluid dynamics: Re‐ labeling symmetry and ion-acoustic wave stability, dissertation presented to the Fac‐ ulty of the Graduate School of The University of Texas at Austin in Partial Fullment

[20] Padhye, N. and P. J. Morrison (1996 a). Relabeling Symmetries in Hydrodynamics

[21] Padhye, N. and P. J. Morrison (1996 b). Fluid Element Relabeling Symmetry, *Physics*

[22] Rai Choudhuri, A. (1998). The Physics of Fluids and Plasmas-an introduction for as‐

[23] Zemansky, M. W. and Dittman, R. H. (1996). Heat and Thermodynamics, Mcgraw-

*Physics: Conference Series*, pp. 1-12, Vol. 169, 012006.

pp. 790-793, Vol. 45, No. 10, September 8.

of the Requirements for the Degree of PhD.

trophysicists, Cambridge University Press.

*Letters* 219A, 287.

138 Hydrodynamics - Concepts and Experiments

Hill College; 7th edition.

and Magnetohydrodynamics, *Plasma Physics Reports* 22, 869.

Tamara Pogosian, Ivan Melikhov, Anastasia Shutova, Alexey Amosov and Sergey Chivilikhin

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/59226

## **1. Introduction**

Stokes steady-state approximation [1] is widely used for calculation of viscous flows with low Reynolds numbers. This approach was first applied to describe the flow with free boundaries in [2]. Many classical problems of high-viscous liquids hydrodynamics were solved in this approximation: the destruction of bubbles [3], the formation of the cusp on the liquid surface [4,5], the coalescence of liquid particles [6] and some problems of nanofluidics [7].

The idea of using inertia-less approximation to calculate wave relaxation on the surface of a fluid with high viscosity was suggested by [8]. The relaxation of harmonic perturbations in this approximation was considered by [9,10]. This approach has also been used in studying thin film flow [11], two-dimensional Stokes flow with free boundary [12-13].

In this work we study an aperiodic relaxation of a small localized perturbation of the planar surface of an incompressible fluid under the influence of gravity and surface tension. The description is given in the Stokes approximation. It is shown that this imposes limitations on the three-dimensional spectrum of the considered perturbations. An equation describing perturbation damping in Fourier space was obtained. It is shown that the volume of the perturbation becomes zero in the time that is small by comparison to the time of perturbation damping. In the short wave limit, the law of evolution of the perturbation permits a simple geometric interpretation.

The second part of the work is focused on theoretical investigation of instability of the surface shape under the influence of diffusion mass flux through the surface. A local eminence on the surface can lead to perturbation growth. Stability and instability of the perturbation will

depend on which mechanism prevails. Diffusion causes distortion of surface shape, gravita‐ tional and capillary forces have the opposite effect.

#### **2. Evolution of half-space free boundary perturbation**

#### **2.1. The general equations**

We consider a half-space filled with a highly viscous incompressible Newton liquid of constant viscosity. At the initial moment of time the planar surface of half-space obtains a small local perturbation. We will investigate further the evolution of this perturbation due to gravitation, capillary forces and external mass flux through the free surface. We introduce the Cartesian coordinate system with axes (*x*1, *x*2) located on the unperturbed surface and the axis *z* directed along the external normal to it. Then the perturbation of the surface can be defined as *z* =*h* (*x*, *t*), where *x* is the vector with components *x*1, *x*2 and *t* is the time. We express the equations of motion in the Stokes approximation [1], the continuity equation, the equation for the renormalized pressure *Π* =*P* + *ρgz* (where *P* is the true pressure) as

$$\begin{cases} \hat{\sigma}\_{xx}^{2}\mathbf{V}\_{x} + \hat{\sigma}\_{\text{f\beta}\emptyset}^{2}\mathbf{V}\_{x} = \frac{1}{\mu}\hat{\sigma}\_{x}\Pi, \\\\ \hat{\sigma}\_{xx}^{2}\mathbf{V}\_{\alpha} + \hat{\sigma}\_{\text{f\beta}\emptyset}^{2}\mathbf{V}\_{\alpha} = \frac{1}{\mu}\hat{\sigma}\_{\alpha}\Pi, \\\\ \hat{\sigma}\_{x}\mathbf{V}\_{x} + \hat{\sigma}\_{\emptyset}\mathbf{V}\_{\emptyset} = \mathbf{0}, \\\\ \hat{\sigma}\_{xx}^{2}\Pi + \hat{\sigma}\_{\emptyset\emptyset}^{2}\Pi = \mathbf{0}. \end{cases} \tag{1}$$

The indices *α*, *β* can be 1 or 2; summation over repeated indices is assumed. The boundary conditions which express the continuity of the normal and tangential stresses on the free surface, in the approximation linear by *h* , have the form

$$\left. \left( -\Gamma \mathbf{I} + 2\mu \hat{\boldsymbol{\sigma}}\_x \mathbf{V}\_x \right) \right|\_{\mathbf{z}=0} = \mathbf{c} \hat{\boldsymbol{\sigma}}\_{\beta\beta}^2 \mathbf{h} - \rho \mathbf{g} \mathbf{h}, \quad \left. \left( \hat{\boldsymbol{\sigma}}\_a \mathbf{V}\_x + \hat{\boldsymbol{\sigma}}\_x \mathbf{V}\_a \right) \right|\_{\mathbf{z}=0} = \mathbf{0}. \tag{2}$$

The following conditions in case of infinite half-space must also be met:

$$
\begin{pmatrix} \mathbf{V}\_x \\ \mathbf{V}\_\alpha \\ \mathbf{I} \end{pmatrix}\_{\mathbf{z}\to\mathbf{-}\mathbf{v}} = \mathbf{0}.\tag{3}
$$

Here *Vz* and *Vα* are the components of velocity, ∂*z* and ∂*α* are the derivatives along the coordinates *z* and *xα* accordingly, *g* is the acceleration due to gravity; *µ*, *ρ* and *σ* are the viscosity, density and surface tension coefficient of the liquid.

Equation of free surface evolution and its initial condition have the form

$$\left. \mathcal{C}\_t \mathbf{h} = \left. \mathbf{V}\_x \right|\_{x=0} + \mathbf{W}, \quad \left. \mathbf{h} \right|\_{t=0} = \mathbf{h}^0 \left( \mathbf{x} \right), \tag{4}$$

where *W* describes the influence of external mass fluxes on the velocity of free boundary.

Applying to (1)-(4) the two-dimensional Fourier transform with respect to the longitudinal coordinates

$$\mathbf{f\_{k}} = \left[ \mathbf{f} \left( \mathbf{x} \right) \mathbf{e}^{-i\mathbf{k}\mathbf{x}} \mathbf{d}^{2} \mathbf{x} \right]$$

we obtain

(1)

(3)

depend on which mechanism prevails. Diffusion causes distortion of surface shape, gravita‐

We consider a half-space filled with a highly viscous incompressible Newton liquid of constant viscosity. At the initial moment of time the planar surface of half-space obtains a small local perturbation. We will investigate further the evolution of this perturbation due to gravitation, capillary forces and external mass flux through the free surface. We introduce the Cartesian coordinate system with axes (*x*1, *x*2) located on the unperturbed surface and the axis *z* directed along the external normal to it. Then the perturbation of the surface can be defined as *z* =*h* (*x*, *t*), where *x* is the vector with components *x*1, *x*2 and *t* is the time. We express the equations of motion in the Stokes approximation [1], the continuity equation, the equation for

tional and capillary forces have the opposite effect.

**2.1. The general equations**

140 Hydrodynamics - Concepts and Experiments

**2. Evolution of half-space free boundary perturbation**

the renormalized pressure *Π* =*P* + *ρgz* (where *P* is the true pressure) as

ì

ï

surface, in the approximation linear by *h* , have the form

density and surface tension coefficient of the liquid.

2 2 zz z z z

<sup>1</sup> VV ,

bb

ï ¶ +¶ = ¶ P <sup>m</sup> <sup>ï</sup>

ïï¶ +¶ = ¶ P <sup>í</sup> <sup>m</sup> <sup>ï</sup> ï ¶ +¶ =

<sup>1</sup> VV ,

a bb a a

V V 0,

b b bb

The indices *α*, *β* can be 1 or 2; summation over repeated indices is assumed. The boundary conditions which express the continuity of the normal and tangential stresses on the free

0.

z z z z z 0 z 0 2 V h gh, V V 0. bb <sup>a</sup> <sup>a</sup> <sup>=</sup> <sup>=</sup> -P + m¶ = s¶ - r ¶ + ¶ = (2)

2 2 zz

> z z 2 2 zz

ï ¶ P+¶ P= î

( ) ( ) <sup>2</sup>

z

V

æ ö ç ÷ <sup>=</sup> ç ÷ ç ÷ è ø P

z

V 0. <sup>a</sup> ®-¥

Here *Vz* and *Vα* are the components of velocity, ∂*z* and ∂*α* are the derivatives along the coordinates *z* and *xα* accordingly, *g* is the acceleration due to gravity; *µ*, *ρ* and *σ* are the viscosity,

The following conditions in case of infinite half-space must also be met:

$$\begin{cases} \hat{\sigma}\_{\rm zx}^{2} \mathbf{V}\_{\rm x\mathbf{k}} - \mathbf{k}^{2} \mathbf{V}\_{\rm x\mathbf{k}} = \frac{1}{\mu} \hat{\sigma}\_{x} \Pi\_{\mathbf{k}} \\\\ \hat{\sigma}\_{\rm zx}^{2} \mathbf{V}\_{\rm o\mathbf{k}} - \mathbf{k}^{2} \mathbf{V}\_{\rm o\mathbf{k}} = \frac{\mathbf{i} \mathbf{k}\_{\alpha}}{\mu} \Pi\_{\mathbf{k}} \\\\ \hat{\sigma}\_{x} \mathbf{V}\_{\rm x\mathbf{k}} + \mathbf{i} \mathbf{k}\_{\beta} \mathbf{V}\_{\beta\mathbf{k}} = 0 \\\\ \hat{\sigma}\_{\rm zx}^{2} \Pi\_{\mathbf{k}} - \mathbf{k}^{2} \Pi\_{\mathbf{k}} = 0 \end{cases} \tag{5}$$

$$\left. \left( -\boldsymbol{\Pi}\_{\mathbf{k}} + 2\mu \boldsymbol{\hat{o}}\_{x} \mathbf{V}\_{x\mathbf{k}} \right) \right|\_{\mathbf{z}=0} = -\left( \sigma \mathbf{k}^{2} + \rho \mathbf{g} \right) \mathbf{h}\_{\mathbf{k}}, \quad \left( \mathrm{ik}\_{\alpha} \mathbf{V}\_{x\mathbf{k}} + \boldsymbol{\hat{o}}\_{x} \mathbf{V}\_{\alpha \mathbf{k}} \right) \Big|\_{\mathbf{z}=0} = \mathbf{0}, \quad \begin{pmatrix} \mathbf{V}\_{x\mathbf{k}} \\ \mathbf{V}\_{\alpha \mathbf{k}} \\ \boldsymbol{\Pi}\_{\mathbf{k}} \end{pmatrix}\_{\mathbf{z} \to -\eta} = \mathbf{0}, \tag{6}$$

$$\left.\hat{\sigma}\_{\mathbf{t}}\mathbf{h}\_{\mathbf{k}} = \mathbf{V}\_{\mathbf{z}\mathbf{k}}\right|\_{\mathbf{x}\_{\perp}=0} + \mathbf{W}\_{\mathbf{k}}.\tag{7}$$

Integrating (5) with boundary conditions (6), we obtain

$$
\begin{pmatrix}
\mathbf{V\_{zk}} \\
\mathbf{V\_{\alpha k}} \\
\mathbf{I}\mathbf{I\_{k}}
\end{pmatrix} = \begin{pmatrix}
\mathbf{k}\mathbf{z} - \mathbf{l} \\
\mathbf{i}\mathbf{k}\_{\alpha}\mathbf{z} \\
2\mu\mathbf{k}
\end{pmatrix} \frac{\sigma\mathbf{k}^{2} + \rho\mathbf{g}}{2\mu\mathbf{k}} \mathbf{h}\_{\mathbf{k}} \mathbf{e}^{\mathbf{k}\mathbf{z}}.\tag{8}
$$

It is apparent that the Fourier component *<sup>h</sup> <sup>k</sup>* has an effect in the layer of depth *<sup>k</sup>* <sup>−</sup><sup>1</sup> . Based on (7), (8), the equation of the free surface motion can be written as

$$\left(\hat{\sigma}\_{\text{t}}\mathbf{h}\_{\text{k}} + \frac{1}{2\mu\mathbf{k}}(\sigma\mathbf{k}^{2} + \rho\mathbf{g})\mathbf{h}\_{\text{k}} = \mathbf{W}\_{\text{k}}\tag{9}$$

#### **2.2. Relaxation of the initial perturbation of the liquid surface, exact solution**

Consider the relaxation of the initial perturbation of the surface due to capillary and gravitation forces in the absence of mass flux at the boundary (*W<sup>k</sup>* =0). In this case, the solution of equation (9) with initial condition *<sup>h</sup> <sup>k</sup>* <sup>|</sup> *<sup>t</sup>*=0 <sup>=</sup>*<sup>h</sup> <sup>k</sup>* 0 becomes

$$\mathbf{h\_{k}} = \mathbf{h\_{k}^{0}e^{-t f^{\Gamma}(k)}},\tag{10}$$

where *<sup>T</sup>* (*k*)= <sup>2</sup>*µk <sup>σ</sup><sup>k</sup>* <sup>2</sup> <sup>+</sup> *<sup>ρ</sup><sup>g</sup>* is the characteristic relaxation time of the corresponding harmonics. For *k* →0 and *k* →*∞* the relaxation time tends to zero assuming non-zero gravity and surface tension, respectively. At *<sup>k</sup>* <sup>=</sup>*k*\*=(*ρ<sup>g</sup>* / *<sup>σ</sup>*)1/2 the relaxation time has its maximal value *T*\*=*<sup>µ</sup>* / (*ρgσ*)1/2 . Small relaxation time for high-frequency harmonics means that the surface tension forces lead to a rapid smoothing of sharp irregularities. On the other hand *<sup>T</sup>* <sup>|</sup> *<sup>k</sup>* =0 =0, therefore *<sup>h</sup> <sup>k</sup>* <sup>|</sup> *<sup>k</sup>* =0 =0 at *t* >0. Noting that

$$\left. \mathbf{h\_{k}} \right|\_{\mathbf{k} = 0} = \left\{ \mathbf{h(x,t)} \mathbf{d^{2}x} \right\}$$

constitutes the "volume of perturbation". According to (10) this value becomes equal to zero due to the gravity force. Immediacy of this process is determined using a quasi-stationary Stokes approximation. Detailed discussion of this topic see in [14].

Switching to the case of arbitrary wave numbers, we introduce the dimensionless variables

$$\mathbf{\dot{\tau}} = \mathbf{t} \mathbf{\dot{\tau}} \mathbf{T}\_\*, \quad \Theta = \mathbf{T} \mathbf{\dot{\tau}} \mathbf{T}\_\*, \quad \mathbf{\dot{\xi}} = \mathbf{k}\_\* \mathbf{x}, \quad \mathbf{\kappa} = \mathbf{k} \,/ \,\, \mathbf{k}\_{\*-} $$

Then (10) takes the form

$$\mathbf{h}\_{\kappa} = \mathbf{h}\_{\kappa}^{0} \mathbf{e}^{-\mathfrak{q}\hbar \mathbf{e}(\kappa)}, \quad \Theta\left(\kappa\right) = \frac{2\kappa}{\kappa^{2} + \mathbf{l}}, \quad \mathbf{h}\_{\kappa} = \int \mathbf{h}\left(\mathfrak{k}\_{\nu}\mathbf{r}\right) \mathbf{e}^{-i\kappa \mathbf{\tilde{x}}}.\tag{11}$$

Applying inverse Fourier transform to (11) we have

$$\mathrm{ch}\left(\boldsymbol{\xi},\boldsymbol{\tau}\right) = \left\{ \mathbf{g}\left(\boldsymbol{\xi} - \boldsymbol{\xi}',\boldsymbol{\tau}\right) \mathbf{h}^{0}\left(\boldsymbol{\xi}'\right) \mathbf{d}^{2}\boldsymbol{\xi}', \quad \mathbf{g}\left(\boldsymbol{\xi},\boldsymbol{\tau}\right) = \left\{ \mathbf{e}^{i\mathbf{x}\boldsymbol{\xi} - \boldsymbol{\eta}'\boldsymbol{\theta}\left(\boldsymbol{\kappa}\right)} \frac{\mathbf{d}^{2}\boldsymbol{\kappa}}{\left(2\pi\right)^{2}},\tag{12}$$

The expression (12) allows calculation of the surface perturbation in any time for an arbitrary initial perturbation.

#### **2.3. Qualitative description of the perturbation relaxation**

**2.2. Relaxation of the initial perturbation of the liquid surface, exact solution**

becomes

<sup>0</sup> tTk( ) h he , -

0

Stokes approximation. Detailed discussion of this topic see in [14].

Applying inverse Fourier transform to (11) we have

(9) with initial condition *<sup>h</sup> <sup>k</sup>* <sup>|</sup> *<sup>t</sup>*=0 <sup>=</sup>*<sup>h</sup> <sup>k</sup>*

142 Hydrodynamics - Concepts and Experiments

respectively. At *<sup>k</sup>* <sup>=</sup>*k*\*=(*ρ<sup>g</sup>* / *<sup>σ</sup>*)1/2

where *<sup>T</sup>* (*k*)= <sup>2</sup>*µk*

*t* >0. Noting that

Then (10) takes the form

initial perturbation.

Consider the relaxation of the initial perturbation of the surface due to capillary and gravitation forces in the absence of mass flux at the boundary (*W<sup>k</sup>* =0). In this case, the solution of equation

*k* →0 and *k* →*∞* the relaxation time tends to zero assuming non-zero gravity and surface tension,

relaxation time for high-frequency harmonics means that the surface tension forces lead to a rapid smoothing of sharp irregularities. On the other hand *<sup>T</sup>* <sup>|</sup> *<sup>k</sup>* =0 =0, therefore *<sup>h</sup> <sup>k</sup>* <sup>|</sup> *<sup>k</sup>* =0 =0 at

constitutes the "volume of perturbation". According to (10) this value becomes equal to zero due to the gravity force. Immediacy of this process is determined using a quasi-stationary

Switching to the case of arbitrary wave numbers, we introduce the dimensionless variables

\* \*\* \* t= q= = = t T, T T, k , /k . **ξ x κk**

( ) ( ) ( ) 0 i 2 <sup>2</sup> h he , ,h h ,e . 1 -t q k <sup>k</sup> - = qk = = t <sup>ò</sup> k +

<sup>2</sup> 0 2 <sup>i</sup>

( ) ( ) ( ) ( ) ( )

<sup>d</sup> h, g ,h d ,g, e .


The expression (12) allows calculation of the surface perturbation in any time for an arbitrary

( ) <sup>2</sup> k 0 h h ,t d x <sup>=</sup> **<sup>k</sup>** = ò **x**

*<sup>σ</sup><sup>k</sup>* <sup>2</sup> <sup>+</sup> *<sup>ρ</sup><sup>g</sup>* is the characteristic relaxation time of the corresponding harmonics. For

the relaxation time has its maximal value *T*\*=*<sup>µ</sup>* / (*ρgσ*)1/2

**κξ**

( )

2

2

**κ κ <sup>κ</sup> ξ** (11)

**κξ ξ ξξ ξ ξ** (12)

**k k** = (10)

. Small

For a qualitative description of the relaxation of the initial perturbation is convenient to use the original expression (11). We have

$$\mathbf{h}\left(\xi,\tau\right) = \left\{ \mathbf{e}^{i\mathbf{k}\cdot\mathbf{\tilde{x}}} \mathbf{h}\_{\mathbf{k}}^{0} \mathbf{e}^{-\tau f^{0}\left(\kappa\right)} \frac{\mathbf{d}^{2}\kappa}{\left(2\pi\right)^{2}}, \quad \Theta\left(\kappa\right) = \frac{2\kappa}{\mathbf{k}^{2} + \mathbf{l}}.\tag{13}$$

From (13) it is clear that the configuration at the moment *τ* is determined by the Fourier transform of the initial perturbation *h <sup>κ</sup>* 0 and the multiplier *e* <sup>−</sup>*τ*/*θ*(*κ*) . The maximum value of this factor is equal to *e* <sup>−</sup>*<sup>τ</sup>* , and corresponds to *τ* =1. Let *κ*1,2 (*τ*) be the values of *κ* at which the multiplier takes a value in *e* times less than the maximum-see Fig. 1.

**Figure 1.** Graph of the function *e* <sup>−</sup>*τ*/*θ*(*κ*) at *τ* =1 and its approximation with a step function.

By solving the equation *<sup>τ</sup> <sup>θ</sup>*(*κ*) <sup>=</sup>*<sup>τ</sup>* <sup>+</sup> 1, we obtain

$$\kappa\_{1,2} = \frac{1}{\pi} \Big(\pi + \mathbf{l} \pm \sqrt{2\pi + \mathbf{l}}\Big). \tag{14}$$

$$\text{At } \pi \to 0, \quad \kappa\_1 = \frac{\pi}{2}, \quad \kappa\_2 = \frac{2}{\pi}. \text{ On the other hand, at } \pi \to \infty, \quad \kappa\_{1,2} = 1 \mp \sqrt{\frac{2}{\pi}}.$$

To estimate perturbation evolution we replace function *e* <sup>−</sup>*τ*/*θ*(*κ*) with the step function *f* (*τ*, *κ*) which is equal to zero for *κ* <*κ*<sup>1</sup> (*τ*) and for *κ* >*κ*<sup>2</sup> (*τ*). In the interval *κ*<sup>1</sup> (*τ*)≤*κ* ≤*κ*<sup>2</sup> (*τ*) this function is equal to *e* <sup>−</sup>*<sup>τ</sup>* -see Fig.1:

$$\mathbf{f}\left(\kappa,\tau\right) = \mathbf{e}^{-\tau}\boldsymbol{\Theta}\left(\kappa - \kappa\_1\left(\tau\right)\right)\boldsymbol{\Theta}\left(\kappa\_2\left(\tau\right) - \kappa\right),\tag{15}$$

where *ϑ*(*x*) is the Heaviside step function.

Returning to the expression (13) we see that the evolution of perturbations can be considered as a result the initial perturbation "cut-off" in Fourier space. The system under consideration behaves as filter which "cuts" in the spectrum of the initial perturbation wave numbers the annular region with characteristic dimensions *κ*<sup>1</sup> (*τ*)≤ |*κ* | ≤*κ*<sup>2</sup> (*τ*). At the same time the amplitude of the Fourier transform of the perturbation decreases as *e* <sup>−</sup>*<sup>τ</sup>* . It can be seen that the width of the annular region decreases with increasing time, rapidly at the beginning-mainly due to the reduction of the outer radius of the ring *κ*2, then slowly-due to symmetric increasing *κ*<sup>1</sup> and decreasing *κ*2. At large times, the ring is localized in the neighborhood of a circle of radius *κ* =1 and has a characteristic width *Δκ* =2 2 *<sup>τ</sup>* . Such filtering process determines the qualitative features of the relaxation of the initial perturbation.

Spatial filtering of the Fourier transform of the initial perturbation in the region of small wave numbers is equivalent to subtracting its part localized in a circle of radius *κ*1. It leads to the subtraction of additional perturbation localized in a circle with radius *ξ*<sup>1</sup> <sup>=</sup> <sup>2</sup>*<sup>π</sup> k*1 from the initial perturbation. When *<sup>τ</sup>* <sup>→</sup>0, *<sup>ξ</sup>*<sup>1</sup> <sup>=</sup> <sup>4</sup>*<sup>π</sup> <sup>τ</sup>* <sup>→</sup>*∞*. This long-range effect of the initial perturbation is determined by the force of gravity. This is clearly seen in the transition to dimensional movements: the radius of the long-range perturbations is *r*<sup>1</sup> = *ξ*1 *k*\* <sup>=</sup> <sup>4</sup>*π<sup>µ</sup> <sup>ρ</sup>gt* . Figuratively speaking, the perturbation is drowning in the ambient fluid, distorting the shape of the surface at long distance.

Localization of the Fourier transform of the perturbation in a circle of radius *κ*2 leads to an expansion of its external borders, in the coordinate representation, to a circle of radius *<sup>ξ</sup>*<sup>2</sup> <sup>=</sup> <sup>2</sup>*<sup>π</sup> κ*2 . In dimensional variables *r*<sup>2</sup> = *ξ*2 *k*\* <sup>=</sup> *πσ<sup>t</sup> <sup>µ</sup>* , i.e. this effect is due to the action of surface tension forces.

To evaluate the value of the perturbation at the origin of coordinates we use of the fact that

$$h\left(0,\ \tau\right) = \int h\_{\mathbf{n}}(\tau) \frac{d^2\kappa}{(2\,\pi)^2}.$$

Then

$$h\left(0,\tau\right) = \int\_{\kappa\_1^{(\epsilon)} \le \left|\sigma\right| \le \kappa\_2^{(\epsilon)}} h\_{\mathfrak{A}}^{\frac{\epsilon}{2}} \frac{d^2 \kappa}{(2\pi)^2}$$

The above statements are based on the approximate replacement of real *e* <sup>−</sup>*τ*/*θ*(*κ*) functions with the step function (15). Therefore, we can only expect an approximate correspondence between the description, based on these statements, and the actual process of the perturbation relaxa‐ tion.

Let us consider the relaxation of perturbations with different spatial scale and shape. To illustrate the qualitative dependencies we will use graphics produced on the basis of the exact solution (12).

#### **2.4. Relaxation long-wave perturbations by gravity**

f, e ( ) ( 1 2 ( )) ( ( ) ), -t k t = J k-k t J k t -k (15)

(*τ*). At the same time the amplitude

*k*1

*<sup>ρ</sup>gt* . Figuratively speaking,

from the initial

*<sup>τ</sup>* . Such filtering process determines the qualitative

*<sup>τ</sup>* <sup>→</sup>*∞*. This long-range effect of the initial perturbation is

*ξ*1 *k*\* <sup>=</sup> <sup>4</sup>*π<sup>µ</sup>*

*<sup>µ</sup>* , i.e. this effect is due to the action of surface tension

. It can be seen that the width of

Returning to the expression (13) we see that the evolution of perturbations can be considered as a result the initial perturbation "cut-off" in Fourier space. The system under consideration behaves as filter which "cuts" in the spectrum of the initial perturbation wave numbers the

the annular region decreases with increasing time, rapidly at the beginning-mainly due to the reduction of the outer radius of the ring *κ*2, then slowly-due to symmetric increasing *κ*<sup>1</sup> and decreasing *κ*2. At large times, the ring is localized in the neighborhood of a circle of radius *κ* =1

Spatial filtering of the Fourier transform of the initial perturbation in the region of small wave numbers is equivalent to subtracting its part localized in a circle of radius *κ*1. It leads to the

determined by the force of gravity. This is clearly seen in the transition to dimensional

the perturbation is drowning in the ambient fluid, distorting the shape of the surface at long

Localization of the Fourier transform of the perturbation in a circle of radius *κ*2 leads to an expansion of its external borders, in the coordinate representation, to a circle of radius

To evaluate the value of the perturbation at the origin of coordinates we use of the fact that

(τ) d2 κ

∫ κ1 (τ)≤|κ|≤κ2 (τ) h κ

(2π)2

> 0 d2 κ

> > (2π)2

2

subtraction of additional perturbation localized in a circle with radius *ξ*<sup>1</sup> <sup>=</sup> <sup>2</sup>*<sup>π</sup>*

*ξ*2 *k*\* <sup>=</sup> *πσ<sup>t</sup>*

> h (0, τ)=∫h κ

h (0, τ)= (*τ*)≤ |*κ* | ≤*κ*<sup>2</sup>

where *ϑ*(*x*) is the Heaviside step function.

144 Hydrodynamics - Concepts and Experiments

annular region with characteristic dimensions *κ*<sup>1</sup>

features of the relaxation of the initial perturbation.

and has a characteristic width *Δκ* =2

perturbation. When *<sup>τ</sup>* <sup>→</sup>0, *<sup>ξ</sup>*<sup>1</sup> <sup>=</sup> <sup>4</sup>*<sup>π</sup>*

. In dimensional variables *r*<sup>2</sup> =

distance.

*<sup>ξ</sup>*<sup>2</sup> <sup>=</sup> <sup>2</sup>*<sup>π</sup> κ*2

forces.

Then

of the Fourier transform of the perturbation decreases as *e* <sup>−</sup>*<sup>τ</sup>*

movements: the radius of the long-range perturbations is *r*<sup>1</sup> =

If the characteristic dimensionless scale of the region occupied by the perturbation *a* > >1, then the Fourier transform of this disturbance is concentrated in the region with a typical size 2*π <sup>a</sup>* <sup>&</sup>lt; <1. At the initial stage the relaxation of such perturbation is mainly determined by filtering out small wave numbers of its Fourier representation, i.e., by gravity. In this case, as shown above, to the initial perturbation is added, taken with the opposite sign, the Fourier inverse transform of the initial perturbation, lying in a circle of radius *κ*<sup>1</sup> (*τ*). In the coordinate repre‐ sentation this additive is concentrated in a circle of radius *ξ*<sup>1</sup> (*τ*)=2*π* / *κ*<sup>1</sup> (*τ*) and has an amplitude

$$\Delta h(x) = \int\_{|\mathfrak{a}| \le \kappa\_1(x)} h\_{\mathfrak{a}}^{\;0} \frac{d^2 \kappa}{(2\pi)^2}$$

At the initial stage of relaxation process *κ*<sup>1</sup> < <2*π* / *a*, respectively, *ξ*<sup>1</sup> > >*a*,

$$\Delta h(\pi) = \frac{1}{4\pi} \kappa\_1^2(\pi) \left. h\_{\mathfrak{A}}^{\;0} \right|\_{\mathfrak{A}^\equiv 0} = \frac{\pi^2}{16\pi} \int h^{\;0}(\mathfrak{F}) \, d^2 \xi$$

Through this long-wavelength addition, the initial perturbation goes down (in case of positive volume) on *Δh* (*τ*), practically without distortion of the shape. Outside the region of the initial perturbation arises an additional perturbation with the amplitude *Δh* (*τ*) and attenuation radius *ξ*1 (*τ*). At large times, the shape of additional perturbation is somewhat distorted, its radius decreases and becomes comparable to the size of the initial perturbation.

By the time *τ*<sup>1</sup> =4*π* / *a*, when *κ*<sup>1</sup> =2*π* / *a*, the process of relaxation of the perturbation is mainly completed. To illustrate these effects on the Fig.2 presented the relaxation of axisymmetric long-wave perturbation shape. On the Fig.3 we can see the increasing of this perturbation amplitude.

If the characteristic dimensions of the initial perturbation in two mutually perpendicular directions are very different, the relaxation process is more complicated. Note that for the perturbation, elongated along an axis *ξ*1 (*a*<sup>1</sup> > >*a*<sup>2</sup> > >1), the region of localization of the Fourier transform has the inverse aspect ratio 2*π* / *a*<sup>1</sup> < <2*π* / *a*2. Removal from the spectrum of the initial perturbation the circle with radius *κ*<sup>1</sup> (*τ*)=*τ* / 2 leads initially to the formation of a symmetric long-wave perturbation of the radius 4*<sup>π</sup>* / *<sup>τ</sup>* and the amplitude *Δ<sup>h</sup>* <sup>=</sup> *<sup>τ</sup>* <sup>2</sup> <sup>16</sup>*<sup>π</sup> <sup>∫</sup><sup>h</sup>* 0(*ξ*)*<sup>d</sup>* <sup>2</sup> *ξ*. At this stage, the difference of the transverse scales the perturbation does not affect practically the nature of its evolution. Perturbation settles as a whole, and its amplitude decreases according to *h* (0, *τ*)=*h* 0(0)−*Δh* (*τ*).

**Figure 2.** The relaxation of axisymmetric long-wave perturbations. The initial shape of perturbation is *h* <sup>0</sup> = *Ae* <sup>−</sup>*<sup>ξ</sup>* <sup>2</sup> /*a* <sup>2</sup> , *a* =20.

**Figure 3.** Increasing of amplitude of long-wave axisymmetric perturbation. The initial shape of perturbation is *h* <sup>0</sup> = *Ae* <sup>−</sup>*<sup>ξ</sup>* <sup>2</sup> /*a* <sup>2</sup> , *a* =20.

From the time moment *τ* =4*π* / *a*1, when *k*<sup>1</sup> (*τ*)=4*π* / *a*1, the area of influence of long-range pertur‐ bations loses circular symmetry. In the direction of the axis *ξ*<sup>1</sup> the scale of the area *ξ*<sup>1</sup> \* becomes, and subsequently maintained, equal to *a*1, and along the axis *ξ*2 this area continues to shrink as *ξ*<sup>2</sup> \* =4*π* / *τ*. The amplitude of the perturbation begins to decrease linearly *<sup>h</sup>* (0, *<sup>τ</sup>*)=*<sup>h</sup>* 0(0)−*Δ<sup>h</sup>* (*τ*), *<sup>Δ</sup><sup>h</sup>* <sup>=</sup> *<sup>τ</sup>* 4*πL* <sup>1</sup> *∫h* 0(*ξ*)*d* <sup>2</sup> *ξ*.

long-wave perturbation of the radius 4*<sup>π</sup>* / *<sup>τ</sup>* and the amplitude *Δ<sup>h</sup>* <sup>=</sup> *<sup>τ</sup>* <sup>2</sup>

*h* (0, *τ*)=*h* 0(0)−*Δh* (*τ*).

146 Hydrodynamics - Concepts and Experiments

*h* <sup>0</sup> = *Ae* <sup>−</sup>*<sup>ξ</sup>* <sup>2</sup>

*h* <sup>0</sup> = *Ae* <sup>−</sup>*<sup>ξ</sup>* <sup>2</sup>

/*a* <sup>2</sup>

, *a* =20.

/*a* <sup>2</sup>

, *a* =20.

the difference of the transverse scales the perturbation does not affect practically the nature of its evolution. Perturbation settles as a whole, and its amplitude decreases according to

**Figure 2.** The relaxation of axisymmetric long-wave perturbations. The initial shape of perturbation is

**Figure 3.** Increasing of amplitude of long-wave axisymmetric perturbation. The initial shape of perturbation is

<sup>16</sup>*<sup>π</sup> <sup>∫</sup><sup>h</sup>* 0(*ξ*)*<sup>d</sup>* <sup>2</sup>

*ξ*. At this stage,

At the moment *τ* =4*π* / *a*2, when *κ*<sup>1</sup> (*τ*)=2*π* / *a*2 the perturbation is practically relaxed. Figure 4 shows the characteristic configuration of the long-wavelength perturbation, the extent of which the two axes are substantially different.

**Figure 4.** The configuration of the long-wave non-axisymmetric perturbations in the time *τ* =0.04. The initial shape of perturbation is *h* <sup>0</sup> = *A*exp( − *ξ*1 2 *a*1 <sup>2</sup> − *ξ*2 2 *a*2 <sup>2</sup> ), *a*<sup>1</sup> =100, *a*<sup>2</sup> =20.

#### **2.5. Relaxation of short-wave perturbations due to capillary forces**

The evolution of perturbation with the scale *a* < <1 is different. The Fourier transform of this disturbance is concentrated in the region with a typical size of 2*π* / *a* > >1. In this case occurs the evolution of the perturbation in the initial stage due to filtration of its Fourier transform through small wave numbers. However, the contribution of this region to the total amplitude of the perturbation is small, and therefore, the effects of subsidence of the perturbation and the formation of the region of the far influence, in this case, is irrelevant.

Relaxation of short-wave perturbations is mainly determined by the action of surface tension forces, i.e. filtration its Fourier transform by the large wave numbers. The relaxation process begins at the moment *τ* =*a* / *π*, when *κ*<sup>2</sup> =2*π* / *a*. The size of the region occupied by the perturbation begins to grow as *ξ*(*τ*)=2*π* / *κ*<sup>2</sup> =*πτ*, and its amplitude decreases at the initial stage, as 1 *πτ* <sup>2</sup> *<sup>∫</sup><sup>h</sup>* 0(*ξ*)*<sup>d</sup>* <sup>2</sup> *ξ*. Characteristic decay time of the perturbation is of order unity. Fig. 5 shows the relaxation of short-wave axisymmetric perturbation.

**Figure 5.** Relaxation of short-wave axisymmetric perturbation. The initial shape of perturbation *h* <sup>0</sup> = *Ae* <sup>−</sup>*<sup>ξ</sup>* <sup>2</sup> /*a* <sup>2</sup> , *a* =0.1.

Finally let's describe the relaxation of small-scale perturbations with significantly different scales along the axes *ξ*1 and *ξ*2 : *a*<sup>1</sup> < <*a*<sup>2</sup> < <1. The evolution of this perturbation begins at the moment *τ* =*a*<sup>1</sup> / *π*. The scale begins to increase like *πτ* in the *ξ*1 axis direction and in *ξ*<sup>2</sup> axis direction it varies slightly. The amplitude of the perturbation decreases according to *<sup>h</sup>* (0, *<sup>τ</sup>*)= <sup>2</sup> *<sup>π</sup>a*2*<sup>τ</sup> <sup>∫</sup><sup>h</sup>* 0(*ξ*)*<sup>d</sup>* <sup>2</sup> *ξ*.

After the moment *τ* =*a*<sup>2</sup> / *π* the magnitudes of the perturbation in the direction of the two axes become comparable. From this moment the perturbation becomes axially symmetric, and its radius increases as *πτ*. Amplitude of perturbation thus decreases faster than during the first stage: *<sup>h</sup>* (0, *<sup>τ</sup>*)= <sup>1</sup> *πτ* <sup>2</sup> *<sup>∫</sup><sup>h</sup>* 0(*ξ*)*<sup>d</sup>* <sup>2</sup> *ξ*.

The relaxation process ends substantially at the moment *τ* =1. Fig.5 shows relaxation of shortwave non-axisymmetric perturbations. One can see that over time the perturbation becomes axisymmetric.

Evolution of Small Perturbations of the Free Surface of Viscous Fluid in the Stokes Approximation http://dx.doi.org/10.5772/59226 149

Figure 6. Configurations of the short‐wave non‐axisymmetric perturbation. The initial shape of perturbation is 2 2 <sup>0</sup> 1 2 2 2 1 2 1 2 h Aexp , a 0.02, a 0.2. a a **Figure 6.** Configurations of the short-wave non-axisymmetric perturbation. The initial shape of perturbation is *h* <sup>0</sup> = *A*exp( − *ξ*1 2 *a*1 <sup>2</sup> − *ξ*2 2 *a*2 <sup>2</sup> ), *a*<sup>1</sup> =0.02, *a*<sup>2</sup> =0.2.

In this section we study surface perturbations of a flat liquid layer with depth H resting on a solid bottom. We use system (1) and conditions on the free boundary (2) to calculate perturbation evolution due to gravity and surface tension.

Again applying the Fourier transform in the longitudinal coordinates and solving the resulting system of ordinary

**k k**

<sup>0</sup> z z0

#### The conditions at infinity (3) are replaced by the no‐slip condition along the solid surface <sup>z</sup> z H z H V 0, V 0. (16) **3. Finite layer over solid bottom and over a half-space**

**3.1 Finite layer over solid bottom**

**3. Finite layer over solid bottom and over a half‐space**

#### differential equations, we obtain **3.1. Finite layer over solid bottom**

Relaxation of short-wave perturbations is mainly determined by the action of surface tension forces, i.e. filtration its Fourier transform by the large wave numbers. The relaxation process begins at the moment *τ* =*a* / *π*, when *κ*<sup>2</sup> =2*π* / *a*. The size of the region occupied by the perturbation begins to grow as *ξ*(*τ*)=2*π* / *κ*<sup>2</sup> =*πτ*, and its amplitude decreases at the initial stage, as

**Figure 5.** Relaxation of short-wave axisymmetric perturbation. The initial shape of perturbation

Finally let's describe the relaxation of small-scale perturbations with significantly different scales along the axes *ξ*1 and *ξ*2 : *a*<sup>1</sup> < <*a*<sup>2</sup> < <1. The evolution of this perturbation begins at the moment *τ* =*a*<sup>1</sup> / *π*. The scale begins to increase like *πτ* in the *ξ*1 axis direction and in *ξ*<sup>2</sup> axis direction it varies slightly. The amplitude of the perturbation decreases according to

After the moment *τ* =*a*<sup>2</sup> / *π* the magnitudes of the perturbation in the direction of the two axes become comparable. From this moment the perturbation becomes axially symmetric, and its radius increases as *πτ*. Amplitude of perturbation thus decreases faster than during the first

The relaxation process ends substantially at the moment *τ* =1. Fig.5 shows relaxation of shortwave non-axisymmetric perturbations. One can see that over time the perturbation becomes

relaxation of short-wave axisymmetric perturbation.

*ξ*. Characteristic decay time of the perturbation is of order unity. Fig. 5 shows the

1 *πτ* <sup>2</sup> *<sup>∫</sup><sup>h</sup>* 0(*ξ*)*<sup>d</sup>* <sup>2</sup>

148 Hydrodynamics - Concepts and Experiments

*h* <sup>0</sup> = *Ae* <sup>−</sup>*<sup>ξ</sup>* <sup>2</sup>

*<sup>h</sup>* (0, *<sup>τ</sup>*)= <sup>2</sup>

stage: *<sup>h</sup>* (0, *<sup>τ</sup>*)= <sup>1</sup>

axisymmetric.

/*a* <sup>2</sup>

, *a* =0.1.

*<sup>π</sup>a*2*<sup>τ</sup> <sup>∫</sup><sup>h</sup>* 0(*ξ*)*<sup>d</sup>* <sup>2</sup>

*ξ*.

*πτ* <sup>2</sup> *<sup>∫</sup><sup>h</sup>* 0(*ξ*)*<sup>d</sup>* <sup>2</sup>

*ξ*.

 0 H <sup>H</sup> sh k z H z H V V , 2 sh kH 2 sh k z H sh kz , sh kH sh kH **k k k** (17) where In this section we study surface perturbations of a flat liquid layer with depth H resting on a solid bottom. We use system (1) and conditions on the free boundary (2) to calculate pertur‐ bation evolution due to gravity and surface tension. The conditions at infinity (3) are replaced by the no-slip condition along the solid surface

$$\left. \mathbf{V}\_{\mathbf{z}} \right|\_{\mathbf{z} = -\mathbf{H}} = \mathbf{0}, \quad \left. \mathbf{V}\_{\alpha} \right|\_{\mathbf{z} = -\mathbf{H}} = \mathbf{0}. \tag{16}$$

**k k**

Again applying the Fourier transform in the longitudinal coordinates and solving the resulting system of ordinary differential equations, we obtain

$$\begin{split} \mathbf{V\_{xk}} &= \left(\mathbf{V\_{xk0}} - \frac{\mathrm{H}\Pi\_{\mathbf{k}0}}{2\mu}\right) \frac{\mathrm{sh}\left(\mathrm{k}\left(\mathrm{z} + \mathrm{H}\right)\right)}{\mathrm{sh}\left(\mathrm{k}\mathrm{H}\right)} + \frac{\left(\mathrm{z} + \mathrm{H}\right)\Pi\_{\mathbf{k}}}{2\mu}, \\ \boldsymbol{\Pi}\_{\mathbf{k}} &= \boldsymbol{\Pi}\_{\mathbf{k}0} \frac{\mathrm{sh}\left(\mathrm{k}\left(\mathrm{z} + \mathrm{H}\right)\right)}{\mathrm{sh}\left(\mathrm{k}\mathrm{H}\right)} - \boldsymbol{\Pi}\_{\mathbf{k}\mathrm{H}} \frac{\mathrm{sh}\left(\mathrm{k}\mathrm{z}\right)}{\mathrm{sh}\left(\mathrm{k}\mathrm{H}\right)}, \end{split} \tag{17}$$

where

$$
\begin{pmatrix}
\mathbf{V\_{k0}} \\
\Pi\_{\mathbf{k}0} \\
\Pi\_{\mathbf{k}\mathbf{H}}
\end{pmatrix} = \begin{pmatrix}
2\mathbf{k}\mathbf{H} - \text{sh}\left(2\mathbf{k}\mathbf{h}\right) \\
2\mu\mathbf{k}\left(\text{ch}\left(2\mathbf{k}\mathbf{H}\right) + \mathbf{l}\right) \\
4\mu\mathbf{k}\left(\text{k}\mathbf{H}\text{sh}\left(\mathbf{k}\mathbf{H}\right) + \text{ch}\left(\mathbf{k}\mathbf{H}\right)\right)
\end{pmatrix} \frac{\sigma\mathbf{k}^2 + \rho\mathbf{g}}{2\mu\mathbf{k}\left(\text{ch}\left(2\mathbf{k}\mathbf{H}\right) + \mathbf{l} + 2\left(\mathbf{k}\mathbf{H}\right)^2\right)}\mathbf{h}\_{\mathbf{k}}.\tag{18}
$$

The evolution equation of the free surface and the initial condition has the form

$$\left\|\mathbf{\tilde{c}}\_{\text{i}}\mathbf{h}\_{\text{k}} + \frac{2\mathbf{k}\mathbf{H} - \text{sh}\left(2\mathbf{k}\mathbf{h}\right)}{\text{ch}\left(2\mathbf{k}\mathbf{H}\right) + \text{l} + 2\left(\mathbf{k}\mathbf{H}\right)^{2}} \frac{\sigma\mathbf{k}^{2} + \rho\mathbf{g}}{2\mu\mathbf{k}}\mathbf{h}\_{\text{k}} = \mathbf{0}, \quad \left\|\mathbf{h}\_{\text{k}}\right\|\_{\text{t}=0} = \mathbf{h}\_{\text{k}}^{0}.\tag{19}$$

The solution of equation (19) is similar to (10)

$$\mathbf{h\_{k}} = \mathbf{h\_{k}^{0}e^{-\mathbf{q}^{\mathrm{tr}}\mathbf{\hat{r}}(k)}},\tag{20}$$

where *<sup>T</sup>* (*k*)= *ch* (2*kH* ) <sup>+</sup> <sup>1</sup> <sup>+</sup> 2(*kH* )<sup>2</sup> *sh* (2*kh* )−2*kH* 2*µk <sup>σ</sup><sup>k</sup>* <sup>2</sup> <sup>+</sup> *<sup>ρ</sup><sup>g</sup>* is the characteristic relaxation time of corresponding harmonic.

At *kH* <sup>&</sup>gt; >1, *<sup>T</sup>* (*k*, *<sup>H</sup>* )= <sup>2</sup>*µk <sup>σ</sup><sup>k</sup>* <sup>2</sup> <sup>+</sup> *<sup>ρ</sup><sup>g</sup>* coincides with (10), i.e. with the case of a fluid of infinite depth. Therefore, for *k* →*∞* the relaxation time *T* (*H* , *k*)→0 due to the surface tension. On the other hand, for *k* →0, *T* (*H* , *k*)→*∞* unlike the case of a fluid of infinite depth. Then, according to (20),

*h <sup>k</sup>* | *<sup>k</sup>* =0 =*h <sup>k</sup>* <sup>0</sup> <sup>|</sup> *<sup>k</sup>* =0 or *∫<sup>h</sup>* (*x*, *<sup>t</sup>*)*<sup>d</sup>* <sup>2</sup> *x* =*∫h* 0(*x*)*d* <sup>2</sup> *x*, i.e. the "volume of perturbation" is kept. Using the same scales as we used for infinite liquid, we present (20) in dimensionless form

$$\mathbf{h}\_{\mathbf{k}} = \mathbf{h}\_{\mathbf{k}}^{0} \mathbf{e}^{-\frac{\mathbf{\tilde{r}} \cdot \mathbf{\tilde{r}} \left(\mathbf{e}, \mathbf{\tilde{H}}\right)}{\mathbf{\tilde{r}} \left(\mathbf{e}, \mathbf{\tilde{H}}\right)}, \quad \boldsymbol{\Theta}\left(\kappa, \mathbf{\tilde{H}}\right) = \frac{\text{ch}\left(2\kappa \mathbf{\tilde{H}}\right) + \text{l} + 2\left(\kappa \mathbf{\tilde{H}}\right)^{2}}{\text{sh}\left(2\kappa \mathbf{\tilde{H}}\right) - 2\kappa \mathbf{\tilde{H}}} \frac{2\kappa}{\kappa^{2} + \mathbf{l}}, \quad \mathbf{h}\_{\mathbf{k}} = \left[\mathbf{h}\left(\mathbf{\tilde{g}}, \mathbf{\tau}\right)\mathbf{e}^{-\mathbf{\tilde{g}} \cdot \mathbf{\tilde{r}}}, \tag{21}$$

where *H*¯ <sup>=</sup>*k*\**H* is the dimensionless thickness of layer.

#### **3.2. Finite layer over a half-space**

Consider a plane layer with thickness *H* , filled with liquid with viscosity *µ*1. Layer is located over the half-space which is filled with liquid viscosity *µ*2. At the interface between the layer and half-space we require continuity of the velocity components and continuity of the convolution of the stress tensor with the normal vector. Taking into account the gravitational forces and capillary forces on the free surface of the layer, using the technique described above, we obtain the equation of evolution of the free surface of the layer:

$$\mathbf{h}\_{\mathbf{k}} = \mathbf{h}\_{\mathbf{k}}^{0} \mathbf{e}^{-\frac{\mathbf{v}\_{\mathbf{f}} \mathbf{h} \mathbf{e} \left(\mathbf{x}, \overline{\mathbf{H}}\right)}{\mathbf{h}}}, \quad \boldsymbol{\Theta}\left(\kappa, \overline{\mathbf{H}}, \mathbf{m}\right) = \frac{2\mathbf{m} \cdot \mathbf{sh}\left(2\kappa \overline{\mathbf{H}}\right) + \left(\mathbf{m}^{2} + 1\right)\mathbf{ch}\left(2\kappa \overline{\mathbf{H}}\right) + \left(1 + 2\kappa^{2} \overline{\mathbf{H}}^{2}\right)\left(\mathbf{m}^{2} - 1\right)}{\left(\mathbf{m}^{2} + 1\right)\mathbf{sh}\left(2\kappa \overline{\mathbf{H}}\right) + 2\mathbf{m} \cdot \mathbf{ch}\left(2\kappa \overline{\mathbf{H}}\right) - 2\kappa \overline{\mathbf{H}}\left(\mathbf{m}^{2} - 1\right)} \frac{2\kappa}{\kappa^{2} + 1},\tag{22}$$

where *m*= *µ*2 *µ*1 . As expected, the characteristic relaxation time *θ*(*κ*, *<sup>H</sup>*¯, *<sup>m</sup>*) at *m*=1 or *H*¯→*<sup>∞</sup>* transferred to the corresponding expression (13) obtained for the half-space, at *m*→*∞*, (22) transferred to (21)-expression for finite layer over solid bottom. At *<sup>κ</sup>* =0, *θ*(*κ*, *<sup>H</sup>*¯, *<sup>m</sup>*)=0. Conse‐ quently, the "volume of perturbation" instantly vanishes, as in the case of the half-space.

#### **3.3. Comparing the cases of finite layer and finite layer over a half-space**

( ( )) ( )

( ( ) ( )) ( ( ) ( ) )

k g 2 k ch 2kH 1 h .

**k k**

H sh k z H z H V V , 2 sh kH 2

> ( ) ( )

( ( )) ( )

( ) ( ( ) )

z 0 2

æ ö æ ö - ç ÷ ç ÷ s +r P= m + ç ÷ ç ÷ m ++ è ø <sup>P</sup> m + è ø

> ( ) ( ) ( )

> > 2*µk*

*x* =*∫h* 0(*x*)*d* <sup>2</sup>

where *H*¯ <sup>=</sup>*k*\**H* is the dimensionless thickness of layer.

scales as we used for infinite liquid, we present (20) in dimensionless form

( ) ( ) ( ) ( )

**k k k**

+ P =P - P

0 H

sh k z H sh kz , sh kH sh kH

0 2

The evolution equation of the free surface and the initial condition has the form

t 2 t 0 2kH sh 2kh k g <sup>h</sup> h 0, h h . ch 2kH 1 2 kH 2 k <sup>=</sup>

<sup>0</sup> tTk( ) h he , -

Therefore, for *k* →*∞* the relaxation time *T* (*H* , *k*)→0 due to the surface tension. On the other hand, for *k* →0, *T* (*H* , *k*)→*∞* unlike the case of a fluid of infinite depth. Then, according to (20),

0 ,H i

Consider a plane layer with thickness *H* , filled with liquid with viscosity *µ*1. Layer is located over the half-space which is filled with liquid viscosity *µ*2. At the interface between the layer

ch 2 H 1 2 H <sup>2</sup> h h e , ,H ,h h ,e , sh 2 H 2 H <sup>1</sup> -t q k - k ++ k <sup>k</sup> <sup>=</sup> qk = = t <sup>ò</sup> k -k k +

2 k ch 2kH 1 2 kH 4 k kHsh kH ch kH

**k k**

æ ö <sup>P</sup> <sup>+</sup> + P <sup>=</sup> ç ÷ - + è ø m m

<sup>0</sup> z z0

**k k**

V 2kH sh 2kh

H

The solution of equation (19) is similar to (10)

*sh* (2*kh* )−2*kH*

**k**

where *<sup>T</sup>* (*k*)= *ch* (2*kH* ) <sup>+</sup> <sup>1</sup> <sup>+</sup> 2(*kH* )<sup>2</sup>

<sup>0</sup> <sup>|</sup> *<sup>k</sup>* =0 or *∫<sup>h</sup>* (*x*, *<sup>t</sup>*)*<sup>d</sup>* <sup>2</sup>

**3.2. Finite layer over a half-space**

At *kH* <sup>&</sup>gt; >1, *<sup>T</sup>* (*k*, *<sup>H</sup>* )= <sup>2</sup>*µk*

harmonic.

*h <sup>k</sup>* | *<sup>k</sup>* =0 =*h <sup>k</sup>*

**k**

150 Hydrodynamics - Concepts and Experiments

where

( )

<sup>2</sup> <sup>0</sup>

*<sup>σ</sup><sup>k</sup>* <sup>2</sup> <sup>+</sup> *<sup>ρ</sup><sup>g</sup>* coincides with (10), i.e. with the case of a fluid of infinite depth.

( ) ( )

2

2

**κ κ <sup>κ</sup> ξ** (21)

**k k** = (20)

*<sup>σ</sup><sup>k</sup>* <sup>2</sup> <sup>+</sup> *<sup>ρ</sup><sup>g</sup>* is the characteristic relaxation time of corresponding

*x*, i.e. the "volume of perturbation" is kept. Using the same

**κξ**


(17)

(18)

Fig.7 shows function *e* <sup>−</sup>*τ*/*<sup>θ</sup>* for the cases of finite layer and finite layer over a half-space. Qualitative differences in the behavior of this function lead to a difference in the perturbation decay for these cases.

**Figure 7.** Graph of the function *e* <sup>−</sup>*τ*/*θ* at *τ* =1 for half-space, finite layer and layer with half-space.

Fig. 8 and Fig. 9 show decrease of the amplitude for the perturbation with the same initial characteristics. The figures show that the perturbation on the surface layer over solid bottom decays faster than the perturbation on the surface layer over the half-space.

**Figure 8.** Relaxation of axisymmetric perturbation of free surface of the layer with thickness *<sup>H</sup>*¯ =1.The initial shape of perturbation *h* <sup>0</sup> = *Ae* <sup>−</sup>*<sup>ξ</sup>* <sup>2</sup> /*a* <sup>2</sup> , *a* =0.1

**Figure 9.** Relaxation of axisymmetric perturbation of free surface of the layer with thickness *<sup>H</sup>*¯ =1 over a half-space. The initial shape of perturbation *h* <sup>0</sup> = *Ae* <sup>−</sup>*<sup>ξ</sup>* <sup>2</sup> /*a* <sup>2</sup> , *a* =0.1, *m*=2.

#### **3.4. Comparing the relaxation of finite layer free surface perturbation in linear and nonlinear approximation**

Fig. 8 and Fig. 9 show decrease of the amplitude for the perturbation with the same initial characteristics. The figures show that the perturbation on the surface layer over solid bottom

**Figure 8.** Relaxation of axisymmetric perturbation of free surface of the layer with thickness *<sup>H</sup>*¯ =1.The initial shape of

**Figure 9.** Relaxation of axisymmetric perturbation of free surface of the layer with thickness *<sup>H</sup>*¯ =1 over a half-space.

, *a* =0.1, *m*=2.

/*a* <sup>2</sup>

perturbation *h* <sup>0</sup> = *Ae* <sup>−</sup>*<sup>ξ</sup>* <sup>2</sup>

152 Hydrodynamics - Concepts and Experiments

/*a* <sup>2</sup>

The initial shape of perturbation *h* <sup>0</sup> = *Ae* <sup>−</sup>*<sup>ξ</sup>* <sup>2</sup>

, *a* =0.1

decays faster than the perturbation on the surface layer over the half-space.

Using finite-element software Comsol Multiphysics we developed a 2D axisymmetric model of the viscous flow to calculate transient relaxation of the perturbation on the free surface and compare it with analytical results reported in the previous sections. We used "Laminar Two-Phase Flow, Moving Mesh" module to enable deformation of the computational domain during the solution. Full Navier-Stokes equations are solved on the moving mesh to describe deformation of the domain due to capillary forces acting along the free surface.

Computational mesh was refined near the perturbation to resolve high curvature and velocity gradients (see Fig. 10).

**Figure 10.** Initial non-deformed finite-element discretization of the domain. The initial shape of perturbation *h* <sup>0</sup> = *Ae* <sup>−</sup>*<sup>ξ</sup>* <sup>2</sup> /*a* <sup>2</sup> , *A*=1, *a* =1. Boundary conditions: 1 – axial symmetry, 2 – free surface, 3 – slip, and 4 – no-slip wall.

Multiple cases were calculated corresponding to small and large non-dimensional perturba‐ tion amplitude and width. The evolution of the perturbation amplitude with time was compared for numerical and analytical models (Fig. 11). It can be seen that the analytics and FEM are in excellent agreement for small amplitude A=0.1 for both small and large width (a=0.1, a=1.0). Some discrepancy is observed in the case of large amplitude A=1.0 due to violation of the analytic model assumptions.

**Figure 11.** Comparison of the analytical model (solid lines) and FEM (points) results

## **4. Evolution of small perturbations of half-space free boundary due to external mass flux, gravitation and capillary forces**

Let's describe the influence of external mass flow on stability of the free surface of half-space. A small perturbation on the free surface leads to a perturbation of the diffusion boundary layer above it. The corresponding mass flow perturbations on the free boundary defines the source term in the right hand side of (9) *W<sup>k</sup>* =*κVcth* (*kHD*), where *HD* is the thickness of the diffusion boundary layer, *V* is the unperturbed velocity of the liquid level rise due to external mass flow. In this case the equation of perturbation evolution has the following shape

$$\left\|\left\|\left\|\mathbf{h}\_{\mathbf{k}} + \mathbf{\boldsymbol{y}}\right\|\_{\mathbf{k}} = 0, \quad \mathbf{h}\_{\mathbf{k}} = \mathbf{h}\_{\mathbf{k}}^{0} \mathbf{c}^{\gamma(\mathbf{k})\mathbf{t}},\tag{23}$$

where *γ*(*k*)= <sup>−</sup> <sup>1</sup> <sup>2</sup>*µk* (*σ<sup>k</sup>* <sup>2</sup> <sup>+</sup> *<sup>ρ</sup>g*) <sup>+</sup> *kVcth* (*kHD*). At *<sup>γ</sup>* >0 correspondent harmonic is unstable, and at *γ* <0 it is stable. On the Fig.12 we can see the regions in the space of parameters, correspondent stable and unstable regimes. On the Fig.13 present the graphs of *γ*(*k*) for these regions. In the case (a) all the harmonics with the wave number exceeding some critical value are unstable. In the case (b) all harmonics are stable. In the case (c) there is a range of wave numbers corresponding to unstable harmonics. In the case (d) the harmonics in a certain range of wave numbers, and with a wave number exceeding a critical value are unstable.

**Figure 12.** Space of parameters with regions, corresponds to regimes a-unstable, b-stable, c-partly unstable, d-mixed.

**Figure 13.** Graphs of *γ*(*k*) for regimes: a-unstable, b-stable, c-partly unstable, d-mixed

**Figure 11.** Comparison of the analytical model (solid lines) and FEM (points) results

**external mass flux, gravitation and capillary forces**

where *γ*(*k*)= <sup>−</sup> <sup>1</sup>

154 Hydrodynamics - Concepts and Experiments

In this case the equation of perturbation evolution has the following shape

**4. Evolution of small perturbations of half-space free boundary due to**

Let's describe the influence of external mass flow on stability of the free surface of half-space. A small perturbation on the free surface leads to a perturbation of the diffusion boundary layer above it. The corresponding mass flow perturbations on the free boundary defines the source term in the right hand side of (9) *W<sup>k</sup>* =*κVcth* (*kHD*), where *HD* is the thickness of the diffusion boundary layer, *V* is the unperturbed velocity of the liquid level rise due to external mass flow.

*γ* <0 it is stable. On the Fig.12 we can see the regions in the space of parameters, correspondent stable and unstable regimes. On the Fig.13 present the graphs of *γ*(*k*) for these regions. In the case (a) all the harmonics with the wave number exceeding some critical value are unstable. In the case (b) all harmonics are stable. In the case (c) there is a range of wave numbers

<sup>0</sup> (k t) th h 0, h h e , <sup>g</sup> ¶ +g = = **k k kk** (23)

<sup>2</sup>*µk* (*σ<sup>k</sup>* <sup>2</sup> <sup>+</sup> *<sup>ρ</sup>g*) <sup>+</sup> *kVcth* (*kHD*). At *<sup>γ</sup>* >0 correspondent harmonic is unstable, and at

### **5. Conclusion**

The relaxation process of a small perturbation of the free planar surface of a viscous incom‐ pressible fluid under the action of capillary and gravitational forces was investigated. The case of half-space, layer over hard bottom and layer over the half-space was analyzed. We propose a method for qualitative analysis of the configuration of the perturbation based on spatial filtering the Fourier transform of the initial perturbation.

It is shown that under the gravity force there is an additional long-wave perturbation arising around the initial perturbation which provides a zero total volume of the perturbation. Capillary forces lead to an increase of the area affected by the perturbation. The influence of amplitude of perturbation on the rate of its relaxation was investigated numerically and compared with the analytics.

The stability of the free surface of the half-space with the mass flux at the boundary was studied.

#### **Acknowledgements**

This work was financially supported by Government of Russian Federation, Grant 074-U01.

#### **Author details**

Tamara Pogosian1 , Ivan Melikhov1 , Anastasia Shutova1 , Alexey Amosov2 and Sergey Chivilikhin1


#### **References**


**5. Conclusion**

156 Hydrodynamics - Concepts and Experiments

compared with the analytics.

**Acknowledgements**

**Author details**

Tamara Pogosian1

Sergey Chivilikhin1

**References**

studied.

The relaxation process of a small perturbation of the free planar surface of a viscous incom‐ pressible fluid under the action of capillary and gravitational forces was investigated. The case of half-space, layer over hard bottom and layer over the half-space was analyzed. We propose a method for qualitative analysis of the configuration of the perturbation based on spatial

It is shown that under the gravity force there is an additional long-wave perturbation arising around the initial perturbation which provides a zero total volume of the perturbation. Capillary forces lead to an increase of the area affected by the perturbation. The influence of amplitude of perturbation on the rate of its relaxation was investigated numerically and

The stability of the free surface of the half-space with the mass flux at the boundary was

This work was financially supported by Government of Russian Federation, Grant 074-U01.

, Anastasia Shutova1

[1] Happel S.J., Brenner H. Low Reynolds Number Hydrodynamics. Englewood Cliffs:

[2] Frenkel, J. Viscous flow of crystalline bodies under the action of surface tension. Jour‐

[3] Tanveer S., Vasconcelos G.L. Bubble Breakup in two-dimensional Stokes flow. Physi‐

, Alexey Amosov2

and

filtering the Fourier transform of the initial perturbation.

, Ivan Melikhov1

2 Corning Scientific Center, Saint-Petersburg, Russia

1 ITMO University, Saint-Petersburg, Russia

Prentice-Hall; 1965.

nal of Physics 1945;9(5) 385-391.

cal Review Letters 1994;73(21) 2845-2848.


## *Edited by Harry Edmar Schulz*

The present book on Hydrodynamics - Concepts and Experiments contains six chapters, each of which devoted to a topic considered interesting to compose this book. Although the tools used in the different chapters are somewhat superposed, a first classification may be presented as: the first two chapters describe numerical models; the third and fourth chapters use experimental tools and results; and the fourth, fifth and sixth chapters present conceptual aspects of fluid dynamics. A more detailed description is furnished in the sequence.

The six topics that compose this book may be resumed as: Flow-Structure Interactions, Hydrodynamic Models and Biological processes in Coastal Waters, Useful Cavitation, New Concepts in Hydraulic Jumps, Lagrangian Formulation of Fluid Dynamics, and Small Perturbations in Viscous Fluids.

Photo by Tim Marshall / unsplash

Hydrodynamics - Concepts and Experiments

Hydrodynamics

Concepts and Experiments

*Edited by Harry Edmar Schulz*