Integrodifferential Hyperbolic Equations and its Application for 2-D - - PowerPoint PPT Presentation

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14, 2006

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Outline

• 1. Integrodifferential hyperbolic models.
• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

2.1. Mathematical model. 2.2. Hyperbolicity conditions. 2.3. Riemann problem. 2.4. The equation of the two-layer flow with constant vorticity. 2.5. Wave of flows interaction. Particles trajectories. 2.6. Discontinuous solutions. 2.7. Conservation laws. General case.

• 3. Exact solutions for the shear shallow water equations.

3.1. Shallow water equations for shear flows. 3.2. Traveling waves. 3.3. The sequences of the smooth solutions converging to the discontinuous. 3.4. Self-similar solutions.

• 4. Conclusions.

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 1. Integrodifferential hyperbolic models.
• 1. Integrodifferential hyperbolic models.

A number of problems of theoretical hydrodynamics can be deduced to the integrodifferential equations which can be written as Ut + A(U)Ux = 0 (1) Here U(t, x, λ) - desired vector and A(U) is a nonlocal matrix operator acting

• ver the variable λ. For example,

1) long-wave models for 2-D rotational (shear) flows of an ideal homogeneous and barotropic fluid with free surface and in channels; 2) Vlasov type 1-D kinetic equations, in particular bubbly flow kinetic models (written in Euler-Lagrangian coordinates) are belong to this class. Theoretical analysis of such systems is based on a generalization of the concept

• f hyperbolicity and characteristics for equations with operator coefficients

suggested by V.M. Teshukov. It was established that integrodifferential models contain both discrete and continuous spectra of the characteristics velocities. Hiperbolicity conditions for this type equations were formulated. The theory of simple waves was developed.

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 1. Integrodifferential hyperbolic models.

However, general theory of integrodifferential equations is under development now, and great number of the problems is not solved including development of numerical schemes, theory of discontinuous solutions, etc. This paper deals with models of type (1) describing in long-wave approximation plane-parallel shear flows of a perfect fluid. There are two main goals:

• 1. To study characteristic properties and to solve the Riemann problem

for the equations describing shear plane-parallel flows of an ideal incompressible liquid in a narrow channel;

• 2. To construct the sequences of the exact smooth solutions converging

to discontinues for the shear shallow water equations.

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.
• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

2.1. Mathematical model. The Euler equations and boundary conditions in dimensionless form (an ideal incompressible fluid in a narrow straight channel): ut + uux + vuy + ρ−1px = 0, ε2(vt + uvx + vvy) + ρ−1py = −g, ux + vy = 0, v(t, x, 0) = 0, v(t, x, h0) = 0. Long-wave approximation ε = 0: ut + uux + vuy + p∗

x = 0,

v = − y ux dy′, h0 ux dy = 0. (2)

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

Note that flow rate Q in the channel does not depend on the cross section Q(t) = h0 u(t, x, y) dy. Using equations for velocities at the walls uit + uiuix + p∗

x = 0,

(u0 = u(t, x, 0), u1 = u(t, x, 1)) we can eliminate p∗(t, x) (pressure at the upper boundary of the channel). As result we have evolutionary equation for relative velocity w = u − u1: wt + (w2/2 + wu1)x + vwy = 0, v = −yu1x− y wx(t, x, y′) dy′, u1 = 1 h0

• Q(t)−

h0 w(t, x, y) dy

• .

We pass to the Euler-Lagrangian coordinates (x, λ) by the substitution y = Φ(t, x, λ): Φt + u(t, x, Φ)Φx = v(t, x, Φ), Φ(0, x, λ) = Φ0(x, λ).

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

As result, we obtain following system which belong to the class (1): wt + (w2/2 + wu1)x = 0, Ht + ((w + u1)H)x = 0, (3) u1 = 1 h0

• Q(t) −

1 wH dλ

• Here w = u − u1, H = Φλ, λ ∈ [0, 1]. The equation ht = 0, which

