arxiv: v1 [math.ap] 26 Mar 2014
|
|
- Claud Fowler
- 6 years ago
- Views:
Transcription
1 Stability and Optimization in Structured Population Models on Graphs Rinaldo M. Colombo 1 Mauro Garavello arxiv: v1 [math.ap] 6 Mar 14 October 9, 17 Abstract We prove existence and uniqueness of solutions, continuous dependence from the initial datum and stability with respect to the boundary condition in a class of initial boundary value problems for systems of balance laws. The particular choice of the boundary condition allows to comprehend models with very different structures. In particular, we consider a juvenileadult model, the problem of the optimal mating ratio and a model for the optimal management of biological resources. The stability result obtained allows to tackle various optimal management/control problems, providing sufficient conditions for the existence of optimal choices/controls. Keywords: renewal equation; balance laws; juvenile-adult model; management of biological resources; optimal mating ratio. 1 MSC: 35L5, 9D5 1 Introduction This paper is devoted to the following initial boundary value problem for a system of balance laws in one space dimension: t u i + x gi t, x u i = di t, x u i t, x R + R + g i t, u i t, + = B i t, u1 t,..., u n t t R + u i, x = u o i x x R + i = 1,..., n. 1.1 Here, i = 1,..., n and t R + is time. The space variable x varies in R + and in the applications of 1.1 will have the meaning of a biological age, or size. The scalar functions g 1,..., g n are growth functions, d 1,..., d n are the death rates and u o 1,..., u o n constitute the initial data. A key role is played by our choice of the birth function B i, for i = 1,..., n, which we assume of the form B i t, u 1,..., u n = α i t, u1 x 1,..., u n x n + β i u 1 x dx,..., u n x dx 1. I 1 I n for suitable functions α i, β i, points x i > and measurable I i R +, for i = 1,..., n. The literature on equations similar to 1.1 is vast. We refer for instance to the exhaustive monograph [14] or to the more recent edition of [5] and to the references therein. Specific features of 1.1 are that it is a system, boundary conditions may contain both a local term, the α i, and a nonlocal term, the β i. 1 INDAM Unit, University of Brescia Department of Mathematics and Applications, University of Milano Bicocca 1
2 From the analytical point of view, in the present treatment we emphasize the role of the total variation, setting the main result in BV. In particular, this allows to consider a function of the type 1. and to prove that the boundary data are attained in the sense of traces, also due to the boundary being non characteristic. In this setting, the stability of solutions with respect to α i and β i is also obtained. Moreover, the techniques used in the sequel can easily be extended to more general source terms as well as to situations where also the space distribution needs to be taken into account. From the modeling point of view, the use of boundary conditions of the type 1. unifies the treatment of rather diverse situations. First, it comprises the standard case always covered in the literature on renewal equations, where the independent variable x varies along a segment or a half line, see Figure 1, left. The dependent variable u represents the population density that at time t j a Figure 1: Biological structures comprised in Left, a standard linear setting and, right, a juvenile-adult situation. is of size or age x. A more complicate structure was recently considered in [1], see Figure 1, right. There, the size/age biological variable varies along a graph consisting of distinct sets, corresponding to the juvenile and to the adult stages in the development of the considered species. Here, we are able to deal also with this situation, as depicted in Figure, right. m j s f r Figure : Biological structures comprised in Left, a framework corresponding, for instance, to sexual reproduction: the two branches correspond to males M and females F. Right, a structure possibly accounting for the exploitation of biological resources: when juveniles reach the adult stage, they are split into a part S which is, say, sold and a part R used for reproduction. The finite propagation speed intrinsic to models of the type 1.1 clearly allows to combine various instances of the graphs above. Other situations of biological interest can be for instance a three stage linear structure or a tree shaped one, see Figure 3. These schemes, as well as many others, all fit into the scope of Theorem.4 below. In this connection, we recall that similar network structures are widely considered in the framework of vehicular traffic modeling, see [9]. In the case of nonlinear systems of balance laws with flow independent from the space variable, the initial boundary value problem has been widely investigated, see for instance [7]. For the relations between the problems with boundaries and with junction see [8, Proposition 4.]. The present treatment is self-contained. Section is devoted to the analytical results. Specific applications are in Section 3, where sample numerical integrations are also provided. All technical details are deferred to Section 4.
3 Figure 3: General graphs for further biological structures comprised in Analytical Results Throughout, we use the standard notation R + = [, + [ and R + = ], + [. When A and B are suitable subsets of R m, C A; B, respectively C,1 A; B, L 1 A; B or L A; B, is the set of continuous, respectively Lipschitz continuous, Lebesgue integrable or essentially bounded, maps defined on A and attaining values in B. For the basic theory of BV functions we refer to [3]. When referring to a function u: R + R + R, the first argument is time, the second is the biological age/size variable. If I R + is an interval, we denote TV ut, N = sup ut, xh ut, x h 1 : N N, xh 1 < x h, h = 1,..., N h=1 TV ut, ; I N = sup ut, xh ut, x h 1 : N N, xh 1 < x h, x h I, h = 1,..., N h=1 TV u, x; I N = sup uth, x ut h 1, x : N N, th 1 < t h, t h I, h = 1,..., N. h=1 Preliminarily, we consider the following initial boundary value problem for a linear scalar balance law, or renewal equation in [14, Chapter 3]: under the following assumptions b b BV loc R + ; R ; t u + x gt, x u = dt, x u t, x R + R + u, x = u o x x R +.1 gt, ut, + = bt t R + g g C 1 R + R + ; [, ĝ] for positive, ĝ and d d C 1 L R + R + ; R, sup t R + TV dt, < +. { sup t R + TV gt, < + sup t R + TV x gt, < + ; The solutions to.1 can be written in terms of the ordinary differential equation ẋ = gt, x. If g satisfies g, we can introduce the globally defined maps { { ẋ = gt, x t = 1 t Xt; t o, x o that solves and x T x; t xt o = x o, x o that solves gt,x. o tx o = t o. Denote γt = Xt;,, its inverse being t = Γx. Note that if x γt then X; t, x [, x] and if x < γt then T ; t, x [, t]. Recall the following definition of solution to.1, see also [4, 6, 11, 14, 18]. 3
4 Definition.1. Let b, g and d hold. Choose an initial datum u o L 1 R + ; R. The function u C R + ; L 1 R + ; R is a solution to.1 if 1. for all ϕ C 1 c R + R + ; R, R +. u = u o R + ut, x t ϕt, x + gt, x ut, x x ϕt, x + dt, x ut, x ϕt, x dt dx = 3. for a.e. t R +, lim x + gt, x ut, x = bt. The following Lemma summarizes various properties of the solution to.1, see also [14]. Here, we stress the role of BV estimates. The proof is deferred to Section 4. Lemma.. Let b, g and d hold. Then, for any u o L 1 BVR + ; R, the map u: R + R + R defined by t u o X; t, x exp dτ, Xτ; t, x x g τ, Xτ; t, x dτ x > γt ut, x = bt ;t,x gt ;t,x, exp t dτ, Xτ; t, x T ;t,x x g τ, Xτ; t, x.3 dτ x < γt solves.1 in the sense of Definition.1. Moreover, there exists a constant C dependent only on g and d, see 4.6, such that the following a priori estimates hold for all t R + : ut L u R + ;R o L R + ;R + 1 b L [,t];r e Ct.4 ut L u 1 R + ;R o L1 R + ;R + 1 b L 1 [,t];r e Ct.5 TV ut [ u o L R + ;R + TVu o + C + b L [,t];r + 1 ] TVb; [, t] e Ct.6 u o TV u, x; [, t] L [X;t,x,x];R + TVu o; [X; t, x, x] e Ct x > γt u o L [,x];r + TVu o; [, x] e Ct.7 1 C t + TVb; [, T ; t, x] + b L [,T ;t,x];r e Ct x < γt. Moreover, for any interval I R +, TV u, x dx ; [, t] C I uτ L I;R + TV uτ, ; I dτ..8 For every t R +, there exists a positive L dependent on, ĝ, C and TVb; [, t], b L [,t];r, such that, for t, t [, t], ut ut L1 R + ;R L t t..9 For u o, u o L 1 BVR + ; R and b, b as in b, the solutions u and u to t u + x gt, x u = dt, x u t u + x gt, x u = dt, x u u, x = u ox gt, ut, + = b t and u, x = u ox gt, ut, + = b t.1 4
5 satisfy the stability and monotonicity estimates u t u t [ u L 1 R + ;R o u o L 1 R + ;R + 1 b b ] L e Ct,.11 1 [,t];r u ox u ox for all x R + b t b t for all t R + u t, x u t, x for all t, x R + R +..1 Recall that in the present case of a linear conservation law, the definition of weak solution at. is equivalent to the definition of Kružkov solution [11, Definition 1]. It is immediate to verify that for u o = and b =, problem.3 admits the solution u =. Hence, the monotonicity property.1 also ensures that non-negative initial and boundary data in lead to non-negative solutions. In order to pass to system 1.1, we need the following notation for norms and total variations of functions attaining values in R n : u L 1 R + ;R n = n i=1 n u i L 1 R + ;R, u L R + ;R n = n u i L R + ;R, TVu = TVu i. As a reference for the usual definition of weak solution to scalar conservation laws, see [4, 11]. Definition.3. Let T >. Consider 1.1 with g 1,..., g n satisfying assumptions g and d 1,..., d n satisfying d. Fix an initial datum u o L 1 BVR + ; R n. A map u C [, T ]; L 1 BVR + ; R n is a solution to if, setting b i t = α i t, u1 t, x 1,..., u n t, x n + β i u 1 t, x dx,..., u n t, x dx I 1 I n for all i = 1,..., n, the i-th component u i is a solution to t u i + x gi t, x u i = di t, x u i t, x R + R + u i, x = u i ox x R +.13 g i t, u i t, + = b i t t R + in the sense of Definition.1. The following result ensures the well posedness of Its proof is presented in Section 4. Theorem.4. Let n N\{, x 1,..., x n R +, g 1,..., g n satisfy g and d 1,..., d n satisfy d. Assume that the maps α α 1,..., α n and β β 1,..., β n satisfy i=1 α C,1 R + R n ; R n, β C,1 R n ; R n and αt, = β =..14 Then, for any u o L 1 BVR + ; R n, the problem 1.1 admits a unique solution in the sense of Definition.3. Moreover, there exists an increasing function K C R + ; R + dependent only on Lipα, Lipβ and on C in 4.6 such that for any initial data u o, u o L 1 BVR + ; R n, the corresponding solutions satisfy u t u t u L Kt 1 R + ;R n o u o L 1 R + ;R m + t u o u o L,.15 R + ;R m u t u t u L Kt R + ;R n o u L 1 R + ;R m + u o u..16 Moreover, if u o =, then the solution is ut, x = for all t, x R + R +. o i=1 o L R + ;R m 5
6 We now state separately the stability of solutions to with respect to the birth function B. This result plays a key role in the optimization problems considered below. Theorem.5. Let both systems t u i + x gi t, x u i = di t, x u i g i t, u i t, + = B i t, u1 t,..., u n t u i, x = u o i x t u i + x gi t, x u i = di t, x u i g i t, u i t, + = B i t, u1 t,..., u n t u i, x = u o i x.17 with B it, u 1,..., u n = α i t, u1 x 1,..., u n x n + β i u 1 x dx,..., u n x dx, I 1 I n B i t, u 1,..., u n = α i t, u1 x 1,..., u n x n + β i u 1 x dx,..., u n x dx I 1 I n satisfy the assumptions of Theorem.4. Then, the corresponding solutions u and u are such that u t u t Ht L1 α α R + ;R n C R + R n ;R n + β β C.18 R n ;R n where H C R + ; R + is such that H =. The proof is deferred to Section 4. In applications of Theorem.4 to systems motivated by, for instance, structured population biology, further assumptions are natural and lead to further reasonable properties. Proposition.6. Under the assumptions of Theorem.4, if the boundary functions and the initial data are such that uj α i for all i, j = 1,..., n, wj β i for all i, j = 1,..., n, u o i u o i for all i = 1,..., n, then, the corresponding solutions satisfy u i t, x u i t, x for all t, x R+ R + and i = 1,..., n. In particular, if u o i for i = 1,..., n, then u i t, x for i = 1,..., n. The proof follows immediately from Theorem.4 and from.1, hence it is omitted. 3 Applications This section is devoted to sample applications of Theorem.4 and Theorem.5 to models inspired by structured population biology. We selected three cases corresponding to three different graphs, namely those in Figure 1, right, and in Figure. First, the well posedness ensured by Theorem.4 provides a ground for the reliability of each model. Then, the stability result in Theorem.5 allows to consider further problems. On the one hand, it ensures the existence of a choice of parameters in the equations that lead to solution that best approximate a given set of data. On the other hand, it allows to tackle the problem of optimal mating ratio in a population with sexual reproduction. Finally, we consider the problem of the optimal management of a biological resource. In the former case, the presentation is based on [1, ] where a sensitivity analysis for a model belonging to the class is proved. In the latter cases, we provide numerical integrations showing further qualitative properties of the models considered. 6
7 3.1 A Nonautonomous Juvenile Adult Model In , setting n = and with reference to the structure in Figure 1, right, u 1 t, x = Jt, x g 1 t, x = 1 α 1 t, u 1, u = u t, x = At, x + x min g t, x = gt, x + x min α t, u 1, u = u 1 x 1 = a max d 1 t, x = νt, x β 1 w 1, w = w x = d t, x = µt, x + x min β w 1, w = 3.1 with moreover I = [, x max x min ], we recover [1, Formula.1] in the case β = 1, which we state here for completeness: t J + a J = νt, a J t, a R + [, a max ] t A + x gt, x A = µt, x A t, x R + [x min, x max ] xmax Jt, = At, x dx t R + x min 3. gt, x min At, x min = Jt, a max t R + J, a = J o a a [, a max ] A, x = A o x x [x min, x max ]. Theorem.4 then applies and ensures the well posedness of 3. under assumptions slightly different from those in [1]. Corollary 3.1. In 3., assume that ν C 1 L R + [, a max ]; R and sup t R + TV νt, < + µ C 1 L R + [x min, x max ]; R and sup t R + TV µt, { < + g C 1 R + sup [x min, x max ]; [, ĝ] and t R + TV gt, < + sup t R + TV x gt, < + J o BV[, a max ]; R + A o BV[x min, x max ]; R +. Then, problem 3. admits a unique solution in the sense of Definition.3, the continuous dependence estimates and the stability estimate.18 apply. For completeness, we remark that the model in [1] contains the following slightly more general boundary inflow: Jt, = xmax x min βt, x At, x dx. As soon as β C 1 R + [x min, x max ]; [ ˇβ, + [ for a suitable ˇβ >, the change of variables At, x = βt, x At, x 3.3 still allows to apply Theorem.4. Indeed, with this variable, the second equation in 3. becomes t A + x gt, x A = t βt, x + gt, x x βt, x µt, x A, which is again of the type.1 and hence Theorem.4 can still be applied. The stability proved above allows to tackle the problem of parameter identification. Indeed, through a Weierstraß argument based on Theorem.5, one can prove the existence of a set of parameters in 3. that minimizes a continuous functional representing the distance between the computed solution and a set of experimental data. For a detailed sensitivity analysis for a juvenile adult model we refer to []. 7
8 3. Optimal Mating Ratio Consider a species consisting of males and females, whose densities at time t and age a are described through the functions M = Mt, a and F = F t, a on a structure as that in Figure, left. A natural model is then t M + a M = κ µ M t F + a F = 1 κ µ F{ m f Mt, + F t, = ν min ϑ Mt, a da, 1 ϑ F t, a da m 1 f η Mt, = 1 η F t, M, a = M o a F, a = F o a. Here, κ µ, respectively 1 κ µ, is the mortality rate of males, respectively females, with µ > and κ [, 1]. The positive parameter η [, 1] defines the ratio of male to female newborns, in the sense that every η M males, 1 η F females are born. The constant ν is the fertility rate. We describe the mating ratio at age a through the parameter ϑ, with ϑ [, 1] as follows. The fertile ages are those in the intervals [m 1, m ] for males and [f 1, f ] for females, where m 1, m, f 1, f are positive constants. According to 3.4, all individuals in their fertile age might contribute to reproduction provided the condition imposed by the presence of the mating ratio ϑ is met. If ϑ m m 1 Mt, a da exceeds 1 ϑ f f 1 F t, a da, then only 1 ϑ f ϑ f 1 F t, a da males contribute to the overall population s fertility. Problem 3.4 fits into setting u 1 = M g 1 = 1 d 1 t, x = κ u = F g = 1 d t, x = 1 κ µ I 1 = [m 1, m ] α 1 = β 1 w 1, w = 1 η ν min{ϑ w 1, 1 ϑw I = [f 1, f ] α = β w 1, w = η ν min{ϑ w 1, 1 ϑw. Corollary 3.. Let µ, ν R + ; η, ϑ, κ [, 1]; m 1, m, f 1, f R + with m 1 < m and f 1 < f. For M o, F o L 1 BVR + ; R, problem 3. has a unique solution in the sense of Definition.3, the continuous dependence estimates and the stability estimate.18 apply. The proof is immediate and, hence, omitted. Here, we note that the presence of C 1 positive weights in the integrands defining the boundary data can be recovered through a change of variables entirely similar to that in 3.3. A first immediate property of the solutions to 3.4 is that a zero initial density in either of the two sexes leads to the extinction of the other at exponential speed. Several different optimization problems can be tackled in the framework of 3.4. It is possible to investigate the relations between the parameters κ identifying relative mortality, η the relative natality and ϑ the mating ratio. Below, we look for the optimal mating ratio for given relative natality and mortality coefficients. To this aim, consider the instantaneous average fertility rate over the fertile population { ν min ϑ m m 1 Mt, a da, 1 ϑ f f 1 F t, a da R = m m 1 Mt, a da + f. 3.5 f 1 F t, a da Remark that the functions M and F in 3.5 are solutions to 3.4, hence they depend on the mating ratio ϑ that enters the boundary condition throughout the time interval [, t]. It is natural to assume that a key role is played by the maximal value of R, which is obtained by the choice f f ϑ = 1 F t, a da m m 1 Mt, a da + f f 1 F t, a da
9 Remarkably, this leads to the maximal fertility rate m f m 1 Mt, a da f 1 R = ν F t, a da m m 1 Mt, a da + f f 1 F t, a da coherently with the classical harmonic mean law, see [1, 13, 15, 16, 17, 19]. On the other hand, the right hand sides in 3.5 and 3.6 are time dependent and it can be hardly accepted that ϑ is instantaneously adjusted to the value that maximizes R. More reasonably, one may imagine that ϑ is optimal 1 over a suitably long time interval. We are thus lead to introduce the utility function { Rϑ; T, M o, F o = 1 T ν min ϑ m m 1 Mt, a da, 1 ϑ f f 1 F t, a da T m m 1 Mt, a da + f dt. 3.7 f 1 F t, a da We thus consider the problem find ϑ that maximizes Rϑ; T, M o, F o. A straightforward corollary of Theorem.5 ensures the existence of one such ϑ. Corollary 3.3. Under the assumptions of Corollary 3., for any T R + and any initial datum M o, F o L 1 BVR + ; R there exists a ϑ ], 1[ such that Rϑ ; T, M o, F o = max ϑ [,1] Rϑ; T, M o, F o. The proof is immediate: thanks to Theorem.5, the function ϑ Rϑ; T, M o, F o is continuous for any choice of T R + and M o, F o L 1 BVR + ; R. By the compactness of [, 1], Weierstraß Theorem ensures the existence of ϑ. Moreover, since R; T, M o, F o = R1; T, M o, F o =, we also have ϑ ], 1[. It can be of interest to note that M and F may well increase exponentially with time, but Rϑ; T, M o, F o [, ν] for all T R + and M o, F o L 1 BVR + ; R. As a specific example, we consider the situation identified by the following choices of functions and parameters in : κ =.6 µ =. m 1 = 18 f 1 = 16 η =.485 ν = 3 m = 6 f = and we consider ϑ as a control parameter in [, 1]. As initial datum we choose M o a = 1 and F o a = 1 for all a. 3.9 The graph of the average fertility rate Rϑ; T, M o, F o as a function of ϑ for T = 5 is in Figure 4. The outcome shows a reasonable qualitative behavior. As ϑ or ϑ 1, the number of newborns goes to ; hence the population extinguish. Near to ϑ =.77 there is an optimal choice for the parameter ϑ with respect to the average fertility rate 3.7, which yields a maximal value of.83, see Figure 5 For completeness, we precise that the numerical integration above was obtained using a Lax Friedrichs algorithm, see [1, 1.5], with space mesh a = Here and in the sequel, optimal is understood in the sense that it is the value that maximizes R. However, an excessive natality rate might turn out to be not optimal from the biological point of view. 9
10 Figure 4: Average fertility rate 3.7 along the solutions to with initial datum 3.9 plotted as a function of ϑ. For ϑ = or ϑ = 1, there is no reproduction and the population extinguishes. For ϑ.77, the average natality rate reaches its maximum value, which is approximately equal to.83. Figure 5: Solutions to The integrals 8 Mt, a da blue and 9 F t, a da green as a function of time. Left, for ϑ =.4647, the population vanishes. Middle, with ϑ =.767, there is an equilibrium and right, for ϑ =.777, the population grows exponentially. The latter value gives the maximal average natality rate, see Figure Management of a Biological Resource In biological resource management, one typically rears/breeds a species up to a suitable stage, then part of the population is sold and part is used for reproduction. The equations comprehend this situation. Indeed, call J = Jt, a the density of the juveniles at time t of age or size a. Juveniles reaching the age/size ā are then selected. The density S = St, a refers to those individuals that are going to be sold, while R = Rt, a stands for the density of those reserved for reproduction purposes. One is thus lead to the following model, defined on the structure in Figure, right: t J + a gj t, a J = d J t, a J t, a R + [, ā] t S + a gs t, a S = d S t, a S t, a R + [ā, + [ t R + a gr t, a R = d R t, a R t, a R + [ā, + [ amax g J t, Jt, = β Rt, x dx t R + ā g S t, ā St, ā = η g J t, ā Jt, ā t R + g R t, ā Rt, ā = 1 η g J t, ā Jt, ā t R + J, a = J o a a [, ā] S, a = S o a a [ā, + [ R, a = R o a a [ā, + [
11 Above, we used the obvious notation for the growth and mortality functions g J, g S, g R and d J, d S, d R. The birth rate is described through the function β. A key role is played by the parameter η [, 1] which quantifies the percentage of juveniles selected for the market. System 3.1 fits into setting u 1 t, s = Jt, x g 1 t, x = g J t, x α 1 t, w 1, w, w 3 = u t, x = St, x + ā g t, x = g S t, x + ā α t, w 1, w, w 3 = η w 1 g 1 t, x 1 u 3 t, x = Rt, x + ā g 3 t, x = g R t, x + ā α 3 t, w 1, w, w 3 = 1 η w 1 g 1 t, x 1 d 1 t, x = d J t, x β 1 w 1, w, w 3 = βw 3 x 1 = ā d t, x = d S t, x + ā β w 1, w, w 3 = I 3 = [ā, a max ] d 3 t, x = d R t, x + ā β 3 w 1, w, w 3 = Corollary 3.4. Let g J, g S, g R satisfy g for suitable, ĝ R + with ĝ > >. Let d J, d S, d R satisfy d. Let β C,1 R + ; R be such that β =. For any η [, 1] and any initial data J o BV[, ā]; R + and S o, R o L 1 BV[ā, + [, R +, system 3.1 admits a unique non negative solution and the stability estimates in Theorem.4 apply. A natural question based on model 3.1 is: find the optimal percentage η of juveniles that have to be chosen for the market. To this aim, we postulate simple, though reasonable, cost and gain functionals [ T ā amax [ Cη; T = C J a Jt, a da + CS a St, a + C R a Rt, a ] ] da dt, ā 3.11 Gη; T = T amax ā Ga St, a da dt. Here, C J a is the unit cost to grow a juvenile at age a, and similarly C R, C S are the costs for the other two groups. The gain obtained selling an adult at age a is Ga. We denoted by J = Jt, a, S = St, a and R = Rt, a the solution to 3.1 with initial datum J o and S o =, R o = and with the selection parameter η. A direct consequence of Theorem.5 is the following corollary. Corollary 3.5. In the same assumptions of Corollary 3.4, for any T R + and any J o, S o, R o L 1 BVR + ; R, there exists an optimal choice η such that Gη ; T Cη ; T = max Gη; T Cη; T. η [,1] The proof relies on Weierstraß Theorem, exactly as that of Corollary 3.3 and is here omitted. As a specific example, we consider the situation identified by the following choices of functions and parameters in : g J t, a = 1 d J t, a = C J t, a = a βw = w g S t, a = 1 d S t, a = a ā C S t, a = Gt, a = 1 g R t, a = 1 d R t, a = a ā C R t, a =.5 [ā, a max ] = [1, ] and we consider η as a control parameter in [, 1]. As initial datum we choose 3.1 J o a = 5, S o a =, R o a = The graph of the cost Gη; T Cη; T see 3.11 for T = 15 with respect to λ is in Figure 6. The outcome shows a reasonable qualitative behavior. As η, nothing is sold, all population members are kept for reproduction, the population increases exponentially as also does the functional Gη; T Cη; T. On the contrary, for η 1, all population members are immediately sold giving a positive gain and the population vanishes as also G C. Near to η.3 there is an optimal balance, given the chosen unitary costs and gain With the chosen parameters, the optimal choice for η is η.3, which yields a gain of about 6.48 at time t = 15. As expected, different choices of η have deep influences on the solutions to , as shown in Figure 7. For completeness, we precise that the numerical integration above was obtained using a Lax Friedrichs algorithm, see [1, 1.5], with space mesh a =.1. 11
12 Cost Cost at time t = G C t η η Figure 6: Cost 3.11 along the solutions to with initial datum Left, as a function of η on the horizontal axis and t on the vertical axis. Right, as a function of η at time t = 15.. Recall that for η =, all individual are kept for reproduction and no one is sold. On the contrary, for η = 1, they are all sold and no one is kept for reproduction. J S t t t a J a S t a J a S t a a R t a R 1 t t t R a a Figure 7: Solutions to Time varies along the vertical axis and a along the horizontal one. Above, η =.3 is near to the optimal choice. Middle, η =.5 and, below, η =.91. 1
13 4 Technical Details Throughout, when BV functions are considered, we refer to a right continuous representative. We now recall the following elementary estimates on BV functions. u BVR + ; R w BVR + ; R TVu w u L R + ;R TVw + TVu w L R + ;R 4.1 f C,1 R; R u BVR + ; R TVf u Lipf TVu 4. u BVR + ; R f BVR + ; [ ˇf u + [ ˇf > TV 1ˇf TVu + 1ˇf f L R + ;R 4.3 u L 1 t R + ; BVR + ; R TV uτ, dτ TV uτ dτ 4.4 u BVR + ; R h L R; R + R + u x + hx ux dx TVu h L R + ;R. 4.5 Inequality 4.1 follows from [3, Formula 3.1]. The definition of total variation directly implies 4., 4.3 and 4.4. For a proof of 4.5 see for instance [6, Lemma.3]. Proof of Lemma.. To verify that.3 solves.1, a standard integration along characteristics is sufficient. The bounds.4 and.5 are an immediate consequence of.3. Passing to the estimates on the total variation, introduce C = max x g L R + R + ;R, tg L R + R + ;R, sup t R + TV gt,, sup t R + TV x gt,, d L R + R + ;R, sup t R + TV dt, which is finite by g and d. Consider now the total variation estimates. Using 4.1, 4., 4.3, 4.4, compute: TV ut = TV ut,, [, γt] + TV ut,, [ γt, + [ 4.6 TV b /g, ; [, t] e Ct + b L [,t];r e Ct + TVu o e Ct + u o L R + ;R ect u o L R + ;R + TVu o C b L [,t];r + TVb; [, t] e Ct proving.6. The bound.6 directly follows from.3, using 4.1, 4. and 4.3. We exploit now 4.4 and [3, Definition 3.4], and in the lines below, for typographical reasons, we denote by J the real interval [ 1, 1]. TV u, x dx ; [, t] I { = sup uτ, x dx t ϕτ dτ : ϕ C 1 c], t[ ; J { sup { = sup I I I uτ, x ψx t ϕτ dx dτ : ϕ C1 c], t[ ; J ψ C 1 c I; J gτ, xuτ, x x ψxϕτ + dτ, xuτ, xψxϕτ dx dτ : ϕ C1 c], t[ ; J ψ C 1 c I; J 13
14 { = sup gτ, x uτ, x x ψx ϕτ dx dτ : ϕ C1 c], t[ ; J I ψ C 1 c I; J { + sup dτ, x uτ, x ψx ϕτ dx dτ : ϕ C1 c], t[ ; J I ψ C 1 c I; J { { sup sup gτ, x uτ, x x ψx dx : ψ C 1 c I; J ϕτ dτ : ϕ C 1 c], t[ ; J { { + sup sup I I dτ, x uτ, x ψx dx : ψ C 1 c I; J ϕτ dτ : ϕ C 1 c], t[ ; J { = sup TV gτ, uτ, ϕτ dτ : ϕ C 1 c], t[ ; J { + sup TV dτ, uτ, ϕτ dτ : ϕ C 1 c], t[ ; J TV gτ, uτ, + TV dτ, uτ, dτ Apply now 4.1 to obtain: TV u, x dx ; [, t] I + C TV gτ, + TV dτ, uτ L R + ;R dτ gτ L R +;R + dτ L R +;R TV uτ dτ uτ L R +;R + TV uτ dτ, completing the proof of.8. Concerning the stability bounds,.3 implies γt + γt ut, x 1 t dx bτ dτ + + ut, x dx uo x dx + γt + γt dτ, x uτ, x dx dτ dτ, x uτ, x dx dτ ut L1 u R + ;R o L 1 R + ;R + 1 b L 1 [,t];r + C uτ dτ. L1 R + ;R An application of Gronwall Lemma yields the desired estimate.11. Finally, the monotonicity property.1 directly follows from.3. To prove.9, fix t, t R + with t < t. Then, ut ut L1 R + ;R = γt γt + + γt + γt ut, x ut, x dx 4.7 ut, x ut, x dx 4.8 ut, x ut, x dx 4.9 and we deal with the three terms separately, using.3 as follows. Begin with 4.7: γt ut, x ut, x dx 14
15 γt γt + b T ; t, x g T ; t, x, exp T ;t,x b T ; t, x g T ; t, x, exp T ;t,x b T ; t, x g T ; t, x, b T ; t, x g T ; t, x, exp γt T ;t,x d τ, Xτ; t, x x g τ, Xτ; t, x dτ d τ, Xτ; t, x x g τ, Xτ; t, x dτ d τ, Xτ; t, x x g τ, Xτ; t, x dτ b T ; t, x g T ; t, x, exp d τ, Xτ; t, x x g τ, Xτ; t, x dτ T ;t,x dx dx exp T ;t,x d τ, Xτ; t, x x g τ, Xτ; t, x dτ ĝ t TVb; [, t ] TV g, ; [, t ] + b L [,t ];R t t C t e +ĝ t b L [,t ];R e C t t t C + C + ĝ t d L [,t ];R + TVd + x g, + TV L [,t ];R xg, L t t for a suitable positive constant L dependent on t,, ĝ, C and TVb; [, t ], b L [,t ];R. Passing to 4.8, use 4.6 and.5: γt γt ut, x ut, x { ut dx ĝ max, L R + ;R ut L t t R + ;R ĝ u o L R + ;R + 1 b L [,t ];R e Ct t t dx Finally, deal with 4.9 using 4.5: + γt + γt + γt + + γt ut, x ut, x dx ut, x u t, Xt ; t, x [ exp ut, x u t, Xt ; t, x dx u t, Xt ; t, x t exp t TVu gt L R + ;R t t + ut L 1 R + ;R t d τ, Xτ; t, x x g τ, Xτ; t, x dτ d τ, Xτ; t, x x g τ, Xτ; t, x dτ exp Ct t 1 ] dx 1 dx 15
16 ĝ TVu + C ut L1 t t R + ;R Completing the proof of.9. The following elementary lemma is of use below. Lemma 4.1. Let H, K R + and assume that the numbers B k R + satisfy B k+1 H + K B k for all k N. Then, B k K k B + 1 Kk 1 K H. Proof of Theorem.4. Fix a time T so that ] [ γt, min x i i=1,...,n 4.1 and define u t, x = u o x for t [, T ]. Recursively, for k 1 let u k u k 1,..., u k n solve t u k i + x g i t, x u k i = d i t, x u k i t t, x R+ R + g i t, u k i t, = bk i t t R+ i = 1,..., n 4.11 u i, x = u o i x x R+ where b k i t = α i t, u k 1 1 t, x 1,..., u k 1 n t, x n +β i u k 1 1 t, x dx,..., I 1 u k 1 n I n t, x dx. 4.1 Note that b is satisfied and Lemma. applies. Indeed, if k = 1, then b 1 i is independent on time. Let k > 1, then by.14 and.8 TV b k i ; [, T ] TV α i, u k 1 j, x j=1,...,n ; [, T ] + TV β i u k 1 j, x dx j=1,...,n ; [, T ] I j Lipα T + TV u k 1 j, x j=1,...,n ; [, T ] + Lipβ Lipα T + TV u k 1 j, x j j=1,...,n ; [, T ] +C Lipβ n i=1 [ u k 1 i τ + TV L R + ;R n TV u k 1 i, x dx ; [, T ] I i i=1 u k 1 i τ ] dτ Lipα T + TV u k 1 j, x j j=1,...,n ; [, T ] u +C Lipβ τ k 1 + TV u τ k 1 dτ L R + ;R n which is finite by induction. Lemma. then ensures existence and uniqueness of a solution to for any k >. By construction, the choice 4.1 ensures that u k t, x = u 1 t, x for all x > γt and k Therefore, also α i t, u k t, x i = α i t, u 1 t, x i for all t [, T ], for all k 1 and all i = 1,..., n. Compute now t u k i t u k+1 i L 1 R + ;R 16
17 + 1 b k+1 i b k i L1 [,t];r α i, u k, x i α i, u k 1, x i L 1 [,t];r i β I uk 1 1, x dx,..., I uk n n, x dx β i I 1 u k 1 1, x dx,..., I n u k 1 n, x dx n Lipβ u k j, x u k 1 j, x dx j=1 I j L 1 [,t];r n Lipβ u k j, x u k 1 j, x dx Lipβ Lipβ j=1 R + n j=1 Adding up all the components, u k j τ u k 1 τ dτ L 1 R + ;R j L 1 [,t];r L 1 [,t];r u k τ u k 1 τ dτ 4.14 L 1 R + ;R n u k+1 u k C [,T ];L 1 R + ;R n T n Lipβ u k u k 1 L1 dτ R + ;R n n Lipβ T u k u k 1 C [,T ];L 1 R + ;R n and, recursively, u k+1 u k C n Lipβ T k u 1 τ u τ [,T ];L 1 R + ;R n C [,T ];L 1 R + ;R n. Choosing now also T < 1/ n Lipβ, the sequence u k is a Cauchy sequence and we obtain the existence of a map u C [, T ]; L 1 R + ; R n which is the limit of the sequence u k, in the sense that lim sup u k t u t = L 1 R + ;R n k + t [,T ] To prove that u solves 1.1, it is sufficient to check that the boundary condition is attained. Indeed, proving that u is a weak solution to the balance law is a standard procedure. Clearly, the initial datum is attained, since u k = u o for all k. Using Lemma.,.3, 4.13 and 4.15, for all large k N we have b k+1 L [,T ];R n Lipα n u k, x i + Lipβ u k C L [,T ];R [,T ];L 1 R + ;R n i=1 Lipα u o L R + ;R n ect + Lipβ u C [,T ];L 1 R + ;R n With the above choice of T, using 4.14 and Lemma., TV b k+1 i ; [, T ] Lipα T + Lipα u o L [,maxj xj];rn + TV u o ; [, max x j ] e CT j T +C Lipβ u o L R + ;R n + 1 b k L e Cτ dτ [,τ];r n 17
18 [ T +C Lipβ u o L R + ;R n + TVu o + 1 [ C + b k L + TVb k ; [, τ]] ] e Cτ dτ [,τ];r n Lipα T + Lipα u o L R + ;R n + TVu o e CT + C Lipβ T u o L R + ;R n + TVu o e Cτ dτ + C T + C Lipβ b k L + TVb k ; [, τ] e Cτ dτ [,τ];r n Lipα T + Lipα u o L R + ;R n + TVu o e CT + Lipβ u o L R + ;R n + TVu o e CT 1 + C T + C Lipβ b k L + TVb k ; [, τ] e Cτ dτ [,τ];r n Lipα T + Lipα u o L R + ;R n + TVu o e CT + Lipβ u o L R + ;R n + TVu o e CT C b Lipβ k + TVb k ; [, T ] e CT 1 L [,T ];R n Inserting now the estimate 4.16 in the latter term above, we can apply Lemma 4.1 to the inequality B k+1 H + K B k, where B k = TV b k ; [, T ] H = n Lipα T + n Lipα e CT + Lipβe CT 1 u o L R + ;R n + TVu o +n + C Lipβ e CT 1 Lipα u o L R +;R n ect + Lipβ u C [,T ];L 1R +;R n + 1 e CT 1 K = n + C Lipβ e CT 1. Hence, as soon as T is so small that n + C Lipβ e CT 1 < 1, we obtain a bound on TV b k ; [, T ] uniform in k. This bound, thanks to Lemma., ensures that also sup t [,T ] sup k N TV u k t < + and, by the lower semicontinuity of the total variation with respect to the L 1 topology, also sup t [,T ] TV u t < +. Therefore, the trace lim x + u t, x exists for all t [, T ]. The uniform bound on TV b k ; [, T ], together with.9 and [6, Theorem.4], ensure that for a.e. x R +, we have lim k + u k t, x = u t, x. Choose one such x and observe that: u t, x = lim k + uk t, x = lim k + b k T ; t, x g T ; t, x, exp = lim k + b k T ; t, x g T ; t, x, exp T ;t,x T ;t,x dτ, Xτ; t, x x g τ, Xτ; t, x dτ dτ, Xτ; t, x x g τ, Xτ; t, x dτ 18
19 = lim k + B T ; t, x, u k 1 [ g T ; t, x, exp = B T ; t, x, u g T ; t, x, exp T ;t,x T ;t,x dτ, Xτ; t, x x g τ, Xτ; t, x ] dτ dτ, Xτ; t, x x g τ, Xτ; t, x dτ where in the last step above we used the convergences: lim α i t, u k j t, x j j=1,...,n = α i t, u t, x j j=1,...,n k + lim β i u k j t, x dx j=1,...,n = β i u t, x dx j=1,...,n k + I j I j by by The time T chosen above depends only on β, min i x i, on d and on g. In particular, it is independent from the initial datum. Hence, the above procedure can be iterated, extending u to a function defined on all R +, i.e. u C R + ; L 1 R + ; R n. Let now u o, u o be two initial data. Define x = min i=1,...,n x i and t = Γ x. Denote Ī = [ x, + [. To prove the stability estimate, with obvious notation, preliminarily compute for t [, t]: u t u t L 1 Ī;Rn u o u o L 1 R + ;R n ect u t u t L u Ī;Rn o u o L R + ;R n ect Moreover, by 4.1 and.14, for t [, t], b b L n Lipα u t u t + n Lipβ [, t];r n L Ī;Rn u t u t L1 R + ;R n n Lipα u o u o L R + ;R n ect n Lipβ u t u t L 1 R + ;R n. b b L1 n Lipα ect 1 [, t];r n u o u o C L R + ;R n +n Lipβ u τ u τ L dτ R + ;R n Therefore, by.9, for t [, t], u t u t L 1 R;R n 1 b b + L 1 [,t];r u o u o n Lipα e Ct 1 e Ct u o u o C + u o u o L1 R + ;R n ect +n Lipβ L 1 R + ;R n L R + ;R n e Ct u τ u τ L 1 R + ;R n dτ ect. An application of Gronwall Lemma yields, for t [, t], u t u t L 1 R;R n n Lipα e Ct 1 u o u o C L R + ;R n + u o u o exp Ct + n Lipβe Ct To iterate beyond time t, using.5, 4.17 and the above bound to estimate u t u t L R;R n L 1 R + ;R n 19
20 1 b b L [,t];r n + u o u o L e Ct R + ;R n n Lipα u o u o L R + ;R n ect + n Lipβ u t u t L1 R + ;R n ect + u o u o L R + ;R n ect n Lipα u e Ct + 1 o u o L R + ;R n + n Lipβ u t u t L 1 R + ;R n n Lipα e Ct n Lipβ n Lipα e Ct 1 exp Ct + n Lipβe Ct C u o u o L R + ;R n ect + n Lipβ exp Ct + n Lipβe Ct u o u o L1 R + ;R n completing the proof. Proof of Theorem.5. Define, for i = 1,..., n, b it = α i t, u 1 t, x 1,..., u nt, x n + β i u 1t, x dx,..., u nt, x dx, I 1 I n b i t = α i t, u 1t, x 1,..., u nt, x n + β i u 1t, x dx,..., u nt, x dx. I 1 I n e Ct Preliminary, using.14, let us estimate the term b b n L1 [,t];r n α i i=1 + n i=1 s, u 1 s, x 1,..., u ns, x n α i s, u 1s, x 1,..., u ns, x n ds β i I j u js, x dx β i t α α C R + R n ;R n + Lip α n i=1 u j s, x dx I j ds u js, x j u j s, x j ds +t β β + Lip C β u s u s R n ;R n L1 R + ;R n ds. Define x = min i=1,...,n x i and t = Γ x. As long as s [, T ], we have u j s, x j = u j s, x j. Hence the above estimate leads to b b L 1 [,t];r n t for all t [, t]. In the same time interval, α α C R+ Rn;Rn + β β CRn;Rn + Lip β u s u s L1 R + ;R n ds u t u t L1 ect b b R + ;R n L1 [,t];r n t ect α α C R+ Rn;Rn + β β CRn;Rn + Lipβ e Ct u s u s L1 R + ;R n ds
21 so that by Gronwall Lemma, for t [, t], u t u t [ t e Ct L1 α α R + ;R n C R+ Rn;Rn + ] β β CRn;Rn [ ] Lipβ t e Ct exp. A repeated application of the estimate above on the intervals [k 1 t, k t] allows to complete the proof. Proof of Corollary 3.1. Introduce u 1, u,... as in table 3.1. Then, extend d 1, d and g to R + R maintaining the required regularity and bounds on the total variation. The resulting system fits into Hence, d, g and.14 hold. Theorem.4 applies, ensuring the well posedness of the Cauchy problem. Finally, the solution to 3. is obtained restricting the solution to to [, a max ] and to [x min, x max ]. Acknowledgment: The authors were supported by the INDAM GNAMPA project Conservation Laws: Theory and Applications. References [1] A. S. Ackleh and K. Deng. A nonautonomous juvenile-adult model: well-posedness and long-time behavior via a comparison principle. SIAM J. Appl. Math., 696: , 9. [] A. S. Ackleh, K. Deng, and X. Yang. Sensitivity analysis for a structured juvenile adult model. Comput. Math. Appl., 643:19, 1. [3] L. Ambrosio, N. Fusco, and D. Pallara. Functions of bounded variation and free discontinuity problems. Oxford Mathematical Monographs. The Clarendon Press Oxford University Press, New York,. [4] C. Bardos, A. Y. le Roux, and J.-C. Nédélec. First order quasilinear equations with boundary conditions. Comm. Partial Differential Equations, 49: , [5] F. Brauer and C. Castillo-Chavez. Mathematical models in population biology and epidemiology, volume 4 of Texts in Applied Mathematics. Springer, New York, second edition, 1. [6] A. Bressan. Hyperbolic systems of conservation laws, volume of Oxford Lecture Series in Mathematics and its Applications. Oxford University Press, Oxford,. The one-dimensional Cauchy problem. [7] R. M. Colombo and G. Guerra. On general balance laws with boundary. J. Differential Equations, 485: , 1. [8] R. M. Colombo, G. Guerra, M. Herty, and V. Schleper. Optimal control in networks of pipes and canals. SIAM J. Control Optim., 483:3 5, 9. [9] M. Garavello and B. Piccoli. Traffic flow on networks, volume 1 of AIMS Series on Applied Mathematics. American Institute of Mathematical Sciences AIMS, Springfield, MO, 6. Conservation laws models. [1] N. Keyfitz. The mathematics of sex and marriage. In Proceedings of the Sixth Berkeley Symposium on Mathematical Statistics and Probability, Volume 4: Biology and Health, pages 89 18, Berkeley, Calif., 197. University of California Press. [11] S. N. Kružhkov. First order quasilinear equations with several independent variables. Mat. Sb. N.S., 81 13:8 55, 197. [1] R. J. LeVeque. Finite volume methods for hyperbolic problems. Cambridge Texts in Applied Mathematics. Cambridge University Press, Cambridge,. [13] P. Manfredi. Logistic effects in the two-sex model with harmonic mean fertility function. Genus, 491-:43 65, [14] B. Perthame. Transport equations in biology. Frontiers in Mathematics. Birkhäuser Verlag, Basel, 7. 1
22 [15] R. Schoen. The harmonic mean as the basis of a realistic two-sex marriage model. Demography, 18:1 16, cited By since [16] R. Schoen. Relationships in a simple harmonic mean two-sex fertility model. Journal of Mathematical Biology, 183:1 11, [17] R. Schoen. Modeling multigroup populations. Springer, [18] D. Serre. Systems of conservation laws.. Cambridge University Press, Cambridge,. Geometric structures, oscillations, and initial-boundary value problems, Translated from the 1996 French original by I. N. Sneddon. [19] A. Sundelof and P. Aberg. Birth functions in stage structured two-sex models. Ecological Modeling, 193: , 6.
L 1 Stability for scalar balance laws. Control of the continuity equation with a non-local flow.
L 1 Stability for scalar balance laws. Control of the continuity equation with a non-local flow. Magali Mercier Institut Camille Jordan, Lyon Beijing, 16th June 2010 Pedestrian traffic We consider tu +
More informationNONLINEAR DIFFERENTIAL INEQUALITY. 1. Introduction. In this paper the following nonlinear differential inequality
M athematical Inequalities & Applications [2407] First Galley Proofs NONLINEAR DIFFERENTIAL INEQUALITY N. S. HOANG AND A. G. RAMM Abstract. A nonlinear differential inequality is formulated in the paper.
More informationLocal semiconvexity of Kantorovich potentials on non-compact manifolds
Local semiconvexity of Kantorovich potentials on non-compact manifolds Alessio Figalli, Nicola Gigli Abstract We prove that any Kantorovich potential for the cost function c = d / on a Riemannian manifold
More informationOn uniqueness of weak solutions to transport equation with non-smooth velocity field
On uniqueness of weak solutions to transport equation with non-smooth velocity field Paolo Bonicatto Abstract Given a bounded, autonomous vector field b: R d R d, we study the uniqueness of bounded solutions
More informationBose-Einstein Condensation and Global Dynamics of Solutions to a Hyperbolic Kompaneets Equation
Bose-Einstein Condensation and Global Dynamics of Solutions to a Hyperbolic Kompaneets Equation Joshua Ballew Abstract In this article, a simplified, hyperbolic model of the non-linear, degenerate parabolic
More informationCoupling conditions for transport problems on networks governed by conservation laws
Coupling conditions for transport problems on networks governed by conservation laws Michael Herty IPAM, LA, April 2009 (RWTH 2009) Transport Eq s on Networks 1 / 41 Outline of the Talk Scope: Boundary
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationRegularity and compactness for the DiPerna Lions flow
Regularity and compactness for the DiPerna Lions flow Gianluca Crippa 1 and Camillo De Lellis 2 1 Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy. g.crippa@sns.it 2 Institut für Mathematik,
More informationOn a Lyapunov Functional Relating Shortening Curves and Viscous Conservation Laws
On a Lyapunov Functional Relating Shortening Curves and Viscous Conservation Laws Stefano Bianchini and Alberto Bressan S.I.S.S.A., Via Beirut 4, Trieste 34014 Italy. E-mail addresses: bianchin@mis.mpg.de,
More informationFunctional Differential Equations with Causal Operators
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.11(211) No.4,pp.499-55 Functional Differential Equations with Causal Operators Vasile Lupulescu Constantin Brancusi
More informationM. HERTY, CH. JÖRRES, AND B. PICCOLI
EXISTENCE OF SOLUTION TO SUPPLY CHAIN MODELS BASED ON PARTIAL DIFFERENTIAL EQUATION WITH DISCONTINUOUS FLUX FUNCTION M. HERTY, CH. JÖRRES, AND B. PICCOLI Abstract. We consider a recently [2] proposed model
More informationOn (h, k) trichotomy for skew-evolution semiflows in Banach spaces
Stud. Univ. Babeş-Bolyai Math. 56(2011), No. 4, 147 156 On (h, k) trichotomy for skew-evolution semiflows in Banach spaces Codruţa Stoica and Mihail Megan Abstract. In this paper we define the notion of
More informationA TWO PARAMETERS AMBROSETTI PRODI PROBLEM*
PORTUGALIAE MATHEMATICA Vol. 53 Fasc. 3 1996 A TWO PARAMETERS AMBROSETTI PRODI PROBLEM* C. De Coster** and P. Habets 1 Introduction The study of the Ambrosetti Prodi problem has started with the paper
More informationIdentification of Parameters in Neutral Functional Differential Equations with State-Dependent Delays
To appear in the proceedings of 44th IEEE Conference on Decision and Control and European Control Conference ECC 5, Seville, Spain. -5 December 5. Identification of Parameters in Neutral Functional Differential
More informationHyperbolic Systems of Conservation Laws
Hyperbolic Systems of Conservation Laws III - Uniqueness and continuous dependence and viscous approximations Alberto Bressan Mathematics Department, Penn State University http://www.math.psu.edu/bressan/
More informationNonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm
More informationLinear Ordinary Differential Equations
MTH.B402; Sect. 1 20180703) 2 Linear Ordinary Differential Equations Preliminaries: Matrix Norms. Denote by M n R) the set of n n matrix with real components, which can be identified the vector space R
More informationSTABILITY ESTIMATES FOR SCALAR CONSERVATION LAWS WITH MOVING FLUX CONSTRAINTS. Maria Laura Delle Monache. Paola Goatin
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X doi:10.3934/xx.xx.xx.xx pp. X XX STABILITY ESTIMATES FO SCALA CONSEVATION LAWS WITH MOVING FLUX CONSTAINTS Maria Laura Delle Monache Department
More informationL 1 stability of conservation laws for a traffic flow model
Electronic Journal of Differential Equations, Vol. 2001(2001), No. 14, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu ftp ejde.math.unt.edu (login:
More informationEXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM. Saeid Shokooh and Ghasem A. Afrouzi. 1. Introduction
MATEMATIČKI VESNIK MATEMATIQKI VESNIK 69 4 (217 271 28 December 217 research paper originalni nauqni rad EXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM Saeid Shokooh and Ghasem A.
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationLebesgue-Stieltjes measures and the play operator
Lebesgue-Stieltjes measures and the play operator Vincenzo Recupero Politecnico di Torino, Dipartimento di Matematica Corso Duca degli Abruzzi, 24, 10129 Torino - Italy E-mail: vincenzo.recupero@polito.it
More informationAN OVERVIEW ON SOME RESULTS CONCERNING THE TRANSPORT EQUATION AND ITS APPLICATIONS TO CONSERVATION LAWS
AN OVERVIEW ON SOME RESULTS CONCERNING THE TRANSPORT EQUATION AND ITS APPLICATIONS TO CONSERVATION LAWS GIANLUCA CRIPPA AND LAURA V. SPINOLO Abstract. We provide an informal overview on the theory of transport
More informationExistence Results for Multivalued Semilinear Functional Differential Equations
E extracta mathematicae Vol. 18, Núm. 1, 1 12 (23) Existence Results for Multivalued Semilinear Functional Differential Equations M. Benchohra, S.K. Ntouyas Department of Mathematics, University of Sidi
More informationA Lieb-Robinson Bound for the Toda System
1 FRG Workshop May 18, 2011 A Lieb-Robinson Bound for the Toda System Umar Islambekov The University of Arizona joint work with Robert Sims 2 Outline: 1. Motivation 2. Toda Lattice 3. Lieb-Robinson bounds
More informationExponential stability of families of linear delay systems
Exponential stability of families of linear delay systems F. Wirth Zentrum für Technomathematik Universität Bremen 28334 Bremen, Germany fabian@math.uni-bremen.de Keywords: Abstract Stability, delay systems,
More informationLimit solutions for control systems
Limit solutions for control systems M. Soledad Aronna Escola de Matemática Aplicada, FGV-Rio Oktobermat XV, 19 e 20 de Outubro de 2017, PUC-Rio 1 Introduction & motivation 2 Limit solutions 3 Commutative
More informationConvergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions
Convergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions Marie Hélène Vignal UMPA, E.N.S. Lyon 46 Allée d Italie 69364 Lyon, Cedex 07, France abstract. We are interested
More informationINITIAL AND BOUNDARY VALUE PROBLEMS FOR NONCONVEX VALUED MULTIVALUED FUNCTIONAL DIFFERENTIAL EQUATIONS
Applied Mathematics and Stochastic Analysis, 6:2 23, 9-2. Printed in the USA c 23 by North Atlantic Science Publishing Company INITIAL AND BOUNDARY VALUE PROBLEMS FOR NONCONVEX VALUED MULTIVALUED FUNCTIONAL
More informationLocal null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain
Local null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain Nicolás Carreño Université Pierre et Marie Curie-Paris 6 UMR 7598 Laboratoire
More informationVANISHING VISCOSITY IN THE PLANE FOR VORTICITY IN BORDERLINE SPACES OF BESOV TYPE
VANISHING VISCOSITY IN THE PLANE FOR VORTICITY IN BORDERLINE SPACES OF BESOV TYPE ELAINE COZZI AND JAMES P. KELLIHER Abstract. The existence and uniqueness of solutions to the Euler equations for initial
More informationThe 2D Magnetohydrodynamic Equations with Partial Dissipation. Oklahoma State University
The 2D Magnetohydrodynamic Equations with Partial Dissipation Jiahong Wu Oklahoma State University IPAM Workshop Mathematical Analysis of Turbulence IPAM, UCLA, September 29-October 3, 2014 1 / 112 Outline
More informationControl of the Continuity Equation with a Non Local Flow
Control of the Continuity Equation with a Non Local Flow Rinaldo M. Colombo Michael Herty Magali Mercier. September 29, 29 Abstract This paper focuses on the analytical properties of the solutions to the
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationRégularité des équations de Hamilton-Jacobi du premier ordre et applications aux jeux à champ moyen
Régularité des équations de Hamilton-Jacobi du premier ordre et applications aux jeux à champ moyen Daniela Tonon en collaboration avec P. Cardaliaguet et A. Porretta CEREMADE, Université Paris-Dauphine,
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationNonlinear Control Lecture 5: Stability Analysis II
Nonlinear Control Lecture 5: Stability Analysis II Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 Farzaneh Abdollahi Nonlinear Control Lecture 5 1/41
More informationIntersection Models and Nash Equilibria for Traffic Flow on Networks
Intersection Models and Nash Equilibria for Traffic Flow on Networks Alberto Bressan Department of Mathematics, Penn State University bressan@math.psu.edu (Los Angeles, November 2015) Alberto Bressan (Penn
More informationA NOTE ON ALMOST PERIODIC VARIATIONAL EQUATIONS
A NOTE ON ALMOST PERIODIC VARIATIONAL EQUATIONS PETER GIESL AND MARTIN RASMUSSEN Abstract. The variational equation of a nonautonomous differential equation ẋ = F t, x) along a solution µ is given by ẋ
More informationAN ALGORITHM FOR DETERMINATION OF AGE-SPECIFIC FERTILITY RATE FROM INITIAL AGE STRUCTURE AND TOTAL POPULATION
J Syst Sci Complex (212) 25: 833 844 AN ALGORITHM FOR DETERMINATION OF AGE-SPECIFIC FERTILITY RATE FROM INITIAL AGE STRUCTURE AND TOTAL POPULATION Zhixue ZHAO Baozhu GUO DOI: 117/s11424-12-139-8 Received:
More informationPeriodic solutions of weakly coupled superlinear systems
Periodic solutions of weakly coupled superlinear systems Alessandro Fonda and Andrea Sfecci Abstract By the use of a higher dimensional version of the Poincaré Birkhoff theorem, we are able to generalize
More informationMildly degenerate Kirchhoff equations with weak dissipation: global existence and time decay
arxiv:93.273v [math.ap] 6 Mar 29 Mildly degenerate Kirchhoff equations with weak dissipation: global existence and time decay Marina Ghisi Università degli Studi di Pisa Dipartimento di Matematica Leonida
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 5. Input-Output Stability DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Input-Output Stability y = Hu H denotes
More informationLeast Squares Based Self-Tuning Control Systems: Supplementary Notes
Least Squares Based Self-Tuning Control Systems: Supplementary Notes S. Garatti Dip. di Elettronica ed Informazione Politecnico di Milano, piazza L. da Vinci 32, 2133, Milan, Italy. Email: simone.garatti@polimi.it
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationarxiv: v1 [math.ap] 5 Nov 2018
STRONG CONTINUITY FOR THE 2D EULER EQUATIONS GIANLUCA CRIPPA, ELIZAVETA SEMENOVA, AND STEFANO SPIRITO arxiv:1811.01553v1 [math.ap] 5 Nov 2018 Abstract. We prove two results of strong continuity with respect
More informationConservation laws and some applications to traffic flows
Conservation laws and some applications to traffic flows Khai T. Nguyen Department of Mathematics, Penn State University ktn2@psu.edu 46th Annual John H. Barrett Memorial Lectures May 16 18, 2016 Khai
More informationPERIODIC PROBLEMS WITH φ-laplacian INVOLVING NON-ORDERED LOWER AND UPPER FUNCTIONS
Fixed Point Theory, Volume 6, No. 1, 25, 99-112 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.htm PERIODIC PROBLEMS WITH φ-laplacian INVOLVING NON-ORDERED LOWER AND UPPER FUNCTIONS IRENA RACHŮNKOVÁ1 AND MILAN
More informationNONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS
Fixed Point Theory, Volume 9, No. 1, 28, 3-16 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.html NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS GIOVANNI ANELLO Department of Mathematics University
More informationON THE FOLIATION OF SPACE-TIME BY CONSTANT MEAN CURVATURE HYPERSURFACES
ON THE FOLIATION OF SPACE-TIME BY CONSTANT MEAN CURVATURE HYPERSURFACES CLAUS GERHARDT Abstract. We prove that the mean curvature τ of the slices given by a constant mean curvature foliation can be used
More informationAW -Convergence and Well-Posedness of Non Convex Functions
Journal of Convex Analysis Volume 10 (2003), No. 2, 351 364 AW -Convergence Well-Posedness of Non Convex Functions Silvia Villa DIMA, Università di Genova, Via Dodecaneso 35, 16146 Genova, Italy villa@dima.unige.it
More informationOSCILLATION AND ASYMPTOTIC STABILITY OF A DELAY DIFFERENTIAL EQUATION WITH RICHARD S NONLINEARITY
2004 Conference on Diff. Eqns. and Appl. in Math. Biology, Nanaimo, BC, Canada. Electronic Journal of Differential Equations, Conference 12, 2005, pp. 21 27. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
More informationExistence and uniqueness of solutions for nonlinear ODEs
Chapter 4 Existence and uniqueness of solutions for nonlinear ODEs In this chapter we consider the existence and uniqueness of solutions for the initial value problem for general nonlinear ODEs. Recall
More informationInstitut für Mathematik
U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Martin Rasmussen, Peter Giesl A Note on Almost Periodic Variational Equations Preprint Nr. 13/28 14. März 28 Institut für Mathematik, Universitätsstraße,
More informationUNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE
Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (2010), 275 284 UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Iuliana Carmen Bărbăcioru Abstract.
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 3. Fundamental properties IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Example Consider the system ẋ = f
More informationDifferentiability with respect to initial data for a scalar conservation law
Differentiability with respect to initial data for a scalar conservation law François BOUCHUT François JAMES Abstract We linearize a scalar conservation law around an entropy initial datum. The resulting
More informationNON-EXTINCTION OF SOLUTIONS TO A FAST DIFFUSION SYSTEM WITH NONLOCAL SOURCES
Electronic Journal of Differential Equations, Vol. 2016 (2016, No. 45, pp. 1 5. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu NON-EXTINCTION OF
More informationHow large is the class of operator equations solvable by a DSM Newton-type method?
This is the author s final, peer-reviewed manuscript as accepted for publication. The publisher-formatted version may be available through the publisher s web site or your institution s library. How large
More informationExtremal Solutions of Differential Inclusions via Baire Category: a Dual Approach
Extremal Solutions of Differential Inclusions via Baire Category: a Dual Approach Alberto Bressan Department of Mathematics, Penn State University University Park, Pa 1682, USA e-mail: bressan@mathpsuedu
More informationGIOVANNI COMI AND MONICA TORRES
ONE-SIDED APPROXIMATION OF SETS OF FINITE PERIMETER GIOVANNI COMI AND MONICA TORRES Abstract. In this note we present a new proof of a one-sided approximation of sets of finite perimeter introduced in
More informationStability Theory for Nonnegative and Compartmental Dynamical Systems with Time Delay
1 Stability Theory for Nonnegative and Compartmental Dynamical Systems with Time Delay Wassim M. Haddad and VijaySekhar Chellaboina School of Aerospace Engineering, Georgia Institute of Technology, Atlanta,
More informationExistence and uniqueness of the weak solution for a contact problem
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 (216), 186 199 Research Article Existence and uniqueness of the weak solution for a contact problem Amar Megrous a, Ammar Derbazi b, Mohamed
More informationNavier Stokes and Euler equations: Cauchy problem and controllability
Navier Stokes and Euler equations: Cauchy problem and controllability ARMEN SHIRIKYAN CNRS UMR 888, Department of Mathematics University of Cergy Pontoise, Site Saint-Martin 2 avenue Adolphe Chauvin 9532
More informationNOTES ON CALCULUS OF VARIATIONS. September 13, 2012
NOTES ON CALCULUS OF VARIATIONS JON JOHNSEN September 13, 212 1. The basic problem In Calculus of Variations one is given a fixed C 2 -function F (t, x, u), where F is defined for t [, t 1 ] and x, u R,
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationDIFFERENTIAL CRITERIA FOR POSITIVE DEFINITENESS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX)- DIFFERENTIAL CRITERIA FOR POSITIVE DEFINITENESS J. A. PALMER Abstract. We show how the Mellin transform can be used to
More informationInternal Stabilizability of Some Diffusive Models
Journal of Mathematical Analysis and Applications 265, 91 12 (22) doi:1.16/jmaa.21.7694, available online at http://www.idealibrary.com on Internal Stabilizability of Some Diffusive Models Bedr Eddine
More informationLinear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems
Linear ODEs p. 1 Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems Linear ODEs Definition (Linear ODE) A linear ODE is a differential equation taking the form
More informationEntropic Schemes for Conservation Laws
CONSTRUCTVE FUNCTON THEORY, Varna 2002 (B. Bojanov, Ed.), DARBA, Sofia, 2002, pp. 1-6. Entropic Schemes for Conservation Laws Bojan Popov A new class of Godunov-type numerical methods (called here entropic)
More informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
More informationSingularity formation for compressible Euler equations
Singularity formation for compressible Euler equations Geng Chen Ronghua Pan Shengguo Zhu Abstract In this paper, for the p-system and full compressible Euler equations in one space dimension, we provide
More information2 Statement of the problem and assumptions
Mathematical Notes, 25, vol. 78, no. 4, pp. 466 48. Existence Theorem for Optimal Control Problems on an Infinite Time Interval A.V. Dmitruk and N.V. Kuz kina We consider an optimal control problem on
More informationExponential Decay: From Semi-Global to Global
Exponential Decay: From Semi-Global to Global Martin Gugat Benasque 2013 Martin Gugat (FAU) From semi-global to global. 1 / 13 How to get solutions that are global in time? 1 Sometimes, the semiglobal
More informationContractive metrics for scalar conservation laws
Contractive metrics for scalar conservation laws François Bolley 1, Yann Brenier 2, Grégoire Loeper 34 Abstract We consider nondecreasing entropy solutions to 1-d scalar conservation laws and show that
More informationEXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM
1 EXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM V. BERTI and M. FABRIZIO Dipartimento di Matematica, Università degli Studi di Bologna, P.zza di Porta S. Donato 5, I-4126, Bologna,
More information1.4 The Jacobian of a map
1.4 The Jacobian of a map Derivative of a differentiable map Let F : M n N m be a differentiable map between two C 1 manifolds. Given a point p M we define the derivative of F at p by df p df (p) : T p
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationFractal Conservation Laws: Global Smooth Solutions and Vanishing Regularization
Progress in Nonlinear Differential Equations and Their Applications, Vol. 63, 217 224 c 2005 Birkhäuser Verlag Basel/Switzerland Fractal Conservation Laws: Global Smooth Solutions and Vanishing Regularization
More informationOn the dynamics of strongly tridiagonal competitive-cooperative system
On the dynamics of strongly tridiagonal competitive-cooperative system Chun Fang University of Helsinki International Conference on Advances on Fractals and Related Topics, Hong Kong December 14, 2012
More informationConvergence Rate of Nonlinear Switched Systems
Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the
More informationVariational and Topological methods : Theory, Applications, Numerical Simulations, and Open Problems 6-9 June 2012, Northern Arizona University
Variational and Topological methods : Theory, Applications, Numerical Simulations, and Open Problems 6-9 June 22, Northern Arizona University Some methods using monotonicity for solving quasilinear parabolic
More informationA NOTE ON THE INITIAL-BOUNDARY VALUE PROBLEM FOR CONTINUITY EQUATIONS WITH ROUGH COEFFICIENTS. Gianluca Crippa. Carlotta Donadello. Laura V.
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume, Number, Xxxx XXXX pp. A NOTE ON THE INITIAL-BOUNDARY VALUE PROBLEM FOR CONTINUITY EQUATIONS WITH ROUGH COEFFICIENTS Gianluca
More informationThe Bang-Bang theorem via Baire category. A Dual Approach
The Bang-Bang theorem via Baire category A Dual Approach Alberto Bressan Marco Mazzola, and Khai T Nguyen (*) Department of Mathematics, Penn State University (**) Université Pierre et Marie Curie, Paris
More informationOn a weighted total variation minimization problem
On a weighted total variation minimization problem Guillaume Carlier CEREMADE Université Paris Dauphine carlier@ceremade.dauphine.fr Myriam Comte Laboratoire Jacques-Louis Lions, Université Pierre et Marie
More informationEXISTENCE RESULTS FOR NONLINEAR FUNCTIONAL INTEGRAL EQUATIONS VIA NONLINEAR ALTERNATIVE OF LERAY-SCHAUDER TYPE
ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA IAŞI Tomul LII, s.i, Matematică, 26, f.1 EXISTENCE RESULTS FOR NONLINEAR FUNCTIONAL INTEGRAL EQUATIONS VIA NONLINEAR ALTERNATIVE OF LERAY-SCHAUDER TYPE
More informationON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR TERM
Internat. J. Math. & Math. Sci. Vol. 22, No. 3 (999 587 595 S 6-72 9922587-2 Electronic Publishing House ON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR
More informationOn the Weak Solutions of the McKendrick Equation: Existence of Demography Cycles
Math. Model. Nat. Phenom. Vol. 1, No. 1, 26, pp. 1-3 On the Weak Solutions of the McKendrick Equation: Existence of Demography Cycles R. Dilão 1 and A. Lakmeche Nonlinear Dynamics Group, Instituto Superior
More informationStability of solutions to abstract evolution equations with delay
Stability of solutions to abstract evolution equations with delay A.G. Ramm Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA ramm@math.ksu.edu Abstract An equation u = A(t)u+B(t)F
More informationExistence And Uniqueness Of Mild Solutions Of Second Order Volterra Integrodifferential Equations With Nonlocal Conditions
Applied Mathematics E-Notes, 9(29), 11-18 c ISSN 167-251 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Existence And Uniqueness Of Mild Solutions Of Second Order Volterra Integrodifferential
More informationA model for a network of conveyor belts with discontinuous speed and capacity
A model for a network of conveyor belts with discontinuous speed and capacity Adriano FESTA Seminario di Modellistica differenziale Numerica - 6.03.2018 work in collaboration with M. Pfirsching, S. Goettlich
More informationOn a Class of Multidimensional Optimal Transportation Problems
Journal of Convex Analysis Volume 10 (2003), No. 2, 517 529 On a Class of Multidimensional Optimal Transportation Problems G. Carlier Université Bordeaux 1, MAB, UMR CNRS 5466, France and Université Bordeaux
More informationObserver design for a general class of triangular systems
1st International Symposium on Mathematical Theory of Networks and Systems July 7-11, 014. Observer design for a general class of triangular systems Dimitris Boskos 1 John Tsinias Abstract The paper deals
More informationGLOBAL EXISTENCE FOR THE ONE-DIMENSIONAL SEMILINEAR TRICOMI-TYPE EQUATIONS
GLOBAL EXISTENCE FOR THE ONE-DIMENSIONAL SEMILINEAR TRICOMI-TYPE EQUATIONS ANAHIT GALSTYAN The Tricomi equation u tt tu xx = 0 is a linear partial differential operator of mixed type. (For t > 0, the Tricomi
More informationHysteresis rarefaction in the Riemann problem
Hysteresis rarefaction in the Riemann problem Pavel Krejčí 1 Institute of Mathematics, Czech Academy of Sciences, Žitná 25, 11567 Praha 1, Czech Republic E-mail: krejci@math.cas.cz Abstract. We consider
More informationExistence and Uniqueness
Chapter 3 Existence and Uniqueness An intellect which at a certain moment would know all forces that set nature in motion, and all positions of all items of which nature is composed, if this intellect
More informationOn the Front-Tracking Algorithm
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 7, 395404 998 ARTICLE NO. AY97575 On the Front-Tracking Algorithm Paolo Baiti S.I.S.S.A., Via Beirut 4, Trieste 3404, Italy and Helge Kristian Jenssen
More informationExistence and Uniqueness of the Weak Solution for a Contact Problem
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. x (215), 1 15 Research Article Existence and Uniqueness of the Weak Solution for a Contact Problem Amar Megrous a, Ammar Derbazi b, Mohamed Dalah
More information