arxiv: v1 [math.ap] 3 Feb 2016

Transcription

1 Interactions between discontinuities for binary miture separation problem and hodograph method arxiv:62.463v [math.ap] 3 Feb 26 M S Elaeva, E V Shiryaeva 2, M Yu Zhukov 3, 4 Faculty of Applied Mathematics and Computer Science, Financial University under the Government of the Russian Federation, 49 Leningradsky prospekt, 25993, Moscow, Russia mselaeva@fa.ru 2 Institute of Mathematics, Mechanics and Computer Science, Southern Federal University, 5/42 Bolshaya Sadovaya Street, 3446, Rostov-on-Don, Russia shir@math.sdedu.ru 3 Institute of Mathematics, Mechanics and Computer Science, Southern Federal University, 5/42 Bolshaya Sadovaya Street, 3446, Rostov-on-Don, Russia myuzhukov@gmail.com 4 South Mathematical Institute, Vladikavkaz Center of RAS, 22 Markus Street, 36227, Vladikavkaz, Russia February 26 Abstract. The Cauchy problem for first-order PDE with the initial data which have a piecewise discontinuities localized in different spatial points is completely solved. The interactions between discontinuities arising after breakup of initial discontinuities are studied with the help of the hodograph method. The solution is constructed in analytical implicit form. To recovery the eplicit form of solution we propose the transformation of the PDEs into some ODEs on the level lines (isochrones) of implicit solution. In particular, this method allows us to solve the Goursat problem with initial data on characteristics. The paper describes a specific problem for zone electrophoresis (method of the miture separation). However, the method proposed allows to solve any system of two first-order quasilinear PDEs for which the second order linear PDE, arising after the hodograph transformation, has the Riemann-Green function in eplicit form. PACS numbers: 2.3.Ik, 2.3.Jr, 82.8.Bg, 87.5.Tt AMS classification scheme numbers: 35L, 35L67, 35L4, 35L45 35L5, 35L65 Keywords: hyperbolic quasilinear equations, hodograph method, zone electrophoresis Submitted to: J. Phys. A: Math. Theor.

2 Interactions between discontinuities and hodograph method 2. Introduction The paper is devoted to the investigation of the discontinuities interactions for the firstorder quasilinear hyperbolic PDEs. We consider the Cauchy problem for systems of first-order hyperbolic quasilinear equations with discontinuous piecewise constant initial data and additional the Rankine Hugoniot conditions on the discontinuities u k t + ϕ k (u, u 2 ) =, u k (, ) = f k (), D[u k ] = [ϕ k ], k =, 2, () where f k () is the piecewise constant function, D is the discontinuity velocity, [ ] is the notation for jump of function across discontinuity. When the discontinuities of the functions f k () are localized only at one point the problem () is well-known as the Riemann problem. In this case, the Cauchy problem solution is a moving sequence of shock waves (strong discontinuities S ) and rarefaction waves (weak discontinuities W ). We discuss a slightly more general problem in the case when f k () have discontinuities at two points and 2. Obviously, in this case, a moving shock wave and rarefaction wave arising at the points and 2 are interacting with each other in the process of the solution evolution. Four types of interactions are possible: (i) a shock wave S with a shock wave S 2, (ii) a shock wave S with a rarefaction wave W 2, (iii) a rarefaction wave W with a shock wave S 2, (iv) a rarefaction wave W with a rarefaction wave W 2. Studies have shown, that the case of interaction of two rarefaction waves is the most difficult to describe. Indeed, the interaction of weak discontinuities there is the Goursat problem with the initial data on the characteristics because the fronts of the rarefaction waves (weak discontinuities) move along the characteristics. In the case of a single discontinuity at the point the different ways of constructing the Goursat problem solution (and the Riemann problem) for two quasilinear equations are presented, for instance, in []. These methods in general case are based on iterative method of the solution constructing. For one particular quasilinear equations, e.g., equations of isotachophoresis and chromatography, it is possible to obtain an analytical solution in the case when there are several discontinuities of the function f(). In the case of two quasilinear equations, the problem of discontinuities interaction (S, S 2 ), (S, W 2 ), (W, S 2 ) for the equations of electrophoresis is completely solved in [2]. To date, the problem of interaction type (W, W 2 ) was an open problem. To solve problem of interaction between discontinuities we propose to use the hodograph method for system of two quasilinear hyperbolic equations written in the Riemann invariants. We recall that in this method the dependent and independent variables reverse roles ((u, u 2 ) (t, )). In this case, the solution is an implicit form t = t(u, u 2 ), = (u, u 2 ). To determine, for eample, the function t(u, u 2 ) there is a system of first-order linear PDEs with variable coefficients (see, e.g., []). The solvability condition allows us to transform the system into a single second-order linear PDE with variable coefficients. For this equation the Riemann Green function is constructed (see, e.g., [3]) and then the Goursat problem solution is obtained (see e.g. [4]).

