Open AccessArticle

Periodic Waves and Ligaments on the Surface of a Viscous Exponentially Stratified Fluid in a Uniform Gravity Field

Yuli D. Chashechkin

^1,* and

Artem A. Ochirov

Laboratory of Fluid Mechanics, Ishlinsky Institute for Problems in Mechanics RAS, 119526 Moscow, Russia

Department of Intelligent Information Radiophysical Systems, Physics Faculty, P.G. Demidov Yaroslavl State University, 150003 Yaroslavl, Russia

Author to whom correspondence should be addressed.

Axioms 2022, 11(8), 402; https://doi.org/10.3390/axioms11080402

Submission received: 11 July 2022 / Revised: 8 August 2022 / Accepted: 10 August 2022 / Published: 15 August 2022

(This article belongs to the Special Issue Singularly Perturbed Problems: Asymptotic Analysis and Approximate Solution)

Download

Browse Figures

Review Reports Versions Notes

Abstract

The theory of singular perturbations in a unified formulation is used, for the first time, to study the propagation of two-dimensional periodic perturbations, including capillary and gravitational surface waves and accompanying ligaments in the

10^{- 4} < ω < 10^{3} s^{- 1}

frequency range, in a viscous continuously stratified fluid. Dispersion relations for flow constituents are given, as well as expressions for phase and group velocities for surface waves and ligaments in physically observable variables. When the wave-length reaches values of the order of the stratification scale, the liquid behaves as homogeneous. As the wave frequency approaches the buoyancy frequency, the energy transfer rate decreases: the group velocity of surface waves tends to zero, while the phase velocity tends to infinity. In limiting cases, the expressions obtained are transformed into known wave dispersion expressions for an ideal stratified or actually homogeneous fluid.

Keywords:

waves; ligaments; gradient flow; singular solution; stratification; viscous fluid

1. Introduction

Waves on the surface of rivers, seas and oceans became one of the main objects of theoretical research as the first fundamental equations of continuum mechanics and closed systems of equations were formulated [1,2,3]. Due to a large difference in the physical properties of the atmosphere and water [4], the theory of waves evolved with the approximation of the homogeneity and immutability of the density of the contacting media. B. Franklin in the 18th century observed the water–olive oil interface motion in a ship’s lighting lamp and pointed out that the variability of density with depth should be taken into account when analyzing the wave phenomena in a liquid [5].

The results of the analysis of the first papers, which considered the variability of the density of the liquid in depth, were presented [6], but for some unknown reason, mathematicians and mechanics did not pay attention to these topics. For example, such a fundamental characteristic as the buoyancy frequency calculated in [7], which is the limiting frequency for propagating internal waves in a continuously stratified fluid, escaped from attention. Independently the buoyancy frequency was re-discovered as the natural oscillation frequency of probe balls drifting in a stratified atmosphere [8] and later it was associated with the local extreme in the spectrum of high-frequency pressure oscillations recorded by a microbarograph [9].

Due to the practical importance of the issue, for a long time, researchers have been concentrating their efforts on studying periodic phenomena in a liquid in the context of the force action of waves on obstacles and the wave resistance of bodies moving in a liquid. The scientists assumed the constancy of density and incompressibility of the liquid [10] to be the most important qualities. Nowadays the constant density approximation is still considered high priority in the description of wave processes [11].

In theoretical works much attention is paid to the study of nonlinear properties of waves. The first heuristic discovery—the theory of periodic nonlinear vortex waves [12]—is significant in wave theory development. The approach was extended in [13,14], then analyzed in [15] and it has recently been supplemented with new results [16,17]. Another large cycle of studies of nonlinear waves, initiated by observations of a solitary wave in a shipping channel [18], is successfully proceeding at the present time. The number of various model equations [19,20] for the description of nonlinear wave properties is continuously increasing [21].

Seminal works [22,23] occupied a special place among the first publications. In these studies, the parameters of infinitesimal waves were determined, the waves of finite amplitude were calculated and showed the existence of the wave transfer of matter (Stokes drift) using the methods of the regular perturbations theory. The mass transfer is caused by nonlinear effects, which are distinctly revealed in the deviation of the waveform from the ideal one, i.e., in a nonlinear wave, the crests become sharper, and the troughs become flatter. The limiting angle between the tangents to the right and left sides of the crest is calculated as

α = 120 °

in [23], and as

α = 90 °

in [13]. It depends on the nonlinearity parameter

ε = ω^{2} ζ_{0} / g

[24] when describing the propagating potential waves of finite amplitude by means of Lambert’s complex functions (

ζ_{0}

and

ω

are the amplitude and frequency of the wave,

g

is the gravity acceleration).

Calculations of the viscous attenuation of waves in a deep liquid and a channel of finite depth were carried out by asymptotic methods in a linear and nonlinear formulation [6,25,26]. Such calculations continue to be conducted at the present time using various approximations [27,28,29,30].

The “boundary layer” ideas (i.e., a mathematical model of the flow adjacent to the surface of the liquid, in which the influence of viscosity is considerably revealed), which were formulated at the beginning of the 20th century, had a significant impact on the development of the theory of waves in a liquid. The approximation of the boundary layer, accompanied by the reformulation of the defining equations [31], stimulated the development of the search for new analytical methods of finding their solutions [32]. The analysis of the wave equations by the theory of regular perturbations methods showed that in a homogeneous liquid with kinematic viscosity

ν

, a surface gravitational wave with a frequency

ω

is accompanied by a boundary layer with a specific thickness

δ_{ω} = \sqrt{2 ν / ω}

. Separate boundary layers are formed on the free surface and on the solid bottom of the channel through which the waves propagate [33].

The boundary layer greatly influences all the parameters of fluid flows, namely pressure distribution, velocity and substance transfer characteristics [34]. Corrections of the theory formulas [33] due to the viscous attenuation of waves were later calculated in the second order of the perturbation theory [35]. The analysis of the dispersion relation for surface gravitational waves and accompanying boundary layers in a viscous liquid, supplemented by calculations of the attenuation coefficient and the scale of spatial attenuation of the wave was carried out in [36].

A more complex model of a “double boundary layer”, inside which there is a periodic boundary layer with a length scale

δ_{ω} = \sqrt{2 ν / ω}

, was proposed to describe standing waves [37].

With the frequency increasing, the surface tension influences more the pattern of waves on the surface of the liquid. The surface tension coefficient

σ

is a parameter that determines the type of dispersion relation for short waves [38] and accompanying fine constituents [39]. The surface tension changes the dependences of the group and phase velocity, the attenuation coefficient and the velocity of matter transfer

U_{ρ}

on frequency

ω

, wave vector

k

or wavelength

λ

. The pattern of waves generated by a short-range localized source in a viscous liquid, taking into account the effects of surface tension, was calculated in [40].

The data of the first systematic experimental studies of surface waves, which were conducted in laboratory basins at the end of the 19th century [41], generally showed satisfactory agreement with the calculations of velocity in the liquid thickness [23]. The visualization of the velocity profile using a “marker” (i.e., a colored wake of a submerging particle) showed the existence of a drift in the near-surface and bottom layers in the direction of wave propagation and slow counterflow in the middle of the liquid layer [42]. Additional cleaning from dust and film of the target liquid surface and consideration of viscous attenuation in the calculations of the near-surface boundary layer significantly reduced the difference between the data of calculations and experiments in a laboratory channel [43]. The observations of rearrangement of the distribution pattern of initially homogeneous suspension in the field of standing gravity waves in a vertically oscillating basin are given in [44].

In the field of propagating gravitational-capillary waves (

5 < ω < 50 Hz

), the substance of the colored drop is distributed unevenly over the surface of the primary cavity [45], the same as in the case of a drop merging with a target liquid at rest [46]. In the evolutionary process a slowly drifting colored primary contact area, a near-surface jet with a vortex head and a sinking ring vortex remain in the distribution pattern of the substance. There is a pronounced fibrous structure in the distribution of the drop substance in the target fluid [45].

Widespread capillary waves, which are observed in various rivers, seas and oceans, are caused by weak wind gusts in calm [47] and are certainly present on the slopes of large gravitational waves [48] in the stormy open ocean [49] and coastal regions [50]. Capillary waves contribute to the formation of bubbles, drops and foam [51] and participate in the processes of generating sound packets when drops fall [52]. Drop impact sounds form rain noise [53], which is one of the main sources determining the acoustic background of the ocean [54].

The capillary ripples determine the roughness and the real area of the contact surface. The estimation of their influence on the transfer of momentum, energy and matter between the atmosphere and the hydrosphere is of scientific and practical interest. A continuously expanding list of scientific tasks supports the interest in studying the dynamics of capillary waves in a wide range of conditions.

The analysis of the additional dissipation of short gravitational and gravitational-capillary waves caused by “parasitic” capillary waves on the slopes of longer waves in a laboratory pool are given in [55]. The presence of a surfactant film dampens short capillary waves, which in turn amplifies decimeter waves [56]. The discussion of modern optical methods for studying capillary waves and the data of detailed experiments is carried out in [57]. A review of papers on the nonlinear interaction of coexisting waves of different frequencies within the framework of the theory of “wave turbulence” is presented in [58]. The basic mechanism of energy transfer in weak turbulence theory is validated experimentally in the gravity (four-wave interactions) and capillary (three-wave interactions) regimes. Advanced experiments enable the achievement of full spatiotemporal reconstruction of the wave field in a weakly or strong nonlinear regime, to infer wave statistics as well as waves’ nonlinear dispersion relationship, to compare with theory [58]. The study in [59] discusses the influence of attenuation on nonlinear interactions of short waves and surface waves and the properties of “wave turbulence”. In all cases, the density of the liquid is assumed to be homogeneous.

However, in real conditions, the influence of differences in atmospheric and hydrosphere temperatures, insolation, radiative heat transfer, variability in the concentration of dissolved and suspended particles and flows ensure the heterogeneity of the density profile throughout the depth of the liquid [60]. The same happens in a thin near-surface layer, where temperature changes rapidly with depth and forms a “cold” or “warm” film [61].

