Conservative cross diffusions and pattern formation through relaxation
|
|
- Beatrice Mason
- 5 years ago
- Views:
Transcription
1 Conservative cross iffusions an pattern formation through relaxation Mostafa Benahmane, Thomas Lepoutre, Americo Marrocco, Benoît Perthame To cite this version: Mostafa Benahmane, Thomas Lepoutre, Americo Marrocco, Benoît Perthame. Conservative cross iffusions an pattern formation through relaxation. Analyse mathématique et numérique e sytèmes cross-iffusion. 29. <inria > HAL I: inria Submitte on 17 Mar 29 HAL is a multi-isciplinary open access archive for the eposit an issemination of scientific research ocuments, whether they are publishe or not. The ocuments may come from teaching an research institutions in France or abroa, or from public or private research centers. L archive ouverte pluriisciplinaire HAL, est estinée au épôt et à la iffusion e ocuments scientifiques e niveau recherche, publiés ou non, émanant es établissements enseignement et e recherche français ou étrangers, es laboratoires publics ou privés.
2 Conservative cross iffusions an pattern formation through relaxation Mostafa Benahmane Thomas Lepoutre Americo Marrocco Benoît Perthame 13th March 29 Abstract This paper is aime at stuying the formation of patches in a cross-iffusion system without reaction terms when the iffusion matrix can be negative but with positive self-iffusion. We prove existence results for small ata an global a priori bouns in space-time Lebesgue spaces for a large class of iffusion matrices. This result inicates that blow-up shoul occur on the graient. One can tackle this issue using a relaxation system with global solutions an prove uniform a priori estimates. Our proofs are base on a uality argument à la M. Pierre which we exten to treat egeneracy an growth of the iffusion matrix. We also analyze the linearize instability of the relaxation system an a Turing type mechanism can occur. This gives the range of parameters an ata for which instability may occur. Numerical simulations show that patterns arise inee in this range an the solutions ten to exhibit patches with stiff graients on boune solutions, in accorance with the theory. 1 Introuction The ynamics of interacting population with cross-iffusion have been wiely investigate by several researchers. The concept of this phenomena was stuie by Levin [18], Levin an Segel, [17], Okubo [27], Mimura an Murray [24], Mimura an Kawasaki [23], Mimura an Yamaguti [25], an many other authors. All these papers base the pattern formation on a reaction term as prey-preator interactions. Spatial patterns can however emerge from pure iffusions without reaction terms nor oriente rift at the iniviual level. This is the case of N populations escribe microscopically by a brownian process which intensity epens upon the macroscopic ensity U = (U 1,..., U N ) of the populations X k (t) = σ k ( U(X, t) ) Wk (t), 1 k N. Departamento e Ingeniería Matemática Faculta e Ciencias Físicas y Matemáticas, Universia e Concepción, Casilla 16-C, Concepción, Chile mostafa benahmane@yahoo.fr INRIA-Paris-Rocquencourt, Team BANG, Domaine e Voluceau, BP15, LeChesnay Ceex, France. s: thomas.lepoutre@inria.fr Université Pierre et Marie Curie-Paris 6, UMR 7598 LJLL, BC187, 4, place Jussieu, F Paris ceex 5 americo.marrocco@inria.fr Institut Universitaire e France. benoit.perthame@upmc.fr 1
3 When set on a boune omain with reflexion on the bounary, the corresponing moels for the population ensity are cross-iffusions { tu A(U) =, in, (1) n A(U) = on, where U = U(x, t) R N, x a smooth boune omain of R, n enotes the outwar normal to. Finally A : R N R N is a nonlinearity relate to the intensity of the interactions by the relation A k (U) = U k a k (U), (2) an a k (U) = 1 2 σ k(u).σk t (U). We also complete the system with an initial ata U(t = ) = U = (U 1,..., U N ) with U k. The properties, an pattern formation capacity of such systems are better escribe by introucing the more general form N t U k iv[d kl (U) U l ] =, (3) l=1 where D kl (U) are the components of a N N matrix, the erivative of A in the case (1). Bounary conitions have to be impose an we consier here the case of Neumann conitions D(U). U.n =, on. For such bounary conitions, mass conservation yiels naturally U(t) = U, t, where U enotes the average U = 1 U(x)x. The Lotka-Volterra competition with cross-iffusion has recently receive great attention. They are many establishe results concerning the global existence of classical solutions (see [33, 19] an the references therein) where most of the proofs rely on Amann s theorem [1, 2]. We point out that stanar parabolic theory is not irectly applicable to our moel ue to the presence of cross-iffusion terms. In opposition with pattern formations, an important issue has been wiely stuie which is to know in which circumstances the solutions exist globally an behave like in the case a single heat equation, i.e., relax to a constant state as t. Typically three kins of special methos have been helpful in this scope. The first metho is to rely on the maximum principle. It can occur on certain combinations of the U i as in [1, 2, 3]. Entropy methos also applies to particular systems an has also been a useful tool because of the relate symmetrization of the system following [15, 8]. It provies a natural metho both for existence an relaxation to steay state as in the recent stuies in [6, 7] of the Shigezaa-Kawasaki prey-preator system [31], or for tumous moels [12]. This metho typically applies in the special case of the square entropy when D is efinite positive meaning that there is ν > such that N ξ k D kl (U)ξ l ν ξ 2, ξ R N. (4) k,l=1 2
4 This strong positivity property gives the energy inequality 1 2 t U(t) 2 x = N k,l=1 U k D kl (U) U l ν U(t) 2 x. The thir metho, by uality, has been use in [5] on a particular upper-iagonal iffusion with Dirichlet bounary conition; we show here that the metho can be extene to very general systems with Neuman conitions. Concerning instabilities, the interplay between iffusion an reaction terms has raise surprising results in the spirit of the Turing instability mechanism [32]. The question to know if cross-iffusion or self-iffusion gives an avantage to competing species is stuie in [2, 21]. Our interest in this paper concerns instability mechanisms that may appear only from the iffusion intensity an the lost of positivity in the secon orer matrix. We stuy in which circumstances the increase of this intensity with higher ensity of the other species can lea to a segregation phenomena. Of course such instability is incompatible with any entropy inequality, an thus (4) cannot hol. This rises several mathematical questions which seem to be new in the omain of cross-iffusions. Can it still be that small solutions exist globally even though the maximum principle oes not hol in general? For large ata, what kin of regularity or blow-up can we expect? Finally, how o regularize systems behave in the instability regime. We will stuy these questions with a moel problem in min relate to Shigezaa-Kawasaki s system an that represents two species with stronger interactions { t U 1 [U 1 (1 + a 11 U p 1 + a 12U p 2 )] =, x, (5) t U 2 [U 2 (1 + a 21 U p 1 + a 22U p 2 )] =, still with Neumann bounary conitions. One can check that as soon as p > 1, the matrix D kl is negative for U large. We approach these questions both theoretically an numerically. In particular we prove existence for small initial ata (section 2) an we give a priori bouns in L p t,x for possible solutions to (5) thus showing that the break-own shoul come from the blow-up of graient estimates rather than usual L p norms (section 3). Our main tool here is a general estimate ue to M. Pierre [29, 28] in the context of semilinear parabolic systems (arising in population ynamics or more generally reaction-iffusion systems) (see also [9]). In section 4, these bouns are extene to a relaxation system that takes into account a local measurement of ensities; we show that the metho is well aapte to general (even not parabolic) cross-iffusions an prove global existence for the relaxation system. For non-parabolic cases, we show that Turing instability occurs in a certain range of ata an for small relaxation parameters. Numerical simulations of this relaxation system are performe in section 5. They show that the oscillatory initial regime reorganizes to create patches where one species ensity ominates the other an interfaces are generate which with is relate to the relaxation length. The technical an general extension of M. Pierre s estimate to boune omains for Neumann bounary conitions is kept for an Appenix as well as another remarkable energy estimate which hols for particular cross-iffusion coefficients in (5). 2 Global solutions for small ata The lack of maximum principle for general iffusion systems is a major ifficulty that arises for systems as (3). Using stronger H 1 estimates, we can show in one imension that for small initial ata there 3
5 is global existence. Such solutions ecay to the constant state for large time an this is incompatible with the patterns formation we are intereste in. This inicates that large initial ata are necessary for pattern formation as expecte in general. We consier the system uner the form (1) (2) an assume that a k C 1 (R N ), (6) a k () ν >. (7) Theorem 2.1 (Global small solutions in 1 imension) In one imension, uner assumptions (6),(7), an for an initial ata satisfying, with α small enough, U L () + U L 2 () α, there is a global solution U(t, x) to the cross-iffusion system (1) (2). It satisfies for all t > with C inepenent of t an U(t) L () Cα, U(t) L 2 () α, U(t) t U. Proof. Firstly, since U(t) is a priori conserve, we notice that U(t) L () U L () + U L 2 () α + U L 2 (). Seconly, we multiply (1) by the vector A (U) an ifferentiate. We obtain t A(U) = ( A (U) t U ) = ( A (U) A(U) ), we multiply by A(U) an integrate A(U) 2 = A(U) ( A (U) A(U) ) = t 2 A(U) A (U) A(U). (8) But from (2) an (7), we have A () kl = a kk ()δ kl, an thus for U L () ε small enough, we have XA (U)X ν 2 X 2, X R N, (9) We now consier T efine by A (U)X 2 ν2 2 X 2, X R N. (1) T := sup{t, U(t) ε}. For α < ε, then T >. Suppose T <, then for t T, we have, from (8) an (9), A(U(t)) 2 x A(U ) 2 x. Therefore, A ( U(t) ) U(t) 2 4 A(U ) 2 x.
6 Using (1), this leas to U(t) 2 x 2 ν 2 A(U ) 2 x 2C ν 2 U 2 x. Now, we choose α such that α ε 2C 3 an α ε ν 3, this ensures, U(t) L () 2ε 3 an thus, T is not maximal. Therefore T =. Finally, the existence of a solution for small times follows from stanar parabolic theory an the a priori boun above shows that these are global solutions. We now prove the time convergence to U. Because (U U ) 2 = U 2 U 2, we compute U ) t (U 2 = U 2 = A(U)U = A (U) U. U. t But as the solution stays in the omain U(t) L () ε, we have for some constant C, U t (U ) 2 ν U 2 2C (U U ) 2, 2 thanks to the Poincaré Wirtinger inequality. We conclue from Gronwall lemma that U U L 2 () e C t U U L 2 (). 3 A priori bouns for large ata For large initial ata an when the conition (4) oes not hol, we cannot expect in general the existence of solutions for the cross-iffusion system (3) (2). For a single equation, the corresponing situation is when A (u) can be negative on some interval I R +. u A(u) =. t The situation is analyze in [3, 26, 11] (see also the survey in [1]) an it is better analyze in term of relaxation systems, an approach which we will follow later. We expect that oscillations or jumps occur at positive times, but a first issue is a priori control in L for possible solutions. This follows from the maximum principle for a single equation (an possibly from entropy constructions for relaxation systems, see [3]). For systems this is an open question an we give here a first a priori boun Theorem 3.1 Smooth solutions to (3) (2) satisfy the a priori bouns ( T N N ) 1 2 A k (U) U k xt C1 () U ( N L 2 () + C 2, Uk ) T. (11) k=1 Therefore if we assume in (2) that k=1 then we also have, with the notation Q T = [, T ], a k (U) ν > k = 1,..., N, (12) ν U L 2 (Q T ) C 1 () U L 2 () + C 2 (, U ) T. (13) k=1 5
7 In the particular case of moel (5), we observe that the larger is p, the best is the boun in (11). In particular, A(U) is always integrable. The L 2 estimate in (13) is much weaker. Proof. Our proof is base on a variant of a general uality argument ue to [28, 9], that is presente in Appenix A. We enote w = N U k. We sum up the equations an we fin which we write k=1 t w N a k (U)U k =, k=1 t w α(t, x)w =, α(t, x) := We can use now (27) in Appenix A an obtain, N k=1 a k(u(t, x))u k (t, x) w(t, x) = α(u(t, x)). α w L 2 (Q T ) C() w L 2 () + 2 w α L 2 (Q T ). (14) In orer to control the right han sie by α w L 2 (Q T ), we use a truncation metho. Since the coefficients a k (U) are continuous, an U k are nonnegative, we may efine for any R >, sup α(u) := M(R) < +. w R Furthermore, we may truncate w away from values less than R, a parameter to be fixe later on, with the inicator function 1I {w R} an rewrite (14) as α w 1I {w R} L 2 (Q T ) C() w L 2 () + 2 w α1i {w R} L 2 (Q T ) + 2 w α 1I {w R} L 2 (Q T ), α w 1I {w R} L 2 (Q T ) C() w L 2 () + 2 w α w R 1I {w R} L 2 (Q T ) + 2 w M(R) L 2 (Q T ). We choose R = 4 w an obtain α w 1I {w R} L 2 (Q T ) 2C() w L 2 () + 4 w T M(R). Since we also know that αw1i {w R} L 2 (Q T ) R M(R) T = 4 w M(R) T, we conclue α w L 2 (Q T ) 2C() w L 2 () + C 2 ( w ) T, with C 2 ( w ) = 8 w M(4 w ). This is exactly the a priori estimate (11). The other statement is a simple an irect consequence. 6
8 4 A relaxation system If we assume that the intensity of the brownian motion epens on the ensity of the populations measure with a space scale δ > an not at the exact location x, then the system (5) can be replace by a cross-iffusion relaxation system { t u k [a k (ũ)u k ] =, x, k = 1,..., N, δ 2 ũ k + ũ k = u k, (15) together with Neumann bounary conitions both on u k an ũ k. Relaxation proceures are usual an several other examples for cross-iffusions can be foun in [14, 4] an for phase transitions see [13, 1, 3]. In terms of the ecological interpretation, it is also more realistic than the initial system (5), because iniviuals are unlikely to be able to access a pointwise ensity, but might estimate their environment from sensing at a smaller scale. We can expect that the system (15) is well-pose, an we first stuy this question. Then we prove uniform bouns inepenent of δ which inicates that instability shoul arise from the blow-up of graients. To tackle the question of instabilities, we show that the system exhibits Turing patterns for δ small, this is our secon goal in this section. We keep in min the example (5) an assume that for some p > one has < ν a k (U) C (1 + U p ), k {1,..., N}. (16) For later purpose, we also introuce the assumption that for some constant K > an some η >, we have a k (U) a η k (U) K U. (17) We have in min coefficients of the form (1 + Ũ p p 1 j ) an then we can take η = p. 4.1 Uniform estimates for p < 2 We first exten the a priori estimate of section 3 to this relaxation system. The coupling inuces a limitation on the possible growth of the nonlinearities a k (U) an we have the Theorem 4.1 Assume that (16) hols for some < p < 2, then, the a priori boun hols for a constant C inepenent of δ u L 2 (Q T ) C( u L 2 (), T ), T >. This is weaker than the a priori boun in (11). The ifficulty in the case at han comes from the epenency of a k (ũ) which we cannot lower boun from u itself. Proof. We enote by ã k the quantity a k (ũ). The estimate (27) of Appenix A gives, for all k {1,..., N} ã k u k L 2 (Q T ) C() u k L 2 () + 2 u k ã k L 2 (Q T ). 7
9 The last term may be estimate as T T ã k 2 L 2 (Q T ) = ã k C (1 + l ũ p l ). Thanks to Holer inequality an irect estimate on the solution to the elliptic equation on ũ l, we have, T ũ p l ( T ) 2 p 2 ũl p 2 p 2 L 2 (Q T C()T 2 ul p ) L 2 (Q T ). Finally, back to the original inequality we arrive at ν uk L 2 (Q T ) C() u k L 2 () + 2 C T u k + 2 u k C ( T ) 2 p p 4 u 2 L 2 (Q T ), which leas to ν u L 2 (Q T ) C(, u L 2 (), T ) + C(, u )T 2 p 2 u p 2 L 2 (Q T ). As p/2 < 1, this proves that u L 2 (Q T ) is a priori boune as by a constant epening only on, T, u L 2 () an the two constants in (16). 4.2 Existence of solutions We now show stronger estimates from which strong compactness of solutions follows. They use funamentally the regularity on ũ in (15) by elliptic regularizing effects. Existence of global solutions follow an the etails are carrie out in [16]. The main result of this section is the following Proposition 4.2 Assume that (16) hols with p > 1, an p < when > 2. Then, the a priori estimate hols ã i u i L 2 (Q T ) C(δ, u L 1 L 2 (), T ). (18) Furthermore, if we assume (17) in imension 1 with any η >, an in imension 2 with < η < 1, then we have for all 1 q <, 2 2 u(t) L q () C(q, δ, u L 1 L q (), T ), t T, (19) T u q/2 2 xt C(q, δ, u L 1 L q (), T ). (2) Proof. We begin with the proof of (18) which improves that of the theorem 4.1. We use again the estimate (27) applie to u i which yiels ν ui L 2 (Q T ) ã i u i L 2 (Q T ) C() u i L 2 () + 2 u i ã i L 2 (Q T ). (21) We use the hypothesis (16) to get T ã i L 2 (Q T ) C T + C ũ p C () T T + 8 ũ p = C ( T T + ũ p p).
10 Thanks to elliptic regularity we also have ũ p C(δ, r) u r, for any r > 1 satisfying also 1 p 1 r 2 (particularly it is true for any r if = 1, 2). Then, using interpolation inequality an choosing r < 2, we fin successively ũ p C(δ, r) u 1 θ(r) 1 u θ(r) L 2 () = C(δ, r) u 1 θ(r) 1 u θ(r) L 2 () θ(r) = 1 1 r = 2(1 1 r ), ã i L 2 (Q T ) C ( T + T ã i L 2 (Q T ) C ( T + C(δ, r) p/2 u (1 θ(r))p/2 1 C(δ, r) p u (1 θ(r))p 1 u pθ(r) L 2 ()), T Now, if we may choose r such that θ(r)p < 2, we get, thanks to Jensen s inequality T which we rewrite as u(t) pθ(r) L 2 () t = T Replacing in (22), we have T An replacing in (21), we obtain u(t) pθ(r) L 2 () t). (22) ( T ) pθ(r)/2 u(t) 2(pθ(r)/2) t T 1 pθ(r)/2 u(t) 2 L 2 () L 2 () t, u(t) pθ(r) L 2 () t T 1 pθ(r)/2 u pθ(r) L 2 (Q T ). ã i L 2 (Q T ) C T + C(δ, r, p, u 1 )T 2 pθ(r) pθ(r) 4 2 u L 2 (Q T ) ν u L 2 (Q T ) C() u L 2 () + C () T + C(δ, r, p, u 1 )T 2 pθ(r) pθ(r) 4 u 2 L 2 (Q T ) (23) This conclues the first inequality when pθ(r)/2 < 1, an it remains to fin the range of p in orer to fulfill the constraints. These can be obtaine choosing r close enough to 1 for = 1, 2. For > 2, we nee the conitions { 1 p 1 r 2, 1 < r < 2, pθ(r) 2 = p(1 1 r ) < 1. We choose to satisfy the secon line 1 r > p 1 p it a posteriori). This leas, in the first line, to the conition p <, but close to equality (which gives 1 < r < 2 as we check 2 2, but close to equality (which imposes r < 2 +2 ). The boun on u L 2 (Q T ) gives then a boun on ã i u i L 2 (Q T ) thus concluing the proof of (18). This estimate leas to the stronger a priori bouns (19) that we prove now. We go back to the equation an write u q + C q ã u q/2 2 = Cq 1 u q/2 u q/2 ã. t 9
11 From this equality, an writing ã u q/2 = (ã 1/2 u q/2 ) u q/2 ã 1/2, we erive irectly the inequality t u q + Cq 3 ã u q/2 2 + Cq 3 (ã 1/2 u q/2 ) 2 Cq 2 u q ã 2 C q u q ũ 2 ã 2η 1. ã (24) We show separately how in imensions 1 an 2 this allows us to control any L q norm for q < +. The case = 1. The proof is easier for = 1 because there exists a constant C (epening only on δ an ) such that ũ L () + ũ L () C u C( u 1 ). Therefore from (24) we euce t u q C.C q We conclue thanks to Gronwall lemma. The case = 2. In imension 2, we ivie the proof into two steps. Step 1, 1 < q < 2. We first focus on small values of the exponent q namely 1 < q < 2 an 2η < 1 + q 2 (the limitation η < 1 comes from choosing q close to 2). From (24) an using successively (17) an Höler inequality, we obtain t u q C 3 q u q. u q ã q/2 ũ 2. Cq 3 u ã q qr ũ 2 2r thanks to elliptic regularity, with Choosing r = 2 q then 1 m = 1 q 4 t C(q, δ, r) u ã q qr u 2 m 1 m = 1 2r = 1 1 2r. an we arrive at u q C(δ, r) u ã q 2 u 2 m C(δ, r) u ã q 2 u 2θ q u 2(1 θ) L 2 () C(δ, r) u ã q 2 u q q u 2 q? L 2 () because, as we have 1 2 < 1 m < 1 q, we may interpolate m between q an 2 with θ = 1 m q 1 2 = 2 q 4 2 q 2q = q 2. We finally obtain by Young s inequality, [ u q C(q, δ, r) u ] ã 2 2 t + u 2 L 2 () u q q, 1
12 an we may then conclue with Gronwall lemma using the estimates on T u ã 2 2 an T u 2 2 in (18). By interpolation, this also gives a priori boun for any L q norm for q [1, 2[. This ens step 1. Step 2 We now focus on L q norms for q 2. We notice that now, controlling any L q norm for q < 2, we control by elliptic regularity any L q norm of ũ except the L norm. We also control the L norm of ũ an therefore of ã. We go back to (24) an conclue t u q + C u q/2 2 C u q ã 2 ã C(δ, q, T ) u q ũ 2 C(δ, q, T ) u q qr ũ 2 2r C(δ, q, r, T ) u q qr. We use now interpolation: for any s > r > 1, we have u qr u 1 θ q u θ qs θ = 1 qr 1 qs 1 q 1 qs = 1 r 1 s 1 1. s Using Young s inequality, we obtain for ε small (to be chosen later) u q + C u q/2 2 C(δ, q, r, T, ε) u q q t + ε u q qs. We have by Poincaré Wirtinger inequality u q qs = u q/2 2 2s 2 u q/2 u q/2 2 2s + 2 uq/2 2 2s C(, s, ν) uq/2 2 + C(, s) u q q/2 C(, s, ν) uq/2 2 + C(, s) [ u q 1 + u q q], from interpolation an Young s inequality. We fix s as above, choose ε small enough an we obtain u q C(, q) [ u q 1 t + q] u q. An we conclue (19) by Gronwall lemma. The last estimate (2) also follows from (24). 4.3 Turing patterns In orer to go further an stuy the instability occurring in the regularize moel, we consier the following particular system: t u ( u(1 + ṽ 2 ) ) =, t v ( v(1 + ũ 2 ) ) =, δ 2 ũ + ũ = u, δ 2 ṽ + ṽ = v, still with Neumann bounary conitions an initial ata u, v. (25) 11
13 It is rather intuitive that for δ large, iffusion is ominant; this is also the case for small initial ata thanks to the argument in section 2. Therefore, the appearance of patterns epens upon a relation between the average ensities of populations u an v an the parameter δ. In orer to stuy this in etails, we begin with some notations the only possible constant steay state of the system is given by u = ũ = u an v = ṽ = v, we enote by (λ >, w) the non-zero solutions to the Neumann eigenproblem w = λw, n w = on we also enote by λ 1 the first eigenvalue for the Laplacian. In orer to investigate when the (in)stability of the constant steay state occurs, we stuy the linearize system: t u (1 + v 2 ) u 2 u v ṽ =, t v (1 + u 2 ) v 2 u v ũ =, δ 2 ũ + ũ = u, δ 2 ṽ + ṽ = v. As usual, we look for solutions of type e µt (a, b, c, )w. Such solutions shoul satisfy c = a 1+δ 2 λ, = b 1+δ 2 λ, µa + λa(1 + v 2 ) + λ2 u v b 1+δ 2 λ =, µb + λb(1 + u 2 ) + λ2 u v a 1+δ 2 λ =, which may be written uner the matrix form 1 + v 2 2 u v 1+δ 2 λ 2 u v 1+δ 2 λ 1 + u 2 We enote by M = M( u, v, δ, λ) this symmetric matrix. ( ) ( ) a = µ a. b λ b The question of the stability of the constant steay state can now be formulate in terms of eigenvalues of the matrix M. It is unstable if µ > an thus if M has negative eigenvalues. In this case, local behavior aroun the equilibrium shoul lea to segregation since the associate eigenvector to µ/λ shoul satisfy a.b <. We have the following Lemma 4.3 If the initial populations u an v are large enough an the relaxation parameter δ is small enough, then the constant steay state is linearly unstable. More precisely, it occurs uner the conitions γ := 4 u 2 v 2 ( 1 + u 2)( 1 + v 2) >, (1 2 u v + u 2)( 1 + v 2) δ 2 < (1 λ 1 + u 2)( 1 + v 2). (26) 12
14 The omain influences instability only through the smallness conition on δ 2 λ 1 () when the initial ata are such that γ >. Notice that the first conition ensures that the limiting system (δ = ) has a negative iffusion matrix D in the setting (3). Proof. As mentione earlier, the constant steay state is unstable if the symmetric matrix M amits negative eigenvalues, i.e., if et(m) <. We calculate et(m) = 4 u 2 v 2 (1 + δ 2 λ) 2 + (1 + u 2 )(1 + v 2 ) γ. As et(m) is a non-ecreasing function of δ, with limit γ as δ we first nee γ > that is our first conition. The secon conition gives the upper boun on δ to satisfy this inequality. 5 Numerical results The theoretical results inicate that solutions of the relaxation system (15) remain boune in L 2. Therefore, we expect that the instability obtaine for large initial ata or a small δ (through Turing mechanism) shoul lea to stiff graients. We present several numerical tests for the particular cubic system (25). They aim at showing that (i) the conitions of Lemma 4.3 are accurate an escribe the numerical transition to instability, (ii) stationary patterns are inee obtaine in this range of ata with stiff graients. These numerical results also show the variety of possible steay state, an interesting phenomena wiely stuie theoretically ([21] an the references therein). We have performe both 1D an 2D simulations in the following omains In interval = [, 1] (1D simulation) In rectangle = [, 2] [,.5] (in both cases = 1) In unit square =], 1[ 2 In 2D, the computations use an unstructure gri an a mixe finite element metho for space an backwar Euler scheme for time. The metho is alreay presente in [22]. We recall the eigenvectors of Laplace operator with Neumann bounary conition: for = [, 1], e n (x) = cos(nπx) an particularly, the first nonzero eigenvalue is π 2, associate to the eigenvector cos(π x). For = [, 2] [,.5], the eigenvectors are given by e n,m (x, y) = cos( nπx 2 ) cos(2mπy), the first nonzero eigenvalue is π 2 /4. We compare the theoretical formula of Lemma 4.3 an the numerical stability of the steay state. In all simulations we take u = 2, v = 1. In this case, instability might occur, since γ = 16 1 = 6 > an the limiting values of δ are given in the table 1. Therefore, in a first series of numerical tests, we choose the parameters δ 2 =.26 < δ 2 an δ2 =.27 > δ 2 in 1D an for δ2 =.1 < δ 2, an δ 2 =.11 > δ 2 in 2D. In both cases, we have obtaine relaxation to constant equilibrium when δ is taken larger than the critical value (for all the initial ata we have teste), an instability of the constant equilibrium when δ is smaller than the critical value. 13
15 Domain =], 1[ =], 2[ ],.5[ =], 1[ 2 Critical value δ Table 1: Critical value of the parameter δ for Turing instability, compute from formula (26). (a) Intital conition (b) Steay state Figure 1: Initial conition (left) an steay state (right) in 2D simulations. The relaxation parameter δ 2 =.1 is small enough to fulfill conition (26). The scales for the solutions are not the same in the two figures. We illustrate the instability case with steay states in figure 3 for 1D simulations an in figure 1 for 2D simulations. For the 1D simulation, we took v 1 an u = I {].1,.6[}. Next we stuy the singularity that occurs on the transients for small relaxation parameter δ. Numerical solutions show that strong oscillations occur. In figure 2 we epict, for the same initial ata, the effect of δ on the solution at a given time. A Appenix: Michel Pierre s estimate Consier the problem { t u [a(t, x)u] =, u(t = ) = u, together with Neumann bounary conition in a boune omain. We enote Q T = (, T ) ). We assume that a(t, x) > is smooth an u is a weak solution. We can also assume withour lack of generality that u. Then we have the a priori estimate Lemma A.1 For any T >, we have a u L 2 (Q T ) C() u L 2 () + 2 u a L 2 (Q T ), (27) where C() is the constant of Poincaré Wirtinger s inequality. 14
16 (a) δ 2 = (b) δ 2 = (c) δ 2 = () δ 2 = Figure 2: Cuts, at a given time, in the y irection an in the mile of the omain = (, 1) 2 in 2D. The piecewise initial conition is also represente in ashe line. As expecte strong oscillations occur with species seggregation. These oscillations are stronger when δ is smaller. 15
17 (a) time t = (b) time t = 2 (c) time t = 5 () time t = 9 Figure 3: Time evolution for a 1D simulation for δ 2 =.25 < δ 2. This figure shows how a small perturbation is amplifie an creates a steay pattern. Because δ is large (close to δ ) there are not strong oscillations as in the case of smaller values. 16
18 Proof. Consier smooth functions F (t, x) an the solutions to the ajoint problem { t v + a(t, x) v = F (t, x), still with Neumann conitions. We have v(t = T ) =, uv = F u, t an thanks to the final conition for the ajoint problem, T u v = F u. (29) Multiplying (28) by v, we get integrating by parts on, we obtain, v 2 t 2 t v v + a v 2 = F v, + a v 2 ( F 2 2a + a 2 v 2 ), which gives after integration in time, using again v(t ) =, T v 2 + a v 2 an by consequence, T F 2 a, (28) v L 2 () F a L 2 (Q T ), (3) a v L 2 (Q T ) F a L 2 (Q T ). (31) We nee aitionally a boun on v that we erive as follows. We use again (28) to fin T T v ( a v F ) = a v F a + a, which gives, thanks to the Cauchy Schwarz inequality an (31), v 2 a L 2 (Q T ) F a L 2 (Q T ). (32) Finally, we get using Poincaré-Wirtinger ineqality, (32) an then (3), u v u (v v ) + u v C() u L 2 () v L 2 () + 2 u a L 2 (Q T ) F a L 2 (Q T ) C() u L 2 () F a L 2 (Q T ) + 2 u a L 2 (Q T ) F a L 2 (Q T ). 17
19 Back to (29), we conclue that T T ( F u = F au C() u a L 2 () + 2 u ) a L 2 (Q T ) F a L 2 (Q T ), which is equivalent to (27). B Energy for a particular cross-iffusion system A particular choice of cross-iffusion terms in (5) permits for an energy inequality even for negative secon orer matrices. This is the case of the system t U (U(1 + V 2 )) =, t V (V (1 + U 2 (33) )) =, still with Neumann bounary conitions an initial ata U, V. For this system, the energy is given by One can easily check that it hols E(x, t) := (1 + U 2 )(1 + V 2 ). t E(x, t) = 2U(1 + V 2 ) (U(1 + V 2 )) + 2(V (1 + U 2 )) (V (1 + U 2 )), which leas immeiately to (1 + U 2 )(1 + V 2 ) = 2 t (U(1 + V 2 )) 2 2 (V (1 + U 2 )) 2. It follows an a priori estimate in the space L t (L 2 x) that completes the L p tx boun prove in section 3. The system (33) is not always elliptic. This is relate to the non-convexity of this energy (still for large ata), an important ifference with the Shigezaa-Kawasaki prey-preator system which comes with a convex entropy functional ([6]). Acknowlegment This research has been supporte by INRIA-CONICYT scientific cooperation program. References [1] H. Amann, Dynamic theory of quasilinear parabolic equations. I. Abstract evolution equations, Nonlinear Anal., 12 (1988), pp [2], Dynamic theory of quasilinear parabolic equations. II. Reaction-iffusion systems, Differential Integral Equations, 3 (199), pp [3], Erratum: Dynamic theory of quasilinear parabolic systems. III. Global existence [Math. Z. 22 (1989), no. 2, ; MR11386 (9i:35125)], Math. Z., 25 (199), p
20 [4] K. Anguige an C. Schmeiser, A one-imensional moel of cell iffusion an aggregation, incorporating volume filling an cell-to-cell ahesion, J. Math. Biol., (29). [5] N. Bouiba an M. Pierre, Global existence for couple reaction-iffusion systems, J. Math. Anal. Appl., 25 (2), pp [6] L. Chen an A. Jüngel, Analysis of a parabolic cross-iffusion population moel without selfiffusion, J. Differential Equations, 224 (26), pp [7], Analysis of a parabolic cross-iffusion semiconuctor moel with electron-hole scattering, Comm. Partial Differential Equations, 32 (27), pp [8] P. Degon, S. Génieys, an A. Jüngel, Symmetrization an entropy inequality for general iffusion equations, C. R. Aca. Sci. Paris Sér. I Math., 325 (1997), pp [9] L. Desvillettes, K. Fellner, M. Pierre, an J. Vovelle, Global existence for quaratic systems of reaction-iffusion, Av. Nonlinear Stu., 7 (27), pp [1] L. C. Evans, A survey of entropy methos for partial ifferential equations, Bull. Amer. Math. Soc. (N.S.), 41 (24), pp (electronic). [11] L. C. Evans an M. Portilheiro, Irreversibility an hysteresis for a forwar-backwar iffusion equation, Math. Moels Methos Appl. Sci., 14 (24), pp [12] L. Hajaj, K. Hamache, an H. D., An existence result to a strongly couple egenerate system arising in tumor moeling, Abstract an Applie Analysis, (28), pp. ID 23987, 19pp. (electronic). [13] M. Iia, M. Mimura, an H. Ninomiya, Diffusion, cross-iffusion an competitive interaction, J. Math. Biol., 53 (26), pp [14] H. Izuhara an M. Mimura, Reaction-iffusion system approximation to the cross-iffusion competition system, Hiroshima Math. J., 38 (28), pp [15] S. Kawashima an Y. Shizuta, On the normal form of the symmetric hyperbolic-parabolic systems associate with the conservation laws, Tohoku Math. J. (2), 4 (1988), pp [16] t. Lepoutre, Partial ifferential equations arising in biology, PhD Thesis Univ. P. et M. Curie, In preparation. [17] S. Levin an L. Segel, Hypothesis for origin of planktonic patchiness, Nature, 259 (1976), p [18] S. A. Levin, A more functional response to preator-prey stability, The American Naturalist, 18 (1977), pp [19] Y. Li an C. Zhao, Global existence of solutions to a cross-iffusion system in higher imensional omains, Discrete Contin. Dyn. Syst., 12 (25), pp [2] Y. Lou an S. Martinez, Evolution of cross-iffusion an self-iffusion, Preprint, (27), pp
21 [21] Y. Lou an W.-M. Ni, Diffusion, self-iffusion an cross-iffusion, J. Differential Equations, 131 (1996), pp [22] A. Marrocco, Numerical simulation of chemotactic bacteria aggregation via mixe finite elements, M2AN Math. Moel. Numer. Anal., 37 (23), pp [23] M. Mimura an K. Kawasaki, Spatial segregation in competitive interaction-iffusion equations, J. Math. Biol., 9 (198), pp [24] M. Mimura an J. D. Murray, On a iffusive prey-preator moel which exhibits patchiness, J. Theoret. Biol., 75 (1978), pp [25] M. Mimura an M. Yamaguti, Pattern formation in interacting an iffusing systems in population biology., Av Biophys, 15 (1982), pp [26] A. Novick-Cohen an R. L. Pego, Stable patterns in a viscous iffusion equation, Trans. Amer. Math. Soc., 324 (1991), pp [27] A. Okubo, Diffusion an ecological problems: mathematical moels, vol. 1 of Biomathematics, Springer-Verlag, Berlin, 198. An extene version of the Japanese eition, ıt Ecology an iffusion, Translate by G. N. Parker. [28] M. Pierre, Weak solutions an supersolutions in L 1 for reaction-iffusion systems, J. Evol. Equ., 3 (23), pp Deicate to Philippe Bénilan. [29] M. Pierre an D. Schmitt, Blowup in reaction-iffusion systems with issipation of mass, SIAM Rev., 42 (2), pp (electronic). [3] P. I. Plotnikov, Passage to the limit with respect to viscosity in an equation with a variable irection of parabolicity, Differentsialnye Uravneniya, 3 (1994), pp , 734. [31] N. Shigesaa, K. Kawasaki, an E. Teramoto, Spatial segregation of interacting species, J. Theoret. Biol., 79 (1979), pp [32] A. M. Turing, The chemical basis of morphogenesis, Phil. Trans. Roy. Soc. B, 237 (1952), pp [33] Y. Wang, The global existence of solutions for a cross-iffusion system, Acta Math. Appl. Sin. Engl. Ser., 21 (25), pp
On some parabolic systems arising from a nuclear reactor model
On some parabolic systems arising from a nuclear reactor moel Kosuke Kita Grauate School of Avance Science an Engineering, Wasea University Introuction NR We stuy the following initial-bounary value problem
More informationBOUNDEDNESS IN A THREE-DIMENSIONAL ATTRACTION-REPULSION CHEMOTAXIS SYSTEM WITH NONLINEAR DIFFUSION AND LOGISTIC SOURCE
Electronic Journal of Differential Equations, Vol. 016 (016, No. 176, pp. 1 1. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu BOUNDEDNESS IN A THREE-DIMENSIONAL ATTRACTION-REPULSION
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationAbstract A nonlinear partial differential equation of the following form is considered:
M P E J Mathematical Physics Electronic Journal ISSN 86-6655 Volume 2, 26 Paper 5 Receive: May 3, 25, Revise: Sep, 26, Accepte: Oct 6, 26 Eitor: C.E. Wayne A Nonlinear Heat Equation with Temperature-Depenent
More informationThe effect of dissipation on solutions of the complex KdV equation
Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,
More informationChaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena
Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos
More informationOn the Cauchy Problem for Von Neumann-Landau Wave Equation
Journal of Applie Mathematics an Physics 4 4-3 Publishe Online December 4 in SciRes http://wwwscirporg/journal/jamp http://xoiorg/436/jamp4343 On the Cauchy Problem for Von Neumann-anau Wave Equation Chuangye
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationA nonlinear inverse problem of the Korteweg-de Vries equation
Bull. Math. Sci. https://oi.org/0.007/s3373-08-025- A nonlinear inverse problem of the Korteweg-e Vries equation Shengqi Lu Miaochao Chen 2 Qilin Liu 3 Receive: 0 March 207 / Revise: 30 April 208 / Accepte:
More informationSINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationLocal well-posedness of nonlocal Burgers equations
Local well-poseness of nonlocal Burgers equations Sylvie Benzoni-Gavage To cite this version: Sylvie Benzoni-Gavage. Local well-poseness of nonlocal Burgers equations. Differential Integral Equations,
More informationGlobal Solutions to the Coupled Chemotaxis-Fluid Equations
Global Solutions to the Couple Chemotaxis-Flui Equations Renjun Duan Johann Raon Institute for Computational an Applie Mathematics Austrian Acaemy of Sciences Altenbergerstrasse 69, A-44 Linz, Austria
More informationNONLINEAR QUARTER-PLANE PROBLEM FOR THE KORTEWEG-DE VRIES EQUATION
Electronic Journal of Differential Equations, Vol. 11 11), No. 113, pp. 1. ISSN: 17-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu NONLINEAR QUARTER-PLANE PROBLEM
More informationError estimates for 1D asymptotic models in coaxial cables with non-homogeneous cross-section
Error estimates for 1D asymptotic moels in coaxial cables with non-homogeneous cross-section ébastien Imperiale, Patrick Joly o cite this version: ébastien Imperiale, Patrick Joly. Error estimates for
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More informationGLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS
Electronic Journal of Differential Equations, Vol. 015 015), No. 99, pp. 1 14. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu GLOBAL SOLUTIONS FOR D COUPLED
More informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More informationWELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL PRESSURE IN SOBOLEV SPACES
Electronic Journal of Differential Equations, Vol. 017 (017), No. 38, pp. 1 7. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL
More information3 The variational formulation of elliptic PDEs
Chapter 3 The variational formulation of elliptic PDEs We now begin the theoretical stuy of elliptic partial ifferential equations an bounary value problems. We will focus on one approach, which is calle
More informationMARKO NEDELJKOV, DANIJELA RAJTER-ĆIRIĆ
GENERALIZED UNIFORMLY CONTINUOUS SEMIGROUPS AND SEMILINEAR HYPERBOLIC SYSTEMS WITH REGULARIZED DERIVATIVES MARKO NEDELJKOV, DANIJELA RAJTER-ĆIRIĆ Abstract. We aopt the theory of uniformly continuous operator
More informationExistence and Uniqueness of Solution for Caginalp Hyperbolic Phase Field System with Polynomial Growth Potential
International Mathematical Forum, Vol. 0, 205, no. 0, 477-486 HIKARI Lt, www.m-hikari.com http://x.oi.org/0.2988/imf.205.5757 Existence an Uniqueness of Solution for Caginalp Hyperbolic Phase Fiel System
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationASYMPTOTICS TOWARD THE PLANAR RAREFACTION WAVE FOR VISCOUS CONSERVATION LAW IN TWO SPACE DIMENSIONS
TANSACTIONS OF THE AMEICAN MATHEMATICAL SOCIETY Volume 35, Number 3, Pages 13 115 S -9947(999-4 Article electronically publishe on September, 1999 ASYMPTOTICS TOWAD THE PLANA AEFACTION WAVE FO VISCOUS
More informationVarious boundary conditions for Navier-Stokes equations in bounded Lipschitz domains
Various bounary conitions for Navier-Stokes equations in boune Lipschitz omains Sylvie Monniaux To cite this version: Sylvie Monniaux. Various bounary conitions for Navier-Stokes equations in boune Lipschitz
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationOn a class of nonlinear viscoelastic Kirchhoff plates: well-posedness and general decay rates
On a class of nonlinear viscoelastic Kirchhoff plates: well-poseness an general ecay rates M. A. Jorge Silva Department of Mathematics, State University of Lonrina, 8657-97 Lonrina, PR, Brazil. J. E. Muñoz
More informationMathematical methods in biology Master of Mathematics, second year
Mathematical methos in biology Master of Mathematics, secon year Benoît Perthame Sorbonne Université, UMR 7598, Laboratoire J.-L. Lions, BC187 4, place Jussieu, F-75252 Paris ceex 5 October 26, 2018 2
More informationarxiv: v1 [math.ap] 6 Jul 2017
Local an global time ecay for parabolic equations with super linear first orer terms arxiv:177.1761v1 [math.ap] 6 Jul 17 Martina Magliocca an Alessio Porretta ABSTRACT. We stuy a class of parabolic equations
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationTHE KELLER-SEGEL MODEL FOR CHEMOTAXIS WITH PREVENTION OF OVERCROWDING: LINEAR VS. NONLINEAR DIFFUSION
THE KELLER-SEGEL MODEL FOR CHEMOTAXIS WITH PREVENTION OF OVERCROWDING: LINEAR VS. NONLINEAR DIFFUSION MARTIN BURGER, MARCO DI FRANCESCO, AND YASMIN DOLAK Abstract. The aim of this paper is to iscuss the
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationII. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More informationA transmission problem for the Timoshenko system
Volume 6, N., pp. 5 34, 7 Copyright 7 SBMAC ISSN -85 www.scielo.br/cam A transmission problem for the Timoshenko system C.A. RAPOSO, W.D. BASTOS an M.L. SANTOS 3 Department of Mathematics, UFSJ, Praça
More informationSecond order differentiation formula on RCD(K, N) spaces
Secon orer ifferentiation formula on RCD(K, N) spaces Nicola Gigli Luca Tamanini February 8, 018 Abstract We prove the secon orer ifferentiation formula along geoesics in finite-imensional RCD(K, N) spaces.
More informationArithmetic Distributions of Convergents Arising from Jacobi-Perron Algorithm
Arithmetic Distributions of Convergents Arising from Jacobi-Perron Algorithm Valerie Berthe, H Nakaa, R Natsui To cite this version: Valerie Berthe, H Nakaa, R Natsui Arithmetic Distributions of Convergents
More informationTWO-SPECIES COMPETITION WITH HIGH DISPERSAL: THE WINNING STRATEGY. Stephen A. Gourley. (Communicated by Yasuhiro Takeuchi)
MATHEMATICAL BIOSCIENCES http://www.mbejournal.org/ AND ENGINEERING Volume, Number, April 5 pp. 345 36 TWO-SPECIES COMPETITION WITH HIGH DISPERSAL: THE WINNING STRATEGY Stephen A. Gourley Department of
More informationA new proof of the sharpness of the phase transition for Bernoulli percolation on Z d
A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationA REMARK ON THE DAMPED WAVE EQUATION. Vittorino Pata. Sergey Zelik. (Communicated by Alain Miranville)
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 5, Number 3, September 2006 pp. 6 66 A REMARK ON THE DAMPED WAVE EQUATION Vittorino Pata Dipartimento i Matematica F.Brioschi,
More informationSuppressing Chemotactic Blow-Up Through a Fast Splitting Scenario on the Plane
Arch. Rational Mech. Anal. Digital Object Ientifier (DOI) https://oi.org/1.17/s25-18-1336-z Suppressing Chemotactic Blow-Up Through a Fast Splitting Scenario on the Plane Siming He & Eitan Tamor Communicate
More informationL p Theory for the Multidimensional Aggregation Equation
L p Theory for the Multiimensional Aggregation Equation ANDREA L. BERTOZZI University of California - Los Angeles THOMAS LAURENT University of California - Los Angeles AND JESÚS ROSADO Universitat Autònoma
More informationHYPOCOERCIVITY WITHOUT CONFINEMENT. 1. Introduction
HYPOCOERCIVITY WITHOUT CONFINEMENT EMERIC BOUIN, JEAN DOLBEAULT, STÉPHANE MISCHLER, CLÉMENT MOUHOT, AND CHRISTIAN SCHMEISER Abstract. Hypocoercivity methos are extene to linear kinetic equations with mass
More informationAn Optimal Algorithm for Bandit and Zero-Order Convex Optimization with Two-Point Feedback
Journal of Machine Learning Research 8 07) - Submitte /6; Publishe 5/7 An Optimal Algorithm for Banit an Zero-Orer Convex Optimization with wo-point Feeback Oha Shamir Department of Computer Science an
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationINVERSE PROBLEM OF A HYPERBOLIC EQUATION WITH AN INTEGRAL OVERDETERMINATION CONDITION
Electronic Journal of Differential Equations, Vol. 216 (216), No. 138, pp. 1 7. ISSN: 172-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu INVERSE PROBLEM OF A HYPERBOLIC EQUATION WITH AN
More informationMATHEMATICAL ANALYSIS OF A PARABOLIC-ELLIPTIC MODEL FOR BRAIN LACTATE KINETICS
MATHEMATICAL ANALYSIS OF A PARABOLIC-ELLIPTIC MODEL FOR BRAIN LACTATE KINETICS Alain Miranville To cite this version: Alain Miranville. MATHEMATICAL ANALYSIS OF A PARABOLIC-ELLIPTIC MODEL FOR BRAIN LACTATE
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationTHE RELAXATION SPEED IN THE CASE THE FLOW SATISFIES EXPONENTIAL DECAY OF CORRELATIONS
HE RELAXAIO SPEED I HE CASE HE FLOW SAISFIES EXPOEIAL DECAY OF CORRELAIOS Brice Franke, hi-hien guyen o cite this version: Brice Franke, hi-hien guyen HE RELAXAIO SPEED I HE CASE HE FLOW SAISFIES EXPOEIAL
More informationON ISENTROPIC APPROXIMATIONS FOR COMPRESSIBLE EULER EQUATIONS
ON ISENTROPIC APPROXIMATIONS FOR COMPRESSILE EULER EQUATIONS JUNXIONG JIA AND RONGHUA PAN Abstract. In this paper, we first generalize the classical results on Cauchy problem for positive symmetric quasilinear
More informationLOCAL SOLVABILITY AND BLOW-UP FOR BENJAMIN-BONA-MAHONY-BURGERS, ROSENAU-BURGERS AND KORTEWEG-DE VRIES-BENJAMIN-BONA-MAHONY EQUATIONS
Electronic Journal of Differential Equations, Vol. 14 (14), No. 69, pp. 1 16. ISSN: 17-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu LOCAL SOLVABILITY AND BLOW-UP
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationResearch Article Global and Blow-Up Solutions for Nonlinear Hyperbolic Equations with Initial-Boundary Conditions
International Differential Equations Volume 24, Article ID 724837, 5 pages http://x.oi.org/.55/24/724837 Research Article Global an Blow-Up Solutions for Nonlinear Hyperbolic Equations with Initial-Bounary
More informationReactive Power Compensation in Mechanical Systems
Reactive Power Compensation in Mechanical Systems Carlos Rengifo, Bassel Kaar, Yannick Aoustin, Christine Chevallereau To cite this version: Carlos Rengifo, Bassel Kaar, Yannick Aoustin, Christine Chevallereau.
More informationAdaptive Control of the Boost DC-AC Converter
Aaptive Control of the Boost DC-AC Converter Carolina Albea-Sanchez, Carlos Canuas De Wit, Francisco Gorillo Alvarez To cite this version: Carolina Albea-Sanchez, Carlos Canuas De Wit, Francisco Gorillo
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationLocal Input-to-State Stabilization of 1-D Linear Reaction-Diffusion Equation with Bounded Feedback
Local Input-to-State Stabilization of -D Linear Reaction-Diffusion Equation with Boune Feeback Aneel Tanwani, Swann Marx, Christophe Prieur To cite this version: Aneel Tanwani, Swann Marx, Christophe Prieur.
More informationON A NONLOCAL SELECTION-MUTATION MODEL WITH A GRADIENT FLOW STRUCTURE
ON A NONLOCAL SELECTION-MUTATION MODEL WITH A GRADIENT FLOW STRUCTURE PIERRE-EMMANUEL JABIN AND HAILIANG LIU Abstract. In this paper, we are intereste in an integro-ifferential moel with a nonlinear competition
More informationDiscrete Operators in Canonical Domains
Discrete Operators in Canonical Domains VLADIMIR VASILYEV Belgoro National Research University Chair of Differential Equations Stuencheskaya 14/1, 308007 Belgoro RUSSIA vlaimir.b.vasilyev@gmail.com Abstract:
More informationAll s Well That Ends Well: Supplementary Proofs
All s Well That Ens Well: Guarantee Resolution of Simultaneous Rigi Boy Impact 1:1 All s Well That Ens Well: Supplementary Proofs This ocument complements the paper All s Well That Ens Well: Guarantee
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationRobust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k
A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationTractability results for weighted Banach spaces of smooth functions
Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March
More informationMath 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions
Math 3 Winter 2 Avance Bounary Value Problems I Bessel s Equation an Bessel Functions Department of Mathematical an Statistical Sciences University of Alberta Bessel s Equation an Bessel Functions We use
More informationWell-posedness of hyperbolic Initial Boundary Value Problems
Well-poseness of hyperbolic Initial Bounary Value Problems Jean-François Coulombel CNRS & Université Lille 1 Laboratoire e mathématiques Paul Painlevé Cité scientifique 59655 VILLENEUVE D ASCQ CEDEX, France
More informationAsymptotic estimates on the time derivative of entropy on a Riemannian manifold
Asymptotic estimates on the time erivative of entropy on a Riemannian manifol Arian P. C. Lim a, Dejun Luo b a Nanyang Technological University, athematics an athematics Eucation, Singapore b Key Lab of
More informationEnergy behaviour of the Boris method for charged-particle dynamics
Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationTHE VLASOV-MAXWELL-BOLTZMANN SYSTEM NEAR MAXWELLIANS IN THE WHOLE SPACE WITH VERY SOFT POTENTIALS
THE VLASOV-MAXWELL-BOLTZMANN SYSTEM NEAR MAXWELLIANS IN THE WHOLE SPACE WITH VERY SOFT POTENTIALS RENJUN DUAN, YUANJIE LEI, TONG YANG, AND HUIJIANG ZHAO Abstract. Since the work [4] by Guo [Invent. Math.
More informationExistence and Uniqueness of Solutions of Mathematical Models of Predator Prey Interactions
Punjab University Journal of Mathematics (ISSN 1016-56) Vol. 49()(017) pp. 75-87 Existence an Uniqueness of Solutions of Mathematical Moels of Preator Prey Interactions Muhamma Shakil 1, Hafiz Abul Wahab
More informationGeneralized Tractability for Multivariate Problems
Generalize Tractability for Multivariate Problems Part II: Linear Tensor Prouct Problems, Linear Information, an Unrestricte Tractability Michael Gnewuch Department of Computer Science, University of Kiel,
More informationGradient flow of the Chapman-Rubinstein-Schatzman model for signed vortices
Graient flow of the Chapman-Rubinstein-Schatzman moel for signe vortices Luigi Ambrosio, Eoaro Mainini an Sylvia Serfaty Deicate to the memory of Michelle Schatzman (1949-2010) Abstract We continue the
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More informationA COMBUSTION MODEL WITH UNBOUNDED THERMAL CONDUCTIVITY AND REACTANT DIFFUSIVITY IN NON-SMOOTH DOMAINS
Electronic Journal of Differential Equations, Vol. 2929, No. 6, pp. 1 14. ISSN: 172-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu A COMBUSTION MODEL WITH UNBOUNDED
More informationOn a limit theorem for non-stationary branching processes.
On a limit theorem for non-stationary branching processes. TETSUYA HATTORI an HIROSHI WATANABE 0. Introuction. The purpose of this paper is to give a limit theorem for a certain class of iscrete-time multi-type
More informationarxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003
Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep
More information2 GUANGYU LI AND FABIO A. MILNER The coefficient a will be assume to be positive, boune, boune away from zero, an inepenent of t; c will be assume con
A MIXED FINITE ELEMENT METHOD FOR A THIRD ORDER PARTIAL DIFFERENTIAL EQUATION G. Li 1 an F. A. Milner 2 A mixe finite element metho is escribe for a thir orer partial ifferential equation. The metho can
More informationMany problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx
Math 53 Notes on turm-liouville equations Many problems in physics, engineering, an chemistry fall in a general class of equations of the form w(x)p(x) u ] + (q(x) λ) u = w(x) on an interval a, b], plus
More informationarxiv: v1 [cond-mat.stat-mech] 9 Jan 2012
arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationARBITRARY NUMBER OF LIMIT CYCLES FOR PLANAR DISCONTINUOUS PIECEWISE LINEAR DIFFERENTIAL SYSTEMS WITH TWO ZONES
Electronic Journal of Differential Equations, Vol. 2015 (2015), No. 228, pp. 1 12. ISSN: 1072-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu ARBITRARY NUMBER OF
More informationCOUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS
VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND
More informationIERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210
IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.210 Discrete an Continuous Dynamics in Nonlinear Monopolies Akio Matsumoto Chuo University Ferenc
More informationMartin Luther Universität Halle Wittenberg Institut für Mathematik
Martin Luther Universität alle Wittenberg Institut für Mathematik Weak solutions of abstract evolutionary integro-ifferential equations in ilbert spaces Rico Zacher Report No. 1 28 Eitors: Professors of
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationBrooke L. Hollingsworth and R. E. Showalter Department of Mathematics The University of Texas at Austin Austin, TX USA
SEMILINEAR DEGENERATE PARABOLIC SYSTEMS AND DISTRIBUTED CAPACITANCE MODELS Brooke L. Hollingsworth an R. E. Showalter Department of Mathematics The University of Texas at Austin Austin, TX 7872 USA Abstract.
More informationSEMILINEAR DEGENERATE PARABOLIC SYSTEMS AND DISTRIBUTED CAPACITANCE MODELS. Brooke L. Hollingsworth and R.E. Showalter
DISCRETE AND CONTINUOUS DYNAMICAL SYSTEMS Volume, Number, January 995 pp. 59 76 SEMILINEAR DEGENERATE PARABOLIC SYSTEMS AND DISTRIBUTED CAPACITANCE MODELS Brooke L. Hollingsworth an R.E. Showalter Department
More informationThe Generalized Incompressible Navier-Stokes Equations in Besov Spaces
Dynamics of PDE, Vol1, No4, 381-400, 2004 The Generalize Incompressible Navier-Stokes Equations in Besov Spaces Jiahong Wu Communicate by Charles Li, receive July 21, 2004 Abstract This paper is concerne
More informationLeChatelier Dynamics
LeChatelier Dynamics Robert Gilmore Physics Department, Drexel University, Philaelphia, Pennsylvania 1914, USA (Date: June 12, 28, Levine Birthay Party: To be submitte.) Dynamics of the relaxation of a
More information