3 Interactions between discontinuities and hodograph method 3 In a recent paper [5] a variant of the hodograph method based on the presence of conservation laws for the Riemann invariants is proposed. This method allows us to obtain effectively an implicit solution in the form t = t(u, u 2 ), = (u, u 2 ) for the case when there is an eplicit form of the Riemann Green function, for eample, for the shallow water equations, chromatography equations (see e.g. [,6]), zone electrophoresis equations (see e.g. [7]), etc. Moreover, the solution of the Cauchy problem can be obtained for arbitrary initial data including discontinuities. In the papers [8] the method [5] was only slightly modified. In addition, the numerically-analytical method proposed allows to reconstruct the eplicit solution u = u (, t), u 2 = u 2 (, t) with the help of its implicit form t = t(u, u 2 ), = (u, u 2 ). This allows us to transform the Cauchy problem for PDE to the Cauchy problem for ODEs on the level lines of an implicit solution. In fact, the calculations are required only for integration of the Cauchy problem for ODEs. In this paper the numerical part of the method is almost not required. Most of the relations describing the solution are obtained in an eplicit analytical form. We note that the problem of the discontinuities interaction has important practical application, in particular, for the method of electrophoresis. Electrophoresis is a separation method of the multicomponent miture into individual components under action of an electric field. A diffusionless variant of the electrophoresis mathematical models is the conservation laws for hyperbolic equations. The discontinuous initial data correspond to the boundaries of the spatial regions occupied by the individual components. The study of the discontinuities behavior and their interactions allows us to understand the comple scenario of components motion for the separation process. The general formulation of the zone electrophoresis problem is given in [7]. In the case of the two-component miture the interactions (S, S 2 ), (S, W 2 ), (W, S 2 ) are studied in [2, 9 ]. The results presented below indicate the analytical method of the solving for the problem of the rarefaction waves interaction. This gives us a complete analytical solution of the problem (together with the results of [2, 9 ]). Note that the method proposed is weak linked with the electrophoresis problem. This method can be used for any system of two hyperbolic quasilinear equations. However, the application of the method requires eplicit analytic form of the Riemann Green function for the linear equation arising from the use of the hodograph method. The paper is organized as follows. In section 2 the basic electrophoresis equations, written in the original variables (concentrations of components) and the Riemann invariants, are presented. In section 3 we formulate the Riemann problem for the case of initial data with discontinuities at different points in space. The Cauchy problem with initial data in the moment of interaction of weak discontinuities is presented. In section 4 the hodograph method is described. In section 5 the method of recovering the eplicit form of the solution from its implicit form is described. In sections 6 8 different variants of the discontinuities interactions are considered. In Appendi A the brief description of numerical-analytical method for the solving of the Cauchy problem is presented.

4 Interactions between discontinuities and hodograph method 4 2. Zone Electrophoresis Equations To describe zone electrophoresis of binary miture in diffusionless approimation we use two first-order quasilinear hyperbolic equations [2, 7, 9 ] u k t + ( ) µ k u k =, k =, 2, s = u + u 2 >, (2) + s where u k are the effective concentrations, µ k are the component mobilities, ( + s) is the conductivity of the miture, E = ( + s) is the intensity of the electric field. System (2) can be written in the Riemann invariants R k t + λ k (R, R 2 )R k =, k =, 2, (3) λ (R, R 2 ) = R R R 2, λ 2 (R, R 2 ) = R 2 R R 2. (4) The concentrations u, u 2 are connected to the Riemann invariants R, R 2 as u = µ2 (R µ )(R 2 µ ) R R 2 (µ µ 2 ) The roots of the polynomial, u 2 = µ (R µ 2 )(R 2 µ 2 ). (5) R R 2 (µ 2 µ ) P (R, R 2 ) ( + u + u 2 )(R) 2 (µ + µ 2 + u µ 2 + u 2 µ )R + µ µ 2 allows us to get inverse relations R = R (u, u 2 ), R 2 = R 2 (u, u 2 ). 3. The analogue of the Riemann problem We study the equations (3), (4) with the discontinuous initial data where R t= = R (), R 2 t= = R 2 (), (6) µ, <, R() = q, < < 2, µ, 2 <, Here, q, q 2 are given constants such that µ 2, <, R() 2 = q 2, < < 2, µ 2, 2 <. < q < µ < µ 2 < q 2. (8) Note that for the classical Riemann problem the initial data usually have the discontinuities at only one point (not and 2 ). We also note that for original equations (2) the initial data (6), (7) correspond to zero concentrations outside the interval (, 2 ). Inside the interval (, 2 ) the initial concentrations are constant u = u, u 2 = u 2. The values u, u 2 are determined with the help q, q 2, and (5). Inequality (8) means that the effective concentrations u > and u 2 <. To solve the Riemann problem (3), (4), (6) (8) we add the Rankine Hugoniot conditions which correspond to conservation laws (2). Taking into account (5) we rewrite these conditions for the Riemann invariants in the following form [ ] (µ k R )(µ k R 2 ) D = [ (µ k R )(µ k R 2 ) ], k =, 2, (9) µ k R R 2 (7)