The given study of the properties of surface waves and accompanying fine components in a viscous stably stratified liquid is based on a complete system of fundamental equations, which was firstly collected in [62] and later was considered in [63,64]. In this paper, a reduced system of the fundamental equations in which the density variability is preserved without considering variations in temperature, salinity or pressure that cause it to change is used. This system of equations is analyzed by methods of the singular perturbations theory [65] taking into account the compatibility condition [66], enabling calculation of both the waves themselves and the family of accompanying fine constituents (i.e., ligaments).

Analysis of the solutions of the fundamental equations system in a continuously stratified fluid has shown that for periodic flows (with a frequency less than the buoyancy frequency), regular solutions describe waves. The wave frequency

ω

is related to the wave number

k

by an algebraic dispersion relation

ω = ω (k, A k)

, which may include the amplitude

A

[66,67].

Singular solutions describe ligaments perturbations—extended in some direction and fine in others—outlining wave beams. The transverse scale of the ligaments

δ_{ω}^{ν}

is determined by the kinematic viscosity of the liquid

ν

and the wave frequency

ω

δ_{ω}^{ν} = \sqrt{ν / ω}

(Stokes scale) or equivalent parameter with the buoyancy frequency

N

δ_{N}^{ν} = \sqrt{ν / N}

. The number and positions of the ligaments in space are determined by the geometry of the source, the amplitude and frequency of its oscillations. Calculations of internal wave beams and accompanying ligaments in stratified liquids and gases are given in [68,69].

In this paper, a dispersion equation in a linear approximation for two-dimensional periodic perturbations on the surface of a viscous exponentially stratified liquid is constructed for the first time. The properties of its complete solutions are analyzed. Approximate solutions are obtained by perturbation theory. Regular solutions describe waves in the range from infra-low frequency to gravitational-capillary and capillary waves. The dispersion relations continuously transform into the well-known formulas of the theory of waves in a homogeneous viscous or ideal liquid. Singular solutions characterize ligaments—fine constituents that complement waves. Graphs of the dependence of the wavelengths and ligaments on the frequency, phase and group velocity of the constituents on the wavelength are given. The results can be used to solve problems of generation and propagation of surface waves with physically justified initial and boundary conditions and comparison with a high-resolution experiment.

2. Periodic Flows in a Viscous Exponentially Stratified Fluid

We consider the propagation of periodic perturbations over the surface of a viscous exponentially stratified fluid in a uniform gravity field. The undisturbed liquid filling the lower half-space is regarded to be incompressible; meanwhile the effects of thermal conductivity and diffusion are neglected.

We perform the study in a Cartesian coordinate system

O x y z

in which the plane

O x y

coincides with the equilibrium position of a free surface of the liquid, and the axis

O z

is directed vertically upwards against the direction of the gravity acceleration

g .

The undisturbed density distribution over depth is exponential

ρ_{0} (z) = ρ_{00} \exp (- z / Λ)

. It is characterized by reference density

ρ_{00} (z_{0})

(i.e., the density value at the equilibrium level

z = z_{0}

), as well as the scale

Λ = {| d \ln ρ / d z |}^{- 1},

frequency

N = \sqrt{g / Λ}

and buoyancy period

T_{b} = 2 π / N

Mostly, changes of real fluids density are usually small and produce a small impact on the inertial properties of the flows. Nevertheless, the conservation of terms describing the stratification effects in the governing equations set is important since gravity acceleration is large. In this regard, it is useful to consider three types of medium: stratified fluids when buoyancy scale

Λ

, frequency

N

and period

T_{b}

which are included in the list of main parameters; then very weakly stratified fluids, when the scale of buoyancy substantially exceeds the values of other length scales of the problem (so called potentially homogeneous fluid) and actually homogeneous liquid whose density is assumed to be constant in the entire space [70]. Using a weak but variable density helps to save the rank of complete non-linear set and order of the linearized set of governing equations [66,71] and analyzes additional solutions which are lost in approximation of homogeneous fluid.

Depending on the magnitude of the density gradient created by the corresponding temperature or salinity distributions, it is customary to distinguish strongly and weakly stratified fluids, as well as potentially and actually homogeneous fluids. The values of the characteristic physical quantities are given in Table 1 [70].

Strong stratification (

T_{b} ~ 10 s

) is created in laboratory installations. Relatively weak stratification (

T_{b} ~ 10 \min

) is observed in natural conditions. Figure 1 shows the density–depth dependencies for exponential stratification (solid line) and for linear stratification (dotted line) for a strongly stratified fluid

N = 1 s^{- 1}

To compare the properties of exponentially and linearly stratified media, we analyze the changes in the density gradient with depth. For exponential stratification, the magnitude of the density gradient depends on the vertical coordinate

z

\frac{d ρ}{d z} = - \frac{ρ_{00}}{Λ} e^{- \frac{z}{Λ}}

(1)

For linear stratification, the density gradient does not depend on the depth:

\frac{d ρ}{d z} = - \frac{ρ_{00}}{Λ}

(2)

For small variations

z ≪ Λ

the dependencies of the density gradient are not distinguishable for exponential and linear stratification. The exponent can be represented as a Taylor series expansion over a small parameter

z / Λ

and expression (1) will be written as follows:

\frac{d ρ}{d z} = - \frac{ρ_{00}}{Λ} e^{- \frac{z}{Λ}} ≃ - \frac{ρ_{00}}{Λ} (1 - \frac{z}{Λ} + \frac{z^{2}}{2 Λ^{2}})

(3)

The numerical values for the stratification scale for a strongly stratified fluid are

Λ = 9.81 \cdot 10^{2} cm

, and for a weakly stratified fluid

Λ = 9.81 \cdot 10^{6} cm

. At depths of less than 10% of the stratification scale, it can be argued with a high degree of accuracy that the calculations in the model of a linearly stratified fluid and exponentially stratified fluid coincide.

We study below two-dimensional periodic flows of the form

A = A_{0} \exp i (k x - ω t)

with a positive definite frequency

ω

and a complex wave number

k

, the imaginary part of which characterizes the spatial attenuation of the flow. The disturbances which are homogeneous in the direction of the transverse horizontal coordinate

y

are selected. They include the displacement of the free surface position

ζ (x, t)

. The velocity

u = u_{x} e_{x} + u_{z} e_{z}

with horizontal and vertical velocity components

u_{x}, u_{z}

in an incompressible fluid (

div u = 0

) is represented by derivatives of the stream function

ψ

u_{x} = \partial_{z} ψ, u_{z} = - \partial_{x} ψ

(4)

The expression for the density of the liquid

ρ = ρ_{00} (r_{0} (z) + s (x, z, t))

replacing the equation of state in [66] includes the initial distribution

r_{0} (z)

and the perturbation caused by the examined periodic flow

s (x, z, t)

Taking into account these assumptions, the system of continuity and Navier–Stokes equations is extremely simplified and takes the traditional form [62,63,69]:

