On Lagrangian solutions for the semi-geostrophic system with singular initial data
|
|
- Everett Fisher
- 5 years ago
- Views:
Transcription
1 On Lagrangian solutions for the semi-geostrophic system with singular initial data Mikhail Feldman Department of Mathematics University of Wisconsin-Madison Madison, WI 5376, USA Adrian Tudorascu Department of Mathematics West Virginia University Morgantown, WV 2656, USA January 8, 213 Abstract We show that weak (Eulerian) solutions for the Semi-Geostrophic system in physical space exhibiting some mild regularity in time cannot yield point masses in dual space. However, such solutions are physically relevant to the model. Thus, we discuss a natural generalization of weak Lagrangian solutions in the physical space to include the possibility of singular measures in dual space. We prove existence of such solutions in the case of discrete measures in dual space. We also prove that weak Lagrangian solutions in physical space determine solutions in the dual space. This implies conservation of geostrophic energy along the Lagrangian trajectories in the physical space. 1 Introduction The Semi-Geostrophic (abbreviated SG in this work) equations have been proposed as simplifications of the primitive equations (Boussinesq) when the rate of change of momentum is much smaller than the Coriolis force (small Rossby number) [16]. The advected quantity momentum is approximated by its geostrophic value, but the trajectories are not. Throughout the entire paper, IR 3 is a given open, bounded set, and T (, ) is fixed. A Key words: SG system, flows of maps, optimal mass transport, Wasserstein metric, optimal maps, absolutely continuous curves. 1
2 2 version of the 3D Semi-Geostrophic system [6, 9] is D t X = J [ X x ], X = P, u =, u ν =, on [, T ), P (, ) = P in, (1) where P is convex defined on. One looks for solutions (P, u) satisfying the Cullen-Purser stability condition (see, e.g., [1]), which amounts to imposing that P t ( ) := P (t, ) be convex for all t [, T ). Here, 1 J = 1. Henceforth, we shall assume without loss of generality that L 3 () = 1 (otherwise all the measures considered will have total mass equal to the volume of ) and denote χ := L 3. If P t# χ =: α t are all absolutely continuous with respect to L 3, then one can use the Legendre- Fenchel transforms Pt of P t to formally rewrite (1) as the so-called SG in dual variables t α + (Uα) = in [, T ) IR 3, (2) P (t, )#χ = α(t, ) for any t [, T ); (3) U(t, X) = J[X P (t, X)], (4) α(, X) = α (X) for a.e. X IR 3. (5) Existence of solutions of the problem recast in dual variables was obtained by J. D. Benamou, Y. Brenier [6], and for some related models by M. Cullen and W. Gangbo [1], and M. Cullen and H. Maroofi [12]. They considered case when α L q () for some q > 1, and the solution satisfies α(t, ) L q for all t. These results were extended to the case q = 1 in [18]. These solutions are not known to be regular enough to be translated into Eulerian solutions of the problem in physical space. Recently, existence of Eulerian solutions for a class of initial data, where the conditions include the requirement that the support of α = P # χ in the dual space is the whole space, was obtained by L. Ambrosio, M. Colombo, G. De Philippis, A. Figalli [2, 3] based on the results of G. De Philippis and M. Figalli [13] on regularity of solutions for the Monge-Ampere equation. Existence of Eulerian solutions in physical space for more general initial data, when the support of α in the dual space may have a non-empty boundary, is presently not known. Another approach is to consider Lagrangian solutions in the physical space. Such solutions were introduced by M. Cullen and M. Feldman [9], and existence of Lagrangian solutions of (1) was shown in [9] for the case α = P # χ L q () for q > 1, on the basis of Ambrosio s theory on transport equations and flows of BV vector fields [1]. These results were extended to the case q = 1 in [14]. The definition of weak Lagrangian solutions in the physical space is following (where we slightly modify the definition given in [9] by relaxing assumptions on P ). Definition 1.1. Let P H 1 () be convex, F : [, T ) be a Borel map such that F C([, T ); L 2 (; IR 3 )), and let P C([, T ); H 1 ()). Assume P (t, ) is convex in for each t [, T ). Then the pair (P, F ) is called a weak Lagrangian solution of (1) in [, T ) if
3 3 i. F (, x) = x, P (, x) = P (x) for a.e. x, ii. for any t > the mapping F t = F (t, ) : is Lebesgue measure preserving, in the sense that F t# χ = χ; iii. There exists a Borel map F : [, T ) such that for every t (, T ) the map F t = F (t, ) : is Lebesgue measure preserving: F t# χ = χ, and satisfies F t F t (x) = x and F t F t (x) = x for a.e. x ; iv. The function is a distributional solution of Z(t, x) = P (t, F t (x)) (6) t Z(t, x) = J [ Z(t, x) F (t, x) ] in [, T ), Z(, x) = P (x) in. Note that the sense in which (7) must be satisfied is [Z(t, x) t ϕ(t, x) + J [ Z(t, x) F (t, x) ] ϕ(t, x)]dxdt + T for any ϕ C 1 c ([, T ) ; IR 3 ). (7) P (x) ϕ(, x)dx = (8) In this paper we consider the case of singular initial data, i.e. when α = P # χ is a singular measure. The dual problem in this case was studied by G. Loeper [17], and L. Ambrosio and W. Gangbo [4]. Note that equation (2) represents the fact that the dual flow t α t is weakly (in the sense of distributions) transported by the dual velocity U defined by (4). The change of variable X = P t (x) is reversible (although there may not be enough smoothness to transport the dual space solutions back to physical space even in this case) if and only if Pt pushes α t forward to χ, one such situation being provided if α t L 3. In general, if Pt is not necessarily the a.e. inverse of P t, then Pt may not necessarily push α t forward to χ, α t may not necessarily be absolutely continuous with respect to L 3 and the equation (4) in the dual-variable system must be generalized to (see [4]) U(t, X) = J [ X γ t (X) ], (9) where γ t is the barycentric projection onto α t of the (unique) optimal Kantorovich plan [2] γ t := ( P t Id ) # χ having α t and χ as first and second marginals, respectively. It is defined by (see [4], [5]) ξ(x) γ t (X)dα t (X) = ξ(x) ydγ t (X, y) (1) IR 3 IR 3 for all continuous ξ : IR 3 IR 3 of at most quadratic growth. Since χ is absolutely continuous with respect to L 3, we deduce ξ(x) γ t (X)dα t (X) = ξ( P t (x)) xdx. (11) IR 3 To justify generalizing (4) to (9), let us assume (1) has a smooth, stable solution (P, u). Set Y := P and compute, for any ζ Cc (IR 3 ), d ζ(x)dα t (X) = t Y t (x) ζ(y t (x))dx dt IR 3 [ = (u(t, x) )Yt (x) ] ζ(y t (x))dx + J [ Y t (x) x ] ζ(y t (x))dx.
4 4 Since u is divergence-free and Y t = 2 P t is symmetric, we deduce that the first integral in the right hand side of the second equation of the above display is [( Yt (x) ) t ] [( u(t, x) ζ(yt (x))dx = Yt (x) ) t ζ(yt (x)) ] u(t, x)dx = [ ζ(y t (x)) ] u(t, x)dx =. As for the second integral, we have, due to J t = J, Y t# χ = α t and (11), J [ Y t (x) x ] ζ(y t (x))dx = Y t (x) (J ζ)(y t (x))dx + x (J ζ)(y t (x))dx = (JX) ζ(x)dα t (X) [J γ t (X)] ζ(x)dα t (X). IR 3 IR 3 Thus, which implies d ζ(x)dα t (X) = J [ X γ t (X) ] ζ(x)dα t (X), (12) dt IR 3 IR 3 is satisfied in the sense of distributions. t α(t, X) + X (J [ X γ(t, X) ] α(t, X) ) = (13) Existence of a solution for the dual problem (2), (3), (9), (5) for any probability measure α with finite second moment was proved in [4]. In this paper we define weak Lagrangian solutions in physical space for any convex initial data P (i.e. any initial probability measure α = P # χ), and prove existence of such solutions in the case of discrete measures. Moreover, we show conservation of geostrophic energy along the Lagrangian trajectories in the physical space. One stand-out feature of Lagrangian solutions in physical space in the case of singular initial data is the absence of the Lebesgue measure preserving property for the flow map F t, given in (ii) of Definition 1.1. Heuristically, since the solution of the dual problem (2), (3), (9), (5) does not provide enough information to separate the particle paths (in physical space) within the subgradient set P t (X) for X supp α t, we define the flow map F t which collapses the set P t (X) to one point, its α t -barycenter γ t (X) P t (X), and thus the χ-measure carried by the set P t (X) should now be concentrated at the point γ t (X). Then the measure preserving property F t# χ = χ is replaced by F t# χ = µ t for t >, where µ t := ( γ t P t ) # χ. (14) Note that in the case when α t is absolutely continuous with respect to the Lebesgue measure, we have ( γ t (X) P t )(x) = ( P t P t )(x) = x for a.e. x, and thus we recover the Lebesgue measure preserving property of F t given in (ii) of Definition 1.1. We justify this definition of F t in the case of singular data, and its generalized measure-preserving property (14), by giving, in Section 4.1 below, an explicit example of a sequence of Lagrangian solutions (P (i), F (i) ) in the sense of Definition 1.1 with absolutely continuous initial measures α i space, such that the sequence (P (i), F (i) ) converges weakly to a Lagrangian solution (P, F ) in the sense of Definition 4.12 with a singular initial measure P # χ = δ z. We show that this limiting solution (P, F ) satisfies the (generalized) measure-preserving property (14), yet does not satisfy (ii) of Definition 1.1. Thus, Definition 4.12 appears to be a natural extension of Definition 1.1 (i) := P # χ in the dual
5 5 to the case of singular initial data if we expect to have stability of weak Lagrangian solutions in physical space. We note that in the case of Lagrangian solutions with absolutely continuous initial measure α = P # χ, a stability property was proved in [14]. For solutions with singular initial data α, stability is not proved at present, although the example described above suggests that some stability properties may be expected in this case too. The rest of the paper is organized as following. In the next section we collect some results on the existence and properties of solutions in dual space. In Section 3 we prove that, under a mild time-regularity assumption, Eulerian solutions do not exist if the dual-space solutions have point masses at non-negligibly many times. This will advocate finding a suitable notion of weak Lagrangian solutions in physical case that makes sense in the general case of measurevalued solutions α in the dual space. This will be achieved in Section 4. There we will also prove that such solutions exist in the case where α t are convex combinations of point masses, or equivalently, when P is the maximum of finitely many affine functions. In Section 5 we will show that weak Lagrangian solutions in physical space can be translated into weak solutions of the problem in the dual space. This will lead to conservation of energy for weak Lagrangian solutions for the Semi-Geostrophic system. 2 Solutions in dual space In this section we collect a number of results on the existence and some properties of solutions in dual space. Our main source is [4]. Before that, let us recall the definitions of some important objects. In the spirit of [5], one defines AC p (, T ; P 2 (IR 3 )) (for 1 p ) as the set of all paths µ : [, T ] t µ t P 2 (IR 3 ) for which there exists β L p (, T ) such that W 2 (µ s, µ t ) t s β(τ)dτ for all s t T, where W 2 is the quadratic Wasserstein distance [2]. The smallest of the functions β satisfying the inequality above is called the metric derivative of the curve µ, it is denoted by µ and it satisfies [5] that µ W 2 (µ s, µ t ) (t) = lim for a.e. t (, T ). s t s t There exists a Borel velocity v : (, T ) IR 3 IR 3 transporting µ in the sense of distributions, i.e. t µ + x (vµ) = in D ((, T ) IR 3 ), (15) such that v(t, ) L 2 (µ t ; IR 3 ) and v t L 2 (µ t;ir 3 ) = µ (t) for a.e. t (, T ). It turns out that this velocity (called of minimal norm ) minimizes v t L 2 (µ t;ir 3 ) among all possible Borel velocities (i.e. satisfying (15)), and it can be uniquely selected (µ t a.e.) for a.e. t (, T ) by requiring that it belong to the closure of Cc (IR 3 ) in L 2 (µ t ; IR 3 ), which is denoted [5] by T µt P 2 (IR 3 ) and called the tangent space to P 2 (IR 3 ) at µ t. Finally, if F : P 2 (IR 3 ) IR is lower semicontinuos with respect to the topology induced by the Wasserstein distance W 2, then we can define the subdifferential of F at some µ P 2 (IR 3 ) as the set of all ξ L 2 (µ; IR 3 ) such that F (ν) F (µ) + ξ(x) (y x)γ(dx, dy) + o ( W 2 (µ, ν) ) IR 3 IR 3 for all ν P 2 (IR 3 ) and all the optimal transport plans [2] γ between µ and ν. This set is denoted by F (µ) and its element of minimal L 2 (µ; IR 3 ) norm is denoted by w F (µ) and called [4] the Wasserstein gradient of F at µ.
6 6 Adapting to our context (by, for example, replacing J by J t = J which still satisfies the required orthogonality property in [4]), a Hamiltonian ODE solution/trajectory for the lower semicontinuous Hamiltonian H : P 2 (IR 3 ) IR is defined as follows. Definition 2.1. A curve µ AC 1 (, T ; P 2 (IR 3 )) with the property that there exists a vector field v : (, T ) IR 3 IR 3 such that and (, T ) t v t L 2 (µ t;ir 3 ) belongs to L1 (, T ) (16) t µ (Jvµ) =, µ = µ as distributions, (17) v t T µt P 2 (IR 3 ) H(µ t ) for a.e. t (, T ), (18) is called a solution of the Hamiltonian ODE associated to the Hamiltonian H, starting at µ. Fix ν P 2 (IR 3 ) and define H ν : P 2 (IR 3 ) IR, H ν (µ) = 1 2 W 2 2 (µ, ν), (19) which is, obviously, continuous. Then H ν (µ) consists of all functions of the type γ Id, where γ is any optimal plan between µ and ν [2], and γ is its barycentric projection onto µ, i.e. given by (see, e.g. [5]) y ξ(x)γ(dx, dy) = IR 3 IR 3 γ(x) ξ(x)µ(dx) IR 3 for all ξ C(IR 3 ; IR 3 ) of at most quadratic growth. In light of the fact that the y marginal of γ is ν, the above definition implies γ L 2 (µ;ir 3 ) Id L 2 (ν;ir 3 ). Thus, W2 2 (µ, ν) = x 2 µ(dx) + IR 3 y 2 ν(dy) 2 IR 3 x γ(x)µ(dx) γ Id 2 IR 3 L 2 (µ;ir 3 ) for all such optimal plans, which implies w H ν (µ) L 2 (µ;ir 3 ) 1 + 2H ν ( µ) =: a( µ) whenever W 2 (µ, µ) 1. (2) Then one combines Lemma 7.6 [4] and Theorem 7.4 [4] to obtain on [, 1/a( µ) =: T ] a solution µ AC (, T ; P 2 (IR 3 )) (Lipschitz curve) of the Hamiltonian ODE starting at µ and satisfying that it is a( µ) Lipschitz and conservative, i.e. [, T ] t H ν (µ t ) is constant. Thus, a(µ T ) = a( µ), which means that (2) is also satisfied with the same bounds if we replace µ by µ T. Thus, in light of the same results from [4], we infer that we can extend the solution to [T, 2T ] (whereas it preserves the Lipschitz bound and conserves the Hamiltonian). By induction, we obtain: Theorem 2.2. For any µ P 2 (IR 3 ) there exists a solution µ AC (, ; P 2 (IR 3 )) of (17) and (18), which starts at µ, conserves [, ) t H ν (µ t ), and is globally a( µ) Lipschitz. Let us now specialize to the case ν = χ (defined in the introduction) and denote H : P 2 (IR 3 ) IR, H(µ) = 1 2 W 2 2 (µ, χ). (21) The absolute continuity of χ with respect to L 3 adds the benefit of the fact that there is a unique optimal plan now between µ and χ for every µ P 2 (IR 3 ), namely γ µ = ( Φ Id) # χ, where Φ is the optimal map pushing χ forward to µ. Thus, we have H(µ) T µ P 2 (IR 3 ) = { γ µ Id} = { w H(µ)} for all µ P 2 (IR 3 ). This means that Theorem 2.2 has the following:
7 7 Corollary 2.3. Let α P 2 (IR 3 ) be given. Then there exists a distributional solution α AC (, ; P 2 (IR 3 )) for (13) with α() = α, where γ t is the barycentric projection of the optimal plan between α t and χ. In other words, the curve α is globally Lipschitz and satisfies the system (2), (3), (9), (5) in the sense of distributions in (, ) IR 3. Also, the Hamiltonian energy t H(α t ) is conserved. Remark 2.4. Since [, ) t α t P 2 (IR 3 ) is continuous, we use Proposition 3.2 [19] to conclude that there is a family P C([, ); H 1 ()) of convex functions P (t, ) such that x P (t, ) =: P t pushes χ forward to α t optimally for all t. 3 Lagrangian vs Eulerian As of this date, the only weak Lagrangian solutions in the physical space have been shown to exist in the case P # χ =: α L 1 (). This is achieved in [14] by improving the L q (q > 1) result in [9]. The solution (P, F ) constructed in these references satisfies P t# χ =: α t L q () for almost all t (, T ). We work with solutions that may not satisfy any of these conditions, i.e. the measures α t (for t [, T )) may be singular. In particular, α t can be Dirac measures δ z(t) as in the example given in [14] (see also Section 4.1 below). An interesting question is existence of Eulerian solutions in physical space. Recently, existence of such solutions for a class of initial data, where the conditions include the requirement that the support of the function α = P # χ in the dual space be the whole space, was obtained by L. Ambrosio, M. Colombo, G. De Philippis, A. Figalli in two-dimensional periodic case [2] and in three-dimensional case with = IR 3 [3], based on the results of G. De Philippis and A. Figalli [13] on the regularity of solutions for the Monge-Ampere equation. Existence of Eulerian solutions in physical space for more general initial data α L q (), when the support of α in the dual space may have a non-empty boundary, or if α is a singular measure, is presently not known. In this section, we show that Eulerian solutions in the physical space exhibiting some mild regularity in time cannot give rise to very irregular solutions in dual space. Let us begin by recalling the definition of a weak solution of (1). Definition 3.1. Let u : [, T ) IR 3 and P : [, T ) IR satisfy u L 1 ([, T ) ; IR 3 ), P L ([, T ) ) C([, T ); L 1 ()), and P (t, ) is convex in for every t [, ). The pair (P, u) is a weak Eulerian solution of (1) if T { P (t, x) [ t φ(t, x) + φ(t, x)u(t, x)] (22) +J[ P (t, x) x] φ(t, x)}dxdt + P (x) φ(, x)dx = for any φ Cc 1 ([, T ) ; IR 3 ), and for any ψ C 1 c ([, T ) ). T u(t, x) ψ(t, x)dxdt = (23) In what follows, we show that one cannot have weak Eulerian solutions exhibiting mild time regularity and spatial flat parts except, possibly, at negligibly many times. We begin with
8 8 Proposition 3.3 below, whose proof is a relatively straight-forward adaptation of Corollary 2.3 and Proposition 2.4 in [19]. The only difference is that in the said reference the set O is a subset of (, T ) of full measure. Before stating the result, we need some preliminary observations. If X H 1 (, T ; L 2 (; IR 3 )), we denote by Ẋ L2 (, T ; L 2 (; IR 3 )) its functional derivative, defined by lim X t+h X t h h Ẋt = for L 1 a.e. t (, T ). L 2 (;IR 3 ) In the next lemma, we extend X to a map in AC 2 (IR; L 2 (; IR 3 )) by setting X t = X + for t and X t = X T for t T. Lemma 3.2. Let X H 1 (, T ; L 2 (; IR 3 )) and T lim h Ẋ be its functional derivative. Then X t+h x X t x h Ẋtx 2 dxdt =. (24) As a consequence, there exist sequences h + k +, h k and a measurable subset A IR such that L 4 ((IR ) \ A) = and for all (t, x) A. lim k X t+h + x X t x k h + = lim k k X t+h x X t x k h = Ẋtx (25) k The proof in [15] needs no modification. The philosophy behind this result is that, in some specified sense, Ẋ can be viewed as almost a classical pointwise time-derivative of X. Also, since X H 1 (, T ; L 2 (; IR 3 )), we have that it admits a Borel representative. Equation (24) shows that Ẋ itself has that property. Throughout the paper we identify both X and Ẋ with their Borel representatives. Proposition 3.3. Let X H 1 (, T ; L 2 (; IR 3 )) be such that X t = P t, where P t is convex for all t (, T ). Let A (, T ) as in Lemma 3.2. Furthermore, let O A be a Borel set with L 4 (O) > and such that L 3 ([X t x]) > for all (t, x) O, where [X t x] := {y : X t y = X t x}. Then, there exists a Borel map w : (, T ) IR 3 IR 3 such that Ẋ(t, x) = w(t, X(t, x)) for L 4 a.e. (t, x) O. (26) Proof: ( Let ) λ denote the L 4 measure restricted to O, Ψ : O (, T ) IR 3 given by Ψ(t, x) = t, X(t, x), and set ϑ := Ψ# λ. Denote by η the vector-measure whose density with respect to λ is Ẋ, then set σ := Ψ # η. Clearly, σ ϑ, which means there exists a Borel vector field w : (, T ) IR 3 IR 3 such that dσ = wdϑ. The disintegration theorem (see, for example, Theorem [5]) applies to the Borel vector field Ψ and the measure λ. Thus, for ϑ a.e. (t, y) (, T ) IR 3, there exists a unique Borel probability measure λ t,y on O such ( that the map (t, y) λ t,y (B) is Borel measurable for each Borel set B O. Furthermore, λ t,y Ψ 1 (t, y) ) = 1 for ϑ a.e. (t, y) (, T ) IR 3 and O f(t, x) dλ(t, x) = T IR 3 ( ) f(t, x) λ t,y (dt, dx) ϑ(dt, dy) Ψ 1 (t,y)
9 9 for every Borel measurable f : (, T ) [, ]. We showed in [19], Theorem 2.2, that w(t, X(t, x)) = Ẋ(t, z) dλ t,x(t,x) (t, z) for λ a.e. (t, x) O. (27) Ψ 1 (t,x(t,x)) Note that (t, z) Ψ 1 (t, X(t, x)) is equivalent to (t, z) O and X(t, z) = X(t, x), so we have Ψ 1 (t, X(t, x)) = O [X t x]. Then we apply Proposition 2.4 in [19] to get that Ẋ(t, x) = Ẋ(t, z) for all (t, z) O [X t x]. According to (27), we get(26). QED. Before proving the main theorem of this section, we need a measurability lemma. Lemma 3.4. Let α AC 2 (, T ; P 2 (IR 3 )) for some T >. Then the set D(α) := { (t, X) (, T ) IR 3 : X is an atom of α t } (28) is Borel. Proof: Denote by C c + the nonnegative cone of C c (IR 3 ). Consider, for every positive integer m and every ξ C c +, the set { Dm ξ := (t, X) (, T ) IR 3 : ξ(y )dα t (Y ) 1 } IR 3 m ξ(x). Note that the absolute continuity in time of the left hand side of the above inequality [5] and the continuity in X of the right hand side imply that the difference is a continuous function of (t, X). Therefore, Dm ξ is the nonnegative set of a continuous function, which makes it a closed subset of (, T ) IR 3. Thus, { D m := (t, X) (, T ) IR 3 : ξ(y )dα t (Y ) 1 } IR 3 m ξ(x) for all ξ C+ c is closed as well, by being an arbitrary intersection of closed sets. Since D(α) = m 1 D m, the proof is concluded. QED. Now we prove the main result of this section: that weak (Eulerian) solutions in the physical space exhibiting some mild regularity in time cannot give rise to very irregular solutions in dual space. Theorem 3.5. Let (P, u) be a weak solution for the Semi-Geostrophic system in the physical space such that P H 1 (, T ; L 2 (; IR 3 )). Then α t := P t# χ is atom-free for L 1 a.e. t (, T ). Proof: Set Ψ(t, x) := (t, X(t, x)). Since Ψ is a Borel map on (, T ), due to Lemma 3.4 we infer that O := Ψ 1( D(α) ) is a Borel subset of (, T ). One can see that O = {(t, x) (, T ) : L 3( [X t x] ) > }.
10 1 By Fubini s Theorem and by the convexity of the potentials whose gradients push χ forward to α t, we infer that there exists a Borel subset T (, T ) such that O = t T ( {t} Ot ), (29) where O t is the union of all (at most countably many) convex subsets of of positive L 3 measure on each of which X t is constant. Assume by contradiction that L 1 (T ) >. We further throw out of T the L 1 negligible set of times at which ζ(x)u(t, x)dx IR 3 (3) for all ζ Cc (; IR 3 ) (via the separability of this space with respect to the sup norm), but we keep the notation T for the remaining subset. Consider A as the Borel subset of (, T ) of full measure defined in Lemma 3.2, i.e. where the time pseudo-derivative of X in the sense of (25) exists. The set A O is a Borel set (which we still denote by O) with L 4 (O) >. According to Proposition 3.3, we infer there exists a Borel map w : (, T ) IR 3 IR 3 such that (26) is satisfied. By taking φ(t, x) = ξ(t)ζ(x) with ξ Cc (, T ) and ζ Cc (; IR 3 ) in (22) we discover that t X(t, x) ζ(x)dx is absolutely continuous and d X(t, x) ζ(x)dx = {X(t, x) [ ζ(x)u(t, x)] + J[X(t, x) x] ζ(x)}dx dt for a.e. t (, T ) and every ζ Cc (; IR 3 ) (via the usual argument involving the separability of this space endowed with the sup norm). Throwing out, if necessary, a negligible set of times, we conclude Ẋ(t, x) ζ(x)dx = {X(t, x) [ ζ(x)u(t, x)] + J[X(t, x) x] ζ(x)}dx for a.e. t (, T ) and every ζ Cc (; IR 3 ). Choose such a t that also lies in T and consider now only test functions ζ Cc (ω ; IR 3 ), where ω is a connected component of O t (which is a convex subset of of positive L 3 measure on which X t is constant). Since X t c IR 3 in ω, we infer w(t, c) ζ(x)dx = c ω ζ(x)u(t, x)dx + ω J[c x] ζ(x)}dx ω for all ζ Cc (ω ; IR 3 ). Due to (3) the first term in the right hand side vanishes and since the equality holds for all ζ Cc (ω ; IR 3 ), we deduce w(t, c) = J[c x] for a.e. x ω, which contradicts the fact that ω has nonempty interior. QED. Remark 3.6. Thus, in order to accommodate singular solutions in dual space, we see the need for defining Lagrangian solutions instead of Eulerian ones in the physical space. Whereas solutions in the dual space may come in any form or shape (from pure Dirac deltas to functions), only the absolutely continuous ones with respect to the Lebesgue measures have been so far known to give rise to Lagrangian solutions in the physical space [9], [14]. In the next section we discuss an extension to this notion and prove some existence results.
11 11 4 Weak Lagrangian solutions in physical space for the case of singular initial data In [9], Lagrangian solutions in the physical space with initial data P # χ = α L p (IR 3 ) (p > 1), were constructed by the following procedure. First, time-stepping approximation from [1] combined with results of [1] yield existence of a solution (α, P ) of the dual space system (2) (4) with initial data (5), and a locally bounded map Φ : (, T ) IR 3 IR 3 satisfying such that Φ = U(, Φ), L 4 a.e. in (, T ) P (), Φ(, X) = X for α a.e. X IR 3, (31) α t = Φ t# α for all t >. There also exists a Borel map Φ such that Φ t preserve L 3 and Φ t Φ t = Id = Φ t Φ t a.e. in IR 3. The physical space flow is defined as F t := P t Φ t P. (32) Then it is shown that (P, F ) satisfies all the requirements of Definition 1.1. In order to see what can be expected in the case of general initial data, we consider an example in which α is a point mass. 4.1 An example: a limit of Lagrangian solutions, with initial measures weakly converging to a point mass We discuss the case α = 4 3 πδ z for some z IR 3 by considering a sequence of approximations α (ε) = 1 ε 3 χ B z (ε), where ε >, where χ A denotes the indicator function of the set A. This is a version of the counterexample in [14]. For notational simplicity, we drop the requirements that α t be probability measures and that L 3 () = 1; thus, α t are measures of total mass equal to L 3 (). We fix the domain in physical space to be the ball = B 1 () IR 3, which means χ = L 3 B1 (). It is easy to see that α = P # χ and α (ε) = P (ε) χ, where # P (x) = z x, P (ε) (x) = z x + ε x 2 2. (33) Following the calculations in [14], we find that for each ε >, one Lagrangian solution (P (ε), F (ε) ) in the sense of Definition 1.1 with initial data P (ε), can be obtained as follows. Let z(t) be determined by solving ż(t) = Jz(t), z() = z ; thus, z(t) = L 1 (t)z, where cos ct sin ct L c (t) = sin ct cos ct for c IR. (34) 1 Then P (ε) t (x) = z(t) x + ε x 2 2, F (ε) t (x) = L c (t)x, with c = 1 1 ε. (35)
12 12 We also note that the Legendre transform (over ) of P (ε) t is X z(t) 2 (P (ε) t ), if X B (X) = 2ε z(t) (ε); X z(t) ε 2, if X / B z(t)(ε), and the flow map in the dual variables is Φ (ε) t (X) = z(t) + M c (t)(x z ), with c = 1 1 ε, where M c = cos ct sin ct sin ct cos ct. (36) Expression (35) was obtained by (32) using (P (ε) t ), Φ (ε) t given above. Now we take limits as ε +. Clearly, where P (ε) t P t uniformly in for all t, (37) P t (x) = z(t) x for x, t. (38) Note that by taking t = in (38) we obtain the function P in (33). On the other hand, as far as F (ε) t is concerned, we see that F (ε) (x) = x for all x, ε >, yet, for any t > the sequence F (ε) t ( ) does not converge even in a weak sense. Then we consider a limit of F (ε) as functions of (t, x), and we note that, as ε +, F (ε) weakly in L 2 ([, T ] ; IR 3 ) for any T >. (39) (39) determines the limit F (ε) only up to the set of points (t, x) of L 4 -measure zero. However, since the initial measure α for the limiting problem is a point mass, it is natural to expect that the flow map for such problem is sensitive to changes on sets of L 4 -measure zero. In particular, we cannot take the flow map of the form F for the limiting problem, since we require F (, ) Id. Thus, in the following lemma we note some properties which are satisfied by any map ˆF : [, T ) which is a weak limit of F (ε). We first note the following: for P t in (38) we have α t := P t# χ = 4 3 πδ z(t), and the Kantorovich optimal plan γ t := ( P t Id ) # χ is defined by γ t(a B) = δ z(t) (A)χ(B) for all Borel A IR 3, B. Then (1) implies that the barycentric projection of γ t onto α t is which defines γ t (X) for α t -a.e. X IR 3. γ t (z(t)) =, (4)
13 13 Lemma 4.1. Let T > and ˆF : [, T ) be a Borel map such that F (ε) ˆF weakly in L 2 ([, T ] ; IR 3 ) for any T > as ɛ +, where F (ε) is defined by (35). Then i. ˆFt# χ = 4 3 πδ for a.e. t >. In particular, the map ˆF t does not preserve L 3 for a.e. t. ii. (P, ˆF ) satisfies property (14) for a.e. t >. Here P is defined by (38). Proof: Assertion (i) follows from (39) and Fubini Theorem. To show (ii), we fix any t > and note that (38), (4) imply γ t ( P t (x)) = γ t (z(t)) = for every x, which implies ( γ t P t ) # χ = 4 3 πδ. Now (ii) follows from (i). QED. Remark 4.2. Lemma 4.1 suggests that for the case of singular initial measure α, the L 3 - measure preserving property of F t# χ = χ of F t should be replaced by (14) if we expect to have stability of weak Lagrangian solutions in physical space. Now we discuss a possible construction of a Lagrangian solution in physical space for the limiting initial data P (x) = z x, i.e. with α = 4 3 πδ z, and compare it with the limit of (P (ε), F (ε) ) obtained above. We use P t in (38) to see that α t := P t# χ = 4 3 πδ z(t). One can readily check that the property ż(t) = Jz(t), z() = z implies that α is a weak solution of the SG in dual space, i.e. of the problem (2), (3), (9), (5). In order to construct a flow map F t in in physical space, let us go back to the definition of F t = Pt Φ t P from [9]. Since in the general case (when α t is a measure) the place of Pt is taken by the barycentric projection γ t (see (13)), it is natural to ask if one can obtain a solution by putting F Id, F t := γ t Φ t P for t >. (41) Of course, for this function to be well-defined, one needs the Lagrangian flow in dual space Φ t to be defined at z for all t [, T ]. Thus, we solve for X = z. Using (4) we rewrite this as Φ(t, X) = J[Φ(t, X) γ t (Φ(t, X))] with Φ(, X) = X (42) Φ(t, z ) = JΦ(t, z ) with Φ(, z ) = z. Since ż(t) = Jz(t), z() = z, then uniqueness for the Cauchy problem associated to this ODE system implies Φ(t, z ) = z(t). Now, F t is well-defined by (41), and explicitly is F (t, x) = F t (x) = { x, if t =, x,, if t >, x. (43)
14 14 From (39), (43) we conclude that, as ε +, We also notice that F (ε) F weakly in L 2 ([, T ] ; IR 3 ), for any T >, F (ε) F Id in, for any ε >. Moreover, we showed in the proof of Lemma 4.1 that (44) F t# χ = 4 3 πδ for all t >. (45) (γ t P t ) # χ = 4 3 πδ for t >. Then (45) implies that F t satisfies the generalized measure-preserving property (14) for all t >. Furthermore, since ż(t) = Jz(t), z() = z, and Z(t, x) := ( P t )(F (t, x)) = z(t), it follows that (P, F ) satisfies equation (7) in the sense (8). Consequently, we have arrived to: Lemma 4.3. Let P be defined by (38). Then F defined by (41) has the explicit form (43), and (P, F ) satisfies the following: (a) (P, F ) is a weak limit (in the sense of (37), (44)) as ε +, of a family of Lagrangian solutions (P (ε), F (ε) ) in the sense of Definition 1.1; (b) (P, F ) satisfies equation (7) in the sense (8); (c) For all t >, (P, F ) satisfies property (14), which takes the form (45). In other words, the limit (P, F ) satisfies the properties of Lagrangian solution (as in Definition 1.1) except for the L 3 -measure preserving property of F t, which is replaced by (14). An important remaining issue, the continuity of t F t ( ), will be clarified in the next section. Remark 4.4. The above computations of the Lagrangian flow in dual space highlight the difficulty in applying the theory of regular Lagrangian flows [1] to the case of singular initial measures α. One can extend γ t, defined by (4), to all IR 3 by using that P t (X) = X z(t) is differentiable for X z(t), defining γ t (X) := P t (X) = X z(t) X z(t) if X z(t). Then one can easily check that U(t, X) := J[X γ t (X)] satisfies the conditions in [1], therefore a unique solution to the ODE (42) exists and is unique for Lebesgue-a.e. X IR 3 such that Φ t preserves L 3. However, when α is singular, the set {X IR 3 : Φ(t, X) exists and is unique} may not contain all (if any) of the support of α. In the case α = δ z, the solution Φ(, z ) can be found explicitly by solving the ODE system, and this defines Φ(, X) for α -a.e. X IR. In general, for singular measures α, it is not clear how to find a suitable Lagrangian flow Φ(t, X) in the dual space that is defined for α -a.e X. 4.2 Definition of Lagrangian solutions with singular initial data, and existence results In order to see the motivation for Definition 4.12 of Lagrangian solutions for the case of singular initial data given below, we first study properties of flows in physical space given by (41) in the
15 15 case when α = P # χ is possibly a singular measure, under the following assumptions: (i) P H 1 () is convex, α = P # χ P 2 (IR 3 ); (ii) There exists a Borel map Φ : [, T ] IR 3 IR 3 such that (42) holds in the integral sense for α -a.e. X IR 3 ; (iii) The family of measures α t := Φ t# α for t [, T ) is a solution in D ((, T ) IR 3 ) for the equation (2) with U defined by (9), (1), where γ t is the unique optimal transport plan between α t and χ. This also defines a convex P t in (for each t) such that γ t := ( P t Id ) # χ; (iv) α AC (, ; P 2 (IR 3 )). Note that assumptions (46) say that there exists a Lagrangian solution of the dual problem with initial data P. In the case when P W 1, () and α = P # χ L p (IR 3 ) for p (1, ] existence of Φ t such that (46) is satisfied is shown in [9]. Below it is convenient to work with a Borel representative of the barycentric projection γ. That is why we prove the following lemma: Lemma 4.5. There exists a Borel measurable function defined on (, T ) IR 3 which for L 1 -a.e. t (, T ) coincides with γ(t, X) for α(t, )-a.e. X IR 3. Proof: Let ϑ P ( (, T ) IR 3 ) be the Borel probability given by T IR 3 ϕ(t, X, y)ϑ(dt, dx, dy) = 1 T T IR 3 ϕ(t, X, y)γ(t, dx, dy)dt for any continuous and bounded ϕ. Since α(t, ) is the X-marginal of γ(t,, ), the (t, X) marginal of ϑ is ϑ P ( (, T ) IR 3) given by T Thus, by disintegrating ϑ we get T IR 3 IR 3 ζ(t, X) ϑ(dt, dx) = 1 T ϕ(t, X, y)ϑ(dt, dx, dy) = 1 T T T IR 3 IR 3 ζ(t, X)α(t, dx)dt. ( ) ϕ(t, X, y)ϑ(t, X; dy) α(t, dx)dt, where ϑ(t, X; ) are Borel probabilities on such that the map (t, X) ϑ(t, X; B) is Borel for any Borel set B. In particular, the maps (t, X) f(y)ϑ(t, X; dy) ( are Borel for all f C b ; IR 3 ). By taking ϕ(t, X, y) = u(t)ξ(x) y and using (1), we conclude that γ(t, X) = yϑ(t, X; dy) for L 1 a.e. t (, T ) and α(t, ) a.e. X IR 3. This finishes the proof. (46) QED.
16 16 Remark 4.6. In light of Lemma 4.5 and equation (9), we see immediately that there exists a Borel measurable function defined on (, T ) IR 3 which for L 1 -a.e. t (, T ) coincides with U(t, X) for α(t, )-a.e. X IR 3. Remark 4.7. In view of the above lemma, from now on we shall use the notation γ to denote the Borel representative of the barycentric projection defined by (1). Likewise, by following the same proof, it is easy to see that the barycentric projection of ( Id P t χ onto its first )# marginal (namely, χ) can also be extended to a Borel map from (, T ) into IR 3. Since this coincides with P t (x) for L 4 -a.e. (t, x) (, T ), we shall assume in the remainder of the paper that P is Borel measurable in both variables. We also note the following: Lemma 4.8. For (t, X) [, T ) IR 3 denote t,x := {x : P t (x) exists and P t (x) = X}. Then In particular, γ t (X) t,x for α(t, )-a.e. X IR 3. (47) ( P t γ t )(X) = X for every t [, T ) and α(t, )-a.e. X IR 3. (48) Proof: Fix t [, T ). Since γ t = ( P t Id) # χ and π 1# γ t = α t, by disintegrating γ t as in Theorem [5], we get γ t = γ t,x dα t (X), IR 3 where γ t,x is a family of probability measures on such that the map X γ t,x (B) is Borel for any Borel set B, and γ t,x (IR 3 \ t,x ) = for α(t, )-a.e. X IR 3. From (1) and the disintegration, we get γ t (X) = y dγ t,x (y) = y dγ t,x (y) for α(t, )-a.e. X IR 3. t,x The convexity of P t ( ) implies that t,x is a convex set (of dimension k(t, X) {, 1, 2, 3}) for every X such that t,x. Thus, γ t (X) t,x for α(t, )-a.e. X IR 3, and for such X we get ( P t γ t )(X) = X. QED. The continuity of F in time was proved in [9], [14] in the case of absolutely continuous α t. In general we cannot expect that, which is clear from the structure (41) of F t, especially looking at the example when P is linear, i.e. α = δ z discussed above: streamlines in the physical space are concentrated on the barycenters γ t (X) of the sets Pt (X) where X supp (α t ). Thus it is natural to expect continuity of F t relative to P t. Indeed, we have the following:
17 17 Proposition 4.9. Assume that (46) hold. Then the map F defined in (41) satisfies lim ξ( P t F t (x)) F t (x)dx = ξ( P t F t (x)) F t (x)dx (49) t t for all t [, ) (which we define as meaning t + if t = ) and all ξ C c (IR 3 ; IR 3 ). and Furthermore, lim ξ( P t F t (x)) F t (x)dx = t t for all t [, ) and all ξ C c (IR 3 ; IR 3 ). ξ( P t F t (x)) F t (x)dx, (5) lim ξ( P (x)) F t (x)dx = ξ( P (x)) x dx (51) t + Proof: The (unique) optimal transport plan between its X-marginal α t and it y-marginal χ is γ t = ( P t Id) # χ. Let t [, ). Since W 2 (α t, α t ) as t t (which we define as meaning t + if t = ), by stability of optimal plans [5], we infer that γ t converges to γ t in P 2 (IR 3 ). We use the definition (41) of F t, the fact that P # χ = α and Φ t# α = α t, and (1) to get for t > ξ(φ t P (x)) F t (x)dx = ξ(x) yγ t (dx, dy). IR 3 IR 3 For t =, we use F = Id, Φ = Id, P # χ = α and (1) to get the equality displayed above. Now, by the continuity of t γ t proved above, we deduce for t [, ) lim ξ(φ t P (x)) F t (x)dx = ξ(x) yγ t (dx, dy) = ξ(φ t P (x)) F t (x)dx. (52) t t IR 3 Next we note that (Φ t P ) # χ = α t. Combining this with (48), we obtain for every t (, T ): P t γ t Φ t P (x) = Φ t P (x) for a.e. x. Then using (41) we obtain for every t (, T ) P t F t (x) = Φ t P (x) for a.e. x. (53) For t =, we deduce (53) from F = Id, Φ = Id. Now (52) implies (49). Furthermore, (31) implies This yields t Φ t (X) = Φ t (X) + U(s, Φ s (X))ds for α a.e. X IR 3 and all t t. t ξ(φ t (X)) ξ(φ t (X)) 2 ξ 2 ( t (t t ) ξ 2 t U(s, Φ s (X)) ds t ) 2 t U(s, Φ s (X)) 2 ds
18 18 for α -a.e. X IR 3. Note that (, T ) IR 3 (t, X) U(t, Φ t (X)) is Borel (as composition of Borel maps). Next we use P # χ = α and Φ s# α = α s to get ξ Φ t P ξ Φ t P L 2 (;IR 3 ) ( t ) 1/2 t t ξ U(s, Y ) 2 dα s (Y )ds. IR 3 Then, using (53) we obtain ξ P t F t ξ P t F t L 2 (;IR 3 ) ( t ) 1/2 t t ξ U(s, Y ) 2 dα s (Y )ds. IR 3 Since U t L 2 (α t;ir 3 ) L (, ) and F is bounded uniformly in time-space, we use (49) and the inequality displayed above to deduce (5). For t = we use the fact that F (x) = x for a.e. x to get (51) from (5). QED. t t Remark 4.1. In the case α t L 3 and 1 < p <, note that the strong L p -continuity in time F t F t L p () as t t is equivalent to weak L p -continuity, since all the maps F t share the same L p -norm (by being measure-preserving). In fact, we can further prove that when α t L 3 for all t, (5) is equivalent to this weak L p -continuity. Thus, our result above simply generalizes the L p -continuity obtained in [9]. Proof: Let 1 < q < be such that 1/p + 1/q = 1. First, one can easily check that ζ L q (α t ; IR 3 ) ζ Φ t P L q (; IR 3 ), ξ L q (; IR 3 ) ξ P Φ t L q (α t ; IR 3 ). (54) Thus, pick ξ L q (; IR 3 ) and a sequence {ξ n } n Cc 1 (IR 3 ; IR 3 ) which converges to ξ P Φ t strongly in L q (α t ; IR 3 ), i.e. ξ n Φ t P converges strongly to ξ in L q (; IR 3 ). We have ξ (F t F t )dx ξ n Φ t P (F t F t )dx + C(p, ) ξ n Φ t P ξ L q (;IR 3 ). We choose n sufficiently large such that the last term in the right hand side is sufficiently small. Then use (5) for ξ = ξ n to conclude that (5) implies the weak L p -continuity for t F t at t. A similar argument can be used to prove the converse. Proposition 4.9 motivates the following QED. Definition Let P : [, ) IR be such that P C([, ); H 1 ()) and P (t, ) is convex in for each t [, ). A map F : [, ) is called weakly P continuous if (5), (51) hold for all t (, ) and all ξ C c (IR 3 ; IR 3 ). Furthermore, we note that considering the example of initial data in the dual space being a linear combination of Dirac masses (see also Proposition 4.14 below), the solution α t at each time t is also a linear combination of Dirac masses concentrated in the time-dependent location, it is not clear whether distinct initial location of Dirac masses should imply that their locations are distinct at all times. Thus we cannot expect existence of the map F as in Definition 1.1(iii). In light of the above, we are ready to generalize Definition 1.1 as follows:
19 19 Definition Let P : [, ) IR be such that P C([, ); H 1 ()) and P (t, ) is convex in for each t [, ). Let F : [, ) be a weakly P -continuous Borel map. Denote by γ t the barycentric projection of the measure γ t := ( P t Id ) χ, defined in (1), and # set α t := P t# χ, µ t := γ t# α t. (55) Then the pair (P, F ) is called a weak Lagrangian solution of (1) in [, T ) if i. F (, x) = x and P (, x) = P (x) for all x ; ii. for any t > the mapping F t = F (t, ) : satisfies F t# χ = µ t and F t# µ = µ t ; iii. The function Z : (, T ) IR 3 defined by Z t = P t F t (56) lies, along with F, in L (, T ; L 2 (; IR 3 )) and is a distributional solution of (7) in the sense of (8). The following proposition gives sufficient conditions for the existence of weak Lagrangian solutions in physical space. Note that these conditions are exactly what was proved in [9] in the case α t L q (IR 3 ). Proposition Assume that (46) hold. Define F by (41). Then the pair (P, F ) is a weak Lagrangian solution in physical space in the sense of Definition Furthermore, if there exists a Borel map Φ : [, T ] IR 3 IR 3 such that Φ # α t = α for all t [, T ), then there exists a Borel mapping F : [, T ] satisfying F t# µ t = µ and F t F t = Id µ -a.e., and F t F t = Id µ t -a.e. (57) Proof: Proposition 4.9 implies that F is weakly P -continuous. From (48) and since P t# χ = α t, we have: P t γ t P t (x) = P t (x) for every t (, T ) and χ a.e. x. (58) Now, in order to prove (i) of Definition 4.12 for F, we find that γ (Φ P γ ) P (x) = γ P (x) for χ-a.e. x, due to (58) for t = and the hypothesis on Φ. Since µ = γ P # χ, then (i) follows. To prove (ii) we check that for t > F t# χ = ( γ t Φ t P ) # χ = ( γ t Φ t ) # α = γ t# α t = µ t, where we used the definitions (55), and that Φ t# α = α t. Similarly, for t > F t# µ = ( γ t Φ t P ) # µ = ( γ t Φ t P γ ) # α = ( γ t Φ t ) # α = γ t# α t = µ t,
20 2 where we used in the second line that ( P γ ) # α = α by (48) for t =. Next we prove (iii) of Definition We first note that Z and F are Borel (therefore, Lebesgue) measurable as compositions of Borel maps (see Lemma 4.5). By the definition of γ, we have ( ) 1/2 ( ) 1/2 ξ(x) γ t (X)dα t (X) ξ(x) IR IR 2 dγ t (X, y) y 2 dγ t (X, y) 3 3 ( 1/2 = ξ L 2 (α t;ir 3 ) y dy) 2. It follows that for all t (, T ) we have γ t L 2 (α t ; IR 3 ) with γ t L 2 (α t;ir 3 ) R (L 3 ()) 1/2, (59) where < R < is large enough such that B(, R ). Thus, there exists C IR independent of t such that F t (x) 2 dx = γ t Φ t 2 dα = IR 3 γ t 2 dα t = γ t 2 IR 3 L 2 (α t;ir 3 ) C. Also, using that F t# χ = µ t as proved above, and also using (48), we have for Z t = P t F t : Z t (x) 2 dx = P t (y) 2 dµ t (y) = P t ( γ t (X)) 2 dα t (X) = X 2 dα t (X) = P t (x) 2 dx C, IR 3 where C < is a constant coming from the fact that P C([, T ]; H 1 ()). Thus, both Z and F belong to L (, T ; L 2 (; IR 3 )). Next we note that assumption (ii) of (46) implies that for for χ-a.e. x, the function Φ(, P (x)) is a weak solution of the problem Φ(t, P (x)) = J[Φ(t, P (x)) γ t (Φ(t, P (x)))] with Φ(, P (x)) = P (x). (6) From (53), for each t [, T ) the equality Z t (x) = Φ t P holds for a.e. x. Thus using the integrability of Z and F proved above, using a function ϕ C 1 ([, T ) IR3 ) in the weak form of (6) and integrating with respect to x, we get (8). Finally, in order to prove (57) under the additional assumption of the existence of Φ, we set F (x) = x, F t (x) := γ Φ t P t for t >. Then property (57) is obvious for t =. Thus we fix t > and and compute F t# µ t = ( γ Φ t P t γ t ) # α t = ( γ Φ t ) # α t = γ # α = µ. (61) To finish proving (57), note that F t F t (y) = y for µ -a.e. y amounts to for χ-a.e. x. Since F t F t γ P (x) = γ P (x) (62) (F t F t γ )( P (x)) = ( γ Φ t P t γ t Φ t P γ )( P (x))
21 21 and P t γ t, P γ, Φ t Φ t are all equal to the identity on the corresponding domains, we deduce (62). QED. In the case of discrete measures we can prove that the construction works. Even though not explicitly present in [4], a simple argument inserted in the proof of Theorem 7.4 from said reference yields that the solutions to the Hamiltonian ODE constructed there are convex combinations of point masses provided that the initial measures are of the same form (coefficients of the convex combinations are time-invariant). Our Hamiltonian in (21) satisfies all the requirements for Theorem 7.4 [4] to apply (see Lemma 7.6 [4]). Proposition Let x = (x 1,..., xn ) IR3n be arbitrary for some integer n 1, and let µ = n c i δ x i (63) i=1 be a convex combination of the Diracs at these points, i.e. nonnegative constants c i satisfy n i=1 c i = 1. Then the solution constructed in Theorem 7.4 [4] for the initial-value problem associated to the Hamiltonian ODE as in Definition 2.1 for the Hamiltonian in (21) is of the form n α t = c i δ x i (t), (64) where [, T ] t x(t) is in W 1, (, T ; IR 3n ) and x() = x. i=1 Proof: Let m be a positive integer and set h = T/m. Then take w m := J H(ᾱ), where H(α) denotes the element of H(α) with least L 2 (α; IR 3 )-norm, and set, for all t [, h], α m t = (Id + tw m ) # ᾱ, ν m t = (Id + tw m ) # (w m ᾱ) and w m t := dνm t dα m t where we used the fact (see Lemma 7.1 [4]) that νt m αt m to get the Radon-Nykodim derivative wt m. On the next time subinterval [h, 2h] one defines αt m and wt m similarly by using αh m and wh m instead of ᾱ and wm, and t h instead of t. In general, the construction can be extended to [kh, (k + 1)h] for k =,..., m 1 by using αkh m and wm kh, t kh instead of t, and repeating the steps above. It is proved in [4] that the paths of measures t αt m are uniformly bounded in P 2 (IR 3 ) and uniformly Lipschitz continuous. For a subsequence m j we have a limiting t α t, which is shown to satisfy the Hamiltonian ODE with α = ᾱ. Since α m j t = [Id + (t kh)w m j kh ] #α m j kh for t [kh, (k + 1)h], we deduce that all probabilities α m j t are convex combinations of n Dirac masses if α m j kh is (with same coefficients). This is true for all k =,..., m j 1, so we deduce that it holds for α m j t all t [, T ]. The uniform bounds on αt m mentioned above translate into the uniform L (, T ; IR 3n ) bounds on t x m j (t) = (x mj,1 (t),..., x mj,n (t)) (where x mj,i (t) are the points in the support of α m j t ). Furthermore, the uniform Lipschitz continuity of the paths t αt m in the Wasserstein space P 2 (IR 3 ) gives a finite constant C > for which W 2 (α m j t, α m j s ) C t s for all t, s [, T ]. Now fix t (, T ). Note that x m are piecewise linear and continuous in time, thus it is clear that C 2 t t 2 W2 2 (α m j t, α m j t ) = n c i x mj,i (t) x mj,i (t ) 2 for all t close enough to t. i=1,
22 22 Thus, the vector functions t x m j (t) are uniformly Lipschitz. By Ascoli-Arzela s theorem, a subsequence converges in the sup norm to a function x W 1, (, T ; IR 3n ), which implies the limiting measures α t found above must have the form (64) for all t [, T ]. This follows from W2 2 (α m j t, β t ) n c i x mj,i (t) x i (t) 2 as j, i=1 where β t is the convex combination of the Dirac masses at x i (t) with coefficients c 1,..., c n. QED. So far we are not aware of any reason why distinct initial x k should give rise to distinct xk (t) at all later times. As a consequence, existence of the map Φ is uncertain (since transport maps from an average of n points masses to one of m point masses exist if and only if n m). Proposition Let n be a positive integer and α t := n c i δ x i (t), for t [, T ) i=1 be the solution of SG in dual space constructed in Proposition 4.14 with initial data α := n c i δ x i, i=1 where x i are arbitrary in IR3, i = 1,..., n, and c i are nonnegative and n i=1 c i = 1. Then the map Φ n : [, T ] IR 3 IR 3 given by Φ n (t, X) = X if X x i, and Φn (t, x i ) = xi (t) for t [, T ], i = 1,..., n satisfies (46), (ii) (iii). Proof: It is obvious that Φ n is Borel and Φ n satisfies Φ n t# α = α t. To show that Φ n solves (42) we start from the fact that α t solves the system (17), (18). In fact, the Hamiltonian (21) has the property that T α P 2 (IR 3 ) H(α) = { γ Id } (see, e.g. [4]), so α t solves (13) in the sense of distributions. This is equivalent to t ξ(x)dα t (X) is absolutely continuous IR 3 and for a.e. t (, T ) we have d ξ(x)dα t (X) = ξ(x) U(t, X)dα t (X) for all ξ Cc 1 (IR 3 ), dt IR 3 IR 3 for U given in (9). According to Proposition 4.14, we have that t x i (t) is in W 1, (, T ; IR 3 ) for all i = 1,..., n, which implies n c i ξ(x i (t)) ẋ i (t) = i=1 n c i ξ(x i (t)) U(t, x i (t)) for all ξ Cc 1 (IR 3 ). i=1 This leads to the desired conclusion. QED. We are now in a position to formulate:
On Lagrangian solutions for the semi-geostrophic system with singular initial data
On Lagrangian solutions for the semi-geostrophic system with singular initial data Mikhail Feldman Department of Mathematics University of Wisconsin-Madison Madison, WI 5376, USA feldman@math.wisc.edu
More informationThe semi-geostrophic equations - a model for large-scale atmospheric flows
The semi-geostrophic equations - a model for large-scale atmospheric flows Beatrice Pelloni, University of Reading with M. Cullen (Met Office), D. Gilbert, T. Kuna INI - MFE Dec 2013 Introduction - Motivation
More informationSobolev regularity for the Monge-Ampère equation, with application to the semigeostrophic equations
Sobolev regularity for the Monge-Ampère equation, with application to the semigeostrophic equations Alessio Figalli Abstract In this note we review some recent results on the Sobolev regularity of solutions
More informationCLASSICAL SOLUTIONS TO SEMI-GEOSTROPHIC SYSTEM WITH VARIABLE CORIOLIS PARAMETER
CLASSICAL SOLUTIONS TO SEMI-GEOSTROPHIC SYSTEM WITH VARIABLE CORIOLIS PARAMETER JINGRUI CHENG, MICHAEL CULLEN, AND MIKHAIL FELDMAN Abstract. We prove short time existence and uniqueness of smooth solutions
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 information3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable
More informationLECTURE 15: COMPLETENESS AND CONVEXITY
LECTURE 15: COMPLETENESS AND CONVEXITY 1. The Hopf-Rinow Theorem Recall that a Riemannian manifold (M, g) is called geodesically complete if the maximal defining interval of any geodesic is R. On the other
More informationOn the velocities of flows consisting of cyclically monotone maps
On the velocities of flows consisting of cyclically monotone maps A. Tudorascu Department of Mathematics West Virginia University Morgantown, WV 2656, USA adriant@math.wvu.edu August 25, 21 Abstract Motivated
More informationIntegral Jensen inequality
Integral Jensen inequality Let us consider a convex set R d, and a convex function f : (, + ]. For any x,..., x n and λ,..., λ n with n λ i =, we have () f( n λ ix i ) n λ if(x i ). For a R d, let δ a
More informationThe optimal partial transport problem
The optimal partial transport problem Alessio Figalli Abstract Given two densities f and g, we consider the problem of transporting a fraction m [0, min{ f L 1, g L 1}] of the mass of f onto g minimizing
More informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More informationSEMI-GEOSTROPHIC SYSTEM WITH VARIABLE CORIOLIS PARAMETER
SEMI-GEOSTROPHIC SYSTEM WITH VARIABLE CORIOLIS PARAMETER JINGRUI CHENG, MICHAEL CULLEN, AND MIKHAIL FELDMAN Abstract. We prove short time existence and uniqueness of smooth solutions ( in C k+2,α with
More information2 (Bonus). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due 9/5). Prove that every countable set A is measurable and µ(a) = 0. 2 (Bonus). Let A consist of points (x, y) such that either x or y is
More information3. (a) What is a simple function? What is an integrable function? How is f dµ defined? Define it first
Math 632/6321: Theory of Functions of a Real Variable Sample Preinary Exam Questions 1. Let (, M, µ) be a measure space. (a) Prove that if µ() < and if 1 p < q
More informationIntegration on Measure Spaces
Chapter 3 Integration on Measure Spaces In this chapter we introduce the general notion of a measure on a space X, define the class of measurable functions, and define the integral, first on a class of
More informationUniqueness of the solution to the Vlasov-Poisson system with bounded density
Uniqueness of the solution to the Vlasov-Poisson system with bounded density Grégoire Loeper December 16, 2005 Abstract In this note, we show uniqueness of weak solutions to the Vlasov- Poisson system
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 informationMetric Spaces and Topology
Chapter 2 Metric Spaces and Topology From an engineering perspective, the most important way to construct a topology on a set is to define the topology in terms of a metric on the set. This approach underlies
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationFrom now on, we will represent a metric space with (X, d). Here are some examples: i=1 (x i y i ) p ) 1 p, p 1.
Chapter 1 Metric spaces 1.1 Metric and convergence We will begin with some basic concepts. Definition 1.1. (Metric space) Metric space is a set X, with a metric satisfying: 1. d(x, y) 0, d(x, y) = 0 x
More informationA LOCALIZATION PROPERTY AT THE BOUNDARY FOR MONGE-AMPERE EQUATION
A LOCALIZATION PROPERTY AT THE BOUNDARY FOR MONGE-AMPERE EQUATION O. SAVIN. Introduction In this paper we study the geometry of the sections for solutions to the Monge- Ampere equation det D 2 u = f, u
More information1 Directional Derivatives and Differentiability
Wednesday, January 18, 2012 1 Directional Derivatives and Differentiability Let E R N, let f : E R and let x 0 E. Given a direction v R N, let L be the line through x 0 in the direction v, that is, L :=
More informationPARTIAL REGULARITY OF BRENIER SOLUTIONS OF THE MONGE-AMPÈRE EQUATION
PARTIAL REGULARITY OF BRENIER SOLUTIONS OF THE MONGE-AMPÈRE EQUATION ALESSIO FIGALLI AND YOUNG-HEON KIM Abstract. Given Ω, Λ R n two bounded open sets, and f and g two probability densities concentrated
More informationPart III. 10 Topological Space Basics. Topological Spaces
Part III 10 Topological Space Basics Topological Spaces Using the metric space results above as motivation we will axiomatize the notion of being an open set to more general settings. Definition 10.1.
More information1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer 11(2) (1989),
Real Analysis 2, Math 651, Spring 2005 April 26, 2005 1 Real Analysis 2, Math 651, Spring 2005 Krzysztof Chris Ciesielski 1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer
More informationDynamic and Stochastic Brenier Transport via Hopf-Lax formulae on Was
Dynamic and Stochastic Brenier Transport via Hopf-Lax formulae on Wasserstein Space With many discussions with Yann Brenier and Wilfrid Gangbo Brenierfest, IHP, January 9-13, 2017 ain points of the
More informationTHEOREMS, ETC., FOR MATH 515
THEOREMS, ETC., FOR MATH 515 Proposition 1 (=comment on page 17). If A is an algebra, then any finite union or finite intersection of sets in A is also in A. Proposition 2 (=Proposition 1.1). For every
More informationREGULARITY OF MONOTONE TRANSPORT MAPS BETWEEN UNBOUNDED DOMAINS
REGULARITY OF MONOTONE TRANSPORT MAPS BETWEEN UNBOUNDED DOMAINS DARIO CORDERO-ERAUSQUIN AND ALESSIO FIGALLI A Luis A. Caffarelli en su 70 años, con amistad y admiración Abstract. The regularity of monotone
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationA generic property of families of Lagrangian systems
Annals of Mathematics, 167 (2008), 1099 1108 A generic property of families of Lagrangian systems By Patrick Bernard and Gonzalo Contreras * Abstract We prove that a generic Lagrangian has finitely many
More informationChapter 2 Metric Spaces
Chapter 2 Metric Spaces The purpose of this chapter is to present a summary of some basic properties of metric and topological spaces that play an important role in the main body of the book. 2.1 Metrics
More informationTools from Lebesgue integration
Tools from Lebesgue integration E.P. van den Ban Fall 2005 Introduction In these notes we describe some of the basic tools from the theory of Lebesgue integration. Definitions and results will be given
More informationLecture 5 Theorems of Fubini-Tonelli and Radon-Nikodym
Lecture 5: Fubini-onelli and Radon-Nikodym 1 of 13 Course: heory of Probability I erm: Fall 2013 Instructor: Gordan Zitkovic Lecture 5 heorems of Fubini-onelli and Radon-Nikodym Products of measure spaces
More information2. Function spaces and approximation
2.1 2. Function spaces and approximation 2.1. The space of test functions. Notation and prerequisites are collected in Appendix A. Let Ω be an open subset of R n. The space C0 (Ω), consisting of the C
More information5 Measure theory II. (or. lim. Prove the proposition. 5. For fixed F A and φ M define the restriction of φ on F by writing.
5 Measure theory II 1. Charges (signed measures). Let (Ω, A) be a σ -algebra. A map φ: A R is called a charge, (or signed measure or σ -additive set function) if φ = φ(a j ) (5.1) A j for any disjoint
More informationSPACES ENDOWED WITH A GRAPH AND APPLICATIONS. Mina Dinarvand. 1. Introduction
MATEMATIČKI VESNIK MATEMATIQKI VESNIK 69, 1 (2017), 23 38 March 2017 research paper originalni nauqni rad FIXED POINT RESULTS FOR (ϕ, ψ)-contractions IN METRIC SPACES ENDOWED WITH A GRAPH AND APPLICATIONS
More informationat time t, in dimension d. The index i varies in a countable set I. We call configuration the family, denoted generically by Φ: U (x i (t) x j (t))
Notations In this chapter we investigate infinite systems of interacting particles subject to Newtonian dynamics Each particle is characterized by its position an velocity x i t, v i t R d R d at time
More informationLocally Lipschitzian Guiding Function Method for ODEs.
Locally Lipschitzian Guiding Function Method for ODEs. Marta Lewicka International School for Advanced Studies, SISSA, via Beirut 2-4, 3414 Trieste, Italy. E-mail: lewicka@sissa.it 1 Introduction Let f
More informationA few words about the MTW tensor
A few words about the Ma-Trudinger-Wang tensor Université Nice - Sophia Antipolis & Institut Universitaire de France Salah Baouendi Memorial Conference (Tunis, March 2014) The Ma-Trudinger-Wang tensor
More informationSome SDEs with distributional drift Part I : General calculus. Flandoli, Franco; Russo, Francesco; Wolf, Jochen
Title Author(s) Some SDEs with distributional drift Part I : General calculus Flandoli, Franco; Russo, Francesco; Wolf, Jochen Citation Osaka Journal of Mathematics. 4() P.493-P.54 Issue Date 3-6 Text
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationIntroduction to Real Analysis Alternative Chapter 1
Christopher Heil Introduction to Real Analysis Alternative Chapter 1 A Primer on Norms and Banach Spaces Last Updated: March 10, 2018 c 2018 by Christopher Heil Chapter 1 A Primer on Norms and Banach Spaces
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationCHAPTER 7. Connectedness
CHAPTER 7 Connectedness 7.1. Connected topological spaces Definition 7.1. A topological space (X, T X ) is said to be connected if there is no continuous surjection f : X {0, 1} where the two point set
More informationTopological properties
CHAPTER 4 Topological properties 1. Connectedness Definitions and examples Basic properties Connected components Connected versus path connected, again 2. Compactness Definition and first examples Topological
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More information+ 2x sin x. f(b i ) f(a i ) < ɛ. i=1. i=1
Appendix To understand weak derivatives and distributional derivatives in the simplest context of functions of a single variable, we describe without proof some results from real analysis (see [7] and
More informationEuler Poisson systems as action-minimizing paths in the Wasserstein space
ARMA manuscript No. (will be inserted by the editor) Euler Poisson systems as action-minimizing paths in the Wasserstein space W. Gangbo, T. Nguyen, A. Tudorascu Abstract This paper uses a variational
More informationSecond order forward-backward dynamical systems for monotone inclusion problems
Second order forward-backward dynamical systems for monotone inclusion problems Radu Ioan Boţ Ernö Robert Csetnek March 6, 25 Abstract. We begin by considering second order dynamical systems of the from
More informationTopological vectorspaces
(July 25, 2011) Topological vectorspaces Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ Natural non-fréchet spaces Topological vector spaces Quotients and linear maps More topological
More informationPart V. 17 Introduction: What are measures and why measurable sets. Lebesgue Integration Theory
Part V 7 Introduction: What are measures and why measurable sets Lebesgue Integration Theory Definition 7. (Preliminary). A measure on a set is a function :2 [ ] such that. () = 2. If { } = is a finite
More informationHESSIAN VALUATIONS ANDREA COLESANTI, MONIKA LUDWIG & FABIAN MUSSNIG
HESSIAN VALUATIONS ANDREA COLESANTI, MONIKA LUDWIG & FABIAN MUSSNIG ABSTRACT. A new class of continuous valuations on the space of convex functions on R n is introduced. On smooth convex functions, they
More informationC 1 regularity of solutions of the Monge-Ampère equation for optimal transport in dimension two
C 1 regularity of solutions of the Monge-Ampère equation for optimal transport in dimension two Alessio Figalli, Grégoire Loeper Abstract We prove C 1 regularity of c-convex weak Alexandrov solutions of
More informationNotions such as convergent sequence and Cauchy sequence make sense for any metric space. Convergent Sequences are Cauchy
Banach Spaces These notes provide an introduction to Banach spaces, which are complete normed vector spaces. For the purposes of these notes, all vector spaces are assumed to be over the real numbers.
More informationThe Arzelà-Ascoli Theorem
John Nachbar Washington University March 27, 2016 The Arzelà-Ascoli Theorem The Arzelà-Ascoli Theorem gives sufficient conditions for compactness in certain function spaces. Among other things, it helps
More informationA description of transport cost for signed measures
A description of transport cost for signed measures Edoardo Mainini Abstract In this paper we develop the analysis of [AMS] about the extension of the optimal transport framework to the space of real measures.
More informationMañé s Conjecture from the control viewpoint
Mañé s Conjecture from the control viewpoint Université de Nice - Sophia Antipolis Setting Let M be a smooth compact manifold of dimension n 2 be fixed. Let H : T M R be a Hamiltonian of class C k, with
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationWiener Measure and Brownian Motion
Chapter 16 Wiener Measure and Brownian Motion Diffusion of particles is a product of their apparently random motion. The density u(t, x) of diffusing particles satisfies the diffusion equation (16.1) u
More informationLecture Notes in Advanced Calculus 1 (80315) Raz Kupferman Institute of Mathematics The Hebrew University
Lecture Notes in Advanced Calculus 1 (80315) Raz Kupferman Institute of Mathematics The Hebrew University February 7, 2007 2 Contents 1 Metric Spaces 1 1.1 Basic definitions...........................
More informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
More informationSobolev spaces. May 18
Sobolev spaces May 18 2015 1 Weak derivatives The purpose of these notes is to give a very basic introduction to Sobolev spaces. More extensive treatments can e.g. be found in the classical references
More informationFULL CHARACTERIZATION OF OPTIMAL TRANSPORT PLANS FOR CONCAVE COSTS
FULL CHARACTERIZATION OF OPTIMAL TRANSPORT PLANS FOR CONCAVE COSTS PAUL PEGON, DAVIDE PIAZZOLI, FILIPPO SANTAMBROGIO Abstract. This paper slightly improves a classical result by Gangbo and McCann (1996)
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 informationMATHS 730 FC Lecture Notes March 5, Introduction
1 INTRODUCTION MATHS 730 FC Lecture Notes March 5, 2014 1 Introduction Definition. If A, B are sets and there exists a bijection A B, they have the same cardinality, which we write as A, #A. If there exists
More informationProduct measures, Tonelli s and Fubini s theorems For use in MAT4410, autumn 2017 Nadia S. Larsen. 17 November 2017.
Product measures, Tonelli s and Fubini s theorems For use in MAT4410, autumn 017 Nadia S. Larsen 17 November 017. 1. Construction of the product measure The purpose of these notes is to prove the main
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 informationHILBERT SPACES AND THE RADON-NIKODYM THEOREM. where the bar in the first equation denotes complex conjugation. In either case, for any x V define
HILBERT SPACES AND THE RADON-NIKODYM THEOREM STEVEN P. LALLEY 1. DEFINITIONS Definition 1. A real inner product space is a real vector space V together with a symmetric, bilinear, positive-definite mapping,
More informationthe convolution of f and g) given by
09:53 /5/2000 TOPIC Characteristic functions, cont d This lecture develops an inversion formula for recovering the density of a smooth random variable X from its characteristic function, and uses that
More informationOptimal Transportation. Nonlinear Partial Differential Equations
Optimal Transportation and Nonlinear Partial Differential Equations Neil S. Trudinger Centre of Mathematics and its Applications Australian National University 26th Brazilian Mathematical Colloquium 2007
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 informationX n D X lim n F n (x) = F (x) for all x C F. lim n F n(u) = F (u) for all u C F. (2)
14:17 11/16/2 TOPIC. Convergence in distribution and related notions. This section studies the notion of the so-called convergence in distribution of real random variables. This is the kind of convergence
More informationGLOBAL LIPSCHITZ CONTINUITY FOR MINIMA OF DEGENERATE PROBLEMS
GLOBAL LIPSCHITZ CONTINUITY FOR MINIMA OF DEGENERATE PROBLEMS PIERRE BOUSQUET AND LORENZO BRASCO Abstract. We consider the problem of minimizing the Lagrangian [F ( u+f u among functions on R N with given
More informationREAL AND COMPLEX ANALYSIS
REAL AND COMPLE ANALYSIS Third Edition Walter Rudin Professor of Mathematics University of Wisconsin, Madison Version 1.1 No rights reserved. Any part of this work can be reproduced or transmitted in any
More informationConvexity in R n. The following lemma will be needed in a while. Lemma 1 Let x E, u R n. If τ I(x, u), τ 0, define. f(x + τu) f(x). τ.
Convexity in R n Let E be a convex subset of R n. A function f : E (, ] is convex iff f(tx + (1 t)y) (1 t)f(x) + tf(y) x, y E, t [0, 1]. A similar definition holds in any vector space. A topology is needed
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 informationStanford Mathematics Department Math 205A Lecture Supplement #4 Borel Regular & Radon Measures
2 1 Borel Regular Measures We now state and prove an important regularity property of Borel regular outer measures: Stanford Mathematics Department Math 205A Lecture Supplement #4 Borel Regular & Radon
More informationA new Hellinger-Kantorovich distance between positive measures and optimal Entropy-Transport problems
A new Hellinger-Kantorovich distance between positive measures and optimal Entropy-Transport problems Giuseppe Savaré http://www.imati.cnr.it/ savare Dipartimento di Matematica, Università di Pavia Nonlocal
More informationJoint work with Nguyen Hoang (Univ. Concepción, Chile) Padova, Italy, May 2018
EXTENDED EULER-LAGRANGE AND HAMILTONIAN CONDITIONS IN OPTIMAL CONTROL OF SWEEPING PROCESSES WITH CONTROLLED MOVING SETS BORIS MORDUKHOVICH Wayne State University Talk given at the conference Optimization,
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 information16 1 Basic Facts from Functional Analysis and Banach Lattices
16 1 Basic Facts from Functional Analysis and Banach Lattices 1.2.3 Banach Steinhaus Theorem Another fundamental theorem of functional analysis is the Banach Steinhaus theorem, or the Uniform Boundedness
More informationOn duality theory of conic linear problems
On duality theory of conic linear problems Alexander Shapiro School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia 3332-25, USA e-mail: ashapiro@isye.gatech.edu
More informationII - REAL ANALYSIS. This property gives us a way to extend the notion of content to finite unions of rectangles: we define
1 Measures 1.1 Jordan content in R N II - REAL ANALYSIS Let I be an interval in R. Then its 1-content is defined as c 1 (I) := b a if I is bounded with endpoints a, b. If I is unbounded, we define c 1
More informationCourse 212: Academic Year Section 1: Metric Spaces
Course 212: Academic Year 1991-2 Section 1: Metric Spaces D. R. Wilkins Contents 1 Metric Spaces 3 1.1 Distance Functions and Metric Spaces............. 3 1.2 Convergence and Continuity in Metric Spaces.........
More informationBehaviour of Lipschitz functions on negligible sets. Non-differentiability in R. Outline
Behaviour of Lipschitz functions on negligible sets G. Alberti 1 M. Csörnyei 2 D. Preiss 3 1 Università di Pisa 2 University College London 3 University of Warwick Lars Ahlfors Centennial Celebration Helsinki,
More informationChapter 1. Measure Spaces. 1.1 Algebras and σ algebras of sets Notation and preliminaries
Chapter 1 Measure Spaces 1.1 Algebras and σ algebras of sets 1.1.1 Notation and preliminaries We shall denote by X a nonempty set, by P(X) the set of all parts (i.e., subsets) of X, and by the empty set.
More information7 Complete metric spaces and function spaces
7 Complete metric spaces and function spaces 7.1 Completeness Let (X, d) be a metric space. Definition 7.1. A sequence (x n ) n N in X is a Cauchy sequence if for any ɛ > 0, there is N N such that n, m
More informationOptimization and Optimal Control in Banach Spaces
Optimization and Optimal Control in Banach Spaces Bernhard Schmitzer October 19, 2017 1 Convex non-smooth optimization with proximal operators Remark 1.1 (Motivation). Convex optimization: easier to solve,
More informationCombinatorics in Banach space theory Lecture 12
Combinatorics in Banach space theory Lecture The next lemma considerably strengthens the assertion of Lemma.6(b). Lemma.9. For every Banach space X and any n N, either all the numbers n b n (X), c n (X)
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 information4th Preparation Sheet - Solutions
Prof. Dr. Rainer Dahlhaus Probability Theory Summer term 017 4th Preparation Sheet - Solutions Remark: Throughout the exercise sheet we use the two equivalent definitions of separability of a metric space
More informationReal Analysis Problems
Real Analysis Problems Cristian E. Gutiérrez September 14, 29 1 1 CONTINUITY 1 Continuity Problem 1.1 Let r n be the sequence of rational numbers and Prove that f(x) = 1. f is continuous on the irrationals.
More informationLectures on. Sobolev Spaces. S. Kesavan The Institute of Mathematical Sciences, Chennai.
Lectures on Sobolev Spaces S. Kesavan The Institute of Mathematical Sciences, Chennai. e-mail: kesh@imsc.res.in 2 1 Distributions In this section we will, very briefly, recall concepts from the theory
More information2. Dual space is essential for the concept of gradient which, in turn, leads to the variational analysis of Lagrange multipliers.
Chapter 3 Duality in Banach Space Modern optimization theory largely centers around the interplay of a normed vector space and its corresponding dual. The notion of duality is important for the following
More informationMeasure and Integration: Solutions of CW2
Measure and Integration: s of CW2 Fall 206 [G. Holzegel] December 9, 206 Problem of Sheet 5 a) Left (f n ) and (g n ) be sequences of integrable functions with f n (x) f (x) and g n (x) g (x) for almost
More informationRecall that if X is a compact metric space, C(X), the space of continuous (real-valued) functions on X, is a Banach space with the norm
Chapter 13 Radon Measures Recall that if X is a compact metric space, C(X), the space of continuous (real-valued) functions on X, is a Banach space with the norm (13.1) f = sup x X f(x). We want to identify
More informationOptimal Transport for Applied Mathematicians
Optimal Transport for Applied Mathematicians Calculus of Variations, PDEs and Modelling Filippo Santambrogio 1 1 Laboratoire de Mathématiques d Orsay, Université Paris Sud, 91405 Orsay cedex, France filippo.santambrogio@math.u-psud.fr,
More informationTopological properties of Z p and Q p and Euclidean models
Topological properties of Z p and Q p and Euclidean models Samuel Trautwein, Esther Röder, Giorgio Barozzi November 3, 20 Topology of Q p vs Topology of R Both R and Q p are normed fields and complete
More informationNewtonian Mechanics. Chapter Classical space-time
Chapter 1 Newtonian Mechanics In these notes classical mechanics will be viewed as a mathematical model for the description of physical systems consisting of a certain (generally finite) number of particles
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 information