arxiv: v1 [math.ap] 4 Sep 2017

Transcription

1 RELATIVE ENTROPY METHOD FOR MEASURE-VALUED SOLUTIONS IN NATURAL SCIENCES TOMASZ DEBIEC, PIOTR GWIAZDA, KAMILA ŁYCZEK, AND AGNIESZKA ŚWIERCZEWSKA-GWIAZDA. arxiv:79.4v [math.ap] 4 Sep 27 ABSTRACT. We describe the applications of the relative entropy framework introduced in []. In particular uniqueness of an entropy solution is proven for a scalar conservation law, using the notion of measure-valued entropy solutions. Further we survey recent results concerning measurevalued-strong uniqueness for a number of physical systems incompressible and compressible Euler equations, compressible Navier-Stokes, polyconvex elastodynamics and general hyperbolic conservation laws, as well as long-time asymptotics of the McKendrick-Von Foerster equation.. INTRODUCTION The origins of the relative entropy method can be traced back to physics. The underlying principle behind it is the simple idea to measure in a certain way how much two evolutions of a given physical system, whose initial states are "close", differ and to investigate how this "distance" evolves in time. This framework, closely related to the second law of thermodynamics, is a useful tool in obtaining a variety of interesting analytical results. For instance it can be used to show uniqueness of solutions to a conservation law in the scalar case, while for many systems of equations it provides the so-called weak-strong uniqueness property, i.e. establishes uniqueness of classical solutions in a wider class of weak solutions. This application, first described by Dafermos [, ], will be highlighted in this article. Other areas where relative entropy method is found useful include stability studies, asymptotic limits and dimension reduction problems (e.g. [9], [23], [7], [2]). The method is also applied to problems arising from biology, cf. [32], [33], [37], [26], known in this context as General Relative Entropy (GRE). It is essentially used for showing asymptotic convergence of solutions to steadystate solutions. On the level of weak solutions for various physical systems, including Navier-Stokes and Euler, the story seems quite complete. However recent years have delivered many new results on the level of measure-valued solutions (e.g. [4], [3], [27], [8]). This shows that even though mv solutions are considered a very weak notion of solution, not carrying much information, they do play an important role in the analysis of physical systems. We begin our discussion on the level of the scalar conservation law t u(x,t)+div x f(u(x,t)) =, u(x,) = u (x), int d R +, int d. (.) in the framework of measure-valued solutions. Here R + =[,+ ), T d =(R/2πZ) d and u is a given initial datum. The main ideas come from Tartar [42] and DiPerna [5], who defined entropy mv solutions in the language of classical Young measures, see also [35]. Definition. (Measure-valued solution). A measurable measure-valued map ν :(x,t) ν (x,t) Prob(R)

2 2 TOMASZ DEBIEC, PIOTR GWIAZDA, KAMILA ŁYCZEK, AND AGNIESZKA ŚWIERCZEWSKA-GWIAZDA. from T d R + to the space of probability measures on R is a measure-valued solution of (.) if t ν (x,t),λ +div x ν (x,t), f(λ) = (.2) in the sense of distributions, that is R + for all ϕ C c (Td R + ). T d { ν(x,t),λ t ϕ+ ν (x,t), f(λ) x ϕ } dx dt = Here, denotes the duality pairing between C (R), the closure with respect to the supremum norm of the space of continuous functions on R with compact support, and M(R), the space of signed Radon measures onr, i.e. µ,g(λ) := g(λ) dµ(λ). R By measurability of a measure ν we mean the weak -measurability of the measures ν (x,t), i.e. measurability of the map (x,t) ν (x,t),g(λ) for each g C (R). Such a measure-valued solution often arises from a weakly convergent approximating sequence. The framework of Young measures is in a sense a way of immersing the initial problem into a wider space in this way one gains linearity at the cost of having to deal with measure spaces rather than function spaces. In other words Young measures allow to deal with the non-commutativity of weak limits with nonlinearities. Indeed, it can be shown that if {u k } is a sequence uniformly bounded in L p, then, along a non-relabelled subsequence, the weak limit of { f(u k )} can be represented by means of a parameterised family of measures, a Young measure. Lemma.2 (Fundamental Lemma of Classical Young Measures). be a sequence uniformly bounded in L (R + ;L p (T d )), i.e. Let {u k } : T d R + R sup u k (,t) L p (T d ) C, for almost every t R +. (.3) k Then there exists a subsequence, still denoted {u k } and a measurable measure-valued mapping ν :T d R + Prob(R), such that for each g C(R d ), satisfying the growth condition g(λ) C(+ λ q ) for q< p, the weak limit of g(u k (x,t)) exists and is represented by ν (x,t),g(λ), i.e. lim g(u k (x,t))ϕ(x,t) dx dt= ν (x,t),g(λ) ϕ(x,t) dx dt, ϕ L (T d R + ). k T d R + T d R + (.4) DiPerna [5] and later Szepessy [4] under different assumptions on the continuous flux showed the so-called averaged contraction principle, which is a crucial estimate in showing uniqueness of mv solutions, and is essentially a form of relative entropy inequality. Importantly the uniqueness result is proven under the assumption that the initial data is a Dirac delta measure. For non-atomic initial data uniqueness might fail in the class of measure-valued solutions, even in the scalar case and even under an entropy inequality, as it provides information only on certain moments of the solution, cf. [2] To show existence of a measure-valued entropy solution to (.) a parabolic approximate problem is considered. This generates a sequence of approximate solutions, which can be shown to be uniformly integrable. Thus one can see that there is going to be no concentration effect. However,

3 RELATIVE ENTROPY METHOD FOR MEASURE-VALUED SOLUTIONS IN NATURAL SCIENCES 3 other approximation schemes can be considered, which will not posses sufficient integrability. In fact we prove uniqueness in a wider class of mv solutions not necessarily corresponding to any approximation scheme. Further, the situation differs substantially in the case of hyperbolic systems. The result introducing measure-valued solutions for the incompressible Euler describes both oscillations and concentrations, cf [6], []. This is because, contrary to the scalar case, for systems there is usually only one entropy-entropy flux pair forming a companion law. The corresponding entropy inequality lacks the symmetry which is present in the scalar case. One then considers the so-called generalized Young measure. Here by a generalized Young measure we mean a triple (ν,m,ν ) describing oscillations, concentrations and concentrationangle respectively. The following representation result, proven in [] is then true Lemma.3. Let {u k } :T d R + R n be a sequence bounded in L (R + ;L (T d )), i.e. sup k u k (,t) L (T d ) C, for almost every t R +. (.5) Then there exists a subsequence, still denoted {u k }, a measurable measure-valued mapping ν : T d R + Prob(R), a non-negative measure m on T d R + and a parameterised probability measure ν L w (Td R +,m,prob(s n )) such that for each Caratheodory function g : T d R + R n R n g(x,t,u k (x,t)) ν (x,t),g(x,t,λ) + ν (x,t),g m (.6) weak in the sense of measures, provided the function g, defined to be g (x,t,β) := lim s is of class C(T d R + S n ). lim (x,t,β ) (x,t,β) g(x,t,sβ ), where β S n. s In this article we want to show in as technically simple a way as possible the main idea of this framework. To this end we provide in Section 2 the full proof of an averaged contraction principle for the scalar conservation law under the artificial assumption that concentration effects cannot be excluded. This highlights an interesting phenomenon that the concentration measure indeed vanishes thus concluding uniqueness on the basis of averaged contraction principle. The same observation transfers to systems, of course on the level of proving only weak(measure-valued)-strong uniqueness, cf. [3], [27]. Interestingly, even more can be claimed, namely that the information about concentration angle does not give any essential information, cf. [8]. Indeed it is enough to know that the concentration measure appearing in the weak formulation is dominated by the concentration measure coming from the energy/entropy inequality. Then in Sections 3-5 we survey the weak-strong uniqueness results obtained using the relative entropy method for equations of fluid dynamics (Section 3), polyconvex elastodynamics (Section 4) and general systems of conservation laws. Finally we display how the relative entropy method is used for the renewal equation of mathematical biology (Section 6). 2. MEASURE-VALUED ENTROPY SOLUTIONS FOR SCALAR CONSERVATION LAWS We consider the Cauchy problem for the scalar conservation law t u+div x f(u)=, u(x,)=u (x), (2.)