z < ζ : {\begin{cases} ρ = ρ_{00} (r (z) + s (x, z, t)) \\ ρ \partial_{t} u + ρ (u \cdot \nabla) u - ρ ν Δ u = ρ g - \nabla P \\ \partial_{t} ρ + div (ρ u) = 0 \end{cases}

(5)

Kinematic and dynamic boundary conditions on a perturbed surface are traditional [62]:

z = ζ : {\begin{cases} \partial_{t} (z - ζ) + u \cdot \nabla (z - ζ) = 0 \\ τ \cdot (n \cdot \nabla u) + n (τ \cdot \nabla u) = 0 \\ P - P_{0} - σ \nabla \cdot n - 2 ρ ν n (n \cdot \nabla u) = 0 \end{cases} n = \frac{\nabla (z - ζ)}{| \nabla (z - ζ) |} = \frac{- \partial_{x} ζ e_{x} + e_{z}}{\sqrt{1 + {(\partial_{x} ζ)}^{2}}}, τ = \frac{e_{x} + \partial_{x} ζ e_{z}}{\sqrt{1 + {(\partial_{x} ζ)}^{2}}},

(6)

where

P

is the hydrodynamic pressure,

σ = γ ρ_{00}

is the total coefficient of the surface tension of the liquid,

n

and

τ

are the vectors of the external normal and tangent to the free surface of the liquid, respectively.

The studied problem includes the following dimensional parameters: density

ρ

and its gradient

\frac{d ρ}{d z}

, dynamic

μ

and density-normalized kinematic viscosity

ν = μ / ρ

, gravity acceleration

g

, surface tension coefficient

σ

and density-normalized surface tension coefficient

γ = σ / ρ

, amplitude

A

, frequency

ω

period

T_{ω} = 2 π / ω

and wavelength

λ

. The wave is also characterized by a wave vector

k

, the modulus of which is related to the wavelength

λ = 2 π / | k |

In the transition to a mathematical description, physical quantities serve as the basis for the introduction of characteristic scales of length, time and speed. These scales are used in the physical interpretation of the results obtained, assessing the degree of influence of various physical factors.

The degree of relative influence of viscosity and gravity on fluid flows is characterized by the scale

δ_{g}^{ν} = \sqrt[3]{ν^{2} / g}

that appears in [71].

Traditionally, stratification is characterized by its own length scale which takes into account or neglect the effects of compressibility, the frequency

N = \sqrt{g / Λ}

, period

T_{b} = 2 π / N

and the inverse value of the buoyancy frequency

T_{N} = T_{b} / 2 π

. The stratification itself is characterized by a combinational viscous wave scale

L_{N}^{g ν} = {(g ν)}^{1 / 3} / N

[72,73], as well as a length scale

δ_{N}^{ν} = \sqrt{ν / N}

– a functional analogue of the dissipative Stokes microscale

δ_{ω}^{ν} = \sqrt{ν / ω}

The dualism of the nature of surface waves, whose properties depend on the gravity acceleration

g

, normalized coefficients of surface tension

γ

and kinematic viscosity

ν

results in the introduction of several scales, including the capillary length

δ_{g}^{γ} = \sqrt{γ / g}

. The inequality

λ < δ_{g}^{γ}

indicates the severity of surface force effects, and

λ > > δ_{g}^{γ}

—indicates the severity of gravitational effects.

The structure of near-surface flows is also characterized by a known proper velocity scale

v_{c} = γ / ν

and an additional time scale of kinematic nature

T_{ν g}^{γ} = v_{c} / g = γ / ν g

Wave propagation is characterized by group

c_{g} = \frac{\partial ω}{\partial k}

and phase velocity

c_{p h} = \frac{ω k}{k^{2}}

Ratios of uniform scales form a set of dimensionless parameters, parts of which are used in further calculations.

2.1. Equations of Periodic Flows and Dispersion Relations for Plane Infinitesimal Waves

To simplify expressions, we accept the initial density distribution to be exponential, as it is customary in most models:

ρ_{0} (z) = ρ_{00} r (z) = ρ_{00} \exp (- z / Λ)

(7)

The change of variables proposed in [74] allows the transformation of the equations of motion for an arbitrary smooth density profile into equations with constant coefficients that determine the preferential choice of the exponential density profile.

Further calculations are carried out in the Boussinesq approximation, considering the smallness of absolute density variations in stratified media (in particular, a small value of the wavelength and buoyancy scale ratio

C = λ / Λ

). The density variations are neglected everywhere, except for the term with a large coefficient (i.e., the gravity acceleration

g

). In this approximation, the system (5) takes the form:

{\begin{cases} \partial_{t z} ψ + \partial_{z} ψ \partial_{x z} ψ - \partial_{x} ψ \partial_{z z} ψ - ν \partial_{z} Δ ψ + \partial_{x} P = 0 \\ g (r + s) - g - \partial_{t x} ψ + \partial_{x} ψ \partial_{x z} ψ - \partial_{z} ψ \partial_{x x} ψ - ν \partial_{x} Δ ψ + \partial_{z} P = 0 \\ \partial_{t} s + \partial_{z} ψ \partial_{x} s - \partial_{x} ψ \partial_{z} (r + s) = 0 \end{cases}

(8)

In a linear approximation the system (8) reduces to:

{\begin{cases} \partial_{t z} ψ - ν \partial_{z} Δ ψ + \partial_{x} P = 0 \\ g (r + s) - g - \partial_{t x} ψ - ν \partial_{x} Δ ψ + \partial_{z} P = 0 \\ \partial_{t} s - \partial_{x} ψ \partial_{z} r = 0 \end{cases}

(9)

Cross-differentiation of spatial coordinates of the upper and middle equations of the system (9) allows getting rid of pressure:

{\begin{cases} \partial_{t} Δ ψ - g \partial_{x} s - ν Δ Δ ψ = 0 \\ \partial_{t} s - \partial_{x} ψ \partial_{z} r = 0 \end{cases}

(10)

Subtraction of the second equation, differentiated by coordinate and multiplied by a coefficient

g

from the first equation of the system (10), differentiated by time, allows getting rid of the function

s (x, z, t)

\partial_{t t} Δ ψ - ν \partial_{t} Δ Δ ψ + \partial_{x x} ψ \partial_{z} r = 0

(11)

Equation (11) takes the next form for an exponentially stratified fluid (7):

\partial_{t t} Δ ψ - ν \partial_{t} Δ Δ ψ + N^{2} \exp (- z / Λ) \partial_{x x} ψ = 0, N = \sqrt{g / Λ}

(12)

The boundary conditions (6) for all variables of the infinitesimal waves

A = A_{0} \exp (i k x - i ω t)

are traditionally carried away from the wave surface

z = ζ

to the unperturbed level

z = 0

and take the form

z = 0 {\begin{cases} \partial_{t} ζ + \partial_{x} ψ = 0 \\ - ρ g ζ + ρ_{00} P + 2 ρ_{00} ν \partial_{x z} ψ + σ \partial_{x x} ζ = 0 \\ \partial_{z z} ψ - \partial_{x x} ψ = 0 \end{cases}

(13)

The small coefficient (kinematic viscosity

ν

), which ensures slow attenuation of the wave as it propagates in the liquid, allows the classification of the system (12) as a system of singularly perturbed equations and the application of the theory of singular perturbations [65] for its analysis. The compatibility condition, which requires the analysis of all the roots of the dispersion equation and the system (12) solutions, should be taken into account [66].

For surface waves with

ω > N

, the solution of Equation (12) is sought in the form of plane waves:

ψ = (A_{+} e^{i k_{x} x - i ω t} + A_{-} e^{- i k_{x} x - i ω t}) (e^{k_{z} z} + β e^{k_{l} z})

(14)

Here

k_{z}

corresponds to the regular solution of the dispersion equation, and

k_{l}

corresponds to the singular solution of the dispersion equation.

We obtain the shape of the deviation of the free surface of the liquid, the condition of the relationship between the amplitudes of the regular and singular component by substituting (14) into Equation (12) and boundary conditions (13):

\begin{array}{l} ζ = \frac{k_{x}}{ω} (A_{+} e^{i k_{x} x - i ω t} - A_{-} e^{- i k_{x} x - i ω t}) (1 + β) \\ β = - \frac{k_{x}^{2} + k_{z}^{2}}{k_{x}^{2} + k_{l}^{2}} \end{array}

(15)

as well as the dispersion relations:

{\begin{cases} - N^{2} k_{x}^{2} + e^{z / Λ} (k_{x}^{2} - k_{z, l}^{2}) ω (i k_{x}^{2} ν - i k_{z, l}^{2} ν + ω) = 0 \\ g k_{x}^{2} - 3 i k_{x}^{2} k_{z, l} ν ω + k_{z, l} (i k_{z, l}^{2} ν - ω) ω + k_{x}^{4} γ = 0, γ = σ / ρ_{00} \end{cases}

(16)

Calculated for the first time the expressions (16) are transformed into the relations for a viscous homogeneous fluid and for an ideal exponentially stratified fluid in limiting transitions

N \to 0

and

ν \to 0

respectively. Similar dispersion relations for internal waves and ligaments in the thickness of a stratified fluid were presented in [75].

2.2. Solution of the Dispersion Equation

It is convenient to analyze the dispersion relations (13) in a dimensionless form. The time scale will be the parameter, which is inverse to the buoyancy frequency

N

, and the viscous wave scale

L_{ν} = {(ν g)}^{1 / 3} / N

is chosen as the length scale. The ratios of the eigenscales of the problem

δ = {(\frac{δ_{g}^{γ}}{δ_{N}^{ν}})}^{2} = \frac{γ}{g} \cdot \frac{N}{v}

and

ε = \frac{L_{v}}{Λ} = \frac{N ν^{1 / 3}}{g^{2 / 3}}

are used to construct new dimensionless parameters

δ, ε

involved in further calculations. Then, expressions (16) could be rewritten in a dimensionless form:

\begin{matrix} {\begin{cases} - k_{* x}^{2} + i e^{z / Λ} {(k_{* x}^{2} - k_{* z, l}^{2})}^{2} ε ω_{*} + e^{z / Λ} (k_{* x}^{2} - k_{* z, l}^{2}) ω_{*}^{2} = 0 \\ k_{* x}^{2} + k_{* x}^{4} δ ε - 3 i k_{* x}^{2} k_{* z, l} ε^{2} ω_{*} + i k_{* z, l}^{3} ε^{2} ω_{*} - k_{* z, l} ε ω_{*}^{2} = 0 \end{cases} \\ δ = N γ / ν g, ε = N ν^{1 / 3} / g^{2 / 3} \end{matrix}

(17)

The upper dispersion relation in (17) has regular solutions, which we have denoted

k_{* z}

, and singular solutions, which are denoted

k_{* l}

. For a large number of real liquids, the parameter

ε ≪ 1

and regular solutions are found by direct decomposition over a small parameter

ε

k_{* z} = k_{* 0 z} + ε k_{* 1 z} + ε^{2} k_{* 2 z} + \dots

(18)

Substituting (18) into the upper expression in (17) we obtain up to the terms of the order

O (ε^{2})

k_{* z} = \pm \frac{\sqrt{ω_{*}^{2} - e^{- z / Λ}}}{ω_{*}} k_{* x} \pm ε \frac{i e^{- 2 z / Λ} k_{* x}^{3}}{2 ω_{*}^{4} \sqrt{ω_{*}^{2} - e^{- z / Λ}}}

(19)

Since the equations in the system (17) have the fourth degree, it is necessary to find two more solutions, which are singular perturbed ones in the form of decomposition [65]:

k_{* l} = ε^{- η} (k_{* 0 l} + ε k_{* 1 l} + ε^{2} k_{* 2 l} + \dots)

(20)

The parameter

η

in (20) is chosen in such a way that at the highest degree the main term of the decomposition persists. Substituting instead of

k_{* l}

in the upper expression in (17)

k_{* l} = ε^{- η} k_{* 0 l}

and equating the exponents

ε

in different terms we get:

\begin{array}{l} 1 - 4 η = - 2 η \\ 1 - 4 η = 1 - 2 η \\ 1 - 4 η = 0 \\ 1 - 2 η = - 2 η \\ 1 - 2 η = 0 \\ - 2 η = 0 \end{array}

(21)

At the highest degree, the main term of the decomposition remains only at

η = 1 / 2

. Substituting the value

η = 1 / 2

in (20) and then in (17) up to the terms of the order

O (ε^{3 / 2})

we get:

k_{* l} = \pm \frac{(1 - i) \sqrt{ω_{*}}}{\sqrt{2 ε}} \pm \frac{(1 + i) (e^{- z / Λ} + ω_{*}^{2}) k_{* x}^{2}}{2 \sqrt{2} ω_{*}^{5 / 2}} \sqrt{ε}

(22)

Leaving only the main terms of

ε

in the lower dispersion relation (17), we obtain the dispersion equation for the regular part of the solution:

k_{* x} + (k_{* x}^{3} δ - ω_{*} \sqrt{ω_{*}^{2} - 1}) ε = 0

(23)

For the singular part of the solution, leaving only the main terms of the expansion of

ε

, we obtain:

1 + k_{* x}^{2} δ ε + \frac{(1 + i) (1 - 2 ω_{*}^{2})}{\sqrt{2 ω_{*}}} ε^{3 / 2} = 0

(24)

The solution of the dispersion Equation (23) for the regular wave part is:

\begin{matrix} k_{* x} = - \frac{2^{1 / 3}}{α} + \frac{α}{3 \cdot 2^{1 / 3} δ ε} \\ k_{* x} = \frac{1 \pm i \sqrt{3}}{2^{2 / 3} α} - \frac{(1 \mp i \sqrt{3}) α}{6 \cdot 2^{1 / 3} δ ε} \\ α = {(27 δ^{2} ε^{3} ω_{*} \sqrt{ω_{*}^{2} - 1} + \sqrt{108 δ^{3} ε^{3} + 729 δ^{4} ε^{6} ω_{*}^{2} (ω_{*}^{2} - 1)})}^{1 / 3} \end{matrix}

(25)

From the condition of the physical realization of the solution (i.e., the damping of the flow with depth) it follows that only roots with

Re (k_{* z}) > 0

possess a physical meaning. Decomposing the solution (25) for the regular wave part into a Taylor series, we obtain:

\begin{array}{l} k_{* x} = ω_{*} ε \sqrt{ω_{*}^{2} - 1} + O (ε^{4}) \\ k_{* x} = \pm \frac{i}{\sqrt{δ ε}} - ω_{*} ε \frac{\sqrt{ω_{*}^{2} - 1}}{2} + O (ε^{5 / 2}) \end{array}

(26)

With consideration of (19) and the condition

Re (k_{* z}) > 0

, it can be seen that we implement only one root of (26) and finally for the regular wave solution we get:

\begin{matrix} k_{* x} = - \frac{2^{1 / 3}}{α} + \frac{α}{3 \cdot 2^{1 / 3} δ ε} \\ k_{* z} = \frac{\sqrt{ω_{*}^{2} - e^{- z / Λ}}}{ω_{*}} (- \frac{2^{1 / 3}}{α} + \frac{α}{3 \cdot 2^{1 / 3} δ ε}) + ε \frac{i e^{- 2 z / Λ} {(- \frac{2^{1 / 3}}{α} + \frac{α}{3 \cdot 2^{1 / 3} δ ε})}^{3}}{2 ω_{*}^{4} \sqrt{ω_{*}^{2} - e^{- z / Λ}}} \end{matrix}

(27)

We perform similar transformations for singular ligament-solution. The solutions (24), which are found taking into account (22) and the conditions of the physical implementation of the roots

Re (k_{* l}) > 0

of the singular ligament, take the form:

\begin{array}{l} k_{* x} = \pm \sqrt{{(δ ε)}^{- 1} (((1 + i) (2 ω_{*}^{2} - 1) / \sqrt{2 ω_{*}}) ε^{3 / 2} - 1)} \\ k_{* l} = \frac{(1 - i) 2 δ ω_{*}^{3} + (1 + i) (e^{- z / Λ} + ω_{*}^{2}) (((1 + i) (2 ω_{*}^{2} - 1) / \sqrt{2 ω_{*}}) ε^{3 / 2} - 1)}{2 \sqrt{2} ω_{*}^{5 / 2} δ \sqrt{ε}} \end{array}

(28)

Expressions (27), (28) describe new dispersion relations for surface waves (27) and associated ligaments (28) in a viscous exponentially stratified fluid.

2.3. Low Frequency Waves

The solution for low frequency waves (

ω < N

) is sought as the sum of the propagation of gravitational waves (waves with wave number components

k_{z}

) and fine perturbation (with wave number component

k_{s}

) as well:

ψ = (A_{+} e^{i k_{x} x - i ω t} + A_{-} e^{- i k_{x} x - i ω t}) (α e^{i k_{z} z} + β e^{- i k_{z} z} + χ e^{k_{s} z})

(29)

Substituting (26) into (9) leads to dispersion relations:

{\begin{cases} - N^{2} k_{x}^{2} + e^{z / Λ} (k_{x}^{2} + k_{z}^{2}) ω (i k_{x}^{2} ν + i k_{z}^{2} ν + ω) = 0 \\ - N^{2} k_{x}^{2} + i e^{z / Λ} (k_{s}^{2} - k_{x}^{2}) ω (k_{s}^{2} ν - k_{x}^{2} ν + i ω) = 0 \end{cases}

(30)

Substituting solution (29) into the boundary conditions (13), we obtain an expression for the shape of the free surface and additional relations. These additional relations determine the relationship between the amplitudes of the various components of periodic motion:

ζ = (A_{+} e^{i k_{x} x - i ω t} - A_{-} e^{- i k_{x} x - i ω t}) (α + β + χ) k_{x} / ω

(31)

(α + β) (k_{x}^{2} - k_{z}^{2}) + χ (k_{x}^{2} + k_{s}^{2}) = 0

(32)

k_{x}^{2} (g + γ k_{x}^{2}) (α + β + χ) + ν ω (k_{z} (α - β) (3 k_{x}^{2} + k_{z}^{2}) + i χ k_{s} (k_{s}^{2} - 3 k_{x}^{2})) - ω^{2} (i k_{z} (α - β) - χ k_{s}) = 0

(33)

In a dimensionless form, the dispersion relations (30) take the form:

{\begin{cases} - k_{* x}^{2} + i e^{z / Λ} {(k_{* x}^{2} - k_{* s}^{2})}^{2} ε ω_{*} + e^{z / Λ} (k_{* x}^{2} - k_{* s}^{2}) ω_{*}^{2} = 0 \\ - k_{* x}^{2} + i e^{z / Λ} {(k_{* x}^{2} + k_{* z, l}^{2})}^{2} ε ω_{*} + e^{z / Λ} (k_{* x}^{2} + k_{* z, l}^{2}) ω_{*}^{2} = 0 \end{cases}

(34)

The expressions for the surface component of periodic motion, up to the notation, coincide with the dispersion equations of surface waves discussed above. For low frequency waves, we will look for regular wave solutions in the form of expansion (18) and up to the terms of the order

O (ε^{2})

k_{* z} = \pm \frac{\sqrt{e^{- z / Λ} - ω_{*}^{2}}}{ω_{*}} k_{* x} \mp ε \frac{i e^{- 2 z / Λ} k_{* x}^{3}}{2 ω_{*}^{4} \sqrt{ω_{*}^{2} - e^{- z / Λ}}}

(35)

Similarly to surface waves case, we find singular solutions of the dispersion equation for low frequency waves using decomposition (20). Just as in the case of high frequency surface waves

ω > N

, the exponent of the degree

η = 1 / 2

is the only one that satisfies the condition of the main term of the decomposition at the highest degree

k_{l}

. Making a substitution up to the terms of the order

O (ε^{3 / 2})

we get:

k_{* l} = \pm \frac{(1 + i) \sqrt{ω_{*}}}{\sqrt{2 ε}} \mp \frac{(1 - i) (e^{- z / Λ} + ω_{*}^{2}) k_{* x}^{2}}{2 \sqrt{2} ω_{*}^{5 / 2}} \sqrt{ε}

(36)

From the relation (32) we obtain that

χ = (α + β) \frac{k_{z}^{2} - k_{x}^{2}}{k_{x}^{2} + k_{s}^{2}}

(37)

α + β + χ = (α + β) \frac{k_{s}^{2} + k_{z}^{2}}{k_{x}^{2} + k_{s}^{2}}

(38)

From the ratio (33) follows:

α + β = - \frac{2 i (k_{s}^{2} + k_{x}^{2}) k z α (3 k_{x}^{2} ν + k_{z}^{2} ν - i ω) ω}{i k_{x}^{2} (k_{s}^{2} + k_{z}^{2}) (g + k_{x}^{2} γ) - (k_{s} + i k_{z}) (- k_{s}^{2} k_{x}^{2} + 3 k_{x}^{4} + 4 i k_{s} k_{x}^{2} k_{z} + (k_{s}^{2} + k_{x}^{2}) k_{z}^{2}) ν ω - (k_{s} - i k_{z}) (i k_{x}^{2} + k_{s} k_{z}) ω^{2}}

(39)

In the low viscosity approximation, we obtain that

α + β ≃ 0, χ ≃ - (α + β) = 0

(40)

The relations (40) correspond to the situation when all the energy is concentrated in gravitational waves.

2.4. Periodic Flows on the Surface of a Viscous Exponentially Stratified Liquid

First of all, let us consider the dependence of the wavelength

λ

on the frequency of the wave motion

ω

. We define the wavelength as follows:

λ = 2 π / \sqrt{Re (k_{x}^{2}) + Im (k_{z}^{2})}

(41)

This method of determination of the wavelength is due to the fact that the imaginary part of the wave number

k_{x}

component and the real part of the wave number

k_{z}

component are responsible for the spatial attenuation of motion and are not impactful to wave propagation. Substituting the dispersion relations (27) in (32), we can construct the desired dependencies. Figure 2 shows the dependence of the wavelength on the frequency of wave motion

λ (ω)

The velocity of movement of the phase front of structures (wave and ligament) and the rate of energy transfer are of particular interest. The phase front moves with the phase velocity of the wave (and its analogue for the singular solution):

c_{p h} = ω k / {| k |}^{2}

(42)

The energy is transferred with a group velocity

c_{g r}

, which is defined as follows:

c_{g r} = {{(\frac{\partial k_{x}}{\partial ω})}^{- 1}, {(\frac{\partial k_{z, l}}{\partial ω})}^{- 1}}

(43)

Figure 3 shows the dependencies for the phase (dashed line) and group (solid line) velocities.

Figure 2 and Figure 3 show that viscosity has a noticeable effect on capillary waves with a wavelength

λ < δ_{g}^{γ}

, and the stratification influences the waves with frequencies close to the buoyancy frequency

N

. With the advent of stratification, a forbidden part of the spectrum appears for surface waves (with frequencies lower than the buoyancy frequency). With the tending to the stratification frequency, the group velocity goes to zero, and the phase velocity tends to infinity. A similar pattern is observed when electromagnetic waves propagate in waveguides. The size of the waveguide sets a certain critical size, beyond which the electromagnetic wave does not propagate in the waveguide channel. In stratified fluids, this mechanism is not due to external geometric constraints, but is determined by the characteristics of the medium (stratification).

3. Reduction to Approximation of Actually Homogeneous Fluid

The fine constituents of periodic flows—ligaments—are a consequence of the dissipative properties of the medium, which exist not only in stratified, but also in homogeneous liquids. They are described by singular solutions of dispersion equations, the appearance of which can be observed experimentally in the structure of flows in an inhomogeneous medium. If the effects associated with buoyancy are neglected, then the density of the liquid can be considered as actually homogeneous:

ρ = ρ_{00} \equiv const

(44)

and the basic equations of motion are simplified. For a viscous homogeneous liquid in a 2D formulation in a linear approximation equations for the stream function are shortened to:

\partial_{t} Δ ψ - ν Δ Δ ψ = 0

(45)

Equation (45) can also be obtained by performing a limiting transition

N \to 0

in (12). The boundary conditions that are removed to the equilibrium surface in the case of a homogeneous liquid are as follows:

z = 0 {\begin{cases} \partial_{t} ζ + \partial_{x} ψ = 0 \\ - g ζ + P + 2 ν \partial_{x z} ψ + γ \partial_{x x} ζ = 0 \\ \partial_{z z} ψ - \partial_{x x} ψ = 0 \end{cases}

(46)

Similarly to the basic equations of motion, the boundary conditions (46) can also be obtained from the boundary conditions (13) using a limit transition

N \to 0

. Substitution of the solution in the form of a propagating wave

ψ = A_{+} \exp (k_{z} z + i k_{x} x - i ω t) + A_{-} \exp (k_{z} z - i k_{x} x - i ω t)

(47)

in (45) leads to dispersion relations between the components of the wave number:

(k_{x}^{2} - k_{z}^{2}) (ν (k_{x}^{2} - k_{z}^{2}) - i ω) = 0

(48)

The relation (48) naturally decomposes into two solutions:

\begin{array}{l} k_{z}^{2} = k_{x}^{2} \\ k_{l}^{2} = k_{x}^{2} + \frac{i ω}{ν} \end{array}

(49)

Here, the reassignment of one solution

k_{z}

into

k_{l}

is introduced. The component

k_{z}

corresponds to the wave solution, and the component

k_{l}

corresponds to the ligament solution. Taking into account (49), the solution of the problem in the form of a propagating wave is transformed into:

ψ = (A_{+} \exp (i k_{x} x - i ω t) + A_{-} \exp (- i k_{x} x - i ω t)) (\exp (k_{z} z) + B \exp (k_{l} z))

(50)

Substituting (50) into the kinematic boundary condition, we obtain the relationship between the amplitudes of the velocity field and the deviation of the free surface from the equilibrium value:

ζ = \frac{k_{x}}{ω} (1 + B) (A_{-} \exp (- i k_{x} x - i ω t) - A_{+} \exp (i k_{x} x - i ω t))

(51)

and substitution of (50) into the dynamic condition for tangential tensions leads to the expression for the amplitude of the ligament component:

B = - \frac{(k_{x}^{2} + k_{z}^{2})}{(k_{x}^{2} + k_{l}^{2})} = - 2 k_{x}^{2} {(2 k_{x}^{2} + \frac{i ω}{ν})}^{- 1}

(52)

Getting rid of the pressure in the dynamic boundary condition and taking normal components, we rewrite it as:

ν \partial_{t z} Δ ψ - \partial_{t t z} ψ + g \partial_{x x} ψ + 2 ν \partial_{t x z z} ψ - γ \partial_{x x x x} ψ = 0

(53)

Substituting the solution (50) into the boundary conditions (53) taking into account (49) leads to the dispersion relation for the wave constituent:

γ k_{z}^{3} - 2 i ν ω k_{z}^{2} + g k_{z} - ω^{2} = 0

(54)

and ligament constituent:

γ {(k_{l}^{2} - \frac{i ω}{ν})}^{2} - 2 i ν ω (k_{l}^{2} - \frac{i ω}{ν}) k_{l} + g (k_{l}^{2} - \frac{i ω}{ν}) + ω k_{l} (i ν k_{l}^{2} - ω) = 0

(55)

It can be noticed that Equations (54) and (55) are also obtained from (16) with a limit transition

N \to 0

. Further we will transfer the dispersion relations (54) and (55) to dimensionless forms in the same way as it has been done in the previous paragraph. We will choose our own viscous scale

δ_{g}^{ν} = \sqrt[3]{ν^{2} / g}

as the length scale, and the ratio as the time scale will be

T_{ν g}^{γ} = γ / ν g

and the decomposition parameter

ε = \sqrt[3]{{{(δ_{g}^{ν})}^{6} / (δ_{g}^{γ})}^{6}} = \sqrt[3]{(ν^{4} g) / γ^{3}}

In a dimensionless form, the conventional dispersion Equation [15] for the wave component of the solution will be rewritten taking into account (54):

k_{* x}^{3} - 2 i ε^{2} ω_{*} k_{* x}^{2} + ε k_{* x} - ε^{3} ω_{*}^{2} = 0

(56)

The solution of Equation (56) is obtained by the standard method and takes the form:

k_{x}^{*} = \frac{2}{3} i ε^{2} ω^{*} - \frac{ε (3 + 4 ε^{3} ω^{*}^{2})}{3 β} + \frac{1}{3} β

(57)

k_{x} = \frac{2}{3} i ε^{2} ω^{*} \pm \frac{1}{6} i (\pm i + \sqrt{3}) β + \frac{(1 \pm i \sqrt{3}) ε (3 + 4 ε^{3} ω^{*}^{2})}{6 \cdot 2^{2 / 3} β} β = {(- 8 i ε^{6} ω^{*}^{3} + \frac{9}{2} ε^{3} ω^{*} (- 2 i + 3 ω^{*}) + \frac{3}{2} \sqrt{3} \sqrt{ε^{3} (4 - 32 i ε^{6} ω^{*}^{5} + ε^{3} ω^{*}^{2} (4 + 9 ω^{*} (- 4 i + 3 ω^{*})))})}^{1 / 3}

(58)

Decomposition (57) into a Taylor series by a small parameter

ε

at least up to the summands

O (ε^{9 / 2})

gives the expression:

k_{* z} = ε^{2} ω_{*}^{2} + O (ε^{9 / 2})

(59)

Decomposition of the roots (58) into a Taylor series by a small parameter

ε

at least up to the summands

O (ε^{9 / 2})

gives the expression:

k_{* z} = \pm \frac{1}{8} i \sqrt{ε} (8 \pm 4 ε^{3 / 2} ω_{*} (2 + i ω_{*}) + ε^{3} ω_{*}^{2} (4 - 4 i ω_{*} + 3 ω_{*}^{2})) + O (ε^{9 / 2})

(60)

As it follows from the condition of the physical realization of the roots

Re (k_{z}) > 0

, considering the positive definite frequency of the wave motion, only one root, which is given by the expression (57), can be physically realized This solution describes the wave part of a periodic motion in a liquid.

To analyze the dispersion equation for a ligament solution (55) and carry out the nondimensionalization procedure, the characteristic scales can be selected the same as in the wave solution. Then in a dimensionless form, the dispersion equation (55) is written as:

k_{* l}^{4} - 2 i ε ω_{*} k_{* l}^{2} - ε^{2} ω_{*}^{2} - 3 i ε^{2} k_{* l}^{3} ω_{*} - 3 ε^{3} k_{* l} ω_{*}^{2} + ε k_{* l}^{2} - i ε^{2} ω_{*} + i ε^{2} ω_{*} k_{* l}^{3} - ε^{3} k_{* l} ω_{*}^{2} = 0

(61)

Equation (61) has four roots. The analysis shows that the spatial attenuation condition is satisfied for one root only. Due to the cumbersomeness of the expression, the calculated solution is not given here.

Let us construct the dispersion dependences of the components of periodic motion in log–log scales in dimensional variables for liquids with glycerin and water parameters. Figure 4 shows the dependence of the wavelength

λ

on the frequency for the wave component (a) and the analog of the wavelength

δ_{i}

on the frequency for the ligament component (b). The letter

(W)

indicates the dependencies for water and the letter

(G l)

indicates the dependencies for glycerin.

Figure 5 shows the dependences of the phase and group velocity on the frequency of periodic motion for the wave component of periodic motion (a) and the ligament associated with the wave component (b). The dependences are constructed for liquids with the parameters of water and glycerin. Figure 6 shows similar dependencies, but on the wavelength. There are several remarkable velocities in a viscous homogeneous liquid. The minimum group velocity and the velocity at which the group and phase velocities are compared. We show that the velocities are compared when the value of the phase velocity is minimal. The extremum condition for the phase velocity will be written as:

\partial_{ω} c_{p h} = \partial_{ω} (\frac{ω}{k}) = \frac{1}{k^{2}} (k - ω \partial_{ω} k) = 0

(62)

The ratio (62) leads us to equality:

c_{p h} = \frac{ω}{k} = {(\partial_{ω} k)}^{- 1} = c_{g r}

(63)

The minimum value of the group velocity for a liquid with water parameters is

c_{g r \min}^{W} = 17.71 cm / s

and is achieved at frequency

ω_{g r \min}^{W} = 40.53 s^{- 1}

and wavelength

λ_{g r \min}^{W} = 4.33 cm

. For a liquid with glycerin parameters, the corresponding values are as follows:

c_{g r \min}^{G l} = 29.39 cm / s

at the frequency

ω_{g r \min}^{W} = 21.05 s^{- 1}

and wavelength

λ_{g r \min}^{G l} = 14.88 cm

. The minimum phase velocity for a liquid with water parameters is

c_{p h \min}^{W} = 23.05 cm / s

at the frequency

ω_{p h \min}^{W} = 84.82 s^{- 1}

and wavelength

λ_{p h \min}^{W} = 1.71 cm

. For a liquid with glycerin parameters, the corresponding values are as follows:

c_{p h \min}^{G l} = 40.85 cm / s

at frequency

ω_{p h \min}^{W} = 38.96 s^{- 1}

and wavelength

λ_{p h \min}^{G l} = 6.59 cm

The viscosity of the liquid affects the capillary wave motion. An increase in viscosity leads to an increase in the wavelength at a constant frequency in the region of capillary waves. As a consequence, the values of the phase and group velocity increase. The characteristic values of fluid velocities for the wave component of periodic motion increase with increasing viscosity and shift to the region of lower frequencies (longer wavelengths). Taking into account the viscosity in the model makes it possible to calculate, in addition to the wave component, the ligament component of periodic motion in a liquid.

4. Reduction to Inviscid Fluid

In an inviscid exponentially stratified fluid, the equations of motion and boundary conditions (12), (13) are reduced

ν \to 0

z < ζ : {\begin{cases} ρ = ρ_{00} (r (z) + s (x, z, t)) \\ ρ \partial_{t} u + ρ (u \cdot \nabla) u = ρ g - \nabla P \\ \partial_{t} ρ + div (ρ u) = 0 \end{cases}

(64)

z < 0 : \partial_{t t} Δ ψ + N^{2} \exp (- z / Λ) \partial_{x x} ψ = 0

(65)

z = 0 : {\begin{cases} \partial_{t} ζ + \partial_{x} ψ = 0 \\ - ρ g ζ + ρ_{00} P + σ \partial_{x x} ζ = 0 \end{cases}

(66)

Dispersion relations in an inviscid liquid allow us to obtain only the wave component of periodic motion. The dispersion relations for the wave component are obtained using the limit transition

ν \to 0

from expressions (16) and (30) for high frequency (compared to the buoyancy frequency) and low frequency oscillations, respectively. For high-frequency wave disturbances of the free surface (

ω > N

) we obtain:

{\begin{cases} - N^{2} k_{x}^{2} + e^{z / Λ} (k_{x}^{2} - k_{z}^{2}) ω^{2} = 0 \\ g k_{x}^{2} - k_{z} ω^{2} + k_{x}^{4} γ = 0 \end{cases}

(67)

For low frequency waves (

ω < N

{\begin{cases} - N^{2} k_{x}^{2} + e^{z / Λ} (k_{x}^{2} + k_{z}^{2}) ω^{2} = 0 \\ - N^{2} k_{x}^{2} - e^{z / Λ} (k_{s}^{2} - k_{x}^{2}) ω^{2} = 0 \end{cases}

(68)

Thus, the approximation of an ideal fluid enables us to find solutions that describe the wave constituent of periodic motion in the region of gravitational waves quite well. In the region of capillary waves, the wave constituent calculated in the model of an ideal liquid underestimates the value of the wavelength at a given frequency. Note that in the model of an ideal fluid, it is impossible to obtain a solution for the ligament constituent of periodic motion along the free surface of the fluid. This makes the solution incomplete.

5. Discussion

For the first time the theory of singular perturbations was used to describe the propagation of two-dimensional periodic perturbations over the surface of a viscous exponentially stratified incompressible fluid occupying the lower half-space. Calculations were performed in a single formulation in a wide frequency range, which includes the propagation of capillary, gravitational and internal waves. The system of equations, which includes the equations of continuity, momentum transfer and density profile with depth, replacing the equation of state, is solved together with physically justified kinematic and dynamic conditions on a free surface, taking into account the effect of surface tension. The theory of singular perturbations with respect to the compatibility condition was applied to study the propagation of infinitesimal harmonic flows with a positive definite frequency and a complex wavenumber, which considers spatial attenuation.

A complete set of solutions of the obtained dispersion relation for infinitesimal periodic perturbations includes two inseparable constituents of periodic flows. The theory of regular perturbations determines the parameters of gravitational-capillary or infra-low-frequency waves with a frequency lower than the buoyancy frequency.

The accompanying fine constituents that are ligaments are calculated by the methods of the theory of singular perturbations. The dispersion equation designed for infinitesimal periodic perturbations describes the relationship between the components of the wave vector in periodic flows in a wide frequency range

10^{- 4} < ω < 10^{3} s^{- 1}

which includes infra-low frequency, gravitational and capillary waves. The result, represented by length- frequency dependence, is convenient for experimental determination. The dependences of phase and group wave velocities on frequency and wavelength have also been obtained.

The dependences of phase and group velocities of regular perturbations (i.e., waves in water and glycerin) on frequency and wavelength have been represented in graphs. Wave properties are significantly modified during the transition through buoyancy frequency. The group velocity of wave propagation tends to zero and the phase one tends to infinity with the approximation of wave and buoyancy frequencies. Viscosity has a significant influence on short waves when their length becomes compared to or less than capillary scale

λ < δ_{g}^{γ} = \sqrt{γ / g}

The wave parameters are calculated in an actual homogeneous liquid in the limit of the buoyancy frequency

N \to 0

as well as in a constant density approximation

ρ \equiv const, T_{b} = 0

taking into account the compatibility conditions. These calculations correspond to the given results.

Singular solutions describe ligaments (i.e., thin currents accompanying waves in a viscous stratified or homogeneous liquid). To compare them we have shown the graphs of the dependence of the wavelength and scale of the ligaments on the wave frequency in a viscous homogeneous liquid. The sizes of the ligaments differ by several orders of magnitude at low frequencies. The dependences of the phase and group velocities of waves and ligaments on the frequency and wavelength are also given.

The application of the theory of singular perturbations makes it possible to design complete solutions of the dispersion equation and the system of fundamental equations without additional hypotheses and limitations of the type of boundary layer approximation [33,34].

The evolving approach admits the extrapolation to the investigation of three-dimensional periodic perturbations. In this case, propagating surface waves accompany several types of ligaments, as well as internal waves in the thickness of a continuously stratified liquid [68,69].

In general cases, both waves and ligaments continuously interact with each other and generate new groups of flow constituents [75]. The influence of forcing together with the effects of nonlinearity and dissipation causes the evolution of the wave structure, which has different shapes at the phases of growth and attenuation.

In the low frequency range, the forcibly oscillating free surface of a stratified ocean can be a source of internal and inertial waves that transfer energy into the ocean. Internal waves, in turn, interact non-linearly with each other [76], generating new wave groups and ligaments (i.e., accompanying flows that cause the observed modulation of the surface waves) [77]).

6. Conclusions

The analysis of the linearized reduced version of the fundamental equations system by the methods of the theory of singular perturbations has been carried out with respect to the compatibility condition. It showed that the complete solutions describe waves propagating over the surface of a viscous stratified incompressible fluid, and small-scale constituents that are ligaments accompanying the waves. In extreme cases that are in limit of viscous homogeneous liquid, ideal stratified and ideal homogeneous liquid, the obtained dispersion relations for waves transfer to the widely known ones.

The experimental studies of the fine structure of surface waves in a continuously stratified liquid with high-resolution instruments that allow recording the influence of all constituents of periodic flows are of great interest. We should pay particular attention to the use of high-resolution optical methods for recording variations in the density gradient, and highly sensitive temperature or electrical conductivity sensors in a liquid with salt stratification.

Author Contributions

Conceptualization, Y.D.C. and A.A.O.; methodology, Y.D.C. and A.A.O.; software, A.A.O.; validation, Y.D.C. and A.A.O.; formal analysis, Y.D.C.; investigation, Y.D.C.; resources, Y.D.C.; data curation, Y.D.C.; writing—original draft preparation, Y.D.C. and A.A.O.; writing—review and editing, Y.D.C.; visualization, A.A.O.; supervision, Y.D.C.; project administration, Y.D.C.; funding acquisition, Y.D.C. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by the Russian Science Foundation (project 19-19-00598” Hydrodynamics and energetics of drops and droplet jets: formation, motion, break-up, interaction with the contact surface” https://rscf.ru/en/project/19-19-00598/. (accessed on 12 August 2022).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

Bernoulli, D. Hydrodynamica, Sive de Viribus et Motibus Fluidorum Commentarii, Opus Academicum; Sumptibus, J.R., Ed.; Argentorati: Dulseckeri, The Netherlands, 1738; 304p. [Google Scholar]
D’Alembert, J.-L.R. Réflexions sur la Cause Générale des Vents; David: Paris, France, 1747; 372p. [Google Scholar]
Euler, L. Principes généraux du mouvement des fluids. Mém. L’académie Des. Sci. Berl. 1757, 11, 274–315. [Google Scholar]
Feistel, R. Thermodynamic properties of seawater, ice and humid air: TEOS-10, before and beyond. Ocean Sci. 2018, 14, 471–502. [Google Scholar] [CrossRef]
Franklin, B. Behavior of oil on water. Letter to J. Pringle. In Experiments and Observations on Electricity; R. Cole: London, UK, 1769; pp. 142–144. [Google Scholar]
Lamb, H. Treatise on the Mathematical Theory of the on the Motion of Fluids, 6th ed.; Cambridge University Press: Cambridge, UK, 1879; 258p. [Google Scholar]
Rayleigh. Investigation of the character of the equilibrium of an incompressible heavy fluid of variable density. Proc. Lond. Math. Soc. 1882, 1–14, 170–177. [Google Scholar] [CrossRef]
Väisälä, V. Uber die Wirkung der Windschwankungen auf die Pilotbeoachtungen. Soc. Sci. Fenn. Commentat. Phys. Math. 1925, 2, 19–37. [Google Scholar]
Brunt, D. The period of simple vertical oscillations in the atmosphere. Q. J. R. Meteorol. Soc. 1927, 53, 30–32. [Google Scholar] [CrossRef]
Darrigol, O. Worlds of Flow: A History of Hydrodynamics from the Bernoullis to Prandtl; Oxford University Press: Oxford, UK, 2005; 356p. [Google Scholar]
Basile, R.; De Serio, F. Flow Field around a Vertical Cylinder in Presence of Long Waves: An Experimental Study. Water 2022, 14, 1945. [Google Scholar] [CrossRef]
Gerstner, F.J. Theorie der Wellen. Abhandlunger der Königlichen Böhmischen Geselschaft der Wissenschaften, Prague. Ann. der Phys. 1809, 32, 412–445. [Google Scholar] [CrossRef]
Rankine, W.J.M. On the exact form of waves near the surface of deep water. Philos. Trans. R. Soc. 1863, 153, 127–138. [Google Scholar]
Rayleigh. On waves. Phil. Mag. Ser. 1876, 5, 257–279. [CrossRef]
Kochin, N.E.; Kibel, I.A.; Rose, N.V. Theoretical Hydromechanics; Part 1; OGIZ. Gostekhizdat: Leningrad-Moscow, USSR, 1948; 728p. (In Russian) [Google Scholar]
Henry, D. On Gerstner’s Water Wave. J. Nonlinear Math. Phys. 2008, 15, 87–95. [Google Scholar] [CrossRef]
Ionescu-Kruse, D. On the particle paths and the stagnation points in small-amplitude deep-water. J. Math. Fluid Mech. 2013, 15, 41–54. [Google Scholar] [CrossRef]
Russell, J.S. Report on Waves. Report of the Fourteenth Meeting of the British Association for the Advancement of Science; Murray, J., Street, A., Eds.; Scientific Research: London, UK, 1845; pp. 311–391. [Google Scholar]
Boussinesq, J. Theorie de I’intumescence Liquid, Appleteonde Solitaire au de Translation, se Propageantdansun Canal Rectangulaire. Les Comptes Rendus De L’académie Des Sci. 1871, 72, 755–759. [Google Scholar]
Korteweg, D.J.; de Vries, G. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1895, 39, 422–443. [Google Scholar] [CrossRef]
Whitham, G.B. Linear and Nonlinear Waves; Wiley Interscience: New York, NY, USA, 1999; 660p. [Google Scholar]
Stokes, G.G. On the theory of oscillatory waves. Trans. Camb. Phil. Soc. 1847, 8, 441–455. [Google Scholar]
Stokes, G.G. Supplement to a paper on the theory of oscillatory waves. In Mathematical and Physical Papers; CUP: Cambridge, UK, 1880; Volume I, pp. 314–326. [Google Scholar]
Kistovich, A.V.; Chashechkin, Y.D. Propagating stationary surface potential waves in a deep ideal fluid. Water Resour. 2018, 45, 719–727. [Google Scholar] [CrossRef]
Basset, B. A Treatise on Hydrodynamics with Numerous Examples; Deighton, Bell and Co.: Cambridge, UK, 1888; Volume II, 368p. [Google Scholar]
Hough, S.S. On the Influence of Viscosity on Waves and Currents. Proc. Lond. Math. Soc. 1896, s1–28, 264–288. [Google Scholar] [CrossRef]
Harrison, W.J. The Influence of Viscosity on the Oscillations of Superposed Fluids. Proc. Lond. Math. Soc. 1908, s2–6, 96–405. [Google Scholar] [CrossRef]
Biesel, F. Calcul de l’amortissement d’une houle dans un liquide visqueux de profondeur finie. La Houille Blanche 1949, 4, 630–634. [Google Scholar] [CrossRef]
Tyvand, P.A. A note on gravity waves in a viscous liquid with surface tension. Z. Angew. Math. Phys. 1984, 35, 592–597. [Google Scholar] [CrossRef]
Antuono, M.; Colagross, A. The damping of viscous gravity waves. Wave Motion 2013, 50, 197–209. [Google Scholar] [CrossRef]
Prandtl, L. Über Flüssigkeitsbewegung bei sehr kleiner Reibung. In Proceedings of the Verhandlungen des Dritten Internationalen Mathematiker-Kongresses, Heidelberg, Germany, 8–13 August 1904; Teubner: Leipzig, Germany, 1905; pp. 485–491. [Google Scholar]
Schlichting, H. Boundary Layer Theory; McGraw Hill Co.: New York, NY, USA, 1955; 812p. [Google Scholar]
Longuet-Higgins, M.S. Mass transport in water waves. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Sci. 1953, 245, 535–581. [Google Scholar] [CrossRef]
Longuet-Higgins, M.S. Mass transport in the boundary layer at a free oscillating surface. J. Fluid Mech. 1960, 8, 293–306. [Google Scholar] [CrossRef]
Liu, A.; Davis, S. Viscous attenuation of mean drift in water waves. J. Fluid Mech. 1977, 81, 63–84. [Google Scholar] [CrossRef]
Robertson, S.; Rousseaux, G. Viscous dissipation of surface waves and its relevance to analogue gravity experiments. arXiv 2018, arXiv:1706.05255. [Google Scholar]
Dore, B. Mass transport in layered fluid systems. J. Fluid Mech. 1970, 40, 113–126. [Google Scholar] [CrossRef]
Thomson, W. Hydrokinetic solutions and observations. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1871, 42, 362–377. [Google Scholar] [CrossRef]
Kistovich, A.V.; Chashechkin, Y.D. Dynamics of gravity-capillary waves on the surface of a nonuniformly heated fluid. Izv. Atmos. Ocean. Phys. 2007, 43, 95–102. [Google Scholar] [CrossRef]
Chen, X.; Duan, W.-Y.; Lu, D.Q. Gravity waves with effect of surface tension and fluid viscosity. J. Hydrodyn. 2006, 18, 171–176. [Google Scholar] [CrossRef]
Caligny, A.F.H. Expériences sur les mouvements des molécules liquides des ondes courantes, considérées dans leur mode d’action sur la marche des navires. Compt. Rendu Acad. Sci. Paris 1878, 87, 1019–1023. [Google Scholar]
Bagnold, R.A. Sand movement by waves: Some small-scale experiments with sand of very low density. J. Inst. Civ. Eng. 1947, 27, 447–469. [Google Scholar] [CrossRef]
Van Dorn, W.G. Boundary dissipation of oscillatory waves. J. Fluid Mech. 1966, 24, 769–779. [Google Scholar] [CrossRef]
Kalinichenko, V.A.; Chashechkin, Y.D. Structurization and restructurization of a homogeneous suspension in a standing wave field. Fluid Dyn. 2012, 47, 778–788. [Google Scholar] [CrossRef]
Chashechkin, Y.D. Transfer of the Substance of a Colored Drop in a Liquid Layer with Travelling Plane Gravity–Capillary Waves. Izv. Atmos. Ocean. Phys. 2022, 58, 188–197. [Google Scholar] [CrossRef]
Chashechkin, Y.D.; Ilinykh, A.Y. Drop decay into individual fibers at the boundary of the contact area with a target fluid. Dokl. Phys. 2021, 66, 101–105. [Google Scholar] [CrossRef]
Toffoli, A.; Bitner-Gregersen, E.M. Types of ocean surface waves, wave classification. In Encyclopedia of Maritime and Offshore Engineering; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2017. [Google Scholar] [CrossRef]
Hung, L.-P.; Tsa, W.-T. The Formation of parasitic capillary ripples on gravity–capillary waves and the underlying vortical structures. J. Phys. Oceanogr. 2009, 39, 263–289. [Google Scholar] [CrossRef]
Kinsman, B. Wind Waves; Prentice-Hall: Englewood Cliffs, NJ, USA, 1965; 676p. [Google Scholar]
Holthuijsen, L.H. Waves in Oceanic and Coastal Waters; Cambridge University Press: Cambridge, UK, 2007; 388p. [Google Scholar]
Chashechkin, Y.D. Packets of capillary and acoustic waves of drop impact. Her. Bauman Mosc. State Technol. Univ. Ser. Nat. Sci. 2021, 1, 73–92. [Google Scholar] [CrossRef]
Chashechkin, Y.D.; Prohorov, V.E. Evolution of the structure of acoustic signals caused by the impact of a falling drop on a liquid. Acoust. Phys. 2020, 66, 362–374. [Google Scholar] [CrossRef]
Pumphrey, H.C.; Crum, L.A.; Bjorno, L. Underwater sound produced by individual drop impacts and rainfall. J. Acoust. Soc. Am. 1989, 85, 1518–1526. [Google Scholar] [CrossRef]
Bjorno, L. Underwater rain noise: Sources, spectra and interpretations. J. De Phys. IV Proc. EDP Sci. 1994, 4, C5/1023–C5/1030. [Google Scholar] [CrossRef]
Zhang, X. Enhanced dissipation of short gravity and gravity capillary waves due to parasitic capillaries. Phys. Fluids 2002, 14, L81. [Google Scholar] [CrossRef]
Sergievskaya, I.A.; Ermakov, S.A.; Ermoshkin, A.V.; Kapustin, I.A.; Molkov, A.A.; Danilicheva, O.A.; Shomina, O.V. Modulation of dual-polarized x-band radar backscatter due to long wind waves. Remote Sens. 2019, 11, 423. [Google Scholar] [CrossRef]
Slavchov, R.; Peychev, B.; Ismail, A.S. Characterization of capillary waves: A review and a new optical method. Phys. Fluids 2021, 33, 101303. [Google Scholar] [CrossRef]
Falcon, E.; Mordant, N. Experiments in surface gravity-capillary wave turbulence. Annu. Rev. Fluid Mech. 2022, 54, 1–27. [Google Scholar] [CrossRef]
Berhanu, M. Impact of the dissipation on the nonlinear interactions and turbulence of gravity-capillary waves. Fluids 2022, 7, 137. [Google Scholar] [CrossRef]
Thorpe, S.A. An Introduction to Ocean Turbulence; Cambridge University Press: Bangor, UK, 2012; 244p. [Google Scholar]
McWilliams, J.C.; François Colas, F.; Molemaker, M.J. Cold filamentary intensification and oceanic surface convergence lines. Geophys. Res. Lett. 2009, 36, L18602. [Google Scholar] [CrossRef]
Landau, L.D.; Lifshitz, E.M. Fluid Mechanics. V.6. Course of Theoretical Physics; Pergamon Press: Oxford, UK, 1987; 560p. [Google Scholar]
Müller, P. The Equations of Oceanic Motions; Cambridge University Press: Cambridge, UK, 2006; 291p. [Google Scholar]
Vallis, G.K. Atmospheric and Oceanic Fluid Dynamics; Cambridge University Press: Cambridge, UK, 2006; p. 745. [Google Scholar]
Nayfeh, A.H. Introduction to Perturbation Technique; John Wiley & Sons: New York, NY, USA, 1993; 536p. [Google Scholar]
Chashechkin, Y.D. Foundations of engineering mathematics applied for fluid flows. Axioms 2021, 10, 286. [Google Scholar] [CrossRef]
Lighthill, J. Waves in Fluids; Cambridge University Press: Cambridge, UK, 1978; 524p. [Google Scholar]
Bardakov, R.N.; Vasil’ev, A.Y.; Chashechkin, Y.D. Calculation and measurement of conical beams of three-dimensional periodic internal waves excited by a vertically oscillating piston. Fluid Dyn. 2007, 42, 612–626. [Google Scholar] [CrossRef]
Kistovich, A.V.; Chashechkin, Y.D. Fine structure of a conical beam of periodical internal waves in a stratified fluid. Izv. Atmos. Ocean. Phys. 2014, 50, 103–110. [Google Scholar] [CrossRef]
Zagumennyi, Y.V.; Chashechkin, Y.D. Pattern of unsteady vortex flow around plate under a zero angle of attack (two-dimensional problem). Fluid Dyn. 2016, 51, 53–70. [Google Scholar] [CrossRef]
Prandtl, L. Führer Durch die Strömungslehre; Verlagskatalog Von Friedr; Vieweg & Sohn in Braunschweig: Berlin, Germany, 1942; 337p. [Google Scholar]
Makarov, S.A.; Neklyudov, V.I.; Chashechkin, Y.D. Spatial structure of twodimensional monochromatic internal-wave beams in an exponentially stratified liquid. Izv. Atmos. Ocean. Phys. 1990, 26, 548–554. [Google Scholar]
Kistovich, Y.V.; Chashechkin, Y.D. An exact solution of a linearized problem of the radiation of monochromatic internal waves in a viscous fluid. J. Appl. Maths Mech. 1999, 63, 587–594. [Google Scholar] [CrossRef]
Kistovich, Y.V.; Chashechkin, Y.D. Linear theory of the propagation of internal wave beams in an arbitrarily stratified liquid. J. Appl. Mech. Technol. Phys. 1998, 39, 729–737. [Google Scholar] [CrossRef]
Chashechkin, Y.D. Conventional partial and new complete solutions of the fundamental equations of fluid mechanics in the problem of periodic internal waves with accompanying ligaments generation. Mathematics 2021, 9, 586. [Google Scholar] [CrossRef]
Veeresha, P.; Baskonus, H.M.; Gao, W. Strong Interacting internal waves in rotating ocean: Novel fractional approach. Axioms 2021, 10, 123. [Google Scholar] [CrossRef]
Lenain, L.; Pizzo, N. Modulation of surface gravity waves by internal waves. J. Phys. Oceanogr. 2021, 51, 2735–2748. [Google Scholar] [CrossRef]

Figure 1. The density–depth dependencies for exponential stratification (solid line) and for linear stratification (dotted line) for strongly stratified fluid (

N = 1 s^{- 1}

ν = 0.01 St

σ = 72 dyn / cm,

ρ_{00} = 1.0 g / {cm}^{3})

Figure 1. The density–depth dependencies for exponential stratification (solid line) and for linear stratification (dotted line) for strongly stratified fluid (

N = 1 s^{- 1}

ν = 0.01 St

σ = 72 dyn / cm,

ρ_{00} = 1.0 g / {cm}^{3})

Figure 2. The dependence of the wavelength on the frequency of wave motion in a viscous exponentially stratified fluid for a periodic solution. The curves indicated by the letter

(W)

are constructed for a liquid with

ν = 0.01 St

σ = 72 dyn / cm,

ρ_{00} = 1 g / {cm}^{3}

, and by the letter

(G l)

—for a liquid with glycerin parameters

(ν = 11.746 St,

ρ_{00} = 1.26 g / {cm}^{3}

σ = 64.7 dyn / cm

). The numbers indicate a different degree of stratification. Index (1) corresponds to a weak pycnocline

N = 0.001 s^{- 1}

, index (2) to a weakly stratified fluid

N = 0.01 s^{- 1}

, and index (3) to a strongly stratified fluid

N = 1 s^{- 1}

for water and glycerin.

Figure 2. The dependence of the wavelength on the frequency of wave motion in a viscous exponentially stratified fluid for a periodic solution. The curves indicated by the letter

(W)

are constructed for a liquid with

ν = 0.01 St

σ = 72 dyn / cm,

ρ_{00} = 1 g / {cm}^{3}

, and by the letter

(G l)

—for a liquid with glycerin parameters

(ν = 11.746 St,

ρ_{00} = 1.26 g / {cm}^{3}

σ = 64.7 dyn / cm

). The numbers indicate a different degree of stratification. Index (1) corresponds to a weak pycnocline

N = 0.001 s^{- 1}

, index (2) to a weakly stratified fluid

N = 0.01 s^{- 1}

, and index (3) to a strongly stratified fluid

N = 1 s^{- 1}

for water and glycerin.

Figure 3. Dependences of the phase (dashed lines) and group (solid lines) velocities on the frequency of wave motion (a) and on the wavelength (b) in a viscous exponentially stratified fluid for a periodic solution. The curves indicated by the letter

(W)

are constructed for a liquid with water parameters, and by the letter

(G l)

—for a liquid with glycerin parameters. The numbers indicate a different degree of stratification. Indexes (1) and (6) correspond to a weak pycnocline

N = 0.001 s^{- 1}

, indexes (2) and (5) correspond to a weakly stratified liquid

N = 0.01 s^{- 1}

and indexes (3) and (4) correspond to a strongly stratified liquid

N = 1 s^{- 1}

(W)

are constructed for a liquid with water parameters, and by the letter

(G l)

—for a liquid with glycerin parameters. The numbers indicate a different degree of stratification. Indexes (1) and (6) correspond to a weak pycnocline

N = 0.001 s^{- 1}

, indexes (2) and (5) correspond to a weakly stratified liquid

N = 0.01 s^{- 1}

and indexes (3) and (4) correspond to a strongly stratified liquid

N = 1 s^{- 1}

Figure 4. Wavelength dependences of the wave solution (a) and the ligament solution (b) on the frequency in a viscous homogeneous fluid. The curves indicated by the letter

(W)

are constructed for a liquid with water parameters, and by the letter

(G l)

—for a liquid with glycerin parameters.

Figure 4. Wavelength dependences of the wave solution (a) and the ligament solution (b) on the frequency in a viscous homogeneous fluid. The curves indicated by the letter

(W)

are constructed for a liquid with water parameters, and by the letter

(G l)

—for a liquid with glycerin parameters.

Figure 5. Dependences of the phase (dashed lines) and group (solid lines) velocities on the frequency of the wave solution (a) and the ligament solution (b) on the frequency in a viscous homogeneous fluid. The curves indicated by the letter

(W)

are constructed for a liquid with water parameters, and by the letter

(G l)

for a liquid with glycerin parameters.

(W)

are constructed for a liquid with water parameters, and by the letter

(G l)

for a liquid with glycerin parameters.

Figure 6. Dependences of the phase (dashed lines) and group (solid lines) velocities on the wavelength of the wave solution (a) and the ligament solution (b) on the frequency in a viscous homogeneous fluid. The curves indicated by the letter

(W)

are constructed for a liquid with water parameters, and by the letter

(G l)

for a liquid with glycerin parameters.

(W)

are constructed for a liquid with water parameters, and by the letter

(G l)

for a liquid with glycerin parameters.

Table 1. Values of characteristic physical quantities.

Parameter	Fluids
	Stratified (SF)		Homogeneous (HF)
	Strongly	Weakly	Potentially	Actually
$Buoyancy frequency N, s^{- 1}$	1	0.01	0.00001	0.0
$Period T_{b}$	6.28 s	10.5 min	7.3 days	$\infty$
$Scale of stratification Λ$	9.8 m	100 km	$10^{8}$ km	$\infty$
$Viscous wave scale L_{N}^{g ν} = \sqrt[3]{g ν} / N$ , cm	2.14	200	$2 \cdot 10^{5}$	$\infty$
Stokes microscale $δ_{N}^{ν} = \sqrt{ν / N}$ , cm	0.1	1.0	30	$\infty$

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Chashechkin, Y.D.; Ochirov, A.A. Periodic Waves and Ligaments on the Surface of a Viscous Exponentially Stratified Fluid in a Uniform Gravity Field. Axioms 2022, 11, 402. https://doi.org/10.3390/axioms11080402

AMA Style

Chashechkin YD, Ochirov AA. Periodic Waves and Ligaments on the Surface of a Viscous Exponentially Stratified Fluid in a Uniform Gravity Field. Axioms. 2022; 11(8):402. https://doi.org/10.3390/axioms11080402

Chicago/Turabian Style

Chashechkin, Yuli D., and Artem A. Ochirov. 2022. "Periodic Waves and Ligaments on the Surface of a Viscous Exponentially Stratified Fluid in a Uniform Gravity Field" Axioms 11, no. 8: 402. https://doi.org/10.3390/axioms11080402

APA Style

Chashechkin, Y. D., & Ochirov, A. A. (2022). Periodic Waves and Ligaments on the Surface of a Viscous Exponentially Stratified Fluid in a Uniform Gravity Field. Axioms, 11(8), 402. https://doi.org/10.3390/axioms11080402

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

Article Menu

Periodic Waves and Ligaments on the Surface of a Viscous Exponentially Stratified Fluid in a Uniform Gravity Field

Abstract

1. Introduction

2. Periodic Flows in a Viscous Exponentially Stratified Fluid

2.1. Equations of Periodic Flows and Dispersion Relations for Plane Infinitesimal Waves

2.2. Solution of the Dispersion Equation

2.3. Low Frequency Waves

2.4. Periodic Flows on the Surface of a Viscous Exponentially Stratified Liquid

3. Reduction to Approximation of Actually Homogeneous Fluid

4. Reduction to Inviscid Fluid

5. Discussion

6. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI