Symmetry Analysis of Ideal Fluid Equations in Terms of Trajectories and Weber’s Potential Текст научной статьи по специальности «Математика»

УДК 532.51

Symmetry Analysis of Ideal Fluid Equations in Terms of Trajectories and Weber's Potential

Victor K. Andreev*

Institute of Computational Modelling SB RAS Akademgorodok, 50/44, Krasnoyarsk, 660036 Institute of Mathematics and Fundamental Informatics

Siberian Federal University Svobodny, 79, Krasnoyarsk, 660041

Russia

Daria A. Krasnova^

Institute of Mathematics and Fundamental Informatics

Siberian Federal University Svobodny, 79, Krasnoyarsk, 660041

Russia

Received 24.11.2018, received in revised form 26.12.2018, accepted 20.01.2019 The 2D perfect fluid motions equations in Lagrangian coordinates are considered. If body forces are potential one, then there is the general integral called Weber's integral and the resulting system includes initial data which in fact make the problem of group-theoretical classification actual. It is established that the basic group becomes infinite-dimensional with respect to the space variable too. The exceptional values of arbitrary initial vorticity are obtained at which we can be observed further extension of the group. Group properties of Euler equations in arbitrary Lagrangian coordinates are also considered and some exact solutions are constructed.

Keywords: Euler equations, symmetry analysis, Weber's transformation, equivalence transformation, group classification.

DOI: 10.17516/1997-1397-2019-12-2-133-144.

Introduction

Group analysis of differential equations (Lie-Ovsyannikov symmetry method [1-3]) became a powerful tool for studying nonlinear equations and boundary value problems. This analysis is especially fruitful in application to the basic equations of mechanics and physics because the invariance principles are already involved in their derivation. Symmetry methods are successfully used to study the mathematical models in hydrodynamics and to construct exact solutions of non-linear problems, e.g. in [4-6]. An admissible group characterizes symmetry properties of differential equations and used for complete integration or construction of certain classes of exact solutions as well as qualitative investigation of the equations.

In describing the motion of an ideal incompressible fluid that has a free boundary the problem, in particular waves problem, is reduced to finding the solution of the Euler equations with

* [email protected]

1 [email protected] © Siberian Federal University. All rights reserved

the fulfilment of the kinematic and the dynamic condition at the free surface. The kinematic condition allows us to convert this problem to another one in which the domain is fixed. It is achieved by the transition to the Lagrangian coordinates which are introduced as the coordinate values of the fluid particles at the initial time [7]. The new system of differential equations arised very is often used for the investigation solvability of this problem [8-11].

The paper purpose is to study group properties of equations which describe 2D perfect liquid flows. It is well known [12] that, there are two integrals of such equations in Lagrangian coordinates. First of them is the Cauchy one which implies that the particle vorticity is preserved under planar motions. The second integral is appeared due to Weber transformation which applies to Euler equations [12]. The differential equations system corresponding Cauchy integral is the closed system for only trajectories. In this case the problem of group-theoretic classification was solved in [13,14]. Another differential equations system is the closed system for trajectories and Weber's potential. The passage from one system to another mentioned above is a nonlocal transformation, therefore there should be no isomorphism between the basic Lie group of the equations under consideration. So, we consider the last differential equations system in terms of trajectories and Weber's potential. It is established that the basic group is infinite-dimensional with respect to the space variables. The exceptional values of initial vorticity at which we observe further extension of the group are obtained. Group properties in arbitrary Lagrangian coordinates are considered and two exact solution are constructed.

1. Governing equations

A nonstationary 3D flow of an ideal incompressible fluid is described by the following system of equations (Euler equations)

ut + uVu + 1 Vp = g(x,t), div u = 0, (1.1)

P

where p = const is the fluid density and g is the vector of body forces, u(x, t) and p(x, t) are the velocity vector and the pressure, respectively.

Let us introduce the Lagrangian coordinates £ by solving the Cauchy problem

§ = u(x^ x|t=0 = (O)

Without loss of generality we can assume that the fluid density equals unity. Then the system (1.1) can be rewritten in the following form [7]

M*(xtt - g(x,t))+ Vp = 0, div M-1xt =0. (1.3)

Here M is the Jacobi matric

z

"c

f*

/ x x n x

M = | * Vv y

V z n z

d(x) d(0

and M* is inverse matrix to M. The second equations of (1.3) is equivalent to equation det M = 1 [7].

If body forces have a potential g = Vxh then from the first eq'n of (1.3) we derive [7]

x

T*

M *xt = V^> + uo(£), (1.4)

ft + P =2 |xt|2 + h + X(t), (1.5)

where u0(£) is the initial value of u(x, t) and x(t) is arbitrary function. The relation (1.5) is the Weber's integral. If flow is potential one then this integral coincides with Cauchy-Lagrange integral [6,7].

Function f(£, t) is the solution of second order elliptic eq'n

div [M-iM*-i(Vf + uo(£))] =0 (1.6)

and can be called Weber's potential.

For 2D motions let us put u0 = (u(£,n),v(£,n)), £ = (€,n), x = (x(£,n,t),y(£,n,t)). Then the equations (1.4), (1.6) imply the closed differential system

xt = yv(fi + u(n)) - yi(fn + v(& n)); (1.7)

yt = -xn(fi + u(£, n)) + Xi(fn + v(£, n)); (1.8)

Xiyv - Xvyi = 1; (1.9)

ui + vn = 0, (1.10)

where x(£, n, t), y(£, n, t) are the trajectories of liquid particles and f(£, n, t) is the Weber potential. The relationship (1.10) is the compatibility condition for velocity vector at initial time. The pressure distribution can be found from Weber's integral (1.5) in the following form

P(tn,t) = 2(x2 + yl) - ft + h(tn,t) (1.11)

and h = -gy(£, n, t) for water waves, where g = const is the gravitational acceleration. Remark 1. Equations (1.7), (1.8) are equivalent to the following

fi = xtxi + ytyi - u, fn = xtxn + yty - v-The compatibility condition of these equations implies the eq'n

xixnt - xnxit + yiynt - ynyit = n) = vi - Un = 0. (1.12)

It is the conservation vorticity law (Cauchy integral [7]). The equations (1.9) and (1.12) are the closed system for functions x(£, n, t), y(£, n, t). The group properties for this system were studied in [13,14]. Due to inequality (1.12) an unsteady flow is rotational one for all time.

2. Symmetry properties of system (1.7)—(1.9)

Firstly, let us compute the equivalence transformation. We can demonstrate that the equivalence transformation is

t = ait, x = ai(x cos 0,2 + y sin ai) + 0,3, y = ai(-x sin 0,2 + y cos ai) + 0,4,

t = ai(£cos a2 + nsin a2) + a3, n = ai(-Csin a2 + ncos a2) + a4, (2.1)

f = ai f + + b2n + d(t);

u(£,fj ) = a5u(Z,n )=0. (2.2)

Here a1 = 0, a2, a3, a4, a5 = 0, b1, b2 are arbitrary constants and d(t) is arbitrary function. Of course compatibility condition (1.10) is invariant with respect to transformation (2.1). At last, we should supplement (2.1) with the discrete transformation

x = -x, y = -y. (2.3)

Generator for system (1.7)-(1.9) is sought in the form

Y=e1 -+ed+ed+v1—+n2d+rf~, (2.4)

dt de dn dx dy By1

where £\ rf depend on all variables t, £, n, x, y, y. To derive the determining equations it is necessary to apply the first prolongation Y(1) of generator Y to system (1.7)-(1.9). After a large member of calculations the determining system can be reduced to the following one

e1 = e1(t), e2 = e2(e,n), e3 = e3(e,n),

n1 = C1x + C2y + n(t), n2 = C\y - C2x + m(t), e| + e3 = 2C,

n3 = (2C1 - e1)y + ntx + mty + h(e, n) + d(t), (2.5)

«(en - e!) - e2«« - e3«v - efv - h« = 0, v(ei - e1) - e2v« - e\ - el« - hv = 0,

here n(t), m(t), d(t) are arbitrary smooth functions, C1, C2 are constants. From the last two equations (2.5) for h we get that e1 = C3t + C4 and

e2w« + e3^^ + C3 w = 0, (2.6)

It is the classifying equation and now we can find the function h in the form

h = J [u(en - C3) - e2«« - e3«n - e|v] de + Kef - C3) - e2v« - e% - el«] dn. (2.7)

Basis for the Lie algebra is

Lo : Y1 = dt, Y2 = ydx - xdy, Yn = n(t)3x + nt(t)xdv,

