Next Article in Journal
Estimations of Covering Functionals of Convex Bodies Based on Relaxation Algorithm
Next Article in Special Issue
Semi-Analytical Methods in the Problem of Deformation of a Fluid Strip
Previous Article in Journal
A Stochastic Convergence Result for the Nelder–Mead Simplex Method
Previous Article in Special Issue
Significance of Nanoparticle Radius and Gravity Modulation on Dynamics of Nanofluid over Stretched Surface via Finite Element Simulation: The Case of Water-Based Copper Nanoparticles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Dynamics of Perturbations of a Plane Potential Fluid Flow: Nonlocal Generalization of the Hopf Equation

by
Nikolay M. Zubarev
1,2
1
Institute of Electrophysics, UB RAS, 620016 Yekaterinburg, Russia
2
Lebedev Physical Institute, RAS, 119991 Moscow, Russia
Mathematics 2023, 11(9), 1999; https://doi.org/10.3390/math11091999
Submission received: 27 March 2023 / Revised: 18 April 2023 / Accepted: 21 April 2023 / Published: 23 April 2023
(This article belongs to the Special Issue Numerical and Analytical Study of Fluid Dynamics)

Abstract

:
In this paper, we analytically study the two-dimensional unsteady irrotational flow of an ideal incompressible fluid in a half-plane whose boundary is assumed to be a linear sink. It is shown that the nonlinear evolution of perturbations of the initial uniform flow is described by a one-dimensional integro-differential equation, which can be considered as a nonlocal generalization of the Hopf equation. This equation can be reduced to a system of ordinary differential equations (ODEs) in the cases of spatially localized or spatially periodic perturbations of the velocity field. In the first case, ODEs describe the motion of a system of interacting virtual point vortex-sinks/sources outside the flow domain. In the second case, ODEs describe the evolution of a finite number of harmonics of the velocity field distribution; this is possible due to the revealed property of the new equation that the interaction of initial harmonics does not lead to generation of new ones. The revealed reductions made it possible to effectively study the nonlinear evolution of the system, in particular, to describe the effect of nonlinearity on the relaxation of velocity field perturbations. It is shown that nonlinearity can significantly reduce the relaxation rate by more than 1.5 times.

1. Introduction

The Hopf (Burgers-Hopf) equation [1,2]
u t + u u x = 0
is the simplest model in fluid mechanics. It has the meaning of the transport equation for the fluid velocity u ( x , t ) , where t is the time, x is the coordinate, u t and u x are the corresponding partial derivatives. This equation arises, for example, when studying the one-dimensional (1D) motion of a gas of noninteracting particles: the 1D Euler equation in the absence of any forces. The nonlinear term of the Hopf equation is a characteristic nonlinearity that occurs in almost any problem of fluid mechanics.
The Hopf equation is easily integrated by the method of characteristics. Its solution is defined implicitly by the formula u = u 0 ( x u t ) where u 0 ( x ) = u ( x , 0 ) . The solution exists only for a finite time, being destroyed as a result of the wave breaking. The Hopf equation corresponds to the dispersionless limit of the Burgers, Korteweg-de Vries, and Benjamin–Ono equations. The presence of linear dispersion terms in these equations suppresses the wave breaking process, ensuring the existence of solutions at all times t > 0 .
In the two-dimensional (2D) case, where there are two x - and y -components of the velocity, u and v , respectively, the complex analogue of the Hopf equation
U t + U U z = 0
can be applied to describe unsteady potential flows of an incompressible fluid with a free surface [3,4,5]. Here U = u i v is the complex velocity, which is an analytic function of the complex variable z = x + i y (see also works [6,7,8] and references therein). Equation (2) turns out to be compatible with the 2D Euler equations
u t + u u x + v u y = p x / ρ
v t + u v x + v v y = p y / ρ
if the pressure p is related to the velocity component v by the formula p = ρ v 2 (here ρ is the fluid density). The free fluid boundary corresponds to the condition p = 0 [5]. The complete integrability of the Hopf Equation (2) made it possible to construct a number of nontrivial examples of 2D flows of a fluid with a free surface in the absence of the gravitational force and capillarity. Solutions of the complex Hopf equation are destroyed in a finite time; the analyticity of the function U is violated in the flow domain. This leads to the formation of various types of singularities (depending on initial conditions) on the free surface such as cusps, droplets, and bubbles [5,9,10]. Note that works [11,12] study plane potential flows with a free boundary described by the modified complex Hopf equation
U t + e i γ U U z = 0 ,
where γ is a real constant. Equation (3) can be reduced to the form (2) by replacing U e i γ U .
In [13,14], when analyzing the 2D inertial motion of a deep fluid with a free surface in the absence of external forces and capillarity, the following integro-differential equation has been derived in the weakly nonlinear approximation:
u t + u u x H [ u ]   H [ u x   ] = 0 ,
where H stands for the Hilbert transform with respect to x (see Section 3). The nonlocal Equation (4) is transformed into the complex Hopf Equation (2) by introducing the complex velocity U = u i H [ u ] .
The equation
u t + u u x H [ u ]   H [ u x   ] = V H [ u x   ]
with the nonlinear part coinciding with that of Equation (4), where V is a constant, arose in the study of the instability of the boundary of a conducting liquid in an external electric field [15], as well as of the Kelvin–Helmholtz instability of the interface between two fluids with strongly different densities [16] (compare with [17]). In terms of the complex velocity, it rewrites as
U t + U U x = i V U x ,
which differs from the Hopf Equation (2) by the linear term on the right-hand side (this term can be eliminated by the simple change of variables U U i V ).
Both Equations (4) and (5), due to their connection with the Hopf equation, describe the formation of singularities at the fluid boundary in a finite time (to be precise, the formation of Moore’s singularities [18]).
In the present paper, when considering a certain 2D flow in the half-plane < x < and y 0 , we arrive at the new nonlinear integro-differential equation
u t + u u x + H [ u ]   H [ u x   ] = V H [ u x   ] .
It contains the same terms as Equation (5) but differs from it in the signs of the second nonlinear term and of the linear term on the right-hand side. As a result, it cannot be reduced to the complex Hopf equation and its modifications. Equation (6) has the same symmetries as the Hopf equation, however, as will be shown, it has fundamentally different properties that are closer to those of the Laplacian growth equation [19,20]. This leads to the need to use different methods of integrating the equation and to other properties of its solutions. The solutions remain smooth at all times t > 0 .
The plan of the paper is as follows: Section 2 contains the formulation of the problem where we write down the equations describing the two-dimensional unsteady flow of an incompressible fluid in a half-plane, the boundary of which corresponds to a linear sink. A physical interpretation of such a flow is given. Section 3 derives the key integro-differential Equation (6), which can be considered as a nonlocal generalization of the Hopf equation. Section 4 and Section 5 show that this equation can be reduced to a finite number of ordinary differential equations (ODEs) for the cases of either spatially localized or spatially periodic perturbations of the initial uniform flow. In the first case, ODEs describe the motion of a system of interacting virtual point vortex-sinks/sources outside the flow domain. In the second case, ODEs describe the evolution of a finite number of harmonics of the velocity field distribution. Such a reduction is possible due to the following revealed property of Equation (6): the interaction of initial harmonics does not lead to the generation of new ones. It is demonstrated that the influence of the nonlinearities of Equation (6) leads to a rather complex dynamics of velocity field perturbations. The character of perturbation relaxation described by exact solutions of the problem and by solutions obtained in the linear approximation is compared. It is shown that nonlinearity significantly decreases the relaxation rate. Finally, Section 6 contains some concluding remarks.

2. Formulation of the Problem and Initial Equations

We will consider the following 2D boundary-value problem:
Φ x x + Φ y y = 0 , 0 < y < ,
Φ t + 1 2   ( Φ x ) 2 + 1 2   ( Φ y ) 2 = 1 2   V 2 , y = 0 ,
Φ y V , y ,
where V > 0 is some constant. These equations describe the irrotational flow of an incompressible liquid in the half-plane < x < and y > 0 with the velocity potential Φ ( x , y , t ) (we have { u , v } = { Φ x , Φ y } for the x - and y -components of the velocity vector). The potential Φ satisfies the Laplace Equation (7), which is compatible with the non-stationary Bernoulli equation
Φ t + 1 2   ( Φ x ) 2 + 1 2   ( Φ y ) 2 = p ρ + Γ ( t ) , 0 < y < ,
where ρ = const is the fluid density, p ( x , y , t ) is the pressure, and Γ is some function of time arising from the spatial integration of the Euler equations. According to Condition (9), the flow is uniform at infinity y and directed along the y axis in the negative direction; V is the absolute value of the velocity of this flow. The fluid flows out of the domain < x < and y > 0 , crossing its linear boundary y = 0 . We suppose that the pressure along this boundary is constant, p | y = 0 = p 0 = const . In view of Equation (10), this boundary condition can be written as
Φ t + 1 2   ( Φ x ) 2 + 1 2   ( Φ y ) 2 = p 0 ρ + Γ ( t ) , y = 0 .
Using some arbitrariness in determining the velocity potential, we choose Γ = p 0 / ρ + V 2 / 2 = const , which reduces (11) to the Form (8). It is clear that the uniform flow { u , v } = { 0 , V } gives the trivial, unperturbed solution to the problem. In terms of the velocity potential, this solution is
Φ = V y .
We will be interested in the dynamics of perturbations of this basic flow.
The formulated problem can be considered as describing an upstream flow before a linear sink on the straight line < x < and y = 0 . For clarity, one can imagine the flow in front of a waterfall. In the upstream region y > 0 , i.e., before the waterfall, we have a plane potential flow of a fluid placed between two horizontal semi-infinite plates; the depth of the fluid layer (it corresponds to the vertical direction in three-dimensional consideration) is much less than the characteristic scale of the flow in the plane { x , y } . The fluid flows across the plateau y > 0 towards the cliff edge y = 0 , where the plates end, and then it falls freely (see, for example, [21,22,23]). If we again assume that the flow scale in the plane { x , y } is dominant, then we can admit that the pressure is unloaded at the boundary y = 0 to its atmospheric value p 0 , which is embedded in Equation (8). Thus, Equations (7)–(9) can be interpreted as a toy model of 2D flow in front of a linear sink/waterfall.
In view of the indicated interpretation of the problem (7)–(9), it must be supplemented with the condition that the liquid cannot flow into the domain y > 0 through its boundary y = 0 ; the y -component of the velocity cannot be positive at y = 0 , i.e., v | y = 0 0 . In terms of the velocity potential Φ , this corresponds to the kinematic condition in the form of the inequality
Φ y 0 , y = 0 .
Indeed, the fluid falls freely after crossing the boundary y = 0 , so that it cannot subsequently return to the region y > 0 . This immediately leads to Condition (13) for the velocity component normal to the flow boundary.

3. Derivation of the Key Equation

Let us reduce the original 2D problem to a single 1D evolution equation. Considering a uniform flow with the velocity potential Φ = V y as unperturbed, we introduce the potential perturbation as φ = Φ + V y . Substituting it into (7)–(9), we obtain
φ x x + φ y y = 0 , 0 < y < ,
φ t + 1 2   ( φ x ) 2 + 1 2   ( φ y ) 2 = V φ y , y = 0 ,
φ y 0 , y .
It is necessary to eliminate the variable y from these equations. Note that the solution of the two-dimensional Laplace Equation (14) at y > 0 in view of the condition (16) is given by
φ ( x , y , t ) = 1 π +   y   ψ ( x , t ) y 2 + ( x x ) 2   d x ,
where ψ ( x , t ) φ | y = 0 is the function that specifies the perturbation of the velocity potential at the boundary y = 0 . Differentiating this expression with respect to y , we find φ y = H [ ψ x   ] at y = 0 , where H is the Hilbert transform with respect to x (the “P.V.” denotes the principal-value integral):
H [ f ( x ) ] = 1 π   P . V .     + f ( x ) x x   d x .
As a result, we obtain from (15) the one-dimensional nonlinear evolution equation for the function ψ ( x , t ) , which will be the subject of our further consideration:
  ψ t + 1 2 ( ψ x ) 2 + 1 2 ( H [ ψ x   ] ) 2 = V H [ ψ x   ] .