4 4 TOMASZ DEBIEC, PIOTR GWIAZDA, KAMILA ŁYCZEK, AND AGNIESZKA ŚWIERCZEWSKA-GWIAZDA. where f :R R d represents the flux of the quantity u :T d R + R. We choose to work in the spatially periodic setting to avoid inessential technical issues. We introduce now the concept of entropy measure-valued solutions, which will then be used to show existence of a unique entropy solution to the problem (2.). The following will be standing assumptions on the flux function f and the initial datum u : f C(R) and u L (T d ), (2.2) f(λ) C(+ λ ), for some C>. (2.3) The definition of a measure-valued entropy solution to (2.) consists of a classical Young measure ν as well as two concentration measures m and m 2. We will assume that m 2 (T d A) C m (T d A) (2.4) for any borel set A R +. We will later see that for a solution arising as a limit of an approximating sequence this is guaranteed by characterization of the corresponding concentration measures. Definition 2. (Measure-valued entropy solution). The triple (ν,m,m 2 ) generated by a sequence which satisfies (.5) is called a measure-valued entropy solution with concentration of conservation law (2.) if t ( ν, λ k +m )+div x ( ν,sgn(λ k)( f(λ) f(k)) +m 2 ) (2.5) in the sense of distributions ont d R + for all k R, and if { T lim ν T + T (x,t), λ u (x) dx dt+ T d T T d [,T] } m (dxdt) =. (2.6) Notice that the entropy inequality (2.5) is sufficient to guarantee that the triple (ν,m,m 2 ) satisfies also the weak formulation. To see this assume, for simplicity, that the measure ν has a bounded support supp ν (x,t) ( A,A). Then taking k< A we have λ k> and (2.5) becomes for any ϕ ν (x,t),λ k t ϕ dx dt+ t ϕ m (dxdt) T d R + T d R + + ν (x,t), f(λ) f(k) ϕ dx dt+ ϕ m 2 (dxdt) T d R + T d R + = ν (x,t),λ t ϕ dx dt+ t ϕ m (dxdt) T d R + T d R + + ν (x,t), f(λ) ϕ dx dt+ ϕ m 2 (dxdt). T d R + T d R + Taking in turn k>a one obtains the reversed inequality. It follows that t ( ν (x,t),λ +m )+div x ( ν (x,t), f(λ) +m 2 )=. Now we prove a uniqueness result concerning mv entropy solutions. Theorem 2.2. Suppose (ν,m,m 2 ) and (σ,m,m 2 ) are Young measure-concentration measure triples satisfying (2.5) and ν (x,) = σ (x,) = δ {u (x)}. Then there exists a function such that w L (R + ;L (T d )) (2.7) ν (x,t) = σ (x,t) = δ {w(x,t)} (2.8) for almost each (x,t) T d R + and m i = m i =, i=,2.