expresses the fact that the upper boundary is fixed, is a consequence of the system. 2.2. Hyperbolicity conditions. The characteristics of the system (1) is determined by the equation x′(t) = k(t, x), where k is eigenvalue of the operator A∗. System (1) is hyperbolic if all the eigenvalues k are real and set of equations for the characteristics (F, Ut + AUx) = 0 (F is eigenfunctional acting on functions of the variable λ) is equivalent to equations (1)

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

Equation defining characteristic velocity takes form: χ(k) = − 1 ω1(u1 − k) + 1 ω0(u0 − k) + 1

• 1

ω

• λ

dλ u − k = 0. Integrodifferential model (3) has continuous characteristic spectrum kλ(t, x) = u(t, x, λ). Hiperbolicity conditions are formulated in terms of the limiting values χ± of the complex function χ(z) from upper and lower half-planes on the real axis. For flows with monotonic velocity depth profile uλ = 0, conditions △arg(χ+(u)/χ−(u)) = 0, χ±(u) = 0 are necessary and sufficient for Eq. (3) to be hyperbolic (functions u and H are sufficiently smooth, ω = uλ/H, and the argument increment is calculated for u varied from u0 to u1).

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

2.3. Riemann problem. As was shown above, the solution of the equation wt + (w2/2 + wu1)x + vwy = 0, (4) v = −yu1x − y wx(t, x, y′) dy′, u1 = − 1 h0 h0 w(t, x, y) dy. describes plane-parallel shear motions of an ideal fluid in a narrow channel (Q = 0). Equation (4) admits particular solutions of the form w = w(y), that correspond to steady shear flows. A natural generalization of the formulation of the Riemann problem in the case of Eq. (4) is the Cauchy problem w|t=0 = ur(y) − ur(h0), x > 0 ul(y) − ul(h0), x < 0

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

2.4. The equation of the two-layer flow with constant vorticity.

2.4. The equation of the two-layer flow with constant vorticity. First, we consider special initial data ur(y) = ω1y + ur

0,

ul(y) = ω2y + ul

0,

which correspond to the interaction of flows with constant vorticities. In this case the solution of the Eqs. (2) takes form u(t, x, y) = Ω1y + u0(t, x), 0 ≤ y ≤ h(t, x) Ω2(y − h0) + u1(t, x), h(t, x) ≤ y ≤ h0 v(t, x, y) = −yu0x, 0 ≤ y ≤ h(t, x) (h − y)u1x − hu0x, h(t, x) ≤ y ≤ h0 Here Ω1 = ω1, Ω2 = ω2 if ωi > 0 (figure ”a”) and Ω1 = ω2, Ω2 = ω1 if ωi < 0 (figure ”b”). The form of the solution and condition that the velocity vector is continuous on the boundary of the layers y = h(t, x) allow us to obtain following

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

2.4. The equation of the two-layer flow with constant vorticity.

scalar conservation law ∂h ∂t + ∂ϕ ∂x = 0, ϕ(h) = Ω1 − Ω2 2h0 h3 +

• Ω2 − Ω1

2

• h2 − Ω2h0

2

• h. (5)

with initial data h(0, x) =

• r

h0, x < 0 h0

• r

0, x > 0 The function ϕ(h) is convex if 2−1 < α0 < 2, (α0 = Ω1/Ω2). In this case continuous solution (simple wave of flows interaction) exists. Since the equations and boundary conditions are invariant with respect to uniform stretching of the variables t and x, we seek self-similar solution h = h(k), k = x/t. Let 0 < ω2 < ω1, 1 < ω1/ω2 < 2 (we should take Ω1 = ω1, Ω2 = ω2). Integrating Eq. (5) we obtain h(k) = −(2ω2 − ω1)h0 3(ω1 − ω2) +

• (2ω2 − ω1)2h2

9(ω1 − ω2)2 + h0(2k + ω2h0) 3(ω1 − ω2)

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

2.5. Wave of flows interaction. Particles trajectories.