To clarify the notation used in this equation, we can rewrite it in the equivalent form
  ψ t + 1 2 ( ψ x ) 2 + 1 2 ( K [ ψ   ] ) 2 = V K [ ψ   ] ,
where K x H is the integral operator with the Fourier transform equal to | k | , i.e., K [ e i k x   ] =     | k | e i k x (see [14,15,16]).
The obvious difficulty in studying Equation (17) is that it is no longer a partial differential equation, similar to the original Equations (7)–(9), but an integro-differential equation. We also recall that we are not interested in all solutions of Equation (17), but only in those that satisfy the Condition (13), which, in terms of the function ψ , can be rewritten as
H [ ψ x   ] V .
By differentiating Equation (17) with respect to x , we arrive at the following equation for the x -component of the velocity at the boundary y = 0 , u | y = 0 = ψ x :
u t + u u x + H [ u ]   H [ u x   ] = V H [ u x   ] .
It can be considered as a nonlocal generalization of the Hopf equation (see Section 1 for details).
In the linear approximation, Equation (17) turns into the compact equation
  ψ t = V H [ ψ x   ] .
The dispersion relation for (20) can be obtained by substituting ψ exp ( i k x i ω t ) , where ω is the frequency and k is the wave number. Taking into account the property H [ e i k x   ] = i   sgn ( k ) e i k x of the Hilbert transform, we find the simple expression
  ω = i V | k | .
According to this expression, perturbations decay with the characteristic time 1 / ( V | k | ) . Such relaxation is obviously due to the fact that perturbations of the velocity field are carried downstream with the velocity V , i.e., the linear stability of the unperturbed flow (12) is inherent in the formulation of the problem itself. The relaxation character is determined by the fact that Equation (20) is actually elliptical, that is, it corresponds to one of the branches of the dispersion law   ω 2 = V 2 k 2 of the equation   ψ t t + V 2 ψ x x = 0 . Note the difference in the character of relaxation processes from those described by the parabolic and hyperbolic heat diffusion equations [24].
The present work is devoted to a large extent to the consideration of how the nonlinear terms of Equation (17) will affect the character of the relaxation of perturbations. First of all, this is important for situations where the perturbations are not small and, therefore, the approximation (20) together with the dispersion law (21) are not applicable. Below we will pay considerable attention to situations where, at the initial moment of time t = 0 at some point, for example, at x = 0 , the y -component of the velocity at the boundary (for convenience we denote it as v b ( x , t ) v | y = 0 = V + H [ ψ x   ] ) is equal to zero: v b ( 0 , 0 ) = 0 . For instance, one can take the initial velocity distribution in the form
  v b ( x , 0 ) = V + 4 a 2 V V 2 x 2 + 4 a 2 ,
for which   v b ( x , 0 ) V at | x | and   v b ( 0 , 0 ) = 0 (here   a is some constant). Such a situation can be considered as the presence of an obstacle in the downflow, which is removed at the moment t = 0 . Then we need to describe the relaxation of the flow to the unperturbed state   { u , v } = { 0 , V } taking into account the influence of nonlinear terms in Equation (17). For the initial velocity distribution (22) (the perturbation of the   y -component of the velocity at x = 0 is equal to the velocity of the basic flow, taken with the opposite sign), the contribution of linear and nonlinear terms will be comparable.
Let us represent the perturbation of the velocity potential at the boundary in the form
ψ = ψ + + ψ ,
where ψ ± are analytic continuations of the function ψ into the upper and lower half-planes, respectively, of the complex variable x . The functions ψ ± can be written as the actions of the projectors P ± = ( 1 i H ) / 2 on the function ψ : ψ ± = P ± [ ψ ] . Due to the property H [ ψ ± ] = ± i ψ ± of the Hilbert transform, Equation (17) takes the form
  ψ t + +   ψ t + 2 ψ x + ψ x = i V ψ x + i V ψ x .
It can be separated into a pair of equations for parts that are analytic in the upper and lower half-planes of the x variable:
  ψ t ± + 2 P ± [ ψ x + ψ x ] = ± i V ψ x ± .
Since the functions ψ ± are complex conjugates, these equations are also complex conjugates. Then it suffices to consider, for instance, only the equation for the function ψ + , namely
  ψ t + + 2 P + [ ψ x + ψ x ] = i V ψ x + .
The nonlinear part of Equation (23) resembles the well-known Laplacian growth equation (LGE),
  Im ( ψ t + ψ x ) = 1 ,
which is used for describing the 2D motion of the boundary between two liquids with noticeably different viscosities (Hele-Shaw flow) [25,26,27], evolution of the free surface of a liquid in a gravitational field in the high-Jacobian approximation [28,29], electrostatic aggregation in two dimensions [30], electrohydrodynamic and quantum Kelvin–Helmholtz instabilities of the free surface of liquid helium [31,32], and so on. It is integrable in the sense of the existence of an infinitely large number of exact particular solutions and relation to the dispersionless limit of the integrable Toda hierarchy [33]. LGE can be written in the form
  ψ t + ψ x ψ t ψ x + = 2 i ,
or, in order to emphasize the similarity with Equation (23), in the form
  ψ t + + P + [ ψ t + ψ x ψ t ψ x + ] = i ψ x + ,
where we have made the changes ψ ± = ψ ± ± i t + x and applied the operator P + . A common feature of Equations (23) and (24) is that their nonlinear parts contain only cross terms with respect to the functions ψ + and ψ (in other words, they are linear in ψ + and ψ ). This will allow us to apply to Equation (23) a number of techniques [26,27,30,34], which were developed for LGE, and, as a result, to construct spatially localized (Section 4) and spatially periodic (Section 5) solutions of the problem.

4. Spatially Localized Solutions

We are looking for a solution of Equation (23) as a sum of N logarithmic singularities,
  ψ + ( x , t ) = j = 1 N s j ln ( x Z j ( t ) ) ,