5 Interactions between discontinuities and hodograph method 5 where D is the discontinuity velocity, [ ] is the notation for jump of function across discontinuity. In the case when the solution has strong discontinuities such solution should additionally satisfy to the La conditions. These conditions are well known and we omit them (see e.g. [, 2, 3] and also [2, 7, 9 ]). For completeness, we present the results obtained in [2, 9 ]). The equations (2) (and (3), (4)) have the hyperbolic type. In particular, it means that the breakup of the initial discontinuities at the points and 2 can be analyzed independently. In the initial moment t = + at the point and also at point 2 arise the rarefaction waves (self-similar solutions, weak discontinuities) and shock waves (discontinuous solutions, strong discontinuities). After a time, strong and weak discontinuities of the solution arising at point interact with strong and weak discontinuities of the solution arising at point 2. We illustrate the evolution of solution using (, t)-plane and zones on this plane. It is convenient to introduce the notation Z k = {R k(, t), R 2 k(, t); a k (t), b k (t)}, where Z k is the zone on (, t)-plane, a k (t), b k (t) are the zone boundaries at the moment t, Rk (, t) and R2 k (, t) are the Riemann invariants in this zone. For initial data (7) the evolution of zones on (, t)-plane is shown in figure. t T 6 T 3 s Z 5 Z 2 ϕ Z 7 Z 6 θ 2 s T int Z 2 r 2 l Z 4 Z 3 r Z 8 l 2 X int X 3 X 6 Figure. The zones in the (, t) plane. Weak discontinuities are solid lines, strong discontinuities dashed lines. The fronts of the rarefaction waves and the shock waves are shown by solid and dashed lines, respectively. Figure corresponds to some finite time interval (, T 3 ). It is obvious that the intersection of the boundaries corresponds to the interaction of discontinuities. There are only four types of interactions: (i) the interaction of two strong discontinuities, (ii, iii) the interaction of strong and weak discontinuities, (iv) the interaction of two weak discontinuities. For instance, the point (X int, T int ) is interaction of two weak discontinuities 2 r and l. We recall that the waves arising at the breakup of initial discontinuity can be classified using the eigenvalues λ k. The wave corresponding to the eigenvalue λ k is

6 Interactions between discontinuities and hodograph method 6 called k-wave. Thus, the line s corresponds to the -shock wave, the lines 2 l and 2 r correspond to the left and right fronts of 2-rarefaction wave, etc. (see figure ). The solution of the problem (3), (4), (6) (8) in time interval (, T int ) is described with the help of zones (see figure ) Z = {µ, µ 2 ;, s(t)}, Z 2 = {q, µ 2 ; s(t), 2 l (t)}, () Z 3 = {q, R 2 a(z 2 ); 2 l (t), 2 r(t)}, Z 6 = {R a(z ), q 2 ; l (t), r(t)}, Z 4 = {q, q 2 ; 2 r(t), l (t)}, Z 7 = {µ, q 2 ; r(t), 2 s(t)}, where Z 8 = {µ, µ 2 ; 2 s(t), + }, R a(z ) = z q 2, z = 2, R 2 t a(z 2 ) = z 2 q, z2 =. () t The motion laws of the left and right rarefaction waves fronts (weak discontinuities) have the following form l (t) = 2 + q q q 2 t, r(t) = 2 + µ µ q 2 t, (2) 2 l (t) = + q µ 2 µ 2 t, 2 r(t) = + q q 2 q 2 t. Strong discontinuities = s(t), = 2 s(t) (shock waves) are moved with the velocities D = q µ µ 2, D 2 = µ µ 2 q 2 as s(t) = + q µ µ 2 t, 2 s(t) = 2 + µ µ 2 q 2 t. (3) As already mentioned, the solution () (3) of problem (3), (4), (6) (8) is valid for the time interval (, T int ). In the time moment t = T int the right front of 2-rarefaction wave = 2 r(t) reaches the left front of -rarefaction wave = l (t). This moment corresponds to interaction of waves. The coordinates of the point (X int, T int ) on (, t)- plane (see figure ) can be calculated as the solution of the equation 2 r(t) = l (t) T int = 2 q q 2 (q 2 q ), X int = q 2 q 2 q q 2. For completeness, the solutions of the original problem (2) for concentrations u k at moments t = T int and t = t < T int are presented in figure 2 in the case of the parameters q = 2, µ = 5, µ 2 = 8, q 2 =, =, 2 =. (4) As a result of interaction between two rarefaction waves the zone Z 4 is disappeared. New zone Z 5 is appeared and the boundaries of the zones Z 3, Z 6 are modified (but the Riemann invariants in these zones are saved the same) Z 5 = {R (, t), R 2 (, t); ϕ(t), θ(t)}, Z 3 = {q, R 2 a(z 2 )); 2 l (t), ϕ(t)}, Z 6 = {R a(z 2 ), q 2, θ(t), r(t)}.

7 Interactions between discontinuities and hodograph method 7 Here ϕ(t), θ(t) are functions that determine new boundaries of the zone. Z Z 6 Z 7 Z 8 u k Z 2 Z 3 Z 4 u k Z 2 Z 3 Z Z 6 Z 7 Z s s 2 l 2 r l r 2 s 2 l 2 r r 2 s Figure 2. Concentrations u k at the moments t =., t = T int =.25 Obviously, the lines = ϕ(t), = θ(t) are weak discontinuities that are identified by corresponding characteristic directions λ and λ 2. On these lines the Riemann invariants R k have the weak discontinuity. To determine, for eample ϕ(t), we have the Cauchy problem dϕ dt = λ (q, Ra), 2 ϕ(t int ) = X int. (5) Integrating (5) we get ( ) ϕ /2 ( ) ( ) = q 3/2 t /2 T /2 int + ( X int ) /2. (6) Similarly, the line of weak discontinuity = θ(t) is determined by formula ( ) θ 2 /2 ( ) ( ) = q 2 3/2 t /2 T /2 int + ( X int 2) /2. (7) Easily to describe the behavior of the lines = ϕ(t), = θ(t) on the (, t)-plane. The zone Z 3 (see figure ) disappears at the moment t = T 3 in the point = X 3 when characteristic = ϕ(t) intersects characteristic = 2 l (t). The point (X 3, T 3 ) is determined by the relation X 3 = ϕ(t 3 ) = 2 l (T 3) T 3 = (q2 q ) 2 (q µ 2 ) 2 T int, X 3 = + q µ 2 µ 2 T 3. The zone Z 6 (see figure ) disappears at the moment t = T 6 in the point = X 6, when characteristic = θ(t) intersects characteristic = r(t). Point (X 6, T 6 ) is determined by the relation X 6 = θ(t 6 ) = r(t 6 ) T 6 = (q2 q ) 2 (q 2 µ ) 2 T int, X 6 = 2 + µ µ q 2 T 6. The results listed in this section are given in one form or another in the papers [2, 9 ]. In particular, these papers indicate that further study of interactions between