5 RELATIVE ENTROPY METHOD FOR MEASURE-VALUED SOLUTIONS IN NATURAL SCIENCES 5 Proof. We begin by mollifying the measures ν, σ, m i and m i, i =,2. Let η be a non-negative, symmetric smooth function ont d R + compactly supported in the open unit ball inr d R + and such that R d R + η(x,t)=, and let η ε (x,t)=ε (d+) η( x ε, t ε ). By νε we denote the parameterised measure satisfying for all g C(R) ν ε (x,t),g = T d R + η ε (x x,t t ) ν (x,t ),g dx dt C (K;M(R)) where K T d R + is a set whose ε-neighbourhood is entirely contained int d R +. Furthermore we observe that by virtue of the Riesz representation theorem, for each (x,t) T d R + there are bounded measures x ν(x,t) ε and tν(x,t) ε s.t. for all g C(R) α ν(x,t) ε,g = α ν(x,t) ε,g, for α {t,x}, and the map (x,t) α ν(x,t) ε,g is continuous. Observe that mollifying measures m i and m i with a regular kernel gives smooth functions m ε,mε C (K;R + ) and m2 ε,mε 2 C (K;R). Observe also that ν and σ have finite first moments. We will now show that the regularized measures satisfy the entropy inequality. Let V T d R + be an open bounded set. For sufficiently small ε > we have, for each non-negative ϕ C c (V) and g C(R), ( ν(x,t),g ) t (ϕ η ε ) dx dt+ t (ϕ η ε ) m (dxdt) T d R + T d R + ( ) = t ν(x,t) Td R+ ε,g + tm ε ϕ(x, t) dx dt, and analogously for the spatial derivative. Hence, choosing in turn g(λ) = λ µ and g(λ) = sgn(λ µ)( f(λ) f(µ)) and using (2.5) ) ( t ν + Td R+ ε(x,t), λ µ t m ε + x ν ε(x,t),sgn(λ µ)( f(λ) f(µ)) + x m2 ε ϕ(x,t) dx dt [ = ν(x,t), λ µ t (ϕ η ε )+ ν (x,t),q(λ,µ) (ϕ η ε ) ] dx dt T d R + t (ϕ η ε ) m (dxdt) (ϕ η ε ) m 2 (dxdt) T d R + T d R +, where we denote q(λ,µ) := sgn(λ µ)( f(λ) f(µ)) for brevity. Therefore (2.9) t ν ε (x,t), λ µ + tm ε + x ν ε (x,t),sgn(λ µ)( f(λ) f(µ)) + xm ε 2 (2.) for all µ R and(x,t) V. Symmetrically it can be seen that the triple(σ ε,m ε,mε 2 ) satisfies (2.5) for all λ R. Next we observe that the function q is continuous on R 2 and it has sublinear growth due to growth conditions (2.3) on f. It follows that the maps µ q(λ,µ) dν(x,t) ε (λ) and λ q(λ,µ) dσ(x,t) ε (µ) are continuous (by virtue of the dominated convergence theorem). Let now ε correspond to the smoothing of(ν,m,m 2 ) and ε 2 correspond to smoothing of(σ,m,m 2 ). We can then compute ) div x ( ν ε (x,t) σ ε 2 (x,t),q(λ,µ) + m ε 2 + mε 2 2 ) = div x( q(λ,µ) dν ε (x,t) dσ ε 2 (x,t) + q(λ,µ) dν ε (x,t) d( xσ ε 2 (x,t) )+ xm ε 2 + xm ε 2 2 R R R R = x ν ε (x,t),q(λ,µ) dσ ε 2 (x,t) + x σ ε 2 (x,t),q(λ,µ) dνε (x,t) + xm ε 2 + xm ε 2 2 R R

6 6 TOMASZ DEBIEC, PIOTR GWIAZDA, KAMILA ŁYCZEK, AND AGNIESZKA ŚWIERCZEWSKA-GWIAZDA. (2.) where we have used Fubini s theorem in the last line. The tensor product ν (x,t) σ (x,t) is defined as the product measure dν (x,t) (λ)dσ (x,t) (µ). Similarly ) t ( ν ε (x,t) σ ε 2 (x,t), λ µ + m ε + mε 2 (2.2) = t ν ε (x,t), λ µ dσ ε 2 (x,t) + t σ ε 2 (x,t), λ µ dνε (x,t) + tm ε + tm ε 2. R Consequently using (2.) ( ) ν Td R+ ε (x,t) σ ε 2 (x,t), λ µ +mε + mε 2 t ϕ dx dt ( ) + ν Td R+ ε (x,t) σ ε 2 (x,t),q(λ,µ) +mε 2 + mε 2 2 ϕ dx dt {( ) = t ν ε (x,t), λ µ + xν ε (x,t),q(λ,µ). ϕ T d R + R ϕ T d R + R R } dσ ε 2 (x,t) tm ε xm ε 2 dx dt } {( ) t σ ε 2 (x,t), λ µ + xσ ε 2 (x,t),q(λ,µ) dν ε (x,t) tm ε 2 xm ε 2 2 (2.3) This holds for any ϕ C c(v), thus establishing the following inequality in D (T d R + ) ( ) t ( ν ε (x,t) σ ε 2 (x,t), λ µ + m ε + mε 2 )+div x ν ε (x,t) σ ε 2 (x,t),q(λ,µ) + m ε 2 + mε 2 2. (2.4) Observe that the function ϕ(x,t) = ψ(x)ϕ(t) with ψ and ϕ C c (R +) is an admissible test function, since T d is a compact manifold without boundary. Testing (2.4) with such a function yields ( ) ν Td R+ ε (x,t) σ ε 2 (x,t), λ µ +mε + mε 2 t ϕ dx dt. (2.5) The argument presented in [4] establishes the limit lim lim ε ε 2 Td R+ ν ε (x,t) σ ε 2 (x,t), λ µ t ϕ dx dt= T d R + ν(x,t) σ (x,t), λ µ dx dt. Clearly we have the convergence m ε T d tϕ dx dt = t (ϕ η ε ) m (dxdt) ε t ϕ m (dxdt). R + T d R + T d R + Similarly m ε 2 tϕ t ϕ m (dxdt). We thus have for any non-negative ϕ C c (R +) ν(x,t) σ (x,t), λ µ t ϕ(t) dx dt+ t ϕ m (dxdt)+ t ϕ m (dxdt). T d R + T d R + T d R + (2.6) Let now A(t) := ν (x,t) σ (x,t), λ µ dx. (2.7) T d Then A is a non-negative locally integrable function on R +. Let τ > be a fixed Lebesgue point of A and define ϕ :R + R + by ( t ) ( ϕ(t)= ε χ (ε,2ε) (t)+ χ (2ε,τ) (t)+ t 2ε + τ+ ε ) χ 2ε (τ ε,τ+ε) (t). dx dt

7 RELATIVE ENTROPY METHOD FOR MEASURE-VALUED SOLUTIONS IN NATURAL SCIENCES 7 We denote by ϕ δ, δ < ε, a mollification of ϕ by a smooth kernel. Likewise m δ and mδ denote mollifications of measures m and m with respect to the time variable. Notice that t ϕ is supported only in the intervals (ε,2ε) and (τ ε,τ+ ε), where it is equal to ε and 2ε respectively. Hence using the following inequalities m (T d (a,b)) liminf Φ m(dxdt) limsup δ Φ m(dxdt) m δ (T d [a,b]), δ + R + δ + R + where Φ=χ (a,b) a.e. with respect to the Lebesgue measure, we obtain using also (2.6) 2ε A(t) dt+ (m (dxdt)+m (dxdt)) ε ε ε T d [ε,2ε] τ+ε A(t) dt (m (dxdt)+m (dxdt)) 2ε τ ε 2ε T d (τ ε,τ+ε) lim A(t) t ϕ δ dt+ t ϕ δ (m (dxdt)+m (dxdt)) δ T d R + T d R +. This implies that τ+ε A(t) dt+ 2ε τ ε 2ε 2ε A(t) dt+ ε ε ε 2ε A(t) dt+ ε ε T d (τ ε,τ+ε) T d [ε,2ε] T d [,2ε] (m (dxdt)+m (dxdt)) (m (dxdt)+m (dxdt)) (m (dxdt)+m (dxdt)). Furthermore A(t) ν (x,t) σ (x,t), λ u (x) dx+ ν (x,t) σ (x,t), µ u (x) dx T d T d = ν (x,t), λ u (x) dx+ σ (x,t), µ u (x) dx. T d T d Consequently by the initial condition (2.6) we have { T lim A(t) dt+ m (dxdt)+ T + T T T d [,T] T T d [,T] } m (dxdt) =. (2.8) (2.9) It follows that the right hand side of (2.9) converges to zero as ε. Therefore, since A and the measures m, m are non-negative, we see that τ+ε lim A(t) dt =, and lim (m (dxdt)+m (dxdt))=. (2.2) ε 2ε ε 2ε τ ε T d (τ ε,τ+ε) Hence, by Lebesgue differentiation theorem, it follows that A(τ)=for a.a. τ R +. This implies that λ µ dν (x,t) (λ)dσ (x,t) (µ)=. (2.2) R R It can be easily seen from (2.2) that the measures ν (x,t) and σ (x,t) have a common support consisting of a point w(x,t) for a.e. (x,t) T d R +. Further, the second limit in (2.2) implies that m = m =. To see this consider an arbitrary time interval [a,b] R +. It can be covered with a finite number of overlapping open intervals of radius ε, denoted B(t i,ε), so that L (B(t i,ε)) 2(b a), i J

8 8 TOMASZ DEBIEC, PIOTR GWIAZDA, KAMILA ŁYCZEK, AND AGNIESZKA ŚWIERCZEWSKA-GWIAZDA. where L denotes the -dimensional Lebesgue measure. From (2.2) we have that m (T d B(τ,ε)) C(ε)L d+ (T d B(τ,ε)) for any ball, where C(ε) as ε. Hence m (T d [a,b]) C(ε) L d+ (T d B(t i,ε)) 2C(ε)L d (T d )(b a) ε. i J Hence m (T d [a,b])=for any a<b. Similarly m =. This in turn implies m 2 = m 2 = as a consequence of (2.4). Finally the claimed regularity of w follows from the convergence g(w k (x,t)) ϕ(x,t) dx dt k ν (x,t),g ϕ(x,t) dx dt T d R + T d R + }{{} g(w(x,t)) and the boundedness of the sequence u k. Remark 2.3. One can also pass to the limit in the divergence term of (2.4) to obtain the following averaged contraction principle t ( ν(x,t) σ (x,t), λ µ + m + m ) + divx ( ν(x,t) σ (x,t),q(λ,µ) + m 2 + m 2 ). (2.22) The theorem is extendible to the problem defined on R d R + rather than T d R +. An approximation argument is needed in that case since the constant unit function no longer belongs to the class C c(r d ). We have thus established uniqueness in the class of mv entropy solutions. We conclude this section with proving existence of a unique entropy solution. First we observe that recession functions for η = λ µ and q=sgn(λ µ)( f(λ) f(µ)) can be easily described as follows: η λv k = lim = (2.23) λ ± λ where v S, and q (v)= { f (), for v=, f ( ), for v=, (2.24) where we define f f(λ) () := lim λ + λ and f f( λ) ( ) := lim λ λ (2.25) Notice that in fact the assumption on existence of the limits in (2.25) together with continuity of the flux implies condition (2.3). Theorem 2.4. There exists a unique entropy solution w to (2.) such that for a.e. t R +. w(,t) L (T d ) u L (T d ) (2.26)

9 RELATIVE ENTROPY METHOD FOR MEASURE-VALUED SOLUTIONS IN NATURAL SCIENCES 9 Proof. We adapt the proof presented in [4]. First we prove existence of a mv entropy solution (ν,m,m 2 ) - to this end we consider the following parabolic regularization t w n + div x f n (w n )= n w n, ont d R +, w n (x,)=u n (x), ont d (2.27) where u n C c (Td ) and u n u in L (T d ) and f n := f η ε n. The sequence ε n is chosen so that for z u L we have sup h ε n f(z+h) f(z) n. (2.28) This is possible because of the assumptions on asymptotic behaviour (2.25) of the flux. Problem (2.27) has a unique smooth solution w n satisfying the bound w n (,t) L (T d ) u L (T d ) for a.e. t. Since the sequence {w n } is bounded in L (R + ;L (T d )), there is, after passing to a subsequence, a triple (ν,m,m 2 ) of associated Young measure and concentration measures. We then have by (2.28) lim sgn(w n k)( f n (w n ) f n (k)) ϕ dx dt n T d R + = lim sgn(w n k)[( f(w n ) f(k))+(f n (w n ) f(w n ))+(f n (k) f(k))] ϕ dx dt n T d R + = lim sgn(w n k)( f(w n ) f(k)) ϕ n T d R + ( ) = ν(x,t),sgn(λ k)( f(λ) f(k)) +m 2 ϕ dx dt, T d R + (2.29) where the last equality follows from Lemma.3. Moreover, by virtue of the same lemma, we have ( ) lim w n k t ϕ dx dt = ν(x,t), λ k +m t ϕ dx dt. (2.3) n T d R + T d R + We observe that from Lemma.3 and (2.23)- (2.25), the following characterization of the concentration measures holds true ( ) m 2 = q (v) dν S (x,t) m. Clearly this implies the bound (2.4). We remark that the measures m and m 2 are common for each choice of entropy-entropy flux pair (i.e. for each choice of k R), because the associated recession functions do not depend on k. Let now sgn δ, θ δ and q δ be regularizations of the functions z sgn(z), z k and q(z,k) respectively, for k R. Multiplying (2.27) by ϕ sgn δ (w n k), integrating in time and space and integrating by parts yields (θ δ (w n,k) t ϕ+ q δ (w n,k) ϕ) dx dt T d R + = w n 2 sgn n T d δ (w n k)ϕ dx dt (2.3) θ R + n δ (w n,k) ϕ dx dt. T d R + The right hand side of the last identity can be bounded from below by C(ϕ) ( ) w n n L (R + ;L (T d ))+. (2.32)

10 TOMASZ DEBIEC, PIOTR GWIAZDA, KAMILA ŁYCZEK, AND AGNIESZKA ŚWIERCZEWSKA-GWIAZDA. Combining (2.3) and (2.32) and passing with δ to zero we see that T d R + w n k t ϕ+ sgn(w n k)( f n (w n ) f n (k)) ϕ dx dt C n. Then passing to infinity with n we get t ( ν(x,t), λ k +m ) + divx ( ν(x,t),q(λ,k)+m 2 ). Therefore the generalized Young measure generated by the sequence {w n } satisfies the entropy inequlity (2.5). Standard methods of analysis for PDEs can be employed to show that the initial condition (2.6) is satisfied as well. Therefore (ν,m,m 2 ) is a mv entropy solution of (2.). By the previous theorem there is a function w such that ν (x,t) = δ {w(x,t)}. This function is then the unique entropy solution. Remark 2.5. Other approximation schemes can be used rather than the viscosity approximation employed here. Indeed from the point of view of numerics other schemes are favourable. In [3], [38], [39] and [5] a kinetic approximation (i.e. an approximation by a Boltzmann type equation) is used, while in [8] a hyperbolic conservation law is realized as a limit of the attracitve zero range process (ZRP). In the latter paper a discontinuous flux is considered. An extension of the averaged contraction principle to the case of discontinuous flux (both in x and u) was considered in a series of papers [5 7, 28]. 3. RELATIVE ENTROPY METHOD FOR EQUATIONS OF FLUID DYNAMICS In this section we survey selected weak-strong uniqueness results for equations of fluid dynamics. A global existence of measure-valued solutions to the incompressible Euler system was proven in [6] for any finite-energy initial data. Later existence was also shown for the compressible Euler and Navier-Stokes systems, cf. [36], [29]. The first weak-strong uniqueness result was proven for incompressible Euler in [3]. Note however that existence of a strong solution is needed otherwise uniqueness for admissible solutions might not hold, cf. [2], [4]. 3.. Euler equations. First we consider the incompressible Euler equations t u+div(u u)+ p=, divu=, (3.) where u :T d R + R d is the velocity of a fluid and p is the scalar pressure. Definition 3.. Let ν be a Young measure, m a matrix-valued measure on T d [,T] satisfying m(dx dt)=m t (dx) dt for some family{m t } t (,T) of uniformly bounded measures onr d. Further let D L (,T) with D such that m t (T d ) CD(t) for some constant C> and almost every t [,T]. The triple(ν,m,d) is called a dissipative measure-valued solution to (3.) with initial datum u if it satisfies the following conditions: () for any divergence-free ϕ C (T d [,T];R d ) the equation T t ϕ(x,t) ν (x,t),id + ϕ(x,t) :( ν (x,t),id id +m t ) dx dt T d = ν (x,τ),id ϕ(x,τ) u (x) ϕ(x,) dx T d holds for almost every τ (, T); (3.2)

11 RELATIVE ENTROPY METHOD FOR MEASURE-VALUED SOLUTIONS IN NATURAL SCIENCES (2) the divergence free condition T d ν (x,t),id ψ(x) dx= holds for every ψ C (T d ) and a.e. t (,T); (3) the admissibility condition E(τ) u (x,τ) 2 dx 2 T d is satisfied for a.e τ (, T), where the measure-valued energy is defined by E(τ) := ν 2 (x,τ), u 2 dx+d(τ). T d We will now show how measuring the relative entropy between a dissipative mvs and a strong solution leads to a uniqueness result. We repeat the proof presented in [43]. Theorem 3.2. Let (ν,m,d) be a dissipative measure-valued solution and U C (T d [,T]) a strong solution to (3.) with the same initial datum u. Then ν (x,τ) = δ {U(x,τ)} for almost every (x,τ) T d (,T), m= and D=. Proof. We begin by defining the relative entropy E rel as E rel (τ) := ν 2 (x,τ), id U(x,τ) 2 dx+d(τ). T d This quantity can be estimated as follows. E rel (τ)= U(x,τ) 2 dx+ ν 2 T d 2 (x,τ), id U 2 dx ν (x,τ),id U(x,τ) dx+d(τ) T d T d u 2 dx+ u 2 dx ν 2 T d 2 (x,τ),id U(x,τ) dx T d T d = u 2 dx u (x) u (x) dx T d T d τ t U(x,t) ν (x,t),id + U(x,t) : ν (x,t),id id dx dt T d τ U(x,t) dm(x,t) T d τ = div(u U)(x,t) ν (x,t),id U(x,t) ν (x,t),id id dx dt T d τ sym U(x,t) dm(x,t) T d τ = ν (x,t),(u(x,t) id) sym U(x,t)(id U(x,t)) dx dt T d τ sym U(x,t) dm(x,t) T d τ sym U(t) L E rel (t) dt. It now follows from Gronwall s inequality that the relative entropy is zero almost everywhere. This in turn implies that ν (x,τ) = δ {U(x,τ)} and D(τ)=. Now consider the isentropic compressible Euler system t ρ+ div(ρu)=, t (ρu)+div(ρu u)+ (ρ γ )=, (3.3)

12 2 TOMASZ DEBIEC, PIOTR GWIAZDA, KAMILA ŁYCZEK, AND AGNIESZKA ŚWIERCZEWSKA-GWIAZDA. where γ > is the adiabatic coefficient. The definition of a measure-valued solution to the above system requires a slight refinement of the Alibert-Bouchitté framework, which we ignore here; see Section 3 in [27] for details. The following notation is used for brevity f(dx dt)= ν (x,t), f dx dt+ ν (x,t), f m(dx dt). Definition 3.3. The triple (ν,m,ν ) is called a measure-valued solution of (3.3) with initial data (ρ,u ) such that ρ and ρ u are integrable if for every τ [,T], ψ C (T d [,T]) and ϕ C (T d [,T];R n ) and T T t ψρ+ ψ ρ dx dt+ ψ(x,)ρ ψ(x,t)ρ(x,t) dx= (3.4) T d T d t ϕ ρu+ ϕ : ρu u+ divϕρ γ dx dt T d + ϕ(x,) ρ u ϕ(x,t) ρu(x,t) dx=. T d We then define the entropy of such a measure-valued solution and E mvs (t) := T d 2 ρ u 2 (x,t)+ γ ργ (x,t) dx E := T d 2 ρ u 2 (x)+ γ ργ (x) dx An admissibility criterion is then posed as follows E mvs E. It was shown in [27] that one can use the relative entropy method to prove the following weakstrong uniqueness result. (3.5) Theorem 3.4. Suppose (P,U) W, (T d [,T]) C (T d [,T]) is a strong solution of (3.3) with initial data (ρ,u ) s.t. ρ c>, ρ L r (T d ), ρ u L (T d ) and P c>. Then if (ν,m,ν ) is an admissible mv solution with the same initial data, then and m=. ν (x,t) = δ ( P(x,t), P(x,t)U(x,t)) for a.e.(x,t) The relative entropy functional used in this case has the form E rel (t) := 2 T d ρ u U 2 + γ ργ γ γ Pγ ρ+ P γ dx.

13 RELATIVE ENTROPY METHOD FOR MEASURE-VALUED SOLUTIONS IN NATURAL SCIENCES Navier-Stokes equations. We now consider the barotropic Navier-Stokes system t ρ+ div x (ρu) =, (3.6) t (ρu)+div x (ρu u)+ x p(ρ) = div x S( x u), (3.7) u =. (3.8) Here is a regular bounded domain in R 2 or R 3 and S is the Newtonian viscous stress. The following definition of a dissipative measure-valued solution and subsequent results are taken from [8]. Definition 3.5. We say that a parametrized measure {ν (x,t) } (x,t) (O,T), ν L weak( (,T);P ( R N [, ) )), ν (x,t) ;s ρ, ν (x,t) ;v u is a dissipative measure-valued solution of the Navier-Stokes system ( ) in (, T), with the initial conditions ν (x,) and dissipation defect D, D L (,T), D, if the following holds. Equation of continuity. There exists a measure r C L (,T ;M()) and χ L (,T) such that for a.a. τ (,T) and every ψ C ( [,T]), r C (τ); x ψ χ(τ)d(τ) ψ C () (3.9) and ν (x,τ) ;s ψ(,τ) dx ν ;s ψ(, ) dx τ [ ] τ = ν (x,t) ;s t ψ+ ν (x,t) ;sv x ψ dx dt+ r C ; x ψ dt. (3.) Momentum equation. u= ν (x,t) ;v L 2 (,T ;W,2 (;R N )), and there exists a measure r M L (,T ;M()) and ξ L (,T) such that for a.a. τ (,T) and every ϕ C ( [,T];R N ), ϕ =, r M (τ); x ϕ ξ(τ)d(τ) ϕ C () (3.) and ν (x,τ) ;sv ϕ(,τ) dx ν ;sv ϕ(, ) dx τ [ ] = ν (x,t) ;sv t ϕ+ ν (x,t) ;s(v v) : x ϕ+ ν (x,t) ; p(s) div x ϕ τ τ S( x u) : x ϕ dx dt+ r M ; x ϕ dt. Energy inequality. ( ) ν (x,τ) ; 2 s v 2 + P(s) dx+ τ S( x u) : x u dx dt+d(τ) ( ) ν (x,) ; 2 s v 2 + P(s) dx dx dt (3.2) (3.3)

14 4 TOMASZ DEBIEC, PIOTR GWIAZDA, KAMILA ŁYCZEK, AND AGNIESZKA ŚWIERCZEWSKA-GWIAZDA. for a.a. τ (, T). In addition, the following version of Poincaré s inequality" holds for a.a. τ (,T): τ ν(x,t) ; v u 2 dx dt c P D(τ). (3.4) One can show (Theorem 2. in [8]) that if the pressure satisfies the following coercivity assumptions p C[, ) C 2 (, ), p()=, p (ρ)> for ρ >, liminf ρ P(ρ) p (ρ)>, liminf ρ p(ρ) >, (3.5) then there exists a dissipative mv solution with a prescribed finite-energy initial data. The following weak-strong uniqueness result can then be proven Theorem 3.6. Let R N, N = 2,3 be a bounded smooth domain. Suppose the pressure p satisfies (3.5). Let {ν (x,t),d} be a dissipative measure-valued solution to the barotropic Navier- Stokes system ( ) in (,T), with the initial state represented by ν (x,), in the sense specified in Definition 3.5. Let [r,u] be a strong solution of ( ) in (,T) belonging to the class r, x r, U, x U C( [,T]), t U L 2 (,T ;C(;R N )), r>, U =. Then there is a constant Λ=Λ(T), depending only on the norms of r, r, U, χ, and ξ in the aforementioned spaces, such that ν (x,τ) ; 2 s v U 2 + P(s) P (r)(s r) P(r) dx τ + Λ(T) x u x U 2 dx dt+d(τ) ν (x,) ; 2 s v U(, ) 2 + P(s) P (r(, ))(s r(, )) P(r(, )) for a.a. τ (, T). In particular, if the initial states coincide, meaning ν (x,) = δ [r(x,),u(x,)] for a.a. x then D =, and ν (x,τ) = δ [r(x,τ),u(x,τ)] for a.a. τ (,T), x. 4. POLYCONVEX ELASTODYNAMICS In this section we consider the system of elasticity 2 y = S( y), (4.) t2 where y : Q R + R 3 stands for the motion, F = y, v= t y, and S stands for the Piola-Kirchoff stress tensor obtained as the gradient of a stored energy function, S = W F. Here we assume that W is polyconvex, that is W(F) = G(Φ(F)) where G : Mat 3 3 Mat 3 3 R [, ) is a strictly convex function and Φ(F)=(F,cofF,det F) Mat 3 3 Mat 3 3 R stands for the vector of null- Lagrangians: F, the cofactor matrix coff and the determinant detf. It is observed in [3] and [4] dx

15 RELATIVE ENTROPY METHOD FOR MEASURE-VALUED SOLUTIONS IN NATURAL SCIENCES 5 that this system can be embedded into the following symmetrizable hyperbolic system in a new dependent variable Ξ=(F,Z,w) taking values in Mat 3 3 Mat 3 3 R v i = ( ) G t x α Ξ A(Ξ) ΦA (F), F iα Ξ A = ( ) (4.2) Φ A t x α (F)v i. F iα This system admits the following entropy-entropy flux pair η(v,f,z,w)= 2 v 2 + G(F,Z,w), G q α = v i Ξ A(Ξ) ΦA (F). F iα A strong solution to (4.) is a function w W,. It automatically satisfies (4.3) t η+ α q α =. (4.4) A weak solution which satisfies (4.4) as an inequality is called an entropy weak solution. The following definition of a dissipative mv solution is taken from [3]. Definition 4.. Let the pair (y,ν) consist of a map y, with distributional time and space derivatives (v,f) L (L 2 ) L (L p ) and a Young measure ν = (ν (x,t) ) x,t QT generated by a sequence satisfying sup η(v ε,f ε,z ε,w ε )dx < ε,t which represents weak limits in the following way: wk- lim f(v ε,f ε,z ε,w ε )= f(λ v,λ Ξ )dν (x,t) (λ v,λ Ξ ) ε continuous f = f(λ v,λ Ξ ) with lim λ v + λ Ξ f(λ v,λ Ξ ) 2 λ v 2 + G(λ Ξ ) = where λ v R 3, λ Ξ =(λ F,λ Z,λ w ) Mat 3 3 Mat 3 3 R=R 9. The Young measure is connected with the map y through the requirements that (almost everywhere) (4.5) F = ν,λ F, v= ν,λ v, Ξ= ν,λ X i. (4.6) The pair (y,ν) is a measure-valued solution to (4.) if for i=,2,3 t v i α ν, G Ξ A(λ Ξ) ΦA F iα (λ F ) = (4.7) and for A=,...,9 in distributions with t Φ A (F) α ( Φ A F iα (F)v i ) = (4.8) Ξ=Φ( ν,λ F )=Φ(F). (4.9) The solution is said to be a dissipative measure-valued solution with concentration if it is a measure-valued solution which verifies in addition: dθ ( ) ν,η +γ dx dt+ θ()η (x)dx, dt for all non-negative functions θ = θ(t) C c ([,T)) with θ. Here η means the entropy η evaluated on the initial data and γ is the non-negative concentration measure.

16 6 TOMASZ DEBIEC, PIOTR GWIAZDA, KAMILA ŁYCZEK, AND AGNIESZKA ŚWIERCZEWSKA-GWIAZDA. Under the following additional growth assumptions on the function G: (H) G C 3 (Mat 3 3 Mat 3 3 R;[, )) is a strictly convex function satisfying for some γ > the bound D 2 G γ >, (H2) G(F,Z,w) c ( F p + Z q + w r + ) c 2 where p (4, ), q,r [2, ), (H3) G(F,Z,w) c( F p + Z q + w r + ), (H4) F G + Z G p p + w G p p 2 o()( F p + Z q + w r +) where o() as Ξ, existence of dissipative mv solutions as well as a weak-strong uniqueness result are proven, cf [4] and [3]. Theorem 4.2. Let G satisfy (H) (H3),(H4) and let(y,ν,γ) be a dissipative measure-valued solution in the sense of definition 4.. If the initial data equal those of a Lipschitz bounded solution ( ˆv, ˆF) W, (Q T ): (v(x,),ξ(x,)) =( ˆv(x,),Φ( ˆF(x,))) then γ is zero, (v,ξ)=( ˆv,Φ( ˆF)) and ν = δ ˆv,Φ( ˆF). 5. WEAK-STRONG UNIQUENESS FOR GENERAL HYPERBOLIC CONSERVATION LAWS In [9] the method of relative entropy is extended to a more general class of problems of hyperbolic and hyperbolic-parabolic type. In this section we describe the results of that paper. Consider the following hyperbolic problem t A(u)+ α F α (u)=, (5.) where u(t,x) takes values in R n, t R +, x R d and A,F α :R n R n are given smooth functions with α =,,d. It is assumed that this system is symmetrizable. The following hypotheses are assumed throughout. (H ) A :R n R n is a C 2 globally invertible map, (H 2 ) existence of an entropy-entropy flux pair (η,q), that is G :R n R n, G=G(u) smooth such that η = G A q α = G F α, α =,...,d, (H 3 ) the symmetric matrix 2 η(u) G(u) 2 A(u) is positive definite, together with the following growth assumptions on the entropy η(u), the functions F α and A β ( u p + ) B η(u) β 2 ( u p + ) for u R n (A ) for some positive constants β, β 2, B and for some p (, ). F α (u) η(u) = o() as u, α =,...,d, (A 2 ) A(u) η(u) = o() as u. (A 3) C ( A(u) q + ) B η(u) C( A(u) q + ), q>, (A 4 ) for some uniform constant C > and B>.

17 RELATIVE ENTROPY METHOD FOR MEASURE-VALUED SOLUTIONS IN NATURAL SCIENCES 7 Definition 5.. A dissipative measure-valued solution (u, ν, γ) with concentration to (5.) consists of u L (L p ), a Young measure ν = (ν (x,t) ) {(x,t) Q T } and a non-negative Radon measure γ M + (Q T ) such that u(x,t)= ν (x,t),λ and ν (x,t),a i (λ) t ϕ i dxdt+ ν (x,t),f i,α (λ) α ϕ i dxdt+ ν (x,),a i ϕ i (x,)dx= (5.2) for i=,...,n and any ϕ C c(q [,T)), and dξ [ ν(x,t),η(λ) dxdt+ γ(dxdt) ] + dt for all ξ = ξ(t) C c ([,T)) with ξ. ξ() [ ν (x,),η dx+γ (dx) ], (5.3) The following theorem provides weak measure-valued versus strong uniqueness in the L p framework for < p<. Theorem 5.2. Suppose that (H) (H3) hold, the growth properties (A ) (A 4 ) are satisfied, and the entropy η(u). Assume that (u,ν,γ) is a dissipative measure-valued solution, u = ν (x,t),λ, and ū W, (Q T ) is a strong solution to (5.). Then, if the initial data satisfy γ = and ν x = δū (x), then ν = δū and u=ū almost everywhere on Q T. Remark 5.3. The work of Christoforou and Tzavaras [9] is a generalization of the results presented in [3], where the case A(u)=u was considered and concentration effects were ignored. 6. GENERAL RELATIVE ENTROPY METHOD IN MATHEMATICAL BIOLOGY In this section we give a short overview of an extension of relative entropy method to linear PDEs which was introduced in the context of biological systems in [34] and further developed in [32] and [33]. It was extended to measure solutions with measure initial data in [26]. Here the relative entropies are a family of renormalizations to the initial linear problem, obtained by multiplying the original equation by a nonlinear function. The notation used in this section is intentionally inconsistent with the remainder of this paper however it is consistent with the notation used in the above mentioned papers on the subject. Following [37] and [26] we consider the McKendrick-Von Foerster equation t n(t,x)+ x n(t,x) = on (R + ) 2, n(t,x=)= n(x,t = )=n (x). B(y)n(t, y)dy, (6.) Here, n(x,t) denotes the population density at time t with age x, and B L (R + ;R + ) is a birth rate such that B(x)dx>. The associated primal and dual eigenvalue problems have the form x N(x)+λ N(x)=, x, N()= N >, B(y)N(y)dy, N(x)dx= (6.2)

18 8 TOMASZ DEBIEC, PIOTR GWIAZDA, KAMILA ŁYCZEK, AND AGNIESZKA ŚWIERCZEWSKA-GWIAZDA. and x ϕ(x)+λ ϕ(x)=ϕ()b(x), x, ϕ, N(x)ϕ(x)dx=, where λ >. Under the above assumptions both these problems posses a unique solution. In fact, it can be seen that the solution of (6.2) is given by N(x)=λ e λ x. Since the death rate is ignored and the birth rate integrates to more than one, an exponential growth of the population is expected. In order to quotient out this growth, we set ñ(t,x)= n(t,x)e λ t. Then (6.) becomes t ñ(t,x)+ x ñ(t,x)+λ ñ(t,x)= on (R + ) 2, ñ(t,x=)= ñ(x,t = )=n (x). B(y)ñ(t, y)dy, It can be shown that (6.3) has a unique solution in the weak sense for any n M + ([, )), see [25]. The following results on long-time asymptotics of this solution are proven in [26]. Theorem 6.. Let n M + ([; )). Then there is y >, σ > and a bounded function η, positive on supp ϕ, such that the solution of the renewal equation satisfies η(x)d ñ(t,x) m N(x)dx e σ(t y ) where m = ϕ(x)dn (x). (6.3) η(x)d ñ (x) m N(x)dx, (6.4) Theorem 6.2. Assume in addition there exists C> such that B(x) Cϕ(x). Let n M + ([; )). Then the solution of the renewal equation satisfies lim t where m = ϕ(x)dn (x). ϕ(x)d ñ(t,x) m N(x)dx =, (6.5) 7. ACKNOWLEDGEMENT This work was partially supported by the Simons Foundation grant 3463 and the Polish Government MNiSW matching fund. T.D acknowledges the support of the National Science Centre, DEC-22/5/E/ST/228. The research was partially supported by the Warsaw Center of Mathematics and Computer Science. P.G and A.Ś-G received support from the National Science Centre (Poland), 25/8/M/ST/75. REFERENCES [] J.-J. Alibert and G. Bouchitté. Non-uniform integrability and generalized young measure. J. Convex Anal., 4:29 48, 997. [2] P. Bella, E. Feireisl, A. Novotný. Dimension reduction for compressible viscous fluids. Acta Appl. Math. 34: 2, 24. [3] Y. Brenier, C. De Lellis and L. Székelyhidi Jr. Weak-strong uniqueness for measure-valued solutions. Comm. Math. Phys., 35(2):35 36, 2. [4] J. Březina, E. Feireisl. Measure-valued solutions to the complete Euler system. arxiv:72.487v. [5] M. Bulíček, P. Gwiazda, A. Świerczewska-Gwiazda. On unified theory for scalar conservation laws with fluxes and sources discontinuous with respect to the unknown. J. Diff. Eq., 262:33 364, 27.

19 RELATIVE ENTROPY METHOD FOR MEASURE-VALUED SOLUTIONS IN NATURAL SCIENCES 9 [6] M. Bulíček, P. Gwiazda, A. Świerczewska-Gwiazda. Multi-dimensional scalar conservation laws with fluxes discontinuous in the unknown and the spatial variable. Math. Models Methods Appl. Sci., 23(3):47 439, 23. [7] M. Bulíček, P. Gwiazda, J. Malek, A. Świerczewska-Gwiazda. On scalar hyperbolic laws with discontinuous flux. Math. Models Methods Appl. Sci., 2(), 2. [8] G-Q. Chen, N. Even, M. Klingenberg. Hyperbolic conservation laws with discontinuous fluxes and hydrodynamic limit for particle systems. J. Diff. Eq, 245: , 28. [9] C. Christoforou, A.E. Tzavaras. Relative entropy for hyperbolic-parabolic systems and application to the constitutive theor of thermoviscoelasticity. arxiv: [] C. Dafermos. The second law of thermodynamics and stability. Arch. Rational Mech. Anal. 7:67 79, 979. [] C. Dafermos. Hyperbolic conservation laws in continuum physics. Third edition. Springer-Verlag, Berlin, 2. [2] C. De Lellis, L. Székelyhidi Jr. On admissibility criteria for weak solutions of the Euler equations. Arch. Rational Mech. Anal. 95():225 26, 2. [3] S. Demoulini, D.M.A. Stuart and A.E. Tzavaras. Weak-strong uniqueness of dissipative measure-valued solutions for polyconvex elastodynamics. Arch. Ration. Mech. Anal., 25(3):927 96, 22. [4] S. Demoulini, D.M.A. Stuart, A.E. Tzavaras, A variational approximation scheme for three-dimensional elastodynamics with polyconvex energy. Arch. Rational Mech. Analysis 57: , 2. [5] R. J. DiPerna. Measure-valued solutions to conservation laws. Arch. Ration. Mech. Anal., 88(3):223 27, 985. [6] R. J. DiPerna and A. J. Majda. Oscillations and concentrations in weak solutions of the incompressible fluid equations. Comm. Math. Phys., 8(4): , 987. [7] E. Feireisl, B.J. Jin and A. Novotný. Relative entropies, suitable weak solutions, and weak-strong uniqueness for the compressible Navier-Stokes system. J. Math. Fluid Mech., 4:72 73, 22. [8] E. Feireisl, P. Gwiazda, A. Świerczewska-Gwiazda and E. Wiedemann. Dissipative measure-valued solutions to the compressible Navier-Stokes system. Calc. Var. PDE, 55(6):Art. 4, 2 pp., 26. [9] E. Feireisl, A. Novotný. Weak-strong uniqueness property for the full Navier-Stokes-Fourier system. Arch. Rational Mech. Anal. 24:683 76, 22. [2] U.S. Fjordholm, R. Käppeli, S. Mishra, E. Tadmor. Construction of approximate entropy measure-valued solutions for hyperbolic systems of conservation laws. Found. Comput. Math. 7(3): , 27. [2] U.S. Fjordholm, S. Lanthaler, S. Mishra. Statistical solutions of hyperbolic conserwation laws I: Foundations. arxiv: [22] U. S. Fjordholm, S. Mishra and E. Tadmor. On the computation of measure-valued solutions. Acta Numerica, 25: , 26. [23] J. Giesselmann and A.E. Tzavaras. Stability properties of the Euler-Korteweg system with monotone pressures. Applicable Analysis, 27. [24] P. Gwiazda. On measure-valued solutions to a two-dimensional gravity-driven avalanche flow model. Math. Methods Appl. Sci., 28(8): , 25. [25] P. Gwiazda, T. Lorenz, A. Marciniak-Czochra. A nonlinear structured population model: Lipschitz continuity of measure valued solutions with respect to model ingredients. J. Diff. Eq. 248: , 2. [26] P. Gwiazda, E. Wiedemann Generalized entropy method for the renewal equation with measure data. Commun. Math. Sci. 5(2): , 27. [27] P. Gwiazda, A. Świerczewska-Gwiazda and E. Wiedemann. Weak-strong uniqueness for measure-valued solutions of some compressible fluid models. Nonlinearity, 28(): , 25. [28] P. Gwiazda, A. Świerczewska-Gwiazda, P. Wittbold, A. Zimmermann. Multi-dimensional scalar balance laws with discontinuous flux. J. Funct. Anal., 267(8): , 24. [29] D. Kröner, W. Zajączkowski. Measure-valued solutions of the Euler system for ideal compressible polytropic fluids. Math. Methods Appl. Sci. 9(3): , 996. [3] S.N. Kružkov. First order quasilinear equations in several independent variables. Math. USSR-Sb :27 243, 97. [3] P.L. Lions, B. Perthame, E. Tadmor. A Kinetic formulation of scalar multidimensional conservation laws. J. AMS 7:69 9, 994. [32] P. Michel, P. Mischler and B. Perthame. General entropy equations for structured population models and scattering. C. R. Acad. Sci. Paris, Ser. I 338: , 24. [33] P. Michel, P. Mischler and B. Perthame. General relative entropy inequality: an illustration on growth models. J. Math. Pures Appl. 84: , 25. [34] S. Mischler, B. Perthame, L. Ryzhik. Stability in a Nonlinear Population Maturation Model. Math. Models Methods Appl. Sci. 2:75 772, 22. lin, Lecture Notes in Math. 48: 6.

20 2 TOMASZ DEBIEC, PIOTR GWIAZDA, KAMILA ŁYCZEK, AND AGNIESZKA ŚWIERCZEWSKA-GWIAZDA. [35] J. Necas, J. Malek, M. Rokyta, M. Ruzicka. Weak and Measure-Valued Solutions to Evolutionary PDEs. Chapman and Hall/CRC, 996. [36] J. Neustupa. Measure-valued solutions of the Euler and Navier-Stokes equations for compressible barotropic fluids. Math. Nachr., 63:27 227, 993. [37] B. Perthame. Transport Equations in Biology. Frontiers in Mathematics, Birkhäuser Basel, 27. [38] B. Perthame, E. Tadmor. A kinetic equation with kinetic entropy functions for scalar conservation laws. Commun. Math. Phys. 36:5 57, 99. [39] B. Perthame, A.E. Tzavaras. Kinetic formulation for systems of two conservation laws and elastodynamics. Arch. Ration. Mech. Anal. 55(): 48, 2. [4] L. Székelyhidi Jr., E. Wiedemann. Young measures generated by ideal incompressible fluid flows. Arch. Rational Mech. Anal. 26(): , 22. [4] A. Szepessy An existence result for scalar conservation laws using measure valued solutions. Comm. PDE, 4():329 35, 989. [42] L. Tartar. Compensated compactness method applied to systems of conservation laws. Systems Nonlinear PDE, Ball, J.M (ed.) NATO ASI Series, C.Reidel Publishing Col., 983. [43] E. Wiedemann. Weak-strong uniqueness in fluid dynamics. arxiv: Tomasz Dębiec: INSTITUTE OF APPLIED MATHEMATICS AND MECHANICS, UNIVERSITY OF WARSAW, BA- NACHA 2, 2-97 WARSZAWA, POLAND address: t.debiec@mimuw.edu.pl Piotr Gwiazda: INSTITUTE OF MATHEMATICS, POLISH ACADEMY OF SCIENCES, ŚNIADECKICH 8, -656 WARSZAWA, POLAND, AND INSTITUTE OF APPLIED MATHEMATICS AND MECHANICS, UNIVERSITY OF WAR- SAW, BANACHA 2, 2-97 WARSZAWA, POLAND address: pgwiazda@mimuw.edu.pl Kamila Łyczek: INSTITUTE OF APPLIED MATHEMATICS AND MECHANICS, UNIVERSITY OF WARSAW, BA- NACHA 2, 2-97 WARSZAWA, POLAND address: k.lyczek@mimuw.edu.pl Agnieszka Świerczewska-Gwiazda: INSTITUTE OF APPLIED MATHEMATICS AND MECHANICS, UNIVERSITY OF WARSAW, BANACHA 2, 2-97 WARSZAWA, POLAND address: aswiercz@mimuw.edu.pl