where s j are complex constants, Z j ( t ) are complex functions specifying positions of singularities in the complex plane x . Let Z j ( t ) = X j ( t ) + i Y j ( t ) , where X j ( t ) = R e Z j ( t ) and Y j ( t ) = Im Z j ( t ) . Since   ψ + must be analytic in the upper half-plane of the variable x , all singularities are located in the lower half-plane, i.e., the conditions Y j ( t ) < 0 are valid for all j = 1 ,   2 , ,   N . If it becomes Y j ( t ) = 0 at some moment, then the solution will be destroyed. Note that the representation (25) corresponds to a system of moving virtual point vortex-sinks/sources in the region y < 0 , i.e., outside the flow domain. Each singularity in (25) corresponds to a separate spatially localized perturbation of the velocity field.
By substituting (25) into (23) and then expanding the resulting expression into simple fractions (pole decomposition), we get the following system of N ODEs for the motion of singular points:
d d t Z j = m = 1 N 2 s m ¯ Z j Z m ¯ i V , j = 1 , 2 , , N .
Thus, the analysis of the integro-differential Equation (23) is reduced to the analysis of the ODE system (26), for which a significant number of exact solutions can be found (note that a similar system of ODEs arising from the LGE can be easily integrated).
First, we consider the simplest case of N = 1 , i.e., we look for a solution in the form
  ψ + ( x , t ) = s ln ( x Z ( t ) ) .
According to (26), the singularity position is described by ODE as
d d t Z = 2 s ¯ Z Z ¯ i V .
Let s = a + i b and Z ( t ) = X ( t ) + i Y ( t ) . Separating (28) into real and imaginary parts, we arrive at a pair of ODEs
d d t X = b Y ,
d d t Y = a Y V
with the initial conditions X ( 0 ) = X 0 and Y ( 0 ) = Y 0 < 0 .
Equation (30) for Y admits the stationary solution
Y = Y st = Y 0 = a / V ,
where a > 0 . The corresponding solution to Equation (29) has the form
X = ( b / a ) V t + X 0 ,
i.e., the singularity moves with constant velocity parallel to the real axis. Such a solution can be interpreted as a soliton solution of Equation (17); however, this exact solution does not make sense in the context of the problem under consideration, since it violates the Condition (18) of non-positivity of the normal velocity at the fluid boundary, v b ( x , t ) V + i ψ x + i ψ x 0 . Indeed, the normal velocity for the representation (27) is given by the expression
  v b ( x , t ) = V + i s x Z ( t ) i s ¯ x Z ¯ ( t ) = V 2 b ( x X ( t ) ) + a Y ( t ) ( x X ( t ) ) 2 + Y 2 ( t ) .
For the soliton solutions (31) and (32), the characteristic distributions of the velocity along the boundary are shown in Figure 1. It can be seen that there is always a region of positive velocities, i.e., formally, the fluid flows into the region y > 0 through the boundary y = 0 , which violates the requirements imposed on possible solutions.
It can be seen from Equation (30) that the singularity always moves away from the real axis for a 0 . For a > 0 , it moves away from the real axis if Y 0 < Y st and approaches it if Y 0 > Y st . In the latter case, the solution is destroyed as a result of violation of the condition of analyticity of the function ψ + when the singularity reaches the real axis at some finite time t c with asymptotic behavior Y 2 a ( t c t ) . However, such a scenario of the solution destruction has no physical meaning, since when the condition Y 0 > Y st holds, the non-positivity condition for the normal velocity at the fluid boundary is obviously violated (see Figure 1 for Y 0 = Y st ).
We consider several simple cases that do not violate the condition   v b 0 . Let a > 0 and b = 0 . According to (33), the distribution of the normal velocity at the boundary is determined by the expression
  v b ( x , t ) = V 2 a Y ( t ) x 2 + Y 2 ( t ) .
In this case, as can be seen from ODE (29), the singularity does not move along the real axis, i.e., X ( t ) = X 0 = const . We choose the following initial conditions: X 0 = 0 and Y 0 = 2 a / V . This gives the initial velocity distribution (22), where the velocity vanishes at a single point at the origin:   v b ( 0 , 0 ) = 0 . As stated in Section 3, such an initial condition corresponds to the situation where there was an obstacle at the boundary y = 0 , which was removed at time   t = 0 .
In the linear approximation (i.e., we consider the linearized Equation (20) instead of the original one (17)), ODEs (29) and (30) take the trivial form
d d t X = 0 , d d t Y = V ,
i.e., the singularity moves away from the real axis with constant velocity V :
X lin ( t ) = X 0 = 0 , Y lin ( t ) = Y 0 V t =   2 a / V V t .
According to (34), the normal velocity at the center of the perturbation evolves as
  v b ( 0 , t ) = V 2 a Y ( t ) .
Then, as a result of the linear relaxation of the perturbation, the velocity reaches the value of   0.5 V (i.e., half of the unperturbed flow velocity) at the time instant t lin = 2 a / V 2 , when Y ( t lin ) =   4 a / V .
Let us consider how nonlinearity will affect the relaxation of the perturbation. The exact solution of the nonlinear Equation (30) is written in implicit form as
Y ( t ) + 2 a / V V a V 2 l n ( a V Y ( t ) a ) = t .
The velocity reaches the value of   0.5 V at the time instant
t nl = 2 a V 2 ( 1 + l n 3 ) .
The ratio of the nonlinear and linear relaxation times is
t nl / t lin = 1 + l n 3 1.55 ,
i.e., the relaxation proceeds more than one and a half times slower than in the linear approximation. This is demonstrated in Figure 2, which shows the distribution of the normal velocity at the boundary ( v b ) at the initial moment   t = 0 and at the moment   t = t nl both for the exact solution and for the linear approximation. It can be seen that the influence of nonlinear terms of Equation (17) significantly affects the process of relaxation of the perturbation of the basic fluid flow.
Now let a = 0 and b > 0 . Then, according to (33), the distribution of the normal velocity at the boundary is determined by the expression
  v b ( x , t ) = V 2 b ( x X ( t ) ) ( x X ( t ) ) 2 + Y 2 ( t ) .