2.5. Wave of flows interaction. Particles trajectories. Liquid particles in the flow move along the trajectories x = x(t) and y = y(t): ˙ x = u(t, x, y), ˙ y = v(t, x, y). It is convenient to use variables (k, y). Then we have equations dk ds = u(k, y)−k, dy ds = V (k, y)

• s = ln t

t0 , V = tv, k = x t

• The singular points (where u(k, y) − k = 0, V (k, y) = 0) of the system are

the points A and C. Let us separate out the linear part of the system in a neighborhood of the point A (C) and calculate the eigenvalues of the coefficient

• matrix. This analysis shows that A is a saddle and C is a node (figure ”a”).

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

2.5. Wave of flows interaction. Particles trajectories.

It is evident that the flow in the region of interaction of the flows has a substantially 2-D nature. Along the line of contact of the two rotating flows, jet flows is formed that is directed to the upper wall or the bottom, depending on the ratio of the vorticities.

Figure: Relative velocity field (u − k, V ). Figure: Velocity field (u, V ).

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

2.6. Discontinuous solutions.

2.6. Discontinuous solutions. If condition 2−1 < α0 < 2, (α0 = Ω1/Ω2) is fell then smooth solution of the equation (5) does not exists. To construct solution of the equation ht + ϕ(h)x = 0 we should combine simple wave and shock.

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

2.6. Discontinuous solutions.

Oleinik conditions ϕ(h) − ϕ(h−) h − h− ≥ D (h− < h < h+) ϕ(h) − ϕ(h+) h − h+ ≤ D (h+ < h < h−) Possible configuration of the simple wave and shock are shown. Let 0 < ω2 < ω1, ω1/ω2 > 2. Discontinuous solution h(t, x) =    0, x < Dt ¯ h(k), Dt ≤ x ≤ ϕ′(h0)t h0, x > ϕ′(h0)t

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

2.6. Discontinuous solutions.

Particles trajectories in the discontinuous solution (0 < ω2 < ω1, ω1/ω2 > 2).

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

2.7. Conservation laws. General case.

2.7. Conservation laws. General case. The Riemann problem for arbitrary monotonic velocity profiles ul(y) and ur(y) can be solve numerically. We propose to use the conservation laws (y = Φ(t, x, λ), H = Φλ) wt + (w2/2 + wu1)x = 0, Ht + (uH)x = 0, u1 = − 1 H dλ −1 1 wH dλ, u = w + u1 Then we divide interval [0, 1] into N intervals 0 = λ0 < λ1 < ... < λN = 1 and introduce the designations yi = Φ(t, x, λi), ηi = yi − yi−1, ui = u(t, x, λi) ωi = (ui − ui−1)/ηi, uci = (ui + ui−1)/2. (such approach was used by Teshukov, Russo, Chesnokov (2004) for shear shallow water equation). Taking into account the equality H dλ and using the piecewise linear approximation for horizontal velocity u = ωi(y − yi−1) + ui−1, y ∈ [yi−1, yi]

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 2. Unsteady interaction of uniformly shear flows in a narrow channel.

2.7. Conservation laws. General case.

we arrive to the differential conservation laws. In this case unknown vector (η1, ..., ηN, ω1η1, ..., ωNηN) satisfies the following 2 × 2N system ∂ηi ∂t + ∂ ∂x(uciηi) = 0, ∂ ∂t(ωiηi) + ∂ ∂x(uciωiηi) = 0, uci = u0 + ωiηi 2 +

i−1

• j=1

ωjηj, u0 = − N

• i=1

ηi −1 ×

N

• i=1

ηi

• ωiηi

2 +

i−1

• j=1

ωjηj

• .

Equation (η1 + ...ηN)t = 0 (which means that depth of the channel is equal to h0 = const) is the sequence of the system of equations.

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 3. Exact solutions for the shear shallow water equations.
• 3. Exact solutions for the shear shallow water equations.