8 Interactions between discontinuities and hodograph method 8 waves requires the solution of the Goursat problem with discontinuous initial data on characteristics. Primarily, it concerns to the interaction between weak discontinuities. An effort to solve a similar problem with the help of the numerical methods for the modified diffusive approimation has shown that the numerical results do not allow us to describe adequately the interaction of weak discontinuities (see []). In the following sections we show that the Goursat problem can be solved using a variant of the hodograph method based on the conservation laws (see [5] and also [8]). Moreover, this method allows us to study any interaction of the discontinuities for the original problem (3), (4), (6), (7) in the case when inequalities (8) are valid. 3.. The Goursat problem To determine the Riemann invariant R k (, t) in zone Z 5 we solve the Goursat problem for the equations (3), (4) with data that have weak discontinuities on the characteristics = ϕ(t), = θ(t) R =ϕ(t) = q, R 2 =ϕ(t) = R 2 a(z 2 (ϕ(t), t)), (8) R =θ(t) = R a(z (θ(t), t)), R 2 =θ(t) = q 2. (9) Note that data (8), (9) are matched at the point (X int, T int ) R t=tint = q, R 2 t=tint = q 2. (2) To avoid confusions we note that the problem (3), (4), (8) (2) is the initialboundary value problem for hyperbolic equations. Analysis of the solution methods of such problems one can find, for eample, in []. This problem has some originality in view of the fact that the boundary conditions (8), (9) are specified on characteristics, but the initial condition (2) is specified only in one point (X int, T int ). 4. Hodograph method We recall that the hodograph method allows us to transform a quasi-linear system of equation in a linear system, using the replacement of dependent and independent variables: (R, R 2 ) (, t). To determine t(r, R 2 ), (R, R 2 ) we have the linear systems of equations R 2 λ t R 2 =, R λ 2 t R =. (2) The compatibility condition of these equations leads to a linear hyperbolic system with variable coefficients. Taking into account (4), for the function t(r, R 2 ) we have t R R R 2 R (t R t R2) =. (22) The equation (22) has the Riemann-Green function (see, e.g. [3]) V (r, r 2 R, R 2 ) = ((R + R 2 )(r + r 2 ) 2(R R 2 + r r 2 ))(r r 2 ) (R R 2 ) 3,

9 Interactions between discontinuities and hodograph method 9 which satisfies to equation (22) for the variables R, R 2 and satisfies to conjugate equation for the variables r, r 2. In addition the condition V (r, r 2 r, r 2 ) = is valid. The eplicit form of the Riemann-Green function allows us to obtain the solution of the Goursat problem for equation (22) as (see, e.g. [4, formula (2.2)]) where t(r, R 2 ) = V (R, R 2, R, R 2 )t(r, R 2 ) + (23) + 2(R R )(R 2 R ) (R R 2 ) 3 H 2 (R 2 ) 2(R R 2 )(R 2 R 2 ) (R R 2 ) 3 H (R ) + (R R 2 ) 2 (R R 2 ) 2 t(r, R 2 ) + (R2 R ) 2 (R R 2 ) 2 t(r, R 2 ), H (R ) = R t(τ, R 2 ) dτ, H 2 (R 2 ) = R2 t(r, τ) dτ. (24) R Here t(r, R), 2 t(r, R 2 ) are the function t(r, R 2 ) on characteristics R 2 = R 2 and R = R, that one can get using conditions (8), (9). The condition (8) means that for characteristic = ϕ(t) the Riemann invariant is R = q. Using the relations (8), (), (6) and ecluding function ϕ(t) we get t = (q2 q ) 2 (R 2 q ) 2 T int = t(r, R 2 ), R = q. (25) Similarly, the condition (9) leads to (it is enough to replace indices 2) t = (q q 2 ) 2 (R q 2 ) 2 T int = t(r, R 2 ), R 2 = q 2. (26) Finaly, substituting (25), (26) in (23), (24) we obtain t(r, R 2 ) = T int V (q, q 2 R, R 2 ) = (27) = (2 )(2R R 2 + 2q q 2 (q + q 2 )(R + R 2 )) q q 2 (R R 2 ) 3. Note that the initial condition (2) is automatically satisfied because the Riemann-Green function has properties: t(q, q 2 ) = T int V (q, q 2 q, q 2 ) = T int. To obtain the function (R, R 2 ) it is enough to pay attention that after the substitutions R i = /K i (see [4]) the equation for (R, R 2 ) K K K 2 K ( K K 2) = formally coincides with the equation (22). Performing the replacements t, R i /K i in the formula (27) leads to (R, R 2 ) = ( 2 )(R R 2 ) 2 R2 + R 2(q + q 2 ) q q 2 (R R 2 ) 3 + (28) + ( (R ) 3 2 (R 2 ) 3 ) + 3R R 2 (R 2 2 R ) (R R 2 ) 3. Thus, the relations (27), (28) give us the implicit solution of the Goursat problem (3), (4), (8) (2). R 2