We take the initial conditions X 0 = 0 and Y 0 = b / V , which correspond to the initial velocity distribution
  v b ( x , 0 ) = V 2 b V 2 x V 2 x 2 + b 2 .
According to this expression, at the initial time   t = 0 , the velocity vanishes at the only point   x = b / V , i.e.,   v b ( b / V , 0 ) = 0 (the red line in Figure 3). Note that this perturbation of the velocity field, in contrast to (22), is not symmetric. This can be associated with a different geometry of the obstacle on the flow path.
The solution of the nonlinear ODEs (29) and (30) describing the singularity motion is
X ( t ) = b V ln ( b + V 2 t b ) , Y ( t ) = b / V V t .
Trivial equations
d d t X = 0 , d d t Y = V
are the linear analog of ODEs (29) and (30). Their solution has the form
X lin ( t ) = 0 , Y lin ( t ) = Y ( t ) = b / V V t .
Figure 3 shows the evolution of the normal velocity at the boundary ( v b ) corresponding to the exact solution (36) as well as to the linear approximation (37). The relaxation rate of the flow velocity perturbation for this case is the same in linear and nonlinear approximations. The singular point moves away from the real axis at the same velocity. However, linear and nonlinear velocity distributions are noticeably different, which is associated with the displacement of the perturbation along the x axis according to the logarithmic dependence in (36).
Let us now discuss the case N 2 . In the general case, the corresponding ODEs are non-integrable, and the dynamics of singular points are extremely complex; however, two relatively simple cases, viz. the cases of (i) independent and (ii) collective motion of singularities can be singled out.
The case of independent motion is realized if the interaction between singularities is weak. This, as can be easily seen from ODEs (26), occurs if the conditions
| Z j Z m ¯ | > > | Z j Z j ¯ | , j m
are satisfied, i.e., singular points are located relative to each other at greater distances than the distances from them to the real axis. In this case, the velocity field perturbations corresponding to the singularities are separated in space and evolve independently of each other according to the solutions described above for N = 1 .
The case of collective motion of singularities is realized when distances between them on the complex plane are much smaller than the distances to the real axis
| Z j Z m | < < | Z j Z j ¯ | .
where j = 1 , 2 , , N and m = 1 , 2 , , N . Note that these inequalities are not opposite to the inequalities (38).
It is convenient to analyze ODEs (26) using the “center-of-mass” system of coordinates, the origin of which is given by the expression
Z ( t ) = 1 S j = 1 N s j Z j ( t ) , S = j = 1 N s j .
We put
Z j ( t ) = Z ( t ) + δ Z j ( t ) , j = 1 , 2 , , N ,
where δ Z j ( t ) are small deviations of the positions of singularities from the position of the “center of mass”. In view of the Definition (40), the following relation
j = 1 N s j δ Z j ( t ) = 0
holds. By substituting (41) into ODEs (26) and then expanding these equations in power series of small quantities ( δ Z j δ Z m ¯ ) / ( Z Z ¯ ) , we get, in the main order, the equation of motion of the “center of mass”:
d d t Z = 2 S ¯ Z Z ¯ i V .
After the replacement S s , it will coincide with the equation of motion of a single singularity (28). Its exact implicit solution is
Y ( t ) Y 0 V a V 2 l n ( a + V Y ( t ) a + V Y 0 ) = t ,
X ( t ) = X 0 0 t b Y ( t ) d t ,
where, as before, Z = X + i Y , S = a + i b , X 0 = X ( 0 ) , and Y 0 = Y ( 0 ) .
In the next order of smallness, we obtain linear equations for the deviations δ Z j ( t ) ,
d d t δ Z j = 1 ( Z Z ¯ ) 2 m = 1 N 2 s m ¯ ( δ Z j δ Z m ¯ ) , j = 1 , 2 , , N ,
or, taking into account the Formulas (40) and (42), the following compact equations:
d d t δ Z j = 2 S ¯ ( Z Z ¯ ) 2 δ Z j = S ¯ 2 ( Im Z ) 2 δ Z j , j = 1 , 2 , , N .
Note that the Condition (42) agrees with these equations:
d d t j = 1 N s j δ Z j = 2 S ¯ ( Z Z ¯ ) 2 j = 1 N s j δ Z j = 0 .
Thus, in the considered approximation, the initial system of Equation (26) is split into independent equations for the quantities δ Z j ( t ) , which drastically simplifies its analysis. We introduce the auxiliary time variable
τ = d t 2 ( I m Z ) 2 .
Then Equation (45) take the simple form
d d τ δ Z j = S ¯ δ Z j , j = 1 , 2 , , N .
Their obvious solutions are
δ Z j ( τ ) = δ Z j ( 0 ) e S ¯ τ , j = 1 , 2 , , N .
Hence it is clear that the character of the collective behavior of singularities is determined by the quantity S ¯ . If R e S ¯ < 0 , then singular points converge with time. If R e S ¯ > 0 , then singular points move away from each other. For R e S ¯ = 0 , they rotate around the “center of mass”, which moves according to Solutions (43) and (44).
It can be expected that, in the intermediate case where inequalities (38) or (39) are not satisfied, the dynamics of singularities will be extremely complex and possibly chaotic.

5. Spatially Periodic Solutions

In this section, we will demonstrate the possibility of constructing exact spatially periodic solutions of the problem under consideration. We will look for solutions with the period 2 π / k ( k > 0 is the wave number of the basic harmonic) in the form of a sum of N + 1 harmonics, including the zero one, with (in the general case) complex amplitudes A j ,
  ψ + ( x , t ) = 1 2 A 0 ( t ) + j = 1 N A j ( t ) e i j k x .