3.1. Shallow water equations for shear flows. The solution of the problem ut + uux + vuy + px = 0, ε2(vt + uvx + vvy) + py = −g, ux + vy = 0; ht + u(t, x, h)hx = v(t, x, h), v(t, x, 0) = 0, p(t, x, h) = p0; (ε = H0/L0 → 0). describes plane-parallel free-boundary flows of an ideal incompressible fluid in a gravity field (dimensionless variables).

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 3. Exact solutions for the shear shallow water equations.

Ignoring terms of order ε2, we can find pressure p = g(h − y) + p0. After transformations, one obtains the equations for unknown u(t, x, y) and h(t, x) describing shear shallow water flows: ut+uux+vuy +ghx = 0, ht+ h u dy

• x

= 0, v = − y ux dy′. We formulate the equations using Euler-Lagrangian coordinates: y = Φ(t, x, λ), (0 ≤ λ ≤ 1) Φt + u(t, x, Φ)Φx = v(t, x, Φ), Φ(0, x, λ) = λh(0, x). Substitution is reversible if Φλ = 0. The functions u(t, x, λ) and H(t, x, λ) = Φλ satisfy the system ut + uux + g 1 Hx dλ = 0, Ht + (uH)x = 0. (6)

• Eqs. (6) belong to the class (1).

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 3. Exact solutions for the shear shallow water equations.

The dimensionless vorticity of the flow is equal to ω = −uy + ε2vx satisfies the equation ωt + uωx + vωy = 0. In the long-wave theory ω = uy and also conserve along the trajectories. Vlasov-like formulation can be obtain for flows with ω = 0: Wt + uWx − ghxWu = 0, h = u1

u0

W du, u1t + u1u1x + ghx = 0, u0t + u0u0x + ghx = 0. (7) Here W = ω−1 (analog of the distribution function) is considered as an unknown function depending on t, x, u; functions u0(t, x) and u1(t, x) are horizontal velocities at the bottom and at free surface, respectively. It follows from equation (7) that the quantity W (t, x, u) is conserved along the trajectories dx dt = u, du dt = −ghx.

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 3. Exact solutions for the shear shallow water equations.

3.2. Traveling waves.

3.2. Traveling waves. A solution of the form W = W (ζ, u), ui = ui(ζ), where ζ = x − Dt, describes a travelling wave propagating at constant velocity D. For this class of solutions, Eq. (7) takes the form (u − D)Wζ − gh′(ζ)Wu = 0, h = u1

u0

W du, (ui − D)u′

i(ζ) + gh′(ζ) = 0,

(i = 0, 1). Previous equations are integrated: W = Φ(η), η = u2 − 2Du + 2gh, u0 ≤ u ≤ u1, u2

i − 2Dui + 2gh = ηi.

Let us consider the Cauchy problem ζ = ζ0 : W = W0(u), u0 = u00, u1 = u10

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 3. Exact solutions for the shear shallow water equations.

3.2. Traveling waves.

The solution takes constant values on the characteristics η = (u − D)2 + 2gh − D2 = const Designations: η0 = 2gh0 − D2, η00 = (u00 − D)2 − 2gh0 − D2, η10 = (u10 − D)2 − 2gh0 − D2; (u00 < D < u10). The solution of the Cauchy problem is uniquely determined from initial data in the domains Ω1 and Ω2. Ω1 = {(u, h) | 0 < h < h0, u0(h) < u < u∗(h)} Ω2 = {(u, h) | 0 < h < h0, u∗(h) < u < u1(h)} Ω3 = {(u, h) | 0 < h < h0, u∗(h) < u < u∗(h)}

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 3. Exact solutions for the shear shallow water equations.

3.2. Traveling waves.

Here u0(h) = D −

• 2g(h0 − h) + (u00 − D)2

u1(h) = D +

• 2g(h0 − h) + (u10 − D)2

u∗(h) = D −

• 2g(h0 − h),

u∗(h) = D +

• 2g(h0 − h).