10 Interactions between discontinuities and hodograph method 5. Construction of the solution on isochrone The functions R (, t), R 2 (, t) that are inverse to the functions t(r, R 2 ), (R, R 2 ) allow us to obtain the Riemann invariants for the zone Z 5. Relations describe the solution to R (, t), R 2 (, t) on the interval T int t min(t 3, T 6 ) only. This should be kept in mind because, from the moment T int, the zone Z 5 disappears as a result of weak discontinuities interaction. We also note that the direct solution of algebraic equations t = t(r, R 2 ), = (R, R 2 ), in practice, is much more difficult than the solution of differential equations which are presented below. To obtain the differential equations for functions R (, t), R 2 (, t) we differentiate functions t = t(r, R 2 ), = (R, R 2 ) with respect to = t R R + t R 2R 2, = R R + R 2R 2. Obviously, these equations are valid on any level line t = t(r, R 2 ) (isochrone), where t is fied. For isochrone this system is the ODEs which can be written in the following form where dr (, t ) d = t R 2(R, R 2 ) (R, R 2 ), dr 2 (, t ) d = t R (R, R 2 ) (R, R 2 ), (29) = t R R 2 t R 2 R = (λ λ 2 )t R t R 2. (3) Here, to eclude R, R, we use (2). For the equations (29), (3) the initial conditions result from the relations (8) R (, t ) = q, R 2 (, t ) = R 2 a(, t ), = ϕ(t ). (3) Integrating the Cauchy problem (29) (3) from ϕ(t ) to θ(t ) we get the solution R (, t ), R 2 (, t ), ϕ(t ) θ(t ) for any t. We note that the derivatives t R, t R 2 in the right part of equations (29), (3) can be calculated in the eplicit form with the help (27) only. The formula (28) for function (R, R 2 ) is not required! For parameters (4) the numerical results for the Cauchy problem (29) (3) are presented in figure 3 u k Z 3 Z Z 2 Z 5 Z 6 Z 7 Z 8 u k Z Z 2 Z 5 Z 6 Z 7 Z l s ϕ θ r 2 s s ϕ = 2 l θ r 2 s Figure 3. Concentrations u k at the moments t =.8, t = T 3 = = /45.

11 Interactions between discontinuities and hodograph method 6. Interaction of discontinuities at t T 3 (and at t T 6 ) Complete scenario of the solution behavior is shown on the figure 4. To investigate the behavior of the zone boundaries after the interaction it should be remembered that the weak discontinuities of the solution do not disappear over time and move along characteristics. It means that the boundaries = w(t) and = ϕ(t) of the new zone Z 9 are the characteristics that have the weak discontinuity of the Riemann invariants R and R 2 correspondingly. Similarly, the boundaries = θ(t) and = 2 w(t) of the new zone Z are the characteristics that have the weak discontinuity of the Riemann invariants R and R 2. Note that the line of weak discontinuity for the Riemann invariant R is denoted as = w(t), and the line of weak discontinuity for the Riemann invariant R 2 is denoted as = 2 w(t). The lines of discontinuities for the Riemann invariants R and R 2 are still denoted as = θ(t) and = ϕ(t). Obviously, for zones Z 9, Z the functions θ(t), ϕ(t) are the new functions that do not coincide with previous. We focus to the solving of the problem only in zone Z 9. For zone Z the solution is obtained in a similar way. Zone Z 9 arises at time t = T 3, when the line = 2 l (t) (R2 weak discontinuity) intersects with line = ϕ(t) (also R 2 weak discontinuity). In the moment of interaction the weak discontinuity of the Reamann invariant R is appeared on line = w(t) and the weak discontinuity of the Reamann invariant R 2 is saved on line = ϕ(t). t T fin T T 9 T 6 T 3 T int Z s 2 r 2 l Z 4 Z Z 5 θ 2 Z 7 Z 8 ϕ Z Z 6 3 l r w Z 9 2 X int X 9 X 3 X 6 X fin X Φ ϕ θ 2 s f Z 2 f 2 w Z Θ Figure 4. Evolution of initial discontinuities Note that the invariant R 2 in the zone Z 9 has the same value as in zone Z 2, because line = w(t) is discontinuity line only for invariant R. Hence, for zone Z 9 we have R 2 = µ 2 (similarly for zone Z we have R = µ ). Thus, in the zone Z 9 varies the invariant R only and one can consider the equation R t + λ (R, µ 2 )R =