A remarkable feature of the nonlinear part of Equation (17) (or of Equation (19)) is that it does not generate new harmonics with numbers greater than N . By substituting (46) into Equation (23), we arrive at the system of N + 1 equations for the temporal evolution of the amplitudes A j :
d d t A j = 2 k 2 m = 0 N j ( j + m ) m A j + m A m ¯ j k V A j , j = 0 ,   1 ,   2 , ,   N .
It is clear that this reduction opens up the possibility of constructing a significant number of exact analytical solutions to the problem.
We restrict ourselves here to a relatively simple special case of N = 2 with real amplitudes A 0 , 1 , 2 . For simplicity, let us take k = 1 (in which case the wavelength equals 2 π ). This corresponds to the following representation for the function ψ + :
  ψ + ( x , t ) = 1 2 A 0 ( t ) + A 1 ( t ) e i x + A 2 ( t ) e 2 i x .
According to (47), such a substitution leads to the system of three nonlinear ODEs,
  d d t A 0 = 2 A 1 2 8 A 2 2 ,
d d t A 1 = 4 A 1 A 2 V A 1 ,
d d t A 2 = 2 V A 2 .
It is clear that, in the linear approximation, these equations will describe the exponential relaxation of perturbations in accordance with the dispersion law (21). Let us consider how nonlinearities affect the dynamics of perturbations. It is immediately evident from Equation (50) that, for   A 2 < 0 , the nonlinearity slows down the relaxation of the basic harmonic   A 1 . The exact solution of ODEs (50) and (51) is
  A 2 ( t ) = A 2 ( 0 ) e 2 V t ,
  A 1 ( t ) = A 1 ( 0 ) exp [ 2 A 2 ( 0 ) V ( e 2 V t 1 ) V t ] ,
or
  A 2 ( t ) = A 2 ( 0 ) e 2 V t ,   A 1 ( t ) = A 1 ( 0 ) e V t
in the linear approximation. We do not present the solution of Equation (49) for the zero harmonic   A 0 ( t ) here since it does not affect the velocity distribution.
The distribution of the normal velocity at the boundary for the representation (48) is given by the expression
  v b ( x , t ) = V 2 A 1 ( t ) cos x 4 A 2 ( t ) cos 2 x .
We take   A 1 ( 0 ) = 0.4 V and A 2 ( 0 ) = 0.05 V as initial conditions. This corresponds to a periodic system of obstacles on the flow path at the boundary y = 0 removed at time t = 0 . We have v b ( 2 π n , 0 ) = 0 , where n is an integer, and v b ( x , 0 ) < 0 for x 2 π n , i.e., the necessary condition v b ( x , 0 ) 0 is satisfied.
The velocity field relaxation dynamics for the exact solutions (52) and (53) is shown in Figure 4. We do not show the solution (54) obtained in the linear approximation here since it differs only slightly from the exact one. So, for instance, the velocity   v b at the point   x = 0 reaches the value of   0.5 V (i.e., half of the unperturbed flow velocity) at the time t nl 0.66 , while in the linear approximation, this happens only slightly faster, at the time t lin 0.60 . The ratio of these times is t nl / t lin 1.1 , which is noticeably smaller than that of the spatially localized solution considered in Section 4, where the analogous ratio exceeded the value of 1.5.

6. Conclusions

