GENERAL EXISTENCE OF SOLUTIONS TO DYNAMIC PROGRAMMING PRINCIPLE. 1. Introduction
|
|
- Virgil Moody
- 5 years ago
- Views:
Transcription
1 GENERAL EXISTENCE OF SOLUTIONS TO DYNAMIC PROGRAMMING PRINCIPLE QING LIU AND ARMIN SCHIKORRA Abstract. We provide an alternative approach to the existence of solutions to dynamic programming equations arising in the discrete gametheoretic interpretations for various nonlinear partial differential equations including the infinity Laplacian, mean curvature flow and Hamilton-Jacobi type. Our general result is similar to Perron s method but adapted to the discrete situation. 1. Introduction Recently, the deterministic or stochastic game-theoretic approach to various nonlinear partial differential equations has attracted a lot of attention; see, for example, a variety of games in [10, 11, 20, 21, 7, 19, 1, 5] and related topics in [13, 18, 17, 22]. The results in the literature so far can be summarized in the following way. Consider a partial differential equation in a bounded domain Ω R n : { F (x, u, 2 u) = 0 in Ω, (1.1) u = g on Ω, where F is a function defined on Ω R n S n and g is a continuous function on Ω. Here S n denotes the set of all n n symmetric matrices. Suppose there exists a unique (viscosity) solution u of the equation. Then (under certain additional conditions) one may construct a deterministic or stochastic discrete game, with step size ε > 0, whose value function u ε converges to u (locally) uniformly as ε 0. The most important step in the proof of these game-theoretic approximations is to establish the so-called dynamic programming principle (2.1) for each particular F. More precisely, one proves that u ε satisfies a relation between its value at the current position and that of the next possible positions, which 2010 Mathematics Subject Classification. 35A35, 49C20, 91A05, 91A15. Key words and phrases. Dynamic programming principle, differential games, nonlinear partial differential equations. A.S. was supported by DAAD fellowship D/12/
2 2 Q. LIU AND A. SCHIKORRA is reminiscent of nonlinear semigroup-properties. A formal derivation for (1.1) from the games often follows by adopting the Taylor expansion on its dynamic programming principle. One example we have in mind is the Dirichlet boundary problem for the infinity Laplacian with continuous boundary data g : Ω R: ( ) u u u := tr 2 u = 0 in Ω, u 2 u = g on Ω, whose unique viscosity solution u can be approximated locally uniformly by solutions u ε to the DPP associated with a tug-of-war game, first proposed in [20], u ε (x) = 1 2 sup u ε + 1 B ε(x) 2 inf u ε for x Ω B ε(x) (1.2) u ε = g on Ω. Let us define sub- and supersolutions for this particular problem. Definition 1.1. A function u : Ω R with sup Ω u < is called a subsolution of (1.2), if it satisfies u(x) 1 2 sup u + 1 B ε(x) 2 inf u for x Ω, B ε(x) u g on Ω. The set of all subsolutions is denoted by S. Conversely, a function u : Ω R with inf Ω u > is called a supersolution of (1.2), if u(x) 1 2 sup u + 1 B ε(x) 2 inf u for x Ω, B ε(x) u g on Ω. The set of all supersolutions is denoted by S. A function u : Ω R is a solution, if it is simultaneously a subsolution and a supersolution. In the literature this is also called an infinity harmonious function, cf. [13]. The set of all solutions is denoted by S. A first easy observation is that S and S are non-empty. Example 1.2. The function u : Ω R given by { inf Ω g x Ω u(x) := g x Ω.
3 EXISTENCE FOR DYNAMIC PROGRAMMING PRINCIPLE 3 is a subsolution, the function v : Ω R given by { sup v(x) := Ω g x Ω g x Ω. is a supersolution. The next observation is that any subsolution is uniformly bounded from above, and any supersolution is uniformly bounded from below. This follows from the following maximum principle: Proposition 1.3 (Maximum principle). Any subsolution u S satisfies sup Ω and any supersolution u S satisfies inf Ω u sup g, Ω u inf Ω g. Proof. We are only going to show the subsolution case. It suffices to prove that for any ν > 0, sup u sup g + ν. (1.3) Ω Ω So fix ν > 0, and set K := diam Ω + 1. For δ := (K + 1) 1 ε, let x ε 0 Ω sup u u(x 0 ) + δ. Ω If x 0 Ω, then (1.3) is proven. If x 0 Ω, we pick a sequence (x i ) K i=0 such that x i B ε (x i 1 ) and x K Ω. Since u is a subsolution, u(x 0 ) 1 2 sup u + 1 B ε(x 0 ) 2 inf u 1 B ε(x 0 ) 2 u(x 0) u(x 1) δ. Consequently, sup u u(x 0 ) + δ u(x 1 ) + 2δ. Ω By induction, sup u u(x i ) + (i + 1)δ, Ω and in particular, since x K Ω, by the choice of δ, sup Ω u u(x K ) + (K + 1)δ sup g + ν. Ω
4 4 Q. LIU AND A. SCHIKORRA Having this boundedness, we can define the largest subsolution and the smallest supersolution u(x) := sup u(x), x Ω, u S u(x) := inf u(x), x Ω. u S Then we have the following result, which extends the classical Perron s method used in viscosity solution theory [8, 6]. Theorem 1.4 (Perron s Method). Both u and u are solutions to (1.2). Proof. We only give a proof for u. It is easy to check that u S. Indeed, for any subsolution u, we have u g on Ω and u(x) 1 2 sup u + 1 B ε(x) 2 inf u for any x Ω. B ε(x) A pointwise supremum yields u g on Ω and u(x) 1 2 sup u + 1 B ε(x) 2 inf u for any x Ω. B ε(x) It suffices to prove that u S. Since u is a subsolution, so is v : Ω R given by 1 sup u + 1 inf u for x Ω 2 2 v(x) := B ε(x) B ε(x) g on Ω. Since u is the pointwise supremum of all subsolutions, in particular 1 2 sup u + 1 B ε(x) 2 inf u = v(x) u(x), B ε(x) for x Ω and g u on Ω, we deduce that u is a supersolution. Consequently, u is a solution. The arguments leading to Theorem 1.4 are very general. Indeed, they allow us to give an elementary approach for finding solutions to DPPs for a very general class of equations in a general metric space X (equipped with a measure when necessary); Namely, for U {u : X R} and a given operator T : U U we find solutions u to u(x) = T u(x) for x X. (1.4)
5 EXISTENCE FOR DYNAMIC PROGRAMMING PRINCIPLE 5 In our case above, T is given by T u(x) = 1 sup u B ε(x) T u(x) = g(x) inf u for x Ω, B ε(x) for x Ω. (1.5) The main properties for T and U are as follows: (H1) The operator T is monotone. (H2) There exists at least one subsolution (or supersolution). (H3) There is a uniform upper bound for all subsolutions (or a uniform lower bound for all supersolutions). (H4) The function space U is closed under the pointwise supremum operations (or infimum operations). We specify the precise meaning of these conditions in Section 2, and discuss examples of known DPPs and there respective operator T in Section 5. Our main result, Theorem 2.2, is that if T satisfies (H1) (H4), then there exists a solution to (1.4). Usually, the game-theoretic approach and Perron s method are seen as separate ways to obtain existence of solutions to fully nonlinear PDEs. Our result indicates that they are connected on the discrete level. It is not unusual to assume the monotonicity (H1). Similar conditions are used in schemes for fully nonlinear second order PDE, see, e.g. [4], which are elliptic, i.e., F (x, p, A) F (x, p, B) for all x Ω = Ω Ω, p R n and n n symmetric matrices A, B with A B. The upper bound for subsolutions (H3) is necessary as to justify the pointwise supremum of subsolutions. It is, however, non-trivial to obtain the boundedness. In this work, we construct a bounded strict supersolutions and prove a comparison principle for any subsolution and strict supersolution, Theorem 3.2. We show the construction of strict supersolutions for each example in Section 5. Concerning (H4), we must choose a suitable function space U to guarantee the closedness. The classical choice for elliptic and parabolic PDEs is the (E)
6 6 Q. LIU AND A. SCHIKORRA space of all upper/lower semicontinuous solutions and an application of semicontinuous envelope is also necessary; see for example [8, 6]. In our case, we generally cannot expect even semicontinuity, see the examples in [16]. In the case of integral DPPs, U has to be a subset of all measurable functions, which makes it more difficult to ensure (H4). The remaining issue is the uniqueness of solutions to DPP, which is equivalent to the usual comparison principle (as shown in Lemma 2.5), since we have Perron s method in hand for existence of solutions. Once the uniqueness is established, we may conclude that the solutions we found via Perron s method coincide with the value function obtained in games. In [16] we obtained uniqueness for a particular problem related to the biased tug-of-war game with nonzero running payoff. In Section 4, we generalize this to dynamic programming equations in the presence of running payoffs. However, the uniqueness problem in general, especially without running payoffs, is still not known: The classical approach of a comparison principle of sub- and supersolutions, in the case of (1.2) see [2], relies on semicontinuity which in general we cannot expect. Another earlier related result for existence and uniqueness is due to [13] but for a quite different DPP; the radius ε of the ball B ε (x) where the extrema are taken depends on x and diminishes as x approaches the boundary of the domain. In our general setting of (1.2), one might hope that a finer analysis of game-trees as in [16] leads to answers related to this question. Acknowledgments. The authors thank Juan J. Manfredi and Marta Lewicka for their interests and valuable suggestions. 2. General Existence We consider a function space U {u : X R}, and an operator T : U U. Then, we look for solutions u U to T u = u in X. (2.1) Let us start defining sub- and supersolutions to (2.1). The following definition is consistent with the usual classical definition of sub- and supersolutions of PDEs with Dirichlet boundary conditions. Definition 2.1. Given T : U U, a function u U is called a subsolution of (2.1) if sup X u < and u T u in X.
7 EXISTENCE FOR DYNAMIC PROGRAMMING PRINCIPLE 7 Similarly, a function u U is called a supersolution of (2.1) if inf X u > and u T u in X. A function u U is called a solution if it is both a subsolution and a supersolution, i.e., if T u = u pointwise in X. In the case of tug-of-war (1.2), this definition is consistent with Definition 1.1, with T defined as in (1.5). Similar to classical Perron s method for the existence of solutions, our strategy is to take either the pointwise supremum of all subsolutions or the pointwise infimum of all supersolutions. To accomplish the former, we assume the following conditions on T and U. (A1) (Monotonicity) If u v in X then T u T v in X. (A2) (Non-emptyness for subsolutions) There exists at least one subsolution of (2.1). (A3) (Boundedness of subsolutions) For a uniform constant C = C(T ) > 0 any subsolution u of (2.1) satisfies that sup X u C. (A4) (Closedness under supremum) If Ũ U, and then sup X sup u <, u Ũ ũ 1 (x) := sup u(x) U. u Ũ If one chooses to consider the infimum of all supersolutions, then we should replace (A2) (A4) above with (B2) (B4) below. (B2) (Non-emptiness for supersolutions) There exists at least one supersolution of (2.1). (B3) (Boundedness of supersolutions) For a uniform constant C = C(T ) > 0 any supersolution u of (2.1) satisfies that inf X u C. (B4) (Closedness under infimum) If Ũ U, and then inf X inf u >, u Ũ ũ 2 (x) := inf u(x) U. u Ũ
8 8 Q. LIU AND A. SCHIKORRA Theorem 2.2 (Existence, I). Let X be a set of points and U {u : X R}, and let T : U U. Assume T and U satisfy (A1), (A2), (A3) and (A4). Let U be the set of subsolutions and Then u U is a solution to (2.1). u(x) := sup u(x). u U Theorem 2.3 (Existence, II). Let X be a set of points and U {u : X R}, and let T : U U. Assume T and U satisfy (A1), (B2), (B3) and (B4). Let U be the set of supersolutions and Then u U is a solution to (2.1). u(x) := inf u(x). u U We only prove Theorem 2.2, since the argument for Theorem 2.3 is analogous. Proof of Theorem 2.2. For any subsolution u U, by the monotonicity (A1), we get T u U. By (A3), all elements u of U are bounded by a uniform constant, thus u < is well defined. By the monotonicity (A1), we have that u T u T u in X, for any u U. Taking pointwise the supremum, this implies u T u in X, and thus u U. As observed above, this implies in particular that T u U, and since u is the pointwise supremum, we have Thus, u is a solution. T u u in X. Remark 2.4. It is straight-forward to extend the proof of Theorem 2.2 to a more general function space U {u : X Y }, where Y is a conditionally complete lattice Y. Theorem 2.2 can be viewed as a general version of Perron s method. The regularity assumptions on u are minimal. It is also clear that if all subsolutions are bounded from above and all supersolutions are bounded from below, then u is the maximal solution of (2.1) and u is the minimal solution of (2.1).
9 EXISTENCE FOR DYNAMIC PROGRAMMING PRINCIPLE 9 It thus suffices to show that u u to conclude the uniqueness of the solutions but in this generalality this remains an open problem. It is unlikely that the usual comparison principle holds in such a general setting, but it would be interesting to find conditions on T when such a principle holds. The comparison principle is however equivalent to uniqueness: Proposition 2.5. Assume that T : U U satisfies (A1) (A4) and (B2) (B4). Then the following are equivalent (i) For any u 1, u 2 Ũ with T u i = u i in X, i = 1, 2, then u 1 = u 2 in X. (ii) For any u 1, u 2 : X R such that u 1 T u 1 and u 2 T u 2 we have u 1 u 2 in X. Proof. (i) (ii): Let u, u : X R be the solutions obtained respectively from Theorem 2.2 and Theorem 2.3. Note that any subsolutions u satisfies u u and any supersolution v satisfies v u. The assumption (i) implies that u 1 u = u u 2. (ii) (i): Since both solutions u 1, u 2 are both, sub- and supersolutions, we have immediately u 1 u 2 u 1, and thus u 1 = u 2. Another important issue is to ensure the uniform boundedness (A3) of subsolutions (or supersolution, respectively). In the Section 3 we discuss this in more detail. 3. The Boundedness We discuss for subsolutions u of (2.1) the boundedness from above. In [16] we obtained boundedness for a quite general type of DPP using an iteration method. Here, we give a different type of proof for our situation: The idea is to construct a strict supersolution of (2.1) and to use a weaker type of comparison principle. Let us first give a definition for strict sub- and supersolutions. Definition 3.1. For any given operator T : U U, a function u U with sup X u < is called a strict subsolution of (2.1) if there exists a constant σ > 0 such that u T u σ in X.
10 10 Q. LIU AND A. SCHIKORRA Similarly, a function u U with inf X u > is called a strict supersolution of (2.1) if there exists a σ > 0 such that u T u + σ in X. It is clear that a strict subsolution (resp., strict supersolution) is a subsolution (resp., supersolution). Again, the above definition in practice contains the boundary condition u = g on X \ Ω. Theorem 3.2 (Strict comparison theorem). Let U be a linear space. Let T : U U be any given operator that satisfies monotonicity (A1) of Theorem 2.2. Assume that for any x X, any c R and any u U, the following two relations holds: T (u + c)(x) (T u)(x) + c, (3.1) Then any subsolution u U of (2.1) and any strict supersolution v U of (2.1) satisfies u v in X. The same relation holds if u is a strict subsolution and v is a supersolution. Proof of Theorem 3.2. By definition, we have and there exists σ > 0 such that Suppose by contradiction that for some m (0, ) u T u in X (3.2) v T v + σ in X. (3.3) sup(u v) = m. X Then we take ˆv := v + m. By (3.1),ˆv is also a strict supersolution. Indeed, v(x) + m (T v)(x) + m + σ T (v + m)(x) + σ. Also, by (3.3), T ˆv(x) T v(x) + m for any x X. (3.4) Moreover, sup(u ˆv) = 0. X Then for any δ > 0 there exists x δ Ω such that u(x δ ) ˆv(x δ ) δ. (3.5) On the other hand, u(x) ˆv(x) for all x X. The monotonicity condition (A1) implies that T u(x δ ) T ˆv(x δ ). (3.6)
11 EXISTENCE FOR DYNAMIC PROGRAMMING PRINCIPLE 11 Combining the inequalities (3.2), (3.4), (3.5) and (3.6), we get v(x δ ) + m = ˆv(x δ ) u(x δ ) + δ T u(x δ ) + δ T ˆv(x δ ) + δ T v(x δ ) + m + δ, and consequently v(x δ ) T v(x δ ) + δ. which contradicts (3.3) if we choose δ < σ. In order to obtain the boundedness of all subsolutions in terms of the given operator T, it is therefore important to build a strict supersolution that is bounded from above. Corollary 3.3 (Boundedness). Suppose that T : U U is an operator that satisfies the conditions of Theorem 2.2. Assume that there exists C R depending on T and a strict supersolution v U of (2.1) satisfying v C. Then u C for all subsolutions u of (2.1). Similarly, if there exists C R depending on T and a strict subsolution v U of (2.1) satisfying v C, then u C for all supersolutions u of (2.1). This result follows immediately from Theorem 3.2. However, there seems to be no universal method to get the existence of a bounded strict supersolution for a general T. One needs to discuss it case by case. 4. Comparison Principle: A Special Case For general discrete dynamic programming equations, the uniqueness problem is challenging. In [16], we showed that in the case of the infinity laplacian as in (1.2), uniqueness follows if one assumes running costs. Technically, there it was shown that a certain discretized flow converges uniformly to the solution, which implies as an immediate corollary the uniqueness. Here we show, that running costs imply uniqueness in our more general context. Uniqueness follows from the comparison result, Theorem 3.2, by approximating a supersolution (or subsolution) with a strict supersolution (or strict subsolution). This method works well especially when the corresponding games bear positive or negative running payoffs. The idea in what follows is inspired by [9]. Theorem 4.1. Let U be a linear space. Let T : U U be any given operator that satisfies monotonicity (A1) of Theorem 2.2. Assume that for any x X, any c R and any u U, (3.1) holds. In addition let T satisfy
12 12 Q. LIU AND A. SCHIKORRA (A5) for any u U and any λ > 1, there exists σ > 0 such that T (λu) λt u σ. Then any subsolution u U of (2.1) and any supersolution v U of (2.1) satisfies u v in X. Proof. Suppose u and v are respectively a subsolution and a supersolution of (2.1). It is easily seen that under the assumption (A5), λv is a strict supersolution for any λ > 1. Indeed, there exists σ > 0 depending on λ and v such that λv λt v T (λv) + σ. Thus, using Theorem 3.2, we have u λv for any λ > 1. We conclude the proof by sending λ 1. In order to use Theorem 4.1, an extra step may be needed in practical use: (A5) is not necessarily satisfied if T assigns nonpositive values on the boundary. In this case, one only needs to shift the boundary value up to make it positive in the definition of T. The uniqueness of solutions then follows immediately. See the examples below for more details. 5. Examples of Dynamic Programming Principle We give several typical examples in R n. Hereafter B ε (x) denotes the ball centered at x with radius ε > 0 while B ε simply means B ε (0). We first consider (degenerate) elliptic equations in a bounded domain Ω R n with Dirichlet boundary condition g C( Ω). In our exposition below, it is enough to assume that g is bounded. To connect the results in the previous sections, we let X = Ω = Ω Ω. Example 5.1. A more general variant of the example discussed in Introduction is related to the so called biased tug-of-war games. More precisely, one considers the problem { u + c u = f(x) in Ω, (5.1) u = g on Ω, (5.2) where c 0 is a fixed constant and f is assumed to be a bounded function on Ω. The associated DPP is discussed in [19] for the case f 0 using games
13 EXISTENCE FOR DYNAMIC PROGRAMMING PRINCIPLE 13 and in [16] for the case min Ω f > 0 with a tree approach. The DPP for the value function is given by u ε (x) = µ sup u ε (y) + (1 µ) inf y B ε(x) y B uε (y) + f(x) ε(x) 2 ε2 in Ω, (5.3) u ε (x) = g(x) on Ω,(5.4) where µ = 1/2 cε/4 with ε > 0 small such that µ (0, 1). We take U to be the set of all functions Ω R, which is clearly closed under the operators of supremum and infimum. Let f(x) µ sup u(y) + (1 µ) inf u(y) + T u(x) := y B ε(x) y B ε(x) 2 ε2 for x Ω g(x) for x Ω, with µ = 1/2 cε/4 with ε > 0 small such that µ (0, 1), It is not difficult to see that T satisfies (A1) in Theorem 2.2 and (3.1) as well. The existence of a subsolution is also easily obtained by taking u inf Ω g. In order to show the boundedness property (A3), we construct the following strict supersolution v. Suppose the diameter of domain Ω is D. Let σ > 0. We divide the domain Ω into subregions according to the distance to the boundary Ω. Set Ω 0 := R n \ Ω, (5.5) Ω k := {x Ω : (k 1)ε < dist (x, Ω) kε} for all k = 1, 2,..., N, where N < is the total number for the partition, as Ω is bounded. It is clear that Ω = N k=1 Ω k. We then define a bounded function v U as follows: v(x) := a(e kε 1) + sup g + 1 Ω when x Ω k for any k = 0, 1, 2,..., N where a > 0 is determined later. Note that v(x) g(x) + 1 = T (x) + 1 for any x Ω and v(x) T v(x) = v(x) µ sup v (1 µ) inf v f(x) B ε(x) B ε(x) 2 ε2 = ae ( (k 1)ε e ε + µe 2ε + (1 µ) ) f(x) ( ) c + 1 ae D ε 2 ε 2 sup f + o(ε 2 ) 2 Ω 2 ε2 for all x Ω when ε > 0 is small. Since c 0, it follows that v is a strict supersolution if a is large enough and therefore all of the subsolutions are
14 14 Q. LIU AND A. SCHIKORRA bounded from above by Proposition 3.3. We thus conclude that there exists a solution to (5.3). The solutions are unique if σ 1 := inf Ω f > 0, as follows from the arguments in [16]. For an alternative argument, we first notice that a constant C in X is a strict subsolution for a sufficiently large C > 0 and therefore u C in X for any solution u. It suffices to show that û = u + C uniquely solves û = ˆT û with µ sup u(y) + (1 µ) ˆT u(x) := y B ε(x) g(x) + C inf y B ε(x) u(y) + f(x) 2 ε2 for x Ω for x Ω, where C is a positive constant satisfying G := inf x Ω g(x) + C > 0. Note that ˆT satisfies (A5) in Theorem 4.1: ˆT (λu)(x) λ ˆT u(x) + 1 λ f(x) λ 2 u(x) (1 λ)σ 1 for x Ω; ˆT (λu)(x) = ˆT u(x) λ ˆT u(x) (1 λ)g for x Ω. We may use Theorem 4.1 to reach the conclusion. Example 5.2 (Stationary mean curvature operator). A typical elliptic PDE involving level set mean curvature operator is as follows: ( ) u u div 1 = 0 in Ω, (5.6) (MCF) u u = 0 on Ω, (5.7) where Ω is a bounded domain in R n. A deterministic game-theoretic interpretation is also available, given in [10]. For simplicity, we only study the case when n = 2. The DPP is then written as u ε (x) = inf sup u ε (x + 2bw) + ε 2 for x Ω, w B ε b=±1 (5.8) u ε (x) = 0 for x Ω. See more details about the related games in [10, 14, 15]. We here apply Theorem 2.2 to seek a solution to (5.8). Let U be the set of all functions: Ω R. The DPP in this case is again written as u = T u for any u U with inf sup u(x + 2bw) + ε 2 for x Ω, T u(x) := w B ε b=±1 0 for x Ω. It is clear that u 0 in Ω is a subsolution of u = T u. It is then easily seen that T and U satisfy the conditions (A1), (A2) and (A4) in Theorem 2.2. It
15 EXISTENCE FOR DYNAMIC PROGRAMMING PRINCIPLE 15 remains to verify (A2). We may take R > 0 such that Ω B R. We then take a partition for the disk B R : O 1 = B 2ε, O k = B 2kε \ O k 1, for k = 2, 3,... Without loss, we assume R 2 = 2K for some K N, which implies that K k=1 O k = B R. We denote O K+1 = R 2 \ B R Now we define v U to be v(x) := 2(K + 1 k)ε 2 if x O k for 1 k K + 1 and claim that v is a bounded strict supersolution in Ω. The boundedness is clear. To show that v is a strict supersolution, we first notice that v(x) 2ε 2 when x Ω. Moreover, for any x Ω, there exists 1 k K such that x O k and therefore v(x) = 2(K + 1 k)ε 2. It is clear that x ± 2w O k+1 for w B ε (x) and orthogonal to x (w can be arbitrary in B ε (x) if x = 0). This yields and therefore v(x) T v(x) + ε 2. T v(x) = 2(K k)ε 2 + ε 2 = v(x) ε 2 The uniqueness of solutions holds as well in this case. The proof, omitted here, is based on Theorem 4.1 and similar to that in Example 5.1. We remark that despite our solution is the same as the value function u ε in games, it still remains as an open question whether the solution converges to a unique solution of (MCF) as ε 0, when Ω is a general non-convex domain [10]. Example 5.3 (Discrete games for the Eikonal equations). Our general DPP even applies to the first order Hamilton-Jacobi equations, but it is a discretized version of those well studied, for instance, in [3, 12]. We take an easy example of the Eikonal equation. { u = f(x) in Ω, (5.9) u = 0 on Ω, (5.10) where f 0 is a given bounded function in Ω. The DPP for the associated optimal control problem is { u ε (x) = inf y B ε(x) uε (y) + εf(x) for x Ω, u ε = 0 on Ω. We let U be the set of all functions Ω R again and T u(x) := { inf u(y) + εf(x) for x Ω, y B ε(x) 0 for x Ω. (5.11)
16 16 Q. LIU AND A. SCHIKORRA We only show that (A3) in Theorem 2.2 is satisfied since (A1), (A2) and (A4) are straightforward. We adopt the same partition of Ω as in (5.5) and take v(x) := 2k(sup f + 1) Ω We then have v(x) 1 = T v(x) + 1 if x Ω and v(x) T v(x) = 2 sup Ω when x Ω k for any k = 1, 2,..., N f + 2 sup f 2 Ω for every x Ω, which implies that v is a strict supersolution. By Theorem 3.2, we get u v in Ω for any subsolution u and therefore (A3) is verified. It follows from Theorem 2.2 that there exists a solution u = T u in Ω for the Eikonal operator (5.11). The solutions are unique since T satisfies (A5) in Theorem 4.1 after a translation of values on the boundary, i.e., we set { inf u(y) + εf(x) y B ˆT for x Ω u(x) := ε(x) C for x Ω, and discuss uniqueness of the solution û = u + C to û = ˆT û by following the argument in Example 5.1. It is possible to study the games for time-dependent problems as well. Consult [17] and [10] respectively for the games related to parabolic p-laplace equations and mean curvature flow equations. Our last example is from [11] about the deterministic games for a general parabolic equation. Example 5.4 (General fully nonlinear parabolic equations). Consider the parabolic equation { ut + F (x, u, 2 u) = 0 in R n (0, ), (5.12) u = g on R n {0}, (5.13) where F is assumed to fulfill the following: (F1) F satisfies the ellipticity condition. (F2) There exists C 1 > 0 such that F (x, p, Γ) F (x, p, Γ ) C 1 ( p p + Γ Γ ) for all x, p, p R n and symmetric matrices Γ, Γ. (F3) There exists C 2 > 0 and q, r 1 such that for all x, p and Γ. F (x, p, Γ) C 2 (1 + p q + Γ r )
17 EXISTENCE FOR DYNAMIC PROGRAMMING PRINCIPLE 17 The dynamic programming for the value function u ε in the associated games in this case is u ε (x, t) = g(x) for 0 t < ε 2 and u ε (x, t) = sup inf p, Γ w [ u ε (x + εw, t ε 2 ) εp w 1 ] 2 ε2 Γw w ε 2 F (x, p, Γ) for any t ε 2, x R n (0, ), where w, p R n with and Γ is an n n symmetric matrix satisfying w ε α, p ε β (5.14) Γ ε γ. (5.15) Here the parameters α, β and γ are all positive constants and satisfy proper relations below: α + β < 1, 2α + γ < 2, max(βq, γr) < 2; γ < 1 α, β(q 1) < α + 1, γ(r 1) < 2α, γr < 1 + α. (5.16) We let U be the set of all functions R n (0, ) R as before and define T u(x, t) := [ sup inf u(x + εw, t ε 2 ) εp w 1 ] p, Γ w 2 ε2 Γw w ε 2 F (x, p, Γ) if t ε 2, g(x) if t < ε 2. (5.17) Let us use Theorem 2.2 to get a solution of u = T u in X = R n [0, τ) for a fixed τ > 0 and a small ε > 0. Denote Ω = R n [ε 2, τ). The condition (A1) and (A4) are clearly fulfilled. It is also not difficult to see (A2) holds, since u 1 (x, t) = C F t + inf X\Ω g is a subsolution: for every x R n ( u1 (x + εw, t ε 2 ) ε 2 F (x, 0, O) ) T u 1 (x, t) inf w C F (t ε 2 ) + inf g C F ε 2 = u 1 (x, t) X\Ω when ε 2 t < τ and u 1 (x, t) g(x) when 0 t < ε 2. Concerning (A3), we construct a strict supersolution v in R n [0, τ). We first define an operator S for any x R n and any bounded function φ on R n : [ S[x, φ] := sup inf φ(x + εw) εp w 1 ] p, Γ w 2 ε2 Γw w ε 2 F (x, p, Γ)
18 18 Q. LIU AND A. SCHIKORRA with w, p, Γ satisfying (5.14) (5.15) and (5.16). Then we aim to find v(x, t) such that v(x, t) S[x, v(, t ε 2 )] + σ for t ε 2. with some σ > 0. We adapt [11, Lemma 4.1] to our setting, presented below for the reader s convenience. Lemma 5.5 (Lemma 4.1 in [11]). Assume that F satisfy (F1) (F3). Let S be defined as above with (5.14) (5.15) and (5.16). Then for any x R n and smooth bounded function φ, S[x, φ] φ = ε 2 F (x, φ, 2 φ) + o(ε 2 ). (5.18) The constant implied in the error term is uniform with respect to x. Now fix σ > 0. We show that v(x, t) := (σ + C 2 )t + sup X\Ω g + σε 2 is a strict supersolution when ε > 0 is small. By Lemma 5.5, for sufficiently small ε > 0, we then get S[x, v(, t ε 2 )] = v(x, t ε 2 ) + ε 2 F (x, 0, O) + o(ε 2 ) v(x, t ε 2 ) + C 2 ε 2 + o(ε 2 ) v(x, t) 1 2 σε2 when ε 2 t < τ and clearly v(x, t) g(x) + σε 2 when 0 t < ε 2. In conclusion, by Theorem 2.2, there exists a function u in R n [0, τ) satisfying u = T u, where T is as in (5.17). This solution coincides with the game value u ε in [11], since the uniqueness of solutions holds, as shown in the comparison principle below. Proposition 5.6 (Comparison principle for parabolic equations). If u 1 and u 2 are respectively a subsolution and a supersolution of u = T u in R n [0, τ), where T is defined in (5.17), then u 1 u 2 in R n [0, τ). Proof. We take h > 0 arbitrarily and define u h 2(x, t) := u 2 (x, t) + ht + h for (x, t) R n [0, τ). It is clear that u h 2 u uniformly as h 0. By direct calculation, we get u h 2(x, t) T u h 2(x, t) + hε 2 for ε 2 t < τ and u h 2(x, t) g(x) + δ for 0 t < ε 2. We are led to u 1 u h 2 in virtue of Theorem 3.2 and conclude by letting h 0. The last two examples of DPP below are related to p-laplacian, for which the associated dynamic programming principle involves integrals. It is natural to include measurability into consideration when choosing the function space U, but it is not known whether one may still obtain the closedness under
19 EXISTENCE FOR DYNAMIC PROGRAMMING PRINCIPLE 19 supremum or infimum, as in (A4) of Theorem 2.2. Extra work on a modification of the notion of extrema, compatible with the measurability, seems to be necessary. We leave it as a future topic. Example 5.7 (p-laplacian, I). We consider the normalized p-laplace equation: { u 2 p div( u p 2 u) = 0 in Ω, (5.19) u = g on Ω. (5.20) There are two dynamic programming principles known to generate the p- Laplace operator. One is based on the so-called asymptotic mean value property [18]: ( 1 2 sup y B ε(x) u ε (y) u ε (x) = α inf y B uε (y) + β u ε (y) dy in Ω, ε(x) B ε(x) u ε (x) = g(x) on Ω, (5.21) where α = p 2 2+n and β =. Here one would tend to take U to be the space p+n p+n of all measurable functions, and T be given by ( ) 1 α T u(x) = 2 sup u(y) + 1 inf y B ε(x) 2 u(y) + β u(y) dy for x Ω, y B ε(x) B ε(x) g(x) for x Ω. Nevertheless, it is not obvious how to show the closedness of T. For α > 0, the boundedness follows from the arguments in [16]. Example 5.8 (p-laplacian, II). The work by Peres and Sheffield [21] gave another game-theoretic approach to describe p-harmonic functions. The authors did not provide the DPP but a suggested version is u = T u with T defined below: T u(x) = ˆ 1 2 sup v B ε g(x) u(x + v + z) dµ v (z) + 1 ˆ 2 inf v B ε ) u(x + v + z)dµ v (z) for x Ω, for x Ω, where µ v is the uniform distribution on the sphere S v orthogonal to v with radius n 1 r = p 1. Here the verification of boundedness and closedness is a challenge.
20 20 Q. LIU AND A. SCHIKORRA References [1] T. Antunović, Y. Peres, S. Sheffield, and S. Somersille. Tug-of-war and infinity Laplace equation with vanishing Neumann boundary condition. Comm. Partial Differential Equations, 37(10): , [2] S. N. Armstrong and C. K. Smart. An easy proof of Jensen s theorem on the uniqueness of infinity harmonic functions. Calc. Var. Partial Differential Equations, 37(3-4): , [3] M. Bardi and I. Capuzzo-Dolcetta. Optimal control and viscosity solutions of Hamilton- Jacobi-Bellman equations. Systems & Control: Foundations & Applications. Birkhäuser Boston Inc., Boston, MA, With appendices by Maurizio Falcone and Pierpaolo Soravia. [4] G. Barles and P. E. Souganidis. Convergence of approximation schemes for fully nonlinear second order equations. Asymptotic Anal., 4(3): , [5] F. Charro, J. García Azorero, and J. D. Rossi. A mixed problem for the infinity Laplacian via tug-of-war games. Calc. Var. Partial Differential Equations, 34(3): , [6] M. G. Crandall, H. Ishii, and P.-L. Lions. User s guide to viscosity solutions of second order partial differential equations. Bull. Amer. Math. Soc. (N.S.), 27(1):1 67, [7] Y. Giga and Q. Liu. A billiard-based game interpretation of the Neumann problem for the curve shortening equation. Adv. Differential Equations, 14(3-4): , [8] H. Ishii. Perron s method for Hamilton-Jacobi equations. Duke Math. J., 55(2): , [9] H. Ishii. A simple, direct proof of uniqueness for solutions of the Hamilton-Jacobi equations of eikonal type. Proc. Amer. Math. Soc., 100(2): , [10] R. V. Kohn and S. Serfaty. A deterministic-control-based approach to motion by curvature. Comm. Pure Appl. Math., 59(3): , [11] R. V. Kohn and S. Serfaty. A deterministic-control-based approach to fully nonlinear parabolic and elliptic equations. Comm. Pure Appl. Math., 63(10): , [12] S. Koike. A beginner s guide to the theory of viscosity solutions, volume 13 of MSJ Memoirs. Mathematical Society of Japan, Tokyo, [13] E. Le Gruyer and J. C. Archer. Harmonious extensions. SIAM J. Math. Anal., 29(1): (electronic), [14] Q. Liu. Fattening and comparison principle for level-set equations of mean curvature type. SIAM J. Control Optim., 49(6): , [15] Q. Liu. Waiting time effect for motion by positive second derivatives and applications. preprint, [16] Q. Liu and A. Schikorra. A game-tree approach to discrete infinity Laplacian with running costs. preprint, [17] J. J. Manfredi, M. Parviainen, and J. D. Rossi. An asymptotic mean value characterization for a class of nonlinear parabolic equations related to tug-of-war games. SIAM J. Math. Anal., 42(5): , [18] J. J. Manfredi, M. Parviainen, and J. D. Rossi. An asymptotic mean value characterization for p-harmonic functions. Proc. Amer. Math. Soc., 138(3): , [19] Y. Peres, G. Pete, and S. Somersille. Biased tug-of-war, the biased infinity Laplacian, and comparison with exponential cones. Calc. Var. Partial Differential Equations, 38(3-4): , 2010.
21 EXISTENCE FOR DYNAMIC PROGRAMMING PRINCIPLE 21 [20] Y. Peres, O. Schramm, S. Sheffield, and D. B. Wilson. Tug-of-war and the infinity Laplacian. J. Amer. Math. Soc., 22(1): , [21] Y. Peres and S. Sheffield. Tug-of-war with noise: a game-theoretic view of the p- Laplacian. Duke Math. J., 145(1):91 120, [22] M. B. Rudd and H. A. Van Dyke. Median values, 1-harmonic functions, and functions of least gradient. Commun. Pure Appl. Anal., 12(2): , Qing Liu, Department of Mathematics, University of Pittsburgh, Pittsburgh, PA 15260, USA, qingliu@pitt.edu Armin Schikorra, Max-Planck Institut MiS Leipzig, Inselstr. 22, Leipzig, Germany, armin.schikorra@mis.mpg.de
AN ASYMPTOTIC MEAN VALUE CHARACTERIZATION FOR p-harmonic FUNCTIONS. To the memory of our friend and colleague Fuensanta Andreu
AN ASYMPTOTIC MEAN VALUE CHARACTERIZATION FOR p-harmonic FUNCTIONS JUAN J. MANFREDI, MIKKO PARVIAINEN, AND JULIO D. ROSSI Abstract. We characterize p-harmonic functions in terms of an asymptotic mean value
More informationAN ASYMPTOTIC MEAN VALUE CHARACTERIZATION FOR p-harmonic FUNCTIONS
AN ASYMPTOTIC MEAN VALUE CHARACTERIZATION FOR p-harmonic FUNCTIONS JUAN J. MANFREDI, MIKKO PARVIAINEN, AND JULIO D. ROSSI Abstract. We characterize p-harmonic functions in terms of an asymptotic mean value
More informationAN INTRODUCTION TO VISCOSITY SOLUTION THEORY. In this note, we study the general second-order fully nonlinear equations arising in various fields:
AN INTRODUCTION TO VISCOSITY SOLUTION THEORY QING LIU AND XIAODAN ZHOU 1. Introduction to Fully Nonlinear Equations In this note, we study the general second-order fully nonlinear equations arising in
More informationMEAN VALUE PROPERTY FOR p-harmonic FUNCTIONS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 40, Number 7, July 0, Pages 453 463 S 000-9939(0)8-X Article electronically published on November, 0 MEAN VALUE PROPERTY FOR p-harmonic FUNCTIONS
More informationOn the infinity Laplace operator
On the infinity Laplace operator Petri Juutinen Köln, July 2008 The infinity Laplace equation Gunnar Aronsson (1960 s): variational problems of the form S(u, Ω) = ess sup H (x, u(x), Du(x)). (1) x Ω The
More informationOn the definition and properties of p-harmonious functions
On the definition and properties of p-harmonious functions University of Pittsburgh, UBA, UAM Workshop on New Connections Between Differential and Random Turn Games, PDE s and Image Processing Pacific
More informationPERRON S METHOD FOR p-harmonious FUNCTIONS
Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 123, pp. 1 12. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu PERRON S METHOD FOR p-harmonious FUNCTIONS
More informationON THE DEFINITION AND PROPERTIES OF p-harmonious FUNCTIONS
ON THE DEFINITION AND PROPERTIES OF p-harmonious FUNCTIONS J. J. MANFREDI, M. PARVIAINEN, AND J. D. ROSSI Abstract. We consider functions that satisfy the identity } u ε(x) = α sup u ε + inf u ε + β u
More informationSOLUTIONS OF NONLINEAR PDES IN THE SENSE OF AVERAGES
SOLUTIONS OF NONLINEAR PDES IN THE SENSE OF AVERAGES BERND KAWOHL, JUAN MANFREDI, AND MIKKO PARVIAINEN Abstract. We characterize p-harmonic functions including p = 1 and p = by using mean value properties
More informationHausdorff Continuous Viscosity Solutions of Hamilton-Jacobi Equations
Hausdorff Continuous Viscosity Solutions of Hamilton-Jacobi Equations R Anguelov 1,2, S Markov 2,, F Minani 3 1 Department of Mathematics and Applied Mathematics, University of Pretoria 2 Institute of
More informationu(y) dy. In fact, as remarked in [MPR], we can relax this condition by requiring that it holds asymptotically
AN ASYMPTOTIC MEAN VALUE CHARACTERIZATION FOR A CLASS OF NONLINEAR PARABOLIC EQUATIONS RELATED TO TUG-OF-WAR GAMES JUAN J. MANFREDI, MIKKO PARVIAINEN, AND JULIO D. ROSSI Key words. Dirichlet boundary conditions,
More informationHomogenization and error estimates of free boundary velocities in periodic media
Homogenization and error estimates of free boundary velocities in periodic media Inwon C. Kim October 7, 2011 Abstract In this note I describe a recent result ([14]-[15]) on homogenization and error estimates
More informationExample 1. Hamilton-Jacobi equation. In particular, the eikonal equation. for some n( x) > 0 in Ω. Here 1 / 2
Oct. 1 0 Viscosity S olutions In this lecture we take a glimpse of the viscosity solution theory for linear and nonlinear PDEs. From our experience we know that even for linear equations, the existence
More informationEverywhere differentiability of infinity harmonic functions
Everywhere differentiability of infinity harmonic functions Lawrence C. Evans and Charles K. Smart Department of Mathematics University of California, Berkeley Abstract We show that an infinity harmonic
More informationSébastien Chaumont a a Institut Élie Cartan, Université Henri Poincaré Nancy I, B. P. 239, Vandoeuvre-lès-Nancy Cedex, France. 1.
A strong comparison result for viscosity solutions to Hamilton-Jacobi-Bellman equations with Dirichlet condition on a non-smooth boundary and application to parabolic problems Sébastien Chaumont a a Institut
More informationThe main motivation of this paper comes from the following, rather surprising, result of Ecker and Huisken [13]: for any initial data u 0 W 1,
Quasilinear parabolic equations, unbounded solutions and geometrical equations II. Uniqueness without growth conditions and applications to the mean curvature flow in IR 2 Guy Barles, Samuel Biton and
More informationCONTINUOUS DEPENDENCE ESTIMATES FOR VISCOSITY SOLUTIONS OF FULLY NONLINEAR DEGENERATE ELLIPTIC EQUATIONS
Electronic Journal of Differential Equations, Vol. 20022002), No. 39, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu login: ftp) CONTINUOUS DEPENDENCE
More informationBoundary value problems for the infinity Laplacian. regularity and geometric results
: regularity and geometric results based on joint works with Graziano Crasta, Roma La Sapienza Calculus of Variations and Its Applications - Lisboa, December 2015 on the occasion of Luísa Mascarenhas 65th
More informationVISCOSITY SOLUTIONS OF ELLIPTIC EQUATIONS
VISCOSITY SOLUTIONS OF ELLIPTIC EQUATIONS LUIS SILVESTRE These are the notes from the summer course given in the Second Chicago Summer School In Analysis, in June 2015. We introduce the notion of viscosity
More informationThe infinity-laplacian and its properties
U.U.D.M. Project Report 2017:40 Julia Landström Examensarbete i matematik, 15 hp Handledare: Kaj Nyström Examinator: Martin Herschend December 2017 Department of Mathematics Uppsala University Department
More informationTHE NEUMANN PROBLEM FOR THE -LAPLACIAN AND THE MONGE-KANTOROVICH MASS TRANSFER PROBLEM
THE NEUMANN PROBLEM FOR THE -LAPLACIAN AND THE MONGE-KANTOROVICH MASS TRANSFER PROBLEM J. GARCÍA-AZORERO, J. J. MANFREDI, I. PERAL AND J. D. ROSSI Abstract. We consider the natural Neumann boundary condition
More informationA GAME-THEORETIC PROOF OF CONVEXITY PRESERVING PROPERTIES FOR MOTION BY CURVATURE. 1. Introduction
A GAME-THEORETIC PROOF OF CONVEXITY PRESERVING PROPERTIES FOR MOTION BY CURVATURE QING LIU, ARMIN SCHIKORRA, AND XIAODAN ZHOU Abstract. In this paper, we revisit the convexity preserving properties for
More informationCOMPARISON PRINCIPLES FOR CONSTRAINED SUBHARMONICS PH.D. COURSE - SPRING 2019 UNIVERSITÀ DI MILANO
COMPARISON PRINCIPLES FOR CONSTRAINED SUBHARMONICS PH.D. COURSE - SPRING 2019 UNIVERSITÀ DI MILANO KEVIN R. PAYNE 1. Introduction Constant coefficient differential inequalities and inclusions, constraint
More informationTUG OF WAR INFINITY LAPLACIAN
TUG OF WAR and the INFINITY LAPLACIAN How to solve degenerate elliptic PDEs and the optimal Lipschitz extension problem by playing games. Yuval Peres, Oded Schramm, Scott Sheffield, and David Wilson Infinity
More informationA NOTE ON ZERO SETS OF FRACTIONAL SOBOLEV FUNCTIONS WITH NEGATIVE POWER OF INTEGRABILITY. 1. Introduction
A NOTE ON ZERO SETS OF FRACTIONAL SOBOLEV FUNCTIONS WITH NEGATIVE POWER OF INTEGRABILITY ARMIN SCHIKORRA Abstract. We extend a Poincaré-type inequality for functions with large zero-sets by Jiang and Lin
More informationSome notes on viscosity solutions
Some notes on viscosity solutions Jeff Calder October 11, 2018 1 2 Contents 1 Introduction 5 1.1 An example............................ 6 1.2 Motivation via dynamic programming............. 8 1.3 Motivation
More informationViscosity Solutions of the Bellman Equation for Perturbed Optimal Control Problems with Exit Times 0
Viscosity Solutions of the Bellman Equation for Perturbed Optimal Control Problems with Exit Times Michael Malisoff Department of Mathematics Louisiana State University Baton Rouge, LA 783-4918 USA malisoff@mathlsuedu
More informationRELATIONSHIP BETWEEN SOLUTIONS TO A QUASILINEAR ELLIPTIC EQUATION IN ORLICZ SPACES
Electronic Journal of Differential Equations, Vol. 2014 (2014), No. 265, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu RELATIONSHIP
More informationTHE OBSTACLE PROBLEM FOR THE p-laplacian VIA OPTIMAL STOPPING OF TUG-OF-WAR GAMES
THE OBSTACLE PROBLEM FOR THE p-laplacian VIA OPTIMAL STOPPING OF TUG-OF-WAR GAMES MARTA LEWICKA AND JUAN J. MANFREDI Abstract. We present a probabilistic approach to the obstacle problem for the p-laplace
More informationarxiv: v2 [math.ap] 10 Mar 2016
Hölder gradient estimates for parabolic homogeneous p-laplacian equations arxiv:1505.05525v2 [math.ap] 10 Mar 2016 Tianling Jin and Luis Silvestre March 11, 2016 Abstract We prove interior Hölder estimates
More informationNEUMANN BOUNDARY CONDITIONS FOR THE INFINITY LAPLACIAN AND THE MONGE-KANTOROVICH MASS TRANSPORT PROBLEM
NEUMANN BOUNDARY CONDITIONS FOR THE INFINITY LAPLACIAN AND THE MONGE-KANTOROVICH MASS TRANSPORT PROBLEM J. GARCIA-AZORERO, J. J. MANFREDI, I. PERAL AND J. D. ROSSI Abstract. In this note we review some
More informationMEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW
MEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW GREGORY DRUGAN AND XUAN HIEN NGUYEN Abstract. We present two initial graphs over the entire R n, n 2 for which the mean curvature flow
More informationOn the Bellman equation for control problems with exit times and unbounded cost functionals 1
On the Bellman equation for control problems with exit times and unbounded cost functionals 1 Michael Malisoff Department of Mathematics, Hill Center-Busch Campus Rutgers University, 11 Frelinghuysen Road
More informationHAMILTON-JACOBI EQUATIONS : APPROXIMATIONS, NUMERICAL ANALYSIS AND APPLICATIONS. CIME Courses-Cetraro August 29-September COURSES
HAMILTON-JACOBI EQUATIONS : APPROXIMATIONS, NUMERICAL ANALYSIS AND APPLICATIONS CIME Courses-Cetraro August 29-September 3 2011 COURSES (1) Models of mean field, Hamilton-Jacobi-Bellman Equations and numerical
More informationFlux-limited solutions for quasi-convex Hamilton-Jacobi equations on networks
Flux-limited solutions for quasi-convex Hamilton-Jacobi equations on networks C. Imbert and R. Monneau June 24, 2014 Abstract We study Hamilton-Jacobi equations on networks in the case where Hamiltonians
More informationA PDE Perspective of The Normalized Infinity Laplacian
A PDE Perspective of The ormalized Infinity Laplacian Guozhen Lu & Peiyong Wang Department of Mathematics Wayne State University 656 W.Kirby, 1150 FAB Detroit, MI 48202 E-mail: gzlu@math.wayne.edu & pywang@math.wayne.edu
More informationEQUAZIONI A DERIVATE PARZIALI. STEKLOV EIGENVALUES FOR THE -LAPLACIAN
EQUAZIONI A DERIVATE PARZIALI. STEKLOV EIGENVALUES FOR THE -LAPLACIAN J. GARCIA-AZORERO, J. J. MANFREDI, I. PERAL AND J. D. ROSSI Abstract. We study the Steklov eigenvalue problem for the - laplacian.
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationA MAXIMUM PRINCIPLE FOR SEMICONTINUOUS FUNCTIONS APPLICABLE TO INTEGRO-PARTIAL DIFFERENTIAL EQUATIONS
Dept. of Math. University of Oslo Pure Mathematics ISBN 82 553 1382 6 No. 18 ISSN 0806 2439 May 2003 A MAXIMUM PRINCIPLE FOR SEMICONTINUOUS FUNCTIONS APPLICABLE TO INTEGRO-PARTIAL DIFFERENTIAL EQUATIONS
More informationOptimal Control and Viscosity Solutions of Hamilton-Jacobi-Bellman Equations
Martino Bardi Italo Capuzzo-Dolcetta Optimal Control and Viscosity Solutions of Hamilton-Jacobi-Bellman Equations Birkhauser Boston Basel Berlin Contents Preface Basic notations xi xv Chapter I. Outline
More informationMINIMAL GRAPHS PART I: EXISTENCE OF LIPSCHITZ WEAK SOLUTIONS TO THE DIRICHLET PROBLEM WITH C 2 BOUNDARY DATA
MINIMAL GRAPHS PART I: EXISTENCE OF LIPSCHITZ WEAK SOLUTIONS TO THE DIRICHLET PROBLEM WITH C 2 BOUNDARY DATA SPENCER HUGHES In these notes we prove that for any given smooth function on the boundary of
More informationVISCOSITY SOLUTIONS. We follow Han and Lin, Elliptic Partial Differential Equations, 5.
VISCOSITY SOLUTIONS PETER HINTZ We follow Han and Lin, Elliptic Partial Differential Equations, 5. 1. Motivation Throughout, we will assume that Ω R n is a bounded and connected domain and that a ij C(Ω)
More informationLipschitz continuity for solutions of Hamilton-Jacobi equation with Ornstein-Uhlenbeck operator
Lipschitz continuity for solutions of Hamilton-Jacobi equation with Ornstein-Uhlenbeck operator Thi Tuyen Nguyen Ph.D student of University of Rennes 1 Joint work with: Prof. E. Chasseigne(University of
More informationAsymptotic behavior of infinity harmonic functions near an isolated singularity
Asymptotic behavior of infinity harmonic functions near an isolated singularity Ovidiu Savin, Changyou Wang, Yifeng Yu Abstract In this paper, we prove if n 2 x 0 is an isolated singularity of a nonegative
More informationON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS
Bendikov, A. and Saloff-Coste, L. Osaka J. Math. 4 (5), 677 7 ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS ALEXANDER BENDIKOV and LAURENT SALOFF-COSTE (Received March 4, 4)
More informationNumerical Solutions of Geometric Partial Differential Equations. Adam Oberman McGill University
Numerical Solutions of Geometric Partial Differential Equations Adam Oberman McGill University Sample Equations and Schemes Fully Nonlinear Pucci Equation! "#+ "#* "#) "#( "#$ "#' "#& "#% "#! "!!!"#$ "
More informationThuong Nguyen. SADCO Internal Review Metting
Asymptotic behavior of singularly perturbed control system: non-periodic setting Thuong Nguyen (Joint work with A. Siconolfi) SADCO Internal Review Metting Rome, Nov 10-12, 2014 Thuong Nguyen (Roma Sapienza)
More informationOn Generalized and Viscosity Solutions of Nonlinear Elliptic Equations
Advanced Nonlinear Studies 4 (2004), 289 306 On Generalized and Viscosity Solutions of Nonlinear Elliptic Equations David Hartenstine, Klaus Schmitt Department of Mathematics, University of Utah, 155 South
More informationTHE DIRICHLET PROBLEM FOR THE CONVEX ENVELOPE
THE DIRICHLET PROBLEM FOR THE CONVEX ENVELOPE LUIS SILVESTRE AND ADAM M. OBERMAN Abstract. This work studies the Dirichlet problem for the Convex Envelope. While the convex envelope is a natural object
More informationEquilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains
Equilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains J. Földes Department of Mathematics, Univerité Libre de Bruxelles 1050 Brussels, Belgium P. Poláčik School
More informationSome aspects of vanishing properties of solutions to nonlinear elliptic equations
RIMS Kôkyûroku, 2014, pp. 1 9 Some aspects of vanishing properties of solutions to nonlinear elliptic equations By Seppo Granlund and Niko Marola Abstract We discuss some aspects of vanishing properties
More informationLandesman-Lazer type results for second order Hamilton-Jacobi-Bellman equations
Author manuscript, published in "Journal of Functional Analysis 258, 12 (2010) 4154-4182" Landesman-Lazer type results for second order Hamilton-Jacobi-Bellman equations Patricio FELMER, Alexander QUAAS,
More informationNon-radial solutions to a bi-harmonic equation with negative exponent
Non-radial solutions to a bi-harmonic equation with negative exponent Ali Hyder Department of Mathematics, University of British Columbia, Vancouver BC V6TZ2, Canada ali.hyder@math.ubc.ca Juncheng Wei
More informationBrunn-Minkowski inequality for the 1-Riesz capacity and level set convexity for the 1/2-Laplacian
Brunn-Minkowski inequality for the 1-Riesz capacity and level set convexity for the 1/2-Laplacian M. Novaga, B. Ruffini January 13, 2014 Abstract We prove that that the 1-Riesz capacity satisfies a Brunn-Minkowski
More informationMATH Final Project Mean Curvature Flows
MATH 581 - Final Project Mean Curvature Flows Olivier Mercier April 30, 2012 1 Introduction The mean curvature flow is part of the bigger family of geometric flows, which are flows on a manifold associated
More informationRisk-Sensitive Control with HARA Utility
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 46, NO. 4, APRIL 2001 563 Risk-Sensitive Control with HARA Utility Andrew E. B. Lim Xun Yu Zhou, Senior Member, IEEE Abstract In this paper, a control methodology
More informationSome lecture notes for Math 6050E: PDEs, Fall 2016
Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.
More informationTitle. Author(s)Giga, Yoshikazu; Hamamuki, Nao. CitationHokkaido University Preprint Series in Mathematics, Issue Date DOI 10.
Title On a dynamic boundary condition for singular degener Author(s)Giga, Yoshikazu; Hamamuki, Nao CitationHokkaido University Preprint Series in Mathematics, Issue Date 2018-04-28 DOI 10.14943/84298 Doc
More informationEQUIVALENCE OF VISCOSITY AND WEAK SOLUTIONS FOR THE p(x)-laplacian
EQUIVALENCE OF VISCOSITY AND WEAK SOLUTIONS FOR THE p(x-laplacian PETRI JUUTINEN, TEEMU LUKKARI, AND MIKKO PARVIAINEN Abstract. We consider different notions of solutions to the p(x-laplace equation div(
More informationBoundary value problems for the infinity Laplacian. regularity and geometric results
: regularity and geometric results based on joint works with Graziano Crasta, Roma La Sapienza Calculus of variations, optimal transportation, and geometric measure theory: from theory to applications
More informationMotivation Power curvature flow Large exponent limit Analogues & applications. Qing Liu. Fukuoka University. Joint work with Prof.
On Large Exponent Behavior of Power Curvature Flow Arising in Image Processing Qing Liu Fukuoka University Joint work with Prof. Naoki Yamada Mathematics and Phenomena in Miyazaki 2017 University of Miyazaki
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationRégularité des équations de Hamilton-Jacobi du premier ordre et applications aux jeux à champ moyen
Régularité des équations de Hamilton-Jacobi du premier ordre et applications aux jeux à champ moyen Daniela Tonon en collaboration avec P. Cardaliaguet et A. Porretta CEREMADE, Université Paris-Dauphine,
More informationAsymptotic Behavior of Infinity Harmonic Functions Near an Isolated Singularity
Savin, O., and C. Wang. (2008) Asymptotic Behavior of Infinity Harmonic Functions, International Mathematics Research Notices, Vol. 2008, Article ID rnm163, 23 pages. doi:10.1093/imrn/rnm163 Asymptotic
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 informationHAIYUN ZHOU, RAVI P. AGARWAL, YEOL JE CHO, AND YONG SOO KIM
Georgian Mathematical Journal Volume 9 (2002), Number 3, 591 600 NONEXPANSIVE MAPPINGS AND ITERATIVE METHODS IN UNIFORMLY CONVEX BANACH SPACES HAIYUN ZHOU, RAVI P. AGARWAL, YEOL JE CHO, AND YONG SOO KIM
More informationViscosity Solutions for Dummies (including Economists)
Viscosity Solutions for Dummies (including Economists) Online Appendix to Income and Wealth Distribution in Macroeconomics: A Continuous-Time Approach written by Benjamin Moll August 13, 2017 1 Viscosity
More informationBrunn-Minkowski inequality for the 1-Riesz capacity and level set convexity for the 1/2-Laplacian
Brunn-Minkowski inequality for the 1-Riesz capacity and level set convexity for the 1/2-Laplacian M. Novaga, B. Ruffini Abstract We prove that that the 1-Riesz capacity satisfies a Brunn-Minkowski inequality,
More informationEXISTENCE RESULTS FOR QUASILINEAR HEMIVARIATIONAL INEQUALITIES AT RESONANCE. Leszek Gasiński
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS SUPPLEMENT 2007 pp. 409 418 EXISTENCE RESULTS FOR QUASILINEAR HEMIVARIATIONAL INEQUALITIES AT RESONANCE Leszek Gasiński Jagiellonian
More informationNECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES
NECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES JUHA LEHRBÄCK Abstract. We establish necessary conditions for domains Ω R n which admit the pointwise (p, β)-hardy inequality u(x) Cd Ω(x)
More informationA VERY BRIEF REVIEW OF MEASURE THEORY
A VERY BRIEF REVIEW OF MEASURE THEORY A brief philosophical discussion. Measure theory, as much as any branch of mathematics, is an area where it is important to be acquainted with the basic notions and
More informationSYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS. M. Grossi S. Kesavan F. Pacella M. Ramaswamy. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 12, 1998, 47 59 SYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS M. Grossi S. Kesavan F. Pacella M. Ramaswamy
More informationOn a mean value property related to the p-laplacian and p-harmonious functions
On a mean value property related to the p-laplacian and p-harmonious functions Ángel Arroyo García Universitat Autònoma de Barcelona Workshop on Complex Analysis and Operator Theory Málaga, June 23, 2016
More informationResearch Article Almost Periodic Viscosity Solutions of Nonlinear Parabolic Equations
Hindawi Publishing Corporation Boundary Value Problems Volume 29, Article ID 873526, 15 pages doi:1.1155/29/873526 Research Article Almost Periodic Viscosity Solutions of Nonlinear Parabolic Equations
More informationSymmetry of nonnegative solutions of elliptic equations via a result of Serrin
Symmetry of nonnegative solutions of elliptic equations via a result of Serrin P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract. We consider the Dirichlet problem
More informationMaths 212: Homework Solutions
Maths 212: Homework Solutions 1. The definition of A ensures that x π for all x A, so π is an upper bound of A. To show it is the least upper bound, suppose x < π and consider two cases. If x < 1, then
More informationConvexity of level sets for solutions to nonlinear elliptic problems in convex rings. Paola Cuoghi and Paolo Salani
Convexity of level sets for solutions to nonlinear elliptic problems in convex rings Paola Cuoghi and Paolo Salani Dip.to di Matematica U. Dini - Firenze - Italy 1 Let u be a solution of a Dirichlet problem
More informationA FINITE DIFFERENCE APPROACH TO THE INFINITY LAPLACE EQUATION AND TUG-OF-WAR GAMES
TRANSACTIONS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 364, Number 2, February 2012, Pages 595 636 S 0002-9947(2011)05289-X Article electronically published on September 14, 2011 A FINITE DIFFERENCE
More informationarxiv: v3 [math.ap] 28 Feb 2017
ON VISCOSITY AND WEA SOLUTIONS FOR NON-HOMOGENEOUS P-LAPLACE EQUATIONS arxiv:1610.09216v3 [math.ap] 28 Feb 2017 Abstract. In this manuscript, we study the relation between viscosity and weak solutions
More informationShape from Shading with discontinuous image brightness
Shape from Shading with discontinuous image brightness Fabio Camilli Dip. di Matematica Pura e Applicata, Università dell Aquila (Italy), camilli@ing.univaq.it Emmanuel Prados UCLA Vision Lab. (USA), eprados@cs.ucla.edu
More informationSOME PROPERTIES ON THE CLOSED SUBSETS IN BANACH SPACES
ARCHIVUM MATHEMATICUM (BRNO) Tomus 42 (2006), 167 174 SOME PROPERTIES ON THE CLOSED SUBSETS IN BANACH SPACES ABDELHAKIM MAADEN AND ABDELKADER STOUTI Abstract. It is shown that under natural assumptions,
More informationExistence of weak solutions for general nonlocal and nonlinear second-order parabolic equations
Existence of weak solutions for general nonlocal and nonlinear second-order parabolic equations Guy Barles, Pierre Cardaliaguet, Olivier Ley, Aurélien Monteillet To cite this version: Guy Barles, Pierre
More informationObstacle Problems Involving The Fractional Laplacian
Obstacle Problems Involving The Fractional Laplacian Donatella Danielli and Sandro Salsa January 27, 2017 1 Introduction Obstacle problems involving a fractional power of the Laplace operator appear in
More informationTHE TWO-PHASE MEMBRANE PROBLEM REGULARITY OF THE FREE BOUNDARIES IN HIGHER DIMENSIONS. 1. Introduction
THE TWO-PHASE MEMBRANE PROBLEM REGULARITY OF THE FREE BOUNDARIES IN HIGHER DIMENSIONS HENRIK SHAHGHOLIAN, NINA URALTSEVA, AND GEORG S. WEISS Abstract. For the two-phase membrane problem u = λ + χ {u>0}
More informationON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 29, 2004, 167 176 ON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES Piotr Haj lasz and Jani Onninen Warsaw University, Institute of Mathematics
More informationEXISTENCE AND UNIQUENESS OF p(x)-harmonic FUNCTIONS FOR BOUNDED AND UNBOUNDED p(x)
UNIVERSITY OF JYVÄSKYLÄ DEPARTMENT OF MATHEMATICS AND STATISTICS REPORT 30 UNIVERSITÄT JYVÄSKYLÄ INSTITUT FÜR MATHEMATIK UND STATISTIK BERICHT 30 EXISTENCE AND UNIQUENESS OF p(x)-harmonic FUNCTIONS FOR
More informationNon-trivial pinning threshold for an evolution equation involving long range interactions
Non-trivial pinning threshold for an evolution equation involving long range interactions Martin Jesenko joint work with Patrick Dondl Workshop: New trends in the variational modeling of failure phenomena
More information1. Introduction. The objective of this paper is to prove that the viscosity solutions of the p-laplace equation
SIAM J. MATH. ANAL. Vol. 33, No. 3, pp. 699 717 c 2001 Society for Industrial and Applied Mathematics ON THE EQUIVALENCE OF VISCOSITY SOLUTIONS AND WEAK SOLUTIONS FOR A QUASI-LINEAR EQUATION PETRI JUUTINEN,
More informationA CONNECTION BETWEEN A GENERAL CLASS OF SUPERPARABOLIC FUNCTIONS AND SUPERSOLUTIONS
A CONNECTION BETWEEN A GENERAL CLASS OF SUPERPARABOLIC FUNCTIONS AND SUPERSOLUTIONS RIIKKA KORTE, TUOMO KUUSI, AND MIKKO PARVIAINEN Abstract. We show to a general class of parabolic equations that every
More informationOn the infinity Laplacian and Hrushovski s fusion. Charles Krug Smart
On the infinity Laplacian and Hrushovski s fusion by Charles Krug Smart A dissertation submitted in partial satisfaction of the requirements for the degree of Doctor of Philosophy in Mathematics in the
More informationA Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations
A Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations Xiaobing Feng Department of Mathematics The University of Tennessee, Knoxville, U.S.A. Linz, November 23, 2016 Collaborators
More informationOn a Fractional Monge Ampère Operator
Ann. PDE (015) 1:4 DOI 10.1007/s40818-015-0005-x On a Fractional Monge Ampère Operator Luis Caffarelli 1 Fernando Charro Received: 16 November 015 / Accepted: 19 November 015 / Published online: 18 December
More informationRegularity estimates for fully non linear elliptic equations which are asymptotically convex
Regularity estimates for fully non linear elliptic equations which are asymptotically convex Luis Silvestre and Eduardo V. Teixeira Abstract In this paper we deliver improved C 1,α regularity estimates
More informationLevel-set convex Hamilton-Jacobi equations on networks
Level-set convex Hamilton-Jacobi equations on networks C. Imbert and R. Monneau January 17, 2014 Abstract The paper deals with Hamilton-Jacobi equations on networks with level-set convex (in the gradient
More informationA SIMPLE, DIRECT PROOF OF UNIQUENESS FOR SOLUTIONS OF THE HAMILTON-JACOBI EQUATIONS OF EIKONAL TYPE
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 100, Number 2, June 1987 A SIMPLE, DIRECT PROOF OF UNIQUENESS FOR SOLUTIONS OF THE HAMILTON-JACOBI EQUATIONS OF EIKONAL TYPE HITOSHI ISHII ABSTRACT.
More informationCompactness in Ginzburg-Landau energy by kinetic averaging
Compactness in Ginzburg-Landau energy by kinetic averaging Pierre-Emmanuel Jabin École Normale Supérieure de Paris AND Benoît Perthame École Normale Supérieure de Paris Abstract We consider a Ginzburg-Landau
More informationGLOBAL EXISTENCE RESULTS AND UNIQUENESS FOR DISLOCATION EQUATIONS
GLOBAL EXISTENCE RESULTS AND UNIQUENESS FOR DISLOCATION EQUATIONS GUY BARLES, PIERRE CARDALIAGUET, OLIVIER LEY & RÉGIS MONNEAU Abstract. We are interested in nonlocal Eikonal Equations arising in the study
More informationHartogs Theorem: separate analyticity implies joint Paul Garrett garrett/
(February 9, 25) Hartogs Theorem: separate analyticity implies joint Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ (The present proof of this old result roughly follows the proof
More informationOPTIMAL REGULARITY FOR A TWO-PHASE FREE BOUNDARY PROBLEM RULED BY THE INFINITY LAPLACIAN DAMIÃO J. ARAÚJO, EDUARDO V. TEIXEIRA AND JOSÉ MIGUEL URBANO
Pré-Publicações do Departamento de Matemática Universidade de Coimbra Preprint Number 18 55 OPTIMAL REGULARITY FOR A TWO-PHASE FREE BOUNDARY PROBLEM RULED BY THE INFINITY LAPLACIAN DAMIÃO J. ARAÚJO, EDUARDO
More information