Ym = m(t)dy + mt(t)ydv, Yd = d(t)3v.

It is convenient to use the new unknowns £2(e,n), £3(e, n) by the formulae

e2 = C1e+e2(e,n), e3 = cm+e3(e,n), (2.9)

then the classifying part of system (2.5) of determining equations can be rewritten as

ee + en = 0, h« = u(C1 - C3 + e3) - (C1e+e2)«« - (cn+e3« - e3v,

hv = v(C1 - C3+e|) - (C1e+e2)v« - (cm+ev - ev«, (2.10)

(C1e+e2)^« + (Cm+fn + = 0.

If w = w0 = const, then C3 = 0 and L0 is prolonged by the generators

Y^ = ^d« - ^«+ hdv, Y3 = ed« + ndn + xdx + ydy + (2y + h)8v (2.11)

(2.8)

with arbitrary function Function h(£,r) is determined from (2.7), where C3 = 0,

£2 = ,2 = фп, £3 = £3 = —ф% with arbitrary function Ф(£,г), u = u(£,r) — ш0г for generator Y^ and C3 =0, £2 = £, £3 = n, u = u(£, r) + шоГ for generator Y3. The functions u(£, r), v(£, r) are the solution of Cauchy-Riemann system: u% + vn = 0, un — v% = 0.

If ш(£, r) = const basis for the Lie algebra extended L0 is given by the generators

Yf = д% — ^F)^ дп, Y4 = Stdt + £д% + rdn + xdx + ydy + (2 — S)^8V. (2.12)

Function F is arbitrary for YF and h(£, r) is found from (2.10), where C\ = C3 = 0, £2 = dF(ш)/dr, £3 = —dF(ш)/д£. For generator Y4 S is a constant and the vorticity ш(£,г) has the form

ш(£,г) = r- f (/), (2.13)

with arbitrary function f (Z), Z = £/r. In this case h = 0.

3. Invariance of the initial conditions

Solving the system of equations (1.7)-(1.9), which is not normal with respect to time, it is necessary to take into account the initial conditions

x = Ç, y = v, t = 0. (3.1)

We require the relations (3.1) are invariant under action of the generator Y from (2.5). The system of determining equations has the simpler form

Çi = C3t, e = C\Ç + C2V + n(0), Ç3 = Cl n - C2Ç + m(0),

n1 = Cix + C2y + n(t), n2 = Ciy - C2X + m(t), n3 = (2Ci - C3)p + ntx + mty + h(Ç, n) + d(t), (3.2)

hç = (Ci - C3)u - Ç2uç - Ç3uv + C2V, hn = (Ci - C3)v - Ç2ve - Ç3vv - C2U; (CiÇ + C2n + n(0))uç + (Cin - C2Ç + m(0))u>n + C3W = 0. (3.3)

The basic algebra L00 includes the generators

Zn = n(t)dx + nt(t)xdv, Zm = m(t)dy + mt(t)ydv, Zd = d(t)dv, n(0) = m(0) = 0.

(3.4)

Here the equivalence transformations have the form (2.1), (2.2). To carry out the group classification let us rewrite the classifying eq'n (3.3) in the form

(A£ + Bn + C)^ + (An — + D)wn + Hw = 0 (3.5)

with some constants A, B, C, D and H. It is easy to verity that eq'n (3.5) in invariant with respect to equivalence transformations (2.1) with the transformations of A, B, C, D, H such as

A = A, B = B, C = cos a2 C — sin a2D — a3A — a4B,

- „ - (3.6)

D = cos a2D — sin a2C — a4A + a4B, H = H. - 137-

Under integration of equation (3.6) we consider the following cases:

I. AB = 0. Then using (3.6) we can set C = D = 0. Eq'n (3.5) has the solution

w = exp[72 arctg(£/n)]/ ((£2 + nY/2 exp[7i arctg(£/n)j) (3.7)

with constants y2 = 0, y2 and arbitrary function /.

II. A = 0, B = 0. Then from (3.6) C = 0 and the solution of eq'n (3.5) has the form

w = (n + Y3 )74 /MM (3.8)

\n + 73/

with constants y3, 74 and arbitrary function /.

III. A = 0, B = 0. Then C = 0 and we obtain the solution of eq'n (3.5) in the form

w = exp

' . f 75 - e V

76 arcsml + (Y5 - e)2 )

fin2 + (75 - e)2) (3.9)

with constants 75, 76 and arbitrary function /.

IV. A = 0, B = 0. Then the solution of eq'n (3.5) can be written as

= f (7se - 7on), (3.10)

w

if C = 0, or

w = e-Y7n f (7se - Y9V), (3.11)

if D = 0 with 77, 7s, 79 are constants and f is arbitrary function.

Substituting solutions (3.7)—(3.11) to classifying equation (3.3) and using formulae (3.2) and (2.2) we obtain the extension of the Lie symmetry algebras spanned by generators Loo. The result of group classification is presented in Tab. 1. There we can find the form of the function w and basis of generator which are possessed by equations (1.7)—(1.9) and initial conditions (3.1).

Table 1. Results of group classification

w Generators Remarks

(3.7), f is arbitrary function L00 Z1

(3.7), f = 1 L00 Z1, Z2, Z3

(3.7), f = C78, 78 = 0 L00 Z4, Z5 Z = (e2 + n2)1/2 exp[Y1 arctg(e/n)]

(3.7), f = C72, 72 = 0, 71 = 1 L00 Z1, Z6 Z = (e2 + n2)1/2 exp[arctg(e/n)]

(3.8), f is arbitrary function L00

(3.8), f = 1 L00 Z7

(3.8), f = (C + 76)74, 74 = 0 L00 Zg z = e/(n+73 )

(3.9), f is arbitrary function L00

(3.9), f = 1 L00 Z9

(3.9), f = (C - 752)77, 77 = 0 L00 Z10, Z11 z = n2 + (75 - e)2

(3.10), f is arbitrary function L00

(3.10), f = 1 L00 Z12

(3.10), f = eYloZ, 710 = 0 L00 Z13 Z = 78e - 79n, + Y2 = 0

(3.10), f = c712 L00 Z14, Z15 Z = Y8e - Y9n,Y9 = 0

(3.10), f = (Z + 711)712, 711 =0 L00 Z16, Z17 Z = Y8e - Y9n,Y9 = 0

In Tab. 1 the following generators are involved in extension of the algebra Loo : Zi = Y2tdt + (e - nd + (e + n)dv + (x - y)dx + (x + y)dy + [(2 - 72)^ + h(£, n)]dv,

he = (1 - Y2)u - (e - n)ui - (e + n)uv - v, hn = (1 - Y2)v - (e - n)ve - (e + n)vn + u;

Z2 = ede + ndn + ydy + [2y + h(e, n)]dv,

he = u - eue - nun, hn = v - eve - nvn;

Z3 = -Y2tdt + nde - edn + ydx - xdy + [72^ + h(e, n)]dv,

he = 72U - nue + eun + v, hv = 72v - nve + evn - u;

Z4 = -Ystdt+ede+ndn + xdx+ydy + [(2 + 7s)y + h(e,n)]dv,

he = (1 + Y8)u - eue - nun, hn = (1 + 78)v - eve - nvn;

Z5 = -(72 + Y8)tdt + nde - edn + ydx - xdy + [(72 + Y8)y + h(e, n)]dv,

he = (72 + 78)u - nue + eun + v, hn = (72 + Y8)v - nve + evn - u;

Z6 = a(0)de - a(0)dn + a(t)dx - a(t)dy + [at(x + y) + h(e, n)]dv, he = -a(0)ue + a(0)un, hn = -a(0)ve + a(0)vn;

Z7 = -Y4tdt + ede + (n + 73)dn + xdx + (y + Y3a(t))dy + [(2 + 74)^ + 73aty + h(e, n)]dv, he = (1 + Y3)u - eue - (n + Y3)un, hn = (1 + Y3)v - eve - (n + Y3)vn;

Z8 = -Y6de - Y6a(t)dx + dn + a(t)dx + (-Y6atx + aty + h(e,n))dv, he = -Y6ue - un, hn = -Y6ve - vn, a(0) = 1;

Z9 = (e - Y5)de + ndn + (x - Y5a(t))dx + ydy + (2y - Y5atx + h(e, n))dv, he = u - (e - Y5)ue - nun, hn = v - (e - Y5)ve - nun;

Zio = -2Y7tdt + (e - Y5)de + ndn + (x - Y5a(t))dx + ydy + [2(y7 + 1)^ - Y5atx + h(e,n)]dv

he = (2Yr + 1)u - (e - Y5)ue - nun, hn = (2y7 + 1)v - (e - Y5)ve - nun; Zii = Y6tdt + nde - (e - Y5)dn + ydx - (x - Y5a(t))dy + (-Y6<f + Y5aty + h(e, n))dv,

he = Y6u - nue + eun + v, hn = Y6v - nve + evn - u, a(0) = 1;

Zi2 = Y7tdt + de + a(t)dx + (-Y7^ + atx)dv,

he = -Y7u - ue, hn = -Y7v - ve, a(0) = 1;

Zi3 = YioY9tdt + dn + a(t)dy + (-YioY9^ + aty + h(e, n))dv, he = -YioY9u - un, hn = -YioY9v - vn, a(0) = 1;

Zi4 = -Yi2tdt + ede + a(t)xdx + [(2 + Yi2y + atx + h(e, n)]dv,

he = (1 + Yi2)u - eue - nun, hn = (1 + Yi2)v - eve - nvn;

Zi5 = de - Y-iY8dn + a(t)dx - Y-iY8a(t)dy + (atx - y-1 Y8aty + h(e, n))dv, he = -ue - Y-iY8un, hn = -ve - Y-iY8vn, a(0) = 1;

Z\6 = -Yi2tdt + £d$ + (n- Y-1 m)dv + xdx + [y - Yc-1Y11a(t)]dy+

+ (<Pdv - Yg1Yiiaty + h(£, n))dv,

= (1 + Y12U - £u$ - (n - Yg1Yii)un, hn = (1 + Y12v - £ve - (n - Yg1Yn)vv, a(0) = 1'■>

Z17 = + a(t)dx + Y-1Y&a(t)dy + Yc-Ssdv + [yCYaty + h(£, n)]dv, = - Yg1Y8Uv, hv = - Y-1Y8vv, a(0) = 1.

Remark 2. For potential flows u = , v = p0n and range of p ^ p + p0 gives u = v = 0. Basic algebra Lie of generators is the following

H(t)3t - Ht(t)pOv, ydx - xdy, n(t)dx + nt(t)3v, m(t)dy + mt(t)3v,

d(t)6v, - dv, £3$ + ndv + xdx + ydy +

with arbitrary smooth functions n(t), n(t), m(t), d(t), ^(£,n).

4. Arbitrary Lagrangian coordinates

The Cauchy problem

dx dy

— = u(x,y,t), — = v(x,y,t),

x\t=0 = £, y\t=0 = n

defines the Lagrangian coordinates £ and n as the Cartesian coordinates of a fluid particle at the initial time. However, instead of £ and n distinguishing between points, we can take any quantities a and b connected with £ and n via a one-to-one correspondence, see Fig. 1,

£ = f (a, b), n = g(a,b), J = fagb - fbQa =0. (4.1)

Fig. 1. Arbitrary Lagrangian coordinates

System (1.7)-(1.9) is equivalent to the system

JXt = Yb($a + Ui) - Y0($6 + Vi); (4.2)

JYt = -Xb($a + Ui)+ X„($b + Vi); (4.3)

XaYb - XbYa = J(a,b) (4.4)

with the initial data

X|t=o = / (a,b), Y|t=o = g(a,b). (4.5) New "velocities" Ui(a,b), Vi(a,b) are connected with u(e,n), v(e,n) by the relations Ui(a, b) = /au(/(a, b), g(a, b)) + gav(/(a, b), g(a, b)),

(4.6)

Vi(a, b) = /bu(/(a, b), g(a, b)) + gbv(/(a, b), g(a, b)).

We substitute X(a,b,t)= x(e,n,t), Y(a,b,t)= y(e,n,t), $(a,b,t)= y(e,n,t) in to (4.2)-(4.4) with e and n defined by (4.1). It is easy to see that the function x, y, y satisfy system (1.7)-(1.9), moreover x = e, y = n under t = 0. Hence, the Lie group of transformations for the systems (1.7)-(1.9) and (4.2)-(4.4) are similar. The generator admitted by system (4.2)-(4.4) has the form

Yi = e i(t)dt + e2(a, b)da + e3(a, b)db + nidx + n2dy + n3d$,

where e 1 (t), ni, n2, n3 are the same as for generator Y. The coordinates e2(a, b), e3(a,b) satisfy the determining equations

(Je2)a + (je3)b = 2cj,

(4.7)

e2(w i /J )a + e3(w i/J)b + C3W i/J = 0, where wi(a,b) = J(a,b)w(/(a,b),g(a,b)) = Via - Ui,b is a new "vorticity".

5. Examples of exact solutions

5.1. Unsteady flow in a layer with one free boundary

Let us take v = v(n), vn = 0, then from (1.10) u(e,n) = -vne + ui(n) with some function ui(n). In such case the initial vorticity has the form w(e,n) = vnne - ui n. Because of ui(n) is arbitrary the system (1.7)-(1.9) admits 3D subalgebra (de,dx,dv). Since the variables t, n and y are invariant, we seek a partially invariant solution of rank 2 and defect 2 in the form x = x(e,n,t), y = y(n,t) and y = y(e,n,t). Inserting it into (1.9), we obtain

rn dn

x = a(n,t)e + b(n,t), y = (5.1)

Jo a(z,t)

with initial data

a(n, 0) = 1, at,(n, 0) = -vn, b(n, 0) = 0, bt(n, 0) = ui(n). (5.2)

The functions a(n,t), b(n,t) and y = y(-,n,t) can be found from equations (1.7) and (1.8) as

fd dr

= f (t)

r dT

1 — v

n lo f 2(t)J

6 = « 1 (n)f (t)J0 f 2 (t)

(5.3)

—2 rn y

y = (aat + vv)+ [(abt)v - uiv]- + J (btbv + y - v) dn,

where /(t) is arbitrary function satisfying the conditions f (0) = 1, ft (0) = 0. Then, from Weber's integral (1.11) and formulae (5.1), (5.3) we find the pressure

k(t)

p = l(n,t) x2(-,n) - gy(n,t), (5.4)

where

l(n,t) =l(t) - f dn, k(t) = /1 (5.5)

and l(t) is arbitrary function.

Solution (5.1)-(5.5) can be interpreted as the motion of plane layer (or a filled-up fluid layer) with a free boundary (or a moving rigid wall). Really, assume at the initial time there is a liquid

layer of thickness y = h0 = const. The lower rigid wall y = 0 is fixed and upper boundary y = h0

n

is a free one. The initial velocity field is given the formulae u = ui(n)- + u2(n), v = - J ui(n) dn.

o

The further motion is described by formulae (5.1)-(5.5). Of course, the outer pressure must be equals pout = l(t) - k(t)x2(-, n)/2 - gy(n,t) and the evolution of free boundary h(t) is given by the expression

1 fh 0 h(t) = ml

ft

1 - vn (n)Jo f^J

-1

dn. (5.6)

5.2. Example of Gerstner's waves

Let us take mapping (4.1) in the form

£ = a + k ekb sin(ka), n = b - k ekb cos(ka), J =1 - e2kb. (5.7)

It is one-to-one under conditions k > 0, b < b0 < 0, but there is not explicit dependence a, b on n. However, the formulas

(5.8)

X = a + 1 ekb sink (a + ct)], Y = b - 1 ekb cos[k(a + ct)],

kk

cekb

$ = —:— {sin[k(a + ct)] — sin(ka)} k

give us the exact solution of system (4.2)-(4.4) (k, c are constant). Physical vorticity is

= 2kcekb ^ = — 1 _ e2kb '

The fluid pressure can be found from integral (1.5)

p = g(bo - b) + 2 c2(e2kb - e2kb0), c2 = g/k. (5.10)

(5.9)

a

Formulas (5.7), (5.8), (5.10) determine the water waves called Gerstner's waves [12]. These waves are stationary in coordinate system which moves with velocity -c along the x-axis. Free boundary is a straight line b = b0 (in variables a, b) or it is trochoid with wave number 2n/k. Such waves are always vortex flow and vorticity has a maximum on free boundary and decreases exponentially with the fluid depth.

If b > 0 then mapping (5.7) is not one-to-one and such solution has not physical meaning. The case b0 = 0 gives us the cycloid equation

x = a + k sin(ka), y = b - cos(ka).

For ka = (2n + 1)n, n G N, this cycloid has a cusps.

In conclusion, Gerstner's waves are the invariant solution of system (4.2)-(4.4) with respect to two-dimensional subalgebra (c-18t -da -dX; d^) when J =1 - e2kb, f = a + k-1ekb sin(ka), g = b - k-1ekb cos(ka). New initial "velocities" have the form

U1 = ce2kb + cekb cos(ka), V1 = cekb sin(ka),

where c = (g/k)1/2. Notice that U1a + Vu = 0.

Conclusion

The equations describing 2D ideal liquid flows under potential force action are considered. Using Lagrangian coordinates this system is reduced to the system for particle trajectories and Weber's potential. Group classification with respect to the initial vorticity is performed. The equivalence group admitted by the governing equations is calculated. This group is infinite one and is used to reduce arbitrary element to the simpler form. As the example two exact solutions are constructed.

References

[1] L.V.Ovsyannikov, Group Properties of Differential Equations, Novosibirsk, Siberian Branch of USSR Ac. Sci., 1962 (in Russian).

[2] L.V.Ovsyannikov, Group Analysis of Differential Equations, Academic Press, New York, 1982.

[3] P.J.Olver, Application of Lie Groups to Differential Equations, Springer-Verlag, New York Inc., 1986.

[4] Andreev V. K., Kaptsov O. V., Pukhnachev V. V., Rodionov A. A., Applications of Group-Theoretical Methods in Hydrodynamics, Springer, Dordrecht, 2010.

[5] V.K.Andreev, V.V.Bublik, V.O.Bytev, Symmetries of Non-classical Hydrodynamics Models, Novosibirsk, Nauka, 2003 (in Russian).

[6] V.K.Andreev, Yu.A.Gaponenko, O.N.Goncharova, V.V.Pukhnachev, Mathematical Models of Convection, Walter de Gruyter GmbH & CO KG, Berlin-Boston, 2012.

[7] V.K.Andreev, Stability of Unsteady Motions of a Fluid with a Free Boundary, Novosibirsk, Nauka, 1992 (in Russian).

[8] V.I.Nalimov, The Cauchy-Poisson Problem, Dinamica Sploshnoy Sredy, 1, 1974 (in Russian).

[9] H.Yosihara, Gravity Waves on the Free Surface of an Incompressible Perfect Fluid of Finite Depth, Publ. RIMS Kyoto Univ., Vol. 18, 1982.

[10] H.Yosihara, Capillary-gravity Waves for an Incompressible Ideal Fluid, J. Math. Kyoto Univ., (1983), 23, 649-694.

[11] V.I.Nalimov, Unsteady Vortex Surface Waves, Siberian Math. J., 37(1996), no. 6, 1189-1198 (in Russian).

[12] H.Lamb, Hydrodynamics, Moscow, Gostekhizdat, 1947 (in Russian).

[13] V.K.Andreev, A.A.Rodionov, Group-Theoretic Classification and Exact Solutions of equations for a Plane and Rotationally-Symmetric Flow of an Ideal Fluid in Lagrangian Coordinates, Differential Equations, 24(1988), no. 9, 1041-1049 (in Russian).

[14] V.K.Andreev, A.A.Rodionov, Group Analysis of Equations for Planar Flows of an Ideal Fluid in Lagrangian Coordinates, Dokl. Akad. Nauk SSSR, 298(1988), no. 6, (in Russian).

Симметрииныи анализ уравнении идеальном жидкости в терминах траектории — потенциал Вебера

Виктор К. Андреев

Институт вычислительного моделирования СО РАН Академгородок, 50/44, Красноярск, 660036 Институт математики и фундаментальной информатики Сибирский федеральный университет Свободный, 79, Красноярск, 660041

Россия

Дарья А. Краснова

Институт математики и фундаментальной информатики Сибирский федеральный университет Свободный, 79, Красноярск, 660041

Россия

Рассматриваются уравнения двумерных движений идеальной жидкости в координатах Лагран-жа. Для потенциальных внешних сил они имеют общий интеграл Вебера, причем новая система включает в себя начальные данные. Это делает актуальной задачу групповой классификации. Установлено, что основная группа непрерывных преобразований является бесконечномерной по пространственным координатам. Найдены специальные зависимости начальной завихренности, при которых происходит 'расширение группы. Кроме того, изучены групповые свойства исходной системы в произвольных лагранжевых координатах и найдены точные решения.

Ключевые слова: уравнения Эйлера, симметрийный анализ, преобразование Вебера, преобразование эквивалентности, групповая классификация.