In the present work, the 2D unsteady flow described by the partial differential Equations (7)–(9) with the additional kinematic Condition (13) has been considered. The problem was reduced to the single integro-differential Equation (19), which can be considered as a nonlocal generalization of the Hopf Equation (1). It has been demonstrated that the new equation (to be precise, Equation (17) related to Equation (19) by the change u = ψ x ) can be reduced to a system of a finite number of ODEs that describe a rather complex evolution of spatially localized or spatially periodic perturbations of the velocity field.
For the problem under consideration, the unperturbed uniform flow (12) is stable in the linear approximation, i.e., when the amplitude of the velocity field perturbations is small compared to the velocity of the unperturbed flow ( V in absolute value). The main attention in the work was paid to studying the influence of nonlinearity on the relaxation of perturbations of the velocity field. The revealed possibility of reducing the original equations to ODEs provided a unique opportunity to analytically compare the evolution of perturbations described by exact solutions and by solutions obtained in the linear approximation. We have demonstrated that nonlinearity, as expected, significantly affects the dynamics of the velocity field in the case when the amplitude of perturbations is comparable with the velocity of the unperturbed flow (see Figure 2 and Figure 3). For all considered examples that satisfy the Condition (13), accounting for nonlinearity did not lead to the loss of stability of the basic flow. Nonlinear effects only led to a significant decrease in the rate of relaxation of perturbations. This suggests that the basic flow is stable not only with respect to small perturbations but also with respect to perturbations of finite amplitude.
The analysis of Equation (17) derived in this work (or of the equivalent Equation (19) which generalizes the Hopf equation) was carried out only in the case of a finite value of the parameter V . This was due to the specifics of the original problem (7)–(9) and its interpretation as describing the upstream flow in front of the linear sink. However, in our opinion, Equation (19) for the case of V = 0 , when it takes the form
u t + u u x + H [ u ]   H [ u x   ] = 0 ,
may also be of independent interest. As a rule, in the absence of dispersion, the influence of nonlinearity leads to wave breaking. Examples of this are the classical Hopf Equation (1), its complex version (2), as well as the related nonlocal Equation (4), which contains the same terms as (56), but with different signs. As can be easily understood from the analysis of Section 4, Equation (56) with two competing quadratic nonlinearities, u u x and H [ u ]   H [ u x   ] , characteristics for the problems of fluid dynamics, exhibits fundamentally different properties. It admits solutions corresponding to an arbitrary number of structurally stable solitary waves. In particular, it admits a solution of the type of the “algebraic” Benjamin–Ono soliton [35], which moves without breaking with the constant velocity c ,
u ( x , t ) = 2 c Y 0 2 ( x X 0 c t ) 2 + Y 0 2 ,
where X 0 , Y 0 are some constants. We note that a preliminary analysis of the interaction of several solitons of the form (57) indicates its rather complex character. The form of Poincaré sections in the numerical solution of ODEs (26) for V = 0 in a number of examples testifies to the chaotic dynamics of N 3 interacting solitary waves. However, a detailed consideration of the properties of Equation (56) is far beyond the scope of this work and will be the subject of a separate study.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Hopf, E. The partial differential equation ut + uux = μxx. Commun. Pure Appl. Math. 1950, 3, 201–230. [Google Scholar] [CrossRef]
  2. Whitham, G.B. Linear and Nonlinear Waves; John Wiley & Sons: New York, NY, USA, 1999. [Google Scholar] [CrossRef]
  3. Karabut, E.A.; Zhuravleva, E.N. Unsteady flows with a zero acceleration on the free boundary. J. Fluid Mech. 2014, 754, 308–331. [Google Scholar] [CrossRef]
  4. Karabut, E.A.; Zhuravleva, E.N. Unsteady flows with a zero acceleration on the free boundary. Dokl. Phys. 2014, 59, 480–483. [Google Scholar] [CrossRef]
  5. Zubarev, N.M.; Karabut, E.A. Exact local solutions for the formation of singularities on the free surface of an ideal fluid. JETP Lett. 2018, 107, 412–417. [Google Scholar] [CrossRef]
  6. Caflisch, R.E.; Ercolani, N.; Hou, T.Y.; Landis, Y. Multi-valued solutions and branch point singularities for nonlinear hyperbolic or elliptic systems. Commun. Pure Appl. Math. 1993, 46, 453–499. [Google Scholar] [CrossRef]
  7. Gao, Y.; Gao, Y.; Liu, J.G. Large time behavior, bi-Hamiltonian structure and kinetic formulation for a complex Burgers equation. Quart. Appl. Math. 2021, LXXIX, 55–102. [Google Scholar] [CrossRef]
  8. Konopelchenko, B.G.; Ortenzi, G. Homogeneous Euler equation: Blow-ups, gradient catastrophes and singularity of mappings. J. Phys. A Math. Theor. 2021, 55, 035203. [Google Scholar] [CrossRef]
  9. Karabut, E.A.; Zhuravleva, E.N.; Zubarev, N.M. Application of transport equations for constructing exact solutions for the problem of motion of a fluid with a free boundary. J. Fluid Mech. 2020, 890, A13. [Google Scholar] [CrossRef]
  10. Karabut, E.A.; Zhuravleva, E.N.; Zubarev, N.M.; Zubareva, O.V. Evolution of nonlinear perturbations for a fluid flow with a free boundary. Exact results. J. Fluid Mech. 2022, 953, A1. [Google Scholar] [CrossRef]
  11. Karabut, E.A.; Zhuravleva, E.N. Reproduction of solutions in the plane problem on motion of a free-boundary fluid. Dokl. Phys. 2016, 61, 347–350. [Google Scholar] [CrossRef]
  12. Zhuravleva, E.N. Some exact solutions that describe an ideal fluid flow with a free boundary. Phys. Fluids 2023, 35, 017104. [Google Scholar] [CrossRef]
  13. Kuznetsov, E.A.; Spector, M.D.; Zakharov, V.E. Surface singularities of ideal fluid. Phys. Lett. A 1993, 182, 387–393. [Google Scholar] [CrossRef]
  14. Kuznetsov, E.A.; Spector, M.D.; Zakharov, V.E. Formation of singularities on the free surface of an ideal fluid. Phys. Rev. E 1994, 49, 1283–1290. [Google Scholar] [CrossRef]
  15. Zubarev, N.M. Formation of root singularities on the free surface of a conducting fluid in an electric field. Phys. Lett. A 1998, 243, 128–131. [Google Scholar] [CrossRef]
  16. Zubarev, N.M.; Kuznetsov, E.A. Singularity formation on a fluid interface during the Kelvin-Helmholtz instability development. J. Exp. Theor. Phys. 2014, 119, 169–178. [Google Scholar] [CrossRef]
  17. Zhdanov, S.K.; Trubnikov, B.A. Quasigas approximation in treating electron-beam bunching in a plasma and the tangential discontinuity in hydrodynamics. Sov. Phys. JETP 1988, 67, 1575–1582. [Google Scholar]
  18. Moore, D.W. The spontaneous appearance of a singularity in the shape of an evolving vortex sheet. Proc. R. Soc. Lond. Ser. A 1979, 365, 105–119. [Google Scholar] [CrossRef]
  19. Polubarinova-Kochina, P.Y. On motion of the contour of an oil layer. Dokl. Akad. Nauk SSSR 1945, 47, 254–257. [Google Scholar]
  20. Galin, L.A. Unsteady filtration with a free surface. Dokl. Akad. Nauk SSSR 1945, 47, 246–249. [Google Scholar]
  21. Clarke, N.S. On two-dimensional inviscid flow in a waterfall. J. Fluid Mech. 1965, 22, 359–369. [Google Scholar] [CrossRef]
  22. Chow, W.L.; Han, T. Inviscid solution for the problem of free overfall. J. Appl. Mech. 1979, 46, 1–5. [Google Scholar] [CrossRef]
  23. Goh, K.H.M.; Tuck, E.O. Thick waterfalls from horizontal slots. J. Eng. Math. 1985, 19, 341–349. [Google Scholar] [CrossRef]
  24. Akbari, S.; Faghiri, S.; Poureslami, P.; Hosseinzadeh, K.; Shafii, M.B. Analytical solution of non-Fourier heat conduction in a 3-D hollow sphere under time-space varying boundary conditions. Heliyon 2022, 8, e12496. [Google Scholar] [CrossRef] [PubMed]
  25. Saffman, P.G.; Taylor, G.I. The penetration of a fluid into a porous medium or Hele-Shaw cell containing a more viscous liquid. Proc. R. Soc. Lond. Ser. A 1958, 245, 312–329. [Google Scholar] [CrossRef]
  26. Howison, S.D. Cusp development in Hele–Shaw flow with a free surface. SIAM J. Appl. Math. 1986, 46, 20–26. [Google Scholar] [CrossRef]
  27. Bensimon, D.; Kadanoff, L.P.; Liang, S.; Shraiman, B.I.; Tang, C. Viscous flows in two dimensions. Rev. Mod. Phys. 1986, 58, 977–999. [Google Scholar] [CrossRef]
  28. Dyachenko, A.I.; Zakharov, V.E.; Kuznetsov, E.A. Nonlinear dynamics of the free surface of an ideal fluid. Plasma Phys. Rep. 1996, 22, 829–840. [Google Scholar]
  29. Zakharov, V.E.; Dyachenko, A.I. High-Jacobian approximation in the free surface dynamics of an ideal fluid. Phys. D Nonlinear Phenom. 1996, 98, 652–654. [Google Scholar] [CrossRef]
  30. Shraiman, B.I.; Bensimon, D. Singularities in nonlocal interface dynamics. Phys. Rev. A 1984, 30, 2840–2842. [Google Scholar] [CrossRef]
  31. Zubarev, N.M. Charged-surface instability development in liquid helium: An exact solution. JETP Lett. 2000, 71, 367–369. [Google Scholar] [CrossRef]
  32. Lushnikov, P.M.; Zubarev, N.M. Exact solutions for nonlinear development of a Kelvin-Helmholtz instability for the counterflow of superfluid and normal components of Helium II. Phys. Rev. Lett. 2018, 120, 204504. [Google Scholar] [CrossRef] [PubMed]
  33. Mineev-Weinstein, M.; Wiegmann, P.B.; Zabrodin, A. Integrable structure of interface dynamics. Phys. Rev. Lett. 2000, 84, 5106–5109. [Google Scholar] [CrossRef] [PubMed]
  34. Mineev-Weinstein, M.B.; Dawson, S.P. Class of nonsingular exact solutions for Laplacian pattern formation. Phys. Rev. E 1994, 50, R24. [Google Scholar] [CrossRef] [PubMed]
  35. Ono, H. Algebraic solitary waves in stratified fluids. J. Phys. Soc. Jpn. 1975, 39, 1082–1091. [Google Scholar] [CrossRef]