12 Interactions between discontinuities and hodograph method 2 only for invariant R instead of the system. Easily to formulate the continuity conditions for invariant R (and R 2 ) on the lines = w(t), = ϕ(t) R = w (t) = q, R 2 = w (t) = µ2, (32) R =ϕ(t) = R 5 =ϕ(t), µ 2 = R 2 5 =ϕ(t), (33) where R 5(, t), R 2 5(, t) are the invariants for zone Z 5 that are obtained in section 5. The conditions (32) mean that on the line = w(t) we have λ (R, R 2 ) = q q µ 2. It allows us to get the function w(t) (characteristic), i.e the motion law for left boundary of zone Z 9, w(t) = X 3 + q q µ 2 (t T 3 ). The right boundary of zone Z 9, i.e. function ϕ, can be represented in parametric form with the help of the eplicit relations (27), (28) for t(r, R 2 ) and (R, R 2 ) and the continuity conditions (33) for the Riemann invariant R 2 ϕ(t) = (ρ, µ 2 ), t = t(ρ, µ 2 ), q ρ µ, (34) where ρ is the parameter (function R on the line = ϕ(t), i.e. R 5(ϕ(t), t), plays the role of parameter). Similarly, parametric representation of the function θ(t) is θ(t) = (µ, ρ 2 ), t = t(µ, ρ 2 ), µ 2 ρ 2 q 2, (35) where ρ 2 is the parameter. To obtain the Riemann invariant R (, t) for zone Z 9 we have the Cauchy problem (classical method of characteristic) dr dt =, d dt = λ (R, µ 2 ), R t=τ = R 5(ϕ(τ), τ), t=τ = ϕ(τ). Obviously, that the solution of problem (36) is d dt = t + λ (R, µ 2 ), (36) R (, t) = R 5(ϕ(τ), τ), τ = τ(, t), (37) where τ(, t) is implicitly determined by the relation = ϕ(τ) + R 5(ϕ(τ), τ)r 5(ϕ(τ), τ)µ 2 (t τ). (38) Thus, the zone Z 9 is completely determined by formula Z 9 = {(R 5(ϕ(τ), τ), µ 2 ), ( w(t), ϕ(t))}, where for any point (, t) parameter τ is given by (38). At first sight the relations obtained are complicated. In actuality, using again the idea of isochrone, it is easy to construct a solution which depends on only one parameter ρ (see figure 5). We fi a time t and consider the isochrone t = t(r, R 2 ). For zone Z 9 we have t = t(r, µ 2 ). We introduce notation for the Riemann invariants on the isochrone:

13 Interactions between discontinuities and hodograph method 3 R (, t ) = ρ. Obviously (see figure 5), that the parameter ρ satisfies to inequality q ρ ρ, where ρ is root of equation t = t(ρ, µ 2 ), q ρ µ. (39) This root ρ is the value of Riemann invariant which corresponds to the intersection of isochrone t = t(r, µ 2 ) and line = ϕ(t) at moment t. The root ρ always eists, because the equation (39) coincides with one equation of system (34), that determine the function ϕ(t). It is easy to check that the equation (39) is a cubic equation. The solution of this equation should be chosen from the interval [q, µ ]. = w w(t) R (, t ) = ρ = ϕ(t) t = t(ρ, µ 2 ) R ( w (t ), t ) = q Z 9 R (ϕ(t ), t ) = ρ t (X 3, T 3 ) Figure 5. Fragment of the zone Z 9 After replacement R 5(ϕ(τ), τ) ρ in the relations (37), (38) we get the oneparametric representation of the Riemann invariant R for any point (, t ) where R (, t ) = ρ, q ρ ρ µ, (4) = ϕ(τ) + ρ ρ µ 2 (t τ), τ = t(ρ, µ 2 ), ϕ(τ) = (ρ, µ 2 ). (4) It is easy to check that the first condition (32) is satisfied at ρ = q. Similarly, one can obtain solution for zone Z. To write this solution it is enough to replace the indices 2 and the functions ϕ θ. For parameters (4) the numerical results are presented in figures 6, 7. u k Z 2 Z Z 9 Z 5 Z 6 Z 7 Z 8 u k Z 2 Z Z 9 Z 5 Z 7 Z ϕ w θ s r 2 s ϕ w s θ = r 2 s Figure 6. Concentrations u k at the moments t =.28, t = T 6 =.32

14 Interactions between discontinuities and hodograph method 4 u k Z 2 Z Z 5 Z 7 Z 8 Z Z 9 u k Z Z 9 Z 5 Z Z 7 Z 8 w s ϕ θ 2 w 2 s w = s ϕ 5 5 θ 2 w 2 s 2 Figure 7. Concentrations u k at the moments t =.36, t = T 9 = = 2/45 7. Interaction between strong and weak discontinuities It should be clarified that the solution presented in section 6 is valid before time t = T 9. In moment t = T 9 the left boundary of the zone Z 9 (line = w(t)) intersects the right boundary of the zone Z 2 (line = s(t)). It corresponds to the intersection of strong and weak discontinuities for the Riemann invariant R (see figure 4). Solving the equation w(t) = s(t) one can obtain the interaction point (X 9, T 9 ) T 9 = µ2 q µ q T 3, X 9 = + q µ µ 2 T 9. Starting at moment t = T 9 left boundary of zone Z 9 is modified. To get new boundary = Φ(t) we use the Rankine Hugoniot conditions (9). The Riemann invariant R 2 as before is a constant (R 2 = µ 2 ). It means that one of the Rankine-Hugoniot conditions (9) is automatically satisfied, and another condition has the form D µ 2 [ R ] = [ R ]. (42) On the new boundary = Φ(t) of the zone Z 9 (see figure 4) we have R (Φ(t), t) = µ, R (Φ(t) +, t) = R 9(Φ(t) +, t), where R 9(, t) is the Riemann invariant R for the zone Z 9 that is obtained with the help of the relations (37), (38). We introduce the notation (compare with (4)) ρ (t) = R (Φ(t) +, t) = R 9(Φ(t) +, t). Using (34), (38), and (42) we get dφ(β) dβ = D = µ µ 2 ρ (β), (43) Φ(β) = ϕ(τ) + µ 2 ρ (β)ρ (β)(β τ), (44) τ = t(ρ (β), µ 2 ), ϕ(τ) = (ρ (β), µ 2 ). (45)