We determine the functions Φ∓(η) = W0(D ∓

• η − 2gh0 + D2),

Using functions Φ±, we obtain the following solution in the domains Ω1 and Ω2 W = Φ−(u2 − 2Du + 2gh), W = Φ+(u2 − 2Du + 2gh). To construct the solution in the domain Ω3, we transform the relation h = u1

u0

W du, W = Φ(η); s = 2gh − D2 to an integral equation for the function Φ(η):

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 3. Exact solutions for the shear shallow water equations.

3.2. Traveling waves.

η0

• s

Φ(η) dη √η − s = F (s), (s = 2gh − D2) F (s) = s + D2 2g − 1 2

η00

• η0

Φ−

0 (η) dη

√η − s − 1 2

η10

• η0

Φ+

0 (η) dη

√η − s . The solution of the Abel equation yields Φ(η) = 1 π

• F (η0)

√η0 − η −

η0

• η

F ′(s) ds √s − η

• 3.3. The sequences of the smooth solutions converging to the discontinuous.

Traveling wave solution involve arbitrary function h = h(x − Dt). Form the sequence smooth functions hn converging to the discontinuous solution. Let we have following initial data W (ζ0, u) = W0(u) = ω−1 =const, u0(ζ0) = u00, u1(ζ0) = u10 (u00 < D < u10), corresponding to flow u = ωy + u00, v = 0, h = (u10 − u00)/ω.

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 3. Exact solutions for the shear shallow water equations.

3.2. Traveling waves.

Then Φ−(η) = Φ+(η) = ω−1, Φ(η) = 1 ω − 1 πω

• arcsin
• η0 − η

η00 − η + arcsin

• η0 − η

η10 − η

• − 1

πg

• η0 − η.

The solution W (h, u) =    ω−1, u0(h) ≤ u ≤ u∗(h) Φ(u2 − 2Du + 2gh), u∗(h) < u < u∗(h) ω−1, u∗(h) ≤ u ≤ u1(h) Here u0,1(h) = D ∓

• 2g(h0 − h) + (u00 − D)2

u∗(h) = D −

• 2g(h0 − h),

u∗(h) = D +

• 2g(h0 − h).

Let us take D = 0, h0 = 0, hmin = 0.55, hn(x) = 2−1(h0 − hmin)(2π−1arctg(nx) + 1) + hmin.

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 3. Exact solutions for the shear shallow water equations.

3.2. Traveling waves.

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 3. Exact solutions for the shear shallow water equations.

3.2. Traveling waves.

Streamlines and velocity field

Figure: Velocity field (u, v) in the travel- ing wave. Figure: Streamlines of the flow.

Flow with critical layer is considered (there is line, where fluid velocity coincides with traveling wave velocity).

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 3. Exact solutions for the shear shallow water equations.

3.2. Traveling waves.

Flows without critical layer (D < u00 < u10 or u00 < u10 < D): Streamlines and velocity field

Figure: Velocity field (u, v). Figure: Streamlines in the traveling wave.

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 3. Exact solutions for the shear shallow water equations.

3.4. Self-similar solutions.

3.4. Self-similar solutions. This algorithm also works for construction self-similar solutions (at least in the case, when critical layer is absent) W = W (ξ, u), u0 = u0(ξ), u1 = u1(ξ), ξ = x/t : (u − ξ)Wξ − gh′Wu = 0, (ui − ξ)u′

i + gh′ = 0,

h = u1

u0

W du The solution takes constant values along the integral curves of the equation du/dξ = −gh′(ξ)/(u − ξ)

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

• 4. Conclusions.
• 4. Conclusions.

The Riemann problem (the problem of decay of an arbitrary discontinuity) for the system of equationd describing shear plane-parallel flows of an ideal incompressible liquid in a channel is solved in a long-wave approximation. Exact solutions of the shallow water equations for shear flows for the class of traveling waves and simple waves are derived. The sequence of smooth solutions converging to discontinuous solution is constructed.