Figure 1. Normal velocity at the boundary v b for soliton solutions (31)–(33) of the key Equation (17). Three characteristic solutions are shown for t = 0 , X 0 = 0 , Y 0 = 1 , V = 1 , a = 1 , and b = 0 ,   0.5 ,   1 .
Figure 1. Normal velocity at the boundary v b for soliton solutions (31)–(33) of the key Equation (17). Three characteristic solutions are shown for t = 0 , X 0 = 0 , Y 0 = 1 , V = 1 , a = 1 , and b = 0 ,   0.5 ,   1 .
Mathematics 11 01999 g001
Figure 2. Relaxation of the perturbation (22) of the normal velocity at the flow boundary for V = 1 , a = 1 , b = 0 , X 0 = 0 , Y 0 = 2 . The red solid line corresponds to the initial time instant t = 0 ; the blue solid line corresponds to the time instant t = t nl 3.1 when the velocity at the center of the perturbation reaches the value 0.5 ; the green dashed line corresponds to the evolution of the velocity field in the linear approximation (the same time instant t 3.1 ).
Figure 2. Relaxation of the perturbation (22) of the normal velocity at the flow boundary for V = 1 , a = 1 , b = 0 , X 0 = 0 , Y 0 = 2 . The red solid line corresponds to the initial time instant t = 0 ; the blue solid line corresponds to the time instant t = t nl 3.1 when the velocity at the center of the perturbation reaches the value 0.5 ; the green dashed line corresponds to the evolution of the velocity field in the linear approximation (the same time instant t 3.1 ).
Mathematics 11 01999 g002
Figure 3. Relaxation of the perturbation (35) of the normal velocity at the boundary for V = 1 , a = 0 , b = 1 , X 0 = 0 , Y 0 = 1 . The red solid line corresponds to the initial time instant t = 0 ; the blue solid line corresponds to the time instant t = 2 ; the green solid line corresponds to the time instant t = 3 . The blue and green dashed lines correspond to the evolution of the velocity field in the linear approximation (time instants t = 2 ,   3 , respectively).
Figure 3. Relaxation of the perturbation (35) of the normal velocity at the boundary for V = 1 , a = 0 , b = 1 , X 0 = 0 , Y 0 = 1 . The red solid line corresponds to the initial time instant t = 0 ; the blue solid line corresponds to the time instant t = 2 ; the green solid line corresponds to the time instant t = 3 . The blue and green dashed lines correspond to the evolution of the velocity field in the linear approximation (time instants t = 2 ,   3 , respectively).
Mathematics 11 01999 g003
Figure 4. Normal velocity at the flow boundary for the periodic solution (52), (53), and (55) of the problem. Velocity distributions are shown at successive times t = 0 ,   0.5 ,   1 ,   1.5 for V = 1 ,   A 1 ( 0 ) = 0.4 , and A 2 ( 0 ) = 0.05 .
Figure 4. Normal velocity at the flow boundary for the periodic solution (52), (53), and (55) of the problem. Velocity distributions are shown at successive times t = 0 ,   0.5 ,   1 ,   1.5 for V = 1 ,   A 1 ( 0 ) = 0.4 , and A 2 ( 0 ) = 0.05 .
Mathematics 11 01999 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zubarev, N.M. Nonlinear Dynamics of Perturbations of a Plane Potential Fluid Flow: Nonlocal Generalization of the Hopf Equation. Mathematics 2023, 11, 1999. https://doi.org/10.3390/math11091999

AMA Style

Zubarev NM. Nonlinear Dynamics of Perturbations of a Plane Potential Fluid Flow: Nonlocal Generalization of the Hopf Equation. Mathematics. 2023; 11(9):1999. https://doi.org/10.3390/math11091999

Chicago/Turabian Style

Zubarev, Nikolay M. 2023. "Nonlinear Dynamics of Perturbations of a Plane Potential Fluid Flow: Nonlocal Generalization of the Hopf Equation" Mathematics 11, no. 9: 1999. https://doi.org/10.3390/math11091999

APA Style

Zubarev, N. M. (2023). Nonlinear Dynamics of Perturbations of a Plane Potential Fluid Flow: Nonlocal Generalization of the Hopf Equation. Mathematics, 11(9), 1999. https://doi.org/10.3390/math11091999

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

Article Metrics

Back to TopTop