15 Interactions between discontinuities and hodograph method 5 Here β is the current time t. We change the notation in order not to confuse the current time t and function t(ρ, µ 2 ). We note that (44), (45) differ from (4) only in that the parameter ρ depends from time. We also note that (43) is not the equation for the function Φ(t) which is determined by (44), (45). It is ODE for the function ρ (t) which is written in the implicit form. Adding the obvious initial conditions ρ (T 9 ) = q, Φ(T 9 ) = X 9, T 9 = µ2 q µ q T 3, X 9 = + q µ µ 2 T 9, (46) we obtain the Cauchy problem (43) (46) for determination of the function ρ (t). Similar way one can study the interaction of discontinuities which occurs at the point (X, T ) T = µ q 2 µ 2 q 2 T 6, X = 2 + q 2 µ µ 2 T, when the lines = 2 w(t) and = 2 s(t) are intersected (see figure 4). To obtain the solution it is enough to replace the indices 2 and the functions Φ Θ. In figures 8, 9 the numerical results are presented for parameters (4). u k Z Z 9 Z 5 Z Z 7 Z 8 u k Z Z 9 Z 5 Z Z Φ ϕ θ 2 w 2 s Φ ϕ θ 2 w = 2 s Figure 8. Concentrations u k at the moments t =.5, t = T =.8 u k Z 5 Z Z 9 Z Z 8 Z Z 9 Z Z 8 u k Φ 2 ϕ θ Θ Φ ϕ = θ 4 5 Θ Figure 9. Concentrations u k at the moments t =., t = T fin = = 2/5 8. The behavior of solutions at t T fin The last interaction of the discontinuities is occurred in the point (X fin, T fin ) when the lines = ϕ(t) and = θ(t) are intersected (see figure 4). Functions ϕ(t), θ(t) are

16 Interactions between discontinuities and hodograph method 6 determined by the relations (34), (35). The new zone Z appears (see figure 4): Z = {(µ, µ 2 ), ( f(t), 2 f(t))}, where the boundaries of the zone are determined by f = X fin + µ µ µ 2 (t T fin ), 2 f = X fin + µ 2 µ µ 2 (t T fin ). The boundaries of the zones Z 9 and Z are also changed Z 9 = {(R (, t), µ 2 ), (Φ(t), f(t))}, Z = {(µ, R 2 (, t)), ( 2 f(t), Θ(t))}, where functions Φ(t), R (, t), Θ(t) and R 2 (, t) can be obtained by the method which is described in section 7. The values T fin, X fin are determined by the formulae (34), (35), where ρ = µ, ρ 2 = µ 2 T fin = t(µ, µ 2 ), X fin = ϕ(t fin) = θ(t fin) = (µ, µ 2 ). We note that component concentrations in the zone Z are zero (see formulae (5), where R = µ and R 2 = µ 2 ). Thus the miture is separated into individual components. In other words, zone Z is an empty zone between the zones Z 9 and Z. It is interesting that if we a priori assume the separation of the miture then T fin can be obtained with the help of the formula (27) without the results of sections 5 7 T fin = T int V (q, q 2 µ, µ 2 ). (47) For parameters (4) the numerical results are presented in figure. u k Z Z 9 Z Z Z Φ f 2 f Θ Figure. Concentrations u k at the moment t = Conclusion The eistence of an analytic eplicit formula for the Riemann function plays a key role for the method proposed. The application of this method is not restricted to the problem of zone electrophoresis. Various equations of mathematical physics which can be studied by this method are listed for eample in [5]. Such equations are the shallow water equations, gas dynamics equations, the Born-Infeld equation, nonlinear hyperbolic heat equation, Euler Poisson Darbou equation, chromatography equations etc. Note that there is an alternative way: this is application of the generalized hodograph [5].

17 Interactions between discontinuities and hodograph method 7 If the initial data of the Cauchy problem allow to obtain an eplicit formulae for the commuting flows, then these flows can be used instead of the RiemannGreen function for construction of the solution implicit form, i. e. t(r, R 2 ), (R, R 2 ). Another important result is the method proposed for recovery of an eplicit solution with the help of its implicit form. In fact, it is the transformation of the Cauchy problem for PDEs into the Cauchy problem for ODEs on the level lines (isochrones). This method allows not only to solve the Goursat problem, but also to solve effectively the Cauchy problem with arbitrary initial data (including discontinuous data, see [8]). Finally, to emphasize the importance of eact solutions, we present in figure the comparison of the eact solution with numerical solution..5 u k..5 a b a b Figure. Calculations (dotted lines) and theory, a) t =.67, b) t =.28 The solid line corresponds to eact solution of the Goursat problem obtained in sections 5, 6, and the dotted line corresponds to the solution constructed by a numerical method (FEM and VOF). Figure illustrates the smoothing of the numerical solution. It indicates to uncertain identification of weak discontinuities and, as a result, the impossibility of the discontinuities interaction checking. Acknowledgments The research was supported by the Base part of the project 23.-/24- (subject 26 year) the Ministry of Education and Science of the Russian Federation, Southern Federal University. The authors are grateful to N.M. Polyakova for proofreading the manuscript. Appendi A. Numerical-analytical method for the Cauchy problem We give a brief description of the numerical-analytical method, which allows to solve the Cauchy problem (3), (6) with any initial data (in detail, see [8]). The implicit solution

18 Interactions between discontinuities and hodograph method 8 of the problem can be written in the following form (compare with (23)) where t = t(a, b) 2(b a) (r r 2 ) (r + r 2 )F (a, b) + 2r r 2 G(a, b) 3 (r r 2 ) 3 (r r 2 ), (A.) 3 r = R (b), r 2 = R 2 (a), F (a, b) = b f(τ) dτ, G(a, b) = b g(τ) dτ, a a f(τ) = R (τ) + R 2 (τ) R (τ)r 2 (τ), g(τ) = R (τ)r 2 (τ). On the (a, b)-plane we choose the point (a, b ) which identifies the isochrone t = t(a, b ). In practice point (a, b ) can be selected with the help of the level lines of the function t(a, b) for some interval of the parameters a, b. We assume that the isochrone is parameterized by functions a = a(µ), b = b(µ), where µ is the parameter. To obtain the functions a = a(µ), b = b(µ) on isochrone we solve the Cauchy problem da dµ = t b, db dµ = t a, df dµ = F µ, dg dµ = G µ, dx dµ = X µ, (A.2) a µ= = a, b µ= = b, F µ= = F, G µ= = G, X µ= = X, where the right parts of differential equations can be calculated in the eplicit form F µ = f(b)t a + f(a)t b, G µ = g(b)t a + g(a)t b, X µ = (λ 2 (r, r 2 ) λ (r, r 2 ))t a t b, Integrating from b = a to b = b the Cauchy problem dy (b) db X (a, b). = (λ 2 (r (b), r 2 (a ))t b (a, b), Y (a ) = a, (A.3) df (a, b) dg(a, b) = f(b), = g(b), F (a, a ) =, G(a, a ) =, db db we obtain values X = Y (b ), F = F (a, b ), G = G(a, b ). For any value of the parameter µ the solution of the original problem (3), (4) with initial data (6) has the form R (, t ) = R (b(µ)), R 2 (, t ) = R 2 (a(µ)), = X(µ). The method described can be used for the solution of any Cauchy problem, if the eplicit relation for the function t(a, b) is given. Strictly speaking, numerical part of the method is required only to solve the Cauchy problems (A.2) and (A.3). In some cases, these problems can be solve eplicitly, as, for instance, for the Born Infeld equation (see e.g. [6]). In a sense, the method is a generalization of the method of characteristics for the case of two equations. In this connection see paper [7], which develops quite similar method for solving equations of chromatography.

19 Interactions between discontinuities and hodograph method 9 References [] Rozdestvenskii B L, Janenko N N 983 Systems of Quasilinear Equations and Their Applications to Gas Dynamics (Translation of Mathematical Monographs vol 55) (Provience, RI: American Mathematical Society) [2] Elaeva M S 2 Interaction between strong and weak discontinuities in the Riemann problem for hyperbolic equations University news. North-Caucasian Region. Natural Sciences series [3] Copson E T 958 On the Riemann-Green Function Arch. Ration. Mech. Anal [4] Bitsadze A V 98 Some Classes of Equations in Partial Derivatives (Moscow: Nauka) p 448 [5] Senashov S I, Yakhno A 22 Conservation laws, hodograph transformation and boundary value problems of plane plasticity SIGMA 8 6 [6] Ferapontov E V, Tsarev S P 99 Systems of hydrodynamic type that arise in gas chromatography. Riemann invariants and eact solutions Mat. Model [7] Babskii V G, Zhukov M Yu, Yudovich V I 989 Mathematical Theory of Electrophoresis (New York: Plenum Publishers Corporation) p 24 [8] Shiryaeva E V, Zhukov M Yu 24 Hodograph Method and Numerical Integration of Two Quasilinear Hyperbolic Equations. Part II. The Zone Electrophoresis Equations arxiv: [9] Elaeva M S 2 Investigation of zone elecrophoresis for two component miture Mat. Model [] Elaeva M S 2 Mathematical modeling of capillary zone electrophoresis. Dissertation of candidate of physico-mathematical sciences: (Rostov-on-Don: SFU) p 6 [] Elaeva M S 22 Separation of two component miture under action an electric field Comp. Math. and Mat. Phys [2] La P D 26 Hyperbolic Partial Differention Equations (Provience, RI: American Mathematical Society) [3] Liu T-P 2 Hyperbolic and Viscous Conservation Laws (CMBS-NSF Regional Conference Series in Applied Mathematics, SIAM) [4] Pavlov M V 987 Hamiltonian formalism of electophoresis equations. Integrable equations of hydrodynamics (Moscow: Preprint of Landau Inst. for Theor. Phys. 7) p 7 [5] Tsarev S P 99 The geometry of hamiltonian systems of hydrodynamic type. The generalized hodograph method USSR Izvestiya Mathematics 54: [6] Shiryaeva E V, Zhukov M Yu 25 Hodograph Method and Numerical Integration of Two Quasilinear Hyperbolic Equations. Part III. Two-Beam Reduction of the Dense Soliton Gas Equations arxiv: arxiv:52.67 [7] Curro C, Fusco D, Manganaro N 25 Eact description of simple wave interactions in multicomponent chromatography J. Phys. A: Math. Theor